发布时间 : 星期三 文章gaussian更新完毕开始阅读
AEA(adiabatic electron affinity)=E(A)-E(A-)
VEA（vertical electron affinity）=E(A) –E(A- at A geometry) VDE （vertical detachment energy）=E(A+ at A geometry) –E（A） 优化关键词
# opt b3lyp/gen iop(2/15=1,5/13=1) pseudo=read (分子说明行后空一格)。基组产生网站：https://bse.pnl.gov/bse/portal
在优化的基础上保存chk，再频率计算
# freq=noram scrf=check guess=read geom=allcheck chkbas genchk test
一个文件指定频率和优化 %chk=a
# method freq ** **
--link1— %chk=a %nosave
#method geom=check guess=read freq=readfc
Frequence at T
0 1
400.0 2.0
混合基组输入
# opt hf/genecp iop（2/15=1,5/13=1）optcycle=100 ***
C H 0 6-31G(D) **** Pt 0 SDD **** Pt 0
SDD
NBO 计算
# b3lyp/6-311+g** po=nboread
$nbo bndidx $end
BOMD计算
由于计算BOMD需要指定分裂路径，如果不指定的话，需要在分子说明完后空两行。 如果自定义输入基组的话，建议这样输
分解路径为0，前空一行，后空1行
关键词
MaxStep=n步长n*0.0001fs
Maxpoints=400 每个轨迹中包含的点数
Stepsize=10000 Sets the dynamic step size to n*0.0001, in the appropriate units. The default step size is 0.25 amu1/2*Bohr except for
NTraj=N 计算N步轨迹
组合ntraj maxpoints
Ntraj和maxstep不能组合
过渡态计算 先小基组qst2或qst计算出过渡态，然后高精度qst确认即可。 可以先TS2，
#p opt=(calcfc,qst2,noeigen，maxcycle=N) 方法基组
也可以初猜一个进行TS
#p fopt=(calcfc,ts) freq=noraman ump2/6-31g(d) iop(1/11=1) iop(5/13=1) optcyc=100
IRC 过渡态确认以后（频率计算），进行反应路径扫描
Forward
# irc=(forward,calcfc,maxpoint=10,maxcycle=400) method/basis set cp
Reverse
# irc=(reverse,calcfc,maxpoint=10,maxcycle=400) method/basis set iop(2/15=1,5/13=1)
扫描
柔性扫描_笛卡尔坐标形式
nosymm)
结果中每一步的SCF要检验收敛性 (最好指定
柔性扫描 opt=z-matrix
hf/3-21g opt(z-matrix) nosymm ** 1 1
N H 1 B1 H 1 B2 2 A1
H 1 B3 2 A2 3 D1 O 1 B4 2 A3 3 D2 H 5 B5 1 A4 2 D3 H 5 B6 1 A5 2 D4 H 5 B7 1 A6 2 D5
B1 1.015627 B2 1.015596 B3 1.015722
B4 4.000000 s 20 -0.1 B5 0.991790 B6 0.976661 B7 0.976674 A1 106.126364 A2 106.159425