[論文レビュー] Overcoming timescale and finite-size limitations to compute condensation rates at physically relevant conditions

Matteo Salvalaglio, Pratyush Tiwary+2 more|arXiv (Cornell University)|Aug 7, 2015

nanoparticles nucleation surface interactions被引用数 3

ひとこと要約

本稿では、アルゴン凝縮の分子動力学的シミュレーションにおける時間スケール的制限および有限サイズ効果を克服する向上されたサンプリング手法を提示している。この手法により、物理的に現実的な過飽和度条件下での核化率の信頼性ある計算が可能となり、実験的時間スケール（数時間）に近い凝縮時間スケールが達成された。

ABSTRACT

Condensation of a liquid droplet from a supersaturated vapour phase is a classical nucleation phenomenon. As such it is challenging to compute its rate from atomistic molecular dynamics simulations. In fact at realistic supersaturation conditions condensation occurs on time scales that far exceed what can be reached by conventional molecular dynamics methods. Another known problem in this context is the distortion of the free energy profile associated to nucleation due to the small, finite size of typical simulation boxes. In this work the problem of time scale is addressed with a recently developed enhanced sampling method while contextually correcting for finite size effects. We demonstrate our approach studying the condensation of argon, and showing that characteristic nucleation times of the order of magnitude of hours can be reliably calculated, approaching realistic supersaturation conditions. We thus bridge the gap between what standard molecular dynamics can do and real physical systems.

研究の動機と目的

• 現実的な過飽和度における原子レベルの分子動力学的シミュレーションにおいて、核化率をシミュレートする際の極めて長い時間スケールの課題に対処すること。
• 小さなシミュレーションボックスサイズに起因する自由エネルギー曲線の歪みを補正すること。
• 物理的に現実的な核化率の計算を可能にすること。実験的時間スケール（数時間）に近づけること。
• 実際の物理的状態に近い過飽和度条件下で、アルゴンをモデル系として用いて手法の有効性を示すこと。

提案手法

• 最近開発された強化サンプリング手法を用い、分子動力学における稀な核化イベントのサンプリングを加速する。
• 小さなシミュレーションボックスに内在する有限サイズ効果を補正するため、自由エネルギー曲線に補正を施す。
• 強化サンプリングと有限サイズ補正を組み合わせることで、正確な核化率計算を実現する。
• 物理的システムに近い条件でアルゴン凝縮をシミュレートするためにこの手法を適用する。
• シミュレートされた核化時間スケールを実験的に関連する値と比較することで、手法の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 強化サンプリング手法は、従来の分子動力学が抱える時間スケール的制限を、凝縮のシミュレーションにおいて克服できるか？
• RQ2 小さなシミュレーションボックス内での核化の自由エネルギー曲線に、有限サイズ効果がどのように歪みをもたらすか？
• RQ3 有限サイズ補正を施すことで、核化率予測の正確性はどの程度向上するか？
• RQ4 組み合わせた手法により、物理的過飽和度条件下で実験的観察と同等の核化時間スケール（例：数時間）を達成できるか？
• RQ5 物理的過飽和度条件下でのアルゴン凝縮は、この手法を用いて十分な精度で計算可能か？

主な発見

• 本手法は、物理的に現実的な過飽和度条件下でのアルゴン凝縮の核化率を信頼性を持って計算できた。
• 数時間オーダーの凝縮時間スケールが安定的に計算され、実験的条件に近づいた。
• 自由エネルギー曲線に対する有限サイズ効果が効果的に補正され、計算精度が向上した。
• 強化サンプリング手法により、従来の分子動力学では到達不可能な稀な核化イベントにアクセス可能となった。
• 結果として、原子レベルのシミュレーションと現実の物理的システムとの間の大きなギャップが、核化現象において顕著に埋められた。