From the Multiple-Well Master Equation to Phenomenological Rate Coefficients: Reactions on a C3H4 Potential Energy Surface

Theoretical and kinetic study of the H-atom abstraction reactions by Ḣ atom from alkyl cyclohexanes

Reaction kinetics of hydrogen atom abstraction from six alkyl cyclohexanes, methyl cyclohexane (MCH), ethyl cyclohexane (ECH), n-propyl cyclohexane (nPCH), iso-propyl cyclohexane (iPCH), sec-butyl cyclohexane (sBCH) and iso-butyl cyclohexane (iBCH), by the Ḣ atom are systematically studied in this work. The M06-2X method combined with the 6-311++G(d,p) basis set is used to perform geometry optimization, frequency analysis and zero-point energy calculations for all species. The intrinsic reaction coordinate (IRC) calculations are performed to confirm the transition states connecting the reactants and products correctly. One-dimensional hindered rotors are used to treat the low frequency torsional models with potentials scanned at the M06-2X/6-31G level of theory. Electronic single-point energy calculations for all reactants, transition states, and products are performed at the QCISD(T)/CBS level of theory. High-pressure limiting rate constants of 39 reaction channels are obtained using conventional transition state theory with asymmetric Eckart tunneling corrections in the temperature range 298.15-2000 K. Reaction rate rules for H-atom abstraction by the Ḣ atom from alkyl cyclohexanes on primary, secondary and tertiary carbon sites on both the side chain and ring are provided. The obtained rate constants are given by the Arrhenius expression in the temperature range 500-2000 K, which can be used for the combustion kinetics model development for alkyl cyclohexanes.

Pressure and temperature dependent kinetics and the reaction mechanism of Criegee intermediates with vinyl alcohol: a theoretical study

Criegee intermediates (CIs), the key intermediates in the ozonolysis of olefins in atmosphere, have received much attention due to their high activity. The reaction mechanism of the most simple Criegee intermediate CH2OO with vinyl alcohol (VA) was investigated by using the HL//M06-2X/def2TZVP method. The temperature and pressure dependent rate constant and product branching ratio were calculated using the master equation method. For CH2OO + syn-VA, 1,4-insertion is the main reaction channel while for the CH2OO + anti-VA, cycloaddition and 1,2-insertion into the O-H bond are more favorable than the 1,4-insertion reaction. The 1,4-insertion or cycloaddition intermediates are stabilized collisionally at 300 K and 760 torr, and the dissociation products involving OH are formed at higher temperature and lower pressure. The rate constants of the CH2OO reaction with syn-VA and anti-VA both show negative temperature effects, and they are 2.95 × 10-11 and 2.07 × 10-13 cm3 molecule-1 s-1 at 300 K, respectively, and the former is agreement with the prediction in the literature.

Modeling-Experiment-Theory Analysis of Reactions Initiated from Cl + Methyl Formate

Methyl formate (MF; CH3OCHO) is the smallest representative of esters, which are common components of biodiesel. The present study characterizes the thermal dissociation kinetics of the radicals formed by H atom abstraction from MF─CH3OCO and CH2OCHO─through a combination of modeling, experiment, and theory. For the experimental effort, excimer laser photolysis of Cl2 was used as a source of Cl atoms to initiate reactions with MF in the gas phase. Time-resolved species profiles of MF, Cl2, HCl, CO2, CH3, CH3Cl, CH2O, and CH2ClOCHO were measured and quantified using photoionization mass spectrometry at temperatures of 400-750 K and 10 Torr. The experimental data were simulated using a kinetic model, which was informed by ab initio-based theoretical kinetics calculations and included chlorine chemistry and secondary reactions of radical decomposition products. We calculated the rate coefficients for the H-abstraction reactions Cl + MF → HCl + CH3OCO (R1a) and Cl + MF → HCl + CH2OCHO (R1b): k1a,theory = 6.71 × 10-15·T1.14·exp(-606/T) cm3/molecule·s; k1b,theory = 4.67 × 10-18·T2.21·exp(-245/T) cm3/molecule·s over T = 200-2000 K. Electronic structure calculations indicate that the barriers to CH3OCO and CH2OCHO dissociation are 13.7 and 31.6 kcal/mol and lead to CH3 + CO2 (R3) and CH2O + HCO (R5), respectively. The master equation-based theoretical rate coefficients are k3,theory (P = ∞) = 2.94 × 109·T1.21·exp(-6209/T) s-1 and k5,theory (P = ∞) = 8.45 × 108·T1.39·exp(-15132/T) s-1 over T = 300-1500 K. The calculated branching fractions into R1a and R1b and the rate coefficient for R5 were validated by modeling of the experimental species time profiles and found to be in excellent agreement with theory. Additionally, we found that the bimolecular reactions CH2OCHO + Cl, CH2OCHO + Cl2, and CH3 + Cl2 were critical to accurately model the experimental data and constrain the kinetics of MF-radicals. Inclusion of the kinetic parameters determined in this study showed a significant impact on combustion simulations of larger methyl esters, which are considered as biodiesel surrogates.

