Atmospheric Kinetics: Bimolecular Reactions of Carbonyl Oxide by a Triple-Level Strategy

Quantitative kinetics reveal that reactions of HO2 are a significant sink for aldehydes in the atmosphere and may initiate the formation of highly oxygenated molecules via autoxidation

Large aldehydes are widespread in the atmosphere and their oxidation leads to secondary organic aerosols. The current understanding of their chemical transformation processes is limited to hydroxyl radical (OH) oxidation during daytime and nitrate radical (NO3) oxidation during nighttime. Here, we report quantitative kinetics calculations of the reactions of hexanal (C5H11CHO), pentanal (C4H9CHO), and butanal (C3H7CHO) with hydroperoxyl radical (HO2) at atmospheric temperatures and pressures. We find that neither tunneling nor multistructural torsion anharmonicity should be neglected in computing these rate constants; strong anharmonicity at the transition states is also important. We find rate constants for the three reactions in the range 3.2-7.7 × 10-14 cm3 molecule-1 s-1 at 298 K and 1 atm, showing that the HO2 reactions can be competitive with OH and NO3 oxidation under some conditions relevant to the atmosphere. Our findings reveal that HO2-initiated oxidation of large aldehydes may be responsible for the formation of highly oxygenated molecules via autoxidation. We augment the theoretic studies with laboratory flow-tube experiments using an iodide-adduct time-of-flight chemical ionization mass spectrometer to confirm the theoretical predictions of peroxy radicals and the autoxidation pathway. We find that the adduct from HO2 + C5H11CHO undergoes a fast unimolecular 1,7-hydrogen shift with a rate constant of 0.45 s-1. We suggest that the HO2 reactions make significant contributions to the sink of aldehydes.

First-principles study of molecular hydrogen binding to heme in competition with O2, NO and CO

Molecular hydrogen shows antioxidant activity and distinct efficacy towards vascular diseases, but the understanding of this is not yet satisfactory at the atomic level. In this work, we study the binding properties of H2 to the heme group in relation with other diatomic molecules (DMs), including O2, NO and CO, and their displacement reactions, using first-principles calculations. We carry out molecular modeling of the heme group, using iron-porphyrin with the imidazole ligand, i.e., FePIm, and smaller models of Fe(CnHn+2N2)2NH3 with n = 3 and 1, and of molecular complexes of heme-DM and -H. Through analysis of optimized geometries and energetics, it is found that the order of binding strength of DMs or H to the Fe of heme is NO > O2 > CO > H > H2 for FePIm-based systems, while it is H > O2 > NO > CO > H2 for model-based systems. We calculate the activation energies for displacement reactions of H2 and H by other DMs, revealing that the H2 displacements occur spontaneously while the H displacements require a large amount of energy. Finally, our calculations corroborate that the rate constants increase with increasing temperature according to the Arrhenius relation.

Pressure and temperature dependent kinetics and the reaction mechanism of Criegee intermediates with vinyl alcohol: a theoretical study

Criegee intermediates (CIs), the key intermediates in the ozonolysis of olefins in atmosphere, have received much attention due to their high activity. The reaction mechanism of the most simple Criegee intermediate CH2OO with vinyl alcohol (VA) was investigated by using the HL//M06-2X/def2TZVP method. The temperature and pressure dependent rate constant and product branching ratio were calculated using the master equation method. For CH2OO + syn-VA, 1,4-insertion is the main reaction channel while for the CH2OO + anti-VA, cycloaddition and 1,2-insertion into the O-H bond are more favorable than the 1,4-insertion reaction. The 1,4-insertion or cycloaddition intermediates are stabilized collisionally at 300 K and 760 torr, and the dissociation products involving OH are formed at higher temperature and lower pressure. The rate constants of the CH2OO reaction with syn-VA and anti-VA both show negative temperature effects, and they are 2.95 × 10-11 and 2.07 × 10-13 cm3 molecule-1 s-1 at 300 K, respectively, and the former is agreement with the prediction in the literature.

Experimental and theoretical study of Criegee intermediate (CH2OO) reactions with n-butyraldehyde and isobutyraldehyde: kinetics, implications and atmospheric fate

