24小时热门版块排行榜    

查看: 2169  |  回复: 4

gyy86_happy

新虫 (初入文坛)

[交流] nve系综下,不加热浴,温度可能持续升高吗? 已有3人参与

刚接触分子动力学,有很多地方不大明白。

我用lammps模拟三维空间的LJ 粒子,边界条件是p,p,p。粒子是均匀分布在空间的,不存在个别粒子间距离过近的问题。

模拟过程中没有对粒子加热浴,采用的是nve系综。

粒子的初始温度设置为0.1,模拟进行1000步后,温度升高到1.08,之后就一直在1.08上下波动。

但是老板认为这样的温度变化有问题:1. 温度从0.1升高到1.08,这样的温度变化太大; 2. 没有加热浴的情况下,温度会一直缓慢增加,不可能处于稳定状态!

可是我认为刚开始的时候温度升高是粒子突然间具有速度,它们之间发生相互碰撞,导致温度升高。一段时间后,体系应该趋于平衡,即温度,能量,压强只是在很小的范围内才会波动。

我现在很困惑,请各位大侠指教,不胜感激!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gyy86_happy

新虫 (初入文坛)

十分感谢!
2楼2014-06-05 13:52:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher


小木虫: 金币+0.5, 给个红包,谢谢回帖
初始势能太高,nve后转化成动能所以温度升高。
说明初始结构不太合理,离稳定结构很远。
3楼2014-06-06 15:42:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

madfarmer

铁虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼长正解,初始的结构决定最终温度。稳定态温度取决于初始结构,NVE是保持能量不变,input里面没有类似MOVE之类的命令加入新能量的话,最后温度是能稳定的。
4楼2014-06-23 17:37:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

453695608

银虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
能不能告诉我一下这么设置升温啊,比如我要每步上升1度
学习更多的知识,为了未来
5楼2015-10-09 20:47:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gyy86_happy 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:39 by 9i7PCSZyMzot
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 00:51 by 9i7PCSZyMzot
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:15 by 9i7PCSZyMzot
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:03 by 9i7PCSZyMzot
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 4pRyVGCVerAA 2026-09-19 4/200 2026-09-20 00:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:51 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:39 by 9i7PCSZyMzot
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-18 7/350 2026-09-19 21:03 by 9i7PCSZyMzot
[教师之家] 现代”学阀”该如何界定 +3 iaeyuan 2026-09-19 3/150 2026-09-19 20:52 by 没昵称呀
[基金申请] 国社科系统bug了,是不是要放榜了? +7 kynobel 2026-09-16 9/450 2026-09-19 14:52 by 接好运婷
[考博] 申博发邮件 +8 Lmengk 2026-09-15 9/450 2026-09-19 09:27 by 中间体超市
[有机交流] 各位大神,目前国内有哪些比较好用的逆合成软件? 20+3 zapen 2026-09-15 9/450 2026-09-18 18:17 by zapen
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
信息提示
请填处理意见