A Gauss-Newton method for iterative optimization of memory kernels for generalized Langevin thermostats in coarse-grained molecular dynamics simulations

V Klippenstein and N Wolf and NFA van der Vegt, JOURNAL OF CHEMICAL PHYSICS, 160, 204115 (2024).

DOI: 10.1063/5.0203832

In molecular dynamics simulations, dynamically consistent coarse-grained (CG) models commonly use stochastic thermostats to model friction and fluctuations that are lost in a CG description. While Markovian, i.e., time-local, formulations of such thermostats allow for an accurate representation of diffusivities/long-time dynamics, a correct description of the dynamics on all time scales generally requires non- Markovian, i.e., non-time-local, thermostats. These thermostats typically take the form of a Generalized Langevin Equation (GLE) determined by a memory kernel. In this work, we use a Markovian embedded formulation of a position-independent GLE thermostat acting independently on each CG degree of freedom. Extracting the memory kernel of this CG model from atomistic reference data requires several approximations. Therefore, this task is best understood as an inverse problem. While our recently proposed approximate Newton scheme allows for the iterative optimization of memory kernels (IOMK), Markovian embedding remained potentially error-prone and computationally expensive. In this work, we present an IOMK-Gauss-Newton scheme (IOMK- GN) based on IOMK that allows for the direct parameterization of a Markovian embedded model. (c) 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/)

Return to Publications page