计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 透彻认识氢键本质、简单可靠地估计氢键强度:一篇2019年JCC上的重要研究文章介绍

  [复制链接 Copy URL]

84

帖子

11

威望

4702

eV
积分
5006

Level 6 (一方通行)

16#
发表于 Post on 2019-10-13 17:14:48 | 只看该作者 Only view this author
谢谢分享,对于理解弱相互作用很有帮助!

84

帖子

11

威望

4702

eV
积分
5006

Level 6 (一方通行)

17#
发表于 Post on 2020-1-10 19:41:15 | 只看该作者 Only view this author
请教卢老师:帖子中,第二幅能量分解的柱状图,是用什么作图软件实现这种层叠的柱状图效果?

6万

帖子

99

威望

5万

eV
积分
120080

管理员

公社社长

18#
 楼主 Author| 发表于 Post on 2020-1-10 22:12:14 | 只看该作者 Only view this author
Yaqi 发表于 2020-1-10 19:41
请教卢老师:帖子中,第二幅能量分解的柱状图,是用什么作图软件实现这种层叠的柱状图效果?

那个图是合作者Saeed画的,应该用的是Origin
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

84

帖子

11

威望

4702

eV
积分
5006

Level 6 (一方通行)

19#
发表于 Post on 2020-1-11 20:04:39 | 只看该作者 Only view this author
sobereva 发表于 2020-1-10 22:12
那个图是合作者Saeed画的,应该用的是Origin

好的,我试一下,谢谢老师

1187

帖子

5

威望

2859

eV
积分
4146

Level 6 (一方通行)

20#
发表于 Post on 2020-2-27 23:47:46 | 只看该作者 Only view this author
请问比如氨基酸水溶液中两性离子,分子内羧基负离子和氨基阳离子之间的氢键,适合用这个BCP处电子密度的方法估计么,两边有明显电荷,用中性氢键体系预测公式是否合理呢?谢谢

6万

帖子

99

威望

5万

eV
积分
120080

管理员

公社社长

21#
 楼主 Author| 发表于 Post on 2020-2-28 03:05:24 | 只看该作者 Only view this author
snljty 发表于 2020-2-27 23:47
请问比如氨基酸水溶液中两性离子,分子内羧基负离子和氨基阳离子之间的氢键,适合用这个BCP处电子密度的方 ...

文中预测中性氢键键能的公式,对于氢键作用位置明确的(不显著牵扯二级作用等情况)的中性氢键很普适
但是你的情况,COO-上负电荷同时明显分布在两个氧上,比较复杂,原理上肯定不可能光靠一个位置的BCP靠谱地预测氢键作用。应当用下文里其它做法计算
计算分子内氢键键能的几种方法
http://sobereva.com/522
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1187

帖子

5

威望

2859

eV
积分
4146

Level 6 (一方通行)

22#
发表于 Post on 2020-2-28 08:25:12 | 只看该作者 Only view this author
sobereva 发表于 2020-2-28 03:05
文中预测中性氢键键能的公式,对于氢键作用位置明确的(不显著牵扯二级作用等情况)的中性氢键很普适
但 ...

了解,谢谢老师!

1043

帖子

0

威望

4174

eV
积分
5217

Level 6 (一方通行)

23#
发表于 Post on 2020-7-30 00:06:47 | 只看该作者 Only view this author
有点觉得B3LYP-D3BJ优化很弱的氢键复合物(如结合能只有-0.5 kcal/mol左右)似乎未必可靠吧(毕竟 B3LYP不胜任弱相互作用,且DFT-D3属于经验方法)?由于体系不是很大,为什么不用CCSD(T)甚或MP2优化个别体系测试一下呢?

6万

帖子

99

威望

5万

eV
积分
120080

管理员

公社社长

24#
 楼主 Author| 发表于 Post on 2020-7-30 03:09:24 | 只看该作者 Only view this author
