Interaction regimes for oppositely charged plates with multivalent counterions

Strong-coupling theory of counterions with hard cores between symmetrically charged walls

By a combination of Monte Carlo simulations and analytical calculations, we investigate the effective interactions between highly charged planar interfaces, neutralized by mobile counterions (salt-free system). While most previous analysis have focused on pointlike counterions, we treat them as charged hard spheres. We thus work out the fate of like-charge attraction when steric effects are at work. The analytical approach partitions counterions in two subpopulations, one for each plate, and integrates out one subpopulation to derive an effective Hamiltonian for the remaining one. The effective Hamiltonian features plaquette four-particle interactions, and it is worked out by computing a Gibbs-Bogoliubov inequality for the free energy. At the root of the treatment is the fact that under strong electrostatic coupling, the system of charges forms an ordered arrangement, that can be affected by steric interactions. Fluctuations around the reference positions are accounted for. To dominant order at high coupling, it is found that steric effects do not significantly affect the interplate effective pressure, apart at small distances where hard-sphere overlap are unavoidable, and thus rule out configurations.

Ionic size effects on the Poisson-Boltzmann theory

In this paper, we develop a simple theory to study the effects of ionic size on ionic distributions around a charged spherical particle. We include a correction to the regular Poisson-Boltzmann equation in order to take into account the size of ions in a mean-field regime. The results are compared with Monte Carlo simulations and a density functional theory based on the fundamental measure approach and a second-order bulk expansion which accounts for electrostatic correlations. The agreement is very good even for multivalent ions. Our results show that the theory can be applied with very good accuracy in the description of ions with highly effective ionic radii and low concentration, interacting with a colloid or a nanoparticle in an electrolyte solution.

Strong-coupling theory of counterions between symmetrically charged walls: from crystal to fluid phases

We study thermal equilibrium of classical pointlike counterions confined between symmetrically charged walls at distance d. At very large couplings when the counterion system is in its crystal phase, a harmonic expansion of particle deviations is made around the bilayer positions, with a free lattice parameter determined from a variational approach. For each of the two walls, the harmonic expansion implies an effective one-body potential at the root of all observables of interest in our Wigner strong-coupling expansion. Analytical results for the particle density profile and the pressure are in good agreement with numerical Monte Carlo data, for small as well as intermediate values of d comparable with the Wigner lattice spacing. While the strong-coupling theory is extended to the fluid regime by using the concept of a correlation hole, the Wigner calculations appear trustworthy for all electrostatic couplings investigated. Our results significantly extend the range of accuracy of analytical equations of state for strongly interacting charged planar interfaces.

Mean-field beyond mean-field: the single particle view for moderately to strongly coupled charged fluids

In a counter-ion only charged fluid, Coulomb coupling is quantified by a single dimensionless parameter. Yet, the theoretical treatment of moderately to strongly coupled charged fluids is a difficult task, central to the understanding of a wealth of soft matter problems, including biological systems. We show that the corresponding coupling regime can be remarkably well described by a single particle treatment, which, at variance with previous works, takes due account of inter-ionic interactions. To this end, the prototypical problem of a planar charged dielectric interface is worked out. Testing our predictions against Monte Carlo simulation data reveals an excellent agreement.

Charge neutrality breakdown in confined aqueous electrolytes: Theory and simulation

We study, using Density Functional theory (DFT) and Monte Carlo simulations, aqueous electrolyte solutions between charged infinite planar surfaces, in contact with a bulk salt reservoir. In agreement with recent experimental observations [Z. Luo et al., Nat. Commun. 6, 6358 (2015)], we find that the confined electrolyte lacks local charge neutrality. We show that a DFT based on a bulk-HNC expansion properly accounts for strong electrostatic correlations and allows us to accurately calculate the ionic density profiles between the charged surfaces, even for electrolytes containing trivalent counterions. The DFT allows us to explore the degree of local charge neutrality violation, as a function of plate separation and bulk electrolyte concentration, and to accurately calculate the interaction force between the charged surfaces.

Screening like charges in one-dimensional Coulomb systems: Exact results

The possibility that like charges can attract each other under the mediation of mobile counterions is by now well documented experimentally, numerically, and analytically. Yet, obtaining exact results is in general impossible, or restricted to some limiting cases. We work out here in detail a one-dimensional model that retains the essence of the phenomena present in higher-dimensional systems. The partition function is obtained explicitly, from which a wealth of relevant quantities follow, such as the effective force between the charges or the counterion profile in their vicinity. Isobaric and canonical ensembles are distinguished. The case of two equal charges screened by an arbitrary number N of counterions is first studied, before the more general asymmetric situation is addressed. It is shown that the parity of N plays a key role in the long-range physics.

Numerical study of density functional theory with mean spherical approximation for ionic condensation in highly charged confined electrolytes

We investigate numerically a density functional theory (DFT) for strongly confined ionic solutions in the canonical ensemble by comparing predictions of ionic concentration profiles and pressure for the double-layer configuration to those obtained with Monte Carlo (MC) simulations and the simpler Poisson-Boltzmann (PB) approach. The DFT consists of a bulk (ion-ion) and an ion-solid part. The bulk part includes nonideal terms accounting for long-range electrostatic and short-range steric correlations between ions and is evaluated with the mean spherical approximation and the local density approximation. The ion-solid part treats the ion-solid interactions at the mean-field level through the solution of a Poisson problem. The main findings are that ionic concentration profiles are generally better described by PB than by DFT, although DFT captures the nonmonotone co-ion profile missed by PB. Instead, DFT yields more accurate pressure predictions than PB, showing in particular that nonideal effects are important to describe highly confined ionic solutions. Finally, we present a numerical methodology capable of handling nonconvex minimization problems so as to explore DFT predictions when the reduced temperature falls below the critical temperature.

Wigner-crystal formulation of strong-coupling theory for counterions near planar charged interfaces

We present a new analytical approach to the strong electrostatic coupling regime (SC) that can be achieved equivalently at low temperatures, high charges, low dielectric permittivity, etc. Two geometries are analyzed in detail: one charged wall first, and then two parallel walls at small distances that can be likely or oppositely charged. In all cases, only one type of mobile counterions is present, and ensures electroneutrality (salt-free case). The method is based on a systematic expansion around the ground state formed by the two-dimensional Wigner crystal(s) of counterions at the plate(s). The leading SC order stems from a single-particle theory, and coincides with the virial SC approach that has been much studied in the last 10 years. The first correction has the functional form of the virial SC prediction, but the prefactor is different. The present theory is free of divergences and the obtained results, both for symmetrically and asymmetrically charged plates, are in excellent agreement with available data of Monte Carlo simulations under strong and intermediate Coulombic couplings. All results obtained represent relevant improvements over the virial SC estimates. The present SC theory starting from the Wigner crystal and therefore coined Wigner SC, sheds light on anomalous phenomena like the counterion mediated like-charge attraction, and the opposite-charge repulsion.