Active Thermochemical Tables: Water and Water Dimer

Composition of the Water Dimer and the Heterodimers of Water with N2 and O2 in Earth's Atmosphere

The mole fractions χ and number concentrations n of the water dimer and the heterodimers H2O-N2 and H2O-O2 in Earth's atmosphere are reported up to 20 km. The water dimer data is obtained from published values of the equilibrium constant based on the water equation of state. The mixed equilibrium constants for the heterodimers are obtained from the respective second virial coefficients using an approach introduced by Stogryn and Hirschfelder that extracts the components pertaining to pairwise interactions producing bound and metastable dimers. From these calculations, χ and n for the water dimer and the (H2O)(N2) and (H2O)(O2) heterodimers at standard sea level are 1.79(6) × 10-5, 4.77(12) × 10-5 and 9.90(5) × 10-6 and 4.55(15) × 1014, 1.23(3) × 1016 and 2.56(1) × 1015, respectively. Analytical expressions are provided for these quantities for altitudes between 0-20 km and temperatures from 200-300 K. Sea level values of χ and n are given for two specific locations.

Ab initio spectroscopy and thermochemistry of the platinum hydride ions, PtH+ and PtH

Rovibrational levels of low-lying electronic states of the gas-phase, diatomic molecules, PtH+ and PtH-, are computed on potential-energy functions obtained by using a hybrid spin-orbit configuration-interaction procedure. PtH- has a well-separated Σ0++1 ground state, while the first two electronic states of PtH+ (Σ0++1 and 3Δ3) are nearly degenerate. Combining the experimental photoelectron (PE) spectra of PtH- with theoretical photodetachment spectroscopy leads to an improved value for the electron affinity of PtH, EA(PtH) = (1.617 ± 0.015) eV. When PtH- is a product of photodissociation of PtHCO2-, its PE spectrum is broad because of rotational excitation. Temperature-dependent thermodynamic functions and thermochemistry of dissociation are computed from the theoretical energy levels. Previously published energetic quantities for PtH+ and PtH- are revised. The ground 1Σ+ term of PtH+ is not well described using single-reference theory.

Kinetics of the Gas-Phase Reactions of syn- and anti-CH3CHOO Criegee Intermediate Conformers with SO2 as a Function of Temperature and Pressure

Kinetics of reactions between SO2 and CH3CHOO Criegee intermediate conformers have been measured at temperatures between 242 and 353 K and pressures between 10 and 600 Torr using laser flash photolysis of CH3CHI2/O2/N2/SO2 gas mixtures coupled with time-resolved broadband UV absorption spectroscopy. The kinetics of syn-CH3CHOO + SO2 are pressure-dependent and exhibit a negative temperature dependence, with the observed pressure dependence reconciling apparent discrepancies between previous measurements performed at ∼298 K. Results indicate a rate coefficient of (4.80 ± 0.46) × 10-11 cm3 s-1 for the reaction of syn-CH3CHOO with SO2 at 298 K and 760 Torr. In contrast to the behavior of the syn-conformer, the kinetics of anti-CH3CHOO + SO2 display no significant dependence on temperature or pressure over the ranges investigated, with a mean rate coefficient of (1.18 ± 0.21) × 10-10 cm3 s-1 over all conditions studied in this work. Results indicate that the reaction of syn-CH3CHOO with SO2 competes with unimolecular decomposition and reaction with water vapor in areas with high SO2 concentration and low humidity, particularly at lower temperatures.

Effect of an Electric Field on the Structure and Stability of Atmospheric Clusters

We study the influence of an applied electric field on the structure and stability of some common bimolecular clusters that are found in the atmosphere. These clusters play an important role in new particle formation (NPF). For low values of the electric field (i.e., |E| ≤ 0.01 V Å-1), we demonstrate that the field response of the clusters can be predicted from simply calculating the dipole moment of the cluster and the dipole moments of the constituent molecules and that the influence on the association energy of the cluster is minimal (i.e., <0.5 kcal mol-1). For higher field strengths |E| > 0.2 V Å-1, there can be more dramatic effects on both structure and energetics, as the induced dipole, charge transfer, and geometric distortion play a larger role. Although such large fields are not very relevant in the atmosphere, they do exist in some situations of experimental interest, such as near interfaces and in intense laser fields.

