Progress toward a better BOCS: Systematic coarse-graining with local density potentials

技术摘要：迈向更佳的 BOCS：利用局部密度势进行系统性粗粒化

问题陈述
自下而上的粗粒化（CG）模型为模拟原子级（AA）模型无法触及的长度和时间尺度提供了计算效率。虽然力匹配（FM）变分原理允许对 CG 势进行参数化，以最优地逼近精确的 CG 势，但标准的实现通常依赖于传统的对势（pair potentials）。这些对势难以捕捉非均匀系统（如气-液界面）的热力学和结构特性，并且往往在本体（bulk）与界面环境之间缺乏可迁移性。此外，体积势虽然在重现本体体系的压力-密度状态方程方面非常有效，但无法模拟界面。因此，需要一个系统性的框架，将依赖于局部密度（LD）和局部密度平方梯度（SG）的势函数纳入其中，以提高复杂软材料 CG 模型的高保真度和可迁移性。

方法论
作者描述了 BOCS 软件包的 5.0 版本以及针对 LAMMPS 的 PKG-BOCS 扩展包。

• 理论框架： 该软件使用力匹配变分原理对 CG 势进行参数化，该原理通过最小化映射系综中 AA 模型与 CG 模型之间力的均方差来实现。总非键势被扩展为包含三个部分：

  1. 传统的中心对势（$U_{pair}$）。
  2. 局部密度势（$U_{LD}$），其取决于每个粒子 $I$ 周围的瞬时局部密度 $\rho_I$。
  3. 平方梯度势（$U_{SG}$），其取决于局部密度的梯度平方 $|\nabla_I \rho_I|^2$。
     局部密度被定义为周围粒子贡献之和，并由权重函数 $w(r)$ 进行加权，该函数实现了五种不同的函数形式（例如 DPD、Lucy、Shell、smooth、sphere）。

• 软件实现：

  • BOCS v5.0： 此版本引入了用于从 AA 轨迹计算局部密度及其梯度的工具（ldcalc）。它采用广义 Yvon-Born-Green (g-YBG) 框架，在不需要键合、对项和 LD 项的显式力信息的情况下执行力匹配。该软件使用基于样条的基函数来表示力函数和系数，从而实现灵活的函数形式。它通过奇异值分解求解生成的正规方程，以确定最优参数。
  • PKG-BOCS (LAMMPS)： 在 LAMMPS 中实现了一个新的 pair style (ldd) 来模拟这些势函数。该实现支持复杂的分子拓扑、多种粒子类型，并考虑了分子内贡献和排除体积。

• 验证与应用： 作者通过三个基准测试验证了实现的有效性：

  1. 一个包含 400 个粒子的线性 LD 势系统。
  2. 一个包含 2 个粒子的二次 LD 势和线性 SG 势系统。
  3. 一个涉及键合、对、LD 和 SG 相互作用的 125 个分支分子（含 4 种位点类型）的复杂系统。
     此外，软件被用于参数化基于 SPC/E AA 模型的单位点 CG 水模型，并将通过本体模拟、平板（slab）模拟以及在本体模拟中增加 LD 势进行参数化的模型进行了对比。

关键结果

• 实现精度： 在所有基准测试中，BOCS v5.0 均定量地恢复了用于 LAMMPS 模拟的已知“真值”势和力。力匹配计算产生的误差极小（$\chi^2 \approx 0$），证明了该软件正确实现了 LD 和 SG 势的变分原理。
• 物理洞察： 研究表明，由 LD 和 SG 势产生的有效对相互作用对权重函数 $w(r)$ 的选择高度敏感。例如，即使底层的 LD 或 SG 函数形式相同，不同的权重函数也会导致显著不同的有效对势和力剖面。
• 水模型性能：
  • 使用本体水参数化的标准仅对势模型（PAIR-BULK）未能稳定液相，导致了汽化现象。
  • 从平板模拟中参数化的仅对势模型（PAIR-SLAB）虽然稳定了液体，但显著高估了密度和压缩率，且无法很好地重现液-气界面的宽度。
  • LD-BULK 模型（由本体模拟参数化的 pair + LD 势）显著提升了性能。它高精度地重现了 SPC/E 的状态方程（密度低估仅 0.1%），并较好地重现了液-气界面宽度和密度分布，尽管其参数化仅基于本体模拟。
  • 分析显示，LD 势根据局部环境调节对作用力，有效地在低密度区域加深了吸引作用（模拟界面行为），而在高密度区域保持斥力。

意义与主张
作者声称 BOCS 5.0 版本及其 PKG-BOCS 扩展为系统化参数化具有局部密度和平方梯度势的 CG 模型提供了鲁棒且准确的工具。论文指出，这些势函数通过提高结构保真度、热力学性质和可迁移性，为底向上 CG 模型（特别是对于液-气界面等非均匀系统）提供了显著改进。这项工作证明了 LD 势可以有效地捕获协同多体相互作用，并在无需在模拟中引入显式多体力项的情况下改善对密度波动的描述，因为环境依赖性是通过局部密度来捕获的。作者将该软件定位为开发复杂软材料 CG 模型现有工具包的有益补充。