🔬 凝聚态物理

Enhanced Coulomb pseudopotentials in high-pressure hydrides from first principles

本研究通过第一性原理计算表明，由于费米能级处存在高态密度，高压氢化物中的库仑伪势显著大于传统的经验范围，这导致以往的预测大幅高估了超导临界温度。

原作者： Kieran Bozier, Chris J. Pickard

发布于 2026-10-12

📖 1 分钟阅读☕ 轻松阅读

原作者： Kieran Bozier, Chris J. Pickard

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

几十年来，科学家们一直在寻找一种特定的材料：一种在远高于深空极寒温度下仍能实现零电阻导电的材料。这种被称为超导现象的现象，有望彻底改变从电网到医学扫描仪的一切领域，但通常需要将材料冷却至接近绝对零度。近年来，一类被称为氢化物（富含氢的化合物）的新型材料脱颖而出，成为打破这一温度壁垒的有力竞争者。在巨大的压力下，这些物质可以在足以融化冰块的温度下实现超导，这一发现使该领域从纯粹的实验转向了计算机引导的预测。然而，在计算机预测与实验测量之间，一直存在着一个持久的差距。缺失的部分是材料内部一种微妙且无形的力：电子之间的排斥力。虽然允许电子配对并无电阻流动的吸引力已被深入理解，但将它们推开的力量在计算中通常被视为一个粗略的猜测，导致预测结果有时比现实高出两倍。

剑桥大学和东北大学的研究团队现在通过从第一性原理出发计算这种排斥力，而非依赖标准的猜测，填补了这一空白。他们专注于一组高压氢化物，其中包括一些最有希望实现高温超导的候选材料。利用他们开发的一种专门计算机程序，研究人员绘制出了这些材料在极端压力下电子如何相互作用的图谱。他们的工作表明，排斥力并不是此前假设的那种恒定的、普适的数值，而是一个根据材料特定电子结构而剧烈变化的变量。在某些情况下，这种力量非常强大，以至于显著降低了材料实现超导的温度，从而修正了之前的过度估计，使理论预测与实验数据保持一致。

这些材料的故事始于超导体内部微妙的平衡。为了让电流无电阻流动，电子必须克服它们相互排斥的自然倾向，转而形成配对。这种配对通常由晶格中的振动介导，这一过程就像是一种“胶水”。然而，电子也会互相推挤，如果这种推力过强，就会破坏配对并摧毁超导状态。过去，试图预测新材料性能的科学家必须估算这种排斥推力的强度。由于直接计算极其困难，他们通常使用源自铝等简单金属的固定平均值作为估算。这对于旧材料来说足够有效，但随着研究人员开始探索高压下的复杂富氢化合物，这种捷径开始失效。标准的猜测导致预测的超导温度往往过于乐观，在理论与实验室之间造成了令人困惑的脱节。

为了解决这个问题，研究人员采用了一种将排斥力视为可计算的物理现实，而非猜测的方法。他们研究了包括硫化氢、十氢化镧以及一种新预测的化合物——六氢化镁铱在内的几种氢化物系统。他们使用一种新的软件工具来模拟这些材料中电子的行为，精确计算了在超导发生的能量水平处排斥力的行为。结果令人震惊。他们发现，对于许多此类氢化物，其排斥力比传统猜测的要强得多。在六氢化镁铱的案例中，计算出的排斥力几乎是大多数研究中使用的常规估值的四倍。这种差异并非细枝末节；它使预测的超导温度改变了40开尔文，这是一个巨大的偏移，解释了为什么之前的模型高估了该材料的潜力。

研究人员发现，这种排斥力的强度与费米能级处的可用电子态密度紧密相关，费米能级是一个作为电子活动阈值的特定能量水平。在这些电子密度较高的材料中，电子排列得如此紧密，以至于它们的相互排斥力成为了一个主导因素。这看似违反直觉，因为高态密度通常被视为超导的利好因素，因为它有助于电子配对的初步形成。然而，研究表明，在这些高压氢化物中，正是这种有助于形成配对的高密度，同时也放大了试图破坏配对的排斥力。研究人员发现，在像六氢化镁铱这样的系统中，这种排斥力如此剧烈，以至于显著抑制了超导温度，从而解决了两个此前存在两倍差异的理论方法之间的长期分歧。

研究还观察了这些力量在材料受到不同压力挤压时是如何变化的。他们发现，虽然排斥力会随着压力的增加而移动，但这种变化主要是由电子态密度的偏移驱动的，而非由排斥力本质属性的变化驱动。这意味着，一旦科学家了解了某种材料在某一压力下的电子结构，他们就可以通过追踪电子密度的变化，对该材料在其他压力下的行为做出相当准确的估计。这一洞察简化了寻找新超导体的过程，表明准确预测的关键不在于对排斥力进行复杂的、多变的猜测，而在于对电子景观进行精确的计算。

