The role of multipoles in counterion-mediated interactions between charged surfaces: strong and weak coupling

Impact of the inner solute structure on the electrostatic mean-field and strong-coupling regimes of macromolecular interactions

The structural diversity of the solute molecules involved in biomolecular processes necessitates the characterization of the forces between charged macromolecules beyond the point-ion description. From the field-theoretic partition function of an electrolyte confined between two anionic membranes, we derive a contact-value identity valid for general intramolecular solute structure and electrostatic coupling strength. In the electrostatic mean-field regime, the inner charge spread of the solute particles is shown to induce the twofold enhancement of the short-range Poisson-Boltzmann level membrane repulsion and a longer-range depletion attraction. Our contact theorem indicates that the twofold repulsion enhancement by solute size is equally present in the opposite strong-coupling regime of linear and spherical solute molecules. Upon the inclusion of the dielectric contrast between the electrolyte and the interacting membranes, the emerging polarization forces substantially amplify the solute specificity of the macromolecular interactions. Namely, the finite size of the dumbbell-like solute particles composed of similar terminal charges weakens the intermembrane repulsion. However, the extended structure of the solute molecules carrying opposite elementary charges such as ionized atoms and zwitterionic molecules enhances the membrane repulsion by several factors. We also show that these polarization forces can extend the range of the solute structure effects up to intermembrane distances exceeding the solute size by an order of magnitude. This radical alteration of the intermembrane interactions by the salt structure identifies the solute specificity as a key ingredient of the thermodynamic stability in colloidal systems.

Interactions between Rigid Polyelectrolytes Mediated by Ordering and Orientation of Multivalent Nonspherical Ions in Salt Solutions

Multivalent ions in solutions with polyelectrolytes (PEs) induce electrostatic correlations that can drastically change ion distributions around the PEs and their mutual interactions. Using coarse-grained molecular dynamics simulations, we show how in addition to valency, ion shape and concentration can be harnessed as tools to control rigid like-charged PE-PE interactions. We demonstrate a correlation between the orientational ordering of aspherical ions and how they mediate the effective PE-PE attraction induced by multivalency. The interaction type, strength, and range can thus be externally controlled in ionic solutions. Our results can be used as generic guidelines to tune the self-assembly of like-charged polyelectrolytes by variation of the characteristics of the ions.

Theory of Inhomogeneous Rod-like Coulomb Fluids

A field theoretic representation of the classical partition function is derived for a system composed of a mixture of anisotropic and isotropic mobile charges that interact via long range Coulomb and short range nematic interactions. The field theory is then solved on a saddle-point approximation level, leading to a coupled system of Poisson–Boltzmann and Maier–Saupe equations. Explicit solutions are finally obtained for a rod-like counterion-only system in proximity to a charged planar wall. The nematic order parameter profile, the counterion density profile and the electrostatic potential profile are interpreted within the framework of a nematic–isotropic wetting phase with a Donnan potential difference.

Like-charge polymer-membrane complexation mediated by multivalent cations: One-loop-dressed strong coupling theory

We probe the electrostatic mechanism driving adsorption of polyelectrolytes onto like-charged membranes upon the addition of tri- and tetravalent counterions to a bathing monovalent salt solution. We develop a one-loop-dressed strong coupling theory that treats the monovalent salt at the electrostatic one-loop level and the multivalent counterions within a strong-coupling approach. It is shown that the adhesive force of the multivalent counterions mediating the like-charge adsorption arises from their strong condensation at the charged membrane. The resulting interfacial counterion excess locally maximizes the screening ability of the electrolyte and minimizes the electrostatic polymer grand potential. This translates into an attractive force that pulls the polymer to the similarly charged membrane. We show that the high counterion valency enables this adsorption transition even at weakly charged membranes. Additionally, strongly charged membranes give rise to monovalent counterion-induced correlations and intensify the interfacial multivalent counterion condensation, strengthening the complexation of the polymer with the like-charged membrane, as well as triggering the orientational transition of the molecule prior to its adsorption. Finally, our theory provides two additional key features as evidenced by previous adsorption experiments: first, the critical counterion concentration for polymer adsorption decreases with the rise of the counterion valency and, second, the addition of monovalent salt enhances the screening of the membrane charges and suppresses monovalent counterion correlations close to the surface. This weakens the interfacial multivalent counterion condensation and results in the desorption of the polymer from the substrate.

