键合与反键合共振：单线态双自由基特征的一个统一框架

《Chemistry – A European Journal》：Bond and Antibond Resonances: A Unified Framework for Singlet Biradical Character

研究人员发现，目前对单线态双自由基(oid)的结构与成键，尤其是其“双自由基特征”程度的化学概念理解尚不充分，且文献中使用的指标各异，难以直接比较，阻碍了对不同双自由基(oid)物种的系统性理解和比较。为解决这一问题，研究者采取整体性方法，考虑了多种双自由基指标（如LUNO占据数、单线态-三重态（ST）能隙、自由基中心间键级等），并应用多种计算化学方法（DFT、CCSD(T)、CASSCF、NEVPT2、MRCI），旨在获得一个涵盖从闭壳层分子到“完美”开壳层双自由基整个双自由基特征范围的统一描述框架。他们提出，对于中低程度双自由基特征的物种，可以利用化学直观的“键”和“反键”来描述自由基电子间的相互作用。这项研究为解决长期存在的关于臭氧O3和二氮化二硫S2N2等分子双自由基特征的争议提供了新见解，并为理解具有双自由基特征的物种提供了一种更统一的、直观的化学描述框架。

　　

在化学的奇妙世界里，有一类分子因其独特的电子结构和反应性而备受关注，它们被称为双自由基（Biradical）或双自由基体（Biradicaloid）。与所有电子都成对存在的常规闭壳层分子不同，这些物种的特征在于两个电子之间几乎没有或只有微弱的相互作用，这导致了开壳层特征，并赋予它们高反应性以及非线性光学、小分子无势垒活化等不寻常的性质。这些分子在光化学过程、分子开关和自旋信息存储技术中扮演着关键角色。然而，尽管如今对双自由基电子结构的数学描述已较为成熟，但科学家们感到，对于其结构、成键，尤其是单线态双自由基体的化学概念理解仍然不够深入。更令人困惑的是，文献中讨论单线态双自由基的“双自由基特征”时，使用了各种各样的指标，这些指标往往无法直接比较，从而阻碍了建立系统性理解以及比较不同双自由基物种的进程。这就好比用不同的尺子去测量同一个物体的长度，结果自然难以统一和比较。为了改善这种状况，一篇发表在《Chemistry – A European Journal》上的研究论文，题为“键与反键共振：单线态双自由基特征的一个统一框架”，旨在为理解这类迷人的分子提供一个更清晰、更直观的蓝图。

为了开展这项研究，研究人员采取了一种整体性的方法。他们精心挑选了一系列模型体系，主要是无机四元环体系，这些体系在形式上与S₂N₂是等瓣相似的。研究涵盖了从文献已知的到假想的共计28个分子，此外还包括了氧气、臭氧、烯丙基阴子等分子，以展示其结论的广泛适用性。研究团队运用了一套强大的计算化学“工具箱”，包括密度泛函理论（DFT）、耦合簇理论（CCSD(T)）、完全活性空间自洽场方法（CASSCF）、N电子价态微扰理论（NEVPT2）以及多参考组态相互作用方法（MRCI）。他们计算了多种与双自由基特征相关的关键指标，包括垂直和绝热单线态-三重态能隙（ST gap，ΔE_ST）、从CASSCF(2,2)波函数导出的双自由基特征β、基于NEVPT2和MRCI自然轨道占据数计算得到的最低未占自然轨道（LUNO）占据数、基于定域轨道和CASSCF密度矩阵计算的自由基中心间的键级（Bond Order, BO）以及重要路易斯共振结构的权重。通过比较这些指标，并结合对双电子-两轨道（2e2o）模型的深入理论分析，研究人员旨在建立一个统一且直观的描述框架，以跨越从闭壳层分子到“完美”开壳层双自由基的整个光谱。

3.1 结构优化

为了获得可靠的分子结构用于高精度计算，研究人员首先对模型体系进行了几何优化。他们主要使用了密度泛函理论（DFT）中的PBE泛函，并将其结果与更高精度的CCSD(T)方法优化的结构进行了基准对比。研究发现，DFT与CCSD(T)优化的结构总体上吻合得很好，平均键长误差仅为0.014 ?，表明PBE这类纯泛函对于研究具有不同程度双自由基特征的分子结构是可靠的。研究还注意到，在处理单线态双自由基体时，需要关注DFT方法可能出现的对称性破缺问题，PBE泛函在大约50%双自由基特征时会发生这种情况。

3.2 单线态-三重态能隙

单线态-三重态能隙是描述双自由基中两个未成对电子相互作用强度的关键参数。研究人员使用PBE-D3、CCSD(T)、NEVPT2和MRCI+Q多种方法计算了垂直ST能隙，并以MRCI结果为基准进行了比较。结果显示，CCSD(T)与MRCI结果最为接近，平均误差很小，而PBE-D3的结果也相当准确，尽管误差分布稍大。这表明即使在处理具有较大取代基的实际分子时，PBE-D3等纯泛函也是一个可行的选择。研究还对绝热ST能隙进行了分析，发现其值可能因三重态与单线态的几何结构差异而与垂直能隙有很大不同，例如P₂C₂^Cl的绝热能隙非常小，这与其自由基中心的几何形变有关。

3.3 双自由基特征：理论考量

