变分激发态计算新突破：基于约束的轨道优化方法(COOX)在准确性、稳定性和适用性上的系统性超越

《Communications Chemistry》：A detailed comparison of ΔSCF methods with the constraint-based orbital-optimized excited state method

本文聚焦量子化学中快速准确计算电子激发态这一核心挑战。为解决线性响应理论(LR-TDDFT)在描述电荷转移(CT)、芯激发及双重激发等方面的局限，以及传统ΔSCF方法存在的收敛不稳定、可能收敛到错误态等难题，研究团队系统比较了已建立的ΔSCF方法与其最新提出的基于约束的轨道优化激发态方法(COOX)。通过将COOX应用于ΔSCF所定义的特定轨道旋转，提出了可直接比较的ΔCOOX方法，并对两者在精度和稳定性上的性能差异进行了深入探讨。研究涵盖价态、芯态、里德伯态、双重激发及电荷转移激发等多种分子体系，结果表明COOX方法在保持可比结果质量的同时，在变分激发态计算的整体稳定性方面显著优于传统ΔSCF方法，为激发态计算提供了更可靠、高效的新工具。

　　

在化学、材料科学和生命科学等领域，深入理解分子的电子激发态对于揭示光化学反应机理、设计新型光电材料乃至解读复杂光谱实验数据至关重要。因此，发展能够快速、准确计算电子激发态的理论方法是现代量子化学的核心挑战之一。长久以来，线性响应含时密度泛函理论(LR-TDDFT)因其在计算成本与精度间较好的平衡而成为最广泛应用的方法。然而，它存在一些众所周知的短板，例如在描述芯激发、长程电荷转移激发以及双重激发等方面存在困难或根本性局限。为了突破这些限制，研究者们将目光投向了轨道优化的激发态方法，这类方法通过自洽场(SCF)过程变分地直接确定激发态，其中ΔSCF及其改进版(如MOM, iMOM)是代表性方案。ΔSCF方法确实在改善上述LR-TDDFT的缺陷方面展现了潜力，甚至能以近似于基态SCF计算的成本处理双重激发，并且无需求解昂贵的Z-vector方程即可计算核梯度。但ΔSCF方法自身也面临严峻挑战：其激发由明确定义的单一轨道旋转(如HOMO→LUMO)给出，难以描述涉及多个占据或虚轨道混合的激发；更重要的是，这些方法常常遭受收敛缓慢、甚至完全无法收敛的困扰，还可能发生“变分崩溃”——即计算退回到初始的基态。尽管近十年来出现了如平方梯度最小化(SGM)等改进算法以提升稳定性，但它们通常计算代价更高，且不一定能收敛到目标态。这些固有缺陷限制了ΔSCF类方法的广泛应用和可靠性。

为了从根本上解决这些问题，由J?rg Kussmann、Yannick Lemke等人组成的研究团队在《Communications Chemistry》上发表了他们的最新研究成果。他们提出并系统评估了一种全新的基于约束的轨道优化激发态方法——COOX。与直接违反Aufbau原理的ΔSCF不同，COOX方法植根于约束密度泛函理论(cDFT)框架。它通过引入一个由(简化的)TDDFT跃迁振幅构建的约束势，在保持Aufbau原理的前提下，通过直接操纵分子轨道来“旋转”特征向量，从而实现目标激发态的变分优化。这种方法不仅保持了ΔSCF能够改进LR-TDDFT短板、计算成本相当等优点，还带来了几个关键优势：显著提升的收敛稳定性和速度、避免变分崩溃的内在特性，以及描述任意(包括混合)轨道旋转激发的能力。为了与传统的ΔSCF进行公平、直接的比较，研究者们特别设计了ΔCOOX方案，即让COOX使用与ΔSCF完全相同的特定轨道旋转来构建约束，从而在相同“起跑线”上比拼性能。

