24小时热门版块排行榜    

查看: 1566  |  回复: 15

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!
已阅   回复此楼   关注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吗?直接算,误差大去了。
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要用赝势基组,还有多重度一系列问题要考虑
Anewday,anewhour,anewminute,anewpeople.
10楼2011-04-25 15:25:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gtolv8688 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 328求调剂,英语六级551,有科研经历 +4 生物工程调剂 2026-03-17 8/400 2026-03-21 02:12 by JourneyLucky
[考研] 279分求调剂 一志愿211 +11 chaojifeixia 2026-03-19 12/600 2026-03-21 01:49 by 星空星月
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 材料专硕英一数二306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 一志愿中南化学(0703)总分337求调剂 +8 niko- 2026-03-19 9/450 2026-03-20 21:57 by luoyongfeng
[考研] 265求调剂 +12 梁梁校校 2026-03-19 13/650 2026-03-20 21:01 by 无际的草原
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 260求调剂 +3 朱芷琳 2026-03-20 3/150 2026-03-20 20:35 by 学员8dgXkO
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
信息提示
请填处理意见