Effective screening length and quasiuniversality for the restricted primitive model of an electrolyte solution

Systematic Incorporation of Ionic Hard-Core Size into the Debye-Hückel Theory via the Cumulant Expansion of the Schwinger-Dyson Equations

The Debye-Hückel (DH) formalism of bulk electrolytes equivalent to the Gaussian-level closure of the electrostatic Schwinger-Dyson identities without the interionic hard-core (HC) coupling is extended via the cumulant treatment of these equations augmented by HC interactions. By comparing the monovalent ion activity and pressure predictions of our cumulant-corrected DH (CCDH) theory with hypernetted-chain results and Monte Carlo simulations from the literature, we show that this rectification extends the accuracy of the DH formalism from submolar to molar salt concentrations. In the case of internal energies or the general case of divalent electrolytes mainly governed by charge correlations, the improved accuracy of the CCDH theory is limited to submolar ion concentrations. Comparison with experimental data from the literature shows that, via the adjustment of the hydrated ion radii, CCDH formalism can equally reproduce the nonuniform effect of salt increment on the ionic activity coefficients up to molar concentrations. The inequality satisfied by these HC sizes coincides with the cationic branch of the Hofmeister series.

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.

Electric-field-induced oscillations in ionic fluids: a unified formulation of modified Poisson-Nernst-Planck models and its relevance to correlation function analysis

We theoretically investigate an electric-field-driven system of charged spheres as a primitive model of concentrated electrolytes under an applied electric field. First, we provide a unified formulation for the stochastic charge and density dynamics of the electric-field-driven primitive model using the stochastic density functional theory (DFT). The stochastic DFT integrates the four frameworks (the equilibrium and dynamic DFTs, the liquid state theory and the field-theoretic approach), which allows us to justify in a unified manner various modifications previously made for the Poisson-Nernst-Planck model. Next, we consider stationary density-density and charge-charge correlation functions of the primitive model with a static electric field. We predict an electric-field-induced synchronization between emergences of density and charge oscillations. We are mainly concerned with the emergence of stripe states formed by segregation bands transverse to the external field, thereby demonstrating the following: (i) the electric-field-induced crossover occurs prior to the conventional Kirkwood crossover without an applied electric field, and (ii) the ion concentration dependence of the decay lengths at the onset of oscillations bears a similarity to the underscreening behavior found by recent simulation and theoretical studies on equilibrium electrolytes. Also, the 2D inverse Fourier transform of the correlation function illustrates the existence of stripe states beyond the electric-field-induced Kirkwood crossover.

Rigorous treatment of pairwise and many-body electrostatic interactions among dielectric spheres at the Debye-Hückel level

Electrostatic interactions among colloidal particles are often described using the venerable (two-particle) Derjaguin-Landau-Verwey-Overbeek (DLVO) approximation and its various modifications. However, until the recent development of a many-body theory exact at the Debye-Hückel level (Yu in Phys Rev E 102:052404, 2020), it was difficult to assess the errors of such approximations and impossible to assess the role of many-body effects. By applying the exact Debye-Hückel level theory, we quantify the errors inherent to DLVO and the additional errors associated with replacing many-particle interactions by the sum of pairwise interactions (even when the latter are calculated exactly). In particular, we show that: (1) the DLVO approximation does not provide sufficient accuracy at shorter distances, especially when there is an asymmetry in charges and/or sizes of interacting dielectric spheres; (2) the pairwise approximation leads to significant errors at shorter distances and at large and moderate Debye lengths and also gets worse with increasing asymmetry in the size of the spheres or magnitude or placement of the charges. We also demonstrate that asymmetric dielectric screening, i.e., the enhanced repulsion between charged dielectric bodies immersed in media with high dielectric constant, is preserved in the presence of free ions in the medium.

A multiple decay-length extension of the Debye-Hückel theory: to achieve high accuracy also for concentrated solutions and explain under-screening in dilute symmetric electrolytes

The Poisson-Boltzmann and Debye-Hückel approximations for the pair distributions and mean electrostatic potential in electrolytes predict that these entities have one single decay mode with a decay length equal to the Debye length 1/κD, that is, they have a characteristic contribution that decays with distance r like e-κDr/r. However, in reality, electrolytes have several decay modes e-κr/r, e-κ'r/r etc. with different decay lengths, 1/κ, 1/κ' etc., that in general are different from the Debye length. As an illustration of the significance of multiple decay modes in electrolytes, the present work uses a very simple extension of the Debye-Hückel approximation with two decay lengths, which predicts oscillatory modes when appropriate. This approach gives very accurate results for radial distribution functions and thermodynamic properties of aqueous solutions of monovalent electrolytes for all concentrations investigated, including high ones. It is designed to satisfy necessary statistical mechanical conditions for the distributions. The effective dielectric permittivity of the electrolyte plays an important role in the theory and each mode has its own value of this entity. Electrolytes with high electrostatic coupling, like those with multivalent ions and/or with solvent of low dielectric constant, have decay lengths in dilute solutions that substantially deviate from the Debye length. It is shown that this is caused by nonlinear ion-ion correlation effects and the origin of under-screening, i.e., 1/κ > 1/κD, in dilute symmetric electrolytes is analyzed. The under-screening is accompanied by an increase in the effective dielectric permittivity that is also caused by these correlations. The theoretical results for the decay length are successfully compared with recent experimental data for simple electrolytes in various solvents. The paper includes background material on electrolyte theory and screening in order to be accessible for nonexperts in the field.

