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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 高斯计算中金属离子与氯苯的相互作用结构优化错误
北京石油化工学院2026年研究生招生接收调剂公告
52/1返回列表
查看: 2670  |  回复: 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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +4 Y-cap 2026-03-29 4/200 2026-03-29 10:16 by lbsjt
[考研] 298求调剂 +4 种圣赐 2026-03-28 4/200 2026-03-29 08:42 by q1092522407
[考研] 289求调剂 +13 新时代材料 2026-03-27 13/650 2026-03-29 01:16 by 544594351
[考研] 生物学学硕，一志愿湖南大学，初试成绩338 +6 YYYYYNNNNN 2026-03-26 7/350 2026-03-28 20:52 by 唐沐儿
[考研] 083000学硕274求调剂 +8 Li李鱼 2026-03-26 8/400 2026-03-28 20:33 by 加油向未来啊
[考研] 286求调剂 +12 PolarBear11 2026-03-26 12/600 2026-03-28 12:14 by zllcz
[考研] 材料277求调剂 +7 min3 2026-03-24 7/350 2026-03-28 11:39 by xuxiang
[考研] 291求调剂 +15 hhhhxn.. 2026-03-23 21/1050 2026-03-28 11:26 by self2008
[考研] 张芳铭-中国农业大学-环境工程专硕-298 +4 手机用户 2026-03-26 4/200 2026-03-28 07:17 by mmm just
[考研] 求调剂 +8 张zz111 2026-03-27 9/450 2026-03-28 03:41 by fmesaito
[考研] 352分 化工与材料 +5 海纳百川Ly 2026-03-27 5/250 2026-03-28 03:39 by fmesaito
[考研] 322求调剂 +4 我真的很想学习 2026-03-23 4/200 2026-03-27 13:51 by 杨杨杨紫
[考研] 材料求调剂 +5 .m.. 2026-03-25 5/250 2026-03-27 11:08 by 不吃魚的貓
[考研] 284求调剂 +11 junqihahaha 2026-03-26 12/600 2026-03-27 04:37 by wxiongid
[考研] 352求调剂 +4 大米饭！ 2026-03-22 4/200 2026-03-26 16:40 by 不吃魚的貓
[考研] 一志愿哈工大，085400，320，求调剂 +4 gdlf9999 2026-03-24 4/200 2026-03-25 23:01 by boxking200
[考研] 求调剂 +3 李李不服输 2026-03-25 3/150 2026-03-25 13:03 by cmz0325
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 求老师收我 +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
信息提示
关闭
请填处理意见
关闭