Convergence of a Numerical Method for Solving Discontinuous Fokker-Planck Equations

In studies of molecular motors, the stochastic motion is modeled using the Langevin equation. If we consider an ensemble of motors, the probability density is governed by the corresponding Fokker–Planck equation. Average quantities, such as average velocity, effective diffusion coefficient, and randomness parameter, can be calculated from the probability density. A numerical method was previously developed to solve Fokker–Planck equations [H. Wang, C. Peskin, and T. Elston, J. Theoret. Biol., 221 (2003), pp. 491–511]. It preserves detailed balance, which ensures that if the system is forced to an equilibrium, the numerical solution will be the same as the Boltzmann distribution. Here we study the convergence of this numerical method when the potential has a finite number of discontinuities at half numerical grid points. We prove that this numerical method is stable and is consistent with the differential equation. Based on the consistency analysis, we propose a modified version of this numerical method to...

[1]  George Oster,et al.  Energy transduction in the F1 motor of ATP synthase , 1998, Nature.

[2]  F. Reif,et al.  Fundamentals of Statistical and Thermal Physics , 1965 .

[3]  Hong Wang,et al.  Mathematical theory of molecular motors and a new approach for uncovering motor mechanism. , 2003, IEE proceedings. Nanobiotechnology.

[4]  H. Berg.  Random Walks in Biology , 2018 .

[5]  R. Cox,et al.  Steady-state currents in sharp stochastic ratchets. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[6]  T. Elston,et al.  Numerical and analytical studies of nonequilibrium fluctuation-induced transport processes , 1996 .

[7]  R. Vale,et al.  The load dependence of kinesin's mechanical cycle. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[8]  G. D. Maso,et al.  Definition and weak stability of nonconservative products , 1995 .

[9]  R. Leighton,et al.  The Feynman Lectures on Physics; Vol. I , 1965 .

[10]  Charles S. Peskin,et al.  The Role of Protein Flexibility in Molecular Motor Function: Coupled Diffusion in a Tilted Periodic Potential , 2000, SIAM J. Appl. Math..

[11]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[12]  Lubomira A. Dontcheva,et al.  Constructive role of noise: Fast fluctuation asymptotics of transport in stochastic ratchets. , 1998, Chaos.

[13]  小寺 武康,et al.  On the theory of the Brownian motion , 1959 .

[14]  George Oster,et al.  Energy transduction in ATP synthase , 1998, Nature.

[15]  Prost,et al.  Asymmetric pumping of particles. , 1994, Physical review letters.

[16]  P. Raviart,et al.  Numerical Approximation of Hyperbolic Systems of Conservation Laws , 1996, Applied Mathematical Sciences.

[17]  Jan Pieter Abrahams,et al.  Structure at 2.8 Â resolution of F1-ATPase from bovine heart mitochondria , 1994, Nature.

[18]  T. Elston,et al.  A robust numerical algorithm for studying biomolecular transport processes. , 2003, Journal of theoretical biology.

[19]  C S Peskin,et al.  Cellular motions and thermal fluctuations: the Brownian ratchet. , 1993, Biophysical journal.

[20]  Kazuhiko Kinosita,et al.  Direct observation of the rotation of F1-ATPase , 1997, Nature.

[21]  R. Astumian.  Thermodynamics and kinetics of a Brownian motor. , 1997, Science.

[22]  David M. Young,et al.  ON THE CRANK-NICOLSON PROCEDURE FOR SOLVING PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS, , 1957 .

[23]  Mark J. Schnitzer,et al.  Single kinesin molecules studied with a molecular force clamp , 1999, Nature.