[论文解读] 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

一句话总结

本文挑战了Palmer等人（2014年）提出的观点，即ST2水模型的分子动力学模拟表现出两种亚稳态液相之间的稳定界面，暗示存在亚稳态双液相临界性。本文认为，所观察到的行为是由于平衡不足和亚稳态约束所致的假象，这些因素将涨落尺寸限制在约2–3纳米，使得此类系统中无法实现真正的相共存。

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.

研究动机与目标

• 挑战Palmer等人（2014年）提出的ST2水模型表现出两种亚稳态液相之间稳定界面的主张。
• 证明模拟中观察到的双相行为是由于平衡不足和亚稳态约束所导致的人工效应。
• 表明过冷水中的涨落长度尺度从根本上受亚稳态寿命限制，无法观察到真正的临界性。
• 指出Palmer等人研究中自由能计算和采样方法存在的技术缺陷。
• 阐明模拟中看似双液相的行为并不意味着真实的相共存或临界性，尤其是在2–3纳米以上的尺度上。

提出的方法

• 分析过冷水中的时间与长度尺度，表明弛豫时间常数（τR）和亚稳态寿命（τMS）随温度降低而增长，但其比值始终低于约10^3。
• 应用不等式 ξ/a < (τMS/τR)^(1/3) 推导关联长度 ξ 的上限，得出在 T ≈ 215 K 时 ξ < 2–3 nm。
• 评估界面自由能 ΔF 随系统尺寸 N 的变化，检验 Δ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纳米，无法更大。
• Palmer 等人数据中的界面自由能 ΔF 无法区分 ΔF ∝ N^(2/3) 与 ΔF ∝ N 的情况，因此不支持存在稳定界面的证据。
• 模拟中 Q6 有序参数的动力学是非平稳的，并表现出滞后现象，表明平衡不足且采样不可逆。
• Palmer 等人自由能计算中显示出路径依赖性和低密度区域的不对称性，意味着误差至少是报告值的三倍。
• ST2 模型在305 K时即表现出玻璃态动力学，显著高于真实水（277 K），使其在模拟中极易产生不可逆的人工效应。
• ST2 模型中看似双液相的行为并非源于真实的相共存，而是冰粗化过程中短暂存在的介观涨落所致，而非临界性。