A unifying mode-coupling theory for transport properties of electrolyte solutions. I. General scheme and limiting laws

On analytical theories for conductivity and self-diffusion in concentrated electrolytes

Describing analytically the transport properties of electrolytes, such as their conductivity or the self-diffusion of the ions, has been a central challenge of chemical physics for almost a century. In recent years, this question has regained some interest in light of Stochastic Density Field Theory (SDFT) - an analytical framework that allows the approximate determination of density correlations in fluctuating systems. In spite of the success of this theory to describe dilute electrolytes, its extension to concentrated solutions raises a number of technical difficulties, and requires simplified descriptions of the short-range repulsion between the ions. In this article, we discuss recent approximations that were proposed to compute the conductivity of electrolytes, in particular truncations of Coulomb interactions at short distances. We extend them to another observable (the self-diffusion coefficient of the ions) and compare them to earlier analytical approaches, such as the mean spherical approximation and mode-coupling theory. We show how the treatment of hydrodynamic effects in SDFT can be improved, that the choice of the modified Coulomb interactions significantly affects the determination of the properties of the electrolytes, and that comparison with other theories provides a guide to extend SDFT approaches in this context.

Electrical noise in electrolytes: a theoretical perspective

Seemingly unrelated experiments such as electrolyte transport through nanotubes, nano-scale electrochemistry, NMR relaxometry and surface force balance measurements, all probe electrical fluctuations: of the electric current, the charge and polarization, the field gradient (for quadrupolar nuclei) and the coupled mass/charge densities. The fluctuations of such various observables arise from the same underlying microscopic dynamics of the ions and solvent molecules. In principle, the relevant length and time scales of these dynamics are encoded in the dynamic structure factors. However, modelling the latter for frequencies and wavevectors spanning many orders of magnitude remains a great challenge to interpret the experiments in terms of physical processes such as solvation dynamics, diffusion, electrostatic and hydrodynamic interactions between ions, interactions with solid surfaces, etc. Here, we highlight the central role of the charge-charge dynamic structure factor in the fluctuations of electrical observables in electrolytes and offer a unifying perspective over a variety of complementary experiments. We further analyze this quantity in the special case of an aqueous NaCl electrolyte, using simulations with explicit ions and an explicit or implicit solvent. We discuss the ability of the standard Poisson-Nernst-Planck theory to capture the simulation results, and how the predictions can be improved. We finally discuss the contributions of ions and water to the total charge fluctuations. This work illustrates an ongoing effort towards a comprehensive understanding of electrical fluctuations in bulk and confined electrolytes, in order to enable experimentalists to decipher the microscopic properties encoded in the measured electrical noise.

Theoretical description of the electrical double layer for a mixture of n ionic species with arbitrary size and charge asymmetries. I. Spherical geometry

In this work, we propose a theoretical finite element description of the ionic profiles of a general mixture of n species of spherical charged particles dissolved in an implicit solvent, with arbitrary size and charge asymmetries, neutralizing a spherical macroion. This approach aims to close the gap between the nano- and the micro-scales in macroion solutions, taking into account the ion correlations and ionic excluded volume effects consistently. When these last two features are neglected, the classical non-linear Poisson-Boltzmann theory for n ionic species-with different ionic closest approach distances to the colloidal surface-is recovered as a limit case. As a proof of concept, we study the electrical double layer of an electroneutral mixture of oppositely charged colloids and small microions, with an asymmetry 1:333 in size and 1:10 in valence, in salt-free and added salt environments. Our theoretical approach displays a good agreement regarding the ionic profiles, the integrated charge, and the mean electrostatic potential obtained from molecular dynamics simulations with explicit-sized microions. Although the non-linear Poisson-Boltzmann colloid-colloid and colloid-microion profiles differ notably from those obtained via molecular dynamics simulations with explicit small-sized ions, the associated mean electrostatic potential agrees well with the corresponding explicit microion simulations.

Frequency and field-dependent response of confined electrolytes from Brownian dynamics simulations

