| 查看: 598 | 回复: 1 | |||
牛很牛新虫 (小有名气)
|
[求助]
水分子温度趋近于零
|
| 目前我遇到的问题就是在模拟过程中，尽管我设置了分子的温度为一个定值，但是随着计算的进行，水分子的温度逐渐趋近于0。这就是目前我遇到的最尴尬，最头疼的问题，目前我自己查找参数的设置是没有问题，但一定是哪里出了问题，希望内行的大佬指点一二。稍后贴出代码 |
回复此楼» 猜你喜欢
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有11人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有5人回复
-
重要来源：本周末出结果
已经有10人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有11人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有6人回复
-
不应该看fileCode
已经有9人回复
-
Filecode 又变了，巨变
已经有4人回复
-
分享一下我之前已中青C的计划书的filecode
已经有5人回复
-
FileCode能看出啥？
已经有26人回复
-
关于Filecode分析方法
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
牛很牛
新虫 (小有名气)
- 应助: 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
