关于用高斯优化Encounter complex问题求助

快乐如我

请教各位老师，我最近在做一个SiO2 带正电的cluster活化H2O分子，这个反应应该有反应能垒的，可我开始优化的第一步Encounter complex就直接生成产物了，不管我怎么调节反应物之间的距离都是这样，我是用ump2/6-31g(d) 优化。以下是输入文件：%chk=SiO2-H2O-ec1.chk
%mem=8000MB
# opt=(cartesian,gdiis,maxcycle=1000) freq ump2/6-31g(d)
iop(5/13=1,1/11=1,2/16=1) scf=maxcycle=600

Title Card Required

1 2
Si                 1.94895663    1.05941539   -0.62206440
O                  3.39739807    0.46344931   -0.62206440
O                  0.39629554    0.85348390   -0.62206440
O                  7.56043969    2.72063872   -2.08249616
H                  6.78994581    2.27171856   -1.72695452
H                  8.30070795    2.11023312   -2.11424384

输出结果如图

然后我又用了势能面扫描来优化cluster和H2O之间的距离，可优化三步就出错了，请问改怎么解决呢？下面是出错信息：
New curvilinear step not converged.
RedCar failed.
Error termination via Lnk1e in /afs/math/software/gaussian/g09/@sys/l103.exe

分子在优化过程中会变成线性，所以我加了虚原子，下面是势能面扫描的输入文件：
%chk=SiO2-H2O-scan3.chk
%mem=8000MB
# opt=(modredundant,maxcycle=1000) ump2/6-31g(d) nosymm
iop(5/13=1,1/11=1,2/16=1) scf=maxcycle=600

Title Card Required

1 2
Si                 2.04510342    1.25555581   -0.67199734
O                  3.49354486    0.65958973   -0.67199734
O                  0.49244233    1.04962432   -0.67199734
O                  7.29661705    3.39139061    2.66989626
H                  6.55998643    2.77582604    2.66244114
H                  8.11018130    2.91537754    2.85189834
X                  1.99984832    2.75257991   -0.12854394

B 2 5 S 45 -0.100000


[size=21.3333px]

sobereva

优化时候用小步长上限，诸如opt(maxstep=10,notrust)，摆成不同的距离，如果优化完了都成产物了，说明就是没势垒的过程。这种反应没势垒不奇怪。


PS：

不要用UMP2，诸如M062X明显快得多也好得多。基组建议用TZVP。

不要把maxcyc设那么大，毫无意义还可能白费很多计算时间，尤其是绝对不要用5/13=1！这里都提到了，文末专门狠狠批判了5/13=1：
解决SCF不收敛问题的方法
http://sobereva.com/61

1/11=1毫无意义，opt根本就不做Hessian本征值测试
2/16=1也没有任何用处。

nosymm加得也完全多余。多余的关键词一律不要写。

快乐如我

sobereva 发表于 2015-12-3 01:15
优化时候用小步长上限，诸如opt(maxstep=10,notrust)，摆成不同的距离，如果优化完了都成产物了，说明就是 ...

谢谢老师，可能是我用的方法问题，我换成Ｍ062X方法就可以优化出来了。因为这段时间刚开始学习量化，的确好多都不懂，需要好好学习一下，谢谢老师指导