QST2方法寻找过渡态出现l9999.exe

ct_2015 · 发表于 Post on 2015-8-1 19:08:10

RT，我正在使用QST2方法寻找过渡态，现出现了l9999.exe错误。输入文件如下：#p b3lyp/6-31g* fopt(qst2,maxcyc=800) iop(5/13=1) iop(2/16=1) iop(1/11=1) freq

输出文件最后几行如下：
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=       -6347.5856 YYYY=       -1327.2480 ZZZZ=          -514.2104 XXXY=             2.7877
XXXZ=          40.4207 YYYX=          -53.7467 YYYZ=          17.9651 ZZZX=          26.9504
ZZZY=          -0.9658 XXYY=       -1255.3903 XXZZ=       -1291.3123 YYZZ=          -320.5586
XXYZ=             2.5513 YYXZ=          13.3472 ZZXY=          -5.7214
N-N= 1.516128350120D+03 E-N=-5.068417167781D+03  KE= 8.693373371463D+02
No NMR shielding tensors so no spin-rotation constants.
Leave Link  601 at Sat Aug  1 09:36:25 2015, MaxMem=  786432000 cpu:    1.2
(Enter /home/yysheng/software/g09/l9999.exe)

NO MATTER HOW THIN YOU SLICE IT, IT'S STILL BALONEY

                                 -- ALFRED E. SMITH
Error termination request processed by link 9999.
Error termination via Lnk1e in /home/yysheng/software/g09/l9999.exe at Sat Aug  1 09:36:25 2015.
Job cpu time:  0 days  1 hours 13 minutes 44.2 seconds.
File lengths (MBytes):  RWF= 132 Int=    0 D2E=    0 Chk=    10 Scr=    2
"1ts2C_qst2.log" 9012L, 566800C                                                                                           9012,2       Bot

之前以为是maxcyc值设置小了，输入文件里面的maxcyc的值一开始设置为150的，然后修改为了500和800，还是出现同样了的问题。我是菜鸟，没招了，请教高手帮忙解答，谢谢！

smutao · 发表于 Post on 2015-8-1 22:55:13

試試 opt=noeigentest

sobereva · 发表于 Post on 2015-8-2 01:21:01

几乎一定是没收敛。增大maxcyc是白费功夫，只会浪费计算量。

解决方法此文有讨论，虽然讨论的是几何优化的，但大多数对于过渡态搜索同样管用
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

这些关键词别写：iop(5/13=1) iop(2/16=1) iop(1/11=1) maxcyc=800
1/11=1等价于noeigen，凡是有现成关键词的都不要以iop来写

ct_2015 · 发表于 Post on 2015-8-3 20:09:51

sobereva 发表于 2015-8-2 01:21
几乎一定是没收敛。增大maxcyc是白费功夫，只会浪费计算量。

解决方法此文有讨论，虽然讨论的是几何优化 ...

谢谢Sob老师给推荐的文章，可是作为一个菜鸟，没怎么看懂，只能根据文章最后一段的方法解决问题：
具体怎么解决优化不收敛得看实际情况，特别是要用gview检查优化过程的轨迹，是否最后开始震荡了。虽然上面讨论了很多，但从实用性出发，对于初学者笔者推荐按照以下顺序去尝试解决震荡问题，如下：
1 对于DFT，特别是明尼苏达系列泛函（如M062X、M06、M06L），首先加上int=ultrafine再试。
2 如果体系小，尝试opt=calcall。如果体系大算得慢，此关键词会令计算量进一步激增，实在没辙再考虑这个。
3 尝试opt=gdiis。
4 尝试opt(gdiis,maxstep=x,notrust)，x为3~5。
5 如果上述做法都不灵，干脆改成loose收敛限，但如果之后做振动分析容易有虚频。

这五条都尝试了，尝试3和4建议运行不久就终止了，1,2和5可以正常运行可是结果跟之前的一样，如输入文件：#p b3lyp/6-31g* fopt(qst2,verytight) freq

