$d$-wave superconductivity in quintuple-layer bulk superconductor La$_6$Ni$_5$O$_{12}$

技术摘要：五层体相超导体 La$_6$Ni$_5$O$_{12}$ 中的 d 波超导电性

问题陈述
近期发现的五层镍氧化物 La$_6$Ni$_5$O$_{12}$ (6512) 体相超导电性，为解决降低阶梯型 Ruddlesden–Popper (RP) 镍氧化物的配对对称性提供了关键契机。虽然无限层镍氧化物 (112) 已显示出存在节点超导性的证据，但由于这些体系的超导性主要在超晶格薄膜中实现，确定的表征工作一直受到阻碍。这种局限性限制了常规体相敏感探测手段的应用，如核磁共振 (NMR)、穿透深度测量以及热力学研究。此外，尽管关于无限层和有限层镍氧化物的理论研究主要倾向于由反铁磁自旋涨落介导的自旋单态 $d_{x^2-y^2}$ 波配对，但对于配对趋势如何随整个降低阶梯型 RP 系列（从无限层到有限层化合物）演变的理解仍不完整。

方法论
作者采用结合密度泛函理论 (DFT) 与变分蒙特卡洛 (VMC) 计算的现实理论框架，来研究 La$_6$Ni$_5$O$_{12}$ 的基态。

1. 模型构建： 基于 DFT 计算，作者构建了一个二维五能带紧束缚 (TB) 模型，重点关注 Ni $d_{x^2-y^2}$ 轨道。该模型处理五个 NiO$_2$ 层（标记为外–内–中心–内–外），并包含一个作为自掺杂库的 La 衍生能带。层间跳跃参数化以匹配多层铜氧化物中的动量依赖性。
2. 相互作用与变分方法： 系统使用具有原位 Hubbard 斥力 $U$ 的五层 Hubbard 哈密顿量进行建模。利用带有 Gutzwiller 投影波函数 ($|\Psi\rangle = P_G |\Phi_{mf}\rangle$) 的 VMC 研究基态。变分参数包括层间跳跃振幅、层依赖化学势以及层内配对振幅。
3. 对比分析： 为了确定受 Favored 的配对对称性，作者比较了三种平均场拟设：$d_{x^2-y^2}$ 波态、扩展 $s$ 波态以及无配对参考态。该框架同样被扩展到三层 La$_4$Ni$_3$$_8$ (438) 和无限层 LaNiO$_2$ (112)，以考察整个家族的趋势。
4. 可观测物理量计算： 利用优化的变分态，作者计算了实验可观测物理量，包括隧穿态密度 (TDOS)、超流密度、奈特移 (Knight shift) 以及核自旋-晶格弛豫率 ($1/T_1T$)，并引入 Gutzwiller 约化因子以考虑相关效应。

核心结果

• 6512 中的配对对称性： VMC 计算表明，自旋单态 $d_{x^2-y^2}$ 波态在能量上优于非超导参考态。相比之下，扩展 $s$ 波配对振幅受到强烈抑制，比 $d$ 波振幅弱约两个数量级。层分辨电子密度在五个层中保持近乎均匀。
• 实验特征： 计算出的 $d$ 波态的可观测物理量为洁净节点超导体提供了一致的特征：
  • TDOS： 表现出特征性的 V 型低能谱以及线性消失的态密度，反映了其节点结构。
  • 超流密度： 在低温下表现出线性耗尽 ($\rho_s(0) - \rho_s(T) \propto T$)，这归因于靠近能隙节点的热激发准粒子。
  • 奈特移： 在 $T_c$ 以下随温度呈现近似线性下降，符合自旋单态配对特征。
  • 弛豫率： 自旋-晶格弛豫率遵循 $1/T_1T \propto T^2$ (等价于 $1/T_1 \propto T^3$)，这是洁净节点超导体的标志。
• 跨家族比较：
  • 无限层 (112)： 该框架证实了稳健的 $d_{x^2-y^2}$ 波配对趋势，与之前的研究一致。
  • 三层 (438)： 在其标称空穴填充度 ($x=1/3$) 下，该体系处于 $d$ 波超导态与非超导态的边界附近。扩展 $s$ 波振幅可以忽略不计。作者认为，电子掺杂（减少空穴浓度）可能会稳定该化合物的超导性。

意义与主张
本文将体相 La$_6$Ni$_5$O$_{12}$ 确立为一个极具前景且实验可及的平台，用于确定镍氧化物超导体的配对对称性。通过直接预测体相敏感探测手段（如穿透深度和 NMR）的响应，本研究为明确表征序参量提供了一条路径，而这在薄膜体系中一直难以实现。

作者主张，其结果表明在降低阶梯型 RP 镍氧化物家族中存在共同的 $d_{x^2-y^2}$ 波配对趋势，且超导强度受载流子浓度和层结构的控制。该研究将体相 6512 定位为连接已建立的无限层 112 超导性与新兴有限层化合物（如 438）超导机制之间的桥梁。这项工作并非提出新的实验装置，而是旨在阐释现有体相探测手段在验证该类新材料中 $d$ 波对称性理论预测方面的潜力。