Physics

JOURNAL OF CHEMICAL PHYSICS

VOLUME 119, NUMBER 8

22 AUGUST 2003

Ab initio calculation of the ground „ 1 A … potential energy surface and theoretical rate constant for the Si¿O2\SiO¿O reaction

F. Dayou and A. Spielfiedel

` Laboratoire d’Etude du Rayonnement et de la Matiere en Astrophysique, UMR 8112 du CNRS, Observatoire de Paris-Meudon, 5 Place Jules Janssen, 92195 Meudon Cedex, France

Received 21 February 2003; accepted 2 June 2003 The entrance channel of the Si O2 →SiO O reaction has been investigated for collinear and perpendicular approach of the silicon atom to the O2 molecule by ab initio electronic structure calculations using the multireference configuration interaction MRCI method and Davidson correction MRCI Q . Results show that the reaction can proceed through the ground singlet ( 1 A ) and first triplet ( 3 A ) electronic states at low temperatures. The ground 1 A three-dimensional potential energy surface PES which correlates the Si( 3 P) O2 (X 3 g ) reactants to the SiO(X 1 ) O( 1 D) products was computed at the MRCI Q level of theory using the Woon and Dunning cc-pVTZ basis sets. The reaction was found barrierless and three minima have been PES with energy ordering: linear OSiO( 1 g ) triangular characterized on the 1 A 1 1 OSiO( A 1 ) linear SiOO( ). About 2500 ab initio data points have been fitted to a many body expansion using the method of Aguado and Paniagua, with a global root-mean-square of 1.49 kcal/mol. The analytical 1 A PES has been used to determine the thermal rate constants in the temperature range 15–300 K by quasiclassical trajectory calculations. Comparison with experimental results shows a quite good agreement for temperature dependence of the rate constants when the spin–orbit structure of the reactants is taken into account. The rate constants are also compared with earlier results of adiabatic capture calculations. The excellent agreement between both theoretical results for temperatures above 50 K points out an...