First-Principles Calculation of Superconducting $T_c$ in Superhard B-C-N Metals

本研究采用第一性原理计算预测，三元超硬金属 B$_2$C$_3$N 和 B$_4$C$_5$N$_3$ 在常压下展现出极具前景的超导转变温度，分别约为 40 K 和 20 K，这主要由其高德拜温度以及与已知合成化合物相当的形成能所驱动。

要点

想象一下，你正试图为电力建造一条超快、无摩擦的高速公路。在物理学世界中，这被称为超导性。通常情况下，电流会遇到阻碍（电阻）并损失能量转化为热量。而超导体就像是神奇的道路，电流可以在上面飞驰，而不损失哪怕一滴能量。

伟大的梦想是找到一种能在“室温”（就像一个温暖的夏日）下实现这种功能的材料，这样我们就能在各处使用它。然而，目前发现的最好的材料只有在承受巨大压力（就像被埋在行星深处那样）时才能发挥作用。这对于你的家庭或汽车来说并不实用。

这篇论文是一场关于寻找新型“神奇道路”的计算机化寻宝之旅。以下是他们是如何做的，以及他们的发现：

寻找“超硬”金属

研究人员研究了一类由三种常见元素组成的材料家族：硼 (B)、碳 (C) 和 氮 (N)。你可以把这些元素看作原子世界的“乐高积木”。

他们专注于两种特定的配方：B₂C₃N 和 B₄C₅N₃。

• 为什么选择这些？ 这些材料被预测是超硬的。想象一下一种坚硬到可以划伤几乎任何其他物质的材料，类似于钻石。
• 其中的联系： 通常，硬质材料的原子紧密锁定在一起，振动得非常快。研究人员怀疑，正因为这些材料如此坚硬且“紧张”，它们即使在没有巨型压力挤压的情况下，也可能成为优秀的无电阻导电体。

计算机模拟（“虚拟实验室”）

由于在真实实验室中制造这些材料非常困难，科学家们使用超级计算机充当了虚拟实验室。他们不仅仅是在猜测；他们使用了“第一性原理”计算。

• 类比： 想象你在预测一个复杂的舞池会如何运作。与其邀请真实的舞者，不如创建一个完美的数字模拟，模拟每一个舞者（原子）、他们如何牵手（化学键）以及他们如何摇摆（振动）。
• 他们模拟了电子（电流）如何在这些原子舞池中移动，以及它们如何与原子的振动（声子）进行相互作用。

重大发现：温暖天气的超导体？

结果令人兴奋。计算机模拟预测，这些超硬金属可以在比此类材料通常高的温度下实现超导：

• B₂C₃N 的超导温度可能约为 -233°C (40 K)。
• B₄C₅N₃ 的超导温度可能约为 -253°C (20 K)。

为什么这很重要？
为了让你有个直观的概念，目前常压超导体的冠军是一种名为 MgB₂（二硼化镁）的材料，它在 20 年前被发现，其工作温度约为 -234°C (40 K)。

• 新材料 B₂C₃N 被预测可以达到这位“冠军”的表现。
• 研究人员发现，这种材料的“硬度”在这里反而成了一种“超能力”。就像走钢丝的人需要紧绷且坚硬的绳索来保持平衡一样，这些超硬材料拥有坚硬的原子“绳索”，能够让电流平稳流动。

“各向异性”的转折

论文还发现了一个有趣的现象，即关于电流是如何流动的。

• 在某些材料中，电流在各个方向上的流动方式是相同的（就像水在圆形的管子里流动）。
• 在这些新材料中，流动过程更为复杂。研究人员必须使用高级数学（Eliashberg 方程）来弄清楚，电流的流动会根据其行进的方向而有所不同，就像足球在草地上滚动与在泥地上滚动的感觉不同一样。
• 他们发现，如果忽略这种复杂性，可能会低估这些材料的优越性。当他们进行正确的数学计算时，结果看起来非常有前景。

我们真的能制造出它吗？

论文谨慎地指出：“我们还没有制造它。”
然而，他们对原料进行了“成本检查”。他们计算了构建这些结构所需的能量，发现其能量水平与其他科学家已在实验室中成功制造的类似材料相当。

• 结论： 人类化学家极有可能利用现有的高科技方法（如高压炉或等离子体设备）创造出这些材料。

总结

研究人员利用超级计算机设计了一种由硼、碳和氮组成的新型“超硬”金属。他们预测，由于这些材料如此坚硬且紧绷，它们可以在约 -233°C 的温度下实现无电阻导电，达到了目前我们拥有的最佳材料水平。虽然他们尚未在真实实验室中制造出它，但数学表明这是可能的，这为寻找不需要在极端压力下被压碎的更好超导体提供了一条新路径。

技术摘要

技术摘要：超硬 B-C-N 金属中超导 $T_c$ 的第一性原理计算

