Long-range transition state theory

Reassessment of the High-Temperature Oxidation of Di Ethyl Ether through Ab Initio Calculations

Recent investigations of diethyl ether (DEE) high-temperature pyrolysis and fuel-rich oxidation have highlighted the failure of existing kinetic models to describe experimental CO production. The DEE high-temperature pyrolysis and oxidation chemistry is thus investigated through ab initio calculations. Geometries, frequencies, and hindered-rotor potentials of reactants, products, and transition states of key reactions (fuel decomposition radical decomposition and H-abstraction reactions) are calculated with the B2PLYP-D3/def2-TZVPD method, whereas final energies are refined using CCSD(T)/aug-cc-pV(D,T)Z. Temperature- and pressure-dependent rate constants are then derived from either canonical transition state theory (CTST) or ME/RRKM analysis with the inclusion of tunneling effect and hindered-rotor corrections and compared to experimental measurements when available as well as to previously suggested values. This new information is then merged with a C0-C3 core chemistry model and a low-temperature chemistry DEE subset from the literature to propose a new kinetic model for the combustion of DEE. This model is tested successfully against a large database related to the high-temperature oxidation and pyrolysis chemistry of DEE, including ignition delay times, shock tube speciation data, time-resolved CO profiles, laminar flame speeds, flame structures, and jet-stirred reactor data.

Insights into the role of the H-abstraction reaction kinetics of amines in understanding their degeneration fates under atmospheric and combustion conditions

Amines, a class of prototypical volatile organic compounds, have garnered considerable interest within the context of atmospheric and combustion chemistry due to their substantial contributions to the formation of hazardous pollutants in the atmosphere. In the current energy landscape, the implementation of carbon-neutral energy and strategic initiatives leads to generation of new amine sources that cannot be overlooked in terms of the emission scale. To reduce the emission level of amines from their sources and mitigate their impact on the formation of harmful substances, a comprehensive understanding of the fundamental reaction kinetics during the degeneration process of amines is imperative. This perspective article first presents an overview of both traditional amine sources and emerging amine sources within the context of carbon peaking and carbon neutrality and then highlights the importance of H-abstraction reactions in understanding the atmospheric and combustion chemistry of amines from the perspective of reaction kinetics. Subsequently, the current experimental and theoretical techniques for investigating the H-abstraction reactions of amines are introduced, and a concise summary of research endeavors made in this field over the past few decades is provided. In order to provide accurate kinetic parameters of the H-abstraction reactions of amines, advanced kinetic calculations are performed using the multi-path canonical variational theory combined with the small-curvature tunneling and specific-reaction parameter methods. By comparing with the literature data, current kinetic calculations are comprehensively evaluated, and these validated data are valuable for the development of the reaction mechanism of amines.

The Interplay between the Temperature and Pressure on the Reaction Pathways of the Prenol Oxidation by Hydroxyl Radicals

Despite prenol emerging as a next-generation biofuel, some questions about its mechanism still need to be adequately proposed to rationalize its consumption and evaluate its efficiency in spark-ignition (SI) engines. Here, we present new insights into the reaction mechanism of prenol (3-methyl-2-buten-1-ol) with OH radicals as a function of temperature and pressure. We have determined that the different temperature and pressure conditions control the preferred products. At combustion temperatures and low pressures, OH-addition adducts are suppressed, increasing the formation of α and δ allylic radicals responsible for the auto-ignition.

New insights into the mechanism and kinetics of the addition reaction of unsaturated Criegee intermediates to CF3COOH and tropospheric implications

In this work, we have investigated the mechanism, thermochemistry and kinetics of the reaction of syn-cis-CH2RzCRyCO+O- (where Rz, Ry = H, CH3-) unsaturated Criegee intermediates (CIs) with CF3COOH using quantum chemical methods. The rate coefficients for the barrierless reactions were calculated using variable reaction coordinate variational transition state theory (VRC-VTST). For the syn-cis-CH2RzCRyCO+O- conformation in which conjugated CC and CO double bonds are aligned with each other, we propose a new pathway for the unidirectional addition of an OC-OH molecule (CF3COOH) to the CC double bond of syn-cis-CH2RzCRyCO+O-. The rate coefficient for the 1,4-CC addition reaction at 298 K is ∼10-10 to 10-11 cm3 s-1, resulting in the formation of CF3C(O)OCH2CRzRyCOOH trifluoroacetate alkyl allyl hydroperoxide (TFAAAH) as a new transitory adduct. It can act as a precursor for the formation of secondary organic aerosols (SOAs). This novel TFAAAH hydroperoxide was identified through a detailed quantum chemical study of the 1,4-addition mechanism and will provide new insights into the significance of the 1,4-addition reaction of unsaturated Cls with trace tropospheric gases on -CRzCH2 vinyl carbon atoms.

Experimental and Theoretical Study of the OH-Initiated Degradation of Piperidine under Simulated Atmospheric Conditions

The OH-initiated photo-oxidation of piperidine and the photolysis of 1-nitrosopiperidine were investigated in a large atmospheric simulation chamber and in theoretical calculations based on CCSD(T*)-F12a/aug-cc-pVTZ//M062X/aug-cc-pVTZ quantum chemistry results and master equation modeling of the pivotal reaction steps. The rate coefficient for the reaction of piperidine with OH radicals was determined by the relative rate method to be kOH-piperidine = (1.19 ± 0.27) × 10-10 cm3 molecule-1 s-1 at 304 ± 2 K and 1014 ± 2 hPa. Product studies show the piperidine + OH reaction to proceed via H-abstraction from both CH2 and NH groups, resulting in the formation of the corresponding imine (2,3,4,5-tetrahydropyridine) as the major product and in the nitramine (1-nitropiperidine) and nitrosamine (1-nitrosopiperidine) as minor products. Analysis of 1-nitrosopiperidine photolysis experiments under natural sunlight conditions gave the relative rates jrel = j1-nitrosoperidine/jNO2 = 0.342 ± 0.007, k3/k4a = 0.53 ± 0.05 and k2/k4a = (7.66 ± 0.18) × 10-8 that were subsequently employed in modeling the piperidine photo-oxidation experiments, from which the initial branchings between H-abstraction from the NH and CH2 groups, kN-H/ktot = 0.38 ± 0.08 and kC2-H/ktot = 0.49 ± 0.19, were derived. All photo-oxidation experiments were accompanied by particle formation that was initiated by the acid-base reaction of piperidine with nitric acid. Primary photo-oxidation products including both 1-nitrosopiperidine and 1-nitropiperidine were detected in the particles formed. Quantum chemistry calculations on the OH initiated atmospheric photo-oxidation of piperidine suggest the branching in the initial H-abstraction routes to be ∼35% N1, ∼50% C2, ∼13% C3, and ∼2% C4. The theoretical study produced an atmospheric photo-oxidation mechanism, according to which H-abstraction from the C2 position predominantly leads to 2,3,4,5-tetrahydropyridine and H-abstraction from the C3 position results in ring opening followed by a complex autoxidation, of which the first few steps are mapped in detail. H-abstraction from the C4 position is shown to result mainly in the formation of piperidin-4-one and 2,3,4,5-tetrahydropyridin-4-ol, whereas H-abstraction from N1 under atmospheric conditions primarily leads to 2,3,4,5-tetrahydropyridine and in minor amounts of 1-nitrosopiperidine and 1-nitropiperidine. The calculated rate coefficient for the piperidine + OH reaction agrees with the experimental value within 35%, and aligning the theoretical numbers to the experimental value results in k(T) = 2.46 × 10-12 × exp(486 K/T) cm3 molecule-1 s-1 (200-400 K).

Generation of 3-aminopropanamide and its cluster formation with nucleation precursors- a theoretical exploration

