Vacancy-free cubic superconducting NbN enabled by quantum anharmonicity

本研究证明，量子非谐效应稳定了一种先前未知的、无空位的化学计量比 NbN 立方相，其超导转变温度为 20 K，这挑战了空位对其结构稳定性至关重要的长期观点。

要点

想象你有一支舞者团队（原子），试图在舞池上组成一个完美、刚性的正方形队形。多年来，科学家们一直认为，对于以超导能力（零电阻导电）闻名的氮化铌（NbN）而言，这种完美的正方形队形无法维持在一起。

旧的故事是这样的：为了防止舞者互相绊倒并导致队形崩塌，你必须在舞池上留一些空位（空位）。你需要移除少数舞者才能使正方形稳定。如果你试图填满每一个位置（1:1 比例），队形就会摇晃并分崩离析。

新发现：“量子抖动”

这篇论文讲述了一个不同的故事。研究人员发现，如果你不再将舞者视为僵硬、冻结的雕像，而是意识到它们实际上是量子粒子，整个图景就会改变。

在量子世界中，原子并非静止不动；即使在绝对零度，它们也持续地颤动和振动。这被称为“零点运动”。此外，将它们结合在一起的力不像简单弹簧那样均匀回拉；它们是“非谐的”，意味着当被强力拉伸时，弹簧会变得怪异且富有弹性。

作者利用超级计算机模拟了这些“量子抖动”和“弹性弹簧”。他们发现，当原子被允许伴随这些量子运动起舞时，它们并不需要空位来保持稳定。相反，它们自然地转变为一种新的、略微扭曲的形状，这种形状实际上比旧的完美正方形更加稳定。

隐喻：颤动的果冻

将旧的“完美正方形”结构想象成一块太僵硬而无法站立的果冻块；它会坍塌。科学家过去认为，你必须给果冻戳出洞（空位）才能让它保持形状。

这篇论文表明，如果你让果冻颤动（量子非谐性），它就不会坍塌。相反，这种颤动导致果冻 settle 到一个略微被压扁、摇晃的形状，这个形状实际上比刚性块更坚固、更舒适。这种新形状就是作者发现的“无空位”立方相。

他们的发现

1. 一种新形状：他们确定了一种特定的、此前未知的原子排列（空间群称为 $P\bar{4}3m$）。这就像舞者们找到了一个新的、略微偏离中心的队形，比完美的正方形更有效。
2. 它更稳定：这种新的、摇晃的形状在能量上更“快乐”（能量更低），优于旧的“完美”正方形形状，即使没有任何舞者缺失。
3. 超导性能：他们计算了这种新形状在无电阻情况下导电的能力。他们发现它在20 开尔文的温度下工作。这与现实中近乎完美（近化学计量比）样品实验观察到的结果非常吻合。
4. 旧数学为何失败：先前的计算机模型假设原子是刚性弹簧（谐性）。那些模型认为完美正方形是不稳定的。当研究人员加入“量子抖动”（非谐性）后，数学终于与现实达成一致：完美正方形可以存在，但它只需要略微扭曲才能保持站立。

核心结论

长期以来，科学家们认为需要缺陷（缺失的原子）才能使立方氮化铌发挥作用。这篇论文论证了其实并不需要。我们在实验中看到的“缺陷”可能只是因为我们未能理解原子的量子舞步。如果我们能够合成这种完美的、无空位的材料，它作为超导体的表现可能实际上比我们目前认为的还要好。

这篇论文表明，与其试图通过添加或移除原子来修复材料，我们可能只需要让原子进行其自然的量子舞蹈，以找到它们最稳定、性能最高的形状。

技术摘要

技术摘要：量子非谐性实现的无空位立方超导 NbN

问题陈述
氮化铌（NbN）是一种技术关键材料，以其高临界磁场和超导转变温度（$T_c$）而闻名，其面心立方（fcc）$\delta$相（$Fm\bar{3}m$）的$T_c$可达 17 K。然而，理论预测表明，理想的 1:1 化学计量比$\delta$-NbN 相在谐性近似下是动力学不稳定的，在布里渊区的大部分区域表现出虚声子模。因此，实验观察一致表明，氮空位对于稳定立方结构是必要的。这些空位诱导的结构畸变虽然稳定了该相，但也改变了其电子性质。挑战在于确定是否存在稳定的、无空位的、1:1 化学计量比的立方相，如果存在，其超导性质是否与空位稳定的$\delta$相有显著差异。传统的计算方法难以解决这一问题，因为模拟无序和空位需要大得令人望而却步的超胞，而标准的谐性近似无法捕捉理想结构的真实稳定性。

