reaxff力场来模拟聚酰亚胺的热裂解，弛豫过程控温失败，体系不平衡

wangqian123

我想利用reaxff力场来模拟聚酰亚胺的热裂解情况，但温度设置在300k的情况下也发生裂解，300k为常温，应该不会发生裂解，看运行日志发现，在弛豫过程中，温度一下子从我设置的5k跳到5000kin 文件部分为
timestep        0.1
velocity all create 5.0 321456 rot yes dist gaussian
thermo 100
thermo_style custom step temp pe ke etotal
dump   1 all atom 100 mini.lammpstrij
fix  1 all nve
unfix 1
undump 1
reset_timestep 0

fix  1 all nvt temp 100 100 100.0
fix     2 all qeq/reax 1 0.0 10.0 1e-6 reax/c
fix     4 all reax/c/species 1 100 100 species.out element C H O N
run 10000
unfix 1
unfix 2
reset_timestep 0

fix  1 all nve
fix  3 all temp/berendsen 300 300 100.0
fix     2 all qeq/reax 1 0.0 10.0 1e-6 reax/c
fix     4 all reax/c/species 1 100 100 species.out element C H O N
thermo 100
thermo_style custom step temp pe ke etotal3
timestep        0.1
dump        1 all custom 100 PE.xyz type x y z #每次输出升高1K
run 150000

弛豫过程日志部分为

wildon

data文件的mass部分发个截图看看

wangqian123

这样

wangqian123

wildon 发表于 2023-4-23 18:15
data文件的mass部分发个截图看看

有什么不对吗

cyx98

第一步是 fix  1 all nve。你要弛豫不就 fix 1 all nvt temp 300 300 10

Kmetsch

wangqian123 发表于 2023-4-23 19:00
这样

“pair coeff * * ffield.reax C H O N”
你的data中写的是NCOH

wildon

Kmetsch正解，in文件中的元素顺序要与data文件的元素顺序保持一致。