Toward an Accurate Black-Box Tool for the Kinetics of Gas-Phase Reactions Involving Barrier-less Elementary Steps

An enhanced computational protocol has been devised for the accurate characterization of gas-phase barrier-less reactions in the framework of the reaction-path (RP) and variable reaction coordinate variational transition-state theory. In particular, the synergistic combination of density functional theory and Monte Carlo sampling to optimize reactive fluxes led to a reliable yet effective computational workflow. A black-box strategy has been developed for selecting the most suited density functional with reference to a high-level one-dimensional reference potential. At the same time, different descriptions of hindered rotations are automatically selected, depending on the corresponding harmonic frequencies along the RP. The performance of the new tool is investigated by means of two prototypical reactions involving different degrees of static and dynamic correlation, namely, H2S + Cl and CH3 + CH3. The remarkable agreement of the computed kinetic parameters with the available experimental data confirms the accuracy and robustness of the proposed approach. Together with their intrinsic interest, these results also pave the way toward systematic investigations of gas-phase reactions involving barrier-less elementary steps by a reliable, user-friendly tool, which can be confidently used also by nonspecialists.

Reaction Profiles and Kinetics for Radical-Radical Hydrogen Abstraction via Multireference Coupled Cluster Theory

Radical-radical abstractions in hydrocarbon oxidation chemistry are disproportionation reactions that are generally exothermic with little or no barrier yet are underappreciated and poorly studied. Such challenging multireference electronic structure problems are tackled here using the recently developed state-specific multireference coupled cluster methods Mk-MRCCSD and Mk-MRCCSD(T), as well as the companion perturbation theory Mk-MRPT2 and the established MRCISD, MRCISD+Q, and CASPT2 approaches. Reaction paths are investigated for five prototypes involving radical-radical hydrogen abstraction: H + BeH → H₂+ Be, H + NH₂ → H₂ + NH, CH₃ + C₂H₅ → CH₄ + C₂H₄, H + C₂H₅ → H₂ + C₂H₄, and H + HCO → H₂ + CO. Full configuration interaction (FCI) benchmark computations for the H + BeH, H + NH₂, and H + HCO reactions prove that Mk-MRCCSD(T) provides superior accuracy for the interaction energies in the entrance channel, with mean absolute errors less than 0.3 kcal mol^-1 and percentage deviations less than 10% over the fragment separations of relevance to kinetics. To facilitate combustion studies, energetics for the CH₃ + C₂H₅, H + C₂H₅, and H + HCO reactions were computed at each level of theory with correlation-consistent basis sets (cc-pVXZ, X = T, Q, 5) and extrapolated to the complete basis set (CBS) limit. These CBS energies were coupled with CASPT2 projected vibrational frequencies along a minimum energy path to obtain rate constants for these three reactions. The rigorous Mk-MRCCSD(T)/CBS results demonstrate unequivocally that these three reactions proceed with no barrier in the entrance channel, contrary to some earlier predictions. Mk-MRCCSD(T) also reveals that the economical CASPT2 method performs well for large interfragment separations but may deteriorate substantially at shorter distances.

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.

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).

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.

Theory, measurements, and modeling of OH and HO2 formation in the reaction of cyclohexyl radicals with O2

The production of OH and HO(2) in Cl-initiated oxidation of cyclohexane has been measured using pulsed-laser photolytic initiation and continuous-laser absorption detection. The experimental data are modeled by master equation calculations that employ new G2(MP2)-like ab initio characterizations of important stationary points on the cyclo-C(6)H(11)O(2) surface. These ab initio calculations are a substantial expansion on previously published characterizations, including explicit consideration of conformational changes (chair-boat, axial-equatorial) and torsional potentials. The rate constants for the decomposition and ring-opening of cyclohexyl radical are also computed with ab initio based transition state theory calculations. Comparison of kinetic simulations based on the master equation results with the present experimental data and with literature determinations of branching fractions suggests adjustment of several transition state energies below their ab initio values. Simulations with the adjusted values agree well with the body of experimental data. The results once again emphasize the importance of both direct and indirect components of the kinetics for the production of both HO(2) and OH in radical + O(2) reactions.

Long-range transition state theory

The implementation of variational transition state theory (VTST) for long-range asymptotic potential forms is considered, with particular emphasis on the energy and total angular momentum resolved (microJ-VTST) implementation. A long-range transition state approximation yields a remarkably simple and universal description of the kinetics of reactions governed by long-range interactions. The resulting (microJ-VTST) implementation is shown to yield capture-rate coefficients that compare favorably with those from trajectory simulations (deviating by less than 10%) for a wide variety of neutral and ionic long-range potential forms. Simple analytic results are derived for many of these cases. A brief comparison with a variety of low-temperature experimental studies illustrates the power of this approach as an analysis tool. The present VTST approach allows for a simple analysis of the applicability conditions for some related theoretical approaches. It also provides an estimate of the temperature or energy at which the "long-range transition state" moves to such short separations that short-range effects, such as chemical bonding, steric repulsion, and electronic state selectivity, must be considered.

