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