[论文解读] Third-order many-body expansion of OSV-MP2 wavefunction for low-order scaling analytical gradient computation

Qiujiang Liang, Jun Yang|arXiv (Cornell University)|Jun 10, 2021

Advanced Chemical Physics Studies参考文献 111被引用 8

一句话总结

本论文提出了一种三阶多体展开（MBE(3)）的轨道特定虚轨道MP2（OSV-MP2）波函数，以实现对大分子的低标度解析梯度计算。通过利用轨道特定聚类和拟合积分及OSV弛豫中的稀疏性，该方法实现了能量计算的线性标度和梯度计算的二次标度，其精度与规范OSV-MP2相当，从而实现了对复杂体系（如亚卟啉和铜轮烷）的高效分子动力学模拟与几何优化。

ABSTRACT

We present a many-body expansion (MBE) formulation and implementation for efficient computation of analytical energy gradients from OSV-MP2 theory based on our earlier work (Zhou et al. J. Chem. Theory Comput. 2020, 16, 196-210). The third-order MBE(3) expansion of OSV-MP2 wavefunction was developed to adopt the orbital-specific clustering and long-range termination schemes, which avoids term-by-term differentiations of the MBE energy bodies. We achieve better efficiency by exploiting the algorithmic sparsity that allows to prune out insignificant fitting integrals and OSV relaxations. With these approximations, the present implementation is benchmarked on a range of molecules that show an economic scaling in the linear and quadratic regimes for computing MBE(3)-OSV-MP2 amplitude and gradient equations, respectively, and yields normal accuracy comparable to the original OSV-MP2 results. The MPI-3-based parallelism through shared memory one-sided communication is further developed for improving parallel scalability and memory accessibility by sorting the MBE(3) orbital clusters into independent tasks that are distributed on multiple processes across many nodes, supporting both global and local data locations in which selected MBE(3)-OSV-MP2 intermediates of different sizes are distinguished and accordingly placed. The accuracy and efficiency level of our MBE(3)-OSV-MP2 analytical gradient implementation is finally illustrated in two applications: we show that the subtle coordination structure differences of mechanically interlocked Cu-catenane complexes can be distinguished when tuning ligand lengths; and the porphycene molecular dynamics reveals the emergence of the vibrational signature arising from softened N-H stretching associated with hydrogen transfer, using an MP2 level of electron correlation and classical nuclei for the first time.

研究动机与目标

• 开发一种基于多体展开和OSV-MP2的低标度解析梯度方法，用于大分子相关体系。
• 通过将MBE(3)与轨道特定虚轨道（OSVs）结合，克服规范MP2梯度的高计算成本。
• 实现对大分子在MP2相关电子水平下的高效分子动力学与几何优化。
• 通过混合MPI-3单边通信与多节点间数据分布，实现高并行可扩展性与内存效率。
• 在基准体系（包括亚卟啉和机械互锁Cu-轮烷）上展示方法的精度与效率。

提出的方法

• 对OSV-MP2振幅和密度矩阵进行三阶多体展开（MBE(3)），将系统分解为局域轨道簇。
• 应用轨道特定聚类与长程截断，避免对MBE能量项逐项求导。
• 利用算法稀疏性，剔除不重要的拟合积分与OSV弛豫矢量，降低计算成本。
• 基于MPI-3实现并行计算，采用共享内存单边通信，将独立的MBE(3)簇分发至多个节点。
• 通过区分簇大小，合理分配全局与本地中间变量数组，优化内存访问。
• 结合RMA与SHM混合并行机制，支持能量与解析梯度的可扩展、高性能计算。

实验结果

研究问题

• RQ1 MBE(3)-OSV-MP2能否在大分子中实现能量的线性标度与梯度的二次标度？
• RQ2 MBE(3)-OSV-MP2在分子几何构型与动力学中能否准确再现规范OSV-MP2梯度？
• RQ3 通过调节配体长度，该方法能否分辨机械互锁分子（如Cu-轮烷）中的细微结构差异？
• RQ4 MBE(3)-OSV-MP2能否通过分子动力学模拟捕捉亚卟啉中氢转移的振动特征？
• RQ5 稀疏性与簇筛选对计算效率与精度有何影响？

主要发现

• MBE(3)-OSV-MP2方法实现了能量的线性标度与梯度的二次标度，支持大体系的高效计算。
• 采用标准MBE(3)簇选择策略，相关能恢复精度超过99.85%（相对于规范RI-MP2）。
• 梯度均方根偏差（RMSD）小于3×10⁻⁴ a.u.，与规范RI-MP2一致，达到原始OSV-MP2的精度水平。
• 中等大小分子（最多5000个拟合函数）的优化几何结构中，短键误差为0.1–0.5 pm，长非键距离误差为1.1–1.5 pm。
• 对40个原子的亚卟啉进行NVE分子动力学模拟（10 ps），成功捕捉到与互变异构态相关的N–H伸缩振动特征，与基准OSV-MP2结果一致。
• 在24个CPU上，C60@catcher（148个原子）与(H2O)190（570个原子）的每步梯度计算分别耗时约2.7小时与约46小时，验证了良好的可扩展性。