输出文件最后几行：
Quadrupole moment (field-independent basis, Debye-Ang):
XX=          -100.3184 YY=          -109.6487 ZZ=          -129.8615
XY=          -6.1587 XZ=             4.1677 YZ=             2.4977
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
XX=          12.9578 YY=             3.6275 ZZ=          -16.5853
XY=          -6.1587 XZ=             4.1677 YZ=             2.4977
Octapole moment (field-independent basis, Debye-Ang**2):
  XXX=          -43.9367  YYY=          23.5706  ZZZ=          -19.5324  XYY=          10.1315
  XXY=          -2.7735  XXZ=          -47.0646  XZZ=             0.2170  YZZ=          -3.3515
  YYZ=          -8.6280  XYZ=             0.9601
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=       -6347.5856 YYYY=       -1327.2480 ZZZZ=          -514.2104 XXXY=             2.7877
XXXZ=          40.4207 YYYX=          -53.7467 YYYZ=          17.9651 ZZZX=          26.9504
ZZZY=          -0.9658 XXYY=       -1255.3903 XXZZ=       -1291.3123 YYZZ=          -320.5586
XXYZ=             2.5513 YYXZ=          13.3472 ZZXY=          -5.7214
N-N= 1.516128350121D+03 E-N=-5.068417167782D+03  KE= 8.693373371463D+02
No NMR shielding tensors so no spin-rotation constants.
Leave Link  601 at Tue Aug  4 00:48:52 2015, MaxMem=  786432000 cpu:    1.7
(Enter /home/yysheng/software/g09/l9999.exe)

THE WORLD OF CHEMICAL REACTIONS IS LIKE A STAGE,
ON WHICH SCENE AFTER SCENE IS CEASELESSLY PLAYED.
THE ACTORS ON IT ARE THE ELEMENTS.
-- CLEMENS WINKLER BER. 30,13(1897)
(DISCOVERER OF GERMANIUM, FEB 6, 1886)
Error termination request processed by link 9999.
Error termination via Lnk1e in /home/yysheng/software/g09/l9999.exe at Tue Aug  4 00:48:52 2015.
Job cpu time:  0 days  1 hours 54 minutes 37.5 seconds.
Faile lengths (MBytes):  RWF= 132 Int=    0 D2E=    0 Chk=    10 Scr=    2

Sob老师看看出现了什么问题，谢谢！

liyuanhe211 · 发表于 Post on 2015-8-3 21:49:11

ct_2015 发表于 2015-8-3 20:09
谢谢Sob老师给推荐的文章，可是作为一个菜鸟，没怎么看懂，只能根据文章最后一段的方法解决问题：
具体 ...

输出文件往上找，在表Non-Optimized Param之上找到形如：
Optimization stopped.
-- Wrong number of Negative eigenvalues: Desired= 1 Actual= 0 # Or something else
-- Flag reset to prevent archiving.
的行s，贴出来。

liyuanhe211 · 发表于 Post on 2015-8-3 21:50:06

ct_2015 发表于 2015-8-3 20:09
谢谢Sob老师给推荐的文章，可是作为一个菜鸟，没怎么看懂，只能根据文章最后一段的方法解决问题：
具体 ...

以及考虑你的反应是不是基元反应，反应物和产物是否已经优化，结构是否合理

sobereva · 发表于 Post on 2015-8-3 23:07:18

ct_2015 发表于 2015-8-3 20:09
谢谢Sob老师给推荐的文章，可是作为一个菜鸟，没怎么看懂，只能根据文章最后一段的方法解决问题：
具体 ...

不要停留在：1 出错 2 没出错这么低级的判断上。
要从原理出发，搞明白每种帮助收敛的做法的原理，而且必须要观看优化过程中结构变化的轨迹、能量/受力的变化曲线，根据实际情况判断怎么处理最合适。

没有万能的解决不收敛的方法！

ct_2015 · 发表于 Post on 2015-8-4 14:05:55

