计算化学公社

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

[建模与可视化] gromacs结果的快照根据分子密度进行染色,得到heatmap

[复制链接 Copy URL]

12

帖子

0

威望

131

eV
积分
143

Level 2 能力者

如题,我想对gromacs运算后得到的某一帧结构进行后处理,我研究的是在NVT系综下,水分子在载体表面的蒸发行为,能够观察到随着模拟时间的增长,载体表面的水分子会有向盒子上方运动的趋势,水分子的密度出现稀疏之分,我在文献中看到,作者根据水分子密度进行颜色映射,进行了密度的可视化,如图所示,其中水分子低的位置用红色染色,水分子密度高的位置用蓝色染色,密度从低到高用颜色渐变体现。我尝试利用VMD的volmap功能来统计并可视化水分子密度,得到结果如下图所示,我发现整体只能体现一种颜色,我通过调整无法实现颜色体现密度高低的模式,想请教一下大家该怎么处理。
我想尝试参考sob大佬提供的ESP染色模式,但我不是很清楚我该去获得什么形式的文件来体现水分子密度的高低,例如ESP染色中根据静电势高低和正负进行染色,ESP可以使用Multiwfn获得,那我需要的这种反应密度的"ESP"该如何获得呢,恳请各位大佬指点一二!

202506022029023741..png (816.12 KB, 下载次数 Times of downloads: 20)

volmap绘制结果

volmap绘制结果

202506022021275301..png (588.53 KB, 下载次数 Times of downloads: 22)

文献中的水分子密度可视化图

文献中的水分子密度可视化图

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

2#
发表于 Post on 2025-6-2 20:47:50 | 只看该作者 Only view this author
用诸如quicksurf方式显示等值面,用空间分布函数的格点数据进行着色
北京科音自然科学研究中心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

12

帖子

0

威望

131

eV
积分
143

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2025-6-2 22:29:51 | 只看该作者 Only view this author
sobereva 发表于 2025-6-2 20:47
用诸如quicksurf方式显示等值面,用空间分布函数的格点数据进行着色

谢谢sob老师,我想问问volmap给的这个数据我能用于着色嘛,如果不行,我去哪里获取空间分布函数的格点数据呢

12

帖子

0

威望

131

eV
积分
143

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2025-6-2 22:30:17 | 只看该作者 Only view this author
这是我的volmap的参数

202506022230069523..png (82.09 KB, 下载次数 Times of downloads: 21)

202506022230069523..png

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

5#
发表于 Post on 2025-6-2 23:30:10 | 只看该作者 Only view this author
杂算者 发表于 2025-6-2 22:29
谢谢sob老师,我想问问volmap给的这个数据我能用于着色嘛,如果不行,我去哪里获取空间分布函数的格点数 ...

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

12

帖子

0

威望

131

eV
积分
143

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2025-6-3 23:13:15 | 只看该作者 Only view this author

sob老师,我使用了volmap生成的格点数据进行染色,格点数据范围从0-3.85,但我发现,比如使用RWB的色彩模式,我设置color scale data range跟格点数据的范围匹配,即0-3.85,此时显示的染色效果就只有蓝色,我将color scale data range换成-3.85到0则只有红色,想请教sob老师我该怎么样调整,能够在格点数据只有正值的情况下绘制出颜色渐变的热图呢

12

帖子

0

威望

131

eV
积分
143

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2025-6-3 23:14:56 | 只看该作者 Only view this author
这是color scale data range从0到3.85的染色效果

202506032314005362..png (584.13 KB, 下载次数 Times of downloads: 23)

202506032314005362..png

202506032314359600..png (630.8 KB, 下载次数 Times of downloads: 19)

202506032314359600..png

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

8#
发表于 Post on 2025-6-3 23:20:26 | 只看该作者 Only view this author
杂算者 发表于 2025-6-3 23:13
sob老师,我使用了volmap生成的格点数据进行染色,格点数据范围从0-3.85,但我发现,比如使用RWB的色彩模 ...

不需要上、下限非得分别跟格点数据最小、最大值一致。颜色区分不出来就减小色彩刻度范围
北京科音自然科学研究中心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

12

帖子

0

威望

131

eV
积分
143

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2025-6-3 23:33:20 | 只看该作者 Only view this author
sobereva 发表于 2025-6-3 23:20
不需要上、下限非得分别跟格点数据最小、最大值一致。颜色区分不出来就减小色彩刻度范围

感谢sob老师,我再尝试一下

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 06:34 , Processed in 0.152775 second(s), 23 queries , Gzip On.

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