Quasi-Classical Trajectory Calculation of Rate Constants Using an Ab Initio Trained Machine Learning Model (aML-MD) with Multifidelity Data

Machine learning (ML) provides a great opportunity for the construction of models with improved accuracy in classical molecular dynamics (MD). However, the accuracy of a ML trained model is limited by the quality and quantity of the training data. Generating large sets of accurate ab initio training data can require significant computational resources. Furthermore, inconsistent or incompatible data with different accuracies obtained using different methods may lead to biased or unreliable ML models that do not accurately represent the underlying physics. Recently, transfer learning showed its potential for avoiding these problems as well as for improving the accuracy, efficiency, and generalization of ML models using multifidelity data. In this work, ab initio trained ML-based MD (aML-MD) models are developed through transfer learning using DFT and multireference data from multiple sources with varying accuracy within the Deep Potential MD framework. The accuracy of the force field is demonstrated by calculating rate constants for the H + HO2 → H2 + 3O2 reaction using quasi-classical trajectories. We show that the aML-MD model with transfer learning can accurately predict the rate constants while reducing the computational cost by more than five times compared to the use of more expensive quantum chemistry training data sets. Hence, the aML-MD model with transfer learning shows great potential in using multifidelity data to reduce the computational cost involved in generating the training set for these potentials.

Isomer-resolved unimolecular dynamics of the hydroperoxyalkyl intermediate (•QOOH) in cyclohexane oxidation

The oxidation of cycloalkanes is important in the combustion of transportation fuels and in atmospheric secondary organic aerosol formation. A transient carbon-centered radical intermediate (•QOOH) in the oxidation of cyclohexane is identified through its infrared fingerprint and time- and energy-resolved unimolecular dissociation dynamics to hydroxyl (OH) radical and bicyclic ether products. Although the cyclohexyl ring structure leads to three nearly degenerate •QOOH isomers (β-, γ-, and δ-QOOH), their transition state (TS) barriers to OH products are predicted to differ considerably. Selective characterization of the β-QOOH isomer is achieved at excitation energies associated with the lowest TS barrier, resulting in rapid unimolecular decay to OH products that are detected. A benchmarking approach is employed for the calculation of high-accuracy stationary point energies, in particular TS barriers, for cyclohexane oxidation (C6H11O2), building on higher-level reference calculations for the smaller ethane oxidation (C2H5O2) system. The isomer-specific characterization of β-QOOH is validated by comparison of experimental OH product appearance rates with computed statistical microcanonical rates, including significant heavy-atom tunneling, at energies in the vicinity of the TS barrier. Master-equation modeling is utilized to extend the results to thermal unimolecular decay rate constants at temperatures and pressures relevant to cyclohexane combustion.

Phenalenyl growth reactions and implications for prenucleation chemistry of aromatics in flames

The energetics and kinetics of phenalene and phenalenyl growth reactions were studied theoretically. Rate constants of phenalene and phenalenyl H-abstraction and C2H2 addition to the formed radicals were evaluated through quantum-chemical and rate-theory calculations. The obtained values, assigned to all π radicals, were tested in deterministic and kinetic Monte Carlo simulations of aromatics growth under conditions of laminar premixed flames. Kekulé and non-Kekulé structures of the polycyclic aromatic hydrocarbons (PAHs) evolving in the stochastic simulations were identified by on-the-fly constrained optimization. The numerical results demonstrated an increased PAH growth and qualitatively reproduced experimental observations of Homann and co-workers of non-decaying PAH concentrations with nearly equal abundances of even and odd carbon-atom PAHs. The analysis revealed that the PAH growth proceeds via alternating and sterically diverse acetylene and methyl HACA additions. The rapid and diverse spreading in the PAH population supports a nucleation model as PAH dimerization, assisted by the non-equilibrium phenomena, forming planar aromatics first and then transitioning to the PAH-PAH stacking with size.

Coherences of Photoinduced Electron Spin Qubit Pair States in Photosystem I

