Department of Chemistry, James Frank Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA.
J Chem Phys. 2018 Mar 14;148(10):102335. doi: 10.1063/1.5010270.
Coarse-grained (CG) models serve as a powerful tool to simulate molecular systems at much longer temporal and spatial scales. Previously, CG models and methods have been built upon classical statistical mechanics. The present paper develops a theory and numerical methodology for coarse-graining in quantum statistical mechanics, by generalizing the multiscale coarse-graining (MS-CG) method to quantum Boltzmann statistics. A rigorous derivation of the sufficient thermodynamic consistency condition is first presented via imaginary time Feynman path integrals. It identifies the optimal choice of CG action functional and effective quantum CG (qCG) force field to generate a quantum MS-CG (qMS-CG) description of the equilibrium system that is consistent with the quantum fine-grained model projected onto the CG variables. A variational principle then provides a class of algorithms for optimally approximating the qMS-CG force fields. Specifically, a variational method based on force matching, which was also adopted in the classical MS-CG theory, is generalized to quantum Boltzmann statistics. The qMS-CG numerical algorithms and practical issues in implementing this variational minimization procedure are also discussed. Then, two numerical examples are presented to demonstrate the method. Finally, as an alternative strategy, a quasi-classical approximation for the thermal density matrix expressed in the CG variables is derived. This approach provides an interesting physical picture for coarse-graining in quantum Boltzmann statistical mechanics in which the consistency with the quantum particle delocalization is obviously manifest, and it opens up an avenue for using path integral centroid-based effective classical force fields in a coarse-graining methodology.
粗粒化 (CG) 模型是一种强大的工具,可用于在更长的时间和空间尺度上模拟分子系统。以前,CG 模型和方法是基于经典统计力学构建的。本文通过将多尺度粗粒化 (MS-CG) 方法推广到量子玻尔兹曼统计,为量子统计力学中的粗粒化建立了理论和数值方法。首先通过虚时间费曼路径积分,给出了充分热力学一致性条件的严格推导。它确定了 CG 作用泛函和有效量子 CG(qCG) 力场的最佳选择,以生成与 CG 变量上投影的量子细粒模型一致的平衡系统的量子 MS-CG(qMS-CG) 描述。然后变分原理提供了一类最优逼近 qMS-CG 力场的算法。具体来说,基于力匹配的变分方法,也被应用于经典 MS-CG 理论,被推广到量子玻尔兹曼统计。还讨论了 qMS-CG 数值算法和实施这种变分最小化过程的实际问题。然后,给出了两个数值示例来说明该方法。最后,作为一种替代策略,导出了 CG 变量中热密度矩阵的准经典近似。这种方法为量子玻尔兹曼统计力学中的粗粒化提供了一个有趣的物理图景,其中与量子粒子离域的一致性显然显而易见,它为在粗粒化方法中使用基于路径积分质心的有效经典力场开辟了一条途径。