Van-der-Waals exchange-correlation functionals and their high pressure and warm dense matter applications

本文评估了多种交换相关泛函，确定 r2SCAN 在重现氢分子性质方面最为准确，旨在确定其在模拟高压下热致致密氢的分子到金属转变过程中的适用性。

要点

想象一下，你正在试图预测一群人在一个巨大的加压房间里的行为。有时他们作为个体独自站立，有时他们两两牵手（分子），有时压力变得如此之高，以至于他们松开了手，变成了一团混乱的流体（金属）。

这篇论文是一次科学性的“口味测试”，旨在弄清楚哪一套规则（称为交换相关泛函）能最好地预测这些氢“人”在被挤压在一起形成暖致密物质（Warm Dense Matter）时的行为——这种物质存在于像土星这样的巨行星内部或聚变实验中。

以下是利用简单类比对他们调查过程的拆解：

1. 问题所在：选择正确的规则书

科学家使用计算机模拟来模拟这种物质。为了做到这一点，他们需要一本“规则书”（数学公式）来描述电子如何相互作用。目前有很多可用的规则书：

• PBE： 标准且可靠的规则书，大家都在使用。
• vdW（范德华）泛函： 专门设计的规则书，用于解释中性物体之间非常微弱、长距离的“幽灵般”吸引力（就像摩擦头发后的气球会粘在墙上一样）。
• r2SCAN & HSE06： 更新、更复杂的规则书。

核心问题是：这些专门处理“幽灵吸引力”的规则书（vdW）是否真的能让高压氢气的预测变得更好？

2. 实验室测试：检查基础情况

在模拟整个人群之前，作者先在简单的、孤立的场景下测试了这些规则书，以观察哪一个最准确。

• 测试 A：握手（键长与能量）
  他们观察了两个氢原子如何握手。

  • 结果： 标准的 PBE 规则书说这次握手有点太松了（原子间距太远）。vdW 规则书则说握手太紧了（原子间距太近）。
  • 获胜者： r2SCAN 规则书将握手的距离和强度几乎做得完美无缺。
  • 转折： HSE06 规则书也非常出色，但 r2SCAN 是这项特定测试中的冠军。

• 测试 B：“幽灵式”拥抱（两个分子之间的相互作用）
  他们将两个氢分子靠近放置，观察它们是否会感受到那种微弱的“幽灵般”吸引力。

  • 结果： 在这里，vdW 规则书失败了。它们预测分子会拥抱得太紧且距离太远。
  • 获胜者： 出人意料的是，标准的 PBE 和 HSE06 规则书比专门的 vdW 规则书更好地预测了这种“幽灵拥抱”。
  • 输家： r2SCAN 规则书完全忽略了这个拥抱；它根本没有察觉到这种吸引力。

3. 大挤压：模拟高压房间

现在，他们在高压和高温（暖致密物质）下模拟了整个人群。

• “幽灵效应”微乎其微： 作者发现，分子之间的“幽灵”吸引力极其微弱（像耳语一样），相比于打破分子所需的能量（像呐喊一样）。
• 结论： 由于“幽灵”力量如此微弱，它实际上并不会改变高压房间里人群的行为。无论你使用 vdW 规则书还是标准的 P表 PBE 规则书，最终的“人群行为”（压力、密度、结构）看起来几乎是一模一样的。
• 为什么 vdW 规则书之前改变了结果？ 作者意识到，以往使用 vdW 规则书的研究之所以预测了不同的“相变”（即人群何时松开手），并不是因为它们更好地观察到了“幽灵拥抱”，而是因为它们高估了握手的强度。它们认为分子抓得太紧了，所以需要更高的压力才能把它们拆开。

4. 最终裁定

该论文得出结论，对于暖致密氢气：

1. “幽灵”力量并不重要： 分子之间微弱的吸引力太小，不足以改变材料在高压下的宏观行为。
2. 握手最重要： 最关键的是要准确掌握分子键的强度。
3. 最佳规则书： 由于 r2SCAN 把握手（键长和能量）做得非常完美，尽管它没能察觉到“幽灵拥抱”，作者建议它是模拟这种物质的最佳选择。它抓住了最关键的部分，而专门的 vdW 规则书却把握手搞错了。

简而言之： 作者试图寻找一种工具，去测量嘈杂房间里的一声微弱耳语（vdW 力）。他们发现，这声耳语太轻了，根本无关紧要。相反，他们意识到标准的工具其实更擅长听清那声响亮的呐喊（分子键），而一种新工具（r2SCAN）甚至能完美地听清那声呐喊。因此，他们建议使用那个能听清呐喊的工具，而不是那个专门用来听耳语的工具。

技术摘要

技术摘要：范德华交换相关泛函及其在高压与热致致密物质中的应用

问题陈述
近期关于高密度氢的模拟高度依赖于密度泛函理论（DFT），以模拟分子向金属的转变以及液-液相变（LLPT）。尽管采用了包含非定域范德华（vdW）修正的高级交换相关（xc）泛函来改善与实验 LLPT 曲线的一致性，但这些泛函在热致致密物质（WDM）机制中的优越性在物理学上仍存在争议。范德华相互作用的能量尺度（数 meV）比分子解离能（~4.5 eV）或电离能（13.6 eV）低了几个数量级。因此，目前尚不清楚引入范德华修正是否显著改变了氢的物态方程（EOS），或者观察到的 LLPT 线偏移是否仅仅是不同泛函在描述基本分子属性（如键长和解离能）方式上的人工产物。

