计算化学公社

标题: 关于gaussian QST2计算过渡态出错的问题 [打印本页]

作者
Author: ddls2017 时间: 2017-9-7 10:02
标题: 关于gaussian QST2计算过渡态出错的问题
各位老师，本人在用qst2算过渡态时，老是一启动就报错：
Old curvilinear step not converged, using linear step:
SCX= 6.31D+00 DXMaxT= 2.13D+00 SCLim= 1.07D+00 Fact= 1.69D-01
RedCar/ORedCr failed for GTrans.
Error termination via Lnk1e in C:\G09W\l101.exe at Thu Sep 07 09:29:55 2017.
Job cpu time:    0 days  0 hours  0 minutes  0.0 seconds.

请大神指点迷津，不胜感谢！

输入文件如下：
%nprocshared=4
%rwf=D:\USER\1qst2.rwf
%nosave
%chk=D:\USER\1qst2.chk
%mem=200MW
# opt=(calcfc,qst2,maxcycle=300) ub3lyp/6-31g(d) scf=qc geom=connectivity
gfinput iop(1/8=2)

1qst2

0 1
C                -1.20451500 -0.97975800 0.00000000
C                -1.13431100 0.44081200 0.00000000
N                0.00000000 -1.65056900 0.00000000
N                0.00000000 1.10181500 0.00000000
C                1.13431100 0.44081200 0.00000000
C                1.20451500 -0.97975800 0.00000000
O                0.00000000 -2.95227300 0.00000000
N                2.35801800 1.24960700 0.00000000
O                3.43770300 0.63062800 0.00000000
O                2.25314900 2.46307800 0.00000000
N                -2.35801800 1.24960700 0.00000000
O                -2.25314900 2.46307800 0.00000000
O                -3.43770300 0.63062800 0.00000000
N                2.27290900 -1.77290400 0.00000000
N                -2.27290900 -1.77290400 0.00000000
H                2.07991800 -2.76957400 0.00000000
H                3.19197100 -1.35357900 0.00000000
H                -3.19197100 -1.35357900 0.00000000
H                -2.07991800 -2.76957400 0.00000000

1 2 1.5 3 1.5 15 1.5
2 4 2.0 11 1.0
3 6 1.5 7 1.0
4 5 2.0
5 6 1.5 8 1.0
6 14 1.5
7
8 9 1.5 10 2.0
9
10
11 12 2.0 13 1.5
12
13
14 16 1.0 17 1.0
15 18 1.0 19 1.0
16
17
18
19

product-1

0 1
C                -1.20364200 0.08239000 -0.01504400
C                -1.12815900 1.46664600 0.00382400
N                0.00000000 -0.60710200 -0.00487200
N                -0.00000100 2.10912000 0.01856700
C                1.12815700 1.46664800 0.00382700
C                1.04894178 0.08247719 0.00799904
O                0.00000100 -1.90335900 0.02640300
N                2.52140548 2.22556188 -0.20331644
O                2.81443734 2.44662901 -1.34908061
O                3.06797733 2.42529201 0.84971774
N                -2.52914273 2.23191408 0.11157065
O                -2.88774980 2.40519667 -1.02372487
O                -3.03857381 2.44554510 1.18038474
N                2.06745952 -0.82706567 -0.03939875
N                -2.31713700 -0.70712200 -0.07593700
H                1.76248693 -1.76442311 0.20620673
H                2.96004178 -0.52288036 0.32171256
H                -3.17035500 -0.31144900 0.29162100
H                -2.11854400 -1.67655600 0.15369400

1 2 1.5 3 1.0 15 1.5
2 4 2.0
3 6 1.0 7 1.0
4 5 2.0
5 6 1.5
6 14 1.5
7
8 9 2.0 10 2.0
9
10
11 13 2.0 12 2.0
12
13
14 16 1.0 17 1.0
15 18 1.0 19 1.0
16
17
18
19

作者
Author: ddls2017 时间: 2017-9-7 10:07
补充说明下，之前用的是G09-B01版本算一个类似的过渡态问题，没出错。这次用了G09-D01，就出错了。而且用D01版算之前那个问题也同样出错了，不知是为何？

作者
Author: sobereva 时间: 2017-9-7 10:26
程序bug
高斯的QST2经常容易遇到莫名其妙的bug
这体系也不大，用opt=TS就完了，完全没必要用QST2

作者
Author: ddls2017 时间: 2017-9-7 10:32

sobereva 发表于 2017-9-7 10:26
程序bug
高斯的QST2经常容易遇到莫名其妙的bug
这体系也不大，用opt=TS就完了，完全没必要用QST2

