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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 求助MP2算弱相互作用
52/1返回列表
查看: 1388  |  回复: 15
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

gtolv8688

铜虫 (小有名气)

[求助] 求助MP2算弱相互作用

我最近在用gaussian09计算含有过渡金属的聚合物吸附小分子的问题，遇到一个很奇怪的问题。为了使用gaussian，我先把聚合物拆成几个小的cluster，对于每个cluster我用x3lyp/6-31g*优化看氢气分子吸附到什么位置上。找到吸附位以后，进一步将cluster拆成更小的cluster，然后用MP2/cc-pvtz计算B.E.（B.E.=EAB-EA-EB，其中A、B分别为cluster和H2）。结果计算得到的吸附能是正的！这意味着氢气分子和cluster之间是相互排斥的？我又试了其他的优化方法，比如PBEPBE/6-31g*，发现氢分子仍然吸附在用X3Lyp优化的位置附近。这意味着优化没有什么问题啊，可是MP2计算得到的B.E.>0是什么原因呢？我搜索了一下过去的帖子，看到有人说MP2算某些过渡金属结果不好，但没有过多的解释为什么不好。在我这个问题中是否也是因为这个原因导致MP2结果出问题了呢？我刚刚接触gaussian没有多久，属于十分菜的菜虫，希望版上各位老师指点迷津啊！~
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
Livewithproblem！
1楼2011-04-25 14:04:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gtolv8688

铜虫 (小有名气)
引用回帖:
Originally posted by qzhaosdu at 2011-04-27 11:12:09:
建议楼主你不要先用大的cluster去做了，要是资源够就直接大的去算，要是不够就用小的去算。
还有，用DFT就行了（带D的），MP2就算了。

我昨天算了一下，含Fe的cluster基态自旋是3，果然是多重度的问题没有考虑~我倒是用DFT-D里的WB97XD方法优化过结构，但是不知道算能量用WB97XD是否会被认可，MP2肯定没问题，因为要拟合力场，所以精度还是很重要的。不知道为什么送不了金币给你，就送一朵大红花吧
赞一下 回复此楼
Livewithproblem！
15楼2011-04-27 12:04:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主

【答案】应助回帖

★ ★
gtolv8688(金币+1): 感谢热心回复~ 2011-04-26 10:21:17
ben_ladeng(金币+2): 谢谢版主 2011-04-27 06:22:33
（B.E.=EAB-EA-EB，其中A、B分别为cluster和H2）。
看到这里我就笑了，难道MP2计算不需要BSSE吗？直接算，误差大去了。
赞一下(1人) 回复此楼
2楼2011-04-25 14:20:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gtolv8688

铜虫 (小有名气)
引用回帖:
Originally posted by ghcacj at 2011-04-25 14:20:11:
（B.E.=EAB-EA-EB，其中A、B分别为cluster和H2）。
看到这里我就笑了，难道MP2计算不需要BSSE吗？直接算，误差大去了。

先别笑啊郭版主，我是考虑了BSSE的，只是写帖子的时候没有说明。现在的问题是不考虑BSSE，B.E.都大于0了，再加上BSSE那就更大了。用MP2计算显示吸附并不稳定，但是用其他方法优化的结果显示应该是吸附上的，而且计算B.E.也是小于0的。MP2计算的结果有明显的硬伤，那就是一直要排斥到4点几个埃才会有一点点弱的吸附作用，很显然是错的。MP2为什么结果不对呢？算苯环的时候结果是非常好的，是不是因为有过渡金属的缘故呢？
赞一下 回复此楼
Livewithproblem！
3楼2011-04-25 14:39:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gtolv8688

铜虫 (小有名气)
引用回帖:
Originally posted by ghcacj at 2011-04-25 14:20:11:
（B.E.=EAB-EA-EB，其中A、B分别为cluster和H2）。
看到这里我就笑了，难道MP2计算不需要BSSE吗？直接算，误差大去了。

而且就我目前的计算来看，BSSE并没有大到影响定性结果的程度，有很多文献在用MP2算吸附能的时候竟然都没考虑BSSE。。。
赞一下 回复此楼
Livewithproblem！
4楼2011-04-25 14:42:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