On the theory of the CO+OH reaction, including H and C kinetic isotope effects

Full-dimensional quantum mechanical calculations for the tunneling behavior of HOCO dissociation to H + CO2

The tunneling behavior during HOCO dissociation to H + CO₂ was investigated by full-dimensional quantum mechanical calculations based on an accurate global potential energy surface. The tunneling lifetimes for the low-lying 1500 vibrational states were calculated using the low-storage filter diagonalization method after a 1 million-step Chebyshev propagation. In the calculated energy range, the lifetimes of different vibrational states with similar energy are found to differ by 3-4 orders of magnitude, and the lower limit for these tunneling lifetimes is consistent with the experimental results reported by Continetti et al. For the given vibrational progressions, the lifetime of the vibrational states decreases with the increasing energy level, which is consistent with the results of 1D simulation calculations reported by Bowman, but the declining curve obviously fluctuates, and the declining slope is significantly different from that obtained by 1D simulation. Due to a difference in the effective barrier width, the mode-specific behavior of the tunneling effect is manifested in that the O-C-O' and H-O-C bend modes lead to the largest enhancement and an inhibitory effect on the tunneling process, respectively.

Simultaneous H/D and 13C/12C Anomalous Kinetic Isotope Effects during the Sonolysis of Water in the Presence of Carbon Monoxide

Splitting of water molecules driven by ultrasound plays a central role in sonochemistry. While studies of sonoluminescence revealed the formation of a plasma inside the cavitation bubble, much less is known about the contribution of plasma chemical processes to the sonochemical mechanisms. Herein, we report for the first time sonochemical processes in water saturated with pure CO. The presence of CO causes a large increase in the H/D kinetic isotope effect (KIE) to α_H = 14.6 ± 1.8 in a 10% H₂O/D₂O mixture under 20 kHz ultrasound. The anomalous H/D KIE is attributed to electron quantum tunneling in the plasma produced by cavitation. In addition, CO₂ formed simultaneously with hydrogen during the sonochemical process is enriched with the ¹³C isotope, which indicates a V-V pumping mechanism typical for non-equilibrium plasma. Both observed KIEs unambiguously point to the contribution of quantum effects in sonochemical mechanisms.

Nonstatistical Reaction Dynamics

Nonstatistical dynamics is important for many chemical reactions. The Rice-Ramsperger-Kassel-Marcus (RRKM) theory of unimolecular kinetics assumes a reactant molecule maintains a statistical microcanonical ensemble of vibrational states during its dissociation so that its unimolecular dynamics are time independent. Such dynamics results when the reactant's atomic motion is chaotic or irregular. Intrinsic non-RRKM dynamics occurs when part of the reactant's phase space consists of quasiperiodic/regular motion and a bottleneck exists, so that the unimolecular rate constant is time dependent. Nonrandom excitation of a molecule may result in short-time apparent non-RRKM dynamics. For rotational activation, the 2J + 1 K levels for a particular J may be highly mixed, making K an active degree of freedom, or K may be a good quantum number and an adiabatic degree of freedom. Nonstatistical dynamics is often important for bimolecular reactions and their intermediates and for product-energy partitioning of bimolecular and unimolecular reactions. Post–transition state dynamics is often highly complex and nonstatistical.

Combined Experimental-Theoretical Study of the OH + CO → H + CO2 Reaction Dynamics

A combined experimental-theoretical study is performed to advance our understanding of the dynamics of the prototypical tetra-atom, complex-forming reaction OH + CO → H + CO₂, which is also of great practical relevance in combustion, Earth's atmosphere, and, potentially, Mars's atmosphere and interstellar chemistry. New crossed molecular beam experiments with mass spectrometric detection are analyzed together with the results from previous experiments and compared with quasi-classical trajectory (QCT) calculations on a new, full-dimensional potential energy surface (PES). Comparisons between experiment and theory are carried out both in the center-of-mass and laboratory frames. Good agreement is found between experiment and theory, both for product angular and translational energy distributions, leading to the conclusion that the new PES is the most accurate at present in elucidating the dynamics of this fundamental reaction. Yet, small deviations between experiment and theory remain and are presumably attributable to the QCT treatment of the scattering dynamics.

Spectroscopy and dynamics of the HOCO radical: insights into the OH + CO → H + CO2 reaction

Photoelectron–photofragment coincidence experiments coupled with quantum chemistry and dynamics calculations have significantly enhanced our understanding of the reactive intermediate HOCO.

Complex organic molecules along the accretion flow in isolated and externally irradiated protoplanetary disks

