[論文レビュー] Elastic Collision Based Dynamic Partitioning Scheme for Hybrid Simulations

Björn Kirchhoff, Elvar Ö. Jónsson+3 more|arXiv (Cornell University)|Mar 30, 2021

Theoretical and Computational Physics参考文献 64被引用数 8

ひとこと要約

本稿では、可動境界における弾性衝突を用いてエネルギーと運動量を保存する、ハイブリッド分子シミュレーションのための新規な動的パーティショニング手法SAFIRESを提案する。この手法は、衝突を検出するために時間ステップを動的に調整し、粒子の領域間通過を排除することで、静的または近似境界手法（例：FIRES）に比べて、制約なしのシミュレーションと顕著なずれのない熱力学的整合性を維持する。

ABSTRACT

The scattering-adapted flexible inner region ensemble separator (SAFIRES) is a partitioning scheme designed to divide a simulation cell into two regions to be treated with different computational methodologies. SAFIRES prevents particles from crossing between regions and resolves boundary events through elastic collisions of the particles mediated by the boundary, conserving energy and momenta. A multiple-time-step propagation algorithm is introduced where the time step is scaled automatically to identify the moment a collision occurs. If the length of the time step is kept constant, the new propagator reduces to a regular algorithm for Langevin dynamics, and to the velocity Verlet algorithm for classical dynamics if the friction coefficient is set to zero. SAFIRES constitutes the exact limit of the premise behind boundary-based methods such as FIRES, BEST, and BCC which take advantage of the indistinguishability of molecules on opposite sides of the separator. It gives correct average ensemble statistics despite the introduction of an ensemble separator. SAFIRES is tested in simulations where the molecules on the two sides are treated in the same way, for a Lennard-Jones (LJ) liquid and a LJ liquid in contact with a surface, as well as for liquid modelling simulations using the TIP4P force field. Simulations using SAFIRES are shown to reproduce the unconstrained reference simulations without significant deviations.

研究の動機と目的

• 内側と外側の領域に異なる計算手法を適用するハイブリッドシミュレーションにおいて、正しいアンサンブル統計を維持する境界に基づくパーティショニング手法を開発すること。
• 領域間の粒子の通過を排除し、適応的結合スキームで一般的な不連続性やエネルギーのずれを回避すること。
• FIRES、BEST、BCCなどの既存の境界手法を改善し、界面で正確な弾性衝突を強制することで運動量とエネルギーを保存すること。
• 動的で時間ステップに適応する伝搬を可能にし、極限状態では標準的なランジュバンまたはヴェルレ力学に還元できること。
• Lennard-Jones液体やTIP4P水を含む複数の系において、制約なしのシミュレーションと一致することを確認すること。

提案手法

• SAFIRESは、内側領域に存在する最も外側の粒子に基づいて、可動境界が拡張または収縮することで、シミュレーションセルを内側領域と外側領域にパーティショニングする。
• 境界を越える粒子間の弾性衝突を強制し、回転に基づく衝突モデルによりエネルギーと運動量を保存する。
• 複数の時間ステップを用いた伝搬法を導入し、衝突が発生する正確な瞬間を検出するために時間ステップを動的に調整する。
• 相対速度ベクトルを法線方向に整列させる回転行列を用い、その後、換算質量に基づく運動量の交換を実行する。
• ランジュバンまたは速度ヴェルレ伝搬法を統合し、摩擦係数が0.0または0に設定された場合に標準的な力学に還元されるように時間ステップをスケーリングする。
• 衝突検出、時間ステップの外挿、ベクトル代数と力の伝搬を用いた衝突解消を含む、疑似コードフレームワークを提供する。

実験結果

リサーチクエスチョン

• RQ1 弾性衝突に基づく動的境界手法は、領域間の粒子移動なしに、ハイブリッドシミュレーションで正しいアンサンブル統計を維持できるか？
• RQ2 SAFIRESにおける動的時間ステップ調整は、固定時間ステップまたはFIRESのような近似境界手法に比べて、どの程度精度を向上させるか？
• RQ3 Lennard-Jones液体およびTIP4P水系において、SAFIRESは制約なしのシミュレーションと同様の径方向分布関数（g(r)）をどの程度再現できるか？
• RQ4 弾性衝突モデルは、FIRESやハードスフィア境界手法で知られる境界での不自然な粒子蓄積を防げるか？
• RQ5 極限状態において、SAFIRESは標準的なランジュバンまたはヴェルレ力学に滑らかに還元可能か？これにより、既存のシミュレーションフレームワークとの互換性が保証されるか？

主な発見

• SAFIRESシミュレーションは、密度ρ = 1.374 g cm⁻³、温度T = 94.4 KのLennard-Jones液体について、制約なしの基準シミュレーションと同一の径方向分布関数（g(r)）を再現し、顕著なずれは認められない。
• 低密度Lennard-Jones液体（ρ = 1.100 g cm⁻³）においても、SAFIRESは正確なg(r)と境界位置確率分布を維持し、基準シミュレーションとよく一致する。
• SAFIRESにおける境界位置確率分布は、FIRESよりも広がっており、人工的な粒子蓄積を避けるより物理的に現実的な動的界面を示している。
• 摩擦係数を0.05に設定した場合、SAFIRESは標準的なランジュバン力学に正しく還元され、摩擦係数が0の場合には速度ヴェルレ力学に還元される。
• 動的時間ステップ調整により、衝突が高精度に検出され、バッファ領域やエネルギー補正スキームの必要性がなくなる。
• 弾性衝突モデルにより、FIRESやハードスフィア境界手法の既知の制限である境界での不自然な粒子蓄積が防止され、即時の運動量とエネルギー保存が確保される。