An Ab Initio Neural Network Potential Energy Surface for the Dimer of Formic Acid and Further Quantum Tunneling Dynamics

We construct a full-dimensional ab initio neural network potential energy surface (PES) for the isomerization system of the formic acid dimer (FAD). This is based upon ab initio calculations using the DLPNO-CCSD(T) approach with the aug-cc-pVTZ basis set, performed at over 14000 symmetry-unique geometries. An accurate fit to the obtained energies is generated using a general neural network fitting procedure combined with the fundamental invariant method, and the overall energy-weighted root-mean-square fitting error is about 6.4 cm^-1. Using this PES, we present a multidimensional quantum dynamics study on tunneling splittings with an efficient theoretical scheme developed by our group. The ground-state tunneling splitting of FAD calculated with a four-mode coupled method is in good agreement with the most recent experimental measurements. The PES can be applied for further dynamics studies. The effectiveness of the present scheme for constructing a high-dimensional PES is demonstrated, and this scheme is expected to be feasible for larger molecular systems.

Efficient Quantum Mechanical Calculations of Mode-Specific Tunneling Splittings upon Fundamental Excitation in the Dimer of Formic Acid

The formic acid dimer (FAD) is an important benchmark system for understanding the double hydrogen transfer process. Most recently, Zhang et al. measured a few tunneling splittings upon fundamental excitation of FAD precisely (Zhang, Y. et al. J. Chem. Phys. 2017, 146, 244306); however, relevant theoretical studies are very limited. Here, we present a multidimensional quantum dynamics study on mode-specific tunneling splittings upon fundamental excitation in FAD with an efficient theoretical scheme developed by our group in which the process-oriented basis function customization strategy is combined with the preconditioned inexact spectral transform method. Various mode-specific tunneling splittings upon fundamental excitation are systematically calculated, and interesting mode-specific excitation effects on tunneling rate are identified. In particular, the calculated tunneling splittings for the ν₂₂ and ν₂₁ states are in good agreement with experiment, and the remarkable mode-specific suppression effects upon excitation should result from that the antisymmetric vibrational modes hinder the concerted double H-transfer. The present work is helpful to acquire a better understanding of the mode-specific excitation effects on tunneling processes.

Hydrogen-Atom Tunneling in Metaphosphorous Acid

Metaphosphorous acid (HOPO), a key intermediate in phosphorus chemistry, has been generated in syn- and anti-conformations in the gas phase by high-vacuum flash pyrolysis (HVFP) of a molecular precursor ethoxyphosphinidene oxide (EtOPO→C₂ H₄ +HOPO) at ca. 1000 K and subsequently trapped in an N₂ -matrix at 2.8 K. Unlike the two conformers of the nitrogen analogue HONO, the anti-conformer of HOPO undergoes spontaneous rotamerization at 2.8 K via hydrogen-atom tunneling (HAT) with noticeable kinetic isotope effects for H/D (>10⁴ for DOPO) and ¹⁶ O/¹⁸ O (1.19 for H¹⁸ OPO and 1.06 for HOP¹⁸ O) in N₂ -matrices.

Riddles of the structure and vibrational dynamics of HO3 resolved near the ab initio limit

