Reactions of O(1D) with atmospheric molecules

Unraveling the interaction between singlet state atomic oxygen O(1D) and water: toward the formation of oxywater and hydrogen peroxide

We performed high-level quantum mechanical calculations to explore the interaction of atomic oxygen in the ground triplet state, O(3P), and the excited singlet state, O(1D), with water. We reported the potential energy curves for a few lowest electronic states when an atomic oxygen approaches the oxygen of a water molecule. Our results predict the formation of a singlet oxywater species as the product of O(1D) and H2O which lies about 149.33 kJ mol-1 below the total energy of a singlet oxygen atom and a water molecule. Our calculations predict that an O(3P) atom interacting with a water molecule forms a triplet oxywater complex with a shallow minimum on the triplet potential energy surfaces. We examined the transition of the singlet state oxywater species to hydrogen peroxide through the unimolecular reaction pathway, a (1,2)-hydrogen shift. We reported the structural properties, vibrational frequencies, and dipole moments of oxywater species, the transition state, and hydrogen peroxide. We also reported the energy barrier for the transition, and we provided an estimate for the respective reaction rate constant. In addition, we investigated the impact of solvents on the reaction pathway using an implicit solvation model of water. We predict that a singlet state oxywater species has a longer lifetime in a water environment than in the gas phase.

Photophysical oxidation of HCHO produces HO2 radicals

Formaldehyde, HCHO, is the highest-volume carbonyl in the atmosphere. It absorbs sunlight at wavelengths shorter than 330 nm and photolyses to form H and HCO radicals, which then react with O2 to form HO2. Here we show HCHO has an additional HO2 formation pathway. At photolysis energies below the energetic threshold for radical formation we directly detect HO2 at low pressures by cavity ring-down spectroscopy and indirectly detect HO2 at 1 bar by Fourier-transform infrared spectroscopy end-product analysis. Supported by electronic structure theory and master equation simulations, we attribute this HO2 to photophysical oxidation (PPO): photoexcited HCHO relaxes non-radiatively to the ground electronic state where the far-from-equilibrium, vibrationally activated HCHO molecules react with thermal O2. PPO is likely to be a general mechanism in tropospheric chemistry and, unlike photolysis, PPO will increase with increasing O2 pressure.

Quantum and quasi-classical dynamics of the C(3P) + O2(3Σ-g) → CO(1Σ+) + O(1D) reaction on its electronic ground state

The dynamics of the C(³P) + O₂(³Σ⁻_g) → CO(¹Σ⁺) + O(¹D) reaction on its electronic ground state is investigated by using time-dependent wave packet propagation (TDWP) and quasi-classical trajectory (QCT) simulations. For the moderate collision energies considered (E_c = 0.001 to 0.4 eV, corresponding to a range from 10 K to 4600 K) the total reaction probabilities from the two different treatments of the nuclear dynamics agree very favourably. The undulations present in P(E) from the quantum mechanical treatment can be related to stabilization of the intermediate CO₂ complex with lifetimes on the 0.05 ps time scale. This is also confirmed from direct analysis of the TDWP simulations and QCT trajectories. Product diatom vibrational and rotational level resolved state-to-state reaction probabilities from TDWP and QCT simulations agree well except for the highest product vibrational states (v′ ≥ 15) and for the lowest product rotational states (j′ ≤ 10). Opening of the product vibrational level CO(v′ = 17) requires ∼0.2 eV from QCT and TDWP simulations with O₂(j = 0) and decreases to 0.04 eV if all initial rotational states are included in the QCT analysis, compared with E_c > 0.04 eV obtained from experiments. It is thus concluded that QCT simulations are suitable for investigating and realistically describe the C(³P) + O₂(³Σ⁻_g) → CO(¹Σ⁺) + O(¹D) reaction down to low collision energies when compared with results from a quantum mechanical treatment using TDWPs.

The dynamics of the C(³P) + O₂(³Σ⁻_g) → CO(¹Σ⁺) + O(¹D) reaction on its electronic ground state is investigated by using time-dependent wave packet propagation (TDWP) and quasi-classical trajectory (QCT) simulations.

The C(3P) + O2(3Σg-) → CO2 ↔ CO(1Σ+) + O(1D)/O(3P) reaction: thermal and vibrational relaxation rates from 15 K to 20 000 K