Aminoamides are formed in the atmospheric environments by the auto-oxidation of the parent diamines. In this work, the oxidation chemistry of diamine (1,3-Diaminopropane, Dap) to the amino amide (3- aminopropanamide, 3-APA) and its new particle formation potential with small atmospheric molecules such as NH3 (A), H2O (W) and H2SO4 (SA) are theoretically investigated using the M062X/6-311++G** theory. The bimolecular rate coefficient of the ·OH initiated H-atom abstraction is computed to be 1.01 × 10-11 cm3 molecule-1 s-1. Further reaction of the peroxy radical intermediate indicates that the pathway involving γ H- shift of the initially formed radical intermediates to be more favourable on kinetic grounds with the effective bimolecular rate coefficient of 3.87 × 10-14 cm3 molecule-1s-1. The thermodynamic barrier associated with the H-shifts involved in this pathway is in the range of 13-20 kcal/mol. The cluster formation of APA with SA is more favourable than the clusters with W and A, wherein the free energy of formation of (APA)(SA) and (APA)(SA)2 are -11.3 and -22.6 kcal/mol, respectively. However, the feasibility of cluster formation with W and A increases with the altitude and becomes spontaneous in the case of water at an altitude of 12 km. The present work indicates that aminoamides like 3-APA can participate in the initial stages of new particle formation events by forming clusters with SA molecules. The scattering parameters and topological analysis of different (Amide)(SA) clusters indicate more scattering properties for the (APA)(SA) cluster, which has an adverse effect on the atmosphere. Furthermore, topological analysis indicates that H-bond formation is more prominent in the (APA)(SA) cluster.

Reactions of Ethynyloxy Radical with Hydroperoxyl Radical: Bridging Theoretical Reaction Dynamics and Chemical Modeling of Combustion

A detailed and accurate combustion reaction mechanism is crucial for understanding the nature of fuel combustion. In this work, a theoretical study of reaction HCCO+HO2 using M06-2X/6-311++G(d,p) for geometry optimization and combined methods based on spin-unrestricted CCSD(T)/CBS level of theory with basis set extrapolation from MP2/aug-cc-pVnZ (n=T and Q) for energy calculations were performed. The temperature- and pressure-dependent rate coefficients at 300-2000 K and 0.01-100 atm, suitable for combustion conditions, were derived using the Rice-Ramsberger-Kassel-Marcus/Master-Equation approach. Furthermore, temperature-dependent thermochemistry data of key species for the HCCO+HO2 system has also been studied. Finally, an updated ketene model is developed by supplementing the most recent theoretical work and the theoretical work in this paper. This updated model was tested to simulate the speciation of ketene oxidation in available experimental research. It is shown that the updated model for predicting ketene oxidation exhibits a high level of agreement with experimental data across a wide range of species profiles. An analysis was conducted to identify the crucial reactions that influence ketene ignition. This paper's research findings are essential for enhancing the combustion mechanism of ketene and other hydrocarbons and oxygenated hydrocarbon fuels.

Effect of temperature on the gas-phase reaction of CH3CN with OH radicals: experimental (T = 11.7-177.5 K) and computational (T = 10-400 K) kinetic study

Acetonitrile (CH3CN) is present in the interstellar medium (ISM) in a variety of environments. However, at the ultracold temperatures of the ISM, radical-molecule reactions are not widely investigated because of the experimental handicap of getting organic molecules in the gas phase by conventional techniques. The CRESU (French acronym for Reaction Kinetics in a Uniform Supersonic Flow) technique solves this problem. For this reason, we present in this work the kinetic study of the gas-phase reaction of CH3CN with one of the most ubiquitous radicals, the hydroxyl (OH) radical, as a function of temperature (11.7-177.5 K). The kinetic technique employed to investigate the CH3CN + OH reaction was the pulsed laser photolysis-laser induced fluorescence. The rate coefficient for this reaction k(T) has been observed to drastically increase from 177.5 K to 107.0 K (about 2 orders of magnitude), while the increase in k(T) from 107.0 K to 11.7 K was milder (around 4 times). The temperature dependent expressions for k(T) are provided in the two distinct T-ranges, excluding the upper limit obtained for k(177.5 K): In addition, the rate coefficients estimated by the canonical competitive unified statistical (CCUS) theory show a similar behaviour to the experimental results, when evaluated within the high-pressure limit. This is consistent with the experimentally observed independence of k(T) with total gas density at selected temperatures. Astrochemical networks, such as the KIDA database or UMIST, do not include the CH3CN + OH reaction as a potential depletion process for acetonitrile in the ISM because the current studies predict very low rate coefficients at IS temperatures. According to the model (T = 10 K), the impact of the titled reaction on the abundances of CH3CN appears to be negligible in dark molecular clouds of the ISM (∼1% of the total depletion reactions included in UMIST network). With respect to the potential formation of the CH2CN radical in those environments, even in the most favourable scenario, where this radical could be formed in a 100% yield from the CH3CN + OH reaction, this route would only contribute around 2% to the current assumed formation routes by the UMIST network.

Gas-Phase Phenyl Radical + O2 Reacts via a Submerged Transition State

In the gas-phase chemistry of the atmosphere and automotive fuel combustion, peroxyl radical intermediates are formed following O2 addition to carbon-centered radicals which then initiate a complex network of radical reactions that govern the oxidative processing of hydrocarbons. The rapid association of the phenyl radical-a fundamental radical related to benzene-with O2 has hitherto been modeled as a barrierless process, a common assumption for peroxyl radical formation. Here, we provide an alternate explanation for the kinetics of this reaction by deploying double-hybrid density functional theory (DFT), at the DSD-PBEP86-D3(BJ)/aug-cc-pVTZ level of theory, and locate a submerged adiabatic transition state connected to a prereaction complex along the reaction entrance pathway. Using this potential energy scheme, experimental rate coefficients k(T) for the addition of O2 to the phenyl radical are accurately reproduced within a microcanonical kinetic model. This work highlights that purportedly barrierless radical oxidation reactions may instead be modeled using stationary points, which in turn provides insight into pressure and temperature dependence.

HIO3-HIO2-Driven Three-Component Nucleation: Screening Model and Cluster Formation Mechanism

Iodine oxoacids (HIO3 and HIO2)-driven nucleation has been suggested to efficiently contribute to new particle formation (NPF) in marine atmospheres. Abundant atmospheric nucleation precursors may further enhance HIO3-HIO2-driven nucleation through various multicomponent nucleation mechanisms. However, the specific enhancing potential (EP) of different precursors remains largely unknown. Herein, the EP-based screening model of precursors and enhancing mechanism of the precursor with the highest EP on HIO3-HIO2 nucleation were investigated. The formation free energies (ΔG), as critical parameters for evaluating EP, were calculated for the dimers of 63 selected precursors with HIO2. Based on the ΔG values, (1) a quantitative structure-activity relationship model was developed for evaluating ΔG of other precursors and (2) atmospheric concentrations of 63 (precursor)1(HIO2)1 dimer clusters were assessed to identify the precursors with the highest EP for HIO3-HIO2-driven nucleation by combining with earlier results for the nucleation with HIO3 as the partner. Methanesulfonic acid (MSA) was found to be one of the precursors with the highest EP. Finally, we found that MSA can effectively enhance HIO3-HIO2 nucleation at atmospheric conditions by studying larger MSA-HIO3-HIO2 clusters. These results augment our current understanding of HIO3-HIO2 and MSA-driven nucleation and may suggest a larger impact of HIO2 in atmospheric aerosol nucleation.

OH Roaming during the Ozonolysis of α-Pinene: A New Route to Highly Oxygenated Molecules?

The formation of low-volatility organic compounds in the ozonolysis of α-pinene, the dominant atmospheric monoterpene, provides an important route to aerosol formation. In this work, we consider a previously unexplored set of pathways for the formation of highly oxygenated molecules in α-pinene ozonolysis. Pioneering, direct experimental observations of Lester and co-workers have demonstrated a significant production of hydroxycarbonyl products in the dissociation of Criegee intermediates. Theoretical analyses indicate that this production arises from OH roaming-induced pathways during the OO fission of the vinylhydroperoxides (VHPs), which in turn come from internal H transfers in the Criegee intermediates. Ab initio kinetics computations are used here to explore the OH roaming-induced channels that arise from the ozonolysis of α-pinene. For computational reasons, the calculations consider a surrogate for α-pinene, where two spectator methyl groups are replaced with H atoms. Multireference electronic structure calculations are used to illustrate a variety of energetically accessible OH roaming pathways for the four VHPs arising from the ozonolysis of this α-pinene surrogate. Ab initio transition-state theory-based master equation calculations indicate that for the dissociation of stabilized VHPs, these OH roaming pathways are kinetically significant with a branching that generally increases from ∼20% at room temperature up to ∼70% at lower temperatures representative of the troposphere. For one of the VHPs, this branching already exceeds 60% at room temperature. For the overall ozonolysis process, these branching ratios would be greatly reduced by a limited branching to the stabilized VHP, although there would also be some modest roaming fraction for the nonthermal VHP dissociation process. The strong exothermicities of the roaming-induced isomerizations/additions and abstractions suggest new routes to fission of the cyclobutane rings. Such ring fissions would facilitate further autoxidation reactions, thereby providing a new route for producing highly oxygenated nonvolatile precursors to aerosol formation.

