24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2743  |  回复: 18

juty136

金虫 (小有名气)


[交流] 出现四个yes 却出现不收敛 9999错误

各位大侠好,
  最近我进行几何优化的时候出现了一个怪现象,就是出现了四个yes 但是还是没有收敛, 而且出现了好多的四个YES ,最后就9999错误结速, 不知道大侠们遇见了这种情况了没? 不知怎么解决!希望大家交流一下!  谢谢 !

GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Internal  Forces:  Max     0.000309573 RMS     0.000061536
Search for a local minimum.
Step number  61 out of a maximum of  100
All quantities printed in internal units (Hartrees-Bohrs-Radians)
Swaping is turned off.
Update second derivatives using D2CorX and points    6    7    8    9   10
                                                     11   12   13   14   15
                                                     16   17   18   19   20
                                                     21   22   23   24   25
                                                     26   27   28   29   30
                                                     31   32   33   34   35
                                                     36   37   38   39   40
                                                     41   42   43   44   45
                                                     46   47   48   49   50
                                                     51   52   53   54   55
                                                     61
Linear search not attempted -- option 19 set.
Steepest descent instead of Quadratic search.
Iteration  1 RMS(Cart)=  0.00002249 RMS(Int)=  0.00000000
Iteration  2 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        4.46682   0.00000   0.00000   0.00000   0.00000   4.46682
    R2        7.45408   0.00020   0.00000   0.00010   0.00010   7.45418
    R3        5.79140   0.00031   0.00000   0.00015   0.00015   5.79155
    R4        5.91873   0.00012   0.00000   0.00006   0.00006   5.91879
    R5        7.43495   0.00018   0.00000   0.00009   0.00009   7.43504
    R6        2.89394   0.00000   0.00000   0.00000   0.00000   2.89394
    R7        2.08241   0.00000   0.00000   0.00000   0.00000   2.08241
    R8        2.07941   0.00000   0.00000   0.00000   0.00000   2.07941
    R9        2.07862   0.00000   0.00000   0.00000   0.00000   2.07862
   R10        2.07854   0.00000   0.00000   0.00000   0.00000   2.07854
   R11        2.08293   0.00000   0.00000   0.00000   0.00000   2.08293
   R12        2.07934   0.00000   0.00000   0.00000   0.00000   2.07934
    A1        0.68820  -0.00001   0.00000  -0.00001  -0.00001   0.68819
    A2        0.68972  -0.00001   0.00000   0.00000   0.00000   0.68971
    A3        1.93114   0.00002   0.00000   0.00001   0.00001   1.93115
    A4        1.95114  -0.00001   0.00000  -0.00001  -0.00001   1.95113
    A5        1.94485   0.00000   0.00000   0.00000   0.00000   1.94485
    A6        1.88275  -0.00001   0.00000   0.00000   0.00000   1.88275
    A7        1.87102   0.00000   0.00000   0.00000   0.00000   1.87102
    A8        1.87967  -0.00001   0.00000   0.00000   0.00000   1.87967
    A9        1.94534   0.00000   0.00000   0.00000   0.00000   1.94534
   A10        1.92919   0.00002   0.00000   0.00001   0.00001   1.92921
   A11        1.95185  -0.00001   0.00000   0.00000   0.00000   1.95185
   A12        1.87141   0.00000   0.00000   0.00000   0.00000   1.87141
   A13        1.88027  -0.00001   0.00000   0.00000   0.00000   1.88027
   A14        1.88252  -0.00001   0.00000   0.00000   0.00000   1.88252
   A15        1.90568   0.00002   0.00000   0.00001   0.00001   1.90569
   A16        2.50085   0.00000   0.00000   0.00000   0.00000   2.50085
   A17        2.52992  -0.00003   0.00000  -0.00001  -0.00001   2.52991
   A18        1.93574  -0.00001   0.00000  -0.00001  -0.00001   1.93573
    D1       -0.28375   0.00001   0.00000   0.00001   0.00001  -0.28374
    D2       -0.50549   0.00003   0.00000   0.00002   0.00002  -0.50548
    D3       -0.48473   0.00000   0.00000   0.00000   0.00000  -0.48473
    D4       -0.28641   0.00002   0.00000   0.00001   0.00001  -0.28640
    D5        1.09169   0.00003   0.00000   0.00001   0.00001   1.09171
    D6       -0.98775   0.00001   0.00000   0.00001   0.00001  -0.98775
    D7       -3.08549   0.00001   0.00000   0.00001   0.00001  -3.08548
    D8       -1.00718   0.00003   0.00000   0.00002   0.00002  -1.00717
    D9       -3.08663   0.00002   0.00000   0.00001   0.00001  -3.08662
   D10        1.09882   0.00002   0.00000   0.00001   0.00001   1.09883
   D11       -3.11158   0.00005   0.00000   0.00002   0.00002  -3.11156
   D12        1.09216   0.00003   0.00000   0.00002   0.00002   1.09217
   D13       -1.00557   0.00003   0.00000   0.00002   0.00002  -1.00556
   D14        0.66111  -0.00001   0.00000  -0.00001  -0.00001   0.66111
   D15        0.97427  -0.00001   0.00000   0.00000   0.00000   0.97426
   D16        2.80039  -0.00002   0.00000  -0.00001  -0.00001   2.80037
   D17        3.11354  -0.00002   0.00000  -0.00001  -0.00001   3.11353
   D18       -1.46286  -0.00003   0.00000  -0.00002  -0.00002  -1.46288
   D19       -1.14971  -0.00003   0.00000  -0.00001  -0.00001  -1.14972
   D20        0.99877  -0.00004   0.00000  -0.00002  -0.00002   0.99875
   D21        0.66886  -0.00003   0.00000  -0.00001  -0.00001   0.66885
   D22       -1.12490  -0.00006   0.00000  -0.00003  -0.00003  -1.12492
   D23       -1.45481  -0.00004   0.00000  -0.00002  -0.00002  -1.45483
   D24        3.13756  -0.00004   0.00000  -0.00002  -0.00002   3.13754
   D25        2.80765  -0.00003   0.00000  -0.00001  -0.00001   2.80764
         Item               Value     Threshold  Converged?
