Colloidal Systems in Concentrated Electrolyte Solutions Exhibit Re-entrant Long-Range Electrostatic Interactions due to Underscreening

Exceptionally Strong Double-Layer Barriers Generated by Polyampholyte Salt

Experiments using the surface force apparatus have found anomalously long-range interactions between charged surfaces in concentrated salt solutions. Ion clustering has been suggested as a possible origin of this behavior. In this work, we demonstrate that if such stable clusters indeed form, they are able to induce remarkably strong free energy barriers under conditions where a corresponding solution of simple salt provides negligible forces. Our cluster model is based on connected ions producing a polyampholyte salt containing a symmetric mixture of monovalent cationic and anionic polyampholytes. Ion distributions and surface interactions are evaluated utilizing statistical-mechanical (classical) polymer density functional theory, cDFT. In the Supporting Information, we briefly investigate a range of different polymer architectures (connectivities), but in the main part of the work, a polyampholyte ion is modeled as a linear chain with alternating charges, in which the ends carry an identical charge (hence, a monovalent net charge). These salts are able to generate repulsions, between similarly charged surfaces, of a remarkable strength, exceeding those from simple salts by orders of magnitude. The underlying mechanism for this is the formation of brush-like layers at the surfaces, i.e., the repulsion is strongly related to excluded volume effects, in a manner similar to the interaction between surfaces carrying grafted polymers. We believe our results are relevant not only to possible mechanisms underlying anomalously long-ranged underscreening in concentrated simple salt solutions but also for the potential use of synthesized polyampholyte salt as extremely efficient stabilizers of colloidal dispersions.

Detecting ion-specific forces between fatty acid colloids and salt crystals in brines using colloidal probe AFM

HYPOTHESIS

Ion-specific forces in concentrated salt solutions play critical roles in many applications, ranging from biology to engineering, e.g., separating water-soluble minerals in brines by flotation using air bubbles. There should be some differences in colloidal forces between the surfactant precipitates in brines and NaCl crystals, and KCl crystals, making their selective aggregation and flotation separation possible.

EXPERIMENTS

Micron-sized spheres of fatty acid colloids were successfully prepared using lead laurate, characterized, and used to fabricate the AFM probes. Using a special AFM cell design and procedure, interaction force spectroscopy and force mapping on NaCl and KCl crystal surfaces using the probes were performed in their brines (7 M) and quantified using numerical solutions of advanced van der Waals and electrical double-layer theories to reveal valuable distributions of attractive, repulsive, and adhesive colloidal forces between the surfactant colloids and salt crystals.

FINDINGS

Attraction and adhesion between the lead laurate colloidal probe and the NaCl crystal surface were much stronger than those measured on the KCl crystal surface, explaining the selective separation between NaCl and KCl crystals by flotation in the brines. Theoretical analysis of the measured forces shows the potential role of ion-specific interactions in predicting selective aggregation and flotation separation. Our work provides an innovative approach to quantifying the intermolecular interactions between surfactant colloids and NaCl and KCl crystals, offering new theories on colloid and surface chemistry regarding ion-specific forces that underpin aggregation and separation in brines and beyond.

Ion-Ion Association in Bulk Mixed Electrolytes Using Global and Local Electroneutrality Constraints

Ion atmospheres play a critical role in modulating the interactions between charged components in solutions. However, a detailed description of the nature of ion atmospheres remains elusive. Here, we use Kirkwood-Buff theory, an exact theory of solution mixtures, together with a series of local and bulk electroneutrality constraints to provide relationships between all the net ion-ion distributions in bulk electrolyte mixtures. The validity of the underlying relationships is then confirmed using classical explicit solvent molecular simulations of a range of electrolyte mixtures. Further analysis indicates the ion distributions can be separated into two contributions, one resulting in charge neutralization, for which each ion contributes in proportion to its ionic strength, and the other accounting for all the solution thermodynamics. The relationships hold for atomic and molecular ions of any size and valency regardless of ionic strength, temperature, or pressure, in any solvent system.

