| 查看: 1543 | 回复: 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» 猜你喜欢
0703 调剂
已经有5人回复
-
上海华东理工大学机械与动力工程学院先进材料与制造研究组招聘公告
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有144人回复
-
304求调剂
已经有6人回复
-
西南某211招27届博士
已经有30人回复
-
大连工业杰青/长江团队-生物质材料-储能电池方向招收2026级博士生
已经有0人回复
-
304求调剂
已经有1人回复
-
一志愿985中药学306，求调剂
已经有0人回复
-
河南一本高校化学，材料，化工专业招收研究生
已经有0人回复
-
大连海事大学船舶洁净能源研究中心2026年博士研究生招生启事
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 2-氨基苄醇中的羟基做成氯或者溴，有什么条件
已经有16人回复
- 各位大侠，我可以优化结构，但只要一计算频率就出现下面的问题？
已经有4人回复
- 计算 Gaussian 03已经 弄坏两台电脑了~~~~~
已经有46人回复
- gaussian计算
已经有4人回复
- 求助一个计算标准偏差的方法
已经有12人回复
- 在Gaussian 09 计算中怎么做核冻结处理？
已经有6人回复
- 重奖征集国内国外做理论计算课题组的信息
已经有44人回复
6楼2018-01-23 11:39:09
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