Orientational transition and complexation of DNA with anionic membranes: Weak and intermediate electrostatic coupling

We characterize the role of charge correlations in the adsorption of a short, rodlike anionic polyelectrolyte onto a similarly charged membrane. Our theory reveals two different mechanisms driving the like-charge polyelectrolyte-membrane complexation: In weakly charged membranes, repulsive polyelectrolyte-membrane interactions lead to the interfacial depletion and a parallel orientation of the polyelectrolyte with respect to the membrane; while in the intermediate membrane charge regime, the interfacial counterion excess gives rise to an attractive "salt-induced" image force. This furthermore results in an orientational transition from a parallel to a perpendicular configuration and a subsequent short-ranged like-charge adsorption of the polyelectrolyte to the substrate. A further increase of the membrane charge engenders a charge inversion, originating from surface-induced ionic correlations, that act as a separate mechanism capable of triggering the like-charge polyelectrolyte-membrane complexation over an extended distance interval from the membrane surface. The emerging picture of this complexation phenomenon identifies the interfacial "salt-induced" image forces as a powerful control mechanism in polyelectrolyte-membrane complexation.

Role of metallic core for the stability of virus-like particles in strongly coupled electrostatics

Electrostatic interactions play important roles in the formation and stability of viruses and virus-like particles (VLPs) through processes that often involve added, or naturally occurring, multivalent ions. Here, we investigate the electrostatic or osmotic pressure acting on the proteinaceous shell of a generic model of VLPs, comprising a charged outer shell and a metallic nanoparticle core, coated by a charged layer and bathed in an aqueous electrolyte solution. Motivated by the recent studies accentuating the role of multivalent ions for the stability of VLPs, we focus on the effects of multivalent cations and anions in an otherwise monovalent ionic solution. We perform extensive Monte-Carlo simulations based on appropriate Coulombic interactions that consistently take into account the effects of salt screening, the dielectric polarization of the metallic core, and the strong-coupling electrostatics due to multivalent ions. We specifically study the intricate roles these factors play in the electrostatic stability of the model VLPs. It is shown that while the insertion of a metallic nanoparticle by itself can produce negative, inward-directed, pressure on the outer shell, addition of only a small amount of multivalent counterions can robustly engender negative pressures, enhancing the VLP stability across a wide range of values for the system parameters.

Charged nanorods at heterogeneously charged surfaces

We study the spatial and orientational distribution of charged nanorods (rodlike counterions) as well as the effective interaction mediated by them between two plane-parallel surfaces that carry fixed (quenched) heterogeneous charge distributions. The nanorods are assumed to have an internal charge distribution, specified by a multivalent monopolar moment and a finite quadrupolar moment, and the quenched surface charge is assumed to be randomly distributed with equal mean and variance on the two surfaces. While equally charged surfaces are known to repel within the traditional mean-field theories, the presence of multivalent counterions has been shown to cause attractive interactions between uniformly charged surfaces due to the prevalence of strong electrostatic couplings that grow rapidly with the counterion valency. We show that the combined effects due to electrostatic correlations (caused by the coupling between the mean surface field and the multivalent, monopolar, charge valency of counterions) as well as the disorder-induced interactions (caused by the coupling between the surface disorder field and the quadrupolar moment of counterions) lead to much stronger attractive interactions between two randomly charged surfaces. The interaction profile turns out to be a nonmonotonic function of the intersurface separation, displaying an attractive minimum at relatively small separations, where the ensuing attraction can exceed the maximum strong-coupling attraction (produced by multivalent monopolar counterions between uniformly charged surfaces) by more than an order of magnitude.

Electrostatics of Nanoparticle-Wall Interactions within Nanochannels: Role of Double-Layer Structure and Ion-Ion Correlations

We perform computational investigations of the electrolyte-mediated interactions of charged nanoparticles with the walls of nanochannels. We investigate the role of discrete ion effects, valence, and electrolyte strength on nanoparticle-wall interactions. We find for some of the multivalent charge regimes that the like-charged nanoparticles and walls can have attractive interactions. We study in detail these interactions and the free-energy profile for the nanoparticle-wall separation. We find there are energy barriers and energy minima giving preferred nanoparticle locations in the channel near the center and at a distance near to but separated from the channel walls. We characterize contributions from surface overcharging, condensed layers, and overlap of ion double layers. We perform our investigations using coarse-grained particle-level simulations with Brownian dynamics, classical density functional theory, and the mean-field Poisson-Boltzmann theory. We discuss the implications of our results for phenomena in nanoscale devices.

