Physics – Nuclear Physics – Nuclear Theory
Scientific paper
2000-06-30
Phys.Rev. E63 (2001) 016703
Physics
Nuclear Physics
Nuclear Theory
19 pages, 2 figures
Scientific paper
10.1103/PhysRevE.63.016703
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.
Chin Siu A.
Forbert Harald A.
No associations
LandOfFree
Fourth Order Algorithms for Solving the Multivariable Langevin Equation and the Kramers Equation does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Fourth Order Algorithms for Solving the Multivariable Langevin Equation and the Kramers Equation, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Fourth Order Algorithms for Solving the Multivariable Langevin Equation and the Kramers Equation will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-528695