应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 高斯计算中金属离子与氯苯的相互作用结构优化错误
52/1返回列表
查看: 3111  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

alexwpch

木虫 (正式写手)

[求助] 高斯计算中金属离子与氯苯的相互作用结构优化错误

尝试使用高斯研究金属离子与氯苯的相互作用，在结构优化的时候始终没有结果，出现各种错误，主要是收敛问题。而且经过优化后有的金属离子和氯苯会离得很远，这是初始构型不正确吗，需要如何改动。
将各种错误的输入和输出文件上传，请各位专业人士帮忙分析下，重新提供些输入方法。本人不是理论计算出身，所以可能会有些低级错误，大家轻喷。其中有个唯一运行正确的文件。[ Last edited by alexwpch on 2012-6-11 at 12:15 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : Fe-right.out
    • 2012-06-11 10:41:07, 764.22 K
  • 附件 2 : Nb-L502.out
    • 2012-06-11 10:41:16, 579.93 K
  • 附件 3 : Nb-L9999.out
    • 2012-06-11 10:41:19, 609.11 K
  • 附件 4 : Ti-L9999.out
    • 2012-06-11 10:41:25, 754.43 K
  • 附件 5 : Nb.gjf
    • 2012-06-11 10:41:32, 1.06 K
  • 附件 6 : Ti.gjf
    • 2012-06-11 10:41:46, 1.02 K
  • 附件 7 : Zr-L103.out
    • 2012-06-11 12:15:21, 755.39 K

» 猜你喜欢
1楼 2012-06-11 10:45:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alexwpch

木虫 (正式写手)
引用回帖:
12楼: Originally posted by meteoric30 at 2012-06-13 09:05:33
帖输入文件 我看下~~报错的时候几个收敛了？...

Item               Value     Threshold  Converged?
Maximum Force            0.040574     0.000450     NO
RMS     Force            0.013112     0.000300     NO
Maximum Displacement     0.064171     0.001800     NO
RMS     Displacement     0.014782     0.001200     NO
Predicted change in Energy=-1.206551D-03
四个都不收敛。没法上传附件。
最后一段是
Isotropic Fermi Contact Couplings
        Atom                 a.u.       MegaHertz       Gauss      10(-4) cm-1
     1  C(13)              0.00000       0.00000       0.00000       0.00000
     2  C(13)              0.00000       0.00000       0.00000       0.00000
     3  C(13)              0.00000       0.00000       0.00000       0.00000
     4  C(13)              0.00000       0.00000       0.00000       0.00000
     5  C(13)              0.00000       0.00000       0.00000       0.00000
     6  C(13)              0.00000       0.00000       0.00000       0.00000
     7  H(1)               0.00000       0.00000       0.00000       0.00000
     8  H(1)               0.00000       0.00000       0.00000       0.00000
     9  H(1)               0.00000       0.00000       0.00000       0.00000
    10  H(1)               0.00000       0.00000       0.00000       0.00000
    11  H(1)               0.00000       0.00000       0.00000       0.00000
    12  Cl(35)             0.00000       0.00000       0.00000       0.00000
    13  Zr(91)             0.00000       0.00000       0.00000       0.00000
--------------------------------------------------------
       Center         ----  Spin Dipole Couplings  ----
                      3XX-RR        3YY-RR        3ZZ-RR
--------------------------------------------------------
     1   Atom        0.000000      0.000000      0.000000
     2   Atom        0.000000      0.000000      0.000000
     3   Atom        0.000000      0.000000      0.000000
     4   Atom        0.000000      0.000000      0.000000
     5   Atom        0.000000      0.000000      0.000000
     6   Atom        0.000000      0.000000      0.000000
     7   Atom        0.000000      0.000000      0.000000
     8   Atom        0.000000      0.000000      0.000000
     9   Atom        0.000000      0.000000      0.000000
    10   Atom        0.000000      0.000000      0.000000
    11   Atom        0.000000      0.000000      0.000000
    12   Atom        0.000000      0.000000      0.000000
    13   Atom        0.000000      0.000000      0.000000
--------------------------------------------------------
                        XY            XZ            YZ
--------------------------------------------------------
     1   Atom        0.000000      0.000000      0.000000
     2   Atom        0.000000      0.000000      0.000000
     3   Atom        0.000000      0.000000      0.000000
     4   Atom        0.000000      0.000000      0.000000
     5   Atom        0.000000      0.000000      0.000000
     6   Atom        0.000000      0.000000      0.000000
     7   Atom        0.000000      0.000000      0.000000
     8   Atom        0.000000      0.000000      0.000000
     9   Atom        0.000000      0.000000      0.000000
    10   Atom        0.000000      0.000000      0.000000
    11   Atom        0.000000      0.000000      0.000000
    12   Atom        0.000000      0.000000      0.000000
    13   Atom        0.000000      0.000000      0.000000
--------------------------------------------------------


---------------------------------------------------------------------------------
              Anisotropic Spin Dipole Couplings in Principal Axis System
---------------------------------------------------------------------------------

       Atom             a.u.   MegaHertz   Gauss  10(-4) cm-1        Axes

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     1 C(13)  Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     2 C(13)  Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     3 C(13)  Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     4 C(13)  Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     5 C(13)  Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     6 C(13)  Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     7 H(1)   Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     8 H(1)   Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
     9 H(1)   Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
    10 H(1)   Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
    11 H(1)   Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
    12 Cl(35) Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000

              Baa     0.0000     0.000     0.000     0.000  1.0000  0.0000  0.0000
    13 Zr(91) Bbb     0.0000     0.000     0.000     0.000  0.0000  1.0000  0.0000
              Bcc     0.0000     0.000     0.000     0.000  0.0000  0.0000  1.0000


---------------------------------------------------------------------------------



YOU ARE IN A MAZE OF TWISTY LITTLE PASSAGES.
Error termination request processed by link 9999.
Error termination via Lnk1e in d:\gaussian09\l9999.exe at Tue Jun 12 23:56:21 2012.
Job cpu time:  0 days  1 hours 44 minutes 18.0 seconds.
File lengths (MBytes):  RWF=     25 Int=      0 D2E=      0 Chk=      2 Scr=      1
赞一下 回复此楼
14楼2012-06-13 09:56:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

Mr_WuSH

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
alexwpch: 金币+20, ★★★很有帮助, 谢谢哈，我试试这个方法看，先送你20金币，再没人回复的话剩下的都给你。 2012-06-12 09:09:13
gmy1990: 金币+2 2012-06-12 09:22:02
我觉得你可以用Fragment来优化试一试。不熟悉先不要用link，一步一步来·····
（下面的坐标是乱写的··可能跟你的不一样······，只是说一下Fragment的格式··）
# ···············  Counterpoise=2

   X   

0 1
C(Fragment=1) +0.00000000 +0.00000000 +0.00000000
C(Fragment=1) +0.00000000 +0.00000000 +1.39516000
C(Fragment=1) +1.20775100 +0.00000000 +2.09269800
C(Fragment=1) +2.41626000 -0.00119900 +1.39504400
C(Fragment=1) +2.41618200 -0.00167800 +0.00021900
C(Fragment=1) +1.20797600 -0.00068200 -0.69738200
H(Fragment=1) -0.95231700 +0.00045000 -0.54975900
H(Fragment=1) -0.95251300 +0.00131500 +1.94466800
H(Fragment=1) +1.20783100 +0.00063400 +3.19237800
H(Fragment=1) +3.36840300 -0.00125800 +1.94524400
H(Fragment=1) +3.36846301 -0.00262292 -0.54990300
H(Fragment=1) +1.20815900 -0.00086200 -1.79698600
H(Fragment=2) +1.28409831 +1.81161709 +0.70347558
Ti(Fragment=2) +1.33116591 +2.91011442 +0.68530890
赞一下 回复此楼
2楼2012-06-12 08:11:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
alexwpch: 金币+5, ★有帮助, 不只是L502，还有L9999等，还有其他金属离子也是一样的错误。我试着加maxcyc=200，这种情况容易出现L9999错误，qc没加过。 2012-06-12 10:32:47
2楼说的是bsse吧~~你如果主要是~收敛（l502）的问题的话，主要还是自洽场不收敛，在含fe的物质中计算 常见错误，加上（maxcycle=300，xqc）或者qc掉 比较合适~~
赞一下 回复此楼
3楼2012-06-12 09:25:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)
★
gmy1990: 金币+1, 3X 2012-06-12 20:14:35
引用回帖:
3楼: Originally posted by meteoric30 at 2012-06-12 09:25:05
2楼说的是bsse吧~~你如果主要是~收敛（l502）的问题的话，主要还是自洽场不收敛，在含fe的物质中计算 常见错误，加上（maxcycle=300，xqc）或者qc掉 比较合适~~

999不是错 最后一个拉出去 继续算 就行了 maxcycle 最好300比较合适
赞一下 回复此楼
4楼2012-06-12 10:35:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +37 医学老男孩 2026-08-13 87/4350 2026-08-18 06:35 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
[基金申请] 时间戳变了，能看出什么问题？ +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
关闭
请填处理意见
关闭