This publication presents the first comprehensive experimental study of electron spin coherences in photosynthetic reaction center proteins, specifically focusing on photosystem I (PSI). The ultrafast electron transfer in PSI generates spin-correlated radical pairs (SCRPs), which are entangled spin pairs formed in well-defined spin states (Bell states). Since their discovery in our group in the 1980s, SCRPs have been extensively used to enhance our understanding of structure-function relationships in photosynthetic proteins. More recently, SCRPs have been utilized as tools for quantum sensing. Electron spin decoherence poses a significant challenge in realizing practical applications of electron spin qubits, particularly the creation of quantum entanglement between multiple electron spins. This work is focused on the systematic characterization of decoherence in SCRPs of PSI. These decoherence times were measured as electron spin echo decay times, termed phase memory times (TM), at various temperatures. Decoherence was recorded on both transient SCRP states P700+A1- and thermalized states. Our study reveals that TM exhibits minimal dependence on the biological species, biochemical treatment, and paramagnetic species. The analysis indicates that nuclear spin diffusion and instantaneous diffusion mechanisms alone cannot explain the observed decoherence. As a plausible explanation we discuss the assumption that the low-temperature dynamics of methyl groups in the protein surrounding the unpaired electron spin centers is the main factor governing the loss of the spin coherence in PSI.

Direct observation of coherence transfer and rotational-to-vibrational energy exchange in optically centrifuged CO2 super-rotors

Optical centrifuges are laser-based molecular traps that can rotationally accelerate molecules to energies rivalling or exceeding molecular bond energies. Here we report time and frequency-resolved ultrafast coherent Raman measurements of optically centrifuged CO₂ at 380 Torr spun to energies beyond its bond dissociation energy of 5.5 eV (J_max = 364, E_rot = 6.14 eV, E_rot/k_B = 71, 200 K). The entire rotational ladder from J = 24 to J = 364 was resolved simultaneously which enabled a more accurate measurement of the centrifugal distortion constants for CO₂. Remarkably, coherence transfer was directly observed, and time-resolved, during the field-free relaxation of the trap as rotational energy flowed into bending-mode vibrational excitation. Vibrationally excited CO₂ (ν₂ > 3) was observed in the time-resolved spectra to populate after 3 mean collision times as a result of rotational-to-vibrational (R-V) energy transfer. Trajectory simulations show an optimal range of J for R-V energy transfer. Dephasing rates for molecules rotating up to 5.5 times during one collision were quantified. Very slow decays of the vibrational hot band rotational coherences suggest that they are sustained by coherence transfer and line mixing.

Radical-Radical Reactions in Molecular Weight Growth: The Phenyl + Propargyl Reaction

The mechanism for hydrocarbon ring growth in sooting environments is still the subject of considerable debate. The reaction of phenyl radical (C₆H₅) with propargyl radical (H₂CCCH) provides an important prototype for radical-radical ring-growth pathways. We studied this reaction experimentally over the temperature range of 300-1000 K and pressure range of 4-10 Torr using time-resolved multiplexed photoionization mass spectrometry. We detect both the C₉H₈ and C₉H₇ + H product channels and report experimental isomer-resolved product branching fractions for the C₉H₈ product. We compare these experiments to theoretical kinetics predictions from a recently published study augmented by new calculations. These ab initio transition state theory-based master equation calculations employ high-quality potential energy surfaces, conventional transition state theory for the tight transition states, and direct CASPT2-based variable reaction coordinate transition state theory (VRC-TST) for the barrierless channels. At 300 K only the direct adducts from radical-radical addition are observed, with good agreement between experimental and theoretical branching fractions, supporting the VRC-TST calculations of the barrierless entrance channel. As the temperature is increased to 1000 K we observe two additional isomers, including indene, a two-ring polycyclic aromatic hydrocarbon, and a small amount of bimolecular products C₉H₇ + H. Our calculated branching fractions for the phenyl + propargyl reaction predict significantly less indene than observed experimentally. We present further calculations and experimental evidence that the most likely cause of this discrepancy is the contribution of H atom reactions, both H + indenyl (C₉H₇) recombination to indene and H-assisted isomerization that converts less stable C₉H₈ isomers into indene. Especially at low pressures typical of laboratory investigations, H-atom-assisted isomerization needs to be considered. Regardless, the experimental observation of indene demonstrates that the title reaction leads, either directly or indirectly, to the formation of the second ring in polycyclic aromatic hydrocarbons.

Formation of phosphorus monoxide through the [Formula: see text] reaction

