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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 如何计算分子间的弱相互作用以及分析氢键
141/2返回列表12下一页
查看: 3938  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

十年磨练

木虫 (正式写手)

[求助] 如何计算分子间的弱相互作用以及分析氢键 已有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资源

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
时刻保持虚心学习的态度
1楼 2015-06-22 16:33:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

kong890109

金虫 (小有名气)
★ ★
gmy1990: 金币+2, 感谢应助 2015-06-24 17:09:32
需要BSSE校正，具体的弱相互作用分析方法可以参看此贴http://muchong.com/html/201408/7827455.html；使用过这里面的RDG方法，感觉挺形象的，你可以试试。
赞一下(2人) 回复此楼
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修正和色散修正？
赞一下(1人) 回复此楼
时刻保持虚心学习的态度
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的回帖
相关版块跳转 我要订阅楼主 十年磨练 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856求调剂 +5 zhn03 2026-03-25 6/300 2026-03-25 22:58 by peike
[考研] 303求调剂 +6 蓝山月 2026-03-25 6/300 2026-03-25 22:47 by 418490947
[考研] 材料与化工304求B区调剂 +3 邱gl 2026-03-25 3/150 2026-03-25 19:03 by Ainin_
[考研] 求b区院校调剂 +4 周56 2026-03-24 5/250 2026-03-25 17:12 by yishunmin
[考研] 282求调剂 +3 wcq131415 2026-03-24 3/150 2026-03-25 12:16 by userper
[考研] 303求调剂 +6 元夕元 2026-03-20 7/350 2026-03-25 12:00 by edmund7
[考研] 求调剂 一志愿 本科 北科大 化学 343 +4 13831862839 2026-03-24 5/250 2026-03-25 09:47 by 无际的草原
[考研] 求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分 +6 @taotao 2026-03-19 6/300 2026-03-25 08:37 by 木托莫露露
[考研] 食品专硕 一志愿双一流 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[考博] 申博26年 +4 八6八68 2026-03-19 4/200 2026-03-24 15:49 by 小Ben呵呵
[考研] 279分求调剂 一志愿211 +18 chaojifeixia 2026-03-19 20/1000 2026-03-24 10:34 by dolphin_ycj
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 352求调剂 +3 大米饭！ 2026-03-22 3/150 2026-03-22 23:28 by king123！
[考研] 材料与化工085600，总分304，本科有两篇sci参与，求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 一志愿南大，0703化学，分数336，求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
信息提示
关闭
请填处理意见
关闭