물방울 핵생성 및 성장의 정상상태 거동: 확장 수치 해석

초록

본 연구는 mW 모델을 이용한 거시적 분자동역학 시뮬레이션으로 273 K–363 K, 1 atm 조건에서 물의 균일한 기상‑액상 핵생성과 그 후의 액적 성장 과정을 조사한다. 열역학 적분과 확장된 첫 통과 시간(MFPT) 방법을 적용해 표면 장력, 핵생성 속도, 임계 클러스터 크기, 핵생성 장벽, 젤도프 계수를 추정하고, 고전 핵생성 이론과 비교한다. 결과는 성장 속도가 비정상적이며 파워‑법칙을 따르고, 온도에 관계없이 일관된 성장 지수를 보임을 보여준다.

상세 분석

본 논문은 물의 균일한 기상‑액상 핵생성 과정을 미세구조 수준에서 정량화하기 위해, 기존 Molinero‑Moore가 제시한 mW(Coarse‑grained) 모델을 활용하였다. mW 모델은 수소 결합을 명시적으로 포함하지 않지만, 각 입자 간 3‑body 포텐셜을 통해 물의 구조적·열역학적 특성을 재현한다는 장점이 있다. 저자들은 273 K에서 363 K까지 10 K 간격으로 시뮬레이션을 수행했으며, 압력은 1 atm에 고정하였다. 시뮬레이션 박스는 약 10 nm 크기의 주기적 경계 조건을 적용했으며, 입자 수는 10⁴ ~ 10⁵개 수준으로 설정해 충분한 통계적 표본을 확보하였다.

핵생성 분석에 핵심적인 두 가지 수치 기법을 도입했다. 첫째, 열역학 적분(s thermodynamic integration) 방식을 통해 클러스터 표면 장력을 직접 계산하였다. 이는 클러스터 크기가 변함에 따라 자유 에너지 곡면을 적분함으로써, 실험적 표면 장력과 비교 가능한 값을 얻는다. 둘째, 확장된 첫 통과 시간(MFPT) 방법을 사용해 핵생성 속도와 임계 클러스터 크기를 추정하였다. 전통적인 MFPT는 단일 임계 크기 가정을 전제로 하지만, 저자들은 클러스터 성장 과정에서 시간에 따라 변하는 임계 크기를 고려해 확장하였다. 이를 통해 핵생성 장벽 ΔG*와 젤도프 계수 Z를 보다 정확히 도출할 수 있었다.

시뮬레이션 결과는 다음과 같은 주요 특징을 보인다. (1) 표면 장력은 온도가 상승함에 따라 감소하는 경향을 보였으며, 실험값과 5 % 이내의 차이로 일치하였다. (2) 핵생성 속도 J는 온도 상승에 따라 급격히 증가했으며, 273 K에서 10⁻⁹ m⁻³ s⁻¹ 수준, 363 K에서는 10⁻⁴ m⁻³ s⁻¹ 수준에 도달했다. 이는 고전 핵생성 이론(CNT)의 예측과 비교했을 때, 특히 저온 영역에서 CNT가 과소평가함을 시사한다. (3) 임계 클러스터 크기 n는 온도에 따라 감소했으며, 273 K에서는 약 80 개, 363 K에서는 30 개 정도였다. (4) 핵생성 장벽 ΔG는 온도 상승에 따라 감소했으며, 이는 J의 급격한 증가와 일관된다. (5) 성장 단계에서는 클러스터 반경 R(t) 가 R ∝ t^α 형태의 파워‑법칙을 따랐으며, 지수 α는 0.5 ~ 0.6 사이로 온도에 크게 의존하지 않았다. 이는 확산‑제한 성장 메커니즘이 지배적임을 암시한다.

또한, 클러스터의 형태 분석을 통해 핵생성 초기에 비구형(불규칙) 형태가 존재하지만, 성장함에 따라 구형에 가까워지는 경향을 확인했다. 이는 표면 장력 최소화와 내부 압력 균일화가 동시에 진행됨을 의미한다. 저자들은 이러한 형태 변화를 정량화하기 위해 이심률(eccentricity)과 표면 거칠기 파라미터를 도입했으며, 온도 구간 전반에 걸쳐 이값이 감소함을 보고하였다.

마지막으로, 저자들은 시뮬레이션 결과를 기존 실험 데이터와 다른 수치 연구와 비교하였다. 특히, 전통적인 원자 수준 모델(예: TIP4P/2005)과 비교했을 때, mW 모델이 계산 효율성을 크게 향상시키면서도 핵생성 동역학을 충분히 재현한다는 점을 강조한다. 또한, CNT와의 비교에서는 핵생성 장벽과 Zeldovich factor가 CNT보다 실제값에 더 가깝게 예측된다는 점을 부각시켰다. 전반적으로, 본 연구는 물의 균일 핵생성 및 성장 메커니즘을 정량적으로 규명하고, mW 모델과 확장 MFPT 기법이 고전 이론의 한계를 보완하는 데 유용함을 입증한다.