Phosphorus is a key and vital element for a diverse set of important biological molecules, being indispensable for life as we know. A deeper comprehension of its role in astrochemistry and atmospheric chemistry may aid in finding answers to how this element became available on Earth. The PO molecule is one of the main reservoirs of phosphorus in the interstellar medium (ISM), and a better understanding of the mechanisms and rate coefficients for its formation in the ISM is important for modelling its abundances. In this work, we perform multireference configuration interaction calculations on the formation of PO via the [Formula: see text] reaction, analyzing its potential energy surface and rate coefficients for the global reaction on both doublet and quartet states. We also perform DFT (M06-2X) and CCSD(T) calculations, in order to compare the results. We found that the OPO system possesses a high multiconfigurational character, making DFT and CCSD methodologies not suitable for its potential energy landscape calculation. The rate coefficients have been calculated using the master equation system solver (MESS) package, and the results compared to recent experimental data. It is shown that the quartet state contributes for temperatures higher than 700K. The computed rate coefficient can be described by a modified Arrhenius equation [[Formula: see text]] with [Formula: see text], [Formula: see text] and [Formula: see text] K.

Identification of the acetaldehyde oxide Criegee intermediate reaction network in the ozone-assisted low-temperature oxidation of trans-2-butene

Uni- and bi-molecular reactions involving Criegee intermediates (CIs) have been the focus of many studies due to the role these molecules play in atmospheric chemistry. The reactivity of CIs is known to strongly depend on their structure. The reaction network of the second simplest CI, acetaldehyde oxide (CH₃CHOO), is investigated in this work in an atmospheric pressure jet-stirred reactor (JSR) during the ozonolysis of trans-2-butene to explore the kinetic pathways relevant to atmospheric chemistry and low-temperature combustion. The mole fraction profiles of reactants, intermediates, and final products are determined by means of molecular-beam mass spectrometry in conjunction with single-photon ionization employing tunable synchrotron-generated vacuum ultraviolet radiation. A network of CI reactions is identified in the temperature region below 600 K, characterized by CI addition to trans-2-butene, water, formaldehyde, formic acid, and methanol. No sequential additions of the CH₃CHOO CI are observed, in contrast with the reactivity of the simplest CI (H₂COO) and the earlier observation of an extensive reaction network with up to four H₂COO sequential additions (Phys. Chem. Chem. Phys., 2019, 21, 7341-7357). Experimental photoionization efficiency scans recorded at 300 K and 425 K and ab initio threshold energy calculations lead to the identification and quantification of previously elusive intermediates, such as ketohydroperoxide and hydroperoxide species. Specifically, the C₄H₈ + O₃ adduct is identified as a ketohydroperoxide (KHP, 3-hydroperoxybutan-2-one, CH₃C(O)CH(CH₃)OOH), while hydroxyacetaldehyde (glycolaldehyde, HCOCH₂OH) formation is attributed to unimolecular isomerization of the CIs. Other hydroperoxide species such as methyl hydroperoxide (CH₃OOH), ethyl hydroperoxide (C₂H₅OOH), butyl hydroperoxide (OOH), hydroperoxyl acetaldehyde (HOOCH₂CHO), hydroxyethyl hydroperoxide (CH₃CH(OH)OOH), but-1-enyl-3-hydroperoxide, and 4-hydroxy-3-methylpentan-2-one (HOCH(CH₃)CH(CH₃)C(O)CH₃) are also identified. Detection of additional oxygenated species such as methanol, ethanol, ketene, and aldehydes suggests multiple active oxidation routes. These results provide additional evidence that CIs are key intermediates of the ozone-unsaturated hydrocarbon reactions providing critical inputs for improved kinetics models.

Parsimonious Potential Energy Surface Expansions Using Dictionary Learning with Multipass Greedy Selection

Potential energy surfaces fit with basis set expansions have been shown to provide accurate representations of electronic energies and have enabled a variety of high-accuracy dynamics, kinetics, and spectroscopy applications. The number of terms in these expansions scales poorly with system size, a drawback that challenges their use for systems with more than ∼10 atoms. A solution is presented here using dictionary learning. Subsets of the full set of conventional basis functions are optimized using a newly developed multipass greedy regression method inspired by forward and backward selection methods from the statistics, signal processing, and machine learning literatures. The optimized representations have accuracies comparable to the full set but are 1 or more orders of magnitude smaller, and notably, the number of terms in the optimized multipass greedy expansions scales approximately linearly with the number of atoms.

Permutationally Invariant Polynomial Expansions with Unrestricted Complexity