The reactions of the simplest Criegee intermediate (CH2OO) with n-butyraldehyde (nBD) and isobutyraldehyde (iBD) were studied at 253-318 K and (50 ± 2) torr, using Cavity Ring-down spectroscopy (CRDS). The rate coefficients obtained at room temperature were (2.63 ± 0.14) × 10-12 and (2.20 ± 0.21) × 10-12 cm3 molecule-1 s-1 for nBD and iBD, respectively. Both the reactions show negative temperature-dependency, following equations, knBD(T = 253-318 K) = (11.51 ± 4.33) × 10-14 × exp{(918.1 ± 107.2)/T} and kiBD(T = 253-318 K) = (6.23 ± 2.29) × 10-14 × exp{(1051.4 ± 105.2)/T} cm3 molecule-1 s-1. High-pressure limit rate coefficients were determined from theoretical calculations at the CCSD(T)-F12/cc-pVTZ-F12//B3LYP/6-311+G(2df, 2p) level of theory, with <40% deviation from the experimental results at room temperature and above. The kinetic simulations were performed using a master equation solver to predict the temperature-dependency of the rate coefficients at the experimental pressure, as well as to predict the contribution of individual pathways. The major products predicted from the theoretical calculations were formaldehyde and formic acid, along with butyric acid from nBD and isobutyric acid from iBD reactions.

Rapid Atmospheric Reactions between Criegee Intermediates and Hypochlorous Acid

Hypochlorous acid (HOCl) is a paramount compound in the atmosphere due to its significant contribution to both tropospheric oxidation capacity and ozone depletion. The main removal routes for HOCl are photolysis and the reaction with OH during the daytime, while these processes are unimportant during the nighttime. Here, we report the rapid reactions of Criegee intermediates (CH2OO and anti/syn-CH3CHOO) with HOCl by using high-level quantum chemical methods as the benchmark; their accuracy is close to coupled cluster theory with single, double, and triple excitations and quasiperturbative connected quadruple excitations with a complete basis limit by extrapolation [CCSDT(Q)/CBS]. Their rate constants have been calculated by using a dual-level strategy; this combines conventional transition state theory calculated at the benchmark level with variational transition state theory with small-curvature tunneling by a validated density functional method. We find that the introduction of the methyl group into Criegee intermediates not only affects their reactivities but also exerts a remarkable influence on anharmonicity. The calculated results uncover that anharmonicity increases the rate constants of CH2OO + HOCl by a factor of 18-5, while it is of minor importance in the anti/syn-CH3CHOO + HOCl reaction at 190-350 K. The present findings reveal that the loose transition state for anti-CH3CHOO and HOCl is a rate-determining step at 190-350 K. We also find that the reaction of Criegee intermediates with HOCl contributes significantly to the sink of HOCl during the nighttime in the atmosphere.

Ammonolysis of Glyoxal at the Air-Water Nanodroplet Interface

The reactions of glyoxal (CHO)2 ) with amines in cloud processes contribute to the formation of brown carbon and oligomer particles in the atmosphere. However, their molecular mechanisms remain unknown. Herein, we investigate the ammonolysis mechanisms of glyoxal with amines at the air-water nanodroplet interface. We identified three and two distinct pathways for the ammonolysis of glyoxal with dimethylamine and methylamine by using metadynamics simulations at the air-water nanodroplet interface, respectively. Notably, the stepwise pathways mediated by the water dimer for the reactions of glyoxal with dimethylamine and methylamine display the lowest free energy barriers of 3.6 and 4.9 kcal ⋅ mol-1 , respectively. These results showed that the air-water nanodroplet ammonolysis reactions of glyoxal with dimethylamine and methylamine were more feasible and occurred at faster rates than the corresponding gas phase ammonolysis, the OH+(CHO)2 reaction, and the aqueous phase reaction of glyoxal, leading to the dominant removal of glyoxal. Our results provide new and important insight into the reactions between carbonyl compounds and amines, which are crucial in forming nitrogen-containing aerosol particles.

