计算化学公社

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

[VMD] 关于统计分子两两之间的质心距离脚本

[复制链接 Copy URL]

170

帖子

0

威望

770

eV
积分
940

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 yangjunfang 于 2022-1-13 16:01 编辑

各位大佬好,我想统计分析一下图中聚集体的邻近分子两两之间的质心距离,应该怎么处理呀?请问有什么脚本可以参考吗?

1642039670(1).png (44.61 KB, 下载次数 Times of downloads: 53)

1642039670(1).png

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

2#
发表于 Post on 2022-1-13 14:23:16 | 只看该作者 Only view this author
稍微改一下就行
在VMD中将距离较近的分子质心连线的脚本
http://sobereva.com/410
北京科音自然科学研究中心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

170

帖子

0

威望

770

eV
积分
940

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-1-13 16:00:00 | 只看该作者 Only view this author
本帖最后由 yangjunfang 于 2022-1-13 16:09 编辑
sobereva 发表于 2022-1-13 14:23
稍微改一下就行
在VMD中将距离较近的分子质心连线的脚本
http://sobereva.com/410

老师,是需要加上输出距离的命令行吗如图?按图中脚本进行统计之后出现了“measure center: bad weight sum, would cause divide by zero”这样的问题,请问脚本应该怎么修改呀

1642059231(1).png (56.52 KB, 下载次数 Times of downloads: 53)

重新编辑的脚本图1

重新编辑的脚本图1

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

4#
发表于 Post on 2022-1-13 16:05:08 | 只看该作者 Only view this author
yangjunfang 发表于 2022-1-13 16:00
老师,是需要加上输出距离的命令行吗?再有就是模拟完之后是有水分子的,我不想统计水分子之间的质心距离 ...


选择语句排除掉水

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

170

帖子

0

威望

770

eV
积分
940

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-1-13 16:12:57 | 只看该作者 Only view this author
sobereva 发表于 2022-1-13 16:05

选择语句排除掉水

好的老师,水的问题我解决了,还是那个脚本的问题加上输出质心距离的语句之后出现了measure center: bad weight sum, would cause divide by zero的问题,麻烦您看一下应该怎么改呀

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

6#
发表于 Post on 2022-1-13 17:00:15 | 只看该作者 Only view this author
yangjunfang 发表于 2022-1-13 16:12
好的老师,水的问题我解决了,还是那个脚本的问题加上输出质心距离的语句之后出现了measure center: bad  ...

加上puts语句来调试
比如可能选区什么都没选中
北京科音自然科学研究中心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

170

帖子

0

威望

770

eV
积分
940

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2022-1-13 18:57:27 | 只看该作者 Only view this author
sobereva 发表于 2022-1-13 17:00
加上puts语句来调试
比如可能选区什么都没选中

老师,您能帮忙看一下这个脚本中哪里有问题吗,我加了puts语句进行了调试,但是还是出现同样的问题

com_distance.tcl

925 Bytes, 下载次数 Times of downloads: 50

53

帖子

0

威望

787

eV
积分
840

Level 4 (黑子)

8#
发表于 Post on 2022-10-12 22:28:42 | 只看该作者 Only view this author
yangjunfang 发表于 2022-1-13 18:57
老师,您能帮忙看一下这个脚本中哪里有问题吗,我加了puts语句进行了调试,但是还是出现同样的问题

我也碰到同样的问题
仔细看脚本,set sel [atomselect top "resid $ires"]这里,残基是从1号开始的,我的残基是从0号开始的,报同样的错,改了脚本已经解决。

15

帖子

0

威望

227

eV
积分
242

Level 3 能力者

9#
发表于 Post on 2023-4-24 11:39:22 | 只看该作者 Only view this author
如果只算两个指定的分子之间的质心距离该怎么改呢?小白不会使用vmd脚本

170

帖子

0

威望

770

eV
积分
940

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2023-5-17 11:58:33 | 只看该作者 Only view this author
llh 发表于 2023-4-24 11:39
如果只算两个指定的分子之间的质心距离该怎么改呢?小白不会使用vmd脚本

可以用gromacs中distance命令直接计算

55

帖子

0

威望

195

eV
积分
250

Level 3 能力者

11#
发表于 Post on 2024-7-9 20:04:59 | 只看该作者 Only view this author
yangjunfang 发表于 2022-1-13 18:57
老师,您能帮忙看一下这个脚本中哪里有问题吗,我加了puts语句进行了调试,但是还是出现同样的问题

请教一下,这个脚本是统计任意分子之间的距离吗,我体系有3种分子只想统计其中一种分子任意两个分子之间的最小距离要怎么修改呀

本版积分规则 Credits rule

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

GMT+8, 2026-2-27 23:16 , Processed in 0.301368 second(s), 23 queries , Gzip On.

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