| 查看: 1986 | 回复: 5 | ||||
[求助]
lammps模拟时为什么原子群会发生集体移动？是模拟有问题么？有没有什么解决办法？ 已有1人参与
|
| 请求群里的各位大神惠赐指导！小弟使用lammps在nvt系综下做模拟，均是周期性边界，初始状态下原子群位于z方向的正中间（上下均有三倍于其长度的真空层），但是在模拟的过程中，发现原子群在z方向不断发生移动，并且能从z方向的一端移出去再从另一端移动进来！可是我的本意是想让这组原子在z方向的中间位置不动的，模拟时可以很顺利跑完的，但这样的移动是出乎我意料的，请问是不是我的模拟哪里出了问题呢？谢谢各位大神！ |
回复此楼
无远不界
铁虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2651.6
- 红花: 6
- 帖子: 481
- 在线: 347.5小时
- 虫号: 2102396
- 注册: 2012-11-02
- 性别: GG
- 专业: 交通工程
2楼2018-02-26 18:15:17
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
whu2015_8: 金币+10, ★有帮助, 非常感谢您的建议！昨天一个老师告诉我在中间很薄的一层进行setforce也可以，您的方法我稍后尝试一下，非常感谢！ 2018-02-27 19:14:36
感谢参与，应助指数 +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
4楼2018-02-27 19:11:42
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
5楼2018-02-28 05:06:28
6楼2018-05-04 09:16:34