A general strategy is presented for constructing and validating permutationally invariant polynomial (PIP) expansions for chemical systems of any stoichiometry. Demonstrations are made for three categories of gas-phase dynamics and kinetics: collisional energy-transfer trajectories for predicting pressure-dependent kinetics, three-body collisions for describing transient van der Waals adducts relevant to atmospheric chemistry, and nonthermal reactivity via quasiclassical trajectories. In total, 30 systems are considered with up to 15 atoms and 39 degrees of freedom. Permutational invariance is enforced in PIP expansions with as many as 13 million terms and 13 permutationally distinct atom types by taking advantage of petascale computational resources. The quality of the PIP expansions is demonstrated through the systematic convergence of in-sample and out-of-sample errors with respect to both the number of training data and the order of the expansion, and these errors are shown to predict errors in the dynamics for both reactive and nonreactive applications. The parallelized code distributed as part of this work enables the automation of PIP generation for complex systems with multiple channels and flexible user-defined symmetry constraints and for automatically removing unphysical unconnected terms from the basis set expansions, all of which are required for simulating complex reactive systems.

Watching a hydroperoxyalkyl radical (•QOOH) dissociate

A prototypical hydroperoxyalkyl radical (•QOOH) intermediate, transiently formed in the oxidation of volatile organic compounds, was directly observed through its infrared fingerprint and energy-dependent unimolecular decay to hydroxyl radical and cyclic ether products. Direct time-domain measurements of •QOOH unimolecular dissociation rates over a wide range of energies were found to be in accord with those predicted theoretically using state-of-the-art electronic structure characterizations of the transition state barrier region. Unimolecular decay was enhanced by substantial heavy-atom tunneling involving O-O elongation and C-C-O angle contraction along the reaction pathway. Master equation modeling yielded a fully a priori prediction of the pressure-dependent thermal unimolecular dissociation rates for the •QOOH intermediate-again increased by heavy-atom tunneling-which are required for global models of atmospheric and combustion chemistry.

On the Rate Constant for NH2+HO2 and Third-Body Collision Efficiencies for NH2+H(+M) and NH2+NH2(+M)

In low-temperature flash photolysis of NH₃/O₂/N₂ mixtures, the NH₂ consumption rate and the product distribution is controlled by the reactions NH₂ + HO₂ → products (R1), NH₂ + H (+M) → NH₃ (+M) (R2), and NH₂ + NH₂ (+M) → N₂H₄ (+M) (R3). In the present work, published flash photolysis experiments by, among others, Cheskis and co-workers, are re-interpreted using recent direct measurements of NH₂ + H (+N₂) and NH₂ + NH₂ (+N₂) from Altinay and Macdonald. To facilitate analysis of the FP data, relative third-body collision efficiencies compared to N₂ for R2 and R3 were calculated for O₂ and NH₃ as well as for other selected molecules. Results were in good agreement with the limited experimental data. Based on reported NH₂ decay rates in flash photolysis of NH₃/O₂/N₂, a rate constant for NH₂ + HO₂ → NH₃ + O₂ (R1a) of k_1a = 1.5(±0.5) × 10¹⁴ cm³ mol^-1 s^-1 at 295 K was derived. This value is higher than earlier determinations based on the FP results but in good agreement with recent theoretical work. Kinetic modeling of reported N₂O yields indicates that NH₂ + HO₂ → H₂NO + O (R1c) is competing with R1a, but perturbation experiments with addition of CH₄ indicate that it is not a dominating channel. Measured HNO profiles indicate that this component is formed directly by NH₂ + HO₂ → HNO + H₂O (R1b), but theoretical work indicates that R1b is only a minor channel. Based on this analysis, we estimate k_1c = 2.5 × 10¹³ cm³ mol^-1 s^-1 and k_1b = 2.5 × 10¹² cm³ mol^-1 s^-1 at 295 K, with significant uncertainty margins.

Extreme Low-Temperature Combustion Chemistry: Ozone-Initiated Oxidation of Methyl Hexanoate

