[논문 리뷰] Surface Phase Diagrams from Nested Sampling

Mingrui Yang, Lívia B. Pártay외 1명|arXiv (Cornell University)|2023. 08. 16.

nanoparticles nucleation surface interactionsEarth and Planetary Sciences인용 수 3

한 줄 요약

이 논문은 유한 온도에서 흡착된 레너드존스 입자에 대한 표면 상도를 계산하기 위해 내재 샘플링(NS) 프레임워크를 도입한다. 이는 분할 함수를 직접 샘플링함으로써 정확한 자유 에너지 계산을 가능하게 한다. 이 방법은 계단진 표면인 LJ(110)와 LJ(311)에서 더 깊은 결합 위치로 인해 엔탈피 구동의 응축과 엔트로피 억제의 배열이 발생하며, 상전이가 재배열이 아닌 응축에 의해 지배됨을 밝혀내며, T ≈ 0.8 kBT/ε에서 표면 밀도에 영향을 받지 않는 응축 온도를 규명한다.

ABSTRACT

Studies in atomic-scale modeling of surface phase equilibria often focus on temperatures near zero Kelvin due to the challenges in calculating the free energy of surfaces at finite temperatures. The Bayesian-inference-based nested sampling (NS) algorithm allows for modeling phase equilibria at arbitrary temperatures by directly and efficiently calculating the partition function, whose relationship with free energy is well known. This work extends NS to calculate adsorbate phase diagrams, incorporating all relevant configurational contributions to the free energy. We apply NS to the adsorption of Lennard-Jones (LJ) gas particles on low-index and vicinal LJ solid surfaces and construct the canonical partition function from these recorded energies to calculate ensemble averages of thermodynamic properties, such as the constant-volume heat capacity and order parameters that characterize the structure of adsorbate phases. Key results include determining the nature of phase transitions of adsorbed LJ particles on flat and stepped LJ surfaces, which typically feature an enthalpy-driven condensation at higher temperatures and an entropy-driven reordering process at lower temperatures, and the effect of surface geometry on the presence of triple points in the phase diagrams. Overall, we demonstrate the ability and potential of NS for surface modeling.

연구 동기 및 목표

• 복잡한 표면 시스템에서 유한 온도에서 자유 에너지를 계산하는 데 있어 전통적인 몬테카를로 및 아비니오 메서드의 한계를 극복하기 위해.
• 직접 분할 함수를 계산함으로써 영도 온도 이상의 표면 상 평형을 모델링하기 위해.
• 특히 계단 가장자리가 상전이의 성격과 순서에 어떤 영향을 미치는지, 삼중점의 존재 및 성격을 포함하여 조사하기 위해.
• 저지수 및 경사면 표면에서 엔탈피 구동의 응축과 엔트로피 구동의 배열 간 상호작용을 정량화하기 위해.
• 사전에 안정 상에 대한 지식이 없이도 열역학적으로 일관된 상도를 구축할 수 있는 강력하고 효율적인 방법을 개발하기 위해.

제안 방법

• 흡착된 레너드존스 입자에 대한 캐논리컬 분할 함수를 계산하기 위해 베이지안 추론 기반 내재 샘플링(NS) 알고리즘을 적용한다.
• 표준 MCMC 방법이 에르고디시티 붕괴 문제를 겪는 깊은 잠재 에너지 우물을 가진 시스템에서도 효율적으로 구성 공간을 샘플링한다.
• 기록된 에너지의 집합 평균에서 등적열용량(CV) 및 순서 매개변수와 같은 열역학적 성질을 계산한다.
• 순서 매개변수를 사용하여 구조 전이를 특성화하고, 무질서한 상과 질서 있는 상을 구분한다.
• 온도와 표면 밀도에 따른 표면 협곡의 점유 수를 분석하여 상전이의 미세한 기구를 탐구한다.
• 다중 독립적인 NS 실행을 통해 결과를 검증하고, 집합 평균의 표준편차를 통해 통계적 불확실성을 계산한다.

실험 결과

연구 질문

• RQ1 특히 계단 가장자리가 흡착된 레너드존스 시스템에서 상전이의 성격과 순서에 어떤 영향을 미치는가?
• RQ2 엔탈피 구동의 응축과 엔트로피 구동의 배열이 서로 공존하거나 경쟁하는 정도는 어느 정도이며, 다양한 표면 면적에서 어떻게 나타나는가?
• RQ3 내재 샘플링는 계단진 표면과 같은 과도한 전이 피크가 겹치는 시스템의 상 경계를 정확히 해석할 수 있는가?
• RQ4 표면 대칭성과 결합 위치의 이방성은 엔트로피 구동의 배열 전이를 억제하는 데 어떤 역할을 하는가?
• RQ5 관측된 상전이 온도에 따르면 계단진 표면에서 응축 전이가 표면 밀도에 영향을 받지 않는가?

주요 결과

• 평탄한 LJ(111) 및 LJ(100) 표면에서는 흡착체 상이 명확한 두 단계 전이를 보이며, 고온에서 엔탈피 구동의 응축이 일어난 후 저온에서 엔트로피 구동의 재배열이 발생한다.
• 계단진 표면인 LJ(110) 및 LJ(311)에서는 깊은 협곡에서의 강한 결합으로 인해 무질서 상이 사라지며, 기체에서 질서 있는 상으로 직접 전이된다.
• LJ(110)에서의 응축 전이는 표면 밀도에 영향을 받지 않으며, 모든 표면 밀도에서 일관된 온도 T ≈ 0.8 kBT/ε에서 발생한다. 이는 상도에서 점선으로 표시된 수직선으로 나타난다.
• LJ(110)에서의 열용량 곡선은 응축에 기인한 넓고 고온의 피크를 보이며, 저온의 배열 전이 피크를 가리켜 상 경계 식별을 어렵게 한다.
• 순서 매개변수 분석 및 협곡 점유 데이터는 T ≈ 0.8 kBT/ε 이하에서 입자가 협곡 내에서 빠르게 점유하고 최대 배치를 이루며 협동적인 배열 과정을 이끈다고 확인한다.
• 점유된 협곡 수는 온도가 감소함에 따라 단조적으로 증가하며, 고밀도에서는 단위 세포당 4개에 도달하고, 저밀도에서는 ⌈N/4⌉에 도달함으로써 강한 국소적 질서를 나타낸다.