One-dimensional counterion gas between charged surfaces: exact results compared with weak- and strong-coupling analyses

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.

Ranked diffusion, delta Bose gas, and Burgers equation

We study the diffusion of N particles in one dimension interacting via a drift proportional to their rank. In the attractive case (self-gravitating gas) a mapping to the Lieb-Liniger quantum model allows one to obtain stationary time correlations, return probabilities, and the decay rate to the stationary state. The rank field obeys a Burgers equation, which we analyze. It allows one to obtain the stationary density at large N in an external potential V(x) (in the repulsive case). In the attractive case the decay rate to the steady state is found to depend on the initial condition if its spatial decay is slow enough. Coulomb gas methods allow one to study the final equilibrium at large N.

Relaxation dynamics of two interacting electrical double-layers in a 1D Coulomb system

We consider an out-of-equilibrium one-dimensional model for two electrical double-layers. With a combination of exact calculations and Brownian dynamics simulations, we compute the relaxation time (τ) for an electroneutral salt-free suspension, made up of two fixed colloids, withNneutralizing mobile counterions. ForNodd, the two double-layers never decouple, irrespective of their separationL; this is the regime of like-charge attraction, whereτexhibits a diffusive scaling inL²for largeL. On the other hand, for evenN,Lno longer is the relevant length scale for setting the relaxation time; this role is played by the Bjerrum length. This leads to distinctly different dynamics: forNeven, thermal effects are detrimental to relaxation, increasingτ, while they accelerate relaxation forNodd. Finally, we also show that the mean-field theory is recovered for largeNand moreover, that it remains an operational treatment down to relatively small values ofN(N> 3).

One-dimensional colloidal model with dielectric inhomogeneity

We consider a one-dimensional model allowing analytical derivation of the effective interactions between two charged colloids. We evaluate exactly the partition function for an electroneutral salt-free suspension with dielectric jumps at the colloids' position. We derive a contact relation with the pressure that shows there is like-charge attraction, whether or not the counterions are confined between the colloids. In contrast to the homogeneous dielectric case, there is the possibility for the colloids to attract despite the number of counterions (N) being even. The results are shown to recover the mean-field prediction in the limit N→∞.

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.

Interactions between charged particles with bathing multivalent counterions: experiments vs. dressed ion theory

We compare the recent experimentally measured forces between charged colloidal particles, as well as their effective surface potentials (surface charge) in the presence of multivalent counterions in a bathing monovalent salt solution, with the predictions of the dressed ion theory of strongly charged colloidal systems. The benchmark for comparison is provided by the DLVO theory and the deviations from its predictions at small separations are taken as an indication of the additional non-DLVO attractions that can be fitted by an additional phenomenological exponential term. The parameters characterizing this non-DLVO exponential term as well as the dependencies of the effective potential on the counterion concentration and valency predicted by the dressed ion theory are well within the experimental values. This suggests that the deviations from the DLVO theory are probably caused by ion correlations as formalized within the dressed ion theory.

Configurational and energy landscape in one-dimensional Coulomb systems

We study a one-dimensional Coulomb system, where two charged colloids are neutralized by a collection of point counterions, with global neutrality. With temperature being given, two situations are addressed: Either the colloids are kept at fixed positions (canonical ensemble) or the force acting on the colloids is fixed (isobaric-isothermal ensemble). The corresponding partition functions are worked out exactly, in view of determining which arrangement of counterions is optimal. How many counterions should be in the confined segment between the colloids? For the remaining ions outside, is there a left-right symmetry breakdown? We evidence a cascade of transitions as system size is varied in the canonical treatment or as pressure is increased in the isobaric formulation.

Surface Tension of Acid Solutions: Fluctuations beyond the Nonlinear Poisson-Boltzmann Theory

We extend our previous study of surface tension of ionic solutions and apply it to acids (and salts) with strong ion-surface interactions, as described by a single adhesivity parameter for the ionic species interacting with the interface. We derive the appropriate nonlinear boundary condition with an effective surface charge due to the adsorption of ions from the bulk onto the interface. The calculation is done using the loop-expansion technique, where the zero loop (mean field) corresponds of the full nonlinear Poisson-Boltzmann equation. The surface tension is obtained analytically to one-loop order, where the mean-field contribution is a modification of the Poisson-Boltzmann surface tension and the one-loop contribution gives a generalization of the Onsager-Samaras result. Adhesivity significantly affects both contributions to the surface tension, as can be seen from the dependence of surface tension on salt concentration for strongly absorbing ions. Comparison with available experimental data on a wide range of different acids and salts allows the fitting of the adhesivity parameter. In addition, it identifies the regime(s) where the hypotheses on which the theory is based are outside their range of validity.

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.