Quantitative kinetics of the atmospheric reaction between isocyanic acid and hydroxyl radicals: post-CCSD(T) contribution, anharmonicity, recrossing effects, torsional anharmonicity, and tunneling

Hydroxyl radicals (OH) are the most important atmospheric oxidant, initiating atmospheric reactions for the chemical transformation of volatile organic compounds. Here, we choose the HNCO + OH reaction as a prototype reaction because it contains the fundamental reaction processes for OH radicals: H-abstraction reaction by OH and OH addition reaction. However, its kinetics are unknown under atmospheric conditions. We investigate the reaction of HNCO with OH by using the GMM(P).L method close to the accuracy of single, double, triple, and quadruple excitations and noniterative quintuple excitations with a complete basis set (CCSDTQ(P)/CBS) as benchmark results and a dual-level strategy for kinetics calculations. The calculated rate constant of HNCO + OH is in good agreement with the experimental data available at the temperatures between 620 and 2500 K. We find that the rate constant cannot be correctly obtained by using experimental data to extrapolate the atmospheric temperature ranges. We find that the post-CCSD(T) contribution is very large for the barrier height with the value of -0.85 kcal mol-1 for the H-abstraction reaction, while the previous investigations were done up to the CCSD(T) level. Moreover, we also find that recrossing effects, tunneling, torsional anharmonicity, and anharmonicity are important for obtaining quantitative kinetics in the OH + HNCO reaction.

Theoretical Kinetics studies of isoprene peroxy radical chemistry: The fate of Z-δ-(4-OH, 1-OO)-ISOPOO radical

The OH radical recycling mechanism in isoprene oxidation is one of the most exciting topics in atmospheric chemistry, and the corresponding studies expand our understanding of oxidation mechanisms of volatile organic compounds in the troposphere and provide reliable evidence to improve and develop conventional atmospheric models. In this work, we performed a detailed theoretical kinetics study on the Z-δ-(4-OH, 1-OO)-ISOPOO radical chemistry, which is proposed as the heart of OH recycling in isoprene oxidation. With the full consideration of its accumulation and consumption channels, we studied and discussed the fate of Z-δ-(4-OH, 1-OO)-ISOPOO radical by solving the energy-resolved master equation over a broad range of conditions, including not only room temperatures but also high temperatures of a forest fire or low temperatures and pressures of the upper troposphere. We found non-negligible pressure dependence of its fate at combustion temperatures (up to two orders of magnitude) and demonstrated the significance of both the multi-structural torsional anharmonicity and tunneling for accurately calculating kinetics of the studied system. More interestingly, the tunneling effect on the phenomenological rate constants of the H-shift reaction channel is also found to be pressure-dependent due to the competition with the O2 loss reaction. In addition, our time evolution calculations revealed a two-stage behavior of critical species in this reaction system and estimated the shortest half-lives for the Z-δ-(4-OH, 1-OO)-ISOPOO radical at various temperatures, pressures and altitudes. This detailed kinetics study of Z-δ-(4-OH, 1-OO)-ISOPOO radical chemistry offers a typical example to deeply understand the core mechanism of OH recycling pathways in isoprene oxidation, and provides valuable insights for promoting the development of relevant atmospheric models.

Significant influence of water molecules on the SO3 + HCl reaction in the gas phase and at the air-water interface

The products resulting from the reactions between atmospheric acids and SO3 have a catalytic effect on the formation of new particles in aerosols. However, the SO3 + HCl reaction in the gas-phase and at the air-water interface has not been considered. Herein, this reaction was explored exhaustively by using high-level quantum chemical calculations and Born Oppenheimer molecular dynamics (BOMD) simulations. The quantum calculations show that the gas-phase reaction of SO3 + HCl is highly unlikely to occur under atmospheric conditions with a high energy barrier of 22.6 kcal mol-1. H2O and (H2O)2 play obvious catalytic roles in reducing the energy barrier of the SO3 + HCl reaction by over 18.2 kcal mol-1. The atmospheric lifetimes of SO3 show that the (H2O)2-assisted reaction dominates over the H2O-assisted reaction within the altitude range of 0-5 km, whereas the H2O-assisted reaction is more favorable within an altitude range of 10-50 km. BOMD simulations show that H2O-induced formation of the ClSO3-⋯H3O+ ion pair and HCl-assisted formation of the HSO4-⋯H3O+ ion pair were identified at the air-water interface. These routes followed a stepwise reaction mechanism and proceeded at a picosecond time scale. Interestingly, the formed ClSO3H in the gas phase has a tendency to aggregate with sulfuric acids, ammonias, and water molecules to form stable clusters within 40 ns simulation time, while the interfacial ClSO3- and H3O+ can attract H2SO4, NH3, and HNO3 for particle formation from the gas phase to the water surface. Thus, this work will not only help in understanding the SO3 + HCl reaction driven by water molecules in the gas-phase and at the air-water interface, but it will also provide some potential routes of aerosol formation from the reaction between SO3 and inorganic acids.