Using Brownian dynamics simulations, we investigate the effects of confinement, adsorption on surfaces, and ion-ion interactions on the response of confined electrolyte solutions to oscillating electric fields in the direction perpendicular to the confining walls. Nonequilibrium simulations allows to characterize the transitions between linear and nonlinear regimes when varying the magnitude and frequency of the applied field, but the linear response, characterized by the frequency-dependent conductivity, is more efficiently predicted from the equilibrium current fluctuations. To that end, we (rederive and) use the Green-Kubo relation appropriate for overdamped dynamics, which differs from the standard one for Newtonian or underdamped Langevin dynamics. This expression highlights the contributions of the underlying Brownian fluctuations and of the interactions of the particles between them and with external potentials. Although already known in the literature, this relation has rarely been used to date, beyond the static limit to determine the effective diffusion coefficient or the DC conductivity. The frequency-dependent conductivity always decays from a bulk-like behavior at high frequency to a vanishing conductivity at low frequency due to the confinement of the charge carriers by the walls. We discuss the characteristic features of the crossover between the two regimes, most importantly how the crossover frequency depends on the confining distance and the salt concentration, and the fact that adsorption on the walls may lead to significant changes both at high and low frequencies. Conversely, our results illustrate the possibility to obtain information on diffusion between walls, charge relaxation, and adsorption by analyzing the frequency-dependent conductivity.

Assessment of the Wolf method using the Stillinger-Lovett sum rules: From strong electrolytes to weakly charged colloidal dispersions

The Ewald method has been the cornerstone in molecular simulations for modeling electrostatic interactions of charge-stabilized many-body systems. In the late 1990s, Wolf and collaborators developed an alternative route to describe the long-range nature of electrostatic interactions; from a computational perspective, this method provides a more efficient and straightforward way to implement long-range electrostatic interactions than the Ewald method. Despite these advantages, the validity of the Wolf potential to account for the electrostatic contribution in charged fluids remains controversial. To alleviate this situation, in this contribution, we implement the Wolf summation method to both electrolyte solutions and charged colloids with moderate size and charge asymmetries in order to assess the accuracy and validity of the method. To this end, we verify that the proper selection of parameters within the Wolf method leads to results that are in good agreement with those obtained through the standard Ewald method and the theory of integral equations of simple liquids within the so-called hypernetted chain approximation. Furthermore, we show that the results obtained with the original Wolf method do satisfy the moment conditions described by the Stillinger-Lovett sum rules, which are directly related to the local electroneutrality condition and the electrostatic screening in the Debye-Hückel regime. Hence, the fact that the solution provided by the Wolf method satisfies the first and second moments of Stillinger-Lovett proves, for the first time, the reliability of the method to correctly incorporate the electrostatic contribution in charge-stabilized fluids. This makes the Wolf method a powerful alternative compared to more demanding computational approaches.

Accurate modeling of a biological nanopore with an extended continuum framework

Despite the broad success of biological nanopores as powerful instruments for the analysis of proteins and nucleic acids at the single-molecule level, a fast simulation methodology to accurately model their nanofluidic properties is currently unavailable. This limits the rational engineering of nanopore traits and makes the unambiguous interpretation of experimental results challenging. Here, we present a continuum approach that can faithfully reproduce the experimentally measured ionic conductance of the biological nanopore Cytolysin A (ClyA) over a wide range of ionic strengths and bias potentials. Our model consists of the extended Poisson-Nernst-Planck and Navier-Stokes (ePNP-NS) equations and a computationally efficient 2D-axisymmetric representation for the geometry and charge distribution of the nanopore. Importantly, the ePNP-NS equations achieve this accuracy by self-consistently considering the finite size of the ions and the influence of both the ionic strength and the nanoscopic scale of the pore on the local properties of the electrolyte. These comprise the mobility and diffusivity of the ions, and the density, viscosity and relative permittivity of the solvent. Crucially, by applying our methodology to ClyA, a biological nanopore used for single-molecule enzymology studies, we could directly quantify several nanofluidic characteristics difficult to determine experimentally. These include the ion selectivity, the ion concentration distributions, the electrostatic potential landscape, the magnitude of the electro-osmotic flow field, and the internal pressure distribution. Hence, this work provides a means to obtain fundamental new insights into the nanofluidic properties of biological nanopores and paves the way towards their rational engineering.

Persistent Anti-Correlations in Brownian Dynamics Simulations of Dense Colloidal Suspensions Revealed by Noise Suppression

