计算化学公社

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

[GROMACS] 计算真空盒子内的液滴里各种成分的径向分布函数图的问题

[复制链接 Copy URL]

107

帖子

0

威望

467

eV
积分
574

Level 4 (黑子)

各位老师好,这是一篇文献中对于真空盒子内的液滴里各种成分的径向分布函数图,请问这个曲线可以怎么得到呢?谢谢
文献地址 pubs.acs.org/doi/10.1021/acs.analchem.6b00650

202212161626083386..png (81.71 KB, 下载次数 Times of downloads: 13)

202212161626083386..png

310

帖子

0

威望

1503

eV
积分
1813

Level 5 (御坂)

2#
发表于 Post on 2022-12-17 12:19:35 | 只看该作者 Only view this author
这个图的做法是先算出中心的位置,然后计算各组分中的原子到中心的距离,分别统计画出直方图并归一化,由于在球坐标下,均匀分布的径向分别是和4πr^2成正比的,所以要再除以这个因子,最后为了图片美观,作者把蛋白质的径向分布函数数值按比例减小了。

107

帖子

0

威望

467

eV
积分
574

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-12-18 10:23:52 | 只看该作者 Only view this author
Frozen-Penguin 发表于 2022-12-17 12:19
这个图的做法是先算出中心的位置,然后计算各组分中的原子到中心的距离,分别统计画出直方图并归一化,由于 ...

谢谢老师,请问怎么算中心的位置再算各组分中的原子到中心的距离呀?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2022-12-18 10:50:31 | 只看该作者 Only view this author
961106 发表于 2022-12-18 10:23
谢谢老师,请问怎么算中心的位置再算各组分中的原子到中心的距离呀?

中心的位置在哪看具体体系,自己最清楚,既可能是固定的也可能不固定

距离中心的距离自己写VMD tcl脚本统计就完了。这是北京科音分子动力学与GROMACS培训班(http://www.keinsci.com/workshop/KGMX_content.html)里的相关例子,课上非常详细讲怎么编写脚本







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

107

帖子

0

威望

467

eV
积分
574

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-12-18 15:26:39 | 只看该作者 Only view this author
sobereva 发表于 2022-12-18 10:50
中心的位置在哪看具体体系,自己最清楚,既可能是固定的也可能不固定

距离中心的距离自己写VMD tcl脚 ...

谢谢老师,我参加过您的培训,刚才翻看了讲义里的这部分内容,我的体系也是气相中的液滴,但是我添加了电场,所以液滴在模拟过程中一直发生移动,请问老师对于这种体系中心移动的情况,如何在脚本文件里修改呢?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2022-12-18 18:22:21 | 只看该作者 Only view this author
961106 发表于 2022-12-18 15:26
谢谢老师,我参加过您的培训,刚才翻看了讲义里的这部分内容,我的体系也是气相中的液滴,但是我添加了电 ...

求体系的几何中心作为被统计的中心,或者trjconv结合-center把体系居中
北京科音自然科学研究中心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

8

帖子

0

威望

155

eV
积分
163

Level 3 能力者

7#
发表于 Post on 2023-3-30 12:18:21 | 只看该作者 Only view this author
sobereva 发表于 2022-12-18 10:50
中心的位置在哪看具体体系,自己最清楚,既可能是固定的也可能不固定

距离中心的距离自己写VMD tcl脚 ...

老师,我用VMD TkConsole导入报错 invalid conmmand name "rangenum"

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

8#
发表于 Post on 2023-4-1 05:19:10 | 只看该作者 Only view this author
Weijie 发表于 2023-3-30 12:18
老师,我用VMD TkConsole导入报错 invalid conmmand name "rangenum"

没有以正确方式运行
不要用“导入”这种含糊其辞的描述
北京科音自然科学研究中心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

8

帖子

0

威望

155

eV
积分
163

Level 3 能力者

9#
发表于 Post on 2023-4-2 11:39:21 | 只看该作者 Only view this author
sobereva 发表于 2023-4-1 05:19
没有以正确方式运行
不要用“导入”这种含糊其辞的描述

好的老师,谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 21:10 , Processed in 0.169112 second(s), 23 queries , Gzip On.

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