Counter-ions at single charged wall: Sum rules

For inhomogeneous classical Coulomb fluids in thermal equilibrium, like the jellium or the two-component Coulomb gas, there exists a variety of exact sum rules which relate the particle one-body and two-body densities. The necessary condition for these sum rules is that the Coulomb fluid possesses good screening properties, i.e. the particle correlation functions or the averaged charge inhomogeneity, say close to a wall, exhibit a short-range (usually exponential) decay. In this work, we study equilibrium statistical mechanics of an electric double layer with counter-ions only, i.e. a globally neutral system of equally charged point-like particles in the vicinity of a plain hard wall carrying a fixed uniform surface charge density of opposite sign. At large distances from the wall, the one-body and two-body counter-ion densities go to zero slowly according to the inverse-power law. In spite of the absence of screening, all known sum rules are shown to hold for two exactly solvable cases of the present system: in the weak-coupling Poisson-Boltzmann limit (in any spatial dimension larger than one) and at a special free-fermion coupling constant in two dimensions. This fact indicates an extended validity of the sum rules and provides a consistency check for reasonable theoretical approaches.

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.

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.

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.

Interaction regimes for oppositely charged plates with multivalent counterions

Within a mean-field treatment of the interaction between two oppositely charged plates in a salt-free solution, the distance at which a transition from an attractive to a repulsive regime appears can be computed analytically. The mean-field description, however, breaks down under strong Coulombic couplings, which can be achieved at room temperature with multivalent counterions and highly charged surfaces. Making use of the contact theorem and simple physical arguments, we propose explicit expressions for the equation of state in several situations at short distances. The possibility of Bjerrum pair formation is addressed and is shown to have profound consequences on the interactions. To complete the picture, we consider the large-distance limit, from which schematic phase diagram discriminating attractive from repulsive regions can be proposed.

Counter-ions at charged walls: two-dimensional systems

We study equilibrium statistical mechanics of classical point counter-ions, formulated on 2D Euclidean space with logarithmic Coulomb interactions (infinite number of particles) or on the cylinder surface (finite particle numbers), in the vicinity of a single uniformly charged line (one single double layer), or between two such lines (interacting double layers). The weak-coupling Poisson-Boltzmann theory, which applies when the coupling constant [Formula: see text] is small, is briefly recapitulated (the coupling constant is defined as [Formula: see text] [Formula: see text] [Formula: see text] e (2) , where [Formula: see text] is the inverse temperature, and e the counter-ion charge). The opposite limit ( [Formula: see text] [Formula: see text] ∞ is treated by using a recent method based on an exact expansion around the ground-state Wigner crystal of counter-ions. These two limiting results are compared at intermediary values of the coupling constant [Formula: see text] = 2[Formula: see text] ([Formula: see text] = 1, 2, 3) , to exact results derived within a 1D lattice representation of 2D Coulomb systems in terms of anti-commuting field variables. The models (density profile, pressure) are solved exactly for any particles numbers N at [Formula: see text] = 2 and up to relatively large finite N at [Formula: see text] = 4 and 6. For the one-line geometry, the decay of the density profile at asymptotic distance from the line undergoes a fundamental change with respect to the mean-field behavior at [Formula: see text] = 6 . The like-charge attraction regime, possible for large [Formula: see text] but precluded at mean-field level, survives for [Formula: see text] = 4 and 6, but disappears at [Formula: see text] = 2 .

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.

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

We present general arguments for the importance, or lack thereof, of structure in the charge distribution of counterions for counterion-mediated interactions between bounding symmetrically charged surfaces. We show that on the mean field or weak coupling level, the charge quadrupole contributes the lowest order modification to the contact value theorem and thus to the intersurface electrostatic interactions. The image effects are non-existent on the mean field level even with multipoles. On the strong coupling level the quadrupoles and higher order multipoles contribute additional terms to the interaction free energy only in the presence of dielectric inhomogeneities. Without them, the monopole is the only multipole that contributes to the strong coupling electrostatics. We explore the consequences of these statements in all their generality.