The birth environment of the Sun will have influenced the physical and chemical structure of the pre-solar nebula, including the attainable chemical complexity reached in the disk, important for prebiotic chemistry. The formation and distribution of complex organic molecules (COMs) in a disk around a T Tauri star is investigated for two scenarios: (i) an isolated disk, and (ii) a disk irradiated externally by a nearby massive star. The chemistry is calculated along the accretion flow from the outer disk inwards using a comprehensive network which includes gas-phase reactions, gas-grain interactions, and thermal grain-surface chemistry. Two simulations are performed, one beginning with complex ices and one with simple ices only. For the isolated disk, COMs are transported without major chemical alteration into the inner disk where they thermally desorb into the gas reaching an abundance representative of the initial assumed ice abundance. For simple ices, COMs can efficiently form on grain surfaces under the conditions in the outer disk. Gas-phase COMs are released into the molecular layer via photodesorption. For the irradiated disk, complex ices are also transported inwards; however, they undergo thermal processing caused by the warmer conditions in the irradiated disk which tends to reduce their abundance along the accretion flow. For simple ices, grain-surface chemistry cannot efficiently synthesise COMs in the outer disk because the necessary grain-surface radicals, which tend to be particularly volatile, are not sufficiently abundant on the grain surfaces. Gas-phase COMs are formed in the inner region of the irradiated disk via gas-phase chemistry induced by the desorption of strongly bound molecules such as methanol; hence, the abundances are not representative of the initial molecular abundances injected into the outer disk. These results suggest that the composition of comets formed in isolated disks may differ from those formed in externally irradiated disks with the latter composed of more simple ices.

Dynamics on the HOCO potential energy surface studied by dissociative photodetachment of HOCO− and DOCO−

An experimental study of the dissociative photodetachment (DPD) dynamics of HOCO(-) and DOCO(-) at a photon energy of 3.21 eV has been carried out to probe the potential energy surface of the HOCO free radical and the dynamics of the OH+CO-->H+CO(2) reaction. These photoelectron-photofragment coincidence experiments allow the identification of photodetachment processes leading to the production of stable HOCO free radicals and both the H+CO(2) and OH+CO dissociation channels on the neutral surface. Isotopic substitution by deuterium in the parent ion is observed to reduce the product branching ratio for the D+CO(2) channel, consistent with tunneling playing a role in this dissociation pathway. Other isotope effects on the detailed partitioning of kinetic energy between photoelectrons and photofragments are also discussed. The results are compared to recent theoretical predictions of this DPD process, and evidence for the involvement of vibrationally excited HOCO(-) anions is discussed.

Methyl Vinyl Ketone+OH and Methacrolein+OH Oxidation Reactions: A Master Equation Analysis of the Pressure- and Temperature-Dependent Rate Constants

High-level electronic structure calculations and master equation analyses were carried out to obtain the pressure- and temperature-dependent rate constants of the methyl vinyl ketone+OH and methacrolein+OH reactions. The balance between the OH addition reactions at the high-pressure limit, the OH addition reactions in the fall-off region, and the pressure-independent hydrogen abstractions involved in these multiwell and multichannel systems, has been shown to be crucial to understand the pressure and temperature dependence of each global reaction. In particular, the fall-off region of the OH addition reactions contributes to the inverse temperature dependence of the rate constants in the Arrhenius plots, leading to pressure-dependent negative activation energies. The pressure dependence of the methyl vinyl ketone+OH reaction is clearly more important than in the case of the methacrolein+OH reaction owing to the weight of the hydrogen abstraction process in this second system. Comparison of the theoretical rate constants and the experimental measurements shows quite good agreement.

On the theory of the reaction rate of vibrationally excited CO molecules with OH radicals

The dependence of the rate of the reaction CO+OH-->H+CO2 on the CO-vibrational excitation is treated here theoretically. Both the Rice-Ramsperger-Kassel-Marcus (RRKM) rate constant kRRKM and a nonstatistical modification knon [W.-C. Chen and R. A. Marcus, J. Chem. Phys. 123, 094307 (2005).] are used in the analysis. The experimentally measured rate constant shows an apparent (large error bars) decrease with increasing CO-vibrational temperature Tv over the range of Tv's studied, 298-1800 K. Both kRRKM(Tv) and knon(Tv) show the same trend over the Tv-range studied, but the knon(Tv) vs Tv plot shows a larger effect. The various trends can be understood in simple terms. The calculated rate constant kv decreases with increasing CO vibrational quantum number v, on going from v=0 to v=1, by factors of 1.5 and 3 in the RRKM and nonstatistical calculations, respectively. It then increases when v is increased further. These results can be regarded as a prediction when v state-selected rate constants become available.