Gaussian计算溶剂环境溶质自由能模板

分布计算需要设置四个文件，以下是这四个文件的关键词内容：

(1)B3LYP/6-311G* em=GD3BJ opt freq scrf(SMD,solvent=ethanol)

对应SMD溶剂模型表现乙醇环境的情况下用B3LYP-D3(BJ)/6-311G*对体系进行优化和振动分析。这个级别不昂贵，而用于opt freq目的足够给出可靠的结果。

(2)B2PLYPD3/def2TZVP geom=allcheck，

使用较好的B2PLYP-D3(BJ)/def2-TZVP级别基于(1)的结构计算真空下较高精度电子能量

(3)M052X/6-31G* geom=allcheck，

基于(1)的结构用M05-2X/6-31G*在真空下计算电子能量

(4)M052X/6-31G* scrf(SMD,solvent=ethanol) geom=allcheck，

基于(1)的结构用M05-2X/6-31G*在SMD溶剂模型表现乙醇环境的情况下计算电子能量

溶质溶解自由能：（单位：Hartree）

ΔE = E_[in_4.out] - E_[in_3.out]

溶质在溶剂中的自由能：（单位：Hartree）

G_corr_[in_1.out] + E_[in_2.out] + ΔE + 1.89/627.51

直接在一个文件中设置多步计算。

【多步计算示例】：

%nprocshared=20

%mem=80GB

%chk=SL-1.chk

#p opt freq B3LYP/6-311G* em=GD3BJ scrf=(iefpcm,solvent=generic,read) pressure=0.3

temperature=443

Title Card Required

0 1

 #坐标

eps=39.9

epsinf=1.9

---link1----

%oldchk=SL-1.chk

%nprocshared=20

%mem=80GB

%chk=SL-1_2.chk

#p B2PLYPD3/def2TZVP geom=allcheck

Title Card Required

0 1

--link1---

%oldchk=SL-1.chk

%nprocshared=20

%mem=80GB

%chk=SL-1_3.chk

#p M052X/6-31G* geom=allcheck

Title Card Required

0 1

---link1---

%oldchk=SL-1.chk

%nprocshared=20

%mem=80GB

%chk=SL-1_4.chk

#p M052X/6-31G* scrf=(iefpcm,solvent=generic,read) geom=allcheck

Title Card Required

0 1

eps=39.9

epsinf=1.9