Atmospheric Chemistry of NH2SO3H in Polluted Areas: An Unexpected Isomerization of NH2SO3H in Acid-Polluted Regions

NH2SO3H is an effective nucleation agent for the formation of atmospheric aerosols and cloud particles. So, the ammonolysis of SO3 to form NH2SO3H without and with neutral (H2O) and basic (NH3) trace gases has been extensively investigated. However, the acidic trace gas X (X = H2SO4 and CH3SO3H)-assisted ammonolysis of SO3 is still up for debate. In this work, a comprehensive theoretical investigation of X-assisted ammonolysis of SO3 and its reverse reaction (the isomerization of NH2SO3H to form SO3-···NH3+) was carried out in the gas phase and at the air-water interface. The gas-phase results show that X-assisted isomerization of NH2SO3H to form SO3-···NH3+ is more energetically and kinetically favorable than its reverse reaction and the isomerization of NH2SO3H in the presence of H2O and NH3. Such unexpected findings revealed that gas-phase NH2SO3H is highly reactive in the presence of acidic trace gas in contrast to the high stability of NH2SO3H in neutral and basic conditions. At the air-water interface, the X-assisted isomerization reaction of NH2SO3H involves multiple water molecules. The loop structure of the reaction center (X···NH2SO3H···3H2O) promotes the transfer of protons in the water molecules to form the SO3-···NH3+ ion pair, which can then interact with several interfacial water molecules to form ammonium bisulfate. These interfacial reaction channels follow a stepwise mechanism and proceed at the picosecond time-scale. The findings of this study will contribute to a better understanding of the atmospheric behavior of NH2SO3H in polluted acidic trace gases.

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.

Kinetics of Sulfur Trioxide Reaction with Water Vapor to Form Atmospheric Sulfuric Acid

Although experimental methods can be used to obtain the quantitative kinetics of atmospheric reactions, experimental data are often limited to a narrow temperature range. The reaction of SO3 with water vapor is important for elucidating the formation of sulfuric acid in the atmosphere; however, the kinetics is uncertain at low temperatures. Here, we calculate rate constants for reactions of sulfur trioxide with two water molecules. We consider two mechanisms: the SO3···H2O + H2O reaction and the SO3 + (H2O)2 reaction. We find that beyond-CCSD(T) contributions to the barrier heights are very large, and multidimensional tunneling, unusually large anharmonicity of high-frequency modes, and torsional anharmonicity are important for obtaining quantitative kinetics. We find that at lower temperatures, the formation of the termolecular precursor complexes, which is often neglected, is rate-limiting compared to passage through the tight transition states. Our calculations show that the SO3···H2O + H2O mechanism is more important than the SO3 + (H2O)2 mechanism at 5-50 km altitudes. We find that the rate ratio between SO3···H2O + H2O and SO3 + (H2O)2 is greater than 20 at altitudes between 10 and 35 km, where the concentration of SO3 is very high.

Unimolecular Reactions of E-Glycolaldehyde Oxide and Its Reactions with One and Two Water Molecules

