24小时热门版块排行榜    

查看: 4460  |  回复: 13

十年磨练

木虫 (正式写手)

[求助] 如何计算分子间的弱相互作用以及分析氢键 已有2人参与

如题 如何计算分子间的弱相互作用以及分析氢键


大家好,最近我在用gaussian软件计算分子间的相互作用,但是对于操作步骤现在还是不太确定是不是正确的,所以向各位大牛详细请教一下。
现在我把我的计算步骤详细说一下
1)比如两个分子,分别称为分子A和分子B,两个分子之间会形成氢键作用   
2)我用gaussian软件用B3LYP/6-311+G*基组分别计算了只有分子A存在的能量E(A)、只有分子B存在的能量E(B)以及同时存在两个分子(即分子A和分子B)的能量E(AB)        3)最后计算 认为E(AB)-E(A)-E(B)的绝对值是分子间的弱相互作用,
这样的计算方法对不对,还是说要加上BSSE修正E(BSSE)就行了  即|E(AB)-E(A)-E(B)|+E(BSSE)

请大家不吝赐教,在此谢谢大家了
回复此楼

» 收录本帖的淘帖专辑推荐

gaussian资源

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

时刻保持虚心学习的态度
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

kong890109

金虫 (小有名气)

★ ★
gmy1990: 金币+2, 感谢应助 2015-06-24 17:09:32
需要BSSE校正,具体的弱相互作用分析方法可以参看此贴http://muchong.com/html/201408/7827455.html;使用过这里面的RDG方法,感觉挺形象的,你可以试试。
2楼2015-06-22 20:06:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

十年磨练

木虫 (正式写手)

引用回帖:
3楼: Originally posted by zhou2009 at 2015-06-23 10:23:29
B3LYP算氢键并不那么合用(参见“乱谈DFT-D”:http://sobereva.com/83),如果起点不高,计算结果很难说。
不知分子体系有多大。最起码要用b97d。最好用DFT-D。如果氢键体系较小,可以用mp4(sdtq)或者CCSD。...

首先,谢谢你的热心解答,我再问一下
一般一个分子有八九个原子,不过有时候我会研究好几个分子之间的相互作用。这样的话,原子个数大概就是30~60个左右,对于这样的体系,选取mp2/aug-cc-pVDZ这样一个基组怎么样?计算时需不需要添加BSSE修正和色散修正?
时刻保持虚心学习的态度
4楼2015-06-23 10:36:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
十年磨练: 金币+5, 有帮助 2015-11-03 09:32:15
B3LYP算氢键并不那么合用(参见“乱谈DFT-D”:http://sobereva.com/83),如果起点不高,计算结果很难说。
不知分子体系有多大。最起码要用b97d。最好用DFT-D。如果氢键体系较小,可以用mp4(sdtq)或者CCSD。
3楼2015-06-23 10:23:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

十年磨练

木虫 (正式写手)

引用回帖:
2楼: Originally posted by kong890109 at 2015-06-22 20:06:50
需要BSSE校正,具体的弱相互作用分析方法可以参看此贴http://muchong.com/html/201408/7827455.html;使用过这里面的RDG方法,感觉挺形象的,你可以试试。

嗯,好的,我先看看帖子
时刻保持虚心学习的态度
5楼2015-06-23 10:37:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

★ ★
gmy1990: 金币+2, 感谢应助 2015-06-24 17:09:45
引用回帖:
4楼: Originally posted by 十年磨练 at 2015-06-23 10:36:26
首先,谢谢你的热心解答,我再问一下
一般一个分子有八九个原子,不过有时候我会研究好几个分子之间的相互作用。这样的话,原子个数大概就是30~60个左右,对于这样的体系,选取mp2/aug-cc-pVDZ这样一个基组怎么样 ...

原子个数不算太大。mp2/aug-cc-pVDZ应是可用的。但最好用Linux操作系统,避免数据溢出。
以我算氢键经验看,如果有蓝移氢键,mp2没有mp4(sdtq)好。
6楼2015-06-23 10:53:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

十年磨练

木虫 (正式写手)

引用回帖:
6楼: Originally posted by zhou2009 at 2015-06-23 10:53:43
原子个数不算太大。mp2/aug-cc-pVDZ应是可用的。但最好用Linux操作系统,避免数据溢出。
以我算氢键经验看,如果有蓝移氢键,mp2没有mp4(sdtq)好。...

哦,我就是用linux下的gaussian算的

另外怎么区分得到的结果是氢键能量和范德华相互作用等
时刻保持虚心学习的态度
7楼2015-06-23 13:58:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
gmy1990: 金币+2, 感谢应助 2015-06-24 17:09:51
十年磨练: 金币+10, ★★★很有帮助 2015-09-02 20:03:18
引用回帖:
7楼: Originally posted by 十年磨练 at 2015-06-23 13:58:47
哦,我就是用linux下的gaussian算的

另外怎么区分得到的结果是氢键能量和范德华相互作用等...

氢键是分子端头处于能级高位的电子(如HOMO)向能级低位的势阱空域偏移电子形成的。其中一方是LUMO之类,是空域。
范德华作用的双方通常都是配对的电子所形成的。(少见有一方是单电子)
8楼2015-06-23 16:00:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

十年磨练

木虫 (正式写手)

引用回帖:
8楼: Originally posted by zhou2009 at 2015-06-23 16:00:18
氢键是分子端头处于能级高位的电子(如HOMO)向能级低位的势阱空域偏移电子形成的。其中一方是LUMO之类,是空域。
范德华作用的双方通常都是配对的电子所形成的。(少见有一方是单电子)...

哦,我的意思是如何分别计算氢键和范德华相互作用的能量,或者算出来的能量如何区分多少是氢键多少是范德华相互作用能量   
对于你上面说的我不是很明白,我看看那方面的资料。对了,你有没有推荐的这方面的资料可以给我发一下吗? 谢谢了
时刻保持虚心学习的态度
9楼2015-06-25 09:25:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

809011916

铜虫 (小有名气)

【答案】应助回帖

引用回帖:
9楼: Originally posted by 十年磨练 at 2015-06-25 09:25:18
哦,我的意思是如何分别计算氢键和范德华相互作用的能量,或者算出来的能量如何区分多少是氢键多少是范德华相互作用能量   
对于你上面说的我不是很明白,我看看那方面的资料。对了,你有没有推荐的这方面的资料可 ...

是不是应该用波函数分析这个软件,比如multiwfn这个,通过rdg分析可以表征弱相互作用和氢键作用,你可以试试。
选择最喜欢的,专注到无情。一个,就够了。
10楼2015-07-24 01:45:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 十年磨练 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 01:30 by o3EMxHUAQDtS
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:12 by o3EMxHUAQDtS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见