| 查看: 962 | 回复: 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 哪位大神指点一下！ |
回复此楼» 猜你喜欢
2025冷门绝学什么时候出结果
已经有5人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
LAMMPS温度控制----如何达到平衡态
已经有4人回复- LAMMPS温度－步长曲线震荡
已经有11人回复
- 【讨论】lammps-怪异的group温度
已经有9人回复
- 【讨论】lammps中NVE和NVT定温有什么区别
已经有7人回复
 |
2楼2013-12-28 10:09:57