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

登录注册

应《网络安全法》要求，自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

» 猜你喜欢

散金求好运已经有315人回复
博士生招生, 材料模拟方向，University of South Carolina 已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有289人回复
Postdoc position, materials simulations/machine learning, U. of South Carolina 已经有0人回复
祈福青基已经有34人回复
牛来！米来！面来！已经有8人回复
2027申博咨询已经有11人回复
寻求CrC2Al的cif文件，感谢已经有0人回复
大湾区大学刘天辉课题组招收2027级博士生（与中山大学物理学院联合培养）已经有55人回复
ORCA 纯新手向入门指南0.7版已经有0人回复
哈尔滨工业大学（深圳）绿色合成化学团队招收27年春季或秋季入学的博士生已经有5人回复

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