[論文レビュー] Theory of complex fluids in the warm-dense-matter regime, and application to an unusual phase-transitions in liquid carbon

M. W. C. Dharma‐wardana|arXiv (Cornell University)|Jul 25, 2016

High-pressure geophysics and materials被引用数 4

ひとこと要約

本稿では、温暖密剤質（WDM）領域における複雑な流体を第一原理的にモデル化するためのニュートラル・擬似原子（NPA）＋ハイパーニットチェーン（HNC）法を導入し、液体炭素における一時的な共有結合を正確に捉える。7 eVおよび1.0 g/cm³付近で、同時に液体-蒸気および金属-半金属の相転移が発生し、平均イオン化度 $ Z $ に急激な不連続性が生じるが、これはDFT+MDおよびPIMCシミュレーションと一致する。

ABSTRACT

Data from recent laser-shock experiments, density-functional theory (DFT) with molecular-dynamics (MD), and path-integral Monte Carlo (PIMC) simulations on carbon are compared with predictions from the neutral-pseudo-atom (NPA)+ hyper-netted-chain (HNC) approach for carbon, a complex liquid in the warm-dense matter regime. The NPA results are in good agreement, not only with high-density regimes that have been studies via PIMC, but even at low densities and low temperatures where transient covalent bonding dominates ionic correlations. Thus the `pre-peak' due to the C-C bond at $\\sim$1.4-1.6 \\AA$\\,$ and other features found in the pair-distribution function from DFT+MD simulations at 0.86 eV and 3.7 g/cm$^3$ etc., are recovered accurately in the NPA+HNC calculations. Such C-C bonding peaks have not been captured via average-atom ion-sphere (IS) models. Evidence for an unusual liquid $\ o$ vapor and metal$\ o$ semi-metal transition occurring simultaneously is presented. Here a strongly correlated metallic-liquid with transient C-C bonds, i.e., carbon at density $\\sim$ 1.0 g/cm$^3$ and mean ionization $Z=4$ transits abruptly to a disordered mono-atomic vapour at 7 eV, with $Z\\simeq$ 3. Other cases where $Z$ drops abruptly are also noted. The nature of $Z$, its discontinuities, and the role of exchange correlation, are reviewed. The limitations of IS models in capturing the physics of transient covalent bonding in warm dense matter are discussed.

研究の動機と目的

• WDM領域における炭素のような複雑な液体を正確にモデル化できる第一原理的手法の開発。従来のモデルでは失敗する領域である。
• 平均原子イオン球（IS）モデルがWDMにおける一時的な共有結合および電子相関を捉えられないという限界を克服すること。
• 強相関金属液体における $ Z $、構造因子、対ポテンシャル、電子スペクトルを計算する統一的理論枠組みの提供。
• 炭素のような低原子番号元素における $ Z $ の不連続転移の性質と、相挙動との関連を解明すること。
• 広範な密度および温度範囲において、高精度なDFT+MDおよびPIMCシミュレーションと比較して、NPA+HNC手法の妥当性を検証すること。

提案手法

• NPAモデルは、電子を1体密度 $ n(r) $ で、イオンを $ \rho(r) $ で扱い、有効な電子-イオンポテンシャルを介してKohn-Sham DFT計算を可能にする。
• ハイパーニットチェーン（HNC）積分方程式を用いて、NPAから得られるポテンシャルから対分布関数（PDF）を計算し、イオン相関を捉える。
• Friedelの和と電荷中性条件を用いて、熱力学的制約と整合性を持つように、自己無撞着に平均イオン化度 $ Z $ を決定する。
• 本手法は、共有結合が支配する低温領域から、DFT+MDでは到達できない高温WDM条件まで滑らかに拡張可能である。
• 有限温度におけるKohn-Sham-Mermin形式を用いて電子構造を計算し、占有数はフェルミ・ディラック統計から導出する。
• 本モデルは、束縛電子と自由電子の寄与を区別し、$ Z = Z_n - n_b $ で表される。ここで $ Z_n = 6 $、$ n_b $ は束縛電子数である。

実験結果

リサーチクエスチョン

• RQ1 NPA+HNC手法は、広範な密度および温度範囲において、液体炭素の対分布関数（PDF）を、一時的なC-C共有結合に起因する前ピークを含めて正確に再現できるか？
• RQ2 7 eVおよび1.0 g/cm³付近で観測された液体-蒸気および金属-半金属の同時相転移の性質は何か？
• RQ3 平均原子イオン球（IS）モデルがWDMにおける共有結合特徴を捉えられない理由は何か？NPAモデルはその限界をどのように克服するか？
• RQ4 7 eV付近での平均イオン化度 $ Z $ の不連続性は、電子構造の変化および炭素における相転移とどのように関連しているか？
• RQ5 NPA+HNC手法は、炭素のような軽元素のWDM研究において、計算コストの高いDFT+MDおよびPIMCシミュレーションにどの程度置き換え可能か？

主な発見

• NPA+HNC手法は、イオン-イオンPDFにおけるC-C結合の前ピーク（~1.4–1.6 Å付近）を正確に再現した。これはISモデルでは存在しないが、DFT+MDおよびPIMCシミュレーションでは確認される特徴である。
• 0.86 eVおよび3.7 g/cm³の条件下で、NPA+HNCのPDFはDFT+MD結果と高い忠実度で一致し、一時的な共有結合の特徴も再現した。
• 7 eVおよび1.0 g/cm³付近で、同時に液体-蒸気および金属-半金属の相転移が発生し、$ Z $ が約4から約3へ急激に低下する。これは電子構造の不連続的変化を示している。
• 平均イオン化度 $ Z $ は7 eV付近で不連続性を示し、T = 7 eVの条件下で2つの解が存在する：$ Z = 4 $（1s軌道が完全に占有、2s軌道が空）と $ Z = 2.78 $（強い相関状態を示す）。
• NPAモデルは、$ Z_{cn} $、$ Z_b $ およびFriedel和の間で、~93%以内の一致を示し、相転移点付近でも一貫性と明確性を持つ $ Z $ 値を提供する。
• 本手法は、一時的なC-C結合を有する強相関金属液体から、不規則な単原子蒸気への遷移を的確に捉えており、圧縮率および導電率の顕著な変化が観察された。