Theoretical Kinetics studies of isoprene peroxy radical chemistry: The fate of Z-δ-(4-OH, 1-OO)-ISOPOO radical

The OH radical recycling mechanism in isoprene oxidation is one of the most exciting topics in atmospheric chemistry, and the corresponding studies expand our understanding of oxidation mechanisms of volatile organic compounds in the troposphere and provide reliable evidence to improve and develop conventional atmospheric models. In this work, we performed a detailed theoretical kinetics study on the Z-δ-(4-OH, 1-OO)-ISOPOO radical chemistry, which is proposed as the heart of OH recycling in isoprene oxidation. With the full consideration of its accumulation and consumption channels, we studied and discussed the fate of Z-δ-(4-OH, 1-OO)-ISOPOO radical by solving the energy-resolved master equation over a broad range of conditions, including not only room temperatures but also high temperatures of a forest fire or low temperatures and pressures of the upper troposphere. We found non-negligible pressure dependence of its fate at combustion temperatures (up to two orders of magnitude) and demonstrated the significance of both the multi-structural torsional anharmonicity and tunneling for accurately calculating kinetics of the studied system. More interestingly, the tunneling effect on the phenomenological rate constants of the H-shift reaction channel is also found to be pressure-dependent due to the competition with the O2 loss reaction. In addition, our time evolution calculations revealed a two-stage behavior of critical species in this reaction system and estimated the shortest half-lives for the Z-δ-(4-OH, 1-OO)-ISOPOO radical at various temperatures, pressures and altitudes. This detailed kinetics study of Z-δ-(4-OH, 1-OO)-ISOPOO radical chemistry offers a typical example to deeply understand the core mechanism of OH recycling pathways in isoprene oxidation, and provides valuable insights for promoting the development of relevant atmospheric models.

First principles reaction discovery: from the Schrodinger equation to experimental prediction for methane pyrolysis

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theory based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.

Bimolecular Peroxy Radical (RO2) Reactions and Their Relevance in Radical Initiated Oxidation of Hydrocarbons

The kinetics of peroxy radical (RO₂) reactions have been of long-standing interest in atmospheric and combustion chemistry. Nevertheless, the lack of kinetic studies at higher temperatures for their reactions with other radicals such as OH has precluded the inclusion of this class of reactions in detailed kinetics models developed for combustion applications. In this work, guided by the limited room-temperature experimental studies on selected alkyl-peroxy radicals and literature theoretical kinetics on the prototypical CH₃O₂ + OH system, we have performed parametric studies on the effect of uncertainties in the rate coefficients and branching ratios to potential product channels for RO₂ + OH reactions at higher temperatures. Literature kinetics models were used to simulate autoignition delays, laminar flame speeds, and speciation profiles in flow and stirred reactors for a variety of common combustion-relevant fuels. Inclusion of RO₂ + OH reactions was found to retard autoignition in fuel-lean (φ = 0.5) mixtures of ethane and dimethyl ether in air. The observed effects were noticeably more pronounced in ozone-enriched combustion of ethane and dimethyl ether. The simulations also examined the influence of ozone doping levels, pressures, and equivalence ratios for both ethane and dimethyl ether oxidation. Sensitivity and flux analyses revealed that the RO₂ + OH reaction is a significant sink of RO₂ radicals at the early stage of autoignition, affecting fuel oxidation through RO₂ ↔ QOOH, RO₂ ↔ alkene + HO₂, or RO₂ + HO₂ ↔ ROOH + O₂. Additionally, the kinetic stability of the trioxide formed from RO₂ + OH reactions was investigated using master equation analyses. Last, we discuss other bimolecular reactions that are missing in literature kinetics models but are relevant to hydrocarbon oxidation initiated by external radical sources (plasma-enhanced, ozone-enriched combustion, etc.). The present simulations provide a strong motivation for better characterizing the bimolecular kinetics of peroxy radicals.

Correcting the Experimental Enthalpies of Formation of Some Members of the Biologically Significant Sulfenic Acids Family

Sulfenic acids are important intermediates in the oxidation of cysteine thiol groups in proteins by reactive oxygen species. The mechanism is influenced heavily by the presence of polar groups, other thiol groups, and solvent, all of which determines the need to compute precisely the energies involved in the process. Surprisingly, very scarce experimental information exists about a very basic property of sulfenic acids, the enthalpies of formation. In this Article, we use high level quantum chemical methods to derive the enthalpy of formation at 298.15 K of methane-, ethene-, ethyne-, and benzenesulfenic acids, the only ones for which some experimental information exists. The methods employed were tested against well-known experimental data of related species and extensive CCSD(T) calculations. Our best results consistently point out to a much lower enthalpy of formation of methanesulfenic acid, CH₃SOH (Δ_fH⁰(298.15K) = -35.1 ± 0.4 kcal mol^-1), than the one reported in the NIST thermochemical data tables. The enthalpies of formation derived for ethynesulfenic acid, HC≡CSOH, +32.9 ± 1.0 kcal/mol, and benzenesulfenic acid, C₆H₅SOH, -2.6 ± 0.6 kcal mol^-1, also differ markedly from the experimental values, while the enthalpy of formation of ethenesulfenic acid CH₂CHSOH, not available experimentally, was calculated as -11.2 ± 0.7 kcal mol^-1.

