24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1587  |  回复: 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的回帖

mileyfan

新虫 (小有名气)

楼主请问你现在解决这个问题了么?
我也是经常出现个别水分子在非周期边界上越跑越远。很困扰
十分感谢!!!!
5楼2016-04-13 10:35:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工学硕294分,求导师收留 +12 yzyzx 2026-04-12 12/600 2026-04-13 00:08 by solbeg
[考研] 考研求调剂 +5 子木呐 2026-04-12 5/250 2026-04-12 22:15 by paopaotu326
[考研] 0854调剂 +9 长弓傲 2026-04-12 11/550 2026-04-12 22:13 by paopaotu326
[考研] 化工调剂求导师收留!一志愿失利,踏实肯干,有植物提取科研经历 +20 yzyzx 2026-04-09 21/1050 2026-04-12 00:12 by 小小小小啦啦啦
[考研] 药学305求调剂 +8 玛卡巴卡boom 2026-04-10 8/400 2026-04-12 00:07 by zhouwenxian
[考研] 农学0904 312求调剂 +3 Say Never 2026-04-11 3/150 2026-04-11 17:22 by daydayup2005
[考研] 本人女孩 +7 吼吼, 2026-04-10 9/450 2026-04-11 14:45 by ACS Nano——
[考研] 0854调剂 +5 音像店听花鼓戏 2026-04-10 5/250 2026-04-11 10:49 by qingpingzhu
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +17 努力奋斗112 2026-04-06 20/1000 2026-04-11 00:31 by wangjihu
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 071000生物学调剂求助 +17 zzzzwww 2026-04-09 20/1000 2026-04-10 15:55 by 求调剂zz
[考研] 材料调剂 +5 hzhahg 2026-04-06 5/250 2026-04-10 10:10 by may_新宇
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考研] 085404,285分求调剂 +12 薇薇考研 2026-04-07 14/700 2026-04-09 23:10 by parmtree
[考研] 材料307分求大佬组收留 +17 Hll胡 2026-04-07 17/850 2026-04-09 10:53 by liuhuiying09
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 293分求调剂,外语为俄语 +7 加一一九 2026-04-07 10/500 2026-04-08 20:14 by yutian743
[考研] 085602调剂 初试总分335 +3 19123253302 2026-04-06 3/150 2026-04-07 18:00 by jp9609
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 考研调剂 +3 Wwwwwww哇 2026-04-06 3/150 2026-04-06 20:55 by lbsjt
信息提示
请填处理意见