本研究主要依赖于量子化学计算与理论分析。核心方法包括：1) 理论方法对比：详细阐述了传统ΔSCF(特别是MOM和iMOM)与新型ΔCOOX/COOX方法的理论基础、共同点与本质差异。2) 系统性基准测试：使用多种密度泛函(如PBE, PBE0, LRC-ωPBE, ωB97X)和基组，对一系列具有代表性的分子体系进行了大规模计算比较，涵盖了苯的低位激发态、铍的价态与里德伯态、标准体系的双重激发、乙烯-四氟乙烯复合物的分子间电荷转移激发、QUEST数据库中的分子内电荷转移激发、跨越元素周期表的K-、L-、M-边芯激发，以及复杂环境(如一氧化碳嵌入C₆₀) 中的激发态。3) 性能评估指标：不仅比较激发能精度，还重点考察了SCF迭代步数(收敛速度)、收敛成功率、是否收敛到预期态等稳定性和鲁棒性指标，并引入了“污染/偏差”量化参数来评估激发态与指定初始猜想的匹配程度。4) 扩展应用：展示了COOX在嵌入方案(e-COOX) 中处理环境效应的便利性，以及作为ΔSCF预处理器的潜在价值。

结果：

理论与概念差异

研究首先从理论上厘清了ΔSCF与ΔCOOX的核心区别。ΔSCF的实现依赖于在SCF过程中手动违反Aufbau原理，改变特定轨道的占据数，其本质是在优化激发态的同时重新优化“特征值”(占据数)。而ΔCOOX则保持Aufbau原则，通过约束势来移动涉及激发的轨道的能级，从而间接地优化“特征向量”(分子轨道)以实现电子从指定占据轨道到虚轨道的转移。这种根本差异被认为是ΔSCF诸多收敛问题的根源，而ΔCOOX则因其不违反电子占据的基本规则而天然地避免了变分崩溃。

苯的低位激发态：稳定性与收敛性的压倒性优势

对苯分子所有单激发HOMO?x → LUMO+y (x, y = 0, …, 6)的计算提供了令人信服的证据。在收敛性方面，ΔCOOX表现极为稳健，所有计算均在8-10步内收敛。相比之下，MOM和iMOM不仅需要多得多的迭代步数(经常超过50步)，还分别有1个和4个态完全无法在200步内收敛。在精度上，大多数态三者结果相近，但ΔCOOX预测的激发能普遍略高。分析表明，激发能差异与“污染/偏差”值密切相关：ΔCOOX的污染/偏差严格为零(由约束保证)，而MOM/iMOM的结果则显示不同程度的偏离，偏离越大，与ΔCOOX的能差也越大。这说明ΔSCF可能收敛到一个发生了额外轨道松弛或混合的态，而ΔCOOX则严格实现了指定的电子跃迁。

铍的价态与里德伯态

对铍原子2s2p和2s3s激发的计算表明，ΔCOOX和iMOM能提供质量相近的激发能。但iMOM的收敛行为依然不稳定且不可预测，所需迭代步数在5到80之间剧烈波动，而ΔCOOX则在5-8步内稳定收敛。研究还发现，当使用高精确交换比例泛函时，iMOM结果中指定虚轨道的占据偏差显著增大，暗示单轨道跃迁图像可能不再完全适用。此时，COOX可以采用与处理芯激发类似的芯-价分离(CVS)方案来构建约束，从而获得更合理的描述。

双重激发

对一组标准双重激发的测试显示，在排除个别未收敛或收敛到错误态的情况后，ΔCOOX与iMOM/MOM给出的激发能总体相近。但ΔCOOX再次展现出无与伦比的收敛稳定性，最多仅需9步即可收敛，而iMOM/MOM则波动很大。这确认了ΔCOOX在计算这类LR-TDDFT无法直接描述的过程时，同样具有显著的计算效率优势。

电荷转移激发

电荷转移态是展示轨道优化方法优势的经典场景。在乙烯-四氟乙烯复合物的分子间长程CT研究中，ΔCOOX(配合适当的α和β约束分离)在整个距离范围内都表现出稳定、正确的渐进行为(-1/R)，且收敛快速。而MOM和iMOM在短距离处面临严重的数值不稳定性：iMOM多次无法收敛，MOM则可能崩溃到局域Frenkel激发。在>5.5 ?后，所有方法才都表现出稳定收敛。对于QUEST数据库的分子内CT基准测试，ΔCOOX(尤其是使用长程校正杂化泛函LRC-ωPBE和ωB97X时)的总体精度略优于或与iMOM相当，平均绝对误差在0.24-0.25 eV。一个有趣的观察是，iMOM对长程分子间CT的激发能表现出惊人的泛函无关性，而ΔCOOX则对泛函更敏感(需要长程校正泛函以获得正确渐进行为)。作者指出，前者是由于非Aufbau解中潜在的误差抵消，其物理根源值得未来深入研究。