The kinetics of Criegee intermediates are important for atmospheric modeling. However, the quantitative kinetics of Criegee intermediates are still very limited, especially for those with hydroxy groups. Here, we calculate rate constants for the unimolecular reaction of E-glycolaldehyde oxide [E-hydroxyethanal oxide, E-(CH₂OH)CHOO], for its reactions with H₂O and (H₂O)₂, and for the reaction of the E-(CH₂OH)CHOO…H₂O complex with H₂O. For the highest level of electronic structure, we use W3X-L//CCSD(T)-F12a/cc-pVDZ-F12 for the unimolecular reaction and the reaction with water and W3X-L//DF-CCSD(T)-F12b/jun-cc-pVDZ for the reaction with 2 water molecules. For the dynamics, we use a dual-level strategy that combines conventional transition state theory with the highest level of electronic structure and multistructural canonical variational transition state theory with small-curvature tunneling with a validated density functional for the electronic structure. This dynamical treatment includes high-frequency anharmonicity, torsional anharmonicity, recrossing effects, and tunneling. We find that the unimolecular reaction of E-(CH₂OH)CHOO depends on both temperature and pressure. The calculated results show that E-(CH₂OH)CHOO…H₂O + H₂O is the dominant entrance channel, while previous investigations only considered Criegee intermediates + (H₂O)₂. In addition, we find that the atmospheric lifetime of E-(CH₂OH)CHOO with respect to 2 water molecules is particularly short with a value of 1.71 × 10^-6 s at 0 km, which is about 2 orders of magnitude shorter than those usually assumed for Criegee intermediate reactions with water dimer. We also find that the OH group in E-(CH₂OH)CHOO enhances its reactivity.

Direct gas-phase formation of formic acid through reaction of Criegee intermediates with formaldehyde

Ozonolysis of isoprene is considered to be an important source of formic acid (HCOOH), but its underlying reaction mechanisms related to HCOOH formation are poorly understood. Here, we report the kinetic and product studies of the reaction between the simplest Criegee intermediate (CH₂OO) and formaldehyde (HCHO), both of which are the primary products formed in ozonolysis of isoprene. By utilizing time-resolved infrared laser spectrometry with the multifunctional dual-comb spectrometers, the rate coefficient k_CH2OO+HCHO is determined to be (4.11 ± 0.25) × 10^-12 cm³ molecule^-1 s^-1 at 296 K and a negative temperature dependence of the rate coefficient is observed and described by an Arrhenius expression with an activation energy of (-1.81 ± 0.04) kcal mol^-1. Moreover, the branching ratios of the reaction products HCOOH + HCHO and CO + H₂O + HCHO are explored. The yield of HCOOH is obtained to be 37-54% over the pressure (15-60 Torr) and temperature (283-313 K) ranges. The atmospheric implications of the reaction CH₂OO + HCHO are also evaluated by incorporating these results into a global chemistry-transport model. In the upper troposphere, the percent loss of CH₂OO by HCHO is found by up to 6% which can subsequently increase HCOOH mixing ratios by up to 2% during December-January-February months.

Comprehensive Theoretical Study on Four Typical Intramolecular Hydrogen Shift Reactions of Peroxy Radicals: Multireference Character, Recommended Model Chemistry, and Kinetics

Intramolecular hydrogen shift reactions in peroxy radicals (RO₂^• → ^•QOOH) play key roles in the low-temperature combustion and in the atmospheric chemistry. In the present study, we found that a mild-to-moderate multireference character of a potential energy surface (PES) is widely present in four typical hydrogen shift reactions of peroxy radicals (RO₂^•, R = ethyl, vinyl, formyl methyl, and acetyl) by a systematic assessment based on the T₁ diagnostic, %TAE diagnostic, M diagnostic, and contribution of the dominant configuration of the reference CASSCF wavefunction (C₀²). To assess the effects of these inherent multireference characters on electronic structure calculations, we compared the PESs of the four reactions calculated by the multireference method CASPT2 in the complete basis set (CBS) limit, single-reference method CCSD(T)-F12, and single-reference-based composite method WMS. The results showed that ignoring the multireference character will introduce a mean unsigned deviation (MUD) of 0.46-1.72 kcal/mol from CASPT2/CBS results by using the CCSD(T)-F12 method or a MUD of 0.49-1.37 kcal/mol by WMS for three RO₂^• reactions (R = vinyl, formyl methyl, and acetyl) with a stronger multireference character. Further tests by single-reference Kohn-Sham (KS) density functional theory methods showed even larger deviations. Therefore, we specifically developed a new hybrid meta-generalized gradient approximation (GGA) functional M06-HS for the four typical H-shift reactions of peroxy radicals based on the WMS results for the ethyl peroxy radical reaction and on the CASPT2/CBS results for the others. The M06-HS method has an averaged MUD of 0.34 kcal/mol over five tested basis sets against the benchmark PESs, performing best in the tested 38 KS functionals. Last, in a temperature range of 200-3000 K, with the new functional, we calculated the high-pressure-limit rate coefficients of these H-shift reactions by the multi-structural variational transition-state theory with the small-curvature tunneling approximation (MS-CVT/SCT) and the thermochemical properties of all of the involved key radicals by the multi-structural torsional (MS-T) anharmonicity approximation method.

