[論文レビュー] 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 ポテンシャルを含む8種類の半経験的水モデルを用いて、分子動力学シミュレーションにより氷 I h の弾性定数を計算した。mW モデルは定量的に最も正確な弾性剛性値を提供するが、特に弾性異方性という重要な定性的な特徴を再現できないことから、水素結合の分子的性質が氷 I h の正確な機械的挙動に不可欠であることが示された。

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$.

研究の動機と目的

• さまざまな半経験的水ポテンシャルを用いた分子動力学シミュレーションにより、氷 I h の弾性定数を計算すること。
• 明示的プロトンモデルおよび粗粒度モデルを含む、さまざまな水モデルの実験的弾性特性を再現する精度を評価すること。
• 水素結合や回転自由度などの分子的詳細が、機械的異方性を決定する役割を調査すること。
• 計算効率に優れるものの、mW 粗粒度モデルが氷 I h の機械的応答を適切に捉えられるかどうかを評価すること。
• 剛体および可動性の明示的プロトンモデルが、定量的および定性的な弾性挙動をどれほど正確に再現できるかを比較すること。

提案手法

• 制御されたひずみ変形における応力応答を計算することで弾性剛性定数を特定するため、分子動力学シミュレーションを用いた。
• 計算の容易さを図るため、ランク4の弾性剛性テンソル $C_{ijkl}$ を $6\times6$ 行列としてボイチェット表記で表現した。
• 定温定圧条件下で有限変形技術を適用し、応力-ひずみ応答から弾性定数を抽出した。
• SPC/E、TIP4P、TIP4P2005、TIP4P/Ice、TIP5P、qTIP4P/Fw、SPC/Fw、mW のすべてのモデルを、さまざまな熱力学的条件下でシミュレートした。
• 定量的精度および弾性異方性などの定性的傾向を含め、実験データと結果を比較した。
• 特に mW モデルにおける異常を示す傾向に注目し、弾性定数の温度および圧力依存性を分析した。

実験結果

リサーチクエスチョン

• RQ1 異なる水モデルを用いて計算された氷 I h の弾性定数は、実験値とどのように比較できるか？
• RQ2 mW 粗粒度モデルは、簡略化された表現にもかかわらず、なぜ実験値とより良い定量的一致を示すのか？
• RQ3 明示的プロトンモデルは、特にヤング率およびせん断率において、氷 I h で観察される弾性異方性をどの程度再現できるか？
• RQ4 mW モデルにおいて温度上昇に伴い $C_{12}$ および $C_{13}$ が異常的に増加する原因は何か？
• RQ5 mW モデルが分子的自由度を欠いている場合、氷 I h の非中心力および角度依存性を正しく記述できるか？

主な発見

• mW モデルは、特にゼロ温度・圧力条件下で、実験値に最も近い弾性剛性定数を生成した。
• すべての明示的プロトンモデル（剛体・可動性を問わず）は、実験値から顕著な定量的ずれを示し、特に TIP5P は剛性を約2倍に過剰評価した。
• mW モデルはヤング率およびせん断率の異方性を再現できず、実験的に観察されたよりもより等方的な機械的応答を予測した。
• 明示的プロトンモデルは一貫して弾性異方性の定性的な傾向を再現しており、この挙動が氷 I h の水素結合の分子的性質に深く関連していることが示された。
• mW モデルは温度上昇に伴い $C_{12}$ および $C_{13}$ が異常的に増加する傾向を示し、ポテンシャルが角度依存性および非中心力の正確な記述に失敗している可能性を示唆した。
• mW モデルが機械的異方性を記述できず、物理的に不自然な温度依存性を示すことは、分子的詳細（特に回転自由度および非対称水素結合）を粗粒度化することで、氷 I h の機械的応答をモデル化する上で顕著な制限が生じることを強調している。