Degradation of Carbonyl Hydroperoxides in the Atmosphere and in Combustion

Identification of Torsional Modes in Complex Molecules Using Redundant Internal Coordinates: The Multistructural Method with Torsional Anharmonicity with a Coupled Torsional Potential and Delocalized Torsions

Identification of internal-rotation modes in the normal-mode analysis of complex molecules is important for accurately describing the thermodynamic properties and kinetics of complex molecules when it is necessary to treat the anharmonicity of torsions and the multiconformer anharmonicity caused by the internal rotations. However, identifying and distinguishing torsional modes are very challenging because they are coupled to one another. In this work, we present a new strategy to automatically identify torsional vibrations and separate them from the other vibrational modes. By combining a redundant-internal-coordinate auto-generation procedure with torsional projection techniques, we automate the procedure of identifying and separating the coupled torsions, and we show that we can obtain robust and consistent results with various reasonable definitions of redundant-internal-coordinate sets. This model has been implemented in a new development version of the MSTor program to reduce the user input needed for multistructural and torsional anharmonicity (MS-T) calculations. The new method is called multistructural and torsional anharmonicity with a coupled torsional potential and delocalized torsions ([MS-T(CD)]. As example applications, we consider MS-T(CD) calculations on three molecules (2-hexyl radical, n-propylbenzene, and 5-hydroperoxy-6-oxohexanoylperoxy radical) that have multiple rotors and that provide challenges to choosing good sets of nonredundant-internal coordinates, and we compare the performance of the new strategy to five other torsion identification methods. The new strategy is demonstrated to be efficient in separating the torsional and nontorsional elements in the Hessian matrix, as well as in providing reasonable projected nontorsional frequencies to be used for calculations of partition function and thermochemistry.

Lowering of Reaction Rates by Energetically Favorable Hydrogen Bonding in the Transition State. Degradation of Biofuel Ketohydroperoxides by OH

Ketohydroperoxides (KHPs) are oxygenates with carbonyl and hydroperoxy functional groups, and they are generated under combustion and atmospheric conditions. Their fate is crucial for secondary organic aerosol formation in the troposphere and for the ignition processes of biofuels in advanced combustion engines. We investigated the thermodynamics and kinetics of nine hydrogen abstraction reactions from four ether KHPs by OH. We find that the rate constants are strongly affected by entropic effects whose estimation requires a consideration of higher-energy conformers of the transition state. A density functional was selected for these reactions by comparison to coupled cluster calculations, and it was used for calculations by multistructural canonical transition-state theory with multidimensional tunneling over the temperature range of 200-2000 K. We find that the effect of multistructural torsional anharmonicity is very large and quite different for the various ether KHP reactions. A leading cause of the structural dependence is the dominance of entropic factors due to the lack of hydrogen bonding in some of the higher-energy conformers of the transition states. Four of the reactions involve abstraction from the α-carbon (the carbon vicinal to the hydroperoxide group); they exhibit nonmonotonic temperature dependence with complex fuel-specific dependence. The rate constants for abstraction from a non-α-carbon of a given KHP can be faster than the ones for abstraction from an α-carbon; in two cases, this is due to entropy, and in one case, the non-α-carbon abstraction has a lower energy barrier. Tunneling and recrossing effects are also found to be important.

Effects of Anharmonicity, Recrossing, Tunneling, and Pressure on the H-Abstractions from Dimethylamine by Triplets O and O2

Rate constants of the H-abstraction reactions from dimethylamine (DMA) by triplets O and O₂ are theoretically determined with the canonical variational transition-state theory (CVT). By comparing the barrier heights and reaction energies obtained from different density-functional theory methods to those computed from the gold-standard method CCSD(T)/CBS(T-Q), we identify the M08-HX/ma-TZVP method as the best with a mean unsigned deviation of 1.0 kcal mol^-1. On the basis of the optimized geometries and frequencies with the selected method, the rate constants are calculated using the CVT method combined with the multistructural torsional anharmonicity and small-curvature tunnelling (MS-CVT/SCT) options in the temperature range 200-2000 K. The calculations show that OH and HO₂ are mainly produced from the direct abstraction from the C-H bond. The multistructural torsional anharmonicity has a large contribution to the rate constants, and the effects of recrossing and tunneling at the N-site are more important than those at C-site. Additionally, given the formation of reactant complex between DMA and triplet O, the H-abstraction channel is not favored at high pressure. Our calculations with both the Polyrate and MESS codes agree with the reported data within the uncertainty.

Develop the low-temperature oxidation mechanism of cyclopentane: an experimental and theoretical study

At present, the reactivity of cyclic alkanes is estimated by analogy comparison with non-cyclic hydrocarbons. Due to the difference in the structure of cycloalkanes and non-cycloalkanes, the thermodynamic data obtained by analogy is not applicable. In this study, a molecular beam sampling vacuum ultraviolet photoionization time-of-flight mass spectrometer (MB-VUV-PI-TOFMS) was applied to study the low-temperature oxidation of cyclopentane (CPT) at a total pressure range from 1-3 atm and low-temperature range between 500 K and 800 K. The low-temperature reaction products including cyclic olefins, cyclic ethers, highly oxygenated intermediates (e.g., ketohydroperoxide KHP, keto-dihydroperoxide KDHP, olefinic hydroperoxides OHP and ketone structure products) were observed. Further investigation of the oxidation of CPT, the electronic structure calculations at the UCCSD(T)-F12a/aug-cc-pVDZ//B3LYP/6-31+ G(d,p) level were carried out to explore the reactivity of O2 molecules adding sequentially to cyclopentyl radicals. Experimental and theoretical observations showed that the dominant product channel in the reaction of CPT radicals with O2 is HO2 elimination, yielding cyclopentene. The pathways of second and third O2 addition, the dissociation of hydroperoxide were further confirmed. The results of this study will develop the low-temperature oxidation mechanism of CPT, which can be used for future research on accurately simulating the combustion process of CPT.

Differences in the torsional anharmonicity between reactant and transition state: the case of 3-butenal + H abstraction reactions

Thermal rate coefficients for the hydrogen-abstraction reactions of 3-butenal by a hydrogen atom were obtained applying multipath canonical variational theory with small-curvature tunneling (MP-CVT/SCT). Torsional anharmonicity due to the hindered rotors was taken into account by calculating the rovibrational partition function using the extended two-dimensional torsional (E2DT) method. For comparison, rovibrational partition functions were also estimated using the multistructural method with torsional anharmonicity based on a coupled torsional potential (MS-T(C)). By contrast, with (E)-2-butenal reactions, the abstraction reactions of 3-butenal proceed via five reaction channels (R1)-(R5). In a conformational search, 45 distinguishable structures of transition states were found, including enantiomers, which were separated into six conformational reaction channels (CRCs). The individual reactive paths were constructed, the recrossing and semiclassical transmission coefficients estimated, and the multipath rate constants were obtained. High torsional barriers between the wells of CRC2/CRC6 indicate a harmonic behavior. Consequently, a difference between the torsional anharmonicity of 3-butenal and the transition states is responsible for the increase in the thermal rate constants for channel (R2). Analysis of the contributions of each conformer of the transition state shows an important contribution of the high-energy rotamers in the total flux of (R1)-(R5). After fitting the rate constants in a four-parameter equation, the activation energy estimation showed a strong temperature dependence.

Multistructural Variational Reaction Kinetics of the Simplest Unsaturated Methyl Ester: H-Abstraction from Methyl Acrylate by H, OH, CH3, and HO2 Radicals

The H-abstraction reaction kinetics of methyl acrylate (MA) + H/OH/CH₃/HO₂ radicals have been investigated theoretically in the present work. For these reactions, the reaction energies and barrier heights are first computed using several density functionals and compared to the coupled cluster CCSD(T)-F12/jun-cc-pVTZ benchmark calculations. The M062X/maug-cc-pVTZ method shows the best performance with the smallest mean unsigned deviation (MUD) of 0.42 kcal mol^-1. Combined with the electronic structure calculations using the M062X/maug-cc-pVTZ method, the multistructural canonical variational transition-state theory (MS-CVT) with small-curvature tunneling (SCT) is employed to calculate the reaction rate constants at 500-2000 K. The variational effect is between 0.56 and 1.0, the multistructural torsional anharmonicity factor ranges from 0.004 to 4.57, and the tunneling coefficient is in the range of 0.5-4.70. Notably, given the existence of reactant complexes (RCs) between reactants and transition states for the reaction systems MA + OH/HO₂, we further compare the rate constants under the low-pressure limit (LPL) kinetic model, which treats the reaction as a single-step process and neglects RCs, and the pre-equilibrium model, which takes RCs into account in the reaction and treats the reaction as a two-step process. The rate constants calculated by these two models are similar within the combustion temperature range, and apparent differences occur at lower temperatures. In addition, we determine the branching ratios as a function of temperature and find that the methyl site (S3) abstractions by OH and H radicals are dominant in the low- and high-temperature ranges, respectively. Moreover, we update the kinetic model with the calculated H-abstraction rate constants to simulate the ignition delay times of MA. The simulations of the updated model are in good agreement with experimental results. The accurate reaction kinetics determined in this work are useful for the understanding and prediction of consumption branching fractions and ignition properties of the unsaturated methyl esters.

Rate coefficients and product branching ratios for (E)-2-butenal + H reactions

Thermal rate constants for the hydrogen abstraction reactions of (E)-2-butenal by hydrogen atoms were calculated, for the first time, using the multipath canonical variational theory with small-curvature tunneling (MP-CVT/SCT). After a torsional potential energy surface exploration, ten conformations of the transition states (including the mirror images) were found and separated into four conformational reaction channels (CRCs). Individual energy paths of each CRC were built, recrossing and quantum tunneling effects estimated, and the thermal rate constants obtained. Due to the hindered rotors, the torsional anharmonicity was incorporated in the rate coefficient through the calculations of the rovibrational partition functions using the extended two-dimensional torsional method (E2DT). For comparison, the one-well (1W-CVT/SCT) and harmonic multipath (MP-CVT/SCT) thermal rate constants were also estimated. In addition, kinetic Monte Carlo (KMC) simulations were performed to predict the product branching ratios. For all kinetic approaches, the formation of products of (R1) is predominant. Compared to the harmonic multipath estimation, the percentage of reaction (R4) increases by approximately 9% when the torsional anharmonicity is taken into account. For the reactions (R2) and (R3), the product branching ratio is slightly decreased when compared with the harmonic simulation.

Multistructural Anharmonicity Controls the Radical Generation Process in Biofuel Combustion

The OH radical plays an important role in combustion, and isopentanol (3-methylbutan-1-ol) is a promising sustainable fuel additive and second-generation biofuel. The abstractions of H atoms from fuel molecules are key initiation steps for chain branching in combustion chemistry. In comparison with the more frequently studied ethanol, isopentanol has a longer carbon chain that allows a greater number of products, and experimental work is unavailable for the branching fractions to the various products. However, the site-dependent kinetics of isopentanol with OH radicals are usually experimentally unavailable. Alcohol oxidation by OH is also important in the atmosphere, and in the present study we calculate the rate constants and branching fractions of the hydrogen abstraction reaction of isopentanol by OH radical in a broad temperature range of 298-2400 K, covering temperatures important for atmospheric chemistry and those important for combustion. The calculations are done by multipath variational transition state theory (MP-VTST). With a combination of electronic structure calculations, we determine previously missing thermochemical data. With MP-VTST, a multidimensional tunneling approximation, multiple-structure anharmonicity, and torsional potential anharmonicity, we carried out more realistic rate constant calculations than can be computed by conventional single-structure harmonic transition state theory or by the empirical relations that are currently used in atmospheric and combustion modeling. The roles of various factors in determining the rates are elucidated, and we show that recrossing, tunneling, and multiple structures are all essential for accurate work. We conclude that the multiple structure anharmonicity is the most important correction to conventional transition state theory for this reaction, although recrossing effects and tunneling are by no means insignificant and the tunneling depends significantly on the path. The thermodynamic and kinetics data determined in this work are indispensable for the gas-phase degradation of alcohols in the atmosphere and for the detailed understanding and prediction of ignition mechanisms of biofuels in combustion.