Ground and electronically excited states of methyl hydroperoxide: Comparison with hydrogen peroxide

Deciphering the Influence of Ground-State Distributions on the Calculation of Photolysis Observables

Nonadiabatic molecular dynamics offers a powerful tool for studying the photochemistry of molecular systems. Key to any nonadiabatic molecular dynamics simulation is the definition of its initial conditions (ICs), ideally representing the initial molecular quantum state of the system of interest. In this work, we provide a detailed analysis of how ICs may influence the calculation of experimental observables by focusing on the photochemistry of methylhydroperoxide (MHP), the simplest and most abundant organic peroxide in our atmosphere. We investigate the outcome of trajectory surface hopping simulations for distinct sets of ICs sampled from different approximate quantum distributions, namely harmonic Wigner functions and ab initio molecular dynamics using a quantum thermostat (QT). Calculating photoabsorption cross-sections, quantum yields, and translational kinetic energy maps from the results of these simulations reveals the significant effect of the ICs, in particular when low-frequency (∼ a few hundred cm-1) normal modes are connected to the photophysics of the molecule. Overall, our results indicate that sampling ICs from ab initio molecular dynamics using a QT is preferable for flexible molecules with photoactive low-frequency modes. From a photochemical perspective, our nonadiabatic dynamics simulations offer an explanation for a low-energy tail observed at high excitation energy in the translational kinetic energy map of MHP.

A comprehensive benchmark ab initio survey of the stationary points and products of the OH· + CH3OH system

Reactions between methanol and the hydroxyl radical are of significant interest for combustion-, atmospheric-, and astrochemistry. While the two primary product channels (the formation of H₂O with either CH₃O· or ·CH₂OH) have been the subject of numerous studies, the possibility of other products has seen little attention. Here, we present a comprehensive thermochemical survey of the stationary points and plausible products of the reaction, featuring 29 geometries optimized at the UCCSD(T)-F12b/aug-cc-pVTZ level, followed by accurate composite ab initio computations for all stationary points (including ·CH₂OH dissociation and isomerization) and five product channels, with a detailed evaluation of basis set convergence and efficiency. The computations reveal that the formation of methanediol and the hydroxymethoxy radical is thermodynamically favorable and the endothermicity of formaldehyde formation is low enough to be a plausible product channel. We also observe unexpectedly large energy deviations between the partially-spin-adapted ROHF-RCCSD(T) method and ROHF-UCCSD(T) as well as between UHF-UCCSDT(Q) and ROHF-UCCSDT(Q) results.

Spectral Signatures of Hydrogen Thioperoxide (HOSH) and Hydrogen Persulfide (HSSH): Possible Molecular Sulfur Sinks in the Dense ISM

For decades, sulfur has remained underdetected in molecular form within the dense interstellar medium (ISM), and somewhere a molecular sulfur sink exists where it may be hiding. With the discovery of hydrogen peroxide (HOOH) in the ISM in 2011, a natural starting point may be found in sulfur-bearing analogs that are chemically similar to HOOH: hydrogen thioperoxide (HOSH) and hydrogen persulfide (HSSH). The present theoretical study couples the accuracy in the anharmonic fundamental vibrational frequencies from the explicitly correlated coupled cluster theory with the accurate rotational constants provided by canonical high-level coupled cluster theory to produce rovibrational spectra for use in the potential observation of HOSH and HSSH. The ν6 mode for HSSH at 886.1 cm−1 is within 0.2 cm−1 of the gas-phase experiment, and the B0 rotational constant for HSSH of 6979.5 MHz is within 9.0 MHz of the experimental benchmarks, implying that the unknown spectral features (such as the first overtones and combination bands) provided herein are similarly accurate. Notably, a previous experimentally-attributed 2ν1 mode, at 7041.8 cm−1, has been reassigned to the ν1+ν5 combination band based on the present work’s ν1+ν5 value at 7034.3 cm−1. The most intense vibrational transitions for each molecule are the torsions, with HOSH having a more intense transition of 72 km/mol compared to HSSH’s intensity of 14 km/mol. Furthermore, HOSH has a larger net dipole moment of 1.60 D compared to HSSH’s 1.15 D. While HOSH may be the more likely candidate of the two for possible astronomical observation via vibrational spectroscopy due to the notable difference in their intensities, both HSSH and HOSH have large enough net dipole moments to be detectable by rotational spectroscopy to discover the role these molecules may have as possible molecular sulfur sinks in the dense ISM.

