High-accuracy theoretical thermochemistry of atmospherically important sulfur-containing molecules

In this study, several sulfur-containing molecules with atmospherical importance were investigated by means of high-accuracy quantum chemical calculations including: HSO, HOS, HOSO2, HSNO, SH, CH2SO, CH2SH, S2COH, and SCSOH. After identifying the stable conformers of the molecules, a coupled-cluster-based composite model chemistry, which includes contributions up to quadruple excitations as well as corrections beyond the nonrelativistic and Born–Oppenheimer approximations, was applied to calculate the corresponding heat of formation (Δ(f)H(0)° and Δ(f)H(298)°) and entropy (S(298)°) values. In most of the cases, this study delivers more reliable estimates for the investigated thermodynamic properties than those reported in previous investigations. Our data also suggest that the experimental heats of formation associated with the HSO molecule are very likely to belong to its structural isomer, HOS. It is also confirmed by the calculated thermodynamic properties including standard reaction entropies, enthalpies, and equilibrium constants that, in the reaction CS2 + OH CS2OH, the SCSOH structural isomer is produced. It is also noted that the currently accepted Δ(f)H(0)°(S(gas)) = 274.73 ± 0.3 kJ/mol value is in need of revision, and based on a recent measurement, which is also confirmed by our computations, it is advised to update it to Δ(f)H(0)°(S(gas)) = 277.25 ± 0.3 kJ/mol.     

Isomers of OCS2: IR absorption spectra of OSCS and O(CS2) in solid Ar

Irradiation of an Ar matrix sample containing O(3) and CS(2) with a KrF excimer laser at 248 nm yielded new lines at 1402.1 (1404.7), 1056.2 (1052.7), and 622.3 (620.5) cm(-1); numbers in parentheses correspond to species in a minor matrix site. Secondary photolysis at 308 nm diminished these lines and produced mainly OCS and SO(2). Annealing of this matrix to 30 K yielded a second set of new lines at 1824.7 and 617.8 cm(-1). The first set of lines are assigned to C=S stretching, O-S stretching, and S-C stretching modes of carbon disulfide S-oxide, OSCS; and the second set of lines are assigned to C=O stretching and OCS bending modes of dithiiranone, O(CS(2)), respectively, based on results of (34)S- and (18)O-isotopic experiments and quantum-chemical calculations. These calculations using density-functional theory (B3LYP/aug-cc-pVTZ) predict four stable isomers of OCS(2): O(CS(2)), SSCO, OSCS, and SOCS, listed in order of increasing energy. According to calculations, O(CS(2)) has a cyclic CS(2) moiety and is the most stable isomer of OCS(2). OSCS is planar, with bond angles angle OSC congruent with 111.9 degrees and angle SCS congruent with 177.3 degrees ; it is less stable than SSCO and O(CS(2)) by approximately 102 and 154 kJ mol(-1), respectively, and more stable than SOCS by approximately 26 kJ mol(-1). Calculated vibrational wave numbers, IR intensities, (34)S- and (18)O-isotopic shifts for OSCS and O(CS(2)) fit satisfactorily with experimental results.     

Atmospheric reaction of the HOSO radical with NO2: a theoretical study

The gas-phase reaction between HOSO and NO(2) was examined using density functional theory. Geometry optimizations and frequency computations were performed at the B3LYP/6-311++G(2df,2pd) level of theory for all minimum species and transition states. The ground-state potential energy surface, including activation energies and enthalpies, were calculated using the ab initio CBS-QB3 composite method. The results suggest that the addition of HOSO and NO(2) leads to two possible intermediates, HOS(O)NO(2) and HOS(O)ONO, without any energy barrier. The HOS(O)NO(2) easily decomposes into HONO + SO(2) through the low energy product complex HONO···SO(2), whereas the HOS(O)ONO dissociates to HOSO(2) + NO products. This latter dissociation is preferred from the isomerization of the HOS(O)ONO to HOS(NO)O(2). Also, HOS(O)NO(2) isomerization to HOS(O)ONO is hindered due to the presence of a large energy barrier. From the thermodynamic aspect, the main products in the title reaction are HONO + SO(2), whereas HOSO(2) + NO are expected as a minor products.     

