[論文レビュー] Evidence for supercritical behavior of high-pressure liquid hydrogen

Bingqing Cheng, Guglielmo Mazzola+1 more|arXiv (Cornell University)|Jun 7, 2019

High-pressure geophysics and materials被引用数 3

ひとこと要約

本研究では、第一原理的量子計算に基づく機械学習ポテンシャルを用いて、極限条件下における高密度水素の相挙動をシミュレートする。液体-液体転移は、急激な一次転移ではなく、臨界点が約350 GPaおよび380 Kに位置する超臨界現象であることが判明し、実験的不一致を解消する。これは、急激な変化ではなく、ウィドム線に類似した漸近的変化を示すことを示している。

ABSTRACT

Hydrogen exhibits unusual behaviors at megabar pressures, with consequences for planetary science, condensed matter physics and materials science. Experiments at such extreme conditions are challenging, often resulting in hard-to-interpret and controversial observations. We present a theoretical study of the phase diagram of dense hydrogen, using machine learning to overcome time and length scale limitations while describing accurately interatomic forces. We reproduce the re-entrant melting behavior and the polymorphism of the solid phase. In simulations based on the machine learning potential we find evidence for continuous metallization in the liquid, as a first-order liquid-liquid transition is pre-empted by freezing. This suggests a smooth transition between insulating and metallic layers in giant gas planets, and reconciles existing discrepancies between experiments as a manifestation of supercritical behavior.

研究の動機と目的

• 高圧水素における液体-液体転移（LLT）の性質に関して長年の実験的・理論的不一致を解消すること。
• 大規模系および長時間スケールのシミュレーションにおいて、従来のDFTおよびQMCシミュレーションの限界を克服すること。
• 高密度水素におけるLLTが一次転移であるか、連続的で超臨界現象であるかを特定すること。
• 実験データおよび惑星モデルの解釈に、熱力学的に整合性のあるフレームワークを提供すること。
• 臨界点行動および多形性が高密度水素の相図に与える役割を評価すること。

提案手法

• ビーラー＝パリネロフレームワークに基づく機械学習ポテンシャル（MLP）を用い、高精度のab initioデータでトレーニングする。
• MLPを用いて、広範な圧力範囲（25–400 GPa）および温度範囲（100–4000 K）において大規模かつ長時間スケールの分子動力学的シミュレーションを実施する。
• 等温等圧系を用いて相挙生をシミュレートし、熱力学的安定性を評価する。
• 小スケール系における直接ab initioシミュレーションとの比較を通じて、MLPの正確性を検証する。
• 液体相における分子的・原子的水素の不混和性を記述するために、多形的溶液モデルを適用する。
• 密度、モル比熱、拡散係数などの熱力学的観測量を分析し、臨界およびウィドム線行動を同定する。

実験結果

リサーチクエスチョン

• RQ1 高圧水素における液体-液体転移は一次転移であるか、連続的で超臨界現象であるか？
• RQ2 メガバー圧力下における水素相図における臨界点の位置と性質は何か？
• RQ3 LLTの実験的観測における不一致はどのように生じるのか。超臨界行動によってそれらを統合できるか？
• RQ4 過去のDFTシミュレーションにおける有限サイズ効果が、誤った一次転移を引き起こす程度はどの程度か？
• RQ5 多形的溶液モデルは、高密度水素における分子的から原子的への転移の熱力学を正確に記述できるか？

主な発見

• 臨界点が約350 GPaおよび380 Kに位置することが同定され、液体-液体転移が連続的かつ超臨界的になる。
• この臨界温度を超えると、系は相分離を示さず、超臨界行動を示し、熱力学的異常は明確なウィドム線に従う。
• 過去のDFT研究で観測された密度の不連続性は、有限サイズ効果および欠陥を有する固体相の形成に起因する。
• 臨界点は融解線の下にあるため、凍結が一次転移を妨げ、滑らかな金属化過程が生じる。
• 比熱および密度は、異なるウィドム線に従う明確な異常を示し、これが異なる実験が異なる転移位置を報告する理由を説明する。
• MLPに基づくシミュレーションは再入る融解および固体多形性を再現し、ab initioベンチマークとの整合性により、手法の正確性を検証した。