24小时热门版块排行榜    

查看: 1362  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 09:48 by G6APbkg8SA6w
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:38 by G6APbkg8SA6w
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 09:28 by G6APbkg8SA6w
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 05:45 by gy1nBQXYQJqL
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见