24小时热门版块排行榜    

查看: 1541  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 310求调剂 +3 baibai1314 2026-03-16 3/150 2026-03-21 03:56 by JourneyLucky
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 278求调剂 +6 烟火先于春 2026-03-17 6/300 2026-03-21 01:57 by JourneyLucky
[考研] 280求调剂 +7 咕噜晓晓 2026-03-18 8/400 2026-03-21 01:27 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕,总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 14/700 2026-03-20 09:30 by kkcoco25
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 328求调剂,英语六级551,有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
信息提示
请填处理意见