计算化学公社

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

[波函数分析交流] 用Hirshfeld电荷和信息增益研究亲电性和亲核性

[复制链接 Copy URL]

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

用Hirshfeld电荷和信息增益研究亲电性和亲核性
Studying electrophilicity and nucleophilicity using Hirshfeld charge and information gain

Sobereva@北京科音   2014-Nov-1


对于化学体系,Shannon熵定义为-∫ρ(r)lnρ(r) dr,其中ρ是体系的电子密度。

划分分子成为原子空间的方法很多,如Hirshfeld、AIM、Becke、ISA等。Nalewajski和Parr在PNAS, 97, 8879 (2000)中指出,如果用Hirshfeld方式划分分子空间来定义原子(其密度记为ρ_A、ρ_B、ρ_C...),那么形成分子后,各个原子相对于它在孤立状态(密度记为ρ0_A、ρ0_B、ρ0_C...)的信息熵变,或者说原子的信息增益,都是最小的(如A原子的增益=-∫ρ_A(r)ln[ρ_A(r)/ρ0_A(r)] dr)。同时分子整体的信息增益(即原子信息增益之和)也是最小的,这称为最小信息增益原理。在这篇文章发表后,和分子、原子的信息熵以及Hirshfeld划分相关的文章大量出现。

在最近的一篇文章Shubin Liu, Chuying Rong, Tian Lu, Information Conservation Principle Determines Electrophilicity, Nucleophilicity, and Regioselectivity, J. Phys. Chem. A 2014, 118, 3698−3704当中,作者指出在一阶近似下,Hirshfeld方式划分导致的分子信息增益恰为0,并称之为信息保守原理。并且,在一阶近似下,原子的信息增益和原子的Hirshfeld电荷是等价的!这将抽象的信息熵概念和计算化学工作者熟知的原子电荷直接联系了起来。

此文通过50个只含碳、氢的分子作为测试集,基于Multiwfn程序计算了信息增益和Hirshfeld电荷,发现它们与Mayr实验亲电/亲核性标度有绝好的线性关系!如下图所示。并且发现Hirshfeld电荷不仅与亲电/亲核反应速率密切相关,还能正确地预测出亲核、亲电反应活性位点。

