Lawrence Berkeley National Laboratory

Despite their great importance in low-temperature geochemistry, self-diffusion coefficients of noble gas isotopes in liquid water (D) have been measured only for the major isotopes of helium, neon, krypton and xenon. Data on the self-diffusion coefficients of minor noble gas isotopes are essentially non-existent and so typically are estimated by a kinetic theory model in which D varies as the inverse square root of the isotopic mass (m): D proportional to m-0.5. To examine the validity of the kinetic theory model, we performed molecular dynamics (MD) simulations of the diffusion of noble gases in ambient liquid water with an accurate set of noble gas-water interaction potentials. Our simulation results agree with available experimental data on the solvation structure and self-diffusion coefficients of the major noble gas isotopes in liquid water and reveal for the first time that the isotopic mass-dependence of all noble gas self-diffusion coefficients has the power-law form D proportional to m-beta with 0 < beta < 0.2. Thus our results call into serious question the widespread assumption that the square root model can be applied to estimate the kinetic fractionation of noble gas isotopes caused by diffusion in ambient liquid water.