Quantum chemical and theoretical kinetics study of the O(3P) + CS2 reaction

The triplet potential energy surface of the O((3)P) + CS(2) reaction is investigated by using various quantum chemical methods including CCSD(T), QCISD(T), CCSD, QCISD, G3B3, MPWB1K, BB1K, MP2, and B3LYP. The thermal rate coefficients for the formation of three major products, CS + SO ((3)Σ(-)), OCS + S ((3)P) and CO + S(2) ((3)Σ(-)(g)) were computed by using transition state and RRKM statistical rate theories over the temperature range of 200-2000 K. The computed k(SO + CS) by using high-level quantum chemical methods is in accordance with the available experimental data. The calculated rate coefficients for the formation of OCS + S ((3)P) and CO + S(2) ((3)Σ(-)(g)) are much lower than k(SO + CS); hence, it is predicted that these two product channels do not contribute significantly to the overall rate coefficient.     

氧原子与二硫化碳反应的机理

用从头计算法、内禀反应坐标和电子密度拓扑分析方法研究了3P态氧原子与二硫化碳的反应.找到了分别形成CS SO,S OCS和S2 CO三个反应通道上的极小点和过渡态.采用UHF/631G进行几何构型优化,并在UMP2/631G水平上进行能量校正.三个反应通道上的稳定点和过渡态都经内禀反应坐标(IRC)跟踪得以确认,并用电子密度拓扑分析方法考察了反应过程中化学键的变化.计算结果表明,反应过程中所有稳定点和过渡态都具有Cs对称性,即对每个反应通道而言在整个反应过程中分子始终保持在同一平面内.在三个反应通道中,第一个反应通道O CS2→CS SO由于具有较小的活化能而更容易发生,与实验结果相一致.文中对反应机理进行了较详细的讨论.     

氧气和CS自由基反应势能面的密度泛函理论研究

应用量子化学从头计算和密度泛函理论(DFT)对O_2和CS自由基的反应进行了研 究。在B3LYP/6-311G~(**)水平上计算出了各物种的优化构型、振动频率和零点振 动能(ZPVE)。各物种的总能量由CCSD（T）/6-311G~(**)//B3LYP/6-311G~(**)给出 ，并对总能量进行了零点能校正。计算结果表明：CS自由基中的C端沿着O_2的双键 中线方向进攻，进行加成反应，反应的第一步放出大量的热量(450 kJ/mol)，推动 反应继续进行，从稳定的中间体4(Cs)出发，反应主要通过O的相邻位置的迁移生成 P_1(CO+SO)和P_3(COS+O)；通过S的相邻位置的迁移生成了重要的反应复合物 (complex 1)，进一步离解为产物P_2(CO_2+S)。计算结果可以很好地解释反应机理 。     

Thermochemical stabilities and structures of the cluster ions OCS+, S2+, H+(OCS), and C2H5+ with OCS molecules in the gas phase

The gas-phase clustering reactions of OCS+, S2+, H+(OCS), and C2H5+ ions with carbonyl sulfide (OCS) molecules were studied using a pulsed electron-beam high-pressure mass spectrometer and applying density functional theory (DFT) calculations. In the cluster ions OCS+(OCS)(n) and H+(OCS)(OCS)(n), a moderately strong, here referred to as "semi-covalent", bond was formed with n = 1. However, the nature of bonding changed from semi-covalent to electrostatic with n = 1 --> 2. The bond energy of S2(+)(OCS) was determined experimentally to be 12.9 +/- 1 kcal/mol, which is significantly smaller than that of the isovalent S2(+)(CS2) complex (30.9 +/- 1.5 kcal/mol). DFT based calculations predicted the presence of several isomeric structures for H+(OCS)(OCS)(n) complexes. The bond energies in the C2H5+(OCS)(n) clusters showed an irregular decrease for n = 1 --> 2 and 7 --> 8. The nonclassical bridge structure for the free C2H5+ isomerized to form a semi-covalent bond with one OCS ligand, [H3CCH2...SCO]+, i.e., reverted to classical structure. However, the nonclassical bridge structure of C2H5+ was preserved in the cluster ions C2H5+(OCS)(n) below 140 K attributable to the lack of thermal energy for the isomerization. DFT calculations revealed that stability orders of the geometric isomers of H+(OCS)(OCS)(n) and C2H5+(OCS)(n) changed with increasing n values.     