Thermodynamics of the van der Waals Dimers of O2, N2 and the Heterodimer (N2)(O2) and Their Presence in Earth's Atmosphere

The dimerization thermodynamics of N2 and O2, the principal components of Earth's atmosphere, have been determined from the respective second virial coefficients of the bound and metastable dimers calculated using the method of Stogryn and Hirschfelder that utilizes the Lennard-Jones (LJ) potential to account for intermolecular interactions. In addition, the thermodynamic properties of the heterodimer (N2)(O2) have been obtained using the same approach, employing combining rules to construct the LJ potential. Thus, Keq, ΔH, and ΔS for the three dimers are reported between 80-120 K. Over this temperature range, the ranking of Keq is (N2)(O2) > (O2)(O2) > (N2)(N2). The same trend is found for the exoethalpicity of dimer formation. For example, at 100 K, the Keq values are, respectively, 0.0406(14), 0.0215(5), and 0.0181(10), and the corresponding ΔH values are -2401(5), -2344(7), and -2279(1) J/mol. The mole fraction composition of the dimers in the atmosphere was calculated for altitudes up to 20 km. These calculations show that in the troposphere and the lower stratosphere (up to 20 km), the three dimers rank fifth to seventh in abundance, between CO2 and Ne. In this region, the average mole fractions of (N2)(N2), (O2)(O2), and (N2)(O2) are calculated to be 3.4(2) × 10-4, 2.80(9) × 10-5, and 1.95(7) × 10-4, respectively.

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.

The simplest Criegee intermediate CH2OO reaction with dimethylamine and trimethylamine: kinetics and atmospheric implications

We have used the OH laser-induced fluorescence (LIF) method to measure the kinetics of the simplest Criegee intermediate (CH2OO) reacting with two abundant amines in the atmosphere: dimethylamine ((CH3)2NH) and trimethylamine ((CH3)3N). Our experiments were conducted under pseudo-first-order approximation conditions. The rate coefficients we report are (2.15 ± 0.28) × 10-11 cm3 molecule-1 s-1 for (CH3)2NH at 298 K and 10 Torr, and (1.56 ± 0.23) × 10-12 cm3 molecule-1 s-1 for (CH3)3N at 298 K and 25 Torr with Ar as the bath gas. Both reactions exhibit a negative temperature dependence. The activation energy and pre-exponential factors derived from the Arrhenius equation were (-2.03 ± 0.26) kcal mol-1 and (6.89 ± 0.90) × 10-13 cm3 molecule-1 s-1 for (CH3)2NH, and (-1.60 ± 0.24) kcal mol-1 and (1.06 ± 0.16) × 10-13 cm3 molecule-1 s-1 for (CH3)3N. We propose that the electronegativity of the atom in the co-reactant attached to the C atom of CH2OO, in addition to the dissociation energy of the fragile covalent bonds with H atoms (H-X bond), plays an important role in the 1,2-insertion reactions. Under certain circumstances, the title reactions can contribute to the sink of amines and Criegee intermediates and to the formation of secondary organic aerosol (SOA).

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.

Analysis of the Electron Density of a Water Molecule Encapsulated by Two Cholic Acid Residues

