24小时热门版块排行榜    

查看: 520  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 牛很牛 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 293求调剂 +12 zjl的号 2026-03-16 17/850 2026-03-22 16:51 by i_cooler
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 280求调剂 +11 咕噜晓晓 2026-03-18 12/600 2026-03-21 22:40 by ACS Nano——
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 306求0703调剂一志愿华中师范 +5 纸鱼ly 2026-03-21 5/250 2026-03-21 17:11 by 学员8dgXkO
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
信息提示
请填处理意见