谢谢Sob老师，用TS的话，过渡态不太好猜啊，有什么好方法吗？

作者
Author: sobereva 时间: 2017-9-7 12:51

ddls2017 发表于 2017-9-7 10:32
谢谢Sob老师，用TS的话，过渡态不太好猜啊，有什么好方法吗？

这种体系这么小，反应牵扯的原子又很少，不会不好猜。比这复杂得多的反应，只要稍微有点直觉，通过opt=TS反复尝试都不难找到过渡态。QST2看起来很美，实则烂得很，而很多经验不丰富的人却将QST2视为默认的找过渡态的方法，结果不断被坑。

作者
Author: ddls2017 时间: 2017-9-7 15:37

sobereva 发表于 2017-9-7 12:51
这种体系这么小，反应牵扯的原子又很少，不会不好猜。比这复杂得多的反应，只要稍微有点直觉，通过opt= ...

好的，谢谢Sob老师

作者
Author: YuLONG 时间: 2018-3-10 10:42
请问这种BUG有解决的办法吗？ @sobereva

作者
Author: sobereva 时间: 2018-3-10 16:20

YuLONG 发表于 2018-3-10 10:42
请问这种BUG有解决的办法吗？ @sobereva

别用QST2就是最好的解决办法

作者
Author: YuLONG 时间: 2018-3-10 17:16
那有没有替代的方法呢？除了用猜测和TS 优化的方法。@sobereva
谢谢您的回复

作者
Author: Frank 时间: 2020-2-4 06:07

sobereva 发表于 2018-3-10 16:20
别用QST2就是最好的解决办法

社长，我用Opt=TS也遇到了这个问题，是做的ONIOM(QM:PM7)的限制性优化,只优化溶质，溶剂环境不动。报
错信息是：
Old curvilinear step not converged, using linear step:
SCX= 2.27D-02 DXMaxT= 2.32D+00 SCLim= 1.16D+00 Fact= 1.00D+00
Error imposing constraints
请问有办法解决吗？

作者
Author: sobereva 时间: 2020-2-4 10:01

YuLONG 发表于 2018-3-10 17:16
那有没有替代的方法呢？除了用猜测和TS 优化的方法。@sobereva
谢谢您的回复

“猜测和TS 优化的方法”语义不明
Gaussian里能找过渡态的方法此文都说了
简谈Gaussian里找过渡态的关键词opt=TS和QST2、3
http://sobereva.com/460（http://bbs.keinsci.com/thread-12093-1-1.html）

作者
Author: sobereva 时间: 2020-2-4 10:02

Frank 发表于 2020-2-4 06:07
社长，我用Opt=TS也遇到了这个问题，是做的ONIOM(QM:PM7)的限制性优化,只优化溶质，溶剂环境不动。报
...

去掉限制，或者尝试修改限制

作者
Author: Frank 时间: 2020-2-4 16:23

sobereva 发表于 2020-2-4 10:02
去掉限制，或者尝试修改限制

多谢社长，去掉限制体系就太大了，或许可以尝试下对溶质分子周围的几个溶剂分子取消限制，我正在试opt=Cartesian

作者
Author: sobereva 时间: 2020-2-6 02:50

Frank 发表于 2020-2-4 16:23
多谢社长，去掉限制体系就太大了，或许可以尝试下对溶质分子周围的几个溶剂分子取消限制，我正在试opt=Ca ...

是否做限制原理上不直接影响耗时和体系大小
直接与被计算分子挨着的肯定得取消限制，否则强烈引入人为虚假因素

作者
Author: Frank 时间: 2020-2-6 16:58

sobereva 发表于 2020-2-6 02:50
是否做限制原理上不直接影响耗时和体系大小
直接与被计算分子挨着的肯定得取消限制，否则强烈引入人为虚 ...

不好意思社长，之前说错了。我限制溶剂不动不是因为体系太大，是因为在MD中通过电场使溶剂定向了，我想用ONIOM去模拟这个定向的溶剂是否可以催化反应，如果不限制的话，溶剂又会回到混乱的状态。

作者
Author: sobereva 时间: 2020-2-8 10:08

Frank 发表于 2020-2-6 16:58
不好意思社长，之前说错了。我限制溶剂不动不是因为体系太大，是因为在MD中通过电场使溶剂定向了，我想用 ...

还不如甭用ONIOM，溶剂就用小点的基组，并且Gaussian计算的时候也用field关键词带着电场

作者
Author: fuhua 时间: 2024-9-15 17:23
您是否有限制原子，我也遇到了这个问题，同sobls说的一样，把限制的原子取消就可以运行了。

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