[논문 리뷰] Molecular dynamic simulation of water vapor interaction with blind pore of dead-end and saccate type

E. G. Nikonov, Miron Pavluš외 1명|arXiv (Cornell University)|2017. 07. 28.

Nanopore and Nanochannel Transport Studies인용 수 4

한 줄 요약

이 연구는 분자 동역학(MD) 시뮬레이션과 매크로스코픽 확산 모델을 결합하여 사선형과 사선형이 아닌 유형의 막힘 포어에서 수증기 이동을 연구한다. MD 시뮬레이션에서 유도된 확산 계수를 3차원 확산 방정식에 적용함으로써, 저자들은 하이브리드 접근법이 건조 및 습윤 동역학을 정확하게 캡처함을 입증하며, 수증기 농도의 시간에 따른 마이크로 및 매크로 스케일 예측 간 강한 일치를 보였다.

ABSTRACT

One of the varieties of pores, often found in natural or artificial building materials, are the so-called blind pores of dead-end or saccate type. Three-dimensional model of such kind of pore has been developed in this work. This model has been used for simulation of water vapor interaction with individual pore by molecular dynamics in combination with the diffusion equation method. Special investigations have been done to find dependencies between thermostats implementations and conservation of thermodynamic and statistical values of water vapor - pore system. The two types of evolution of water-pore system have been investigated: drying and wetting of the pore. Full research of diffusion coefficient, diffusion velocity and other diffusion parameters has been made.

연구 동기 및 목표

• 분자 동역학를 사용하여 사선형 및 사선형이 아닌 유형의 막힘 포어에서 수증기 상호작용을 모델링하기.
• 수증기-포어 시스템에서 열역학적 안정성에 영향을 미치는 열제어 매개변수의 영향을 조사하기.
• 건조 및 습윤 과정 중 확산 계수와 속도를 결정하기.
• MD에서 유도된 확산 계수를 사용하여 매크로스코픽 확산 모델을 개발하고 검증하기.
• 수증기 농도에 대한 수치적 예측을 비교하기 위해 마이크로 스케일 MD 결과와 매크로 스케일 확산 모델 예측을 비교하기.

제안 방법

• MD 시뮬레이션을 위한 3차원 막힘 포어 모델(500×50×50 nm)을 구축하였다.
• 물 분자 상호작용을 기술하기 위해 σ=3.17 Å 및 ε=6.74×10⁻³ eV를 사용한 레너드존 포텐셜을 사용하였다.
• 속도 버틀 방법과 베렌센 열제어(τB ∈ [0.5,2] ps)를 사용하여 25 °C에서 온도를 유지하였다.
• 65.3 ns 동안 시뮬레이션을 수행하였으며, 시간 간격 Δt=0.016 ps로, 포화 내 1000–1800개의 H₂O 분자와 외부 영역 내 1800–9000개의 분자를 추적하였다.
• 건조 및 습윤 과정에서의 MD 데이터로부터 확산 계수를 계산하였다.
• 포어 개구부에서의 수증기 이동을 고려한 경계 조건을 포함한 3차원 매크로스코픽 확산 모델을 수립하였다. 방정식은 ∂wv/∂t = D(∂²wv/∂x² + ∂²wv/∂y² + ∂²wv/∂z²)로 표현된다.

실험 결과

연구 질문

• RQ1 다양한 열제어 구현 방식이 수증기-포어 시스템에서 열역학적 및 통계적 성질의 보존에 미치는 영향은 무엇인가?
• RQ2 막힘 포어에서 건조 및 습윤 과정 동안 시간에 따라 변화하는 수증기의 확산 계수는 무엇인가?
• RQ3 MD에서 유도된 매크로스코픽 확산 모델의 예측과 마이크로 스케일 MD 결과는 어떻게 비교되는가?
• RQ4 포어 기하학적 형상과 경계 조건이 수증기 이동 동역학에 미치는 역할은 무엇인가?
• RQ5 하이브리드 MD–확산 모델은 막힘 포어에서 수증기 농도의 임시 프로파일을 정확하게 기술할 수 있는가?

주요 결과

• MD 시뮬레이션 결과, 베렌센 열제어가 65.3 ns 동안 시스템을 25 °C에서 안정적으로 유지하며 열역학적 성질을 잘 보존하는 것으로 나타났다.
• 건조 과정의 확산 계수는 약 1779.1 nm²/ps였고, 습윤 과정에서는 536.33 nm²/ps였다.
• 매크로스코픽 확산 모델은 변수 분리 방법을 통해 시간에 따른 수증기 농도에 대한 정확한 해석적 해를 제공하였다.
• 매크로 모델의 공간 평균 수증기 농도는 Fig. 3에서 확인한 바와 같이 MD 시뮬레이션 결과와 강한 일치를 보였다.
• MD에서 유도된 D 값을 사용한 매크로스코픽 확산 모델에 기반한 하이브리드 접근법은 수증기-포어 상호작용을 시뮬레이션하는 데 신뢰할 수 있는 방법으로 검증되었다.
• 모델은 농도가 처음에 높은 건조 과정과 처음에 낮은 습윤 과정의 동역학을 성공적으로 캡처하였으며, 시간 척도에 걸쳐 농도 프로파일이 일치하였다.