24小时热门版块排行榜    

查看: 1430  |  回复: 8

huosuhong

铜虫 (初入文坛)

[求助] 用CCSD(T)/6-311++g(d,p)算能量出错

各位量化同仁请帮忙看看这个问题怎么解决
我的jgf文件
# CCSD(T)/6-311++g(d,p) scf=qc

Title Card Required

0 1
C                  0.77480300    1.79573300   -0.08277500
H                  1.26161700    2.43467400    0.66327900
H                  1.08416100    2.13297800   -1.07598800
C                 -0.77480200    1.79573400    0.08277800
H                 -1.26161700    2.43467600   -0.66327500
H                 -1.08416000    2.13297500    1.07599200
C                  2.79226600   -1.06523300   -0.30515400
H                  3.07668000   -0.74644200   -1.31131600
H                  3.59664100   -0.81032200    0.39050700
H                  2.64812500   -2.14833200   -0.29331600
C                 -2.79226400   -1.06523400    0.30515800
H                 -3.07667000   -0.74645100    1.31132500
H                 -3.59664400   -0.81031800   -0.39049500
H                 -2.64812300   -2.14833400    0.29331100
Ge                -1.12895000   -0.16392300   -0.25559700
Ge                 1.12894900   -0.16392400    0.25559500

out文件如下[ Last edited by huosuhong on 2012-4-2 at 18:25 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : Ge2a.out
  • 2012-04-02 18:18:04, 16.92 K

» 猜你喜欢

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

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

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★
gmy1990: 金币+2 2012-04-02 22:26:29
CCSD(T)方法不可能算这么大的分子
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2012-04-02 19:56:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sculhf

荣誉版主 (知名作家)

断翼天使

优秀版主优秀版主

【答案】应助回帖

感谢参与,应助指数 +1
为什么要scf=qc这个关键词?
本来体系已经够大的了
去掉试一下
然后再适当降低一下scf收敛标准
谋事在人成事在天
3楼2012-04-03 12:07:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+2 2012-04-03 22:16:09
huosuhong: 金币+5 2012-04-05 08:47:55
你的分子中Ge是3价?或要让二Ge成双键?这是否是不收敛的原因?

G03应不能算CCSD(T)。
G09W也不能算CCSD(T),只能算到CCSD。输入CCSD(T)也不会报错,只是输出一个对应基组的RHF计算给你作交待。
你现在RHF计算都不收敛。

[ Last edited by zhou2009 on 2012-4-3 at 20:20 ]
4楼2012-04-03 17:44:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huosuhong

铜虫 (初入文坛)

引用回帖:
3楼: Originally posted by sculhf at 2012-04-03 12:07:27:
为什么要scf=qc这个关键词?
本来体系已经够大的了
去掉试一下
然后再适当降低一下scf收敛标准

非常感谢,我尝试写了个scf=(qc,conver=6)然后就正常结束了。有一点我不明白,就是scf=qc这个关键词不是强制收敛的吗?这个对于体系大的不适用吗?楼主是名量子化学研究生现读研二,许多知识还不成熟,望不吝赐教,非常感谢!
加油
5楼2012-04-05 08:44:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huosuhong

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by beefly at 2012-04-02 19:56:53:
CCSD(T)方法不可能算这么大的分子

beefly您好,非常感谢您的帮助,请问CCSD(T)方法适用于多大的分子啊
加油
6楼2012-04-05 08:47:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huosuhong

铜虫 (初入文坛)

引用回帖:
4楼: Originally posted by zhou2009 at 2012-04-03 17:44:44:
你的分子中Ge是3价?或要让二Ge成双键?这是否是不收敛的原因?

G03应不能算CCSD(T)。
G09W也不能算CCSD(T),只能算到CCSD。输入CCSD(T)也不会报错,只是输出一个对应基组的RHF计算给你作交待。
你现在RHF计 ...

非常感谢zhou2009的帮助,我一直用g03算CCSD(T)能量。楼主的这个问题我通过写scf=(qc,conver=6)解决了,正常结束了,我想把这个的输出结果让您看看,以您的意思,是不是最后那个CCSD(T)的能量不对啊?
加油
7楼2012-04-05 08:55:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huosuhong

铜虫 (初入文坛)

这是楼主通过加上scf=(qc,conver=6)后正常结束了的out文件

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : Ge2a.out
  • 2012-04-05 09:01:29, 31.25 K
加油
8楼2012-04-05 09:01:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

看了你的out文件。感到整个out只是RHF的结果。
虽然在得到“SCF Done:  E(RHF) =”之后,进行了CCSD(T)计算,得到了最后的“E(CORR)=”。但此后的分析却完全是基于RHF的结果,即“Population analysis using the SCF density”。后面的MO也基于此。CCSD(T)计算并没有真正进行下去。

如果是CCSD之类,在最后的“E(CORR)=”收敛之后,它后面的分析是基于“Population analysis using the QCI/CC density”,已经完全不是RHF的MO了。
这时没有MO的能级值,因为在CCSD之类,MO已经只是分数电子的占有了。不是RHF的每个MO占有两个电子。

所以我认为你现在的结果是RHF的,收敛也只是RHF的,不是CCSD(T)的。
作CCSD之类计算,需要同时加关键词:density。

[ Last edited by zhou2009 on 2012-4-5 at 10:01 ]
9楼2012-04-05 09:57:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huosuhong 的主题更新
信息提示
请填处理意见