24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者DOONWOON将赠送您 5 个金币

DOONWOON

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 25.3
帖子: 12
在线: 13.9小时
虫号: 24672078
注册: 2020-12-01
专业: 热物性与热物理测试技术

[求助] 请大家帮帮我，我模拟的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客户端

回复此楼