第一次用高斯，优化过渡态的时候出现结构不收敛，求助各位大佬

zyhhh

优化到一半报错，Error termination request processed by link 9999.
Error termination via Lnk1e in /opt/gaussian/g16/l9999.exe at Mon Sep 29 17:06:36 2025.
查看之前的帖子发现可能是优化不收敛，并在输出文件中找到了Optimization stopped.
    -- Number of steps exceeded,  NStep= 124。
请问接下来改怎么更改？
下面是我的输入文件
%chk=/home/user01/gaosi/jisuan
# opt=(calcfc,qst2,noeigentest) freq b3lyp/6-31g(d) geom=connectivity

Title Card Required

-1 2
C                  2.23473670    0.32969557    0.58582803
O                  2.92224638    0.88171333    1.48370185
O                  1.72632925    1.10444255   -0.50329500
H                  2.32274028    1.83515044   -0.68206826
C                  1.94089442   -1.18019355    0.65993312
H                  1.83767458   -1.57248060   -0.33019643
H                  2.74748217   -1.67650334    1.15792187
C                  0.63472824   -1.41116759    1.44235301
C                  0.34088596   -2.92105671    1.51645810
H                  1.14747371   -3.41736649    2.01444685
H                  0.23766612   -3.31334376    0.52632855
C                 -0.96528023   -3.15203074    2.29887799
O                 -1.30722524   -2.34740896    3.20399012
O                 -1.78957496   -4.28084957    1.99687424
H                 -2.70696499   -4.06532635    2.18002794
C                 -0.52615507   -0.69685257    0.72562154
O                 -0.49880361   -0.54958251   -0.52383194
O                 -1.63519933   -0.20091238    1.47991866
H                 -2.00248256    0.56933948    1.04012883
O                  0.77267625   -0.88689610    2.76561026
H                  1.49634376   -1.33218339    3.21240391
S                 -4.10525384    3.56531382    0.66467171
O                 -5.77525384    3.56533435    0.66467171
O                 -3.54860653    1.99081587    0.66467171
O                 -3.54857749    4.35255252    2.02822100
O                 -3.54857749    4.35255252   -0.69887758

1 2 2.0 3 1.0 5 1.0
2
3 4 1.0
4
5 6 1.0 7 1.0 8 1.0
6
7
8 9 1.0 16 1.0 20 1.0
9 10 1.0 11 1.0 12 1.0
10
11
12 13 2.0 14 1.0
13
14 15 1.0
15
16 17 2.0 18 1.0
17
18 19 1.0
19
20 21 1.0
21
22 23 1.0 24 1.0 25 1.0 26 1.0
23
24
25
26

Title Card Required

-1 2
C                 -0.03696858   -0.24953789    0.00000000
O                  0.65047874    0.30255572    0.89787493
O                 -0.54533361    0.52512458   -1.08920294
H                  0.05107049    1.25583423   -1.26799208
C                 -0.33078022   -1.75942881    0.07418979
H                 -0.43393627   -2.15178444   -0.91591924
H                  0.47579140   -2.25568674    0.57225634
C                 -1.63698437   -1.99038059    0.85655288
C                 -1.93079601   -3.50027151    0.93074267
H                 -1.12422439   -3.99652944    1.42880921
H                 -2.03395206   -3.89262714   -0.05936636
C                 -3.23700016   -3.73122329    1.71310575
O                 -3.57901374   -2.92654886    2.61814516
O                 -4.06125221   -4.86008127    1.41113189
H                 -4.97865732   -4.64456692    1.59422043
C                 -2.79784446   -1.27614020    0.13970944
O                 -2.77042728   -1.12895305   -1.10975237
O                 -3.90694185   -0.78017511    0.89391206
H                 -4.67205103    1.29745429    0.40602177
O                 -1.49912161   -1.46601746    2.17978270
H                 -0.77546857   -1.91125822    2.62664614
S                 -6.29738233    2.72461823    0.08852964
O                 -7.96042266    2.58480403    0.14893537
O                 -5.60668529    1.20865809    0.20558714
O                 -5.77221031    3.67509384    1.35726542
O                 -5.85021106    3.42991697   -1.35766936

1 2 2.0 3 1.0 5 1.0
2
3 4 1.0
4
5 6 1.0 7 1.0 8 1.0
6
7
8 9 1.0 16 1.0 20 1.0
9 10 1.0 11 1.0 12 1.0
10
11
12 13 2.0 14 1.0
13
14 15 1.0
15
16 17 2.0 18 1.0
17
18
19 24 1.0
20 21 1.0
21
22 23 1.0 24 1.0 25 1.0 26 1.0
23
24
25
26

