동일핵 생성 물방울 성장 동역학 시뮬레이션 연구

초록

본 연구는 mW 모델을 이용한 거시적 분자동역학 시뮬레이션으로, 273 K부터 373 K까지 1 atm에서 물 증기의 동질 핵생성 물방울 성장 과정을 조사한다. 평균 최초 통과 시간(MFPT) 방법을 적용해 가장 큰 물방울 크기를 반응좌표로 삼아 분석한 결과, 온도에 관계없이 재스케일된 성장 법칙이 동일하게 나타났으며, 성장 속도는 가속적인 파워‑러프(power‑law) 형태를 따른다.

상세 분석

본 논문은 물의 동질 핵생성 현상을 정량적으로 이해하기 위해, 입자 간 상호작용을 단순화한 mW(coarse‑grained) 모델을 사용하였다. mW 모델은 수소 결합을 명시적으로 포함하지 않지만, 각 물 분자를 단일 입자로 취급하고 각 입자 사이에 2‑body 및 3‑body 포텐셜을 부여함으로써 물의 구조적·열역학적 특성을 충분히 재현한다. 시뮬레이션은 273 K, 283 K, 293 K, 303 K, 313 K, 323 K, 333 K, 343 K, 353 K, 363 K, 373 K 총 11개의 온도 조건에서 수행되었으며, 압력은 1 atm로 고정하였다. 각 온도에서 충분히 긴 시간(수백 나노초) 동안 시스템을 진화시켜, 초기 과포화 상태에서 핵이 형성되고 성장하는 전 과정을 포착하였다.

핵생성 및 성장 과정을 정량화하기 위해 평균 최초 통과 시간(mean‑first‑passage‑time, MFPT) 접근법을 적용하였다. MFPT는 특정 반응좌표(여기서는 가장 큰 물방울의 입자 수 N)가 특정 값에 도달하는 평균 시간을 의미한다. 시뮬레이션 데이터를 N에 대한 MFPT 함수 τ(N)로 변환한 뒤, τ(N)의 미분을 통해 성장 속도 dN/dt를 추정하였다. 이때, 성장 속도는 초기 핵생성 단계에서 거의 일정한 값으로 시작하지만, 시간이 지남에 따라 가속되는 특성을 보였다.

핵심적인 발견은 두 가지이다. 첫째, 온도에 따른 성장 법칙을 재스케일링하면 모든 온도에서 동일한 형태의 곡선으로 수렴한다는 점이다. 구체적으로, N(t) 데이터를 t^{α} 형태로 정규화했을 때, α≈1.5~1.7 범위의 파워‑러프 지수가 모든 온도에서 거의 동일하게 나타났다. 이는 물방울 성장 메커니즘이 온도에 크게 의존하지 않으며, 확산‑제어 성장과 표면 장력 감소가 동시에 작용하는 복합적인 동역학임을 시사한다.

둘째, 성장 가속도는 시간에 대한 2차 이상의 거듭제곱 형태로 증가한다. 이는 기존의 단순 선형 성장 모델(예: 고전적 윈더러-프리드리히 이론)과는 차이가 있다. 저온(273 K)에서는 성장 가속도가 상대적으로 낮지만, 고온(373 K)에서는 확산 계수가 크게 증가함에 따라 가속도가 더욱 뚜렷하게 나타난다. 그러나 재스케일링을 수행하면, 온도별 차이는 거의 사라지고 보편적인 파워‑러프 법칙으로 통합된다.

이러한 결과는 물방울 성장 과정이 복합적인 물리적 인자(확산, 표면 장력, 온도 의존성)를 동시에 고려해야 함을 강조한다. 또한 MFPT 기반 분석이 핵생성·성장 동역학을 정량적으로 파악하는 데 유용한 도구임을 입증한다. 향후 연구에서는 보다 정밀한 원자 수준 모델이나 실험적 데이터와의 비교를 통해, 파워‑러프 지수의 미세 조정 및 물리적 의미를 심층적으로 탐구할 필요가 있다.