Gaussian計算溶劑環(huán)境溶質自由能模板
分布計算需要設置四個文件,以下是這四個文件的關鍵詞內容:
(1)B3LYP/6-311G* em=GD3BJ opt freq scrf(SMD,solvent=ethanol)
對應SMD溶劑模型表現乙醇環(huán)境的情況下用B3LYP-D3(BJ)/6-311G*對體系進行優(yōu)化和振動分析。這個級別不昂貴,而用于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溶劑模型表現乙醇環(huán)境的情況下計算電子能量
溶質溶解自由能:(單位: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