《Journal of Chemical Theory and Computation》：Modeling Xanthophyll Excited States via Cost-Effective Quantum Chemistry methods and Property-Based Diabatization
编辑推荐：
黄体素类胡萝卜素在光合作用光捕获与生物系统光保护中发挥重要作用，然而由于其扩展的π共轭体系与多参考态特性，以可接受的计算成本准确描述其激发态仍具挑战性。本研究中，研究人员比较了DFT/MRCI（密度泛函理论/多参考组态相互作用）、FOMO–CI（浮动占据分子轨
黄体素类胡萝卜素在光合作用光捕获与生物系统光保护中发挥重要作用，然而由于其扩展的π共轭体系与多参考态特性，以可接受的计算成本准确描述其激发态仍具挑战性。本研究中，研究人员比较了DFT/MRCI（密度泛函理论/多参考组态相互作用）、FOMO–CI（浮动占据分子轨道组态相互作用）以及MRSF-TDDFT（混合参考态自旋翻转含时密度泛函理论）对一组代表性黄体素激发态的描述。通过采用基于多重性质的非绝热态构建（MPD）方案对各方法的电子态进行一致性表征，分析了垂直激发能及沿键长交替坐标的势能面。研究发现，三种方法对低能激发态及其势能曲线给出了显著一致的描述。DFT/MRCI与半经验FOMO–CI在各类黄体素中表现出高度一致的行为，而MRSF-TDDFT的性能则依赖于交换-相关泛函的选择。总体而言，计算结果表明在Franck–Condon点处，共价的2Ag–态位于离子性的1Bu+态之下。该分析为黄体素的电子结构提供了普适性认识，包括结构畸变与环境效应诱导的共价/离子构型相互作用。MPD框架实现了态性质的直接比较，凸显了方法间的相似性并阐明了残余差异的起源。本工作为理解黄体素光物理学建立了统一的理论框架，并为采用成本效益方法建模其激发态提供了实践指导。
本研究旨在解决扩展共轭体系激发态计算的准确性-成本权衡难题，聚焦于黄体素类胡萝卜素这一兼具强多参考特性、几何驱动态混杂及环境敏感性的复杂生物分子体系。该类分子在光合作用中承担光捕获、能量传递及非光化学淬灭（NPQ）等关键功能，其低能激发态（特别是暗态2Ag–与亮态1Bu+）的能量排序、耦合关系及驰豫路径直接决定其生物学功能，但精确描述这些态需处理共价-离子构型相互作用、几何依赖的态混杂及环境诱导极化等多重挑战。传统多参考方法如CASSCF（完全活性空间自洽场）结合CASPT2或NEVPT2等微扰理论虽能严格处理此类问题，但需极大活性空间导致计算成本剧增；线性响应TDDFT因无法描述双激发特征的2Ag–态而失效；MRSF-TDDFT虽通过扩展激发空间缓解了自旋污染问题，但泛函依赖性仍是限制；DFT/MRCI作为精度与效率的折中方案，其最新R2022参数化的泛化能力尚待验证；半经验的FOMO–CI方法虽有计算优势，但其参数的可转移性缺乏系统评估。为此，本研究首次将多重性质非绝热态构建（MPD）框架应用于该类体系，实现了跨方法的电子态特征一致性追踪。
研究发表于《Journal of Chemical Theory and Computation》。所得主要结论如下：DFT/MRCI与FOMO–CI对黄体素激发态描述高度一致，均支持Franck–Condon点处2Ag–低于1Bu+的能级排序；MRSF-TDDFT结果强烈依赖泛函选择，BH&HLYP泛函系统性高估暗态能量导致定性错误的态排序，而DTCAM-XIV泛函给出与DFT/MRCI最接近的结果；FOMO–CI以最低计算成本实现了与参考数据（DSRG-MRPT2）偏差小于0.12 eV的精度，但远离平衡区时势能面曲率被高估；环境效应与分子对称性破缺共同诱导S1/S2态混杂并产生分子内电荷转移（ICT）特征，MPD框架可有效分离本征态性质与混合效应。该研究为复杂生物环境中类胡萝卜素光物理模拟确立了具有物理可靠性的成本效益策略。
主要技术方法包括：（1）多重性质非绝热态构建（MPD）方案，以Mulliken原子布居为观测集，通过最小化Frobenius范数残差获得绝热-非绝热变换矩阵，实现电子态特征的平滑追踪；（2）DFT/MRCI（R2022参数化/BH&HLYP/def2-SVP），基于Kohn–Sham轨道的多参考组态相互作用方法；（3）FOMO–CISD(6,9)方法，采用针对黄体素再参数化的AM1半经验哈密顿量，浮动占据宽度0.1 Hartree；（4）MRSF-TDDFT方法，测试BH&HLYP及DTCAM系列（AEE、VEE、VAEE、XI、XIV、STG）共七种泛函（均使用def2-SVP基组）；（5）QM/MM（量子力学/分子力学）静电嵌入方案，以AMBER ff14SB力场描述蛋白质环境，采用ONIOM方案进行几何优化；（6）B3LYP/6–31G(d)水平的几何优化与键长交替（BLA）约束扫描（步长0.025 ?）。样本队列包括五种结构多样性黄体素：紫黄质（Vio）、黄体素（Lut）、玉米黄质（Zea）、角黄质（Can）和深红菌黄素（Sal），其对称/非对称取代模式覆盖羟基、酮基及环氧官能团。
研究结果部分：
3.1 激发能比较
通过三种方法计算了五种黄体素在B3LYP/6–31G(d)基态几何上的垂直激发能（图3）。各方法对给定态在不同黄体素间的能级趋势（Vio > Lut > Zea > Can > Sal）给出一致预测，Bu+态顺序与实验测量一致。FOMO–CISD与DFT/MRCI结果高度吻合；MRSF-TDDFT/BH&HLYP系统性高估所有黄体素激发能，且错误预测2Ag–/1Bu+态反转（2Ag–高于1Bu+约0.4 eV）。DTCAM-XIV泛函表现最优，能再现定性趋势及定量激发能，但仍存在微小态序反转（图S10）。Sal的2Ag–跃迁偶极矩异常偏大，表明在B3LYP几何上2Ag–与1Bu+强烈混杂，该混杂在MRSF-TDDFT/DTCAM-XIV优化几何上显著减弱（图S9–S12）。
3.1.1 键长交替方向的势能曲线
沿BLA坐标的基态约束优化扫描（图4）显示，三种方法均预测：随BLA减小，2Ag–态能量陡峭下降并于负BLA区与基态形成强回避交叉；两个Bu态在低于Franck–Condon点的BLA值处发生交叉。MPD非绝热表示（以S0势能面最小值为参考点）实现了态特征的连续追踪。方法间最大差异出现在基态势能面及其最小值位置：FOMO–CISD预测S0最小值BLA≈0.056 ?，显著低于DFT/MRCI与MRSF-TDDFT的0.076–0.106 ?范围。DFT/MRCI与MRSF-TDDFT在正负BLA区域均表现出良好一致性，仅在负BLA区存在已知偏差。FOMO–CISD在远离Franck–Condon区时势能面曲率偏大，但不改变交叉点位置。MRSF-TDDFT/BH&HLYP的非绝热能曲线与DFT/MRCI形状相近，但暗态能量系统性高估；对2Ag–和1Bu–分别施加–0.45 eV和–0.55 eV能量修正后，可恢复与其余方法一致的交叉行为（图5）。
3.1.2 与从头算多参考计算比较
基于CAM-B3LYP/cc-pVDZ几何的DSRG-MRPT2参考数据（Table 1a）表明，该几何下亮态（最大跃迁偶极矩）对Vio和Lut为S2，对Zea为S1，且三体系S1–S2间隙约0.02 eV。FOMO–CISD与DFT/MRCI预测能隙在0.12 eV参考误差范围内，与参考吻合；MRSF-TDDFT/DTCAM-XIV将暗态置于亮态上方约0.2 eV，偏差稍大。在B3LYP几何上（Table 1b），仅MRSF-TDDFT预测S1为亮态（实际几乎简并）。由于B3LYP几何BLA更小（~0.08 vs ~0.10 ?），暗态能量相对更低，支持Franck–Condon区2Ag– < 1Bu+的能级排序预期。
3.2 激发态性质
3.2.1 化学结构效应
对比非对称黄体素（Lut）与对称黄体素（Can）的态偶极矩与跃迁偶极矩（图6）：Lut基态偶极矩近乎为零，但在S1/S2交叉区出现显著尖峰，表明ICT特征的出现；该峰在非绝热表示中消失，证实ICT源于态混杂而非本征态性质。Can因C2h对称性约束，各态偶极矩保持为零，ICT被对称性禁阻。
3.2.2 环境扰动效应
以Orange Carotenoid Protein（OCP）中的Can为模型（图7），通过三种条件（真空优化/真空计算、OCP环境优化/真空计算、OCP环境优化/OCP环境计算）对比分析。结构畸变使基态获得非零偶极矩，S1/S2交叉区形成带电S1态；加入环境静电后，各态偶极矩整体增强，特别是1Bu+态在非绝热表示中保留显著偶极矩，表明环境极化对离子态更强。对称情况下S1/S2跃迁偶极矩交换急剧发生，而对称破缺条件下交换更为渐进，多组态混杂点呈现中间跃迁偶极矩值。
讨论总结：
研究确立了三种成本效益方法描述黄体素激发态的一致性与局限性。DFT/MRCI与FOMF–CI的高度一致性支持后者作为高效的计算替代方案；FOMO–CI势能面曲率偏大的特征解释了其在超快动力学模拟中的成功应用，尽管可能需针对非平衡几何优化参数。MRSF-TDDFT的泛函依赖性要求谨慎选择，DTCAM-XIV因其同时优化价层激发能与电离势的设计理念，在正BLA区域表现优异。MPD框架的核心价值在于消除几何依赖的态标记混淆，使方法比较建立在物理本质而非表观能级标签之上。环境因素研究揭示了蛋白质环境通过两种途径调制ICT：结构畸变打破对称性启用S1/S2混杂，静电极化则进一步稳定离子特征并加剧态耦合，这对理解NPQ等生物学功能至关重要。
研究结论：
本研究比较了DFT/MRCI、FOMO–CISD和MRSF-TDDFT三种方法对代表性黄体素低能激发态的描述，聚焦于Franck–Condon点垂直激发能及键长交替坐标上的势能曲线。借助统一的多重性质非绝热态构建方案，实现了跨方法的电子态特征一致性比较。
结果表明方法间具有显著一致性，尤其是DFT/MRCI与半经验FOMO–CISD。两者对暗态S1/2Ag–与亮态S2/1Bu+态在Franck–Condon点附近的能量排序及整个黄体素系列的垂直激发能给出一致预测。值得指出的是，FOMO–CISD使用共同的活性空间及原始黄体素参数化，对所有体系提供了平衡描述，表明其在该类胡萝卜素中具有良好的可转移性，尽管紫黄质呈现较高偏差。其主要局限在于远离Franck–Condon区时系统性高估基态与激发态势能面曲率，但这并不显著改变交叉区域位置，这解释了FOMO–CISD成功模拟类胡萝卜素超快激发态动力学的原因。
DFT/MRCI与FOMO–CISD激发能均与CAM-B3LYP基态几何上的多参考微扰理论参考值（偏差在0.12 eV内）吻合良好，这使两种方法对2Ag–和1Bu+态能量的不确定性低于0.2 eV。
对于MRSF-TDDFT，结果展示了强烈的泛函依赖性。BH&HLYP泛函系统性高估暗态激发能并导致错误态序；相比之下，DTCAM-XIV泛函给出与DFT/MRCI接近的激发能和势能曲线，尤其在正BLA的Franck–Condon相关区域。负BLA区偏差更为显著，此处MRSF-TDDFT可靠性较低。非绝热表示揭示BH&HLYP的主要偏差源于系统性的态特定能量偏移，在非绝热 representation 中制定修正策略可能有助于改进势能面。
所采用的非绝热态构建方案为跨方法和体系的态特征与交叉点提供了统一描述，实现了势能面的严格比较。利用该方案，研究展示了各方法描述类胡萝卜素激发态的高度相似性。更 generally，该工作表明非绝热方案可辅助不同几何和理论水平下的电子结构分析。
所有考虑的方法均能一致预测类胡萝卜素的激发态性质。特别是，它们可靠地捕获了反映绝热态间特性交换的跃迁偶极矩趋势。对称破缺与环境效应诱导的态混杂分析，为S1和S2态中分子内电荷转移（ICT）特征的出现提供了清晰机理：该特征产生于态强耦合区域，而非本征态性质。
上述发现为解释黄体素光物理学提供了统一基础，并为类胡萝卜素光物理模拟中的方法选择建立了实践指导原则。
