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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 高斯计算中金属离子与氯苯的相互作用结构优化错误
52/1返回列表
查看: 2456  |  回复: 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的回帖

zdgclf

银虫 (初入文坛)

【答案】应助回帖

四楼的说法是比较准确的，不过你可以多金属离子采用lanl2dz 那样计算的结果会更好！
%chk=F:\fe.chk
# opt rhf/genecp geom=connectivity pop=mk

Title Card Required

0 1
C                  0.84805700    0.97511100    0.33143200
C                  2.14219300    1.49348500    0.18852800
C                  3.20653100    0.64133500   -0.13542700
C                  2.97673300   -0.72918900   -0.31647800
C                  1.68259700   -1.24756300   -0.17357200
C                  0.61825900   -0.39541300    0.15038300
H                  0.03541100    1.62574500    0.57878200
H                  2.31764900    2.53990900    0.32676700
H                  4.19463300    1.03712500   -0.24453500
H                  3.78937900   -1.37982400   -0.56382300
H                  1.50714200   -2.29398900   -0.31180500
Cl                -1.00703000   -1.04643300    0.32985600
Fe                -2.44503600    0.45514600   -0.21777400

1 2 1.5 6 1.5 7 1.0
2 3 1.5 8 1.0
3 4 1.5 9 1.0
4 5 1.5 10 1.0
5 6 1.5 11 1.0
6 12 1.0
7
8
9
10
11
12 13 1.0
13

C H Cl 0
6-31g
****
Fe 0
lanl2dz
****
这是输入文件
赞一下 回复此楼
17楼2012-06-14 15:16:20
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