Reactions of SO3 with the O/H radical pool under combustion conditions

The reactions of SO3 with H, O, and OH radicals have been investigated by ab initio calculations. For the SO3 + H reaction (1), the lowest energy pathway involves initial formation of HSO3 and rearrangement to HOSO2, followed by dissociation to OH + SO2. The reaction is fast, with k(1) = 8.4 x 10(9)T(1.22) exp(-13.9 kJ mol(-1)/RT) cm(3) mol(-1) s(-1) (700-2000 K). The SO3 + O --> SO2 + O2 reaction (2) may proceed on both the triplet and singlet surfaces, but due to a high barrier the reaction is predicted to be slow. The rate constant can be described as k(2) = 2.8 x 10(4)T(2.57) exp(-122.3 kJ mol(-1)/RT) cm(3) mol(-1) s(-1) for T > 1000 K. The SO3 + OH reaction to form SO2 + HO2 (3) proceeds by direct abstraction but is comparatively slow, with k(3) = 4.8 x 10(4)T(2.46) exp(-114.1 kJ mol(-) 1/RT) cm(3) mol(-1) s(-1) (800-2000 K). The revised rate constants and detailed reaction mechanism are consistent with experimental data from batch reactors, flow reactors, and laminar flames on oxidation of SO2 to SO3. The SO3 + O reaction is found to be insignificant during most conditions of interest; even in lean flames, SO3 + H is the major consumption reaction for SO3.     

An ab initio MO study on the structures and electronic states of hydrogen-bonded O3- HF and SO2-HF complexes

