Strong-coupling anisotropic superconductivity in hexagonal HfRuAs from anisotropic Migdal-Eliashberg theory

本研究采用各向异性 Migdal-Eliashberg 理论并结合第一性原理计算，证明六方 HfRuAs 是一种强耦合、声子介导的超导体，其特征为单一各向异性$s$波能隙，且临界温度与实验观测结果一致。

要点

想象一个电力毫无阻力流动的世界。这就是超导性的魔力。长期以来，科学家们一直试图确切理解某些材料究竟是如何施展这一魔法的。

本文深入探讨了一种名为六方 HfRuAs（一种由铪、钌和砷组成的晶体）的特定材料。研究人员利用强大的计算机模拟，弄清了这种材料为何会成为超导体，以及它如何表现。

以下是他们发现的简要概述，分解为简单概念：

1. “舞池”与“音乐”

在这种材料中，电子就像拥挤舞池中的舞者。通常情况下，它们会相互碰撞并损失能量（即电阻）。但当材料冷却到足够低的温度时，它们开始配对，完美同步地起舞。

• 音乐（声子）： 本文解释说，促使这些电子起舞的“音乐”实际上是原子本身的振动。将原子想象成在蹦床上跳跃的人。当他们跳跃时，会形成波浪。
• 强连接： 研究人员发现，跳舞的电子与跳跃的原子之间的连接极其紧密。这并非轻描淡写的轻拍，而是一次坚定的握手。在科学术语中，这被称为"强耦合"。这种连接的强度测量值约为1.56，远高于普通超导体中的数值。

2. “重”舞者与“轻”舞者

该材料拥有不同“层”或电子层（称为费米面）。研究发现，“音乐”并非在所有地方都同等播放：

• 低音： 最重要的振动是缓慢的低频振动。这些振动主要由沉重的铪和钌原子晃动引起。
• 各向异性（失衡的舞蹈）： 舞蹈并非在所有方向上都相同。在电子“舞池”的某些部分，与“音乐”的连接非常强，而在其他部分则较弱。这就像一个舞池，中心区域音乐响亮清晰，而边缘则变得沉闷。这种不均匀性被称为各向异性。

3. 能量中的“能隙”

要成为超导体，电子需要在能级上打开一个“能隙”——一个防止它们受到干扰的保护屏障。

• 单一且波动的护盾： 本文发现，这种材料拥有一个主护盾（单一能隙），而非多个不同的能隙。然而，由于前述的“失衡”舞蹈，这个护盾并非一个完美均匀的圆形。它更像是一个略微被压扁或波动的圆形。
• 无孔洞： 至关重要的是，该护盾是完全闭合的。护盾本身没有孔洞或缝隙。这意味着超导性非常稳定，并遵循经典的"s 波”模式（一种标准且安全的超导类型）。

4. 温度谜题

研究人员计算出，这种材料应在约16 开尔文的温度下成为超导体（非常冷，但并非那么冷）。

• 差异： 现实世界的实验显示，该材料在较低温度下（4 K 到 7 K 之间）才表现出超导性。
• 为何存在差异？ 本文指出，计算机模型代表的是没有缺陷的“完美”晶体。现实世界的样品可能含有微小的杂质、缺陷或混合相，它们就像“减速带”一样，减缓了超导性，并降低了其发生的温度。

5. 主要结论

主要结论是：六方 HfRuAs 是一种“强耦合”超导体。

• 类比： 如果弱耦合超导体就像两个人手牵手轻轻行走，那么强耦合超导体就像两个人紧紧相拥，作为一个整体移动。
• 证据： 能隙与温度的比率远高于弱超导体的标准极限，证明了电子与振动原子之间的“拥抱”非常紧密。

总之： 本文利用高级数学证明，HfRuAs 是一种由自身原子强烈振动驱动的稳健超导体。虽然现实世界的样品不如计算机模型预测的那样完美，但基础物理揭示了一种电子与原子以惊人强度共舞的材料。

技术摘要

技术摘要：基于各向异性 Migdal-Eliashberg 理论的六方 HfRuAs 强耦合各向异性超导电性

问题与动机
三元等原子金属间化合物 $TT'X$（其中$T, T'$为过渡金属，$X$为磷族元素或硅）常表现出超导电性，特别是在六方$ZrNiAl$型（$h$相）结构中。尽管该家族成员如$HfRuP$和$ZrRu(As, Si, P)$的超导转变温度（$T_c$）超过 10 K，但近期关于六方$HfRuAs$（$h$-$HfRuAs$）的实验报道显示$T_c$存在显著差异（根据合成条件不同，范围在 4.3 K 至 7.9 K 之间）。此前针对相关化合物的理论研究也得出了截然不同的$T_c$预测值。此外，虽然实验比热数据表明存在弱耦合、全能隙的$s$波态，且比热跃变减小（$\Delta C/\gamma T_c \approx 1.03$），但针对$h$-$HfRuAs$ 中动量依赖的电子 - 声子耦合（EPC）及各向异性超导（SC）能隙结构的全面理论理解仍然缺失。具体而言，需要全各向异性 Migdal–Eliashberg（ME）计算来解析多带超导特性及费米面（FS）依赖的配对相互作用，这些是标准各向同性近似可能遗漏的。

