First-principles calculations of internal conversion processes in spin defects

本文引入了一种预测性的第一性原理框架，该框架结合了多构型 TDDFT 与解析非绝热耦合，用以精确计算光学活性自旋缺陷中的内转换速率，并成功解决了钻石 NV$^-$ 色心与 SiC 双空位与实验数据之间长期存在的差异。

要点

想象一下，在晶体（如钻石或碳化硅）内部有一个微小的、发光的缺陷，它就像一个微型量子计算机。这些缺陷就像微小的舞台演员。当我们用激光照射它们时，它们会被激发并跳到一个更高的能量层（即“舞台”）。为了回到它们的静止状态，它们必须选择一条路径：要么明亮地发光（辐射衰减），要么悄无声息地滑落回原位（非辐射衰减）。

长期以来，试图预测这些演员“滑落回原位”速度（一个被称为内转换的过程）的科学家们一直使用的是一张非常粗略的地图。他们的计算就像是通过只观察单行道上的一辆车来预测交通流量。他们一直在猜测速度极其缓慢，但实际上，交通流量移动得很快。他们的预测偏差巨大——有时甚至慢了上千倍。

这篇论文引入了一套全新的、高分辨率的 GPS 系统来修正这些预测。以下是作者如何实现这一目标的，使用了简单的类比：

1. “多体”问题：看见整个管弦乐队

以往的方法将缺陷中的电子视为演奏单一音符的独奏音乐家。但实际上，这些电子是一个复杂的爵士乐队，它们同时在进行即兴演奏并相互反应。

• 旧方法： 忽略乐队之间的相互作用，将电子视为仅仅是一个人。
• 新方法： 作者使用了一种先进的方法（结合了杂化泛函的 TDDFT）来聆听整个管弦乐队。通过考虑所有电子是如何共同起舞的（多构型效应），他们终于能够听到能量层真实的复杂性。

2. “振动”问题：计算每一步

当一个电子下降能量层时，它并不只是落下；它必须将其多余的能量传递给晶体的原子，使它们发生振动。你可以把晶体想象成一个由数百万个弹簧组成的巨大蹦床。

• 旧方法： 科学家们过去常常假装这个蹦床只有一个弹簧，或者只有几个“主要”弹簧，以节省时间。他们根据仅有的那几个弹簧来计算能量释放。
• 新方法： 作者意识到，蹦床上的每一个弹簧都对这次落下有所贡献，而不仅仅是靠近缺陷的那几个。他们开发了一种方法，可以同时计算与所有振动原子的相互作用，而不只是那些靠近缺陷的原子。他们通过解析地计算“非绝热耦合”（一种衡量电子如何推动原子的复杂方式），这就像是拥有了一个关于“推力”的数学公式，而不是通过试错法来猜测。

结果：修复地图

作者用这两个著名的“演员”测试了他们的新 GPS：

1. 钻石演员（NV- 中心）：

   • 谜团： 科学家知道这个演员在特定的激发态下寿命非常短，但旧的计算认为它应该存活得久得多。
   • 修复方案： 新方法计算了“滑落”速度，发现它极其迅速（大约每秒 1000 亿次）。这与最近的超快实验测量完美吻合。它证实了这种“滑落”正是该演员无法长时间保持激发状态的主要原因。

2. 碳化硅演员（空位对中心）：

   • 谜团： 对于这个演员，旧的计算认为它应该激发约 37 纳秒（基于仅发光的假设）。但实验显示它仅持续约 15 纳秒。有些东西丢失了。
   • 修复方案： 新方法发现了一个科学家之前忽略的“隐藏门”。他们发现了一个显著的、此前被忽视的“滑落”路径（非辐射通道），该路径加速了衰减。当他们把这条隐藏路径加入到数学模型中时，预测结果终于与实验（15 纳秒）相匹配。

为什么这很重要

这篇论文不仅解决了一个数学问题，它还提供了一个通用工具包。

• 它证明了忽略“整个管弦乐队”（电子相互作用）或“所有的弹簧”（振动）会导致极其错误的答案。
• 它使科学家能够精确预测这些量子缺陷的行为，而无需先进行猜测或进行昂贵的实验。
• 它为设计更好的量子计算机奠定了基础，因为我们可以准确了解这些微小的“量子比特”（缺陷的磁态）在失去能量前能维持多久。

简而言之，作者制造了一台显微镜，既能看到电子复杂的舞蹈，也能看到每一个原子的振动，终于让我们能够准确预测这些量子缺陷何时“关闭”。

技术摘要

技术摘要：自旋缺陷中内转换过程的第一性原理计算