The accelerating chemical effect of ozone addition on the oxidation chemistry of methyl hexanoate [CH₃(CH₂)₄C(═O)OCH₃] was investigated over a temperature range from 460 to 940 K. Using an externally heated jet-stirred reactor at p = 700 Torr (residence time τ = 1.3 s, stoichiometry φ = 0.5, 80% argon dilution), we explored the relevant chemical pathways by employing molecular-beam mass spectrometry with electron and single-photon ionization to trace the temperature dependencies of key intermediates, including many hydroperoxides. In the absence of ozone, reactivity is observed in the so-called low-temperature chemistry (LTC) regime between 550 and 700 K, which is governed by hydroperoxides formed from sequential O₂ addition and isomerization reactions. At temperatures above 700 K, we observed the negative temperature coefficient (NTC) regime, in which the reactivity decreases with increasing temperatures, until near 800 K, where the reactivity increases again. Upon addition of ozone (1000 ppm), the overall reactivity of the system is dramatically changed due to the time scale of ozone decomposition in comparison to fuel oxidation time scales of the mixtures at different temperatures. While the LTC regime seems to be only slightly affected by the addition of ozone with respect to the identity and quantity of the observed intermediates, we observed an increased reactivity in the intermediate NTC temperature range. Furthermore, we observed experimental evidence for an additional oxidation regime in the range near 500 K, herein referred to as the extreme low-temperature chemistry (ELTC) regime. Experimental evidence and theoretical rate constant calculations indicate that this ELTC regime is likely to be initiated by H abstraction from methyl hexanoate via O atoms, which originate from thermal O₃ decomposition. The theoretical calculations show that the rate constants for methyl ester initiation via abstraction by O atoms increase dramatically with the size of the methyl ester, suggesting that ELTC is likely not important for the smaller methyl esters. Experimental evidence is provided indicating that, similar to the LTC regime, the chemistry in the ELTC regime is dominated by hydroperoxide chemistry. However, mass spectra recorded at various reactor temperatures and at different photon energies provide experimental evidence of some differences in chemical species between the ELTC and the LTC temperature ranges.

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.

Isomer-Selective Detection of Keto-Hydroperoxides in the Low-Temperature Oxidation of Tetrahydrofuran

Keto-hydroperoxides (KHPs) are reactive, partially oxidized intermediates that play a central role in chain-branching reactions during the gas-phase low-temperature oxidation of hydrocarbons and oxygenated species. Although multiple isomeric forms of the KHP intermediate are possible in complex oxidation environments when multiple reactant radicals exist that contain nonequivalent O₂ addition sites, isomer-resolved data of KHPs have not been reported. In this work, we provide partially isomer-resolved detection and quantification of the KHPs that form during the low-temperature oxidation of tetrahydrofuran (THF, cycl.-O-CH₂CH₂CH₂CH₂-). We describe how these short-lived KHPs were detected, identified, and quantified using integrated experimental and theoretical approaches. The experimental approaches were based on direct molecular-beam sampling from a jet-stirred reactor operated at near-atmospheric pressure and at temperatures between 500 and 700 K, followed by mass spectrometry with single-photon ionization via tunable synchrotron-generated vacuum-ultraviolet radiation, and the identification of fragmentation patterns. The interpretation of the experiments was guided by theoretical calculations of ionization thresholds, fragment appearance energies, and photoionization cross sections. On the basis of the experimentally observed and theoretically calculated ionization and fragment appearance energies, KHP isomers could be distinguished as originating from H-abstraction reactions from either the α-C adjacent to the O atom or the β-C atoms. Temperature-dependent concentration profiles of the partially resolved isomeric KHP intermediates were determined in the range of 500-700 K, and the results indicate that the observed KHP isomers are formed overwhelmingly (∼99%) from the α-C THF radical. Comparisons of the partially isomer-resolved quantification of the KHPs to up-to-date kinetic modeling results reveal new opportunities for the development of a next-generation THF oxidation mechanism.

Nonthermal rate constants for CH4 * + X → CH3 + HX, X = H, O, OH, and O2

Quasiclassical trajectories are used to compute nonthermal rate constants, k^*, for abstraction reactions involving highly-excited methane CH₄ ^* and the radicals H, O, OH, and O₂. Several temperatures and internal energies of methane, E_vib, are considered, and significant nonthermal rate enhancements for large E_vib are found. Specifically, when CH₄ ^* is internally excited close to its dissociation threshold (E_vib ≈ D₀ = 104 kcal/mol), its reactivity with H, O, and OH is shown to be collision-rate-limited and to approach that of comparably-sized radicals, such as CH₃, with k^* > 10^-10 cm³ molecule^-1 s^-1. Rate constants this large are more typically associated with barrierless reactions, and at 1000 K, this represents a nonthermal rate enhancement, k^*/k, of more than two orders of magnitude relative to thermal rate constants k. We show that large nonthermal rate constants persist even after significant internal cooling, with k^*/k > 10 down to E_vib ≈ D₀/4. The competition between collisional cooling and nonthermal reactivity is studied using a simple model, and nonthermal reactions are shown to account for up to 35%-50% of the fate of the products of H + CH₃ = CH₄ ^* under conditions of practical relevance to combustion. Finally, the accuracy of an effective temperature model for estimating k^* from k is quantified.

