Comparison of effective and stable Langevin dynamics integrators
arXiv:2304.09030 · doi:10.1016/j.cpc.2024.109152
Abstract
Langevin and Brownian simulations play a prominent role in computational research, and state of the art integration algorithms provide trajectories with different stability ranges and accuracy in reproducing statistical averages. The practical usability of integrators is an important aspect to allow choosing large time steps while ensuring numerical stability and overall computational efficiency. In this work, different use cases and practical features are selected in order to perform a cumulative comparison of integrators with a focus on evaluating the derived velocity and position autocorrelation functions, a comparison that is often disregarded in the literature. A standard industrial open-source software methodology is suggested to compare systematically the different algorithms.
34 pages, 12 figures
References in corpus (9)
- A simple and effective Verlet-type algorithm for simulating Langevin dynamics
- Robust and efficient configurational molecular sampling via Langevin Dynamics
- Design of quasi-symplectic propagators for Langevin dynamics
- Multiscale coupling of molecular dynamics and hydrodynamics: application to DNA translocation through a nanopore
- Rayleigh's dissipation function at work
- Comparison of structure and transport properties of concentrated hard and soft sphere fluids
- Mesoscopic simulations at the physics-chemistry-biology interface
- Fourth Order Algorithms for Solving the Multivariable Langevin Equation and the Kramers Equation
- Trotter Derivation of Algorithms for Brownian and Dissipative Particle Dynamics