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

Image charge effects under metal and dielectric boundary conditions

The image charge (IC) effect is a fundamental problem in electrostatics. However, proper treatment at the continuum level for many-ion systems, such as electrolyte solutions or ionic liquids, remains an open theoretical question. Here, we demonstrate and systematically compare the IC effects under metal and dielectric boundary conditions (BCs), based on a renormalized Gaussian-fluctuation theory. Our calculations for a simple 1:1 symmetric electrolyte in the point-charge approximation show that the double-layer structure, capacitance, and interaction forces between like-charged plates depend strongly on the types of boundaries, even in the weak-coupling regime. Like-charge attraction is predicted for both metal and dielectric BCs. Finally, we comment on the effects of a dielectrically saturated solvent layer on the metal surface. We provide these results to serve as a baseline for comparison with more realistic molecular dynamics simulations and experiments.

Strong-coupling effective-field theory for asymmetrically charged plates with counterions only

We are interested in rationalizing the phenomenon of like-charge attraction between charged bodies, such as a pair of colloids, in the strong coupling regime. The two colloids are modelled as uniformly charged parallel plates, neutralized by mobile counterions. In an earlier work [Palaia et al., J. Phys. Chem. B 126, 3143 (2022)1520-610610.1021/acs.jpcb.2c00028], we developed an effective-field theory for symmetric plates, stemming from the ground-state description that holds at infinite couplings. Here, we generalize the approach to the asymmetric case, where the plates bear charges of the same sign, but of different values. In the symmetric situation, the mobile ions, which are localized in the vicinity of the two plates, share equally between both of them. Here, the sharing is nontrivial, depending both on the coupling parameter and the distance between the plates. We thus introduce a counterion occupation parameter that is determined variationally to ensure minimum of the free energy. The analytical results for the pressure as a function of the plate-plate distance d agree well with our Monte Carlo data, in a large interval of strong and intermediate coupling constants Ξ. We show in particular that within this description there exists a range of large distances at which the attractive pressure features a 1/d^{2} behavior.

Electroosmotic Flow Grows with Electrostatic Coupling in Confining Charged Dielectric Surfaces

In this work, the effects of polarization of confining charged planar dielectric surfaces on induced electroosmotic flow are investigated. To this end, we use dissipative particle dynamics to model solvent and ionic particles together with a modified Ewald sum method to model electrostatic interactions and surfaces polarization. A relevant difference between counterions number density profiles, velocity profiles, and volumetric flow rates obtained with and without surface polarization for moderate and high electrostatic coupling parameters is observed. For low coupling parameters, the effect is negligible. An increase of almost 500% in volumetric flow rate for moderate/high electrostatic coupling and surface separation is found when polarizable surfaces are considered. The most important result is that the increase in electrostatic coupling substantially increases the electroosmotic flow in all studied range of separations when the dielectric constant of electrolytes is much higher than the dielectric constant of confining walls. For the higher separation simulated, an increase of around 340% in volumetric flow rate when the electrostatic coupling is increased by a factor of two orders of magnitude is obtained.

Imaging Arrangements of Discrete Ions at Liquid-Solid Interfaces

The individual and collective behavior of ions near electrically charged interfaces is foundational to a variety of electrochemical phenomena encountered in biology, energy, and the environment. While many theories have been developed to predict the interfacial arrangements of counterions, direct experimental observations and validations have remained elusive. Utilizing cryo-electron microscopy, here we directly visualize individual counterions and reveal their discrete interfacial layering. Comparison with simulations suggests the strong effects of finite ionic size and electrostatic interactions. We also uncover correlated ionic structures under extreme confinement, with the channel widths approaching the ionic diameter (∼1 nm). Our work reveals the roles of ionic size, valency, and confinement in determining the structures of liquid-solid interfaces and opens up new opportunities to study such systems at the single-ion level.

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.

Multivalent counterions induced attraction between DNA polyelectrolytes

In this paper we study the electrostatic attraction between two parallel rodlike DNA polyelectrolytes induced by neutralizing multivalent counterions at the zero temperature limit. The counterions crystallize on the charged surfaces of DNA so that we can handle the system by using the Wigner crystal lattice model. We derived the 3D ground state configuration of counterions with minimized energy by use of the gradient descent method, and calculated the interaction between two DNA cylinders with divalent or trivalent counterions when they approach. The results show that the complex ground state configuration of counterions plays a key role in determining the caused attraction. The counterions form three-dimensional Wigner crystals on each cylinder at large separation. When the cylinders are brought together, some counterion lines will move towards the inner region and lead to strong attraction. The calculated interaction from our model is in good agreement with the simulation result, however, the single particle approximation considerably overestimates the attraction.