Identification of the Criegee intermediate reaction network in ethylene ozonolysis: impact on energy conversion strategies and atmospheric chemistry

The reaction network of the simplest Criegee intermediate (CI) CH₂OO has been studied experimentally during the ozonolysis of ethylene.

Low-Temperature Oxidation of Ethylene by Ozone in a Jet-Stirred Reactor

Ethylene oxidation initiated by ozone addition (ozonolysis) is carried out in a jet-stirred reactor from 300 to 1000 K to explore the kinetic pathways relevant to low-temperature oxidation. The temperature dependencies of species' mole fractions are quantified using molecular-beam mass spectrometry with electron ionization and single-photon ionization employing tunable synchrotron-generated vacuum-ultraviolet radiation. Upon ozone addition, significant ethylene oxidation is found in the low-temperature regime from 300 to 600 K. Here, we provide new insights into the ethylene ozonolysis reaction network via identification and quantification of previously elusive intermediates by combining experimental photoionization energy scans and ab initio threshold energy calculations for isomer identification. Specifically, the C₂H₄ + O₃ adduct C₂H₄O₃ is identified as a keto-hydroperoxide (hydroperoxy-acetaldehyde, HOOCH₂CHO) based on the calculated and experimentally observed ionization energy of 9.80 (±0.05) eV. Quantification using a photoionization cross-section of 5 Mb at 10.5 eV results in 5 ppm at atmospheric conditions, which decreases monotonically with temperature until 550 K. Other hydroperoxide species that contribute in larger amounts to the low-temperature oxidation of C₂H₄, like H₂O₂, CH₃OOH, and C₂H₅OOH, are identified and their temperature-dependent mole fractions are reported. The experimental evidence for additional oxygenated species such as methanol, ketene, acetaldehyde, and hydroxy-acetaldehyde suggest multiple active oxidation routes. This experimental investigation closes the gap between ozonolysis at atmospheric and elevated temperature conditions and provides a database for future modeling.

Theoretical Study of the Ti–Cl Bond Cleavage Reaction in TiCl4

In this work the kinetics of the TiCl4 ⇌ TiCl3 + Cl reaction is studied theoretically. A variable-reaction coordinate transition-state theory (VRC-TST) is used to calculate the high-pressure limit rate coefficients. The interaction energy surface for the VRC-TST step is sampled directly at the CASPT2(6e,4o)/cc-pVDZ level of theory including an approximate treatment of the spin-orbit coupling. The pressure-dependence of the reaction in an argon bath gas is explored using the master equation in conjunction with the optimised VRC-TST transition-state number of states. The collisional energy transfer parameters for the TiCl4–Ar system are estimated via a “one-dimensional minimisation” method and classical trajectories. The Ti–Cl bond dissociation energy is computed using a complete basis set extrapolation technique with cc-pVQZ and cc-pV5Z basis sets. Good quantitative agreement between the estimated rate constants and available literature data is observed. However, the fall-off behaviour of the model results is not seen in the current experimental data. Sensitivity analysis shows that the fall-off effect is insensitive to the choice of model parameters and methods. More experimental work and development of higher-level theoretical methods are needed to further investigate this discrepancy.

Low Temperature Kinetics of the First Steps of Water Cluster Formation

We present a combined experimental and theoretical low temperature kinetic study of water cluster formation. Water cluster growth takes place in low temperature (23-69 K) supersonic flows. The observed kinetics of formation of water clusters are reproduced with a kinetic model based on theoretical predictions for the first steps of clusterization. The temperature- and pressure-dependent association and dissociation rate coefficients are predicted with an ab initio transition state theory based master equation approach over a wide range of temperatures (20-100 K) and pressures (10^{-6}-10  bar).

Pressure-dependent rate constants for PAH growth: formation of indene and its conversion to naphthalene