liyuanhe211 发表于 2015-8-3 21:49
输出文件往上找，在表Non-Optimized Param之上找到形如：
Optimization stopped.
-- Wrong number ...

表Non-Optimized Param之上找到如下：
199       D99       0.00318  -0.01522 0.000001000.00000
200       D100    -0.00119  -0.00388 0.000001000.00000
201       D101    0.00298  -0.00971 0.000001000.00000
202       D102    -0.00139 0.00163 0.000001000.00000
RFO could not converge Lambda in  999 iterations.
Linear search not attempted -- option 19 set.
Iteration  1 RMS(Cart)=  0.00000000 RMS(Int)=  0.00000000
Variable    Old X -DE/DX Delta X Delta X Delta X    New X
                              (Linear) (Quad) (Total)
R1       2.65563  -0.00052 0.00000 0.00000 0.00000 2.65563
R2       2.63161 0.00066 0.00000 0.00000 0.00000 2.63161
R3       2.05123 0.00009 0.00000 0.00000 0.00000 2.05123
R4       2.63630 0.00068 0.00000 0.00000 0.00000 2.63630
R5       2.05383  -0.00007 0.00000 0.00000 0.00000 2.05383
R6       2.64106  -0.00121 0.00000 0.00000 0.00000 2.64106
R7       2.05049 0.00007 0.00000 0.00000 0.00000 2.05049
R8       2.66132  -0.00097 0.00000 0.00000 0.00000 2.66132
R9       2.63967 0.00470 0.00000 0.00000 0.00000 2.63967
R10       2.63501  -0.00031 0.00000 0.00000 0.00000 2.63501
R11       2.77485 0.00271 0.00000 0.00000 0.00000 2.77485
R12       2.05248  -0.00005 0.00000 0.00000 0.00000 2.05248
R13       2.69215  -0.00886 0.00000 0.00000 0.00000 2.69215
R14       2.62948  -0.00177 0.00000 0.00000 0.00000 2.62948
R15       2.85268  -0.00106 0.00000 0.00000 0.00000 2.85268
R16       2.70983  -0.00155 0.00000 0.00000 0.00000 2.70983
R17       3.93457  -0.03642 0.00000 0.00000 0.00000 3.93457
R18       2.30962 0.00284 0.00000 0.00000 0.00000 2.30962
R19       2.56361 0.00220 0.00000 0.00000 0.00000 2.56361
R20       2.04844  -0.00029 0.00000 0.00000 0.00000 2.04844
R21       2.72222  -0.00262 0.00000 0.00000 0.00000 2.72222
R22       2.04936 0.00002 0.00000 0.00000 0.00000 2.04936
R23       2.63365 0.00235 0.00000 0.00000 0.00000 2.63365
R24       2.72312  -0.00009 0.00000 0.00000 0.00000 2.72312
R25       2.59991 0.00023 0.00000 0.00000 0.00000 2.59991
R26       2.69454 0.00084 0.00000 0.00000 0.00000 2.69454
R27       2.65479 0.00012 0.00000 0.00000 0.00000 2.65479
R28       2.61155  -0.00003 0.00000 0.00000 0.00000 2.61155
R29       2.64229 0.00015 0.00000 0.00000 0.00000 2.64229
R30       1.91121 0.00012 0.00000 0.00000 0.00000 1.91121
R31       2.62380  -0.00010 0.00000 0.00000 0.00000 2.62380
R32       2.05381  -0.00002 0.00000 0.00000 0.00000 2.05381
R33       2.66495 0.00022 0.00000 0.00000 0.00000 2.66495
R34       2.05286  -0.00003 0.00000 0.00000 0.00000 2.05286
R35       2.62687  -0.00008 0.00000 0.00000 0.00000 2.62687
R36       2.05302 0.00000 0.00000 0.00000 0.00000 2.05302
R37       2.05299  -0.00003 0.00000 0.00000 0.00000 2.05299
A1       2.09326 0.00031 0.00000 0.00000 0.00000 2.09326
A2       2.08908  -0.00012 0.00000 0.00000 0.00000 2.08908
A3       2.10081  -0.00018 0.00000 0.00000 0.00000 2.10081
A4       2.13090  -0.00048 0.00000 0.00000 0.00000 2.13090
A5       2.07997 0.00027 0.00000 0.00000 0.00000 2.07997
A6       2.07231 0.00020 0.00000 0.00000 0.00000 2.07231
A7       2.05208  -0.00035 0.00000 0.00000 0.00000 2.05208
A8       2.10630  -0.00003 0.00000 0.00000 0.00000 2.10630
A9       2.12478 0.00038 0.00000 0.00000 0.00000 2.12478
A10       2.10594 0.00121 0.00000 0.00000 0.00000 2.10594
A11       2.25044 0.00037 0.00000 0.00000 0.00000 2.25044
A12       1.92616  -0.00155 0.00000 0.00000 0.00000 1.92616
A13       2.11369  -0.00053 0.00000 0.00000 0.00000 2.11369
A14       1.88516  -0.00077 0.00000 0.00000 0.00000 1.88516
A15       2.27731 0.00150 0.00000 0.00000 0.00000 2.27731
A16       2.07010  -0.00016 0.00000 0.00000 0.00000 2.07010
A17       2.12655  -0.00004 0.00000 0.00000 0.00000 2.12655
A18       2.08654 0.00020 0.00000 0.00000 0.00000 2.08654
A19       1.90448 0.00146 0.00000 0.00000 0.00000 1.90448
A20       2.23411  -0.00251 0.00000 0.00000 0.00000 2.23411
A21       2.12345 0.00130 0.00000 0.00000 0.00000 2.12345
A22       1.86430 0.00177 0.00000 0.00000 0.00000 1.86430
A23       1.98538 0.00101 0.00000 0.00000 0.00000 1.98538
A24       1.82309  -0.00075 0.00000 0.00000 0.00000 1.82309
A25       2.14124 0.00108 0.00000 0.00000 0.00000 2.14124
A26       1.72009  -0.00485 0.00000 0.00000 0.00000 1.72009
A27       1.88714 0.00052 0.00000 0.00000 0.00000 1.88714
A28       1.82258  -0.00003 0.00000 0.00000 0.00000 1.82258
A29       2.27187  -0.00140 0.00000 0.00000 0.00000 2.27187
A30       2.18553 0.00186 0.00000 0.00000 0.00000 2.18553
A31       2.10565 0.00022 0.00000 0.00000 0.00000 2.10565
A32       2.02280 0.00002 0.00000 0.00000 0.00000 2.02280
A33       2.15325  -0.00020 0.00000 0.00000 0.00000 2.15325
A34       2.07608  -0.00071 0.00000 0.00000 0.00000 2.07608
A35       2.08713 0.00016 0.00000 0.00000 0.00000 2.08713
A36       2.11866 0.00062 0.00000 0.00000 0.00000 2.11866
A37       2.07178  -0.00138 0.00000 0.00000 0.00000 2.07178
A38       2.33123 0.00162 0.00000 0.00000 0.00000 2.33123
A39       1.85565  -0.00007 0.00000 0.00000 0.00000 1.85565
A40       2.12774 0.00039 0.00000 0.00000 0.00000 2.12774
A41       2.18042 0.00148 0.00000 0.00000 0.00000 2.18042
A42       1.91857  -0.00016 0.00000 0.00000 0.00000 1.91857
A43       1.86070 0.00002 0.00000 0.00000 0.00000 1.86070
A44       2.34466  -0.00008 0.00000 0.00000 0.00000 2.34466
A45       2.07743 0.00009 0.00000 0.00000 0.00000 2.07743
A46       1.88627 0.00041 0.00000 0.00000 0.00000 1.88627
A47       2.12937  -0.00022 0.00000 0.00000 0.00000 2.12937
A48       2.26755  -0.00019 0.00000 0.00000 0.00000 2.26755
A49       1.89850 0.00018 0.00000 0.00000 0.00000 1.89850
A50       2.14999  -0.00035 0.00000 0.00000 0.00000 2.14999
A51       2.23466 0.00016 0.00000 0.00000 0.00000 2.23466
A52       2.07764  -0.00007 0.00000 0.00000 0.00000 2.07764
A53       2.10303  -0.00002 0.00000 0.00000 0.00000 2.10303
A54       2.10251 0.00009 0.00000 0.00000 0.00000 2.10251
A55       2.11205 0.00003 0.00000 0.00000 0.00000 2.11205
A56       2.08866 0.00003 0.00000 0.00000 0.00000 2.08866
A57       2.08247  -0.00006 0.00000 0.00000 0.00000 2.08247
A58       2.11812 0.00022 0.00000 0.00000 0.00000 2.11812
A59       2.08346  -0.00011 0.00000 0.00000 0.00000 2.08346
A60       2.08161  -0.00011 0.00000 0.00000 0.00000 2.08161
A61       2.05168  -0.00004 0.00000 0.00000 0.00000 2.05168
A62       2.11516 0.00002 0.00000 0.00000 0.00000 2.11516
A63       2.11634 0.00002 0.00000 0.00000 0.00000 2.11634
D1    -0.00573 0.00014 0.00000 0.00000 0.00000  -0.00573
D2       3.13080 0.00007 0.00000 0.00000 0.00000 3.13080
D3    -3.13932  -0.00013 0.00000 0.00000 0.00000  -3.13932
D4    -0.00279  -0.00020 0.00000 0.00000 0.00000  -0.00279
D5       0.02500  -0.00020 0.00000 0.00000 0.00000 0.02500
D6    -3.11589  -0.00067 0.00000 0.00000 0.00000  -3.11589
D7    -3.12465 0.00006 0.00000 0.00000 0.00000  -3.12465
D8       0.01764  -0.00041 0.00000 0.00000 0.00000 0.01764
D9    -0.00650 0.00016 0.00000 0.00000 0.00000  -0.00650
D10       3.12778  -0.00029 0.00000 0.00000 0.00000 3.12778
D11       3.14014 0.00023 0.00000 0.00000 0.00000 3.14014
D12    -0.00877  -0.00022 0.00000 0.00000 0.00000  -0.00877
D13    -0.00086  -0.00039 0.00000 0.00000 0.00000  -0.00086
D14    -3.09985  -0.00145 0.00000 0.00000 0.00000  -3.09985
D15    -3.13505 0.00007 0.00000 0.00000 0.00000  -3.13505
D16       0.04914  -0.00099 0.00000 0.00000 0.00000 0.04914
D17       0.02072 0.00034 0.00000 0.00000 0.00000 0.02072
D18       3.04785 0.00204 0.00000 0.00000 0.00000 3.04785
D19       3.12698 0.00124 0.00000 0.00000 0.00000 3.12698
D20    -0.12909 0.00295 0.00000 0.00000 0.00000  -0.12909
D21       3.10824 0.00061 0.00000 0.00000 0.00000 3.10824
D22    -0.25291 0.00220 0.00000 0.00000 0.00000  -0.25291
D23       0.00573  -0.00041 0.00000 0.00000 0.00000 0.00573
D24       2.92776 0.00119 0.00000 0.00000 0.00000 2.92776
D25    -0.03257  -0.00001 0.00000 0.00000 0.00000  -0.03257
D26       3.10834 0.00045 0.00000 0.00000 0.00000 3.10834
D27    -3.03087  -0.00194 0.00000 0.00000 0.00000  -3.03087
D28       0.11004  -0.00148 0.00000 0.00000 0.00000 0.11004
D29       0.19099  -0.00404 0.00000 0.00000 0.00000 0.19099
D30    -3.03471 0.00162 0.00000 0.00000 0.00000  -3.03471
D31    -3.07953  -0.00232 0.00000 0.00000 0.00000  -3.07953
D32    -0.02205 0.00334 0.00000 0.00000 0.00000  -0.02205
D33       0.11582  -0.00195 0.00000 0.00000 0.00000 0.11582
D34       2.53370 0.00237 0.00000 0.00000 0.00000 2.53370
D35    -1.69567 0.00306 0.00000 0.00000 0.00000  -1.69567
D36    -2.82275  -0.00285 0.00000 0.00000 0.00000  -2.82275
D37    -0.40487 0.00147 0.00000 0.00000 0.00000  -0.40487
D38       1.64895 0.00216 0.00000 0.00000 0.00000 1.64895
D39    -2.66718  -0.00184 0.00000 0.00000 0.00000  -2.66718
D40       0.53105  -0.00259 0.00000 0.00000 0.00000 0.53105
D41       0.23021  -0.00014 0.00000 0.00000 0.00000 0.23021
D42    -2.85475  -0.00089 0.00000 0.00000 0.00000  -2.85475
D43    -0.18637 0.00353 0.00000 0.00000 0.00000  -0.18637
D44       3.03388  -0.00157 0.00000 0.00000 0.00000 3.03388
D45    -2.52403  -0.00129 0.00000 0.00000 0.00000  -2.52403
D46       0.69623  -0.00638 0.00000 0.00000 0.00000 0.69623
D47       1.69875 0.00140 0.00000 0.00000 0.00000 1.69875
D48    -1.36419  -0.00370 0.00000 0.00000 0.00000  -1.36419
D49       0.26616  -0.00091 0.00000 0.00000 0.00000 0.26616
D50    -2.48451  -0.00657 0.00000 0.00000 0.00000  -2.48451
D51       2.55499 0.00450 0.00000 0.00000 0.00000 2.55499
D52    -0.19568  -0.00116 0.00000 0.00000 0.00000  -0.19568
D53    -1.75028  -0.00091 0.00000 0.00000 0.00000  -1.75028
D54       1.78224  -0.00657 0.00000 0.00000 0.00000 1.78224
D55       0.10503  -0.00151 0.00000 0.00000 0.00000 0.10503
D56    -3.09102  -0.00012 0.00000 0.00000 0.00000  -3.09102
D57    -3.09758  -0.00069 0.00000 0.00000 0.00000  -3.09758
D58    -0.01044 0.00070 0.00000 0.00000 0.00000  -0.01044
D59    -0.23570 0.00264 0.00000 0.00000 0.00000  -0.23570
D60       2.63342 0.00349 0.00000 0.00000 0.00000 2.63342
D61       2.96137 0.00123 0.00000 0.00000 0.00000 2.96137
D62    -0.45270 0.00209 0.00000 0.00000 0.00000  -0.45270
D63       0.03882  -0.00071 0.00000 0.00000 0.00000 0.03882
D64       2.84168 0.00451 0.00000 0.00000 0.00000 2.84168
D65    -2.89817  -0.00167 0.00000 0.00000 0.00000  -2.89817
D66    -0.09531 0.00355 0.00000 0.00000 0.00000  -0.09531
D67    -2.83094  -0.00298 0.00000 0.00000 0.00000  -2.83094
D68       0.27897  -0.00174 0.00000 0.00000 0.00000 0.27897
D69       0.06215  -0.00240 0.00000 0.00000 0.00000 0.06215
D70    -3.11113  -0.00115 0.00000 0.00000 0.00000  -3.11113
D71       2.88285 0.00176 0.00000 0.00000 0.00000 2.88285
D72    -0.24981 0.00220 0.00000 0.00000 0.00000  -0.24981
D73       0.09249  -0.00338 0.00000 0.00000 0.00000 0.09249
D74    -3.04017  -0.00295 0.00000 0.00000 0.00000  -3.04017
D75    -0.00809 0.00041 0.00000 0.00000 0.00000  -0.00809
D76       3.13118 0.00075 0.00000 0.00000 0.00000 3.13118
D77    -3.12377  -0.00061 0.00000 0.00000 0.00000  -3.12377
D78       0.01550  -0.00026 0.00000 0.00000 0.00000 0.01550
D79    -3.11809  -0.00085 0.00000 0.00000 0.00000  -3.11809
D80       0.02624  -0.00085 0.00000 0.00000 0.00000 0.02624
D81    -0.01122 0.00051 0.00000 0.00000 0.00000  -0.01122
D82       3.13311 0.00051 0.00000 0.00000 0.00000 3.13311
D83    -0.05051 0.00174 0.00000 0.00000 0.00000  -0.05051
D84       3.08159 0.00127 0.00000 0.00000 0.00000 3.08159
D85       3.09365 0.00136 0.00000 0.00000 0.00000 3.09365
D54       1.78224  -0.00657 0.00000 0.00000 0.00000 1.78224
D55       0.10503  -0.00151 0.00000 0.00000 0.00000 0.10503
D56    -3.09102  -0.00012 0.00000 0.00000 0.00000  -3.09102
D57    -3.09758  -0.00069 0.00000 0.00000 0.00000  -3.09758
D58    -0.01044 0.00070 0.00000 0.00000 0.00000  -0.01044
D59    -0.23570 0.00264 0.00000 0.00000 0.00000  -0.23570
D60       2.63342 0.00349 0.00000 0.00000 0.00000 2.63342
D61       2.96137 0.00123 0.00000 0.00000 0.00000 2.96137
D62    -0.45270 0.00209 0.00000 0.00000 0.00000  -0.45270
D63       0.03882  -0.00071 0.00000 0.00000 0.00000 0.03882
D64       2.84168 0.00451 0.00000 0.00000 0.00000 2.84168
D65    -2.89817  -0.00167 0.00000 0.00000 0.00000  -2.89817
D66    -0.09531 0.00355 0.00000 0.00000 0.00000  -0.09531
D67    -2.83094  -0.00298 0.00000 0.00000 0.00000  -2.83094
D68       0.27897  -0.00174 0.00000 0.00000 0.00000 0.27897
D69       0.06215  -0.00240 0.00000 0.00000 0.00000 0.06215
D70    -3.11113  -0.00115 0.00000 0.00000 0.00000  -3.11113
D71       2.88285 0.00176 0.00000 0.00000 0.00000 2.88285
D72    -0.24981 0.00220 0.00000 0.00000 0.00000  -0.24981
D73       0.09249  -0.00338 0.00000 0.00000 0.00000 0.09249
D74    -3.04017  -0.00295 0.00000 0.00000 0.00000  -3.04017
D75    -0.00809 0.00041 0.00000 0.00000 0.00000  -0.00809
D76       3.13118 0.00075 0.00000 0.00000 0.00000 3.13118
D77    -3.12377  -0.00061 0.00000 0.00000 0.00000  -3.12377
D78       0.01550  -0.00026 0.00000 0.00000 0.00000 0.01550
D79    -3.11809  -0.00085 0.00000 0.00000 0.00000  -3.11809
D80       0.02624  -0.00085 0.00000 0.00000 0.00000 0.02624
D81    -0.01122 0.00051 0.00000 0.00000 0.00000  -0.01122
D82       3.13311 0.00051 0.00000 0.00000 0.00000 3.13311
D83    -0.05051 0.00174 0.00000 0.00000 0.00000  -0.05051
D84       3.08159 0.00127 0.00000 0.00000 0.00000 3.08159
D85       3.09365 0.00136 0.00000 0.00000 0.00000 3.09365
D86    -0.05743 0.00089 0.00000 0.00000 0.00000  -0.05743
D87    -0.01082  -0.00013 0.00000 0.00000 0.00000  -0.01082
D88       3.13576  -0.00015 0.00000 0.00000 0.00000 3.13576
D89       3.12789 0.00030 0.00000 0.00000 0.00000 3.12789
D90    -0.00871 0.00028 0.00000 0.00000 0.00000  -0.00871
D91       0.00308  -0.00037 0.00000 0.00000 0.00000 0.00308
D92    -3.13741  -0.00014 0.00000 0.00000 0.00000  -3.13741
D93    -3.14125  -0.00037 0.00000 0.00000 0.00000  -3.14125
D94       0.00144  -0.00014 0.00000 0.00000 0.00000 0.00144
D95       0.00154  -0.00003 0.00000 0.00000 0.00000 0.00154
D96    -3.13936 0.00020 0.00000 0.00000 0.00000  -3.13936
D97    -3.14116  -0.00026 0.00000 0.00000 0.00000  -3.14116
D98       0.00113  -0.00003 0.00000 0.00000 0.00000 0.00113
D99       0.00227 0.00028 0.00000 0.00000 0.00000 0.00227
D100    3.13886 0.00030 0.00000 0.00000 0.00000 3.13886
D101    -3.14002 0.00005 0.00000 0.00000 0.00000  -3.14002
D102    -0.00343 0.00007 0.00000 0.00000 0.00000  -0.00343
      Item             Value    Threshold  Converged?
