In search of the next Holy Grail of polyoxide chemistry: Explicitly correlated ab initio full quartic force fields for HOOH, HOOOH, HOOOOH, and their isotopologues

Spectroscopy and Photochemistry of OAlNO and Implications for New Metal Chemistry in the Atmosphere

A new aluminum-bearing species, OAlNO, which has the potential to impact the chemistry of the Earth's upper atmosphere, is characterized via high-level, ab initio, spectroscopic methods. Meteor-ablated aluminum atoms are quickly oxidized to aluminum oxide (AlO) in the mesosphere and lower thermosphere (MLT), where a steady-state layer of AlO then builds up. Concurrent formation of nitric oxide (NO) in the same region of the atmosphere will lead to the bimolecular formation of the OAlNO molecule. Molecular orbital analysis provides fundamental insights into the chemical bonding and energetic arrangement of the triplet (1 3A″) ground state and singlet (1 1A') excited-state species of OAlNO. Additionally, unpaired electrons on the terminal oxygen atom of triplet (1 3A″) OAlNO cause it to be reactive to atmospheric species, potentially impacting climate science and high-altitude chemistry. The triplet (1 3A″) ground-state species exhibits a large permanent dipole moment useful for rotational spectroscopic detection; however, similar rotational constants to the singlet (1 1A') excited-state species will hamper differentiation in a spectrum. Strong infrared intensities will assist in detection and discrimination of the different spin states and isomers. Repulsive electronic excited states of OAlNO will lead to photolysis of the Al-N bond and formation of various electronic states of AlO + NO through nonadiabatic pathways. Reaction through the OAlNO intermediate represents a means for the production of electronically excited AlO, leading to new chemistry in the atmosphere. Excitation to higher-lying electronic states will lead to fluorescence with a minor Stokes shift, useful for laboratory investigation. Such physical properties of this molecule will allow for new, unexplored chemical pathways in the MLT to be considered.

Permutation-Invariant-Polynomial Neural-Network-Based Δ-Machine Learning Approach: A Case for the HO2 Self-Reaction and Its Dynamics Study

Δ-machine learning, or the hierarchical construction scheme, is a highly cost-effective method, as only a small number of high-level ab initio energies are required to improve a potential energy surface (PES) fit to a large number of low-level points. However, there is no efficient and systematic way to select as few points as possible from the low-level data set. We here propose a permutation-invariant-polynomial neural-network (PIP-NN)-based Δ-machine learning approach to construct full-dimensional accurate PESs of complicated reactions efficiently. Particularly, the high flexibility of the NN is exploited to efficiently sample points from the low-level data set. This approach is applied to the challenging case of a HO₂ self-reaction with a large configuration space. Only 14% of the DFT data set is used to successfully bring a newly fitted DFT PES to the UCCSD(T)-F12a/AVTZ quality. Then, the quasiclassical trajectory (QCT) calculations are performed to study its dynamics, particularly the mode specificity.

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.

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.

The atmospheric degradation pathways of BrCH2O2: computational calculation on mechanisms of the reaction with HO2

Mechanisms for the atmospheric degradation reaction of BrCH2O2+HO2 were investigated using quantum chemistry methods. The result indicates that the dominant product is BrCH2OOH+O2((3)Σ). While CH2O+HBr+O3, BrCHO+OH+HO2 and CH2O+Br+HO3 will be competitive to a certain extent in the atmosphere. Meanwhile, the nascent product - BrCH2OOH reacts easily with OH radicals leading to BrCH2O2 again under the atmospheric conditions. Moreover, OH radicals could act as a catalyst in the net reaction of BrCH2OOH→BrCHO+H2O. Thus the proposed product BrCHO+H2O+O2 in the experiment might be generated from the subsequent reaction of BrCH2OOH with extra OH radicals. Comparisons indicate that halogen substitution effect makes minor contributions to the XCH2O2 (X=H, F, Cl and Br)+HO2 reactions in the atmosphere.

Analytic gradients for density cumulant functional theory: the DCFT-06 model

Density cumulant functional theory (DCFT) is one of a number of nascent electron correlation methods that are derived from reduced density matrices and cumulants thereof, instead of the wavefunction. Deriving properties from the density cumulant naturally yields methods that are size extensive and size consistent. In this work, we derive expressions for the analytic gradient, with respect to an external perturbation, for the DCFT-06 variant of density cumulant functional theory. Despite the fact that the DCFT-06 energy functional is stationary with respect to the density cumulant, the analytic gradients of the energy require the solution of perturbation-independent equations for both orbital and cumulant response. These two sets of linear response equations are coupled in nature and are solved iteratively with the solution of orbital and cumulant response equations each macroiteration, exhibiting rapid convergence. The gradients are implemented and benchmarked against coupled cluster theory with single and double excitations (CCSD) and CCSD with perturbative triple excitations [CCSD(T)], as well as accurate empirically corrected experimental data, for a test set comprising 15 small molecules. For most of the test cases, results from DCFT-06 are closer to CCSD(T) and empirical data than those from CCSD. Although the total energy and analytic gradient have the same asymptotic scaling, the present experience shows that the computational cost of the gradient is significantly lower.

Arbitrary order El'yashevich-Wilson B tensor formulas for the most frequently used internal coordinates in molecular vibrational analyses

In recent years, internal coordinates have become the preferred means of expressing potential energy surfaces. The ability to transform quantities from chemically significant internal coordinates to primitive Cartesian coordinates and spectroscopically relevant normal coordinates is thus critical to the further development of computational chemistry. In the present work, general nth order formulas are presented for the Cartesian derivatives of the five most commonly used internal coordinates--bond stretching, bond angle, torsion, out-of-plane angle, and linear bending. To compose such formulas in a reasonably understandable fashion, a new notation is developed that is a generalization of that which has been used previously for similar purposes. The notation developed leads to easily programmable and reasonably understandable arbitrary order formulas, yet it is powerful enough to express the arbitrary order B tensor of a general, N-point internal coordinate, as is done herein. The techniques employed in the derivation of such formulas are relatively straightforward, and could presumably be applied to a number of other internal coordinates as needed.