Chapter 4 The Kinetics of Pressure-Dependent Reactions

Kinetic study of isoprene hydroxy hydroperoxide radicals reacting with sulphur dioxide and their global-scale impact on sulphate formation

Isoprene is the most relevant volatile organic compound emitted during the biosynthesis of metabolism processes. The oxidation of isoprene by a hydroxy radical (OH) is one of the main consumption schemes that generate six isomers of isoprene hydroxy hydroperoxide radicals (ISOPOOs). In this study, the rate constants of ISOPOOs + sulphur dioxide (SO2) reactions that eventually generate sulphur trioxide (SO3), the precursor of sulphate aerosol (SO42-(p)), are determined using microcanonical kinetic theories coupled with molecular structures and energies estimated by quantum chemical calculations. The results show that the reaction rates range from 10-27 to 10-20 cm3 molecule-1 s-1, depending on the atmospheric temperature and structure of the six ISOPOO isomers. The effect of SO3 formation from SO2 oxidation by ISOPOOs on the atmosphere is evaluated by a global chemical transport model, along with the rate constants obtained from microcanonical kinetic theories. The results show that SO3 formation is enhanced in regions with high SO2 or low nitrogen oxide (NO), such as China, the Middle East, and Amazon rainforests. However, the production rates of SO3 formation by ISOPOOs + SO2 reactions are eight orders of magnitude lower than that from the OH + SO2 reaction. This is indicative of SO42-(p) formation from the direct oxidation of SO2 by ISOPOOs, which is almost negligible in the atmosphere. The results of this study entail a detailed analysis of SO3 formation from gas-phase reactions of isoprene-derived products.

Theoretical Study on Ethylamine Dissociation Reactions Using VRC-VTST and SS-QRRK Methods

Barrierless bond dissociation reactions play an important role in fuel combustion. In this work, the pressure-dependent dissociation rate constants of ethylamine (EA) are accurately determined using variable-reaction-coordinate variational transition-state theory combined with the system-specific quantum Rice-Ramsperger-Kassel method. Before the kinetics calculations, the performances of four density functional theory methods in describing the bond dissociation of EA are evaluated against the benchmark method, FIC-MRCISD(T)+Q/cc-pVTZ, and the MN15-L/cc-pVTZ method is the best choice. By comparison of the Gibbs free energies and the rate constants for the bond dissociation reactions of EA, ethanol, and propane, the influence of functional groups on the reaction kinetics is discussed. The kinetics calculations show that the dissociation rate constants of EA are sensitive to pressure at low pressures and high temperatures, and the dominant channel is the reaction that yields C2H5 and NH2 radicals. A literature combustion model of EA is updated with our calculations, and the satisfactory agreement between the model predictions and reported ignition delay times of EA suggests the reliability of our calculations.

Mechanistic analysis by NMR spectroscopy: A users guide

A 'principles and practice' tutorial-style review of the application of solution-phase NMR in the analysis of the mechanisms of homogeneous organic and organometallic reactions and processes. This review of 345 references summarises why solution-phase NMR spectroscopy is uniquely effective in such studies, allowing non-destructive, quantitative analysis of a wide range of nuclei common to organic and organometallic reactions, providing exquisite structural detail, and using instrumentation that is routinely available in most chemistry research facilities. The review is in two parts. The first comprises an introduction to general techniques and equipment, and guidelines for their selection and application. Topics include practical aspects of the reaction itself, reaction monitoring techniques, NMR data acquisition and processing, analysis of temporal concentration data, NMR titrations, DOSY, and the use of isotopes. The second part comprises a series of 15 Case Studies, each selected to illustrate specific techniques and approaches discussed in the first part, including in situ NMR (^1/2H, ^10/11B, ¹³C, ¹⁵N, ¹⁹F, ²⁹Si, ³¹P), kinetic and equilibrium isotope effects, isotope entrainment, isotope shifts, isotopes at natural abundance, scalar coupling, kinetic analysis (VTNA, RPKA, simulation, steady-state), stopped-flow NMR, flow NMR, rapid injection NMR, pure shift NMR, dynamic nuclear polarisation, ¹H/¹⁹F DOSY NMR, and in situ illumination NMR.

