[논문 리뷰] The Boson peak in supercooled water

Pradeep Kumar, Kjartan Thor Wikfeldt외 3명|arXiv (Cornell University)|2013. 05. 19.

Material Dynamics and Properties참고 문헌 54인용 수 4

한 줄 요약

이 연구는 TIP4P/2005 수분 모델의 분자 동역학 시뮬레이션을 통해 초냉각된 물에서 보존 피크가 위도 막선 온도 $T_W \approx 230$ K 이하에서 나타나며, 이는 저밀도 유사 액체로의 구조 전이와 동시에 발생함을 보여준다. 피크는 원형 얼음 Ih의 횡파 진행성 진동 모드에서 기인하며, 이는 보존 피크가 구조적 상전이를 겪는 테트라헤드랄 액체의 일반적 특성임을 시사하며, 국소화에 국한된 것이 아님을 나타낸다.

ABSTRACT

We perform extensive molecular dynamics simulations of the TIP4P/2005 model of water to investigate the origin of the Boson peak reported in experiments on supercooled water in nanoconfined pores, and in hydration water around proteins. We find that the onset of the Boson peak in supercooled bulk water coincides with the crossover to a predominantly low-density-like liquid below the Widom line $T_W$. The frequency and onset temperature of the Boson peak in our simulations of bulk water agree well with the results from experiments on nanoconfined water. Our results suggest that the Boson peak in water is not an exclusive effect of confinement. We further find that, similar to other glass-forming liquids, the vibrational modes corresponding to the Boson peak are spatially extended and are related to transverse phonons found in the parent crystal, here ice Ih.

연구 동기 및 목표

• 나노구속 조건과 수화층에서 실험적으로 관측된 초냉각된 물의 보존 피크 기원을 조사하기 위해
• 초냉각된 물의 보존 피크가 구속 효과에 기인한 잡음인지 아니면 구조적 및 동적 전이와 관련된 막대한 성질인지 판단하기 위해
• 보존 피크와 원형 결정인 얼음 Ih의 진동 모드 간의 관계를 탐색하기 위해
• 보존 피크가 액체-액체 상전이와 위도 막선과 관련이 있는지 테스트하기 위해

제안 방법

• 대기압 하에서 210 K에서 260 K 사이의 온도 범위에서 TIP4P/2005 수분 모델의 분자 동역학 시뮬레이션을 수행하며, 약 230 K에서 위도 막선을 포함한다.
• 균형 상태 확보를 위해 NPT 및 NVT 캐노니컬 집합을 사용하고, 동적 성질 계산을 위해 NVE 집합을 사용하며, 시간 단계는 0.2–1.0 fs 범위로 설정한다.
• 비탄성 동적 구조 인자 $S_S(k,\omega)$ 및 축소된 진동 밀도 상태(VDOS)를 계산하여 보존 피크를 식별한다.
• 횡파 및 猶량파 상관 함수 분석을 통해 보존 피크 강도에 기여하는 진동 모드를 평가한다.
• 균형 상태 구조를 냉각하여 고유 구조로 전환하고, 헤시안 행렬의 대각화를 통해 진동 모드 및 참여 비율을 계산한다.
• TIP5P 수분 모델을 사용하여 240–270 K 범위에서 결과의 타당성을 검증함으로써 다양한 수분 모델 간 결과의 강건성을 확인한다.

실험 결과

연구 질문

• RQ1 초냉각된 물의 보존 피크가 구속 효과에서 기인하는지 아니면 내재된 막대한 성질에서 기인하는가?
• RQ2 보존 피크의 발생 시점이 위도 막선과 저밀도 유사 액체로의 구조 전이와 관련이 있는가?
• RQ3 보존 피크를 일으키는 진동 모드가 원형 얼음 Ih 결정의 횡파 진행성 진동 모드와 관련이 있는가?
• RQ4 막대한 초냉각된 물에서의 보존 피크 주파수는 나노구속된 물에서의 실험 관측치와 일치하는가?
• RQ5 물의 보존 피크는 액체-액체 상전이를 겪는 테트라헤드랄 액체의 일반적 특성인가?

주요 결과

• 초냉각 막대한 물에서의 보존 피크는 정확히 위도 막선 온도 $T_W \approx 230$ K에서 나타나며, 이는 저밀도 유사 액체로의 구조 전이와 동시에 발생한다.
• TIP4P/2005 수분 모델에서 저온 영역에서의 보존 피크 주파수는 $\omega \approx 45$ cm$^{-1}$에 도달하며, 나노구속된 물의 실험 값과 일치한다.
• 보존 피크는 원형 얼음 Ih 결정의 횡파 진행성 진동 모드에 의해 주로 유도되며, 상관 함수 분석 및 모드 분해를 통해 확인되었다.
• LDA 유리 상에서 피크 주파수는 약 $45$ cm$^{-1}$에 가까워지며, 삼각형 구조의 얼음에서의 횡파 진행성 피크와 실험 및 시뮬레이션 결과와 일치한다.
• 계산 크기의 영향은 미미하며, 최대 $N = 45,000$ 개 분자로 수행된 시뮬레이션에서도 시스템 크기 변화에 관계없이 일관된 결과를 보였다.
• TIP5P 수분 모델 역시 그 위도 막선 이하에서 보존 피크를 나타내어, 이 현상이 수분 모델에 의존하지 않고 일반적인 성질임을 확인하였다.