[論文レビュー] Illusions of phase coexistence: Comments on "Metastable liquid-liquid transition ..." by J. C. Palmer et al., Nature 510, 385 (2014)

David Chandler|arXiv (Cornell University)|Jul 25, 2014

Material Dynamics and Properties参考文献 2被引用数 3

ひとこと要約

この論文は、パルマーら（2014年）が提唱した、ST2水モデルの分子動力学的シミュレーションにおいて、2つの準安定液体相の間に安定した界面が存在し、準安定二液相臨界性を示すという主張に反論している。著者らは、観察された振る舞いが、不十分な平衡化と準安定性制約に起因する錯覚であると主張しており、これによりフラクチュエーションのサイズが約2–3 nmに制限され、このような系では真の相分離が不可能であると結論づけている。

ABSTRACT

The recent paper cited above claims that a molecular simulation of one specific model of supercooled water establishes a stable interface separating two metastable liquid phases, which would imply the existence of metastable two-liquid criticality for that model. Here, we note that this claim conflicts with fundamental principles and with earlier work published in the extit{Journal of Chemical Physics}, and we show that the claim is unjustified by the data put forward to support the conclusion. Other technical problems are also noted.

研究の動機と目的

• パルマーら（2014年）が提唱した、ST2水モデルが2つの準安定液体相の間に安定した界面を示すという主張に反論すること。
• シミュレーションで観察された二相行動が、不十分な平衡化と準安定性制約に起因するアーティファクトであることを示すこと。
• 過冷却水中のフラクチュエーションの長さスケールが、準安定寿命によって本質的に制限されることを示し、真の臨界性の観測を不可能にする理由を説明すること。
• パルマーらの研究で用いられた自由エネルギー計算とサンプリング手法における技術的欠陥を特定すること。
• シミュレーションで見られる二液相に類似した振る舞いが、真の相分離や臨界性を示すものではないこと、特に2–3 nmより大きなスケールではそうであると明確にすること。

提案手法

• 過冷却水中の時間スケールと長さスケールを分析し、緩和時間定数（τR）と準安定寿命（τMS）が温度が低下するにつれて増加するが、その比が約10^3未満に保たれることを示している。
• 不等式 ξ/a < (τMS/τR)^(1/3) を適用して相関長さ ξ の上限を導出し、T ≈ 215 K 付近で ξ < 2–3 nm が得られた。
• 系のサイズ N に応じて界面自由エネルギー ΔF を評価し、安定界面では ΔF ∝ N^(2/3)、非共存相では ΔF ∝ N であると予想されるため、両者の違いを検証した。
• Q6秩序パラメータのダイナミクスの定常性と可逆性を評価し、平衡サンプリングを保証するためには緩和時間が必要であり、Q6自己相関時間の10^2倍以上である必要があると明言した。
• ρ–Q6 ウィンドウ全体にわたる可逆性テストを行い、ヒステリシスや経路依存性を検出し、非可逆的サンプリングと過大な誤差推定を示した。
• 対応状態解析を用いてST2モデルを他の水モデルと比較し、ST2が305 Kでガラス的動的挙動を示すのに対し、現実の水では277 Kで同様の挙動を示すことを明らかにした。これにより、アーティファクトのリスクが高まると結論づけた。

実験結果

リサーチクエスチョン

• RQ1 過冷却水中の二つの準安定液体相の間に安定した界面を分子シミュレーションで信頼性高く観測できるか？
• RQ2 準安定性制約がフラクチュエーションのサイズをどの程度制限し、真の臨界性の観測を妨げるか？
• RQ3 他のモデルや以前の研究とは異なり、ST2水モデルのシミュレーションではなぜ二液相に類似した振る舞いが観察されるのか？
• RQ4 不十分な平衡化時間と非定常的ダイナミクスが、自由エネルギー計算や相分離の解釈にどのように影響を与えるか？
• RQ5 ST2モデルの見かけ上の二液相的振る舞いの物理的根拠は何か？なぜ他の水モデルとは異なるのか？

主な発見

• 過冷却水中の最大相関長さ ξ は、準安定寿命 τMS によって制限され、約2–3 nmに抑えられる。
• パルマーらのデータにおける界面自由エネルギー ΔF は、ΔF ∝ N^(2/3) と ΔF ∝ N の区別がつかないため、安定界面の証拠は支持されない。
• シミュレーションにおける Q6 秩序パラメータのダイナミクスは非定常的であり、ヒステリシスを示しており、不十分な平衡化と非可逆的サンプリングを示している。
• パルマーらの自由エネルギー計算では、低密度領域に経路依存性と非対称性が認められ、誤差が報告値の少なくとも3倍以上大きい可能性がある。
• ST2モデルは305 Kでガラス的動的挙動を示し、現実の水（277 K）と比較して顕著に高い温度で発現するため、シミュレーションで非可逆的アーティファクトを生じやすい。
• ST2モデルにおける見かけ上の二液相的振る舞いは、真の相分離によるものではなく、氷の粗大化過程における一時的で中間スケールのフラクチュエーションに起因するものであり、臨界性とは無関係である。