但此文只考虑了含碳体系,为了检验是否这样好的关系对于其它元素也有效,在 周夏禹, 荣春英, 卢天, 刘述斌, 用Hirshfeld电荷定量标度亲电性和亲核性:含氮体系, 物理化学学报, 30(11), 2055 (2014) (http://www.whxb.pku.edu.cn/CN/abstract/abstract28901.shtml)一文当中,作者又使用了5 个不同类别的含氮体系共计40 个分子进行了研究。结果发现对所有五个含氮体系其Hirshfeld电荷与实验亲电亲核性标度之间仍然存在着很强的线性关联。但同时发现,这些相关性依赖于氮元素的化合价类型和键合环境,线性关系只能在同一类型中成立。

那么,Hirshfeld原子电荷,或者说一阶近似下的原子信息增益有什么用?通过上面的讨论,可以知道,起码对于含碳、氮的体系,只要计算出了Hirshfeld原子电荷,利用文中给出的关系式,就能得出对应的亲电/亲核性标度,这与反应速率直接相关。并且还可以预测亲核/亲电反应优先发生在哪个位点。


本文提到的三篇文章在下面下载
Information theory, atoms in molecules, and molecular similarity.pdf (74.49 KB, 下载次数 Times of downloads: 113)
Information Conservation Principle Determines Electrophilicity, Nucleophilicity,.pdf (2.87 MB, 下载次数 Times of downloads: 181)
用Hirshfeld 电荷定量标度亲电性和亲核性-含氮体系.pdf (789.15 KB, 下载次数 Times of downloads: 339)

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
aqhuangry + 4

查看全部评分 View all ratings

北京科音自然科学研究中心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
积分
120110

管理员

公社社长

21#
 楼主 Author| 发表于 Post on 2021-6-24 01:09:25 | 只看该作者 Only view this author
945129408 发表于 2021-6-23 22:59
用静电势考察亲核试剂的亲核位点是不是也是糟糕选择,而应该用Hirshfeld电荷

往往是,但不一定。参考

Comparative study on the methods for predicting the reactive site of nucleophili.pdf (1.89 MB, 下载次数 Times of downloads: 9)
北京科音自然科学研究中心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

4

帖子

0

威望

23

eV
积分
27

Level 2 能力者

20#
发表于 Post on 2021-6-23 22:59:18 | 只看该作者 Only view this author
用静电势考察亲核试剂的亲核位点是不是也是糟糕选择,而应该用Hirshfeld电荷

4

帖子

0

威望

23

eV
积分
27

Level 2 能力者

19#
发表于 Post on 2021-6-23 22:56:21 | 只看该作者 Only view this author
本帖最后由 945129408 于 2021-6-23 22:57 编辑

静电势去考察亲核试剂的亲核位点是不是也是糟糕的选择?

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

18#
 楼主 Author| 发表于 Post on 2015-7-1 11:04:02 | 只看该作者 Only view this author
xieyi 发表于 2015-7-1 09:35
请问Sob老师,我的分子带有C H O,是不是遵循帖子所说的线性关系?亲核亲电能不能用Hirshifeld来预测?我 ...

不光看元素,也看成键类型。帖子提到,同样是氮,在不同化学环境下线性关系不同。你可以先试试。
这类问题最鼓励的就是用Hirshfeld电荷,这类问题不鼓励用NPA电荷。
北京科音自然科学研究中心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

19

帖子

0

威望

226

eV
积分
245

Level 3 能力者

17#
发表于 Post on 2015-7-1 09:35:32 | 只看该作者 Only view this author
sobereva 发表于 2015-5-26 22:32
文中只是全部计算时用scrf模型表现这种溶剂。
这种低极性溶剂对反应位点判断是不会有什么影响的,定量数 ...

请问Sob老师,我的分子带有C H O,是不是遵循帖子所说的线性关系?亲核亲电能不能用Hirshifeld来预测?我目前用的NPA

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

16#
 楼主 Author| 发表于 Post on 2015-5-28 14:04:05 | 只看该作者 Only view this author
此文专门讨论这些问题
Phys.Chem.Chem.Phys.,2015, 17, 4977
你若要用shape function来定义稍微改下代码就行。

此问题我并不怎么关心
北京科音自然科学研究中心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

52

帖子

0

威望

1660

eV
积分
1712

Level 5 (御坂)

15#
发表于 Post on 2015-5-28 13:29:41 | 只看该作者 Only view this author
基于shape function得到的才是正在的Shannon entropy,因为都是归一化的,而基于电子密度的应该给另外一个名字,因为电子密度的积分是等于电子总数,所以文献中有些以讹传讹的概念。

52

帖子

0

威望

1660

eV
积分
1712

Level 5 (御坂)

14#
发表于 Post on 2015-5-28 13:22:45 | 只看该作者 Only view this author
本帖最后由 ORCA_in_TCC 于 2015-5-28 13:27 编辑

这就是“N-derivative of Shannon entropy of shape function for atoms” Chem. Physi. Lett. 410 (2005) 70–76一文。从这篇文章中的表2给出的数据可以发现,Ne之后的基于电子密度的Shannon entropy都是负值,而表1中的基于Shape function的Shannon entropy都是正值。


QQ截图20150528131829.jpg (70.71 KB, 下载次数 Times of downloads: 179)

QQ截图20150528131829.jpg

52

帖子

0

威望

1660

eV
积分
1712

Level 5 (御坂)

13#
发表于 Post on 2015-5-28 13:12:16 | 只看该作者 Only view this author
但是这样计算出来的原子Shannon entropy对于原子序数大于10的原子会导致负熵的出现,这应该怎么去理解?

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

12#
 楼主 Author| 发表于 Post on 2015-5-27 23:06:32 | 只看该作者 Only view this author
ORCA_in_TCC 发表于 2015-5-27 15:58
sob老师,multiwfn里面Shannon entropy是基于电子密度还是shape function(电子密度除以电子总数)?

电子密度。

手册2.6节有公式(编号50)
北京科音自然科学研究中心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

52

帖子

0

威望

1660

eV
积分
1712

Level 5 (御坂)

11#
发表于 Post on 2015-5-27 15:58:16 | 只看该作者 Only view this author
sob老师,multiwfn里面Shannon entropy是基于电子密度还是shape function(电子密度除以电子总数)?

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

10#
 楼主 Author| 发表于 Post on 2015-5-26 22:32:32 | 只看该作者 Only view this author
tjuchan 发表于 2015-5-26 16:08
看来《Information Conservation Principle Determines Electrophilicity,Nucleophilicity, and Regiosel ...

文中只是全部计算时用scrf模型表现这种溶剂。
这种低极性溶剂对反应位点判断是不会有什么影响的,定量数值也不会有多少改变,没有可讨论的。
北京科音自然科学研究中心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

105

帖子

0

威望

867

eV
积分
974

Level 4 (黑子)

9#
发表于 Post on 2015-5-26 16:08:17 | 只看该作者 Only view this author
sobereva 发表于 2015-5-18 21:53
怎么考虑得看情况。
隐式溶剂模型只能表现一个平均场。如果加了碱会导致反应发生了显著变化,那就应该 ...

看来《Information Conservation Principle Determines Electrophilicity,Nucleophilicity, and Regioselectivity》这篇文章,文中提及的采用CH2CL2作为溶剂,但是好像没有看见相关内容呢?

39

帖子

0

威望

166

eV
积分
205

Level 3 能力者

8#
发表于 Post on 2015-5-19 09:15:31 | 只看该作者 Only view this author
sobereva 发表于 2015-5-18 21:53
怎么考虑得看情况。
隐式溶剂模型只能表现一个平均场。如果加了碱会导致反应发生了显著变化,那就应该 ...

谢谢

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

GMT+8, 2025-8-14 18:48 , Processed in 0.990398 second(s), 31 queries , Gzip On.

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