Ferroelectric quantum critical point in superconducting hydrides: The case of H$_3$S

通过利用路径积分分子动力学结合机器学习势函数，本研究揭示了 H$_3$S 中的高温超导性发生在一个由巨大的核量子涨落所主导的顺电区域，该区域位于约 134 GPa 处的铁电量子临界点之上，而该临界点属于 4D Ising 普适类。

要点

想象你有一个由硫原子和氢原子组成的微小、超高密度的球体。在极高的压力下，这个球体变成了一种超导体——一种电阻为零的材料。长期以来，科学家们一直对为什么这种名为 H3S 的特定材料会发生这种现象感到困惑。他们知道它在特定的压力（约 155 GPa）下表现最好，但关于原子如何行为的地图却一直缺失。

这篇论文就像是在绘制那张缺失的地图。研究人员使用了一种强大的计算机模拟，来追踪原子是如何“跳舞”的，不仅是将它们视为固体的球，而是视为“概率云”（一种量子效应）。以下是他们的发现，用简单的语言解释如下：

1. “模糊”的原子与神奇的压力

在微观原子的世界里，事物并非是实心的，它们在不停地抖动和跳动。研究人员发现，在特定压力（约 134 GPa）下，H3S 中的氢原子达到了一个“临界点”。

• 类比： 想象一个放在碗里的球。如果碗很深，球就会留在中间。如果你摇晃这个碗（加热）或挤压它（加压），球可能会开始到处滚动。
• 发现： 在这个被称为**量子临界点（QCP）**的临界点上，原子处于一种极度混乱的状态。它们既没有稳定在一个位置，也没有完全随机。它们正在进行剧烈的“涨落”，就像一群试图决定转向哪边的人群一样。

2. 相变：从“对称”到“不对称”

这种材料可以以两种主要形状（相）存在：

• “完美平衡”相（顺电相）： 氢原子位于硫原子的正中间。它是对称的，就像一个完美平衡的跷跷板。
• “不对称”相（铁电相）： 氢原子被推向了一侧。跷跷板倾斜了。

论文显示，从“平衡”到“不对称”的转变并不正好发生在超导性最强的压力处。相反，超导性的峰值出现在“平衡”区域，但紧邻着原子抖动最剧烈的那个临界点。

3. “超导甜点区”

这里有一个巨大的惊喜：

• 旧观点： 科学家们认为超导峰值是因为材料正在从平衡状态切换到不对称状态。
• 新发现： 论文表明，峰值实际上发生在平衡区域，但紧挨着混沌边缘。
• 类比： 想象一名冲浪者。最好的浪潮既不是平静的水面，也不是混乱的破碎浪花。最好的浪潮正是海洋开始变得汹涌的地方。这种“粗糙度”（量子涨落）有助于电子配对并无阻力地流动。论文表明，临界点附近氢原子的剧烈抖动起到了增强超导性的作用。

4. “4D Ising” 规则手册

研究人员分析了这一临界点背后的数学逻辑，发现它遵循一个被称为 4D Ising 普适类 的特定规则手册。

• 类比： 想象不同的游戏（如国际象棋、跳棋或围棋）。尽管它们看起来不同，但可能都遵循相同的底层逻辑，即棋子的移动方式。研究人员发现，这些原子的行为方式遵循与一种用于描述事物状态变化的复杂四维数学模型相同的“逻辑”。这证实了他们的发现是一个基本的物理定律，而非偶然现象。

5. 为什么这对于“地图”至关重要

在这项研究之前，H3S 的地图是模糊的。研究人员使用了一种新型的“AI 大脑”（机器学习势函数）来运行那些用旧方法无法完成的高成本模拟。

• 他们发现，如果你忽略原子的量子“模糊性”（将它们视为坚硬的台球），你得到的地图就是错误的。量子抖动会将转变压力移动巨大的数值（约 50 GPa）。
• 通过包含这些抖动，他们终于定位了“临界点”（QCP），并展示了超导峰值就位于强量子涨落的区域，即临界点之上。

总结

论文揭示了 H3S 中的超导魔力并非仅仅由材料形状的变化引起。相反，它发生在该材料徘徊在“量子临界点”附近时，此时原子正在剧烈振动。这些剧烈的振动就像一种催化剂，帮助电流完美地流动。研究人员现在精确地绘制了发生这种情况的位置，并证明了它遵循一个特定的、普适的数学规则。

技术摘要

