[论文解读] 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

一句话总结

本文提出了一种基于中性伪原子（NPA）与超网状链（HNC）方法的第一性原理模型，用于模拟温稠密物质（WDM）区域的复杂流体，能够准确捕捉液态碳中的瞬态共价键合行为。该方法揭示了在约 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 模型通过单体电子密度 $ n(r) $ 和离子密度 $ \rho(r) $ 描述电子与离子，支持基于有效电子-离子势的 Kohn-Sham DFT 计算。
• 采用超网状链（HNC）积分方程，从 NPA 推导的势能中计算对分布函数（PDFs），以捕捉离子关联效应。
• 通过 Friedel 总和与电荷中性条件，自洽地确定平均电离度 $ Z $，确保与热力学约束的一致性。
• 该方法可无缝扩展至低温区域（共价键主导）至 DFT+MD 无法触及的高温 WDM 条件。
• 采用 Kohn-Sham-Mermin 形式在有限温度下计算电子结构，占据数基于 Fermi-Dirac 统计分布。
• 模型区分束缚电子与自由电子贡献，其中 $ Z = Z_n - n_b $，$ Z_n = 6 $，$ n_b $ 为束缚电子数。

实验结果

研究问题

• RQ1 NPA+HNC 方法能否在广泛密度与温度范围内，准确再现液态碳的对分布函数（PDFs），包括来自瞬态 C-C 共价键的前峰？
• RQ2 在约 7 eV 和 1.0 g/cm³ 条件下，观察到的液相到气相与金属到半导体的同步相变的本质是什么？
• RQ3 为何平均原子离子球（IS）模型无法捕捉 WDM 中的共价键特征？NPA 模型如何克服这一局限？
• RQ4 平均电离度 $ Z $ 在约 7 eV 附近的不连续性，如何与电子结构变化及碳的相变相关联？
• RQ5 在轻元素如碳的 WDM 研究中，NPA+HNC 方法在多大程度上可替代计算成本高昂的 DFT+MD 和 PIMC 模拟？

主要发现

• NPA+HNC 方法在离子-离子 PDF 中准确恢复了约 1.4–1.6 Å 处的 C-C 键合前峰，该特征在 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 时存在两个解：一个为 $ Z = 4 $（1s 全充满，2s 全空），另一个为 $ Z = 2.78 $，表明存在强关联态。
• NPA 模型在相变点附近仍能提供一致且明确的 $ Z $ 值，$ Z_{cn} $、$ Z_b $ 与 Friedel 总和结果的一致性达到约 93%。
• 该方法成功捕捉了从具有瞬态 C-C 键的强关联金属液体到无序单原子气相的转变，表现为可压缩性与电导率的显著变化。