24小时热门版块排行榜    

查看: 596  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 牛很牛 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 不应该看fileCode +6 且听虎啸 2026-08-12 8/400 2026-08-13 12:07 by 小木虫6752397
[基金申请] FileCode能看出啥? +8 要乐观耀哥 2026-08-10 26/1300 2026-08-13 12:01 by Tide man
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 11/550 2026-08-13 10:40 by ZH远山
[基金申请] 重要来源:本周末出结果 +9 瞬息宇宙 2026-08-12 9/450 2026-08-13 10:12 by Vivilian
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
信息提示
请填处理意见