Gaussian中激发态的结构优化

激发态结构优化是指在量子化学计算中，对分子在电子激发态（如S1、T1等）的几何结构进行优化，找到其最低能量构型。这不同于垂直激发（固定基态几何），优化后可计算绝热激发能（Adiabatic Excitation Energy，基态最小点到激发态最小点的能量差）、荧光发射能（从S1优化几何到S0的垂直能量差）或磷光等。在光化学、OLED材料、光伏设计等领域至关重要，能模拟实际光谱峰位（考虑核弛豫）。

Gaussian软件支持激发态优化，主要用TD-DFT（时间相关密度泛函理论），效率高，适用于中等分子。下面一步步详细解释操作、输入输出和技巧，做一个干货分享。

步骤1: 准备初始结构

从基态优化后的几何开始，因为激发态往往从基态跃迁而来。

选择计算水平：泛函如CAM-B3LYP（长程校正，避免CT激发问题）、ωB97XD；基组6-31+G(d,p)或def2-TZVP。TD-DFT是主力，避免CIS（太不准）。

溶剂：用SCRF=PCM模拟环境，激发态常在溶液中。

干货：激发态优化计算量大（~基态的2-5倍），小分子测试用B3LYP/6-31G(d)，大分子加TDA加速。

步骤2: 激发态几何优化

关键词opt td(root=N)指定优化第N个激发态（root=1为S1）。Gaussian会用TD-DFT计算激发态梯度，进行优化。

输入文件示例（名为gjf，从基态chk读取）：

text

%chk=excited.chk

%oldchk=ground.chk  # 从基态检查点读取初始几何

# opt td(root=1, nstates=5, singlets) cam-b3lyp/6-31+g(d,p) geom=checkpoint  # opt=优化，td(root=1)=优化S1，nstates=5确保根稳定

Excited State S1 Optimization  # 标题

0 1  # 电荷0，自旋1（闭壳层单重态）

运行：g16 excited_opt.gjf。

关键选项：

root=N：指定哪个激发态（如1为最低单重态S1，2为S2）。

nstates=M：计算前M个激发态（M>root，确保能量排序正确，通常M=root+3~5）。

singlets/triplets：单重态/三重态（磷光用三重态T1，root=1）。

TDA：加此加速，牺牲少许精度。

opt=calcfc：初次优化用，计算力常数加速收敛。

开壳层：如三重态，用自旋多重度3，# uCAM-B3LYP。

步骤3: 频率计算验证

优化后，加freq确认是真实最小点（无虚频）。

输入示例（gjf）：

text

%chk=freq.chk

%oldchk=excited.chk

# freq td(root=1, nstates=5) cam-b3lyp/6-31+g(d,p) geom=checkpoint

Excited State Frequency Calculation

0 1

输出检查：搜索"No imaginary frequencies"，虚频表示鞍点，重优化。

步骤4: 输出解读与应用

优化输出类似基态，但有TD部分。

示例片段：

text

Optimization completed.

-- Stationary point found.

Excited State   1:      Singlet-A    3.4567 eV  358.90 nm  f=0.0789

Total Energy, E(TD-HF/TD-DFT) =  -123.456789 Hartree

优化几何：.log末尾提取新坐标。

激发能：优化后TD能量是激发态最小点能量。

绝热激发能：ΔE_adiabatic = E_excited_opt - E_ground_opt（单位Hartree，转eV: ×27.211）。

荧光发射：在S1优化几何上做TD单点（root=1），但计算到S0的垂直能量差（即S1能量 - S0能量，在S1几何上）。

计算公式：荧光能 = E_S1(at S1 geom) - E_S0(at S1 geom)（垂直发射）。

注意事项与优化技巧

根翻转问题：常见于类似能量激发态。解决：增大nstates，监控每步TD输出； CIS先粗优，再TD精优等。

准确性：TD-DFT对锥形交叉（Conical Intersection）不准，用CASSCF/MCSCF高级方法（Gaussian支持# cas）。

计算资源：%mem=16GB, %nproc=16。失败时检查SCF收敛（加scf=maxcycle=200）。

扩展应用：

磷光：优化T1（td(triplets, root=1)），计算磷光发射。

非辐射弛豫：优化后比对基态/激发态RMSD，分析结构变化。

光谱模拟：用优化几何计算振子强度，生成发射谱。

常见错误：

"TDDFT failed to converge"：换泛函或基组。

负激发能：表示激发态低于基态，模型失效，重选方法。

总结干货

激发态优化是“TD-DFT + opt + root指定”的核心组合，从基态几何起步，输出优化结构和能量，用于绝热/发射计算。入门用CAM-B3LYP/6-31+G(d,p)，注意根稳定和频率验证。比垂直激发复杂，但更贴近实验（如荧光波长红移）。多实践小分子如苯、荧光染料，就能上手。如果有具体输入或问题，随时补充！