Tsubasa Otake, Antonio C. Lasaga, Hiroshi Ohmoto
CHEMICAL GEOLOGY 249 (3-4) 357 - 376 0009-2541 2008/04
[Refereed][Not invited] Ab initio calculations of vibrational frequencies (v), reduced partition function ratios (0) of various gaseous and aqueous sulfur-bearing compounds, and isotopic fractionation factors (a) between them were performed on the multiple isotope systems using the Hartree-Fock (HF) and the Becke three-parameter Lee-Yang-Parr (B3LYP) density functional methods. Based on the comparison of the geometry, vibrational frequencies, and fractionation factors with those determined experimentally, the HF/6-31G(d) was chosen as the optimal method for the sulfur isotope systems. The results of computation for multiple sulfur fractionation factors indicate that the mass-dependent relationships, i.e., ((33)alpha - 1)/((34)alpha - 1) and ((36)alpha - 1)/((34)alpha - 1) values, converge only at high temperature (e.g., T>500 degrees C) to 0.515 and 1.89, respectively, which are the values previously derived from the Bigeleisen-Mayer approximations. However, they significantly deviate from these values with decreasing temperature, and depending on the sulfur species pairs that are involved in equilibrium isotope exchange reactions. The ((33)alpha- 1)/((34)alpha - 1) values and the ((36)alpha - 1)/ ((34)alpha- 1) values range from 0.505 to 0.517 and from 1.88 to 1.96, respectively, at 0 degrees C.
An isotopic crossover relationship, in which the fractionation factor changes from > 1.0 to < 1.0 (or from < 1.0 to > 1.0) about a specific temperature (i.e., crossover temperature), was observed in several species couples (e.g., S-2-H2S, S-8-H2S, CS2-H2S, SO3-SO4(aq)2-). Since the crossover temperatures vary among different isotope molecules, large changes in the ((33)alpha -1)/((34)alpha - 1) and ((36)alpha - 1)/((34)alpha - 1) values (e.g., > 1 to <0) also occur near the crossover temperatures.
The effects of solvation on the fractionation factors were investigated using the PCM models. The dielectric constant (P) of solvent, which is the dominant solvent characteristic, affects the fractionation factors between gaseous and aqueous species (e.g., SO2(g) and SO2(aq)), but the effect levels out for epsilon > 20. Among the solvent models, the IEF-PCM model yields the greatest fractionation factors between gaseous and aqueous species.
Overall, temperature, harmonic frequency, and mass of atoms bonded to sulfur (M are the three important parameters affecting the reduced partition function ratios, fractionation factors, and mass-dependent relationships. By decreasing temperature and increasing frequency and mass, the 1000((34)beta - 1) values increase and the ((33)beta - 1)/((34)beta - 1) values decrease (i.e., deviate from "expected" 0.515). The ((33)beta - 1)/((34)beta - 1) values converge to 0.516 only at higher temperature (e.g., T>600 degrees C), lower frequency (omega(e) <500 cm(-1)), and smaller mass (M<7 amu) and become as low as 0.503 at T=50 degrees C, omega(e) = 2500 cm(-1), and M=56 amu. (C) 2008 Elsevier B.V All rights reserved.