查看: 550  |  回复: 11

红精灵

新虫 (初入文坛)

[求助] 全原子力场与RDF的分析是否有关系?如何对原子基团或者分子进行RDF分析?

吉布斯系统Monte Carlo模拟的新手,想问大侠们一个问题。选择全原子力场NERD进行的轻烃类气液相平衡的模拟,想对平衡的液相进行RDF分析,不知是分析分子(或者原子基团)之间的RDF还是原子之间的RDF?这与所选的力场是否有关系,选择全原子力场就必须对原子基团之间进行分析还是可以对单一原子分析?
还有一个很头疼的问题,选择原子作为SETS比较容易,但是选择分子或者基团的话,通过find patterns 可以选择分子,但不能选择基团,而且 这种sets 在RDF 或者是PCF中不能分析啊, 如何选择原子基团或者是分子并进行RDF分析呢?
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

飞刀
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yahoohoo

铁杆木虫 (著名写手)

【答案】应助回帖

红精灵(金币+6): 给出关键建议 2011-06-29 07:04:44
引用回帖:
Originally posted by 红精灵 at 2011-06-28 22:25:26:
吉布斯系统Monte Carlo模拟的新手,想问大侠们一个问题。选择全原子力场NERD进行的轻烃类气液相平衡的模拟,想对平衡的液相进行RDF分析,不知是分析分子(或者原子基团)之间的RDF还是原子之间的RDF?这与所选的力 ...

分析分子还是原子基团的RDF完全取决于你需要得到什么样的信息。与所用的力场是有关系的,不过合理的力场应该能重现实验结果。
2楼2011-06-29 01:32:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

红精灵

新虫 (初入文坛)

引用回帖:
Originally posted by yahoohoo at 2011-06-29 01:32:28:
分析分子还是原子基团的RDF完全取决于你需要得到什么样的信息。与所用的力场是有关系的,不过合理的力场应该能重现实验结果。


嗯,就是说用如果用NERD(全原子力场)是可以分析原子之间的RDF的,我想分析基团之间的RDF,但是不会选中原子基团作为SETS,使用的是MS4.0中Forcite/Anaysisi/Radial distribution function,用edit/select atoms还是edit/find patterns来选择基团? 具体怎么怎么选? 大侠能指点一下吗?
还想问大侠一个问题。下图是我分析原子之间的一个RDF,怎么调整都是很粗糙,我看文献上的都是很平滑的漂亮的曲线,由于毕业限制需要发文章,所以很急,能糊弄过去就行啊,希望您能提出宝贵意见.
提这么多的问题,真是麻烦大侠了,大侠肯相助,我可以把所以的金币都赠送给你
不胜感激!





飞刀
3楼2011-06-29 07:40:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yahoohoo

铁杆木虫 (著名写手)

★ ★ ★
御剑江湖(金币+3): 谢谢 2011-06-29 18:55:46
引用回帖:
Originally posted by 红精灵 at 2011-06-29 07:40:52:

嗯,就是说用如果用NERD(全原子力场)是可以分析原子之间的RDF的,我想分析基团之间的RDF,但是不会选中原子基团作为SETS,使用的是MS4.0中Forcite/Anaysisi/Radial distribution function,用edit/se ...

我没用过MS,不知道该选择哪个,不过你倒是可以看看MS的使用说明啊。

你所计算的RDF看上去是有些粗糙,你模拟的时间是否够长呢,或者说样本有多少?
4楼2011-06-29 13:10:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

红精灵

新虫 (初入文坛)

哦,你说的我都考虑过呢,样本已经是512个分子了,时间总步是1000 000步了,就是以上图的结果了。
谢谢你了
飞刀
5楼2011-06-29 16:07:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yahoohoo

铁杆木虫 (著名写手)

★ ★ ★
御剑江湖(金币+3): 谢谢 2011-06-30 10:03:35
引用回帖:
Originally posted by 红精灵 at 2011-06-29 16:07:20:
哦,你说的我都考虑过呢,样本已经是512个分子了,时间总步是1000 000步了,就是以上图的结果了。
谢谢你了

也就是说你计算RDF时统计了512个分子?

时间步长1M,那么独立的样本数是多少?也就是每个多少步你抽取一个样本呢?
6楼2011-06-29 20:43:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yahoohoo

铁杆木虫 (著名写手)

★ ★ ★
御剑江湖(金币+3): 谢谢 2011-06-30 10:03:43
引用回帖:
Originally posted by 红精灵 at 2011-06-29 16:07:20:
哦,你说的我都考虑过呢,样本已经是512个分子了,时间总步是1000 000步了,就是以上图的结果了。
谢谢你了

你可以考虑使用更小的距离间隔统计RDF,这样你的图中就有更多的数据点。

检查你模拟的体系是否已经平衡了,如果是,那可以增大统计样本的数目以提高统计精度。
7楼2011-06-29 20:53:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

红精灵

新虫 (初入文坛)

引用回帖:
Originally posted by yahoohoo at 2011-06-29 20:53:29:
你可以考虑使用更小的距离间隔统计RDF,这样你的图中就有更多的数据点。

检查你模拟的体系是否已经平衡了,如果是,那可以增大统计样本的数目以提高统计精度。

