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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 求助MP2算弱相互作用
161/2返回列表12下一页
查看: 1510  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 4/200 2026-02-08 08:46 by vs90ilomwc
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 5/250 2026-02-08 08:32 by vs90ilomwc
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:26 by vs90ilomwc
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:07 by vs90ilomwc
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 08:06 by vs90ilomwc
[教师之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:52 by vs90ilomwc
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:46 by vs90ilomwc
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:32 by vs90ilomwc
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:27 by vs90ilomwc
[教师之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:26 by vs90ilomwc
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:07 by vs90ilomwc
[硕博家园] 博士延得我，科研能力直往上蹿 +8 偏振片 2026-02-02 8/400 2026-02-08 06:52 by liyeqik
[教师之家] 有院领导为了换新车，用横向课题经费买了俩车 +7 瞬息宇宙 2026-02-04 7/350 2026-02-07 21:47 by tfang
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[基金申请] 同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗 +4 hitsdu 2026-02-06 4/200 2026-02-07 13:07 by jurkat.1640
[基金申请] 有时候真觉得大城市人没有县城人甚至个体户幸福 +9 苏东坡二世 2026-02-04 10/500 2026-02-07 12:37 by 小毛球
[考博] 天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人） +4 a793625982 2026-02-05 5/250 2026-02-07 10:57 by a793625982
[公派出国] CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优 +4 雨念 2026-02-01 6/300 2026-02-06 23:32 by MelissaPon
[基金申请] 面上项目申报 +3 Tide man 2026-02-01 3/150 2026-02-05 22:56 by god_tian
[教师之家] 遇见不省心的家人很难过 +18 otani 2026-02-03 22/1100 2026-02-04 11:06 by tangmnt
信息提示
关闭
请填处理意见
关闭