计算化学公社

标题: 用G4计算反应的能垒，得到的过渡态单点能和振动绝热基态能量过高 [打印本页]

作者
Author: 阿豪0605 时间: 2026-1-13 17:44
标题: 用G4计算反应的能垒，得到的过渡态单点能和振动绝热基态能量过高
请兄弟姐妹们帮我看看是哪里出了问题，在此谢谢大家看帖和给出指导
这是我验证的反应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

作者
Author: sobereva 时间: 2026-1-14 05:39
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)。

作者
Author: 阿豪0605 时间: 2026-1-14 10:37

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

作者
Author: sobereva 时间: 2026-1-14 23:56

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

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)