请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 1465|回复: 6

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

[复制链接]

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

发表于 2017-12-17 20:41:45 | 显示全部楼层 |阅读模式
绘制有填色效果的用于弱相互作用分析的RDG散点图的方法

文/Sobereva @北京科音  2017-Dec-17

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

本文使用Multiwfn 3.5(dev)版,这里用通过苯酚二聚体来示例,相应的波函数文件是Multiwfn目录下的PhenolDimer.wfn。

启动Multiwfn,依次输入以下命令,让Multiwfn把此体系的RDG vs sign(lambda2)rho的散点数据导出到当前目录下的output.txt中。
100
1
15,13
-10
0
2

得到的output.txt文件有5列,最后两列是被作为散点图的x,y数据。故用Ultraedit等程序,通过列模式把所有行的前三列都删掉。如果在Linux中,也可以将下面这行复制到命令行窗口运行,则当前目录下的output.txt就只剩最后两列数据了。
cut -c 34-66 output.txt > t.txt;mv t.txt output.txt -f

gnuplot是个免费的轻量级的基于命令行的数据作图程序,各种系统都支持,可以在这里下载:http://www.gnuplot.info。本文用的是gnuplot 5.2.2 Windows版。将处理好的output.txt放到gnuplot目录下的bin目录下,然后将RDGmap.gnu(下载: RDGmap.gnu (706 Bytes, 下载次数: 7)

评分

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

查看全部评分

北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

50

帖子

1

威望

968

eV
积分
1038

Level 4 (黑子)

发表于 2017-12-19 11:02:45 | 显示全部楼层
这个真漂亮

166

帖子

3

威望

2465

eV
积分
2691

Level 5 (御坂)

发表于 2017-12-19 11:20:25 | 显示全部楼层
我都是手工用PS通过渐变画出来的
Monte Carlo

10

帖子

0

威望

160

eV
积分
170

Level 3 能力者

发表于 2017-12-19 21:22:54 | 显示全部楼层
gnuplot网站为什么一直500错误啊

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

 楼主| 发表于 2017-12-19 21:53:56 | 显示全部楼层
w33 发表于 2017-12-19 21:22
gnuplot网站为什么一直500错误啊

你的网有问题
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

1

帖子

0

威望

19

eV
积分
20

Level 1 能力者

发表于 2017-12-24 20:39:56 | 显示全部楼层
你好,请问有填充颜色的散点图会画不?

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

 楼主| 发表于 2017-12-24 20:53:56 | 显示全部楼层
hljiang531 发表于 2017-12-24 20:39
你好,请问有填充颜色的散点图会画不?

病句,意义不明
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-9-26 13:46 , Processed in 0.264217 second(s), 28 queries .

快速回复 返回顶部 返回列表