An extended ab initio theory of the V$_{\text{B}}^-$ center in hBN: excited states, Jahn-Teller distortion, and pressure dependence

本文采用高级 CASSCF-NEVPT2 计算，全面模拟了六方氮化硼中 V$_{\text{B}}^-$ 缺陷中心的激发态、结构畸变及应力依赖性质，从而阐明了其复杂的光学循环，并为先进的二维量子传感应用奠定了理论基础。

要点

想象一下，在六方氮化硼（hBN） sheets 内部隐藏着一个微小、不可见的“量子灯泡”，六方氮化硼本质上是一种超薄的原子级平整材料层。这个灯泡是一种特定的缺陷，称为带负电的硼空位（$V^-_B$）。科学家们对它感到兴奋，因为它可以作为磁场和其他微小力的传感器，在室温下工作，甚至能集成到超薄的二维器件中。

然而，长期以来，科学家们并未完全理解这个“灯泡”是如何工作的。他们知道它会发光并对磁场产生反应，但其内部机制一直是个谜，因为涉及的电子是“强关联”的——这是一种委婉的说法，意指它们以复杂、混沌的方式共同“舞蹈”，而标准的计算机模型难以轻易预测。

本文就像一本高分辨率手册，利用先进的计算机模拟终于阐明了这个量子灯泡的内部运作机制。以下是他们研究发现的分解，辅以简单的类比：

1. 变形之舞（Jahn-Teller 畸变）

当灯泡被激光（如绿光）激发时，它不会静止不动。想象一个由三个氮原子组成的完美圆形等边三角形。当电子被激发时，这个三角形会突然在一个方向上被“拉伸”，变成一个不对称的形状。

• 论文主张： 这种拉伸被称为Jahn-Teller 畸变。它不是微小的颤动，而是重大的结构变化。三角形的畸变如此剧烈，以至于在能量景观中形成了一个“三尖帽”形状（想象一顶有三个明显山谷的帽子）。
• 后果： 在低温下（低于 200 K），三角形会“卡”在这三个山谷中的某一个里（静态状态）。但在室温下，它具有足够的能量在三个山谷之间快速跳跃（动态状态）。这种跳跃改变了灯泡的行为方式，以及它如何分裂其磁信号。

2. 缺失原子的“幽灵”

该缺陷的产生是因为缺失了一个硼原子。这使得邻近的氮原子上留下了六个“悬空”的电子轨道。

• 论文主张： 作者绘制了这些电子的能级图。他们发现，灯泡吸收绿光（约 2.3 eV）以被激发。然而，当它弛豫回基态时，并不会发出单一锐利的颜色。相反，它会发出一种宽阔、模糊的光（“声子边带”），因为形状变化如此剧烈，以至于每发射一个光子，就会踢出大约五个“声波”（声子）。
• 结果： 光的“纯净”颜色（零声子线）非常微弱（仅占总光量的 0.4%），几乎不可见，被宽阔、模糊的光晕所掩埋。这解释了为什么实验难以观察到锐利的颜色峰值。

3. 秘密隧道（系间窜越）

这个灯泡用于传感的魔力在于它能够切换不同的“自旋”状态（将其想象为微小内部指南针的不同取向）。

• 论文主张： 作者发现，电子切换自旋所采取的路径高度依赖于其取向（$m_S = 0$ 与 $m_S = \pm 1$）。
  • 一条路径快速且直接。
  • 另一条路径涉及一个“准简并”态，其中单重态（一种自旋类型）和三重态（另一种自旋类型）在能量上如此接近，以至于几乎接触。
• 类比： 想象两条平行的铁轨，它们靠得如此之近，以至于如果铁轨震动，火车就能轻易地在它们之间跳跃。这种“跳跃”（系间窜越）使得设备能够被光学读取。论文表明，这种跳跃对温度和压力高度敏感。

4. 挤压灯泡（压力与应变）

研究人员还测试了如果挤压材料（施加压力）会发生什么。

• 论文主张：
  • 从顶部挤压（垂直压力）： 这使得材料的层与层之间靠得更近。它显著加快了“自旋跳跃”过程，使灯泡变暗，寿命缩短。
  • 从侧面挤压（水平压力）： 这改变了基态的磁“分裂”（D 参数）。
• 结论： 灯泡是一个非常灵敏的应变仪。它对压力的反应取决于你从哪个方向挤压它。论文证实，压力下磁信号的变化是由于原子晶格的物理压缩造成的。

5. 论文未提及的内容

重要的是要注意这篇论文没有声称：

• 它没有声称已经制造出工作的商业传感器。
• 它没有声称解决了所有谜团。作者承认，从“零自旋”态到单重态的跃迁对于他们目前的模型来说仍然过于复杂，无法完美计算。他们建议未来的工作需要更先进的模拟方法来完全理解这一特定的“跳跃”。
• 它没有讨论临床用途或医疗应用。

总结

简而言之，这篇论文利用超先进的计算机建模，绘制了$V^-_B$ 中心的详细地图。它解释说，这种量子缺陷是一个变形者，在受激发时会扭曲其自身的原子结构，从而创造出复杂的能量景观。这种畸变决定了它如何发光、如何切换磁自旋，以及如何对受挤压做出反应。这一理论地图为将该缺陷转化为可靠的纳米级量子传感工具奠定了基础。

技术摘要

技术摘要：hBN 中 V⁻B 中心的扩展从头算理论

问题陈述
六方氮化硼（hBN）中的带负电硼空位（V⁻B）中心已成为极具前景的二维自旋量子比特，适用于纳米级量子传感，其在表面邻近性和空间分辨率方面优于体相金刚石系统。然而，对其光探测磁共振（ODMR）信号的全面理论描述仍难以企及。这一挑战源于光学循环中涉及的电子态具有强关联和多参考特性。先前的理论工作难以完全捕捉激发态精细结构、Jahn-Teller（JT）畸变以及系间窜越（ISC）路径，从而限制了优化 V⁻B 中心以用于集成量子传感器的能力。

