Long-range transition state theory

A theoretical study of the addition of CH2OO to hydroxymethyl hydroperoxide and its implications on SO3 formation in the atmosphere

The reaction of hydroxymethyl hydroperoxide (HMHP, HOCH₂OOH) with the simplest Criegee intermediate, CH₂OO, has been examined using quantum chemical methods with transition state theory. Geometry optimization and IRC calculations were performed using the M06-2X, MN15-L, and B2PLYP-D3 functionals in conjunction with the aug-cc-pVTZ basis set. Single point energy calculations using QCISD(T) and BD(T) with the same basis set have been performed to determine the energy of reactants, reactive complexes, transition states, and products. Rate coefficients have been obtained using variational transition state theory. The addition of CH₂OO on the three different oxygen atoms in HMHP has been considered and the ether oxide forming channel, CH₂OO + HOCH₂OOH → HOCH₂O(O)CH₂OOH (channel 2), is the most favorable. The best computed standard enthalpy of reaction (ΔH) and zero-point corrected barrier height are -20.02 and -6.33 kcal mol^-1, respectively. The reaction barrier is negative and our results suggest that both the inner and outer transition states contribute to the corresponding overall reactive flux in the tropospheric temperature range (220 K to 320 K). A two-transition state model has been used to obtain reliable rate coefficients at the high-pressure limit. The pressure-dependent rate coefficient calculations using the SS-QRRK theory have shown that this channel is pressure-dependent. Moreover, our investigation has shown that the ether oxide formed may rapidly react with SO₂ at 298 K to form SO₃, which can, in turn, react with water to form atmospheric H₂SO₄. A similar calculation has been conducted for the reaction of HMHP with OH, suggesting that the titled reaction may be a significant sink of HMHP. Therefore, the reaction between CH₂OO and HOCH₂OOH could be an indirect source for generating atmospheric H₂SO₄, which is crucial to the formation of clouds, and it might relieve global warming.

Theoretical Insight into the Reaction Mechanism and Kinetics for the Criegee Intermediate of anti-PhCHOO with SO2

In this study, the density functional theory (DFT) and CCSD(T) method have been performed to gain insight into the possible products and detailed reaction mechanism of the Criegee intermediate (CI) of anti-PhCHOO with SO2 for the first time. The potential energy surfaces (PESs) have been depicted at the UCCSD(T)/6-311++G(d,p)//UB3LYP/6-311++G(d,p) levels of theory with ZPE correction. Two different five-membered ring adducts, viz., endo PhCHOOS(O)O (IM1) and exo PhCHOOS(O)O (IM2) have been found in the entrance of reaction channels. Both direct and indirect reaction pathways from IM1 and IM2 have been considered for the title reaction. Our calculations show that the formation of PhCHO+SO3 (P1) via indirect reaction pathways from IM1 is predominant in all the pathways, and the production of P1 via direct dissociation pathway of IM1 and indirect reaction pathways of IM2 cannot be neglected. Moreover, PhCOOH+SO2 (P2) initiated from IM2 is identified as the minor product. According to the kinetic calculation, the total rate constant for the anti-PhCHOO+SO2 reaction is estimated to be 6.98 × 10−10 cm3·molecule−1·s−1 at 298 K.

Atmospheric oxidation mechanism and kinetics of isoprene initiated by chlorine radicals: A computational study

The reaction with chlorine radicals (·Cl) has been considered to be one of indispensable sinks for isoprene. However, the mechanism of ·Cl initiated isoprene reaction was not fully understood. Herein, the reaction of isoprene with ·Cl, and ensuing reactions of the resulting isoprene relevant radicals were investigated by combined quantum chemistry calculations and kinetics modeling. The results indicate that ·Cl addition to two terminal C-atoms of two double bonds of isoprene, forming IM_1-1 and IM_1-4, are more favorable than H-abstractions from isoprene. Interestingly, the predicted reaction rate constant for the direct H-abstraction pathway is much lower than that of the indirect one, clarifying a direct H-abstraction mechanism for previously experimental observation. The IM_1-1 and IM_1-4 have distinct fate in their subsequent transformation. The reaction of IM_1-1 ends after the one-time O₂ addition. However, IM_1-4 can follow auto-oxidation mechanism with two times O₂ addition to finally form highly oxidized multi-functional molecules (HOMs), C₅H₇ClO₃ and ·OH. More importantly, the estimated contribution of ·Cl on HOMs (monomer only) formation from isoprene is lower than that of ·OH in addition pathway, implying overall HOMs yield from atmospheric isoprene oxidation could be overestimated if the role of ·Cl in transforming isoprene is ignored.

