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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 求助MP2算弱相互作用
161/2返回列表12下一页
查看: 1607  |  回复: 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的回帖

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的回帖

ghcacj

荣誉版主 (著名写手)

阿超

优秀版主

【答案】应助回帖

BSSE未必都是正的，首先确定BSSE后计算结果是否合理。
其次再研究是否是过渡金属带来的影响。
赞一下 回复此楼
5楼2011-04-25 14:44:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gtolv8688

铜虫 (小有名气)
引用回帖:
Originally posted by ghcacj at 2011-04-25 14:44:42:
BSSE未必都是正的，首先确定BSSE后计算结果是否合理。
其次再研究是否是过渡金属带来的影响。

我的结果其实都是考虑了BSSE之后的，因为我觉得不是BSSE的问题，所以在写帖子的时候就没有写。嗯，我修改一下原帖~
赞一下 回复此楼
Livewithproblem！
6楼2011-04-25 14:48:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gtolv8688

铜虫 (小有名气)
竟然不能修改原帖了。。。
补充一下，我计算吸附能的时候是考虑了BSSE修正的，但即使如此，发现由MP2计算得到的B.E.还是非常大的正值。而用其他方法，如PBEPBE/6-311+g*，或者WB97XD/6-311+g*算出来的都是合理的负值。请问mp2计算为什么会有问题呢？按理说MP2算弱相互作用应该更准才对啊？
赞一下 回复此楼
Livewithproblem！
7楼2011-04-25 14:55:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qzhaosdu

金虫 (著名写手)
你拆的那些cluster都优化没有？基组有没有问题？是什么金属？
赞一下 回复此楼
Anewday,anewhour,anewminute,anewpeople.
8楼2011-04-25 14:56:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gtolv8688

铜虫 (小有名气)
引用回帖:
Originally posted by qzhaosdu at 2011-04-25 14:56:59:
你拆的那些cluster都优化没有？基组有没有问题？是什么金属？

因为拆出来的小cluster一优化就散了，所以我把骨架原子给固定了，边缘用氢饱和，这些氢原子和放在cluster周围的氢分子可以动。主要是因为mp2无法算大cluster，所以只能这样处理了，我看到文献上也是这样做的，我说的文献跟我做的体系是不同的，我只是参考了里面的方法。我用的机组是cc-pvtz，金属原子是Fe。请问这样做有问题吗？
赞一下 回复此楼
Livewithproblem！
9楼2011-04-25 15:11:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qzhaosdu

金虫 (著名写手)
★ ★
ben_ladeng(金币+2): 谢谢交流 2011-04-27 06:21:47
引用回帖:
Originally posted by gtolv8688 at 2011-04-25 15:11:27:
因为拆出来的小cluster一优化就散了，所以我把骨架原子给固定了，边缘用氢饱和，这些氢原子和放在cluster周围的氢分子可以动。主要是因为mp2无法算大cluster，所以只能这样处理了，我看到文献上也是这样做的， ...

用H饱和的culster是否优化过？小的复合物是否优化过？Fe要用赝势基组，还有多重度一系列问题要考虑
赞一下(1人) 回复此楼
Anewday,anewhour,anewminute,anewpeople.
10楼2011-04-25 15:25:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gtolv8688 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 284求调剂 +9 archer.. 2026-04-10 9/450 2026-04-10 16:37 by 6781022
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 298求调剂 +8 钉叮咚冬瓜 2026-04-07 10/500 2026-04-10 11:25 by 高维春
[考研] 0702物理学学硕299求调剂 +6 祁柒连 2026-04-06 6/300 2026-04-10 11:10 by Roomoo
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 材料调剂 +11 一样YWY 2026-04-05 11/550 2026-04-10 09:32 by 钟洲2011
[考研] 275求调剂 +6 1624447980 2026-04-08 7/350 2026-04-10 08:06 by 探123
[考研] 一志愿211电子信息347求调剂 +3 554916 2026-04-03 4/200 2026-04-09 22:24 by 554916
[考研] 求调剂 +7 chenxrlkx 2026-04-05 9/450 2026-04-09 09:04 by wj165256
[考研] 一志愿哈工大，初试329，求环境科学与工程调剂！ +11 余未辛 2026-04-06 11/550 2026-04-08 15:21 by screening
[考研] 344求调剂 +11 魏子per 2026-04-07 11/550 2026-04-07 23:01 by JourneyLucky
[考研] 071000生物学，一志愿深圳大学296分，求调剂 +12 TIckLw 2026-04-06 13/650 2026-04-07 20:34 by lijunpoly
[考研] 专硕085403，291分，有两篇专利，一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
[考研] 材料求调剂 +18 一样YWY 2026-04-05 18/900 2026-04-07 15:49 by dxlg
[考研] 287分求调剂 有专利国奖一志愿哈工大085406 +6 白易辰 2026-04-06 7/350 2026-04-06 22:46 by 875465
[考研] 华南理工0703化学，总分336求调剂 +10 手机用户 2026-04-03 11/550 2026-04-06 22:20 by michael2011
[考研] 考研调剂 +3 Wwwwwww哇 2026-04-06 3/150 2026-04-06 20:55 by lbsjt
[考研] 332求调剂 +17 小小孟... 2026-04-05 18/900 2026-04-06 09:51 by 蓝云思雨
[考研] 08专硕275调剂 +5 AaAa7420 2026-04-05 5/250 2026-04-05 18:01 by jkddd
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
信息提示
关闭
请填处理意见
关闭