| 查看: 1703 | 回复: 4 | ||
[求助]
fix nvt 体系为何会越来越大？ 已有1人参与
|
|
做lammps模拟，大概是金属板上面一层水，取的是p p fm边界 用nvt跑平衡的时候，系统的体积会越来越大，总会有几个水分子跑得很远很远，这是正常情况吗？ 照理说nvt的话体积应该恒定啊~~ 我看的文献里并没有提到这种现象，是我的边界条件错误了还是漏了些什么东西呢？ 参考的文献是Martin Fitzner, Gabriele C. Sosso, Stephen J. Cox, et al. JACS.137,13658(2015) in文件： echo screen units real newton on dimension 3 boundary p p fm atom_style atomic region box block -30.00 30.00 -30.00 30.00 0.00 70.00 units box create_box 2 box lattice bcc 2.8665 region metel block INF INF INF INF 0.00 4.00 units lattice create_atoms 1 region metel lattice sc 3.15 region water block INF INF INF INF 13.50 48.50 units box create_atoms 2 region water mass 1 55.845 mass 2 18.01528 pair_style hybrid lj/cut 7.53 sw pair_coeff * * sw mW.sw NULL HOO pair_coeff 1 2 lj/cut 1.00 2.34 pair_coeff 1 1 none group M type 1 group W type 2 timestep 10 fix 1 M setforce 0.00 0.00 0.00 dump min all atom 100 min.lammpstrj minimize 1.0e-5 1.0e-6 10000 100000 undump min thermo 100 dump tmp all atom 10000 tmp1.lammpstrj fix 2 all nvt temp 290.0 290.0 500 tchain 10 run 1000000 undump tmp1 |
回复此楼» 猜你喜欢
关于代码变化问题，想知道的进来
已经有5人回复
-
关于filecode
已经有4人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
» 本主题相关商家推荐: (我也要在这里推广)
alundilong
木虫 (正式写手)
- 应助: 129 (高中生)
- 金币: 2019.1
- 红花: 18
- 帖子: 309
- 在线: 56小时
- 虫号: 595468
- 注册: 2008-09-08
- 性别: GG
- 专业: 传热传质学
2楼2016-03-18 23:51:43
3楼2016-03-21 09:18:27
alundilong
木虫 (正式写手)
- 应助: 129 (高中生)
- 金币: 2019.1
- 红花: 18
- 帖子: 309
- 在线: 56小时
- 虫号: 595468
- 注册: 2008-09-08
- 性别: GG
- 专业: 传热传质学
4楼2016-03-21 11:42:22
mileyfan
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 221.8
- 帖子: 53
- 在线: 42.4小时
- 虫号: 3333003
- 注册: 2014-07-21
- 专业: 气体放电与放电等离子体技
5楼2016-04-13 10:35:38
100