Asymmetric Coulomb fluids at randomly charged dielectric interfaces: Anti-fragility, overcharging and charge inversion

Structure of electrolyte solutions at nonuniformly charged surfaces on a variety of length scales

The structures of dilute electrolyte solutions close to nonuniformly charged planar substrates are systematically studied within the entire spectrum of microscopic to macroscopic length scales by means of a unified classical density functional theory approach. This is in contrast to previous investigations, which are applicable either to short or to long length scales. It turns out that interactions with microscopic ranges, e.g., due to the hard cores of the fluid molecules and ions, have a negligible influence on the formation of nonuniform lateral structures of the electrolyte solutions. This partly justifies the Debye-Hückel approximation schemes applied in previous studies of that system. In general, a coupling between the lateral and the normal fluid structures leads to the phenomenology that, upon increasing the distance from the substrate, fewer details of the lateral nonuniformities contribute to the fluid structure, such that ultimately only large-scale surface features remain relevant. It can be expected that this picture also applies to other fluids characterized by several length scales.

Electrostatic potential between charged particles at an oil-water interface

Electrostatic interactions between point charges embedded into interfaces separating dielectric media are omnipresent in soft matter systems and often control their stability. Such interactions are typically complicated and do not resemble their bulk counterparts. For instance, the electrostatic potential of a point charge at an air-water interface falls off as r^{-3}, where r is the distance from the charge, exhibiting a dipolar behavior. This behavior is often assumed to be generic, and is widely referred to when interpreting experimental results. Here we explicitly calculate the in-plane potential of a point charge at an interface between two electrolyte solutions with different, finite dielectric permittivities and Debye screening lengths, such as oil and water. We show that the asymptotic behavior of this potential is neither a dipole, which characterizes the potential at air-water interfaces, nor a screened monopole, which describes the bulk behavior in a single electrolyte solution. By considering the same problem in arbitrary dimensions, we find that the physics behind this difference can be traced to the asymmetric propagation of the interaction in the two media. Our results should be relevant to understand the effective potential acting between interfacial proteins in biofilms, and the self-assembly of charged colloids at droplet surfaces in oil-water emulsions.

Heterogeneous surface charge confining an electrolyte solution

The structure of dilute electrolyte solutions close to a surface carrying a spatially inhomogeneous surface charge distribution is investigated by means of classical density functional theory within the approach of fundamental measure theory. For electrolyte solutions, the influence of these inhomogeneities is particularly strong because the corresponding characteristic length scale is the Debye length, which is large compared to molecular sizes. Here, a fully three-dimensional investigation is performed, which accounts explicitly for the solvent particles, and thus provides insight into effects caused by ion-solvent coupling. The present study introduces a versatile framework to analyze a broad range of types of surface charge heterogeneities even beyond the linear response regime. This reveals a sensitive dependence of the number density profiles of the fluid components and of the electrostatic potential on the magnitude of the charge as well as on the details of the surface charge patterns at small scales.

Long-range forces and charge inversions in model charged colloidal dispersions at finite concentration

In charged colloidal dispersion systems the interest is in finding their stability conditions, phase transitions, and transport properties, either in bulk or confinement, among other physicochemical quantities, for which the knowledge of the dispersions' molecular structure and the associated macroion-macroion forces is crucial. To investigate these phenomena simple models have been proposed. Most of the theoretical and simulation studies on charged particles suspensions are at infinite dilution conditions. Hence, these studies have been focused on the electrolyte structure around one or two isolated central particle(s), where phenomena as charge reversal, charge inversion and surface charge amplification have been shown to be relevant. However, experimental studies at finite volume fraction exhibit interesting phenomenology which imply very long-range correlations. A simple, yet useful, model is the Colloidal Primitive Model, in which the colloidal dispersion is modeled as a mixture of size (and charge) asymmetrical hard spheres, at finite volume fraction. In this paper we review recent integral equations solutions for this model, where very long-range attractive-repulsive forces, as well as new long-range, giant charge inversions are reported. The calculated macroions radial distribution functions, charge distributions, and macroion-macroion forces are qualitatively consistent with existing experimental results, and Monte Carlo and molecular dynamics simulations.

A Monte Carlo study of the electrical double layer of a shape-asymmetric electrolyte around a spherical colloid

In this paper, we present a Monte Carlo simulation study on the structure of the electrical double layer around a spherical colloid surrounded by a binary electrolyte composed of spherical and non-spherical ions. Results are provided for the radial distribution functions between the colloid and ions, the orientation correlations between the colloid and non-spherical particles, and the integrated charge. Work is reported mainly for non-spherical particles modeled as spherocylinders, although a particular comparison is made between spherocylindrical particles and dimers. For the conditions investigated here, spherocylinders and dimers produce essentially the same structural information. Additionally, it is shown that spherocylinders mostly orient tangentially to the colloid at its surface; this preferred orientation disappears for larger distances. We also evidence that, near the colloid, the integrated charge attenuates monotonically when the macroparticle is highly charged, whereas for intermediate and low charged states of the colloid, the integrated charge can display charge reversal, overcharging, or both, with magnitudes that are sensitive to the salt concentration and to the localization of charge inside the spherocylinders.

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.