wal

摆好过渡态初猜 用TS
opt=(calcfc,TS,noeigentest)

Stardust0831

确认一下当前的自旋多重度和净电荷是否可靠。
当前体系真的有单电子/自由基么？

zyhhh

Stardust0831 发表于 2025-9-29 22:14
确认一下当前的自旋多重度和净电荷是否可靠。
当前体系真的有单电子/自由基么？

我的体系里是有一个硫酸根自由基的，所以我这么设置的电荷和spin

zyhhh

wal 发表于 2025-9-29 21:26
摆好过渡态初猜 用TS
opt=(calcfc,TS,noeigentest)

用TS的话是不是我就要换成过渡态的结构而不是像这样的初末态，请问一下这两者计算有什么区别

wal

zyhhh 发表于 2025-9-30 00:17
用TS的话是不是我就要换成过渡态的结构而不是像这样的初末态，请问一下这两者计算有什么区别

是
区别在于qst2非常麻烦而且对于你当前要算的过渡态效率远不如ts。如果你执意要用qst2，你要自己想办法把初末结构基本对齐。你当前输入用connection检查一下能看到对应关系是有问题的

zyhhh

wal 发表于 2025-9-30 00:47
是
区别在于qst2非常麻烦而且对于你当前要算的过渡态效率远不如ts。如果你执意要用qst2，你要自己想办法 ...

好的谢谢，那用TS的话结构是不是就只能自己去猜呢，有没有什么方法能去验证我这个猜的对不对

sobereva

仔细把此文看了
简谈Gaussian里找过渡态的关键词opt=TS和QST2、3
http://sobereva.com/460（http://bbs.keinsci.com/thread-12093-1-1.html）

建议有机会参加北京科音初级量子化学培训班（http://www.keinsci.com/KEQC）（下一届在11月底），里面极为详细讲了过渡态怎么搜索，手把手教，给了近20个例子和练习，学一遍并举一反三就会知道怎么又好又快地找过渡态了。

zyhhh

sobereva 发表于 2025-9-30 06:26
仔细把此文看了
简谈Gaussian里找过渡态的关键词opt=TS和QST2、3
http://sobereva.com/460（http://bbs.k ...

好的谢谢大佬，还想请教一下。如果我用TS计算的话，体系可能会有一个或三个自由基，这个时候自旋度是不是2或6呢，还是说我需要每个都算一下确定能量最低的那个

wal

zyhhh 发表于 2025-9-30 10:34
好的谢谢大佬，还想请教一下。如果我用TS计算的话，体系可能会有一个或三个自由基，这个时候自旋度是不是 ...

哪里有三个自由基

zyhhh

wal 发表于 2025-9-30 10:38
哪里有三个自由基

就是硫酸根自由基，氢自由基和二氧化碳自由基。我看得文献过渡态是这么描述的，可能不太对

wal

zyhhh 发表于 2025-9-30 10:44
就是硫酸根自由基，氢自由基和二氧化碳自由基。我看得文献过渡态是这么描述的，可能不太对

反应前是硫酸根自由基 反应后是羧基自由基，没有三个，自旋按2就行

zyhhh

wal 发表于 2025-9-30 14:24
反应前是硫酸根自由基 反应后是羧基自由基，没有三个，自旋按2就行

想问一下我这个貌似还是没有收敛，需不要添加隐式溶剂模型呢

wal

zyhhh 发表于 2025-9-30 19:44
想问一下我这个貌似还是没有收敛，需不要添加隐式溶剂模型呢

这个输入看着还行 但是你应该贴输出文件
如果你实际体系里有溶剂，那是要加的

zyhhh

wal 发表于 2025-9-30 19:50
这个输入看着还行 但是你应该贴输出文件
如果你实际体系里有溶剂，那是要加的

输出文件太大了，显示无法上传，这是最后一段 I HAVE NOT FAILED, I HAVE ONLY DISCOVERED 10,000 WAYS THAT DIDN'T WORK.  -- THOMAS A. EDISON Error termination request processed by link 9999. Error termination via Lnk1e in /opt/gaussian/g16/l9999.exe at Tue Sep 30 18:37:40 2025. Job cpu time:       0 days  7 hours 53 minutes 23.4 seconds. Elapsed time:       0 days  7 hours 54 minutes 45.3 seconds.
File lengths (MBytes):  RWF=    556 Int=      0 D2E=      0 Chk=     11 Scr=      1
我的溶剂是水。我看到一篇跟我类似的文献用了SMD然后他还用到了M06-2X这个泛函和def2-SV基组，但在我的gview里没有找到这个函数,我就按网上的复制粘贴了一下，不知道可对.