24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 539  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 牛很牛 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学070300 求调剂 +21 哈哈哈^_^ 2026-04-12 21/1050 2026-04-13 21:36 by nxybio2007
[考研] 复试调剂 +17 积极向上; 2026-04-10 19/950 2026-04-13 19:53 by 闲人终南山
[考研] 人工智能320调剂08工类还有机会吗 +15 振—TZ 2026-04-10 16/800 2026-04-13 14:11 by 张zhihao
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +8 慕绝cc 2026-04-09 8/400 2026-04-13 14:08 by 张zhihao
[考研] 339求调剂 +4 hanwudada 2026-04-12 4/200 2026-04-13 12:03 by 蓝云思雨
[考研] 0854调剂 +10 长弓傲 2026-04-11 11/550 2026-04-13 10:38 by wp06
[考研] 求调剂 +9 吃口冰激凌 2026-04-07 9/450 2026-04-13 08:46 by 紫曦紫棋
[考研] 368化学求调剂 +14 wwwwabcde 2026-04-07 15/750 2026-04-13 08:36 by lhj2009
[考研] 296求调剂 +14 汪!?! 2026-04-10 16/800 2026-04-12 10:48 by zhouyuwinner
[考研] 291求调剂 +8 关忆北. 2026-04-11 8/400 2026-04-12 09:32 by 逆水乘风
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 材料与化工300求调剂 +39 肖开文 2026-04-09 43/2150 2026-04-12 01:30 by 秋豆菜芽
[考研] 求调剂 +18 MAX怅惘 2026-04-09 20/1000 2026-04-11 23:31 by zhen~
[考研] 302分求调剂 +9 凡语祈愿 2026-04-08 10/500 2026-04-10 23:26 by 314126402
[考研] 调剂 +12 卷卷卷心菜_ 2026-04-09 13/650 2026-04-10 22:36 by Ftglcn90
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考博] 博士自荐 +7 可可小胖 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[论文投稿] 求助文献原文 10+3 18500821399 2026-04-08 3/150 2026-04-09 16:56 by 北京莱茵润色
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +10 努力奋斗112 2026-04-07 10/500 2026-04-08 15:01 by screening
[考研] 325 调剂 +6 QQ小虾 2026-04-07 6/300 2026-04-07 15:17 by Ccclqqq
信息提示
请填处理意见