Cholic acid is a trihydroxy bile acid with a nice peculiarity: the average distance between the oxygen atoms (O7 and O12) of the hydroxy groups located at C7 and C12 carbon atoms is 4.5 Å, a value which perfectly matches with the O/O tetrahedral edge distance in Ih ice. In the solid phase, they are involved in the formation of hydrogen bonds with other cholic acid units and solvents. This fact was satisfactorily used for designing a cholic dimer which encapsulates one single water molecule between two cholic residues, its oxygen atom (Ow) being exactly located at the centroid of a distorted tetrahedron formed by the four steroid hydroxy groups. The water molecule participates in four hydrogen bonds, with the water simultaneously being an acceptor from the 2 O12 (hydrogen lengths are 2.177 Å and 2.114 Å) and a donor towards the 2 O7 (hydrogen bond lengths are 1.866 Å and 1.920 Å). These facts suggest that this system can be a nice model for the theoretical study of the formation of ice-like structures. These are frequently proposed to describe the water structure found in a plethora of systems (water interfaces, metal complexes, solubilized hydrophobic species, proteins, and confined carbon nanotubes). The above tetrahedral structure is proposed as a reference model for those systems, and the results obtained from the application of the atoms in molecules theory are presented here. Furthermore, the structure of the whole system allows a division into two interesting subsystems in which water is the acceptor of one hydrogen bond and the donor of another. The analysis of the calculated electron density is performed through its gradient vector and the Laplacian. The calculation of the complexation energy used correction of the basis set superposition error (BSSE) with the counterpoise method. As expected, four critical points located in the H…O bond paths were identified. All calculated parameters obey the proposed criteria for hydrogen bonds. The total energy for the interaction in the tetrahedral structure is 54.29 kJ/mol, while the summation obtained of the two independent subsystems and the one between the alkyl rings without water is only 2.5 kJ/mol higher. This concordance, together with the calculated values for the electron density, the Laplacian of the electron density, and the lengths of the oxygen atom and the hydrogen atom (involved in the formation of each hydrogen bond) to the hydrogen bond critical point, suggests that each pair of hydrogen bonds can be considered independent of each other.

Active Thermochemical Tables: Enthalpies of Formation of Bromo- and Iodo-Methanes, Ethenes and Ethynes

The thermochemistry of halocarbon species containing iodine and bromine is examined through an extensive interplay between new Feller-Peterson-Dixon (FPD) style composite methods and a detailed analysis of all available experimental and theoretical determinations using the thermochemical network that underlies the Active Thermochemical Tables (ATcT). From the computational viewpoint, a slower convergence of the components of composite thermochemistry methods is observed relative to species that solely contain first row elements, leading to a higher computational expense for achieving comparable levels of accuracy. Potential systematic sources of computational uncertainty are investigated, and, not surprisingly, spin-orbit coupling is found to be a critical component, particularly for iodine containing molecular species. The ATcT analysis of available experimental and theoretical determinations indicates that prior theoretical determinations have significantly larger uncertainties than originally reported, particularly in cases where molecular spin-orbit effects were ignored. Accurate and reliable heats of formation are reported for 38 halogen containing systems, based on combining the current computations with previous experimental and theoretical work via the ATcT approach.

Equilibrium Thermodynamics of the Dimers and Trimers of H2 and D2 and Their Heterodimer (H2)(D2)

The equilibrium thermochemical properties of the dimers and trimers of H₂ and D₂ are obtained from the equations of state (EOSs) of normal H₂ and D₂. The standard dimer and trimer equilibrium constants K _D ^o and K _T ^o and ΔH _D, T ^o and ΔS _D, T ^o are reported for these weakly bound van der Waals molecules between 25 and 45 K. Statistical thermodynamics (ST) calculations of H₂ and D₂ dimerization using Morse pair potentials to account for intermolecular interactions, obtained from recent experimental work, are in qualitative agreement with the EOS results. The entropies of the H₂ and D₂ dimers and trimers are calculated from the EOS ΔS ^o values and ST calculations of the monomer entropies.

Vibrational Spectroscopy of the Water Dimer at Jet-Cooled and Atmospheric Temperatures

The vibrational spectroscopy of the water dimer provides an understanding of basic hydrogen bonding in water clusters, and with about one water dimer for every 1,000 water molecules, it plays a critical role in atmospheric science. Here, we review how the experimental and theoretical progress of the past decades has improved our understanding of water dimer vibrational spectroscopy under both cold and warm conditions. We focus on the intramolecular OH-stretching transitions of the donor unit, because these are the ones mostly affected by dimer formation and because their assignment has proven a challenge. We review cold experimental results from early matrix isolation to recent mass-selected jet expansion techniques and, in parallel, the improvements in the theoretical anharmonic models. We discuss and illustrate changes in the vibrational spectra of complexes upon increasing temperature, and the difficulties in recording and calculating these spectra. In the atmosphere, water dimer spectra at ambient temperature are crucial. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Thermodynamics of Helium-4 Dimerization and Trimerization

