Crossed-Beam and Theoretical Studies of the O(P, D) + Benzene Reactions: Primary Products, Branching Fractions, and Role of Intersystem Crossing.

Reliable modeling of hydrocarbon oxidation relies critically on knowledge of the branching fractions (BFs) as a function of temperature () and pressure () for the products of the reaction of the hydrocarbon with atomic oxygen in its ground state, O(P). During the past decade, we have performed in-depth investigations of the reactions of O(P) with a variety of small unsaturated hydrocarbons using the crossed molecular beam (CMB) technique with mass spectrometric (MS) detection and time-of-flight (TOF) analysis, combined with synergistic theoretical calculations of the relevant potential energy surfaces (PESs) and statistical computations of product BFs, including intersystem crossing (ISC). This has allowed us to determine the primary products, their BFs, and extent of ISC to ultimately provide theoretical channel-specific rate constants as a function of and . In this work, we have extended this approach to the oxidation of one of the most important species involved in the combustion of aromatics: the benzene (CH) molecule. Despite extensive experimental and theoretical studies on the kinetics and dynamics of the O(P) + CH reaction, the relative importance of the CHO (phenoxy) + H open-shell products and of the spin-forbidden CH (cyclopentadiene) + CO and phenol adduct closed-shell products are still open issues, which have hampered the development of reliable benzene combustion models. With the CMB technique, we have investigated the reaction dynamics of O(P) + benzene at a collision energy () of 8.2 kcal/mol, focusing on the occurrence of the phenoxy + H and spin-forbidden CH + CO and phenol channels in order to shed further light on the dynamics of this complex and important reaction, including the role of ISC. Concurrently, we have also investigated the reaction dynamics of O(D) + benzene at the same . Synergistic high-level electronic structure calculations of the underlying triplet/singlet PESs, including nonadiabatic couplings, have been performed to complement and assist the interpretation of the experimental results. Statistical (RRKM)/master equation (ME) computations of the product distribution and BFs on these PESs, with inclusion of ISC, have been performed and compared to experiment. In light of the reasonable agreement between the CMB experiment, literature kinetic experimental results, and theoretical predictions for the O(P) + benzene reaction, the so-validated computational methodology has been used to predict (i) the BF between the CHO + H and CH + CO channels as a function of collision energy and temperature (at 0.1 and 1 bar), showing that their increase progressively favors radical (phenoxy + H)-forming over molecule (CH + CO and phenol stabilization)-forming channels, and (ii) channel-specific rate constants as a function of and , which are expected to be useful for improved combustion models.