[論文レビュー] Ab initio melting temperatures of bcc and hcp iron under the Earth's inner core condition

Yang Sun, Mikhail I. Mendelev+6 more|arXiv (Cornell University)|May 4, 2022

High-pressure geophysics and materials被引用数 4

ひとこと要約

本研究では、古典的分子動力学（MD）とab initio計算を組み合わせたハイブリッドな熱力学統合（TI）手法を用いて、地球の内核条件下における体心立方（bcc）および立方最密充填（hcp）鉄のab initio融解温度を計算した。有限サイズ効果と不確実性を滑らかな結合経路によって最小化することで、この手法はhcp鉄が熱力学的に安定な相であることを確認した。融解付近では、hcpとbccの自由エネルギー差は約10s meV/原子にとどまり、高圧・高温下での二相のほぼ degeneracy（簡約状態）を示している。

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 の範囲で測定されている実験的およびシミュレーションベースの融解温度測定値に大きな不確実性が存在するという問題に対処すること。
• 有限サイズ効果を低減し、ab initio融点計算の精度を向上させる熱力学的に整合性のある手法を開発・適用すること。
• 融解温度付近での自由エネルギー差を計算することで、hcp鉄とbcc鉄の熱力学的安定性を特定すること。

提案手法

• 古典的系（半経験的ポテンシャルを用いる）からab initio系（密度汎関数理論を用いる）への熱力学統合（TI）を用いて、自由エネルギー差を計算する。
• TI経路は、古典的ポテンシャルとab initio力の間を滑らかに補間する結合パラメータλを介して構築され、自由エネルギーの滑らかな評価が可能になる。
• ギブズ＝ヘルムホルツの式を用いて古典的系の融解温度および自由エネルギー差を計算し、その後TI補正項および体積力働（PV仕事）項による補正を加える。
• 高い精度のDFT設定を用いて、高効率なDFT軌道から再評価することで、自由エネルギー摂動（FEP）法を適用し、精度を向上させる。
• 式 (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/原子のオーダーであり、二相のほぼ同等の安定性（near-degeneracy）を示している。
• 従来のab initio研究に比べ、融解温度の不確実性が著しく低減されて計算された。従来の不確実性は約100–500 Kであったが、本研究ではその範囲が大幅に縮小された。
• 古典的系とab initio系の間を滑らかに結合する経路を用いることで、有限サイズ効果が効果的に低減された。
• 自由エネルギー摂動（FEP）法により、ab initioエネルギー差の信頼性が確認された。そのフラクチュエーションは3–8 meV/原子であり、6000 KにおけるkBT（約500 meV/原子）に比べてはるかに小さい。
• 本研究は、熱力学統合を用いて内核条件下におけるbcc鉄のab initio融解温度を初めて計算し、先行研究における重要な空白を解消した。