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

在量子化学尤其是激发态计算里，单重态（Singlet）和三重态（Triplet）是两个最核心的概念，很多OLED/TADF材料、磷光材料、光催化等研究都绕不开它们。下面我直接给你最实用的干货，重点讲清楚Gaussian里如何用一个任务同时算出3个单重态和3个三重态。

单重态 vs 三重态

单重态（S，Singlet）：总自旋S=0，多重度2S+1=1。两个电子自旋相反（↑↓），总自旋为0。最常见的基态（S₀）就是闭壳层单重态。激发态里S₁、S₂…也是单重态。

三重态（T，Triplet）：总自旋S=1，多重度2S+1=3。两个电子自旋平行（↑↑），总自旋为1（有3个简并分量：Tₓ、Tᵧ、T_z）。磷光、系间窜越（ISC）都发生在三重态。

为什么重要？

TADF材料需要S₁和T₁能量很接近（∆E_ST < 0.1 eV），才能通过热激活实现反向系间窜越（RISC）。

磷光材料需要强的旋轨耦合（SOC），把T₁ → S₀的跃迁偷出来。

计算SOC矩阵元时，通常要求单重态和三重态数量完全一样。

Gaussian里怎么算激发态？（TD-DFT/CIS基础）

Gaussian用TD（Time-Dependent）关键词做激发态计算（支持TD-DFT、TD-HF、CIS等）。默认只算单重态（Singlets），因为大多数有机分子基态是闭壳层单重态（Charge=0, Multiplicity=1）。

关键选项：

NStates=N：总共算N个激发态（默认N=3）。

Singlets：只算单重态（默认，闭壳层）。

Triplets：只算三重态。

50-50：同时算N个单重态 + N个三重态（闭壳层基态专用）！如果指定50-50，NStates的值就是每种态的数量（默认就是3个单重态+3个三重态）。

最常用写法（推荐）：

text

# TD(50-50,NStates=3) 泛函/基组

→ 一次性得到 S₁~S₃ 和 T₁~T₃ 的激发能、振子强度、跃迁偶极矩、NTO等全部信息！

完整实用输入文件示例（直接复制改分子就行）

gaussian

%mem=16GB

%nproc=16

%chk=td_50-50.chk

# wB97XD/def2-TZVP TD(50-50,NStates=5) 6D 10F nosymm GFInput

Title: TD-DFT 同时算5个S + 5个T（推荐多算2个防态混叠）

0 1

...（你的分子坐标）

小贴士：

为什么NStates=5而不是刚好3？因为TD-DFT是按能量从低到高排序的，多算1~2个能避免“感兴趣的态被挤到后面”。

wB97XD或CAM-B3LYP这类长程校正泛函对电荷转移激发态和∆E_ST比较准。

6D 10F + nosymm + GFInput 是PySOC等后续脚本的硬性要求。

基态必须是0 1（闭壳层单重态），不要改成0 3！改成0 3就是直接算三重态基态UDFT，不是TD激发态。

输出文件里怎么看结果？

输出文件里会先列出Singlet Excited States（S₁~S₅），然后是Triplet Excited States（T₁~T₅）。

典型输出片段：

text

Excited State 1: Singlet-A 3.1234 eV ... f=0.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₁)。

常见tips

只适用于闭壳层基态。开壳层基态（双自由基等）不能用50-50，只能靠<S²>自己判断是单重还是三重。

计算量 ≈ 单重态计算的2倍（因为要解两个独立的响应方程）。

想优化S₁/T₁结构？加opt + root=1（或root=6对应T₁）。

想算SOC（旋轨耦合）？必须用TD(50-50,NStates=xx) + 保留.rwf文件，再跑PySOC。

TDA近似更快：TDA(50-50,NStates=5)（Tamm-Dancoff Approximation），对大部分有机体系误差很小。