Superconductivity in RbH$_{12}$ at low pressures: an \emph{ab initio} study

这项基于第一性原理的研究预测，当纳入晶格量子非谐效应时，RbH$_{12}$体系可在极低的压力（低至 10 GPa）下实现高温超导，其临界温度高达 111 K。

要点

想象一下，你正在试图建造一座无论外界多么寒冷都不会散热的房子。在物理学界，这被称为超导性：一种电流以零电阻流动的状态。几十年来，科学家们一直在努力寻找能够在“室温”（或者至少是在无需昂贵液氮即可轻松达到的温度）下实现这一状态的物质。

问题在于，迄今为止发现的最佳候选材料就像冰雕：只有当你用相当于整座山的重量（极端压力）去挤压它们时，它们才能发挥作用。一旦释放这种压力，它们就会崩塌并停止工作。

本文是一项计算研究（一种超先进的计算机模拟），旨在探讨：我们能否找到一种像超导体一样运作，却不需要头顶压着一座山就能保持稳定的材料？ 具体来说，研究人员考察了铷（一种软金属）和氢（最轻的元素）的混合物。

以下是他们研究发现的分解，使用了简单的类比：

1. “量子抖动”问题

在常规物理学中，我们想象原子静止地排列在整齐的网格中。但在原子层面，尤其是对于像氢这样的轻原子，由于量子效应，它们会不断地摇晃和振动。不要把这些原子想象成坚硬的弹珠，而要想象成富有弹性、不停抖动的果冻豆。

之前的研究将这些“果冻豆”当作坚硬的弹珠来处理。本文的研究人员意识到，要得到正确的答案，必须考虑到这些“果冻豆”正在剧烈摇晃的事实。他们使用了一种名为SSCHA（随机自洽谐波近似）的特殊数学工具，来模拟这种“摇晃”以及它如何改变材料的形状。

2. 寻找“金发姑娘”结构

研究人员在不同压力下模拟了铷 - 氢混合物（从 0 到 100 吉帕，这相当于最深海沟底部的压力，但还要高得多得多）。

他们发现了原子排列自身的五种不同方式（五种不同的“结构”）。

• 旧观点：如果不考虑“摇晃”，计算机显示只有两种结构是稳定的，而且仅在极高压力下才稳定。
• 新观点（包含摇晃）：当他们将“量子抖动”加入考量后，规则改变了。“摇晃”实际上有助于稳定这些结构。
  • 一种结构（Immm）在低至 25 GPa 的压力下变得稳定。
  • 另一种结构（P63/mmc）在低至10 GPa的压力下变得稳定。

为什么 10 GPa 是个大新闻？ 这就像发现了一座只需背着一个沉重的背包就能屹立不倒的房子，而不需要整座山。这是迄今为止预测到的此类二元超氢化物所需的最低压力。

3. “超导派对”

一旦确认这些结构可以存在，他们便问道：它们能完美导电吗？

• 答案：是的！他们发现的所有稳定结构都是金属性的（即能导电）。
• 温度：“派对”（超导性）在46 K 到 111 K（约 -227°C 到 -162°C）的温度下开始。
  • 虽然这还不是“室温”，但比这些材料通常所需的 -200°C 到 -270°C 要温暖得多。
  • 至关重要的是，研究人员发现氢原子的“摇晃”实际上有助于电子配对（超导的机制），就像一名指挥家，帮助电子更容易地协同起舞。

4. 如何识别它们（指纹）

由于这些材料难以制造，研究人员为实验学家（那些在实验室里实际制造这些东西的人）提供了一份“指纹”指南。

• X 射线衍射：他们模拟了 X 射线如何从这些结构上反弹。这就像用手电筒照射晶体；光线的图案能确切地告诉你原子处于什么形状。他们展示了不同的结构具有独特的图案，因此科学家不会将它们混淆。
• 拉曼光谱：他们还预测了如果用激光照射这些材料，它们会如何振动。这就像聆听材料的“嗡嗡声”来识别它。

核心结论

本文是一份路线图。它告诉实验科学家：“如果你将铷和氢在约 10 到 25 GPa 的压力下挤压在一起，并且考虑到氢原子是抖动的这一事实，你可能会发现一种在相对较低压力和较高温度下工作的超导体。”

它并不承诺明天就会出现新的电网，但它指明了一条通往未来的道路，在那里，我们或许不再需要庞大而昂贵的机器来维持超导体的存在。

技术摘要

技术摘要：低压下 RbH12 中的超导性：一项从头算研究

