用G4计算反应的能垒，得到的过渡态单点能和振动绝热基态能量过高

阿豪0605

请兄弟姐妹们帮我看看是哪里出了问题，在此谢谢大家看帖和给出指导
这是我验证的反应cyc-C6H12 + NO2 → C6H11 + cis-HONO

下面是该反应的过渡态G4输出文件的能量值，也有正常终止，坐标用的是文献附件给的
Temperature=              298.150000 Pressure=                      1.000000
E(ZPE)=                     0.176579 E(Thermal)=                    0.187662
E(QCISD(T))=             -439.665635 E(Empiric)=                   -0.187769
E(CCSD(T))=                -0.187769 E(Empiric)=                   -0.025401
DE(Plus)=                  -0.368915 DE(2DF)=                      -0.573841
E(Delta-G3XP)=             -0.046009 DE(HF)=                     -440.690991
G4(0 K)=                 -440.679908 G4 Energy=                  -440.678963
G4 Enthalpy=             -440.734306 G4 Free Energy=   
下面是过渡态的输入文件，然后我看了G4输出文件的结构直接干成产物了，频率从-1768.18变成了10，键长变长了两倍，从文献的1.25变成了2.794
%chk=gaiC6H12_TS1_G4_SP.chk
%mem=100GB
%nproc=32
# G4 scf=qc scf=maxcyc=500 int=ultrafine ccsd=maxcyc=500

SPECIES NAME/COMMENTS: C 20/05/2015 17:53:49

0 2
6        2.018878000      1.259465000     -0.073126000
6        0.479414000      1.261585000     -0.127539000
6        0.021946000      0.000000000     -0.810331000
6        0.479414000     -1.261585000     -0.127538000
6        2.018878000     -1.259465000     -0.073126000
6        2.529622000      0.000000000      0.623840000
1       -1.306637000      0.000000000     -0.811295000
1        0.084161000      1.281388000      0.896801000
1        0.112749000      2.153970000     -0.645917000
1        2.418256000      1.298038000     -1.096846000
1        2.369297000      2.161686000      0.441540000
1        0.084161000     -1.281387000      0.896802000
1        0.112749000     -2.153970000     -0.645916000
1        2.369297000     -2.161686000      0.441541000
1        2.418256000     -1.298039000     -1.096846000
1        2.182325000      0.000000000      1.667403000
1        3.625614000      0.000000000      0.653346000
1        0.184184000      0.000000000     -1.896257000
7       -3.037677000      0.000000000      0.300896000
8       -2.276890000      0.000001000      1.210676000
8       -2.558059000      0.000000000     -0.883640000

--Link1--
%chk=gaiC6H12_TS1_G4_SP.chk
# Geom=AllCheck ccsd=(T,T1diag)

SPECIES NAME/COMMENTS: C T1 Diagnostic Calculation 20/05/2015 17:53:49

0 2

sobereva

阿豪0605 发表于 2026-1-14 10:37
谢谢sobereva老师指导，刚入门拿着师哥的输入文件改的，云里雾里，才造成简单的错误
那像CBS-QB3这种复 ...

所有的热力学组合方法的用法都一样
不懂的关键词一律别写、不懂的用法一律别用，绝对不要随便照抄别人的关键词，否则只会自取其辱

建议好好看看
谈谈学量子化学如何正确地入门
http://sobereva.com/355（http://bbs.keinsci.com/thread-4447-1-1.html）

阿豪0605

sobereva 发表于 2026-1-14 05:39
G16 A.03版的G4输出有bug，换其它版本

G4若要用来算过渡态必须同时写上找过渡态的关键词，诸如opt(TS,ca ...

谢谢sobereva老师指导，刚入门拿着师哥的输入文件改的，云里雾里，才造成简单的错误
那像CBS-QB3这种复合方法也是要先加上寻找过渡态的关键词嘛


%chk=C6H12_TS1_CBS.chk
%mem=50GB
%nproc=40
# CBS-QB3=SP

SPECIES NAME/COMMENTS: C 30/10/2019 23:43:23

0 2
6        1.781272000     -1.000414000     -0.565377000
6        2.172992000      0.463371000     -0.782084000
6        0.985391000      1.396433000     -0.533520000
6        0.420330000      1.204833000      0.877213000
6        0.106558000     -0.245421000      1.159291000
6        1.215622000     -1.223349000      0.841710000
1        1.279175000      2.438707000     -0.677162000
1        2.987954000      0.726636000     -0.096393000
1        2.555108000      0.602090000     -1.796438000
1        1.020818000     -1.282883000     -1.301770000
1        2.641075000     -1.654930000     -0.726196000
1        1.175947000      1.526233000      1.611169000
1       -0.465035000      1.822143000      1.036620000
1       -0.420800000     -0.428252000      2.097021000
1       -0.813522000     -0.588037000      0.288537000
1        2.015104000     -1.069599000      1.583217000
1        0.868580000     -2.252144000      0.965626000
1        0.197723000      1.183212000     -1.265853000
8       -1.601578000     -0.988439000     -0.610072000
7       -2.655204000     -0.250427000     -0.610177000
8       -2.717009000      0.633077000      0.181255000

--Link1--
%chk=C6H12_TS1_CBS.chk
# Geom=AllCheck ccsd=(T,T1diag)

SPECIES NAME/COMMENTS: C T1 Diagnostic Calculation 30/10/2019 23:43:23

0 2

sobereva

G16 A.03版的G4输出有bug，换其它版本

G4若要用来算过渡态必须同时写上找过渡态的关键词，诸如opt(TS,calcfc,noeigen)，否则自然给你优化到极小点


int=ultrafine本来就是G16默认的，根本不需要写
慎用scf=qc，容易收敛到不稳定波函数，下文说了
解决SCF不收敛问题的方法
http://sobereva.com/61

--link1--后面都不知道想干嘛，若要做T1诊断做CCSD就够了根本不需要CCSD(T)。