计算化学公社

标题: qst2寻找聚乙烯热解过渡态求助 [打印本页]

作者
Author: Adolphin 时间: 2019-4-23 16:04
标题: qst2寻找聚乙烯热解过渡态求助
新人救助，想要通过聚合度为10的聚乙烯在400k温度下热解生成甲烷和C9长链，09启动说文件错误求指教，输入文件内容如下。
# opt=(calcfc,qst2,noeigen) freq b3lyp/6-31++g(d,p) geom=connectivity temperature=400K

Title Card Required 1

0 1
C                -6.93192150 0.49079754 0.09074382
H                -6.57526708 -0.51801246 0.09074382
H                -6.57524866 0.99519573 -0.78290768
H                -8.00192150 0.49081072 0.09074382
C                -6.41857928 1.21675381 1.34814879
H                -5.34857930 1.21657124 1.34824625
H                -6.77541154 0.71246801 2.22180007
C                -6.93166562 2.66876703 1.34800834
H                -6.57382468 3.17335298 0.47494312
H                -8.00166489 2.66930776 1.34688265
C                -6.41928926 3.39394687 2.60625500
H                -5.34929001 3.39338526 2.60739240
H                -6.77714903 2.88937450 3.47932035
C                -6.93186155 4.84614164 2.60630283
H                -8.00186082 4.84670235 2.60518518
H                -6.57401832 5.35070842 1.73322747
C                -6.41946136 5.57132989 3.86453494
H                -5.34946208 5.57077666 3.86564724
H                -6.77729670 5.06675777 4.73761045
C                -6.93204501 7.02352065 3.86458971
H                -8.00204430 7.02407298 3.86348247
H                -6.57421422 7.52809141 2.99151155
C                -6.41963825 7.74871122 5.12281781
H                -5.34963896 7.74816078 5.12392292
H                -6.77746641 7.24413885 5.99589611
C                -6.93222569 9.20090064 5.12287489
H                -8.00222498 9.20145019 5.12177171
H                -6.57439953 9.70547255 4.24979550
C                -6.41981605 9.92609223 6.38110123
H                -5.34981676 9.92554272 6.38220335
H                -6.77596534 10.93508068 6.38114142
H                -6.77764133 9.42151995 7.25418076

1 2 1.0 3 1.0 4 1.0 5 1.0
2
3
4
5 6 1.0 7 1.0 8 1.0
6
7
8 9 1.0 10 1.0 11 1.0
9
10
11 12 1.0 13 1.0 14 1.0
12
13
14 15 1.0 16 1.0 17 1.0
15
16
17 18 1.0 19 1.0 20 1.0
18
19
20 21 1.0 22 1.0 23 1.0
21
22
23 24 1.0 25 1.0 26 1.0
24
25
26 27 1.0 28 1.0 29 1.0
27
28
29 30 1.0 31 1.0 32 1.0
30
31
32

Title Card Required 2

0 1
C                -5.56933804 11.04734698 7.35926904
H                -5.21266520 11.55174517 8.23292055
H                -5.21266520 11.55174517 6.48561754
H                -6.63933804 11.04736017 7.35926904
C                -5.05602232 9.59541483 7.35926904
H                -5.41269313 9.09101736 8.23292179
H                -3.98602232 9.59540165 7.35926656
C                -5.56936746 8.86945752 6.10186587
H                -6.63936746 8.86946739 6.10187061
H                -5.21270006 9.37385777 5.22821333
C                -5.05604676 7.41752714 6.10186281
H                -5.41271368 6.91312502 6.97551446
H                -3.98604676 7.41751573 6.10185659
C                -5.56939508 6.69157481 4.84445806
H                -5.21432458 7.19710479 3.97080803
H                -6.63939331 6.68989333 4.84544065
C                -5.05378187 5.24045765 4.84304250
H                -3.98378424 5.24213971 4.84154088
H                -5.40842838 4.73507843 5.71695193
C                -5.56773965 4.51428840 3.58601202
H                -5.21310173 5.01967398 2.71210277
H                -6.63773725 4.51259760 3.58751910
C                -5.05211418 3.06317560 3.58458893
H                -3.98211660 3.06486671 3.58306982
C                -5.56608638 2.33700123 2.32756731
H                -5.21144982 2.84237713 1.45365192
H                -6.63608399 2.33531960 2.32908138
C                -7.24369896 8.21552947 2.52506936
H                -6.88704454 7.20671947 2.52506936
H                -6.88702612 8.71992766 3.39872087
H                -6.88702612 8.71992766 1.65141786
H                -8.31369896 8.21554266 2.52506936

1 2 1.0 3 1.0 4 1.0 5 1.0
2
3
4
5 6 1.0 7 1.0 8 1.0
6
7
8 9 1.0 10 1.0 11 1.0
9
10
11 12 1.0 13 1.0 14 1.0
12
13
14 15 1.0 16 1.0 17 1.0
15
16
17 18 1.0 19 1.0 20 1.0
18
19
20 21 1.0 22 1.0 23 1.0
21
22
23 24 1.0 25 2.0
24
25 26 1.0 27 1.0
26
27
28 29 1.0 30 1.0 31 1.0 32 1.0
29
30
31
32

作者
Author: Frank 时间: 2019-4-23 16:36
产物和反应物的原子顺序必须一致

作者
Author: Adolphin 时间: 2019-4-23 17:06
本帖最后由 Adolphin 于 2019-4-23 17:10 编辑

感谢您的解答，这个我已经调整过了但是出现了这个错误请问怎么改？是我选择的基组不正确，还是关键词不正确？本人初步接触这个软件，问的有点多而且基础，谢谢。

作者
Author: lmy 时间: 2019-4-23 19:37
建议把geom=connectivity以及与其相关联的“1 2 1.0 3 1.0 4 1.0 5 ***********”都去掉，然后再试试

作者
Author: sobereva 时间: 2019-4-24 02:25
新人尽量甭用QST2，效果差，用着麻烦，简直是自找不快。

这点专门在这里说了
简谈Gaussian里找过渡态的关键词opt=TS和QST2、3
http://sobereva.com/460（http://bbs.keinsci.com/thread-12093-1-1.html）

当前用的计算级别也忒烂，才2-zeta还加一堆没用的弥散
仔细看此文，省得回头被要求重算
坚持使用B3LYP算有机体系的人的下场
http://bbs.keinsci.com/thread-12773-1-1.html

作者
Author: Adolphin 时间: 2019-4-24 15:00

sobereva 发表于 2019-4-24 02:25
新人尽量甭用QST2，效果差，用着麻烦，简直是自找不快。

这点专门在这里说了

sob老师谢谢你的解答，由于本人是初步接触这个软件，对用ts方法寻找过渡态的方法有点迷茫，请问寻找过渡态是靠自己的知识推测或者猜测出来的一种反应过程中的结构吗？还有一个问题是用TS优化完后如何确定这个物质就是该反应的过渡态呢？最后一点sob老师我提问的输入文件是有错误的吗？望sob老师能在百忙之中给予解答，谢谢。

作者
Author: Frank 时间: 2019-4-24 17:43

Adolphin 发表于 2019-4-24 15:00
sob老师谢谢你的解答，由于本人是初步接触这个软件，对用ts方法寻找过渡态的方法有点迷茫，请问寻找过渡 ...

你问的问题sob的帖子里面都有，建议可以仔细看完再问。

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