Heats of Formation of the H1,2OmSn (m, n = 0−3) Molecules from Electronic Structure Calculations

Correcting the Experimental Enthalpies of Formation of Some Members of the Biologically Significant Sulfenic Acids Family

Sulfenic acids are important intermediates in the oxidation of cysteine thiol groups in proteins by reactive oxygen species. The mechanism is influenced heavily by the presence of polar groups, other thiol groups, and solvent, all of which determines the need to compute precisely the energies involved in the process. Surprisingly, very scarce experimental information exists about a very basic property of sulfenic acids, the enthalpies of formation. In this Article, we use high level quantum chemical methods to derive the enthalpy of formation at 298.15 K of methane-, ethene-, ethyne-, and benzenesulfenic acids, the only ones for which some experimental information exists. The methods employed were tested against well-known experimental data of related species and extensive CCSD(T) calculations. Our best results consistently point out to a much lower enthalpy of formation of methanesulfenic acid, CH₃SOH (Δ_fH⁰(298.15K) = -35.1 ± 0.4 kcal mol^-1), than the one reported in the NIST thermochemical data tables. The enthalpies of formation derived for ethynesulfenic acid, HC≡CSOH, +32.9 ± 1.0 kcal/mol, and benzenesulfenic acid, C₆H₅SOH, -2.6 ± 0.6 kcal mol^-1, also differ markedly from the experimental values, while the enthalpy of formation of ethenesulfenic acid CH₂CHSOH, not available experimentally, was calculated as -11.2 ± 0.7 kcal mol^-1.

Anharmonic fundamental vibrational frequencies and spectroscopic constants of the potential HSO2 radical astromolecule

The recent report that HSO₂ is likely kinetically favored over the HOSO thermodynamic product in hydrogen addition to sulfur dioxide in simulated Venusian atmospheric conditions has led to the need for reference rotational, vibrational, and rovibrational spectral data for this molecule. While matrix-isolation spectroscopy has been able to produce vibrational frequencies for some of the vibrational modes, the full infrared to microwave spectrum of 1 ²A' HSO₂ is yet to be generated. High-level quantum chemical computations show in this work that the >2.5 D dipole moment of this radical makes it a notable target for possible radioastronomical observation. Additionally, the high intensity antisymmetric S-O stretch is computed here to be 1298.3 cm^-1, a 13.9 cm^-1 blueshift up from H₂ matrix analysis. In any case, the full set of rotational and spectroscopic constants and anharmonic fundamental vibrational frequencies is provided in this work in order to help characterize HSO₂ and probe its kinetic favorability.

An accurate full-dimensional potential energy surface for the reaction OH + SO → H + SO2

An accurate full-dimensional PES for the OH + SO ↔ H + SO₂ reaction is developed by the permutation invariant polynomial-neural network approach.

A New Mechanism of Acid Rain Generation from HOSO at the Air-Water Interface

The photochemistry of SO₂ at the air-water interface of water droplets leads to the formation of HOSO radicals. Using first-principles simulations, we show that HOSO displays an unforeseen strong acidity (pK_a = -1) comparable with that of nitric acid and is fully dissociated at the air-water interface. Accordingly, this radical might play an important role in acid rain formation. Potential implications are discussed.

Generalized Valence Bond Description of Chalcogen-Nitrogen Compounds. III. Why the NO-OH and NS-OH Bonds Are So Different

Crabtree et al. recently reported the microwave spectrum of nitrosyl-O-hydroxide (trans-NOOH), an isomer of nitrous acid, and found that this molecule has the longest O-O bond ever observed: 1.9149 Å ± 0.0005 Å. This is in marked contrast to the structure of the valence isoelectronic trans-NSOH molecule, which has a normal NS-OH bond length and strength. Generalized valence bond calculations show that the long bond in trans-NOOH is the result of a weak through-pair interaction that singlet couples the spins of the electrons in singly occupied orbitals on the hydroxyl radical and nitrogen atom, an interaction that is enhanced by the intervening lone pair of the oxygen atom in NO. The NS-OH bond is the result of the formation of a stable recoupled pair bond dyad, which accounts for both its length and strength.

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.

