Friction Coefficients of Ions in Aqueous Solution at 25 °C

Abstract

The diffusion coefficients of ions and of uncharged solutes in aqueous solution at 25 °C and at infinite dilution are studied by computer simulation using the SPC/E model for water and solute−water potentials employed in previous work (Koneshan, S.; et al. J. Phys. Chem. 1998, 102, 4193−4204). The mobilities of the ions calculated from the diffusion coefficients pass through a maximum as a function of ion size, with distinct curves and maximums for positive and negative ions in qualitative agreement with experiment. We aim to understand this at a microscopic level in terms of theoretical studies of the friction coefficient ζ, which is related to diffusion coefficient by the Stokes−Einstein relation. This provides one method of calculating ζ, but it can also be obtained in principle from the random force autocorrelation function which is the starting point of molecular theories of the friction. Molecular and continuum theories divide the friction into hydrodynamic and dielectric components calculated in different ways on the basis of certain approximations that are tested in this paper. The two methods of determining ζ give consistent results in simulations of large ions or uncharged solutes but differ by a factor of nearly 1.5 for smaller ions. This is attributed to the assumption, in our simulations and in some molecular theories, that the random force autocorrelation function of a moving ion can be approximated by the total force autocorrelation function of a fixed ion but the observed trends in ζ with ion size are unchanged. Three different separations of the force autocorrelation function are studied; namely, partitioning this into (a) electrostatic and Lennard-Jones forces (b) hard and soft forces, and (c) forces arising from the first shell and more distant forces. The cross-terms are found to be significant in all cases, and the contributions of the sum of the soft term and the cross-terms, which are of opposite sign, to the total friction in the hard−soft separation, is small for all the ions (large and small). This suggests that dielectric friction calculated using this separation, with the neglect of cross-terms, is less accurate than previously supposed and the success of these theories is due to a cancellation of errors in the approximations. This is supported by recent a theoretical study of Chong and Hirata (Chong, C.; Hirata, F. J. Chem. Phys. 1998, 108, 739) which evaluates the cross-terms. A phenomenological theory due to Chen and Adelman (Chen, J. H.; Adelman, S. A. J. Chem. Phys. 1980, 72, 2819) calculates the friction in terms of effective hydrodynamic and dielectric radii for the ions (cf. solventberg picture) and predicts low dielectric friction for large ions and small strongly solvated ions. This also agrees qualitatively with our simulations but a complete molecular theory, applicable to positive and negative ions in hydrogen-bonded solvents such as water, has yet to be developed.

Diffusion constant of K+ inside Gramicidin A: A comparative study of four computational methods

“…According to the second fluctuation dissipation theorem (SFDT), the memory function is related to the random force autocorrelation function (FACF) [20,32] according to (5) Using Einstein's relation and the connection between friction constant and memory function , then:…”

“…Therefore, one must rely on simulations to predict diffusion constants, and, indeed, several theoretical methods have been developed for this purpose. They are widely used for calculating diffusion properties of ions and molecules in bulk phases [20][21][22], but the applicability of some of these methods to narrow ion channels (e.g., Gramicidin A) is questionable. Currently, there is no consensus in the biophysics literature about the magnitude of diffusion constants of ions inside narrow channels [6,[23][24][25][26][27][28].…”

Diffusion coefficient of ionic solvation shell molecules

“…[1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16] Generally, their motion is found to be substantially slowed down with respect to bulk solvent molecules but still somewhat faster than that of the ion. This conclusion is usually founded on the computation of the diffusion coefficient for the subset of first solvation shell molecules: its value is larger than that of the ion and lower than that of the bulk.…”

Influence of methanol as a buffer additive on the mobilities of organic cations in capillary electrophoresis

“…Thus, the sign of the ion charge is apparently responsible for the variation in dielectric friction. Differences in behavior between anions and cations have been previously reported in the literature [22,[60][61][62][63][64]. The Hubbard-Onsager theory fails at predicting these differences because it is a continuum model, and as such, it does not consider ion-solvent interactions, solvent structure, or the nature of the solvent motion [65,66].…”