问题陈述
对室温超导体的搜索一直由高压氢化物（如 $H_3S$、$LaH_{10}$）所驱动，但其对兆巴级压力的需求限制了实际应用。虽然轻元素系统（Be、B、C、N）提供了替代方案，但许多超硬材料（如金刚石）由于存在巨大的绝缘能隙而无法实现超导。相反，理论预测金属超硬化合物 $BC_5$ 能够表现出超导特性，但实验验证仍处于空白。本研究旨在通过调查三元硼-碳-氮（B-C-N）金属，寻找具有高转变温度（$T_c$）且可在常压下合成的超导体。具体而言，本文评估了预测的超硬结构 $B_2C_3N$ 和 $B_4C_5N_3$ 的超导性能，并将其与基准材料 $BC_5$ 以及常压纪录保持者 $MgB_2$ 进行对比。

方法论
作者利用密度泛函理论（DFT）和密度泛函微扰理论（DFPT）进行了第一性原理电子-声子计算。

• 软件： 计算使用了用于 DFT 的 Quantum ESPRESSO (v7.2) 和用于电子-声子耦合的 EPW 代码 (v5.7)。
• 结构： 本研究分析了 $BC_5$、$B_2C_3N$ 和 $B_4C_5N_3$（三方对称性，空间群 $P3m1$）以及 $MgB_2$（六方对称性，空间群 $P6/mmm$）。结构来源于先前的机器学习和进化搜索 [34] 以及 Materials Project [43]。
• 计算参数： 采用了规范赝势（PseudoDojo 库）和 PBE 交换相关泛函。通过 100 Ry 的波函数截断和密集的 k 点网格（例如 B-C-N 系统使用 $70 \times 70 \times 14$）实现收敛。
• 超导性分析： 研究通过 Wannier 插值求解了完整的 Eliashberg 能隙方程。应用了各向同性和各向异性 Migdal-Eliashberg (ME) 理论来确定超导转变温度（$T_c$）。库仑伪势（$\mu^*$）设定为 0.1，该值旨在重现 $MgB_2$ 的实验 $T_c$（约 40 K）。

关键结果

1. 电子结构： $BC_5$、$B_2C_3N$ 和 $B_4C_5N_3$ 表现出类金刚石的局部共价键，但实际上是空穴掺杂的，其费米能级（$E_F$）位于由 p 轨道主导的价带内。由于基于 B/N 价态比的不同空穴掺杂水平，$B_2C_3N$ 拥有比 $B_4C_5N_3$ 更高的费米能级态密度（$N_F$）。
2. 声子动力学： 所有材料在常压下均具有动力学稳定性。B-C-N 化合物表现出超硬材料特有的高频光学声子（$\sim 80$ meV），这显著贡献了电子-声子耦合（$\lambda$）。强耦合源于 B 原子与 C 原子之间的面内键伸缩运动。
3. 电子-声子耦合（$\lambda$）：
   • $BC_5$: $\lambda \approx 0.76$
   • $B_2C_3N$: $\lambda \approx 0.69$
   • $B_4C_5N_3$: $\lambda \approx 0.57$
   • $MgB_2$: $\lambda \approx 0.69$
     这些超硬材料的晶格刚性支撑了强耦合作用，即使在高声子能量下也能补偿某些情况下较低的 $N_F$。
4. 转变温度（$T_c$）：
   • 各向同性 vs. 各向异性： Allen-Dynes 方程（各向同性近似）低估了 $MgB_2$ 的 $T_c$（22.5 K vs. 实验值 40 K），凸显了各向异性计算的必要性。
   • 各向异性 ME 结果： 通过求解完整的各向异性 Eliashberg 方程得出：
     • $BC_5$: $T_c \approx 41.2$ K
     • $B_2C_3N$: $T_c \approx 38.2$ K
     • $B_4C_5N_3$: $T_c \approx 22.9$ K
     • $MgB_2$: $T_c \approx 40.7$ K
   • 各向异性： $BC_5$ 和 $B_2C_3N$ 表现为具有均匀能隙的各向同性超导体。相比之下，$B_4C_5N_3$ 和 $MgB_2$ 表现出强烈的各向异性以及多能隙结构（双能隙结构），如果使用各向同性近似，会导致 $T_c$ 被显著低估。

意义与主张
论文声称，超硬金属 B-C-N 化合物，特别是 $B_2C_3N$ 和 $BC_5$，是高 $T_c$ 常压超导体的有力候选材料。相对较高的 $T_c$ 值（高达 $\sim 40$ K）归功于高德拜温度（与超硬特性相关）以及由高频光学声子介导的强电子-声子耦合。

至关重要的是，作者指出 $B_2C_3N$ 和 $B_4C_5N_3$ 的理论形成能与近期合成的三元超硬材料（如 $BC_2N$、$BC_{10}N$）相当，这表明这些化合物有望通过高压高温（HPHT）或微波等离子体化学气相沉积（MPCVD）技术进行合成。研究结论认为，研究超硬金属是发现新型常压声子介导超导体的可行途径，前提是必须采用准确的各向异性 Eliashberg 计算来确定 $T_c$。