24小时热门版块排行榜    

查看: 3108  |  回复: 17

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

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的回帖

alexwpch

木虫 (正式写手)

引用回帖:
2楼: Originally posted by Mr_WuSH at 2012-06-12 08:11:05
我觉得你可以用Fragment来优化试一试。不熟悉先不要用link,一步一步来·····
(下面的坐标是乱写的··可能跟你的不一样······,只是说一下Fragment的格式··)
# ···············  C ...

试了你的方法,还是不收敛,L502错误。
输入文件如下:
# opt=(cartesian) uhf/genecp Counterpoise=2

Zr

4,1 0,1 4,1
C(Fragment=1)                     1.13913  -0.12117  -0.3112
C(Fragment=1)                     2.45988  -0.4199    0.04978
C(Fragment=1)                     3.34451   0.613     0.38812
C(Fragment=1)                     2.90839   1.94462   0.36548
C(Fragment=1)                     1.58764   2.24334   0.0045
C(Fragment=1)                     0.70301   1.21044  -0.33384
H(Fragment=1)                     0.4637   -0.90981  -0.56953
H(Fragment=1)                     2.79287  -1.43662   0.06706
H(Fragment=1)                     4.35294   0.38492   0.66373
H(Fragment=1)                     3.58383   2.73326   0.6238
H(Fragment=1)                     1.25466   3.26006  -0.01279
Cl(Fragment=1)                   -0.9557    1.58561  -0.7872
Zr(Fragment=2)                   -2.41086  -0.26799  -0.15452

Zr 0
Lanl2DZ
****
C H 0
3-21G
****
Cl 0
6-311G**(d,p)
****

Zr 0
Lanl2DZ
5楼2012-06-12 18:43:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alexwpch

木虫 (正式写手)

引用回帖:
4楼: Originally posted by meteoric30 at 2012-06-12 10:35:30
999不是错 最后一个拉出去 继续算 就行了 maxcycle 最好300比较合适...

嗯,我试下你的方法。
6楼2012-06-12 18:43:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
alexwpch: 金币+5, 有帮助 2012-06-12 20:38:11
引用回帖:
6楼: Originally posted by alexwpch at 2012-06-12 18:43:52
嗯,我试下你的方法。...

2楼给出其实是 bsse...我不懂 你试这个干什么。。。999不是错  你有空多看看 小卒版主的faq吧。。。
7楼2012-06-12 18:50:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alexwpch

木虫 (正式写手)

引用回帖:
7楼: Originally posted by meteoric30 at 2012-06-12 18:50:19
2楼给出其实是 bsse...我不懂 你试这个干什么。。。999不是错  你有空多看看 小卒版主的faq吧。。。...

不懂BSSE,我基本是自学的,从最简单的烷基苯的电子云密度开始学起,现在想进一步。小卒的faq看过,但是没全搞懂。不懂得只能先放一边,什么时候用了,什么时候再看。而且,我是有机合成的,实验为主,大部分时间需要做实验和看相关文献。
PS:终于知道为什么一些较有名的有机课题组都需要专门的做量化计算的人了。一个人同时学,做,看两个方向很是痛苦啊,而且量化还是没人指导,只能自己摸索。
再PS:我前面发过一个求助贴,请人帮我用G2或G3算一个小分子的生成热,不知道你能不能帮下忙。
http://muchong.com/bbs/viewthread.php?tid=4576382
8楼2012-06-12 20:30:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

多半是自旋多重度不合适,计算了一些的高激发态
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
9楼2012-06-13 00:34:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alexwpch

木虫 (正式写手)

引用回帖:
7楼: Originally posted by meteoric30 at 2012-06-12 18:50:19
2楼给出其实是 bsse...我不懂 你试这个干什么。。。999不是错  你有空多看看 小卒版主的faq吧。。。...

加了maxcyc=300,xqc,还是L9999错误,并且金属离子离得很远很远到了18A。
10楼2012-06-13 08:43:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 alexwpch 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳变了,能看出什么问题? +3 基诺咪客 2026-08-17 3/150 2026-08-17 15:25 by xinshengwu
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] filecode=后面第一个是大写字母 +7 wangze12014 2026-08-14 9/450 2026-08-17 14:59 by wangze12014
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 13:49 by EYK67XJMj64E
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 11:46 by 4GBAYCdVQoK3
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6GojgJvkDudM 2026-08-16 3/150 2026-08-17 09:49 by xLVPIRuSvCUe
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了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 启萌科技
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +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
[基金申请] 重要来源:本周末出结果 +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
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
信息提示
请填处理意见