Cluster Formation Induced by Local Dielectric Saturation in Restricted Primitive Model Electrolytes

Experiments using the Surface Force Apparatus (SFA) have found anomalously long-ranged charge-charge underscreening in concentrated salt solutions. Meanwhile, theory and simulations have suggested ion clustering to be a possible origin of this behavior. The popular Restricted Primitive Model of electrolyte solutions, in which the solvent is represented by a uniform relative dielectric constant, εr, is unable to resolve the anomalous underscreening seen in experiments. In this work, we modify the Restricted Primitive Model to account for local dielectric saturation within the ion hydration shell. The dielectric "constant" in our model locally decreases from the bulk value to a lower saturated value at the ionic surface. The parameters for the model are deduced so that typical salt solubilities are obtained. Our simulations for both bulk and slit geometries show that our model displays strong cluster formation and these give rise to long-ranged density correlations between charged surfaces, at distances similar to what has been observed in SFA measurements. An electrolyte model wherein the dielectric constant remains uniform does not display similar clusters, even with εr equal to the low saturated value at ion contact. Hence, the observed behaviors are not simply due to an enhanced Coulomb interaction. In the latter case, cluster growth is counteracted by long-ranged repulsions between like-charged ions within clusters; this is an effect that is considerably reduced when the dielectric response drop is local. Our results imply that long-ranged interactions in these systems are mainly due to cluster-cluster correlations, rather than large electrostatic screening lengths.

An efficient method to establish electrostatic screening lengths of restricted primitive model electrolytes

We present a novel, and computationally cheap, way to estimate electrostatic screening lengths from simulations of restricted primitive model (RPM) electrolytes. We demonstrate that the method is accurate by comparisons with simulated long-ranged parts of the charge density, at various Bjerrum lengths, salt concentrations and ion diameters. We find substantial underscreening in low dielectric solvent, but with an "aqueous" solvent, there is instead overscreening, the degree of which increases with ion size. Our method also offers a possible path to (future) more accurate classical density functional treatments of ionic fluids.

Ion specific tuning of nanoparticle dispersion in an ionic liquid: a structural, thermoelectric and thermo-diffusive investigation

Dispersions of charged maghemite nanoparticles (NPs) in EAN (ethylammonium nitrate) a reference Ionic Liquid (IL) are studied here using a number of static and dynamical experimental techniques; small angle scattering (SAS) of X-rays and of neutrons, dynamical light scattering and forced Rayleigh scattering. Particular insight is provided regarding the importance of tuning the ionic species present at the NP/IL interface. In this work we compare the effect of Li+, Na+ or Rb+ ions. Here, the nature of these species has a clear influence on the short-range spatial organisation of the ions at the interface and thus on the colloidal stability of the dispersions, governing both the NP/NP and NP/IL interactions, which are both evaluated here. The overall NP/NP interaction is either attractive or repulsive. It is characterised by determining, thanks to the SAS techniques, the second virial coefficient A2, which is found to be independent of temperature. The NP/IL interaction is featured by the dynamical effective charge ξeff0 of the NPs and by their entropy of transfer ŜNP (or equivalently their heat of transport ) determined here thanks to thermoelectric and thermodiffusive measurements. For repulsive systems, an activated process rules the temperature dependence of these two latter quantities.

Scaling perspectives of underscreening in concentrated electrolyte solutions

We present a scaling view of underscreening observed in salt solutions in the range of concentrations greater than about 1 M, in which the screening length increases with concentration. The system consists of hydrated clusters of positive and negative ions with a single unpaired ion as suggested by recent simulations. The environment of this ion is more hydrated than average which leads to a self-similar situation in which the size of this environment scales with the screening length. The prefactor involves the local dielectric constant and the cluster density. The scaling arguments as well as the cluster model lead to scaling of the screening length with the ion concentration, in agreement with observations.

A screening of results on the decay length in concentrated electrolytes

