版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

alexwpch

木虫 (正式写手)

应助: 15 (小学生)
金币: 2566.2
散金: 134
红花: 3
帖子: 451
在线: 705.4小时
虫号: 672273
注册: 2008-12-11
性别: GG
专业: 应用有机化学

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

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

» 猜你喜欢

26年秋季博士申请已经有0人回复
推荐给英语教学者的一本单词书《金鱼单词讲义：从26个拉丁字母到106万个英语单词》已经有59人回复
物理化学论文润色/翻译怎么收费? 已经有117人回复
推荐给教师的一本单词书《金鱼单词讲义：从26个拉丁字母到106万个英语单词》已经有32人回复
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生已经有11人回复
核磁分析软件MestReNova打开文件时报错已经有0人回复
在职博后不能申请博后基金了，那么在职博后意义何在？已经有2人回复
青岛大学化学化工学院分子测量学研究院2026年招收博士研究生已经有0人回复
香港科技大学（广州）诚招电催化方向博士生（2026秋入学）已经有0人回复

1楼 2012-06-11 10:45:09

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

alexwpch

木虫 (正式写手)

应助: 15 (小学生)
金币: 2566.2
散金: 134
红花: 3
帖子: 451
在线: 705.4小时
虫号: 672273
注册: 2008-12-11
性别: GG
专业: 应用有机化学

引用回帖:

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

金虫 (小有名气)

应助: 9 (幼儿园)
金币: 96.2
散金: 290
帖子: 121
在线: 74.3小时
虫号: 1850241
注册: 2012-06-06
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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

木虫 (正式写手)

应助: 43 (小学生)
金币: 1719.5
散金: 12
红花: 2
帖子: 436
在线: 320.7小时
虫号: 1046846
注册: 2010-06-24
性别: GG
专业: 催化化学

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +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

木虫 (正式写手)

应助: 43 (小学生)
金币: 1719.5
散金: 12
红花: 2
帖子: 436
在线: 320.7小时
虫号: 1046846
注册: 2010-06-24
性别: GG
专业: 催化化学

★
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 个回答