Abstract

We compare two numerical algorithms for the computer modeling of the thermal decay process (Kramers problem) solving the stochastic (Langevin) equations for the generalized coordinate and its conjugate momentum. These are known in the literature (i) the ALGO algorithm including in the fluctuation part the terms up to τ3/2 (t is the time step of the numerical modeling) and (ii) the Euler-Maruyama algorithm including in the fluctuation part only the terms proportional to τ1/2. We concentrate on the quasistationary decay rate and transient time. The ALGO algorithm appears to be more efficient, however, the optimal size of the time step has a strong and non-trivial dependence upon the dissipation strength and governing parameter.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.