Ab initio molecular orbital (MO) calculations have been carried out for base-hydrogen fluoride (HF) complexes (base = O3 and SO2) in order to elucidate the structures and energetics of the complexes. The ab initio calculations were performed up to the QCISD(T)/6-311++G(d,p) level of theory. In both complexes, hydrogen-bonded structures where the hydrogen of HF orients toward one of the oxygen atoms of bases were obtained as stable forms. The calculations showed that cis and trans isomers exist in both complexes. All calculations for the SO2-HF complex indicated that the cis form is more stable in energy than the trans form. On the other hand, in O3-HF complexes, the stable structures are changed by the ab initio levels of theory used, and the energies of the cis and trans forms are close to each other. From the most sophisticated calculations (QCISD(T)/6-311++G(d,p)//QCISD/6-311+G(d) level), it was predicted that the complex formation energies for cis SO2-HF, trans SO2-HF, cis O3-HF, and trans O3-HF are 6.1, 5.7, 3.4, and 3.6 kcal/mol, respectively, indicating that the binding energy of HF to SO2 is larger than that of O3. The harmonic vibrational frequencies calculated for cis O3-HF and cis SO2-HF complexes were in good agreement with the experimental values measured by Andrews et al. Also, the calculated rotation constants for cis SO2-HF agreed with the experiment.     

Properties of fluorosulfate-based ionic liquids and geometries of (FO2SOH)OSO2F- and (FO2SOH)2O2SOF-

A room temperature ionic liquid (IL) based on the fluorosulfate anion (SO(3)F(-)) has been synthesized by the reaction of 1-ethyl-3-methylimidazolium (EMIm(+)) chloride and fluorosulfuric acid (HOSO(2)F). The viscosity, ionic conductivity, and electrochemical window of EMImSO(3)F at 25 °C are 46.6 mPa s, 10.8 mS cm(-1), and 4.3 V, respectively. According to a solvatochromic measurement using ILs, there is a trend in the donor ability of fluoro- and oxofluoroanions, PF(6)(-) < BF(4)(-) < N(SO(2)CF(3))(2)(-) < SO(3)CF(3)(-) < SO(3)F(-) < PO(2)F(2)(-), which is explained by the atomic charges obtained from quantum mechanical calculations. The 1 : 2 and 1 : 3 stoichiometric reactions of EMImCl and HOSO(2)F give EMIm(FO(2)SOH)OSO(2)F and EMIm(FO(2)SOH)(2)O(2)SOF, respectively. Both the salts are liquid at room temperature without a HOSO(2)F dissociation pressure (< 1 Pa) and have low viscosity and high ionic conductivity (9.2 mPa s and 30.8 mS cm(-1) for EMIm(FO(2)SOH)OSO(2)F and 5.1 mPa s and 43.2 mS cm(-1) for EMIm(FO(2)SOH)(2)O(2)SOF). The vibrational modes and bonding properties of these anionic species are discussed with the aid of quantum mechanical calculations. The (FO(2)SOH)OSO(2)F(-) anion in EMIm(FO(2)SOH)OSO(2)F does not have an inversion centre, which stands in contrast to the one with an inversion centre (e.g. observed in solid Cs(FO(2)SOH)OSO(2)F). The (FO(2)SOH)(2)O(2)SOF(-) anion in EMIm(FO(2)SOH)(2)O(2)SOF is characterized by vibrational spectroscopy under C(s) symmetry.     

Oxidation of CO by SO2: a theoretical study

The elementary reaction SO(2) + CO --> CO(2) + SO((3)Sigma) (1) and the subsequent reaction SO((3)Sigma) + CO --> CO(2) + S((3)P) (2) have been studied by the application of the Gaussian-3//B3LYP quantum chemical approach to characterize the potential energy surfaces and transition state kinetic analysis to derive rate coefficients. Reaction 1 is found to take place via two transition states (TS), a cis-OSOCO TS and a trans-OSOCO TS. Reaction via the cis-TS is concerted and takes place on a singlet surface. Intersystem crossing to the final products occurs after passage through the barrier on the singlet surface. The trans-TS leads to a very weakly bound singlet OSOCO intermediate that then passes through a second TS (on the triplet surface) to form the products. Reaction 2 takes place on triplet surfaces. There is a concerted reaction through a cis-SOCO TS and a weakly bound trans-SOCO has also been identified. Reaction 2 is analogous to the reaction CO + O(2)((3)Sigma) --> CO(2) + O((3)P) (3), and this reaction has been reinvestigated at a similar level of theory and the rate coefficient derived by quantum chemistry is compared with experiment. The sensitive effects of trace impurities such as H(2), H(2)O, and hydrocarbons on the accurate experimental determination of the rate coefficient of reaction 3 is discussed. Using rate coefficients for reactions 1 and 2 obtained via quantum chemical calculations, we have been unable to model the extent of decomposition of SO(2) measured in a shock tube study of reaction between SO(2) and CO [Bauer, S. H.; Jeffers, P.; Lifshitz, A.; Yadava, B. P. Proc. Combust. Inst. 1971, 13, 417]. In light of the known sensitivity of reaction 3 to trace impurities, we have incorporated trace amounts of H(2), CH(4), or H(2)O, together with our rate coefficients for (1) and (2), in a kinetic model of Alzueta et al. [Combust. Flame 2001, 127, 2234], which is then shown to be able to substantially model the SO(2) data of Bauer et al. In the course of this modeling study we also computed heats of formation for a number of sulfur-containing small molecules: HS, HSO, HSOH, HOSO, HS(2), HSO(2), HOSO(2), HOSOH, and HOSHO.     

Reactions of group 3 transition metal atoms with CS2 and OCS: Matrix isolation infrared spectra and density-functional calculations of SMCS, SM-(eta2-CS), SMCO, and SM-(eta2-CO) in solid argon

Laser-ablated scandium, yttrium, and lanthanum atoms were reacted with CS2 and OCS molecules in an argon matrix. Products of the type SMCX and S-M(eta2-CX) (X = S or O) were formed on sample deposition. Photolysis favored the S-M(eta(2)-CX) complex, while annealing increased the more stable SMCX isomer. Product absorptions are identified by density-functional frequency calculations and isotopic substitutions. This work reports the first vibrational spectroscopic characterization of Sc, Y, and La reaction products with CS2 and OCS and the subsequent interconversion between SMCX and S-M(eta2-CX) structural isomers.     

Mechanism of heterogeneous reaction of carbonyl sulfide on magnesium oxide

Heterogeneous reaction of carbonyl sulfide (OCS) on magnesium oxide (MgO) under ambient conditions was investigated by in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), quadrupole mass spectrometer (QMS), and density functional theory (DFT) calculations. It reveals that OCS can be catalytically hydrolyzed by surface hydroxyl on MgO to produce carbon dioxide (CO2) and hydrogen sulfide (H2S), and then H2S can be further catalytically oxidized by surface oxygen or gaseous oxygen on MgO to form sulfite (SO3(2-)) and sulfate (SO4(2-)). Hydrogen thiocarbonate (HSCO2-) was found to be the crucial intermediate. Surface hydrogen sulfide (HS), sulfur dioxide (SO2), and surface sulfite (SO3(2-)) were also found to be intermediates for the formation of sulfate. Furthermore, the surface hydroxyl contributes not only to the formation of HSCO2- but also to HSCO2- decomposition. On the basis of experimental results, the heterogeneous reaction mechanism of OCS on MgO was discussed.     