Unraveling the mechanisms for growth of polycyclic aromatic hydrocarbons (PAHs) requires accurate temperature- and pressure-dependent rate coefficients for a great variety of feasible pathways. Even the pathways for the formation of the simplest PAHs, indene and naphthalene, are fairly complex. These pathways provide important prototypes for modeling larger PAH growth. In this work we employ the ab initio RRKM theory-based master equation approach to predict the rate constants involved in the formation of indene and its conversion to naphthalene. The reactions eventually leading to indene involve C₉H_x (x = 8–11) potential energy surfaces (PESs) and include C₆H₅ + C₃H₄ (allene and propyne), C₆H₆ + C₃H₃, benzyl + C₂H₂, C₆H₅ + C₃H₆, C₆H₆ + C₃H₅ and C₆H₅ + C₃H₅. These predictions allow us to make a number of valuable observations on the role of various mechanisms. For instance, we demonstrate that reactions which can significantly contribute to the formation of indene include phenyl + allene and H-assisted isomerization to indene of its major product, 3-phenylpropyne, benzyl + acetylene, and the reactions of the phenyl radical with propene and the allyl radical, both proceeding via the 3-phenylpropene intermediate. 3-Phenylpropene can be activated to a 1-phenylallyl radical, which in turn rapidly decomposes to indene. Next, indene can be converted to benzofulvene or naphthalene under typical combustion conditions, via its activation by H atom abstraction and methyl substitution on the five-membered ring followed by isomerization and decomposition of the resulting 1-methylindenyl radical, C₁₀H₉ → C₁₀H₈ + H. Alternatively, the same region of the C₁₀H₉ PES can be accessed through the reaction of benzyl with propargyl, C₇H₇ + C₃H₃ → C₁₀H₁₀ → C₁₀H₉ + H, which therefore can also contribute to the formation of benzofulvene or naphthalene. Benzofulvene easily transforms to naphthalene by H-assisted isomerization. An analysis of the effect of pressure on the reaction outcome and relative product yields is given, and modified Arrhenius fits of the rate constants are reported for the majority of the considered reactions. Ultimately, the implementation of such expressions in detailed kinetic models will help quantify the role of these reactions for PAH growth in various environments.

Kinetics of Propargyl Radical Dissociation

Due to the prominent role of the propargyl radical for hydrocarbon growth within combustion environments, it is important to understand the kinetics of its formation and loss. The ab initio transition state theory-based master equation method is used to obtain theoretical kinetic predictions for the temperature and pressure dependence of the thermal decomposition of propargyl, which may be its primary loss channel under some conditions. The potential energy surface for the decomposition of propargyl is first mapped at a high level of theory with a combination of coupled cluster and multireference perturbation calculations. Variational transition state theory is then used to predict the microcanonical rate coefficients, which are subsequently implemented within the multiple-well multiple-channel master equation. A variety of energy transfer parameters are considered, and the sensitivity of the thermal rate predictions to these parameters is explored. The predictions for the thermal decomposition rate coefficient are found to be in good agreement with the limited experimental data. Modified Arrhenius representations of the rate constants are reported for utility in combustion modeling.

Detection and Identification of the Keto-Hydroperoxide (HOOCH2OCHO) and Other Intermediates during Low-Temperature Oxidation of Dimethyl Ether

In this paper we report the detection and identification of the keto-hydroperoxide (hydroperoxymethyl formate, HPMF, HOOCH2OCHO) and other partially oxidized intermediate species arising from the low-temperature (540 K) oxidation of dimethyl ether (DME). These observations were made possible by coupling a jet-stirred reactor with molecular-beam sampling capabilities, operated near atmospheric pressure, to a reflectron time-of-flight mass spectrometer that employs single-photon ionization via tunable synchrotron-generated vacuum-ultraviolet radiation. On the basis of experimentally observed ionization thresholds and fragmentation appearance energies, interpreted with the aid of ab initio calculations, we have identified HPMF and its conceivable decomposition products HC(O)O(O)CH (formic acid anhydride), HC(O)OOH (performic acid), and HOC(O)OH (carbonic acid). Other intermediates that were detected and identified include HC(O)OCH3 (methyl formate), cycl-CH2-O-CH2-O- (1,3-dioxetane), CH3OOH (methyl hydroperoxide), HC(O)OH (formic acid), and H2O2 (hydrogen peroxide). We show that the theoretical characterization of multiple conformeric structures of some intermediates is required when interpreting the experimentally observed ionization thresholds, and a simple method is presented for estimating the importance of multiple conformers at the estimated temperature (∼100 K) of the present molecular beam. We also discuss possible formation pathways of the detected species: for example, supported by potential energy surface calculations, we show that performic acid may be a minor channel of the O2 + ĊH2OCH2OOH reaction, resulting from the decomposition of the HOOCH2OĊHOOH intermediate, which predominantly leads to the HPMF.

Ab initio Variational Transition State Theory and Master Equation Study of the Reaction (OH)3SiOCH2 + CH3 ⇌ (OH)3SiOC2H5

