24小时热门版块排行榜    

查看: 2302  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

些小

铜虫 (初入文坛)

[求助] 过渡态QST2总是失败,请帮我看看

用QST2搜索过渡态总是出错,请帮我看看。总是提示l9999的错误,在out文件中有这样一句话:Item               Value     Threshold  Converged?
Maximum Force            0.000901     0.000450     NO
RMS     Force            0.000242     0.000300     YES
Maximum Displacement     0.353423     0.001800     NO
RMS     Displacement     0.091115     0.001200     NO
Predicted change in Energy=-2.242666D-04
Optimization stopped.
    -- Number of steps exceeded,  NStep= 100
    -- Flag reset to prevent archiving.
用gaussian view打开out文件,两个反应物变得更远了,不知道是不是输入文件的结构有问题,请帮我看看。新人没有金币可以悬赏,对不起。
输入文件如下:
%Chk=Untitled-1.chk
%nproc=3
# opt=(calcfc,qst2,maxcycle=5000,maxstep=500) freq b3lyp geom=connectivity
scf=maxcycle=1024 gen pseudo=read test

reactants

0 1
C
C                  1              B1
Br                 2              B2    1              A1
Pd                 3              B3    2              A2    1              D1
P                  4              B4    3              A3    2              D2
P                  4              B5    3              A4    2              D3
H                  1              B6    2              A5    3              D4
H                  1              B7    2              A6    3              D5
H                  2              B8    1              A7    3              D6
H                  5              B9    4              A8    3              D7
H                  5             B10    4              A9    3              D8
H                  5             B11    4             A10    3              D9
H                  6             B12    4             A11    3             D10
H                  6             B13    4             A12    3             D11
H                  6             B14    4             A13    3             D12

   B1             1.32660010
   B2             1.90792296
   B3             4.50098403
   B4             2.30438200
   B5             2.30438200
   B6             1.08498400
   B7             1.08711990
   B8             1.08352758
   B9             1.41873393
   B10            1.41873388
   B11            1.41873392
   B12            1.41873392
   B13            1.41873388
   B14            1.41873393
   A1           123.05128037
   A2            55.99949871
   A3            63.43782277
   A4           116.56217723
   A5           122.68208037
   A6           119.63246425
   A7           124.52962183
   A8           121.06541806
   A9           121.06541929
   A10          121.06541815
   A11          121.06541815
   A12          121.06541929
   A13          121.06541806
   D1            80.55375590
   D2           125.23835973
   D3           -54.76164027
   D4            -0.00448079
   D5           179.98821515
   D6           179.99102849
   D7            35.02710965
   D8           155.02710828
   D9           -84.97289346
   D10          144.90848711
   D11           24.90848886
   D12          -95.09150978

1 2 2.0 7 1.0 8 1.0
2 3 1.0 9 1.0
3
4 6 1.0 5 1.0
5 12 1.0 10 1.0 11 1.0
6 13 1.0 15 1.0 14 1.0
7
8
9
10
11
12
13
14
15

Title Card Required

0 1
C
C                  1              B1
Br                 2              B2    1              A1
Pd                 2              B3    1              A2    3              D1
P                  4              B4    2              A3    1              D2
P                  4              B5    2              A4    1              D3
H                  1              B6    2              A5    4              D4
H                  1              B7    2              A6    4              D5
H                  2              B8    1              A7    4              D6
H                  5              B9    4              A8    2              D7
H                  5             B10    4              A9    2              D8
H                  5             B11    4             A10    2              D9
H                  6             B12    4             A11    2             D10
H                  6             B13    4             A12    2             D11
H                  6             B14    4             A13    2             D12

   B1             1.33276933
   B2             3.18020961
   B3             2.02951742
   B4             2.45539153
   B5             2.29749429
   B6             1.08859253
   B7             1.09130219
   B8             1.08983443
   B9             1.40958445
   B10            1.40951409
   B11            1.42199436
   B12            1.41738686
   B13            1.40993619
   B14            1.41028783
   A1           125.41452910
   A2           126.48940219
   A3           171.94989407
   A4            83.39239473
   A5           123.20279898
   A6           120.98034046
   A7           121.55974105
   A8           115.44065467
   A9           115.19693220
   A10          127.10412369
   A11          121.99745370
   A12          117.48939346
   A13          117.05247664
   D1            65.74766237
   D2          -112.67267069
   D3            71.63740640
   D4            -1.91950175
   D5           175.25733871
   D6          -176.31591604
   D7           -54.32107456
   D8            59.57135703
   D9          -177.49552219
   D10          177.45958098
   D11           56.22808693
   D12          -61.10612297

1 2 2.0 7 1.0 8 1.0
2 4 1.0 9 1.0
3 4 1.0
4 5 1.0 6 1.0
5 10 1.0 11 1.0 12 1.0
6 13 1.0 14 1.0 15 1.0
7
8
9
10
11
12
13
14
15

P H C Br 0
6-31G*
****
Pd 0
LANL2DZ
****

Pd 0
LANL2DZ
回复此楼

» 收录本帖的淘帖专辑推荐

Gaussian

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldzh0531

木虫 (小有名气)

【答案】应助回帖

引用回帖:
6楼: Originally posted by 些小 at 2012-05-23 19:17:49:
还想请教一下IRC的问题。
TS计算完成后,IRC的输入文件是怎么生成的???
看论坛说要找到.chk文件,那然后呢???
就是说比如我算了几何优化,想再算freq,就会把几何优化算完的out文件用view打开,再存为GJP文

我习惯自己手动修改路径部分!给你一个我的输入文件吧:# irc=(rcfc,maxpoints=N,stepsize=N,maxcyc=200) b3lyp/6-31g(d,p) geom=connectivity
7楼2012-05-23 22:41:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

chembetsey

木虫 (小有名气)

★ ★
gmy1990: 金币+2 2012-05-22 21:57:51
我的理解是,QST2找过渡态的算法不够“聪明”。如果你能猜到过渡态的结构,最好还是自己优化过渡态(用ts=opt)。
2楼2012-05-22 01:06:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

些小

铜虫 (初入文坛)

我也用TS直接优惠了过渡态,但结果出来确实有有且只有一个虚频,但从结构上看和产物相差不大,所以不知道是否可用。才想用QST2再找一下,但是就总是出上述错误。
3楼2012-05-22 13:11:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldzh0531

木虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
gmy1990: 金币+1 2012-05-23 15:53:20
个人觉得qst2就不要加calcfc了,那个是先计算过渡态的初始猜测的力常数,以方便为过渡态的优化找一个大致方向的。你用qst2了还加这个干什么。另外,循环的步长没必要那么大,200就足够了,再多也没多大用。
4楼2012-05-22 23:30:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见