激发电子态变分计算的冻结释放直接优化方法

《Journal of Chemical Theory and Computation》：Freeze-and-Release Direct Optimization Method for Variational Calculations of Excited Electronic States

　　

激发电子态的含时不依赖密度泛函计算中的轨道的变分优化面临重大挑战，因为激发态通常对应于电子能量曲面上的鞍点。如果激发涉及电子密度的显著重排（如电荷转移激发），优化可能尤为困难。本研究提出了一种用于激发态变分轨道优化的简单策略。该方法通过冻结直接参与激发的轨道的同时最小化能量，然后进行完全无约束的鞍点优化。该冻结释放策略的两个步骤均使用与基态计算相同计算标度的直接优化算法进行。研究人员利用广义梯度近似泛函，在有机分子和分子二聚体的分子内与分子间电荷转移激发态计算中广泛评估了该方法的性能。研究发现，对于基于最大重叠法的传统算法失效的情况（即收敛到虚假的电荷离域解），冻结释放直接优化方法能够避免变分塌缩。对于分子间电荷转移，轨道优化计算提供了能量对给体-受体分离距离的正确依赖性，而无需长程精确交换，这是常见含时密度泛函理论方法未能实现的。

# 冻结释放直接优化法在激发电子态变分计算中的研究解读

## 研究背景与问题

准确模拟大分子发色团和光活性材料中的电子激发对于理解光合作用、视觉等自然现象以及开发高效可持续的太阳能转换器件至关重要。光诱导过程通常涉及多个具有不同电子特性的态，其中电荷转移激发态常起关键作用。然而，当前激发态计算方法在可负担的计算成本与可靠描述不同电子特性态的能力之间存在矛盾。传统方法（如含时密度泛函理论，TDDFT）在处理电荷转移激发时存在显著缺陷：低估激发能，且无法正确给出能量随给体-受体距离的1/R变化规律。这些失败源于标准Kohn-Sham（KS）泛函缺乏导数不连续性以及绝热交换相关（XC）核在无穷远距离的发散问题。尽管范围分离杂化泛函和非绝热近似有所改进，但系统特异性或发展阶段限制使其难以广泛应用。

轨道优化（OO）方法（亦称Δ自洽场，ΔSCF）通过求解非Aufbau占据的KS方程获得激发态，具有与基态等价的变分处理优势，在双激发、核激发和电荷转移激发中优于TDDFT。然而，KS激发态解通常对应电子能量曲面上的鞍点（saddle point），优化过程中易沿负曲率方向发生变分塌缩，收敛到比目标态能量更低的解。传统策略如最大重叠法（MOM）虽能缓解塌缩，但无法完全消除，尤其在电荷转移激发中常导致振荡收敛或塌缩至电荷离域化虚假解。因此，亟需一种既稳健又易于实现的激发态轨道优化策略。

## 研究内容与结论

本研究提出并系统评估了一种称为“冻结释放直接优化”（FR-DO）的两步策略。该方法首先在约束优化中冻结直接参与激发的轨道（即空穴和激发电子轨道），仅松弛其余轨道，然后将部分松弛的轨道作为初猜，进行完全无约束的鞍点优化。研究人员采用PBE泛函（广义梯度近似，GGA），对包含15个有机分子27个分子内电荷转移态的数据集以及四氟乙烯-乙烯和氨-氟两个分子二聚体在不同分子间距离下的分子间电荷转移态进行了广泛测试。结果表明：FR-DO在所有情况下均成功收敛到目标电荷转移解，无变分塌缩；而传统MOM与直接优化（DO-MOM）或SCF-MOM（基于对角化的自洽场MOM）结合的方法存在塌缩或收敛失败。FR-DO虽平均迭代次数略高于DO-MOM，但避免了塌缩，且计算标度与基态计算相同。更重要的是，对于分子间电荷转移态，FR-DO结合GGA泛函正确预测了激发能随给体-受体距离的近似1/R依赖关系，与高精度EOM-CCSD(T)（运动方程耦合簇）参考结果吻合，而线性响应TDDFT在同一泛函下给出了定性错误的能量曲线。该研究发表在《Journal of Chemical Theory and Computation》。

## 关键技术方法

