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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 求助MP2算弱相互作用
161/2返回列表12下一页
查看: 1604  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +4 璃茉一定上岸 2026-04-10 4/200 2026-04-10 13:26 by wp06
[考研] 344求调剂 +7 丶风雪夜归人丶 2026-04-09 7/350 2026-04-10 12:05 by pengliang8036
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 化工调剂求导师收留！一志愿失利，踏实肯干，有植物提取科研经历 +17 yzyzx 2026-04-09 17/850 2026-04-10 10:29 by wp06
[考研] 考研调剂 +22 硕星赴 2026-04-09 23/1150 2026-04-10 09:43 by liuhuiying09
[考研] 297求调剂 +27 GENJIOW 2026-04-07 30/1500 2026-04-09 23:20 by wolf97
[考研] 初试分332，一志愿报考西北工业大学， +11 故人?? 2026-04-09 11/550 2026-04-09 21:54 by JineShine
[考研] 材料专硕调剂 +16 哈哈哈吼吼吼哈 2026-04-07 17/850 2026-04-09 21:16 by wutongshun
[考研] 材料考研求调剂总分280 +30 mkjlz1 2026-04-06 35/1750 2026-04-08 21:25 by cyh—315
[考研] 336求调剂，一志愿中科大 +9 墨彧 yuyu 2026-04-06 9/450 2026-04-08 11:24 by 想读书的菌菌
[考研] 265求调剂 +19 小木虫085600 2026-04-06 21/1050 2026-04-08 10:38 by 逆水乘风
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[考研] 284求调剂 +7 徐同学_001 2026-04-04 13/650 2026-04-05 17:19 by yulian1987
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-04 5/250 2026-04-05 14:06 by imissbao
[考研] 材料化工306分找合适调剂 +14 沧海轻舟e 2026-04-04 14/700 2026-04-05 09:53 by 朱云虎202
[考研] 341求调剂 +3 学无止境，冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
[考研] 材料调剂 +7 dxy调剂 2026-04-04 7/350 2026-04-05 09:15 by 陌秋26
[考研] 一志愿北京化工大学，初试成绩350求调剂 +9 沿岸?贝壳 2026-04-04 14/700 2026-04-05 01:09 by 沿岸?贝壳
[考研] 考研调剂 +5 小sun要好运 2026-04-03 5/250 2026-04-03 21:43 by 啵啵啵0119
信息提示
关闭
请填处理意见
关闭