[論文レビュー] A Reversible Unwrapping Algorithm for Constant Pressure Molecular Dynamics Simulations

Martin Kulke, Josh V. Vermaas|arXiv (Cornell University)|Nov 23, 2021

Nanopore and Nanochannel Transport Studies参考文献 41被引用数 8

ひとこと要約

本論文は、定圧分子動力学シミュレーションにおけるハイブリッド可逆アンラップアルゴリズムを提案し、周期的境界条件の体積変化を考慮することで、従来のヒューリスティック法および変位アンラップ法に生じるアーティファクトを是正する。この手法は、ボックスサイズの変化に基づく追加補正項を含む修正された変位スキームを用い、分子の幾何学的構造と正しい拡散係数を保つことで、長時間シミュレーションにおける精度を顕著に向上させる。

ABSTRACT

Molecular simulation technologies have afforded researchers a unique look into the nanoscale interactions driving physical processes. However, a limitation for molecular dynamics (MD) simulations is that they must be performed on finite-sized systems in order to map onto computational resources. To minimize artifacts arising from finite-sized simulation systems, it is common practice for MD simulations to be performed with periodic boundary conditions (PBC). However, in order to calculate specific physical properties, such as mean square displacements to calculate diffusion coefficients, continuous particle trajectories where the atomic movements are continuous and do not jump between cell faces are required. In these cases, modifying atomic coordinates through unwrapping schemes are an essential post-processing tool to remove these jumps. Here, two established trajectory unwrapping schemes are applied to 1us wrapped trajectories for a small water box. The existing schemes can result in spurious diffusion coefficients, long bonds within unwrapped molecules, and inconsistent atomic coordinates when coordinates are rewrapped after unwrapping. We determine that prior unwrapping schemes do not account for changing periodic box dimensions, and introduce an additional correction term to the existing displacement unwrapping scheme by von Bulow et al. to correct for these artifacts. We also demonstrate that the resulting algorithm is a hybrid between the existing heuristic and displacement unwrapping schemes. After treatment with this new unwrapping scheme, molecular geometries are correct even after long simulations. In anticipation for longer molecular dynamics trajectories, we develop implementations for this new scheme in multiple PBC handling tools.

研究の動機と目的

• NPT分子動力学シミュレーションにおける既存の軌道アンラップ手法の不正確さを是正すること。
• ボックス寸法が時間とともに変動する長時間シミュレーションにおいて、ヒューリスティック法および変位アンラップ法が失敗することを特定すること。
• 周期的ボックス体積の変化に対しても、正確な分子幾何学的構造と拡散挙動を維持できる新しいアンラップアルゴリズムを開発すること。
• 広く使用されているPBCツールへの統合を実現し、分子シミュレーションコミュニティにおける広範な採用を促進すること。

提案手法

• 提案手法は、ヒューリスティック法と変位アンラップ法の要素を統合し、主に変位ベクトルをアンラップの基盤として用いる。
• 時間ステップ間での周期的ボックス寸法の変化を補正するため、補正項 m(Li+1 − Li) を導入する。
• 前回のアンラップ座標、ラップされた変位、およびボックスサイズに応じた補正項を用いて、アンラップ座標を計算する。
• 補正項により、原子位置の補正が細胞の境界を越えた移動だけでなく、ボックスの膨張・収縮に対しても行われる。
• VMDのpbctools、fastpbc、qwrapにおいて実装され、一般的なMD解析ワークフローと互換性を持つ。
• この手法は可逆であり、幾何学的歪みを生じさせることなく、正確な再ラップが可能である。

実験結果

リサーチクエスチョン

• RQ1 ヒューリスティックアンラップ法は、ボックスサイズが変動する長時間のNPTシミュレーションにおいて、正確な拡散係数を生成するか？
• RQ2 変位アンラップ法は、ボックス寸法が変化する長時間のシミュレーションにおいても、正確な分子幾何学的構造を維持できるか？
• RQ3 ボックスサイズの変動は、定圧分子動力学における軌道アンラップの精度にどのように影響するか？
• RQ4 ボックスサイズの時間的変化を考慮するハイブリッドアンラップスキームは、幾何学的忠実性と拡散係数の正確性の両方を向上させられるか？
• RQ5 新しいアルゴリズムは可逆的であり、標準的なMD解析ツールへの統合に適しているか？

主な発見

• ヒューリスティック法は、ボックスサイズの変動に起因して、100ナノ秒程度で不正な拡散係数を生成し、その後失敗する。
• 変位法は、再ラップ後に結合長が不正確になり、特に長時間のシミュレーションで分子間距離が一貫性を欠く。
• 提案されたハイブリッドスキーム（式12）は、1マイクロ秒のシミュレーション後でも水分子の幾何学的構造が正確に維持され、結合長が約0.1 nmの近くを保つ。
• ハイブリッド法は、単精度浮動小数点演算の精度内に原子位置を保持し、幾何学的歪みではなく、大きな原子座標の変化に起因するRMSDの増加のみを示す。
• 新しいアルゴリズムは、m(Li+1 − Li) 項を組み込むことでボックスサイズの変化を適切に補正し、ヒューリスティック法による拡散の過大評価を解消する。
• この手法はVMDのpbctools、fastpbc、qwrapに実装され、分子動力学解析分野における広範な採用が可能となった。