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

1. 单重态 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矩阵元时，通常要求单重态和三重态数量完全一样。

2. 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等全部信息！

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

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激发态。

4. 输出文件里怎么看结果？

输出文件里会先列出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₁)。

5. 常见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），对大部分有机体系误差很小。