| 查看: 1561 | 回复: 5 | |||||
beefly专家顾问 (职业作家)
地沟油冶炼专家
|
[交流]
GAMESS的gMCQDPT做RASPT2计算可能存在BUG 已有4人参与
|
|
对于大活性空间的RASPT2型计算，GAMESS的gMCQDPT模块（$GMCPT group & MRPT=GMCPT in $MRMP group）可能存在BUG。 最早怀疑程序有问题，是三年前算CuCl2的激发能。活性空间是21电子占据17轨道（Cu:3d4s4d + Cl:3p + Cl:3p）。实验值，以及Molpro的RASPT2计算结果，见 J. Chem. Phys. 130, 154313, 2009，表3。用同样的活性空间，GAMESS的RASSCF（$GMCPT group & CISTEP=GMCCI in $MCSCF group）可以得到与Molpro非常接近的参考态绝对能量（未在表3列出），但是加上相关能的RASPT2（gMCQDPT）结果非常差。 最近发现的例子是FeO-。最新一期的JCTC发表了RASPT2级别的计算结果（J. Chem. Theory Comput., 7, 3961, 2011，表9）。其实早在一年前，我用Molpro得到了类似的RASPT2结果，但是和GAMESS的gMCQDPT相差太远。 程序： GAMESS 12.JAN.2009(R3)； Molpro 2009.1作为参考 二者都测试了主页提供的二进制程序，和自己从源码编译的程序，结果没有差别 计算方法： 活性空间：13电子占据17轨道（Fe:3d4s4d + O:2p3p）。一般来说，O:3p不是必须的，但是对FeO-来说，Fe:4d和O:3p的能量非常接近以至于无法分开，只好都放入活性空间。在J. Phys. Chem. A, 113, 10588 2009这篇文章中，就是因为没有加上O:3p，导致Fe:4d的某些分量被O:3p替代，因此对FeO-的强关联描述不够准确，结果和实验不符。 组态限制：Fe:4d + O:3p只允许0-2个电子占据 基组：DKH2标量相对论ANO-RCC（收缩到VQZ） 排除入侵态： Molpro的RASPT2：ipea=0.25,shift=0.2 GAMESS的gMCQDPT：EDSHFT=0.03 结果： 这里只列出键长＝1.64 Ang.的单点能 RASSCF MOLPRO GAMESS 6Sigma+ -1346.4359457 -1346.4359467 Hartree 4Delta -1346.4300975 -1346.4300976 Hartree RASSCF结果完全一致。 RASPT2 MOLPRO GAMESS 6Sigma+ -1347.1556336 -1347.1554413 Hartree 4Delta -1347.1547722 -1347.1682683 Hartree Tv 0.023 -0.349 eV MOLPRO遇到入侵态会自动停止计算；GAMESS的参考态权重都在89%-90%之间。可见以上结果不存在入侵态问题。 6Sigma+态的总能量比较接近，但是4Delta态总能量相差较大。gMCQDPT模块严重高估了4Delta态的相关能，导致垂直激发能从0.023eV一下降到-0.349eV。可供参考的值为文献JCTC, 7, 3961, 2011中给出的RASPT2激发能0.02 eV（用MOLCAS7开发版做的，键长略有不同），和MOLPRO结果一致。 不同的RASPT2方法得到不同结果是正常的，但是gMCQDPT误差太大，连定性结果都得不到。注意GAMESS的RASSCF计算也使用gMCQDPT模块，可见用这个模块做RASSCF计算还是可靠的，但是在BUG排除之前，不建议用于RASPT2级别的计算。 [ Last edited by beefly on 2011-12-14 at 20:59 ] |
回复此楼» 本帖已获得的红花（最新10朵）
calfwoo» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有188人回复
-
调剂招生
已经有0人回复
-
莆田学院环境与生物工程学院资源与环境专硕（环境工程方向） 开放调剂招生
已经有12人回复
-
浙江海洋大学石油化工与环境学院招收工科调剂生（学硕、专硕均可，线上面试）
已经有0人回复
-
浙江海洋大学石油化工学院接收工科调剂（学硕、专硕均可，近期线上面试）
已经有8人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有1人回复
-
调剂招生
已经有0人回复
-
浙江海洋大学石油与天然气工程专业接收工科调剂（13日早上截止，线上面试）
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 2-氨基苄醇中的羟基做成氯或者溴，有什么条件
已经有16人回复
- 各位大侠，我可以优化结构，但只要一计算频率就出现下面的问题？
已经有4人回复
- 计算 Gaussian 03已经 弄坏两台电脑了~~~~~
已经有46人回复
- gaussian计算
已经有4人回复
- 求助一个计算标准偏差的方法
已经有12人回复
- 在Gaussian 09 计算中怎么做核冻结处理？
已经有6人回复
- 重奖征集国内国外做理论计算课题组的信息
已经有44人回复
calfwoo
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1812.1
- 红花: 2
- 帖子: 171
- 在线: 276.9小时
- 虫号: 97355
- 注册: 2005-11-08
- 性别: GG
- 专业: 理论和计算化学
2楼2011-12-14 16:22:10
Jasminer
铁杆木虫 (著名写手)
- 应助: 52 (初中生)
- 金币: 7407.8
- 散金: 1845
- 红花: 27
- 帖子: 1002
- 在线: 728.8小时
- 虫号: 620156
- 注册: 2008-10-08
- 性别: GG
- 专业: 应用有机化学
3楼2011-12-14 20:46:27
calfwoo
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1812.1
- 红花: 2
- 帖子: 171
- 在线: 276.9小时
- 虫号: 97355
- 注册: 2005-11-08
- 性别: GG
- 专业: 理论和计算化学
4楼2011-12-15 14:38:44
dongwen176
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1831.3
- 红花: 1
- 帖子: 179
- 在线: 195.6小时
- 虫号: 91819
- 注册: 2005-09-07
- 性别: GG
- 专业: 理论和计算化学
5楼2011-12-17 10:13:01
6楼2018-01-23 11:39:09
5