这项工作的意义不仅限于单一材料。研究人员证明，对于广泛的氢化物而言，标准的一刀切式电子排斥值是不够的。在具有高电子密度的材料中，使用旧的标准值会导致高估超导温度，从而可能误导研究人员关于哪些材料值得追求。通过直接计算排斥力，团队表明，几种领先候选材料的预测温度都下降到了更接近实验观测值的水平。例如，在硫化氢中，新的计算将预测温度降低了约40度，使其与实验室测得的数据相吻合。这表明，发现新型高温超导体的路径需要一种更严谨、针对特定材料的电子排斥力计算方法。

作者强调，他们的发现是基于经过已知金属和常见超导体验证的详细模拟，这使他们对结果充满信心。他们指出，尽管其方法计算量巨大，但足以作为未来寻找新材料的标准工具。该研究并不声称已经完全解开了超导之谜，但它消除了一个多年来困扰该领域的重大误差来源。通过用精确的计算取代粗略的猜测，研究人员为在复杂的高压氢化物世界中导航提供了一张更清晰的地图。其结果是对这些材料能力的更诚实的描绘，确保了寻找室温超导体的过程是由准确的物理学而非乐观的假设所引导。

技术摘要：基于第一性原理的高压氢化物增强型库仑伪势

问题陈述
在超导性的 Eliashberg 理论中，超导临界温度 ( $T_c$ ) 取决于声子介导的吸引作用与屏蔽库仑排斥作用之间的竞争。虽然利用密度泛函微扰理论 (DFPT) 可以常规地从第一性原理计算出电子-声子耦合强度 ( $\lambda$ )，但库仑排斥作用通常采用经验性的库仑伪势 ( $\mu^*$ ) 来处理。大多数文献依赖于来源于 Tolmachev-Morel-Anderson 结果的固定唯象值 $\mu^*$ （范围在 0.1–0.15 之间）。然而，这种方法假设电子态密度 (DOS) 和屏蔽库仑相互作用是恒定的，这一假设在复杂材料（尤其是高压氢化物）中可能并不成立。

依赖固定 $\mu^*$ 值导致了 $T_c$ 预测的显著差异。一个显著的例子是氢化物 Mg $_2$ IrH $_6$ ，两项独立的理论研究分别报告了 160 K 和 71 K 的 $T_c$ 值。作者将这一两倍的差异主要归因于对库仑相互作用的处理，指出标准的 $\mu^*$ 值可能会在具有高费米能级电子态密度 ( $N_F$ ) 的系统中显著高估 $T_c$ 。

方法论
作者提出了一个第一性原理框架，用于计算各种氢化物系统和压力下的库仑参数 ( $\mu$ ) 以及重整化伪势 ( $\mu^*_{ME}$ )。

计算工具： 本研究使用了 Mantel 代码，该代码实现了 Sanna 和 Pellegrini 的方法论。该代码通过接口与 Quantum ESPRESSO 进行电子结构计算（使用 PBE 伪势）以及与 Yambo 进行静态逆介电矩阵的评估。
理论框架：
- 在 随机相位近似 (RPA) 水平下评估 Bloch 态之间的屏蔽库仑相互作用 ( $W$ )。
- 无量纲库仑参数定义为 $\mu = N_F W(\varepsilon_F, \varepsilon_F)$ ，其中 $N_F$ 是费米能级的 DOS。
- 为了获得重整化伪势 $\mu^*_{ME}$ ，作者应用了包含 Matsubara 切断频率 ( $\omega_c$ ) 和电子能量尺度 ( $E_{el}$ ) 的标准重整化公式。
- 使用 可变态密度 (vDOS) 公式求解 各向同性 Eliashberg 方程 来计算超导 $T_c$ 。这种方法显式地保留了对 DOS 的能量积分，使其适用于在费米能级处存在尖锐峰值的系统，不同于 Allen-Dynes 公式中常用的恒定 DOS 近似。
验证： Mantel 代码针对金属元素（如 Al, Nb, Pd）和常见超导体进行了基准测试，显示出与文献值的一致性良好。研究还测试了 $\mu$ 对不同交换相关泛函（LDA, PBE, PBEsol, r2SCAN）的敏感性，结果表明 $\mu$ 对泛函选择的选择具有高度不敏感性，其变化主要由电子 DOS 的变化驱动。

关键结果
研究计算了一系列包括 LaH $_{10}$ , H $_3$ S, CaH $_6$ , YSbH $_6$ , Mg $_2$ IrH $_6$ , LaSc $_2$ H $_{24}$ 在内的氢化物在不同压力下的 $\mu$ 和 $\mu^*_{ME}$ 。