Dipolar 1,3-cycloaddition of thioformaldehyde S-methylide (CH2 SCH2 ) to ethylene and acetylene. A comparison with (valence) isoelectronic O3 , SO2 , CH2 OO and CH2 SO

Methods rooted in the density functional theory and in the coupled cluster ansatz were employed to investigate the cycloaddition reactions to ethylene and acetylene of 1,3-dipolar species including ozone and the derivatives issued from replacement of the central oxygen atom by the valence-isoelectronic sulfur atom, and/or of one or both terminal oxygen atoms by the isoelectronic CH₂ group. This gives rise to five different 1,3-dipolar compounds, namely ozone itself (O₃ ), sulfur dioxide (SO₂ ), the simplest Criegee intermediate (CH₂ OO), sulfine (CH₂ SO), and thioformaldehyde S-methylide (CH₂ SCH₂ , TSM). The experimental and accurate theoretical data available for some of those molecules were employed to assess the accuracy of two last-generation composite methods employing conventional or explicitly correlated post-Hartree-Fock contributions (jun-Cheap and SVECV-f12, respectively), which were then applied to investigate the reactivity of TSM. The energy barriers provided by both composite methods are very close (the average values for the two composite methods are 7.1 and 8.3 kcal mol^-1 for the addition to ethylene and acetylene, respectively) and comparable to those ruling the corresponding additions of ozone (4.0 and 7.7 kcal mol^-1 , respectively). These and other evidences strongly suggest that, at least in the case of cycloadditions, the reactivity of TSM is similar to that of O₃ and very different from that of SO₂ .

NAST: Nonadiabatic Statistical Theory Package for Predicting Kinetics of Spin-Dependent Processes

We present a nonadiabatic statistical theory (NAST) package for predicting kinetics of spin-dependent processes, such as intersystem crossings, spin-forbidden unimolecular reactions, and spin crossovers. The NAST package can calculate the probabilities and rates of transitions between the electronic states of different spin multiplicities. Both the microcanonical (energy-dependent) and canonical (temperature-dependent) rate constants can be obtained. Quantum effects, including tunneling, zero-point vibrational energy, and reaction path interference, can be accounted for. In the limit of an adiabatic unimolecular reaction proceeding on a single electronic state, NAST reduces to the traditional transition state theory. Because NAST requires molecular properties at only a few points on potential energy surfaces, it can be applied to large molecular systems, used with accurate high-level electronic structure methods, and employed to study slow nonadiabatic processes. The essential NAST input data include the nuclear Hessian at the reactant minimum, as well as the nuclear Hessians, energy gradients, and spin-orbit coupling at the minimum energy crossing point (MECP) between two states. The additional computational tools included in the NAST package can be used to extract the required input data from the output files of electronic structure packages, calculate the effective Hessian at the MECP, and fit the reaction coordinate for more advanced NAST calculations. We describe the theory, its implementation, and three examples of application to different molecular systems.

Master Equation Study of Hydrogen Abstraction from HCHO by OH Via a Chemically Activated Intermediate

The abstraction reaction of hydrogen from formaldehyde by OH radical plays an important role in formaldehyde oxidation. The reaction involves a bimolecular association to form a chemically activated hydrogen-bonded reaction...

Unimolecular isomerisation of 1,5-hexadiyne observed by threshold photoelectron photoion coincidence spectroscopy

The unimolecular isomerisation of the prompt propargyl + propargyl "head-to-head" adduct, 1,5- hexadiyne, to fulvene and benzene by the 3,4-dimethylenecyclobut-1-ene (DMCB) intermediate (all C6H6) was studied in the high-pressure limit...

Spiers Memorial Lecture: theory of unimolecular reactions

One hundred years ago, at an earlier Faraday Discussion meeting, Lindemann presented a mechanism that provides the foundation for contemplating the pressure dependence of unimolecular reactions. Since that time, our...

Initiation reactions in the high temperature decomposition of styrene

