计算长碳链烯酸跟臭氧反应时为啥过渡态很难搜索到？而且容易报错，IRC只往一个方向跑

liwenhao0701

在计算长碳链烯酸跟臭氧反应时，为啥过渡态很难搜索到？而且容易报错，IRC只往一个方向跑

zjxitcc

可以先检查：
（1）如果体系有可能是开壳层单重态，注意在计算时要构造对称破缺初猜，检验波函数稳定性，确保获得稳定波函数。如果明知是开壳层单重态，向导师和审稿人隐瞒此信息，强行使用闭壳层RKS计算，结果往往定性错误。
（2）如果体系不是单重态（无论闭壳层单重态/开壳层单重态），注意在SCF收敛后检验波函数稳定性，确保获得稳定波函数。
（3）在跑IRC时应当将forward/reverse分成两个任务提交，且两个任务均读取过渡态chk文件中的波函数。不应该偷懒塞在一个文件里跑两个方向。具体操作时，应当将过渡态chk文件复制两份，例如xxx_irc_f.chk, xxx_irc_r.chk，两个任务分别读取和更新各自的chk文件，不会互相覆盖。

sobereva

过渡态可以用柔性扫描等方式帮助给出合适的过渡态初猜。同时楼上说的问题必须注意
如果你之前已经找到了过渡态，记得用stable结合guess=read关键词读取并测试之前优化过渡态产生的波函数稳定性

liwenhao0701

zjxitcc 发表于 2026-5-15 17:14
可以先检查：
（1）如果体系有可能是开壳层单重态，注意在计算时要构造对称破缺初猜，检验波函数稳定性， ...

你好，这是我高斯数输入文件，可以帮我看下哪里有问题吗？出书来的IRC就两个点

%nprocshared=48
%mem=32GB
%oldchk=ts1_freq.chk
%chk=ts1_irc_forward_gs2.chk
#p irc=(calcfc, RCFC, GS2, forward, MaxPoints=100, StepSize=3, tight) m062x/6-31g(d) empiricaldispersion=gd3 int=ultrafine scf=xqc geom=allcheck

IRC forward (GS2 tight)

0 1

zjxitcc

liwenhao0701 发表于 2026-5-27 21:02
你好，这是我高斯数输入文件，可以帮我看下哪里有问题吗？出书来的IRC就两个点

%nprocshared=48

您这么问，看起来是没有理解我在2L回复的3点内容，这说明您可能非常不熟悉 当体系可能是开壳层单重态时的DFT计算（包括结构优化、找过渡态、频率分析、走IRC），因此计算可能有明显错误。此时不应该从IRC输入文件开始检查，而应该从找过渡态的输入/输出文件开始检查，因此请你展示找过渡态的输入文件，先检查过渡态是否合理。一定要从gjf文件中复制关键词、粘贴到论坛上；如果你还想贴坐标，内容超过50行，建议直接上传gjf文件。
万一过渡态的计算步骤或计算结果不合理，则不能讨论IRC。

liwenhao0701

好的，目前我刚开始做，确实有点不熟悉，以下我粘贴相关的ts的输入文件如下：
%nprocshared=48
%mem=32GB
%chk=yousuanjiazhi+o3ts1.chk
#p opt=(calcfc,ts,noeigen,recalc=3,maxcycles=100) freq 6-31g(d)
empiricaldispersion=gd3 int=ultrafine m062x scf=xqc

TS1 search (M062X/6-311Gdp, SMD, GD3): Methyl Oleate + O3 -> POZ

0 1
C                  9.62576200   -2.63684100   -0.18453800
O                  8.20345200   -2.67007400   -0.28708000
C                  7.59262600   -1.49695800   -0.04034100
C                  6.09283000   -1.62246700   -0.17051600
O                  8.18945900   -0.48923900    0.24647300
C                  5.37570700   -0.30558100    0.09644300
C                  3.85997100   -0.43491000   -0.02961100
C                 -2.89401300    3.32926400    0.31455800
C                 -1.34965600    3.23697100    0.55116100
C                 -3.57969100    2.22290700   -0.47529900
C                 -0.63023200    1.93708900    0.23896500
C                  0.88807400    2.07641600    0.36601100
C                  1.61379700    0.75739600    0.11204100
C                  3.13096000    0.88106600    0.23217200
C                 -6.78811600   -4.70485400   -0.47720200
C                 -6.08739500   -3.62641100    0.34685500
C                 -7.15827800   -5.93067000    0.35375200
C                 -5.71187400   -2.39333600   -0.47175100
C                 -5.00945100   -1.31656000    0.35243100
C                 -4.63595300   -0.08383300   -0.46744900
C                 -3.92786100    0.98968900    0.35663200
H                  9.96415300   -3.64703700   -0.40962600
H                 10.04323300   -1.92288900   -0.89799400
H                  9.92903600   -2.34463100    0.82327900
H                  5.76063900   -2.40619000    0.52150800
H                  5.87348100   -2.00101900   -1.17642900
H                  5.74706800    0.45426700   -0.60095400
H                  5.64036500    0.05133100    1.09854100
H                  3.49653600   -1.19802000    0.67311400
H                  3.60412700   -0.79773900   -1.03506100
H                 -3.39593700    3.47684400    1.28016300
H                 -1.13368700    3.54965800    1.58139500
H                 -2.96349100    1.95536300   -1.34113500
H                 -4.50632200    2.64864000   -0.87985200
H                 -0.88808200    1.60492800   -0.77299100
H                 -0.98150800    1.16737900    0.93589000
H                  1.13661500    2.44530000    1.37060100
H                  1.24240000    2.83602000   -0.34033900
H                  1.35723500    0.38873700   -0.89092400
H                  1.25291300   -0.00181800    0.82038600
H                  3.49071700    1.64265400   -0.47289400
H                  3.38866400    1.24600700    1.23588800
H                 -6.13622700   -5.00360700   -1.30821500
H                 -7.69154500   -4.27969100   -0.93304700
H                 -5.18186400   -4.05004700    0.80378400
H                 -6.73839600   -3.32606200    1.17998400
H                 -7.65864200   -6.69188400   -0.25150800
H                 -6.26582500   -6.38607400    0.79598700
H                 -7.83122400   -5.65735300    1.17338400
H                 -6.61766500   -1.96919900   -0.92704800
H                 -5.06212000   -2.69416600   -1.30548100
H                 -4.10263100   -1.74032300    0.80640200
H                 -5.65828300   -1.01579300    1.18680900
H                 -5.54202300    0.34375900   -0.91834900
H                 -3.98938200   -0.38424500   -1.30351300
H                 -4.57242200    1.28972400    1.19429400
H                 -3.02503200    0.55640200    0.80448100
O                 -1.82751700    5.23838000   -0.21512600
O                 -3.01963300    4.51434300   -0.47597400
O                 -0.84927600    4.23040500   -0.33659400