系统依赖的库仑排斥： 计算出的 $\mu$ 值表现出广泛的范围，从 0.14 (CaH $_6$ ) 到 0.62 (Mg $_2$ IrH $_6$ )。
- Mg $_2$ IrH $_6$ ： 由于极高的 $N_F$ 以及导致屏蔽相互作用出现峰值的能隙，该系统表现出异常大的 $\mu = 0.62$ 。由此产生的 $\mu^*_{ME} \approx 0.39$ 远高于常规范围。使用该值相比于标准的 $\mu^*_{ME} = 0.16$ ，会将预测的 $T_c$ 降低约 40 K，从而解决了之前理论估计（160 K vs. 71 K）之间的差异。
- H $_3$ S： 该系统表现出屏蔽相互作用的平滑衰减。虽然乘积 $\mu = N_F W(\varepsilon_F, \varepsilon_F)$ 仍保持在适中的 $\mu = 0.30$ ，但由此产生的 $T_c$ 与标准 $\mu^*$ 值相比被 有所抑制。在 140–240 GPa 的压力范围内，第一性原理计算得到的 $T_c$ 值比使用固定 $\mu^* = 0.1$ 得到的 $T_c$ 值 低 35 到 50 K。
- YSbH $_6$ ： 该系统遵循大 $N_F$ 的趋势，产生较大的 $\mu = 0.33$ 。强烈的库仑相互作用因此将 $T_c$ 从文献 [56] 报道的 118 K 降低到 100 K 左右。
- LaH $_{10}$ 和 CaH $_6$ ： 这些系统的 $\mu$ 值较低（分别为 0.18 和 0.14），证实了标准 $\mu^*$ 范围适用于这些系统。
态密度 ( $N_F$ ) 的作用： 观察到一种强相关性，即 $\mu$ 紧密追踪 $N_F$ 。虽然较大的 $N_F$ 会增强极化率和屏蔽作用（从而降低裸相互作用 $W$ ），但定义 $\mu$ 中的线性前因子 $N_F$ 起到了主导作用，导致库仑参数净增加。因此，那些拥有高 $N_F$ 的材料——通常由于具有强电子-声子耦合而被视为高- $T_c$ 超导性的最有希望的候选者——也恰恰是标准 $\mu^*$ 值最不可靠的系统。
压力依赖性： 研究发现屏蔽库仑相互作用 $W$ 对压力基本不敏感。 $\mu$ 随压力的变化主要由 $N_F$ 的变化驱动。对于 H $_3$ S，随着压力从 140 升至 240 GPa， $\mu$ 从 0.24 适度增加到 0.30，导致在使用第一性原理相互作用时，相比于固定 $\mu^*$ ，其 $T_c$ 会受到额外的抑制。
非谐性： 当在 H $_3$ S 中引入非谐声子频率时，由于电子-声子耦合减弱，预测的 $T_c$ 显著下降（下降约 70 K）。非谐性与第一性原理库仑相互作用的结合使得理论预测与实验数据吻合良好，尽管在这一特定案例中，完全能量依赖的 $W$ 处理方法与费米面平均 $\mu$ 方法相比略微低估了 $T_c$ 。

意义与主张
本文认为，在许多高压氢化物中，使用固定的、经验性的 $\mu^*$ (0.1–0.15) 对于准确预测 $T_c$ 是不足够的。作者证明了：

材料特异性： 库仑相互作用具有很强的材料依赖性，特别是对于具有大 $N_F$ 或特定轨道特征（例如局域化的 d 和 f 态）的系统。
解决分歧： 第一性原理计算 $\mu$ 可以解决 $T_c$ 预测中的巨大差异，正如 Mg $_2$ IrH $_6$ 的案例所示。
可行性： 尽管计算 $\mu$ 的计算成本随系统规模缩放，但它对于高通量筛选是可行的，且对于简单系统而言，其成本与典型的声子计算相当甚至更低。
建议： 作者建议，在 ab initio 超导性预测中，应将纳入第一性原理库仑相互作用作为一种标准做法，特别是对于那些电子结构特征（高 $N_F$ 、费米能级附近的能隙）会导致增强库仑排斥作用的氢化物，而标准模型无法捕捉这些特征。

该研究承认了一些局限性，例如使用静态 RPA 介电函数（忽略了交换相关顶点修正和动态频率依赖性），这可能会导致在低密度系统中低估 $\mu$ 。然而，在 RPA 框架内，其结果为预测超导特性提供了一个比唯象方法更严谨、更准确的基准。

技术摘要：基于第一性原理的高压氢化物增强型库仑伪势

类似论文