24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2111  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 红精灵 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 复试调剂 +8 春日来信- 2026-04-03 8/400 2026-04-05 18:58 by 蓝云思雨
[考研] 材料调剂 +11 壹贰贰亿 2026-04-04 11/550 2026-04-05 18:47 by 蓝云思雨
[考研] 280求调剂 +4 李rien 2026-04-04 4/200 2026-04-05 18:44 by imissbao
[考研] 调剂 +4 不逢春 2026-04-05 4/200 2026-04-05 18:07 by ms629
[考研] 332求调剂 +11 小小孟... 2026-04-05 11/550 2026-04-05 18:04 by 猪会飞
[考研] 一志愿南航,数一英一学硕317求调剂!! +5 Acaciad 2026-04-04 5/250 2026-04-05 12:31 by 搏击518
[考研] 271分求调剂学校 +12 zph158488! 2026-04-02 13/650 2026-04-05 10:13 by lqwchd
[考研] 085602调剂 初试总分335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
[考研] 一志愿郑大0705求调剂 +3 橘十一 2026-04-02 4/200 2026-04-05 00:05 by chongya
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 281求调剂 +10 aaawhy 2026-04-03 10/500 2026-04-03 21:42 by lbsjt
[考研] 322求调剂 +4 FZAC123 2026-04-03 4/200 2026-04-03 20:55 by zhq0425
[考研] 一志愿武汉理工0856,初试334 +3 26考研材料 2026-04-02 3/150 2026-04-02 21:22 by dongzh2009
[考研] 261求B区调剂 +5 明仔· 2026-04-01 7/350 2026-04-02 11:17 by 邹尉尉
[考研] 0710生物学,325求调剂 +3 mkkkkkl 2026-04-01 3/150 2026-04-02 09:48 by Jaylen.
[考研] 土木304求调剂 +6 兔突突突, 2026-03-31 7/350 2026-04-02 09:06 by coolminer
[考研] 304求调剂 +12 素年祭语 2026-03-31 15/750 2026-04-01 22:41 by peike
[考研] 301求调剂 +8 axibli 2026-04-01 8/400 2026-04-01 09:51 by 我的船我的海
[考研] 调剂 +4 GK72 2026-03-30 4/200 2026-03-30 20:32 by dick_runner
[考研] 085602 化学工程专硕 340分求调剂 +4 qianbai11 2026-03-29 4/200 2026-03-30 11:34 by 唐沐儿
信息提示
请填处理意见