24小时热门版块排行榜    

查看: 1384  |  回复: 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-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:39 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 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文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:03 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.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
[基金申请] 国社科系统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 水士口
信息提示
请填处理意见