使用逾渗方法揭示水的反常现象

《ACS Omega》：Unravelling Water Anomalies Using a Percolation Approach

　　

计算研究表明，水在277 K时的最大密度以及在1 bar下228 K时与实验测量压缩率的发散，源自低密度和高密度分子群体之间的竞争。然而，这些异常现象（包括203 K时的最小密度）之间的关系仍存在争议。研究人员使用TIP4P/2005模型的分子动力学模拟，沿等温线和等压线研究了亚稳水（metastable water）的压力-温度行为，以证明在两种情况下，密度和压缩率异常均源于低密度和高密度逾渗团簇（percolating clusters）的出现、共存和消失。研究人员将这些行为与亚稳水相分离的不同阶段联系起来，这些阶段发生在样品密度急剧增加的区域。研究结果与液-液临界点（liquid–liquid critical point, LLCP）场景一致，同时建立了一个基于逾渗的现象学框架（percolation-based framework），系统地组织了不同温度和压力下的水异常现象。该框架可能解释在氧化物和金属玻璃中观察到的逾渗转变，为描述纯物质中的非晶转变提供了一种视角。

**论文解读：《ACS Omega》发表——基于逾渗方法揭示水的反常现象**

**研究背景与问题**
水的反常现象（如大气压下277 K出现密度最大值、228 K附近等温压缩率发散，以及203 K出现密度最小值）对地球环境与生命维持至关重要。现有计算研究表明，这些异常源于低密度（low-density, LD）和高密度（high-density, HD）分子群体之间的竞争。然而，这些异常现象之间的内在关联，特别是密度最小值与压缩率发散的关系，仍存在争议。此外，液-液临界点（LLCP）假设虽被广泛讨论，但受限于深过冷区的快速结晶（“无人区”），实验与模拟均难以直接验证。为统一解释不同热力学条件下的水异常，研究人员引入逾渗（percolation）概念，探索其作为几何与现象学框架的可行性，以建立连通性驱动相变与异常现象之间的系统联系。

**研究方法与关键发现**
研究人员使用TIP4P/2005力场的分子动力学（MD）模拟，在NPT系综中对包含8192、2825、896和336个H₂O分子的样品进行广泛温度（T）和压力（p）范围的模拟。采用leapfrog积分器（2 fs步长）、随机速度重标热浴和Parrinello–Rahman恒压器，键和角约束通过六阶LINCS算法实现，范德华与实空间库仑截断为12 ?（最小样品为9 ?），长程静电使用粒子网格Ewald（PME）方案，傅里叶网格间距1.6 ?（最小样品1 ?）。样品先在300 K平衡，然后以1 K/ns速率等压冷却至80 K（896分子样品在1 bar下冷却至10 K），或等温压缩至15 kbar（0.1 kbar/ns）。冷却与压缩速率基于前人TIP4P/2005研究选定，确保结果收敛。分子分类基于3.5 ?截断距离内的氧原子邻居数Z：Z=4为LD，Z=5–7为HD，Z≥8为甚高密度（VHD）。逾渗序参量P_∞定义为属于逾渗团簇的分子分数，关联长度ξ通过连通性分析计算。玻璃转变温度T_g由焓-温度曲线拟合得到。

**研究结果**

**3.1 冷却过程中水的逾渗**
在1 bar等压冷却下，初始为HD连通网络，随温度降低，LD团簇以分形方式从多个成核点生长。P_∞(LD)从0连续增至1，标志着HDL与LDL相分离的起始及逾渗共存区（白色区域）。在~228–233 K，?(LD)=?(HD)，κ_T达到峰值，与实验观测的228 K压缩率发散温度一致。进一步冷却，LD逾渗团簇粗化，HD团簇破碎，P_∞(HD)降至0，ξ(HD)在密度最小值附近达到峰值，标志HD逾渗极限。团簇尺寸分布n(s)在逾渗转变点服从幂律s^–τ，τ≈2.1，接近普通逾渗值（2.18）。三个关键温度点：277 K时ξ(LD)与ρ同时达到最大值；~230 K时κ_T峰值对应?(LD/HD)相等；215 K附近ξ(HD)峰值对应HD逾渗团簇灭绝。

**3.2 压力驱动的逾渗**
等温压缩（T > T_g）触发LDL→HDL逾渗转变，T < T_g时对应LDA→HDA转变。随压力升高，P_∞(HD)从0升至1，κ_T增强，表明LD/HD逾渗共存与相分离。压缩下，κ_T峰值在T ? T_c时与ξ(LD)峰值对齐，作为自旋线状极限的判据。在逾渗共存区，HD团簇粗化，LD团簇碎裂。在更高压力下（T < T_g），HDA→VHDA转变通过压力-密度导数中的凹陷识别，且VHD团簇（Z≥8）出现并逾渗，与HD团簇共存。