An overlooked oxidation mechanism of toluene: computational predictions and experimental validations

Secondary organic aerosols (SOAs) influence the Earth's climate and threaten human health. Aromatic hydrocarbons (AHs) are major precursors for SOA formation in the urban atmosphere. However, the revealed oxidation mechanism dramatically underestimates the contribution of AHs to SOA formation, strongly suggesting the importance of seeking additional oxidation pathways for SOA formation. Using toluene, the most abundant AHs, as a model system and the combination of quantum chemical method and field observations based on advanced mass spectrometry, we herein demonstrate that the second-generation oxidation of AHs can form novel epoxides (TEPOX) with high yield. Such TEPOX can further react with H2SO4 or HNO3 in the aerosol phase to form less-volatile compounds including novel non-aromatic and ring-retaining organosulfates or organonitrates through reactive uptakes, providing new candidates of AH-derived organosulfates or organonitrates for future ambient observation. With the newly revealed mechanism, the chemistry-aerosol box modeling revealed that the SOA yield of toluene oxidation can reach up to 0.35, much higher than 0.088 based on the original mechanism under the conditions of pH = 2 and 0.1 ppbv NO. This study opens a route for the formation of reactive uptake SOA precursors from AHs and significantly fills the current knowledge gap for SOA formation in the urban atmosphere.

Theoretical Kinetics Predictions for Reactions on the NH2O Potential Energy Surface

Recent modeling studies of ammonia oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies among the existing literature mechanisms. In this study, high-level theoretical kinetics predictions have been obtained for reactions on the NH2O potential energy surface, including the NH2 + O, HNO + H, and NH + OH reactions. These reactions have previously been highlighted as important reactions in NH3 oxidation with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations, including large basis sets and high-level corrections to yield high-accuracy (∼0.2 kcal/mol 2σ uncertainty) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature- and pressure-dependent kinetics. For radical-radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constants as well as the branching ratios. We find that the NO + H2 channel contributes 10% of the total NH2 + O flux at combustion temperatures, although this channel is not included in modern NH3 oxidation mechanisms. Modeling is used to illustrate the ramifications of these rate predictions on the kinetics of NH3 oxidation and NOx formation. The present results for NH2 + O are important for predicting the chain branching and formation of NO in the oxidation of NH3 and thermal DeNOx.

Mechanistic and Kinetic Insights into OH-Initiated Atmospheric Oxidation of Hymexazol: A Computational Study

Hymexazol is a volatile fungicide widely used in agriculture, causing its abundance in the atmosphere; thus, its atmospheric fate and conversion are of great importance when assessing its environmental impacts. Herein, we report a theoretical kinetic mechanism for the oxidation of hymexazol by OH radicals, as well as the subsequent reactions of its main products with O2 and then with NO by using the Rice-Ramsperger-Kassel-Marcus-based Master equation kinetic model on the potential energy surface explored at the ROCBS-QB3//M06-2X/aug-cc-pVTZ level. The predicted total rate constants ktotal(T, P) for the reaction between hymexazol and OH radicals show excellent agreement with scarcely available experimental values (e.g., 3.6 × 10-12 vs (4.4 ± 0.8) × 10-12 cm3/molecule/s at T = 300 K and P = 760 Torr); thus, the calculated kinetic parameters can be confidently used for modeling/simulation of N-heterocycle-related applications under atmospheric and even combustion conditions. The model shows that 3,4-dihydroxy-5-methyl-4,5-dihydro-1,2-oxazol-5-yl (IM2), 3,5-dihydroxy-5-methyl-4,5-dihydro-1,2-oxazol-4-yl (IM3), and (3-hydroxy-1,2-oxazol-5-yl)methyl (P8) are the main primary intermediates, which form the main secondary species of (3,4-dihydroxy-5-methyl-4,5-dihydro-1,2-oxazol-5-yl)dioxidanyl (IM4), (3,5-dihydroxy-5-methyl-4,5-dihydro-1,2-oxazol-4-yl)dioxidanyl (IM7), and ([(3-hydroxy-1,2-oxazol-5-yl)methyl]dioxidanyl (IM11), respectively, through the reactions with O2. The main secondary species then can react with NO to form the main tertiary species, namely, (3,4-dihydroxy-5-methyl-4,5-dihydro-1,2-oxazol-5-yl)oxidanyl (P19), (3,5-dihydroxy-5-methyl-4,5-dihydro-1,2-oxazol-4-yl)oxidanyl (P21), and [(3-hydroxy-1,2-oxazol-5-yl)methyl]oxidanyl (P23), respectively, together with NO2. Besides, hymexazol could be a persistent organic pollutant in the troposphere due to its calculated half-life τ1/2 of 13.7-68.1 h, depending on the altitude.

Kinetics of the Simplest Criegee Intermediate Reaction with Water Vapor: Revisit and Isotope Effect

The kinetics of the simplest Criegee intermediate (CH2OO) reaction with water vapor was revisited. By improving the signal-to-noise ratio and the precision of water concentration, we found that the kinetics of CH2OO involves not only two water molecules but also one and three water molecules. Our experimental results suggest that the decay of CH2OO can be described as d[CH2OO]/dt = -kobs[CH2OO]; kobs = k0 + k1[water] + k2[water]2 + k3[water]3; k1 = (4.22 ± 0.48) × 10-16 cm3 s-1, k2 = (10.66 ± 0.83) × 10-33 cm6 s-1, k3 = (1.48 ± 0.17) × 10-50 cm9 s-1 at 298 K and 300 Torr with the respective Arrhenius activation energies of Ea1 = 1.8 ± 1.1 kcal mol-1, Ea2 = -11.1 ± 2.1 kcal mol-1, Ea3 = -17.4 ± 3.9 kcal mol-1. The contribution of the k3[water]3 term becomes less significant at higher temperatures around 345 K, but it is not ignorable at 298 K and lower temperatures. By quantifying the concentrations of H2O and D2O with a Coriolis-type direct mass flow sensor, the kinetic isotope effect (KIE) was investigated at 298 K and 300 Torr and KIE(k1) = k1(H2O)/k1(D2O) = 1.30 ± 0.32; similarly, KIE(k2) = 2.25 ± 0.44 and KIE(k3) = 0.99 ± 0.13. These mild KIE values are consistent with theoretical calculations based on the variational transition state theory, confirming that the title reaction has a broad and low barrier, and the reaction coordinate involves not only the motion of a hydrogen atom but also that of an oxygen atom. Comparing the results recorded under 300 Torr (N2 buffer gas) with those under 600 Torr, a weak pressure effect of k3 was found. From quantum chemistry calculations, we found that the CH2OO + 3H2O reaction is dominated by the reaction pathways involving a ring structure consisting of two water molecules, which facilitate the hydrogen atom transfer, while the third water molecule is hydrogen-bonded outside the ring. Furthermore, analysis based on dipole capture rates showed that the CH2OO(H2O) + (H2O)2 and CH2OO(H2O)2 + H2O pathways will dominate in the three water reaction.

Gas-phase formation of glycolonitrile in the interstellar medium

