摘要
Efficient numerical algorithm for stochastic differential equation has been an important object in the research of statistical physics and mathematics for a long time.In this work we study the highly accurate numerical algorithm for the overdamped Langevin equation.In particular,our interest is in the behaviour of the numerical schemes for solving the overdamped Langevin equation in the harmonic system.Based on the large friction limit of the underdamped Langevin dynamic scheme,three algorithms for overdamped Langevin equation are obtained.We derive the explicit expression of the stationary distribution of each algorithm by analysing the discrete time trajectory for both one-dimensional case and multi-dimensional case.The accuracy of the stationary distribution of each algorithm is illustrated by comparing with the exact Boltzmann distribution.Our results demonstrate that the“BAOA-limit”algorithm generates an accurate distribution of the harmonic system in a canonical ensemble,within a stable range of time interval.The other algorithms do not produce the exact distribution of the harmonic system.
作者
李德彰
杨小宝
De-Zhang Li;Xiao-Bao Yang(Department of Physics,South China University of Technology,Guangzhou 510640,China)
基金
Project supported by the Basic and Applied Basic Research Foundation of Guangdong Province,China(Grant No.2021A1515010328)
the Key-Area Research and Development Program of Guangdong Province,China(Grant No.2020B010183001)
the National Natural Science Foundation of China(Grant No.12074126)。