芯激发

对一系列分子K-、L-、M-边激发的重新评估证实了COOX方法的卓越性能。尽管所有方法都能给出相对误差1%以内的较好精度，但COOX(使用CVS-TDA差分密度)和ΔCOOX在数值稳健性上更胜一筹。iMOM在计算MoS₄^2-和WCl₆时失败，而COOX方案则能稳定收敛，且所需迭代步数更少。

复杂环境中的激发态

通过e-COOX嵌入方案，研究者展示了将COOX应用于复杂量子环境(如蛋白质、溶剂化壳层)的便捷性。以C₆₀富勒烯包裹的一氧化碳分子为例，在完整体系中，目标激发态对应于第455个简化的TDA根，或HOMO?28 → LUMO+6的跃迁，难以直接指定。而e-COOX和e-ΔCOOX方案允许从孤立子体系(CO)轻松获得约束(分别为第一TDA根或HOMO→LUMO)，然后嵌入到完整体系进行优化，无需繁琐的手动轨道指认，即可黑箱式地获得正确的激发态和电子密度差。

作为预处理器

将ΔCOOX作为MOM/iMOM的预处理器被证明是一个行之有效的策略。以硝基苯的HOMO→LUMO激发为例，直接使用iMOM需49步收敛，MOM需73步且收敛到了错误的双重激发态。而先用ΔCOOX(即使以较松标准)进行预优化，再用其轨道作为MOM/iMOM的初始猜测，可将后续迭代步数减少至27步，并确保MOM收敛到正确的目标态。这为在已有ΔSCF代码中快速引入COOX的优越性提供了一条捷径。

结论与讨论：

本研究通过详尽的理论分析和系统的数值比较，清晰地表明基于约束的轨道优化激发态方法(COOX)，特别是其ΔCOOX变体，在绝大多数情况下明显优于传统的ΔSCF方法。其核心优势源于方法论的根本不同：ΔCOOX在保持Aufbau原理的前提下，通过约束势优化轨道，从而根除了导致ΔSCF不稳定的“变分崩溃”风险，并确保了计算始终朝向预设的激发目标。

结论可归纳为以下几点：

1. 1.
   卓越的稳定性与效率：ΔCOOX在各类激发(价态、芯态、CT态、双重激发)的计算中，都展现出远优于MOM/iMOM的收敛速度和成功率，通常能在10步以内稳健收敛，而后者则常需数十步甚至失败。
2. 2.
   可靠的精度：在能够直接比较的情况下，ΔCOOX提供的激发能精度与ΔSCF方法相当，甚至在某些方面(如与高污染/偏差的ΔSCF结果相比)可能更接近真实的变分态。
3. 3.
   功能与应用的扩展性：标准的COOX方法利用TDDFT振幅构建约束，能够天然描述任意混合的激发，突破了ΔSCF对单轨道跃迁的限制。其嵌入方案(e-COOX)使得在复杂环境中进行黑箱式激发态计算成为可能。
4. 4.
   易于实现：COOX的实现与广泛使用的MOM/iMOM一样简便，尤其对于已具备cDFT功能的代码，仅需在构建和对角化Fock矩阵之间对现有SCF流程进行一次修改，并且与DIIS等标准收敛加速技术完全兼容。

尽管本研究未与SGM等更复杂的ΔSCF算法进行直接对比，但基于文献报道和本文展示的COOX原理优势，作者预期ΔCOOX在收敛速度和单步迭代成本上依然具备竞争力。未来，将ΔCOOX与限制性开壳层Kohn-Sham形式结合，以及发展基于COOX的非正交组态相互作用方法，将是富有前景的研究方向。

综上所述，这项研究不仅证实了COOX/ΔCOOX在处理挑战性激发态问题上的强大能力，更重要的是，它为解决长期困扰ΔSCF应用的稳定性和可靠性难题提供了一个优雅而强大的解决方案。这项工作有望推动轨道优化激发态方法在光化学、材料设计、光谱学模拟等领域的更广泛和更可靠的应用，标志着变分激发态计算领域的一个重要进展。