（不超过250字）主要关键技术方法包括：（1）直接轨道优化（Direct Optimization, DO），通过指数变换参数化酉变换矩阵，利用准牛顿算法（L-BFGS和L-SR1）直接优化轨道使KS能量驻定。（2）约束优化，在优化中冻结参与激发的轨道（空穴和激发电子轨道），仅松弛其余轨道，防止变分塌缩。（3）释放约束后的无约束优化，采用对称秩1（L-SR1）拟牛顿算法，处理含负特征值的Hessian矩阵，收敛到鞍点。（4）通过数值对角化电子Hessian矩阵诊断鞍点阶数（负特征值个数）。（5）使用最大重叠法（MOM）作为对照，基于轨道投影权重分配占据数。所有计算使用PBE泛函和PAW赝势，基组为aug-cc-pVDZ+sz（分子内）或cc-pVDZ+sz（分子间）。分子几何来自文献：分子内数据集源于Loos等（CCSD(T)/CC3优化），二聚体几何分别来自CC2和MC-QCISD/3优化。

## 研究结果

### 3.1 分子内电荷转移态

研究人员测试了27个分子内电荷转移态，比较FR-DO、DO-MOM和SCF-MOM的收敛性。SCF-MOM无变分塌缩但2例不收敛；DO-MOM全部收敛但4例塌缩；FR-DO全部收敛无塌缩，平均迭代次数（17.6）略高于DO-MOM（11.7），但低于SCF-MOM（19.2）。FR-DO的改进在于约束优化后，电子Hessian的负特征值数量更接近目标解的真实鞍点阶数（平均偏差0.7），为后续无约束步提供了更准确的负曲率方向。

#### 3.1.1 估计的鞍点阶数
通过比较初猜（基态轨道更改占据）和约束优化后的鞍点阶数发现：初猜时，数值Hessian和解析对角近似分别高估和低估负特征值数量（尤其大电荷转移距离态）；约束优化后两者趋于一致，仅需排除小负特征值即可精确估计目标解阶数。

#### 3.1.2 扭曲的N-苯基吡咯
以扭曲N-苯基吡咯（PP）的A1电荷转移态为例：MOM-DO塌缩到电荷离域解（激发能4.61 eV，电荷转移距离2.06 ?），而FR-DO收敛到电荷局域解（5.56 eV，2.41 ?），接近理论最佳估计（5.65 eV）和CAM-B3LYP TDDFT结果。按轨道旋转角κπphπpy扫描能量曲线，初猜时准牛顿模型给出正曲率，导致DO-MOM步走向最小；约束稳定局域轨道后，模型给出负曲率，释放约束后正确走向鞍点。MOM在旋转角约45°时无法区分，因重叠投影接近0.7。

### 3.2 分子间电荷转移态

对于四氟乙烯-乙烯和氨-氟二聚体，FR-DO在所有距离下收敛到电荷局域解，而DO-MOM在短距离（3.5 ?）严重塌缩（能量低>3 eV），在部分距离不收敛。FR-DO解的能量与EOM-CCSD(T)吻合良好（最大偏差约0.3 eV），且激发能随距离变化呈现正确的1/R趋势，而TDDFT曲线定性错误。电荷转移距离也与参考值一致。

## 讨论与结论总结

讨论部分指出：电荷转移激发的电子能量曲面在基态轨道初猜下高度崎岖，电子Hessian有大量负特征值，远高于目标解鞍点阶数，导致对角化SCF算法收敛困难。直接优化方法虽更稳健，但解析对角近似严重低估负特征值数量，使计算易塌缩至低阶鞍点的电荷离域虚假解。MOM无法防止此类塌缩，因离域涉及混合角度<45°。FR-DO通过约束优化得到部分松弛的初猜，其Hessian负特征值接近目标解阶数，从而引导无约束优化正确收敛至鞍点。

研究结论（翻译原文Section 4）：本工作揭示了在轨道优化密度泛函计算中使用GGA泛函处理电荷转移态时面临的一些挑战。对于有机分子中大量电荷转移激发，基态轨道非Aufbau占据初猜下的电子Hessian具有大量负特征值，远高于目标解鞍点阶数，表明能量曲面高度崎岖。基于直接轨道优化的方法通常比基于哈密顿矩阵特征分解的SCF方法更稳健，但常用的初始Hessian对角近似无法捕捉电子能量曲面的结构，显著低估负特征值数量，导致计算塌缩至低鞍点阶数的电荷离域解。最大重叠法无法防止此类塌缩。本研究提出的冻结释放直接优化（FR-DO）方法通过两步策略稳健收敛PBE计算中具有挑战性的分子内和分子间电荷转移态。对于分子间电荷转移，PBE泛函的FR-DO计算结果给出激发能随给体-受体距离近似1/R的变化关系，与高精度耦合簇计算一致，而线性响应TDDFT给出定性错误的能量曲线。该方法为高效、稳健的激发态变分计算提供了有效路径。