Correlation Length in Concentrated Electrolytes: Insights from All-Atom Molecular Dynamics Simulations

We study the correlation length of the charge-charge pair correlations in concentrated electrolyte solutions by means of all-atom, explicit-solvent molecular dynamics simulations. We investigate LiCl and NaI in water, which constitute highly soluble, prototypical salts for experiments, as well as two more complex, molecular electrolyte systems of lithium bis(trifluoromethane)sulfonimide (LiTFSI), a salt commonly employed in electrochemical storage systems, in water, and in an organic solvent mixture of dimethoxyethane and dioxolane. Our simulations support the recent experimental observations as well as theoretical predictions of a nonmonotonic behavior of the correlation length with increasing salt concentration. We observe a Debye-Hückel like regime at low concentration, followed by a minimum reached when d/λ_D ≃ 1, where λ_D is the Debye correlation length and d is the effective ionic diameter, and an increasing correlation length with salt concentration in very concentrated electrolytes. As in the experiments, we find that the screening length in the concentrated regime follows a universal scaling law as a function d/λ_D for all studied salts. However, the scaling exponent is significantly lower than the experimentally measured one and lies in the range of the theoretical predictions based on much simpler electrolyte models.

Screening length for finite-size ions in concentrated electrolytes

The classical Debye-Hückel (DH) theory clearly accounts for the origin of screening in electrolyte solutions and works rather well for dilute electrolyte solutions. While the Debye screening length decreases with the ion concentration and is independent of ion size, recent surface-force measurements imply that for concentrated solutions, the screening length exhibits an opposite trend; it increases with ion concentration and depends on the ionic size. The screening length is usually defined by the response of the electrolyte solution to a test charge but can equivalently be derived from the charge-charge correlation function. By going beyond DH theory, we predict the effects of ion size on the charge-charge correlation function. A simple modification of the Coulomb interaction kernel to account for the excluded volume of neighboring ions yields a nonmonotonic dependence of the screening length (correlation length) on the ionic concentration, as well as damped charge oscillations for high concentrations.

Failure of Debye-Hückel Screening in Low-Charge Colloidal Suspensions

Derjaguin–Landau–Verwey–Overbeek (DLVO) theory remains the cornerstone of colloid stability. Electrostatic interactions dominate van der Waals attractions at large colloid-colloid separations h, unless strongly screened. Under these conditions, the potential U ( h ) between charged colloids is expected to be exponentially screened, U ( h ) ∼ exp ( − κ h ) / h , with κ − 1 = λ D where λ D is the classical Debye-Hückel screening length. By measuring the force between individual charged particles at dilute electrolyte concentrations (<mM) using optical tweezers, we tested experimentally the prediction κ − 1 = λ D in a nonpolar solvent. At low salt concentrations, we found close agreement between the directly-measured decay length κ − 1 and Debye-Hückel predictions. However, above a critical electrolyte concentration (≈450 μ M), we obtained significant discrepancies between measured and predicted screening lengths, with κ − 1 ≫ λ D . In marked contrast to expectations, we found that the measured screening length κ − 1 appears to grow as the ionic strength of the solution is increased. The origin of this discrepancy is discussed and the importance of considering the surface is highlighted.

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.

Variational approach for electrolyte solutions: from dielectric interfaces to charged nanopores

A variational theory is developed to study electrolyte solutions, composed of interacting pointlike ions in a solvent, in the presence of dielectric discontinuities and charges at the boundaries. Three important and nonlinear electrostatic effects induced by these interfaces are taken into account: surface charge induced electrostatic field, solvation energies due to the ionic cloud, and image-charge repulsion. Our variational equations thus go beyond the mean-field theory, or weak coupling limit, where thermal fluctuations overcome electrostatic correlations, and allows one to reach the opposite strong coupling limit, where electrostatic interactions induced by interfaces dominate. The influence of salt concentration, ion valency, dielectric jumps, and surface charge is studied in two geometries. (i) A single neutral dielectric interface (e.g., air-water or electrolyte-membrane) with an asymmetric electrolyte. A charge separation and thus an electrostatic field get established due to the different image-charge repulsions for coions and counterions. Both charge distributions and surface tension are computed and compared to previous approximate calculations. For symmetric electrolyte solutions close to a charged surface, two zones are characterized. In the first one, in contact with the surface and with size proportional to the logarithm of the coupling parameter, strong image forces and strong coupling impose a total ion exclusion, while in the second zone the mean-field approach applies. (ii) A symmetric electrolyte confined between two dielectric interfaces as a simple model of ion rejection from nanopores in membranes. The competition between image-charge repulsion and attraction of counterions by the membrane charge is studied. For small surface charge, the counterion partition coefficient decreases with increasing pore size up to a critical pore size, contrary to neutral membranes. For larger pore sizes, the whole system behaves like a neutral pore. For strong coupling and small pore size, coion exclusion is total and the counterion partition coefficient is solely determined by global electroneutrality. A quantitative comparison is made with a previous approach, where image and surface charge effects were smeared out in the pore. It is shown that the variational method allows one to go beyond the constant Donnan potential approximation, with deviations stronger at high ion concentrations or small pore sizes. The prediction of the variational method is also compared with MC simulations and good agreement is observed.