24小时热门版块排行榜    

查看: 1573  |  回复: 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

木虫 (正式写手)

【答案】应助回帖

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

那就再仔细阅读文献吧。或者给作者发邮件。
4楼2016-03-21 11:42:22
已阅   回复此楼   关注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的回帖

mileyfan

新虫 (小有名气)

楼主请问你现在解决这个问题了么?
我也是经常出现个别水分子在非周期边界上越跑越远。很困扰
十分感谢!!!!
5楼2016-04-13 10:35:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +6 十三加油 2026-03-21 6/300 2026-03-22 17:00 by i_cooler
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[考研] 资源与环境 调剂申请(333分) +5 holy J 2026-03-21 5/250 2026-03-21 22:42 by Catalysis25
[考研] 化学调剂 +5 yzysaa 2026-03-21 5/250 2026-03-21 22:12 by peike
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 085700资源与环境308求调剂 +12 墨墨漠 2026-03-18 13/650 2026-03-21 01:42 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕,总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
信息提示
请填处理意见