Association of Cl with C2H2 by unified variable-reaction-coordinate and reaction-path variational transition-state theory

Barrierless unimolecular association reactions are prominent in atmospheric and combustion mechanisms but are challenging for both experiment and kinetics theory. A key datum for understanding the pressure dependence of association and dissociation reactions is the high-pressure limit, but this is often available experimentally only by extrapolation. Here we calculate the high-pressure limit for the addition of a chlorine atom to acetylene molecule (Cl + C₂H₂→C₂H₂Cl). This reaction has outer and inner transition states in series; the outer transition state is barrierless, and it is necessary to use different theoretical frameworks to treat the two kinds of transition state. Here we study the reaction in the high-pressure limit using multifaceted variable-reaction-coordinate variational transition-state theory (VRC-VTST) at the outer transition state and reaction-path variational transition state theory (RP-VTST) at the inner turning point; then we combine the results with the canonical unified statistical (CUS) theory. The calculations are based on a density functional validated against the W3X-L method, which is based on coupled cluster theory with single, double, and triple excitations and a quasiperturbative treatment of connected quadruple excitations [CCSDT(Q)], and the computed rate constants are in good agreement with some of the experimental results. The chlorovinyl (C₂H₂Cl) adduct has two isomers that are equilibrium structures of a double-well C≡C-H bending potential. Two procedures are used to calculate the vibrational partition function of chlorovinyl; one treats the two isomers separately and the other solves the anharmonic energy levels of the double well. We use these results to calculate the standard-state free energy and equilibrium constant of the reaction.

Methanesulfonic Acid-driven New Particle Formation Enhanced by Monoethanolamine: A Computational Study

Amines are recognized as significant enhancing species on methanesulfonic acid (MSA)-driven new particle formation (NPF). Monoethanolamine (MEA) has been detected in the atmosphere, and its concentration could be significantly increased once MEA-based postcombustion CO₂ capture technology is widely implemented. Here, we evaluated the enhancing potential of MEA on MSA-driven NPF by examining the formation of MEA-MSA clusters using a combination of quantum chemical calculations and kinetics modeling. The results indicate that the -OH group of MEA can form at least one hydrogen bond with MSA or MEA in all MEA-containing clusters. The enhancing potential of MEA is higher than that of the strongest enhancing agent known so far, methylamine (MA), for MSA-driven NPF. Such high enhancing potential can be ascribed to not only the higher gas-phase basicity but also the role of the additional -OH group of MEA in increasing the binding free energy by forming additional hydrogen bonds. This clarifies the importance of hydrogen-bonding capacity from the nonamino group of amines in enhancing MSA-driven NPF. The main growth pathway for MEA-MSA clusters proceeds via the initial formation of the (MEA)₁(MSA)₁ cluster, followed by alternately adding one MSA and one MEA molecule, differing from the case of MA-MSA clusters.

A Barrierless Pathway Accessing the C9H9 and C9H8 Potential Energy Surfaces via the Elementary Reaction of Benzene with 1-Propynyl

The crossed molecular beams reactions of the 1-propynyl radical (CH₃CC; X²A₁) with benzene (C₆H₆; X¹A_1g) and D6-benzene (C₆D₆; X¹A_1g) were conducted to explore the formation of C₉H₈ isomers under single-collision conditions. The underlying reaction mechanisms were unravelled through the combination of the experimental data with electronic structure and statistical RRKM calculations. These data suggest the formation of 1-phenyl-1-propyne (C₆H₅CCCH₃) via the barrierless addition of 1-propynyl to benzene forming a low-lying doublet C₉H₉ intermediate that dissociates by hydrogen atom emission via a tight transition state. In accordance with our experiments, RRKM calculations predict that the thermodynamically most stable isomer – the polycyclic aromatic hydrocarbon (PAH) indene – is not formed via this reaction. With all barriers lying below the energy of the reactants, this reaction is viable in the cold interstellar medium where several methyl-substituted molecules have been detected. Its underlying mechanism therefore advances our understanding of how methyl-substituted hydrocarbons can be formed under extreme conditions such as those found in the molecular cloud TMC-1. Implications for the chemistry of the 1-propynyl radical in astrophysical environments are also discussed.

