| 查看: 486 | 回复: 1 | |||
牛很牛新虫 (小有名气)
|
[求助]
水分子温度趋近于零
|
| 目前我遇到的问题就是在模拟过程中，尽管我设置了分子的温度为一个定值，但是随着计算的进行，水分子的温度逐渐趋近于0。这就是目前我遇到的最尴尬，最头疼的问题，目前我自己查找参数的设置是没有问题，但一定是哪里出了问题，希望内行的大佬指点一二。稍后贴出代码 |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
牛很牛
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 265.7
- 散金: 10
- 沙发: 1
- 帖子: 158
- 在线: 21.6小时
- 虫号: 8566555
- 注册: 2018-04-19
- 性别: GG
- 专业: 工程热物理相关交叉领域
|
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