Model Chemistry Recommendations for Scaled Harmonic Frequency Calculations: A Benchmark Study

Despite the widespread popularity of scaled harmonic frequency calculations to predict experimental fundamental frequencies in chemistry, sparse benchmarking is available to guide users on the appropriate level of theory and basis set choices (model chemistry) or deep understanding of expected errors. An updated assessment of the best approach for scaling to minimize errors is also overdue. Here, we assess the performance of over 600 popular, contemporary, and robust model chemistries in the calculation of scaled harmonic frequencies, evaluating different scaling factor types and their implications in the scaled harmonic frequencies and model chemistry performance. We can summarize our results into three main findings: (1) Using model-chemistry-specific scaling factors optimized for three different frequency regions (low (<1,000 cm^-1), mid (1,000-2,000 cm^-1), and high (>2,000 cm^-1)) results in substantial improvements in the agreement between the scaled harmonic and experimental frequencies compared to other choices. (2) Larger basis sets and more robust levels of theory generally lead to superior performance; however, the particular model chemistry choice matters and poor choices lead to significantly reduced accuracies. (3) Outliers are expected in routine calculations regardless of the model chemistry choice. Our benchmarking results here do not consider the intensity of vibrational transitions; however, we draw upon previous benchmarking results for dipole moments that highlight the importance of diffuse functions (i.e., augmented basis sets) in high-quality intensity predictions. In terms of specific recommendations, overall, the highest accuracy model chemistries are double-hybrid density functional approximations with a non-Pople augmented triple-ζ basis set, which can produce median frequency errors down to 7.6 cm^-1 (DSD-PBEP86/def2-TZVPD), which is very close to the error in the harmonic approximation, i.e., the anharmonicity error. Double-ζ basis sets should not be used with double-hybrid functionals as there is no improvement compared to hybrid functional results (unlike for double-hybrid triple-ζ model chemistries). Note that 6-311G* and 6-311+G* basis sets perform like a double-ζ basis set for vibrational frequencies. After scaling, all studied hybrid functionals with non-Pople triple-ζ basis sets will produce median errors of less than 15 cm^-1, with the best result of 9.9 cm^-1 with B97-1/def2-TZVPD. Appropriate matching of double-ζ basis sets with hybrid functionals can produce high-quality results, but the precise choice of functional and basis set is more important. The B97-1, TPSS0-D3(BJ), or ωB97X-D hybrid density functionals with 6-31G*, pc-1, or pcseg-1 are recommended for fast routine calculations, all delivering median errors of 11-12 cm^-1. Note that dispersion corrections are not easily available for B97-1; given its strong performance here, we recommend these be added to major programs in coming updates.

New Reactions for the Formation of Organic Nitrate in the Atmosphere

