[논문 리뷰] A general relation between the largest nucleus and all nuclei distributions for free energy calculations

Joël Puibasset|arXiv (Cornell University)|2022. 07. 28.

nanoparticles nucleation surface interactions참고 문헌 24인용 수 3

한 줄 요약

이 논문은 분자 시뮬레이션에서 모든 핵의 확률분포(pa)와 가장 큰 핵의 분포(pl) 사이의 정확한 해석적 관계를 유도하여, 널리 사용되는 s0-방법에서 유래하는 편향 없이 정확한 자유에너지 장벽 계산을 가능하게 한다. 이 방법은 pl과 pa 사이의 엄밀한 관계를 통해 핵성장 자유에너지 프로파일을 편향된 시뮬레이션으로부터 정확히 재구성할 수 있으며, 가장 큰 핵만 모니터링할 수 있는 계산적 제약 조건 하에서도 자유에너지 장벽을 정확히 계산할 수 있다.

ABSTRACT

Prediction of nucleation rates in first order phase transitions requires the knowledge of the barrier associated to the free energy profile $W$. Molecular simulations offer a direct route through $W = -kT \ln p_a$, where $k$ is Boltzmann's constant, $T$ is temperature, and $p_a$ the probability distribution of the size of any nucleus. But in practice, the extremely scarce spontaneous occurrence of large nuclei impedes the full determination of $p_a$, and a numerical bias must be introduced, e.g. on the size of the largest nucleus in the system, leading to the probability size distribution of the largest nucleus $p_l$. Although $p_l$ is known to be system size dependent, unlike $p_a$, it has been extensively used as an approximation for $p_a$. This paper proposes an exact relation between $p_a$ and $p_l$, which cures this approximation and allows an exact calculation of free energy barriers from biased simulations.

연구 동기 및 목표

• 모든 핵 분포(pa)를 가장 큰 핵 분포(pl)로 근사함으로써 발생하는 핵성장 자유에너지 계산의 오랜 편향 문제를 해결하기 위해.
• pa와 pl 사이의 정확한 수학적 관계를 유도하여 임의의 截단값 또는 경험적 조정이 필요 없도록 하기 위해.
• 모든 1차 상전이(기포 형성, 응축, 결정화 포함)에 일반적으로 적용 가능한, 시스템 크기와 무관한 방법을 제공하기 위해.
• 계산 제약 조건으로 인해 가장 큰 핵만 모니터링 가능한 편향된 시뮬레이션에서도 정확한 자유에너지 장벽 계산을 가능하게 하기 위해.

제안 방법

• 순서 통계학과 극값 이론을 사용하여 모든 핵의 확률분포 pa(s)와 가장 큰 핵의 분포 pl(s) 사이의 정확한 해석적 관계를 유도한다.
• s > 0 에서 pl(s) / [1 - F_l(s)] 라는 식을 적용하며, 여기서 F_l(s)는 pl의 누적분포함수로, 이는 기본 통계역학과 일관성을 유지한다.
• 메타안정 상태에서 Lennard-Jones 유체의 기포 형성에 대한 분자 동역학 시뮬레이션을 통해 방법을 검증한다.
• 국소 밀도 기준과 군집 분석을 이용한 박막 기반 탐지 방법을 사용하여 기포(핵)를 식별하고 핵 크기를 정확히 추적한다.
• 이산 박막 부피에 대해 핵 크기 히스토그램(pa 및 pl)을 정규화하여 연속 확률밀도함수와 일관성을 확보한다.
• pl(s)로부터 보정된 pa(s)가 편향 없는 시뮬레이션 결과와 구분 불가능한 자유에너지 프로파일을 제공함을 보여주며, 조정 가능한 매개변수가 없음을 입증한다.

실험 결과

연구 질문

• RQ1 분자 시뮬레이션에서 모든 핵의 분포(pa)와 가장 큰 핵의 분포(pl) 사이에 정확한 해석적 관계를 도출할 수 있는가?
• RQ2 제안된 관계는 자유에너지 계산에서 s0-방법이 유도하는 체계적 오류를 어떻게 제거하는가?
• RQ3 새로운 방법은 임의의 截단값 또는 경험적 조정에 의존하지 않고 진정한 자유에너지 장벽을 어느 정도 정확히 복원하는가?
• RQ4 이 방법은 기포 형성, 응축, 결정화와 같은 다양한 핵성장 현상에 일반적으로 적용 가능한가?
• RQ5 pl(s)로부터 유도된 보정된 pa(s)는 편향 없는 시뮬레이션의 자유에너지 프로파일을 높은 정확도로 재현할 수 있는가?

주요 결과

• 유도된 관계 p_a(s) = p_l(s) / [1 - F_l(s)] 는 가장 큰 핵 분포 pl에서 진정한 모든 핵 분포 pa로의 정확한 매핑을 가능하게 하여 s0-방법에서 발생하는 편향을 완전히 제거한다.
• 보정된 자유에너지 프로파일 −kT ln p_a(s) 는 조정 가능한 매개변수가 없이도 편향 없는 시뮬레이션 결과와 정확히 일치하며, N = 442, 3375, 8000 원자에서 입증되었다.
• 이 방법은 s0-방법에서 발생하는 일정한 이동 오차를 성공적으로 제거하여, s0에 대한 임의의 비물리적 의존성을 제거함으로써 자유에너지 장벽을 왜곡시키는 것을 방지한다.
• 이 접근법은 기포 형성, 응축, 결정화를 포함한 모든 1차 상전이에 일반적으로 적용 가능하며, 시스템 크기나 핵성장 메커니즘과 무관하다.
• 편향된 시뮬레이션에서 가장 큰 핵만 모니터링 가능한 조건 하에서도 정확한 자유에너지 계산을 가능하게 하여 계산 효율성과 신뢰성을 크게 향상시킨다.
• 검증 결과, pl(s)로부터 유도된 보정된 pa(s)는 편향 없는 샘플링으로 확보된 진정한 pa(s)와 구분 불가능한 자유에너지 프로파일을 제공함으로써, 이 방법의 정확성을 확인한다.