24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1083  |  回复: 16

kent1022

金虫 (著名写手)

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

gaussian03中的ccsd(t)计算完成后,最后的能量输出汇总只显示到mp3能量,到09了还是这样,有什么方便的方法可以将每个单点能计算后的CCSD(T)= .......  提取出来?
回复此楼
已阅   回复此楼   关注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 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见