Organic nitrates make an important contribution to the formation of secondary organic aerosols, but the formation mechanisms of organic nitrates are not fully understood at the molecular level. In the present work, we explore a new route for the formation of organic nitrates in the reaction of formaldehyde (HCHO) with nitric acid (HNO₃) catalyzed by water (H₂O), ammonia (NH₃), and dimethylamine ((CH₃)₂NH) using theoretical methods. The present results using CCSD(T)-F12a/cc-pVTZ-F12//M06-2X/MG3S unravel that dimethylamine has a stronger catalytic ability in the reaction of HCHO with HNO₃, reducing the barrier by 21.97 kcal/mol, while water and ammonia only decrease the energy barrier by 7.35 and 13.56 kcal/mol, respectively. In addition, the calculated kinetics combined with the corresponding concentrations of these species show that the HCHO + HNO₃ + (CH₃)₂NH reaction can compete well with the naked HCHO + HNO₃ reaction at 200-240 K, which may make certain contributions to the formation of organic nitrates under some atmospheric conditions.

Quantitative Kinetics of HO2 Reactions with Aldehydes in the Atmosphere: High-Order Dynamic Correlation, Anharmonicity, and Falloff Effects Are All Important

Kinetics provides the fundamental parameters for elucidating sources and sinks of key atmospheric species and for atmospheric modeling more generally. Obtaining quantitative kinetics in the laboratory for the full range of atmospheric temperatures and pressures is quite difficult. Here, we use computational chemistry to obtain quantitative rate constants for the reactions of HO₂ with HCHO, CH₃CHO, and CF₃CHO. First, we calculate the high-pressure-limit rate constants by using a dual-level strategy that combines conventional transition state theory using a high level of electronic structure wave function theory with canonical variational transition state theory including small-curvature tunneling using density functional theory. The wave-function level is beyond-CCSD(T) for HCHO and CCSD(T)-F12a (Level-A) for XCHO (X = CH₃, CF₃), and the density functional (Level-B) is specifically validated for these reactions. Then, we calculate the pressure-dependent rate constants by using system-specific quantum RRK theory (SS-QRRK) and also by an energy-grained master equation. The two treatments of the pressure dependence agree well. We find that the Level-A//Level-B method gives good agreement with CCSDTQ(P)/CBS. We also find that anharmonicity is an important factor that increases the rate constants of all three reactions. We find that the HO₂ + HCHO reaction has a significant dependence on pressure, but the HO₂ + CF₃CHO reaction is almost independent of pressure. Our findings show that the HO₂ + HCHO reaction makes important contribution to the sink for HCHO, and the HO₂ + CF₃CHO reaction is the dominant sink for CF₃CHO in the atmosphere.

The reaction of Criegee intermediates with formamide and its implication to atmospheric aerosols

The reactions of Criegee intermediates (CIs) play an important role in the formation of secondary organic aerosols that have negative effect on visibility, human health, and global climate. New particle formation (NPF) can contribute to more than half of the aerosols in terms of their number concentration. Here, the reactions of CIs with formamide (FA) in the gas-phase and at the air/water interface were investigated using quantum chemistry calculation and Born-Oppenheimer molecular dynamic simulations. The results show that the reaction mechanism of CIs with FA is similar to that with formic acid, and the formation of hydroperoxymethyl formimidate (P4) is the most favorable pathway both in the gas-phase and at the air/water interface. Moreover, the potential contribution of the products to NPF was also evaluated by means of the molecular dynamic simulations. The results indicate that the product (P4) can participate in the SA-based NPF and water molecules are beneficial to enhance the NPF. The exploration will provide insight into the reaction of CIs with amide and the effect of the Criegee chemistry on the atmospheric aerosols.

Temperature-dependent kinetics of the atmospheric reaction between CH2OO and acetone

Criegee intermediates are important oxidants produced in the ozonolysis of alkenes in the atmosphere. Quantitative kinetics of the reactions of Criegee intermediates are required for atmospheric modeling. However, the experimental studies do not cover the full relevant range of temperature and pressure. Here we report the quantitative kinetics of CH₂OO + CH₃C(O)CH₃ by using our recently developed dual strategy that combines coupled cluster theory with high excitation levels for conventional transition state theory and well validated levels of density functional theory for direct dynamics calculations using canonical variational transition theory including tunneling. We find that the W3X-L//DF-CCSD(T)-F12b/jun-cc-pVDZ electronic structure method can be used to obtain quantitative kinetics of the CH₂OO + CH₃C(O)CH₃ reaction. Whereas previous investigations considered a one-step mechanistic pathway, we find that the CH₂OO + CH₃C(O)CH₃ reaction occurs in a stepwise manner. This has implications for the modeling of Criegee-intermediate reactions with other ketones and with aldehydes. In the kinetics calculations, we show that recrossing effects of the conventional transition state are negligible for determining the rate constant of CH₂OO + CH₃C(O)CH₃. The present findings reveal that the rate ratio between CH₂OO + CH₃C(O)CH₃ and OH + CH₃C(O)CH₃ has a significant negative dependence on temperature such that the CH₂OO + CH₃C(O)CH₃ reaction can contribute as a significant sink for atmospheric CH₃C(O)CH₃ at low temperature. The present findings should have broad implications in understanding the reactions of Criegee intermediates with carbonyl compounds and ketones in the atmosphere.

