High-Accuracy Theoretical Thermochemistry of Atmospherically Important Sulfur-Containing Molecules

Probing the Chemistry of Sulfurous Pollutants: Accurate Thermochemistry Determination of Extensive Sulfur-Containing Species

Sulfur-containing fuels, such as petroleum fuels, natural gas, and biofuels, produce SO2, SO3, and other highly toxic gases upon combustion, which are harmful to human health and the environment, making it essential to understand their thermochemical properties. This study used high-level quantum chemistry calculations to determine thermodynamic parameters, including entropy, enthalpy, and specific heat capacity for an extensive set of sulfur-containing species. The B3LYP/cc-pVTZ level of theory was used for geometry optimization, vibration frequency, and dihedral scan calculations. To determine an appropriate ab initio method for energy calculation, the Bland-Altman diagram, a statistical analysis method, was employed to visualize the 298 K enthalpy value between experimental data and three sets of ab initio methods: G3, CBS-QB3, and the average of G3 plus CBS-QB3. The CBS-QB3 method exhibited the highest accuracy and was eventually selected for the energy calculation in this study. Thermochemical property parameters were then calculated with the MultiWell program suite for all these sulfur-containing species, and the results were in good agreement with the thermochemical data of organic compounds and the National Institute of Standards and Technology Chemistry WebBook databases. The thermochemical property database established in this study is essential to studying sulfur-containing species in desulfurization.

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.

Dynamical signatures from competing, nonadiabatic fragmentation pathways of S-nitrosothiophenol

A joint experiment-theory study of the UV photolysis of S-nitrosothiophenol reveals competing photodissociation pathways that produce NO in its spin–orbit ground state and thiophenoxy radical in either its ground or excited electronic state.

Quantifying rival bond fission probabilities following photoexcitation: C-S bond fission in t-butylmethylsulfide

We illustrate a new, collision-free experimental strategy that allows determination of the absolute probabilities of rival bond fission processes in a photoexcited molecule - here t-butylmethylsulfide (BSM). The method combines single photon ('universal') ionization laser probe methods, simultaneous imaging of all probed fragments (multi-mass ion imaging) and the use of an appropriate internal calibrant (here dimethylsulfide). Image analysis allows quantification of the dynamics of the rival B-SM and BS-M bond fission processes following ultraviolet (UV) excitation of BSM and shows the former to be twice as probable, despite the only modest (∼2%) differences in the respective ground state equilibrium C-S bond lengths or bond strengths. Rationalising this finding should provide a stringent test of the two close-lying, coupled excited states of ¹A'' symmetry accessed by UV excitation in BSM and related thioethers, of the respective transition dipole moment surfaces, and of the geometry dependent non-adiabatic couplings that enable the rival C-S bond fissions.

Computational Structural Biology of S-nitrosylation of Cancer Targets

Nitric oxide (NO) plays an essential role in redox signaling in normal and pathological cellular conditions. In particular, it is well known to react in vivo with cysteines by the so-called S-nitrosylation reaction. S-nitrosylation is a selective and reversible post-translational modification that exerts a myriad of different effects, such as the modulation of protein conformation, activity, stability, and biological interaction networks. We have appreciated, over the last years, the role of S-nitrosylation in normal and disease conditions. In this context, structural and computational studies can help to dissect the complex and multifaceted role of this redox post-translational modification. In this review article, we summarized the current state-of-the-art on the mechanism of S-nitrosylation, along with the structural and computational studies that have helped to unveil its effects and biological roles. We also discussed the need to move new steps forward especially in the direction of employing computational structural biology to address the molecular and atomistic details of S-nitrosylation. Indeed, this redox modification has been so far an underappreciated redox post-translational modification by the computational biochemistry community. In our review, we primarily focus on S-nitrosylated proteins that are attractive cancer targets due to the emerging relevance of this redox modification in a cancer setting.

High Accuracy Quantum Chemical and Thermochemical Network Data for the Heats of Formation of Fluorinated and Chlorinated Methanes and Ethanes

Reliable heats of formation are reported for numerous fluorinated and chlorinated methane and ethane derivatives by means of an accurate thermochemical protocol, which involves explicitly correlated coupled-cluster calculations augmented with anharmonic, scalar relativistic, and diagonal Born-Oppenheimer corrections. The theoretical results, along with additional experimental data, are further enhanced with the help of the thermochemical network approach. For 28 species, out of 50, this study presents the best estimates, and discrepancies with previous reports are also highlighted. Furthermore, the effects of the less accurate theoretical data on the results yielded by thermochemical networks are discussed.

