[논문 리뷰] Elastic constants of Ice I$_h$ as described by semi-empirical water models

P.A.F.P. Moreira, Roberto Gomes de Aguiar Veiga외 1명|arXiv (Cornell University)|2018. 11. 23.

nanoparticles nucleation surface interactions인용 수 6

한 줄 요약

이 연구는 강체 및 유연한 직접 수소 원자 모델과 조절된 입자(mW) 잠재력이 포함된 여덟 가지 반경험적 물질 모델을 사용하여 분자 동역학 시뮬레이션을 통해 얼음 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$ 행렬로 표현하기 위해 Voigt 표기법을 사용하였다.
• 등온 및 등압 조건 하에서 유한 변형 기법을 적용하여 응력-변형률 반응에서 탄성 상수를 추출하였다.
• SPC/E, TIP4P, TIP4P2005, TIP4P/Ice, TIP5P, qTIP4P/Fw, SPC/Fw, mW를 포함한 모든 모델을 다양한 열역학 조건에서 시뮬레이션하였다.
• 정량적 정확도와 탄성 이방성과 같은 정성적 경향을 실험 데이터와 비교하였다.
• 특히 mW 모델에서의 이례적 행동을 중심으로 탄성 상수의 온도 및 압력 의존성을 분석하였다.

실험 결과

연구 질문

• RQ1 다양한 물질 모델을 사용해 계산한 얼음 I h의 탄성 상수는 실험 값과 어떻게 비교되는가?
• RQ2 단순화된 표현에도 불구하고 mW 조절된 입자 모델이 실험 결과와 더 나은 정량적 일치를 보이는 이유는 무엇인가?
• RQ3 직접 수소 원자 모델은 얼음 I h에서 관찰된 탄성 이방성, 특히 영우 모odulus와 전단 모odulus에서 어느 정도 잘 재현하는가?
• RQ4 mW 모델에서 온도가 상승함에 따라 $C_{12}$와 $C_{13}$이 이례적으로 강성 증가를 보이는 원인는 무엇인가?
• RQ5 mW 모델에서 분자 자유도가 부재함으로써 얼음 I h의 비중앙 힘과 각도 의존성의 기술 능력에 어떤 영향을 미치는가?

주요 결과

• mW 모델은 영온도 및 영압 조건에서 실험 값에 가장 가까운 정량적 탄성 스티프니스 상수를 생성한다.
• 모든 직접 수소 원자 모델—강체 또는 유연한지에 관계없이—실험 값에서 상당한 정량적 편차를 보이며, 특히 TIP5P는 강성을 약 2배 과대 평가한다.
• mW 모델은 영우 모odulus와 전단 모odulus의 이방성을 재현하지 못하여 실험적으로 관찰된 것보다 더 등방성에 가까운 기계적 반응을 예측한다.
• 직접 수소 원자 모델은 일관되게 탄성 이방성의 정성적 경향을 반영하며, 이는 이러한 거동이 얼음 I h의 수소 결합의 분자적 성격과 깊이 연관되어 있음을 시사한다.
• mW 모델은 온도 상승에 따라 $C_{12}$와 $C_{13}$이 이례적으로 증가하는 경향을 보이며, 이는 잠재력에서 각도 및 비중앙 힘을 적절히 기술하지 못하고 있음을 시사한다.
• mW 모델이 기계적 이방성을 기술하지 못하고 온도 의존성이 물리적으로 타당치 않은 점은, 분자적 세부 사항—특히 회전 자유도와 비대칭 수소 결합—을 간소화함으로써 얼음 I h의 기계적 반응을 모델링하는 데 심각한 제약가 존재함을 강조한다.