Maximum Force            0.000310     0.000450     YES
RMS     Force            0.000062     0.000300     YES
Maximum Displacement     0.000122     0.001800     YES
RMS     Displacement     0.000022     0.001200     YES
Predicted change in Energy=-1.040478D-07
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

No NMR shielding tensors so no spin-rotation constants.
Leave Link  601 at Mon Jul 30 17:06:53 2012, MaxMem=  117964800 cpu:       1.0
(Enter /home/juty/g09/l9999.exe)


Everything starts as somebody's daydream.
                           -- Larry Niven
Error termination request processed by link 9999.
Error termination via Lnk1e in /home/juty/g09/l9999.exe at Mon Jul 30 17:06:53 2012.
Job cpu time:  0 days  4 hours 38 minutes  1.4 seconds.
File lengths (MBytes):  RWF=     25 Int=      0 D2E=      0 Chk=      5 Scr=      1
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

meteoric30

木虫 (正式写手)


★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
juty136: 金币+1, 谢谢交流 2012-07-31 09:36:58
拉出来继续优化呗,应该很快了~~
2楼2012-07-30 18:39:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nyzhaoyin

新虫 (正式写手)


★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
juty136: 金币+1, 谢谢交流 2012-07-31 09:38:15
最后一步的坐标拿出来接着优化
6楼2012-07-30 20:26:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

t13340033021

金虫 (正式写手)


输入文件看下?
8楼2012-07-30 22:03:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juty136

金虫 (小有名气)


引用回帖:
2楼: Originally posted by meteoric30 at 2012-07-30 18:39:07
拉出来继续优化呗,应该很快了~~

我也做过 但是一直没有收敛,
9楼2012-07-31 09:37:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juty136

金虫 (小有名气)


引用回帖:
6楼: Originally posted by nyzhaoyin at 2012-07-30 20:26:31
最后一步的坐标拿出来接着优化

一样的 我试过好多次 还是9999
10楼2012-07-31 09:38:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juty136

金虫 (小有名气)


引用回帖:
8楼: Originally posted by t13340033021 at 2012-07-30 22:03:20
输入文件看下?

