计算化学公社

标题: 关于用高斯优化Encounter complex问题求助 [打印本页]

作者
Author: 快乐如我 时间: 2015-12-2 23:28
标题: 关于用高斯优化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]

作者
Author: sobereva 时间: 2015-12-3 01:15
优化时候用小步长上限，诸如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加得也完全多余。多余的关键词一律不要写。

作者
Author: 快乐如我 时间: 2015-12-3 17:22

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)