方法论
为应对这些局限，作者采用了一种基于波函数的高级别ab initio方法，具体为完全活性空间自洽场（CASSCF）辅以n电子价态微扰理论（NEVPT2）。选择该方法是因为其能够准确描述静态和动态电子关联效应，而这些效应对于 V⁻B 中心的多参考特性至关重要。

• 模型系统：计算在代表 V⁻B 缺陷的 B₆₀N₆₀H₂₇ 团簇模型上进行。
• 电子结构：活性空间包含局域于第一近邻氮原子上的六个缺陷轨道（三个面内$a'_1$和$e'$态，以及三个面外$a''_2$和$e''$态）。
• 动力学与速率：结构弛豫在 CASSCF 水平上计算。光致发光（PL）和系间窜越（ISC）的跃迁速率利用一维有效声子模型进行近似。零场分裂（ZFS）参数通过准简并微扰理论（QDPT）推导得出。
• 外部微扰：通过变形分子团簇几何结构（调整层间距离和面内晶格常数）并重新评估电子性质，来模拟静水压和单向应变的影响。

主要贡献与结果

1. 电子结构与能级排序：
   本研究确立了 V⁻B 中心的能谱，确定基态为$^3\bar{A}'_2$（$D_{3h}$对称性）。第一个光允许跃迁发生在$^3\bar{E}'$态。激发后，系统迅速发生内转换至较低能量的三重态。关键在于，$^3\bar{E}''$态经历显著的伪 Jahn-Teller（pJT）畸变，将其对称性降低至$C_{2v}$，并分裂为较低能量的$^3\bar{A}_2$态和较高能量的$^3\bar{B}_2$态。$^3\bar{A}_2$态被确定为能量最低的三重态激发态。

2. Jahn-Teller 畸变与势能面：
   $^3\bar{E}'' \to ^3\bar{A}_2$弛豫由强烈的 pJT 效应驱动，导致 355 meV 的稳定化能，并产生结构畸变，其中内部氮三角形发生拉伸（N-N 距离：2.49 Å 和 2.65 Å）。势能面被描述为具有三个极小值的“三叉帽”形状，极小值之间由约$\delta_{JT} \approx 308$ meV 的势垒分隔。

   • 温度依赖性：在约 200 K 以下，系统在$C_{2v}$构型下表现为静态 JT 系统。在 200 K 以上，在激发态寿命内，三个 JT 极小值之间发生跃迁，导致动态 JT 机制。这种动态行为解释了横向零场分裂（ZFS）参数（$E$）的温度依赖性。

3. 光学性质与零声子线（ZPL）：
   主导的辐射跃迁路径是$^3\bar{A}_2 \to ^3\bar{A}'_2$。计算得出的 ZPL 能量为 1.88 eV。然而，计算得到的 Huang-Rhys 因子为$S = 5.4$，表明存在显著的电子 - 声子耦合和大的声子边带（PSB）。德拜 - 沃勒因子极低（$f_{DW} \approx 0.4\%$），表明 ZPL 被 PSB 掩盖，这与实验观察到的宽发射光谱一致。计算得出的 PSB 峰值（1.60 eV）与实验数据吻合良好。

4. 系间窜越（ISC）与自旋极化：
   本研究详细阐述了$m_S = 0$和$m_S = \pm 1$自旋子能级的不同 ISC 路径，这与金刚石中的 NV 中心形成对比。

   • $m_S = \pm 1$态通过$^3\bar{A}_2 \to ^1\bar{E}'$通道快速衰变至单重态区域（$\tau \approx 3.2$ ns）。
   • $m_S = 0$通道涉及准简并的单重态 - 三重态（$^3\bar{A}_2$和$^1\bar{A}_2$），使得跃迁速率对温度和应变高度敏感。
   • 系统布居到一个寿命约为 340 ns 的亚稳态单重态 shelving 态（$^1\bar{E}'$），随后衰变回基态，从而促进自旋初始化和读出。

5. 应变与压力依赖性：
   作者模拟了 ODMR 参数对静水压和单向应变的响应。

   • 基态：基态（$^3\bar{A}'_2$）的$D$张量随压力呈现正移，这主要由面内（水平）压缩而非层间（垂直）压缩驱动。计算得出的斜率（~38 MHz/GPa）与实验值定性一致。
   • 激发态动力学：$m_S = \pm 1$通道的 ISC 速率对压力高度敏感，在 2.7 GPa 静水压下增加约两倍。这归因于自旋轨道耦合矩阵元对垂直压缩下能隙减小的敏感性。这一发现从理论上解释了实验观察到的压力下光致发光强度降低和激发态寿命缩短的现象。

意义与主张
本文声称提供了首个针对 V⁻B 中心的扩展、高级别理论框架，成功调和了复杂的多参考电子效应与实验 ODMR 观测结果。通过阐明伪 Jahn-Teller 效应的作用、激发态从静态到动态的转变以及独特的自旋依赖性衰变路径，该工作为理解 V⁻B 中心的基本行为奠定了理论基础。

作者谦逊地总结道，虽然其结果与实验测得的 ZFS 和光致发光参数吻合良好，但对非辐射衰变（特别是$m_S=0$通道）的完整定量描述，需要进一步发展超越玻恩 - 奥本海默近似的多表面动力学模拟方法。尽管如此，该研究为部署 V⁻B 中心作为集成二维量子传感器提供了关键的理论支持，特别是在解释应变和压力传感能力方面。