Eley–Rideal reaction dynamics between O atoms on β-cristobalite (1 0 0) surface: A new interpolated potential energy surface and classical trajectory study

We present a theoretical study of the collisions of atomic oxygen with O-precovered β-cristobalite (1 0 0) surface. We have constructed a multidimensional potential energy surface for the O₂/β-cristobalite (1 0 0) system based mainly on a dense grid of density functional theory points by using the interpolation corrugation-reducing procedure. Classical trajectories have been computed for quasithermal (100–1500 K) and state-specific (e.g., collision energies between 0.01 and 4 eV) conditions of reactants for different O incident angles (θ_v). Atomic sticking and O_2(adsorbed) formation are the main processes, although atomic reflection and Eley–Rideal (ER) reaction (i.e., O₂ gas) are also significant, depending their reaction probabilities on the O incident angle. ER reaction is enhanced by temperature increase, with an activation energy derived from the atomic recombination coefficient (γ_O(θ_v = 0°, T)) equal to 0.24 ± 0.02 eV within the 500–1500 K range, in close agreement with experimental data. Calculated γ_O(θ_v = 0°, T) values compare quite well with available experimental γ_O(T) although a more accurate calculation is proposed. Chemical energy accommodation coefficient β_O(T) is also discussed as a function of ER and other competitive contributions.