[論文レビュー] Implementing Reactivity in Molecular Dynamics Simulations with Harmonic Force Fields

Jordan J. Winetrout, Krishan Kanhaiya+8 more|arXiv (Cornell University)|Jul 30, 2021

Polymer crystallization and properties被引用数 7

ひとこと要約

この論文は、反応的分子動力学シミュレーションを可能にするエネルギー保存型モースポテンシャルに置き換えた古典的調和結合ポテンシャルを用いる、反応性インターフェース力場（IFF-R）を導入する。各結合種別ごとに解釈可能な3つのモースパラメータを用い、解離時のエネルギーをゼロに設定することで、ポリマー、タンパク質、金属を含む多様な材料における結合の切断と形成を可能にするとともに、精度を維持し、従来の反応的手法と比較して約30倍の高速化を達成する。

ABSTRACT

The simulation of chemical reactions and mechanical properties including failure from atoms to the micrometer scale remains a longstanding challenge in chemistry and materials science. Bottlenecks include computational feasibility, reliability, and cost. We introduce a method for reactive molecular dynamics simulations using a clean replacement of non-reactive classical harmonic bond potentials with reactive, energy-conserving Morse potentials, called the Reactive INTERFACE Force Field (IFF-R). IFF-R is compatible with force fields for organic and inorganic compounds such as IFF, CHARMM, PCFF, OPLS-AA, and AMBER. Bond dissociation is enabled by three interpretable Morse parameters per bond type and zero energy upon disconnect. Use cases for bond breaking in molecules, failure of polymers, carbon nanostructures, proteins, composite materials, and metals are shown. The simulation of bond forming reactions was included via template-based methods. IFF-R maintains the accuracy of the corresponding non-reactive force fields and is about 30 times faster than prior reactive simulation methods.

研究の動機と目的

• 原子的からミクロンスケールでの化学反応および機械的破壊のシミュレーションという長年の課題に取り組む。
• 特に計算コストが高く、信頼性に欠ける反応的分子動力学における計算のボトル neck を克服する。
• 既存の非反応的力場の精度を保ちつつ、結合の解離と再形成を可能にする手法を開発する。
• IFF、CHARMM、OPLS-AA、AMBER、PCFF などの広く使用されている力場との互換性を確保する。
• 物理的妥当性を損なわず、反応的シミュレーションにおける顕著な性能向上を達成する。

提案手法

• 古典的調和結合ポテンシャルをエネルギー保存型モースポテンシャルに置き換えることで、結合の切断と形成を可能にする。
• 結合強度、平衡長、解離エネルギーを制御するため、各結合種別ごとに3つの解釈可能なモースパラメータを定義する。
• 無限距離でのポテンシャルエネルギーをゼロ（切断時のエネルギーをゼロ）に設定することで、エネルギー保存を確保する。
• 既存の力場（例：IFF、CHARMM、OPLS-AA、AMBER、PCFF）と統合することで、有機および無機材料へのスムーズな適用を可能にする。
• テンプレートベースの手法を用いて結合形成反応をシミュレートし、動的反応経路を実現する。
• 修正されたポテンシャルエネルギー関数を用いた標準的な分子動力学アルゴリズムを採用することで、計算効率を維持する。

実験結果

リサーチクエスチョン

• RQ1 非反応的調和力場の精度を保ちつつ、結合の切断を可能にする反応的力場を構築することは可能か？
• RQ2 モースポテンシャルを用いて、物理的に整合的かつ計算的に効率的な方法で結合解離をモデル化するにはどうすればよいか？
• RQ3 新規手法は、ポリマー、カーボンナノ構造、複合材料における複雑な破壊機構をどの程度正確にシミュレートできるか？
• RQ4 タンパク質や金属を含む多様な化学系に一般化可能か。再パラメータ化を最小限に抑えることができるか？
• RQ5 既存の反応的分子動力学手法と比較して、どの程度の性能向上が達成できるか？

主な発見

• IFF-Rは、調和ポテンシャルのクリーンな置き換えとして、エネルギー保存型モースポテンシャルを用いることで、分子動力学シミュレーションにおける結合の切断と形成を可能にする。
• 従来の反応的シミュレーション手法と比較して、約30倍の高速化を達成し、計算の現実可能性を顕著に向上させる。
• 無限距離でのポテンシャルエネルギーをゼロに設定することで、結合解離がエネルギー保存を満たす形で正確にモデル化される。
• ポリマー、タンパク質、金属を含む多様な系において、元の非反応的力場の精度を保持する。
• 実用例では、ポリマーおよびカーボンナノ構造における機械的破壊のシミュレーション、およびテンプレートベース手法による動的結合形成の成功が確認された。
• IFF、CHARMM、OPLS-AA、AMBER、PCFF などの主要な力場と互換性があり、広範な適用可能性を有する。