24小时热门版块排行榜    

查看: 1440  |  回复: 3

leoluode

银虫 (小有名气)

[求助] 初始速度设置编程请教 已有1人参与

分子动力学新手,向各位求教。
想请问下,本人采用Lammps进行金属拉伸性能的分子动力学模拟,看文献说设置初始速度符合初始温度下的Maxwell-Boltzman分布,那么请问在程序中是采用哪个命令设置符合这样分布的初始速度呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

未曾彷徨

新虫 (初入文坛)

velocity create dist gaussian

发自小木虫IOS客户端

» 本帖已获得的红花(最新10朵)

2楼2015-12-15 16:55:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leoluode

银虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by 未曾彷徨 at 2015-12-15 16:55:06
velocity create dist gaussian

谢谢您的回答!
Lammps新手,我还是有些疑问,望解答
velocity create dist gaussian命令,
(1)以上分布代表速度符合高斯分布,而非Maxwell Boltzman分布。查阅了相关文献,发现速度在X,Y,Z轴方向的分量均符合高斯分布(正态分布),而速度(即三个分量的矢量和)是符合Maxwell Boltzman分布的。那么有疑问的是,这样设置的高斯分布是否与所要求的Maxwell Boltzman分布等同呢。
(2)例如velocity all create 400.0 4928459 dist gaussian,词句意思是所有原子在400K时的速度均符合高斯分布(或Maxwell),对于中间seed取值为4928459是如何来的,是随便取值的吗,有取值要求么?
期待您的热心解答
3楼2015-12-15 21:52:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiaqz

新虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
3楼: Originally posted by leoluode at 2015-12-15 08:52:05
谢谢您的回答!
Lammps新手,我还是有些疑问,望解答
velocity create dist gaussian命令,
(1)以上分布代表速度符合高斯分布,而非Maxwell Boltzman分布。查阅了相关文献,发现速度在X,Y,Z轴方向的分量均符合 ...

(1)是的
(2)是的,seed取正整数即可。
4楼2015-12-18 04:21:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leoluode 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 6/300 2026-11-12 02:25 by 9fOTqj074tIg
[基金申请] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 02:14 by 9fOTqj074tIg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 4/200 2026-11-12 02:12 by 9fOTqj074tIg
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 02:12 by 9fOTqj074tIg
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 5/250 2026-11-12 02:09 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 02:01 by 9fOTqj074tIg
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 01:50 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 5/250 2026-11-12 01:39 by 9fOTqj074tIg
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 C6CRn66RYYXS 2026-11-10 3/150 2026-11-12 01:28 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 C6CRn66RYYXS 2026-11-10 5/250 2026-11-12 01:25 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 00:55 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-12 00:44 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 4/200 2026-11-11 18:20 by TgvkbZZl9tqB
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-11 17:25 by TgvkbZZl9tqB
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 16:47 by TgvkbZZl9tqB
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 3/150 2026-11-11 11:41 by gRND4NVAZ8VG
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
信息提示
请填处理意见