问题陈述
硫化氢（H$_3$S）是一种重要的超高温超导体，在约 155 GPa 的压力下，其临界温度（$T_c$）可达 203 K。虽然其超导穹顶表现出类似于铜氧化物和掺杂 SrTiO$_3$ 的峰值特征，但其相图的精确性质，特别是超导最大值与结构相变之间的关系，仍未得到解决。先前的理论研究表明，由氢位移驱动的从体心立方 $Im\bar{3}m$ 相到三方 $R3m$ 相的铁电相变存在。然而，这些理论将该转变压力预测得显著偏低（低了超过 40 GPa），未能与实验观测到的 $T_c$ 峰值相匹配。此外，目前仍缺乏对 H$_3$S 在扩展温度和压力范围内的全面理解，包括热效应与核量子效应（NQE）之间的相互作用。

方法论
为了填补这些空白，作者采用了路径积分分子动力学（PIMD）来精确处理核量子效应。鉴于将 PIMD 与 ab initio 电子结构方法（如密度泛函理论 DFT）结合用于所需系统规模和温度时所面临的极高计算成本，作者利用了机器学习原子间势（MLIP）。具体而言，他们基于 DFT-BLYP 配置训练了一个 MACE（高阶等变消息传递神经网络）势函数，配置的压力范围为 40 至 220 GPa，温度为 100 K 和 200 K。

模拟使用 i-PI 软件包进行，采用 NVT 系综，系统规模为 $L \times L \times L$ 原胞（$L \in \{2, 3, 4\}$），以进行有限尺寸缩放分析。研究涵盖了 50 K 到 300 K 的温度范围。关键观测指标包括：

• 序参数： 全局质子位移（$\Delta$）及其绝对值（$\Delta_{abs}$），用于区分顺电相与铁电相。
• 局部观测量： 局部质子位移分布和局部偶极矩。
• 动力学性质： 用于提取非谐性声子频率并研究延迟效应的虚时声子格林函数（$G_{ij}(\tau)$）。
• 缩放分析： 临界点附近序参数的有限尺寸缩放，以确定普适类。

主要结果

1. 相图与量子临界点（QCP）： 研究揭示了一条随温度降低而向低压方向移动的铁电转变线。通过外推至热力学极限和 $T=0$ K，作者确定了一个铁电量子临界点（QCP），其位置为 $p_{QCP} \approx 134 \pm 2$ GPa。
2. 普适类： 临界点附近序参数的有限尺寸缩放分析得出临界指数 $\nu \approx 0.483 \pm 0.014$。假设洛伦兹不变性（$z=1$），该值与 4D Ising 普适类（其中 $\nu = 1/2$）一致。
3. 与实验 $T_c$ 的差异： 计算出的铁电转变线（例如在 200 K 时为 124 GPa）显著低于实验观测到的超导峰值 155 GPa。因此，实验 $T_c$ 最大值落在顺电 $Im\bar{3}m$ 相内，而非铁电相。
4. 量子涨落的作用： QCP 周围以及延伸至顺电相的区域具有显著的核量子涨落。作者通过量化虚时 $\beta/2$ 处的局部格林函数方差与 $\tau=0$ 处方差的比值（$\sigma^2_g(\beta/2)/\sigma^2_g(0)$）来评估这一现象。该比值在靠近 QCP 的顺电相中达到峰值，表明存在强烈的延迟效应和波动的偶极矩，而这些效应在质子冻结的铁电相中会被抑制。
5. 与经典动力学的比较： 经典分子动力学（MD）模拟显示，铁电转变发生的压力远高于量子情况（在 200 K 时偏移了约 50 GPa），这凸显了核量子效应（NQE）在稳定顺电相和降低转变压力方面的关键作用。
6. 声子软化： 驱动转变的光学声子模在系统从顺电侧接近转变时表现出软化，随后在进入铁电相时出现频率骤跳，这与二阶相变特征一致。

意义与主张
本文确立了 H$_3$S 中的高温超导性并非发生在铁电相变本身，而是发生在由强量子涨落和延迟效应主导的 QCP 附近的顺电区域。作者认为，这些涨落（可能由接近铁电 QCP 所介导）可能会增强超导性。这一情景意味着标准的 Migdal-Eliashberg 理论（通常忽略完整的频率依赖性和顶点修正）可能不足以描述 H$_3$S。铁电 QCP 的存在和强光学声子软化支持了超越线性阶的电子配对机制的可能性，类似于针对 SrTiO$_3$ 等稀薄金属提出的机制。这项工作提供了一个详细且具有量子力学准确性的相图，将超导穹顶的位置与其潜在的结构不稳定性联系起来，强调了在未来对超导氢化物进行理论处理时，包含核量子效应和非微扰电子-声子耦合评估的必要性。