Our automated reaction discovery program, AutoMeKin, has been utilized to investigate the formation of glycolonitrile (HOCH2CN) in the gas phase under the low temperatures of the interstellar medium (ISM). The feasibility of a proposed pathway depends on the absence of barriers above the energy of reactants and the availability of the suggested precursors in the ISM. Based on these criteria, several radical-radical reactions and a radical-molecule reaction have been identified as viable formation routes in the ISM. Among the radical-radical reactions, OH + CH2CN appears to be the most relevant, considering the energy of the radicals and its ability to produce glycolonitrile in a single step. However, our analysis reveals that this reaction produces hydrogen isocyanide (HNC) and formaldehyde (CH2O), with rate coefficients ranging from (7.3-11.5) × 10-10 cm3 molecule-1 s-1 across the temperature range of 10-150 K. Furthermore, the identification of this remarkably efficient pathway for HNC elimination from glycolonitrile significantly broadens the possibilities for any radical-radical mechanism proposed in our research to be considered as a feasible pathway for the formation of HNC in the ISM. This finding is particularly interesing given the persistently unexplained overabundance of hydrogen isocyanide in the ISM. Among the radical-molecule reactions investigated, the most promising one is OH + CH2CHNH, which forms glycolonitrile and atomic hydrogen with rate coefficients in the range (0.3-6.6) × 10-10 cm3 molecule-1 s-1 within the 10-150 K temperature range. Our calculations indicate that the formation of both hydrogen isocyanide and glycolonitrile is efficient under the harsh conditions of the ISM.

Enhancement of Atmospheric Nucleation Precursors on Iodic Acid-Induced Nucleation: Predictive Model and Mechanism

Iodic acid (IA) has recently been recognized as a key driver for new particle formation (NPF) in marine atmospheres. However, the knowledge of which atmospheric vapors can enhance IA-induced NPF remains limited. The unique halogen bond (XB)-forming capacity of IA makes it difficult to evaluate the enhancing potential (EP) of target compounds on IA-induced NPF based on widely studied sulfuric acid systems. Herein, we employed a three-step procedure to evaluate the EP of potential atmospheric nucleation precursors on IA-induced NPF. First, we evaluated the EP of 63 precursors by simulating the formation free energies (ΔG) of the IA-containing dimer clusters. Among all dimer clusters, 44 contained XBs, demonstrating that XBs are frequently formed. Based on the calculated ΔG values, a quantitative structure-activity relationship model was developed for evaluating the EP of other precursors. Second, amines and O/S-atom-containing acids were found to have high EP, with diethylamine (DEA) yielding the highest potential to enhance IA-induced nucleation by combining both the calculated ΔG and atmospheric concentration of considered 63 precursors. Finally, by studying larger (IA)_1-3(DEA)_1-3 clusters, we found that the IA-DEA system with merely 0.1 ppt (2.5×10⁶ cm^-3) DEA yields comparable nucleation rates to that of the IA-iodous acid system.

An Experimental and Master Equation Investigation of Kinetics of the CH2OO + RCN Reactions (R = H, CH3, C2H5) and Their Atmospheric Relevance

We have performed direct kinetic measurements of the CH₂OO + RCN reactions (R = H, CH₃, C₂H₅) in the temperature range 233-360 K and pressure range 10-250 Torr using time-resolved UV-absorption spectroscopy. We have utilized a new photolytic precursor, chloroiodomethane (CH₂ICl), whose photolysis at 193 nm in the presence of O₂ produces CH₂OO. Observed bimolecular rate coefficients for CH₂OO + HCN, CH₂OO + CH₃CN, and CH₂OO + C₂H₅CN reactions at 296 K are (2.22 ± 0.65) × 10^-14 cm³ molecule^-1 s^-1, (1.02 ± 0.10) × 10^-14 cm³ molecule^-1 s^-1, and (2.55 ± 0.13) × 10^-14 cm³ molecule^-1 s^-1, respectively, suggesting that reaction with CH₂OO is a potential atmospheric degradation pathway for nitriles. All the reactions have negligible temperature and pressure dependence in the studied regions. Quantum chemical calculations (ωB97X-D/aug-cc-pVTZ optimization with CCSD(T)-F12a/VDZ-F12 electronic energy correction) of the CH₂OO + RCN reactions indicate that the barrierless lowest-energy reaction path leads to a ring closure, resulting in the formation of a 1,2,4-dioxazole compound. Master equation modeling results suggest that following the ring closure, chemical activation in the case of CH₂OO + HCN and CH₂OO + CH₃CN reactions leads to a rapid decomposition of 1,2,4-dioxazole into a CH₂O + RNCO pair, or by a rearrangement, into a formyl amide (RC(O)NHC(O)H), followed by decomposition into CO and an imidic acid (RC(NH)OH). The 1,2,4-dioxazole, the CH₂O + RNCO pair, and the CO + RC(NH)OH pair are atmospherically significant end products to varying degrees.

Fragmentation of interstellar methanol by collisions with He˙+: an experimental and computational study

Methanol is a key species in astrochemistry as its presence and reactivity provides a primary route to the synthesis of more complex interstellar organic molecules (iCOMs) that may eventually be incorporated in newly formed planetary systems. In the interstellar medium, methanol is formed by hydrogenation of CO ices on grains, and its fate upon collisions with interstellar ions should be accounted for to correctly model iCOM abundances in objects at various stages of stellar evolution. The absolute cross sections (CSs) and branching ratios (BRs) for the collisions of He˙⁺ ions with CH₃OH are measured, as a function of the collision energy, using a Guided Ion Beam Mass Spectrometer (GIB-MS). Insights into the dissociative electron (charge) exchange mechanism have been obtained by computing the entrance and exit multidimensional Potential Energy Surfaces (PESs) and by modelling the non-adiabatic transitions using an improved Landau-Zener-Stückelberg approach. Notably, the dynamical treatment reproducing the experimental findings includes a strong orientation effect of the system formed by the small He˙⁺ ion and the highly polar CH₃OH molecule, in the electric field gradient associated to the strongly anisotropic intermolecular interaction. This is a stereodynamical effect that plays a fundamental role in collision events occurring under a variety of conditions, with kinetic energy confined within intervals ranging from the sub-thermal to the hyper-thermal regime.

Reaction of OH radicals with CH3NH2 in the gas phase: experimental (11.7-177.5 K) and computed rate coefficients (10-1000 K)

Nitrogen-bearing molecules, like methylamine (CH₃NH₂), can be the building blocks of amino acids in the interstellar medium (ISM). At the ultralow temperatures of the ISM, it is important to know its gas-phase reactivity towards interstellar radicals and the products formed. In this work, the kinetics of the OH + CH₃NH₂ reaction was experimentally and theoretically investigated at low- and high-pressure limits (LPL and HPL) between 10 and 1000 K. Moreover, the CH₂NH₂ and CH₃NH yields were computed in the same temperature range for both pressure regimes. A pulsed CRESU (French acronym for Reaction Kinetics in a Uniform Supersonic Flow) apparatus was employed to determine the rate coefficient, k(T), in the 11.7-177.5 K range. A drastic increase of k(T) when the temperature is lowered was observed in agreement with theoretical calculations, evaluated by the competitive canonical unified statistical (CCUS) theory, below 300 K in the LPL regime. The same trend was observed in the HPL regime below 350 K, but the theoretical k(T) values were higher than the experimental ones. Above 200 K, the calculated rate coefficients are improved with respect to previous computational studies and are in excellent agreement with the experimental literature data. In the LPL, the formation of CH₃NH becomes largely dominant below ca. 100 K. Conversely, in the HPL regime, CH₂NH₂ is the only product below 100 K, whereas CH₃NH becomes dominant at 298 K with a branching ratio similar to the one found in the LPL regime (≈70%). At T > 300 K, both reaction channels are competitive independently of the pressure regime.

Reaction Kinetics of Criegee Intermediates with Nitric Acid

This work investigated the reaction kinetics of HNO₃ with four Criegee intermediates (CIs): CH₂OO, (CH₃)₂COO, methyl vinyl ketone oxide (MVKO), and methacrolein oxide (MACRO). Our results show that these reactions are extremely fast with rate coefficients of (1.51 ± 0.45) × 10^-10, (3.54 ± 1.06) × 10^-10, (3.93 ± 1.18) × 10^-10, and (3.0 ± 1.0) × 10^-10 cm³ s^-1 for reactions of HNO₃ with CH₂OO, (CH₃)₂COO, syn-MVKO, and anti-MACRO, respectively. This is consistent with previous results for the reactions between CIs and carboxylic acids, but the rate coefficient of CH₂OO + HNO₃ in the literature [Foreman Angew. Chem. 2016, 128, 10575] was found to be overestimated by a factor of 3.6. In addition, we did not observe any significant pressure dependence in the HNO₃ reactions with CH₂OO and (CH₃)₂COO under 100-400 Torr. Our results indicate that in a dry area with severe NO_x pollution, the reactions of CIs with HNO₃ and their products may be worthy of attention, but these reactions may be insignificant under high-humidity conditions. However, CI reactions with HNO₃ may not play an important role in the atmospheric removal processes of HNO₃ because of the low concentrations of CIs.

