Numerical Integration of the Extended Variable Generalized Langevin Equation with a Positive Prony Representable Memory Kernel
arXiv:1304.6446 · doi:10.1063/1.4815917
Abstract
Generalized Langevin dynamics (GLD) arise in the modeling of a number of systems, ranging from structured fluids that exhibit a viscoelastic mechanical response, to biological systems, and other media that exhibit anomalous diffusive phenomena. Molecular dynamics (MD) simulations that include GLD in conjunction with external and/or pairwise forces require the development of numerical integrators that are efficient, stable, and have known convergence properties. In this article, we derive a family of extended variable integrators for the Generalized Langevin equation (GLE) with a positive Prony series memory kernel. Using stability and error analysis, we identify a superlative choice of parameters and implement the corresponding numerical algorithm in the LAMMPS MD software package. Salient features of the algorithm include exact conservation of the first and second moments of the equilibrium velocity distribution in some important cases, stable behavior in the limit of conventional Langevin dynamics, and the use of a convolution-free formalism that obviates the need for explicit storage of the time history of particle velocities. Capability is demonstrated with respect to accuracy in numerous canonical examples, stability in certain limits, and an exemplary application in which the effect of a harmonic confining potential is mapped onto a memory kernel.
23 pages, 8 figures
References in corpus (3)
Cited by in corpus (27)
- Data-driven parameterization of the generalized Langevin equation
- Quantum Dynamics of Vibrational Polariton Chemistry
- Active mechanics reveal molecular-scale force kinetics in living oocytes
- Coarse-Grained Modelling Out of Equilibrium
- The Generalized Langevin Equation: An efficient approach to non-equilibrium molecular dynamics of open systems
- Generalized Langevin dynamics: Construction and numerical integration of non-Markovian particle-based models
- Likelihood-based non-Markovian models from molecular dynamics
- The Derivation and Approximation of Coarse-grained Dynamics from Langevin Dynamics
- Applications of the Generalised Langevin Equation: towards a realistic description of the baths
- Parametrizing linear generalized Langevin dynamics from explicit molecular dynamics simulations
- Long-time persistence of hydrodynamic memory boosts microparticle transport
- Construction of coarse-grained molecular dynamics with many-body non-Markovian memory
- How good is the generalized Langevin equation to describe the dynamics of photo-induced electron transfer in fluid solution?
- Model reduction techniques for the computation of extended Markov parameterizations for generalized Langevin equations
- A large deviation theory perspective on nanoscale transport phenomena
- Data-driven construction of stochastic reduced dynamics encoded with non-Markovian features
- Coarse-graining Langevin dynamics using reduced-order techniques
- Nonequilibrium processes from Generalised Langevin Equations: realistic nanoscale systems connected to two thermal baths
- Machine learning memory kernels as closure for non-Markovian stochastic processes
- Accurate and robust splitting methods for the generalized Langevin equation with a positive Prony series memory kernel
- Surmounting potential barriers: hydrodynamic memory hedges against thermal fluctuations in particle transport
- Uncertainty quantification for generalized Langevin dynamics
- Non-Markovian systems out of equilibrium: Exact results for two routes of coarse graining
- Internal noise driven generalized Langevin equation from a nonlocal continuum model
- On the Asymptotic Behavior of the Kernel Function in the Generalized Langevin Equation: A One-dimensional lattice model
- Non-Markovianity, entropy production, and Jarzynski equality
- GROMACS Stochastic Dynamics and BAOAB are equivalent configurational sampling algorithms