Molecular Collision Diameters and Electronic Polarizabilities: Inherent Relationship and Fast Evaluation

The collision diameter σ for a large set of molecular species is related to the static electronic polarizability α_el. A remarkable correlation between these quantities conceptually similar to the analogous one previously identified for atoms is revealed. Our recommended model is the function σ(α_el) = p₁ + p₂α_el^1/3, where p₁ = 0.768 Å and p₂ = 2.168 are the fitting parameters providing the best overall match to the reference data for collision diameter (181 data points). The obtained correlation allows one to easily find the collision diameter of molecules from the known polarizability and vice versa. These findings can be useful for many applications, where there is a need for inexpensive assessments of the collision diameters or electronic polarizabilities, for example, when developing the transport property databases for modeling of chemically reacting flows.

CRAHCN-O: A Consistent Reduced Atmospheric Hybrid Chemical Network Oxygen Extension for Hydrogen Cyanide and Formaldehyde Chemistry in CO2-, N2-, H2O-, CH4-, and H2-Dominated Atmospheres

Hydrogen cyanide (HCN) and formaldehyde (H₂CO) are key precursors to biomolecules such as nucleobases and amino acids in planetary atmospheres. However, many reactions which produce and destroy these species in atmospheres containing CO₂ and H₂O are still missing from the literature. We use a quantum chemistry approach to find these missing reactions and calculate their rate coefficients using canonical variational transition state theory and Rice-Ramsperger-Kassel-Marcus/master equation theory at the BHandHLYP/aug-cc-pVDZ level of theory. We calculate the rate coefficients for 126 total reactions and validate our calculations by comparing with experimental data in the 39% of available cases. Our calculated rate coefficients are most frequently within a factor of 2 of experimental values and generally always within an order of magnitude of these values. We discover 45 previously unknown reactions and identify 6 from this list that are most likely to dominate H₂CO and HCN production and destruction in planetary atmospheres. We highlight ¹O + CH₃ → H₂CO + H as a new key source and H₂CO + ¹O → HCO + OH as a new key sink, for H₂CO in upper planetary atmospheres. In this effort, we develop an oxygen extension to our consistent reduced atmospheric hybrid chemical network (CRAHCN-O), building off our previously developed network for HCN production in N₂-, CH₄-, and H₂-dominated atmospheres (CRAHCN). This extension can be used to simulate both HCN and H₂CO production in atmospheres dominated by any of CO₂, N₂, H₂O, CH₄, and H₂.

Automated Chemical Kinetic Modeling via Hybrid Reactive Molecular Dynamics and Quantum Chemistry Simulations

An automated scheme for obtaining chemical kinetic models from scratch using reactive molecular dynamics and quantum chemistry simulations is presented. This methodology combines the phase space sampling of reactive molecular dynamics with the thermochemistry and kinetics prediction capabilities of quantum mechanics. This scheme provides the NASA polynomial and modified Arrhenius equation parameters for all species and reactions that are observed during the simulation and supplies them in the ChemKin format. The ab initio level of theory for predictions is easily exchangeable, and the presently used G3MP2 level of theory is found to reliably reproduce hydrogen and methane oxidation thermochemistry and kinetics data. Chemical kinetic models obtained with this approach are ready to use for, e.g., ignition delay time simulations, as shown for hydrogen combustion. The presented extension of the ChemTraYzer approach can be used as a basis for methodological advancement of chemical kinetic modeling schemes and as a black-box approach to generate chemical kinetic models.

Predicting pressure-dependent unimolecular rate constants using variational transition state theory with multidimensional tunneling combined with system-specific quantum RRK theory: a definitive test for fluoroform dissociation

