计算化学公社

 找回密码 Forget password
 注册 Register
Views: 16313|回复 Reply: 20
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 关于是否形成氢键的讨论

[复制链接 Copy URL]

154

帖子

0

威望

1173

eV
积分
1327

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
老师您好,我计算的醇类+HO2的反应势能面当中,在-OH位的氢提取反应以及alpha位的氢提取反应当中,从反应物出来过后在进入TS之前,都会有一个比较深的well生成,我就想是不是因为形成了氢键,所以想请教一下有没有什么方法可以计算出氢键的存在与否?谢谢老师

1043

帖子

0

威望

4174

eV
积分
5217

Level 6 (一方通行)

2#
发表于 Post on 2018-9-10 22:22:18 来自手机 | 只看该作者 Only view this author
简单的话用nbo看下有无donor-acceptor作用,也就是占据的Lewis NBO与空的non-Lewis NBO之间的作用。可以看下相应的二阶微扰近似作用能

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

3#
发表于 Post on 2018-9-11 05:11:28 | 只看该作者 Only view this author
NBO分析氢键作用是严重误导性的,极为不可取。应该用的分析方法见下
Multiwfn支持的弱相互作用的分析方法概览
http://sobereva.com/252
末尾专门批判了用NBO讨论这种问题的做法,很多文章里都在瞎用滥用NBO,明显体现出作者对NBO一知半解或者受到其它文章误导而跟风

分析氢键存在性首先应当考虑用Multiwfn做上面博文中提到的AIM、RDG、IGM等分析,十分方便,特别是RDG、IGM图像非常直观
北京科音自然科学研究中心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

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

4#
发表于 Post on 2018-9-11 05:14:58 | 只看该作者 Only view this author
granvia 发表于 2018-9-10 22:22
简单的话用nbo看下有无donor-acceptor作用,也就是占据的Lewis NBO与空的non-Lewis NBO之间的作用。可以看 ...

这是极为不妥的。NBO根本都没有把氢键作用的本质正确表现出来,甚至就连看似能被E2所描述的CT作用都没法恰当描述,被Stone在文章里大批一番,作者撰文试图反击,结果又被Stone批回去了
北京科音自然科学研究中心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

154

帖子

0

威望

1173

eV
积分
1327

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2018-9-11 08:23:50 | 只看该作者 Only view this author
谢谢

154

帖子

0

威望

1173

eV
积分
1327

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2018-9-11 08:24:20 | 只看该作者 Only view this author
sobereva 发表于 2018-9-11 05:14
这是极为不妥的。NBO根本都没有把氢键作用的本质正确表现出来,甚至就连看似能被E2所描述的CT作用都没法 ...

谢谢老师的解答,因为我还从来没用过Multiwfn,现在又急着计算分析一下结果,请问老师Multiwfn的手册里面或者是您的博文当中有没有介绍IGM分析的例子让我可以模仿着算一下然后分析一下结果?谢谢老师

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

7#
发表于 Post on 2018-9-11 09:09:52 | 只看该作者 Only view this author
XZDZ 发表于 2018-9-11 08:24
谢谢老师的解答,因为我还从来没用过Multiwfn,现在又急着计算分析一下结果,请问老师Multiwfn的手册里面 ...

仔细看我博文里第2条,写得清清楚楚,相关博文都给了,博文里的操作和原理介绍写得没法更明白、易懂。一张图RDG或IGM图撑死了一分钟就画出来了
北京科音自然科学研究中心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 (一方通行)

8#
发表于 Post on 2018-9-11 09:33:59 | 只看该作者 Only view this author
sobereva 发表于 2018-9-11 05:14
这是极为不妥的。NBO根本都没有把氢键作用的本质正确表现出来,甚至就连看似能被E2所描述的CT作用都没法 ...

我记得Weinhold的书《Valency and Bonding-a NBO donor-acceptor perspetive》中专门将了氢键,而且一开始他也利用了点电荷的静电模型计算了水二聚体的氢键能,结果跟实验值是在一个数量级上的,说明Weinhold并不否认氢键的主要成因是静电作用。  但接下来,他完全从donor-acceptor的perspetive去分析一些典型氢键,我个人觉得这种轨道作用角度的分析makes sense,至少能解释氢键的取向性,能够解释不同体系的氢键相对强度(对于NBO的E2高估氢键能的批评,我认为那是因为NBO的轨道作用能没考虑到分子片段的结构构型弛豫以及分子片段的轨道弛豫,E2有点类似EDA说的instantaneous作用能。另外,E2而且也没有考虑Pauli排斥能,因为在正交化NAO时已经去掉电子对间的Pauli排斥了)。 我想问的是,有没有可能NBO的donor-acceptor的perspetive在某种程序上包含了静电作用呢?

当然,我还没下Stone的文章。只扫了摘要,他似乎说NBO的问题出在BSSE?

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

9#
发表于 Post on 2018-9-11 09:56:50 | 只看该作者 Only view this author
granvia 发表于 2018-9-11 09:33
我记得Weinhold的书《Valency and Bonding-a NBO donor-acceptor perspetive》中专门将了氢键,而且一开 ...

