计算化学公社

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

[VMD] 在VMD中使用GaussView的元素着色的方法

[复制链接 Copy URL]

5万

帖子

99

威望

5万

eV
积分
112544

管理员

公社社长

在VMD中使用GaussView的元素着色的方法
Using GaussView element coloring scheme in VMD

文/Sobereva@北京科音   2022-Sep-16


GaussView的元素配色总体上来说不错,而且对所有元素都定义了唯一的颜色,在File - Preference - Colors - Element Colors中可以看到所有元素的颜色,如下所示



VMD是免费、灵活、强大而且超级流行的化学体系可视化工具。在VMD里也可以按照元素着色,见《在VMD程序里对不同元素的原子用不同颜色显示的方法》(http://sobereva.com/624)。然而,至少对于笔者撰文时的最新的正式版1.9.3版来说,VMD自带的元素着色定义很少,对绝大多数元素都是统一用的褐色,导致很多情况下没法区分不同元素,而且颜色也不美观。这里介绍在VMD里使用GaussView里元素颜色定义的做法。

http://sobereva.com/attach/652/gview_color.tcl下载gview_color.tcl文件,将之放到VMD目录下(即VMD启动后在文本窗口里输入pwd命令后显示的目录),然后在VMD的vmd.rc文件末尾加入一行proc gview {} {source gview_color.tcl}。如果你不了解vmd.rc文件的话,看《VMD初始化文件(vmd.rc)我的推荐设置》(http://sobereva.com/545)。

启动VMD后,在文本窗口里输入gview,就可以把默认的元素着色方案替换成和GaussView相同的了。之后载入个含有元素信息的结构文件(比如pdb、xyz等),然后在Graphics - Representation里把Coloring Method设为Element即可看到按元素着色的效果。如果想和GaussView显示的图像特征尽可能接近,将Drawing Method设为CPK,Bond Radius设为0.2,Material设为AOShiny。对顺铂体系显示效果如下,左边是VMD的,右边是GaussView的



上图还可以看到在键的着色方面有差异,VMD的键的两边的颜色分别对应两个原子的,而GaussView都是白色的。如果想让VMD在这点和GaussView统一,可以把当前的Representation里的Bond Radius设为0使得键不显示,然后再建立一个Representation,用CPK,把Sphere Scale设0,Bond Radius设0.15,Coloring Method选ColorID并指定为白色,之后效果如下所示,可见和GaussView显示的很接近了



再把原子半径、光源方向、材质微调一下,就更像GaussView的效果了。不过,GaussView里对多重键的显示效果是VMD里怎么设都模仿不来的。

评分 Rate

参与人数
Participants 12
eV +56 收起 理由
Reason
cynthia12 + 3 谢谢分享
王采岑 + 5 好物!
ReviewReview + 5 精品内容
Novice + 5 赞!
huanghj + 4 好物!
AmonB + 4 赞!
wangyueda + 5 好物!
alonewolfyang + 5 精品内容
zsu007 + 5 赞!
niobium + 5 好物!
kurohimx + 5 GJ!
丁越 + 5 赞!

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口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!

5万

帖子

99

威望

5万

eV
积分
112544

管理员

公社社长

11#
 楼主 Author| 发表于 Post on 2024-6-30 05:26:00 | 只看该作者 Only view this author
不是哪块小饼干 发表于 2024-6-29 18:06
请问使用上述方法进行着色后的图片如何导出?按照之前的方法导出的是未进行修改的样子

平时图像怎么获得还怎么获得,做法跟此文没有直接关系
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口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!

3

帖子

0

威望

59

eV
积分
62

Level 2 能力者

10#
发表于 Post on 2024-6-29 18:06:28 | 只看该作者 Only view this author
请问使用上述方法进行着色后的图片如何导出?按照之前的方法导出的是未进行修改的样子

5万

帖子

99

威望

5万

eV
积分
112544

管理员

公社社长

9#
 楼主 Author| 发表于 Post on 2023-6-6 09:21:44 | 只看该作者 Only view this author
ReviewReview 发表于 2023-6-5 14:37
是的,我刚才试了下也发现了,Bond Radius没法输入0.15只能设为0.1或者0.2,设为0.1时依然是VMD的键的两 ...

用命令行可以输入0.15并生效,只不过界面上只能显示两位小数,会四舍五入
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口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!

113

帖子

0

威望

258

eV
积分
371

Level 3 能力者

8#
发表于 Post on 2023-6-5 14:37:56 | 只看该作者 Only view this author
neocc 发表于 2022-10-1 21:04
按照社长说的“在Graphics - Representation里把Coloring Method设为Element即可看到按元素着色的效果“
g ...

是的,我刚才试了下也发现了,Bond Radius没法输入0.15只能设为0.1或者0.2,设为0.1时依然是VMD的键的两边的颜色分别对应两个原子的;设为0.2时VMD的键就变成白色啦。这个试验结果似乎与0.15四舍五入约等于0.2这一数学公理相一致。
轻舟一过万重山。

113

帖子

0

威望

258

eV
积分
371

Level 3 能力者

7#
发表于 Post on 2023-6-5 11:28:03 | 只看该作者 Only view this author
ReviewReview 发表于 2023-6-4 19:26
Lu老师,我按照您上边的方法"从http://sobereva.com/attach/652/gview_color.tcl ...

已经解决啦老师。是按照pal这位网友在这个帖子"http://bbs.keinsci.com/forum.php?mod=viewthread&tid=37744&page=1#pid250770"里让我调整了下操作顺序弄出来的。万分感谢!
轻舟一过万重山。

113

帖子

0

威望

258

eV
积分
371

Level 3 能力者

6#
发表于 Post on 2023-6-4 19:26:20 | 只看该作者 Only view this author
本帖最后由 ReviewReview 于 2023-6-4 19:58 编辑

Lu老师,我按照您上边的方法"http://sobereva.com/attach/652/gview_color.tcl下载gview_color.tcl文件,将之放到VMD目录下(即VMD启动后在文本窗口里输入pwd命令后显示的目录),然后在VMD的vmd.rc文件末尾加入一行proc gview {} {source gview_color.tcl}。如果你不了解vmd.rc文件的话,看《VMD初始化文件(vmd.rc)我的推荐设置》(http://sobereva.com/545)。

启动VMD后,在文本窗口里输入gview,就可以把默认的元素着色方案替换成和GaussView相同的了。之后载入个含有元素信息的结构文件(比如pdb、xyz等),然后在Graphics - Representation里把Coloring Method设为Element即可看到按元素着色的效果。如果想和GaussView显示的图像特征尽可能接近,将Drawing Method设为CPK,Bond Radius设为0.2,Material设为AOShiny。"改了一下我用Multiwfn+VMD得到的PVDF的ESP图里的原子颜色,最后得到的更改效果是所有原子颜色都变成蓝色的了。。。这是啥原因呢?
轻舟一过万重山。

145

帖子

0

威望

3119

eV
积分
3264

Level 5 (御坂)

5#
发表于 Post on 2022-10-2 13:48:50 | 只看该作者 Only view this author
sobereva 发表于 2022-10-2 04:58


改光源位置:Mouse - Move Light,选择光源序号,然后在图形窗口里按鼠标左键拖动

谢谢社长

5万

帖子

99

威望

5万

eV
积分
112544

管理员

公社社长

4#
 楼主 Author| 发表于 Post on 2022-10-2 04:58:07 | 只看该作者 Only view this author
neocc 发表于 2022-10-1 21:04
按照社长说的“在Graphics - Representation里把Coloring Method设为Element即可看到按元素着色的效果“
g ...
请问光源方向和 Bond Radius设0.15 调整是通过命令行窗口输入代码吗?



改光源位置:Mouse - Move Light,选择光源序号,然后在图形窗口里按鼠标左键拖动
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口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!

145

帖子

0

威望

3119

eV
积分
3264

Level 5 (御坂)

3#
发表于 Post on 2022-10-1 21:04:35 | 只看该作者 Only view this author
本帖最后由 neocc 于 2022-10-1 21:07 编辑

按照社长说的“在Graphics - Representation里把Coloring Method设为Element即可看到按元素着色的效果“
gview在iso前后执行没有影响。
抱歉社长,我看漏了这一行。想当然的用着NAME           /@_@\

不过ESP等值面选透明类型的材质才行,但是 EddyGlass 对 键的着色 效果不是二者的叠加,要么是双色要么整个分子都是白色。

分子的材质选BrushedMetal并关掉Light2,3,看起来光源和gview接近,不过颜色有点淡。
个人推荐Diffuse关掉Light3 ;以及AOEddy,AOChalky和RTChrome,关掉Light1,3,效果不错。

另外,请问光源方向和 Bond Radius设0.15 调整是通过命令行窗口输入代码吗?GUI里面光源方向没找到如何设置,Bond Radius只能设为0.1或者0.2

145

帖子

0

威望

3119

eV
积分
3264

Level 5 (御坂)

2#
发表于 Post on 2022-9-30 21:22:40 | 只看该作者 Only view this author
本帖最后由 neocc 于 2022-9-30 21:32 编辑

请问社长
在VMD中绘制ESP等值面的时候,按照《在VMD中使用GaussView的元素着色的方法》修改vmd.rc的办法,在绘制完ESP等值面后,输入gview,结果原子配色没有变化。
如果按照《在VMD程序里对不同元素的原子用不同颜色显示的方法》修改每一个NAME的RGB百分值,依次输入vmd命令行窗口可以修改原子颜色。
请问绘制ESP等值面后如何用gview整体修改为gaussview的元素着色?

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

GMT+8, 2024-11-27 19:18 , Processed in 0.249040 second(s), 26 queries , Gzip On.

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