24小时热门版块排行榜    

查看: 4428  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 十年磨练 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
请填处理意见