Thermal rates for the C(3P) + O2(3Σg-) ↔ CO(1Σ+)+ O(1D)/O(3P) reaction are investigated over a wide temperature range based on quasi classical trajectory (QCT) simulations on 3-dimensional, reactive potential energy surfaces (PESs) for the 1A', (2)1A', 1A'', 3A' and 3A'' states. These five states are the energetically low-lying states of CO2 and their PESs are computed at the MRCISD+Q/aug-cc-pVTZ level of theory using a state-average CASSCF reference wave function. Analysis of the different electronic states for the CO2 → CO + O dissociation channel rationalizes the topography of this region of the PESs. The forward rates from QCT simulations match measurements between 15 K and 295 K whereas the equilibrium constant determined from the forward and reverse rates is consistent with that derived from statistical mechanics at high temperature. Vibrational relaxation, O + CO(ν = 1,2) → O + CO(ν = 0), is found to involve both, non-reactive and reactive processes. The contact time required for vibrational relaxation to take place is τ ≥ 150 fs for non-reacting and τ ≥ 330 fs for reacting (oxygen atom exchange) trajectories and the two processes are shown to probe different parts of the global potential energy surface. In agreement with experiments, low collision energy reactions for the C(3P) + O2(3Σg-, ν = 0) → CO(1Σ+) + O(1D) lead to CO(1Σ+, ν' = 17) with an onset at Ec ∼ 0.15 eV, dominated by the 1A' surface with contributions from the 3A' surface. Finally, the barrier for the COA(1Σ+) + OB(3P) → COB(1Σ+) + OA(3P) atom exchange reaction on the 3A' PES yields a barrier of ∼7 kcal mol-1 (0.300 eV), consistent with an experimentally reported value of 6.9 kcal mol-1 (0.299 eV).

Distribution of Internal States of CO from O (1D) + CO Determined with Time-Resolved Fourier Transform Spectroscopy

Following collisions of O (1D) with CO, rotationally resolved emission spectra of CO (1 < or = v < or = 6) in the spectral region 1800-2350 cm(-1) were detected with a step-scan Fourier transform spectrometer. O (1D) was produced by photolysis of O3 with light from a KrF excimer laser at 248 nm. Upon irradiation of a flowing mixture of O3 (0.016 Torr) and CO (0.058 Torr), emission of CO (v < or = 6) increases with time, reaches a maximum approximately 10 micros. At the earliest applicable period (2-3 micros), the rotational distribution of CO is not Boltzmann; it may be approximately described with a bimodal distribution corresponding to temperatures approximately 8000 and approximately 500 K, with the proportion of these two components varying with the vibrational level. A short extrapolation from data in the period 2-6 micros leads to a nascent rotational temperature of approximately 10170 +/- 600 K for v = 1 and approximately 1400 +/- 40 K for v = 6, with an average rotational energy of 33 +/- 6 kJ mol(-1). Absorption by CO (v = 0) in the system interfered with population of low J levels of CO (v = 1). The observed vibrational distribution of (v = 2):(v = 3):(v = 4):(v = 5):(v = 6) = 1.00:0.64:0.51:0.32:0.16 corresponds to a vibrational temperature of 6850 +/- 750 K. An average vibrational energy of 40 +/- 4 kJ mol(-1) is derived based on the observed population of CO (2 < or = v < or = 6) and estimates of the population of CO (v = 0, 1, and 7) by extrapolation. The observed rotational distributions of CO (1 < or = v < or = 3) are consistent with results of previous experiments and trajectory calculations; data for CO (4 < or = v < or = 6) are new.

Plasma Chemistry of NO in Complex Gas Mixtures Excited with a Surfatron Launcher

The plasma chemistry of NO has been investigated in gas mixtures with oxygen and/or hydrocarbon and Ar as carrier gas. Surface wave discharges operating at microwave frequencies have been used for this study. The different plasma reactions have been analyzed for a pressure range between 30 and 75 Torr. Differences in product concentration and/or reaction yields smaller than 10% were found as a function of this parameter. The following gas mixtures have been considered for investigation: Ar/NO, Ar/NO/O2, Ar/NO/CH4, Ar/CH4/O2, Ar/NO/CH4/O2. It is found that NO decomposes into N2 and O2, whereas other products such as CO, H2, and H2O are also formed when CH4 and O2 are present in the reaction mixture. Depending on the working conditions, other minority products such as HCN, CO2, and C2 or higher hydrocarbons have been also detected. The reaction of an Ar/NO plasma with deposits of solid carbon has also been studied. The experiments have provided useful information with respect to the possible removal of soot particles by this type of plasma. It has been shown that carbon deposits are progressively burned off by interaction with the plasma, and practically 100% decomposition of NO was found. Plasma intermediate species have been studied by optical emission spectroscopy (OES). Bands and/or peaks due to N2*, NO*, OH*, C2*, CN*, CH*, or H* were detected with different relative intensities depending on the gas mixture. From the analysis of both the reaction products and efficiency and the type of intermediate species detected by OES, different plasma reactions and processes are proposed to describe the plasma chemistry of NO in each particular mixture of gases. The results obtained provide interesting insights about the plasma removal of NO in real gas exhausts.