DCM-IFC中的反卡莎发射：对第三类分离波函数假说的计算评估

《Physical Chemistry Chemical Physics》：Anti-Kasha emission in DCM-IFC: computational evaluation of the type III separated wavefunction hypothesis

　　

研究人员对DCM-IFC分子的反卡莎（anti-Kasha, AK）发射现象进行了系统的理论计算评估。卡莎规则指出，给定多重度的发光电子能级是该多重度的最低激发能级，其例外情况通常由大的S1–S2能隙导致，如薁（azulene）。近年来，报道了一些具有中等S1–S2能隙的大柔性分子表现出AK发射，并被解释为两态波函数高度空间分离的结果。近期，由二氰亚甲基-4H-吡喃（DCM）与集成荧光素-色烯（IFC）构建的荧光团中，螺内酯开/关开关被证明可调控AK行为，其开型分子表现出双发射，分别归属于S1与S2态。研究人员在隐式溶剂条件下，计算了Franck–Condon区域及通向S2/S1交叉区域的插值路径，采用含时密度泛函理论（time-dependent density functional theory, TD-DFT）及多种泛函和高阶波函数理论，发现S1态主要为局域激发特征且振子强度更高，而S2态为电荷转移（charge-transfer, CT）特征且较弱。从头算结果表明两态吸收光谱重叠，并且从S2势能面到S2/S1交叉区域存在无势垒的内转换通道，这意味着S2态的快速去布居会抑制AK发射。此外，各激发态最小值处的S1与S2能量差也不利于在发射时间尺度内实现热平衡。综上，计算结果显示实验观测到的异常发射需要重新解释。

该研究发表于《Physical Chemistry Chemical Physics》，针对DCM-IFC开型分子的反卡莎（AK）发射假说展开理论验证。卡莎规则作为光化学的核心原则，认为给定多重度的发光能级是最低激发态，但近年发现的第三类AK分子因S₁与S₂波函数高度空间分离而被认为可能突破该规则。然而，已有量子动力学模拟显示此类分子的S₂态寿命极短，与实验提出的AK机制存在矛盾。本研究旨在通过高精度量子化学计算，解析DCM-IFC的激发态结构与动力学，验证其AK发射的可能性。

研究人员采用的关键技术方法包括：在隐式溶剂模型下，使用M06-2X、CAM-B3LYP、ωB97X-D等多种泛函进行TD-DFT几何优化与频率分析；结合单参考二阶耦合簇（CC2）与多参考N-电子价态微扰理论（NEVPT2）提高激发态能量精度；通过量子热浴从头算分子动力学（QT-AIMD）采样基态构象并模拟吸收光谱；沿内坐标线性插值（LIIC）构建S₂至S₂/S₁锥形交叉区域的势能面，并采用校正线性响应（cLR）方法修正溶剂效应。

研究结果如下：

3.1 480 nm泵浦是否激发S₂态：垂直激发能计算表明，无论中性还是阴离子形式，S₁均为亮态（振子强度约2.0），S₂为暗态（振子强度约0.2），且480 nm泵浦能量更接近S₁态。吸收光谱模拟显示，经经验红移校正后，480 nm可同时覆盖S₁和部分S₂态吸收带，但S₁吸收截面显著更高。

3.2 是否存在通向S₂/S₁交叉区域的势垒：LIIC路径计算表明，从中性分子和阴离子的Franck–Condon区域到近简并几何结构均存在无势垒通道，且在CC2与NEVPT2水平上交叉点更接近初始结构。这表明S₂→S₁内转换可在亚皮秒尺度完成，与实验推测的~45 ps S₂寿命不符。

3.3 激发态势间热平衡的可能性：cLR校正后的发射能显示，S₁与S₂在各自最小值的能量差仍高于室温热能，且两态几何结构差异显著，难以在发射时间尺度内实现热平衡。QT-AIMD采样进一步表明，即使在S₂最小值附近，态特征切换的概率也极低（≤11%）。

3.4 计算方法的局限性：研究未直接进行非绝热分子动力学模拟，且溶剂处理仅采用隐式模型，可能低估特定溶质-溶剂相互作用；锥形交叉区域仅基于线性插值路径表征，未完全探索多维交叉缝的拓扑结构。

讨论与结论部分指出，DCM-IFC开型分子的S₁态为亮态且S₂态为暗态，且存在快速的无势垒S₂→S₁衰减通道，因此实验观测的双发射难以用AK机制解释。研究人员强调，仅凭波函数的空间分离不足以支撑AK发射，相关荧光探针的异常信号可能源于其他机制。这一结论对基于DCM-IFC motif的光谱传感体系设计具有重要警示意义，建议未来结合多维实验与动态模拟进一步验证。