Rodlike Counterions near Charged Cylinders: Counterion Condensation and Intercylinder Interaction

We study a system composed of like-charged cylinders and dumbbell-like counterions, with the focus laid on the role of the internal structure of counterions, using Monte Carlo simulations. The dumbbell ions are found to exhibit novel counterion condensation behavior governed by their length. Effective electrostatic interactions mediated between charged parallel cylinders also turn out significantly different from the case of pointlike ions, as a result of the complex interplay between the spatially separated charge distribution in the dumbbell counterions, their orientation, and the curvature of the charged cylinder. We show that at a weak-to-moderate electrostatic coupling strength, where effective like-charge interactions are usually found to be repulsive, the intercylinder interaction can become attractive and display a distinctive sensitivity to the cylinder curvature and dumbbell size, proving the significant effect of ion structure.

Mean-field theory of active electrolytes: Dynamic adsorption and overscreening

We investigate active electrolytes within the mean-field level of description. The focus is on how the double-layer structure of passive, thermalized charges is affected by active dynamics of constituting ions. One feature of active dynamics is that particles adhere to hard surfaces, regardless of chemical properties of a surface and specifically in complete absence of any chemisorption or physisorption. To carry out the mean-field analysis of the system that is out of equilibrium, we develop the "mean-field simulation" technique, where the simulated system consists of charged parallel sheets moving on a line and obeying active dynamics, with the interaction strength rescaled by the number of sheets. The mean-field limit becomes exact in the limit of an infinite number of movable sheets.

Adsorption of highly charged Gaussian polyelectrolytes onto oppositely charged surfaces

In many biological processes highly charged biopolymers are adsorbed onto oppositely charged surfaces of macroions and membranes. They form strongly correlated structures close to the surface which cannot be explained by the conventional Poisson-Boltzmann theory. In this work strong coupling theory is used to study the adsorption of highly charged Gaussian polyelectrolytes. Two cases of adsorptions are considered, when the Gaussian polyelectrolytes are confined (a) by one charged wall, and (b) between two charged walls. The effects of salt and the geometry of the polymers on their adsorption-depletion transitions in the strong coupling regime are discussed.

Strong-coupling electrostatic theory of polymer counterions close to planar charges

Strong-coupling phenomena, such as like-charge macroion attraction, opposite-charged macroion repulsion, charge renormalization, and charge inversion, are known to be mediated by multivalent counterions. Most theories treat the counterions as point charges and describe the system by a single coupling parameter that measures the strength of the Coulomb interactions. In many biological systems, the counterions are highly charged and have finite sizes and can be well-described by polyelectrolytes. The shapes and orientations of these polymer counterions play a major role in the thermodynamics of these systems. In this work we apply a field-theoretic description in the strong-coupling regime to the polymer counterions in the presence of a fixed charge distribution. We work out the special cases of rodlike polymer counterions confined by one, and two charged walls, respectively. The effects of the geometry of the rodlike counterions and the excluded volume of the walls on the density, pressure, and free energy of the rodlike counterions are discussed.

Multiscale Multiphysics and Multidomain Models I: Basic Theory

