Wang Shu, Ma Zhan, Pan Wenxiao
Department of Mechanical Engineering, University of Wisconsin-Madison, Madison, WI 53706, USA.
Soft Matter. 2021 Jul 7;17(26):6404-6412. doi: 10.1039/d1sm00413a.
Modeling a high-dimensional Hamiltonian system in reduced dimensions with respect to coarse-grained (CG) variables can greatly reduce computational cost and enable efficient bottom-up prediction of main features of the system for many applications. However, it usually experiences significantly altered dynamics due to loss of degrees of freedom upon coarse-graining. To establish CG models that can faithfully preserve dynamics, previous efforts mainly focused on equilibrium systems. In contrast, various soft matter systems are known to be out of equilibrium. Therefore, the present work concerns non-equilibrium systems and enables accurate and efficient CG modeling that preserves non-equilibrium dynamics and is generally applicable to any non-equilibrium process and any observable of interest. To this end, the dynamic equation of a CG variable is built in the form of the non-stationary generalized Langevin equation (nsGLE), where the two-time memory kernel is determined from the data of the auto-correlation function of the observable of interest. By embedding the nsGLE in an extended dynamics framework, the nsGLE can be solved efficiently to predict the non-equilibrium dynamics of the CG variable. To prove and exploit the equivalence of the nsGLE and extended dynamics, the memory kernel is parameterized in a two-time exponential expansion. A data-driven hybrid optimization process is proposed for the parameterization, which integrates the differential-evolution method with the Levenberg-Marquardt algorithm to efficiently tackle a non-convex and high-dimensional optimization problem.
针对粗粒化(CG)变量在低维空间中对高维哈密顿系统进行建模,可以大幅降低计算成本,并能够对许多应用中系统的主要特征进行高效的自下而上预测。然而,由于粗粒化过程中自由度的损失,该模型通常会经历显著改变的动力学。为了建立能够忠实保留动力学的CG模型,以往的工作主要集中在平衡系统上。相比之下,各种软物质系统已知处于非平衡状态。因此,本工作关注非平衡系统,并实现了准确且高效的CG建模,该建模保留了非平衡动力学,并且普遍适用于任何非平衡过程和任何感兴趣的可观测量。为此,以非平稳广义朗之万方程(nsGLE)的形式构建CG变量的动力学方程,其中双时记忆核由感兴趣可观测量的自相关函数数据确定。通过将nsGLE嵌入扩展动力学框架,可以有效地求解nsGLE以预测CG变量的非平衡动力学。为了证明并利用nsGLE与扩展动力学的等价性,记忆核采用双时指数展开进行参数化。针对参数化提出了一种数据驱动的混合优化过程,该过程将差分进化方法与列文伯格 - 马夸尔特算法相结合,以有效解决非凸和高维优化问题。