计算化学公社

 找回密码 Forget password
 注册 Register
Views: 11055|回复 Reply: 4

[Multiwfn资源与经验] 利用约化密度梯度(RDG)考察AIM临界点的位置

[复制链接 Copy URL]

4万

帖子

99

威望

4万

eV
积分
89863

管理员

公社社长+计算化学玩家

发表于 Post on 2014-12-14 07:17:42 | 显示全部楼层 Show all |阅读模式 Reading model
利用约化密度梯度(RDG)考察AIM临界点的位置

文/Sobereva @北京科音
First release: 2014-Dec-14  Last update: 2022-Aug-17

约化密度梯度(RDG)定义如下
1.png

RDG把电子密度梯度转化为了无量纲的形式,用于GGA泛函中。RDG的其它用处想必很多人已经知道了,就是可以考察弱相互作用区域,这在《使用Multiwfn图形化研究弱相互作用》(http://sobereva.com/68)、《使用Multiwfn研究分子动力学中的弱相互作用》(http://sobereva.com/186)中都已经介绍过了。RDG还有一个不错的用处是可以帮助确定AIM临界点位置,或者解释为什么临界点存在或不存在。

AIM临界点就是指电子密度梯度(具体来说是它的模)为0的位置,介绍见《使用Multiwfn做拓扑分析以及计算孤对电子角度》(http://sobereva.com/108)。因此,如果图形化展现电子密度梯度模,就可以直观地看出临界点位置了,而用不着必须做临界点搜索。然而,电子密度梯度的数量级跨度太大,不便于分析,而RDG就没这个问题,数值范围较小,绘图比较方便。RDG=0的位置和电子密度梯度=0的位置是相同的。

Multiwfn (http://sobereva.com/multiwfn)可以极为方便地绘制RDG,本文使用3.3.6版。默认情况下,电子密度>0.05的区域的RDG会被自动设为100,原因在《使用Multiwfn图形化研究弱相互作用》中有介绍。为了让所有区域的RDG都能显示,我们先把settings.ini里的RDG_maxrho改为0。

我们这里用C2H5F作为例子,启动Multiwfn,输入
examples\C2H5F.wfn
5
13
3
2   //如果想直接用Multiwfn自带的图形界面观看等值面这里就选-1,不过速度不是很快
在当前目录下得到了RDG.cub,是RDG的cube文件。将之载入到VMD观看等值面,将等值面数值设为一个较小的值,比如0.07,得到下图

2.jpg


每个等值面里面都是数值很小,<0.07的区域。其中就有RDG=0的点,即AIM临界点位置。因此,可以直接观看数值很小的RDG等值面来考察哪里有AIM临界点。但注意,虽然原子核也是AIM临界点,但在上图中没有出现等值面,上图中等值面只展现了键临界点(BCP)的位置,这是因为原子核附近的RDG数值远大于BCP附近的,当前RDG数值设得太小,与原子核对应的等值面太小而看不到。


通过等值面图考察不算很精确。如果只需要考察一个平面,绘制RDG填色图是更好的方法。比如这里绘制尿嘧啶的图。启动Multiwfn依次输入
examples\uracil.wfn
4
13
1
[ENTER]
0
2
1
0

3.jpg

从图中可见,在每个原子之间,以及环中心都有颜色非常深的一小块区域表现出相应位置RDG非常小,其中颜色最深的那个点对应RDG=0,正是AIM临界点位置。从图中也可以看出有的原子核的位置上也有个深蓝的小点,对应于核临界点位置(没显示出深蓝小点的原子纯粹是因为那个点太细微了,实际上是都有的)。


利用RDG图形,对于解释AIM临界点的存在与否问题很有用。比如前一阵子有位Multiwfn用户问我,为什么下图中Os1和C9之间的键临界点以及1-5-9三元环中的环临界点没有被Multiwfn搜索出来,表面看上去按理说应该是有临界点才对
4.jpg


他怀疑是临界点搜索方法不对或者搜索参数没用对。实际上,使用RDG图形分析,就可以方便地搞明白这个问题。先对1-5-9三元环绘制RDG填色图,色彩刻度是0~1。

5.jpg


图中可见,1-5-9三元环之间确实没有RDG很小的区域,所以不会有环临界点。但乍一看,奇怪的是,Os1和C9之间明明存在RDG很小区域,为何拓扑分析没找到键临界点?于是猜测到,Os1-C9之间虽然有RDG很小的区域,但其中却没有RDG小到为0的点。为了看清楚,把色彩刻度改为0~0.1再绘制一次,相当于把0~0.1数值范围内的特征进行放大。结果如下所示
6.jpg


这回Os1-C9之间并没有深蓝色,即RDG为0的点,因此可证实实际上并不存在对应的键临界点。


可见,当碰到搜索不出本该有的AIM临界点,或者想搞明白某些区域内到底有没有、有哪些临界点的时候,都可按照本文的方法对RDG作图来分析。
特别值得说明的是,笔者提出的IRI和IGMH分析可以对将各种相互作用都直观展现出来,包括没有AIM临界点因此无法通过AIM理论考察的相互作用,见《使用IRI方法图形化考察化学体系中的化学键和弱相互作用》(http://sobereva.com/598)和《使用Multiwfn做IGMH分析非常清晰直观地展现化学体系中的相互作用》(http://sobereva.com/621)。为什么IRI能有这种能力,见IRI原文里深入详细的分析:Chemistry—Methods, 1, 231 (2021) DOI: 10.1002/cmtd.202100007。

评分 Rate

参与人数
Participants 5
eV +34 收起 理由
Reason
QuantumHan + 5 你太可爱
aqhuangry + 5
ter20 + 10 好物!
captain + 10 好物!
bignever + 4 赞!

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

107

帖子

0

威望

864

eV
积分
971

Level 4 (黑子)

发表于 Post on 2014-12-14 12:06:46 | 显示全部楼层 Show all
学习!

258

帖子

0

威望

1932

eV
积分
2190

Level 5 (御坂)

发表于 Post on 2015-2-22 14:34:59 | 显示全部楼层 Show all
如果用Lanl2Tz(f)对Cu或Ag或Au,6-311g**对C,N,S,H,PBE0泛函计算分子式为C21H15M3N6S3的配合物(M=Cu/Ag/Au)的RDG,是否可以接受?主要想研究金属与硫原子间有没有弱作用。

4万

帖子

99

威望

4万

eV
积分
89863

管理员

公社社长+计算化学玩家

 楼主 Author| 发表于 Post on 2015-2-22 17:32:46 | 显示全部楼层 Show all
虎王 发表于 2015-2-22 14:34
如果用Lanl2Tz(f)对Cu或Ag或Au,6-311g**对C,N,S,H,PBE0泛函计算分子式为C21H15M3N6S3的配合物(M=Cu/Ag/Au ...

足够了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

2

帖子

0

威望

4

eV
积分
6

Level 1 能力者

发表于 Post on 2015-5-27 04:25:30 | 显示全部楼层 Show all
我终于可以发帖了吗?

本版积分规则 Credits rule

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

GMT+8, 2023-2-1 23:04 , Processed in 0.432762 second(s), 31 queries .

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