Ion chemistry of NOO+

The kinetics for the reactions of NOO+ ions with neutral molecules having ionization potentials (IPs) from 9.27 to 15.58 eV was measured in a selected ion flow tube at 298 K. The NOO+ ions are produced from the reaction of N3+ + O2 and have been reacted with the following: NO, C6F6, CS2, CF3I, C3F6, OCS, C2H6, Xe, SO2, O3, N2O, CO2, Kr, CO, D2, and N2. Numerous types of reactions were observed with the various neutral reagents, including production of NO+ (which may involve loss of an O from the ion or addition of O to the neutral reactant, although the two channels could not be distinguished here), charge transfer, isomerization of NOO+ to ONO+, and hydride abstraction. High level theoretical calculations of the structures and energetics of the various isomers, electronic states, and transition states of NOO and NOO+ were performed to better understand the observed reactivity. All neutral species with an IP< or =11.18 eV were observed to react with NOO+ in part by charge transfer. Detailed calculations showed that the recommended adiabatic and vertical IPs of NOO are 10.4 and 11.7 eV, respectively, at the MRCISDQ/AVQZ level of theory. The observed experimental limit for charge transfer of 11.18 eV agreed well with the energetics of the final products obtained from theory if dissociation of the neutral metastable product occurred, i.e., the products were X+ +[O(3P) + NO(2Pi)], where [O(3P)+NO(2Pi)] formed via dissociation of metastable NOO. Charge exchange with neutral reagent X would, therefore, be exothermic if IP(X)<[IPad(NOO)-DeltaE(O+NO)-NOO]= approximately 11.1 eV, where IPad(NOO) is the adiabatic IP. The potential energy surface for the reaction of NOO+ with C2H6 was also calculated, indicating that two pathways for formation of HNO2 + C2H5 (+) exist.     

Kinetics of the sulfur oxidation on palladium: a combined in situ x-ray photoelectron spectroscopy and density-functional study

We studied the reaction kinetics of sulfur oxidation on the Pd(100) surface by in situ high resolution x-ray photoelectron spectroscopy and ab initio density functional calculations. Isothermal oxidation experiments were performed between 400 and 500 K for small amounts (~0.02 ML) of preadsorbed sulfur, with oxygen in large excess. The main stable reaction intermediate found on the surface is SO(4), with SO(2) and SO(3) being only present in minor amounts. Density-functional calculations depict a reaction energy profile, which explains the sequential formation of SO(2), SO(3), and eventually SO(4), also highlighting that the in-plane formation of SO from S and O adatoms is the rate limiting step. From the experiments we determined the activation energy of the rate limiting step to be 85 ± 6 kJ mol(-1) by Arrhenius analysis, matching the calculated endothermicity of the SO formation.     

Al2O3上羰基硫常温催化水解的氧中毒机理

