Gaussian中溶剂模型的设置及计算

Gaussian主要支持隐式溶剂模型，把溶剂看成连续介质（介电常数ε），溶质放在一个空腔里。

主流模型对比：

IEFPCM / PCM（默认）：积分方程形式，腔体由重叠球构成，适合几何优化、频率计算、激发态。默认使用UFF半径×1.1。

SMD（强烈推荐）：Truhlar小组开发，包含非静电项（色散、氢键、空化），在多数体系中对ΔG预测最可靠。

CPCM：导体型变体，通常与IEFPCM精度接近，但数值略有系统差异。

Onsager：球形腔，仅适合简单情况。

什么时候用哪个？

想算ΔG溶剂化能 → 优先SMD

做优化+频率 → IEFPCM或SMD

激发态 →

👉 TD + IEFPCM + NonEq（非平衡溶剂）更标准

基本设置方法

在路由段加 SCRF 关键字即可

text

# B3LYP/6-31G(d) Opt Freq SCRF=(SMD,Solvent=Water)

常见溶剂名称（直接用英文，Gaussian内置100+种）：

Water（水，ε=78.3553）

Methanol（甲醇）

Ethanol（乙醇）

Acetonitrile（乙腈）

DMSO（二甲基亚砜）

Chloroform（氯仿）

Benzene（苯）

Acetone（丙酮）

Dichloromethane（二氯甲烷）

完整列表可去gaussian.com/scrf查（ε值会自动匹配）。

完整输入文件示例（水溶液中优化乙醇分子）：

text

%chk=ethanol.chk

# B3LYP/6-31G(d) Opt Freq SCRF=(SMD,Solvent=Water)

Ethanol in water

0 1

C  0.000000  0.000000  0.000000

C  0.000000  0.000000  1.500000

O  0.000000  1.200000  2.100000

H  0.900000  0.000000  -0.500000

...（完整坐标）

（末尾空行即可）

运行后，输出会显示“Solvent: Water” 和溶剂化自由能信息。

不同计算类型的设置技巧

几何优化 + 频率（最常用）

text

# Opt Freq SCRF=(IEFPCM,Solvent=Ethanol)

或者用SMD更准。

单点能量（算溶剂化能ΔG）

先跑气相：# B3LYP/6-311++G(d,p) SCRF=No

再跑溶液：# B3LYP/6-311++G(d,p) SCRF=(SMD,Solvent=Water)

a）ΔGsolv ≈ G_solution - G_gas（SMD直接给出G值）b）直接用SMD输出Free Energy of Solvation

过渡态搜索

text

# Opt=(TS,CalcFC) Freq SCRF=(SMD,Solvent=Water)

读取旧检查点继续计算（节省时间）

text

# Geom=AllCheck SCRF=(SMD,Solvent=Water,Read)

自定义溶剂

内置没有的溶剂？用 SCRF=(SMD,Solvent=Generic,Read) + 额外输入段。

SMD需要7个参数，部分参数可省（Gaussian会补默认值）：

text

Eps=          （静态介电常数）

EpsInf=       （光学折射率平方 n²）

HbondAcidity= （氢键酸性 α）

HbondBasicity=（氢键碱性 β）

SurfaceTensionAtInterface= （表面张力参数）

CarbonAromaticity= （芳香碳比例）

ElectronegativeHalogenicity= （电负性卤素比例）

示例（自定义溴苯）：

text

# B3LYP/6-31G(d) SCRF=(SMD,Solvent=Generic,Read)

...分子坐标...

Eps=5.3954

EpsInf=2.432664

HbondAcidity=0.00

HbondBasicity=0.09

SurfaceTensionAtInterface=50.72

CarbonAromaticity=0.857

ElectronegativeHalogenicity=0.143

（空行结束）

普通PCM只需Eps和EpsInf即可。

结果解读

输出中搜索 “Free Energy of Solvation” 或 “ΔG” 就是溶剂化自由能。

SMD会额外给出非静电贡献（Dispersion, Cavitation等）。

常见错误：

忘了加SCRF → 还是气相

溶剂名称拼错 → Gaussian会报错

高精度计算用6-311++G(d,p) + SMD

不同模型不能直接比较ΔG，优化和单点模型要一致

小Tips：

混合溶剂：目前不支持，可用参数线性加权（参考文献）。

计算资源：溶剂模型比气相慢20-50%，建议用Def2-TZVP基组平衡精度。

Gaussian 09 vs 16：写法完全一样，16支持更多外部迭代。

实战案例：计算苯胺在水中的溶剂化能

气相单点 → 得到G_gas

水溶液SMD单点 → 得到G_water

ΔGsolv = G_water - G_gas（实验值约-5.5 kcal/mol，可对比验证）

总结

Gaussian中溶剂模型的使用非常简洁，通过一行SCRF关键字即可实现大多数计算需求。在众多模型中，SMD因同时考虑电静与非电静贡献，已成为当前计算溶剂化自由能的首选方案。对于结构优化与反应路径研究，IEFPCM与SMD均可胜任。合理选择模型并结合计算目标，是获得可靠溶剂效应结果的关键。