In this paper we use variable reaction coordinate variational transition state theory (VRC-TST) to calculate the reaction rate constants for the two reactions, R1: (OH)3SiOCH2 + CH3 ⇌ (OH)3SiOC2H5, and R2: CH2OH + CH3 ⇌ C2H5OH. The first reaction is an important channel during the thermal decomposition of tetraethoxysilane (TEOS), and its rate coefficient is the main focus of this work. The second reaction is analogous to the first and is used as a basis for comparison. The interaction energies are obtained on-the-fly at the CASPT2(2e,2o)/cc-pVDZ level of theory. A one-dimensional correction to the sampled energies was introduced to account for the energetic effects of geometry relaxation along the reaction path. The computed, high-pressure rate coefficients were calculated to be, R1: k 1 = 2.406 × 10−10 T −0.301 exp (− 271.4/T) cm3 molecule –1 s –1 and R2: k 2 = 1.316 × 10−10 T −0.189 exp (− 256.5/T) cm3 molecule –1 s –1. These rates differ from each other by only 10% – 30% over the temperature range 300–2000 K. A comparison of the computed rates with experimental data shows good agreement and an improvement over previous results. The pressure dependency of the reaction R1 is explored by solving a master equation using helium as a bath gas. The results obtained show that the reaction is only weakly pressure dependent over the temperature range 300–1700 K, with the predicted rate constant being within 50% of its high-pressure limit at atmospheric pressure.

Predictive a priori pressure-dependent kinetics

The ability to predict the pressure dependence of chemical reaction rates would be a great boon to kinetic modeling of processes such as combustion and atmospheric chemistry. This pressure dependence is intimately related to the rate of collision-induced transitions in energy E and angular momentum J. We present a scheme for predicting this pressure dependence based on coupling trajectory-based determinations of moments of the E,J-resolved collisional transfer rates with the two-dimensional master equation. This completely a priori procedure provides a means for proceeding beyond the empiricism of prior work. The requisite microcanonical dissociation rates are obtained from ab initio transition state theory. Predictions for the CH4 = CH3 + H and C2H3 = C2H2 + H reaction systems are in excellent agreement with experiment.

Separability of tight and roaming pathways to molecular decomposition

Recent studies have questioned the separability of the tight and roaming mechanisms to molecular decomposition. We explore this issue for a variety of reactions including MgH(2) → Mg + H(2), NCN → CNN, H(2)CO → H(2) + CO, CH(3)CHO → CH(4) + CO, and HNNOH → N(2) + H(2)O. Our analysis focuses on the role of second-order saddle points in defining global dividing surfaces that encompass both tight and roaming first-order saddle points. The second-order saddle points define an energetic criterion for separability of the two mechanisms. Furthermore, plots of the differential contribution to the reactive flux along paths connecting the first- and second-order saddle points provide a dynamic criterion for separability. The minimum in the differential reactive flux in the neighborhood of the second-order saddle point plays the role of a mechanism divider, with the presence of a strong minimum indicating that the roaming and tight mechanisms are dynamically distinct. We show that the mechanism divider is often, but not always, associated with a second-order saddle point. For the formaldehyde and acetaldehyde reactions, we find that the minimum energy geometry on a conical intersection is associated with the mechanism divider for the tight and roaming processes. For HNNOH, we again find that the roaming and tight processes are dynamically separable but we find no intrinsic feature of the potential energy surface associated with the mechanism divider. Overall, our calculations suggest that roaming and tight mechanisms are generally separable over broad ranges of energy covering most kinetically relevant regimes.

Identification of Tetrahydrofuran Reaction Pathways in Premixed Flames

Premixed low-pressure tetrahydrofuran/oxygen/argon flames are investigated by photoionization molecular-beam mass spectrometry using vacuum-ultraviolet synchrotron radiation. For two equivalence ratios (φ = 1.00 and 1.75), mole fractions are measured as a function of distance from the burner for almost 60 intermediates with molar masses ranging from 2 (H2) to 88 (C4H6O2), providing a broad database for flame modeling studies. The isomeric composition is resolved by comparisons between experimental photoionization efficiency data and theoretical simulations, based on calculated ionization energies and Franck-Condon factors. Special emphasis is put on the resolution of the first reaction steps in the fuel destruction. The photoionization experiments are complemented by electron-ionization molecular-beam mass-spectrometry measurements that provide data with high mass resolution. For three additional flames with intermediate equivalence ratios (φ = 1.20, 1.40 and 1.60), mole fractions of major species and photoionization efficiency spectra of intermediate species are reported, extending the database for the development of chemical kinetic models.