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.

Charge Selectivity of an Ionic Transistor

The charge selective properties of a long planar nanochannel with an embedded finite uniformly charged section in the middle are studied. The probability flux of a single test ion initially confined to the inlet reservoir is determined by integrating the Smoluchowski equation using a previously published series solution for the Debye-Hückel potential in this geometry. The charge selective properties are characterized by a dimensionless quantity that we call the "fractional blockage". We study how the fractional blockage depends on the dimensionless parameters that characterize the charge state and channel geometry. In the limit of strongly overlapped wall Debye layers, analytical expressions for the fractional blockage are presented that are found to be in good agreement with numerically computed values in the appropriate asymptotic regimes. These results may be helpful in the design of nanofluidic devices that have a variety of applications.

Debye-Hückel Free Energy of an Electric Double Layer with Discrete Charges Located at a Dielectric Interface

Poisson-Boltzmann theory provides an established framework to calculate properties and free energies of an electric double layer, especially for simple geometries and interfaces that carry continuous charge densities. At sufficiently small length scales, however, the discreteness of the surface charges cannot be neglected. We consider a planar dielectric interface that separates a salt-containing aqueous phase from a medium of low dielectric constant and carries discrete surface charges of fixed density. Within the linear Debye-Hückel limit of Poisson-Boltzmann theory, we calculate the surface potential inside a Wigner-Seitz cell that is produced by all surface charges outside the cell using a Fourier-Bessel series and a Hankel transformation. From the surface potential, we obtain the Debye-Hückel free energy of the electric double layer, which we compare with the corresponding expression in the continuum limit. Differences arise for sufficiently small charge densities, where we show that the dominating interaction is dipolar, arising from the dipoles formed by the surface charges and associated counterions. This interaction propagates through the medium of a low dielectric constant and alters the continuum power of two dependence of the free energy on the surface charge density to a power of 2.5 law.

Electrolytes in regimes of strong confinement: surface charge modulations, osmotic equilibrium and electroneutrality

In the present work, we study an electrolyte solution confined between planar surfaces with nanopatterned charged domains, which has been connected to a bulk ionic reservoir. The system is investigated through an improved Monte Carlo (MC) simulation method, suitable for simulation of electrolytes in the presence of modulated surface charge distributions. We also employ a linear approach in the spirit of the classical Debye-Hückel approximation, which allows one to obtain explicit expressions for the averaged potentials, ionic profiles, effective surface interactions and the net ionic charge confined between the walls. Emphasis is placed on the limit of strongly confined electrolytes, in which case local electroneutrality in the inter-surface space might not be fulfilled. In order to access the effects of such a lack of local charge neutrality on the ion-induced interactions between surfaces with modulated charge domains, we consider two distinct model systems for the confined electrolyte: one in which a salt reservoir is explicitly taken into account via the osmotic equilibrium with an electrolyte of fixed bulk concentration, and a second one in which the equilibrium with a charge neutral ionic reservoir is implicitly considered. While in the former case the osmotic ion exchange might lead to non-vanishing net charges, in the latter model charge neutrality is enforced through the appearance of an implicit Donnan potential across the charged interfaces. A strong dependence of the ion-induced surface interactions on the employed model system is observed at all surface separations. These findings strongly suggest that due care is to be taken while choosing among different scenarios to describe the ion exchange in electrolytes confined between charged surfaces, even in cases when the monopole (non zero net charge) surface contributions are absent.

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.

Exclusion-Enrichment Effect in Ionic Transistors

A simple model of a nanofluidic transistor consisting of a uniformly charged central section between a pair of plane parallel walls is considered. The linearized Poisson-Boltzmann equation corresponding to weak surface charge is solved exactly, and the solution is presented as an infinite series. The problem is characterized by three dimensionless parameters, namely the normalized surface charge, the ratio of the channel width to the Debye length, and the length-to-width aspect ratio of the charged section. The first of these parameters is presumed small, but the other two are arbitrary. The dependence of the exclusion-enrichment effect on these three parameters is discussed.

Effect of Nonzero Solid Permittivity on the Electrical Repulsion between Charged Surfaces

The electrical repulsion between two charged solid surfaces separated by an electrolyte is studied as a function of the permittivity ϵ_s of the solid in the limit in which potentials are small, and the gap between the plane solid surfaces is small compared to the Debye length κ^-1 within the electrolyte. The solid surfaces are uniformly charged in a central region |x|< L outside which they are uncharged. When ϵ_s = 0, ions from the charge cloud between the charged surfaces spill out into regions of length O(κ^-1) beyond x = ± L, thereby reducing the pressure between the surfaces from that predicted by Derjaguin-Landau-Verwey-Overbeek theory for infinite, uniformly charged surfaces. When ϵ_s>0, ions spill out over much larger O(L) regions, thereby reducing still further both the electrical potential between the solid surfaces and the repulsive force between them. However, this reduction becomes smaller as κL becomes large.

Experimental Evidence for Algebraic Double-Layer Forces

According to conventional wisdom, electric double-layer forces normally decay exponentially with separation distance. Here, we present experimental evidence of algebraically decaying double-layer interactions. We show that algebraic interactions arise in both strongly overlapping as well as counterion-only regimes, albeit the evidence is less clear for the former regime. In both of these cases, the disjoining pressure profile assumes an inverse square distance dependence. At small separation distances, another algebraic regime is recovered. In this regime, the pressure decays as the inverse of separation distance.

The effect of the finite size of ions and Debye layer overspill on the screened Coulomb interactions between charged flat plates

Screened repulsion between uniformly charged plates with an intervening electrolyte is analyzed for strongly overlapped electrical double layers (EDL), accounting for the steric effect of ions and their expulsion from EDL edges into the surrounding solution. As a generalization of a study by Philipse et al. which does not account for these effects, an analytical expression is derived for the repulsion pressure in the limit of infinitely long plates with a zero-field assumption, which agrees closely with the corresponding numerical solution at low inter-plate separations. Our results show an augmented repulsive pressure for finite-sized ions at strong EDL overlaps. For plates with a finite lateral size, we demonstrate a further extended domain of low inter-plate gaps where the repulsion pressure increases with ion size due to a strong interplay between the steric interaction of ions and the EDL overspill phenomenon, considered earlier in a study by Ghosal & Sherwood limited to the linear Debye-Hückel regime (which cannot account for the steric effect of ions). This investigation on a simple model should enhance our understanding of the interaction between charged particles in electrophoresis, nanoscale self-assembly, active particles, and various other electrokinetic systems.

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.