A simple time coarse graining method for molecular dynamics simulations of liquids

あなたは、ある賑やかな都市の映画を見ようとしていると想像してください。しかし、カメラがあまりにも高精細であるため、誰かが瞬きをする瞬間や、歩行者の指がわずかに動く様子までもがすべて捉えられてしまいます。物語を完全に理解するためには、すべてのフレームを見なければなりません。問題は、その映画があまりにも長いため、世界最速のスーパコンピュータであっても、わずか数分間の映像をレンダリングするのに数年かかるという点です。これは、液体中の原子がどのように動くかを追跡しようとする「分子動力学（Molecular Dynamics）」シミュレーションを実行する科学者たちが日々直面している苦闘です。

通常、これらのシミュレーションは、文字通り「壁」に突き当たります。原子は、同じ空間を占有することができない跳ね返るボールのようなものです。もしシミュレーションが、あまりにも小さなステップで原子を動かしすぎると、別の原子に衝突してしまい、「無限エネルギー」エラーが発生してプログラム全体がクラッシュしてしまいます。これを避けるために、科学者たちは非常に小さく慎重なステップを踏まなければならず、それがシミュレーションを信じられないほど遅くさせているのです。

大きなアイデア：デコボコした道の平滑化

この論文の中で、研究者たちは賢い問いを投げかけました。「もし、クラッシュすることなく、より大きなステップを踏めるとしたらどうだろうか？」彼らは、コンピュータは無限に小さなステップを踏むことはできないため、シミュレーションはすでに完璧な決定論的な映画というよりは、一種の「推測ゲーム」のようなものになっていることに気づきました。彼らは、シミュレーションを確率的な手法（stochastic method）として扱い、時間の「粗視化（coarse graining）」を試みることにしました。

原子の動きを、山に登ろうとするハイカーに例えて考えてみましょう。元のポテンシャルエネルギーの景観は、底に恐ろしい垂直の絶壁（「ハードコア」の壁）がある山のようなものです。もしハイカーが少しでも崖の方へ滑り落ちれば、転落してしまいます。安全を確保するために、ハイカーは非常に小さく、神経質なステップを踏まなければなりません。

研究者たちはシンプルな解決策を提案しました。その恐ろしい垂直の絶壁を、緩やかで滑らかな曲線に置き換えたらどうでしょうか？ 彼らはこれを「二次関数則（quadratic law）」と呼びました。 「止まれ、さもなくば死ぬ」という硬い壁の代わりに、「おっと、スピードを落として。でも、まだ転がり続けることはできるよ」と告げる、柔らかく弾力のある丘を置いたのです。

どのようにテストしたか

このトリックが機能するかどうかを確認するために、チームは単に推測するのではなく、ダンベル型の分子（2つの原子がピーナッツのように結合したもの）からなる液体のシミュレーションを実行しました。彼らは標準的な温度である500K（モデル液体の融点より100度低い温度）を使用しました。

彼らは、恐ろしい短距離の壁が滑らかな二次曲線に置き換えられた、新しい「有効」ポテンシャルを作成しました。彼らは、Rと呼ばれる特定のカットオフ距離を変えることで、これをテストしました。

• Rが小さい場合（2.65 Å未満）、何も変わりませんでした。原子は依然として通常通りに振る舞い、シミュレーション結果は元の遅い手法と同一でした。
• Rが2.65 Åより大きくなると、「絶壁」は「スロープ」へと変わりました。

結果：原子のためのタイムマシン

ここからがエキサイティングなところです。壁を滑らかにした結果、彼らは通常の4倍のタイムステップを取ってもシミュレーションがクラッシュしないことを発見しました。しかし、魔法はそれだけではありませんでした。

新しい高速シミュレーションを古い低速シミュレーションと比較することで、彼らは驚くべき発見をしました。新しいシミュレーションは、元のシミュレーションと全く同じ動きのパターンを、より速く生成していたのです。

• 彼らは、原子が互いに離れていく様子を、「明確なファン・ホーブ相関関数（distinct Van Hove correlation function）」を用いて測定しました。
• 滑らかな壁を持つ高速シミュレーションは、時間シフトを伴いつつも、低速シミュレーションの結果と一致することを見出しました。
• 特定のテストでは、50倍の加速因子を達成しました。これは、以前は50時間かかっていたシミュレーションが、今では1時間で完了でき、なおかつ原子は同じ調べに合わせて踊り続けることができるということを意味します。

何が変わり、（何が変わらなかったのか）

論文では、何が起こらなかったのかについても注意深く指摘しています。長距離相互作用（遠く離れた原子同士のやり取り）は、まったく同じまま残されました。彼らはシステム全体を過度に単純化する「平均場（mean field）」のアプローチは用いていません。彼らは短距離の壁のみに手を加えたのです。

しかし、小さな副作用がありました。「スロープ」によって原子が動きやすくなったため、シミュレーションにおける「協同運動（cooperative motions）」がわずかに減少しました。元の液体では、原子はしばしばグループとして動きます（群衆が「ウェーブ」をするようなものです）。滑らかなシミュレーションでは、これらのグループによる波は少し不明瞭になりました。「非ガウスパラメータ（Non-Gaussian parameter）」（グループの動きが単純なランダムウォークからどれだけ逸脱しているかを測定するもの）は、カットオフ R の増加に伴って減少しました。これは、「時間の粗視化」が液体の複雑な集団的震動を滑らかにしたことを示唆しています。

結論

著者らは、この手法が粘性液体や、液体が固体のようになる「ガラス転移」に近づいている系のシミュレーションを高速化するための強力なツールであることを示唆しています。恐ろしい短距離の壁を滑らかな二次関数則に置き換えることで、より大きなタイムステップを取ることが可能になり、元のポテンシャルと同じダイナミクスを得ながら、大幅なスピードアップを実現できることを見出しました。

彼らは、この手法がタイムステップ中の相互作用の方向の変化（これは一部の複雑なシステムでは重要になる可能性がある）を考慮していないことも述べていますが、テストされた過冷却液体については、見事に機能しました。彼らはさらに、この平滑化効果が「時間－温度超位性（time-temperature superposition）」と呼ばれる物理学の有名な原理に関連している可能性さえ示唆しており、タイムステップを変えることは数学的に温度を変えることと似ていることを示唆しています。しかし現時点では、これは彼らのシミュレーションに基づく示唆であり、宇宙の法則として証明されたものではない、と強調しています。

要約すると、恐ろしい絶壁を緩やかなスロープに変えることで、科学者たちは、筋書きを失うことなく、原子の世界の映画を早送りで見る方法を見つけたのです。