方法论
作者采用了一种多方面的计算方法，将最先进的第一性原理计算与机器学习相结合，以克服标准密度泛函理论（DFT）和谐性近似的局限性。

1. 量子非谐性：为了考虑量子离子运动（零点运动）和离子势的非抛物线性质，研究采用了随机自洽谐性近似（SSCHA）。该方法生成一个显式包含非谐效应的有效谐性描述。
2. 机器学习势函数：鉴于 SSCHA 和分子动力学（MD）模拟的计算成本（需要数十万次 DFT 计算），作者训练了高保真度的机器学习原子间势函数（具体为矩张量势 MTP 和瞬时数据衍生势 EDDP）。这些势函数保留了 DFT 级别的能量和力精度，实现了约$10^4$–$10^5$倍的计算加速。
3. 互补验证：使用三种不同的方法交叉验证了所识别相的稳定性：
   • SSCHA：用于弛豫结构并计算非谐声子色散。
   • 分子动力学谱能量密度（MD-SED）：直接从平衡 MD 模拟计算振动性质，捕捉所有阶次的非谐性。
   • 从头算随机结构搜索（AIRSS）：利用“抖动”超胞和对称性约束搜索，无偏地独立搜索低能结构。
4. 超导性质：使用各向同性 Migdal-Eliashberg（ME）形式（通过 IsoME 代码）和各向异性 ME 计算（通过 EPW）计算电子 - 声子耦合和超导转变温度。屏蔽库仑相互作用（$\mu^*$）使用 GW 近似进行评估。

主要贡献与结果

• 新稳定相的发现：研究识别出一种此前未报道的动力学稳定的 NbN 立方相，空间群为$P\bar{4}3m$。当理想的 1:1 化学计量比$\delta$相在量子非谐效应下完全弛豫时，便出现了这种结构。
• 热力学稳定性：发现$P\bar{4}3m$相的自由能比理想的$Fm\bar{3}m$ $\delta$相低 65 meV/atom。虽然其能量仍高于六方$\epsilon$-NbN 相（预测的$x=1$时的基态），但$P\bar{4}3m$相在热力学上比传统的立方$\delta$相更有利。
• 稳定化机制：$P\bar{4}3m$相的稳定化主要由量子离子运动（零点涨落）驱动，而非高阶非谐效应。结构畸变消除了$Fm\bar{3}m$相的高简并度，特别是分裂了 X 点和$\Gamma$点的能带。
• 电子结构与耦合：与理想的$\delta$相相比，$P\bar{4}3m$相中的晶格畸变降低了费米能级处的电子态密度（DOS）。这种降低导致电子 - 声子耦合常数（$\lambda$）显著减弱，与那些通过大电子展宽人为稳定$\delta$相的计算结果相比尤为明显。
• 超导转变温度（$T_c$）：
  • 使用$P\bar{4}3m$相的 SSCHA 声子色散，计算出的$T_c$为20 K。
  • 该值与近化学计量比块体 NbN 的实验报道$T_c$值（约 16 K）高度吻合。
  • 相比之下，通过大电子展宽人为稳定$Fm\bar{3}m$相的计算将$T_c$高估至 27 K。
• NbN 的特异性：作者指出，将类似的 SSCHA 弛豫程序应用于相关化合物 TiN 和 NbC 并未导致对称性降低的相变，表明这种现象是 NbN 量子非谐景观所特有的。

意义与主张
本文挑战了长期以来认为氮空位对于稳定 NbN 立方相是绝对必要的假设。相反，它提出量子非谐性可以在理想的 1:1 化学计量比下稳定一种独特的、无空位的、畸变的立方相（$P\bar{4}3m$）。

作者声称，这一发现为 NbN 的合成提供了新视角。虽然之前的工作集中在调节空位浓度以定制材料性能，但研究结果表明，合成理想的 1:1 化学计量比立方相（特别是$P\bar{4}3m$结构）可能为获得增强的超导性能提供更直接的途径。无空位相计算出的 20 K 的$T_c$表明，实现完美的化学计量比有可能带来比目前在空位稳定样品中观察到的更高的临界温度。

作者最后鼓励进行进一步的实验研究，特别是利用高分辨率衍射技术对单晶或控制良好的粉末进行分析，以验证$P\bar{4}3m$相的存在并评估其超导应用潜力。他们还建议未来的研究应探讨应变和化学取代对该相稳定性的影响，并调查量子非谐性在其他超导材料中的作用。