24小时热门版块排行榜    

查看: 640  |  回复: 1

牛很牛

新虫 (小有名气)

[求助] 水分子温度趋近于零

目前我遇到的问题就是在模拟过程中,尽管我设置了分子的温度为一个定值,但是随着计算的进行,水分子的温度逐渐趋近于0。这就是目前我遇到的最尴尬,最头疼的问题,目前我自己查找参数的设置是没有问题,但一定是哪里出了问题,希望内行的大佬指点一二。稍后贴出代码
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛很牛

新虫 (小有名气)

dimension       3
boundary        p p p
units           real
#atom_style        full
neighbor        2.0 bin
neigh_modify  every 2 delay 2  
lattice fcc 1

region 1  block  -52.0 252.0 -52.0 252.0 -52.0 252.0   #all


pair_style  lj/cut/coul/long  15.0
kspace_style pppm 1.0e-5

atom_style        full
bond_style   harmonic
angle_style  harmonic
special_bonds lj 0.0 0.0 0.0 coul 0.0 0.0 0.0

read_data    new.data

这个是前面的设置
fix 10 suan rigid molecule
velocity all  create 400 49284
#velocity Cl  create 600 49284
#velocity suan create 600 49284
timestep     0.1

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

compute energy all group/group suan

thermo_style  custom step temp c_1 c_2 c_3
#press vol etotal c_3 c_1 c_4 c_2  c_ake c_gke c_energy
thermo       1000
thermo_modify temp 1

fix      1  all      nvt temp 450.0 450.0 100.0
#fix  rec all recenter INIT INIT INIT

dump     1   all  custom   1000 ch0110.water id type x y z
dump     2   NaCl custom   1000 ch0111.water id type x y z
dump     3   suan custom   1000 ch0112.water id type x y z


run 10000
#unfix 1

#fix       1    all   nvt temp 250.0 250.0 5.0
#fix       2    all    nve
#fix       3    jia    nve

#compute  cluster all cluster/atom  5.5
#dump     4   all custom   1000 ch0113.water id type x y z
#run 500000
write_data restart.water
这是后面的设置
2楼2019-04-08 15:30:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 牛很牛 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 WqtG0Aem8Ufl 2026-10-24 3/150 2026-10-25 20:24 by cLOojpOxkX0G
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 3/150 2026-10-25 20:11 by cLOojpOxkX0G
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 4/200 2026-10-25 20:11 by cLOojpOxkX0G
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 5/250 2026-10-25 20:07 by cLOojpOxkX0G
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 14:20 by XiYP04It3Evb
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 14:07 by XiYP04It3Evb
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 13:54 by XiYP04It3Evb
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-25 13:28 by XiYP04It3Evb
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-25 13:18 by XiYP04It3Evb
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-23 9/450 2026-10-25 13:02 by XiYP04It3Evb
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-25 12:59 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:51 by XiYP04It3Evb
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:49 by XiYP04It3Evb
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 7xHjyB3bAiwg 2026-10-23 4/200 2026-10-25 12:32 by XiYP04It3Evb
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GVGRDgNfv7P4 2026-10-23 9/450 2026-10-25 12:10 by XiYP04It3Evb
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 GVGRDgNfv7P4 2026-10-23 4/200 2026-10-25 12:10 by XiYP04It3Evb
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 11/550 2026-10-24 23:35 by XcOYG53d3x4p
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-24 22:56 by XcOYG53d3x4p
[文学芳草园] 亲密关系 +7 宋昙 2026-10-21 13/650 2026-10-23 17:07 by 小牧11
[考博] 郑州大学申博 +5 stars~ 2026-10-19 8/400 2026-10-22 15:06 by stars~
信息提示
请填处理意见