Large Pressure Effects Caused by Internal Rotation in the s-cis-syn-Acrolein Stabilized Criegee Intermediate at Tropospheric Temperature and Pressure

Criegee intermediates are important atmospheric oxidants, and quantitative kinetics for stabilized Criegee intermediates are key parameters for atmospheric modeling but are still limited. Here we report barriers and rate constants for unimolecular reactions of s-cis-syn-acrolein oxide (scsAO), in which the vinyl group makes it a prototype for Criegee intermediates produced in the ozonolysis of isoprene. We find that the MN15-L and M06-2X density functionals have CCSD(T)/CBS accuracy for the unimolecular cyclization and stereoisomerization of scsAO. We calculated high-pressure-limit rate constants by the dual-level strategy that combines (a) high-level wave function-based conventional transition-state theory (which includes coupled-cluster calculations with quasiperturbative inclusion of quadruple excitations because of the strongly multiconfigurational character of the electronic wave function) and (b) canonical variational transition-state theory with small-curvature tunneling based on a validated density functional. We calculated pressure-dependent rate constants both by system-specific quantum Rice-Ramsperger-Kassel theory and by solving the master equation. We report rate constants for unimolecular reactions of scsAO over the full range of atmospheric temperature and pressure. We found that the unimolecular reaction rates of this larger-than-previously studied Criegee intermediate depend significantly on pressure. Particularly, we found that falloff effects decrease the effective unimolecular cyclization rate constant of scsAO by about a factor of 3, but the unimolecular reaction is still the dominant atmospheric sink for scsAO at low altitudes. The large falloff caused by the inclusion of the stereoisomerization channel in the master equation calculations has broad implications for mechanistic analysis of reactions with competitive internal rotations that can produce stable rotamers.

CARNOT: a Fragment-Based Direct Molecular Dynamics and Virtual-Reality Simulation Package for Reactive Systems

Traditionally, the study of reaction mechanisms of complex reaction systems such as combustion has been performed on an individual basis by optimizations of transition structure and minimum energy path or by reaction dynamics trajectory calculations for one elementary reaction at a time. It is effective, but time-consuming, whereas important and unexpected processes could have been missed. In this article, we present a direct molecular dynamics (DMD) approach and a virtual-reality simulation program, CARNOT, in which plausible chemical reactions are simulated simultaneously at finite temperature and pressure conditions. A key concept of the present ab initio molecular dynamics method is to partition a large, chemically reactive system into molecular fragments that can be adjusted on the fly of a DMD simulation. The theory represents an extension of the explicit polarization method to reactive events, called ReX-Pol. We propose a highest-and-lowest adapted-spin approximation to define the local spins of individual fragments, rather than treating the entire system by a delocalized wave function. Consequently, the present ab initio DMD can be applied to reactive systems consisting of an arbitrarily varying number of closed and open-shell fragments such as free radicals, zwitterions, and separate ions found in combustion and other reactions. A graph-data structure algorithm was incorporated in CARNOT for the analysis of reaction networks, suitable for reaction mechanism reduction. Employing the PW91 density functional theory and the 6-31+G(d) basis set, the capabilities of the CARNOT program were illustrated by a combustion reaction, consisting of 28 650 atoms, and by reaction network analysis that revealed a range of mechanistic and dynamical events. The method may be useful for applications to other types of complex reactions.

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.