E2能解释相对强度,有相关性,达到这一点是很容易的,但并没有正确体现其物理本质,因此从这种角度分析并不恰当。诸如Mayer键级甚至都能和某一系列体系的pi-pi堆积强度找到那么点正相关,但明显讨论角度是错的(众所周知pi-pi堆积本质来自于色散作用),因此这样的相关性只不过是“伴生”的现象,是把次要因素给过度放大了,令读者对作用本质产生了错误理解。既然有更适合的分析氢键的方法,明显就没必要用NBO(甚至宁可通过氢键距离、键角这样几何角度,也不应该用对当前不恰当、不适用的电子结构分析方法)。

E2并没有丝毫体现氢键的静电作用本质,它基于的是轨道作用模型,只能与氢键作用能中占很小比例的CT作用能部分有联系(一个单体的占据轨道与另一个单体的空轨道混合导致CT的发生)。
北京科音自然科学研究中心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 (一方通行)

10#
发表于 Post on 2018-9-11 17:10:47 来自手机 | 只看该作者 Only view this author
sobereva 发表于 2018-9-11 09:56
E2能解释相对强度,有相关性,达到这一点是很容易的,但并没有正确体现其物理本质,因此从这种角度分析并 ...

关于氢键本质究竟静电占多少比重,我以为也分具体的氢键体系吧。比如,史上最强氢键是FHF^-离子,这个氢键共价成分很高,是典型的3c-4e键;而H2O与CH4之间的氢键电荷转移非常少,不能被视为纯静电作用吧,至少donor-acceptor轨道作用有一定比例,最起码这种轨道作用能够完美解释这两个分子儿之间的几何取向。

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

11#
发表于 Post on 2018-9-11 20:22:31 | 只看该作者 Only view this author
granvia 发表于 2018-9-11 17:10
关于氢键本质究竟静电占多少比重,我以为也分具体的氢键体系吧。比如,史上最强氢键是FHF^-离子,这个氢 ...

肯定得分具体体系啊,没说具体前提下都是默认说的一般强度氢键
静电作用主导,一直都是说的是常规氢键。FHF-那种氢键显然不能拿静电作用模型说事
弱氢键众所周知色散作用占很大成份,不考虑定量完全没法用
不管什么情况,都轮不到E2说事
北京科音自然科学研究中心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 (一方通行)

12#
发表于 Post on 2018-9-12 00:28:12 来自手机 | 只看该作者 Only view this author
sobereva 发表于 2018-9-11 20:22
肯定得分具体体系啊,没说具体前提下都是默认说的一般强度氢键
静电作用主导,一直都是说的是常规氢键。 ...

粗略看了一下Stone的第一篇文章,个人感觉他们似乎没弄懂NBO的思路。。 待我有时间看看后面的来回

1043

帖子

0

威望

4174

eV
积分
5217

Level 6 (一方通行)

13#
发表于 Post on 2018-9-12 00:35:27 来自手机 | 只看该作者 Only view this author
sobereva 发表于 2018-9-11 20:22
肯定得分具体体系啊,没说具体前提下都是默认说的一般强度氢键
静电作用主导,一直都是说的是常规氢键。 ...

另外,Stone文章用的方法是PBE0 functional with asymptotic correction,这能很好地描述dispersion贡献?

6万

帖子

99

威望

5万

eV
积分
120179

管理员

公社社长

14#
发表于 Post on 2018-9-12 05:33:17 | 只看该作者 Only view this author
granvia 发表于 2018-9-12 00:35
另外,Stone文章用的方法是PBE0 functional with asymptotic correction,这能很好地描述dispersion贡献 ...

不加色散校正光靠PBE0是肯定没法描述色散作用的
Stone也是个颇有经验的专家,我记好像80年代还和Weinhold合作发过文章。从他写的文章上看,我认为已经认识到NBO的根上了

有些看似NBO也能解释不错的弱相互作用问题,比如weinhold用E2也讨论了卤键的角度问题,其实用物理意义更严格、确切的分析角度讨论,比如静电势,SAPT等,可以从本质上解释得更好(比如Politzer通过静电势讨论卤键的可能近百篇都有了,Stone在J. Am. Chem. Soc., 135, 7005通过扫描卤键键角从SAPT能量成分上进行更深入解释)。一个问题,可能表面上看似A角度和B角度都能说得通,但如果A能体现问题的95%实质,B只能表现本质的5%,那肯定倾向于用A角度分析。
北京科音自然科学研究中心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

50

帖子

0

威望

2550

eV
积分
2600

Level 5 (御坂)

15#
发表于 Post on 2018-9-13 05:18:50 | 只看该作者 Only view this author
granvia 发表于 2018-9-12 00:35
另外,Stone文章用的方法是PBE0 functional with asymptotic correction,这能很好地描述dispersion贡献 ...

Stone文章的PBE0 with asymptotic correction是用来做SAPT(DFT)中的monomer calculation的,分子间的dispersion是由SAPT(DFT)本身计算得到,不依赖于这个functional

本版积分规则 Credits rule

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

GMT+8, 2025-8-17 01:42 , Processed in 0.162462 second(s), 20 queries , Gzip On.

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