Molecular dynamics simulations of kinetic isotope fractionation during the diffusion of ionic species in liquid water

Interpretation of isotope ratios, a powerful tool in geochemical investigations of fluid-rock systems, requires an understanding of all relevant processes that fractionate isotopes. One such process, diffusion in liquid water, has remained problematic despite its potential significance as a major cause of kinetic isotope fractionation. Recent laboratory experiments published by [Richter, F. M., Mendybaev, R. A., Christensen, J. N., Hutcheon, I. D., Williams, R. W., Sturchio, N. C., and Beloso Jr., A. D. (2006) Kinetic isotopic fractionation during diffusion of ionic species in water. Geochim. Cosmochim. Acta 70, 277-289.] have shown clearly for the first time that lithium and chloride isotopes are fractionated by diffusion in liquid water, whereas magnesium isotopes are not. In the present paper, we present the results of molecular dynamics simulations of lithium, chloride, and magnesium diffusion in liquid water that were designed to provide molecular-scale insight into the experimental findings of Richter et al. (2006). Our results indicate that the self-diffusion coefficients of lithium, chloride, and magnesium isotopes follow an inverse power-law dependence on ion mass (D_i ∝ m_i^{- β}, where D_i is the self-diffusion coefficient of a solute with isotopic mass m_i). The power-law exponents (β) deduced for lithium, chloride, and magnesium from the diffusivity data of Richter et al. (2006) are consistent with the mass dependencies found in our simulations. Further analysis of our simulation results showed that the experimental β-values are inversely related to the residence times of water molecules in the first solvation shells of the diffusing ions, as expected from mode-coupling and renormalized kinetic theories.