Mechanism and kinetics of the reaction of the Criegee intermediate CH2OO with acetic acid studied using a step-scan Fourier-transform IR spectrometer

Acetic acid, CH₃C(O)OH, plays an important role in the acidity of the troposphere. The reactions of Criegee intermediates with CH₃C(O)OH have been proposed to be a potential source of secondary organic aerosol in the atmosphere. We investigated the detailed mechanism and kinetics of the reaction of the Criegee intermediate CH₂OO with CH₃C(O)OH. The time-resolved infrared absorption spectra of transient species produced upon irradiation at 308 nm of a flowing mixture of CH₂I₂/O₂/CH₃C(O)OH at 298 K were recorded using a step-scan Fourier-transform infrared spectrometer. The decrease in the intensity of the bands of CH₂OO was accompanied by the appearance of bands near 886, 971, 1021, 1078, 1160, 1225, 1377, 1402, 1434, and 1777 cm^-1, assigned to the absorption of hydroperoxymethyl acetate [CH₃C(O)OCH₂OOH, HPMA], the hydrogen-transferred adduct of CH₂OO and CH₃C(O)OH. Two types of conformers of HPMA, an open form and an intramolecularly hydrogen-bonded form, were identified. At a later reaction period, bands of the open-form HPMA became diminished, and new bands appeared at 930, 1045, 1200, 1378, 1792, and 1810 cm^-1, assigned to formic acetic anhydride [CH₃C(O)OC(O)H, FAA], a dehydrated product of HPMA. The intramolecularly hydrogen-bonded HPMA is more stable. From the temporal profiles of HPMA and FAA, we derived a rate coefficient k = (1.3 ± 0.3) × 10^-10 cm³ molecule^-1 s^-1 for the reaction CH₂OO + CH₃C(O)OH to form HPMA and a rate coefficient k = 980 ± 40 s^-1 for the dehydration of the open-form HPMA to form FAA. Theoretical calculations were performed to elucidate the CH₂OO + CH₃C(O)OH reaction pathway and to understand the distinct reactivity of these two forms of HPMA.

Gas-Phase Peroxyl Radical Recombination Reactions: A Computational Study of Formation and Decomposition of Tetroxides

The recombination (“dimerization”) of peroxyl radicals (RO₂•) is one of the pathways suggested in the literature for the formation of peroxides (ROOR′, often referred to as dimers or accretion products in the literature) in the atmosphere. It is generally accepted that these dimers play a major role in the first steps of the formation of submicron aerosol particles. However, the precise reaction pathways and energetics of RO₂• + R′O₂• reactions are still unknown. In this work, we have studied the formation of tetroxide intermediates (RO₄R′): their formation from two peroxyl radicals and their decomposition to triplet molecular oxygen (³O₂) and a triplet pair of alkoxyl radicals (RO•). We demonstrate this mechanism for several atmospherically relevant primary and secondary peroxyl radicals. The potential energy surface corresponds to an overall singlet state. The subsequent reaction channels of the alkoxyl radicals include, but are not limited to, their dimerization into ROOR′. Our work considers the multiconfigurational character of the tetroxides and the intermediate phases of the reaction, leading to reliable mechanistic insights for the formation and decomposition of the tetroxides. Despite substantial uncertainties in the computed energetics, our results demonstrate that the barrier heights along the reaction path are invariably small for these systems. This suggests that the reaction mechanism, previously validated at a multireference level only for methyl peroxyl radicals, is a plausible pathway for the formation of aerosol-relevant larger peroxides in the atmosphere.

Capture Theory Models: An overview of their development, experimental verification, and applications to ion-molecule reactions

Since Arrhenius first proposed an equation to account for the behaviour of thermally activated reactions in 1889, significant progress has been made in our understanding of chemical reactivity. A number of capture theory models have been developed over the past several decades to predict the rate coefficients for reactions between ions and molecules-ranging from the Langevin equation (for reactions between ions and non-polar molecules) to more recent fully quantum theories (for reactions at ultra-cold temperatures). A number of different capture theory methods are discussed, with the key assumptions underpinning each approach clearly set out. The strengths and limitations of these capture theory methods are examined through detailed comparisons between low-temperature experimental measurements and capture theory predictions. Guidance is provided on the selection of an appropriate capture theory method for a given class of ion-molecule reaction and set of experimental conditions-identifying when a capture-based model is likely to provide an accurate prediction. Finally, the impact of capture theories on fields such as astrochemical modelling is noted, with some potential future directions of capture-based approaches outlined.

Theoretical study of the NO3 radical reaction with CH2ClBr, CH2ICl, CH2BrI, CHCl2Br, and CHClBr2

The potential reaction of the nitrate radical (NO₃), the main nighttime atmospheric oxidant, with five alkyl halides, halons (CH₂ClBr, CH₂ICl, CH₂BrI, CHCl₂Br, and CHClBr₂) has been studied theoretically. The most favorable reaction corresponds to a hydrogen atom transfer. The stationary points on the potential energy surfaces of these reactions have been characterized. The reactions can be classified into two groups based on the number of hydrogen atoms in the halon molecules (1 or 2). The reactions with halons with only one hydrogen atom show more exothermic profiles than those with two hydrogen atoms. In addition, the kinetics of the reaction of NO₃ + CH₂BrI was studied in much higher detail using a multi-well Master Equation solver as a representative example of the nitrate radical reactivity against these halocarbons. These results indicate that the chemical lifetime of the alkyl halides would not be substantially affected by nitrate radical reactions, even in the case of NO₃-polluted atmospheric conditions.

Atmospheric Chemistry of N-Methylmethanimine (CH3N═CH2): A Theoretical and Experimental Study

The OH-initiated photo-oxidation of N-methylmethanimine, CH₃N=CH₂, was investigated in the 200 m³ EUPHORE atmospheric simulation chamber and in a 240 L stainless steel photochemical reactor employing time-resolved online FTIR and high-resolution PTR-ToF-MS instrumentation and in theoretical calculations based on quantum chemistry results and master equation modeling of the pivotal reaction steps. The quantum chemistry calculations forecast the OH reaction to primarily proceed via H-abstraction from the =CH₂ group and π-system C-addition, whereas H-abstraction from the −CH₃ group is a minor route and forecast that N-addition can be disregarded under atmospheric conditions. Theoretical studies of CH₃N=CH₂ photolysis and the CH₃N=CH₂ + O₃ reaction show that these removal processes are too slow to be important in the troposphere. A detailed mechanism for OH-initiated atmospheric degradation of CH₃N=CH₂ was obtained as part of the theoretical study. The photo-oxidation experiments, obstructed in part by the CH₃N=CH₂ monomer–trimer equilibrium, surface reactions, and particle formation, find CH₂=NCHO and CH₃N=CHOH/CH₂=NCH₂OH as the major primary products in a ratio 18:82 ± 3 (3σ-limit). Alignment of the theoretical results to the experimental product distribution results in a rate coefficient, showing a minor pressure dependency under tropospheric conditions and that can be parametrized k(T) = 5.70 × 10^–14 × (T/298 K)^3.18 × exp(1245 K/T) cm³ molecule^–1 s^–1 with k₂₉₈ = 3.7 × 10^–12 cm³ molecule^–1 s^–1. The atmospheric fate of CH₃N=CH₂ is discussed, and it is concluded that, on a global scale, hydrolysis in the atmospheric aqueous phase to give CH₃NH₂ + CH₂O will constitute a dominant loss process. N₂O will not be formed in the atmospheric gas phase degradation, and there are no indications of nitrosamines and nitramines formed as primary products.

Computational Study on the Reaction of β-Hydroxyethylperoxy Radical with HO2 and Effects of Water Vapor