The thermal decomposition of styrene was investigated in a combined experimental, theory and modeling study with particular emphasis placed on the initial dissociation reactions. Two sets of shock tube/time-of-flight mass spectrometry (TOF-MS) experiments were performed to identify reaction products and their order of appearance. One set of experiments was conducted with a miniature high repetition rate shock tube at the Advanced Light Source at Lawrence Berkeley National Laboratory using synchrotron vacuum ultraviolet photoionization. The other set of experiments was performed in a diaphragmless shock tube (DFST) using electron impact ionization. The datasets span 1660–2260 K and 0.5–12 atm. The results show a marked transition from aromatic products at low temperatures to polyacetylenes, up to C₈H₂, at high temperatures. The TOF-MS experiments were complemented by DFST/LS (laser schlieren densitometry) experiments covering 1800–2250 K and 60–240 Torr. These were particularly sensitive to the initial dissociation reactions. These reactions were investigated theoretically and revealed the dissociation of styrene to be a complex multichannel process with strong pressure and temperature dependencies that were evaluated with multi-well master equation simulations. Simulations of the LS data with a mechanism developed in this work are in excellent agreement with the experimental data. From these simulations, rate coefficients for the dissociation of styrene were obtained that are in good agreement with the theoretical predictions. The simulation results also provide fair predictions of the temperature and pressure dependencies of the products observed in the TOF-MS studies. Prior experimental studies of styrene pyrolysis concluded that the main products were benzene and acetylene. In contrast, this study finds that the majority of styrene dissociates to create five styryl radical isomers. Of these, α-styryl accounts for about 50% with the other isomers consuming approximately 20%. It was also found that C–C bond scission to phenyl and vinyl radicals consumes up to 25% of styrene. Finally the dissociation of styrene to benzene and vinylidene accounts for roughly 5% of styrene consumption. Comments are made on the apparent differences between the results of this work and prior literature.

A combined theoretical and experimental study showing styrene primarily decomposes to styryl radicals + H.

Gas-phase synthesis of benzene via the propargyl radical self-reaction

Polycyclic aromatic hydrocarbons (PAHs) have been invoked in fundamental molecular mass growth processes in our galaxy. We provide compelling evidence of the formation of the very first ringed aromatic and building block of PAHs-benzene-via the self-recombination of two resonantly stabilized propargyl (C₃H₃) radicals in dilute environments using isomer-selective synchrotron-based mass spectrometry coupled to theoretical calculations. Along with benzene, three other structural isomers (1,5-hexadiyne, fulvene, and 2-ethynyl-1,3-butadiene) and o-benzyne are detected, and their branching ratios are quantified experimentally and verified with the aid of computational fluid dynamics and kinetic simulations. These results uncover molecular growth pathways not only in interstellar, circumstellar, and solar systems environments but also in combustion systems, which help us gain a better understanding of the hydrocarbon chemistry of our universe.

Substitution Reactions in the Pyrolysis of Acetone Revealed through a Modeling, Experiment, Theory Paradigm

The development of high-fidelity mechanisms for chemically reactive systems is a challenging process that requires the compilation of rate descriptions for a large and somewhat ill-defined set of reactions. The present unified combination of modeling, experiment, and theory provides a paradigm for improving such mechanism development efforts. Here we combine broadband rotational spectroscopy with detailed chemical modeling based on rate constants obtained from automated ab initio transition state theory-based master equation calculations and high-level thermochemical parametrizations. Broadband rotational spectroscopy offers quantitative and isomer-specific detection by which branching ratios of polar reaction products may be obtained. Using this technique, we observe and characterize products arising from H atom substitution reactions in the flash pyrolysis of acetone (CH₃C(O)CH₃) at a nominal temperature of 1800 K. The major product observed is ketene (CH₂CO). Minor products identified include acetaldehyde (CH₃CHO), propyne (CH₃CCH), propene (CH₂CHCH₃), and water (HDO). Literature mechanisms for the pyrolysis of acetone do not adequately describe the minor products. The inclusion of a variety of substitution reactions, with rate constants and thermochemistry obtained from automated ab initio kinetics predictions and Active Thermochemical Tables analyses, demonstrates an important role for such processes. The pathway to acetaldehyde is shown to be a direct result of substitution of acetone's methyl group by a free H atom, while propene formation arises from OH substitution in the enol form of acetone by a free H atom.

Discovery of the propargyl radical (CH2CCH) in TMC-1: one of the most abundant radicals ever found and a key species for cyclization to benzene in cold dark clouds

We present the first identification in interstellar space of the propargyl radical (CH2CCH). This species was observed in the cold dark cloud TMC-1 using the Yebes 40m telescope. The six strongest hyperfine components of the 20,2-10,1 rotational transition, lying at 37.46 GHz, were detected with signal-to-noise ratios in the range 4.6-12.3 σ. We derive a column density of 8.7 × 1013 cm-2 for CH2CCH, which translates to a fractional abundance relative to H2 of 8.7 × 10-9. This radical has a similar abundance to methyl acetylene, with an abundance ratio CH2CCH/CH3CCH close to one. The propargyl radical is thus one of the most abundant radicals detected in TMC-1, and it is probably the most abundant organic radical with a certain chemical complexity ever found in a cold dark cloud. We constructed a gas-phase chemical model and find calculated abundances that agree with, or fall two orders of magnitude below, the observed value depending on the poorly constrained low-temperature reactivity of CH2CCH with neutral atoms. According to the chemical model, the propargyl radical is essentially formed by the C + C2H4 reaction and by the dissociative recombination of C3Hn + ions with n = 4-6. The propargyl radical is believed to control the synthesis of the first aromatic ring in combustion processes, and it probably plays a key role in the synthesis of large organic molecules and cyclization processes to benzene in cold dark clouds.

