Department of Chemistry, University of California, Berkeley, California 94720, USA.
J Chem Phys. 2013 Jun 7;138(21):214504. doi: 10.1063/1.4807479.
This paper extends our earlier studies of free energy functions of density and crystalline order parameters for models of supercooled water, which allows us to examine the possibility of two distinct metastable liquid phases [D. T. Limmer and D. Chandler, J. Chem. Phys. 135, 134503 (2011) and preprint arXiv:1107.0337 (2011)]. Low-temperature reversible free energy surfaces of several different atomistic models are computed: mW water, TIP4P/2005 water, Stillinger-Weber silicon, and ST2 water, the last of these comparing three different treatments of long-ranged forces. In each case, we show that there is one stable or metastable liquid phase, and there is an ice-like crystal phase. The time scales for crystallization in these systems far exceed those of structural relaxation in the supercooled metastable liquid. We show how this wide separation in time scales produces an illusion of a low-temperature liquid-liquid transition. The phenomenon suggesting metastability of two distinct liquid phases is actually coarsening of the ordered ice-like phase, which we elucidate using both analytical theory and computer simulation. For the latter, we describe robust methods for computing reversible free energy surfaces, and we consider effects of electrostatic boundary conditions. We show that sensible alterations of models and boundary conditions produce no qualitative changes in low-temperature phase behaviors of these systems, only marginal changes in equations of state. On the other hand, we show that altering sampling time scales can produce large and qualitative non-equilibrium effects. Recent reports of evidence of a liquid-liquid critical point in computer simulations of supercooled water are considered in this light.
本文扩展了我们早期对过冷水模型的密度和晶体有序参数自由能函数的研究,这使我们能够检验两种不同亚稳液相的可能性[D. T. Limmer 和 D. Chandler,J. Chem. Phys. 135, 134503(2011)和预印本 arXiv:1107.0337(2011)]。计算了几种不同原子模型的低温可逆自由能表面:mW 水、TIP4P/2005 水、Stillinger-Weber 硅和 ST2 水,最后一种比较了三种不同的长程力处理方法。在每种情况下,我们都表明存在一种稳定或亚稳液相,以及一种冰状晶体相。这些系统中的结晶时间尺度远远超过过冷亚稳液相中的结构弛豫时间尺度。我们展示了这种时间尺度的巨大分离如何产生低温液-液转变的错觉。暗示两种不同液相亚稳性的现象实际上是有序冰状相的粗化,我们使用分析理论和计算机模拟来阐明这一点。对于后者,我们描述了计算可逆自由能表面的稳健方法,并考虑了静电边界条件的影响。我们表明,模型和边界条件的合理改变不会对这些系统的低温相行为产生定性变化,只会对状态方程产生微小变化。另一方面,我们表明改变采样时间尺度会产生大的和定性的非平衡效应。最近关于在过冷水中的计算机模拟中存在液-液临界点的报告就是在这种情况下考虑的。