方法论
作者利用以下工作流程对 $h$-$HfRuAs$ 进行了全面的从头算研究：

1. 电子与晶格性质：在 Quantum ESPRESSO 包内，利用密度泛函理论（DFT）和密度泛函微扰理论（DFPT）计算了电子结构、声子色散及电子 - 声子耦合矩阵元。采用了 PBEsol 交换关联泛函及优化的范守制范德华赝势。
2. Wannier 插值：为实现 EPC 计算所需的密集采样，利用 EPW 代码通过最大局域化 Wannier 函数，将粗布里渊区网格上计算的电子 - 声子矩阵元插值到精细的 $k$ 点和 $q$ 点网格（32 $\times$ 32 $\times$ 48）上。
3. 各向异性 Migdal-Eliashberg 形式：通过在虚频轴上求解全各向异性 ME 方程来研究超导特性。计算利用了各向异性 EPC 核及屏蔽库仑赝势（$\mu^*_c = 0.15$），该选择基于费米能级处 $Hf$-$d$ 和 $Ru$-$d$ 轨道的主导地位。通过解析延拓获得实频轴上的超导能隙函数，以计算准粒子态密度（DOS）。

主要结果

• 强电子 - 声子耦合：计算得到的 Eliashberg 谱函数（$\alpha^2F(\omega)$）显示总 EPC 常数 $\lambda \approx 1.56$。该数值将 $h$-$HfRuAs$置于强耦合机制中，显著高于相关化合物如$h$-$ZrRuAs$（$\lambda \approx 0.79\text{--}1.32$）和$h$-$HfRuP$（$\lambda \approx 0.87$）。耦合主要由低频声子模式（$< 10$ meV）主导，这些模式主要与 $Hf$和$Ru$ 的振动相关，贡献了约 82% 的总 EPC。
• 各向异性多带超导性：费米面由三个片层组成：两个类空穴带（$H1, H2$）和一个类电子带（$E1$）。动量分辨的 EPC（$\lambda_{nk}$）和超导能隙（$\Delta_{nk}$）在这些片层上表现出显著的各向异性，其中最大的变化发生在类空穴 $H2$ 带上。
• 能隙结构与对称性：尽管存在显著的各向异性，超导态的特征是具有单一各向异性能隙，整体呈现 $s$ 波对称性。能隙在所有费米面片层上完全打开（无节点），这由 U 型的低能准粒子态密度所证实。能隙幅值集中在 $\Delta \approx 2.9$ meV 左右，展宽约为 $\sim 0.8$ meV。
• 转变温度与能隙比：计算得到的 $T_c$ 约为 16.06 K（基于平均能隙的拟合 BCS 型温度依赖性），高于大多数实验报道，但在数量级上与先前理论预测的上限一致。能隙比 $2\Delta(0)/k_B T_c \approx 4.25$ 显著超过了弱耦合 BCS 极限值 3.53，进一步证实了超导性的强耦合本质。
• 库仑赝势敏感性：研究系统地将 $\mu^*_c$ 从 0.11 变化至 0.30。虽然增加 $\mu^*_c$ 会抑制 $T_c$（从 13.9 K 降至 6.5 K）并减小能隙幅值，但能隙在动量空间的分布及其各向异性特征依然稳健，表明关于强耦合和各向异性的定性结论不依赖于 $\mu^*_c$ 的具体选择。

意义与主张
本文确立了 $h$-$HfRuAs$为一种声子介导的强耦合各向异性超导体。其主要贡献在于证明了该材料的超导性是由电子与$Hf$和$Ru$ 的低频晶格振动之间的强耦合驱动的，从而导致较大的 EPC 常数（$\lambda \approx 1.56$）。作者阐明，尽管该系统在 EPC 和超导能隙方面均表现出显著的多带各向异性，但它仍保持全开能隙的单能隙 $s$ 波对称性。计算得到的 $T_c$（16 K）与实验值（4.3–7.25 K）之间的差异归因于实验样品中的外在因素（无序、相共存、化学计量比偏差），而非本征理论模型的失效。该工作提供了关于动量依赖的电子 - 声子相互作用在决定此类化合物超导特性中所起作用的详细微观描述，为未来针对高质量单晶的研究提供了一致的理论基准。