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

24小时热门版块排行榜    

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

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 的主题更新
172/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +11 Kittylucky 2026-08-27 12/600 2026-08-28 16:26 by lyxiaomc
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 基金未中，这种答复是模板吗？ +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
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 出来了 +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
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +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
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+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
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
关闭
请填处理意见
关闭