需要注意的是,用qst3方法,过渡态和反应物、产物的原子顺序必须一致,有些麻烦。所以用ts方法比较省事,但做完以后必须要用来IRC验证过渡态的准确性。 |
我认为是你的过渡态结构有点问题,你可以把最后这个结构提取出来,看看是否是你预想的接近,然后稍微调整一下结构参数,用opt=(ts,calcfc)重新优化,或者用qst3方法,并且减小优化步长(用 iop(1/8)来控制,如iop(1/8=5)或者更小iop(1/8=3))。 |
学习了 |
liyuanhe211 发表于 2015-8-4 20:05 ![]() |
ct_2015 发表于 2015-8-4 14:05 “结构是否合理”指的是看一下你的过渡态已经变成什么样子了,不是让你把反应式贴出来。我没有遇到过Maximum Displacement是0,然后Gaussian不知道该往哪儿走的情况,不过我猜可能是现在得到的结构已经不合理了,拎出来看一下;如果发现过渡态结构合理,可以直接扔掉反应物和产物的构象,用当前进行不下去的过渡态构象做opt=(ts,calcfc)。(其实我个人一直感觉QST2成功率不高) |
liyuanhe211 发表于 2015-8-3 21:49 表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 发表于 2015-8-3 20:09 不要停留在:1 出错 2 没出错 这么低级的判断上。 要从原理出发,搞明白每种帮助收敛的做法的原理,而且必须要观看优化过程中结构变化的轨迹、能量/受力的变化曲线,根据实际情况判断怎么处理最合适。 没有万能的解决不收敛的方法! |
ct_2015 发表于 2015-8-3 20:09 以及考虑你的反应是不是基元反应,反应物和产物是否已经优化,结构是否合理 |
ct_2015 发表于 2015-8-3 20:09 输出文件往上找,在表Non-Optimized Param之上找到形如: Optimization stopped. -- Wrong number of Negative eigenvalues: Desired= 1 Actual= 0 # Or something else -- Flag reset to prevent archiving. 的行s,贴出来。 |
sobereva 发表于 2015-8-2 01:21 谢谢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老师看看出现了什么问题,谢谢! |
几乎一定是没收敛。增大maxcyc是白费功夫,只会浪费计算量。 解决方法此文有讨论,虽然讨论的是几何优化的,但大多数对于过渡态搜索同样管用 量子化学计算中帮助几何优化收敛的常用方法 http://sobereva.com/164 这些关键词别写:iop(5/13=1) iop(2/16=1) iop(1/11=1) maxcyc=800 1/11=1等价于noeigen,凡是有现成关键词的都不要以iop来写 |
試試 opt=noeigentest |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-14 13:13 , Processed in 0.278656 second(s), 25 queries , Gzip On.