计算化学公社

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

[综合交流] 如何作径向分布函数图及其积分图的叠图?

[复制链接 Copy URL]

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

如图所示的图是怎么做出来的呢?
单独的径向分布函数图及其积分图都会在origin做,但是两者的叠图是怎么做出来的呢?谢谢!

径向分布函数图及积分图的叠图.png (82.81 KB, 下载次数 Times of downloads: 81)

径向分布函数图及积分图的叠图.png

5万

帖子

99

威望

5万

eV
积分
112353

管理员

公社社长

2#
发表于 Post on 2017-3-6 16:24:56 | 只看该作者 Only view this author
基本的origin操作问题而已。设两个层,每个层用不同的数据集来作图。不会的话google
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

3#
 楼主 Author| 发表于 Post on 2017-3-9 15:43:17 | 只看该作者 Only view this author
sobereva 发表于 2017-3-6 16:24
基本的origin操作问题而已。设两个层,每个层用不同的数据集来作图。不会的话google

如图1所示红色的是径向分布,黑色的是积分图,横坐标的单位是Å。我先做的径向图然后再对径向图进行积分(analysis-calculus-integrate)得到单独的积分图及其数据,然后把两图的数据放一块做成如图所示的图,为什么我的图不能做成别人那样的?是我作图方法错了还是我的径向分布图有问题呢?谢谢!

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

图1.png

5万

帖子

99

威望

5万

eV
积分
112353

管理员

公社社长

4#
发表于 Post on 2017-3-9 17:52:53 | 只看该作者 Only view this author
15whyhq 发表于 2017-3-9 15:43
如图1所示红色的是径向分布,黑色的是积分图,横坐标的单位是Å。我先做的径向图然后再对径向图进行 ...


你的rdf看起来有点怪,或者统计的方式就不对。积分曲线相对于rdf曲线来看倒没什么错误
你可以先拿纯水体系的O-O rdf练手
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

206

帖子

3

威望

2814

eV
积分
3080

Level 5 (御坂)

5#
发表于 Post on 2017-3-10 11:02:09 | 只看该作者 Only view this author
sobereva 发表于 2017-3-9 17:52
你的rdf看起来有点怪,或者统计的方式就不对。积分曲线相对于rdf曲线来看倒没什么错误
你可以先拿纯水 ...

sob老师,我算出来的rdf(蛋白与溶剂)的好像和上面的也差不多(S型曲线)。。是不是也有点怪?
主攻: 蛋白-蛋白对接,蛋白de novo设计、蛋白结构建模,抗体设计等方向。Rosetta/PyRosetta

5万

帖子

99

威望

5万

eV
积分
112353

管理员

公社社长

6#
发表于 Post on 2017-3-10 11:07:43 | 只看该作者 Only view this author
kunkun 发表于 2017-3-10 11:02
sob老师,我算出来的rdf(蛋白与溶剂)的好像和上面的也差不多(S型曲线)。。是不是也有点怪?


对这种大分子的全部原子直接做rdf图其实根本没有什么意义,所有化学环境的原子都混淆在一起,什么也考查不了,rdf也表现不出什么结构特征来。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

7#
 楼主 Author| 发表于 Post on 2017-3-20 17:47:34 | 只看该作者 Only view this author
求动力学模拟的中心分子周围合理的溶剂分子个数  是像文献中的图那样用动力学的rdf图求积分图后再把rdf图及积分图叠在一块,从rdf图第一个峰后的最低点的横坐标值对应的积分图的纵坐标就是配位数吗?  这个操作过程对吗?

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

8#
 楼主 Author| 发表于 Post on 2017-3-20 21:52:47 | 只看该作者 Only view this author
如图2 所示是我做的径向分布函数及其积分的叠图,横坐标的单位是nm,如果要把横纵标的单位改为Å,是不是纵坐标也得*10?

图2.png (12.14 KB, 下载次数 Times of downloads: 68)

图2.png

5万

帖子

99

威望

5万

eV
积分
112353

管理员

公社社长

9#
发表于 Post on 2017-3-21 01:30:57 | 只看该作者 Only view this author
积分值乘以体系中被计算的粒子平均数密度才是配位数,从rdf的基本定义上就知道
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

10#
 楼主 Author| 发表于 Post on 2017-3-21 15:30:31 | 只看该作者 Only view this author
本帖最后由 15whyhq 于 2017-3-21 15:35 编辑

这篇文献里的求配位数的方法对吗?明明是确定中心QM部分周围合理的离子液体分子个数,为什么做的rdf图是阳离子与阴离子的呢?为什么不是QM的质心与环境中的离子液体的质心呢?

图.png (82.7 KB, 下载次数 Times of downloads: 68)

图.png

QMMM MD CO2 absorption-计算部分.pdf

1.62 MB, 下载次数 Times of downloads: 486

supporting-QMMM MD CO2 absorption-计算部分.pdf

4.55 MB, 下载次数 Times of downloads: 1713

5万

帖子

99

威望

5万

eV
积分
112353

管理员

公社社长

11#
发表于 Post on 2017-3-21 16:01:13 | 只看该作者 Only view this author
15whyhq 发表于 2017-3-21 15:30
这篇文献里的求配位数的方法对吗?明明是确定中心QM部分周围合理的离子液体分子个数,为什么做的rdf图是阳 ...

我没时间具体读文章,但看样子人家就是要考察体系中阴阳离子质心之间的rdf,就这样做并没有问题
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

586

帖子

0

威望

4191

eV
积分
4777

Level 6 (一方通行)

12#
 楼主 Author| 发表于 Post on 2017-3-21 19:04:55 | 只看该作者 Only view this author
本帖最后由 15whyhq 于 2017-3-21 19:06 编辑
sobereva 发表于 2017-3-21 16:01
我没时间具体读文章,但看样子人家就是要考察体系中阴阳离子质心之间的rdf,就这样做并没有问题

哦,那可以帮忙解释下为什么rdf的第一个峰顶后的峰谷做垂线与rdf的积分图交点所对应的纵坐标就是配位数呢?给我一些相关的资料看也行的,非常感谢!

5万

帖子

99

威望

5万

eV
积分
112353

管理员

公社社长

13#
发表于 Post on 2017-3-21 19:06:41 | 只看该作者 Only view this author
15whyhq 发表于 2017-3-21 19:04
哦,那可以帮忙解释下为什么rdf的第一个峰顶后的峰谷做垂线与rdf的积分图交点所对应的纵坐标就是配位数呢 ...

随便找点rdf的资料,一看rdf定义立马就明白了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 19:56 , Processed in 0.206593 second(s), 24 queries , Gzip On.

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