Experimental and Theoretical Investigations of the Radical-Radical Reaction: N2H3 + NO2

The N₂H₃ + NO₂ reaction plays a key role during the early stages of hypergolic ignition between N₂H₄ and N₂O₄. Here for the first time, the reaction kinetics of N₂H₃ in excess NO₂ was studied in 2.0 Torr of N₂ and in the narrow temperature range 298-348 K in a pulsed photolysis flow-tube reactor coupled to a mass spectrometer. The temporal profile of the product, HONO, was determined by direct detection of the m/z +47 amu ion signal. For each chosen [NO₂], the observed [HONO] trace was fitted to a biexponential kinetics expression, which yielded a value for the pseudo-first-order rate coefficient, k', for the reaction of N₂H₃ with NO₂. The slope of the plot of k' versus [NO₂] yielded a value for the observed bimolecular rate coefficient, k_obs, which could be fitted to an Arrhenius expression of (2.36 ± 0.47) × 10^-12 exp((520 ± 350)/T) cm³ molecule^-1 s^-1. The errors are 1σ and include estimated uncertainties in the NO₂ concentration. The potential energy surface of N₂H₃ + NO₂ was investigated by advanced ab initio quantum chemistry theories. It was found that the reaction occurs via a complex reaction mechanism, and all of the reaction channels have transition state energies below that of the entrance asymptote. The radical-radical addition forms the N₂H₃NO₂ adducts, while roaming-mediated isomerization reactions yield the N₂H₃ONO isomers, which undergo rapid dissociation reactions to several sets of distinct products. The RRKM multiwell master equation simulations revealed that the major product channel involves the formation of trans-HONO and trans-N₂H₂ below 500 K and the formation of NO + NH₂NHO above 500 K, which is nearly pressure independent. The pressure-dependent rate coefficients of the product channels were computed over a wide pressure-temperature range, which encompassed the experimental data.

Low-Pressure Limit of Competitive Unimolecular Reactions

Barker and Ortiz found unusual falloff effects in the flux coefficients of the competitive unimolecular reactions of 2-methylhexyl radicals, and they concluded that this might have important effects on the rate constants of reactions with higher thresholds. To study this effect, we carried out master equation calculations of the same reaction system to learn whether this effect shows up in measurable rate constants, and the answer is yes. We also studied specially designed mechanisms to reveal that the various reactive pathways connecting the reagents can have a large effect on the rate constants, causing them to be quite different than if the reactions proceeded independently, and that reactions with significantly higher barriers may nevertheless have larger rate constants. This provides a new perspective for interpreting and predicting the kinetics of competitive unimolecular reactions.

An experimental and theoretical study of the high temperature reactions of the four butyl radical isomers

The high temperature gas phase chemistry of the four butyl radical isomers (n-butyl, sec-butyl, iso-butyl, and tert-butyl) was investigated in a combined experimental and theoretical study. Organic nitrites were used as convenient and clean sources of each of the butyl radical isomers. Rate coefficients for dissociation of each nitrite were obtained experimentally and are at, or close to, the high pressure limit. Low pressure experiments were performed in a diaphragmless shock tube with laser schlieren densitometry at post-shock pressures of 65, 130, and 260 Torr and post-shock temperatures of 700-1000 K. Additional experiments were conducted with iso-butyl radicals at 805 K and 8.7 bar to elucidate changes in mechanism at higher pressures. These experiments were performed in a miniature shock tube with synchrotron-based photoionization mass spectrometry. The mass spectra confirmed that scission of the O-NO bond is the primary channel by which the precursors dissociate, but they also provided evidence of a minor channel (<7.7%) through HNO loss and formation of an aldehyde. These high pressure experiments were also used to determine the disproportionation/recombination ratio for iso-butyl radicals as 0.3. Reanalysis of the lower-temperature literature and the present data yielded rate constants for the disproportionation reaction, iso-butyl + iso-butyl = iso-butene + iso-butane. A chemical kinetics model was developed for the reactions of the butyl isomers that included new paths for highly energized adducts. These adducts are formed by the addition of H, CH₃ or C₂H₅ to the butyl radicals. Accompanying theoretical investigations show that chemically activated pathways are competitive with stabilization of the adduct by collision under the conditions of the laser schlieren experiments. These calculations also show that at 10 bar and T < 1000 K stabilization is the only important reaction, but at higher temperatures, even at 10 bar, chemically activated product channels should also be considered. Branching fractions and rate coefficients are presented for these reactions. This study also highlights the importance of the radical structure for determining branching ratios for disproportionation and recombination of alkyl radicals, and these were facilitated by theoretical calculations of recombination rate coefficients for the four butyl radical isomers. The results reveal previously unknown features of butyl radical chemistry under conditions that are relevant to a wide range of applications and reaction mechanisms are presented that incorporate pressure dependent rate coefficients for the key steps.

