Reactivity trends of the MoVOx mixed metal oxide catalyst from density functional modeling

Abstract

We modeled the oxidative dehydrogenation of hydrocarbons at various oxygen sites in the basal plane of the M1 (orthorhombic) phase of the MoVO_x catalyst material, determining hydrogen adsorption energies by applying a hybrid DFT method to three-layer slab models. H adsorption is favored at bridging O centers over terminal MoO sites. Thereby, the unpaired electron, originating from the H atom, was calculated to localize preferentially at V centers, except if there are no V⁵⁺ centers near the H adsorption site; in the latter case reduction of a Mo⁶⁺ center is predicted. A reducing electron may transfer from a metal center in the surface layer to a sub-surface V center, with an associated stabilization of ∼10–90 kJ mol⁻¹. We also probed the initial step for hydrolyzing the MoVO_x surface. The adsorption of a water molecule is calculated to be more favorable, by 30 kJ mol⁻¹, at a surface Mo⁶⁺ center than at a surface V⁵⁺ center (by 22 kJ mol⁻¹ over a surface V⁴⁺ center). Remarkably, the trends in the reaction free energies of hydrogenation and hydrolysis are quite similar for various oxygen sites and the position of the newly created polaron in the surface model. The specific energetics, however, depends notably on the latter factors.