I am trying to run simulation about the reaction of alkoxysilanes using ReaxFF.

I’m faced with problems during simulations.

When I’m performing ReaxFF MD using attached input script, the simulation stopped without any error messages. (There are 3328 atoms in the input structure.)

However, when I’m running the simulation with a 1/2 scale structure, this simulation proceeds normally.

I used lammps-10Aug15 version and 28 processes. (SMP cores)

How can I solve this problem?

​your LAMMPS version is over 2 years old. you should upgrade. there have been many bugfixes and improvements since. nobody will look into fixing problems without confirmation that the problem hasn't been fixed already.

also, LAMMPS doesn't just stop without any kind of error message. it may​ just be, that you didn't look in the right place.

also, your input deck is missing the force field file your are using, so nobody can reproduce the run.