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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 求助MP2算弱相互作用
161/2返回列表12下一页
查看: 2093  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 两块石头 +7 阿美_Lml888 2026-09-03 9/450 2026-09-06 17:48 by myrtle
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！ +5 SCSIOyxguo 2026-09-03 8/400 2026-09-06 15:06 by 可爱小花朵
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:39 by E1iiBLMU2XnD
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:19 by E1iiBLMU2XnD
[基金申请] 面上没中，邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[基金申请] 学科评审组评审是指会评吗？ +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
关闭
请填处理意见
关闭