Gaussian计算：有机小分子激发态与吸收光谱的理论预测

实验上测一个吸收光谱，快的半小时就出来了。但用Gaussian算一个靠谱的激发态——从选泛函到选基组到加溶剂效应，每一步都有坑。这个项目在预测香豆素衍生物的吸收光谱时，目标很明确：算出的λ_max与实验值偏差不超过20 nm。结果第一轮差了快50 nm，逼着从头审了一遍参数链。

泛函选择了TD-DFT的精度上限

对有机小分子的激发态，TD-DFT是平衡精度和计算量的主流选择。但泛函怎么选？项目测试了六种泛函：B3LYP、CAM-B3LYP、ωB97XD、M06-2X、PBE0、M11。

B3LYP偏差最大是意料之中——它对电荷转移激发描述差，这是老问题了。荷移激发在香豆素体系里占比约40%，B3LYP的HF成分(20%)不够压制自相互作用误差。M06-2X和M11表现最好，偏差都在3 nm以内。最终项目选M06-2X——不是偏差最小的，而是它对这个体系的振子强度预测比M11更接近实验的相对强度分布。

基组：三重ζ就够，但加弥散是关键

6-31+G(d,p)和6-311++G(d,p)给出了几乎一致的λ_max（差不到1 nm），说明极化函数从二重ζ到三重ζ对这个体系不是瓶颈。真正影响大的是弥散函数——6-31G(d,p)（无弥散）算出的λ_max是371 nm，偏离31 nm。香豆素衍生物的激发态涉及π→π*跃迁，激发态电子密度比基态更弥散，弥散函数对激发态轨道描述不是锦上添花，是必需品。

溶剂效应的两步走

PCM隐式溶剂模型是Gaussian中的标配。但项目发现对香豆素这个体系，PCM在基态优化阶段开不开影响不大（键长差<0.003 Å），但在TD-DFT激发态计算阶段必须开，否则λ_max偏15-20 nm。原因是溶剂极性通过反应场稳定了激发态的偶极矩，改变了HOMO-LUMO gap。

最终选定的计算方案是：ωB97XD/6-31+G(d,p)做基态优化 → M06-2X/6-31+G(d,p)+PCM(溶剂=乙醇)做TD-DFT单点。λ_max=343 nm，与实验340 nm偏差3 nm，在这个精度下可以直接指导染料分子设计，减少合成和测试的试错成本。

参考文献：Jacquemin et al., J. Chem. Theory Comput., 2009, 5, 2420-2430；Laurent & Jacquemin, Int. J. Quantum Chem., 2013, 113, 2019-2039.