This work extends our earlier two-domain formulation of a differential geometry based multiscale paradigm into a multidomain theory, which endows us the ability to simultaneously accommodate multiphysical descriptions of aqueous chemical, physical and biological systems, such as fuel cells, solar cells, nanofluidics, ion channels, viruses, RNA polymerases, molecular motors and large macromolecular complexes. The essential idea is to make use of the differential geometry theory of surfaces as a natural means to geometrically separate the macroscopic domain of solvent from the microscopic domain of solute, and dynamically couple continuum and discrete descriptions. Our main strategy is to construct energy functionals to put on an equal footing of multiphysics, including polar (i.e., electrostatic) solvation, nonpolar solvation, chemical potential, quantum mechanics, fluid mechanics, molecular mechanics, coarse grained dynamics and elastic dynamics. The variational principle is applied to the energy functionals to derive desirable governing equations, such as multidomain Laplace-Beltrami (LB) equations for macromolecular morphologies, multidomain Poisson-Boltzmann (PB) equation or Poisson equation for electrostatic potential, generalized Nernst-Planck (NP) equations for the dynamics of charged solvent species, generalized Navier-Stokes (NS) equation for fluid dynamics, generalized Newton's equations for molecular dynamics (MD) or coarse-grained dynamics and equation of motion for elastic dynamics. Unlike the classical PB equation, our PB equation is an integral-differential equation due to solvent-solute interactions. To illustrate the proposed formalism, we have explicitly constructed three models, a multidomain solvation model, a multidomain charge transport model and a multidomain chemo-electro-fluid-MD-elastic model. Each solute domain is equipped with distinct surface tension, pressure, dielectric function, and charge density distribution. In addition to long-range Coulombic interactions, various non-electrostatic solvent-solute interactions are considered in the present modeling. We demonstrate the consistency between the non-equilibrium charge transport model and the equilibrium solvation model by showing the systematical reduction of the former to the latter at equilibrium. This paper also offers a brief review of the field.

Microscopic formulation of nonlocal electrostatics in polar liquids embedding polarizable ions

Nonlocal electrostatic interactions associated with finite solvent size and ion polarizability are investigated within the mean-field linear response theory. To this end, we introduce a field-theoretic model of a polar liquid composed of linear multipole solvent molecules and embedding polarizable ions modeled as Drude oscillators. Unlike previous dipolar Poisson-Boltzmann formulations treating the solvent molecules as point dipoles, our model is able to qualitatively reproduce the non-local dielectric response behavior of polar liquids observed in molecular dynamics simulations and atomic force microscope experiments for water solvent at charged interfaces. The present theory explains the formation of the associated interfacial hydration layers in terms of a cooperative dipolar response mechanism driven by the reaction of the solvent molecules to their own polarization field. We also incorporate into the theory the relative multipole moments of water molecules obtained from quantum mechanical calculations and show that the multipolar contributions to the dielectric permittivity are largely dominated by the dipolar one. We find that this stems from the mutual cancellation of the first two interfacial hydration layers of opposite net charge for multipolar liquids. Within the same nonlocal dielectric response theory, we show that the induced ion polarizability reverses the interfacial ion density trends predicted by the Poisson-Boltzmann theory, resulting in a surface affinity of coions and exclusion of counterions. The results indicate that the consideration of the discrete charge composition of solvent molecules and ions is the key step towards a microscopic understanding of nonlocal electrostatic effects in polar solvents.

Electrostatic interactions mediated by polarizable counterions: weak and strong coupling limits

We investigate the statistical mechanics of an inhomogeneous Coulomb fluid composed of charged particles with static polarizability. We derive the weak- and the strong-coupling approximations and evaluate the partition function in a planar dielectric slab geometry with charged boundaries. We investigate the density profiles and the disjoining pressure for both approximations. Comparison to the case of non-polarizable counterions shows that polarizability brings important differences in the counterion density distribution as well as the counterion mediated electrostatic interactions between charged dielectric interfaces.

Novel kinetic trapping in charged colloidal clusters due to self-induced surface charge organization

Colloidal clusters are an unusual state of matter where tunable interactions enable a sufficient reduction in their degrees of freedom that their energy landscapes can become tractable - they form a playground for statistical mechanics and promise unprecedented control of structure on the submicron lengthscale. We study colloidal clusters in a system where a short-ranged polymer-induced attraction drives clustering, while a weak, long-ranged electrostatic repulsion prevents extensive aggregation. We compare experimental yields of cluster structures with theory which assumes simple addition of competing isotropic interactions between the colloids. Here we show that for clusters of size 4 ≤ m ≤ 7, the yield of minimum energy clusters is much less than expected. We attribute this to an anisotropic self-organized surface charge distribution which leads to unexpected kinetic trapping. We introduce a model for the coupling between counterions and binding sites on the colloid surface with which we interpret our findings.

The one-dimensional Coulomb lattice fluid capacitor

