24小时热门版块排行榜    

查看: 1683  |  回复: 11

alexwpch

木虫 (正式写手)

[求助] 求助ONIOM出错原因,并请教解决方法

尝试做ONIOM,最后优化出问题了,请教下原因,并寻求解决方法。

结尾错误如下:
Search for a local minimum.
Step number  19 out of a maximum of   67
All quantities printed in internal units (Hartrees-Bohrs-Radians)
                          R1        R2        R3        R4        R5
           -DE/DX    **************************************************
                          R6        R7        R8        R9        R10
           -DE/DX    **************************************************
                          R11       R12       R13       A1        A2
           -DE/DX    **************************************************
                          A3        A4        A5        A6        A7
           -DE/DX    **************************************************
                          A8        A9        A10       A11       A12
           -DE/DX    **************************************************
                          A13       A14       A15       A16       A17
           -DE/DX    **************************************************
                          A18       A19       A20       A21       D1
           -DE/DX    **************************************************
                          D2        D3        D4        D5        D6
           -DE/DX    **************************************************
                          D7        D8        D9        D10       D11
           -DE/DX    **************************************************
                          D12       D13       D14       D15       D16
           -DE/DX    **************************************************
                          D17       D18       D19       D20       D21
           -DE/DX    **************************************************
                          D22       D23
           -DE/DX    ********************
Optimization aborted.
    -- Gradient out of range.
    -- Maximum allowed force =    0.100E+21
Error termination via Lnk1e in D:\G09W\l103.exe at Sun Jul 07 16:03:45 2013.
Job cpu time:  0 days  0 hours  8 minutes 52.0 seconds.
File lengths (MBytes):  RWF=      9 Int=      0 D2E=      0 Chk=      2 Scr=      1


输入文件如下:
#p opt oniom(b3lyp/6-31g:uff) geom=connectivity

ZV-P

0 1 0 2 0 2
Zr-              0   -1.16809461    0.04912786   -0.24112541 L
O-O_3            0   -3.21978300   -0.00804037   -0.92735552 L
O-O_3            0   -2.38753622    1.40752047    0.91992916 L
O-O_3            0   -1.07563141    0.57732286   -2.33773079 L
O-O_3            0   -0.25630218    1.99392294   -0.50329773 L
V-V_3+5          0   -2.79265729    1.11929080   -2.26380189 L
V-V_3+5          0   -1.84461172    2.76732379   -0.13759530 L H-H_ 10
O-O_3            0   -2.87992423    2.86007451   -1.66384313 L
O-O_2            0   -3.66704094    1.00326272   -3.64260751 L
O-O_2            0   -1.80506682    4.25728015    0.53885166 H
O-O_3            0   -0.05205552    0.29123985    1.53311767 L
H-H_             0    0.36061688   -0.54311932    1.76795011 L
O-O_3            0   -1.83478572   -1.35785045    1.18296553 L
H-H_             0   -1.21784436   -2.09242472    1.22016767 L
O-O_3            0   -1.33894851   -1.84135852   -1.16252120 L
H-H_             0   -0.77306675   -2.47057813   -0.70924368 L
O-O_3            0    0.78660995   -0.64592204   -0.62600945 L
H-H_             0    1.02575155   -1.30689306    0.02784727 L
C-C_3            0   -2.86426160    6.97609297   -1.64972010 H
C-C_3            0   -1.44660043    7.38164337   -1.96159637 H
H-H_             0   -3.45118732    7.85543154   -1.29040782 H
H-H_             0   -3.36527509    6.56532784   -2.55950900 H
C-C_3            0   -0.07450889    6.97953477   -2.43815594 H
H-H_             0   -1.45139284    8.18892386   -2.74131635 H
H-H_             0   -0.97070968    7.81543624   -1.04234363 H
H-H_             0    0.61171042    7.86037149   -2.43867988 H
H-H_             0    0.35905664    6.19473045   -1.77212001 H
H-H_             0   -0.12293355    6.56923678   -3.47585272 H
O-O_3            0   -2.56690762    5.74527051   -0.72861262 H
H-H_             0   -2.90524174    5.29547201    0.04908354 H
回复此楼

» 猜你喜欢

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

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

alexwpch

木虫 (正式写手)

没人帮忙解决下问题吗?请各位高手帮下哇。
2楼2013-07-31 10:32:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

568286346

木虫 (小有名气)

【答案】应助回帖

引用回帖:
2楼: Originally posted by alexwpch at 2013-07-31 10:32:38
没人帮忙解决下问题吗?请各位高手帮下哇。

能上传看下输入文件吗?
坚持才能实现自己的理想!
3楼2013-07-31 16:20:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alexwpch

木虫 (正式写手)

引用回帖:
3楼: Originally posted by 568286346 at 2013-07-31 16:20:31
能上传看下输入文件吗?...

原帖有输入文件。。。。。
4楼2013-07-31 19:00:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

568286346

木虫 (小有名气)

【答案】应助回帖

我把你输入文件里面的geom=connectivity删除后可以进行算了,你自己尝试看看。
坚持才能实现自己的理想!
5楼2013-07-31 19:57:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alexwpch

木虫 (正式写手)

引用回帖:
5楼: Originally posted by 568286346 at 2013-07-31 19:57:41
我把你输入文件里面的geom=connectivity删除后可以进行算了,你自己尝试看看。

好的,我试了给你答复。
6楼2013-07-31 21:52:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alexwpch

木虫 (正式写手)

引用回帖:
5楼: Originally posted by 568286346 at 2013-07-31 19:57:41
我把你输入文件里面的geom=connectivity删除后可以进行算了,你自己尝试看看。

你好,删除了geom=connectivity以后,还是出现同样的错误,可能不是这个的原因吧。
7楼2013-07-31 22:30:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

568286346

木虫 (小有名气)

【答案】应助回帖

引用回帖:
7楼: Originally posted by alexwpch at 2013-07-31 22:30:05
你好,删除了geom=connectivity以后,还是出现同样的错误,可能不是这个的原因吧。...

那应该就是你选的基组的问题了,因为你计算的体系中含有过渡金属元素,对于过渡金属元素6-31g是不适用的,因此最好是采用赝势基组计算,具体的基组可以根据自己的要求来选择,一般计算的话lanl2dz应该就可以。
坚持才能实现自己的理想!
8楼2013-08-02 08:21:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alexwpch

木虫 (正式写手)

引用回帖:
8楼: Originally posted by 568286346 at 2013-08-02 08:21:18
那应该就是你选的基组的问题了,因为你计算的体系中含有过渡金属元素,对于过渡金属元素6-31g是不适用的,因此最好是采用赝势基组计算,具体的基组可以根据自己的要求来选择,一般计算的话lanl2dz应该就可以。...

如果是基组问题的话,是不是应该在一开始就会提示错误。仔细看输出文件的
“Search for a local minimum.
Step number  19 out of a maximum of   67”这句话,是在计算到一半的时候出现的错误,错误原因是“Gradient out of range”。
是否可以这么理解?
9楼2013-08-02 08:42:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

568286346

木虫 (小有名气)

【答案】应助回帖

我自己也算了下,后面的确会出现错误,但始终没有出现你说的错误啊。你看下附件里我的输出文件。

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 1.gjf
  • 2013-08-02 10:06:31, 1.97 K
  • 附件 2 : 1.LOG
  • 2013-08-02 10:06:32, 1.04 M
坚持才能实现自己的理想!
10楼2013-08-02 10:06:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 alexwpch 的主题更新
信息提示
请填处理意见