Transport coefficients of electrolyte solutions from Smart Brownian dynamics simulations

We present results of Brownian dynamics simulations of aqueous 1-1 electrolyte solutions in the 1-molar concentration range. The electrical conductivity and the self-diffusion coefficients obtained from the simulations are compared to experimental data. The interaction potential between the ions is modeled by pairwise repulsive 1/rn soft-core interactions (n=9 or n=12) and Coulomb forces. We take into account hydrodynamic interactions and integrate the stochastic equations of motion with large time steps of about 100 femtoseconds, combined with an acceptance criterion known from the Smart Monte Carlo method. In this way, details of the dynamics of particles in close contact are not considered and the short-ranged repulsive forces act effectively as constraint forces preventing overlap configurations. The lengths of the performed simulations (about 10 nanoseconds) and the number of ions (216) allow to obtain single particle as well as collective transport coefficients with sufficient precision. For this purpose we use Kubo expressions which can be applied on the mesoscopic time scale of Brownian dynamics simulations. It is shown that hydrodynamic interactions must be taken into account to obtain agreement with the experimental data. They lower the electrical conductivity, as expected, but increase the self-diffusion coefficients, confirming a recent finding for colloids.We present results of Brownian dynamics simulations of aqueous 1-1 electrolyte solutions in the 1-molar concentration range. The electrical conductivity and the self-diffusion coefficients obtained from the simulations are compared to experimental data. The interaction potential between the ions is modeled by pairwise repulsive 1/rn soft-core interactions (n=9 or n=12) and Coulomb forces. We take into account hydrodynamic interactions and integrate the stochastic equations of motion with large time steps of about 100 femtoseconds, combined with an acceptance criterion known from the Smart Monte Carlo method. In this way, details of the dynamics of particles in close contact are not considered and the short-ranged repulsive forces act effectively as constraint forces preventing overlap configurations. The lengths of the performed simulations (about 10 nanoseconds) and the number of ions (216) allow to obtain single particle as well as collective transport coefficients with sufficient precision. For this pu...

[1]  L. Onsager,et al.  THEORIES AND PROBLEMS OF LIQUID DIFFUSION , 1945, Annals of the New York Academy of Sciences.

[2]  Theodore Shedlovsky,et al.  THE ELECTROLYTIC CONDUCTIVITY OF SOME UNI-UNIVALENT ELECTROLYTES IN WATER AT 25° , 1932 .

[3]  P. Turq,et al.  Acoustophoresis Revisited. 1. Electrolyte Solutions , 1995 .

[4]  Sedimentation of Clusters of Spheres , 1996 .

[5]  Louis J. Durlofsky,et al.  Dynamic simulation of hydrodynamically interacting particles , 1987, Journal of Fluid Mechanics.

[6]  H. Friedman,et al.  Theory of conductance and related isothermal transport coefficients in electrolytes , 1983 .

[7]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

[8]  Thomas S. Thacher,et al.  Theory of Onsager phenomenological coefficients for isothermal linear transport processes in electrolyte solutions , 1984 .

[9]  Lars Onsager,et al.  Conductance of Unassociated Electrolytes. , 1957 .

[10]  Anthony J. C. Ladd,et al.  Hydrodynamic transport coefficients of random dispersions of hard spheres , 1990 .

[11]  H. Friedman,et al.  Self-diffusion and distinct diffusion of ions in solution , 1988 .

[12]  D. Ermak,et al.  Brownian dynamics with hydrodynamic interactions , 1978 .

[13]  J. Perram,et al.  Simulation of electrostatic systems in periodic boundary conditions. I. Lattice sums and dielectric constants , 1980, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[14]  Hans Christian Öttinger,et al.  Stochastic Processes in Polymeric Fluids , 1996 .

[15]  J. García de la Torre,et al.  Hydrodynamic properties of complex, rigid, biological macromolecules: theory and applications , 1981, Quarterly Reviews of Biophysics.

[16]  W. Ebeling,et al.  Generalizations of Onsagers Semiphenomenological Theory of Electrolytic Conductance , 1978 .

[17]  Gerald R. Kneller,et al.  nMOLDYN: A program package for a neutron scattering oriented analysis of Molecular Dynamics simulations , 1995 .

[18]  Harold L. Friedman,et al.  Brownian dynamics: Its application to ionic solutions , 1977 .

[19]  D. Ermak A computer simulation of charged particles in solution. I. Technique and equilibrium properties , 1975 .

[20]  J. Padró,et al.  Langevin dynamics simulation of ions in solution: Influence of the solvent structure , 1989 .

[21]  E. R. Smith,et al.  Hydrodynamic interactions in Brownian dynamics I. Periodic boundary conditions for computer simulations , 1987 .

[22]  Konrad Hinsen,et al.  Friction and mobility of many spheres in Stokes flow , 1994 .

[23]  R. Mills The Self-diffusion of Chloride Ion in Aqueous Alkali Chloride Solutions at 25° , 1957 .

[24]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[25]  S. Prager,et al.  Variational Treatment of Hydrodynamic Interaction in Polymers , 1969 .

[26]  B. Cichocki The generalized Smoluchowski equation for interacting Brownian particles with hard cores , 1987 .

[27]  M. Chemla,et al.  Coefficients d’autodiffusion dans les solutions aqueuses de KCl et de LiCl , 1971 .

[28]  B. U. Felderhof,et al.  Linear response theory of sedimentation and diffusion in a suspension of spherical particles , 1983 .

[29]  H. Friedman,et al.  Self‐diffusion coefficients of ions in electrolyte solutions by nonequilibrium Brownian dynamics , 1990 .

[30]  R. Stokes,et al.  Conductances of Concentrated Aqueous Sodium and Potassium Chloride Solutions at 25 , 1956 .

[31]  C. Mou,et al.  A statistical mechanical theory of the self‐diffusion coefficient of simple ions in electrolyte solutions , 1983 .

[32]  W. Kunz,et al.  Self-diffusion in electrolyte solutions using the mean spherical approximation , 1992 .

[33]  J. Happel,et al.  Low Reynolds number hydrodynamics , 1965 .

[34]  J. Padró,et al.  Langevin dynamics study of NaCl electrolyte solutions at different concentrations , 1990 .

[35]  A. Altenberger Quasihydrodynamic theory of ionic conductance in spatially homogeneous solutions of strong electrolytes , 1981 .

[36]  D. Heyes Mean-field-hydrodynamics Brownian dynamics simulations of viscosity and self-diffusion of near-hard-sphere colloidal liquids , 1995 .

[37]  Peter Baur,et al.  Long-time dynamics of charged colloidal suspensions: hydrodynamic interaction effects , 1997 .

[38]  Klaus Zahn,et al.  Hydrodynamic Interactions May Enhance the Self-Diffusion of Colloidal Particles , 1997 .

[39]  J. D. Doll,et al.  Brownian dynamics as smart Monte Carlo simulation , 1978 .

[40]  Olivier Bernard,et al.  Conductance in electrolyte solutions using the mean spherical approximation , 1992 .

[41]  H. Risken Fokker-Planck Equation , 1984 .

[42]  H. Friedman,et al.  Brownian Dynamics Simulation of Models for Aqueous 1M NaCl and 1/2M CuSO4 , 1987 .

[43]  Harold L. Friedman,et al.  Study of a Refined Model for Aqueous 1‐1 Electrolytes , 1971 .

[44]  Lars Onsager,et al.  Irreversible Processes in Electrolytes. Diffusion, Conductance and Viscous Flow in Arbitrary Mixtures of Strong Electrolytes , 1931 .