Accurate and robust splitting methods for the generalized Langevin equation with a positive Prony series memory kernel
arXiv:2109.07879 · doi:10.1016/j.jcp.2022.111332
Abstract
We study numerical methods for the generalized Langevin equation (GLE) with a positive Prony series memory kernel, in which case the GLE can be written in an extended variable Markovian formalism. We propose a new splitting method that is easy to implement and is able to substantially improve the accuracy and robustness of GLE simulations in a wide range of the parameters. An error analysis is performed in the case of a one-dimensional harmonic oscillator, revealing that several averages are exact for the newly proposed method. Various numerical experiments in both equilibrium and nonequilibrium simulations are also conducted to compare the method with popular alternatives in interacting multi-particle systems.
References in corpus (24)
- Nuclear quantum effects in solids using a colored-noise thermostat
- Colored-noise thermostats à la carte
- Robust and efficient configurational molecular sampling via Langevin Dynamics
- Data-driven parameterization of the generalized Langevin equation
- The computation of averages from equilibrium and nonequilibrium Langevin molecular dynamics
- On the numerical treatment of dissipative particle dynamics and related systems
- Asymptotic analysis for the generalized langevin equation
- The Generalized Langevin Equation: An efficient approach to non-equilibrium molecular dynamics of open systems
- Numerical Integration of the Extended Variable Generalized Langevin Equation with a Positive Prony Representable Memory Kernel
- GENERIC formalism of a Vlasov-Fokker-Planck equation and connection to large-deviation principles
- Generalized Langevin dynamics: Construction and numerical integration of non-Markovian particle-based models
- Data-driven molecular modeling with the generalized Langevin equation
- Efficient numerical integrators for stochastic models
- High-Order Langevin Diffusion Yields an Accelerated MCMC Algorithm
- Trotter Derivation of Algorithms for Brownian and Dissipative Particle Dynamics
- Adaptive Thermostats for Noisy Gradient Systems
- Fluctuation symmetry leads to GENERIC equations with non-quadratic dissipation
- Structure-preserving integrators for dissipative systems based on reversible-irreversible splitting
- Pairwise adaptive thermostats for improved accuracy and stability in dissipative particle dynamics
- Assessing numerical methods for molecular and particle simulation
- Uncertainty quantification for generalized Langevin dynamics
- Time correlation functions of equilibrium and nonequilibrium Langevin dynamics: Derivations and numerics using random numbers
- On the Generalised Langevin Equation for Simulated Annealing
- Accurate and efficient splitting methods for dissipative particle dynamics