24小时热门版块排行榜    

查看: 1092  |  回复: 16

zpppanda111

金虫 (小有名气)

引用回帖:
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中可以得到能量。
11楼2012-05-19 23:07:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zpppanda111

金虫 (小有名气)

引用回帖:
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中可以得到能量。
12楼2012-05-19 23:08:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
综上帖子:

在G03、G09,计算能量、平衡距离可以用CCSD(T)。应是常规计算。
但这个计算只能给出能量,不能给出波函数(包括衍生信息,如电子密度)。

由于CCSD(T)方法不支持生成密度信息,和density关键词一起用会出错。这时,可以用MP4SDQ或者CCSD来得到密度。
如需要作密度差时,加关键词density ,可用out=wfn产生wfn文件。

(CCSD(T)是在CCSD的基础上通过微扰的方法获得比CCSD更高一级能量的方法,只能得到能量而不能提供波函数、密度信息。理论本身就不能给出波函数,并不是因为Gaussian的限制。)

[ Last edited by zhou2009 on 2012-5-20 at 20:54 ]
13楼2012-05-20 10:08:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zpppanda111

金虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
13楼: Originally posted by zhou2009 at 2012-05-20 10:08:45:
综上帖子:

在G03、G09,计算能量、平衡距离可以用CCSD(T)。应是常规计算。
但这个计算只能给出能量,不能给出波函数(包括衍生信息,如电子密度)。

由于CCSD(T)方法不支持生成密度信息,和density关键词 ...

这个问题值得讨论,考虑到CCSD(T)与密度信息的关系,我做了下测试比较,一是使用MP2下优化得到的fchk进行AIM分析,
另一个就是在MP2下得到的优化结构基础上,进行CCSD(T)单点计算,然后用CCSD(T)下的fchk进行AIM分析,
两种情况下得到的结果不同,也可以说是电子密度p不同,

关于版主说的“CCSD不加关键词density,只能作RHF。”
是否意味着无论是MP2还是在其基础上的CCSD(T)得到的电子密度都是RHF水平?
14楼2012-05-20 23:44:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
是的,CCSD(T)下的chk、fchk,里面的记录都是RHF的,得到的电子密度都是RHF的。CCSD(T)没有自己的电子密度输出,它在方法上就没有波函数,当然也没有电子密度。
因此不能对CCSD(T)的结果作电子密度分析,这时分析的只是RHF的。

[ Last edited by zhou2009 on 2012-5-21 at 08:05 ]
15楼2012-05-21 07:55:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zpppanda111

金虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
15楼: Originally posted by zhou2009 at 2012-05-21 07:55:29:
是的,CCSD(T)下的chk、fchk,里面的记录都是RHF的,得到的电子密度都是RHF的。CCSD(T)没有自己的电子密度输出,它在方法上就没有波函数,当然也没有电子密度。
因此不能对CCSD(T)的结果作电子密度分析,这时分析 ...

请看这两篇文献:

J. Chem. Phys. 130, 164102 (2009)  
page 164102-8 最后一句话:We also thank Andrew Tale for providing CCSD(T) electron densities for the Ne atom.

Chemical Physics Letters 298 1998. 113–119.  
page 118, Once again, only the CCSD(T) method provides electron densities of sufficient quality [24].
16楼2012-05-21 10:23:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
16楼: Originally posted by zpppanda111 at 2012-05-21 10:23:58:
请看这两篇文献:

J. Chem. Phys. 130, 164102 (2009)  
page 164102-8 最后一句话:We also thank Andrew Tale for providing CCSD(T) electron densities for the Ne atom.

Chemical Physics Letters 29

我请教了一下网友:
从一般的微扰计算,比如MP2,获得MP2的波函数比起获得MP2的能量要有更多的计算过程、费更多计算时间。计算MP2的能量只需要有HF级别的波函数就能直接通过微扰公式得到。
想来也是类似的,一般做的CCSD(T)只是在CCSD基础上通过微扰方式计算能量而已。CCSD(T)的波函数虽然也能获得,但是获得的话比起计算CCSD(T)的能量要花费明显更多的计算过程、更多的时间。
JCP文中最后那句话特意指出感谢Tale提供CCSD(T)密度,这也表明CCSD(T)级别的密度是不容易得到的,得到Ne这么小体系的密度还得感谢别人提供,而作者自己难以算出来。
看来,应该说是Gaussian没有包括CCSD(T)的波函数计算而已。
17楼2012-05-21 16:40:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kent1022 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见