A theoretical investigation on the atmospheric degradation of the radical: reactions with NO, NO2, and NO3

The radical is the key intermediate in the atmospheric oxidation of benzaldehyde, and its further chemistry contributes to local air pollution. The reaction mechanisms of the radical with NO, NO₂, and NO₃ were studied by quantum chemistry calculations at the CCSD(T)/CBS//M06-2X/def2-TZVP level of theory. The explicit potential energy curves were provided in order to reveal the atmospheric fate of the radical comprehensively. The main products of the reaction of with NO are predicted to be , CO₂ and NO₂. The reaction of with NO₂ is reversible, and its main product would be C₆H₅C(O)O₂NO₂ which was predicted to be more stable than PAN (peroxyacetyl nitrate) at room temperature. The decomposition of C₆H₅C(O)O₂NO₂ at different ambient temperatures would be a potential long-range transport source of NO_x in the atmosphere. The predominant products of the reaction are predicted to be C₆H₅C(O)O₂H, C₆H₅C(O)OH, O₂ and O₃, while HO˙ is of minor importance. So, the reaction of with would be an important source of ozone and carboxylic acids in the local atmosphere, and has less contribution to the regeneration of HO˙ radicals. The reaction of with NO₃ should mainly produce , CO₂, O₂ and NO₂, which might play an important role in atmospheric chemistry of peroxy radicals at night, but has less contribution to the night-time conversion of ( and RO˙) to ( and HO˙) in the local atmosphere. The results above are in good accordance with the reported experimental observations.

Direct Trace Fitting of Experimental Data Using the Master Equation: Testing Theory and Experiments on the OH + C2H4 Reaction

