[論文レビュー] Ab-initio free energies of liquid metal alloys: application to the phase diagrams of Li-Na and K-Na

Yang Huang, Michael Widom+1 more|arXiv (Cornell University)|May 18, 2021

Thermodynamic and Structural Properties of Metals and Alloys被引用数 4

ひとこと要約

本稿では、第一原理的アプローチを用いて液体金属合金の絶対自由エネルギーを計算する手法を提示する。アービトラリー密度汎関数理論（AIMD）を用い、模擬されたエンタルピーと密度および対相関関数から導かれるエントロピーを組み合わせることで実現する。この手法は、Li-Na（臨界点と液体-液体相分離を示す）およびK-Na（33% Naと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.

研究の動機と目的

• 液体金属合金の絶対自由エネルギーを第一原理から計算する手法を開発し、相や組成間の熱力学的安定性を直接比較可能にする。
• ヒストグラムやサンプリング手法に依存しない、密度および対相関関数から導出されるエントロピーを用いて、エントロピーを絶対スケールに位置づける課題に取り組む。
• Li-NaおよびK-Naという、臨界点や共晶を示す異なる相挙動を示すアルカリ金属系にこの手法を適用し、相図を予測する。
• 実験データとの比較を通じて、液体金属のエンタルピー、エントロピー、自由エネルギー予測の正確性を検証する。
• 有限サイズ効果やエントロピーの二体項のみに制限された展開（截断）が、統計的解析を慎重に行うことで、巨視的相挙動の予測精度を損なわないことを示す。

提案手法

• 有限温度におけるアービトラリー密度汎関数理論（AIMD）シミュレーションを用い、DFTで計算されたエネルギーを基準点として、絶対エンタルピーを計算する。
• 対相関関数と密度のシミュレーション結果から、二体相関に基づく近似を用いてエントロピーを導出し、ヒストограмやサンプリング手法に依存しない絶対エントロピーを実現する。
• シミュレーション温度にスケーリングされたTを用いて、全自由エネルギーG = H - TSを、シミュレーションで得られたエンタルピーHと導出されたエントロピーSから計算する。
• 化学ポテンシャルの計算および自由エネルギーの一貫性の検証には、変更を加えたウィドム粒子挿入法が用いられるが、主たるアプローチはアンサンブル平均による直接的熱力学的統合に依存する。
• 特に相分離過程のための平衡化を加速するために、ハイブリッドモンテカルロ/分子動力学（MCMD）が採用され、配置サンプリングの向上を図るための離散的成分交換を許容する。
• 統計誤差は、シミュレーション実行をセグメントに分割し、標準誤差を計算することで推定され、Hと-TSの負の相関がGの不確実性を低減する。

実験結果

リサーチクエスチョン

• RQ1 熱力学的統合や摂動理論に依存せずに、第一原理的AIMDシミュレーションから液体金属合金の絶対自由エネルギーを直接計算できるか？
• RQ2 高次相関を含まない状態で、密度と対相関関数からのみ、液体金属のエントロピーをどれほど正確に近似できるか？
• RQ3 本手法は、実験と一致するように、Li-Naの臨界点および液体-液体相分離を正しく予測できるか？
• RQ4 本手法は、第一原理から出発して、K-Naの深共晶の組成と融点（33% Naと260 K）を正確に再現できるか？
• RQ5 有限サイズ効果およびエントロピー展開の二体項への制限（多体項との比較）が、液体合金の相境界予測にどのように影響するか？

主な発見

• 本手法は、エンタルピー、密度、対相関関数のAIMDシミュレーション結果のみを用いて、熱力学的統合や基準状態に依存せずに、液体金属合金の絶対自由エネルギーを正確に計算できることを示した。
• Li-Naに対しては、予測された相図に約578 Kでの臨界点と液体-液体相分離が確認され、実験観察と整合的であるが、シミュレーションサイズの影響により臨界温度がやや高估されている。
• K-Naに対しては、33% Naと260 Kでの深共晶が正しく予測され、実験データと一致しており、強い混合異常を示す系に対しても本手法の正確性が確認された。
• 純粋な液体金属（Li, Na, K）の計算された絶対エントロピーは、NIST-JANAF表とよく一致しており、エントロピー導出手法の妥当性が裏付けられた。
• K-Na（473 K）における自由エネルギーの統計誤差は約0.4 meV/原子と推定され、有限サンプリングおよび系サイズの影響にもかかわらず、高い精度を示している。
• 有限サイズ効果により、Li-Naでは臨界点を超えると、特に小さい系（N=300）では早期に相分離が発生するが、より大きな系（N=500）を用いることで、全体的な相挙動は実験と定性的かつ定量的に一致する。