Mechanisms of reactions of sulfur hydride hydroxide: tautomerism, condensations, and C-sulfenylation and O-sulfenylation of 2,4-pentanedione

The conformations, equilibrium structures, hydrogen bonds, and non-covalent interactions involved in the mechanisms of tautomerization, condensations, and C-sulfenylation and O-sulfenylation of 2,4-pentanedione by sulfur hydride hydroxide (hydrogen thioperoxide, oxadisulfane, H-SOH) have been studied using BD(T), CCSD(T), and QCISD(T) with the cc-pVTZ basis set and using B3LYP, B3PW91, CAM-B3LYP, PBE1PBE, PBEh1PBE, LC-ωPBE, M06-2X, and ωB97XD with the 6-311+G(d,p) basis set. All levels of theory predict the sulfenyl (H-SOH) tautomer of hydrogen thioperoxide to be lower in energy than the sulfinyl (H2S═O) tautomer. Four reasonable mechanisms were considered for the tautomerization of the sulfenyl tautomer of hydrogen thioperoxide to the sulfinyl tautomer: a cyclic three-membered water-free transition state (TS, CCSD(T) activation energy barrier E(⧧) = 65.1 kcal/mol), a cyclic five-membered transition state with one water molecule (TSH2O, E(⧧) = 31.1 kcal/mol), a cyclic seven-membered transition state with two water molecules (TS2H2O, E(⧧) = 14.5 kcal/mol), and a cyclic nine-membered transition state with three water molecules (TS3H2O, E(⧧) = 5.6 kcal/mol). The mechanisms involve hydrogen-bonded reactant complexes and hydrogen-bonded product complexes. The CCSD(T)-predicted energy barriers for the condensation of hydrogen thioperoxide to form thiosulfinic acid through transition states with zero, one, and two waters are E(⧧) = 42.0, 18.3, and 0 kcal/mol, respectively. Mixed condensation reactions are predicted to afford organosulfur products and compounds containing sulfur-selenium bonds. Hydrogen thioperoxide is predicted to add to 2,4-pentanedione to form C-sulfenylated (sulfide, thioether) and O-sulfenylated (sulfenate ester) products. Similar mechanistic trends and reaction pathways are observed in the tautomerism, condensations, and C-sulfenylation and O-sulfenylation reactions of hydrogen thioperoxide. The water molecules set up proton relay networks (bridges) that reduce ring strain, generate favorable conformations for reactivity, lower energy barriers, and increase the numbers of stabilizing hydrogen bonds and nonbonding interactions.

Is near-"spectroscopic accuracy" possible for heavy atoms and coupled cluster theory? An investigation of the first ionization potentials of the atoms Ga-Kr

Recent developments in ab initio coupled cluster (CC) theory and correlation consistent basis sets have ushered in an era of unprecedented accuracy when studying the spectroscopy and thermodynamics of molecules containing main group elements. These same developments have recently seen application to heavier inorganic or transition metal-containing species. The present work benchmarks conventional single reference coupled cluster theory (up to full configuration interaction for valence electron correlation and coupled cluster with up to full pentuple excitations (CCSDTQP) for core-valence correlation) and explicitly correlated coupled cluster methods [CC with single, double, and perturbative triple substitutions (CCSD(T)-F12)] for the atomic ionization potentials of the six 4p elements (Ga-Kr), a property with experimental error bars no greater than a few cm(-1). When second-order spin orbit coupling effects are included, a composite methodology based on CCSD(T) calculations yielded a mean signed error of just -0.039 kcal mol(-1) and a mean unsigned error of 0.043 kcal mol(-1). Inclusion of post-CCSD(T) correlation corrections reduced both of these values to -0.008 kcal mol(-1) and 0.025 kcal mol(-1), respectively, with the latter corresponding to an average error of just 9 cm(-1). The maximum signed error in the latter scheme was just -0.043 kcal mol(-1) (15 cm(-1)).

