固态点缺陷的分子视角

《ACS Central Science》：Solid-State Point Defects through a Molecular Lens

　　

固体中可光学寻址的点缺陷，通常称为色心（color centers），已成为量子传感（quantum sensing）、量子通信（quantum communication）和量子计算（quantum computing）中有前景的组分。金刚石中带负电的氮空位中心（negatively charged nitrogen-vacancy center, NCVC–）最为著名，但许多其他点缺陷可被实现并可能为特定应用定制。缺陷的表征与理性设计可借助高精度计算机模拟，但许多缺陷的电子结构对传统方法具有挑战性，推动了理论化学与物理学界的大量研究。缺陷计算困难在于：价电子可占据若干近简并、局域化的缺陷轨道，需要多组态电子结构方法，如完整活性空间方法；而宿主晶格原子数巨大，需用大超胞与周期边界条件以消除有限尺寸误差。完整活性空间自洽场（complete active space self-consistent field, CASSCF）加微扰理论等无法扩展到合适超胞的原子数，也很少支持周期边界条件。研究人员在《ACS Central Science》提出嵌入型关联波函数（embedded correlated-wavefunction, CW）框架：碳空位、邻近替代原子（硼、氮或氧）及若干近邻碳原子被表示为“空位中心分子（vacancy-centered molecules）”，嵌入由量子力学导出、近似周围金刚石晶格作用的一电子势中。随后采用多组态CW方法，特别是CASSCF加n电子价态微扰理论（n-electron valence state perturbation theory, NEVPT2），计算基态与激发态性质及各自旋态相对能量，识别出BCVC–与OCVC0可作为NCVC–的候选替代。框架先用周期密度泛函理论（periodic density-functional theory, DFT）给出缺陷几何与电子环境，提取含掺杂剂、空位周围低配位碳原子及近邻原子的局部簇并化学封端，将所得“分子”置于优化嵌入势中，使该分子DFT计算复现周期DFT结果，再将近似DFT升级为CASSCF+NEVPT2高精度计算。这种分工适合固态点缺陷：周期DFT提供几何与晶格电子影响，嵌入CW捕捉局域缺陷态的多组态特征，可计算自旋态热力学、光学跃迁能量与偶极矩及辐射寿命并与实验比较。嵌入分子描述提供类配体场理论的化学直觉：前沿轨道主要来自空位周围低配位碳，掺杂剂调节其能量、占据与性质，对称性解释优选自旋态与Jahn-Teller畸变，并给出“空位提供轨道框架、掺杂剂调节电子结构”的设计原则。已建立的NCVC–是基准，其计算自旋态与光学性质与实验吻合；OCVC0是NCVC–最近电子类似物，但亮跃迁由绿移蓝；BCVC–的更亮类NV跃迁不是最低三重态激发而是第二重，位于紫外；两者均预测比NCVC–辐射寿命更短且存在可能参与 shelving 路径的近单重态。要将电子结构联系到实验光循环，还需垂直激发能之外的结构与动力学信息：零声子线、振子光谱、非辐射衰减率依赖缺陷与晶格弛豫，需要CW层面对基态与激发态几何优化及晶格响应一致处理；自旋轨道与电子声子耦合可将态能量转为辐射与非辐射率，原子结构与电荷态稳定性决定光循环能否实现。该嵌入CW框架虽使高精度计算可行，但可能对高通量材料发现过贵，更适合作为廉价方法基准或分级筛选最后一步；其对缺陷态较离域、与宿主能带强杂化体系（如硅中某些缺陷）的适用性待考察。硼、氧空位中心能否用于量子信息由实验决定，但嵌入分子视角提供了甄别候选者的有力方式。

研究背景方面，固体中可光学寻址的点缺陷被称为色心（color centers，固态中局域电子态导致光学活性的点缺陷），有望用于量子传感、量子通信与量子计算。金刚石中带负电的氮空位中心（N_CV_C^–）最知名，但其他点缺陷可被实现并定制。缺陷表征与理性设计依赖高精度模拟，然而许多缺陷电子结构对传统方法困难：价电子占据近简并、局域化缺陷轨道，需要多组态电子结构描述；宿主晶格原子数巨大，需用大超胞与周期边界条件消除有限尺寸误差，而完整活性空间自洽场（CASSCF）加微扰理论难以扩展到合适超胞原子数，也少有周期边界条件支持。因此，研究人员需要一种兼顾晶格环境与局域多组态电子相关的方法。

