算irc后两端都是产物是否意味着没有能垒

dbxhtw

求助，计算了乙二醇和氯自由基反应的过渡态，对过渡态进行irc发现两侧都是产物的形式，是否代表这个反应没有能垒呢

zjxitcc

有一些建议：
（1）中性体系没必要使用弥散函数，用6-311+g(d,p)没必要，6-311G(d,p)足矣。
（2）对于TS结构，应当检验波函数稳定性，确保波函数稳定。如果不稳定，则频率计算没有意义，跑IRC也没有意义。
（3）IRC计算始终应当读取TS结构算好的轨道和力常数，否则变成从零开始计算，不仅浪费机时，还无法保证与TS结构收敛到同一个解。
（4）对于电子结构稍微复杂一点的体系，尤其是自由基反应和过渡金属反应，IRC（始终）应当两个方向分开跑，提交两个任务，而非偷懒杂糅在一个任务里，这种偷懒容易造成 后跑的那个方向的电子结构有误，从而造成后跑的那个方向有系统性的错误。

那么如何两个方向分开跑，并且读取轨道和力常数？举一个例子，假设我在UPBE0-D3(BJ)/6-311G(d,p)下找到了一个TS，检验过波函数稳定性，也做过freq计算（对应的chk文件称为test_freq.chk)，我先复制两份chk文件，在Linux下就是运行
cp test_freq.chk test_freq_f.chk
cp test_freq.chk test_freq_r.chk

接着我写两个gjf文件，关键词示例如下
(i) forward方向
%chk=test_freq_f.chk
%mem=128GB
%nprocshared=64
#p irc(rcfc,forward,LQA,maxpoints=40) UPBE1PBE chkbasis em=gd3bj nosymm guess=read geom=allcheck

(ii) reverse方向
%chk=test_freq_r.chk
%mem=128GB
%nprocshared=64
#p irc(rcfc,reverse,LQA,maxpoints=40) UPBE1PBE chkbasis em=gd3bj nosymm guess=read geom=allcheck

底下只需要放2个空行，不需要写电荷、自旋多重度、坐标、基组、赝势，因为上述关键词会从chk文件里读取。至于是否写nosymm，应当与我之前TS计算有没写nosymm保持一致。这样做不仅输入文件简单、内容少，不易出错，而且可以从chk文件中读取大量已有信息，节约机时。更重要的是保证了与之前TS收敛到一样的SCF解，forward/reverse两个方向的计算也不会互相干扰。

dbxhtw

zjxitcc 发表于 2026-1-6 10:55
有一些建议：
（1）中性体系没必要使用弥散函数，用6-311+g(d,p)没必要，6-311G(d,p)足矣。
（2）对于TS ...

谢谢老师的建议！我仔细学习一下！