[論文レビュー] A first-principles global multiphase equation of state for hydrogen

Alfredo A. Correa, Lorin X. Benedict+4 more|arXiv (Cornell University)|Jun 4, 2018

Spectroscopy and Quantum Chemical Studies被引用数 3

ひとこと要約

本稿では、量子陽子と電子を含むab initioシミュレーションを用いて、1–10⁹ Kおよび10⁻⁹–1 m³/molの範囲をカバーする、最初からの全相態状態方程式（EOS）を提示する。物理的に妥当な自由エネルギー・モデル、多パラメータ・マルチドレキシブ適合法、熱力学的に整合性のある解析的表現を用い、補間に起因する不連続性を回避して、分子固体、流体、原子液体、相転移を正確に記述する。

ABSTRACT

We present and discuss a wide-range hydrogen equation of state model based on a consistent set of ab initio simulations including quantum protons and electrons. Both the process of constructing this model and its predictions are discussed in detail. The cornerstones of this work are the specification of simple physically motivated free energy models, a general multiparameter/multiderivative fitting method, and the use of the most accurate simulation methods to date. The resulting equation of state aims for a global range of validity ($T = 1-10^9 K$ and $V_m = 10^{-9}-1 m^3/mol$), as the models are specifically constructed to reproduce exact thermodynamic and mechanical limits. Our model is for the most part analytic or semianalytic and is thermodynamically consistent by construction; the problem of interpolating between distinctly different models -often a cause for thermodynamic inconsistencies and spurious discontinuities- is avoided entirely.

研究の動機と目的

• 極端な熱力学的条件下、つまり超冷却分子固体から高温高密度プラズマに至るまで、水素のグローバルに有効な状態方程式を開発すること。
• 別々の相モデル間の補間を避けることで、既存モデルの熱力学的不整合を解消すること。
• 特に低温・低密度領域において顕著な量子核効果（陽子の零点運動）を、すべての相に明示的に組み込むこと。
• ab initioデータから導出した解析的自由エネルギー式を用いて、すべての相で熱力学的整合性を確保すること。
• 慣性核融合および惑星科学の応用分野における水素および重水素・トリチウ化物混合系の統合的・予測可能なモデルを提供すること。

提案手法

• 電子および量子陽子を高精度で取り扱う第一原理的量子分子動力学および密度汎関数理論（DFT）シミュレーションを用いて、グローバルな状態方程式を構築する。
• 分子固体（Mie-Grüneisenおよびコールドカーブ近似を介して）、原子流体（理想気体極限）、混合相（分子-原子流体混合物）に対して、物理的に妥当な自由エネルギー・モデルを開発する。
• 全相にわたるシミュレーションデータに多パラメータ・マルチドレキシブ適合法を適用し、滑らかな遷移と熱力学的整合性を確保する。
• 臨界点付近の分子-原子遷移を記述するため、自己無撞撃的平均場相互作用（結合パラメータJを介して）を用いた非理想混合モデルを適用する。
• 統計力学を用いて、配置エントロピーと非理想相互作用を組み込んだ解析的自由エネルギー式を混合相について導出する。
• 経路積分モンテカルロ（PIMC）データ、衝撃Hugoniot実験、および既知の熱力学的極限（例：理想気体、融解線）との比較を通じて、モデルの妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 第一原理的シミュレーションを用いて、分子固体から原子プラズマに至る全相にわたる、熱力学的に整合性のあるグローバルな水素状態方程式をどのように構築できるか？
• RQ2 量子核効果（零点運動）は、低密度および低温水素相の安定性および熱力学的性質を決定づける役割を果たすか？
• RQ3 別々の相モデル間の補間に起因するアーティファクトを回避して、相転移（例：融解、分子-原子解離）をどのように正確にモデル化できるか？
• RQ4 原子および分子間の非理想相互作用は、高密度水素における臨界定常および相共存にどの程度影響を与えるか？
• RQ5 本モデルは、高圧・高温領域における実験的衝撃HugoniotデータおよびPIMCシミュレーションと、定量的にどの程度一致するか？

主な発見

• 本モデルは、1–10⁹ Kおよび10⁻⁹–1 m³/molの範囲でグローバルに有効であり、分子固体、液体、プラズマを含む、既知のすべての水素相をカバーする。
• 量子陽子の取り入れにより、特に融解線付近の低温分子相の記述が著しく改善され、古典的モデルでは失敗する領域でも良好な記述が可能になる。
• 主な衝撃Hugoniot曲線は、実験データと良好に一致しており、特に10–100 GPaの範囲で、Purgatorioなどの古くからのイオニック・スフィアモデルを上回る性能を示す。
• 自己無撞撃的平均場相互作用（J）を用いた非理想混合モデルは、液体-蒸気ドームおよび臨界領域をうまく再現し、物理的に不自然な不連続性を回避する。
• モデルは、PIMCシミュレーションと整合性のある分子-原子分率の遷移を予測し、流体相全域にわたって滑らかで連続的な変化を示す。
• 熱力学的整合性は構築段階で確保されている：圧力、エントロピー、比熱といったすべての導出物理量が正確な熱力学的関係を満たし、補間に起因するアーティファクトが排除される。