Laser flash photolysis coupled with laser-induced fluorescence observation of OH has been used to observe the equilibration of OH + C₂H₄ ↔ HOC₂H₄ over the temperature range 563-723 K and pressures of bath gas (N₂) from 58 to 250 Torr. The time-resolved OH traces have been directly and globally fitted with a master equation in order to extract Δ_RH₀⁰, the binding energy of the HOC₂H₄ adduct, with respect to reagents. The global approach allows the role that OH abstraction plays at higher temperatures to be identified. The resultant value ofΔ_RH₀⁰, 111.8 kJ mol^-1, is determined to be better than 2 kJ mol^-1 and is in agreement with our ab initio calculations (carried out at the CCSD(T)/CBS//M06-2X/aug-cc-pVTZ level), 111.4 kJ mol^-1, and other state of the art calculations. Parameters for the abstraction channel are also in good agreement with previous experimental studies. To effect this analysis, the MESMER master equation code was extended to directly incorporate secondary chemistry: diffusional loss from the observation region and reaction with the photolytic precursor. These extensions, which, among other things, resolve issues related to the merging of chemically significant and internal energy relaxation eigenvalues, are presented.

Ab Initio Kinetics of Methylamine Radical Thermal Decomposition and H-Abstraction from Monomethylhydrazine by H-Atom

Methylamine radicals (CH₃NH) and amino radicals (NH₂) are major products in the early pyrolysis/ignition of monomethylhydrazine (CH₃NHNH₂). Ab initio kinetics of thermal decomposition of CH₃NH radicals was analyzed by RRKM master equation simulations. It was found that β-scission of the methyl H-atom from CH₃NH radicals is predominant and fast enough to induce subsequent H-abstraction reactions in CH₃NHNH₂ to trigger ignition. Consequently, the kinetics of H-abstraction reactions from CH₃NHNH₂ by H-atoms was further investigated. It was found that the energy barriers for abstraction of the central amine H-atom, two terminal amine H-atoms, and methyl H-atoms are 4.16, 2.95, 5.98, and 8.50 kcal mol^-1, respectively. In units of cm³ molecule^-1 s^-1, the corresponding rate coefficients were found to be k₈ = 9.63 × 10^-20T^2.596 exp(-154.2/T), k₉ = 2.04 × 10^-18T^2.154 exp(104.1/T), k₁₀ = 1.13 × 10^-20T^2.866 exp(-416.3/T), and k₁₁ = 2.41 × 10^-23T^3.650 exp(-870.5/T), respectively, in the 290-2500 K temperature range. The results reveal that abstraction of the terminal amine H-atom to form trans-CH₃NHNH radicals is the dominant channel among the different abstraction channels. At 298 K, the total theoretical H-abstraction rate coefficient, calculated with no adjustable parameters, is 8.16 × 10^-13 cm³ molecule^-1 s^-1, which is in excellent agreement with Vaghjiani's experimental observation of (7.60 ± 1.14) × 10^-13 cm³ molecule^-1 s^-1 ( J. Phys. Chem. A 1997, 101, 4167-4171, DOI: 10.1021/jp964044z).

Ab Initio/Transition-State Theory Study of the Reactions of Ċ5H9 Species of Relevance to 1,3-Pentadiene, Part I: Potential Energy Surfaces, Thermochemistry, and High-Pressure Limiting Rate Constants

In this study, the reactions of Ċ₅H₉ radicals are theoretically investigated, with a particular emphasis on hydrogen atom addition reactions to 1,3-pentadiene (C₅H₈) to form Ċ₅H₉ radicals, although the subsequent isomerization and decomposition reactions of the Ċ₅H₉ radicals are also of direct relevance to the radicals formed from the pyrolysis and oxidation of species including pentene and cyclopentane. Moreover, H-atom abstraction reactions by hydrogen atoms from 1,3-pentadiene are also investigated. The geometries and frequencies of 63 potential energy surface (PES) minima and 88 transition states are optimized at the ωB97XD/aug-cc-pVTZ level of theory. Spin-unrestricted open-shell single-point energies for all the species are calculated at the CCSD(T)/aug-cc-pVTZ level of theory with basis set corrections from MP2/aug-cc-pVXZ (where X = T and Q). A one-dimensional hindered rotor treatment is employed for torsional modes, with the M06-2X/6-311++G(D,P) method used to compute the potential energy as a function of the dihedral angle. The high-pressure limiting rate constants and the thermochemical properties for C₅ species are calculated using the Master Equation System Solver (MESS) with conventional transition-state theory and comparisons made with existing available literature data. A hydrogen atom can add to the terminal carbon atom of 1,3-pentadiene to form the 2,4-Ċ₅H₉ radical and/or the internal carbon atoms to form 2,5-Ċ₅H₉, 1,4-Ċ₅H₉, and 1,3-Ċ₅H₉ radicals. Among the four entrance channels for Ḣ atom addition reactions, the formation of 2,4-Ċ₅H₉ and 1,3-Ċ₅H₉ radicals is more exothermic in comparison to the other Ċ₅H₉ isomers (2,5-Ċ₅H₉, 1,4-Ċ₅H₉) because of the resonantly stabilized allylic structure. Consequently, the formation of the former is generally dominant in terms of barrier heights. Ḣ atom addition reactions to 1,3-pentadiene are compared to available C₃-C₅ alkenes and dienes, with external addition calculated to be kinetically favored over internal addition. However, the correlation between heats of formation and energy barriers for Ḣ atom addition to 1,2-dienes is different from that for 1,3- and 1,4-dienes. Hydrogen atom addition and abstraction rate constants are also compared for 1,3-pentadiene, with addition found to be dominant. The subsequent unimolecular reactions on the Ċ₅H₉ PES are found to be highly complex with reactions taking place on a multiple-well multiple-channel PES. For clarity, the reaction mechanism and kinetics of each Ċ₅H₉ radical are discussed individually in terms of the computed enthalpy of the reaction and activation, the transition-state structure/reaction class, and also in terms of the combustion species for which the reactions are of potential importance. The reactions on the Ċ₅H₉ PES are divided into three reaction classes (H-shift isomerization, cycloaddition, and β-scission reactions), and the reactivity-structure-based estimation rules for energy barriers are derived for these three reaction classes and compared to literature results for alkyl radicals.

Mechanism and kinetics of the oxidation of dimethyl carbonate by hydroxyl radical in the atmosphere

The mechanism and kinetics for the reaction of dimethyl carbonate (DMC) with OH radical have been studied by using quantum chemical methods. Four reaction pathways were identified for the initial reaction. In the first two pathways, hydrogen atom abstraction is taking place and alkyl radical intermediate is formed with the energy barrier of 6.4 and 7.9 kcal/mol. In the third pathway, OH addition reaction to the carbonyl carbon (C2) atom of DMC and intermediate, I2, is formed with an energy barrier of 11.9 kcal/mol. In the fourth pathway, along with CH₃O^●, methyl hydrogen carbonate is formed. For this C-O bond breaking and O-H addition reaction, the energy barrier is 27 kcal/mol. The calculated enthalpy and Gibbs energy values show that the studied initial reactions are exothermic and exoergic except the OH addition reaction. For the initial reactions, the rate constants were calculated by using canonical variational transition state theory (CVT) with small curvature tunneling (SCT) correction over the temperature range of 278-1200 K. At 298 K, the calculated rate coefficient for the in-plane and out-of-plane hydrogen atom abstraction reaction pathway is 2.30 × 10^-13 and 0.02 × 10^-13 cm³ molecule^-1 s^-1. Further, the reaction between alkyl radical intermediate formed from the first pathway and O₂ is studied. The reaction of alkyl peroxy radical intermediate with atmospheric oxidants, HO₂, NO, and NO₂ is also studied. It was found that the formic (methyl carbonic) anhydride is the end product formed from the atmospheric oxidation and secondary reactions of DMC.

Theoretical investigation on the reaction mechanism and kinetics of a Criegee intermediate with ethylene and acetylene

The competition among the possible pathways, the branching ratios of the adduct and the decomposition products at different temperatures and pressures have been evaluated.

Decomposition kinetics for HONO and HNO2

This work presents a detailed investigation into the isomerization and decomposition of HONO and HNO₂.

Modeling of aromatics formation in fuel-rich methane oxy-combustion with an automatically generated pressure-dependent mechanism

An automatic generated mechanism for methane-rich combustion captures the chemistry from small molecules to three-ring aromatic species.

Potential Energy Surface for Large Barrierless Reaction Systems: Application to the Kinetic Calculations of the Dissociation of Alkanes and the Reverse Recombination Reactions

The isodesmic reaction method is applied to calculate the potential energy surface (PES) along the reaction coordinates and the rate constants of the barrierless reactions for unimolecular dissociation reactions of alkanes to form two alkyl radicals and their reverse recombination reactions. The reaction class is divided into 10 subclasses depending upon the type of carbon atoms in the reaction centers. A correction scheme based on isodesmic reaction theory is proposed to correct the PESs at UB3LYP/6-31+G(d,p) level. To validate the accuracy of this scheme, a comparison of the PESs at B3LYP level and the corrected PESs with the PESs at CASPT2/aug-cc-pVTZ level is performed for 13 representative reactions, and it is found that the deviations of the PESs at B3LYP level are up to 35.18 kcal/mol and are reduced to within 2 kcal/mol after correction, indicating that the PESs for barrierless reactions in a subclass can be calculated meaningfully accurately at a low level of ab initio method using our correction scheme. High-pressure limit rate constants and pressure dependent rate constants of these reactions are calculated based on their corrected PESs and the results show the pressure dependence of the rate constants cannot be ignored, especially at high temperatures. Furthermore, the impact of molecular size on the pressure-dependent rate constants of decomposition reactions of alkanes and their reverse reactions has been studied. The present work provides an effective method to generate meaningfully accurate PESs for large molecular system.

Reaction kinetics of OH + HNO3 under conditions relevant to the upper troposphere/lower stratosphere

Key upper atmosphere reaction of HNO₃ + OH studied over extended pressure and temperature range using new alternative detection method.

Full rate constant matrix contraction method for obtaining branching ratio of unimolecular decomposition

The branching ratio of unimolecular decomposition can be evaluated by solving the rate equations. Recent advances in automated reaction path search methods have enabled efficient construction of the rate equations based on quantum chemical calculations. However, it is still difficult to solve the rate equations composed of hundreds or more elementary steps. This problem is especially serious when elementary steps that occur in highly different timescales coexist. In this article, we introduce an efficient approach to obtain the branching ratio from a given set of rate equations. It has been derived from a recently proposed rate constant matrix contraction (RCMC) method, and termed full-RCMC (f-RCMC). The f-RCMC gives the branching ratio without solving the rate equations. Its performance was tested numerically for unimolecular decomposition of C₃ H₅ and C₄ H₅ . Branching ratios obtained by the f-RCMC precisely reproduced the values obtained by numerically solving the rate equations. It took about 95 h to solve the rate equations of C₄ H₅ consisting of 234 elementary steps. In contrast, the f-RCMC gave the branching ratio in less than 1 s. The f-RCMC would thus be an efficient alternative of the conventional kinetic simulation approach. © 2016 Wiley Periodicals, Inc.

Spectroscopy and Photochemistry of Triplet 1,3-Dimethylpropynylidene (MeC3Me)

Photolysis (λ > 472 nm) of 2-diazo-3-pentyne (11) affords triplet 1,3-dimethylpropynylidene (MeC3Me, (3)3), which was characterized spectroscopically in cryogenic matrices. The infrared, electronic absorption, and electron paramagnetic resonance spectra of MeC3Me ((3)3) are compared with those of the parent system (HC3H) to ascertain the effect of alkyl substituents on delocalized carbon chains of this type. Quantum chemical calculations (CCSD(T)/ANO1) predict an unsymmetrical equilibrium structure for triplet MeC3Me ((3)3), but they also reveal a very shallow potential energy surface. The experimental IR spectrum of triplet MeC3Me ((3)3) is best interpreted in terms of a quasilinear, axially symmetric structure. EPR spectra yield zero-field splitting parameters that are typical for triplet carbenes with axial symmetry (|D/hc| = 0.63 cm(-1), |E/hc| = ∼ 0 cm(-1)), while theoretical analysis suggests that the methyl substituents confer significant spin polarization to the carbon chain. Upon irradiation into the near-UV electronic absorption (λmax 350 nm), MeC3Me ((3)3) undergoes 1,2-hydrogen migration to yield pent-1-en-3-yne (4), a photochemical reaction that is typical of carbenes bearing a methyl substituent. This facile process apparently precludes photoisomerization to other interesting C5H6 isomers, in contrast to the rich photochemistry of the parent C3H2 system.