Maximum Force          0.036416    0.000002    NO
RMS    Force          0.003107    0.000001    NO
Maximum Displacement    0.000000    0.000006    YES
RMS    Displacement    0.000000    0.000004    YES
Predicted change in Energy=-1.224573D-16
Optimization aborted.
-- No acceptable step.
                        ----------------------------
                        ! Non-Optimized Parameters !
                        ! (Angstroms and Degrees)  !

老师您看看，谢谢！

ct_2015 · 发表于 Post on 2015-8-4 14:18:23

liyuanhe211 发表于 2015-8-3 21:50
以及考虑你的反应是不是基元反应，反应物和产物是否已经优化，结构是否合理

反应如图，反应过程不难，应该是基元反应吧，反应物和产物已优化了的。

liyuanhe211 · 发表于 Post on 2015-8-4 20:05:46

ct_2015 发表于 2015-8-4 14:05
表Non-Optimized Param之上找到如下：
199 D99 0.00318 -0.01522 0.000001000.00000 ...

“结构是否合理”指的是看一下你的过渡态已经变成什么样子了，不是让你把反应式贴出来。我没有遇到过Maximum Displacement是0，然后Gaussian不知道该往哪儿走的情况，不过我猜可能是现在得到的结构已经不合理了，拎出来看一下；如果发现过渡态结构合理，可以直接扔掉反应物和产物的构象，用当前进行不下去的过渡态构象做opt=(ts,calcfc)。（其实我个人一直感觉QST2成功率不高）

ct_2015 · 发表于 Post on 2015-8-5 14:39:41

liyuanhe211 发表于 2015-8-4 20:05
“结构是否合理”指的是看一下你的过渡态已经变成什么样子了，不是让你把反应式贴出来。我没有遇到过Maxi ...

，呵呵，菜鸟菜鸟。按照你的方法我试试，谢谢！

小书童 · 发表于 Post on 2015-8-6 10:51:24

学习了

happyrabbit · 发表于 Post on 2015-8-13 17:04:29

我认为是你的过渡态结构有点问题，你可以把最后这个结构提取出来，看看是否是你预想的接近，然后稍微调整一下结构参数，用opt=(ts,calcfc）重新优化，或者用qst3方法，并且减小优化步长（用 iop(1/8）来控制，如iop(1/8=5)或者更小iop(1/8=3））。

happyrabbit · 发表于 Post on 2015-8-13 17:09:17

需要注意的是，用qst3方法，过渡态和反应物、产物的原子顺序必须一致，有些麻烦。所以用ts方法比较省事，但做完以后必须要用来IRC验证过渡态的准确性。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] QST2方法寻找过渡态出现l9999.exe

浏览过的版块