关于gaussian QST2计算过渡态出错的问题

ddls2017

各位老师，本人在用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

ddls2017

补充说明下，之前用的是G09-B01版本算一个类似的过渡态问题，没出错。这次用了G09-D01，就出错了。而且用D01版算之前那个问题也同样出错了，不知是为何？

sobereva

程序bug
高斯的QST2经常容易遇到莫名其妙的bug
这体系也不大，用opt=TS就完了，完全没必要用QST2

ddls2017

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

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

sobereva

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

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

ddls2017

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

好的，谢谢Sob老师

YuLONG

请问这种BUG有解决的办法吗？ @sobereva

sobereva

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

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

YuLONG

那有没有替代的方法呢？ 除了用猜测和TS 优化的方法。@sobereva
谢谢您的回复

Frank

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
请问有办法解决吗？

sobereva

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）

sobereva

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

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

Frank

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

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

sobereva

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

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

Frank

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

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