Counterions in the ground state partially condense in the intervening region of two approaching DNA cylinders and result in attraction.

Electric double layers with surface charge modulations: Exact Poisson-Boltzmann solutions

Poisson-Boltzmann theory is the cornerstone for soft matter electrostatics. We provide exact analytical solutions to this nonlinear mean-field approach for the diffuse layer of ions in the vicinity of a planar or a cylindrical macroion. While previously known solutions are for homogeneously charged objects, the cases worked out exhibit a modulated surface charge-or equivalently, surface potential-on the macroion (wall) surface. In addition to asymptotic features at large distances from the wall, attention is paid to the fate of the contact theorem, relating the contact density of ions to the local wall charge density. For salt-free systems (counterions only), we make use of results pertaining to the two-dimensional Liouville equation, supplemented by an inverse approach. When salt is present, we invoke the exact two-soliton solution to the 2D sinh-Gordon equation. This leads to inhomogeneous charge patterns, that are either localized or periodic in space. Without salt, the electrostatic signature of a charge pattern on the macroion fades exponentially with distance for a planar macroion, while it decays as an inverse power law for a cylindrical macroion. With salt, our study is limited to the planar geometry and reveals that pattern screening is exponential.

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.

A review of immune amplification via ligand clustering by self-assembled liquid-crystalline DNA complexes

We examine how the interferon production of plasmacytoid dendritic cells is amplified by the self-assembly of liquid-crystalline antimicrobial peptide/DNA complexes. These specialized dendritic cells are important for host defense because they quickly release large quantities of type I interferons in response to infection. However, their aberrant activation is also correlated with autoimmune diseases such as psoriasis and lupus. In this review, we will describe how polyelectrolyte self-assembly and the statistical mechanics of multivalent interactions contribute to this process. In a more general compass, we provide an interesting conceptual corrective to the common notion in molecular biology of a dichotomy between specific interactions and non-specific interactions, and show examples where one can construct exquisitely specific interactions using non-specific interactions.

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.

Phase diagram and critical properties of Yukawa bilayers

We study the ground-state Wigner bilayers of pointlike particles with Yukawa pairwise interactions, confined to the surface of two parallel hard walls at dimensionless distance η. The model involves as limiting cases the unscreened Coulomb potential and hard spheres. The phase diagram of Yukawa particles, studied numerically by Messina and Löwen [Phys. Rev. Lett. 91, 146101 (2003)], exhibits five different staggered phases as η varies from 0 to intermediate values. We present a lattice summation method using the generalized Misra functions which permits us to calculate the energy per particle of the phases with a precision much higher than usual in computer simulations. This allows us to address some tiny details of the phase diagram. Going from the hexagonal phase I to phase II is shown to occur at η=0. All second-order phase transitions are proved to be of mean-field type. We also derive the asymptotic shape of critical lines close to the Coulomb and hard-spheres limits. In and close to the hard-spheres limit, the dependence of the internal parameters of the present phases on η is determined exactly.

Counterion density profile around a charged disk: From the weak to the strong association regime

We present a comprehensive study of the two-dimensional one-component plasma in the cell model with charged boundaries. Starting from weak couplings through a convenient approximation of the interacting potential we were able to obtain an analytic formulation to the problem deriving the partition function, density profile, contact densities, and integrated profiles that compared well with the numerical data from Monte Carlo simulations. Additionally, we derived the exact solution for the special cases of Ξ=1,2,3,⋯, finding a correspondence between those from weak couplings and the latter. Furthermore, we investigated the strong-coupling regime taking into consideration the Wigner formulation. Elaborating on this, we obtained the profile to leading order, computed the contact density values as compared to those derived in an earlier work on the contact theorem. We formulated adequately the strong-coupling regime for this system that differed from previous formulations. Ultimately, we calculated the first-order corrections and compared them against numerical results from our simulations with very good agreement; these results compared equally well in the planar limit, whose results are well known.

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.