Screening of electrostatic interactions in room-temperature ionic liquids and concentrated electrolytes has recently attracted much attention as surface force balance experiments have suggested the emergence of unanticipated anomalously large screening lengths at high ion concentrations. Termed underscreening, this effect was ascribed to the bulk properties of concentrated ionic systems. However, underscreening under experimentally relevant conditions is not predicted by classical theories and challenges our understanding of electrostatic correlations. Despite the enormous effort in performing large-scale simulations and new theoretical investigations, the origin of the anomalously long-range screening length remains elusive. This contribution briefly summarises the experimental, analytical and simulation results on ionic screening and the scaling behaviour of screening lengths. We then present an atomistic simulation approach that accounts for the solvent and ion exchange with a reservoir. We find that classical density functional theory (DFT) for concentrated electrolytes under confinement reproduces ion adsorption at charged interfaces surprisingly well. With DFT, we study confined electrolytes using implicit and explicit solvent models and the dependence on the solvent's dielectric properties. Our results demonstrate how the absence vs. presence of solvent particles and their discrete nature affect the short and long-range screening in concentrated ionic systems.

Underscreening in concentrated electrolytes: re-entrant swelling in polyelectrolyte brushes

Hypersaline environments are ubiquitous in nature and are found in myriad technological processes. Recent empirical studies have revealed a significant discrepancy between predicted and observed screening lengths at high salt concentrations, a phenomenon referred to as underscreening. Herein we investigate underscreening using a cationic polyelectrolyte brush as an exemplar. Poly(2-(methacryloyloxy)ethyl)trimethylammonium (PMETAC) brushes were synthesised and their internal structural changes and swelling response was monitored with neutron reflectometry and spectroscopic ellipsometry. Both techniques revealed a monotonic brush collapse as the concentration of symmetric monovalent electrolyte increased. However, a non-monotonic change in brush thickness was observed in all multivalent electrolytes at higher concentrations, known as re-entrant swelling; indicative of underscreening. For all electrolytes, numerical self-consistent field theory predictions align with experimental studies in the low-to-moderate salt concentration regions. Analysis suggests that the classical theory of electrolytes is insufficient to describe the screening lengths observed at high salt concentrations and that the re-entrant polyelectrolyte brush swelling seen herein is consistent with the so-called regular underscreening phenomenon.

Solvent Effects on Structure and Screening in Confined Electrolytes

Using classical density functional theory, we investigate the influence of solvent on the structure and ionic screening of electrolytes under slit confinement and in contact with a reservoir. We consider a symmetric electrolyte with implicit and explicit solvent models and find that spatially resolving solvent molecules is essential for the ion structure at confining walls, excess ion adsorption, and the pressure exerted on the walls. Despite this, we observe only moderate differences in the period of oscillations of the pressure with the slit width and virtually coinciding decay lengths as functions of the scaling variable σ_{ion}/λ_{D}, where σ_{ion} is the ion diameter and λ_{D} the Debye length. Moreover, in the electrostatic-dominated regime, this scaling behavior is practically independent of the relative permittivity and its dependence on the ion concentration. In contrast, the crossover to the hard-core-dominated regime depends sensitively on all three factors.

Arbitrary-Shape Dielectric Particles Interacting in the Linearized Poisson-Boltzmann Framework: An Analytical Treatment

This work considers the interaction of two dielectric particles of arbitrary shape immersed into a solvent containing a dissociated salt and assuming that the linearized Poisson-Boltzmann equation holds. We establish a new general spherical re-expansion result which relies neither on the conventional condition that particle radii are small with respect to the characteristic separating distance between particles nor on any symmetry assumption. This is instrumental in calculating suitable expansion coefficients for the electrostatic potential inside and outside the objects and in constructing small-parameter asymptotic expansions for the potential, the total electrostatic energy, and forces in ascending order of Debye screening. This generalizes a recent result for the case of dielectric spheres to particles of arbitrary shape and builds for the first time a rigorous (exact at the Debye-Hückel level) analytical theory of electrostatic interactions of such particles at arbitrary distances. Numerical tests confirm that the proposed theory may also become especially useful in developing a new class of grid-free, fast, highly scalable solvers.