IF50.5｜准弹性中子散射直接观测到塑性冰VII

水在高压高温下存在一种分子位置固定但可快速旋转的塑性冰VII相。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

水是宇宙中第三丰富的分子，在冰巨星、天王星和海王星类系外行星内部广泛存在。由于其独特的分子结构和灵活的氢键，水在高压高温下形成多种晶态和非晶态相。在2 GPa以上，主要存在质子无序的冰VII、质子有序的冰VIII和原子相冰X，它们共享立方氧晶格。理论模拟预测沿冰VII熔化曲线存在塑性相，其中水分子占据固定晶格位置但能像液体一样旋转。塑性晶体相介于液体和晶体之间，具有独特的力学性质，已在分子固体、水混合物和离子固体中发现，但水冰中是否存在塑性相尚未得到实验证实。

科学问题

卡在哪里： 尽管分子动力学模拟预测了塑性冰VII的存在，但实验上从未直接观测到这一相。冰VII熔化曲线的连续斜率与经典模拟预测的不连续熔化相矛盾，引发了对塑性相是否存在的质疑。此外，塑性相变主要涉及氢动力学变化而非结构转变，使得实验观测极具挑战。

本文要回答： 本文旨在通过准弹性中子散射直接探测冰VII高温熔化曲线附近水分子的动力学，验证塑性冰VII是否存在，并表征其旋转机制。

技术路线

作者主要用了：准弹性中子散射、分子动力学模拟、中子衍射、X射线衍射。

作者利用巴黎-爱丁堡压机在3-6 GPa、450-600 K条件下进行准弹性中子散射实验，探测水分子的自扩散和旋转运动。通过分析散射信号随能量和波矢的变化，区分液体、塑性晶体和普通冰VII的动力学特征。同时结合分子动力学模拟，研究旋转机制和相变性质。图1展示了实验相图、散射原理和探测器原始数据，为后续分析奠定基础。

图 1　Fig. 1。图1展示了实验相图、散射原理和探测器原始数据。相图中绿色圆点标记塑性冰VII，白色虚线为模拟预测的稳定区域，与实验熔化线吻合。散射示意图说明氢的高非相干截面使信号主要来自氢的自运动。探测器图像显示高压装置限制了可用角度范围。

核心亮点

亮点 1｜塑性冰VII的直接观测

在5.5 GPa、518 K等七个热力学点，QENS信号显示准弹性展宽不随Q变化，表明水分子失去平移扩散但保留皮秒级旋转动力学，同时X射线和中子衍射证实晶体结构保持冰VII的体心立方排列。这直接证实了塑性冰VII的存在。

图 2　Fig. 2。图2对比了液体、塑性冰VII和普通冰VII在5.5 GPa下的S(Q,ω)和固定Q=1.2 Å⁻¹的QENS谱。液体显示平移和旋转两种展宽，塑性相仅显示旋转展宽且宽度不随Q变化，普通冰VII无准弹性展宽。结构快照直观展示了三种相中分子排列和动力学的差异。

亮点 2｜旋转机制为跳跃旋转而非自由转子

通过分析准弹性强度随Q的演化并与分子动力学模拟对比，发现塑性冰VII中水分子并非自由旋转，而是在晶体场限制下在优先取向间快速跳跃。马尔可夫链分析显示最频繁的旋转跃迁是绕立方轴的四重旋转（C4），旋转时间约0.5-0.7 ps，与液态水相当。

图 3　Fig. 3。图3展示了塑性冰VII在5.5 GPa、518 K下用三个洛伦兹函数拟合的QENS谱，以及各准弹性成分强度随Q的变化。拟合采用C4跳跃旋转模型，三个洛伦兹峰的宽度和相对强度受模型约束，强度因子为自由参数。强度随Q的变化与模型预测一致，支持C4旋转机制。

亮点 3｜氢原子振动幅度显著增大

拟合得到塑性相中氢的平均振动幅度⟨u²⟩约为0.80±0.1 Ų，比普通冰VII高约四倍，与TIP4P/2005模拟结果一致。这表明旋转与晶格畸变耦合，解释了旋转时间对温度和压力的弱依赖性。

图 4　Fig. 4。图4展示了分子动力学模拟得到的氢原子位置概率分布、旋转中间散射函数和密度随温度的变化。塑性相中氢原子概率分布比普通冰VII更弥散，但仍保留优先取向，证实非各向同性旋转。旋转中间散射函数的拟合给出旋转时间。密度在塑性-普通相变处出现约2%的跳变，提示一级相变。

亮点 4｜塑性相变可能为一级相变

分子动力学模拟显示从塑性冰VII到普通冰VII的密度变化对应约2%的体积变化，提示一级相变。实验室X射线衍射在5.7 GPa、约500 K观察到体积跳变，与模拟一致，但单晶形成阻碍了进一步确认。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 准弹性中子散射、分子动力学模拟、中子衍射、X射线衍射
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事水相变、行星内部物理或中子散射研究的科研人员。

带走一句

塑性冰VII的发现提示在高压高温下，氢键网络可以允许分子旋转而保持晶格完整，这种跳跃旋转机制可能在其他氢键材料中普遍存在。