I want to make an oxidation reaction on FeCr alloy using ReaxFF potential, where should I start learning?
You have to find the ReaxFF potential file for Fe-Cr-O systems from the literature or to derive the ReaxFF parametrization on your own which is a quite advanced procedure.
Once you have the file, implementing the commands in LAMMPS is quite straightforward.
Hi, Thank you so much
Btw I have found the ReaxFF potential file for Fe-Cr-O systems, but I need O2 molecule file, where can I get that?
You write it yourself. See the documentation of the “molecule” command.
For ReaxFF you don’t even need to have a Bonds section and correspondingly the Special or Shake related sections.