Outsized Amplitude-Modulated Structure of Very-Long-Range Charge Inversions in Model Colloidal Dispersions

Most theoretical and simulation studies on charged suspensions are at infinite dilution and are focused on the electrolyte structure around one or two isolated particles. Some classic experimental studies with latex particle solutions exhibit interesting phenomenology which imply very-long-range correlations. Here, we apply an integral equation theory to a model charged macroion suspension, at finite volume fraction, and find an amplitude-modulated charge inversion structure, with outsized amplitudes and of very-long-range extension. These inversions are different from the standard charge inversions in that they occur at finite macroions' volume fraction, far away from the central macroion, are outsized, and increase, not decrease, with increasing particle charge and distance to the central particle, which is indicative of long-range correlations. We find our results to be in agreement with our Monte Carlo simulations and qualitatively consistent with existing experimental results.

Electrolyte solutions at heterogeneously charged substrates

The influence of a chemically or electrically heterogeneous distribution of interaction sites at a planar substrate on the number density of an adjacent fluid is studied by means of classical density functional theory (DFT). In the case of electrolyte solutions the effect of this heterogeneity is particularly long ranged, because the corresponding relevant length scale is set by the Debye length which is large compared to molecular sizes. The DFT used here takes the solvent particles explicitly into account and thus captures phenomena, inter alia, due to ion-solvent coupling. The present approach provides closed analytic expressions describing the influence of chemically and electrically nonuniform walls. The analysis of isolated δ-like interactions, isolated interaction patches, and hexagonal periodic distributions of interaction sites reveals a sensitive dependence of the fluid density profiles on the type of the interaction, as well as on the size and the lateral distribution of the interaction sites.

On the physics of both surface overcharging and charge reversal at heterophase interfaces

A series of Monte Carlo simulations are employed to reveal the physics of both surface overcharging and charge reversal at a negatively charged dielectric interface exposed to a bulk solution containing a +2:−1 electrolyte in the absence and presence of a monovalent salt.

Electrostatics of patchy surfaces

In the study of colloidal, biological and electrochemical systems, it is customary to treat surfaces, macromolecules and electrodes as homogeneously charged. This simplified approach is proven successful in most cases, but fails to describe a wide range of heterogeneously charged surfaces commonly used in experiments. For example, recent experiments have revealed a long-range attraction between overall neutral surfaces, locally charged in a mosaic-like structure of positively and negatively charged domains ("patches"). Here, we review experimental and theoretical studies addressing the stability of heterogeneously charged surfaces, their effect on ionic profiles in solution, and the interaction between two such surfaces. We focus on electrostatics, and highlight the important new physical parameters appearing in the heterogeneous case, such as the largest patch size and inter-surface charge correlations.

DNA like-charge attraction and overcharging by divalent counterions in the presence of divalent co-ions

Strongly correlated electrostatics of DNA systems has drawn the interest of many groups, especially the condensation and overcharging of DNA by multivalent counterions. By adding counterions of different valencies and shapes, one can enhance or reduce DNA overcharging. In this paper, we focus on the effect of multivalent co-ions, specifically divalent co-ions such as SO[Formula: see text]. A computational experiment of DNA condensation using Monte Carlo simulation in grand canonical ensemble is carried out where the DNA system is in equilibrium with a bulk solution containing a mixture of salt of different valency of co-ions. Compared to systems with purely monovalent co-ions, the influence of divalent co-ions shows up in multiple aspects. Divalent co-ions lead to an increase of monovalent salt in the DNA condensate. Because monovalent salts mostly participate in linear screening of electrostatic interactions in the system, more monovalent salt molecules enter the condensate leads to screening out of short-range DNA-DNA like charge attraction and weaker DNA condensation free energy. The overcharging of DNA by multivalent counterions is also reduced in the presence of divalent co-ions. Strong repulsions between DNA and divalent co-ions and among divalent co-ions themselves lead to a depletion of negative ions near the DNA surface as compared to the case without divalent co-ions. At large distances, the DNA-DNA repulsive interaction is stronger in the presence of divalent co-ions, suggesting that divalent co-ions' role is not only that of simple stronger linear screening.

The non-dominance of counterions in charge-asymmetric electrolytes: non-monotonic precedence of electrostatic screening and local inversion of the electric field by multivalent coions