Dehydrative cyclocondensation mechanisms of hydrogen thioperoxide and of alkanesulfenic acids

Structural features of hydrogen thioperoxide (oxadisulfane, H-S-O-H) and of alkanesulfenic acids (R-S-O-H; R = CH3, CH2CH3, CH2CH2CH3, CH(CH3)2, C(CH3)3, CF3, CCl3) and the mechanisms of their dehydrative cyclocondensation to the respective sulfinothioic acid (H-(S═O)-S-H) and alkyl alkanethiosulfinates (R-(S═O)-S-R) have been studied using coupled cluster theory with single and double and perturbative triple excitations [CCSD(T)] and quadratic configuration interaction with single and double and perturbative triple excitations [QCISD(T)] with the cc-pVDZ basis set and also using second-order Møller-Plesset perturbation theory (MP2) and the hybrid density functionals B3LYP, B3PW91, and PBE1PBE with the 6-311+G(d,p) basis set. The concerted cyclodehydration mechanisms include cyclic five-center transition states with relatively long distance sulfur-sulfur bonding interactions. Attractive and repulsive nonbonding interactions are predicted in the sulfenic acids, transition states, and thiosulfinates. In the alkyl alkanethiosulfinates attractive cyclic C-H----O═S nonbonding interactions are predicted. CCSD(T) and QCISD(T) predict similar values for the relative energies and CCSD(T) predicts the barrier to the cyclocondensation of H-S-O-H to sulfinothioic acid (H-(S═O)-S-H) to be 41.8 kcal/mol, and barriers in the range of 37.5 to 39.6 kcal/mol are predicted for the condensation of alkanesulfenic acids to alkyl alkanethiosulfinates.

Theoretical investigation of carbon-sulfur triple bonds

By means of first principle calculations we have investigated a set of molecules that are presumed to contain carbon-sulfur triple bonds, namely HCSOH, H(3)SCH, cis-FCSF, F(3)CCSF(3), and F(5)SCSF(3). For HCSOH, FCSF, and H(3)SCH we used the CCSD(T) methodology and the correlation-consistent basis sets. On the other hand, F(3)CCSF(3) and F(5)SCSF(3) were studied at the B3LYP, M06-2X, MP2, and G3 levels of theory. We found that none of these molecules display a carbon-sulfur adiabatic bond dissociation energy (ABDE) as strong as diatomic CS (170.5 kcal mol(-1)), or a diabatic bond dissociation energy (DBDE) larger than the one found in SCO (212.0 kcal mol(-1)), although the DBDE of FCSF comes quite close at 208.3 kcal mol(-1). The CS ABDEs of F(3)CCSF(3), F(5)SCSF(3), and H(3)SCH are comparable to that of a single C-S bond. In contrast with the experimental results, F(3)CCSF(3) and F(5)SCSF(3) are predicted to be linear with C(3v) and C(s) symmetry, respectively, at the B3LYP/6-311+G(3df,2p) level. MP2/6-311+G(2df,2p) calculations support the C(3v) symmetry for F(3)CCSF(3), despite F(5)SCSF(3) not having a perfect linear structure; the CSC angle is 174.6°, which is nearly 20° larger than the experimental value. The analysis of the carbene structures of HCSOH and H(3)SCH revealed that they are not significant, because the triplet state is dissociative in these cases. However, for F(3)CCSF(3) and F(5)SCSF(3) , the carbene triplet states lie 0.81 and 0.77 eV above the singlet state, respectively. In the same vein, our investigation supports the presence of a strong double bond for HCSOH. The conflicting evidence available for F(3)CCSF(3) and F(5)SCSF(3) makes it very difficult to determine the nature of the CS bonds. However, the bond dissociation energies and the singlet-triplet splittings clearly suggest that these compounds should be considered as masked sulfinylcarbenes. The analysis of the bond dissociation energies challenges the existence of a triple bond in these five molecules, but from a strictly thermodynamic standpoint, cis-FCSF is found to be the candidate most likely to exhibit triple-bond character.

