Palmer Jeremy C, Car Roberto, Debenedetti Pablo G
Faraday Discuss. 2013;167:77-94. doi: 10.1039/c3fd00074e.
We investigate the metastable phase behaviour of the ST2 water model under deeply supercooled conditions. The phase behaviour is examined using umbrella sampling (US) and well-tempered metadynamics (WT-MetaD) simulations to compute the reversible free energy surface parameterized by density and bond-orientation order. We find that free energy surfaces computed with both techniques clearly show two liquid phases in coexistence, in agreement with our earlier US and grand canonical Monte Carlo calculations [Y. Liu, J. C. Palmer, A. Z. Panagiotopoulos and P. G. Debenedetti, J Chem Phys, 2012, 137, 214505; Y. Liu, A. Z. Panagiotopoulos and P. G. Debenedetti, J Chem Phys, 2009, 131, 104508]. While we demonstrate that US and WT-MetaD produce consistent results, the latter technique is estimated to be more computationally efficient by an order of magnitude. As a result, we show that WT-MetaD can be used to study the finite-size scaling behaviour of the free energy barrier separating the two liquids for systems containing 192, 300 and 400 ST2 molecules. Although our results are consistent with the expected N(2/3) scaling law, we conclude that larger systems must be examined to provide conclusive evidence of a first-order phase transition and associated second critical point.
我们研究了ST2水模型在深度过冷条件下的亚稳相行为。使用伞形抽样(US)和加权自适应元动力学(WT-MetaD)模拟来研究相行为,以计算由密度和键取向序参数化的可逆自由能表面。我们发现,用这两种技术计算得到的自由能表面都清楚地显示出两种共存的液相,这与我们早期的US和巨正则蒙特卡罗计算结果一致[Y. Liu, J. C. Palmer, A. Z. Panagiotopoulos和P. G. Debenedetti, J Chem Phys, 2012, 137, 214505; Y. Liu, A. Z. Panagiotopoulos和P. G. Debenedetti, J Chem Phys, 2009, 131, 104508]。虽然我们证明了US和WT-MetaD产生了一致的结果,但估计后一种技术的计算效率要高一个数量级。因此,我们表明WT-MetaD可用于研究包含192、300和400个ST2分子的系统中分隔两种液体的自由能垒的有限尺寸标度行为。尽管我们的结果与预期的N(2/3)标度律一致,但我们得出结论,必须研究更大的系统,以提供一级相变和相关二级临界点的决定性证据。