方法论
作者通过使用维也纳从头算模拟软件包（VASP）进行 DFT 和 DFT-分子动力学（DFT-MD）模拟，开展了一项系统性的基准测试研究，对比了标准与高级交换相关泛函。研究采用了用于热力学 DFT 的 Mermin 公式化形式和裸库仑伪势。

测试的泛函包括：

• 标准型：LDA, PBE。
• 非定域 vdW 型：VDW-DF1, VDW-DF2, VDW-DF3 (optPBE-vdW)。
• Meta-GGA 型：r2SCAN。
• 杂化型：HSE06。

为了建立高精度基准，作者采用了以下手段：

1. 量子蒙特卡洛（QMC）： 特别是使用 HANDE-QMC 软件进行的起始全构型相互作用量子蒙特卡洛（i-FCIQMC），用于计算孤立氢分子二聚体的近乎精确的相互作用能。
2. 耦合簇理论（CCSD(T)）： 用于为 QMC 结果生成完全基组（CBS）修正。

研究涵盖了三个截然不同的机制：

1. 理想真空： 孤立的氢分子（键长、解离能）以及孤立的分子-分子相互作用（vdW 能量极小值）在各种几何构型下的表现。
2. 热致致密/高压： 在不同密度（$r_s$）和温度（100 K 至 1500 K）范围内对氢流体进行 DFT-MD 模拟，跨越了分子、相互作用分子及压力解离机制。
3. 结构与电子分析： 计算对相关函数、静态与动态离子结构因子、有效键长以及电子态密度（DOS）。

主要贡献与结果

• 分子属性：

  • 键长： r2SCAN 泛函提供了对实验键长最好的重现（偏差为 0.001 Å），紧随其后的是 HSE06。标准 PBE 和 vdW 泛函表现出较大的偏差（PBE 高估了 0.1 Å；vdW 泛函则低估或轻微高估）。
  • 解离能： 与实验相比，PBE 低估了约 0.21 eV 的解离能。vd vdW 泛函则高估了约 0.25 eV。r2SCAN 表现出最佳的一致性，偏差仅为 ~0.07 eV。
  • 泛函结论： 在基本分子属性方面，r2SCAN 和 HSE06 的表现优于 PBE 和 vdW 泛函。

• 范德华相互作用：

  • 在理想的两分子相互作用场景中，vdW 能量极小值（深度 ~5 meV）被 HSE06 得到最好的重现，且与 QMC 结果几乎无法区分。PBE 位居第二。
  • 至关重要的是，vdW 泛函（VDW-DF1, VDW-DF3）预测的极小值比 QMC 预测的更深且距离更远。
  • r2SCAN 完全未能重现 vdW 能量极小值，缺乏对这些长程效应的适当描述。因此，r2SCAN 在本研究后续的 WDM 研究中被剔除。

• 物态方程（EOS）与 LLPT：

  • 研究发现，在 LLPT 区域，vdW 泛函预测的压力比 PBE 更高。作者将这种偏移归因于 vdW 泛函预测了更高的解离能，而非直接由于引入了 vdW 力。
  • 较高的解离能意味着需要更多的压力功来破坏分子，从而将分子向金属的转变推向更高压力。
  • 在 LLPT 相关温度下，vdW 相互作用的能量尺度相对于热能是微不足道的（例如，1500 K $\approx$ 130 meV，而 vdW 能量 $\sim$ 10 meV）。

• 结构与动力学属性：

  • 对相关与结构因子： 在分子机制中，观察到分子峰的精细结构存在差异，这与不同泛函预测的键长相关。然而，在中间到长程范围（1–10 Å），各泛函预测的对相关函数和静态结构因子几乎没有差异。
  • 动态结构： 探测集体模式和声速的动态离子结构因子显示，PBE 与 vdW 泛函之间不存在显著差异。集体行为受中程和长程关联支配，而所有测试的泛函对这些关联的描述是相似的。
  • 电子结构： 电子态密度（DOS）被发现几乎完全取决于离子结构（快照几何构型），而非生成该 DOS 的特定交换相关泛函。对于给定的离子构型，无论使用 PBE、vdW 还是其他泛型来计算电子性质，带隙都保持稳定。

意义与主张
论文得出结论：对于热致致密氢，引入非定域 vdW 修正对于描述该系统的物理特性并不比标准泛函（如 PBE）具有根本性的优势。使用 vdW 泛函时观察到的 LLPT 线偏移，主要是因为它们倾向于高估分子解离能，而非直接源于对改进的 vdW 力描述。

作者断言：

1. PBE 仍然是描述热致致密氢的稳健选择，特别是考虑到其相对于杂化泛函（如 HSE06）的计算效率，以及它在描述 vdW 相互作用（在分子间吸引方面）方面几乎与 HSE06 一样出色。
2. 如果决策依据仅为分子键长和解离能的准确性，那么 r2SCAN 是计算 EOS 和 LLPT 的更优泛函，尽管它无法描述 vdW 极小值。
3. 在相关的 WDM 温度和压力环境下，vdW 效应对于氢的 EOS 是可以忽略不计的，因为极化和中性-中性相互作用的能量尺度远小于解离能和电离能。

该研究表明，未来理解压力解离的过程应侧重于电子-电子相关函数和新的量子蒙特卡洛方法，而非选择特定的 vdW 修正交换相关泛函。