A new equation of state (EOS) for helium-4 is used to obtain the equilibrium thermochemical properties of helium-4 dimerization (2⁴He ⇌ ⁴He₂) and trimerization (3⁴He ⇌ ⁴He₃) between 3.0 and 10.0 K. It is shown that at sufficiently low temperatures there are appreciable populations of dimer and trimer. The calculations account only for monomer, dimer, and trimer. At 3.0 K, the respective K_P^o values for dimerization and trimerization are 0.4832 and 0.4876, respectively. The standard enthalpy changes at 3.0 K are -54.53 and -110.0 J/mol, and standard entropy changes are -24.22 and -42.97 J/mol-K. Statistical thermodynamic calculations provide results that are qualitatively consistent with those obtained from the EOS calculations.

Molybdenum Disorder in Hydrated Sedovite, Ideally U(MoO4)2·nH2O, a Microporous Nanocrystalline Mineral Characterized by Three-Dimensional Electron Diffraction, Density Functional Theory Computations, and Complexity Analysis

Sedovite, U⁴⁺(Mo⁶⁺O₄)₂·nH₂O, is reported as being one of the earliest supergene minerals formed of the secondary zone. The difficulty of isolating enough pure material limits studies to techniques that can access the nanoscale combined with theoretical analyses. The crystal structure of sedovite has been solved and refined using the dynamical approach from three-dimensional electron diffraction data collected on natural nanocrystals found among iriginite. At 100 K, sedovite is monoclinic a ≈ 6.96 Å, b ≈ 9.07 Å, c ≈ 12.27 Å, and V ≈ 775 ų with space group C2/c. The microporous structure presents a characteristic framework built from uranium polyhedra and disordered Mo pyramids creating pore hosting water molecules. To confirm the formula U⁴⁺(Mo⁶⁺O₄)₂·nH₂O, the possible presence of a hydroxyl group that would promote Mo⁵⁺ was tested with density functional theory (DFT) computations at the ambient temperature. DFT predicts that sedovite is a ferromagnetic insulator with a fundamental bandgap of E_g ∼ 1.7 eV with its chemical and physical properties dominated by U⁴⁺ rather than Mo⁶⁺. The structural complexity, I_G,tot, of sedovite was evaluated in order to get indirect information about the missing formation conditions.

Nonlinear rotational spectroscopy reveals many-body interactions in water molecules

Because of their central importance in chemistry and biology, water molecules have been the subject of decades of intense spectroscopic investigations. Rotational spectroscopy of water vapor has yielded detailed information about the structure and dynamics of isolated water molecules, as well as water dimers and clusters. Nonlinear rotational spectroscopy in the terahertz regime has been developed recently to investigate the rotational dynamics of linear and symmetric-top molecules whose rotational energy levels are regularly spaced. However, it has not been applied to water or other lower-symmetry molecules with irregularly spaced levels. We report the use of recently developed two-dimensional (2D) terahertz rotational spectroscopy to observe high-order rotational coherences and correlations between rotational transitions that were previously unobservable. The results include two-quantum (2Q) peaks at frequencies that are shifted slightly from the sums of distinct rotational transitions on two different molecules. These results directly reveal the presence of previously unseen metastable water complexes with lifetimes of 100 ps or longer. Several such peaks observed at distinct 2Q frequencies indicate that the complexes have multiple preferred bimolecular geometries. Our results demonstrate the sensitivity of rotational correlations measured in 2D terahertz spectroscopy to molecular interactions and complexation in the gas phase.

Kinetics of the gas phase reaction of the Criegee intermediate CH2OO with SO2 as a function of temperature