Experimental and theoretical investigation on the OH + CH3C(O)CH3 reaction at interstellar temperatures (T=11.7-64.4 K)

The rate coefficient, k(T), for the gas-phase reaction between OH radicals and acetone CH3C(O)CH3, has been measured using the pulsed CRESU (French acronym for Reaction Kinetics in a Uniform Supersonic Flow) technique (T = 11.7-64.4 K). The temperature dependence of k(T = 10-300 K) has also been computed using a RRKM-Master equation analysis after partial revision of the potential energy surface. In agreement with previous studies we found that the reaction proceeds via initial formation of two pre-reactive complexes both leading to H2O + CH3C(O)CH2 by H-abstraction tunneling. The experimental k(T) was found to increase as temperature was lowered. The measured values have been found to be several orders of magnitude higher than k(300 K). This trend is reproduced by calculations, with a special good agreement with experiments below 25 K. The effect of total gas density on k(T) has been explored. Experimentally, no pressure dependence of k(20 K) and k(64 K) was observed, while k(50 K) at the largest gas density 4.47×1017 cm-3 is twice higher than the average values found at lower densities. The computed k(T) is also reported for 103 cm-3 of He (representative of the interstellar medium). The predicted rate coefficients at 10 K surround the experimental value which appears to be very close to the low pressure regime prevailing in the interstellar medium. For gas-phase model chemistry of interstellar molecular clouds, we suggest using the calculated value of 1.8×10-10 cm3 molecule-1 s-1 at 10 K and the reaction products are water and CH3C(O)CH2 radicals.

Impact of water on the BrO + HO2 gas-phase reaction: mechanism, kinetics and products

The BrO + HO2 reaction, which participates in the cycle of ozone removal via BrOH formation, was explored both in the absence and in the presence of water using ab initio calculations. Two main sets of products, (i) HBr + O3 and (ii) BrOH + O2, are formed regardless of the presence of water, following a hydrogen abstraction mechanism. The HBr + O3 products are formed from the intermediate BrOOOH adduct, whereas BrOH + O2 are formed either from the intermediate OBrOOH adduct or via a barrierless hydrogen transfer from HO2 to BrO. Owing to the formation of molecular oxygen that can bear different spin configurations, the formation of BrOH + O2 products was examined both on the singlet and the triplet surfaces. Under relevant atmospheric temperatures and pressure, the formation of products (i) is energetically and kinetically less favorable than that of products (ii). The rate coefficient at 298 K for the HBr + O3 formation was determined to be 2.00 × 10-20 cm3 molecule-1 s-1, and found to decrease by 1-2 orders of magnitude when one or both reactants are clustered with water. For the formation of BrOH + O2, a rate coefficient of 2.21 × 10-11 cm3 molecule-1 s-1 is determined on both singlet and triplet surfaces in the absence of water. Though this rate coefficient slightly decreases for the hydrated reactions, the fractions of the reactants that are effectively complexed with water are not high enough to shift the overall BrOH + O2 formation rate. The current study further indicates that humidity plays a negligible role in ozone removal via the BrO + HO2 reaction.

Ozonolysis of 3-carene in the atmosphere. Formation mechanism of hydroxyl radical and secondary ozonides

The gas-phase ozonolysis mechanism of 3-carene is investigated using high level quantum chemistry and kinetic calculations. The reaction follows the Criegee mechanism with an initial addition of O₃ to the [double bond splayed left]C[double bond, length as m-dash]C[double bond splayed right] bond, followed by a chain of unimolecular isomerizations, as 3-carene + O₃→ POZs (primary ozonides) → CIs (Criegee intermediates, 4 conformers) → Ps (products). In the course of the reaction, a large excess of energy retained in the POZs* lead to the prompt unimolecular processes in POZs*, CIs*, and Ps*, and only ∼4% of CIs* could be stabilized by collision at 298 K and 760 Torr. From RRKM-ME calculations, the VHPs* could further dissociate to vinoxy-type radical and OH radical, the SOZs* could isomerize to 3-caronic acid, and DIOs* could be stabilized via collision. The fractional yield of OH radical, in the range of 0.56 to 0.59, agrees reasonably well with the previously measured value of 1.06 (with an uncertainty factor of 1.5).