**3.3 压力-温度相图**
对于T > T_c，冷却与压缩数据的P_∞(LD/HD)交叉点与κ_T峰值定义的Widom线（操作定义）对齐，表现为Kertész类几何逾渗线（无热力学奇点）。等压冷却/加热循环显示1 bar和1.5 kbar下LDL?HDL转变可逆，但压缩/减压循环在200 K以下出现滞后，HD状态在减压分支持续至更低压力，直至ξ′(HD)峰值指示逾渗极限。T < T_c时，压缩分支的κ_T峰值偏离假设的液-液转变线P_LL(T)，而冷却分支的P_∞(LD/HD)交叉点仍接近P_LL(T)，与Gartner III等报告的S(0)极大值轨迹一致。滞后环对称性表明压缩与减压分支分别对应LD网络失稳的自旋线状边界。

**3.4 逾渗的临界群体分数**
LD和HD团簇的临界分数?_c随配位数Z增加而减小，与3D规则晶格键逾渗阈值（Z=4时0.39，Z=6时0.25，Z=8时0.19）接近。冷却与压缩给出相似的HD逾渗阈值，但LD阈值在压缩下展宽，尤其在T ? T_c时。通过冷却/加热与压缩/减压循环，P_∞对?的依赖表现出强塌缩，仅残留弱滞后，表明逾渗转变主要取决于临界群体分数而非热力学路径。

**3.5 非晶冰与晶态冰的类比**
从构型能量景观角度，非晶态与晶态多晶型共享能量巨盆。LDA密度与冰Ih一致（~0.92 g/cm³），主要由四面体环境组成。HDA冰（密度~1.15–1.20 g/cm³）保留完全氢键四面体网络，但第二壳层塌陷，Z=5–7对应中等塌陷四面体。VHDA冰（密度~1.26 g/cm³）具有更高间隙占有率，Z≥8，与冰IV、VI、XII等致密晶型一致。

**讨论与结论总结**
讨论部分指出，T > T_c时，LD/HD团簇的逾渗与竞争导致密度涨落增强，κ_T在P_∞与?相似时达到峰值，与Kertész类几何逾渗一致。T < T_c时，压缩下κ_T峰值转向与ξ(LD)峰值对齐，转变为自旋线状行为。这一观察与刚度逾渗（rigidity percolation）理论定性关联，但水结构的异质性使体模量与密度涨落的反比关系不成立。逾渗概念已在金属玻璃和SiO₂玻璃中得到验证，暗示不同非晶材料间存在统一逾渗行为。3.5 ?截断距离的合理性在于非键相互作用调节局部密度与配位。本研究采用的氢键约束比以往更严格，但发现氢键网络在低温和中压下持续存在，而Z-团簇由非键相互作用驱动。初步有限尺寸标度分析显示，冷却数据中γ/ν=1.6±0.2，偏离Ising模型（1.97）和普通逾渗（2.05），但受限于模拟协议（有限尺寸、统计样本、弛豫时间增长），不能作为判定平衡普适类的依据。然而，P_∞交叉点与Widom线及Kertész线拓扑的对应关系支持LLCP场景。

**结论翻译**
在本工作中，研究人员采用几何和现象学逾渗模型作为统一框架，描述宽压力-温度范围内水的亚稳转变。在该框架下，低密度（LD）和高密度（HD）逾渗团簇的出现、共存和消失可解释为不同热力学条件下相分离的各个阶段。高于T_c时，逾渗转变与Widom线及水的密度和压缩率主要异常相关，系统达到亚稳平衡，与Ising类相关逾渗场景高度相似。低于T_c时，同一逾渗框架捕捉到非晶水的非平衡现象：LD关联长度最大值与压缩下的自旋线状边界吻合，滞后环与LDA–HDA转变的一级特征一致，HD–VHD逾渗转变与向甚高密度非晶冰的致密化一致。此外，从300 K等压冷却后再压缩的系统中获得的P_∞(LD/HD)交叉点轨迹与Gartner等数值报告的S(0)极大值轨迹在P_c附近吻合良好，尽管尚需更系统的有限尺寸分析。这些发现与液-液临界点（LLCP）场景一致，并提供了逾渗与亚稳相分离之间的现象学联系。然而，LLCP附近的逾渗LD或HD团簇可通过Hill型准则改进以再现Ising普适类。无论结果如何，由于逾渗发现与水的不同状态相分离内在相关，该方法也可容纳无LLCP的场景，并为理解SiO₂和金属玻璃中观察到的逾渗行为提供现象学视角。因此，比较逾渗的普适类可揭示水与其他纯物质异常之间的相似性，有望为描述非晶-非晶转变（发生于不同能量巨盆之间而非局部能量极小值）提供有用见解。