| 查看: 1221 | 回复: 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 哪位大神指点一下！ |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
范进中举一文的中心思想
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
-
filecode，4个jtjc了
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- LAMMPS温度控制----如何达到平衡态
已经有4人回复
- LAMMPS温度－步长曲线震荡
已经有11人回复
- 【讨论】lammps-怪异的group温度
已经有9人回复
- 【讨论】lammps中NVE和NVT定温有什么区别
已经有7人回复
 |
2楼2013-12-28 10:09:57
40