Transport properties of a hard-sphere colloidal fluid are investigated by Brownian dynamics simulations. We implement a novel algorithm for the time-dependent velocity-autocorrelation function (VACF) essentially eliminating the noise of the bare random motion. The measured VACF reveals persistent anti-correlations manifested by a negative algebraic power-law tail t^{-5/2} at all densities. At small packing fractions the simulations fully agree with the analytic low-density prediction, yet the amplitude of the tail becomes dramatically suppressed as the packing fraction is increased. The mode-coupling theory of the glass transition provides a qualitative explanation for the strong variation in terms of the static compressibility as well as the slowing down of the structural relaxation.

Anomalies in the equilibrium and nonequilibrium properties of correlated ions in complex molecular environments

Emergent statistical attributes, and therefore the equations of state, of an assembly of interacting charge carriers embedded within a complex molecular environment frequently exhibit a variety of anomalies, particularly in the high-density (equivalently, the concentration) regime, which are not well understood, because they do not fall under the low-concentration phenomenologies of Debye-Hückel-Onsager and Poisson-Nernst-Planck, including their variants. To go beyond, we here use physical concepts and mathematical tools from quantum scattering theory, transport theory with the Stosszahlansatz of Boltzmann, and classical electrodynamics (Lorentz gauge) and obtain analytical expressions both for the average and the frequency-wave vector-dependent longitudinal and transverse current densities, diffusion coefficient, and the charge density, and therefore the analytical expressions for (a) the chemical potential, activity coefficient, and the equivalent conductivity for strong electrolytes and (b) the current-voltage characteristics for ion-transport processes in complex molecular environments. Using a method analogous to the notion of Debye length and thence the electrical double layer, we here identify a pair of characteristic length scales (longitudinal and the transverse), which, being wave vector and frequency dependent, manifestly exhibit nontrivial fluctuations in space-time. As a unifying theme, we advance a quantity (inverse length dimension), g_{scat}^{(a)}, which embodies all dynamical interactions, through various quantum scattering lengths, relevant to molecular species a, and the analytical behavior which helps us to rationalize the properties of strong electrolytes, including anomalies, in all concentration regimes. As an example, the behavior of g_{scat}^{(a)} in the high-concentration regime explains the anomalous increase of the Debye length with concentration, as seen in a recent experiment on electrolyte solutions. We also put forth an extension of the standard diffusion equation, which manifestly incorporates the effects arising from the underlying microscopic collisions among constituent molecular species. Furthermore, we show a nontrivial connection between the current-voltage characteristics of electrolyte solutions and the Landauer's approach to electrical conduction in mesoscopic solids and thereby establish a definite conceptual bridge between the two disjoint subjects. For numerical insight, we present results on the aqueous solution of KCl as an example of strong electrolyte, and the transport (conduction as well as diffusion) of K^{+} ions in water, as an example of ion transport across the voltage-gated channels in biological cells.

An electric-field induced dynamical state in dispersions of charged colloidal rods

The response of concentrated dispersions of charged colloids to low-frequency electric fields is governed by field-induced inter-colloidal interactions resulting from the polarization of electric double layers and the layer of condensed ions, association and dissociation of condensed ions, as well as hydrodynamic interactions through field-induced electro-osmotic flow. The phases and states that can be formed by such field-induced interactions are an essentially unexplored field of research. Experiments on concentrated suspensions of rod-like colloids (fd-virus particles), within the isotropic-nematic phase coexistence region, showed that a number of phases/states are induced, depending on the field amplitude and frequency [Soft Matter, 2010, 6, 273]. In particular, a dynamical state is found where nematic domains form and melt on a time scale of the order of seconds. We discuss the microscopic origin of this dynamical state, which is attributed to the cyclic, electric-field induced dissociation and association of condensed ions. A semi-quantitative theory is presented for the dynamics of melting and formation of nematic domains, including a model for the field-induced dissociation/association of condensed ions. The resulting equation of motion for the orientational order parameter is solved numerically for parameters complying with the fd-virus system. A limit-cycle is found, with a cycling-time that diverges at the transition line in the field-amplitude versus frequency plane where the dynamical state first appears, in accord with experimental findings.