The hydridotrioxygen (HO₃) radical has been investigated in many previous theoretical and experimental studies over several decades, originally because of its possible relevance to the tropospheric HO_x cycle but more recently because of its fascinating chemical bonding, geometric structure, and vibrational dynamics. We have executed new, comprehensive research on this vexing molecule via focal point analyses (FPA) to approach the ab initio limit of optimized geometric structures, relative energies, complete quartic force fields, and the entire reaction path for cis-trans isomerization. High-order coupled cluster theory was applied through the CCSDT(Q) and even CCSDTQ(P) levels, and CBS extrapolations were performed using cc-pVXZ (X = 2-6) basis sets. The cis isomer proves to be higher than trans by 0.52 kcal mol^-1, but this energetic ordering is achieved only after the CCSDT(Q) milestone is reached; the barrier for cis → trans isomerization is a minute 0.27 kcal mol^-1. The FPA central r_e(O-O) bond length of trans-HO₃ is astonishingly long (1.670 Å), consistent with the semiexperimental r_e distance we extracted from microwave rotational constants of 10 isotopologues using FPA vibration-rotation interaction constants (α_i). The D₀(HO-O₂) dissociation energy converges to a mere 2.80 ± 0.25 kcal mol^-1. Contrary to expectation for such a weakly bound system, vibrational perturbation theory performs remarkably well with the FPA anharmonic force fields, even for the torsional fundamental near 130 cm^-1. Exact numerical procedures are applied to the potential energy function for the torsional reaction path to obtain energy levels, tunneling rates, and radiative lifetimes. The cis → trans isomerization occurs via tunneling with an inherent half-life of 1.4 × 10^-11 s and 8.6 × 10^-10 s for HO₃ and DO₃, respectively, thus resolving the mystery of why the cis species has not been observed in previous experiments executed in dissipative environments that allow collisional cooling of the trans-HO₃ product. In contrast, the pure ground eigenstate of the cis species in a vacuum is predicted to have a spontaneous radiative lifetime of about 1 h and 5 days for HO₃ and DO₃, respectively.

Anab initiobased full-dimensional potential energy surface for OH + O2⇄ HO3and low-lying vibrational levels of HO3

A full-dimensional potential energy surface for HO₃, including the HO + O₂dissociation asymptote, is developed and rigorous quantum dynamics calculations based on this PES have been carried out to compute the vibrational energy levels of HO₃.

An Evaluation of Gas Phase Enthalpies of Formation for Hydrogen-Oxygen (HxOy) Species

We have compiled gas phase enthalpies of formation for nine hydrogen-oxygen species (HxOy) and selected recommended values for H, O, OH, H2O, HO2, H2O2, O3, HO3, and H2O3. The compilation consists of values derived from experimental measurements, quantum chemical calculations, and prior evaluations. This work updates the recommended values in the NIST-JANAF (1985) and Gurvich et al. (1989) thermochemical tables for seven species. For two species, HO3 and H2O3 (important in atmospheric chemistry) and not found in prior thermochemical evaluations, we also provide supplementary data consisting of molecular geometries, vibrational frequencies, and torsional potentials which can be used to compute thermochemical functions. For all species, we also provide supplementary data consisting of zero point energies, vibrational frequencies, and ion reaction energetics.

Weakly Bound Clusters in Astrochemistry? Millimeter and Submillimeter Spectroscopy of trans-HO3 and Comparison to Astronomical Observations

The emergence of chemical complexity during star and planet formation is largely guided by the chemistry of unstable molecules that are reaction intermediates in terrestrial chemistry. Our knowledge of these intermediates is limited by both the lack of laboratory studies and the difficulty in their astronomical detection. In this work, we focus on the weakly bound cluster HO3 as an example of the connection between laboratory spectroscopic study and astronomical observations. Here, we present a fast-sweep spectroscopic technique in the millimeter and submillimeter range to facilitate the laboratory search for trans-HO3 and DO3 transitions in a discharge supersonic jet and report their rotational spectra from 70 to 450 GHz. These new measurements enable full determination of the molecular constants of HO3 and DO3. We also present a preliminary search for trans-HO3 in 32 star-forming regions using this new spectroscopic information. HO3 is not detected, and column density upper limits are reported. This work provides additional benchmark information for computational studies of this intriguing radical, as well as a reliable set of molecular constants for extrapolation of the transition frequencies of HO3 for future astronomical observations.

Is HO3− multiple-minimum and floppy? Covalent to van der Waals isomerization and bond rupture of a peculiar anion

The HO₃⁻ anion is multiple-minimum and floppy: the two main isomers and isomerization barrier all lie quite below dissociation.