The asymptotic convergence of the thermodynamic and structural properties of unequally-sized charge-symmetric ions in strong electric fields was postulated more than thirty years ago by Valleau and Torrie as the dominance of counterions via the non-linear Poisson-Boltzmann theory [Valleau and Torrie, J. Chem. Phys., 1982, 76, 4623]. According to this mean field prescription, the properties of the electrical double layer near a highly charged electrode immersed in a size-asymmetric binary electrolyte converge to those of a size-symmetric electrolyte if the properties of counterions are the same in both instances. On the other hand, some of the present authors have shown that, in fact, counterions do not dominate the electrical properties of a spherical macroion in the presence of unequally-sized ions, symmetric in valence, if ion correlations and ionic excluded volume effects are taken into account consistently. These ingredients are neglected in the classical Poisson-Boltzmann picture. In the present work, we show the occurrence of the non-dominance of counterions in the opposite scenario, that is, when ions are equally-sized but asymmetric in valence. This is performed in the presence of highly charged colloidal surfaces of spherical and planar geometries for different ionic volume fractions. In addition to the phenomenon of non-dominance of counterions, our simulations and theoretical data also exhibit a non-monotonic order or precedence in the mean electrostatic potential, or electrostatic screening, at the Helmholtz plane of a charged colloid. This interesting behaviour is analyzed as a function of the coion's valence, the ionic volume fraction, and the charge and size of the colloidal particle. All these phenomena are explained in terms of the decay of the electric field near the colloidal surface, and by the appearance of a local inversion of both the electric field and the integrated surface charge density of the colloidal particle in the presence of monovalent counterions and multivalent coions.

Charge reversal at a planar boundary between two dielectrics

Despite the ubiquitous character and relevance of the electric double layer in the entire realm of interface and colloid science, very little is known of the effect that surface heterogeneity exerts on the underlying mechanisms of ion adsorption. Herein, computer simulations offer a perspective that, in sharp contrast to the homogeneously charged surface, discrete groups promote multivalent counterion binding, leading to charge reversal but possibly having not a sign change of the electrophoretic mobility. Counterintuitively, the introduction of dielectric images yields a significantly greater accumulation of counterions, which further facilitates the magnitude of charge reversal. The reported results are very sensitive to both the degree of ion hydration and the representation of surface charges. Our findings shed light on the mechanism for charge reversal over a broad range of coupling regimes operating the adsorption of counterions through surface group bridging attraction with their own images and provide opportunities for experimental studies and theoretical development.

Ion-mediated interactions between net-neutral slabs: Weak and strong disorder effects

We investigate the effective interaction between two randomly charged but otherwise net-neutral, planar dielectric slabs immersed in an asymmetric Coulomb fluid containing a mixture of mobile monovalent and multivalent ions. The presence of charge disorder on the apposed bounding surfaces of the slabs leads to substantial qualitative changes in the way they interact, as compared with the standard picture provided by the van der Waals and image-induced, ion-depletion interactions. While, the latter predict purely attractive interactions between strictly neutral slabs, we show that the combined effects from surface charge disorder, image depletion, Debye (or salt) screening, and also, in particular, their coupling with multivalent ions, give rise to a more diverse behavior for the effective interaction between net-neutral slabs at nano-scale separations. Disorder effects show large variation depending on the properly quantified strength of disorder, leading either to non-monotonic effective interaction with both repulsive and attractive branches when the surface charges are weakly disordered (small disorder variance) or to a dominating attractive interaction that is larger both in its range and magnitude than what is predicted from the van der Waals and image-induced, ion-depletion interactions, when the surfaces are strongly disordered (large disorder variance).

Strong coupling electrostatics for randomly charged surfaces: antifragility and effective interactions

We study the effective interaction mediated by strongly coupled Coulomb fluids between dielectric surfaces carrying quenched, random monopolar charges with equal mean and variance, both when the Coulomb fluid consists only of mobile multivalent counterions and when it consists of an asymmetric ionic mixture containing multivalent and monovalent (salt) ions in equilibrium with an aqueous bulk reservoir. We analyze the consequences that follow from the interplay between surface charge disorder, dielectric and salt image effects, and the strong electrostatic coupling that results from multivalent counterions on the distribution of these ions and the effective interaction pressure they mediate between the surfaces. In a dielectrically homogeneous system, we show that the multivalent counterions are attracted towards the surfaces with a singular, disorder-induced potential that diverges logarithmically on approach to the surfaces, creating a singular but integrable counterion density profile that exhibits an algebraic divergence at the surfaces with an exponent that depends on the surface charge (disorder) variance. This effect drives the system towards a state of lower thermal 'disorder', one that can be described by a renormalized temperature, exhibiting thus a remarkable antifragility. In the presence of an interfacial dielectric discontinuity, the singular behavior of counterion density at the surfaces is removed but multivalent counterions are still accumulated much more strongly close to randomly charged surfaces as compared with uniformly charged ones. The interaction pressure acting on the surfaces displays in general a highly non-monotonic behavior as a function of the inter-surface separation with a prominent regime of attraction at small to intermediate separations. This attraction is caused directly by the combined effects from charge disorder and strong coupling electrostatics of multivalent counterions, which dominate the surface-surface repulsion due to the (equal) mean charges on the two surfaces and the osmotic pressure of monovalent ions residing between them. These effects can be quite significant even with a small degree of surface charge disorder relative to the mean surface charge. The strong coupling, disorder-induced attraction is typically much stronger than the van der Waals interaction between the surfaces, especially within a range of several nanometers for the inter-surface separation, where such effects are predicted to be most pronounced.