[논문 리뷰] Overcoming timescale and finite-size limitations to compute condensation rates at physically relevant conditions

Matteo Salvalaglio, Pratyush Tiwary외 2명|arXiv (Cornell University)|2015. 08. 07.

nanoparticles nucleation surface interactions인용 수 3

한 줄 요약

이 논문은 아르곤 응축의 분자 동역학 시뮬레이션에서 시간스케일 및 유한한 크기의 제약을 극복하는 개선된 샘플링 방법을 제시한다. 이는 물리적으로 관련 있는 과포화 수준에서 안정적인 핵형성 비율 계산을 가능하게 하며, 시간 단위로 수시간에 이르는 응축 시간스케일을 달성함으로써 표준 시뮬레이션과 실제 조건 사이의 격차를 메운다.

ABSTRACT

Condensation of a liquid droplet from a supersaturated vapour phase is a classical nucleation phenomenon. As such it is challenging to compute its rate from atomistic molecular dynamics simulations. In fact at realistic supersaturation conditions condensation occurs on time scales that far exceed what can be reached by conventional molecular dynamics methods. Another known problem in this context is the distortion of the free energy profile associated to nucleation due to the small, finite size of typical simulation boxes. In this work the problem of time scale is addressed with a recently developed enhanced sampling method while contextually correcting for finite size effects. We demonstrate our approach studying the condensation of argon, and showing that characteristic nucleation times of the order of magnitude of hours can be reliably calculated, approaching realistic supersaturation conditions. We thus bridge the gap between what standard molecular dynamics can do and real physical systems.

연구 동기 및 목표

• 실제 과포화 조건에서 초래되는 너무 긴 시간스케일로 인해 원자 수준의 분자 동역학에서 응축 속도를 시뮬레이션하는 데 도전하는 문제를 해결한다.
• 작은 시뮬레이션 박스 크기로 인해 발생하는 자유 에너지 프로파일의 왜곡을 보정한다.
• 물리적으로 관련 있는 핵형성 비율 계산을 가능하게 하여 실험적 시간스케일(수시간)에 가까운 결과를 도출한다.
• 실제 과포화 조건에서 모델 시스템으로 아르곤을 사용하여 방법의 효과성을 입증한다.

제안 방법

• 최근 개발된 강화 샘플링 방법을 활용하여 분자 동역학에서 희귀 핵형성 사건의 샘플링을 가속화한다.
• 작은 시뮬레이션 박스에 내재된 유한한 크기 효과를 고려해 자유 에너지 프로파일에 보정을 적용한다.
• 강화 샘플링과 유한한 크기 보정을 결합하여 정확한 핵형성 비율 계산을 달성한다.
• 실제 물리 시스템에서 관찰되는 조건에 가까운 조건에서 아르곤 응축을 시뮬레이션한다.
• 모의된 핵형성 시간스케일을 실험적으로 관련 있는 값과 비교하여 접근법의 타당성을 검증한다.

실험 결과

연구 질문

• RQ1 강화 샘플링 방법은 기존 분자 동역학의 시간스케일 제약을 극복하여 응축을 시뮬레이션할 수 있는가?
• RQ2 작은 시뮬레이션 박스에서의 유한한 크기 효과는 핵형성의 자유 에너지 프로파일을 어떻게 왜곡하는가?
• RQ3 유한한 크기 보정은 얼마나 정확도 향상에 기여하는가?
• RQ4 결합된 접근법을 통해 실제 과포화 수준에서 실험 관측치(예: 수시간)에 가까운 핵형성 시간스케일을 도출할 수 있는가?
• RQ5 이 방법을 사용해 물리적 과포화 조건에서 아르곤 응축을 충분한 정확도로 계산할 수 있는가?

주요 결과

• 이 방법은 물리적으로 관련 있는 과포화 수준에서 아르곤 응축의 핵형성 비율을 성공적으로 계산하였다.
• 수시간 수준의 응축 시간스케일이 안정적으로 계산되어 실험 조건에 가까워졌다.
• 자유 에너지 프로파일에 미치는 유한한 크기 효과가 효과적으로 보정되어 정확도가 향상되었다.
• 강화 샘플링 접근법을 통해 기존 표준 분자 동역학로는 도달할 수 없는 희귀 핵형성 사건에 접근할 수 있었다.
• 결과적으로 원자 수준의 시뮬레이션과 실제 물리 시스템 간의 핵형성 현상에서의 격차를 크게 줄였다.