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

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Klippenstein, V., Wolf, N., van der Vegt, N. F. A.
Format: Preprint
Veröffentlicht: 2024
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866917699936845824
author Klippenstein, V.
Wolf, N.
van der Vegt, N. F. A.
author_facet Klippenstein, V.
Wolf, N.
van der Vegt, N. F. A.
contents 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 are typically in 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, Iterative Optimization of memory kernels (IOMK), allows for the iterative optimization of a memory kernel, Markovian embedding remained potentially error-prone and computationally expensive. In this work, we present a IOMK-Gauss-Newton scheme (IOMK-GN) based on IOMK, that allows for the direct parameterization of a Markovian embedded model.
format Preprint
id arxiv_https___arxiv_org_abs_2402_10652
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle A Gauss-Newton method for iterative optimization of memory kernels for generalized Langevin thermostats in coarse-grained molecular dynamics simulations
Klippenstein, V.
Wolf, N.
van der Vegt, N. F. A.
Statistical Mechanics
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 are typically in 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, Iterative Optimization of memory kernels (IOMK), allows for the iterative optimization of a memory kernel, Markovian embedding remained potentially error-prone and computationally expensive. In this work, we present a IOMK-Gauss-Newton scheme (IOMK-GN) based on IOMK, that allows for the direct parameterization of a Markovian embedded model.
title A Gauss-Newton method for iterative optimization of memory kernels for generalized Langevin thermostats in coarse-grained molecular dynamics simulations
topic Statistical Mechanics
url https://arxiv.org/abs/2402.10652