24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1548  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料0856 英一数二 323 求调剂 +14 袁sy 2026-04-01 14/700 2026-04-05 18:18 by cql1109
[考研] 358求调剂 +7 秋gk 2026-04-04 7/350 2026-04-05 13:29 by huangmoli
[考研] +5 化工专硕323分 2026-04-04 5/250 2026-04-05 08:02 by 544594351
[考研] 324求调剂 +9 想上学求调 2026-04-03 9/450 2026-04-04 23:57 by 果冻大王
[考研] 材料调剂 +10 懒羊羊轻置玉臀 2026-04-02 11/550 2026-04-04 21:56 by laoshidan
[考研] 085601,一志愿厦大334复试被刷求调剂 +13 曾仰之 2026-04-03 15/750 2026-04-04 20:13 by dongzh2009
[考研] 0703求调剂383分 +8 W55j 2026-04-03 8/400 2026-04-04 20:09 by xhai2011
[考研] +5 雾与海 2026-04-02 6/300 2026-04-04 19:53 by 蓝云思雨
[考研] 一志愿085404,总分291,四级已过,求调剂 +5 阿俊阿俊阿俊 2026-04-04 7/350 2026-04-04 13:23 by 莲菜就是藕吧
[考研] 334求调剂 +8 曾仰之 2026-04-03 8/400 2026-04-04 11:16 by w_xuqing
[考研] 292求调剂 +21 是妍子也是研子 2026-03-30 22/1100 2026-04-03 21:44 by qlm5820
[考研] 316求调剂 +14 舟自梗 2026-04-01 18/900 2026-04-03 10:28 by linyelide
[考研] 326求调剂 +10 崽崽仔 2026-04-02 10/500 2026-04-03 09:08 by 帕尔马拉特
[考研] 362求调剂 +14 西南交材料专硕3 2026-03-31 14/700 2026-04-02 17:50 by yunlongyang
[考研] 321求调剂 一志愿 浙江工业大学生物医药 +5 嘿嘿HC 2026-04-01 6/300 2026-04-02 15:23 by sophie2180
[考研] 0710生物学,325求调剂 +3 mkkkkkl 2026-04-01 3/150 2026-04-02 09:48 by Jaylen.
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
[硕博家园] 博一被送出联培感觉不适应怎么办 +3 全村的狗 2026-03-31 3/150 2026-04-01 10:44 by 328838485
[考研] 生物考研337分求调剂 +4 cgxin 2026-03-30 6/300 2026-03-31 14:18 by 记事本2026
[考研] 调剂 +4 GK72 2026-03-30 4/200 2026-03-30 20:32 by dick_runner
信息提示
请填处理意见