In numerous studies,
the application of the molecular dynamics
scheme based on the reactive force field (ReaxFF) method has been
proven effective in modeling the catalytic behavior of metal–organic
compounds. Recently, this method has been successfully applied for
M
x
O
y
(M =
Cu, Fe, Mn, Ni) transition-metal oxides. Yet, bimetallic metal oxides
of the type MnMO
x
(M = Cu, Fe, Ni) were
also present in the experimental system but could not be modeled since
not all of the force field parameters were available at the time.
To bridge this gap, the force field for modeling bimetallic metal
oxides had to be developed. Here, we establish the needed force field
parameter sets (namely, Cu/Mn/O, Fe/Mn/O, and Ni/Mn/O) and apply them
to the problem of toluene adsorption on bimetallic oxide catalyst
surfaces to verify their validity. Each training set consisted of
at least 10 crystal structures containing at least Cu–Mn–O,
Fe–Mn–O, or Ni–Mn–O atoms in contact obtained
from the available structure databases. The parameter training has
been done using the in-home-compiled version of the ReaxFF code. After
training the force fields for geometry reproduction, the parameters
were refined using the optimization by atom charges, comparing the
ReaxFF values to those obtained for the respective structures using
periodic crystal density functional theory (DFT) codes. The as-developed
force fields were then applied to the process of toluene adsorption/degradation
on MnMO
x
catalysts. Results obtained show
agreement with previous experimental expectations, although some remarks
are given since the initially presumed crystal structure of bimetallic
oxide Mn
1–
x
M
x
O
y
crystallites may still have
an impact on theoretical predictions. The presented are, to the best
of the authors’ knowledge, the first applications of the ReaxFF
approach to the Mn–(Cu|Fe|Ni)–O–C–H interaction.