利用原位漫反射傅立叶变换红外光谱、X光衍射、BET、离子色谱(IC)等手段, 对Al2O3常温催化水解羰基硫(OCS)的氧中毒机理进行了研究. 实验表明, 表面—OH在OCS的催化水解反应中起关键作用, 表面HSCO-2物种是OCS催化水解反应的中间体. 有氧条件下, 利用原位红外光谱和离子色谱检测到了催化剂表面SO2-4的生成. SO2-4在催化剂表面积累是Al2O3上OCS常温催化水解氧中毒的主要原因.<     

Structure and reactivity of benzoylnitrene radical anion in the gas phase

The open-shell benzoylnitrene radical anion, readily generated by electron ionization of benzoylazide, undergoes unique chemical reactivity with radical reagents and Lewis acids in the gas phase. Reaction with nitric oxide, NO, proceeds by loss of N2 and formation of benzoate ion. This novel reaction is also observed to occur with phenylnitrene anion, forming phenoxide. Similar reactivity was observed in the reaction between benzoylnitrene radical anion and NO2, forming benzoate ion and nitrous oxide. Electronic structure calculations indicate that the reaction has a high-energy barrier that is overcome by the energy released by bond formation. Benzoylnitrene radical anion also transfers oxygen anion to NO and NO2 as well as to CS2 and SO2. In contrast, phenylnitrene anion reacts with carbon disulfide by C+ or CS+ abstraction, forming S- or S2-. Electronic structure calculations indicate that benzoylnitrene in the ground state resembles a slightly polarized benzoate anion, but with a free radical localized on the nitrogen.     

Infrared spectroscopy of SO2 aqueous solutions

The chemistry of SO2 solutions was studied with infrared transmission and total internal reflection (ATR) spectroscopy. The cross-sections of the SO2 stretching vibrations v1, v3 and the combination band of v1 + v3 were obtained and found to be slightly different from their gas phase values. Six features we associate with sulfur-containing ions were observed. We present the first infrared evidence for the bisulfite ion HOSO2- through detection of its OH stretching vibration. Other features may be associated with HSO3-, another isomer of the bisulfite ion, and a new compound we propose to be H2S2O6(2-). We found no evidence for sulfurous acid H2SO3.     

Selective two-step oxidation of μ2-S ligands in trigonal prismatic unit {Re3(μ6-C)(μ2-S)3Re3} of the bioctahedral cluster anion [Re12CS17(CN)6]6–

An oxidation of cluster anion [Re(12)CS(17)(CN)(6)](6-) by H(2)O(2) in water has been investigated. It was shown that selective two-step oxidation of bridging μ(2)-S-ligands in trigonal prismatic unit {Re(3)(μ(6)-C)(μ(2)-S)(3)Re(3)} takes place. The first stage runs rapidly, whereas the speed of the second stage depends on intensity of ultraviolet irradiation of the reaction mixture. Each stage of the reaction is accompanied by a change in the solution's color. In the first stage of the oxidation, the cluster anion [Re(12)CS(14)(SO(2))(3)(CN)(6)](6-) is produced, in which all bridging S-ligands are turned into bridging SO(2)-ligands. The second stage of the oxidation leads to formation of the anion [Re(12)CS(14)(SO(2))(2)(SO(3))(CN)(6)](6-), in which one of the SO(2)-ligands underwent further oxidation forming the bridging SO(3)-ligand. Seven compounds containing these anions were synthesized and characterized by a set of different methods, elemental analyses, IR and UV/vis spectroscopy, and quantum-chemical calculations. Structures of some compounds based on similar cluster anions, [Cu(NH(3))(5)](3)[Re(12)CS(14)(SO(2))(3)(CN)(6)]·9.5H(2)O, [Ni(NH(3))(6)](3)[Re(12)CS(14)(SO(2))(3)(CN)(6)]·4H(2)O, and [Cu(NH(3))(5)](2.6)[Re(12)CS(14)(SO(2))(3)(CN)(6)](0.6)[{Re(12)CS(14)(SO(2))(2)(SO(3))(CN)(5)(μ-CN)}{Cu(NH(3))(4)}](0.4)·5H(2)O, were investigated by X-ray analysis of single crystals.