Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA.
J Chem Phys. 2010 Sep 28;133(12):124107. doi: 10.1063/1.3481574.
Coarse-grained (CG) models provide a computationally efficient means for investigating phenomena that remain beyond the scope of atomically detailed models. Although CG models are often parametrized to reproduce the results of atomistic simulations, it is highly desirable to determine accurate CG models from experimental data. Recently, we have introduced a generalized Yvon-Born-Green (g-YBG) theory for directly (i.e., noniteratively) determining variationally optimized CG potentials from structural correlation functions. In principle, these correlation functions can be determined from experiment. In the present work, we introduce a reference state potential into the g-YBG framework. The reference state defines a fixed contribution to the CG potential. The remaining terms in the potential are then determined, such that the combined potential provides an optimal approximation to the many-body potential of mean force. By specifying a fixed contribution to the potential, the reference state significantly reduces the computational complexity and structural information necessary for determining the remaining potentials. We also validate the quantitative accuracy of the proposed method and numerically demonstrate that the reference state provides a convenient framework for transferring CG potentials from neat liquids to more complex systems. The resulting CG model provides a surprisingly accurate description of the two- and three-particle solvation structures of a hydrophobic solute in methanol. This work represents a significant step in developing the g-YBG theory as a useful computational framework for determining accurate CG models from limited experimental data.
粗粒化 (CG) 模型为研究原子细节模型无法涵盖的现象提供了一种计算效率高的方法。虽然 CG 模型通常是通过参数化来重现原子模拟的结果,但从实验数据中确定准确的 CG 模型是非常可取的。最近,我们引入了一种广义 Yvon-Born-Green (g-YBG) 理论,用于直接(即非迭代)从结构相关函数确定变分优化的 CG 势。原则上,这些相关函数可以从实验中确定。在本工作中,我们将参考状态势引入 g-YBG 框架中。参考状态定义了 CG 势的固定贡献。然后确定势的其余项,使得组合势提供对平均力多体势的最佳逼近。通过指定势的固定贡献,参考状态大大降低了确定剩余势所需的计算复杂性和结构信息。我们还验证了所提出方法的定量准确性,并数值证明了参考状态为将 CG 势从纯液体转移到更复杂系统提供了一个方便的框架。所得到的 CG 模型对疏水溶质在甲醇中的二聚体和三聚体溶剂化结构提供了惊人准确的描述。这项工作代表了将 g-YBG 理论发展为从有限的实验数据中确定准确 CG 模型的有用计算框架的重要一步。