24小时热门版块排行榜    

查看: 3224  |  回复: 19

牛很牛

新虫 (小有名气)

[交流] 水分子模拟 已有5人参与

我现在在做水分子的模拟,但是现在设置初始温度为150,初始速度为400,但是运行不几步之后,水分子的温度就为零了,然后整个系统中的分子都是在非常缓慢的移动,几乎都不动。不知道问题所在,求大神指点

发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

王精英

新虫 (小有名气)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
也不知道用的什么软件

发自小木虫IOS客户端
DFT计算
2楼2019-03-20 09:53:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

arsennnic

金虫 (小有名气)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
几乎不动的话就说明速度为0了呀,当分子速度为0时统计出的温度也就是0了。建议你把代码贴出来看看。

发自小木虫IOS客户端
3楼2019-03-20 10:04:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛很牛

新虫 (小有名气)

引用回帖:
2楼: Originally posted by 王精英 at 2019-03-20 09:53:57
也不知道用的什么软件

lammps
4楼2019-03-20 10:06:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛很牛

新虫 (小有名气)

引用回帖:
3楼: Originally posted by arsennnic at 2019-03-20 10:04:38
几乎不动的话就说明速度为0了呀,当分子速度为0时统计出的温度也就是0了。建议你把代码贴出来看看。

好,我把温度这是这块截图发出来。稍等
5楼2019-03-20 10:06:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛很牛

新虫 (小有名气)

###########initial setting ###########
#fix    fixshake1 suan shake 0.0001 5 0 b 1 2 3 4 a 1 2 3  #这一行有问题,要问
fix 10 suan rigid molecule
velocity all  create 400 49284
timestep     1.0

###############  compute  #############
compute 1 Na  temp
compute 2 Cl  temp
compute 3 suan temp
compute 4 water temp
compute 5 all pe
compute ake all ke

compute energy all group/group suan

thermo_style  custom step temp c_1 c_2 c_3
thermo       1000
thermo_modify temp 1

fix      1  all      nvt temp 300.0 300.0 100.0
fix  rec all recenter INIT INIT INIT

dump     1   all custom   1000 ch0110.water id type x y z

run 1000000

write_data restart.water
6楼2019-03-20 10:08:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛很牛

新虫 (小有名气)

看是要沉下去了,不要啊
7楼2019-03-20 15:07:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)

8楼2019-03-20 22:00:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Nater1ver

木虫 (正式写手)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
你用的什么单位,timestep设成1.0,而且你的fix nvp的tdamp又是100.0,那频率怕不是1/100年

发自小木虫Android客户端
9楼2019-03-21 00:49:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛很牛

新虫 (小有名气)

引用回帖:
8楼: Originally posted by chuntao118 at 2019-03-20 22:00:03
自己写个程序呗!

这不就是我的程序

发自小木虫Android客户端
10楼2019-03-21 08:46:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 牛很牛 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 厉娜 +4 张丹枫1991 2026-11-09 6/300 2026-11-15 00:48 by Equinoxhua
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 YHiw1gou1Eda 2026-11-10 8/400 2026-11-14 22:20 by Equinoxhua
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-14 13:35 by KXypaGgBBlw7
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-13 23:37 by 6PrO1XyZ5URF
[考博] 考博求助 10+3 社会实践对别 2026-11-08 3/150 2026-11-13 14:04 by 沐小虫252
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 FChKX8meu1ds 2026-11-09 10/500 2026-11-13 11:08 by Equinoxhua
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 C6CRn66RYYXS 2026-11-10 6/300 2026-11-13 04:21 by Equinoxhua
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 13:02 by ECLTXnJnRSXf
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:48 by ECLTXnJnRSXf
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:11 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 8/400 2026-11-12 10:16 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 cqv1zKnzd7hr 2026-11-09 12/600 2026-11-12 09:32 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 08:15 by HajghKOMynmz
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 5/250 2026-11-12 08:15 by HajghKOMynmz
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cdqGD6ftHmbl 2026-11-11 3/150 2026-11-12 07:53 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 gRND4NVAZ8VG 2026-11-11 3/150 2026-11-12 07:42 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 05:30 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-12 05:08 by 9fOTqj074tIg
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 0TEMgxdT1eIw 2026-11-09 8/400 2026-11-12 05:00 by 9fOTqj074tIg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 17:03 by TgvkbZZl9tqB
信息提示
请填处理意见