pair_coeff * * ffield.reax.021.CHONSSi C H S
fix 5 all reax/c/species 1 100 2000 species.out
The order in the data file is C, H, S element, and only the C, H, S element in the C/H/O/N/S/Si force field file is used, but the output in the species.out file is the molecule of C, H, O. What is wrong?