.. math:: k = df/dx
and a pulling velocity
-
+unfolindg
.. math:: v = dx/dt
we have a loading rate
def bell_rate(sawsim_runner, num_domains=1, unfolding_rate=1,
unfolding_distance=1, temperature=1/kB, spring_constant=1,
- velocity=1, N=100):
+ velocity=1, N=200):
loading_rate = float(spring_constant * velocity)
rho = kB * temperature / unfolding_distance
alpha = rho * log(loading_rate / (unfolding_rate * rho))