Three-dimensional Ising superconductors designed via inversion-symmetry breaking in intercalated NbSe$_2$ and NbTe$_2$

技术摘要：通过插层 NbSe₂ 和 NbTe₂ 中的反演对称性破缺设计三维伊辛超导体

问题陈述
伊辛超导性（Ising superconductivity）的特征是面内上临界磁场（$H_{c2,\parallel}$）远超泡利顺磁极限（$H_P$），其依赖于强自旋-轨道耦合（SOC）与反演对称性破缺的共存。虽然这种现象在 1H-NbSe₂ 等二维（2D）过渡金属二硫化物（TMDs）单层中已得到充分证实，但在三维（3D）体相晶体中实现这一目标仍然具有挑战性。在大多数体相 TMDs 中，平衡层堆叠（例如 2H-NbSe₂ 中的 AB 堆叠）会恢复反演对称性，从而消除实现自旋-谷锁定所需的反对称自旋-轨道耦合（SOC）。尽管天然的非中心对称多型体（如 4H$_a$-NbSe$_2$）和特定的插层系统已展现出潜力，但仍缺乏一种针对广泛化学空间设计 3D 伊辛超导体的系统性策略。此外，除了单纯的 SOC 量级之外，决定伊辛保护效率的具体因素仍需定量澄清。

方法论
作者采用系统的第一性原理计算方法来筛选和设计候选材料：

• 结构设计： 使用四种插层元素（In, Sn, Pb, Bi）将它们插入 NbSe$_2$ 和 NbTe$_2$ 基体中，构建了十六种插层化合物。每种成分都以两种不同的多型体进行建模：非中心对称的 $P\bar{6}m2$（破缺反演对称性）和中心对称的 $P6_3/mmc$（恢复反演对称性），以便进行直接的对称性控制对比。
• 稳定性与电子结构： 使用带有 SOC 和范德华力校正的密度泛函理论（DFT）计算（VASP, PBE-GGA），通过声子谱（使用 Quantum ESPRESSO）确定结构稳定性。排除了动力学不稳定的系统。
• 超导性分析： 利用 EPW 代码计算电子-声子耦合（EPC），并通过 Allen-Dynes 修订的 McMillan 公式确定超导转变温度（$T_c$）。
• 临界磁场建模： 为了预测上临界磁场，作者基于最大局域化 Wannier 函数（MLWF）构建了 Bogoliubov-de Gennes (BdG) 哈密顿量。这使得计算磁场下的准粒子态密度（DOS）以及确定 $H_{c2,\parallel}$ 和 $H_{c2,\perp}$ 成为可能。同时求解了各向异性 Migdal-Eliashberg 方程，以分析动量分辨的超导能隙。

关键结果

1. 候选物识别： 在 16 个初始候选物中，有 11 个在动力学上是稳定的。其中，四个非中心对称（$P\bar{6}m2$）化合物被确定为可行的 3D 伊辛超导体：InNbSe$_2$, SnNbSe$_2$, PbNbSe$_2$, 以及 PbNbTe$_2$。其中心对称的反面化合物及 Bi 插层系统要么不具备超导性，要么表现出显著较低的 $T_c$ 且不具备伊辛保护。
2. 电子与拓扑性质：
   • 这四种候选物在费米能级附近表现出显著的 SOC 能带分裂（$\Delta_{SOC}$），量级为 80–100 meV，这主要由 Nb-4d 态驱动。
   • 自旋纹理分析确认了在布里渊区（$k_z=0$ 和 $k_z=1/2$）的 $K$ 和 $K'$ 谷处存在稳健的面外自旋极化，这是伊辛超导性的标志。
   • $P\bar{6}m2$ 相由镜像对称性保护的拓扑节点线（topological nodal lines）以及伴随的鼓状表面态（drumhead-like surface states）组成，类似于 PbTaSe$_2$。
3. 超导参数：
   • 计算出的 $T_c$ 值范围为 2.6 K 至 5.4 K。
   • BdG 计算预测面内上临界磁场（$H_{c2,\parallel}$）可达到泡利极限（$H_P$）的 4 到 7 倍。
   • 对于 $P\bar{6}m2$ 相，各向异性比率 $H_{c2,\parallel}/H_{c2,\perp}$ 非常显著（5.5–10.5），而中心对称的 $P6_3/mmc$ 相则表现出近乎各向同性的行为。
4. 保护机制： 一个关键发现是，伊辛保护的效率并非仅由 $\Delta_{SOC}$ 的大小决定。相反，它受控于费米面上面外自旋极化的纯度。具有较大 $\Delta_{SOC}$ 的化合物（例如 Pb 插层系统）由于存在较大的残余面内自旋分量（$S_{x,y}$），导致其 $H_{c2,\parallel}/H_P$ 比率较低，这些分量在面内磁场下提供了寄生性的退相干通道。
5. 能隙各向异性： 各向异性 Eliashberg 计算揭示了高度各向异性的超导性。SnNbSe$_2$ 和 PbNbSe$_2$ 表现出明显的双能隙特征，较大的能隙出现在 $M-K$ 口袋（强 SOC），较小的能隙出现在 $\Gamma$ 口袋。

意义与主张
本文确立了插层法是一种在 TMDs 中实现 3D 伊辛超导体的实用且可调的对称性工程策略。通过系统地筛选插层化学成分和堆叠多型体，作者提供了一个设计能够在体相材料中克服平衡堆叠对称性限制的框架。

该研究声称提供了：

• 一组特定的四种候选材料（InNbSe$_2$, SnNbSe$_2$, PbNbSe$_2$, PbNbTe$_2$），并预测了其 $T_c$ 和稳健的 $H_{c2,\parallel}$。
• 一项定量见解，即自旋纯度而非仅仅是 SOC 量级，才是决定伊辛保护效率的关键因素。
• 证明了 3D 伊辛超导性可以在体相 TMDs 中实现，其临界磁场显著超过泡利极限，弥合了脆弱的 2D 单层与稳健的体相晶体之间的差距。

作者总结认为，这些发现既提供了候选材料，也提供了一个筛选框架，用以指导未来实现场强稳健型 3D 伊辛超导性的实验工作。