Kinetics of ozone destruction due to the recombination of oxygen atoms produced by pulsed 266 nm laser photolysis of O3/M (M = CO2 and/or N2) mixtures was studied using the absorption and emission spectroscopy to follow time evolutions of O3 and electronically excited molecules O2* formed in the recombination process 2O(3P) + M → O2* + M. An unexpected high ozone destruction rate was observed when O2* was present in the system. The kinetic model developed for the oxygen nightglow on the terrestrial planets was adapted to interpret the detected temporal profiles of the ozone number density and the O2* emission intensities. It was deduced that the vibrationally excited singlet delta oxygen molecule O2(a1Δ, υ) formed in the secondary processes reacts efficiently with ozone in the process O2(a1Δ, υ ≥ 3) + O3 → 2O2 + O, and the rate constant of this process was estimated to be 3 × 10−11 cm3 s−1. Ab initio calculations at the CASPT2(14, 12)/cc-pVTZ/UωB97XD/cc-pVTZ level of theory were applied to find the reaction pathway from the reactants to products on the O5 potential energy surface. These calculations revealed that the O2(a1Δ) + O3 reaction is likely to proceed via singlet–triplet intersystem crossing exhibiting an energy barrier of 9.6 kcal/mol, which lies between two and three quanta of vibrational excitation of O2(a1Δ), and hence, O2(a1Δ, υ) with υ ≥ 3 could rapidly react with ozone.