A global ab initio potential energy surface for the X 2Aʹ ground state of the Si + OH → SiO + H reaction

We report the first global potential energy surface (PES) for the X(2)A' ground electronic state of the Si((3)P) + OH(X(2)Π) → SiO(X(1)Σg(+)) + H((2)S) reaction. The PES is based on a large number of ab initio energies obtained from multireference configuration interaction calculations plus Davidson correction (MRCI+Q) using basis sets of quadruple zeta quality. Corrections were applied to the ab initio energies in the reactant channel allowing a proper description of long-range interactions between Si((3)P) and OH(X(2)Π). An analytical representation of the global PES has been developed by means of the reproducing kernel Hilbert space method. The reaction is found barrierless. Two minima, corresponding to the SiOH and HSiO isomers, and six saddle points, among which the isomerization transition state, have been characterized on the PES. The vibrational spectra of the SiOH/HSiO radicals have been computed from second-order perturbation theory and quantum dynamics methods. The structural, energetic, and spectroscopic properties of the two isomers are in good agreement with experimental data and previous high quality calculations.