Reactivity properties of the HOSO and HSO2 isomers in liquid medium: a sequential Monte Carlo/quantum mechanics study

CONTEXT

The rationalization of acid rain formation steps is fundamental for mitigating its effects. It is believed the hydroxysulfinyl radical is an intermediate species for the production of atmospheric sulfuric acid. Two stable configurations HOSO and HSO₂ have been reported for such a radical in the gas phase. This work aims at studying these isomers in the aqueous medium. The electrical and reactivity quantities - electronic chemical potential ([Formula: see text]), chemical hardness ([Formula: see text]), and electrophilicity ([Formula: see text]) - are here calculated and compared. Considering first solvation shells (15 H₂O for HSO₂ and 23 H₂O for HOSO), an increase of 25% in the dipole moment of HSO₂ was obtained, while the dipole moment of HOSO decreases in 11%. Both solvated isomers grow softer ([Formula: see text] decreases) contrasted to the gas phase.

METHODS

HOSO and HSO₂ are studied through a sequential Monte Carlo/quantum mechanics approach. Lennard-Jones plus the Coulomb potentials were used to represent intermolecular potential interaction in the frame of the DICE package. Molecular structure calculations were performed at CASPT2/aug - cc - pV(T + d)Z level of theory using the MOLPRO suite of programs.

Unraveling sulfur chemistry in interstellar carbon oxide ices

Formyl radical (HCO•) and hydroxycarbonyl radical (HOCO•) are versatile building blocks in the formation of biorelevant complex organic molecules (COMs) in interstellar medium. Understanding the chemical pathways for the formation of HCO• and HOCO• starting with primordial substances (e.g., CO and CO2) is of vital importance in building the complex network of prebiotic chemistry. Here, we report the efficient formation of HCO• and HOCO• in the photochemistry of hydroxidooxidosulfur radical (HOSO•)–a key intermediate in SO2 photochemistry–in interstellar analogous ices of CO and CO2 at 16 K through hydrogen atom transfer (HAT) reactions. Specifically, 266 nm laser photolysis of HOSO• embedded in solid CO ice yields the elusive hydrogen‑bonded complexes HCO•···SO2 and HOCO•···SO, and the latter undergoes subsequent HAT to furnish CO2···HOS• under the irradiation conditions. Similar photo-induced HAT of HOSO• in solid CO2 ice leads to the formation of HOCO•···SO2. The HAT reactions of HOSO• in astronomical CO and CO2 ices by forming reactive acyl radicals may contribute to understanding the interplay between the sulfur and carbon ice-grain chemistry in cold molecular clouds and also in the planetary atmospheric chemistry.

Validating experiments for the reaction H2 + NH2- by dynamical calculations on an accurate full-dimensional potential energy surface

Ion-molecule reactions play key roles in the field of ion related chemistry. As a prototypical multi-channel ion-molecule reaction, the reaction H₂ + NH₂^- → NH₃ + H^- has been studied for decades. In this work, we develop a new globally accurate potential energy surface (PES) for the title system based on hundreds of thousands of sampled points over a wide dynamically relevant region that covers long-range interacting configuration space. The permutational invariant polynomial-neural network (PIP-NN) method is used for fitting and the resulting total root mean squared error (RMSE) is extremely small, 0.026 kcal mol^-1. Extensive dynamical and kinetic calculations are carried out on this PIP-NN PES. Impressively, a unique phenomenon of significant reactivity suppression by exciting the rotational mode of H₂ is reported, supported by both the quasi-classical trajectory (QCT) and quantum dynamics (QD) calculations. Further analysis uncovers that exciting the H₂ rotational mode would prevent the formation of the reactant complex and thus suppress the reactivity. The calculated rate coefficients for H₂/D₂ + NH₂^- agree well with the experimental results, which show an inverse temperature dependence from 50 to 300 K, consistent with the capture nature of this barrierless reaction. The significant kinetic isotope effect observed by experiments is well reproduced by the QCT computations as well.

Spectroscopic Characterization of HSO2• and HOSO• Intermediates Involved in SO2 Geoengineering

Sulfur-containing radicals HSO₂^• and HOSO^• are key intermediates involved in stratospheric sulfur geoengineering by SO₂ injection. The spectroscopic characterization and photochemistry of both radicals are crucial to understanding the chemical impact of SO₂ chemistry in the stratosphere. On the basis of the efficient generation of HOSO^• by flash pyrolysis of gaseous sulfinic acid, CHF₂S(O)OH, a strong absorption is observed at 270 nm along with a shoulder up to 350 nm for HOSO^• isolated in low-temperature noble gas matrixes (Ar and Ne). These mainly arise from the excitations from the ground state (X²A) to the C²A/D²A and A²A/B²A states, respectively. Upon a 266 nm laser irradiation, the broad absorption band in the range 320-500 nm for HSO₂^• appears, and it corresponds to the combination of three excitations from the X²A state to the first (A²A), second (B²A), and third (C²A) excited states. Assignment of the UV-vis spectra is consistent with the photochemistry of HOSO^• and HSO₂^• as observed by matrix-isolation IR spectroscopy and also by the agreement with high-level ab initio calculations.

Application of Reaction Path Search Calculations to Potential Energy Surface Fits

There has been significant progress in recent years in the use of machine learning techniques to model high-dimensional reactive potential energy surfaces using large-scale data obtained from ab initio electronic structure calculations. In these methods, the strategy used to gather data becomes a key issue as the molecular size increases. In this work, we examine the applicability of the reaction path search algorithm implemented in the Global Reaction Route Mapping (GRRM) code as a data-gathering approach. The electronic energies and gradients sampled by using the GRRM calculation are directly used in potential energy surface fitting to a permutationally invariant polynomial function. This simple approach was applied to the HNS and HCNO reaction systems, and we found that the fitted potential energy surfaces reasonably reproduce the features of the electronic structure calculations used in the GRRM calculations. This suggests that the GRRM sampling scheme can be used to construct an initial potential energy surface.

HCl–H2O dimer: an accurate full-dimensional potential energy surface and fully coupled quantum calculations of intra- and intermolecular vibrational states and frequency shifts

The present work reports a new full-dimensional potential energy surface (PES) of the HCl–H₂O dimer, and the first fully coupled 9D quantum calculations of the intra- and intermolecular vibrational states of the complex, utilizing this PES.