The kinetics of the gas phase reaction of the Criegee intermediate CH₂OO with SO₂ have been studied as a function of temperature in the range 223-344 K at 85 Torr using flash photolysis of CH₂I₂/O₂/SO₂/N₂ mixtures at 248 nm coupled to time-resolved broadband UV absorption spectroscopy. Measurements were performed under pseudo-first-order conditions with respect to SO₂, revealing a negative temperature dependence. Analysis of experimental results using the Master Equation Solver for Multi-Energy well Reactions (MESMER) indicates that the observed temperature dependence, combined with the reported lack of a pressure dependence in the range 1.5-760 Torr, can be described by a reaction mechanism consisting of the formation of a pre-reaction complex leading to a cyclic secondary ozonide which subsequently decomposes to produce HCHO + SO₃. The temperature dependence can be characterised by k_CH2OO+SO2 = (3.72 ± 0.13) × 10^-11 (T/298)^{(-2.05±0.38)} cm³ molecule^-1 s^-1. The observed negative temperature dependence for the title reaction in conjunction with the decrease in water dimer (the main competitor for the Criegee intermediate) concentration at lower temperatures means that Criegee intermediate chemistry can play an enhanced role in SO₂ oxidation in the atmosphere at lower temperatures.

Micro- vs Macrosolvation in Reichardt's Dyes

