期刊名称：Proceedings of the National Academy of Sciences
印刷版ISSN：0027-8424
电子版ISSN：1091-6490
出版年度：2016
卷号：113
期号：50
页码：14183-14188
DOI：10.1073/pnas.1609587113
语种：English
出版社：The National Academy of Sciences of the United States of America
摘要：SignificanceThe generalized Langevin equation (GLE) provides a precise description of coarse-grained variable dynamics in reduced dimension models. However, computation of the memory kernel poses a major challenge to the practical use of the GLE. This paper presents a data-driven approach to compute the memory kernel, relying on a hierarchy of parameterized rational approximations in terms of the Laplace transform, which can be expanded to arbitrarily high order as necessary. This parameterization makes it convenient to represent the GLE via an extended stochastic model where the memory term is eliminated by properly introducing auxiliary variables. The present method is well-suited for constructing reduced models for nonequilibrium properties of complex systems such as biomolecules, chemical reaction networks, and climate simulations. We present a data-driven approach to determine the memory kernel and random noise in generalized Langevin equations. To facilitate practical implementations, we parameterize the kernel function in the Laplace domain by a rational function, with coefficients directly linked to the equilibrium statistics of the coarse-grain variables. We show that such an approximation can be constructed to arbitrarily high order and the resulting generalized Langevin dynamics can be embedded in an extended stochastic model without explicit memory. We demonstrate how to introduce the stochastic noise so that the second fluctuation-dissipation theorem is exactly satisfied. Results from several numerical tests are presented to demonstrate the effectiveness of the proposed method.