Photodissociation Dynamics of Methyl Hydroperoxide at 193 nm: A Trajectory Surface-Hopping Study

The photodissociation of methyl hydroperoxide (CH₃OOH) at 193 nm has been studied using a direct dynamics trajectory surface-hopping (TSH) method. The potential energies, energy gradients, and nonadiabatic couplings are calculated on the fly at the MRCIS(6,7)/aug-cc-pVDZ level of theory. The hopping of a trajectory from one electronic state to another is decided on the basis of Tully's fewest switches algorithm. An analysis of the trajectories reveals that the cleavage of the weakest O-O bond leads to major products CH₃O(²E) + OH(²Π), contributing about 72.7% of the overall product formation. This OH elimination was completed in the ground degenerate product state where both the ground singlet (S₀) and first excited singlet (S₁) states become degenerate. The O-H bond dissociation of CH₃OOH is a minor channel contributing about 27.3% to product formation, resulting in products CH₃OO + H. An inspection of the trajectories indicates that unlike the major channel OH elimination, the H-atom elimination channel makes a significant contribution (∼3% of the overall product formation) through the nonadiabatic pathway via conical intersection S₁/S₀ leading to ground-state products CH₃OO(X ²A″) + H(²S) in addition to adiabatic dissociation in the first excited singlet state, S₁, correlating to products CH₃OO(1 ²A') + H(²S). The computed translational energy of the majority of the OH products is found to be high, distributed in the range of 70 to 100 kcal/mol, indicating that the dissociation takes place on a strong repulsive potential energy surface. This finding is consistent with the nature of the experimentally derived translational energy distribution of OH with an average translational energy of 67 kcal/mol after the excitation of CH₃OOH at 193 nm.

On the Theoretical Determination of Photolysis Properties for Atmospheric Volatile Organic Compounds

Volatile organic compounds (VOCs) are ubiquitous atmospheric molecules that generate a complex network of chemical reactions in the troposphere, often triggered by the absorption of sunlight. Understanding the VOC composition of the atmosphere relies on our ability to characterize all of their possible reaction pathways. When considering reactions of (transient) VOCs with sunlight, the availability of photolysis rate constants, utilized in general atmospheric models, is often out of experimental reach due to the unstable nature of these molecules. Here, we show how recent advances in computational photochemistry allow us to calculate in silico the different ingredients of a photolysis rate constant, namely, the photoabsorption cross-section and wavelength-dependent quantum yields. The rich photochemistry of tert-butyl hydroperoxide, for which experimental data are available, is employed to test our protocol and highlight the strengths and weaknesses of different levels of electronic structure and nonadiabatic molecular dynamics to study the photochemistry of (transient) VOCs.

Impacts of cloud water droplets on the OH production rate from peroxide photolysis

Understanding the difference between observed and modeled concentrations of HO_x radicals in the troposphere is a current major issue in atmospheric chemistry. It is widely believed that existing atmospheric models miss a source of such radicals and several potential new sources have been proposed. In recent years, interest has increased on the role played by cloud droplets and organic aerosols. Computer modeling of ozone photolysis, for instance, has shown that atmospheric aqueous interfaces accelerate the associated OH production rate by as much as 3-4 orders of magnitude. Since methylhydroperoxide is a main source and sink of HO_x radicals, especially at low NO_x concentrations, it is fundamental to assess what is the influence of clouds on its chemistry and photochemistry. In this study, computer simulations for the photolysis of methylhydroperoxide at the air-water interface have been carried out showing that the OH production rate is severely enhanced, reaching a comparable level to ozone photolysis.

Gas phase structures of peroxides: experiments and computational problems