The reaction of β-hydroxyethylperoxy radical (β-HEP) and HO₂ with and without water was studied using quantum chemistry and kinetic calculations. The main products are HOCH₂CH₂OOH and ³O₂ for the reaction with and without water, while all other reaction channels can be neglected. The rate coefficients of the reaction follow negative temperature dependence. The pseudo-second-order rate coefficients are 2-4 orders of magnitude smaller for the reaction with saturated water vapor, indicating the negligible contribution of water in this reaction. This is probably also true for other peroxy radicals (except for HO₂), indicating that a large part of previous results on the water enhancement of reaction rate coefficients might have overestimated the influence of water.

Reaction of Methylidyne with Ethane: The C-C Insertion Is Unimportant

High-accuracy coupled-cluster calculations in combination with the E,J-resolved master-equation analysis are used to study the reaction mechanism and kinetics of methylidyne with ethane. This reaction plays an important role in the combustion of hydrocarbon fuels and in interstellar chemistry. Two distinct mechanisms, the C-C and the C-H insertions of CH in C₂H₆, are characterized. The C-C insertion pathway is identified to have a large barrier of 34.5 kcal mol^-1 and hence plays no significant role in kinetics. The C-H insertion pathway is found to have no barrier, leading to a highly vibrationally excited n-C₃H₇ radical, which rapidly dissociates (within 50 ps) to yield CH₃ + C₂H₄ and H + C₃H₆ in a roughly 7:3 ratio. These findings are in good agreement with an experimental result that indicates that about 20% of the reaction goes to H + C₃H₆. The reaction of the electronically excited quartet state of the CH radical with C₂H₆ is examined for the first time and found to proceed as a direct H-abstraction via a small barrier of 0.4 kcal mol^-1 to yield triplet CH₂ and C₂H₅. The reaction on the quartet state surface is negligibly slow at low temperatures characteristic of interstellar environments but becomes important at high combustion temperatures.

Theoretical spectroscopic study of acetyl (CH 3CO), vinoxy (CH 2CHO), and 1-methylvinoxy (CH 3COCH 2) radicals. Barrierless formation processes of acetone in the gas phase

Background: Acetone is present in the earth´s atmosphere and extra-terrestrially. The knowledge of its chemical history in these environments represents a challenge with important implications for global tropospheric chemistry and astrochemistry. The results of a search for efficient barrierless pathways producing acetone from radicals in the gas phase are described in this paper. The spectroscopic properties of radicals needed for their experimental detection are provided.   Methods: The reactants were acetone fragments of low stability and small species containing C, O and H atoms. Two exergonic bimolecular addition reactions involving the radicals CH 3, CH 3CO, and CH 3COCH 2, were found to be competitive according to the kinetic rates calculated at different temperatures. An extensive spectroscopic study of the radicals CH 3COCH 2 and CH 3CO, as well as the CH 2CHO isomer, was performed. Rovibrational parameters, anharmonic vibrational transitions, and excitations to the low-lying excited states are provided. For this purpose, RCCSD(T)-F12 and MRCI/CASSCF calculations were performed. In addition, since all the species presented non-rigid properties, a variational procedure of reduced dimensionality was employed to explore the far infrared region. Results: The internal rotation barriers were determined to be V 3=143.7 cm -1 (CH 3CO), V 2=3838.7 cm -1 (CH 2CHO) and V 3=161.4 cm -1 and V 2=2727.5 cm -1 (CH 3COCH 2).The splitting of the ground vibrational state due to the torsional barrier have been computed to be 2.997 cm -1, 0.0 cm -1, and 0.320 cm -1, for CH 3CO, CH 2CHO, and CH 3COCH 2, respectively. Conclusions: Two addition reactions, H+CH 3COCH 2 and CH 3+CH 3CO, could be considered barrierless formation processes of acetone after considering all the possible formation routes, starting from 58 selected reactants, which are fragments of the molecule. The spectroscopic study of the radicals involved in the formation processes present non-rigidity. The interconversion of their equilibrium geometries has important spectroscopic effects on CH 3CO and CH 3COCH 2, but is negligible for CH 2CHO.

Proton leap: shuttling of protons onto benzonitrile

The detailed description of chemical transformations in the interstellar medium allows deciphering the origin of a number of small and medium - sized organic molecules. We present density functional theory analysis of proton transfer from the trihydrogen cation and the ethenium cation to benzonitrile, a recently discovered species in the Taurus Molecular Cloud 1. Detailed energy transformations along the reaction paths were analyzed using the interacting quantum atoms methodology, which elucidated how the proton carrier influences the lightness to deliver the proton to benzonitrile's nitrogen atom. The proton carriers' deformation energy represents the largest destabilizing effect, whereas a proton's promotion energy, the benzonitrile-proton Coulomb attraction, as well as non-classical benzonitrile-proton and carrier-proton interaction are the dominant stabilizing energy components. As two ion-molecule reactions proceed without energy barriers, rate constants were estimated using the classical capture theory and were found to be an order of magnitude larger for the reaction with the trihydrogen cation compared to that with the ethenium cation (∼10^-8 and 10^-9 cm³ s^-1, respectively). These results were obtained both with quantum chemical and ab initio molecular dynamics simulations (the latter performed at 10 K and 100 K), confirming that up to 100 K both systems choose energetically undemanding routes by tracking the corresponding minimum energy paths. A concept of a turning point is introduced, which is an equivalent to the transition state in barrierless reactions.

Master Equation Modelling of Non-equilibrium Chemistry in Stellar Outflows

A current challenge in astrochemistry is to explain the formation of Fe-Mg silicate dust around evolved stars. The dust is observed to form within 2 to 3 stellar radii of...

Improved Computational Modeling of the Kinetics of the Acetylperoxy + HO2 Reaction

The acetylperoxy + HO2 reaction has multiple impacts on the troposphere, with a triplet pathway leading to peracetic acid + O2 (reaction 1a) competing with singlet pathways leading to acetic...

Mai T, Assali M, Nguyen TTD, Nguyen H, Szőri M, Huynh LK, Fittschen C, Farooq A. Reaction Kinetics of 1,4-Cyclohexadiene with OH radicals : An Experimental and Theoretical Study

This work presents OH-initiated oxidation kinetics of 1,4-cyclochexadiene (1,4-CHD). Temperature dependence of the reaction was investigated by utilizing laser flash photolysis flow reactor and laser-induced fluorescence (LPFR/LIF) technique over the...

Dramatic Conformer-Dependent Reactivity of the Acetaldehyde Oxide Criegee Intermediate with Dimethylamine Via a 1,2-Insertion Mechanism

