垂直激发能（Vertical Excitation Energy，简称VEE）是指分子在基态（Ground State，通常记为S0）的几何结构下，从基态跃迁到激发态（Excited State，如S1、S2等）的能量差。根据Franck-Condon原理，电子跃迁发生在极短时间内（~10^-15秒），分子核坐标几乎不变，因此计算时不需优化激发态几何，而是直接在基态几何上计算能量差。这在光谱学、荧光、光伏材料等领域非常关键，能模拟UV-Vis吸收光谱，帮助理解分子光物理行为。
Gaussian软件（Gaussian 09/16等版本）是量子化学计算的常用工具，支持多种方法计算垂直激发能，最常用的是时间相关密度泛函理论（TD-DFT）。它计算效率高，适用于中等大小分子。下面我们来详细解释如何用Gaussian计算垂直激发能，结合实际操作、输入示例和输出解读，做一个纯干货分享。假设你有Gaussian安装和基本命令行知识。
步骤1: 准备分子结构
- 获取初始结构：用ChemDraw、Avogadro等软件绘制分子，导出为Gaussian输入文件（.gjf格式）。坐标可以用Cartesian（x,y,z）或Z-matrix。
- 选择计算水平：
- 泛函：B3LYP、CAM-B3LYP（长程校正，更适合激发态）、ωB97XD等。TD-DFT常用B3LYP作为入门。
- 基组：6-31G(d)、6-311+G(d,p)等。带扩散函数（如+）对激发态更好，但计算量大。
- 溶剂效应：如果分子在溶液中，用PCM模型（如SCRF=PCM, Solvent=Water）模拟。
提示：小分子用6-31G(d)测试，大分子考虑def2-SVP。避免HF方法，它对激发态不准。
步骤2: 基态几何优化
垂直激发基于基态几何，先优化S0结构。
- 输入文件示例（名为ground_opt.gjf，以甲烷为例，实际替换你的分子）：
text
%chk=ground.chk # 保存检查点文件，便于后续使用
# opt b3lyp/6-31g(d) # opt=优化，b3lyp/6-31g(d)=计算水平
Ground State Optimization # 标题
0 1 # 电荷0，自旋多重度1（闭壳层单重态）
C 0.000000 0.000000 0.000000
H 0.629941 0.629941 0.629941
H -0.629941 -0.629941 0.629941
H -0.629941 0.629941 -0.629941
H 0.629941 -0.629941 -0.629941
- 运行命令：g16 ground_opt.gjf（或g09，根据版本）。
- 输出检查：看"Optimization completed"，提取优化后的坐标（在.log文件末尾）。
步骤3: 计算垂直激发能（TD计算）
在优化后的基态几何上，进行TD计算。TD关键词请求时间相关计算，输出激发态能量、振子强度（Oscillator Strength，f，表示跃迁概率）、轨道贡献等。
- 输入文件示例（名为td_calc.gjf）：
text
%chk=td.chk
%oldchk=ground.chk # 从基态检查点读取几何
# td(nstates=10, root=1) b3lyp/6-31g(d) geom=checkpoint # td=TD-DFT，nstates=计算前10个激发态，root=1指定第一个激发态（S1）
Vertical Excitation Energy Calculation
0 1 # 同基态
- 运行：g16 td_calc.gjf。
- 关键选项：
- nstates=N：计算前N个激发态（默认6，建议10-20，避免遗漏）。
- singlets：只算单重态（默认），triplets算三重态。
- 50-50：混合HF和DFT，改善电荷转移激发。
- TDA：Tamm-Dancoff近似，加速计算，适合大分子。
- 溶剂：添加scrf=(pcm,solvent=ethanol)模拟环境。
高级干货：如果用CIS，替换td为cis，但CIS不如TD-DFT准，只适合小分子测试。对于开壳层系统（如自由基），用UDFT，自旋多重度调整为3等。
步骤4: 输出解读
TD计算输出在.log文件，搜索"Excited State"部分。
示例输出片段：
text
Excited State 1: Singlet-A 4.1234 eV 300.78 nm f=0.0456 <S**2>=0.000
29 -> 30 0.70123 # 从HOMO (29) 到 LUMO (30) 的贡献
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-DFT) = -40.123456789
Excited State 2: Singlet-A 5.6789 eV 218.34 nm f=0.1234 <S**2>=0.000
28 -> 30 0.56789
- 垂直激发能：以eV或nm单位给出。4.1234 eV就是S0到S1的垂直激发能。
- nm：波长，模拟吸收峰位置。
- f：振子强度，f>0.01表示允许跃迁。
- 轨道贡献：如HOMO→LUMO，帮你理解跃迁类型（π-π、n-π等）。
- 总能量：E(TD-HF/TD-DFT)是激发态总能，基态能量从优化输出取，差值即激发能（单位Hartree，转换：1 Hartree ≈ 27.211 eV）。
计算公式：垂直激发能 ΔE = E_excited - E_ground（在相同几何上）。
注意事项与优化技巧
- 准确性：TD-DFT对电荷转移（CT）激发低估，用长程校正泛函如CAM-B3LYP。验证时比较实验光谱。
- 计算资源：nstates大时内存高，用%mem=8GB。并行：%nproc=8。
- 错误排查：
- "Convergence failure"：增加迭代次数（iop(9/40=4)）。
- 激发能负值：基组太小或泛函不合适。
- 扩展应用：
- 荧光发射：优化S1几何（# opt td(root=1)），然后计算垂直发射能。
- 光谱模拟：用GaussView可视化输出，生成UV-Vis谱图。
- 大分子：用ONIOM分层，或TDA近似减时。
- 单位转换：eV to kcal/mol: ×23.06；nm to eV: 1240/nm。
- 参考文献：Gaussian官网TD关键词文档；论文如Casida (1998) TD-DFT理论。
总结
使用Gaussian计算垂直激发能（Vertical Excitation Energy）的标准流程如下：
- 先进行基态几何优化（# opt B3LYP/6-31G(d) 或 CAM-B3LYP 等，建议加PCM溶剂模型）。
- 在优化后的基态几何上做TD-DFT单点计算（# td(nstates=10) B3LYP/6-31G(d) geom=checkpoint）。
- 从输出中提取“Excited State”部分的能量（eV或nm），即为从基态到对应激发态的垂直激发能，同时查看振子强度f和主要轨道跃迁（如HOMO→LUMO）。
这套“基态优化 + TD单点”两步法基于Franck-Condon原理，计算成本适中，最接近实验吸收光谱峰位。精度提升可换长程校正泛函（如CAM-B3LYP、ωB97XD），避免CT激发低估。TD输出直接用于高斯展宽即可生成模拟UV-Vis谱图。