The one-dimensional Coulomb lattice fluid in a capacitor configuration is studied. The model is formally exactly soluble via a transfer operator method within a field theoretic representation of the model. The only interactions present in the model are the one-dimensional Coulomb interaction between cations and anions and the steric interaction imposed by restricting the maximal occupancy at any lattice site to one particle. Despite the simplicity of the model, a wide range of intriguing physical phenomena arise, some of which are strongly reminiscent of those seen in experiments and numerical simulations of three-dimensional ionic liquid based capacitors. Notably, we find regimes where over-screening and density oscillations are seen near the capacitor plates. The capacitance is also shown to exhibit strong oscillations as a function of applied voltage. It is also shown that the corresponding mean-field theory misses most of these effects. The analytical results are confirmed by extensive numerical simulations.

Nonlinear Poisson equation for heterogeneous media

The Poisson equation is a widely accepted model for electrostatic analysis. However, the Poisson equation is derived based on electric polarizations in a linear, isotropic, and homogeneous dielectric medium. This article introduces a nonlinear Poisson equation to take into consideration of hyperpolarization effects due to intensive charges and possible nonlinear, anisotropic, and heterogeneous media. Variational principle is utilized to derive the nonlinear Poisson model from an electrostatic energy functional. To apply the proposed nonlinear Poisson equation for the solvation analysis, we also construct a nonpolar solvation energy functional based on the nonlinear Poisson equation by using the geometric measure theory. At a fixed temperature, the proposed nonlinear Poisson theory is extensively validated by the electrostatic analysis of the Kirkwood model and a set of 20 proteins, and the solvation analysis of a set of 17 small molecules whose experimental measurements are also available for a comparison. Moreover, the nonlinear Poisson equation is further applied to the solvation analysis of 21 compounds at different temperatures. Numerical results are compared to theoretical prediction, experimental measurements, and those obtained from other theoretical methods in the literature. A good agreement between our results and experimental data as well as theoretical results suggests that the proposed nonlinear Poisson model is a potentially useful model for electrostatic analysis involving hyperpolarization effects.

Variational multiscale models for charge transport

This work presents a few variational multiscale models for charge transport in complex physical, chemical and biological systems and engineering devices, such as fuel cells, solar cells, battery cells, nanofluidics, transistors and ion channels. An essential ingredient of the present models, introduced in an earlier paper (Bulletin of Mathematical Biology, 72, 1562-1622, 2010), is the use of differential geometry theory of surfaces as a natural means to geometrically separate the macroscopic domain from the microscopic domain, meanwhile, dynamically couple discrete and continuum descriptions. Our main strategy is to construct the total energy functional of a charge transport system to encompass the polar and nonpolar free energies of solvation, and chemical potential related energy. By using the Euler-Lagrange variation, coupled Laplace-Beltrami and Poisson-Nernst-Planck (LB-PNP) equations are derived. The solution of the LB-PNP equations leads to the minimization of the total free energy, and explicit profiles of electrostatic potential and densities of charge species. To further reduce the computational complexity, the Boltzmann distribution obtained from the Poisson-Boltzmann (PB) equation is utilized to represent the densities of certain charge species so as to avoid the computationally expensive solution of some Nernst-Planck (NP) equations. Consequently, the coupled Laplace-Beltrami and Poisson-Boltzmann-Nernst-Planck (LB-PBNP) equations are proposed for charge transport in heterogeneous systems. A major emphasis of the present formulation is the consistency between equilibrium LB-PB theory and non-equilibrium LB-PNP theory at equilibrium. Another major emphasis is the capability of the reduced LB-PBNP model to fully recover the prediction of the LB-PNP model at non-equilibrium settings. To account for the fluid impact on the charge transport, we derive coupled Laplace-Beltrami, Poisson-Nernst-Planck and Navier-Stokes equations from the variational principle for chemo-electro-fluid systems. A number of computational algorithms is developed to implement the proposed new variational multiscale models in an efficient manner. A set of ten protein molecules and a realistic ion channel, Gramicidin A, are employed to confirm the consistency and verify the capability. Extensive numerical experiment is designed to validate the proposed variational multiscale models. A good quantitative agreement between our model prediction and the experimental measurement of current-voltage curves is observed for the Gramicidin A channel transport. This paper also provides a brief review of the field.

Polarizable Poisson-Boltzmann equation: the study of polarizability effects on the structure of a double layer

