[論文レビュー] 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³（f = 正規モード数）まで短縮。

提案手法

• iGVPT2 は、VPT2 およびその変種（VPT2+K、DCPT2、HDCPT2）を用いて非調和補正を計算する。
• ORCA、Gaussian、Molpro、FireFly、GAMESS-US、DFTB+ などの複数の量子化学パッケージと連携し、ポテンシャルエネルギー、双極子モーメント、力定数の微分を計算。
• ハイブリッド計算をサポート：調和振動数は高水準の方法（例：B3LYP/def2-TZVPP）で計算し、三次および四次力定数は低水準の方法（例：MMFF94 や半経験的手法）で計算。
• 調和振動数の数値微分は、ユーザーが指定したステップサイズを用いた有限差分法で計算され、原子および正規モードごとに並列化。
• ハイブリッド QM//MM アプローチは、QM 領域に DFT、残りの領域に MMFF94 を使用し、大規模バイオ分子における非調和振動数予測を高速かつ高精度に実現。
• コードは並列化されており、実時間計算を最大で ~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⁻¹ 未満であり、高い一貫性と信頼性を示した。
• 本ソフトウェアは、並列処理およびハイブリッド手法を活用することで、バイオ分子を含む大規模系においても高精度な非調和振動数予測を可能にした。