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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » ccsd(t)输出
171/2返回列表12下一页
查看: 1235  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

kent1022

金虫 (著名写手)

[交流] ccsd(t)输出 已有4人参与

gaussian03中的ccsd(t)计算完成后，最后的能量输出汇总只显示到mp3能量，到09了还是这样，有什么方便的方法可以将每个单点能计算后的CCSD(T)= .......  提取出来？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-18 11:24:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
zzy870720z: 金币+2, 谢谢指教 2012-05-18 15:45:06
到目前的gaussian09还不能算CCSD(T)的电子密度。
但可以算到CCSD的水平。CCSD加上关键词density，才能进行真正的CCSD计算，输出单点能E(CORR)=.......  （可提取出来）。

[ Last edited by zhou2009 on 2012-5-21 at 07:56 ]
赞一下 回复此楼
2楼2012-05-18 14:19:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanrding

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
啥个意思啊？ CCSD(T)在gaussian中不是常规计算吗？
每个单点能计算后的CCSD(T)= 能量不是直接就给出来了吗？

估计是我理解有问题了。
赞一下 回复此楼
华北电力大学（北京）数理学院，计算化学，凝聚态物理。欢迎合作！招硕士/博士研究生，招老师。
3楼2012-05-18 21:24:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kent1022

金虫 (著名写手)
难道怎样算的不是CCSD(T)能量？
糊涂了
赞一下 回复此楼
4楼2012-05-18 22:27:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
gmy1990: 金币+3 2012-05-19 08:58:47
CCSD(T)在gaussian09-W-A.02中，虽然没有报错，但只给出RHF的同基组计算结果，即能量表示为:E(RHF)=...... 并没有真正有CCSD(T)计算的波函数、电子密度。

但gaussian09可以算到CCSD的水平。
CCSD必须同时加上关键词density，才能进行真正的CCSD计算，输出单点能E(CORR)=.......
这时MO的表示已经完全不是RHF的表示了。
在CCSD，已不再是每个MO占有两个电子，而是分数占有了。
由于MO这时是分数占有，每个MO只标出分数值，不再能标出MO能量值了。

[ Last edited by zhou2009 on 2012-5-21 at 07:58 ]
赞一下 回复此楼
5楼2012-05-19 08:47:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
gmy1990: 金币+2 2012-05-19 22:37:01
（接上）

从Population analysis看：
RHF是：
Population analysis using the SCF density.

而CCSD（加关键词density）是：
Population analysis using the QCI/CC density.
赞一下 回复此楼
6楼2012-05-19 11:36:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
gmy1990: 金币+2 2012-05-19 22:37:07
（再续上）

CCSD不加关键词density，只能输出RHF的波函数，但有CCSD的能量。

CCSD(T)如果加关键词density，会出错，不能运行。出错信息指明只有下列方法才能用关键词density：
MP2, MP3, MP4SDQ, CI, CCD, and QCI。（实际上可以作CCSD）

[ Last edited by zhou2009 on 2012-5-21 at 08:01 ]
赞一下 回复此楼
7楼2012-05-19 12:01:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kent1022

金虫 (著名写手)
引用回帖:
7楼: Originally posted by zhou2009 at 2012-05-19 12:01:40:
（再续上）

CCSD不加关键词density，只能作RHF。

CCSD(T)如果加关键词density，会出错，不能运行。出错信息指明只有下列方法才能用关键词density：
MP2, MP3, MP4SDQ, CI, CCD, and QCI。（实际上可以作CC ...

但是gaussian说明书CCSD关键字里的原话是：
The Coupled Cluster energy appears in the output as follows (following the final correlation iteration):
DE(CORR)=  -.54979226D-01        E(CORR)=     -.75019641794D+02
…
CCSD(T)=  -.75019717665D+02
The CCSD energy is labeled E(CORR), and the energy including the non-iterative triples contribution is given in the final line.

这样看来最后列出的CCSD(T)就是包含非迭代三重激发的CCSD能量？

而在 CCSD(T)之前给出的T4 ， T5 又是什么能量？
如
T4(CCSD)= -0.24124543D-01
T5(CCSD)=  0.54380474D-02
CCSD(T)= -0.44781728665D+03
赞一下 回复此楼
8楼2012-05-19 17:07:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
zhou2009: 金币+3 2012-05-19 19:41:24
....
3=-1.0711161|MP4D=-1.0761802|MP4DQ=-1.0752858|MP4SDQ=-1.0758967|CCSD=-
1.0801377|CCSD(T)=-1.0801377|RMSD=1.950e-010|PG=D*H [C*(H1.H1)]||@

其中的CCSD(T)=后面的就是CCSD(T)能量。
赞一下 回复此楼
9楼2012-05-19 19:04:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zpppanda111

金虫 (小有名气)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
zhou2009: 金币+2 2012-05-20 08:34:37
zhou2009: 编辑内容 2012-05-21 08:02
引用回帖:
9楼: Originally posted by sobereva at 2012-05-19 19:04:26:
....
3=-1.0711161|MP4D=-1.0761802|MP4DQ=-1.0752858|MP4SDQ=-1.0758967|CCSD=-
1.0801377|CCSD(T)=-1.0801377|RMSD=1.950e-010|PG=D*H ||@

其中的CCSD(T)=后面的就是CCSD(T)能量。

正解。CCSD(T)在gaussian09，A.02中可以得到能量。

[ Last edited by zhou2009 on 2012-5-21 at 08:02 ]
赞一下 回复此楼
10楼2012-05-19 23:07:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kent1022 的主题更新
171/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +7 wangmingqi 2026-08-28 8/400 2026-08-28 18:42 by zhuzg0628
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +11 Kittylucky 2026-08-27 12/600 2026-08-28 16:26 by lyxiaomc
[基金申请] 基金不中，共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
关闭
请填处理意见
关闭