DOE OSTI2016
We present theoretical calculations for all three isotope ratios of sulfur ( 33 S/ 32 S, 34 S/ 32 S, 36 S/ 32 S) at the B3LYP/6-31+G(d,p) level of theory for aqueous sulfur compounds modeled in 30–40H 2 O clusters spanning the range of sulfur oxidation state (S n , n = -2 to +6) for estimating equilibrium fractionation factors in aqueous systems. Computed 34 β values based on major isotope ( 34 S/ 32 S) reduced partition function ratios (RPFRs) scale to a first order with sulfur oxidation state and coordination, where 34 β generally increase with higher oxidation state and increasing coordination of the sulfur atom. Exponents defining mass dependent relationships based on β values ( x/34 κ = ln( x β)/ln( 34 β), x = 33 or 36) conform to tight ranges over a wide range of temperature for all aqueous compounds ( 33/34 κ ≈ 0.5148–0.5159, 36/34 κ ≈ 1.89–1.90 from T ≽ 0 °C). The exponents converge near a singular value for all compounds at the high temperature limit ( 33/34 κT→∞= 0.51587 ± 0.00003 and 36/34 κT→∞= 1.8905 ± 0.0002; 1 s.d. of all computed compounds), and typically follow trends based on oxidation state and coordination similar to those seen in 34 β values at lower temperatures. Theoretical equilibrium fractionation factors computed from these β-values are compared to experimental constraints for HSO 3 - T(aq)/SO 2(g, aq) , SO 2(aq) /SO 2(g) , H 2 S (aq) /H 2 S (g) , H 2 S (aq) /HS - (aq) , SO 4 2- (aq) /H 2 S T(aq) , S 2 O 3 2- (aq) (intramolecular), and S 2 O 3 2- (aq) /H 2 S T(aq) , and generally agree within a reasonable estimation of uncertainties. We make predictions of fractionation factors where other constraints are unavailable. Isotope partitioning of the isomers of protonated compounds in the sulfite and sulfoxylate systems depend strongly on whether protons are bound to either sulfur or oxygen atoms. The magnitude of the HSO 3 - T /SO 3 2- major isotope ( 34 S/ 32 S) fractionation factor is predicted to increase with temperature from 0 to 70 °C due to the combined effects of the large magnitude (HS)O 3 - /SO 3 2- fractionation factor (1000ln 34 α (HS)bisulfite-sulfite = 19.9‰, 25 °C) relative to the (HO)SO 2 - /SO 3 2- fractionation factor (1000ln 34 α (HO)bisulfite–sulfite = -2.2‰, 25 °C), and the increased stability of the (HS)O 3 - isomer with increasing temperature. Finally, we argue that isomerization phenomenon should be considered in models of the sulfur cycle, including models that describe the overall sulfur isotope fractionations associated with microbial metabolism (e.g., microbial sulfate reduction).
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