$rungauss
%chk=IM1
%nproc=4
%mem=2000mb
#P UBPW91/genecp SCF(Vshift=200,MaxCycle=2000,xqc)
optcyc=800 opt(noeigen,steep,maxcycle=2000) freq guess=mix
下面的就是坐标 和基组 因为这个中间体是非常弱的键  必须加STEEP才不跑开, 对于过渡态金属用的的 lanl2tz   其他主族元素是用的 6-311++G**
11楼2012-07-31 09:44:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cxyuan09

金虫 (正式写手)


★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
juty136: 金币+2, 谢谢 2012-07-31 10:52:54
这个可能和机器有关系,如果换一台机器就能算完
那么你的优化实际上已经完成了
好像如果使用opt+freq同时算,也只算opt部分。
如果是这样,你不妨宠重装一下Gaussian甚至连系统也重装一下,可能解决问题
12楼2012-07-31 10:41:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juty136

金虫 (小有名气)


引用回帖:
12楼: Originally posted by cxyuan09 at 2012-07-31 10:41:11
这个可能和机器有关系,如果换一台机器就能算完
那么你的优化实际上已经完成了
好像如果使用opt+freq同时算,也只算opt部分。
如果是这样,你不妨宠重装一下Gaussian甚至连系统也重装一下,可能解决问题

这个还可以与机器有关啊,系统重装就很麻烦了!
我先试一试先优化 然后在算频率,
还有一个问题 想请教一下,
我在做stable的时候,遇到不稳定,出现508错误, 不知道怎样调试才能stable  希望大侠多多指点!谢谢了!
13楼2012-07-31 10:52:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
juty136: 金币+2, 谢谢 2012-07-31 16:52:23
引用回帖:
13楼: Originally posted by juty136 at 2012-07-31 10:52:42
这个还可以与机器有关啊,系统重装就很麻烦了!
我先试一试先优化 然后在算频率,
还有一个问题 想请教一下,
我在做stable的时候,遇到不稳定,出现508错误, 不知道怎样调试才能stable  希望大侠多多指点!谢谢 ...

怕重装系统,可以先重装一下高斯
l508出错的信息需要贴出来,光说1508错误,不知道什么原因啊。

一般,我遇到波函数不稳定时,对应的结构我就不要了,但这样做不一定使用你的研究。
14楼2012-07-31 14:43:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tswangan

银虫 (职业作家)


★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
juty136: 金币+1, 谢谢交流 2012-07-31 16:53:44
实际上完成了,只是999程序调用错误,优化结果是最后一次测坐标。
15楼2012-07-31 14:50:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juty136

金虫 (小有名气)


引用回帖:
14楼: Originally posted by hairan at 2012-07-31 14:43:04
怕重装系统,可以先重装一下高斯
l508出错的信息需要贴出来,光说1508错误,不知道什么原因啊。

一般,我遇到波函数不稳定时,对应的结构我就不要了,但这样做不一定使用你的研究。...

谢谢提醒,我已经在另外个求助贴上贴出来了
希望大侠能指导一下! 谢谢 !
http://muchong.com/bbs/viewthread.php?tid=4783280
16楼2012-07-31 16:52:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juty136

金虫 (小有名气)


引用回帖:
15楼: Originally posted by tswangan at 2012-07-31 14:50:27
实际上完成了,只是999程序调用错误,优化结果是最后一次测坐标。

我要算频率,9999错误了就不能接着往下算了!
17楼2012-07-31 16:53:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sch

铜虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
No NMR shielding tensors so no spin-rotation constants.
18楼2012-09-28 13:39:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oyljw

至尊木虫 (职业作家)



小木虫: 金币+0.5, 给个红包,谢谢回帖
把结构存下来(gv存成gjf),重新改下输入 在计算,输入行改成      opt(maxcyc=1024),scf(maxcyc=1024)freq或其他
19楼2012-10-17 19:12:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
young5173楼
2012-07-30 19:10   回复  
jiji1884楼
2012-07-30 19:28   回复  
郑克宁5楼
2012-07-30 20:11   回复  
yxu19867楼
2012-07-30 20:56   回复  
相关版块跳转 我要订阅楼主 juty136 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见