计算化学公社

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

[GROMACS] gmx clustsize分析结果中的尺寸分布怎么计算的

[复制链接 Copy URL]

125

帖子

0

威望

2930

eV
积分
3055

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
如题,这个尺寸分布是怎么计算的,有没有公式可循?可否进行类似不同体系进行对比?

6万

帖子

99

威望

6万

eV
积分
125151

管理员

公社社长

2#
发表于 Post on 2018-9-29 15:55:22 | 只看该作者 Only view this author
只要原子间距离小于3.5埃(数值可以用-cut设定),就被认为处于同一个簇中。
可以不同体系间对比
北京科音自然科学研究中心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

125

帖子

0

威望

2930

eV
积分
3055

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2018-9-29 16:04:56 | 只看该作者 Only view this author
sobereva 发表于 2018-9-29 15:55
只要原子间距离小于3.5埃(数值可以用-cut设定),就被认为处于同一个簇中。
可以不同体系间对比

谢谢sob老师,可能我没说清楚,我指的是生成的文件中,histo-clust.xvg文件,里面的每个聚集体后面的数字是怎么算的?

6万

帖子

99

威望

6万

eV
积分
125151

管理员

公社社长

4#
发表于 Post on 2018-9-29 21:43:53 | 只看该作者 Only view this author
苏玖染 发表于 2018-9-29 16:04
谢谢sob老师,可能我没说清楚,我指的是生成的文件中,histo-clust.xvg文件,里面的每个聚集体后面的数字 ...

诸如3后面的数字,就是指模拟过程中含有3个分子的簇的平均数目
北京科音自然科学研究中心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

125

帖子

0

威望

2930

eV
积分
3055

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2018-9-30 09:23:30 | 只看该作者 Only view this author
sobereva 发表于 2018-9-29 21:43
诸如3后面的数字,就是指模拟过程中含有3个分子的簇的平均数目

好的,十分感谢sob老师

125

帖子

0

威望

2930

eV
积分
3055

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2019-11-7 19:49:16 | 只看该作者 Only view this author
sobereva 发表于 2018-9-29 15:55
只要原子间距离小于3.5埃(数值可以用-cut设定),就被认为处于同一个簇中。
可以不同体系间对比

再请问下sob老师,这个cutoff选为3.5埃是怎么定的,或者说这个3.5埃是否是某个有特定物理意义的值。
我刚投了一篇文章,审稿人质疑这个cutoff。

我做的主要是有机小分子在水溶液中的团聚行为

6万

帖子

99

威望

6万

eV
积分
125151

管理员

公社社长

7#
发表于 Post on 2019-11-11 11:15:06 | 只看该作者 Only view this author
苏玖染 发表于 2019-11-7 19:49
再请问下sob老师,这个cutoff选为3.5埃是怎么定的,或者说这个3.5埃是否是某个有特定物理意义的值。
我 ...

具体怎么来的开发者才知道
但从碳的Bondi范德华半径为1.7埃来看,这个3.5埃的标准差不多是倆碳的范德华半径和,因此有一定道理
北京科音自然科学研究中心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

125

帖子

0

威望

2930

eV
积分
3055

Level 5 (御坂)

8#
 楼主 Author| 发表于 Post on 2019-11-11 11:50:31 | 只看该作者 Only view this author
sobereva 发表于 2019-11-11 11:15
具体怎么来的开发者才知道
但从碳的Bondi范德华半径为1.7埃来看,这个3.5埃的标准差不多是倆碳的范德华 ...

谢谢sob老师,我问了其它老师,他说,这个3.5A是SPC水之间rdf峰的距离,所以具体问题还是看实际情况选择cutoff值。我感觉我这里应该用3.4A好一些

6万

帖子

99

威望

6万

eV
积分
125151

管理员

公社社长

9#
发表于 Post on 2019-11-11 12:07:08 | 只看该作者 Only view this author
苏玖染 发表于 2019-11-11 11:50
谢谢sob老师,我问了其它老师,他说,这个3.5A是SPC水之间rdf峰的距离,所以具体问题还是看实际情况选择c ...

这是我用SPC/E跑的,没看出有3.5埃


如果是指O-O距离,拿这个当判据也意义不大,毕竟不同体系相互作用特征明显不同。


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

125

帖子

0

威望

2930

eV
积分
3055

Level 5 (御坂)

10#
 楼主 Author| 发表于 Post on 2019-11-11 20:58:18 | 只看该作者 Only view this author
本帖最后由 苏玖染 于 2019-11-11 21:02 编辑
sobereva 发表于 2019-11-11 12:07
这是我用SPC/E跑的,没看出有3.5埃

嗯嗯,还是要考虑不同体系不一样吧,我猜测应该是对应于表面算rdf吧,就是在算rdf的时候加上-surf mol,因为毕竟single linkage判断方法是基于分子表面原子的

22

帖子

0

威望

89

eV
积分
111

Level 2 能力者

11#
发表于 Post on 2020-1-6 17:21:52 | 只看该作者 Only view this author
苏玖染 发表于 2019-11-11 20:58
嗯嗯,还是要考虑不同体系不一样吧,我猜测应该是对应于表面算rdf吧,就是在算rdf的时候加上-surf mol, ...

请问-cut的值不是径向分布函数的第一个极小值处x轴的坐标值吗?

125

帖子

0

威望

2930

eV
积分
3055

Level 5 (御坂)

12#
 楼主 Author| 发表于 Post on 2020-1-7 09:20:10 | 只看该作者 Only view this author
LuckyLynn 发表于 2020-1-6 17:21
请问-cut的值不是径向分布函数的第一个极小值处x轴的坐标值吗?

上述只是我的猜测,好像没有定论

22

帖子

0

威望

89

eV
积分
111

Level 2 能力者

13#
发表于 Post on 2020-1-7 10:45:41 | 只看该作者 Only view this author
苏玖染 发表于 2020-1-7 09:20
上述只是我的猜测,好像没有定论

我还有问题想请教你,就是-cut的值你最终是怎么确定的?它的最大值是0.35吗?如果要调那个值的话只能往小调吗?

125

帖子

0

威望

2930

eV
积分
3055

Level 5 (御坂)

14#
 楼主 Author| 发表于 Post on 2020-1-7 16:28:37 | 只看该作者 Only view this author
LuckyLynn 发表于 2020-1-7 10:45
我还有问题想请教你,就是-cut的值你最终是怎么确定的?它的最大值是0.35吗?如果要调那个值的话只能往小 ...

根据你的体系定,你可以用不同的cutoff值都算算看,然后看结果哪个更符合你的体系

22

帖子

0

威望

89

eV
积分
111

Level 2 能力者

15#
发表于 Post on 2020-1-7 17:05:48 | 只看该作者 Only view this author
苏玖染 发表于 2020-1-7 16:28
根据你的体系定,你可以用不同的cutoff值都算算看,然后看结果哪个更符合你的体系

好的,谢谢

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 06:51 , Processed in 0.258817 second(s), 22 queries , Gzip On.

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