Development of a ReaxFF Reactive Force Field for Fe/Cr/O/S and Application to Oxidation of Butane over a Pyrite-Covered Cr₂O₃ Catalyst

Abstract

We developed a ReaxFF force field for Fe/Cr/O/S, which is parametrized against data from quantum mechanical (QM) calculations. Using this force field, we studied the Cr-oxide catalyzed oxidation reaction of butane at 1600 K. Our simulation results demonstrate that the active oxygen species on the oxide surface play an important role in the conversion of butane. Dehydrogenation of butane, which is found to be catalyzed by oxygen species on the oxide surface, initiates the reaction and generates butane radicals and surface OH groups. The radical intermediates are associated with the oxygen atoms to form C–O bonds or make double bonds when neighboring carbon atoms are dehydrogenated, forming light alkenes. On the clean Cr-oxide, the major oxidation product is CH2O. The presence of iron pyrite (FeS2), a common inorganic component in coal-derived fuels and a major slagging component, on Cr-oxide accelerates the complete oxidation of butane forming CO2 and CO. Surface reconstruction by iron pyrite is probably responsible for the change of the catalytic behavior. Reoxidation of the reduced oxide surface can occur through removal of surface H2O and adsorption of gaseous molecular oxygen at the vacancy sites on the clean Cr-oxide. On the other hand, on the modified Cr-oxide, it is found that a considerable amount of SOH molecules are released from the surface. These results can provide the detailed mechanisms for the catalytic oxidation of alkane and product distributions in Cr-oxide catalyst and give, for the first time, atomistic-scale insight in the complex surface chemistry of these catalysts under realistic operating conditions.

Statistical Analysis of Tri-Cresyl Phosphate Conversion on an Iron Oxide Surface Using Reactive Molecular Dynamics Simulations

“…20 Parameters for Fe/O/C/H interactions were taken from Ref. 21 (which used Fe/C parameters from 22 ), and P/O/C/H interactions were taken from Ref. 23 However, to describe the Fe-P interactions, it was necessary to develop a new parameter set.…”

Early and transient stages of Cu oxidation: Atomistic insights from theoretical simulations and in situ experiments

“…Theoretical studies of alloy system is either limited in the small scale at the framework of DFT calculations [146], or at the level of grain boundaries with less atomistic information [147,148]. Force field simulations of alloy models is still at its early stage [149]. It is crucial that the computational methods correlate closely with experimental techniques in addressing the reaction details at different stages and scales.…”

Benchmarking the Performance of the ReaxFF Reactive Force Field on Hydrogen Combustion Systems

“…ReaxFF potentials have been developed and applied to combustion reactions, [51][52][53][54][55][56][57][58][59][60][61][62][63][64][65][66][67] including hydrocarbons, hydrogen, 54,63 and syngas. 67 The ReaxFF method has also been extensively applied to material science, [68][69][70][71][72][73][74][75][76][77][78][79][80][81][82][83] catalysis, [84][85][86][87][88] and other chemical systems. 39,[89][90][91][92][93][94][95][96]…”