The reactivity of carbonyl oxides has previously been shown to exhibit strong conformer and substituent dependencies. Through a combination of synchrotron-multiplexed photoionization mass spectrometry experiments (298 K and 4 Torr) and high-level theory [CCSD(T)-F12/cc-pVTZ-F12//B2PLYP-D3/cc-pVTZ with an added CCSDT(Q) correction], we explore the conformer dependence of the reaction of acetaldehyde oxide (CH₃CHOO) with dimethylamine (DMA). The experimental data support the theoretically predicted 1,2-insertion mechanism and the formation of an amine-functionalized hydroperoxide reaction product. Tunable-vacuum ultraviolet photoionization probing of anti- or anti- + syn-CH₃CHOO reveals a strong conformer dependence of the title reaction. The rate coefficient of DMA with anti-CH₃CHOO is predicted to exceed that for the reaction with syn-CH₃CHOO by a factor of ∼34,000, which is attributed to submerged barrier (syn) versus barrierless (anti) mechanisms for energetically downhill reactions.

Predictive Combustion Kinetics of OH Radical Reactions with a C5 Unsaturated Alcohol: The Competitive H-Abstraction and OH-Addition Reactions of 2-Methyl-3-buten-2-ol

2-Methyl-3-buten-2-ol (MBO232) is a potential biofuel and renewable fuel additive. In a combustion environment, the consumption of MBO232 is mainly through the reaction with a OH radical, one of the most important oxidants. Here, we predict the intricate reactions of MBO232 and OH radicals under a broad range of combustion conditions, that is, 230-2500 K and 0.01-1000 atm. The potential energy surfaces of H-abstraction and OH-addition have been investigated at the CCSD(T)/CBS//M06-2X/def2-TZVP level, and the rate constants were calculated via Rice-Ramsperger-Kassel-Marcus/master equation (RRKM/ME) theory. The decomposition reactions of the critical intermediates from the OH-addition reactions have also been studied. Our results show that OH-addition reactions are dominant below 850 K, while H-abstraction reactions, especially the channel-abstracting H atoms from the methyl groups, are more competitive at higher temperatures. We found that it is necessary to discriminate H atoms attached to the same C atom, as their abstraction rates can differ by up to 1 order of magnitude. The calculated results show good agreement with the reported experimental data. We have provided the modified Arrhenius expressions for rate constants of the dominant channels. The kinetic data determined in this work are of much value for constructing the combustion models of MBO232 and understanding the combustion kinetics and mechanism of other unsaturated alcohols.

Atmospheric Chemistry of 2-Amino-2-methyl-1-propanol: A Theoretical and Experimental Study of the OH-Initiated Degradation under Simulated Atmospheric Conditions

The OH-initiated degradation of 2-amino-2-methyl-1-propanol [CH₃C(NH₂)(CH₃)CH₂OH, AMP] was investigated in a large atmospheric simulation chamber, employing time-resolved online high-resolution proton-transfer reaction-time-of-flight mass spectrometry (PTR-ToF-MS) and chemical analysis of aerosol online PTR-ToF-MS (CHARON-PTR-ToF-MS) instrumentation, and by theoretical calculations based on M06-2X/aug-cc-pVTZ quantum chemistry results and master equation modeling of the pivotal reaction steps. The quantum chemistry calculations reproduce the experimental rate coefficient of the AMP + OH reaction, aligning k(T) = 5.2 × 10^–12 × exp (505/T) cm³ molecule^–1 s^–1 to the experimental value k_exp,300K = 2.8 × 10^–11 cm³ molecule^–1 s^–1. The theoretical calculations predict that the AMP + OH reaction proceeds via hydrogen abstraction from the −CH₃ groups (5–10%), −CH₂– group, (>70%) and −NH₂ group (5–20%), whereas hydrogen abstraction from the −OH group can be disregarded under atmospheric conditions. A detailed mechanism for atmospheric AMP degradation was obtained as part of the theoretical study. The photo-oxidation experiments show 2-amino-2-methylpropanal [CH₃C(NH₂)(CH₃)CHO] as the major gas-phase product and propan-2-imine [(CH₃)₂C=NH], 2-iminopropanol [(CH₃)(CH₂OH)C=NH], acetamide [CH₃C(O)NH₂], formaldehyde (CH₂O), and nitramine 2-methyl-2-(nitroamino)-1-propanol [AMPNO₂, CH₃C(CH₃)(NHNO₂)CH₂OH] as minor primary products; there is no experimental evidence of nitrosamine formation. The branching in the initial H abstraction by OH radicals was derived in analyses of the temporal gas-phase product profiles to be B_CH3/B_CH2/B_NH2 = 6:70:24. Secondary photo-oxidation products and products resulting from particle and surface processing of the primary gas-phase products were also observed and quantified. All the photo-oxidation experiments were accompanied by extensive particle formation that was initiated by the reaction of AMP with nitric acid and that mainly consisted of this salt. Minor amounts of the gas-phase photo-oxidation products, including AMPNO₂, were detected in the particles by CHARON-PTR-ToF-MS and GC×GC-NCD. Volatility measurements of laboratory-generated AMP nitrate nanoparticles gave Δ_vapH = 80 ± 16 kJ mol^–1 and an estimated vapor pressure of (1.3 ± 0.3) × 10^–5 Pa at 298 K. The atmospheric chemistry of AMP is evaluated and a validated chemistry model for implementation in dispersion models is presented.

Substituent Effect in the Reactions between Criegee Intermediates and 3-Aminopropanol

Via intramolecular H atom transfer, 3-aminopropanol is more reactive toward Criegee intermediates, in comparison with amines or alcohols. Here we accessed the substituent effect of Criegee intermediates in their reactions with 3-aminopropanol. Through real-time monitoring the concentrations of two Criegee intermediates with their strong UV absorption at 340 nm, the experimental rate coefficients at 298 K (100-300 Torr) were determined to be (1.52 ± 0.08) × 10^-11 and (1.44 ± 0.22) × 10^-13 cm³ s^-1 for the reactions of 3-aminopropanol with (CH₃)₂COO (acetone oxide) and CH₂CHC(CH₃)OO (methyl vinyl ketone oxide), respectively. Compared to our previous experimental value for the reaction with syn-CH₃CHOO, (1.24 ± 0.13) × 10^-11 cm³ s^-1, we can see that the methyl substitution at the anti position has little effect on the reactivity while the vinyl substitution causes a drastic decrease in the reactivity. Our theoretical calculations based on CCSD(T)-F12 energies reproduce this 2-order-of-magnitude decrease in the rate coefficient caused by the vinyl substitution. Using the activation strain model, we found that the interaction of Criegee intermediates with 3-aminopropanol is weaker for the case of vinyl substitution. This effect can be further rationalized by the delocalization of the lowest unoccupied molecular orbital for the vinyl-substituted Criegee intermediates. These results would help us better estimate the impact of similar reactions like the reactions of Criegee intermediates with water vapor, some of which could be difficult to measure experimentally but can be important in the atmosphere. We also found that the B2PLYP-D3BJ/aug-cc-pVTZ calculation can reproduce the CCSD(T)-F12 reaction barrier energies within ca. 1 kcal mol^-1, indicating that the use of the B2PLYP-D3BJ method is promising for future predictions of the reactions of larger Criegee intermediates.

Theoretical and Experimental Study of 3-Pentanol Autoignition: Ab Initio Calculation, Shock Tube Experiments, and Kinetic Modeling

3-Pentanol is a potential alternative fuel or a green fuel additive for modern engines. The H-abstraction reactions from 3-pentanol by H, CH₃, HO₂, and OH radicals are significant in the 3-pentanol oxidation process. However, corresponding rate constants are forced to rely on either analogy from sec-butanol or estimation from alkanes due to a lack of direct experimental and theoretical study. In this work, stationary points on the potential energy surfaces (PESs) were calculated with the high-level DLPNO-CCSD(T)/CBS(T-Q)//M06-2X/cc-pVTZ method, which is further used to benchmark against the CBS-QB3 method. Then, the high-pressure limit rate constants for target reactions, over a broad range of temperature (400-2000 K), were calculated with the phase-space theory and conventional transition state theory. A comparison was made between the calculated rate constants and the values available in Carbonnier et al. [ Proc. Combust. Inst. 2019, 37(1), 477-484]. The rate constants for the above H-abstraction reactions in the Carbonnier model were updated with the calculated results, followed by a modification based on the computed results of 3-pentanol + HO₂ to obtain the revised model. Validation against the shock tube (ST) and the jet-stirred reactor (JSR) measurements from the literature proved the revised model an optimal one. Furthermore, using an ST, ignition delay times (IDTs) for the 3-pentanol/air mixtures were measured spanning a temperature range of 920-1450 K, pressures of 6, 10, and 20 bar, and equivalence ratios of 0.5, 1.0, and 1.5. Generally, IDTs decrease with increasing temperature and reflected shock pressure. Improved predictions to present experimental data were obtained by using the revised model as compared with the Carbonnier model. Finally, sensitivity analysis was conducted using the revised model to gain an in-depth comprehension of the 3-pentanol autoignition.

New Mechanistic Insights into Atmospheric Oxidation of Aniline Initiated by OH Radicals

This study theoretically reports the comprehensive kinetic mechanism of the aniline + OH reaction in the range of 200-2000 K and 0.76-7600 Torr. The temperature- and pressure-dependent behaviors, including time-resolved species profiles and rate coefficients, were studied within the stochastic RRKM-based master equation framework with the reaction energy profile, together with molecular properties of the species involved, characterized at the M06-2X/aug-cc-pVTZ level. Hindered internal rotation and Eckart tunneling treatments were included. The H-abstraction from the -NH₂ moiety (to form C₆H₅NH (P1)) is found to prevail over the OH-addition on the C atom at the ortho site of aniline (to form 6-hydroxy-1-methylcyclohexa-2,4-dien-1-yl (I2)) with the atmospheric rate expressions (in cm³/molecule/s) as k^abstraction(P1) = 3.41 × 10¹ × T^-4.56 × exp (-255.2 K/T) for 200-2000 K and k^addition(I2) = 3.68 × 10⁹ × T^-7.39 × exp (-1163.9 K/T) for 200-800 K. The U-shaped temperature-dependent characteristics and weakly positive pressure dependence at low temperatures (e.g., T ≤ 800 K and P = 760 Torr) of k_total(T) are also observed. The disagreement in k_total(T) between the previous calculations and experimental studies is also resolved, and atmospheric aniline is found to be primarily removed by OH radicals (τ_OH ∼ 1.1 h) in the daytime. Also, via TD-DFT simulations, it is recommended to include P1 and I2 in any atmospheric photolysis-related model.

Quasi-Classical Trajectory Study of the CN + NH3 Reaction Based on a Global Potential Energy Surface

Based on a combination of valence-bond and molecular mechanics functions which were fitted to high-level ab initio calculations, we constructed an analytical full-dimensional potential energy surface, named PES-2020, for the hydrogen abstraction title reaction for the first time. This surface is symmetrical with respect to the permutation of the three hydrogens in ammonia, it presents numerical gradients and it improves the description presented by previous theoretical studies. In order to analyze its quality and accuracy, stringent tests were performed, exhaustive kinetics and dynamics studies were carried out using quasi-classical trajectory calculations, and the results were compared with the available experimental evidence. Firstly, the properties (geometry, vibrational frequency and energy) of all stationary points were found to reasonably reproduce the ab initio information used as input; due to the complicated topology with deep wells in the entrance and exit channels and a “submerged” transition state, the description of the intermediate complexes was poorer, although it was adequate to reasonably simulate the kinetics and dynamics of the title reaction. Secondly, in the kinetics study, the rate constants simulated the experimental data in the wide temperature range of 25–700 K, improving the description presented by previous theoretical studies. In addition, while previous studies failed in the description of the kinetic isotope effects, our results reproduced the experimental information. Finally, in the dynamics study, we analyzed the role of the vibrational and rotational excitation of the CN(v,j) reactant and product angular scattering distribution. We found that vibrational excitation by one quantum slightly increased reactivity, thus reproducing the only experimental measurement, while rotational excitation strongly decreased reactivity. The scattering distribution presented a forward-backward shape, associated with the presence of deep wells along the reaction path. These last two findings await experimental confirmation.

Comment on “Atmospheric chemistry of oxazole: the mechanism and kinetic studies on oxidation reaction initiated by OH radicals” by A. Shiroudi, M. A. Abdel-Rahman, A. M. El-Nahas and M. Altarawneh, New J. Chem., 2021, , 2237

The large kinetic discrepancy between computational and experimental studies is resolved using the rigorous stochastic RRKM-based master-equation rate model. Detailed mechanistic insights are also revealed to advance its related applications.

Ab initio kinetic mechanism of OH-initiated atmospheric oxidation of pyrrole

The comprehensive kinetic mechanism of the OH-initiated gas-phase oxidation of pyrrole is first theoretically reported in a broad range of conditions (T = 200-2000 K &P = 1-7600 Torr). On the potential energy surface constructed at the M06-2X/aug-cc-pVTZ level, the temperature- and pressure-dependent behaviors of the title reaction were characterized using the stochastic Rice-Ramsperger-Kassel-Marcus based Master Equation (RRKM-ME) rate model. The corrections of the hindered internal rotation and quantum tunneling treatments were included. The calculated results reveal the competition between the two distinct pathways: OH-addition and direct H-abstraction. The former channels are found favorable at low-temperature and high-pressure range (e.g., T < 900 K and P = 760 Torr) where non-Arrhenius and positive pressure-dependent behaviors of the rate constants are noticeably observed, while the latter predominate at temperatures higher than 900 K at atmospheric pressure and no pressure dependence on the rate constant is found. The predicted global rate constants are in excellent agreement with laboratory values; thus, the derived kinetic parameters are recommended for modeling/simulation of N-heterocycle-related applications in atmospheric and even in combustion conditions. Besides, pyrrole should not be considered as a persistent organic pollutant owing to its short atmospheric lifetime (∼1 h) towards OH radicals. The secondary mechanisms of the subsequent reactions of two OH-pyrrole adducts (namely, I1 and I2) with two abundant species, O₂/NO, which are relevant to the atmospheric degradation process, were also investigated. It is also revealed by TD-DFT calculations that two OH-pyrrole adducts (I1 &I2), nine intermediates, Ii (i = 3-11) and four products (P1, P2, P3 and P6) can undergo photodissociation under the sunlight.

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.

Detailed kinetics of hydrogen abstraction from trans-decalin by OH radicals: the role of hindered internal rotation treatment

In this study, the detailed kinetic mechanism of the trans-decalin + OH reaction is firstly investigated for a wide range of conditions (i.e., T = 200-2000 K & P = 0.76-76 000 Torr) using the M06-2X/aug-cc-pVTZ level and stochastic Rice-Ramsperger-Kassel-Marcus based master equation (RRKM-ME) rate model, which includes corrections of the hindered internal rotor (HIR) and tunneling effects. Our predicted global rate constant excellently matches with the scarce experimental measurement (R. Atkinson, et al. Int. J. Chem. Kinet., 1983, 15, 37-50). The H-abstraction channel from Cα of trans-decalin is found to be dominant at low temperatures. A U-shaped temperature-dependent behavior and slightly positive pressure-dependence at low temperatures (e.g., T ≤ 400 K & P = 760 Torr) of the total rate constants are also observed. Detailed analysis reveals that the HIR treatment is essential to capture the kinetic behavior while the tunneling correction only plays a minor role.

Atmospheric Chemistry of Methyl Isocyanide-An Experimental and Theoretical Study

The reaction of CH3NC with OH radicals was studied in smog chamber experiments employing PTR-ToF-MS and long-path FTIR detection. The rate coefficient was determined to be kCH3NC+OH = (7.9 ± 0.6) × 10-11 cm3 molecule-1 s-1 at 298 ± 3 K and 1013 ± 10 hPa; methyl isocyanate was the sole observed product of the reaction. The experimental results are supported by CCSD(T*)-F12a/aug-cc-pVTZ//M06-2X/aug-cc-pVTZ quantum chemistry calculations showing the reaction to proceed primarily via electrophilic addition to the isocyanide carbon atom. On the basis of the quantum chemical data, the kinetics of the OH reaction was simulated using a master equation model revealing the rate coefficient to be nearly independent of pressure at tropospheric conditions and having a negative temperature dependence with kOH = 4.2 × 10-11 cm3 molecule-1 s-1 at 298 K. Additional quantum chemistry calculations on the CH3NC reactions with O3 and NO3 show that these reactions are of little importance under atmospheric conditions. The atmospheric fate of methyl isocyanide is discussed.

Formation Mechanisms of Iodine-Ammonia Clusters in Polluted Coastal Areas Unveiled by Thermodynamics and Kinetic Simulations

It has been revealed that iodine species play important roles in atmospheric new particle formations (NPFs) in pristine coastal areas. However, it is unclear whether other atmospheric species, such as NH₃, for which the levels in coastal areas of China are >2.5 × 10¹⁰ molecules·cm^-3 are involved in the NPFs of iodine species, although NH₃ has been proved to promote particle formation of H₂SO₄. Via high-level quantum chemical calculations and atmospheric cluster dynamic code simulations, this study unveiled new mechanisms of nucleation, in which NH₃ mediates the formation of iodine particles by assisting hydrolysis of I₂O₅ or reacting with HIO₃. The simulated formation rates of iodine-ammonia clusters via the new mechanisms are much higher than those simulated via sequential addition of HIO₃ with subsequent release of H₂O, under the condition that NH₃ concentrations are higher than 10¹⁰ molecules·cm^-3. The new mechanisms can well explain the observed cluster formation rates at a coastal site in Zhejiang of China. The findings not only expand the current understandings of the role of NH₃ in NPFs but also highlight the importance of monitoring and evaluating NPFs via the iodine-ammonia cluster pathway in the coastal areas of China and other regions worldwide.