24小时热门版块排行榜    

查看: 849  |  回复: 1
【悬赏金币】回答本帖问题,作者DOONWOON将赠送您 5 个金币

DOONWOON

新虫 (初入文坛)

[求助] 请大家帮帮我,我模拟的msd有问题

大家好,我模拟的msd为什么别文献差距很大呢,别人的一百万步和我的一百万步差了几十倍,in文件和data文件是一样的,我只是改了一下在nvt下面模拟,究竟是哪里出了问题呢,下面是我的in文件和我模拟出来的数据,别人都是1000左右,我跑一百万步才25
fix                             npt all npt temp $t $t 100 aniso 1 1 100
run                           200000
unfix                        npt

fix                             nvt all nvt temp $t $t 100
run                           200000
unfix                        nvt

reset_timestep       0

fix                             nvt all nvt temp $t $t 100


group                      na_group type 3
group                      k_group type 2
group                      co_group type 1 4
group                      mol molecule 1

compute                3 na_group msd
fix                            3 na_group ave/time   10 100 1000 c_3[4] c_3[1] c_3[2] c_3[3] file na_group.msd

compute                 4 co_group msd
fix                            4 co_group ave/time   10 100 1000 c_4[4] c_4[1] c_4[2] c_4[3] file co_group.msd

compute                 5 mol msd
fix                             5 mol ave/time   10 100 1000 c_5[4] c_5[1] c_5[2] c_5[3] file mol_group.msd

compute                6 k_group msd
fix                            6  k_group ave/time   10 100 1000 c_6[4] c_6[1] c_6[2] c_6[3] file k_group.msd

run                                  1000000

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

DOONWOON

新虫 (初入文坛)

2楼2021-07-14 18:06:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 DOONWOON 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 ufgqpyvpi6 2026-07-18 5/250 2026-07-19 20:39 by mxyes0ckns
[找工作] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 ufgqpyvpi6 2026-07-18 4/200 2026-07-19 20:23 by mxyes0ckns
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 ufgqpyvpi6 2026-07-18 5/250 2026-07-19 20:19 by mxyes0ckns
[基金申请] 5个%2F,非固定段22个字母 +11 jaasxulu 2026-07-15 12/600 2026-07-19 20:15 by 恋城1988
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 ufgqpyvpi6 2026-07-18 5/250 2026-07-19 20:10 by mxyes0ckns
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 ufgqpyvpi6 2026-07-18 7/350 2026-07-19 20:10 by mxyes0ckns
[基金申请] 面上的水很深 +17 布布和一二 2026-07-18 22/1100 2026-07-19 20:03 by botar
[硕博家园] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 ufgqpyvpi6 2026-07-18 5/250 2026-07-19 20:03 by mxyes0ckns
[基金申请] 小木虫没落了,除了祈祷帖子,几乎看不到有价值的帖子 +15 botar 2026-07-14 16/800 2026-07-19 18:23 by wowsunflower
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 ufgqpyvpi6 2026-07-18 6/300 2026-07-19 16:03 by 0u3sh5hhfi
[基金申请] 生命口面上会评结束了吗大家有知道结果吗?现在还没有结果正常不? +6 碧水00 2026-07-16 9/450 2026-07-17 15:21 by jaasxulu
[论文投稿] 跨出版社商投稿 20+5 倪好520 2026-07-13 5/250 2026-07-16 17:24 by 无视军团
[基金申请] 生生命学部会评 +4 wuke100666 2026-07-15 4/200 2026-07-16 00:12 by 啦啦小妮子
[基金申请] H口,面上,时间戳7月7日。 +3 plumage5 2026-07-15 3/150 2026-07-15 21:06 by 6543yes
[基金申请] 会评 +5 无名者登山 2026-07-15 5/250 2026-07-15 18:15 by 杠就是你对
[基金申请] E口会评 +11 布布和一二 2026-07-13 13/650 2026-07-15 16:34 by Kittylucky
[基金申请] E口青C会评结束了吗? 40+3 小小书虫c 2026-07-14 9/450 2026-07-15 16:04 by 安静的知了
[论文投稿] 农业机械学报投稿周期 +3 G0DLIN 2026-07-14 7/350 2026-07-15 09:26 by xs74101122
[基金申请] 国自然面上D口祈祷 +4 monkeynana 2026-07-14 4/200 2026-07-15 09:14 by 白水煮
[基金申请] 没收到消息,再次陪跑 +3 天空星辰fly 2026-07-14 3/150 2026-07-14 19:14 by 6543yes
信息提示
请填处理意见