计算化学公社

标题: 请教如何解决用TS找过渡态无法收敛的问题 [打印本页]

作者
Author: scott__ 时间: 2024-7-30 20:46
标题: 请教如何解决用TS找过渡态无法收敛的问题

各位前辈好，我用TS方法计算过渡态，出现报错，原因是 Optimization stopped. -- Number of steps exceeded, NStep= 288-- Flag reset to prevent archiving.根据站上的指导，优化时能量的变化曲线，发现震荡不大
(, 下载次数 Times of downloads: 26)
这是一个花生四烯酸和单线态氧反应过程的过渡态计算
后来我尝试把上边最后的结构保存，用calcall计算opt=(calcall,ts,noeigen) freq b3lyp/6-31g(d,p)，可以正常结束但是得到的能量曲线如下，并且没有虚频，无法正确找到过渡态
(, 下载次数 Times of downloads: 27) (, 下载次数 Times of downloads: 34)
不好意思，量子化学不是很懂，刚刚开始学习，希望能够得到指点。

作者
Author: sobereva 时间: 2024-7-30 21:35
99%的原因是过渡态初猜不合适。没具体输入文件也没法做更多判断
参考
简谈Gaussian里找过渡态的关键词opt=TS和QST2、3
http://sobereva.com/460（http://bbs.keinsci.com/thread-12093-1-1.html）

“震荡不大”明显是误解。你把末期曲线放大了看，绝对震荡明显

另外，当前分子内色散作用对构象影响明显，最好带上D3校正
DFT-D色散校正的使用
http://sobereva.com/210
谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413

作者
Author: scott__ 时间: 2024-7-30 22:06

sobereva 发表于 2024-7-30 21:35
99%的原因是过渡态初猜不合适。没具体输入文件也没法做更多判断
参考
简谈Gaussian里找过渡态的关键词opt ...

(, 下载次数 Times of downloads: 25) [attach]96132[/attach][attach]96132[/attach]
谢谢老师回复，这个是后边放大之后的曲线，然后下边是输入文件
%chk=C:\Users\H\Desktop\C20H32O2.chk
# opt=(calcfc,ts,noeigen) freq=noraman b3lyp/6-31g(d,p)

Title Card Required

0 2
O                -5.73949200 4.44270200 -1.64318800
O                -5.75147200 4.26119600 0.57747600
C                5.35108000 1.41612800 0.72893600
C                6.82107800 1.73954700 0.41626900
C                5.05884100 -0.08615700 0.94229900
C                7.07271300 3.24936100 0.27857100
C                5.41325100 -0.95195400 -0.24323600
C                8.53846300 3.57620800 -0.03387400
C                -5.20717200 1.44597500 0.53676600
C                -4.55621800 0.07083000 0.34392700
C                4.70110600 -1.97439100 -0.73186300
C                3.38414800 -2.50098400 -0.20652500
C                -4.90399600 2.36306900 -0.65617900
C                -4.71549600 -0.88516600 1.50183200
C                -2.64440200 -2.27750200 1.02994800
C                2.41444600 -2.81019200 -1.32864200
C                -3.84082300 -1.85346500 1.81141900
C                0.24770100 -2.32658800 -0.08433600
C                1.08075100 -2.72348400 -1.28236800
C                -2.10377300 -1.33166100 0.00518100
C                -0.85502300 -1.36708900 -0.47165000
C                -5.49590100 3.75895100 -0.48798400
H                -5.22582800 4.06162900 -2.36954800
H                4.72056200 1.78736800 -0.08710600
H                5.06382800 1.95189100 1.64071700
H                7.11903600 1.24957800 -0.51688900
H                7.45684900 1.34146800 1.21709300
H                3.99652700 -0.20239700 1.17667600
H                5.62409700 -0.42892200 1.81963500
H                6.78107300 3.74780900 1.21019300
H                6.43499700 3.64644500 -0.51955300
H                6.35092000 -0.71732500 -0.74158100
H                8.69311300 4.65456400 -0.12765700
H                9.19862600 3.21023300 0.75834300
H                8.85229200 3.10939800 -0.97250100
H                -4.82501800 1.91844600 1.44469600
H                -6.29027400 1.34056200 0.65128300
H                -3.48436800 0.14808300 0.24426300
H                -4.94952800 -0.40084500 -0.56633600
H                5.09485100 -2.49784800 -1.60204000
H                2.94196000 -1.78571200 0.49095400
H                3.57462900 -3.42292800 0.36312100
H                -5.27588500 1.92009400 -1.58646700
H                -3.81685800 2.47655100 -0.75847200
H                -5.59412700 -0.76836100 2.13296600
H                -1.92894600 -2.89077900 1.57000400
H                2.86689700 -3.12988100 -2.26555100
H                -4.02844400 -2.44472500 2.70324700
H                0.88041400 -1.86358300 0.67998400
H                -0.18752500 -3.22655500 0.36062800
H                0.51771600 -2.96526300 -2.18129700
H                -2.79348700 -0.59361000 -0.36669100
H                -0.58485000 -0.64039800 -1.23671200
O                -2.88173400 -3.73374200 -1.00283100
O                -3.34578800 -3.44185600 0.19650000

作者
Author: sobereva 时间: 2024-7-30 22:15
你那个根本不是花生四烯酸，C15少了氢。如果是自由基需要说清楚

甭写毫无意义的noraman，仔细看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http://bbs.keinsci.com/thread-3460-1-1.html）

放大得还不够

作者
Author: scott__ 时间: 2024-7-31 09:16
本帖最后由 scott__ 于 2024-7-31 09:26 编辑

sobereva 发表于 2024-7-30 22:15
你那个根本不是花生四烯酸，C15少了氢。如果是自由基需要说清楚

甭写毫无意义的noraman，仔细看

嗯嗯，谢谢老师，我去学习一下,改正一下我的关键词使用。

老师，这个我是初学者，好像不太理解怎样是震荡明显，我看着好像后边五个点的值一样，然后跟前一个点差0.000002。
(, 下载次数 Times of downloads: 27)

作者
Author: sobereva 时间: 2024-8-1 01:20

scott__ 发表于 2024-7-31 09:16
嗯嗯，谢谢老师，我去学习一下,改正一下我的关键词使用。

老师，这个我是初学者，好像不太理解怎样是 ...

看out文件里面的四个优化收敛的判据的变化情况
总能量不是收敛判据

作者
Author: scott__ 时间: 2024-8-1 10:34

sobereva 发表于 2024-8-1 01:20
看out文件里面的四个优化收敛的判据的变化情况
总能量不是收敛判据

嗯嗯，谢谢老师

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