Temperature and Pressure Dependent Rate Coefficients for the Reaction of Ketene with Hydroxyl Radical

The reaction of ketene with hydroxyl radical is drawing growing attention, for it is found to constitute an important step during the combustion of hydrocarbon and oxygenated hydrocarbon fuels, e.g., acetylene, propyne, allene, acetone, gasoline, diesel, jet fuels, and biofuels. We studied the potential energy surface (PES) of this reaction using B2PLYP-D3/cc-PVTZ for geometry optimization and composite methods based on CCSD(T)-F12/cc-PVTZ-F12 for energy calculations. From this PES, temperature- and pressure-dependent rate coefficients and branching ratios at 200-3000 K and 0.01-100 atm were derived using the RRKM/ME approach. The reaction is dominated by four product channels: (i) OH addition on the olefinic carbon of ketene to form CH₂OH + CO, which is the most dominant under all conditions; (ii) H abstraction producing HCCO + H₂O, which is favored at high temperatures; (iii) OH addition on the carbonyl carbon to form CH₃ + CO₂, which is favored at low pressures and high temperatures; and (iv) collisional stabilization of CH₂COOH, which is favored at high pressures and low temperatures. With increasing temperatures, the overall rate constant k_overall exhibit first negative but then positive temperature dependency, with its switching point (also the minimum point) at ∼400 K. Both product channel CH₂OH + CO and HCCO + H₂O are independent of pressure, whereas formation of CH₃ + CO₂ and collisional stabilization of CH₂COOH are highly pressure dependent. Fitted modified Arrhenius expressions of the calculated rate constants are provided for the purpose of combustion modeling.

Comment on “Atmospheric oxidation reactions of imidazole initiated by hydroxyl radicals: kinetics and mechanism of reactions and atmospheric implications” by Safaei et al., Phys. Chem. Chem. Phys., 2019, 21, 8445

This communication resolves the disagreement in the chemical kinetics of the imidazole + OH reaction using a more rigorous master equation/Rice–Ramsperger–Kassel–Marcus rate model on a more comprehensive potential energy surface.

Benchmarking dual-level MS-Tor and DLPNO-CCSD(T) methods for H-abstraction from methyl pentanoate by an OH radical

Methyl pentanoate (MP) was recently proposed as a potential biodiesel surrogate due to its negative temperature coefficient region.

Ab initio dynamics of hydrogen abstraction from N2H4 by OH radicals: an RRKM-based master equation study

The detailed kinetic mechanism of the N₂H₄ + OH reaction is comprehensively reported for a wide condition range of conditions (i.e., 200–3000 K & 1–7600 Torr) using the CCSD(T)/CBS//M06-2X/6-311++G(3df,2p) level and the RRKM-based master equation rate model.

Theoretical Study on the Dynamics and Kinetics of the Reaction of CH2OH with OH

The stochastic one-dimensional chemical master equation (CME) simulation method was used to investigate the dynamics of the reaction of CH₂OH with OH. A multiwell multichannel potential energy surface (PES) was constructed at the CCSD(T)/Aug-cc-pVTZ//CBS-QB3 and QCISD(T)/Aug-cc-pVTZ//CBS-QB3 levels of theory. The constructed PES consisted of three chemically activated intermediates and two van der Waals complexes. The fractional population analysis unraveled the role of the energized intermediates and van der Waals complexes in the early stages of this complex reaction. The CME calculations provided the phenomenological rate constants through analysis of the eigenvalues and eigenvectors of collision matrices while Leonard-Jones potential was used to model the collisions. The CME results indicated that CH₂O and H₂O were the major products, in accordance with the literature. Also, the findings declared the temperature and pressure independence of the reaction over a wide range of temperature (250 to 2400 K) and pressure (0.1 to 7 atm). Furthermore, the efficiency of tunneling on the hydrogen transfer isomerization reaction of trans-HCOH to CH₂O was confirmed over the temperature range of 250 to 3000 K. The rate constants for different reaction channels are reported.

Atmospheric Oxidation Mechanism of Sabinene Initiated by the Hydroxyl Radicals

