揭示自旋-轨道耦合的结构起源：SOPEL方法

《Physical Chemistry Chemical Physics》：Revealing the structural origin of spin–orbit coupling with SOPELOpen Access

　　

自旋-轨道耦合（spin–orbit coupling，SOC）是许多光化学过程的关键要素，但它通常被报告为一个不透明的标量值，缺乏化学可解释性。在此，研究人员引入了SOPEL（Spin–Orbit Population ELement）的概念，这是一种对自旋-轨道耦合矩阵元（spin–orbit coupling matrix elements，SOCMEs）进行原子分辨分解的方法，可以揭示其在分子中的结构起源。研究人员进一步提出了基于实原子轨道的非零耦合实用选择定则，为局域轨道特征与自旋-轨道耦合之间提供了直观的联系。总的来说，该概念将SOCME从一个黑箱式的电子结构输出转化为一种化学可解释的量。通过定位自旋-轨道耦合的结构起源，SOPEL不仅为从结构和运动角度理解自旋禁阻过程提供了一个实用框架，也为化学引导的分子设计提供了依据。该方法在碘取代的BODIPY、甲醛、二苯甲酮和三芳胺-对苯二甲腈上的应用表明，它能够：(i) 识别产生自旋-轨道耦合的原子和片段，(ii) 表征自旋-振动机制，(iii) 捕捉耦合起源在分子骨架上的再分布，以及 (iv) 解释为何仅凭重原子的参与并不能保证强耦合。

论文解读文章

研究背景与问题

系间窜越（intersystem crossing，ISC）是光催化、有机发光二极管（OLED）、光动力疗法等领域众多光化学过程的核心。从理论上看，ISC受控于自旋-轨道耦合（spin–orbit coupling，SOC）与振动态密度的相互作用。然而，当前的理论计算虽然能提供SOC矩阵元（SOCMEs），但这些数值如同一个“黑箱”，难以转化为具有化学洞察力的信息。具体而言，研究人员无法直观地知道是分子中哪些原子区域和局域轨道特征控制了耦合，以及结构变化如何重新分配这种耦合。现有的定性分析方法，如分子轨道（MO）或自然跃迁轨道（NTO）可视化，虽能提供一些信息，但无法提供SOCME的定量原子级贡献。因此，亟需一种能将SOCME与空间可分辨的区域联系起来的新方法，从而指导分子设计。

研究内容与结论

为了解决上述问题，研究人员引入了一个全新的概念框架——SOPEL（Spin–Orbit Population ELement，自旋-轨道布居元素）。该框架由两部分组成：一是将SOCME分解为单个原子贡献的SOPEL；二是基于实原子轨道（AO）语言表述的非零SOCME选择定则。通过这一框架，研究人员成功地将SOC从一个抽象的数值转化为了一个具有化学可解释性的量。他们将该方法应用于碘代BODIPY、甲醛、二苯甲酮和三芳胺-对苯二甲腈（TAA–TPN）四个代表性体系，验证了其有效性。研究表明，SOPEL不仅能精确识别出产生SOC的关键原子或片段，还能揭示自旋-振动耦合的微观机制，并发现SOC的结构起源可以在分子内迁移。更重要的是，该工作解释了为什么仅仅在分子中引入重原子并不总能保证强的SOC，关键在于重原子上是否存在合适的局域轨道通道。这项工作为理解和调控自旋禁阻的光化学过程提供了一个强大的新工具，并为基于化学直觉的分子设计开辟了新途径。

关键技术方法

本研究采用的主要技术方法包括：

1. 1.
   SOPEL分析方法：研究人员定义了SOPEL，即原子A对总SOCME的复数贡献，通过将单电子自旋-轨道积分（SO integrals）按原子中心基函数进行划分（类比于Mulliken和L?wdin布居分析），并与无自旋跃迁密度矩阵相结合，实现了SOCME的原子级分解。
2. 2.
   实原子轨道选择定则：研究人员推导并归纳了一套通用的选择定则，用以判断不同实原子轨道之间的SO积分是否非零。该定则以一种简洁的图表形式呈现，无需复杂的MO或NTO可视化即可快速应用。
3. 3.
   计算与可视化程序：研究人员开发了与ORCA程序的CIS/TDDFT输出文件接口的Python代码，用于执行SOPEL分析并将原子贡献以SOPEL范数的形式进行可视化。

研究结果

• •
  BODIPY：尽管碘持续参与，但振动激活了SOCC
  研究人员发现，在一个碘代BODIPY衍生物中，一种远离碘原子的氢原子面外振动模式能显著增强T₂和S₁态之间的自旋-轨道耦合常数（SOCC）。SOPEL分析表明，整个SOCC几乎完全由碘原子贡献。进一步分析碘原子的原子轨道（AO）贡献后发现，在优化几何构型下，两个主导MO均只含p_z轨道，而p_z-p_z的SO积分为零。振动导致p_x和p_y轨道掺入，使得不同p轨道间的非零SO积分成为可能，从而激活了耦合。这表明，重原子的存在本身不足以保证强SOC，关键在于振动能否打开其上的局域轨道通道。
• •
  甲醛：氧中心SOC通道的丧失
  研究人员研究了甲醛沿对称平面H-C-H剪式振动的T₁和S₂态SOC。SOPEL分析显示，SOCC的单调递减源于氧原子SOPEL贡献的消失，而碳原子的贡献基本不变。通过分析氧原子的AO贡献发现，随着振动，氧原子在LUMO+1中的p轨道特征逐渐丧失，转而变为s轨道特征，而s轨道的SO积分为零。这揭示了SOCC降低并非由于耦合的分散，而是单一主导原子通道的关闭。
• •
  二苯甲酮：SOCC几乎恒定下的环间SOPEL再分布
  研究人员分析了二苯甲酮沿最低频振动模式的T₅和S₁态SOC。虽然总的SOCC变化很小，但SOPEL分析揭示了一个有趣的现象：耦合的结构起源在两个苯环之间来回振荡，一些原子失去贡献的同时另一些原子获得贡献。这说明即使总的SOC强度保持不变，其空间起源也可能发生巨大改变，揭示了一个隐藏的自旋-轨道耦合自由度。
• •
  TAA–TPN：用于分子设计的片段级指导
  研究人员对一种典型的热活化延迟荧光（TADF）分子TAA–TPN衍生物进行了SOPEL分析。结果显示，S₁和T₁态之间的SOCC主要来源于二氰基苯基团及其相邻的苯基。基于此，研究人员提出并验证了一个修改方案：移除一个不贡献SOC的苯基，发现修改后分子的SOCC几乎没有变化（从0.56 cm^-1变为0.53 cm^-1），证明了SOPEL可以为简化分子设计、保持性能提供直接的片段级指导。

总结与讨论

本研究提出的SOPEL概念成功地将SOCME从一个不透明的标量转化为一个具有化学可解释性的量，能够精确揭示SOC在分子中的原子级结构起源。通过对四个代表性案例的研究，SOPEL展示了其强大的分析能力：(i) 识别出携带耦合的关键原子或片段；(ii) 阐明自旋-振动耦合的机制，区分了“开启新通道”与“关闭旧通道”两种不同的调制方式；(iii) 捕捉到耦合起源在分子骨架上的动态迁移；(iv) 解释了重原子效应并非总是有效的原因，并能为分子设计提供直接指导。总之，SOPEL提供了一个将SOC与分子结构和运动直接联系起来的实用框架，对于深入理解自旋禁阻过程以及进行化学引导的分子设计具有重要意义。