| 查看: 1001 | 回复: 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 哪位大神指点一下！ |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- LAMMPS温度控制----如何达到平衡态
已经有4人回复
- LAMMPS温度－步长曲线震荡
已经有11人回复
- 【讨论】lammps-怪异的group温度
已经有9人回复
- 【讨论】lammps中NVE和NVT定温有什么区别
已经有7人回复
 |
2楼2013-12-28 10:09:57