We incorporate ion polarizabilities into the Poisson-Boltzmann equation by modifying the effective dielectric constant and the Boltzmann distribution of ions. The extent of the polarizability effects is controlled by two parameters, γ(1) and γ(2); γ(1) determines the polarization effects in a dilute system and γ(2) regulates the dependence of the polarizability effects on the concentration of ions. For a polarizable ion in an aqueous solution γ(1) ≈ 0.01 and the polarizability effects are negligible. The conditions where γ(1) and/or γ(2) are large and the polarizability is relevant involve the low dielectric constant media, high surface charge, and/or large ionic concentrations.

Protein-DNA interactions determine the shapes of DNA toroids condensed in virus capsids

DNA toroids that form inside the bacteriophage capsid present different shapes according to whether they are formed by the addition of spermine or polyethylene glycol to the bathing solution. Spermine-DNA toroids present a convex, faceted section with no or minor distortions of the DNA interstrand spacing with respect to those observed in the bulk, whereas polyethylene glycol-induced toroids are flattened to the capsid inner surface and show a crescent-like, nonconvex shape. By modeling the energetics of the DNA toroid using a free-energy functional composed of energy contributions related to the elasticity of the wound DNA, exposed surface DNA energy, and adhesion between the DNA and the capsid, we established that the crescent shape of the toroidal DNA section comes from attractive interactions between DNA and the capsid. Such attractive interactions seem to be specific to the PEG condensation process and are not observed in the case of spermine-induced DNA condensation.

A submicron device to rectify a square-wave angular velocity

We study a system composed of two thick dielectric disks separated by a thin layer of an electrolyte solution. Initially both plates have the same surface charge distribution. The surface charge distribution has no rotational symmetry. We show that the top plate experiences a torque [Formula: see text]([Formula: see text]) if it rotates about its axis by an angle [Formula: see text] . The torque can be controlled by varying the electrolyte concentration, the separation and the surface charge density of the plates. For a specific example of charged rods attached to the plates, we find [Formula: see text]([Formula: see text]) [Formula: see text] sin(4[Formula: see text]) . We also study the dynamics of the system. We consider the case where the angular velocity of the bottom disk is a square-wave signal. We find that the average angular velocity of the top disk is not zero.

Charge renormalization of helical macromolecules

Some time ago a theory of electrostatic interaction between helical macromolecules was proposed (Kornyshev and Leikin 1997 J. Chem. Phys. 107 3656): the Kornyshev-Leikin (KL) theory. We place this theory on a more rigorous statistical mechanical grounding, starting from the free energy that can be derived from a grand partition function. We see that the long range behaviour of the force is indeed given by the KL theory, no matter whether the distributions of 'condensed' ionic charge are at the surface of the macromolecule or extend away from it. Thus, for the limiting behaviour, we need only self-consistently calculate the distribution of the condensed fraction of ions for a single macro-ion. This distribution can be related back to interaction parameters: KL parameters. Furthermore, we are able to see within the formalism where corrections due to the hard core radius of the ion enter. For the adjustment of the 'condensed' ions, we show an expression for the leading order contribution, as well as relevant decay lengths. As a demonstration of the theoretical 'machinery', as well as a study of qualitative effects, we calculate the KL parameters in one instance. We use a DNA-like surface charge distribution, where a fraction of the ions are assumed to be bound in the grooves at the surface of a DNA molecule, whereas the rest of the charge distribution is calculated self-consistently. Also, the electrostatic contribution to the counter-ion binding potentials that ions experience within the grooves can be calculated.

Dressed counterions: strong electrostatic coupling in the presence of salt

We reformulate the theory of strong electrostatic coupling in order to describe an asymmetric electrolyte solution of monovalent salt ions and polyvalent counterions using field-theoretical techniques and Monte Carlo simulations. The theory is based on an asymmetric treatment of the different components of the electrolyte solution. The weak coupling Debye-Hückel approach is used in order to describe the monovalent salt ions while a strong coupling approach is used to tackle the polyvalent counterions. This combined weak-strong coupling approach effectively leads to dressed interactions between polyvalent counterions and thus directly affects the correlation attraction mediated by polyvalent counterions between like-charged objects. The general theory is specifically applied to a system composed of two uniformly charged plane-parallel surfaces in the presence of salt and polyvalent counterions. In the strong coupling limit for polyvalent counterions, the comparison with Monte Carlo simulations shows good agreement for large enough values of the electrostatic coupling parameter. We delineate two limiting laws that in fact encompass all the Monte Carlo data.

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.