Model reduction techniques for the computation of extended Markov parameterizations for generalized Langevin equations

Abstract

The generalized Langevin equation is a model for the motion of coarse-grained particles where dissipative forces are represented by a memory term. The numerical realization of such a model requires the implementation of a stochastic delay-differential equation and the estimation of a corresponding memory kernel. Here we develop a new approach for computing a data-driven Markov model for the motion of the particles, given equidistant samples of their velocity autocorrelation function. Our method bypasses the determination of the underlying memory kernel by representing it via up to about twenty auxiliary variables. The algorithm is based on a sophisticated variant of the Prony method for exponential interpolation and employs the Positive Real Lemma from model reduction theory to extract the associated Markov model. We demonstrate the potential of this approach for the test case of anomalous diffusion, where data are given analytically, and then apply our method to velocity autocorrelation data of molecular dynamics simulations of a colloid in a Lennard-Jones fluid. In both cases, the VACF and the memory kernel can be reproduced very accurately. Moreover, we show that the algorithm can also handle input data with large statistical noise. We anticipate that it will be a very useful tool in future studies that involve dynamic coarse-graining of complex soft matter systems.

References in corpus (16)

• Anomalous transport in the crowded world of biological cells
• Fractional Langevin equation
• Nuclear quantum effects in solids using a colored-noise thermostat
• Colored-noise thermostats à la carte
• Langevin equation with colored noise for constant-temperature molecular dynamics simulations
• Data-driven parameterization of the generalized Langevin equation
• Iterative Reconstruction of Memory Kernels
• Brownian motion from molecular dynamics
• 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
• On the "generalized Generalized Langevin Equation"
• Generalized Langevin dynamics: Construction and numerical integration of non-Markovian particle-based models
• The Derivation and Approximation of Coarse-grained Dynamics from Langevin Dynamics
• Computation of the Memory Functions in the Generalized Langevin Models for Collective Dynamics of Macromolecules
• From Generalized Langevin Equations to Brownian Dynamics and Embedded Brownian Dynamics
• Uncertainty quantification for generalized Langevin dynamics

Cited by in corpus (5)

• Coarse-Grained Modelling Out of Equilibrium
• Likelihood-based non-Markovian models from molecular dynamics
• Understanding Dynamics in Coarse-Grained Models: I. Universal Excess Entropy Scaling Relationship
• Understanding Dynamics in Coarse-Grained Models: II. Coarse-Grained Diffusion Modeled Using Hard Sphere Theory
• A Gauss-Newton method for iterative optimization of memory kernels for generalized Langevin thermostats in coarse-grained molecular dynamics simulations