请教如何解决用TS找过渡态无法收敛的问题

scott__

各位前辈好，我用TS方法计算过渡态，出现报错，原因是 Optimization stopped. -- Number of steps exceeded,  NStep= 288-- Flag reset to prevent archiving.根据站上的指导，优化时能量的变化曲线，发现震荡不大

这是一个花生四烯酸和单线态氧反应过程的过渡态计算
后来我尝试把上边最后的结构保存，用calcall计算opt=(calcall,ts,noeigen) freq b3lyp/6-31g(d,p)，可以正常结束但是得到的能量曲线如下，并且没有虚频，无法正确找到过渡态

不好意思，量子化学不是很懂，刚刚开始学习，希望能够得到指点。

sobereva

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

scott__

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

谢谢老师回复，这个是后边放大之后的曲线，然后下边是输入文件
%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

sobereva

你那个根本不是花生四烯酸，C15少了氢。如果是自由基需要说清楚

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

放大得还不够

scott__

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

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

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

老师，这个我是初学者，好像不太理解怎样是震荡明显，我看着好像后边五个点的值一样，然后跟前一个点差0.000002。

sobereva

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

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

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

scott__

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

嗯嗯，谢谢老师