[논문 리뷰] Ab-initio free energies of liquid metal alloys: application to the phase diagrams of Li-Na and K-Na

Yang Huang, Michael Widom외 1명|arXiv (Cornell University)|2021. 05. 18.

Thermodynamic and Structural Properties of Metals and Alloys인용 수 4

한 줄 요약

이 논문은 밀도와 쌍상관함수로부터 유도된 엔트로피와 시뮬레이션된 엔탈피를 조합하여 액체 금속 합금의 절대 자유에너지 계산을 위한 최초 원리적 방법을 제시한다. 아보-아이디어 모래(AB-Initio Molecular Dynamics, AIMD)를 사용하여, 리튬-나트륨(Li-Na)의 경우 임계점과 액체-액체 상분리가 나타나고, 칼륨-나트륨(K-Na)의 경우 33% 나트륨 및 260 K에서 깊은 공융점을 보이며, 실험 데이터와 강한 일치를 보인다. 이는 복잡한 액체 시스템에서 절대 자유에너지 계산을 위한 강력하고 계산 효율적인 길을 제시한다.

ABSTRACT

Comparison of free energies between different phases and different compositions underlies the prediction of alloy phase diagrams. To allow direct comparison, consistent reference points for the energies or enthalpies are required, and the entropy must be placed on an absolute scale, yielding absolute free energies. Here we derive absolute free energies of liquids from ab-initio molecular dynamics (AIMD) by combining the directly simulated enthalpies with an entropy derived from simulated densities and pair correlation functions. As an example of the power of this method we calculate the phase diagrams of two binary alkali metal alloys, Li-Na and K-Na, revealing a critical point and liquid-liquid phase separation in the former case, and a deep eutectic in the latter. Good agreement with experimental data demonstrates the power of this simple method.

연구 동기 및 목표

• 액체 금속 합금의 절대 자유에너지 계산 방법을 최초 원리로부터 개발하여, 다양한 상과 조성 간의 열역학적 안정성 비교를 직접 가능하게 한다.
• 기본 상태나 수정 이론에 의존하지 않고, 시뮬레이션된 밀도와 쌍상관함수로부터 엔트로피를 유도하여 엔트로피를 절대 척도에 놓는 문제를 해결한다.
• 리튬-나트륨과 칼륨-나트륨과 같은 이종 상행동을 보이는 알칼리 금속 시스템에 대해 이 방법을 적용하여 상도표를 예측한다. 이들에는 임계점과 공융점이 포함된다.
• 실험 데이터와의 비교를 통해 자유에너지, 엔탈피, 엔트로피 예측의 정확성을 확보한다.
• 유한한 크기 효과와 엔트로피 전개에서의 절단(쌍체 수준에 국한)이 통계적 분석을 철저히 하면 거대계의 상행동 예측 능력을 손상시키지 않음을 보여준다.

제안 방법

• 유한 온도에서 DFT로 계산된 에너지를 기준점으로 삼아, 아보-아이디어 분자 동역학(AIMD) 시뮬레이션을 통해 절대 엔탈피를 계산한다.
• 쌍상관함수와 밀도를 기반으로 한 근사법을 사용하여 시뮬레이션된 궤도분포함수에서 엔트로피를 유도함으로써, 히스토그램이나 샘플링 방법에 의존하지 않는 절대 엔트로피 계산이 가능하다.
• 시뮬레이션 온도에 비례하는 T를 사용하여 시뮬레이션된 엔탈피 H와 유도된 엔트로피 S를 바탕으로 총 자유에너지 G = H - TS를 계산한다.
• 화학적 포텐셜을 계산하고 자유에너지 일관성을 검증하기 위해 수정된 위도움 입자 삽입 방법을 사용하지만, 주요 방법은 앙상블 평균을 통한 직접 열역학적 통합에 기반한다.
• 특히 상분리 과정에서의 평형을 가속화하기 위해 하이브리드 몬테카를로/분자 동역학(MCMD) 방법을 사용하여 이산적 종류 교환을 통해 구성공간 샘플링을 향상시킨다.
• 통계적 오차는 시뮬레이션 실행을 세그먼트로 나누어 표준 오차를 계산함으로써 추정하며, H와 -TS 간의 반상관관계가 G의 불확도를 감소시킨다.

실험 결과

연구 질문

• RQ1 열역학적 통합이나 수정 이론에 의존하지 않고, 아보-아이디어 분자 동역학(AIMD)에서 직접적으로 액체 금속 합금의 절대 자유에너지 계산이 가능한가?
• RQ2 고차 상관관계 없이 쌍상관함수와 밀도만으로 액체 금속의 엔트로피를 얼마나 정확히 근사할 수 있는가?
• RQ3 이 방법은 실험과 일치하는 바, 리튬-나트륨(Li-Na)의 임계점과 액체-액체 상분리를 정확히 예측하는가?
• RQ4 이 방법은 첫 번째 원리로부터 칼륨-나트륨(K-Na)의 깊은 공융점의 조성과 융해 온도를 정확히 재현할 수 있는가?
• RQ5 유한한 크기 효과와 엔트로피 전개의 절단(쌍체 대 다체)이 액체 합금의 상경계 예측에 어떤 영향을 미치는가?

주요 결과

• 이 방법은 엔탈피, 밀도, 쌍상관함수의 AIMD 시뮬레이션 자료만으로도 액체 금속 합금의 절대 자유에너지 계산이 가능하며, 열역학적 통합이나 기준 상태가 필요로 하지 않는다.
• 리튬-나트륨(Li-Na)의 경우 예측된 상도표는 약 578 K에서 임계점과 액체-액체 상분리를 보이며, 실험 관측과 일치하지만, 시뮬레이션 크기 효과로 인해 임계 온도가 과대평가된다.
• 칼륨-나트륨(K-Na)의 경우 33% 나트륨 및 260 K에서 깊은 공융점이 정확히 예측되며, 실험 데이터와 일치하여 강한 혼합 이방성 현상을 보이는 시스템에 대한 방법의 정확성을 확인한다.
• 순수한 액체 금속(Li, Na, K)의 계산된 절대 엔트로피는 NIST-JANAF 표와 잘 일치하여 엔트로피 유도 방법의 타당성을 검증한다.
• K-Na의 473 K에서 자유에너지의 통계적 오차는 약 0.4 meV/atom로 추정되며, 이는 유한한 샘플링과 시스템 크기에도 불구하고 높은 정밀도를 나타낸다.
• 유한한 크기 효과로 인해 리튬-나트륨(Li-Na)은 임계점 이상에서 조기 상분리를 보이지만, 특히 N=300인 작은 세포에서는 그렇지만, 더 큰 시스템(N=500)을 사용할 경우 전체 상행동은 실험과 정량적·정성적으로 일치한다.