24小时热门版块排行榜    

查看: 1191  |  回复: 9

wj925806

铜虫 (小有名气)

[求助] 探究颗粒扩散时,颗粒是否给定热浴 已有1人参与

大家好。我最近在模拟颗粒在polymer中的扩散,想问下大家,此时颗粒是否要控制温度呢?如果控制颗粒的温度,是否会影响了它的扩散速度?  这里多谢大家!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lx_PICO

铁虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
wj925806(abinitio代发): 金币+1, 鼓励交流 2015-04-13 10:58:20
没温度粒子有速度么?
2楼2015-04-12 15:38:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wj925806

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by lx_PICO at 2015-04-12 15:38:27
没温度粒子有速度么?

它周围的polymer对它有作用力,所以它会有速度啊。因为观察的是扩散,所以我在想颗粒要不要给热浴
3楼2015-04-12 23:08:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
wj925806(abinitio代发): 金币+5, 鼓励交流 2015-04-13 10:58:35
wj925806: 金币+50, ★★★很有帮助 2015-04-13 12:08:55
一般研究扩散都是在某个温度下进行,因为温度对扩散系数的影响是非常大(扩散系数和温度的倒数呈指数型关系)。所以在研究扩散的模拟中都要使用热浴。通常的做法是在NPT或NVT下平衡后,在NVE下进行采样计算MSD或VACF来获得扩散系数(如果不是特别追求过程上的严谨也可以直接在NVT下采样计算扩散系数)。
4楼2015-04-12 23:17:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wj925806

铜虫 (小有名气)

引用回帖:
4楼: Originally posted by lsloneil at 2015-04-12 23:17:50
一般研究扩散都是在某个温度下进行,因为温度对扩散系数的影响是非常大(扩散系数和温度的倒数呈指数型关系)。所以在研究扩散的模拟中都要使用热浴。通常的做法是在NPT或NVT下平衡后,在NVE下进行采样计算MSD或VAC ...

嗯,你说的很对。在nvt下计算MSD时,我知道系统是要给定温度的,我观察的是颗粒在polymer中的扩散,颗粒是否也控制温度呢?我是担心对颗粒控制温度后,影响了它的扩散速度?不知道这样的理解对不对,还请多多指教啊
5楼2015-04-13 09:42:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

引用回帖:
5楼: Originally posted by wj925806 at 2015-04-12 13:42:19
嗯,你说的很对。在nvt下计算MSD时,我知道系统是要给定温度的,我观察的是颗粒在polymer中的扩散,颗粒是否也控制温度呢?我是担心对颗粒控制温度后,影响了它的扩散速度?不知道这样的理解对不对,还请多多指教啊...

当你的体系已经在NVT下平衡后,颗粒的温度也已经达到了你指定的温度,不信的话你可以算算颗粒的平均动能看是不是在3/2*k*T附近。

如果你在做NVT平衡的时候只对polymer而没有对颗粒施加热浴,这是不正确的。
6楼2015-04-13 10:23:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wj925806

铜虫 (小有名气)

引用回帖:
6楼: Originally posted by lsloneil at 2015-04-13 10:23:12
当你的体系已经在NVT下平衡后,颗粒的温度也已经达到了你指定的温度,不信的话你可以算算颗粒的平均动能看是不是在3/2*k*T附近。

如果你在做NVT平衡的时候只对polymer而没有对颗粒施加热浴,这是不正确的。...

我想再问一下,如果我用NVT系统采样计算扩散常数时,这个阶段颗粒是不是也施加热浴?
7楼2015-04-13 12:12:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

引用回帖:
7楼: Originally posted by wj925806 at 2015-04-12 16:12:00
我想再问一下,如果我用NVT系统采样计算扩散常数时,这个阶段颗粒是不是也施加热浴?...

是的,NVT热浴是对你的整个体系施加的,包括颗粒和polymer。

如我前面所述,更严谨的计算扩散系数的方法是先把体系用NVT热浴平衡,平衡之后在NVE系综下采样(也就是把热浴关掉)。但对于有些体系由于某些原因(比如力场),NVE下能量漂移会比较大,导致体系整体温度升高,这时最好还是用NVT采样。
8楼2015-04-13 12:42:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wj925806

铜虫 (小有名气)

引用回帖:
8楼: Originally posted by lsloneil at 2015-04-13 12:42:57
是的,NVT热浴是对你的整个体系施加的,包括颗粒和polymer。

如我前面所述,更严谨的计算扩散系数的方法是先把体系用NVT热浴平衡,平衡之后在NVE系综下采样(也就是把热浴关掉)。但对于有些体系由于某些原因( ...

好的,非常感谢了!
9楼2015-04-13 14:35:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chinaboy798

银虫 (小有名气)

【答案】应助回帖

好想学习下扩散的建模   以及in文件的设置     不知道大神们有联系方式么  ?简单说下即可1
10楼2015-04-24 15:58:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wj925806 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 04:45 by 4wMiSEwB6436
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7lpolszZVXgi 2026-08-14 5/250 2026-08-15 04:41 by 4wMiSEwB6436
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 4/200 2026-08-15 04:17 by 4wMiSEwB6436
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 7lpolszZVXgi 2026-08-14 4/200 2026-08-15 01:08 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 7/350 2026-08-15 00:40 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +23 医学老男孩 2026-08-13 49/2450 2026-08-14 20:02 by zhaifei
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见