The atmospheric oxidation mechanism of sabinene initiated by the OH radical has been studied using quantum chemistry calculations at the CBS-QB3 level and reaction kinetic calculations using transition state theory and unimolecular rate theory coupled with collisional energy transfer. The oxidation is initiated by OH radical additions to the CH₂═C< bond with a branching ratio of ∼(92-96)%, while all the hydrogen atom abstractions count for ∼(4-8)% of branching ratio, which was estimated by comparing the rate coefficients of the reactions of sabinene and sabinaketon with the OH radical. Addition of OH to the ═C< carbon forms radical adduct Ra, while addition of OH to the terminal CH₂═ carbon forms radical adduct Rb, which would break the three-membered ring promptly and almost completely to radical Re. RRKM-ME calculations obtained fractional yields of 0.40, 0.09, and 0.51 for radicals syn-Ra, anti-Ra, and Re, respectively, at 298 K and 760 Torr. In the atmosphere, the syn/ anti-Ra radical would ultimately transform to sabinaketone in the presence of ppbv levels of NO, while in the transformation of the Re radical, both bimolecular reactions and unimolecular H-migrations could occur competitively for the peroxy radicals formed. The H-migrations in peroxy radicals result in the formation of unsaturated multifunctional compounds containing >C═O, -OH, and/or -OOH groups. Formation of sabinaketone from syn- and anti-Ra and formation of acetone from Re are predicted with yields of ∼0.37 and ∼0.38 in the presence of high NO, being larger than while in reasonable agreement with the experimental values of 0.19-0.23 and of 0.21-0.27, respectively.

Atmospheric Oxidation of Piperazine Initiated by ·Cl: Unexpected High Nitrosamine Yield

Chlorine radicals (·Cl) initiated amine oxidation plays an important role for the formation of carcinogenic nitrosamine in the atmosphere. Piperazine (PZ) is considered as a potential atmospheric pollutant since it is an alternative solvent to monoethanolamine (MEA), a benchmark solvent in a leading CO₂ capture technology. Here, we employed quantum chemical methods and kinetics modeling to investigate ·Cl-initiated atmospheric oxidation of PZ, particularly concerning the potential of PZ to form nitrosamine compared to MEA. Results showed that the ·Cl-initiated PZ reaction exclusively leads to N-center radicals (PZ-N) that mainly react with NO to produce nitrosamine in their further reaction with O₂/NO. Together with the PZ + ·OH reaction, the PZ-N yield from PZ oxidation is still lower than that of the corresponding MEA reactions. However, the nitrosamine yield of PZ is higher than the reported value for MEA when [NO] is <5 ppb, a concentration commonly encountered in a polluted urban atmosphere. The unexpected high nitrosamine yield from PZ compared to MEA results from a more favorable reaction of N-center radicals with NO compared to O₂. These findings show that the yield of N-center radicals cannot directly be used as a metric for the yield of the corresponding carcinogenic nitrosamine.

Rate Coefficient Measurements and Theoretical Analysis of the OH + ( E)-CF3CH═CHCF3 Reaction

Rate coefficients, k, for the gas-phase reaction of the OH radical with ( E)-CF₃CH═CHCF₃ (( E)-1,1,1,4,4,4-hexafluoro-2-butene, HFO-1336mzz(E)) were measured over a range of temperatures (211-374 K) and bath gas pressures (20-300 Torr; He, N₂) using a pulsed laser photolysis-laser-induced fluorescence (PLP-LIF) technique. k₁( T) was independent of pressure over this range of conditions with k₁(296 K) = (1.31 ± 0.15) × 10^-13 cm³ molecule^-1 s^-1 and k₁( T) = (6.94 ± 0.80) × 10^-13exp[-(496 ± 10)/ T] cm³ molecule^-1 s^-1, where the uncertainties are 2σ, and the pre-exponential term includes estimated systematic error. Rate coefficients for the OD reaction were also determined over a range of temperatures (262-374 K) at 100 Torr (He). The OD rate coefficients were ∼15% greater than the OH values and showed similar temperature dependent behavior with k₂( T) = (7.52 ± 0.44) × 10^-13exp[-(476 ± 20)/ T] and k₂(296 K) = (1.53 ± 0.15) × 10^-13 cm³ molecule^-1 s^-1. The rate coefficients for reaction 1 were also measured using a relative rate technique between 296 and 375 K with k₁(296 K) measured to be (1.22 ± 0.1) × 10^-13 cm³ molecule^-1 s^-1, in agreement with the PLP-LIF results. In addition, the 296 K rate coefficient for the O₃ + ( E)-CF₃CH═CHCF₃ reaction was determined to be <5.2 × 10^-22 cm³ molecule^-1 s^-1. A theoretical computational analysis is presented to interpret the observed positive temperature dependence for the addition reaction and the significant decrease in OH reactivity compared to the ( Z)-CF₃CH═CHCF₃ stereoisomer reaction. The estimated atmospheric lifetime of ( E)-CF₃CH═CHCF₃, due to loss by reaction with OH, is estimated to be ∼90 days, while the actual lifetime will depend on the location and season of its emission. Infrared absorption spectra of ( E)-CF₃CH═CHCF₃ were measured and used to estimate the 100 year time horizon global warming potentials (GWP) of 32 (atmospherically well-mixed) and 14 (lifetime-adjusted).

