| 查看: 1731 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有10人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
为什么到现在没收到通知？
已经有3人回复
-
有没有仍没收到信息的
已经有5人回复
-
面上意见出来了
已经有10人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
投票:
30