| 查看: 1224 | 回复: 1 | |||
[求助]
本人使用lammps软件进行反应模拟时，温度出现飞温现象
|
|
本人使用lammps软件reaxff进行反应模拟时，体系温度出现飞温现象，以下是我的in文件： # ReaxFF potential for MOF units real boundary p p p atom_style charge read_data data.txt mass 1 12.0107 mass 3 1.00794 mass 4 15.9994 mass 5 15.9994 mass 6 65.38 mass 7 15.9994 mass 8 1.00794 pair_style reax 6.0 1 1 1.0e-6 pair_coeff * * ffield.reax 1 1 2 3 3 4 3 2 compute reax all pair reax variable eb equal c_reax[1] variable ea equal c_reax[2] variable elp equal c_reax[3] variable emol equal c_reax[4] variable ev equal c_reax[5] variable epen equal c_reax[6] variable ecoa equal c_reax[7] variable ehb equal c_reax[8] variable et equal c_reax[9] variable eco equal c_reax[10] variable ew equal c_reax[11] variable ep equal c_reax[12] variable efi equal c_reax[13] variable eqeq equal c_reax[14] neighbor 2.5 bin neigh_modify every 10 delay 0 check no group 1 type 1 2 3 4 5 6 group 2 type 7 8 #minimize 1.0e-4 1.0e-6 100 1000 velocity all create 300.0 4928459 units box fix 1 all reax/bonds 1000 bonds.reax fix 2 all npt temp 300.0 300.0 10.0 iso 1.0 1.0 500.0 dump 1 all custom 1000 trj.1water id type q x y z #dump_modify 1 scale no dump_modify 1 sort 1 restart 100 restart1 restart2 thermo 1000 thermo_style custom step cpu temp evdwl ecoul epair etotal press vol enthalpy & v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa v_ehb & v_et v_eco v_ew v_ep v_efi v_eqeq thermo_modify flush yes timestep 1 run 2000000 哪位大神指点一下！ |
回复此楼» 猜你喜欢
朋友圈看到的
已经有9人回复
-
能否退出参与的面上项目解除限项
已经有13人回复
-
让我中一个面上吧！
已经有16人回复
-
放榜前的不淡定
已经有17人回复
-
科研孤儿太难了
已经有19人回复
-
明天应该可查了！？
已经有4人回复
-
什么时候开奖？
已经有11人回复
-
范进中举一文的中心思想
已经有6人回复
-
建议基金发布提前给出明确的时间点
已经有14人回复
-
2026国自然函评费到账
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- LAMMPS温度控制----如何达到平衡态
已经有4人回复
- LAMMPS温度－步长曲线震荡
已经有11人回复
- 【讨论】lammps-怪异的group温度
已经有9人回复
- 【讨论】lammps中NVE和NVT定温有什么区别
已经有7人回复
 |
2楼2013-12-28 10:09:57