我是先对体系跑的NPT动力学,再进行的RDF分析,上图是动力学里的设置,下图是在RDF分析时的数据,你说的样本是上图中的frame output every:5000?这个更小的间距是在这里设定的吗?
还有上次说错了,不是跑了1000 000步 仅是10000步,体系的分子数是512
请赐教



飞刀
8楼2011-06-29 21:47:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yahoohoo

铁杆木虫 (著名写手)

★ ★ ★
御剑江湖(金币+3): 谢谢 2011-06-30 10:03:53
引用回帖:
Originally posted by 红精灵 at 2011-06-29 21:47:18:
我是先对体系跑的NPT动力学,再进行的RDF分析,上图是动力学里的设置,下图是在RDF分析时的数据,你说的样本是上图中的frame output every:5000?这个更小的间距是在这里设定的吗?
还有上次说错了,不是跑了 ...

没用过MS模拟,你自己应该清楚各项对应什么,不然计算出的东西没有任何意义。

就字面意思,第二个图中的interval对应距离的增量,可以设置小点。你可以增加模拟的步长再统计RDF看看。
9楼2011-06-30 02:57:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

红精灵

新虫 (初入文坛)


御剑江湖(金币+1): 鼓励 2011-06-30 10:04:10
引用回帖:
Originally posted by yahoohoo at 2011-06-30 02:57:11:
没用过MS模拟,你自己应该清楚各项对应什么,不然计算出的东西没有任何意义。

就字面意思,第二个图中的interval对应距离的增量,可以设置小点。你可以增加模拟的步长再统计RDF看看。

嗯,总之谢谢大侠,我再增加以下模拟步长试试吧,第二图中的interva对应的量只是对已经输出的图进行的调整,不会增加动力学步骤,之前发的三个图就是已经用第二图中的cutoff 和 interval 调整过的几乎是最好看一些的结果了
飞刀
10楼2011-06-30 07:27:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 红精灵 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 看书可以让我消除一切不良情绪 +5 wait2022 2021-10-21 10/500 2021-10-22 00:24 by 仰望云海
[论文投稿] 作图 +5 suohongyan 2021-10-21 5/250 2021-10-21 23:12 by 刘冲~
[找工作] 大家觉得这两个offer,哪个更好?暨南大学珠海校区与中科院深圳先进院 +17 UGentgreat 2021-10-16 38/1900 2021-10-21 21:58 by UGentgreat
[硕博家园] 换一种方式,科研并不孤独 +6 原是一条鱼 2021-10-21 9/450 2021-10-21 20:11 by 也许爱
[找工作] 内卷这么厉害吗 +16 cslglgf 2021-10-16 48/2400 2021-10-21 13:24 by 葡萄啊
[考博] 导师回复建议多联系几个学校,是不是凉了 +12 957786038 2021-10-20 21/1050 2021-10-21 10:17 by 957786038
[基金申请] 如何回避恶意基金评审人-2022 +7 hxy_001 2021-10-20 9/450 2021-10-21 09:28 by 快乐的水瓶
[论文投稿] IEEE JSAC 一审全部小修 二审接收概率? +6 lty_sjtu 2021-10-20 11/550 2021-10-21 09:23 by lyf1850
[博后之家] 国内博士后出站后可以去北理工做预聘助理教授,或副预聘副教授吗? +3 我爱_学术 2021-10-20 3/150 2021-10-20 22:45 by H二十八画生
[论文投稿] 因为语言被拒有没有可能挽回 20+4 珣Meng 2021-10-16 14/700 2021-10-20 18:17 by blade0070
[基金申请] 任务计划书状态等待基金委审核 +3 jnwek123 2021-10-20 6/300 2021-10-20 17:27 by jnwek123
[论文投稿] 一审:两个接收,一个修,一个拒 +8 hui尘 2021-10-19 11/550 2021-10-20 15:04 by lty_sjtu
[有机交流] 求助反应机理谢谢 +3 kiss20085537 2021-10-18 10/500 2021-10-20 12:33 by hyacinth1944
[硕博家园] 论坛里男女征友贴 +10 不屈的心 2021-10-16 15/750 2021-10-19 17:30 by 毕加思索
[基金申请] 请问关于面上项目预算申请书的填写问题,谢谢 +3 wittyboy 2021-10-19 3/150 2021-10-19 15:13 by rhhua11111
[论文投稿] 关于农学类生物炭研究与应用类OA期刊 +3 xyx2021 2021-10-19 4/200 2021-10-19 12:27 by xyx2021
[找工作] 硕士有什么兼职?国内期刊审稿挣钱吗? +3 996247712 2021-10-16 6/300 2021-10-18 12:00 by kataya123
[论文投稿] 投稿状态由Pending Editor Decision,直接变成Sent to Reviewers啥情况? +3 linyi0714 2021-10-17 3/150 2021-10-18 07:58 by 国际科学编辑
[考博] 南京土壤所和南京地湖所 哪个比较强一点 哪个方向更有前途 +4 727395684 2021-10-17 4/200 2021-10-17 11:28 by 矛dei到
[基金申请] 青基状态没变 +5 sicencess 2021-10-16 6/300 2021-10-17 08:32 by jnwek123
信息提示
请填处理意见