Gaussian同时计算多个单重态和三重态

在量子化学尤其是激发态计算里单重态Singlet和三重态Triplet是两个最核心的概念很多OLED/TADF材料、磷光材料、光催化等研究都绕不开它们。下面我直接给你最实用的干货重点讲清楚Gaussian里如何用一个任务同时算出3个单重态和3个三重态。1. 单重态 vs 三重态单重态SSinglet总自旋S0多重度2S11。两个电子自旋相反↑↓总自旋为0。最常见的基态S₀就是闭壳层单重态。激发态里S₁、S₂…也是单重态。三重态TTriplet总自旋S1多重度2S13。两个电子自旋平行↑↑总自旋为1有3个简并分量Tₓ、Tᵧ、T_z。磷光、系间窜越ISC都发生在三重态。为什么重要TADF材料需要S₁和T₁能量很接近∆E_ST 0.1 eV才能通过热激活实现反向系间窜越RISC。磷光材料需要强的旋轨耦合SOC把T₁ → S₀的跃迁偷出来。计算SOC矩阵元时通常要求单重态和三重态数量完全一样。2. Gaussian里怎么算激发态TD-DFT/CIS基础Gaussian用TDTime-Dependent关键词做激发态计算支持TD-DFT、TD-HF、CIS等。 默认只算单重态Singlets因为大多数有机分子基态是闭壳层单重态Charge0, Multiplicity1。关键选项NStatesN总共算N个激发态默认N3。Singlets只算单重态默认闭壳层。Triplets只算三重态。50-50同时算N个单重态 N个三重态闭壳层基态专用如果指定50-50NStates的值就是每种态的数量默认就是3个单重态3个三重态。最常用写法推荐text# TD(50-50,NStates3) 泛函/基组→ 一次性得到 S₁~S₃ 和 T₁~T₃ 的激发能、振子强度、跃迁偶极矩、NTO等全部信息3. 完整实用输入文件示例直接复制改分子就行gaussian%mem16GB%nproc16%chktd_50-50.chk# wB97XD/def2-TZVP TD(50-50,NStates5) 6D 10F nosymm GFInputTitle: TD-DFT 同时算5个S 5个T推荐多算2个防态混叠0 1...你的分子坐标小贴士为什么NStates5而不是刚好3因为TD-DFT是按能量从低到高排序的多算1~2个能避免“感兴趣的态被挤到后面”。wB97XD或CAM-B3LYP这类长程校正泛函对电荷转移激发态和∆E_ST比较准。6D 10F nosymm GFInput 是PySOC等后续脚本的硬性要求。基态必须是0 1闭壳层单重态不要改成0 3改成0 3就是直接算三重态基态UDFT不是TD激发态。4. 输出文件里怎么看结果输出文件里会先列出Singlet Excited StatesS₁~S₅然后是Triplet Excited StatesT₁~T₅。典型输出片段textExcited State 1: Singlet-A 3.1234 eV ... f0.85...Excited State 6: Triplet-A 2.9876 eV ...S₁ 就是第一个Singlet Excited State。T₁ 就是第一个Triplet Excited State。直接用grep Excited State或GaussView看就能一目了然。如果要做TADF∆E_ST E(S₁) - E(T₁)。5. 常见tips只适用于闭壳层基态。开壳层基态双自由基等不能用50-50只能靠S²自己判断是单重还是三重。计算量 ≈ 单重态计算的2倍因为要解两个独立的响应方程。想优化S₁/T₁结构加opt root1或root6对应T₁。想算SOC旋轨耦合必须用TD(50-50,NStatesxx) 保留.rwf文件再跑PySOC。TDA近似更快TDA(50-50,NStates5)Tamm-Dancoff Approximation对大部分有机体系误差很小。