Reversal of the relative stability of the isomeric radicals HSO and HOS upon hydration and their reactions with ozone

The radical HSO is an oxidation product of pollutants such as H(2)S and CH(3)SH in Earth's atmosphere. For the first time, the interaction of HSO and its tautomer HOS with single water molecules to yield the hydrates HSO.nH(2)O and HOS.nH(2)O was studied for n = 1-3, applying the high-level G3X(MP2) theory. A large number of structures corresponding to local minima on the potential energy surfaces has been identified. While gaseous HSO is more stable than HOS, the enthalpy diffference between HSO.nH(2)O and HOS.nH(2)O decreases with increasing degree of hydration and becomes practically zero for n = 3. Thus, in aqueous solution as well as in fog and rain droplets, HOS is expected to compete with HSO. The barrier for the tautomerization of HSO to HOS is dramatically lowered by the presence of water molecules since a cyclic transition state allows a concerted proton shift within the system of neighboring hydrogen bonds. The corresponding activation enthalpy of only 73.5 kJ mol(-1) predicted for the transformation of HSO.2H(2)O into HOS.2H(2)O may be compared to the 202 kJ mol(-1) reported for the tautomerization of the unhydrated gaseous HSO/HOS molecules. The impact of water of hydration on the fundamental vibrational modes of HSO and HOS has also been studied. Furthermore, HOS is predicted to dimerize at low temperatures to give two van der Waals molecules with singlet (symmetry C(2)) or triplet configuration (symmetry C(2h)), the latter being more stable than the singlet isomer. The disproportionation of 2HSO to H(2)S and SO(2) is predicted to be exothermic by -263.5 kJ mol(-1). The reaction of HSO with ozone to HSO(2) and O(2) is also strongly exothermic by -274.0 kJ mol(-1) and seems to proceed without any barrier. HOS forms a 1:1 van der Waals complex with O(3); the redox reaction of its two components is calculated as exothermic by -410.9 kJ mol(-1) and results in a rather stable adduct between HOSO and O(2) with the structure of a peroxo isomer of HOSO(3). This unprecedented hydrogen peroxosulfite radical might open a novel route to atmospheric sulfate without the intermediate formation of SO(2) and SO(3).

Method and basis set dependence of anharmonic ground state nuclear wave functions and zero-point energies: application to SSSH

One of the largest remaining errors in thermochemical calculations is the determination of the zero-point energy (ZPE). The fully coupled, anharmonic ZPE and ground state nuclear wave function of the SSSH radical are calculated using quantum diffusion Monte Carlo on interpolated potential energy surfaces (PESs) constructed using a variety of method and basis set combinations. The ZPE of SSSH, which is approximately 29 kJ mol(-1) at the CCSD(T)/6-31G* level of theory, has a 4 kJ mol(-1) dependence on the treatment of electron correlation. The anharmonic ZPEs are consistently 0.3 kJ mol(-1) lower in energy than the harmonic ZPEs calculated at the Hartree-Fock and MP2 levels of theory, and 0.7 kJ mol(-1) lower in energy at the CCSD(T)/6-31G* level of theory. Ideally, for sub-kJ mol(-1) thermochemical accuracy, ZPEs should be calculated using correlated methods with as big a basis set as practicable. The ground state nuclear wave function of SSSH also has significant method and basis set dependence. The analysis of the nuclear wave function indicates that SSSH is localized to a single symmetry equivalent global minimum, despite having sufficient ZPE to be delocalized over both minima. As part of this work, modifications to the interpolated PES construction scheme of Collins and co-workers are presented.