granvia 发表于 2020-7-30 00:06
有点觉得B3LYP-D3BJ优化很弱的氢键复合物(如结合能只有-0.5 kcal/mol左右)似乎未必可靠吧(毕竟 B3LYP不 ...

B3LYP-D3(BJ)算弱相互作用能精度已经不错,而优化对于计算级别要求明显更低,所以这级别就够了。我和合作者大量的计算经验、前人的文章都可以证明这个级别对于优化完全够用了。
另外更关键的是,当前文章给出的拟合关系是希望其他人能拿来解决实际问题、用于预测他们自己体系的作用能。若比如当前文章用CCSD/def2-TZVP做优化并且得到波函数信息,虽然名义上结构和电子密度质量会更准确一点点,但这导致拟合公式缺乏实用性,毕竟稍微大一点的体系拿这个级别都根本算不动。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1043

帖子

0

威望

4174

eV
积分
5217

Level 6 (一方通行)

25#
发表于 Post on 2020-7-30 12:29:47 | 只看该作者 Only view this author
sobereva 发表于 2020-7-30 03:09
B3LYP-D3(BJ)算弱相互作用能精度已经不错,而优化对于计算级别要求明显更低,所以这级别就够了。我和合作 ...

我只是觉得对于特别弱的氢键体系,如结合能只有-0.6kcal/mol的H2O-CH4体系,DFT再精确也无法达到这样的精度吧。要知道,CCSD(T)的单点计算都是基于DFT势能面的基础上的,所以应该对于这类个别体系用CCSD(T)验证一下平衡构型是否可靠,毕竟小于1kcal的能量差都可视为柔性体系了,与分子轻微的相对转动能量相当了吧

6万

帖子

99

威望

5万

eV
积分
120080

管理员

公社社长

26#
 楼主 Author| 发表于 Post on 2020-7-31 01:23:48 | 只看该作者 Only view this author
granvia 发表于 2020-7-30 12:29
我只是觉得对于特别弱的氢键体系,如结合能只有-0.6kcal/mol的H2O-CH4体系,DFT再精确也无法达到这样的精 ...

B3LYP-D3(BJ)优化这个不会有什么问题
就算因为势能面平缓,可能对这个体系B3LYP-D3(BJ)结构误差比其它的稍微大一丁点,但由于本身势能面平缓使得高级别算的作用能对结构不敏感,所以用当前计算级别也没什么问题
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1043

帖子

0

威望

4174

eV
积分
5217

Level 6 (一方通行)

27#
发表于 Post on 2020-7-31 12:11:34 来自手机 | 只看该作者 Only view this author
sobereva 发表于 2020-7-31 01:23
B3LYP-D3(BJ)优化这个不会有什么问题
就算因为势能面平缓,可能对这个体系B3LYP-D3(BJ)结构误差比其它的 ...

明白了,多谢解答

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

28#
发表于 Post on 2022-4-7 21:18:45 | 只看该作者 Only view this author
本帖最后由 牧生 于 2022-4-7 21:50 编辑

请教这样一个让我困惑的疑问,以具体实例进行问题描述:

1、以冰晶胞和水,少量NaCl和CO2,低温高压得到水合物以后,取某个水分子及附近的几个水分子为A体系,该体系净电荷为0

2、取CL原子附近的几个分子为B体系,该体系净电荷为-1

3、删掉虚原子,整理好原子顺序,使用orca得到波函数信息(r2scan-3c,极为简单的体系,应该足够)

4、A体系和B体系的IGMH图分别如下
   

5、 从上图中看出,A体系中,四个等值面看起来都一样的。而在B体系中,一个氯原子与周围5个水分子形成了5个氢键,其中三个等值面都是蓝色,为典型的氢键作用;另外两个氢键以绿色为主。那么有个疑问,这两个绿色的部分,由于分子间距离和角度符合氢键的特征,被归为氢键,但等值面又是绿色的,也可以说是范德华相互作用。
问题1:如果非要进行一个总结,这两个绿色等值面的地方,到底应该归结为范德华相互作用,还是氢键呢?

问题2:我该如何定量去描述这几个氢键(a~i)的强度,如何进行对比。

自我发挥的方法如下:
如果都按照体系为电中性计算,参考此文http://sobereva.com/522 套用公式,得到的结果如下
A (neutral)
B (charged)
a
b
c
d
e
f
g
h
i
-9.36
-6.61
-6.99
-10.1
-6.52
-6.77
-3.98
-2.26
-6.05
a~d的键能有所区别,是因为水分子有部分参与形成五边形,部分参与形成六边形,所以键能不一样,这个我能理解。但是对于e~i,虽然数据值看起来合理,但是由于体系带电,上表中直接套用不带电的公式得到的数值应该是不可靠的??
如果套用带电的那个公式,那么得到的值则明显低的多,表明键的强度更大,如下表
B (charged)
e
f
g
h
i
-11.88970222
-12.25455046
-8.102392287
-5.531298989
-11.18637682

这就出现了我不能理解的矛盾。等值面颜色来看,蓝色的强度>绿色,但套用带电的那个公式,数值明显不靠谱。
如果用CVB,得到的结果是
neutral
charged
a
b
c
d
e
f
g
h
i
-0.075187
-0.033891
-0.032233
-0.059152
-0.068066
-0.069687
-0.016576
0.028026
-0.059147


归根结底,我该如何用数值,定量的去区分这几个氢键的强度呢?什么时候用CVB,什么时候用结合能呢?



又菜又爱玩

6万

帖子

99

威望

5万

eV
积分
120080

管理员

公社社长

29#
 楼主 Author| 发表于 Post on 2022-4-8 16:02:25 | 只看该作者 Only view this author
牧生 发表于 2022-4-7 21:18
请教这样一个让我困惑的疑问,以具体实例进行问题描述:

1、以冰晶胞和水,少量NaCl和CO2,低温高压得到 ...

g和i应当算氢键,通过结构和常识就能判断
几个氢键相对强弱可以用BCP的密度讨论
当前涉及到带电和不带电的情况对比,具体定量强度应当用高精度量子化学方法算相互作用能来得到和对比
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1478

帖子

0

威望

4537

eV
积分
6015

Level 6 (一方通行)

30#
发表于 Post on 2022-4-10 15:37:44 | 只看该作者 Only view this author
本帖最后由 牧生 于 2022-4-10 15:48 编辑
sobereva 发表于 2022-4-8 16:02
g和i应当算氢键,通过结构和常识就能判断
几个氢键相对强弱可以用BCP的密度讨论
当前涉及到带电和不带 ...

重新进行了计算,使用orca得到输入文件。为了表明都是同一个级别,且为了照顾阴离子体系B,所以,都选用了-2加了弥散函数,直接选择了7b得到orca的输入文件,得到的结果好一些,具体如下:
neutral   直接套用不带电公式
charged 直接套用带电公式
a
b
c
d
e
f
g
h
i
-11.80013081
-7.454337147
-8.103973693
-9.889485485
-11.83344969
-12.17984292
-11.14618356
-8.002261914
-5.418346706

数值上看起来比之前合理多了,但IGMH图形上出现了我不好解释的地方,以体系A为例,如图所示


同一个等值面,从外向里看都是蓝色,从里往外看有很明显的红圈,这个该如何解释呢?


以下为计算5个水分子的inp文件,请帮忙看一下,精度是否足够了。

! wB97X-2 D3 ma-def2-QZVPP autoaux RIJCOSX tightSCF noautostart miniprint nopop
%maxcore  1000
%pal nprocs   4 end
# DFT-D3(BJ) parameter for wB97X-2, see SI of PCCP, 20, 23175 (2018)
%method
D3S6 0.547
D3A1 3.520
D3S8 0.0
D3A2 7.795
end
* xyz   0   1
O     13.45000000    9.69000000   36.79000000
H     13.69000000    8.86000000   37.20000000
H     13.87000000   10.36000000   37.34000000
O     16.35000000    7.55000000   34.63000000
H     16.29000000    6.96000000   33.87000000
H     17.26000000    7.85000000   34.62000000
O     14.76000000    9.53000000   34.55000000
H     14.25000000    9.39000000   35.35000000
H     15.47000000    8.89000000   34.61000000
O     16.24000000   11.82000000   34.48000000
H     15.83000000   12.34000000   33.78000000
H     15.69000000   11.05000000   34.56000000
O     13.45000000    9.48000000   32.21000000
H     13.81000000    9.48000000   33.09000000
H     12.52000000    9.29000000   32.32000000
*


又菜又爱玩

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-13 03:06 , Processed in 0.261439 second(s), 22 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list