计算化学公社

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

[Multiwfn资源与经验] 绘制有填色效果的用于弱相互作用分析的RDG散点图的方法

  [复制链接 Copy URL]

6万

帖子

99

威望

5万

eV
积分
120086

管理员

公社社长

绘制有填色效果的用于弱相互作用分析的RDG散点图的方法
Method to draw RDG scatter map with coloring effect for weak interaction analysis

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


2010年提出的用于分析弱相互作用的RDG方法(文献中也普遍叫NCI方法)已被广泛用于考察各种分子间/分子内弱相互作用了,笔者也写过不少相关文章,不了解此方法的者看《使用Multiwfn图形化研究弱相互作用》(http://sobereva.com/68)以及里面提及的相关文章,特别是推荐看《一篇最全面介绍各种弱相互作用可视化分析方法的文章已发表!》(http://sobereva.com/667)和《Angew. Chem.上发表了全面介绍各种共价和非共价相互作用可视化分析方法的综述》(http://sobereva.com/746)里介绍的笔者的综述。能做RDG分析的程序不少,Multiwfn(http://sobereva.com/multiwfn)是其中最为流行、强大、好用的。在上面的博文里已经介绍了怎么在Multiwfn中直接绘制RDG vs sign(lambda2)rho的散点图来考察弱相互作用。一些文章里的这种散点图还加上了填色效果,可以使得对应不同横坐标的spike的颜色一目了然,便于与Multiwfn+VMD绘制的RDG填色等值面图相对应来讨论问题。其实这种图稍有photoshop使用技能的人都可以不太困难地作出来,就是把VMD的色彩刻度条在ps里拉伸成与散点图作图范围相同的大小,垫在Multiwfn给出的散点图下方的图层,然后再把散点图的图层当中作图区域的黑色部分以色彩范围选择方式选中,删除,透出来下层的色彩刻度层即可。不过肯定有不少人嫌这种做法麻烦,此文介绍一种利用gnuplot程序的简单快捷的方法绘制这种填色RDG散点图。

读者请使用2019-Aug-24及之后更新的Multiwfn。这里用通过苯酚二聚体来示例,相应的波函数文件是Multiwfn文件包里examples目录下的PhenolDimer.wfn。本文的操作在《使用Multiwfn做NCI分析展现分子内和分子间弱相互作用》(https://www.bilibili.com/video/av71561024)里也有视频演示。

启动Multiwfn,依次输入以下命令,让Multiwfn把此体系的RDG vs sign(lambda2)rho的散点数据导出到当前目录下的output.txt中。
examples\PhenolDimer.wfn
20  //弱相互作用图形化分析
1  //NCI分析
3  //高质量格点
2  //导出散点数据

gnuplot是个免费的轻量级的基于命令行的数据作图程序,各种系统都支持,可以在这里下载:http://www.gnuplot.info。本文用的是gnuplot 5.4 Windows版。将output.txt放到gnuplot目录下的bin子目录下,然后将Multiwfn目录下的examples\scripts\RDGscatter.gnu这个绘图脚本也拷到此目录下。开启操作系统的命令行模式(例如Windows下的cmd环境)并进入此目录,运行命令gnuplot RDGscatter.gnu(对于Windows用户,这一步不知道怎么弄的话直接把RDGscatter.gnu拖到gnuplot.exe图标上也行),就会在当前目录下产生RDGscatter.ps,这就是填色散点图的postscript格式的文件了。这是一种矢量图形格式,可无损缩放,很多程序都可以查看。比如可以直接用acrobat打开,打开后可以无损缩放。也可以用photoshop打开,打开的时候可以选择产生像素为多大的图片。如果机子里装了ghostscript程序,也可以用小巧且强大的看图程序irfanview观看。如果你懒得装单机程序,也可以用免费的在线程序https://cloudconvert.com/image-converter把ps格式转成常见图像格式。此例效果如下:



RDGscatter.gnu脚本里有很多参数可以设定,比较关键的参数就是X,Y轴以及色彩刻度轴的上下限(x/y/cbrange后面的值)、标签的数值范围和步长(x/y/cbtic后面的值)、散点的大小(pointsize后面的值),以及色彩刻度的定义。笔者习惯在VMD中用-0.035~0.02来对RDG等值面着色,色彩刻度是默认的蓝-绿-红,因此脚本中可以看到这样的设定
set palette defined (-0.035 "blue",-0.0075 "green", 0.02 "red")

如果要把填色的散点图与VMD绘制的填色的RDG等值面图相对照,则二者色彩刻度设定必须严格一致。比如在Multiwfn目录下的examples\RDGfill.vmd文件就是VMD里绘制填色等值面图的脚本,这里面mol scaleminmax top 1那一行后面的值应该设为-0.035 0.02才能与上图来对照(默认就是如此)。在这种色彩刻度下绘制的苯酚二聚体的RDG填色图如下所示,很明显散点图上各个spike位置和RDG填色图上的等值面通过颜色很容易进行一一对应。



笔者在《使用IRI方法图形化考察化学体系中的化学键和弱相互作用》(http://sobereva.com/598)中介绍的我提出的IRI方法比RDG方法明显更强大,不仅可以展现弱相互作用区域,还可以展现化学键作用区域,因此强烈建议用IRI取代RDG。Multiwfn里IRI分析的界面和RDG分析如出一辙,可以导出IRI vs sign(lambda2)rho的散点数据到output.txt。之后用examples\scripts\IRIscatter.gnu(2022-Jul-16及以后更新的Multiwfn才有)代替上文的RDGscatter.gnu就可以绘制出填色的散点图。如果把脚本里的横坐标范围设大,比如设到-0.05到0.03范围,还可以使散点图把化学键作用区域的spike展现出来。例子看Multiwfn手册4.20.4节。

笔者提出的IGMH方法现在也特别流行,远比RDG更适合专门考察片段间的弱相互作用,在《使用Multiwfn做IGMH分析非常清晰直观地展现化学体系中的相互作用》(http://sobereva.com/621)博文中有详细介绍。IGMH方法也可以绘制填色散点图,这在此博文里提到的IGMH官方教程里有很具体的例子。

评分 Rate

参与人数
Participants 15
eV +71 收起 理由
Reason
wch990122 + 4 谢谢
ghifi37 + 5 好物!
yoya + 5 精品内容
LJF0220 + 4 谢谢~~
不懂计算 + 5 精品内容
w33 + 4
teller3531 + 5
aqhuangry + 5
yjr + 5 精品内容
jctyp + 4 谢谢分享
longzhixiu + 5 精品内容
978142355 + 5 GJ!
三石草祭 + 5 谢谢
我本是个娃娃 + 5 你太可爱
winterzen + 5 好物!

查看全部评分 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

82

帖子

0

威望

607

eV
积分
689

Level 4 (黑子)

38#
发表于 Post on 2024-7-23 15:37:50 | 只看该作者 Only view this author
sobereva 发表于 2024-7-23 03:36
1 拿不准就一律带着色散校正,下文明确说了
谈谈“计算时是否需要加DFT-D3色散校正?”
http://sobere ...

好的,感谢卢老师指导

6万

帖子

99

威望

5万

eV
积分
120086

管理员

公社社长

37#
 楼主 Author| 发表于 Post on 2024-7-23 03:36:16 | 只看该作者 Only view this author
jingetiema6112 发表于 2024-7-19 20:46
卢老师您好,请教您2个问题,(1)用DFT-B3LYP/6-311G**优化邻硝基苯胺的结构,并且基于这个优化结构算了一 ...

1 拿不准就一律带着色散校正,下文明确说了
谈谈“计算时是否需要加DFT-D3色散校正?”
http://sobereva.com/413

2 严格来说是
北京科音自然科学研究中心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

82

帖子

0

威望

607

eV
积分
689

Level 4 (黑子)

36#
发表于 Post on 2024-7-19 20:46:52 | 只看该作者 Only view this author
卢老师您好,请教您2个问题,(1)用DFT-B3LYP/6-311G**优化邻硝基苯胺的结构,并且基于这个优化结构算了一系列性质。但是当研究邻硝基苯胺中分子内氢键这种弱相互作用时,看您的博文说应该考虑色散。像这种情况还能不能用DFT-B3LYP/6-311G**优化的这个结构啊?还是必须考虑色散,比如在B3LYP-D3(BJ)/6-311G**水平上重新优化,再基于这个重新优化的结构研究弱相互作用?
(2)仍然是用DFT-B3LYP/6-311G**优化邻硝基苯胺的结构,但是预测某些性质时,用到的公式是基于DFT-B3PW91/6-31G*得到的经验公式。像这种情况是不是也需要把邻硝基苯胺在DFT-B3PW91/6-31G*下重新优化,然后再用这个公式进行计算?

6万

帖子

99

威望

5万

eV
积分
120086

管理员

公社社长

35#
 楼主 Author| 发表于 Post on 2022-7-23 14:33:51 | 只看该作者 Only view this author
好运来 发表于 2022-7-23 11:34
老师您好,我完全跟着您的方法用gnuplot生成的.ps文件,ghostscipt和irfanview装好之后里面自带的例子可 ...

没具体信息没法说
至少得提供你的output.txt,能让别人重复你的情况
北京科音自然科学研究中心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

2

帖子

0

威望

17

eV
积分
19

Level 1 能力者

34#
发表于 Post on 2022-7-23 11:34:14 | 只看该作者 Only view this author
sobereva 发表于 2022-6-30 00:48
ghostscript没恰当安装,或者没被irfanview识别,或者你的ps文件本身就有毛病

老师您好,我完全跟着您的方法用gnuplot生成的.ps文件,ghostscipt和irfanview装好之后里面自带的例子可以打开,但我自己生成的ps文件就打不开,我也觉得是文件的问题,该怎么解决呢

2

帖子

0

威望

17

eV
积分
19

Level 1 能力者

33#
发表于 Post on 2022-7-23 09:55:45 | 只看该作者 Only view this author
Manyu 发表于 2022-6-29 11:38
请问为何我下载了ghostscript,还是无法用irfanview打开RDGmap.ps

您好,我也出现了这个问题,请问您解决了吗

6万

帖子

99

威望

5万

eV
积分
120086

管理员

公社社长

32#
 楼主 Author| 发表于 Post on 2022-7-18 11:44:19 | 只看该作者 Only view this author
20mllu1 发表于 2022-7-18 10:02
老师您好,这是我算IRI散点图的输出文件(图片1)

根据老师您放B站的教程视频,我也是用UE删了前面三列 ...

目前版本Multiwfn根本就不需要自己改txt文件,仔细看1L

最新版Multiwfn的examples\scripts目录里直接给了IRIscatter.gnu,直接用就完了
北京科音自然科学研究中心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

11

帖子

0

威望

45

eV
积分
56

Level 2 能力者

31#
发表于 Post on 2022-7-18 10:02:36 | 只看该作者 Only view this author
老师您好,这是我算IRI散点图的输出文件(图片1)

根据老师您放B站的教程视频,我也是用UE删了前面三列(图片2)

但是画出的图是这样子的(图片3),横坐标是从-0.05到0.05,而Multiwfn直接绘图是这样的(图4),横坐标从-0.4到0.1.

然后我用记事本改了RDGscatter.gnu中x轴范围图(5),画出的散点图依旧没有任何变化。

请问这种情况该如何解决呢?

202207180953507288..png (71.98 KB, 下载次数 Times of downloads: 126)

202207180953507288..png

202207180955436862..png (156.16 KB, 下载次数 Times of downloads: 132)

202207180955436862..png

202207180956082155..png (175.47 KB, 下载次数 Times of downloads: 133)

202207180956082155..png

202207180958478416..png (49.33 KB, 下载次数 Times of downloads: 128)

202207180958478416..png

202207181001263022..png (121.6 KB, 下载次数 Times of downloads: 122)

202207181001263022..png

6万

帖子

99

威望

5万

eV
积分
120086

管理员

公社社长

30#
 楼主 Author| 发表于 Post on 2022-6-30 00:48:43 | 只看该作者 Only view this author
Manyu 发表于 2022-6-29 11:38
请问为何我下载了ghostscript,还是无法用irfanview打开RDGmap.ps

ghostscript没恰当安装,或者没被irfanview识别,或者你的ps文件本身就有毛病
北京科音自然科学研究中心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

1

帖子

0

威望

11

eV
积分
12

Level 1 能力者

29#
发表于 Post on 2022-6-29 11:38:35 | 只看该作者 Only view this author
请问为何我下载了ghostscript,还是无法用irfanview打开RDGmap.ps

3

帖子

0

威望

65

eV
积分
68

Level 2 能力者

28#
发表于 Post on 2021-5-7 13:49:34 | 只看该作者 Only view this author
sobereva 发表于 2021-5-7 02:58
photoshop的事
用acrobat或者irfanview+ghostscript都能正常看到字符

问题已解决,谢谢老师!

6万

帖子

99

威望

5万

eV
积分
120086

管理员

公社社长

27#
 楼主 Author| 发表于 Post on 2021-5-7 02:58:19 | 只看该作者 Only view this author
不可理喻 发表于 2021-5-6 08:47
电脑系统是win10的,Multiwfn用的是3.7版本,gnuplot用的是5.4.1版本,图是我用PS2020打开ps文件后导出的 ...

photoshop的事
用acrobat或者irfanview+ghostscript都能正常看到字符
北京科音自然科学研究中心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

3

帖子

0

威望

65

eV
积分
68

Level 2 能力者

26#
发表于 Post on 2021-5-6 08:47:05 | 只看该作者 Only view this author
本帖最后由 不可理喻 于 2021-5-6 08:48 编辑
sobereva 发表于 2021-5-6 02:10
截图,并且说清楚用的是Multiwfn哪个版本,什么系统的

电脑系统是win10的,Multiwfn用的是3.7版本,gnuplot用的是5.4.1版本,图是我用PS2020打开ps文件后导出的图片 D:\gnuplot\gnuplot\bin\RDGscatter.png

6万

帖子

99

威望

5万

eV
积分
120086

管理员

公社社长

25#
 楼主 Author| 发表于 Post on 2021-5-6 02:10:22 | 只看该作者 Only view this author
不可理喻 发表于 2021-5-5 21:33
老师你好,我想请问一下为什么我做出的散点图的横坐标是sign(l_2)r (a.u.),而不是跟视频里一样是sign(λ_2 ...

截图,并且说清楚用的是Multiwfn哪个版本,什么系统的
北京科音自然科学研究中心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

3

帖子

0

威望

65

eV
积分
68

Level 2 能力者

24#
发表于 Post on 2021-5-5 21:33:47 | 只看该作者 Only view this author
本帖最后由 不可理喻 于 2021-5-5 21:38 编辑

老师你好,我想请问一下为什么我做出的散点图的横坐标是sign(l_2)r (a.u.),而不是跟视频里一样是sign(λ_2)ρ (a.u.)

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

GMT+8, 2025-8-13 13:03 , Processed in 0.206670 second(s), 31 queries , Gzip On.

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