多体展开势能的分层截断方法

《Journal of Chemical Theory and Computation》：Hierarchical Truncations for Many-Body Expansion Potentials

　　

本研究提出了一种新的多体展开（Many-Body Expansion，MBE）高阶项截断策略，称为分层多体展开（Hierarchical Many-Body Expansion，HMBE）。该方法基于系统的分层划分，将体系分解为多层级的片段，原则上可应用于任意大型分子体系。研究人员通过对一系列（H2O）64结构的数值测试表明，与传统MBE相比，HMBE在计算显著更少的高阶项的情况下，仍能获得与全簇计算吻合良好的簇间相对能量及单个簇的结合能。此外，通过在分区界面处引入被称为“申根项（Schengen terms）”的特定多体项，可进一步提升计算精度。这项工作确立了HMBE作为一种有前景的框架，可用于模拟具有天然层级结构的大型体系（如蛋白质）。

研究背景与意义

利用高精度从头算（ab initio）方法模拟包含数千原子的大体系目前面临计算瓶颈，因为此类方法的计算标度随基函数数量（N_bas）呈陡峭增长，从局域密度泛函理论的N_bas³到耦合簇方法的N_bas⁷不等。为解决这一问题，碎片法（Fragmentation methods）通过将体系分解为小片段并利用现代计算机的高度并行架构来加速计算，成为过去几十年的研究热点。碎片法通常建立在多体展开（MBE）或多体分析（Many-Body Analysis，MBA）框架下，即将体系总能量表示为各单体能量与各片段间多体相互作用能量之和。尽管完整的MBE是精确的，但其高阶项数量极其庞大，n阶项的数量为组合数C(m, n)，其中m为片段总数。传统截断方法通常直接忽略所有超过n个片段的相互作用（即MBE-n），虽然随着截断阶数提高会收敛至全MBE结果，但收敛速度往往较慢，且计算成本依然高昂。因此，开发一种能系统性减少高阶项计算量的高效截断策略具有重要学术价值。本论文发表于《Journal of Chemical Theory and Computation》。

关键技术方法

研究人员采用了约束k-means聚类算法对（H₂O）₆₄团簇的9个不同结构进行分层划分，构建了三种双层（2,3）-HMBE分区（4×16、16×4、8×8）和一种三层（2,3,4）-HMBE分区（4×4×4）。所有计算均在Hartree-Fock（HF）理论水平结合aug-cc-pVDZ基组下进行，并使用Psi4软件包完成。研究人员分别实现了自上而下和自下而上两种HMBE构建算法。为提升精度，引入了基于片段分子轨道（Fragment Molecular Orbital，FMO）理论截断阈值筛选的“申根项（Schengen terms）”进行补充计算。计算结果与全簇计算及传统MBE-3结果进行了对比验证。

研究结果

1. 分层截断：总体考量

HMBE的核心在于设定各层级的截断阶数T₁、T₂、T₃（分别对应第一、二、三层），且要求T₁≤ T₂≤ T₃。当三者相等时，HMBE退化为传统MBE。这种设定符合直觉，即层级越低片段越小，高阶相互作用通常越重要。

2. 分层截断：自上而下法

该策略从顶层向底层逐级展开能量表达式。以（2,3,4）-HMBE为例，顶层（第一层）仅计算单体和两体相互作用；第二层计算单体及最多三体相互作用；第三层计算单体及最多四体相互作用。通过层级间的嵌套定义，大幅减少了跨层级的远程高阶项计算。

3. 分层截断：自下而上法

该策略首先写出完整的传统MBE（如MBE-4），然后根据HMBE的截断规则（T₁=2, T₂=3）剔除不符合要求的项。例如，若三个片段属于三个不同的第一层分区，则其对应的三体项将被删除。这种方法对现有MBE代码改动极小，只需将剔除项设为零即可。

4. 分层截断：可扩展性

理论分析表明，HMBE的有效计算标度n_eff介于T₁和T₃之间。对于（2,3,4）-HMBE，其四体项的数量远低于传统MBE-4。计算成本相对于传统MBE降低了约M²M_I倍（M为第一层片段数，M_I为第二层片段数）。第一层分区越多，计算量节省越显著。

5. 分层截断：增强申根项

针对分区界面处可能存在的重要相互作用，研究人员引入了“申根项”。这些项涉及来自不同上层分区的下层片段间的多体相互作用。采用FMO理论的截断阈值（基于氧原子间、氢原子间及氧氢原子间的约化距离R'）来决定是否纳入这些项，从而在几乎不增加计算量的情况下显著提升界面描述精度。

6. 计算结果与讨论

在（H₂O）₆₄测试中，所有HMBE分区均能定性复现全簇计算的势能面。双层HMBE的精度随第二层片段数增加而提升，（2,3）-4×16分区表现最佳（每片段平均绝对误差MUE/f为0.036 kcal/mol）。三层（2,3,4）-4×4×4分区精度与之相当。加入申根项后，精度进一步提升，双层分区的MUE/f降至0.013-0.023 kcal/mol，与MBE-3（0.019 kcal/mol）相当，但计算项数量显著少于MBE-3。在结合能计算中，HMBE同样表现出良好的一致性，三层分区甚至优于MBE-3。

结论总结

本研究成功建立了HMBE这一新的多体展开截断框架。通过利用体系的天然或人工层级结构，HMBE能在保证计算精度的前提下，大幅降低高阶项的计算数量，其有效标度介于各层级截断阶数之间。结合申根项的补充，HMBE能以远低于传统MBE的计算成本达到可比甚至更高的精度。该方法特别适用于具有天然层级结构的大体系（如蛋白质）的高精度模拟，未来结合静电嵌入（Electrostatically Embedded Many-Body，EEMB）等策略有望进一步加速收敛并拓展应用。