Reactivity of metal oxide clusters with hydrogen peroxide and water - a DFT study evaluating the performance of different exchange-correlation functionals

We have performed a density functional theory (DFT) investigation of the interactions of H2O2, H2O and HO radicals with clusters of ZrO2, TiO2 and Y2O3. Different modes of H2O adsorption onto the clusters were studied. In almost all the cases the dissociative adsorption is more exothermic than molecular adsorption. At the surfaces where H2O has undergone dissociative adsorption, the adsorption of H2O2 and the transition state for its decomposition are mediated by hydrogen bonding with the surface HO groups. Using the functionals B3LYP, B3LYP-D and M06 with clusters of 26 and 8 units of ZrO2, the M06 functional performed better than B3LYP in describing the reaction of decomposition of H2O2 and the adsorption of H2O. Additionally, we investigated clusters of the type (ZrO2)(2), (TiO2)(2) and (Y2O3) and the performance of the functionals B3LYP, B3LYP-D, B3LYP*, M06, M06-L, PBE0, PBE and PWPW91 in describing H2O2, H2O and HO center dot adsorption and the energy barrier for decomposition of H2O2. The trends obtained for HO center dot adsorption onto the clusters are discussed in terms of the ionization energy of the metal cation present in the oxide. In order to correctly account for the existence of an energy barrier for the decomposition of H2O2, the functional used must include Hartree-Fock exchange. Using minimal cluster models, the best performance in describing the energy barrier for H2O2 decomposition was obtained with the M06 and PBE0 functionals - the average absolute deviations from experiments are 6 kJ mol(-1) and 5 kJ mol(-1) respectively. With the M06 functional and a larger monoclinic (ZrO2)(8) cluster model, the performance is in excellent agreement with experimental data. For the different oxides, PBE0 was found to be the most effective functional in terms of performance and computational time cost.