问题陈述
高压多氢化物已成为室温超导性的主要候选材料，但其实际应用因合成和稳定所需的百万巴级高压而受到严重限制。该领域当前的前沿是识别能够在环境或低压下保持动力学和热力学稳定的氢化物材料。尽管近期的预测主要集中在三元氢化物上，但缺乏关于能在低压下实现稳定性的二元氢化物的预测。此前关于铷氢化物（RbH$_x$）的计算研究表明其在相对较低的压力（例如 50 GPa）下具有动力学稳定性，但这些研究完全在密度泛函理论（DFT）框架内进行，忽略了量子零点运动和显著的非谐性影响，而众所周知，这些因素会极大地改变富氢系统的稳定性图景。

方法论
本研究采用第一性原理计算方法，考察了 RbH$_{12}$体系在 0 至 100 GPa 压力范围内的结构和超导性质。

• 结构预测： 使用 CrySpy 代码进行晶体结构预测，确定了五个竞争相：$Immm$、$C2/m$、$Cmcm$、$R\bar{3}m$和$P6_3/mmc$。
• 电子与振动计算： 使用 Quantum Espresso 软件配合 PBE 泛函进行标准 DFT 计算。为了克服谐波近似的局限性，作者采用了随机自洽谐波近似（SSCHA）。该方法通过最小化总吉布斯自由能，将离子动力学中的量子非谐效应纳入考量。
• 稳定性分析： 通过计算声子能带结构来评估动力学稳定性，分别使用了 SSCHA 辅助动力学矩阵和自由能的 Hessian 矩阵。研究具体考察了温度（高达 100 K）和量子效应对声子不稳定的影响。
• 超导性： 通过求解各向同性 Migdal-Eliashberg 方程估算超导临界温度（$T_c$）。电子 - 声子耦合常数源自密度泛函微扰理论（DFPT）。至关重要的是，$Immm$相在 25 GPa 下的重整化库仑相互作用参数（$\mu^*$）直接从第一性原理计算得出，结果为 0.118，随后将其应用于其他相。

关键结果

• 相稳定性： 在 DFT 层面，$Cmcm$相在 50 GPa 下最稳定。然而，通过 SSCHA 纳入量子非谐效应显著改变了能量图景。虽然相对焓的层级基本保持不变，但这些效应促进了特定相在更低压力下的动力学稳定性。研究发现，$Immm$相在低至 25 GPa 时具有动力学稳定性，而$P6_3/mmc$相在低至 10 GPa 时仍保持稳定。
• 动力学稳定性机制： 在谐波近似下，多个相在低压下表现出虚声子模（不稳定性）。SSCHA 计算表明，非谐性和量子效应（特别是在 100 K 时）稳定了这些模式，将$Immm$和$P6_3/mmc$相的动力学稳定性推向了二元氢化物前所未有的低压记录。
• 电子性质： 所有竞争相均为金属态。费米能级处的电子态主要呈现氢特征（>90%），这是高温超导性的关键指标。$Immm$相在费米能级附近的电子态密度（DOS）相对恒定，验证了使用标准近似进行$T_c$计算的合理性。
• 超导临界温度： 该研究预测动力学稳定相具有高温超导性。
  • 在 50 GPa 下，$T_c$值范围从 59 K（$P6_3/mmc$）到 111 K（$Immm$）。
  • 在较低压力下，25 GPa 处的$Immm$相表现出约 98 K 的$T_c$，而 10 GPa 处的$P6_3/mmc$相显示出 46 K 的$T_c$。
  • 电子 - 声子耦合常数（$\lambda$）被发现相对于其他超氢化物而言较低（25 GPa 处的$Immm$相 barely 超过 1），但高声子频率的补偿作用使其产生了高$T_c$值。
• 实验特征： 作者提供了模拟的 X 射线衍射（XRD）图谱和拉曼活性声子谱函数。他们指出，虽然$Immm$和$C2/m$难以通过 XRD 区分，但其他相（$Cmcm$、$R\bar{3}m$、$P6_3/mmc$）具有独特的特征。拉曼光谱显示出高频 H$_2$振动模的特征分裂，这有助于实验识别。

意义与主张
本文声称发现了一条潜在途径，可将二元超氢化物稳定在远低于此前认为可能的压力下。通过严格纳入量子非谐效应，该研究证明 RbH$_{12}$的$Immm$和$P6_3/mmc$相在低至 10–25 GPa 的压力下可能是亚稳态的。作者得出结论，这些发现为未来在低压条件下探索 Rb-H 二元化合物中的高温超导性提供了理论基础，有望弥合高压物理与实际应用之间的差距。该工作强调，忽略非谐性可能导致对稳定性的错误预测，而这些相的“创纪录低压”稳定化是量子效应的直接结果。