Gas-phase structures of several organic and inorganic peroxides X-O-O-X and X-O-O-X', which have been determined experimentally by gas electron diffraction and/or microwave spectroscopy, are discussed. The OO bond length in these peroxides varies from 1.481(8) Å in Me3 SiOOSiMe3 to 1.214(2) Å in FOOF and the dihedral angle ϕ(XO-OX) between 0° in HC(O)O-OH and near 180° in Bu(t) O-OBu(t) . Some of the peroxides cause problems for quantum chemistry, since several computational methods fail to reproduce the experimental structures. Extreme examples are MeO-OMe and FO-OF. In the case of MeO-OMe only about half of the more than 100 computational methods reported in the literature reproduce the experimentally determined double-minimum shape of the torsional potential around the OO bond correctly. For FO-OF only a small number of close to 200 computational methods reproduce the OO and OF bond lengths better than ±0.02 Å.

UV photochemistry of peroxyformic acid (HC(O)OOH): an experimental and computational study investigating 355 nm photolysis

The photochemistry of peroxyformic acid (PFA), a molecule of atmospheric interest exhibiting internal hydrogen bonding, is examined by exciting the molecule at 355 nm and detecting the nascent OH fragments using laser-induced fluorescence. The OH radicals are found to be formed in their ground electronic state with the vast majority of available energy appearing in fragment translation. The OH fragments are vibrationally cold (v" = 0) with only modest rotational excitation. The average rotational energy is determined to be 0.35 kcal/mol. Further, the degree of OH rotational excitation from PFA is found to be significantly less than that arising from the dissociation of H2O2 as well as other hydroperoxides over the same wavelength. Ab initio calculation at the EOM-CCSD level is used to investigate the first few electronic excited states of PFA. Differences in the computed torsional potential between PFA and H2O2 help rationalize the observed variation in their respective OH fragment rotational excitation. The calculations also establish that the electronic excited state of PFA accessed in the near UV is of (1)A" symmetry and involves a σ*(O-O) ← n(O) excitation. Additionally, the UV absorption cross section of PFA at 355 and 282 nm is estimated by comparing the yield of OH from PFA at these wavelengths to that from hydrogen peroxide for which the absorption cross sections is known.

The HOOH UV spectrum: importance of the transition dipole moment and torsional motion from semiclassical calculations on an ab initio potential energy surface

The absorption cross section of HOOH, a starting point for larger ROOH, was calculated using the "Wigner method." Calculations use the Wigner transform of ground state wave functions and classical approximations for excited state wave functions. Potential energy and transition dipole moment surfaces were calculated using the equation-of-motion coupled-cluster singles and doubles method over an extended Franck-Condon region. The first two O-O stretches and the first five HOOH torsional levels are included. This study also addresses two fundamental questions about ROOH photodissociation. The long wavelength A(1)A:B(1)B excited state preference has been measured from dynamics experiments, but a Franck-Condon overlap explanation has not been directly verified. A moderate barrier to HOOH torsional motion and excited state dynamics affect the temperature dependence in the UV spectrum. Based on these initial findings for HOOH, photodissociation of large ROOH cannot be eliminated as an important factor for ozone and particulate matter production seen in both ambient and laboratory studies.

Ultrafast photochemistry of methyl hydroperoxide on ice particles

Simulations show that photodissociation of methyl hydroperoxide, CH(3)OOH, on water clusters produces a surprisingly wide range of products on a subpicosecond time scale, pointing to the possibility of complex photodegradation pathways for organic peroxides on aerosols and water droplets. Dynamics are computed at several excitation energies at 50 K using a semiempirical PM3 potential surface. CH(3)OOH is found to prefer the exterior of the cluster, with the CH(3)O group sticking out and the OH group immersed within the cluster. At atmospherically relevant photodissociation wavelengths the OH and CH(3)O photofragments remain at the surface of the cluster or embedded within it. However, none of the 25 completed trajectories carried out at the atmospherically relevant photodissociation energies led to recombination of OH and CH(3)O to form CH(3)OOH. Within the limited statistics of the available trajectories the predicted yield for the recombination is zero. Instead, various reactions involving the initial fragments and water promptly form a wide range of stable molecular products such as CH(2)O, H(2)O, H(2), CO, CH(3)OH, and H(2)O(2).