[论文解读] Elastic constants of Ice I$_h$ as described by semi-empirical water models

P.A.F.P. Moreira, Roberto Gomes de Aguiar Veiga+1 more|arXiv (Cornell University)|Nov 23, 2018

nanoparticles nucleation surface interactions被引用 6

一句话总结

本研究通过分子动力学模拟，对八种半经验水模型（包括刚性和柔性显式质子模型以及粗粒化mW势）的冰Ih弹性常数进行了计算。尽管mW模型在弹性刚度值的定量准确性方面表现最佳，但其无法再现关键的定性特征——尤其是弹性各向异性——表明氢键的分子特性对于准确描述冰Ih的力学行为至关重要。

ABSTRACT

Using molecular dynamics simulations we compute the elastic constants of ice I$_h$ for a set of 8 frequently used semi-empirical potentials for water, namely the rigid-molecule SPC/E, TIP4P, TIP4P2005, TIP4P/Ice and TIP5P models, the flexible-molecule qTIP4P/Fw and SPC/Fw models and the coarse-grained atomic mW potential. In quantitative terms, the mW description gives values for the individual stiffness constants that are closest to experiment, whereas the explicit-proton models display substantial discrepancies. On the other hand, in contrast to all explicit-proton potentials, the mW model is unable to reproduce central qualitative trends such as the anisotropy in Young's modulus and the shear modulus. This suggests that the elastic behavior of ice I$_h$ is closely related to its molecular nature, which has been coarse-grained out in the mW model. These observations are consistent with other recent manifestations concerning the limitations of the mW model in the description of mechanical properties of ice I$_h$.

研究动机与目标

• 通过一系列半经验水势的分子动力学模拟，计算冰Ih的弹性常数。
• 评估不同水模型（显式质子与粗粒化模型）在再现实验弹性性质方面的准确性。
• 研究分子细节（如氢键和转动自由度）在决定机械各向异性中的作用。
• 评估mW粗粒化模型尽管计算效率高，是否能充分捕捉冰Ih的机械响应。
• 比较刚性和柔性显式质子模型在再现定量与定性弹性行为方面的表现。

提出的方法

• 采用分子动力学模拟，通过控制应变变形下的应力响应来确定弹性刚度常数。
• 使用Voigt记法将四阶弹性刚度张量 $C_{ijkl}$ 表示为 $6\times6$ 矩阵，以提高计算可行性。
• 在恒温恒压条件下应用有限变形技术，从应力-应变响应中提取弹性常数。
• 在不同热力学条件下模拟所有模型（包括SPC/E、TIP4P、TIP4P2005、TIP4P/Ice、TIP5P、qTIP4P/Fw、SPC/Fw和mW），以评估其行为。
• 将结果与实验数据对比，评估定量准确性及弹性各向异性等定性趋势。
• 分析弹性常数随温度和压力的变化，特别关注mW模型中的异常现象。

实验结果

研究问题

• RQ1 不同水模型计算的冰Ih弹性常数与实验值相比如何？
• RQ2 为何mW粗粒化模型尽管表示简化，却能与实验结果实现更好的定量吻合？
• RQ3 显式质子模型在多大程度上再现了冰Ih中观察到的弹性各向异性，特别是杨氏模量和剪切模量？
• RQ4 mW模型中$C_{12}$和$C_{13}$随温度升高而异常增强的原因是什么？
• RQ5 mW模型中缺乏分子自由度如何影响其对冰Ih中非中心力和角度依赖性的描述能力？

主要发现

• mW模型产生的弹性刚度常数在定量上最接近实验值，尤其是在零温零压条件下。
• 所有显式质子模型——无论刚性或柔性——在定量上均表现出显著偏差，其中TIP5P的刚度高估约两倍。
• mW模型无法再现杨氏模量和剪切模量的各向异性，预测的机械响应比实验观测更趋各向同性。
• 显式质子模型始终能捕捉弹性各向异性的定性趋势，表明该行为与冰Ih中氢键的分子本质密切相关。
• mW模型表现出$C_{12}$和$C_{13}$随温度升高而异常增大的现象，表明其在势函数中未能正确描述角度依赖性和非中心力。
• mW模型无法描述机械各向异性和其非物理的温度依赖性，凸显了粗粒化过程中忽略分子细节（尤其是转动自由度和不对称氢键）会显著限制对冰Ih机械响应的建模能力。