计算化学公社

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

[VMD] 使用VMD计算径向分布函数,最终稳定的数值的意义是什么

[复制链接 Copy URL]

8

帖子

0

威望

108

eV
积分
116

Level 2 能力者

请问径向分布函数最终在什么情况下是趋于1,什么情况是趋于其他数值,比如说0.
我使用的软件是CPMD,体系是一个褐煤大分子和一些水分子。计算之后使用VMD软件的g(r)功能进行分析,选中了褐煤分子中羧基中的一个O,计算这个O原子与体系中其他所有的H原子的径向分布函数,发现最终趋向于0,想请问一下老师这是为什么,我看了一些关于径向分布函数的解释还是不明白。

12

帖子

0

威望

134

eV
积分
146

Level 2 能力者

2#
发表于 Post on 2018-4-6 10:38:48 | 只看该作者 Only view this author
你没有进行RDF的归一化,归一化后RDF应该都趋近于1

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

3#
发表于 Post on 2018-4-6 15:35:02 | 只看该作者 Only view this author
径向分布函数是每个球层位置的被统计的组的原子的数密度与平均数密度的比值,在距离煤很远的地方,水的分布情况和水的体相一样,因此除以平均密度,结果应当为1。
不满足这个情况时,根据体系实际情况具体情况具体分析。也许是在计算时rdf没有考虑PBC?这样的话远处区域已经超过了盒子,rdf肯定为0。
北京科音自然科学研究中心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

威望

108

eV
积分
116

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2018-4-7 22:01:54 | 只看该作者 Only view this author
木易金豆 发表于 2018-4-6 10:38
你没有进行RDF的归一化,归一化后RDF应该都趋近于1

你好,我不是很明白。能否告诉一下什么是归一化,如何进行归一化啊

8

帖子

0

威望

108

eV
积分
116

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2018-4-7 22:13:28 | 只看该作者 Only view this author
sobereva 发表于 2018-4-6 15:35
径向分布函数是每个球层位置的被统计的组的原子的数密度与平均数密度的比值,在距离煤很远的地方,水的分布 ...

你好,sob老师,很感谢你能回复我的问题。不过我设置了盒子的大小,然后计算的时候选择了pbc。如前两个图所示,结果还是最终趋于0了。

HAF68F}CGV_0G%1%2JZG$ES.png (5.99 KB, 下载次数 Times of downloads: 98)

HAF68F}CGV_0G%1%2JZG$ES.png

2.png (10.51 KB, 下载次数 Times of downloads: 82)

2.png

3.png (6.29 KB, 下载次数 Times of downloads: 87)

3.png

1171

帖子

7

威望

6858

eV
积分
8169

Level 6 (一方通行)

6#
发表于 Post on 2018-4-8 01:00:27 | 只看该作者 Only view this author
VMD算的是pair rdf,你要算的是simple rdf,不太一样

不过你结果等于0的原因是你的max r等于10

8

帖子

0

威望

108

eV
积分
116

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2018-4-16 15:26:19 | 只看该作者 Only view this author
fhh2626 发表于 2018-4-8 01:00
VMD算的是pair rdf,你要算的是simple rdf,不太一样

不过你结果等于0的原因是你的max r等于10

你好,感谢你能回答我的问题。不过为什么max r等于10最后的结果就是0啊,能否给简单的解释一下

1171

帖子

7

威望

6858

eV
积分
8169

Level 6 (一方通行)

8#
发表于 Post on 2018-4-16 16:14:22 | 只看该作者 Only view this author
mazhi19950907 发表于 2018-4-16 15:26
你好,感谢你能回答我的问题。不过为什么max r等于10最后的结果就是0啊,能否给简单的解释一下

因为你只统计到10啊。。

72

帖子

0

威望

402

eV
积分
474

Level 3 能力者

9#
发表于 Post on 2018-5-23 09:40:46 | 只看该作者 Only view this author
你好,看你用vmd做的径向分布函数,我的也是趋于零,你找到解决的原因了吗,可以回答下吗,谢谢,那个所谓的归一化你是怎么解决的,谢谢

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

10#
发表于 Post on 2018-5-23 18:47:07 | 只看该作者 Only view this author
chw985550192 发表于 2018-5-23 09:40
你好,看你用vmd做的径向分布函数,我的也是趋于零,你找到解决的原因了吗,可以回答下吗,谢谢,那个所谓 ...

确认你输入pbc box命令的时候能正常显示出盒子边框,否则当前体系并没法利用周期边界条件
北京科音自然科学研究中心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

5

帖子

0

威望

49

eV
积分
54

Level 2 能力者

11#
发表于 Post on 2019-6-5 20:00:54 | 只看该作者 Only view this author
sobereva 发表于 2018-5-23 18:47
确认你输入pbc box命令的时候能正常显示出盒子边框,否则当前体系并没法利用周期边界条件

请问VMD考虑pbc时只能算正交格矢的吗,如何算六方的

5

帖子

0

威望

49

eV
积分
54

Level 2 能力者

12#
发表于 Post on 2019-6-5 20:01:29 | 只看该作者 Only view this author
mazhi19950907 发表于 2018-4-7 22:13
你好,sob老师,很感谢你能回复我的问题。不过我设置了盒子的大小,然后计算的时候选择了pbc。如前两个图所 ...

请问第一张图是在哪里可以设置

513

帖子

0

威望

4384

eV
积分
4897

Level 6 (一方通行)

truffle

13#
发表于 Post on 2019-10-17 09:36:34 | 只看该作者 Only view this author
fhh2626 发表于 2018-4-8 01:00
VMD算的是pair rdf,你要算的是simple rdf,不太一样

不过你结果等于0的原因是你的max r等于10

想请教一下您pair rdf和simple rdf有什么区别呢 谢谢~
No problem is insoluble in all conceivable circumstances.

1171

帖子

7

威望

6858

eV
积分
8169

Level 6 (一方通行)

14#
发表于 Post on 2019-10-17 10:57:56 | 只看该作者 Only view this author
naoki 发表于 2019-10-17 09:36
想请教一下您pair rdf和simple rdf有什么区别呢 谢谢~

pair rdf就是设定两种类型原子,然后统计出任意两个这两种类型原子的距离,结果是这个距离的分布
simple rdf就是给一个点或者一个原子,算围绕这个点的其他原子的分布

评分 Rate

参与人数
Participants 2
eV +4 收起 理由
Reason
ocbrother + 2 谢谢
少年爱吃地三鲜 + 2 赞!

查看全部评分 View all ratings

513

帖子

0

威望

4384

eV
积分
4897

Level 6 (一方通行)

truffle

15#
发表于 Post on 2019-10-17 14:24:09 | 只看该作者 Only view this author
fhh2626 发表于 2019-10-17 10:57
pair rdf就是设定两种类型原子,然后统计出任意两个这两种类型原子的距离,结果是这个距离的分布
simple ...

谢谢啦!
No problem is insoluble in all conceivable circumstances.

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 01:20 , Processed in 0.185124 second(s), 24 queries , Gzip On.

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