Horstmann R, Vogel M
Institute of Condensed Matter Physics, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany.
J Chem Phys. 2021 Feb 7;154(5):054502. doi: 10.1063/5.0037080.
We use molecular dynamics simulations to study relations between thermodymamic, structural, and dynamical properties of TIP4P/2005 water models with systematically reduced partial charges and, thus, weaker hydrogen bonds. Observing a crossing of isochores in the P-T diagram, we show that these water-like models have a readily accessible liquid-liquid critical point (LLCP) associated with a transition between high-density liquid (HDL) and low-density liquid (LDL) forms and determine the dependence of the critical temperature T, pressure P, and density ρ on the charge-scaling factor from fits to a two-structure equation of states. The results indicate that the water-like models exhibit liquid polyamorphism in a wide range of interaction parameters. Considering elongated systems, we observe a decomposition into extended and stable HDL-like and LDL-like regions at appropriate pressures and low temperatures and analyze the respective structural and dynamical properties. We show that the diverse local order results in very different correlation times of local dynamics, while the fragility is hardly changed. The results yield insights into the origin of a dynamical crossover, which is observed when lowering the temperature along isobars and was previously interpreted in terms of a fragile-to-strong transition. Our findings imply that the effect does not involve two liquid phases with an exceptionally large difference of the fragility but rather a high temperature dependence near the LLCP results from a rapid conversion from HDL-like environments with faster dynamics to LDL-like ones with slower dynamics.
我们使用分子动力学模拟来研究TIP4P/2005水模型的热力学、结构和动力学性质之间的关系,这些模型具有系统降低的部分电荷,因此氢键较弱。通过观察P-T图中等容线的交叉,我们表明这些类水模型具有易于达到的液-液临界点(LLCP),该临界点与高密度液体(HDL)和低密度液体(LDL)形式之间的转变相关,并通过拟合两结构状态方程确定临界温度T、压力P和密度ρ对电荷缩放因子的依赖性。结果表明,类水模型在广泛的相互作用参数范围内表现出液体多晶型现象。考虑到拉长的系统,我们观察到在适当的压力和低温下分解为扩展的和稳定的HDL类和LDL类区域,并分析各自的结构和动力学性质。我们表明,不同的局部有序导致局部动力学的相关时间非常不同,而脆性几乎不变。结果揭示了动态交叉的起源,当沿着等压线降低温度时会观察到动态交叉,并且之前是根据从脆性到强的转变来解释的。我们的发现意味着这种效应并不涉及脆性差异极大的两个液相,而是在LLCP附近的高温依赖性源于从具有较快动力学的HDL类环境到具有较慢动力学的LDL类环境的快速转变。