Ab initio kinetics of the HOSO2 + 3O2 → SO3 + HO2 reaction

The detailed kinetic mechanism of the HOSO₂ + ³O₂ reaction, which plays a pivotal role in the atmospheric oxidation of SO₂, was investigated using accurate electronic structure calculations and novel statistical thermodynamic/kinetic models. Explored using the accurate composite method W1U, the detailed potential energy surface (PES) revealed that the addition of O₂ to a HOSO₂ radical to form the adduct (HOSO₄) proceeds via a transition state with a slightly positive barrier (i.e., 0.7 kcal mol^-1 at 0 K). Such a finding compromises a long-term hypothesis about this channel of being a barrierless process. Moreover, the overall reaction was found to be slightly exothermic by 1.7 kcal mol^-1 at 0 K, which is in good agreement with recent studies. On the newly-constructed PES, the temperature- and pressure-dependent behaviors of the title reaction were characterized in a wide range of conditions (T = 200-1000 K & P = 10-760 Torr) using the integrated deterministic and stochastic master equation/Rice-Ramsperger-Kassel-Marcus (ME/RRKM) rate model in which corrections for hindered internal rotation (HIR) and tunneling treatments were included. The calculated numbers were found to be in excellent agreement with literature data. The sensitivity analyses on the derived rate coefficients with respect to the ab initio input parameters (i.e., barrier height and energy transfer) were also performed to further understand the kinetic behaviors of the title reaction. The detailed kinetic mechanism, consisting of thermodynamic and kinetic data (in NASA polynomial and modified Arrhenius formats, respectively), was also provided at different T & P for further use in the modeling/simulation of any related systems.

Theoretical study of the reactions of the hydroselenyl radical (HSe●) with the selenenic radical (HSeO●)

The formation of selenium species in some biological processes involves the generation of ionic and radical intermediates such as RSe^●, RSe^-, RSeO^●, and RSeO^-, among others. We performed a theoretical study of the possible mechanisms for the reaction of the two simplest Se radicals-the hydroselenyl (HSe^●) and selenenic (HSeO^●) radicals, in which the possible products, intermediates, and transition-state structures were investigated. Density functional theory (DFT) was applied at the B3LYP/6-311++G(3df,3pd) level and the Ahlrichs Coulomb fitting basis sets were employed with an effective core potential (ECP) for both Se atoms. The same procedure was used to calculate the electronic density. All calculations were also performed using the M06-2X functional, which describes weaker bonds better than B3LYP does. In the reaction of interest, the so-called CR complex (HSe····SeOH) is formed initially. After passing through the transition state TS1, cis-HSeSeOH is obtained as a product. If a low barrier is then overcome (passing through the transition state TS32), the trans-HSeSeOH species is obtained. The CR complex can also rearrange into the intermediate INT after overcoming the barrier presented by the transition state TS2. Additionally, the decomposition of INT to H₂O and ¹Se₂ is possible through another transition state. This reaction is not included in this study. We also observed a second possible route for the conversion of INT to one of the HSeSeOH species; this route occurs through two pathways (with transition states TS31 and TS32). A comparison of some of the results with those obtained for sulfur analogs along the same pathways is also presented in this work. Graphical abstract Electronic envelopes for HSeO^● and HSe^● radicals.

Accurate Theoretical Thermochemistry for Fluoroethyl Radicals

An accurate coupled-cluster (CC) based model chemistry was applied to calculate reliable thermochemical quantities for hydrofluorocarbon derivatives including radicals 1-fluoroethyl (CH₃-CHF), 1,1-difluoroethyl (CH₃-CF₂), 2-fluoroethyl (CH₂F-CH₂), 1,2-difluoroethyl (CH₂F-CHF), 2,2-difluoroethyl (CHF₂-CH₂), 2,2,2-trifluoroethyl (CF₃-CH₂), 1,2,2,2-tetrafluoroethyl (CF₃-CHF), and pentafluoroethyl (CF₃-CF₂). The model chemistry used contains iterative triple and perturbative quadruple excitations in CC theory, as well as scalar relativistic and diagonal Born-Oppenheimer corrections. To obtain heat of formation values with better than chemical accuracy perturbative quadruple excitations and scalar relativistic corrections were inevitable. Their contributions to the heats of formation steadily increase with the number of fluorine atoms in the radical reaching 10 kJ/mol for CF₃-CF₂. When discrepancies were found between the experimental and our values it was always possible to resolve the issue by recalculating the experimental result with currently recommended auxiliary data. For each radical studied here this study delivers the best heat of formation as well as entropy data.

