计算化学公社

 找回密码 Forget password
 注册 Register
Views: 3010|回复 Reply: 2
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 高斯qst2搜寻过渡态出现错误l101

[复制链接 Copy URL]

1

帖子

0

威望

57

eV
积分
58

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主

错误代码
Inconsistent reactant and product variables.
Error termination via Lnk1e in E:\G09W\l101.exe at Thu Jul 14 18:17:53 2022.
Job cpu time:       0 days  0 hours  0 minutes  1.0 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      1

输入文件
%mem=1000MB
%chk=D:\DESK\QST2.chk
# opt=(calcfc,qst2,modredundant) freq b3lyp/6-31++g
scrf=(smd,solvent=chlorobenzene) geom=connectivity int=ultrafine

QST2

0 1
C                 -2.01221207    0.19824126    2.74692032
C                 -1.32705607   -0.76428138    1.98302124
C                 -0.70546440   -0.36659911    0.79643289
C                 -0.76378376    0.97615632    0.37195626
C                 -1.44408938    1.93239437    1.13061324
C                 -2.06999001    1.53216653    2.32534822
C                  0.50204156    0.43514398   -1.51252771
C                  0.56071188   -0.90742794   -1.08775847
C                  1.24396209   -1.86312844   -1.84460966
C                  1.87210640   -1.46241542   -3.03788575
C                  1.81378584   -0.12846945   -3.45995684
C                  1.12609521    0.83345357   -2.69782367
O                 -0.16363240    1.43476698   -0.80159900
O                 -0.04120456   -1.36653915    0.08468506
Cl                 1.04154356    2.45388316   -3.18743563
Cl                 2.57871405    0.33866528   -4.90520067
Cl                 2.70772803   -2.61227085   -3.97149627
Cl                 1.30044712   -3.47135188   -1.31255275
Cl                -1.23840226   -2.38359826    2.47560651
Cl                -2.77047147   -0.26827361    4.19578725
Cl                -1.49935513    3.54103276    0.59970148
H                 -0.77190358    1.64631001    3.65848275
Cl                -2.90100591    2.68309395    3.26170088
H                 -0.59341729    2.30270610    4.00593886

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

B 11 16 F
B 11 12 F
B 10 11 F
B 10 17 F
B 10 9 F
B 9 18 F
B 9 8 F
B 12 15 F
B 12 7 F
B 8 7 F
B 8 14 F
B 14 3 F
B 7 13 F
B 13 4 F
B 4 3 F
B 3 2 F
B 2 19 F
B 2 1 F
B 1 20 F
B 1 6 F
B 6 5 F
B 5 21 F
B 5 4 F

Title Card Required

0 1
C                  1.65325616   -1.12311929   -2.76694766
C                  0.88783197    0.01100823   -2.43826249
C                  0.22773727    0.03381621   -1.20633940
C                  0.32573968   -1.05735655   -0.32132887
C                  1.09253749   -2.17501338   -0.66479441
C                  1.76304229   -2.20791559   -1.89530869
C                 -1.04102989    0.05490844    1.20770764
C                 -1.13941340    1.14849640    0.32217401
C                 -1.88787881    2.27678187    0.66630175
C                 -2.54303645    2.30497079    1.91109517
C                 -2.44558349    1.22041743    2.79049915
C                 -1.69140687    0.08440955    2.44450695
O                 -0.31288868   -1.09734986    0.92240344
O                 -0.51528978    1.18144651   -0.92554881
Cl                -1.55880300   -1.24066028    3.49388068
Cl                -3.24339225    1.27030400    4.29228688
Cl                -3.45916246    3.67123370    2.34485589
Cl                -1.99063368    3.57824169   -0.41536444
Cl                 0.75664996    1.33612320   -3.48816135
Cl                 2.45951321   -1.17234879   -4.27297206
Cl                 1.19166711   -3.48350956    0.41946466
H                  2.37102416   -3.07747622   -2.16332703
Cl                 3.89065973   -0.73859694   -0.60080284
H                  4.90970397   -0.13971903   -0.06213941

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



4289

帖子

4

威望

9541

eV
积分
13910

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2022-7-15 10:57:50 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2022-7-15 11:01 编辑

看了下反应物和产物的结构,二者结构相差过大,QST2是没法产生过渡态初始结构的。你需要手动构造一个更合理的初始结构,然后用opt=ts或qst3。此外还有些问题:
(1)当前计算没有意义,当本科生作业都不行,基组有点可怕,仔细阅读《谈谈量子化学中基组的选择》http://sobereva.com/336
(2)用B3LYP,却没有加色散。看
《谈谈“计算时是否需要加DFT-D3色散校正?”》http://sobereva.com/413
《坚持使用B3LYP算有机体系能量的人的下场》http://bbs.keinsci.com/thread-12773-1-1.html
(3)只给1GB内存,太少了,容易内存不足报错,或算得超级慢。
自动做多参考态计算的程序MOKIT

211

帖子

0

威望

1466

eV
积分
1677

Level 5 (御坂)

3#
发表于 Post on 2022-7-15 17:30:54 | 只看该作者 Only view this author
modredundant去掉,101还像是命令行出错了

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 07:11 , Processed in 0.181435 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list