[논문 리뷰] Ab initio melting temperatures of bcc and hcp iron under the Earth's inner core condition

Yang Sun, Mikhail I. Mendelev외 6명|arXiv (Cornell University)|2022. 05. 04.

High-pressure geophysics and materials인용 수 4

한 줄 요약

이 연구는 고전적 분자 동역학와 밀도함수이론을 조합한 하이브리드 열역역학 통합(TI) 방법을 사용하여 지구 내핵 조건에서 체적 중심 구조(bcc)와 허브리드 구조(hcp) 철의 ab initio 융해 온도를 계산한다. 유한한 크기 효과와 불확실성을 부드러운 결합 경로를 통해 최소화함으로써, 이 방법은 hcp 철이 열역학적으로 안정된 상임을 확인하였으며, 융해 근처에서 hcp와 bcc 간 자유 에너지 차이가 약 ~10 meV/atom 수준임을 보여주어 고압 및 고온 조건에서 두 상이 거의 일치하는 자유 에너지를 가짐을 시사한다.

ABSTRACT

There has been a long debate on the stable phase of iron under the Earth's inner core conditions. Because of the solid-liquid coexistence at the inner core boundary, the thermodynamic stability of solid phases directly relates to their melting temperatures, which remains considerable uncertainty. In the present study, we utilized a semi-empirical potential fitted to high-temperature ab initio data to perform a thermodynamic integration from classical systems described by this potential to ab initio systems. This method provides a smooth path for thermodynamic integration and significantly reduces the uncertainty caused by the finite-size effect. Our results suggest the hcp phase is the stable phase of pure iron under the inner core conditions, while the free energy difference between the hcp and bcc phases is tiny, on the order of 10s meV/atom near the melting temperature.

연구 동기 및 목표

• 극한의 지구 내핵 조건에서 bcc 및 hcp 철의 융해 온도에 대한 오랫동안 지속된 불확실성을 해결하기 위해.
• 현재 4,850 K에서 7,600 K 사이로 다양하게 측정되고 있는 실험 및 시뮬레이션 기반 융해 온도 측정의 상당한 불확실성을 해결하기 위해.
• 유한한 크기 효과를 줄이고 융해 온도 계산의 정확도를 향상시키는 열역학적으로 일관된 방법을 개발하고 적용하기 위해.
• 융해 온도 근처에서 자유 에너지 차이를 계산하여 hcp 철과 bcc 철의 열역학적 안정성을 규명하기 위해.

제안 방법

• 자유 에너지 차이를 계산하기 위해 고전적 시스템(반경량 잠재력 사용)에서 ab initio 시스템(DFT 사용)으로의 열역학적 통합(TI)을 적용한다.
• TI 경로는 고전적 잠재력과 ab initio 힘 사이를 보간하는 결합 매개변수 λ를 통해 구성되며, 부드러운 자유 에너지 평가를 가능하게 한다.
• 고전적 융해 온도와 자유 에너지 차이를 계산하기 위해 기브스-헬름홀츠 방정식을 사용하며, 이는 TI 및 PV 일항을 통해 보정된다.
• 정확도를 향상시키기 위해 자유 에너지 변동(FEP) 방법을 적용하여 고효율 DFT 궤적에서 유도된 고정밀 DFT 설정을 사용해 ab initio 에너지를 재평가한다.
• 식 (S3)–(S11)은 ab initio 자유 에너지 차이를 고전적 융해 기여도, TI 보정항, 압력-부피 일항으로 분해하는 데에 사용된다.
• 고전적 시스템과 ab initio 시스템 간 부피 및 압력의 차이를 고려하여 열역학적 통합의 일관성을 확보한다.

실험 결과

연구 질문

• RQ1 지구 내핵 압력 및 온도 조건에서 순수 철의 bcc 또는 hcp 상 중 어느 것이 더 높은 융해 온도를 가지는가?
• RQ2 유한한 크기 효과와 잠재력의 부정확성이 최소화되었을 때 ab initio 융해 온도는 얼마나 정확하게 계산될 수 있는가?
• RQ3 내핵 조건에서 융해 온도 근처에서 hcp와 bcc 철 간의 자유 에너지 차이는 얼마인가?
• RQ4 두 상 간 융해 온도의 차이가 작을 경우 열역학적 안정성이 신뢰성 있게 평가될 수 있는가?
• RQ5 고전적 및 ab initio 시뮬레이션의 조합은 융해 온도 예측의 불확실성을 어떻게 줄이는가?

주요 결과

• hcp 철은 지구 내핵 조건에서 열역학적으로 안정된 상이다.
• 융해 온도 근처에서 hcp와 bcc 철 간 자유 에너지 차이는 약 10s meV/atom 수준이며, 이는 두 상이 고압 및 고온 조건에서 거의 자유 에너지가 일치함을 시사한다.
• 이전 ab initio 연구에서 약 100–500 K의 불확실성이 있었던 것과는 달리, hcp 철의 융해 온도는 상당히 감소된 불확실성으로 계산되었다.
• 고전적 시스템과 ab initio 시스템 사이의 부드러운 결합 경로를 사용함으로써 이 방법은 유한한 크기 효과를 성공적으로 줄였다.
• 자유 에너지 변동(FEP) 방법은 ab initio 에너지 차이의 신뢰성을 확인하였으며, 6000 K에서 kBT (~500 meV/atom)에 비해 훨씬 작은 변동(3–8 meV/atom)을 보였다.
• 이 연구는 열역학적 통합을 통해 내핵 조건에서 bcc 철의 첫 번째 ab initio 융해 온도를 제공하여 이전 연구에서의 핵심적 공백을 해결하였다.