[论文解读] iGVPT2 : an interface to computational chemistry packages for anharmonic corrections to vibrational frequencies

Loïc Barnes, Baptiste Schindler+2 more|arXiv (Cornell University)|Apr 7, 2017

Advanced Chemical Physics Studies参考文献 9被引用 3

一句话总结

iGVPT2 是一款免费开源软件，通过二阶振动微扰理论（VPT2）及其变体（VPT2+K、DCPT2、HDCPT2）计算非谐振动频率校正，并与多种量子化学程序接口对接。它支持混合量子力学/分子力学（QM//MM）计算以及力常数的并行计算，将实际计算时间减少约 ~f³（f 为正常模式数），从而实现对大分子（包括生物分子）的高效非谐红外光谱预测。

ABSTRACT

iGVPT2 is a program for computing anharmonic corrections to vibration frequencies, based on force field expansion of the potential energy surface in normal mode coordinates. It includes second order vibrational perturbation theory (VPT2) algorithm and its derived methods (VPT2+K, DCPT2, HDCPT2). iGVPT2 is interfaced with several computation chemistry packages to compute the potential energies and dipoles derivatives. The second, third and quartic derivatives can be computed at the same level of theory but they can be also computed using different methods via one or two computational packages. iGPVT2 includes also a very fast hybrid QM//MM approach for biomolecules. It is provided free-of-charge for non-commercial research (see \url{https://sites.google.com/site/allouchear/igvpt2}).

研究动机与目标

• 开发一种灵活且高效的接口，用于使用 VPT2 及其变体计算非谐振动频率校正。
• 实现混合级量子化学计算，其中谐振频率在高阶理论水平计算，而非谐校正在较低水平进行。
• 通过类似插件的接口支持多种计算化学程序（如 ORCA、Gaussian、Molpro、OpenMOPac），提升互操作性。
• 实现一种快速的混合 QM//MM 方法，用于大生物分子，显著降低计算成本，同时保持精度。
• 提供高度并行化的框架，用于计算三阶和四阶导数，将实际计算时间减少约 ~f³。

提出的方法

• iGVPT2 使用二阶振动微扰理论（VPT2）计算非谐校正，包括 VPT2+K、DCPT2 和 HDCPT2 等变体。
• 它与多种量子化学程序（如 ORCA、Gaussian、Molpro、FireFly、GAMESS-US、DFTB+）接口，用于计算势能、偶极矩和力的导数。
• 该软件支持混合计算：谐振频率在高阶方法（如 B3LYP/def2-TZVPP）下计算，而三阶和四阶力常数在较低阶方法（如 MMFF94 或半经验方法）下计算。
• 谐振频率的数值导数通过有限差分法计算，用户可自定义步长，且该过程在原子和正常模式之间实现并行化。
• 采用 DFT 计算 QM 区域，其余部分使用 MMFF94 的混合 QM//MM 方法，实现对大生物分子的快速且准确的非谐频率预测。
• 代码通过并行化实现，将实际计算时间减少约 ~f³（f = 正常模式数），显著加速三阶和四阶导数的计算。

实验结果

研究问题

• RQ1 是否能够通过单一接口在多种量子化学程序中高效计算各向异性振动频率校正？
• RQ2 iGVPT2 中的混合 QM//MM 方法是否在大幅降低大生物分子计算成本的同时保持了准确性？
• RQ3 通过并行化是否可将三阶和四阶导数的实际计算时间减少约 ~f³？
• RQ4 与实验数据及其他软件实现相比，iGVPT2 计算的非谐频率精度如何？
• RQ5 混合级计算（如高阶方法计算谐振频率，低阶方法计算非谐校正）是否可在不进行所有导数的高阶计算的前提下获得可靠结果？

主要发现

• 对于质子化的二肽 GlyGlyH+，使用全 DFT（B3LYP/def2-TZVPP）计算非谐频率耗时 109,888 个 CPU 小时，而使用 2000 个核心并行计算后，实际计算时间缩短至 56 小时。
• 采用混合 QM//MM 方法（谐振频率使用 B3LYP/def2-TZVPP，非谐校正使用 MMFF94）后，相同计算仅耗时 0.5 小时（三阶和四阶导数计算仅 6 秒）。
• 混合 QM//MM 方法在先前研究中已验证可高精度重现实验非谐频率。
• 在测试分子中，最大偏差（MAX）低于 1224 cm⁻¹，与参考计算相比平均绝对偏差（MAD）低于 2.5 cm⁻¹。
• 当使用一致的步长时，非谐频率的均方根偏差（RMSD）低于 4.8 cm⁻¹，表明结果具有高度一致性和可靠性。
• 该软件通过并行化和混合方法显著加速计算，实现了对大体系（包括生物分子）的高精度非谐频率预测。