Production and Photodissociation of the Methyl Perthiyl Radical

The photodissociation dynamics of the methyl perthiyl (CH3SS) radical are investigated via molecular beam photofragment translational spectroscopy, using "soft" electron ionization to detect the radicals and their photofragments. With this new capability, we have shown that CH3SS can be generated from flash pyrolysis of dimethyl trisulfide. Utilizing this source of radicals and the advantages afforded by soft electron ionization, we have reinvestigated the photodissociation dynamics of CH3SS at 248 nm, finding CH3S + S to be the dominant dissociation channel with CH3 + SS as a minor process. These results differ from previous work reported in our laboratory in which we found CH3 + SS and CH2S + SH as the main dissociation channels. The difference in results is discussed in light of our new capabilities for characterization of radical production.

Photodissociation dynamics of the methyl perthiyl radical at 248 and 193 nm using fast-beam photofragment translational spectroscopy

The photodissociation dynamics of the methyl perthiyl radical (CH3SS) have been investigated using fast-beam coincidence translational spectroscopy. Methyl perthiyl radicals were produced by photodetachment of the CH3SS(-) anion followed by photodissociation at 248 nm (5.0 eV) and 193 nm (6.4 eV). Photofragment mass distributions and translational energy distributions were measured at each dissociation wavelength. Experimental results show S atom loss as the dominant (96%) dissociation channel at 248 nm with a near parallel, anisotropic angular distribution and translational energy peaking near the maximal energy available to ground state CH3S and S fragments, indicating that the dissociation occurs along a repulsive excited state. At 193 nm, S atom loss remains the major fragmentation channel, although S2 loss becomes more competitive and constitutes 32% of the fragmentation. The translational energy distributions for both channels are very broad at this wavelength, suggesting the formation of the S2 and S atom products in several excited electronic states.

Reaction of silylene with sulfur dioxide: some gas-phase kinetic and theoretical studies

Time-resolved kinetic studies of the reaction of silylene, SiH2, with SO2 have been carried out in the gas phase over the temperature range 297-609 K, using laser flash photolysis to generate and monitor SiH2. The second order rate coefficients at 1.3 kPa (SF6 bath gas) fitted the Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-10.10 ± 0.06) + (3.46 ± 0.45 kJ mol(-1))/RT ln 10 where the uncertainties are single standard deviations. The collisional efficiency is 71% at 298 K, and in kinetic terms the reaction most resembles those of SiH2 with CH3CHO and (CH3)2CO. Quantum chemical calculations at the G3 level suggest a mechanism occurring via addition of SiH2 to one of the S=O double bonds leading to formation of the three-membered ring, thione-siloxirane which has a low energy barrier to ring expansion to yield the four-membered ring, 3-thia-2,4-dioxasiletane, the lowest energy adduct found on the potential energy (PE) surface. RRKM calculations, however, show that, if formed, this molecule would only be partially stabilised under the reaction conditions and the rate coefficients would be pressure dependent, in contrast with experimental findings. The G3 calculations reveal the complexity of possible intermediates and end products and taken together with the RRKM calculations indicate the most likely end products to be H2SiO + SO ((3)Σ(-)). The reaction is compared and contrasted with that of SiH2 + CO2.

Benchmarking Experimental and Computational Thermochemical Data: A Case Study of the Butane Conformers

Due to its crucial importance, numerous studies have been conducted to determine the enthalpy difference between the conformers of butane. However, it is shown here that the most reliable experimental values are biased due to the statistical model utilized during the evaluation of the raw experimental data. In this study, using the appropriate statistical model, both the experimental expectation values and the associated uncertainties are revised. For the 133-196 and 223-297 K temperature ranges, 668 ± 20 and 653 ± 125 cal mol(-1), respectively, are recommended as reference values. Furthermore, to show that present-day quantum chemistry is a favorable alternative to experimental techniques in the determination of enthalpy differences of conformers, a focal-point analysis, based on coupled-cluster electronic structure computations, has been performed that included contributions of up to perturbative quadruple excitations as well as small correction terms beyond the Born-Oppenheimer and nonrelativistic approximations. For the 133-196 and 223-297 K temperature ranges, in exceptional agreement with the corresponding revised experimental data, our computations yielded 668 ± 3 and 650 ± 6 cal mol(-1), respectively. The most reliable enthalpy difference values for 0 and 298.15 K are also provided by the computational approach, 680.9 ± 2.5 and 647.4 ± 7.0 cal mol(-1), respectively.