Understanding the falloff in rate constants of gas-phase unimolecular reaction rate constants as the pressure is lowered is a fundamental problem in chemical kinetics, with practical importance for combustion, atmospheric chemistry, and essentially all gas-phase reaction mechanisms. In the present work, we use our recently developed system-specific quantum RRK theory, calibrated by canonical variational transition state theory with small-curvature tunneling, combined with the Lindemann-Hinshelwood mechanism, to model the dissociation reaction of fluoroform (CHF3), which provides a definitive test for falloff modeling. Our predicted pressure-dependent thermal rate constants are in excellent agreement with experimental values over a wide range of pressures and temperatures. The present validation of our methodology, which is able to include variational transition state effects, multidimensional tunneling based on the directly calculated potential energy surface along the tunneling path, and torsional and other vibrational anharmonicity, together with state-of-the-art reaction-path-based direct dynamics calculations, is important because the method is less empirical than models routinely used for generating full mechanisms, while also being simpler in key respects than full master equation treatments and the full reduced falloff curve and modified strong collision methods of Troe.

Variational transition state theory: theoretical framework and recent developments

This article reviews the fundamentals of variational transition state theory (VTST), its recent theoretical development, and some modern applications.

Barrierless association of CF2 and dissociation of C2F4 by variational transition-state theory and system-specific quantum Rice-Ramsperger-Kassel theory

Bond dissociation is a fundamental chemical reaction, and the first principles modeling of the kinetics of dissociation reactions with a monotonically increasing potential energy along the dissociation coordinate presents a challenge not only for modern electronic structure methods but also for kinetics theory. In this work, we use multifaceted variable-reaction-coordinate variational transition-state theory (VRC-VTST) to compute the high-pressure limit dissociation rate constant of tetrafluoroethylene (C₂F₄), in which the potential energies are computed by direct dynamics with the M08-HX exchange correlation functional. To treat the pressure dependence of the unimolecular rate constants, we use the recently developed system-specific quantum Rice-Ramsperger-Kassel theory. The calculations are carried out by direct dynamics using an exchange correlation functional validated against calculations that go beyond coupled-cluster theory with single, double, and triple excitations. Our computed dissociation rate constants agree well with the recent experimental measurements.

Performance of First-Principles-Based Reaction Class Transition State Theory

Performance of the Reaction Class Transition State Theory (RC-TST) for prediction of rates constants of elementary reactions is examined using data from its previous applications to a number of different reaction classes. The RC-TST theory is taking advantage of the common structure denominator of all reactions in a given family combined with structure activity relationships to provide a rigorous theoretical framework to obtain rate expression of any reaction within a reaction class in a simple and cost-effective manner. This opens the possibility for integrating this methodology with an automated mechanism generator for "on-the-fly" generation of accurate kinetic models of complex reacting systems.

Silane-initiated nucleation in chemically active plasmas: validation of density functionals, mechanisms, and pressure-dependent variational transition state calculations

Pressure-dependent rate constants for nucleation in nanodusty plasmas are calculated by variational transition state theory with system-specific quantum RRK theory.

New pathways for formation of acids and carbonyl products in low-temperature oxidation: the Korcek decomposition of γ-ketohydroperoxides

We present new reaction pathways relevant to low-temperature oxidation in gaseous and condensed phases. The new pathways originate from γ-ketohydroperoxides (KHP), which are well-known products in low-temperature oxidation and are assumed to react only via homolytic O-O dissociation in existing kinetic models. Our ab initio calculations identify new exothermic reactions of KHP forming a cyclic peroxide isomer, which decomposes via novel concerted reactions into carbonyl and carboxylic acid products. Geometries and frequencies of all stationary points are obtained using the M06-2X/MG3S DFT model chemistry, and energies are refined using RCCSD(T)-F12a/cc-pVTZ-F12 single-point calculations. Thermal rate coefficients are computed using variational transition-state theory (VTST) calculations with multidimensional tunneling contributions based on small-curvature tunneling (SCT). These are combined with multistructural partition functions (Q(MS-T)) to obtain direct dynamics multipath (MP-VTST/SCT) gas-phase rate coefficients. For comparison with liquid-phase measurements, solvent effects are included using continuum dielectric solvation models. The predicted rate coefficients are found to be in excellent agreement with experiment when due consideration is made for acid-catalyzed isomerization. This work provides theoretical confirmation of the 30-year-old hypothesis of Korcek and co-workers that KHPs are precursors to carboxylic acid formation, resolving an open problem in the kinetics of liquid-phase autoxidation. The significance of the new pathways in atmospheric chemistry, low-temperature combustion, and oxidation of biological lipids are discussed.