A detailed test study of barrier heights for the HO2 + H2O + O3 reaction with various forms of multireference perturbation theory

We report an ab initio multireference perturbation theory investigation of the HO(2) + H(2)O + O(3) reaction, with particular emphasis on the barrier heights for two possible reaction mechanisms: oxygen abstraction and hydrogen abstraction, which are identified by two distinct saddle points. These saddle points and the corresponding pre-reactive complexes were optimized at the CASSCF(11,11) level while the single point energies were calculated with three different MRPT2 theories: MRMP, CASPT2, and SC-NEVPT2. Special attention has been drawn on the "intruder state" problem and the effect of its corrections on the relative energies. The results were then compared with single reference coupled-cluster methods and also with our recently obtained Kohn-Sham density functional theory (KS-DFT) calculations [L. P. Viegas and A. J. C. Varandas, Chem. Phys., (2011)]. It is found that the relative energies of the pre-reactive complexes have a very good agreement while the MRPT2 classical barrier heights are considerably higher than the KS-DFT ones, with the SC-NEVPT2 calculations having the highest energies between the MRPT2 methods. Possible explanations have been given to account for these differences.