ST能隙能否直接衡量双自由基特征？研究人员通过严格的数学推导给出了否定的答案。在2e2o模型中，垂直ST能隙ΔE_ST^v不仅取决于与双自由基特征β相关的参数（ΔE_ab/K_ab，其中ΔE_ab可视为HOMO-LUMO能隙，K_ab是交换积分），还额外依赖于一组库仑积分（J_aa, J_bb, J_ab）。因此，ST能隙与双自由基特征之间没有简单的数学对应关系。对于不交盖双自由基（如拉伸的H₂），其前线轨道是相同原子轨道的成键与反键组合，库仑项C = (J_aa+J_bb)/2 – J_ab通常较小；而对于交盖双自由基（如O₂），其前线轨道空间分布不同，库仑项C较大，这解释了为什么交盖双自由基通常有较大的正ST能隙。此外，交换积分K_ab的大小调节了β对ΔE_ST^v的影响强度。研究通过参数化绘图清晰地表明，ΔE_ST^v与β之间的关系是模糊的，一个量不能简单地由另一个推导出来。

3.4 双自由基特征：计算

为了将理论与高精度计算对比，研究人员计算了所有模型分子的单线态双自由基特征。他们评估了不同理论级别（CASSCF(2,2)、NEVPT2、MRCI）下的LUNO占据数。结果显示，对于大多数分子，包含动态相关性的NEVPT2和MRCI计算得到的LUNO占据数略高于仅考虑静态相关的CASSCF结果。研究列出了所有分子的具体数值，例如S₂N₂的MRCI LUNO占据数为0.14，O₃为0.31，而F₂Ge₂²⁺则高达0.77，接近“完美”双自由基。这些数据与ST能隙数据相结合，为后续的关联分析和化学解释提供了坚实基础。

3.5 统一描述：键、反键与双自由基特征

这是本研究的核心发现。研究人员提出，对于具有中低程度双自由基特征的物种，其电子结构可以通过路易斯图像中的“单键”和“单反键”作为核心成键图案来直观描述。这里的“反键”是一个化学概念，用于描述两个定域在相近空间区域但自旋反平行的电子之间的弱排斥相互作用。研究通过定域分子轨道分析和计算键级证实了这一点。他们发现，随着双自由基特征的增强，连接两个自由基中心的“键”逐渐减弱，而“反键”特征增强。在完美的开壳层双自由基极限下，两个自由基中心之间既没有键也没有反键。这一框架成功地将复杂的多参考波函数信息转化为了化学家熟悉的键和反键语言，从而可以像理解经典路易斯结构一样理解这些物种。应用这一框架，研究人员澄清了关于O₃和S₂N₂等分子双自由基特征的长期争论。例如，臭氧（O₃）可以描述为具有一个弱的O-O σ键和一个强的O-O σ反键，这与其适中的双自由基特征（β ≈ 0.3）和负的ST能隙相符。

3.6 共振结构与双自由基描述符

研究进一步将定域轨道波函数转换为价键（VB）型波函数，从而得到路易斯共振结构。他们计算了“共价”结构和“离子”结构的权重。分析表明，仅凭单个电子分别定域在两个中心的结构（即所谓的“双自由基”共振式）并不能说明双自由基特征，因为即使在典型的共价键（如H₂）中，这种构型也占有50%的概率。双自由基特征必须通过比较“共价”和“离子”贡献的权重来共同判断。在闭壳层极限，两者权重相等；在完美双自由基极限，“共价”结构占主导（权重为1），“离子”结构权重为0。因此，一个分子的双自由基特征可以通过其波函数偏离闭壳层（共价与离子各半）的程度来量化。

3.7 与ST能隙的关联

尽管理论上ST能隙与双自由基特征没有简单关系，但研究人员将计算得到的所有分子的ΔE_ST^v（MRCI）与β（CASSCF(2,2)）进行了绘图关联。结果显示，数据点大致分布在理论预测的区域内。不交盖的双自由基（如研究的四元环体系）主要分布在负ST能隙、中低β区域，而典型的交盖双自由基O₂则位于正ST能隙、高β区域。这个关联图直观地展示了不同双自由基物种在电子结构谱图上的位置，并验证了理论分析的合理性。

结论与讨论

本研究通过结合多种计算方法和理论分析，为单线态双自由基(oid)的电子结构提供了一个统一、直观的描述框架。研究的主要结论和重要意义在于：

1. 1.
   提出了新的化学直观描述：研究首次提出，对于中低程度双自由基特征的物种，可以利用“键”和“反键”的概念来描述自由基电子间的相互作用。这为化学家理解这些复杂的多参考体系提供了熟悉的语言，将抽象的波函数系数与直观的化学图像联系起来。
2. 2.
   澄清了关键概念：研究明确了单线态-三重态能隙（ΔE_ST）与双自由基特征（β）之间的复杂关系，指出两者并非简单对应，并给出了其依赖的物理参数（库仑积分和交换积分）。这纠正了可能存在的误解，即认为小的或负的ST能隙直接意味着高的双自由基特征。
3. 3.
   解决了长期争议：应用新框架，研究澄清了关于臭氧（O₃）和二氮化二硫（S₂N₂）等分子双自由基特征的争论。研究表明，臭氧具有中等的双自由基特征，而S₂N₂的双自由基特征较低，这为其化学反应性提供了合理解释。
4. 4.
   提供了实用的计算指南：研究评估了多种计算方法的性能，表明CCSD(T)在获取可靠结构和能量方面表现优异，而PBE-D3等纯泛函在结构优化和ST能隙估算上也能提供合理结果，这为研究更大、更真实的双自由基分子提供了方法学参考。
5. 5.
   建立了统一描述符关联：通过系统计算多种双自由基指标（ST能隙、β、LUNO占据数、键级、共振结构权重）并展示其关联，研究促进了不同文献中使用的各种描述符之间的比较和理解，有助于该领域研究走向更标准化和系统化。

总之，这项研究不仅深化了对双自由基基本电子结构的理解，而且架起了一座连通高深计算化学与直观化学思维的桥梁。它提出的“键与反键”框架有望成为教育和研究中描述和理解这类重要反应性物种的有力工具，并可能激发在材料科学、催化反应和分子电子学等领域的新应用探索。