[논문 리뷰] Crystallization of Simple Fluids: Relative Stability of f.c.c. and b.c.c Structures

Swarn Lata Singh, Yashwant Singh|arXiv (Cornell University)|2009. 02. 27.

Material Dynamics and Properties인용 수 11

한 줄 요약

이 논문은 오르누슈-져니케 방정식과 로저-영 클로처를 사용한 대칭 보존 및 대칭 위반 성분을 포함하는 직접 쌍 상관 함수(DPCF)를 갖는 자유 에너지 함수를 개발한다. 이는 대칭 위반 성분에 대한 섭동 전개를 포함하며, n=12인 역제곱력 유체는 f.c.c.로 결정화되며, 더 연성 있는 잠재력(n≤6)은 b.c.c.를 선호함을 정확히 예측한다. 이는 밀도 함수 이론에서 오랫동안 해결되지 않은 과제를 해결한다.

ABSTRACT

A free-energy functional for a crystal that contains both the symmetry conserved and symmetry broken parts of the direct pair correlation function is developed. The free-energy functional is used to investigate the crystallization of fluids interacting via the inverse power potential ; $u(r)=ε{(σ/r)}^n$. In agreement with simulation results we find that for $n=12$ the freezing is into close packed f.c.c structure while for soft repulsions $(n\leq 6)$ b.c.c phase is more stable.

연구 동기 및 목표

• 유체의 결정화 과정에서 f.c.c.와 b.c.c. 상의 상대적 안정성을 정확히 기술하는 제1원리 통계역학 이론을 개발하기 위해.
• 기존의 밀도 함수 이론이 직접 쌍 상관 함수(DPCF)의 대칭 위반 기여를 忽略함으로써 유체-고체 전이를 설명하지 못하는 데서 기인하는 실패를 해결하기 위해.
• 결정상에서 나타나며 융해 시에 사라지는, 양상적으로 새로운 대칭 위반 성분의 DPCF를 포함하기 위해.
• 다양한 연성(n=12, 6, 4)을 가진 역제곱력 잠재력에 대해 몬테카를로 시뮬레이션 데이터와의 예측 능력을 시험하기 위해.

제안 방법

• 자유 에너지 함수 A[ρ] = A_id[ρ] + A_ex[ρ]를 설정하며, 여기서 A_ex는 DPCF의 대칭 보존 c^(0) 및 대칭 위반 c^(b) 성분을 모두 포함한다.
• 대칭 보존 c^(0)(r)을 오르누슈-져니케 방정식과 로저-영 클로처를 사용해 계산하여 밀도 전역에서 열역학적 일관성을 확보한다.
• 대칭 위반 c^(b)(r₁,r₂;[ρ])를 기능적 테일러 전개로 모델링하며, 주요 항은 등방성 유체의 삼체 DPCF와 밀도 변동 ρ_b(r)에 대한 적분으로 표현된다.
• 대칭 위반 성분 c^(b)를 역격자 벡터 G를 사용한 푸리에 공간에서 표현하며, 계수 c^(G)(r)는 차이 벡터 r과 결정 구조에 따라 달라진다.
• 밀도 및 질서 매개변수 장에 대한 기능적 통합을 이중 매개변수화(λ 및 ξ)를 통해 수행하여 결정 상의 초과 자유 에너지 A_ex^(b)를 계산한다.
• 유체-고체 공존점을 찾기 위해 광역 열역학적 잠재력 ΔW = W_l - W를 사용하며, 전이점에서 ΔW/N = 0이 된다.

실험 결과

연구 질문

• RQ1 대칭 보존 및 대칭 위반 DPCF 성분을 모두 포함하는 자유 에너지 함수가 역제곱력 유체에서 f.c.c.와 b.c.c. 상의 상대적 안정성을 정확히 예측하는가?
• RQ2 특히 연성 있는 상호작용에 대해, 대칭 위반 DPCF가 유체와 고체 상 사이의 자유 에너지 차이에 얼마나 기여하는가?
• RQ3 등방성 유체의 삼체 DPCF를 통해 대칭 위반 성분 c^(b)를 섭동적으로 처리하는 방법이 고체 상 안정성에 신뢰할 수 있는 추정치를 제공하는가?
• RQ4 왜 기존의 DFT 접근법은 연성 있는 반발력 잠재력에 대해 올바른 결정화 경로를 예측하지 못하는가?
• RQ5 c^(b)를 포함함으로써 몬테카를로 시뮬레이션 결과와의 일치도가 n=12, 6, 4에서 얼마나 향상되는가?

주요 결과

• n=12일 때, 유리한 대칭 보존 및 대칭 위반 기여로 인해 f.c.c. 구조가 더 안정하다고 예측되며, 몬테카를로 시뮬레이션과 일치한다.
• n=6 및 n=4일 때, b.c.c. 구조가 f.c.c.보다 더 안정해지며, 시뮬레이션 결과를 정확히 재현한다. 이는 대칭 위반 DPCF가 결정 안정성에 결정적인 역할을 함을 시사한다.
• 자유 에너지 차이 ΔW/N에서 대칭 위반 기여는 무시할 수 없으며, n=12에서는 대칭 보존 기여의 약 25%를 차지하고, n=4에서는 거의 50%에 이른다.
• 대칭 위반 항이 없을 경우, 특히 연성 있는 잠재력에 대해 유체의 안정성이 과대평가되어 잘못된 예측을 초래한다.
• 삼체 DPCF를 통한 c^(b)의 섭동 표현은 더 연성 있는 잠재력(n=6, 4)에 대해 그 크기를 과대평가하여, 시뮬레이션보다 더 낮은 γ 값에서 전이가 일어나도록 예측한다.
• f.c.c.와 b.c.c. 구조 간에 c^(G)(r)의 기능적 형태는 크기와 r에 대한 의존성 면에서 크게 다름을 보이며, 이는 두 상의 상이한 안정성에 대한 설명을 제공한다.