论文发表在《ACS Central Science》，研究人员提出嵌入型关联波函数（embedded correlated-wavefunction, CW）框架，把金刚石色心当作嵌入晶格势中的空位中心分子来处理，并用多组态CW方法计算基态、激发态、自旋态与光学性质，识别除N_CV_C^–外的可能替代色心。

主要关键技术方法方面，研究人员首先用周期密度泛函理论（periodic DFT）获得弛豫缺陷几何与晶格电子环境；随后提取以空位为中心、含掺杂剂与近邻低配位碳原子的局部簇并进行化学封端，构造可复现周期DFT结果的优化嵌入一电子势；最后在该嵌入模型中执行CASSCF与n电子价态微扰理论（NEVPT2）多组态关联波函数计算。该方法无生物样本队列，属于计算材料与量子缺陷电子结构研究。

Why are defect calculations challenging?

研究人员指出，缺陷计算难点来自两方面：局域缺陷轨道近简并，需要完整活性空间等多组态方法；宿主晶格需大超胞与周期边界条件，而高精度多组态方法难以扩展且少支持周期性。因此传统单一层次方法不能同时兼顾晶格尺度与缺陷态电子相关。

An embedded correlated-wavefunction framework

研究人员以周期DFT提供几何与电子环境，将空位中心区域提取、化学封端后放入优化嵌入势，再进行CASSCF+NEVPT2计算。这样周期DFT负责晶格环境，嵌入CW负责局域缺陷多组态相关，可计算自旋态能量、光学跃迁能量、偶极矩与辐射寿命，并直接与实验比较。

A chemically intuitive picture

研究人员将空位周围低配位碳视为前沿轨道主要来源，掺杂剂调节轨道能量、占据与特征；对称性解释优选自旋态与Jahn-Teller畸变。由此形成化学设计原则：空位提供轨道框架，掺杂剂调节电子结构。该图像类似化学家熟悉的配体场理论。

Benchmarks and candidate color centers

以N_CV_C^–为基准，其计算自旋态与光学性质和实验吻合，证明框架可靠。O_CV_C⁰是N_CV_C^–最近电子类似物，亮跃迁由绿移蓝；B_CV_C^–更亮的类NV跃迁不是最低三重态激发而是第二重，位于紫外。两者均预测比N_CV_C^–辐射寿命更短，并存在可能参与shelving路径的近单重态。

Beyond vertical excitations

研究人员强调，零声子线、振子光谱、非辐射衰减率依赖激发后缺陷与晶格重组，需要CW层面对基态与激发态几何优化及晶格响应一致处理；自旋轨道耦合与电子声子耦合是把态能量转为辐射与非辐射速率所必需的，原子结构与电荷态稳定性决定预测光循环能否实验实现。

Cost and scope of the framework

研究人员认为嵌入CW框架精度高但计算贵，不适合高通量材料发现常规使用，更适合作为廉价方法基准，或用于分级筛选流程中、先由低价方法缩小候选空间后的最终高精度步骤；当缺陷态离域并与宿主能带强杂化时，嵌入簇描述表现仍需检验。

讨论与结论部分总结及结论翻译：研究人员认为该框架为金刚石色心建立了准确电子结构基础，但要把电子结构连接到实验相关光循环，还需超越垂直激发能的结构与动力学信息。嵌入CW计算虽可行，但可能对高通量发现过贵，宜作基准或分级筛选末端高精度环节；其对缺陷态离域、与宿主能带强杂化体系的外推性有待验证。硼空位与氧空位中心能否最终用于量子信息应用仍有待实验判断，但嵌入分子视角提供了一种有说服力的方式，用于识别值得进一步理论与实验研究的候选体系。

结论原文意译：尽管嵌入CW框架已使高精度计算可行，它可能对高通量材料发现而言过于昂贵；这类计算可作为更廉价方法有价值的基准，或作为分层筛选流程中、较廉价计算缩小候选空间后的最终高精度步骤。其超出当前体系的适用性，取决于当缺陷态较不离域且与宿主能带杂化更强时嵌入簇描述的表现，例如硅等较小带隙半导体中的某些缺陷。硼、氧空位中心是否最终在量子信息应用中实用，仍需实验判定；但嵌入分子视角为识别值得进一步理论与实验研究的候选者提供了引人注目的途径。