Theoretical and Experimental Study on the Reaction of tert-Butylamine with OH Radicals in the Atmosphere

The OH-initiated atmospheric degradation of tert-butylamine (tBA), (CH₃)₃CNH₂, was investigated in a detailed quantum chemistry study and in laboratory experiments at the European Photoreactor (EUPHORE) in Spain. The reaction was found to mainly proceed via hydrogen abstraction from the amino group, which in the presence of nitrogen oxides (NO _x), generates tert-butylnitramine, (CH₃)₃CNHNO₂, and acetone as the main reaction products. Acetone is formed via the reaction of tert-butylnitrosamine, (CH₃)₃CNHNO, and/or its isomer tert-butylhydroxydiazene, (CH₃)₃CN═NOH, with OH radicals, which yield nitrous oxide (N₂O) and the (CH₃)₃Ċ radical. The latter is converted to acetone and formaldehyde. Minor predicted and observed reaction products include formaldehyde, 2-methylpropene, acetamide and propan-2-imine. The reaction in the EUPHORE chamber was accompanied by strong particle formation which was induced by an acid-base reaction between photochemically formed nitric acid and the reagent amine. The tert-butylaminium nitrate salt was found to be of low volatility, with a vapor pressure of 5.1 × 10^-6 Pa at 298 K. The rate of reaction between tert-butylamine and OH radicals was measured to be 8.4 (±1.7) × 10^-12 cm³ molecule^-1 s^-1 at 305 ± 2 K and 1015 ± 1 hPa.

Directed gas phase formation of silicon dioxide and implications for the formation of interstellar silicates

Interstellar silicates play a key role in star formation and in the origin of solar systems, but their synthetic routes have remained largely elusive so far. Here we demonstrate in a combined crossed molecular beam and computational study that silicon dioxide (SiO2) along with silicon monoxide (SiO) can be synthesized via the reaction of the silylidyne radical (SiH) with molecular oxygen (O2) under single collision conditions. This mechanism may provide a low-temperature path-in addition to high-temperature routes to silicon oxides in circumstellar envelopes-possibly enabling the formation and growth of silicates in the interstellar medium necessary to offset the fast silicate destruction.

Under what conditions does (SiO)N nucleation occur? A bottom-up kinetic modelling evaluation

Silicon monoxide (SiO) is a structurally complex compound exhibiting differentiated oxide-rich and silicon-rich nano-phases at length scales covering nanoclusters to the bulk. Although nano-sized and nano-segregated SiO has great technological potential (e.g. nano-silicon for optical applications) and is of enormous astronomical interest (e.g. formation of silicate cosmic dust) an accurate general description of SiO nucleation is lacking. Avoiding the deficiencies of a bulk-averaged approach typified by classical nucleation theory (CNT) we employ a bottom-up kinetic model which fully takes into account the atomistic details involved in segregation. Specifically, we derive a new low energy benchmark set of segregated (SiO)_N cluster ground state candidates for N ≤ 20 and use the accurately calculated properties of these isomers to calculate SiO nucleation rates. We thus provide a state-of-the art evaluation of the range of pressure and temperature conditions for which formation of SiO will or will not proceed. Our results, which match with available experiment, reveal significant deficiencies with CNT approaches. We employ our model to shed light on controversial issue of circumstellar silicate dust formation showing that, at variance with the predictions from CNT-based calculations, pure SiO nucleation under such conditions is not viable.

