Fourth Order Algorithms for Solving the Multivariable Langevin Equation and the Kramers Equation

Abstract

We develop a fourth order simulation algorithm for solving the stochastic Langevin equation. The method consists of identifying solvable operators in the Fokker-Planck equation, factorizing the evolution operator for small time steps to fourth order and implementing the factorization process numerically. A key contribution of this work is to show how certain double commutators in the factorization process can be simulated in practice. The method is general, applicable to the multivariable case, and systematic, with known procedures for doing fourth order factorizations. The fourth order convergence of the resulting algorithm allowed very large time steps to be used. In simulating the Brownian dynamics of 121 Yukawa particles in two dimensions, the converged result of a first order algorithm can be obtained by using time steps 50 times as large. To further demostrate the versatility of our method, we derive two new classes of fourth order algorithms for solving the simpler Kramers equation without requiring the derivative of the force. The convergence of many fourth order algorithms for solving this equation are compared.

19 pages, 2 figures

Cited by in corpus (17)

• Accurate sampling using Langevin dynamics
• Lattice Boltzmann simulations of soft matter systems
• Higher Order Decompositions of Ordered Operator Exponentials
• Design of quasi-symplectic propagators for Langevin dynamics
• On the construction of high-order force gradient algorithms for integration of motion in classical and quantum systems
• Extrapolated High-Order Propagators for Path Integral Monte Carlo Simulations
• Fourth-Order Algorithms for Solving the Imaginary Time Gross-Pitaevskii Equation in a Rotating Anisotropic Trap
• Quasi symplectic integrators for stochastic differential equations
• Quantum Statistical Calculations and Symplectic Corrector Algorithms
• Efficient numerical integrators for stochastic models
• High-order Path Integral Monte Carlo methods for solving quantum dot problems
• Event Driven Langevin simulations of Hard Spheres
• Higher harmonics of the magnetoplasmon in strongly coupled Coulomb and Yukawa systems
• Forward Symplectic Integrators and the Long Time Phase Error in Periodic Motions
• The Complete Characterization of Fourth-Order Symplectic Integrators with Extended-Linear Coefficients
• Comparison of effective and stable Langevin dynamics integrators
• Forward and non-forward symplectic integrators in solving classical dynamics problems