[论文解读] Mass-density and Phonon-frequency Relaxation Dynamics of Water and Ice at Cooling

Changqing Sun|arXiv (Cornell University)|Oct 5, 2012

Material Dynamics and Properties参考文献 19被引用 78

一句话总结

本文提出，水分子中 O:H-O 氢键内成键电子对与非成键电子对之间的库仑排斥力，结合 O:H 与 H-O 段之间比热差异，共同主导了冷却过程中异常的密度变化与声子频率弛豫动力学。O:H 键在液态和固态中均充当‘主导’角色，冷却时收缩并变刚，而 H-O 键作为‘从属’角色则伸长并变软；这种在水-冰相变过程中动态角色的反转，解释了冰的低密度及漂浮现象。

ABSTRACT

Coulomb repulsion between the bonding electron pair in the H-O covalent bond (denoted by "-") and the nonbonding electron pair of O (":") and the specific-heat disparity between the O:H and the H-O segments of the entire hydrogen bond (O:H-O) are shown to determine the O:H-O bond angle-length-stiffness relaxation dynamics and the density anomalies of water and ice. The bonding part with relatively lower specific-heat is more easily activated by cooling, which serves as the "master" and contracts, while forcing the "slave" with higher specific-heat to elongate (via Coulomb repulsion) by different amounts. In the liquid and solid phases, the O:H van der Waals bond serves as the master and becomes significantly shorter and stiffer while the H-O bond becomes slightly longer and softer (phonon frequency is a measure of bond stiffness), resulting in an O:H-O cooling contraction and the seemingly "regular" process of cooling densification. In the water-ice transition phase, the master and the slave swap roles, thus resulting in an O:H-O elongation and volume expansion during freezing. In ice, the O--O distance is longer than it is in water, resulting in a lower density, so that ice floats.

研究动机与目标

• 解决水在 277 K 出现密度极大值以及冰在冷却后密度更低的长期谜题。
• 解释液态与固态水中看似规律的冷却致密化现象，与冻结过程中的异常膨胀形成对比。
• 阐明氢键角度-长度-刚度弛豫在连接声子频率位移与宏观密度异常之间的作用。
• 基于电子对库仑排斥力与 O:H-O 键段特异性比热容差异，建立统一理论框架。
• 将测量得到的拉曼与红外声子位移与水及冰所有相态中的结构弛豫动力学相关联。

提出的方法

• 提出氢键的双段模型：O:H（范德华）与 H-O（共价），二者具有不同的比热响应。
• 应用德拜模型估算各段特异性比热容（η_L 表示 O:H，η_H 表示 H-O），表明 η_L 随温度升高比 η_H 更快。
• 利用分子动力学（MD）模拟，研究从 300 K 冷却至 4 K 时 O:H-O 键的弛豫动力学。
• 引入成键（H-O）与非成键（O:）电子对之间的库仑排斥力（f_q）作为 O:H-O 氢键非对称弛豫的驱动力。
• 通过拉曼光谱数据与 MD 预测，将声子频率位移（ω_L 表示 O:H，ω_H 表示 H-O）与键长及刚度变化相关联。
• 分析三个温度区域：T > 273 K（I），273–258 K（II），T < 258 K（III），识别出不同的弛豫行为。

实验结果

研究问题

• RQ1 为何水在 277 K 出现密度极大值，且在冻结时膨胀，与大多数物质相反？
• RQ2 在液态与固态水中冷却过程中，H-O 声子频率（ω_H）为何出现红移，而冻结时又出现蓝移？
• RQ3 O:H 与 H-O 段之间比热差异如何影响键长与刚度弛豫的方向与大小？
• RQ4 成键与非成键电子对之间的库仑排斥力在驱动 O:H-O 氢键非对称弛豫中起何作用？
• RQ5 为何在水-冰相变过程中，O:H 与 H-O 键的主从角色发生反转，导致体积膨胀？

主要发现

• O:H 键（范德华）在液态与固态中均充当‘主导’角色，冷却时收缩并变刚，其收缩程度超过 H-O 键的伸长，导致冷却致密化。
• H-O 键（共价）充当‘从属’角色，因受 O:H 键的库仑排斥而伸长并变软，导致 ω_H 在液态中从约 3200 cm⁻¹ 红移到约 3150 cm⁻¹。
• 在冻结过程中（273–258 K），主从角色发生反转：O:H 键伸长并变软（ω_L 从 220 cm⁻¹ 红移到 215 cm⁻¹），而 H-O 键则收缩并变刚（ω_H 从 3150 cm⁻¹ 蓝移到 3170 cm⁻¹），导致体积膨胀。
• 冰中 O–O 距离大于液态水，导致密度降低 9%，从而解释冰为何能漂浮。
• 冰中 ω_H 声子频率从 209 K 时的 3253 cm⁻¹ 降至 30 K 时的 3218 cm⁻¹，与纳米尺度冰簇的红外光谱一致。
• MD 模拟动画证实了氢键的动态切换与液态水中持续存在的库仑排斥，支持了具有热涨落的四面体配位结构的存在。