Uncovering the Fundamental Chemistry of Alkyl + O2Reactions via Measurements of Product Formation

The reactions of alkyl radicals (R) with molecular oxygen (O(2)) are critical components in chemical models of tropospheric chemistry, hydrocarbon flames, and autoignition phenomena. The fundamental kinetics of the R + O(2) reactions is governed by a rich interplay of elementary physical chemistry processes. At low temperatures and moderate pressures, the reactions form stabilized alkylperoxy radicals (RO(2)), which are key chain carriers in the atmospheric oxidation of hydrocarbons. At higher temperatures, thermal dissociation of the alkylperoxy radicals becomes more rapid and the formation of hydroperoxyl radicals (HO(2)) and the conjugate alkenes begins to dominate the reaction. Internal isomerization of the RO(2) radicals to produce hydroperoxyalkyl radicals, often denoted by QOOH, leads to the production of OH and cyclic ether products. More crucially for combustion chemistry, reactions of the ephemeral QOOH species are also thought to be the key to chain branching in autoignition chemistry. Over the past decade, the understanding of these important reactions has changed greatly. A recognition, arising from classical kinetics experiments but firmly established by recent high-level theoretical studies, that HO(2) elimination occurs directly from an alkylperoxy radical without intervening isomerization has helped resolve tenacious controversies regarding HO(2) formation in these reactions. Second, the importance of including formally direct chemical activation pathways, especially for the formation of products but also for the formation of the QOOH species, in kinetic modeling of R + O(2) chemistry has been demonstrated. In addition, it appears that the crucial rate coefficient for the isomerization of RO(2) radicals to QOOH may be significantly larger than previously thought. These reinterpretations of this class of reactions have been supported by comparison of detailed theoretical calculations to new experimental results that monitor the formation of products of hydrocarbon radical oxidation following a pulsed-photolytic initiation. In this article, these recent experiments are discussed and their contributions to improving general models of alkyl + O(2) reactions are highlighted. Finally, several prospects are discussed for extending the experimental investigations to the pivotal questions of QOOH radical chemistry.

Predictive theory for the combination kinetics of two alkyl radicals

An ab initio transition state theory based procedure for accurately predicting the combination kinetics of two alkyl radicals is described. This procedure employs direct evaluations of the orientation dependent interaction energies at the CASPT2/cc-pvdz level within variable reaction coordinate transition state theory (VRC-TST). One-dimensional corrections to these energies are obtained from CAS+1+2/aug-cc-pvtz calculations for CH3 + CH3 along its combination reaction path. Direct CAS+1+2/aug-cc-pvtz calculations demonstrate that, at least for the purpose of predicting the kinetics, the corrected CASPT2/cc-pvdz potential energy surface is an accurate approximation to the CAS+1+2/aug-cc-pvtz surface. Furthermore, direct trajectory simulations, performed at the B3LYP/6-31G* level, indicate that there is little local recrossing of the optimal VRC transition state dividing surface. The corrected CASPT2/cc-pvdz potential is employed in obtaining direct VRC-TST kinetic predictions for the self and cross combinations of methyl, ethyl, iso-propyl, and tert-butyl radicals. Comparisons with experiment suggest that the present dynamically corrected VRC-TST approach provides quantitatively accurate predictions for the capture rate. Each additional methyl substituent adjacent to a radical site is found to reduce the rate coefficient by about a factor of two. In each instance, the rate coefficients are predicted to decrease quite substantially with increasing temperature, with the more sterically hindered reactants having a more rapid decrease. The simple geometric mean rule, relating the capture rate for the cross reaction to those for the self-reactions, is in remarkably good agreement with the more detailed predictions. With suitable generalizations the present approach should be applicable to a wide array of radical-radical combination reactions.

A Saturated LIF Study on the High Pressure Limiting Rate Constant of the Reaction CN + NO + M → NCNO + M between 200 and 600 K

The reaction CN + NO + M ↔ NCNO + M was studied in the bath gas helium at temperatures between 200 and 600 K and in the pressure range between 1 and 100 bar. CN radicals were generated by laser flash photolysis of BrCN at 193 nm in the presence of NO and high helium pressures. The concentration of CN radicals was detected by recording their non-resonant fluorescence yield at 420 nm after delayed excitation in the (0, 1)-band of the X