求助：原子丢失

kidult

体系密度0.1，约3000个原子，在nvt系综300k下弛豫10ps，然后在nvt、1000k下弛豫10ps，接着把温度提升到3000k恒温燃烧。
出现”原子丢失“报错。试过很多方法（neighbor增大，减小时间步长）都不行
下面是我的in文件
==================== 基本模拟设置 ====================
units        real
atom_style   charge
boundary     p p p
read_data    C16H34.data
neighbor     3 bin
neigh_modify every 1 delay 0 check no
timestep     0.1

# ==================== 力场设置 ====================
pair_style      reaxff NULL      
pair_coeff      * * ffield.reax.cho C C H C C C

fix          02 all qeq/reax 1 0.0 10.0 1e-6 reaxff
# ==================== 能量最小化阶段 ====================
thermo       100
thermo_style custom step press pe temp ke
dump         1 all atom 100 minimize.xyz
minimize     1e-6 1e-6 5000 10000
undump       1
reset_timestep 0

# ==================== 第一阶段：300K NVT弛豫 (10ps) ====================
thermo       1000
thermo_style custom step press pe temp ke
dump         2 all atom 1000 relax_300K.xyz  # 修改输出文件名
fix          nvt all nvt temp 300 300 100
fix          2 all temp/berendsen 300 300 100
run          10000                      # 100ps
unfix        nvt
undump       2
reset_timestep 0
# ==================== 第一阶段：1000K NVT弛豫 (10ps) ====================
thermo       1000
thermo_style custom step press pe temp ke
dump         2 all atom 1000 relax_300K.xyz  # 修改输出文件名
velocity     all create 1500 4928459  
fix          nvt all nvt temp 1000 1000 100
fix          2 all temp/berendsen 1000 1000 100
run          100000                      # 10ps
unfix        nvt
undump       2
reset_timestep 0
# ==================== 燃烧模拟阶段 ====================
velocity all create 3000.0 4928459 rot yes dist gaussian
fix          nvt all nvt temp 3000 3000 50
fix          3 all temp/berendsen 3000 3000 50
fix          species all reaxff/species 50 100 5000 species2.log element C C H C C C
fix          reaxff_bonds all reaxff/bonds  10000 bonds.log
dump         4 all atom 10000 combustion.xyz
thermo       1000
thermo_style custom step temp press pe ke density
run          30000000

# ==================== 结果输出 ====================
write_data combustion_final.data
write_restart combustion.restart