问题陈述
光学活性自旋缺陷，例如金刚石中的负电荷氮空位（$NV^-$）中心和碳化硅（SiC）中的中性双空位（$VV^0$）中心，是量子技术的关键平台。虽然这些系统的辐射跃迁和系际穿梭（ISC）速率已在第一性原理框架下得到了广泛研究，但内转换（IC）过程——即由电子-声子耦合驱动的相同自旋电子态之间的非辐射跃迁——仍缺乏准确描述。常见的近似方法（如基于 Kohn-Sham 轨道的单模或接受模方法）在历史上往往会低估 IC 速率几个数量级。这种差异阻碍了对光学循环的完整第一性原理描述，限制了预测激发态寿命、温度效应以及光检测磁共振（ODMR）光谱的能力。

方法论
作者提出了一个广泛且具有预测性的第一性原理框架，用于计算 IC 过程的非辐射跃迁速率（NRTR）。该方法集成了两项关键进展：

1. 多体电子结构： 该框架利用线性响应时间依赖密度泛函理论（LR-TDDFT）结合混合泛函（特别采用了 DDH 泛函）。这种方法能够捕捉激发态的多构型特性，而这在单轨道 Kohn-Sham 描述中往往是缺失的。作者采用自旋翻转 TDDFT 框架来处理特定的激发态。
2. 解析非绝热耦合（NACs）与全声子模式： 作者没有依赖有限差分法计算波函数导数，而是基于 WEST 代码中扩展拉格朗日方法实现了一种解析形式的 NAC 矢量（$d_{IF,A}$）。这使得高效计算大型超胞（包含数百个原子）中所有声子模式的电子-声子矩阵元（$M^k_{IF}$）成为可能，从而避免了以往方法的瓶颈。

NRTR，$\Gamma_{NR,IF}(T)$，通过费米黄金定则进行计算，并对所有 $N_{modes}$ 个声子模式的贡献进行求和。计算采用生成函数法来评估声子弛豫贡献（$X^k_{IF}$）和电子-声子谱密度函数 $D(E)$，并考虑了多声子过程。

关键结果
该框架被应用于金刚石 $NV^-$ 中心和 4H-SiC $VV^0$ 中心的光学循环：

• $NV^-$ 中心（金刚石）：

  • 单线态（$^1A_1 \to ^1E$）： 计算出的 NRTR 数量级约为 $\approx 10^1$ GHz。结果与飞秒瞬态吸收实验表现出极佳的定量一致性（例如，在 78 K 时 $\tau_{^1A_1} \approx 117.4$ ps，对比实验值约为 $\sim 92-102$ ps）。这证实了 IC 是该高能单线态的主要衰减机制。
  • 三线态（$^3E \to ^3A_2$）： 计算出的 NRTR 极小（$\lesssim 5$ MHz），相对于辐射速率（$\approx 83.3$ MHz）而言可以忽略不计，这与该跃迁主要为辐射跃迁的理解一致。
  • 模式分析： 研究表明，尽管一个位于 $\sim 170$ meV 的局域振动模式具有很强的非绝热特性，但它仅占总速率的大约一半。包含所有声子模式对于保证准确性至关重要；单模近似无法捕捉到完整的谱密度。

• $VV^0$ 中心（4H-SiC）：

  • 三线态（$^3E \to ^3A_2$）： 此前仅考虑辐射贡献的第一性原理估算高估了实验寿命（$\sim 15$ ns）。通过引入 IC 过程，作者解决了这一长期存在的差异，表明非辐射衰减显著降低了寿命，使其与实验观察结果相匹配。
  • 单线态（$^1A_1 \to ^1E$）： 在缺乏实验数据的情况下，该框架预测其 NRTR 约为 $\approx 10^1$ GHz，表明该跃迁同样由 IC 过程主导。

意义与主张
本文声称提供了首个能够计算 IC 速率并与实验实现定量一致的高效且准确的 ab-initio 框架。作者强调了两个主要贡献：

1. 解决差异： 该框架通过正确捕捉电子态的多体性质以及所有声子模式的贡献，解决了此前理论估算与实验寿命之间存在的数量级差异。
2. 方法论进步： 通过实现解析 NACs，这项工作消除了有限差分导数的计算瓶颈，从而能够在包含数百个原子的系统中纳入所有声子模式。

作者得出结论，有效的单模和接受模近似会系统性地低估非辐射速率几个数量级。他们的方法能够实现对金刚石和 SiC 中自旋缺陷完整光学循环的无参数预测。此外，解析 NAC 的实现被视为向扩展体系中 ab initio 非绝热分子动力学以及极化子和激子动力学研究迈进的重要一步。