24小时热门版块排行榜    

查看: 2141  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 只有每年这种时候来逛逛小木虫 +23 yaoyewhu2008 2026-08-20 25/1250 2026-08-21 12:19 by Poppy1104
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 时间戳又变了 +12 wuchongjun 2026-08-20 18/900 2026-08-21 11:53 by wuchongjun
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 7/350 2026-08-21 11:45 by 6543yes
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 科研孤儿太难了 +15 我4大白菜 2026-08-20 16/800 2026-08-21 11:07 by yming1319
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 让我中一个面上吧! +10 大萍1987 2026-08-20 12/600 2026-08-21 09:31 by 020100099
[基金申请] 时间戳今天,20号变了 +4 archvillain 2026-08-20 4/200 2026-08-21 08:10 by gatelove
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +15 医学老男孩 2026-08-20 17/850 2026-08-20 20:35 by CaoLG
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
信息提示
请填处理意见