24小时热门版块排行榜    

查看: 1931  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

whu2015_8

新虫 (初入文坛)

[求助] lammps模拟时为什么原子群会发生集体移动?是模拟有问题么?有没有什么解决办法?已有1人参与

请求群里的各位大神惠赐指导!小弟使用lammps在nvt系综下做模拟,均是周期性边界,初始状态下原子群位于z方向的正中间(上下均有三倍于其长度的真空层),但是在模拟的过程中,发现原子群在z方向不断发生移动,并且能从z方向的一端移出去再从另一端移动进来!可是我的本意是想让这组原子在z方向的中间位置不动的,模拟时可以很顺利跑完的,但这样的移动是出乎我意料的,请问是不是我的模拟哪里出了问题呢?谢谢各位大神!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by whu2015_8 at 2018-02-27 19:11:42
已经找到解决方法了,就是需要在材料的中间某一层使用setforce进行固定,消除掉z方向的速度就可以了

如果确定固定部分原子不影响所研究的目的,此法也可。
5楼2018-02-28 05:06:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

无远不界

铁虫 (正式写手)

2楼2018-02-26 18:15:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
whu2015_8: 金币+10, 有帮助, 非常感谢您的建议!昨天一个老师告诉我在中间很薄的一层进行setforce也可以,您的方法我稍后尝试一下,非常感谢! 2018-02-27 19:14:36
尝试加上例如这个fix:
fix xxx all momentum 100 linear 0 0 1

参考:
http://lammps.sandia.gov/doc/fix_momentum.html
3楼2018-02-26 23:54:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whu2015_8

新虫 (初入文坛)

已经找到解决方法了,就是需要在材料的中间某一层使用setforce进行固定,消除掉z方向的速度就可以了
4楼2018-02-27 19:11:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见