Pressure-dependent kinetics of methyl formate reactions with OH at combustion, atmospheric and interstellar temperatures

Pressure dependence occurs in bimolecular hydrogen abstraction reactions at combustion, atmospheric and interstellar temperatures.

Gas-phase ozonolysis of furans, methylfurans, and dimethylfurans in the atmosphere

Ozonolysis of methylfurans contributes significantly to their atmospheric degradation.

Unveiling the Ionization Energy of the CN Radical

The cyano radical is a ubiquitous molecule and was, for instance, one of the first species detected in astrophysical media such as comets or diffuse clouds. In photodissociation regions, the reaction rate of CN⁺ + CO → CN + CO⁺ is one of the critical parameters defining nitrile chemistry. The enthalpy of this charge transfer reaction is defined as the difference of ionization energies (E_I) between CN and CO. Although E_I(CO) is known accurately, the E_I(CN) values are more dispersed and deduced indirectly from thermodynamic thresholds only, all above E_I(CO), leading to the assumption that the reaction is fast even at low temperature. Using a combination of synchrotron radiation, electron/ion imaging coincidence techniques, and supporting ab initio calculations, we directly determine the first adiabatic ionization energy of CN at 13.956(7) eV, and we demonstrate that E_I(CN) < E_I(CO). The findings suggest a very slow reaction in the cold regions of interstellar media.

First-Principles Chemical Kinetic Modeling of Methyl trans-3-Hexenoate Epoxidation by HO2

The design of innovative combustion processes relies on a comprehensive understanding of biodiesel oxidation kinetics. The present study aims at unraveling the reaction mechanism involved in the epoxidation of a realistic biodiesel surrogate, methyl trans-3-hexenoate, by hydroperoxy radicals using a bottom-up theoretical kinetics methodology. The obtained rate constants are in good agreement with experimental data for alkene epoxidation by HO₂. The impact of temperature and pressure on epoxidation pathways involving H-bonded and non-H-bonded conformers was assessed. The obtained rate constant was finally implemented into a state-of-the-art detailed combustion mechanism, resulting in fairly good agreement with engine experiments.

Atmospheric chemical reaction mechanism and kinetics of 1,2-bis(2,4,6-tribromophenoxy)ethane initiated by OH radical: a computational study

The mechanism and kinetics of OH-initiated oxidation of BTBPE, an alternative of PBDEs, were investigated.

Variational transition state theory: theoretical framework and recent developments

This article reviews the fundamentals of variational transition state theory (VTST), its recent theoretical development, and some modern applications.

Silicon chemistry in the mesosphere and lower thermosphere

Silicon is one of the most abundant elements in cosmic dust, and meteoric ablation injects a significant amount of Si into the atmosphere above 80 km. In this study, a new model for silicon chemistry in the mesosphere/lower thermosphere is described, based on recent laboratory kinetic studies of Si, SiO, SiO2, and Si+. Electronic structure calculations and statistical rate theory are used to show that the likely fate of SiO2 is a two-step hydration to silicic acid (Si(OH)4), which then polymerizes with metal oxides and hydroxides to form meteoric smoke particles. This chemistry is then incorporated into a whole atmosphere chemistry-climate model. The vertical profiles of Si+ and the Si+/Fe+ ratio are shown to be in good agreement with rocket-borne mass spectrometric measurements between 90 and 110 km. Si+ has consistently been observed to be the major meteoric ion around 110 km; this implies that the relative injection rate of Si from meteoric ablation, compared to metals such as Fe and Mg, is significantly larger than expected based on their relative chondritic abundances. Finally, the global abundances of SiO and Si(OH)4 show clear evidence of the seasonal meteoric input function, which is much less pronounced in the case of other meteoric species.

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

Reaction Kinetics of Meteoric Sodium Reservoirs in the Upper Atmosphere

