应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 过渡态计算出错
61/1返回列表
查看: 1611  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

rock0927

铁虫 (正式写手)

[求助] 过渡态计算出错

用QST2计算过渡态，算到最后出现的错误是：
Error termination request processed by link 9999.
请问这是什么原因造成的

输入的文件：
%chk=qst2_formts
#T RHF/6-31G(d) Opt=QST2 Freq Test

Reactant

1 1
C
H                  1              B1
H                  1              B2    2              A1
C                  1              B3    3              A2    2              D1
C                  4              B4    1              A3    3              D2
H                  5              B5    4              A4    1              D3
H                  5              B6    4              A5    1              D4
H                  5              B7    4              A6    1              D5
H                  4              B8    1              A7    5              D6
C                  1              B9    4              A8    5              D7
H                 10             B10    1              A9    4              D8
H                 10             B11    1             A10    4              D9
H                 10             B12    1             A11    4             D10
C                 10             B13    1             A12    4             D11
H                 14             B14   10             A13    1             D12
C                 14             B15   10             A14    1             D13
H                 16             B16   14             A15   10             D14
H                 16             B17   14             A16   10             D15
H                 16             B18   14             A17   10             D16

   B1             1.07756154
   B2             1.07666980
   B3             1.32422058
   B4             1.50258845
   B5             1.08675771
   B6             1.08551966
   B7             1.08386303
   B8             1.07924633
   B9             3.80620601
   B10            1.11030507
   B11            1.08081693
   B12            1.08258036
   B13            1.44521114
   B14            1.07969671
   B15            1.45798967
   B16            1.08011023
   B17            1.08454965
   B18            1.09898194
   A1           116.66300961
   A2           121.45368704
   A3           125.36474644
   A4           111.29458745
   A5           110.11905479
   A6           111.70967743
   A7           118.69056636
   A8            82.53613015
   A9             7.53216862
   A10          106.47172236
   A11          112.71371115
   A12           95.42223037
   A13          117.58700260
   A14          125.23299795
   A15          113.72128073
   A16          111.78104926
   A17          103.80909975
   D1          -175.42774064
   D2           175.87654011
   D3           125.10435457
   D4          -116.74457681
   D5             3.42220268
   D6          -177.48584992
   D7           -91.36634534
   D8           -13.66377278
   D9            86.06705187
   D10          -38.82456258
   D11         -156.79810249
   D12           86.71042084
   D13          -91.89596580
   D14           12.14033946
   D15          140.93999071
   D16         -105.72199284

Product

1 1
C
H                  1              B1
H                  1              B2    2              A1
C                  1              B3    2              A2    3              D1
C                  4              B4    1              A3    2              D2
H                  5              B5    4              A4    1              D3
H                  5              B6    4              A5    1              D4
H                  5              B7    4              A6    1              D5
H                  4              B8    1              A7    5              D6
C                  1              B9    4              A8    5              D7
H                 10             B10    1              A9    4              D8
H                 10             B11    1             A10    4              D9
H                 10             B12    1             A11    4             D10
C                 10             B13    1             A12    4             D11
H                 14             B14   10             A13    1             D12
C                 14             B15   10             A14    1             D13
H                 16             B16   14             A15   10             D14
H                 16             B17   14             A16   10             D15
H                 16             B18   14             A17   10             D16

   B1             1.08122565
   B2             1.08518309
   B3             1.44501077
   B4             1.45881961
   B5             1.08409858
   B6             1.09921072
   B7             1.07978484
   B8             1.08000859
   B9             2.59976424
   B10            1.08667558
   B11            1.08532914
   B12            1.08306312
   B13            1.52970077
   B14            1.08480500
   B15            1.53066904
   B16            1.08449522
   B17            1.08351778
   B18            1.08375726
   A1           111.22720955
   A2           112.54086380
   A3           126.19076504
   A4           112.20675285
   A5           103.39575381
   A6           113.73650450
   A7           116.97757416
   A8            91.11356688
   A9            95.75175400
   A10           90.43891554
   A11          144.15696757
   A12           34.81848455
   A13          109.63167941
   A14          111.81848154
   A15          111.23552095
   A16          111.78314431
   A17          109.16594207
   D1           123.30120193
   D2            13.62126913
   D3           143.32944868
   D4          -103.23240381
   D5            14.09736950
   D6          -179.96976201
   D7           -77.82214346
   D8             1.72139477
   D9           110.19015255
   D10         -128.01389820
   D11         -121.06008242
   D12          119.22860908
   D13         -120.20406450
   D14           60.52480344
   D15         -177.86337370
   D16          -58.51363537
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-06-17 20:56:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤狼

