Abstract:We present a numerical method to compute the approximation of the memory functions in the generalized Langevin models for collective dynamics of macromolecules. We first derive the exact expressions of the memory functions, obtained from projection to subspaces that correspond to the selection of coarse-grain variables. In particular, the memory functions are expressed in the forms of matrix functions, which will then be approximated by Krylov-subspace methods. It will also be demonstrated that the random noise can be approximated under the same framework, and the fluctuation-dissipation theorem is automatically satisfied. The accuracy of the method is examined through several numerical examples.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to efficiently calculate the memory function in the generalized Langevin model in the coarse - grained molecular model. Specifically, directly calculating these memory functions from molecular dynamics simulations requires a large amount of computational resources, and existing methods have difficulties in dealing with large - scale systems. Therefore, the authors propose a new algorithm based on the Krylov subspace method. This method can efficiently calculate the memory function without direct molecular dynamics simulations and can simultaneously approximate the stochastic noise term, thereby automatically satisfying the fluctuation - dissipation theorem. ### Main contributions of the paper: 1. **Proposing a new algorithm**: Developed a numerical algorithm based on the Krylov subspace method for calculating the memory function in the generalized Langevin model. 2. **Efficient calculation**: This method can efficiently calculate the memory function without direct molecular dynamics simulations, greatly reducing the computational cost. 3. **Stochastic noise approximation**: It can not only calculate the memory function, but also approximate the stochastic noise term within the same framework, ensuring the automatic satisfaction of the fluctuation - dissipation theorem. 4. **Verifying accuracy**: Verified the effectiveness and accuracy of this method through multiple numerical examples. ### Specific content of the paper: - **Introduction**: Introduced the basic principles and challenges of molecular dynamics simulations, especially the modeling difficulties for multi - scale biological processes. - **Theoretical derivation**: Derived in detail the expressions of the memory function and the stochastic noise term of the generalized Langevin model, using the Mori - Zwanzig projection form. - **Krylov subspace method**: Introduced how to use the Krylov subspace method to approximate the memory function and the stochastic noise term, solving the difficulties of direct calculation. - **Numerical experiments**: Through specific numerical examples, showed the application effects of this method in different coarse - grained subspaces and verified its effectiveness and accuracy. - **Discussion and outlook**: Discussed the advantages and potential application prospects of this method and pointed out the future research directions. ### Conclusion: The method proposed in this paper provides a new way for efficiently calculating the memory function in the generalized Langevin model, which not only reduces the computational cost but also ensures the accuracy and physical consistency of the model. This method has important application value in the modeling of multi - scale biological processes.

Computation of the Memory Functions in the Generalized Langevin Models for Collective Dynamics of Macromolecules

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

Memory corrections to Markovian Langevin dynamics

Explicit analytical form for memory kernel in the generalized Langevin equation for end-to-end vector of Rouse chains.

Accurate Memory Kernel Extraction from Discretized Time Series Data

Accurate estimates of dynamical statistics using memory

Numerical Integration of the Extended Variable Generalized Langevin Equation with a Positive Prony Representable Memory Kernel

Effective mass approach to memory in non-Markovian systems

Construction of coarse-grained molecular dynamics with many-body non-Markovian memory

On the generalization ability of coarse-grained molecular dynamics models for non-equilibrium processes

Learning Memory Kernels in Generalized Langevin Equations

Accurate Memory Kernel Extraction from Discretized Time-Series Data

Data-driven learning of the generalized Langevin equation with state-dependent memory

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

The Generalized Langevin Equation with a power-law memory in a nonlinear potential well

Accurate and robust splitting methods for the generalized Langevin equation with a positive Prony series memory kernel

Computation of stationary distributions in stochastic models of cellular processes with molecular memory

The Interplay between Memory and Potentials of Mean Force: A Discussion on the Structure of Equations of Motion for Coarse Grained Observables