Solvation is a complex phenomenon involving electrostatic and van der Waals forces as well as chemically more specific effects such as hydrogen bonding. To disentangle global solvent effects (macrosolvation) from local solvent effects (microsolvation), we studied the UV-vis and IR spectra of a solvatochromic pyridinium-N-phenolate dye (a derivative of Reichardt's dye) in rare gas matrices, in mixtures of argon and water, and in water ice. The π-π* transition of the betaine dye in the visible region and its C-O stretching vibration in the IR region are highly sensitive to solvent effects. By annealing argon matrices of the betaine dye doped with low concentrations of water, we were able to synthesize 1:1 water-dye complexes. Formation of hydrogen-bonded complexes leads to small shifts of the π-π* transition only, as long as the global polarity of the matrix environment does not change. In contrast, changes of the global polarity result in large spectral band shifts. Hydrogen-bonded complexes of the betaine dye are more sensitive to global polarity changes than the dye itself, explaining why E_T values determined with Reichardt's dyes are very different for protic and nonprotic solvents, even if the relative permittivities of these solvents are similar.

The favorable routes for the hydrolysis of CH2OO with (H2O)n (n = 1-4) investigated by global minimum searching combined with quantum chemical methods

The hydrolysis reaction of CH₂OO with water and water clusters is believed to be a dominant sink for the CH₂OO intermediate in the atmosphere. However, the favorable route for the hydrolysis of CH₂OO with water clusters is still unclear. Here global minimum searching using the Tsinghua Global Minimum program has been introduced to find the most stable geometry of the CH₂OO(H₂O)_n (n = 1-4) complex firstly. Then, based on these stable complexes, favorable hydrolysis of CH₂OO with (H₂O)_n (n = 1-4) has been investigated using the quantum chemical method of CCSD(T)-F12a/cc-pVDZ-F12//B3LYP/6-311+G(2d,2p) and canonical variational transition state theory with small curvature tunneling. The calculated results have revealed that, although the contribution of CH₂OO + (H₂O)₂ is the most obvious in the hydrolysis of CH₂OO with (H₂O)_n (n = 1-4), the hydrolysis of CH₂OO with (H₂O)₃ is not negligible in atmospheric gas-phase chemistry as its rate is close to the rate of the CH₂OO + H₂O reaction. The calculated results also show that, in a clean atmosphere, the CH₂OO + (H₂O)_n (n = 1-2) reaction competes well with the CH₂OO + SO₂ reaction at 298 K when the concentrations of (H₂O)_n (n = 1-2) range from 20% relative humidity (RH) to 100% RH, and SO₂ is 2.46 × 10¹¹ molecules per cm³. Meanwhile, when the RH is higher than 40%, it is a new prediction that the CH₂OO + (H₂O)₃ reaction can also compete well with the CH₂OO + SO₂ reaction at 298 K. Besides, Born-Oppenheimer molecular dynamics simulation results show that all the favorable channels of the CH₂OO + (H₂O)_n (n = 1-3) reaction cannot react on a time scale of 100 ps in the NVT simulation. However, the NVE simulation results show that the CH₂OO + (H₂O)₃ reaction can be finished well at 8.5 ps, indicating that the gas phase reaction of CH₂OO + (H₂O)₃ is not negligible in the atmosphere. Overall, the present results have provided a definitive example of how the favorable hydrolysis of important atmospheric species with (H₂O)_n (n = 1-4) takes place, which will stimulate one to consider the favorable hydrolysis of water and water clusters with other Criegee intermediates and other important atmospheric species.

Kinetics of Unimolecular Decay of Methyl Vinyl Ketone Oxide, an Isoprene-Derived Criegee Intermediate, under Atmospherically Relevant Conditions

Isoprene is the most abundant unsaturated hydrocarbon in the atmosphere. Ozonolysis of isoprene produces methyl vinyl ketone oxide (MVKO), which may react with atmospheric SO₂, formic acid, and other important species at substantial levels. In this study, we utilized ultraviolet absorption to monitor the unimolecular decay kinetics of syn-MVKO in real time at 278-319 K and 100-503 Torr. After removing the contributions of radical reactions and wall loss, the unimolecular decay rate coefficient of syn-MVKO was measured to be k_uni = 70 ± 15 s^-1 (1σ uncertainty) at 298 K with negligible pressure dependence. In addition, k_uni increases from ca. 30 s^-1 at 278 K to ca. 175 s^-1 at 319 K with an effective Arrhenius activation energy of 8.3 ± 2.5 kcal mol^-1, k_uni(T) = (9.3 × 10⁷)exp(-4200/T) s^-1. Our results indicate that unimolecular decay is the major sink of MVKO in the troposphere. The data would improve the estimation for the steady-state concentrations of MVKO and thus its oxidizing ability.

Nuclear quantum effects on the hydrogen bond donor-acceptor exchange in water-water and water-methanol dimers

We present results from path integral molecular dynamics simulations that describe effects from the explicit incorporation of nuclear quantum fluctuations on the topology of the free energy associated with the geared exchange of hydrogen bonds in the water-water dimer. Compared to the classical treatment, our results reveal important reductions in the free energy barriers and changes at a qualitative level in the overall profile. Most notable are those manifested by a plateau behavior, ascribed to nuclear tunneling, which bridges reactant and product states, contrasting with the usual symmetric double-well profile. The characteristics of the proton localizations along the pathway are examined. An imaginary time analysis of the rotational degrees of freedom of the partners in the dimer at the vicinities of transition states shows a clear "anticorrelation" between intermolecular interactions coupling beads localized in connective and dangling basins of attractions. As such, the transfer is operated by gradual concerted inter-basin migrations in opposite directions, at practically no energy costs. Modifications operated by partial deuteration and by the asymmetries in the hydrogen bonding characteristics prevailing in water-methanol heterodimers are also examined.

Density-functional-theory-predicted symmetry lowering from cubic to tetragonal in nickel hexacyanoferrate

Density-functional-theory (DFT) computations on a Prussian blue analogue (PBA), nickel hexacyanoferrate, Ni2+ 3[Fe3+(CN)6]2·nH2O, predict the existence of a tetragonal (P 4 m2) crystal structure that is energetically degenerate with the previously reported cubic (F 43m) structure for this PBA. The proposed tetragonal structure satisfies observations, such as X-ray diffraction and magnetic measurements, that have been reported previously. A van der Waals corrected exchange-correlation functional is used in the DFT+U computations for an improved description of hydrogen bonding. The results provide strong support for a revised and simplified crystallographic description of Ni2+ 3[Fe3+(CN)6]2·nH2O, and show how H2O molecules stabilize the crystal structure and affect its magnetic and electronic properties. The symmetry lowering in nickel hexacyanoferrate is attributed to the hydration shell of the interstitial nickel cation. Calculations strongly suggest a maximum of n = 7 interstitial H2O molecules per formula unit for nickel hexacyanoferrate at room temperature, and a higher water content at temperatures below T ≃ 200 K. Since the symmetry lowering relies on the presence of interstitial H2O molecules, this revised crystallographic description may be applicable more generally to the large class of F 43m-structured PBAs.

Effects of water vapor on the reaction of CH2OO with NH3

A strong synergic effect of water and ammonia molecules may enhance the formation of H₂NCH₂OOH.

Effect of unsaturated substituents in the reaction of Criegee intermediates with water vapor

Criegee intermediates (CIs), formed in the reactions of unsaturated hydrocarbons with ozone, are very reactive carbonyl oxides and have recently been suggested as important oxidants in the atmosphere. In this work, we studied the substituent effect on the water monomer and dimer reaction with CIs which include up to three carbon atoms at the QCISD(T)/CBS//B3LYP/6-311+G(2d,2p) level. Our calculation showed that for saturated CIs with a hydrogen atom on the same side as the terminal oxygen atom, the reaction with water vapor would likely dominate the removal processes of these CIs in the atmosphere. On the other hand, for unsaturated CIs, the reactivity toward water vapor decreases compared to the saturated species allowing them to survive in humid atmospheric environments. We also evaluated the kinetic isotope effect in the reaction between CI and water vapor by performing calculations with deuterated water. We found that tunneling is not important and the kinetic isotope effect mainly comes from the difference in the zero point energy between water and deuterated water.

The reaction of Criegee intermediate CH2OO with water dimer: primary products and atmospheric impact

The rapid reaction of the smallest Criegee intermediate, CH₂OO, with water dimers is the dominant removal mechanism for CH₂OO in the Earth's atmosphere, but its products are not well understood. This reaction was recently suggested as a significant source of the most abundant tropospheric organic acid, formic acid (HCOOH), which is consistently underpredicted by atmospheric models. However, using time-resolved measurements of reaction kinetics by UV absorption and product analysis by photoionization mass spectrometry, we show that the primary products of this reaction are formaldehyde and hydroxymethyl hydroperoxide (HMHP), with direct HCOOH yields of less than 10%. Incorporating our results into a global chemistry-transport model further reduces HCOOH levels by 10-90%, relative to previous modeling assumptions, which indicates that the reaction CH₂OO + water dimer by itself cannot resolve the discrepancy between the measured and predicted HCOOH levels.

Impact of the water dimer on the atmospheric reactivity of carbonyl oxides

The reactions of twelve carbonyl oxides or Criegee intermediates with the water monomer and with the water dimer have been investigated employing high level theoretical methods. The study includes all possible carbonyl oxides arising from the isoprene ozonolysis and the methyl and dimethyl carbonyl oxides that originated from the reaction of ozone with several hydrocarbons. These reactions have great significance in the chemistry of the atmosphere because Criegee intermediates have recently been identified as important oxidants in the troposphere and as precursors of secondary organic aerosols. Moreover, water vapor is one of the most abundant trace gases in the atmosphere and the water dimer can trigger the atmospheric decomposition of Criegee intermediates. Our calculations show that the nature and position of the substituents in carbonyl oxides play a very important role in the reactivity of these species with both the water monomer and the water dimer. This fact results in differences in rate constants of up to six orders of magnitude depending on the carbonyl oxide. In this work we have defined an effective rate constant (keff) for the atmospheric reaction of carbonyl oxides with water vapor, which depends on the temperature and on the relative humidity as well. With this keff we show that the water dimer, despite its low tropospheric concentration, enhances the atmospheric reactivity of Criegee intermediates, but its effect changes with the nature of carbonyl oxide, ranging between 59 and 295 times in the most favorable case (syn-methyl carbonyl oxide), and between 1.4 and 3 times only in the most unfavorable case.

Mechanisms of photoexcitation and photoionization in small water clusters

The S₀ → S₁ excitation leads to strong polarization and formation of [(H₂O)₂]⁺˙ from which both photoexcited and photoionized products are generated.

Thermal Decomposition of Potential Ester Biofuels. Part I: Methyl Acetate and Methyl Butanoate

Two methyl esters were examined as models for the pyrolysis of biofuels. Dilute samples (0.06-0.13%) of methyl acetate (CH₃COOCH₃) and methyl butanoate (CH₃CH₂CH₂COOCH₃) were entrained in (He, Ar) carrier gas and decomposed in a set of flash-pyrolysis microreactors. The pyrolysis products resulting from the methyl esters were detected and identified by vacuum ultraviolet photoionization mass spectrometry. Complementary product identification was provided by matrix infrared absorption spectroscopy. Pyrolysis pressures in the pulsed microreactor were about 20 Torr and residence times through the reactors were roughly 25-150 μs. Reactor temperatures of 300-1600 K were explored. Decomposition of CH₃COOCH₃ commences at 1000 K, and the initial products are (CH₂═C═O and CH₃OH). As the microreactor is heated to 1300 K, a mixture of CH₂═C═O and CH₃OH, CH₃, CH₂═O, H, CO, and CO₂ appears. The thermal cracking of CH₃CH₂CH₂COOCH₃ begins at 800 K with the formation of CH₃CH₂CH═C═O and CH₃OH. By 1300 K, the pyrolysis of methyl butanoate yields a complex mixture of CH₃CH₂CH═C═O, CH₃OH, CH₃, CH₂═O, CO, CO₂, CH₃CH═CH₂, CH₂CHCH₂, CH₂═C═CH₂, HCCCH₂, CH₂═C═C═O, CH₂═CH₂, HC≡CH, and CH₂═C═O. On the basis of the results from the thermal cracking of methyl acetate and methyl butanoate, we predict several important decomposition channels for the pyrolysis of fatty acid methyl esters, R-CH₂-COOCH₃. The lowest-energy fragmentation will be a 4-center elimination of methanol to form the ketene RCH═C═O. At higher temperatures, concerted fragmentation to radicals will ensue to produce a mixture of species: (RCH₂ + CO₂ + CH₃) and (RCH₂ + CO + CH₂═O + H). Thermal cracking of the β C-C bond of the methyl ester will generate the radicals (R and H) as well as CH₂═C═O + CH₂═O. The thermochemistry of methyl acetate and its fragmentation products were obtained via the Active Thermochemical Tables (ATcT) approach, resulting in Δ_fH₂₉₈(CH₃COOCH₃) = -98.7 ± 0.2 kcal mol^-1, Δ_fH₂₉₈(CH₃CO₂) = -45.7 ± 0.3 kcal mol^-1, and Δ_fH₂₉₈(COOCH₃) = -38.3 ± 0.4 kcal mol^-1.

A vacuum ultraviolet laser pulsed field ionization-photoion study of methane (CH4): determination of the appearance energy of methylium from methane with unprecedented precision and the resulting impact on the bond dissociation energies of CH4and CH4+

High-resolution VUV laser PFI-PI detection method for the study of quantum-state-selected unimolecular ion dissociation.

A kinetic study of the CH2OO Criegee intermediate reaction with SO2, (H2O)2, CH2I2 and I atoms using OH laser induced fluorescence

The OH laser induced fluorescence method was used to study the kinetics of CH₂OO reacting with SO₂, (H₂O)₂, CH₂I₂ and I atoms.

Communication: Interaction of BrO radical with the surface of water

Solvation of a BrO radical in a slab of water is investigated using adaptive buffered force quantum mechanics/molecular mechanics (QM/MM) dynamics simulations. The simulation results show that the BrO radical exhibits preference towards the water surface with respect to the interior region of the water slab, despite BrO's high affinity to water. Another important finding is the weakening of (BrO)Br⋯O(water) interaction at the water surface due to competitive interactions between (BrO)Br⋯O(water) and (water)H⋯O(water). As such, the BrO-water slab interaction is dominated by (BrO)O⋯H(water) interaction, contrary to that in the gas phase, suggesting that the reactive site for the BrO radical at the air/water surface is more likely the Br site. The conclusion from this study can offer deeper insight into the reactivity of the BrO radical at the air/water interface, with regard to atmospheric implications.