24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1588  |  回复: 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
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alundilong

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
体系变大,是因为上边界是m边界。
改成f的话,我猜你的原子会丢失。
2楼2016-03-18 23:51:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月下的淡然

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by alundilong at 2016-03-18 23:51:43
体系变大,是因为上边界是m边界。
改成f的话,我猜你的原子会丢失。

是的啊~这个知道的。。。怎么解决这类问题呢?文献中压根没提到。。。我再上边界加了个反射墙,不知道这样子会不会影响结果、
3楼2016-03-21 09:18:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alundilong

木虫 (正式写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 月下的淡然 at 2016-03-20 19:18:27
是的啊~这个知道的。。。怎么解决这类问题呢?文献中压根没提到。。。我再上边界加了个反射墙,不知道这样子会不会影响结果、...

那就再仔细阅读文献吧。或者给作者发邮件。
4楼2016-03-21 11:42:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mileyfan

新虫 (小有名气)

楼主请问你现在解决这个问题了么?
我也是经常出现个别水分子在非周期边界上越跑越远。很困扰
十分感谢!!!!
5楼2016-04-13 10:35:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 月下的淡然 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料考研调剂 +27 云木达达 2026-04-11 29/1450 2026-04-12 22:48 by kyle12138
[考研] 339求调剂 +3 hanwudada 2026-04-12 3/150 2026-04-12 21:43 by lbsjt
[考研] 本科南方医科大学 一志愿985 药学学硕284分 求调剂 +5 弱水听文 2026-04-09 5/250 2026-04-12 13:16 by shengxi123
[考研] 22专硕求调剂 +6 haoyun上岸 2026-04-11 8/400 2026-04-11 23:21 by labixiaoqiao
[考研] 346,工科求调剂 +3 moser233 2026-04-09 3/150 2026-04-11 10:04 by zhq0425
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 285求调剂 +9 AZMK 2026-04-07 11/550 2026-04-10 15:24 by AZMK
[考研] 机械还有还有名额吗?太难了 +6 笑笑袁 2026-04-10 6/300 2026-04-10 11:54 by 高维春
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 296求调剂 +6 汪!?! 2026-04-08 6/300 2026-04-10 11:02 by mattzhming
[考研] 0703化学求调剂 +21 不知名的小卅 2026-04-08 21/1050 2026-04-09 18:55 by l_paradox
[考研] 材料工程调剂 +12 小刘同学吖吖 2026-04-06 13/650 2026-04-09 17:07 by luoyongfeng
[考研] 调剂 +12 月@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
[考研] 334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 09:28 by wdyheheeh
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 求调剂 +13 柒luck 2026-04-07 13/650 2026-04-08 22:46 by 猪会飞
[考研] 277求调剂 +4 考研调剂lxh 2026-04-06 6/300 2026-04-08 10:40 by 逆水乘风
[考研] 一志愿南科大生物学297分,求调剂推荐 +8 Y-yyusx 2026-04-06 9/450 2026-04-07 19:38 by biomichael
[考研] 372分材料与化工(085600)英二数二求调剂 +4 蓝笺片 2026-04-06 4/200 2026-04-07 12:30 by dongzh2009
[考研] 材料调剂 +5 小刘同学吖吖 2026-04-06 5/250 2026-04-06 18:34 by sherry_1901
信息提示
请填处理意见