On the aqueous solvation of AsO(OH)3vs. As(OH)3. Born-Oppenheimer molecular dynamics density functional theory cluster studiesElectronic supplementary information (ESI) available: Graphs showing the convergence of the vibrational spectra, images illustrating the frames of reference to define the lower and upper hemispheres used in the analysis of the hydration patterns, graphs with the cumulative numbers of oxygens, ∑NO(r), and hydrogens ∑NH(r), and graphs where the As-O distances are depicted.

While arsenous acid, As(OH) 3 , has been the subject of a plethora of studies due to its worldwide ubiquity and its toxicity, pentavalent As in the form of arsenic acid, AsO(OH) 3 , has recently been found in rivers in central Mexico as the most abundant naturally occurring arsenic species. To better understand the solvation patterns of both toxic acids at the molecular level, we report the results of Born-Oppenheimer molecular dynamics simulations on the aqueous solvation of the AsO(OH) 3 and As(OH) 3 molecules at room temperature using the cluster microsolvation approach including 30 water molecules at the B3LYP/6-31G** level of theory. We found that the average per-molecule water binding energy is ca. 1 kcal mol −1 larger for the As( v ) species as compared to the As( iii ) one. To account for the asymmetry of both molecules, the hydration patterns were studied separately for a "lower" hemisphere, defined by the initially protonated oxygens, and for the opposite "upper" hemisphere. Similar lower hydration patterns were found for both As( iii ) and As( v ), with the same coordination number CN = 7. The upper pattern for As( iii ) was found to be of a hydrophobic type, whereas that for As( v ) showed the fourth oxygen to be hydrogen-bonded to the water network, yielding CN = 3.7; moreover, a proton "hopped" from the lower to the upper side, through the Grotthuss mechanism. Theoretical EXAFS spectra were obtained that showed good agreement with experimental data for As( iii ) and As( v ) in liquid water, albeit with somewhat longer As-O distances due to the level of theory employed. Proton transfer processes were also addressed; we found that the singly deprotonated H 2 AsO 3 − species largely dominated (99% of the simulation) for the As( iii ) case, and that the deprotonated H 2 AsO 4 − and HAsO 4 2− species were almost equally present (45% and 55%, respectively) for the As( v ) case, which is in line with the experimental data p K a1 = 2.24 and p K a2 = 6.96. Through vibrational analysis the features of the Eigen and Zundel ions were found in the spectra of the microsolvated As( iii ) and As( v ) species, in good agreement with experimental data in aqueous solutions. BOMD simulations were used to reveal the hydration features of As(OH) 3 and (for the first time) AsO(OH) 3 in aqueous solution.