冰 Ih 的路径积分模拟：压力的影响

通过路径积分分子动力学模拟，在 50 至 300 K 的温度下研究了压力对冰 Ih 的结构和热力学性质的影响。原子间相互作用采用柔性水的有效 q-TIP4P/F 势进行建模。考虑了正压（压缩）和负压（拉伸），这使我们能够接近该固态水相机械稳定性的极限。我们研究了晶体体积、体积模量、原子间距、原子离域和动能的压力依赖性。从体积模量消失推导出正压和负压下的旋节线点。对于 $P < 0$，在 50 至 300 K 范围内，旋节线压力从 -1.38 GPa 变化到 -0.73 GPa。在正压下，旋节线与冰的非晶化相关，在低温下位于 1.1 至 1.3 GPa 之间。量子核效应导致冰 Ih 的亚稳区域减小。