The gas-phase reactions of a selection of sodium-containing species with atmospheric constituents, relevant to the chemistry of meteor-ablated Na in the upper atmosphere, were studied in a fast flow tube using multiphoton ionization time-of-flight mass spectrometry. For the first time, unambiguous observations of NaO and NaOH in the gas phase under atmospheric conditions have been achieved. This enabled the direct measurement of the rate constants for the reactions of NaO with H2, H2O, and CO, and of NaOH with CO2, which at 300-310 K were found to be (at 2σ confidence level): k(NaO + H2O) = (2.4 ± 0.6) × 10(-10) cm(3) molecule (-1) s(-1), k(NaO + H2) = (4.9 ± 1.2) × 10(-12) cm(3) molecule (-1) s(-1), k(NaO + CO) = (9 ± 4) × 10(-11) cm(3) molecule (-1) s(-1), and k(NaOH + CO2 + M) = (7.6 ± 1.6) × 10(-29) cm(6) molecule (-2) s(-1) (P = 1-4 Torr). The NaO + H2 reaction was found to make NaOH with a branching ratio ≥ 99%. A combination of quantum chemistry and statistical rate theory calculations are used to interpret the reaction kinetics and extrapolate the atmospherically relevant experimental results to mesospheric temperatures and pressures. The NaO + H2O and NaOH + CO2 reactions act sequentially to provide the major atmospheric sink of meteoric Na and therefore have a significant impact on the underside of the Na layer in the terrestrial mesosphere: the newly determined rate constants shift the modeled peak to about 93 km, i.e., 2 km higher than observed by ground-based lidars. This highlights further uncertainties in the Na chemistry cycle such as the unknown rate constant of the NaOH + H reaction. The fast Na-recycling reaction between NaO and CO and a re-evaluated rate constant of the NaO + CO2 sink should be now considered in chemical models of the Martian Na layer.

A theoretical study of the atmospherically important radical–radical reaction BrO + HO2; the product channel O2(a1Δg) + HOBr is formed with the highest rate

A theoretical study has been made of the BrO + HO₂ reaction, a radical–radical reaction which contributes to ozone depletion in the atmosphere via production of HOBr.

Quantum Chemical Study on ·Cl-Initiated Atmospheric Degradation of Monoethanolamine

Recent findings on the formation of ·Cl in continental urban areas necessitate the consideration of ·Cl initiated degradation when assessing the fate of volatile organic pollutants. Monoethanolamine (MEA) is considered as a potential atmospheric pollutant since it is a benchmark and widely utilized solvent in a leading CO2 capture technology. Especially, ·Cl may have specific interactions with the N atom of MEA, which could make the MEA + ·Cl reaction have different pathways and products from those of the MEA + ·OH reaction. Hence, ·Cl initiated reactions with MEA were investigated by a quantum chemical method [CCSD(T)/aug-cc-pVTZ//MP2/6-31+G(3df,2p)] and kinetics modeling. Results show that the overall rate constant for ·Cl initiated H-abstraction of MEA is 5 times faster than that initiated by ·OH, and the tropospheric lifetimes of MEA will be overestimated by 6-46% when assuming that [·Cl]/[·OH] = 1-10% if the role of ·Cl is ignored. The MEA + ·Cl reaction exclusively produces MEA-N that finally transforms into several products including mutagenic nitramine and carcinogenic nitrosamine via further reactions with O2/NOx, and the contribution of ·Cl to their formation is about 25-250% of that of ·OH. Thus, it is necessary to consider ·Cl initiated tropospheric degradation of MEA for its risk assessment.

Ring-Polymer Molecular Dynamics for the Prediction of Low-Temperature Rates: An Investigation of the C((1)D) + H2 Reaction

Quantum mechanical calculations are important tools for predicting the rates of elementary reactions, particularly for those involving hydrogen and at low temperatures where quantum effects become increasingly important. These approaches are computationally expensive, however, particularly when applied to complex polyatomic systems or processes characterized by deep potential wells. While several approximate techniques exist, many of these have issues with reliability. The ring-polymer molecular dynamics method was recently proposed as an accurate and efficient alternative. Here, we test this technique at low temperatures (300-50 K) by analyzing the behavior of the barrierless C((1)D) + H2 reaction over the two lowest singlet potential energy surfaces. To validate the theory, rate coefficients were measured using a supersonic flow reactor down to 50 K. The experimental and theoretical rates are in excellent agreement, supporting the future application of this method for determining the kinetics and dynamics of a wide range of low-temperature reactions.