Calculation of the structures, stabilities, and vibrational spectra of arsenites, thioarsenites and thioarsenates in aqueous solution

Structures, stabilities and vibrational spectra have been calculated using molecular quantum mechanical methods for As(OH)₃, AsO(OH)₃, As(SH)₃, AsS(SH)₃ and their conjugate bases and for several species with partial substitution of S for O. Properties for the neutral gas-phase molecules are calculated with state-of-the-art methods which yield AsL distances within 0. 01 Å and AsL stretching frequencies within 10 cm⁻¹ of experiment. Similar accuracy is obtained for neutral molecules in solution using a polarizable continuum model (PCM). For monoanions such as and frequencies can be calculated to within 20 cm⁻¹ of experiment using the polarizable continuum model. Multiply charged anions remain a challenge for accurate frequency calculations, but we have obtained results within the PCM model which at least semiquantitatively reproduce the available data. This allows us to assign the controversial features D, E and F in the Raman data of (Wood S. A., Tait C. D. and Janecky D. R. (2002) A Raman spectroscopic study of arsenite and thioarsenite species in aqueous solution at 25 °C. Geochem. Trans. 3, 31-39).To help in the assignment of the arsenic sulfide spectra we have also calculated energetics for the oxidation of As(III) to As(V) compounds by polysulfides, disproportionation of As(III) compounds and for the dissociation of the oxo- and thio-acids. We have determined that As(III) oxyacids can be transformed to thioacids which can in turn be oxidized to As(V) sulfides by polysulfides and that the pK_a1s of the acids involved can be ordered as follows: AsS(SH)₃ < As(SH)₃ < AsO(OH)₃ < As(OH)₃ in order of increasing pK_a1. We have also established from the calculated energies that the most stable form of the As(III) oxyacid in acidic aqueous solution is indeed As(OH)₃, consistent with previous assignments.