金虫 (正式写手)

独行者

【答案】应助回帖

★ ★
rock0927(金币+2): 谢谢 2011-06-17 23:37:02
gmy1990(金币+2): 2011-06-18 10:07:05
rock0927(金币+2): 谢谢 2011-06-18 14:05:52
优化步数不够 加上opt=300试试~~
赞一下(1人) 回复此楼
轮回的过往，注定的未来；独孤的行者，孤独的路途……
2楼2011-06-17 22:42:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rock0927

铁虫 (正式写手)
引用回帖:
Originally posted by 独孤狼 at 2011-06-17 22:42:34:
优化步数不够 加上opt=300试试~~

未加限制之前算的最大优化步数是114，加上Maxcycle=300之后为什么也只算了114步，最后出现的错误是一样的，link9999。
最后的收敛结果是：
Maximum Force            0.000227     0.000450     YES
RMS     Force            0.000042     0.000300     YES
Maximum Displacement     0.007633     0.001800     NO
RMS     Displacement     0.002613     0.001200     NO
赞一下 回复此楼
3楼2011-06-18 14:09:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤狼

金虫 (正式写手)

独行者

【答案】应助回帖

rock0927(金币+2): 谢谢 2011-06-20 08:49:51
引用回帖:
Originally posted by rock0927 at 2011-06-18 14:09:58:
未加限制之前算的最大优化步数是114，加上Maxcycle=300之后为什么也只算了114步，最后出现的错误是一样的，link9999。
最后的收敛结果是：
Maximum Force            0.000227     0.000450     YES
RMS   ...

那读chk文件接着算，你这个应该快收敛了。
赞一下 回复此楼
轮回的过往，注定的未来；独孤的行者，孤独的路途……
4楼2011-06-18 20:24:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤狼

金虫 (正式写手)

独行者

【答案】应助回帖

rock0927(金币+1): 谢谢 2011-06-20 08:53:41
引用回帖:
Originally posted by rock0927 at 2011-06-18 14:09:58:
未加限制之前算的最大优化步数是114，加上Maxcycle=300之后为什么也只算了114步，最后出现的错误是一样的，link9999。
最后的收敛结果是：
Maximum Force            0.000227     0.000450     YES
RMS   ...

最大优化步数的最大上限是高斯根据你的原子的数量决定的。但是高斯默认的上限一般不是最大上限，通过maxcyc指定的上限是不能超过最大上限的。
赞一下 回复此楼
轮回的过往，注定的未来；独孤的行者，孤独的路途……
5楼2011-06-18 20:26:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aimonkey

金虫 (正式写手)
看了你画的两个结构，它们之间的差异太大，所以在最大优化步数内没能优化成功。这种情况可能可以通过加大优化步数（opt=300）来解决，但最关键的还是你初始构型给得不合理。建议你自己猜测过渡态构型，将C1-C14键稍拉远（1.8-2埃左右吧），其他变化的键长也相应调节到反应物和产物的键长之间的距离，用TS的方法计算过渡态结构.QST2方法是懒人方法，对一些简单的变化的过渡态（如迁移时形成的三元环过渡态）比较适合。
赞一下 回复此楼
6楼2011-06-19 12:01:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 rock0927 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