Ionic condensation and charge renormalization in colloidal suspensions

Effective interactions, structure, and pressure in charge-stabilized colloidal suspensions: Critical assessment of charge renormalization methods

Charge-stabilized colloidal suspensions display a rich variety of microstructural and thermodynamic properties, which are determined by electro-steric interactions between all ionic species. The large size asymmetry between molecular-scale microions and colloidal macroions allows the microion degrees of freedom to be integrated out, leading to an effective one-component model of microion-dressed colloidal quasi-particles. For highly charged colloids with strong macroion-microion correlations, nonlinear effects can be incorporated into effective interactions by means of charge renormalization methods. Here, we compare and partially extend several practical mean-field methods of calculating renormalized colloidal interaction parameters, including effective charges and screening constants, as functions of concentration and ionic strength. Within the one-component description, we compute structural and thermodynamic properties from the effective interactions and assess the accuracy of the different methods by comparing predictions with elaborate primitive-model simulations [P. Linse, J. Chem. Phys. 113, 4359 (2000)]. We also compare various prescriptions for the osmotic pressure of suspensions in Donnan equilibrium with a salt ion reservoir and analyze instances where the macroion effective charge becomes larger than the bare one. The methods assessed include single-center cell, jellium, and multi-center mean-field theories. The strengths and weaknesses of the various methods are critically assessed, with the aim of guiding optimal and accurate implementations.

N,N,N-Trimethyl chitosan as a permeation enhancer for inhalation drug delivery: Interaction with a model pulmonary surfactant

N,N,N-Trimethyl chitosan (TMC), a biocompatible and biodegradable derivative of chitosan, is currently used as a permeation enhancer to increase the translocation of drugs to the bloodstream in the lungs. This article discusses the effect of TMC on a mimetic pulmonary surfactant, Curosurf®, a low-viscosity lipid formulation administered to preterm infants with acute respiratory distress syndrome. Curosurf® exhibits a strong interaction with TMC, resulting in the formation of aggregates at electrostatic charge stoichiometry. At nanoscale, Curosurf® undergoes a profound reorganization of its lipid vesicles in terms of size and lamellarity. The initial micron-sized vesicles (average size 4.8 μm) give way to a froth-like network of unilamellar vesicles about 300 nm in size. Under such conditions, neutralization of the cationic charges by pulmonary surfactant may inhibit TMC permeation enhancer capacity, especially as electrostatic charge complexation is found at low TMC content. The permeation properties of pulmonary surfactant-neutralized TMC should then be evaluated for its applicability as a permeation enhancer for inhalation in the alveolar region.

Electrostatics of soft (bio)interfaces: Corrections of mean-field Poisson-Boltzmann theory for ion size, dielectric decrement and ion-ion correlation

HYPOTHESIS

Electrostatics of soft (ion-permeable) (bio)particles (e.g. microorganisms, core/shell colloids) in aqueous electrolytes is commonly formulated by the mean-field Poisson-Boltzmann theory and integration of the charge contributions from electrolyte ions and soft material. However, the effects connected to the size of the electrolyte ions and that of the structural charges carried by the particle, to dielectric decrement and ion-ion correlations on soft interface electrostatics have been so far considered at the margin, despite the limits of the Gouy theory for condensed and/or multivalent electrolytes.

EXPERIMENTS

Accordingly, we modify herein the Poisson-Boltzmann theory for core/shell (bio)interfaces to include the aforementioned molecular effects considered separately or concomitantly. The formalism is applicable for poorly to highly charged particles in the thin electric double layer regime and to unsymmetrical multivalent electrolytes.

FINDINGS

Computational examples of practical interests are discussed with emphasis on how each considered molecular effect or combination thereof affects the interfacial potential distribution depending on size and valence of cations and anions, size of particle charges, length scale of ionic correlations and shell-to-Debye layer thickness ratio. The origins of here-evidenced pseudo-harmonic potential profile and ion size-dependent screening of core/shell particle charges are detailed. In addition, the existence and magnitude of the Donnan potential when reached in the shell layer are shown to depend on the excluded volumes of the electrolyte ions.

Charge-Selective Aggregation Behavior of Thermoresponsive Polyelectrolytes Having Low Charge Density in Aqueous Solutions of Organic Counterions

The aggregation behavior of thermoresponsive polyelectrolytes with low charge density in aqueous solutions of organic counterions was investigated. We synthesized two thermoresponsive polyelectrolytes: anionic poly(N-isopropylacrylamide-co-(3-sulfopropyl)acrylamide potassium) (P-NIP-SPAK) and cationic poly(N-isopropylacrylamide-co-(3-acrylamidepropyl)trimethylammonium chloride) (P-NIP-AAPTAC). The polyelectrolytes remained soluble in their aqueous solutions even above the lower critical soluble temperature of P-NIP owing to the strong hydration property of the ionic groups. The aggregation occurred when organic counterions were added to the solution. In these solution systems, the concentration of counterions exceeds those of ionic groups introduced into the polyelectrolytes. The aggregation behavior is attributed to the salting-out effect of counterions accommodated near the polyelectrolyte surface by electrostatic interaction. This aggregation behavior was utilized for the charge-selective recognition of amino acids. P-NIP-SPAK aggregated only when basic amino acids were added under acidic conditions, whereas P-NIP-AAPTAC aggregated only when acidic amino acids were added under basic conditions. The results herein demonstrate that P-NIP-SPAK and P-NIP-AAPTAC have the potential to be used as charge-selective polymer sensors for amino acids without having to strictly control the experimental conditions.

Far-Field Electrostatic Signatures of Macromolecular 3D Conformation

In solution as in vacuum, the electrostatic field distribution in the vicinity of a charged object carries information on its three-dimensional geometry. We report on an experimental study exploring the effect of molecular shape on long-range electrostatic interactions in solution. Working with DNA nanostructures carrying approximately equal amounts of total charge but each in a different three-dimensional conformation, we demonstrate that the geometry of the distribution of charge in a molecule has substantial impact on its electrical interactions. For instance, a tetrahedral structure, which is the most compact distribution of charge we tested, can create a far-field effect that is effectively identical to that of a rod-shaped molecule carrying half the amount of total structural charge. Our experiments demonstrate that escape-time electrometry (ETe) furnishes a rapid and facile method to screen and identify 3D conformations of charged biomolecules or molecular complexes in solution.

Modeling charge separation in charged nanochannels for single-molecule electrometry

We model the transport of electrically charged solute molecules by a laminar flow within a nanoslit microfluidic channel with electrostatic surface potential. We derive the governing convection–diffusion equation, solve it numerically, and compare it with a Taylor–Aris-like approximation, which gives excellent results for small Péclet numbers. We discuss our results in light of designing an assay that can measure simultaneously the hydrodynamic size and electric charge of single molecules by tracking their motion in such nanoslit channels with electrostatic surface potential.

Where in the world are condensed counterions?

A scaling model of the concentration profiles of both condensed and free counterions is presented for solutions of spherical and cylindrical charged nanoparticles of different charge valences, nanoparticle sizes, and salt concentrations. The distribution of counterions for both spherical and cylindrical charged particles in salt-free solutions is determined by the condensation parameter γ₀ defined as the ratio of nanoparticle valence Z₀ to the number of Bjerrum lengths l_B = e²/(εkT) per nanoparticle size (γ₀ = Z₀l_B/(2r₀) for spherical nanoparticles with radii r₀ or γ₀ = Z₀l_B/L for cylindrical particles with length L), where ε is solution dielectric permittivity, e is elementary charge and kT is thermal energy. Depending on the magnitudes of the condensation parameter γ₀ and nanoparticle volume fraction ϕ, we find three qualitatively different regimes for the counterion distribution near charged particles: (i) weakly charged particles with no condensed counterions, (ii) regime of weak counterion condensation with less than half of the counterions condensed, and (iii) regime of strong counterion condensation with the majority of counterions condensed. The magnitude of electrostatic energy of a condensed counterion with respect to solution locations with zero electric field is larger than thermal energy kT, and the fraction of condensed counterions increases from less than half in the weak condensation regime to the majority of all counterions in the strong condensation regime. The condensed counterions are not bound to the nanoparticle surface but instead are localized within the condensed counterion zone near the charged particle. The thickness of the condensed counterion zone varies with the condensation parameter γ₀, the nanoparticle shape and volume fraction ϕ, and the salt concentration and can be as narrow as Bjerrum length (∼nm) or as large as the particle size (∼L the length of charged cylinder).

Dendritic polyelectrolytes with monovalent and divalent counterions: the charge regulation effect and counterion release

The charge regulation and the release of counterions are extremely important and substantial in determining the charge state of polyelectrolytes and the interaction between polyelectrolytes and proteins. Going beyond monovalent to multivalent cations, it is well-known that the effects of ions are qualitatively different. Therefore, the well-accepted descriptions of the charge regulation and the counterion release based on monovalent ions do not immediately apply to systems with multivalent ions. Here, we study the key structural and electrostatic features of charged dendrimers at hand of the pharmaceutically important dendritic polyglycerol sulfate (dPGS) macromolecule equilibrated with monovalent and divalent salts by molecular dynamics (MD) simulations. Following a simple but accurate scheme to determine its effective radius, the counterion condensed layer of the dPGS is determined with high accuracy and we observe the sequential replacement of condensed monovalent cations (MCs) to divalent cations (DCs) rendering a smaller dPGS effective charge versus the DC concentration. We resolve and track the release of counterions on the dPGS along its binding pathway with the plasma protein Human Serum Albumin (HSA). We find that the release of MCs remains favorable for the complexation leading to a considerable amount of release entropy as the driving force for complexation. The release of DCs only occurs above a certain DC concentration with a comparably smaller number of released ions than MCs. Its contribution to the binding free energy is small indicating a subtle cancellation between the entropy gain in releasing DCs and the enthalpy penalty from dissociating DCs from the dendrimer.

Charge regulation of a surface immersed in an electrolyte solution

In this paper, we investigate theoretically a model of charge regulation of a single charged planar surface immersed in an aqueous electrolyte solution. Assuming that the adsorbed ions are mobile in the charged plane, we formulate a field theory of charge regulation where the numbers of adsorbed ions can be determined consistently by equating the chemical potentials of the adsorbed ions to that of the ions in the bulk. We analyze the mean-field treatment of the model for electrolyte of arbitrary valences, and then beyond, where correlation effects are systematically taken into account in a loop expansion. In particular, we compute exactly various one-loop quantities, including electrostatic potentials, ion distributions, and chemical potentials, not only for symmetric (1, 1) electrolyte but also for asymmetric (2, 1) electrolyte, and make use of these quantities to address charge regulation at the one-loop level. We find that correlation effects give rise to various phase transitions in the adsorption of ions, and present phase diagrams for (1, 1) and (2, 1) electrolytes, whose distinct behaviors suggest that charge regulation, at the one-loop level, is no longer universal but depends crucially on the valency of the ions.

Non-monotonic concentration dependence of the electro-phoretic mobility of charged spheres in realistic salt free suspensions

Using super-heterodyne Doppler velocimetry with multiple scattering correction, we extend the optically accessible range of concentrations in experiments on colloidal electro-kinetics. Here, we measured the electro-phoretic mobility and the DC conductivity of aqueous charged sphere suspensions covering about three orders of magnitude in particle concentrations and transmissions as low as 40%. The extended concentration range for the first time allows the demonstration of a non-monotonic concentration dependence of the mobility for a single particle species. Our observations reconcile previous experimental observations made on other species over restricted concentration ranges. We compare our results to the state-of-the-art theoretical calculations using a constant particle charge and the carefully determined experimental boundary conditions as input. In particular, we consider the so-called realistic salt free conditions, i.e., we respect the release of counterions by the particles, the solvent hydrolysis, and the formation of carbonic acid from dissolved neutral CO₂. We also compare our results to previous results obtained under similarly well-defined conditions. This allows identification of three distinct regions of differing density dependence. There is an ascent during the build-up of double layer overlap, which is not expected by theory, an extended plateau region in quantitative agreement with theoretical expectation based on a constant effective charge and a sudden decrease, which occurs way before the expected gradual decrease. Our observations suggest a relation of the non-monotonic behavior to a decrease in particle charge, and we tentatively discuss possibly underlying mechanisms.

Effect of Nanoparticles on the Bulk Shear Viscosity of a Lung Surfactant Fluid

Inhaled nanoparticles (<100 nm) reaching the deep lung region first interact with the pulmonary surfactant, a thin lipid film lining the alveolar epithelium. To date, most biophysical studies have focused on particle-induced modifications of the film interfacial properties. In comparison, there is less work on the surfactant bulk properties and on their changes upon particle exposure. Here we study the viscoelastic properties of a biomimetic pulmonary surfactant in the presence of various engineered nanoparticles. The microrheology technique used is based on the remote actuation of micron-sized wires via the application of a rotating magnetic field and on time-lapse optical microscopy. It is found that particles strongly interacting with lipid vesicles, such as cationic silica (SiO₂, 42 nm) and alumina (Al₂O₃, 40 nm) induce profound modifications of the surfactant flow properties, even at low concentrations. In particular, we find that silica causes fluidification, while alumina induces a liquid-to-soft solid transition. Both phenomena are described quantitatively and accounted for in the context of colloidal physics models. It is finally suggested that the structure and viscosity changes could impair the fluid reorganization and recirculation occurring during breathing.

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.

Investigation of Surfactant AOT Mediated Charging of PS Particles Dispersed in Aqueous Solutions

Nano/submicron particles can be activated by surfactants and aggregate at the air-water interface to generate and stabilize foams. Such systems have been applied extensively in the food, medicine, and cosmetic industries. Studying particle charging behavior in a particle/surfactant/water system is a fundamental way to understand the activation of the particle surface. This paper presents an investigation of the charging behavior of polystyrene (PS) particles dispersed in aqueous solutions of the surfactant sodium di-2-ethylhexylsulfosuccinate (AOT). The results showed that zeta potential of PS was related to the AOT concentration with two different concentration regions. Below the critical micelle concentration (CMC), the charging of PS particles was effected by AOT ions; while above the CMC, it came from both AOT ions and AOT micelles. This behavior was different from that observed for PS in aqueous salt solutions. Additionally, the particle concentration and size were found to affect the zeta potential differently in the two AOT concentration regions. By analyzing these results, the charging mechanism of the PS/AOT/water system was revealed to be preferential adsorption. In summary, the study disclosed the internal connection between the PS charging in aqueous AOT solution and the activation of PS particles, as well as their influence to foam formation and stability.

Transport dynamics of charged colloidal particles during directional drying of suspensions in a confined microchannel

Directional drying of colloidal suspensions, experimentally observed to exhibit mechanical instabilities, is a nonequilibrium procedure that is susceptible to geometric confinement and the properties of colloidal particles. Here, we develop an advection-diffusion model to characterize the transport dynamics for unidirectional drying of a suspension consisting of charged particles in a confined Hele-Shaw cell. We consider the electrostatic interactions by means of the Poisson-Boltzmann cell approach with the viscous flow confined to the cell. By solving the nonequilibrium transport equations, we clarify how the multiple parameters, such as drying rate, confinement ratio, and the monovalent slat concentration, affect the transport dynamics of charged colloidal particles. We find that the drying front recedes into the cell with linear behavior, while the liquid-solid transition front recedes with power law behaviors. The faster evaporation rate creates a rapid formation of the drying front and produces a thinner transition layer. We show that confinement is equivalent to raising the effective concentration in the cell, and, accordingly, the drying front appears earlier and grows more rapidly. Under geometric confinement, a longer fully dried film is created while the total drying time is shortened. Moreover, we have theoretically illustrated that low salt loadings cause a large collective diffusivity of charged colloidal particles, which results in a colloidal network by aggregation. Thus, the drying behavior alters dramatically as salt loadings decrease, since the resulting compacted clusters of charged particles eventually convert the suspension into a gel-like material instead of a simple fluid. Our model is consistent with the current experiments and provides a simple insight for applications in directional solidification and microfluidics.

Jellium and cell model for titratable colloids with continuous size distribution

A good understanding and determination of colloidal interactions is paramount to comprehend and model the thermodynamic and structural properties of colloidal suspensions. In concentrated aqueous suspensions of colloids with a titratable surface charge, this determination is, however, complicated by the density dependence of the effective pair potential due to both the many-body interactions and the charge regulation of the colloids. In addition, colloids generally present a size distribution which results in a virtually infinite combination of colloid pairs. In this paper, we develop two methods and describe the corresponding algorithms to solve this problem for arbitrary size distributions. An implementation in Nim is also provided. The methods, inspired by the seminal work of Torres et al., [J. Chem. Phys. 128, 154906 (2008)] are based on a generalization of the cell and renormalized jellium models to polydisperse suspensions of spherical colloids with a charge regulating boundary condition. The latter is described by the one-pK-Stern model. The predictions of the models are confronted to the equations of state of various commercially available silica dispersions. The renormalized Yukawa parameters (effective charges and screening lengths) are also calculated. The importance of size and charge polydispersity as well as the validity of these two models is discussed in light of the results.

Asymmetric electrolytes near structured dielectric interfaces

The ion distribution of electrolytes near interfaces with dielectric contrast has important consequences for electrochemical processes and many other applications. To date, most studies of such systems have focused on geometrically simple interfaces, for which dielectric effects are analytically solvable or computationally tractable. However, all real surfaces display nontrivial structure at the nanoscale and have, in particular, a nonuniform local curvature. Using a recently developed, highly efficient computational method, we investigate the effect of surface geometry on ion distribution and interface polarization. We consider an asymmetric 2:1 electrolyte bounded by a sinusoidally deformed solid surface. We demonstrate that even when the surface is neutral, the electrolyte acquires a nonuniform ion density profile near the surface. This profile is asymmetric and leads to an effective charging of the surface. We furthermore show that the induced charge is modulated by the local curvature. The effective charge is opposite in sign to the multivalent ions and is larger in concave regions of the surface.

A simple model for electrical charge in globular macromolecules and linear polyelectrolytes in solution

We present a model for calculating the net and effective electrical charge of globular macromolecules and linear polyelectrolytes such as proteins and DNA, given the concentration of monovalent salt and pH in solution. The calculation is based on a numerical solution of the non-linear Poisson-Boltzmann equation using a finite element discretized continuum approach. The model simultaneously addresses the phenomena of charge regulation and renormalization, both of which underpin the electrostatics of biomolecules in solution. We show that while charge regulation addresses the true electrical charge of a molecule arising from the acid-base equilibria of its ionizable groups, charge renormalization finds relevance in the context of a molecule’s interaction with another charged entity. Writing this electrostatic interaction free energy in terms of a local electrical potential, we obtain an “interaction charge” for the molecule which we demonstrate agrees closely with the “effective charge” discussed in charge renormalization and counterion-condensation theories. The predictions of this model agree well with direct high-precision measurements of effective electrical charge of polyelectrolytes such as nucleic acids and disordered proteins in solution, without tunable parameters. Including the effective interior dielectric constant for compactly folded molecules as a tunable parameter, the model captures measurements of effective charge as well as published trends of pKa shifts in globular proteins. Our results suggest a straightforward general framework to model electrostatics in biomolecules in solution. In offering a platform that directly links theory and experiment, these calculations could foster a systematic understanding of the interrelationship between molecular 3D structure and conformation, electrical charge and electrostatic interactions in solution. The model could find particular relevance in situations where molecular crystal structures are not available or rapid, reliable predictions are desired.

Charge and hydration structure of dendritic polyelectrolytes: molecular simulations of polyglycerol sulphate

Macromolecules based on dendritic or hyperbranched polyelectrolytes have been emerging as high potential candidates for biomedical applications. Here we study the charge and solvation structure of dendritic polyglycerol sulphate (dPGS) of generations 0 to 3 in aqueous sodium chloride solution by explicit-solvent molecular dynamics computer simulations. We characterize dPGS by calculating several important properties such as relevant dPGS radii, molecular distributions, the solvent accessible surface area, and the partial molecular volume. In particular, as the dPGS exhibits high charge renormalization effects, we address the challenges of how to obtain a well-defined effective charge and surface potential of dPGS for practical applications. We compare implicit- and explicit-solvent approaches in our all-atom simulations with the coarse-grained simulations from our previous work. We find consistent values for the effective electrostatic size (i.e., the location of the effective charge of a Debye-Hückel sphere) within all the approaches, deviating at most by the size of a water molecule. Finally, the excess chemical potential of water insertion into dPGS and its thermodynamic signature are presented and rationalized.

Surface and extrapolated point charge renormalizations for charge-stabilized colloidal spheres

The Derjaguin-Landau-Verwey-Overbeek (DLVO) theory is widely used to model interactions between weakly charged spheres in dilute suspensions. For particles bearing a higher charge, the linearized electrostatics underlying the DLVO theory is no longer valid but it is possible to map the real colloidal system to an auxiliary one that still obeys linear electrostatics but which involves a different, effective pair potential. This procedure, termed renormalization, can be performed in various ways, the most widely used being surface charge renormalization (SCR) based on the cell model. SCR is still limited to dilute suspensions since the auxiliary system is made of spheres interacting through a DLVO-like pair potential. The recent extrapolated point charge (EPC) renormalization overcomes this limitation by using point charges in the auxiliary system and has indeed been shown to produce better results than the SCR in dense suspensions. Here, we recall that the DLVO-like potential used in the SCR can be modified to account for many-body ion-colloid core exclusion effects (a model termed SCRX here); we show that the accuracy of the EPC and SCRX renormalizations is virtually identical, and conclude by explaining why the EPC method is still the most attractive option of the two in many cases.

Direct determination of forces between charged nanogels through coarse-grained simulations

In this work, electrostatic forces between charged nanogels are explored through coarse-grained simulations. These simulations allow us to explicitly consider the complex topology of these nanoparticles and provide reliable force values to examine highly charged nanogels of a few tens of nanometers. The results obtained here clearly reveal that the electrostatic interactions between these nanoparticles are not governed by the net charge of the nanogel, which includes not only the charge of the polymer network but also the charge of ions inside. Thus two theoretical procedures for predicting effective charges are also proposed and investigated. Both provide predictions of the same order and capture the behavior found for the effective charge obtained from simulations.

Ability of the Poisson-Boltzmann equation to capture molecular dynamics predicted ion distribution around polyelectrolytes

Here, we examine polyelectrolyte (PE) and ion chemistry specificity in ion condensation via all-atom molecular dynamics (MD) simulations and assess the ability of the Poisson-Boltzmann (PB) equation to describe the ion distribution predicted by the MD simulations. The PB model enables the extraction of parameters characterizing ion condensation. We find that the modified PB equation which contains the effective PE radius and the energy of the ion-specific interaction as empirical fitting parameters describes ion distribution accurately at large distances but close to the PE, especially when strongly localized charge or specific ion binding sites are present, the mean field description of PB fails. However, the PB model captures the MD predicted ion condensation in terms of the Manning radius and fraction of condensed counterions for all the examined PEs and ion species. We show that the condensed ion layer thickness in our MD simulations collapses on a single master curve for all the examined simple, monovalent ions (Na⁺, Br⁺, Cs⁺, Cl^-, and Br^-) and PEs when plotted against the Manning parameter (and consequently the PE line charge density). The significance of this finding is that, contrary to the Manning radius extracted from the mean field PB model, the condensed layer thickness in the all atom detail MD modelling does not depend on the PE chemistry or counterion type. Furthermore, the fraction of condensed counterions in the MD simulations exceeds the PB theory prediction. The findings contribute toward understanding and modelling ion distribution around PEs and other charged macromolecules in aqueous solutions, such as DNA, functionalized nanotubes, and viruses.

Fast, Robust Evaluation of the Equation of State of Suspensions of Charge-Stabilized Colloidal Spheres

Increasing demand is appearing for the fast, robust prediction of the equation of state of colloidal suspensions, notably with a view to using it as input data to calculate transport coefficients in complex flow solvers. This is also of interest in rheological studies, industrial screening tests of new formulations, and the real-time interpretation of osmotic compression experiments, for example. For charge-stabilized spherical particles, the osmotic pressure can be computed with standard liquid theories. However, this calculation can sometimes be lengthy and/or unstable under some physicochemical conditions, a drawback that precludes its use in multiscale flow simulators. As a simple, fast, and robust replacement, the literature reports estimations of the osmotic pressure that have been built by adding the Carnahan-Starling and the cell model pressures (CSCM model). The first contribution is intended to account for colloid-colloid contacts, and the second, for electrostatic effects. This approximation has not yet been thoroughly tested. In this work, the CSCM is evaluated by comparison with data from experiments on silica particles, Monte Carlo simulations, and solutions of the accurate Rogers-Young integral equation scheme with a hard-sphere Yukawa potential obtained from the extrapolated point-charge renormalization method for a wide range of volume fractions, surface charge densities, and interaction ranges. We find that the CSCM is indeed perfectly adequate in the electrostatically concentrated regime, where it can be used from vanishingly small to high surface charge because there is error cancellation between the Carnahan-Starling and cell model contributions at intermediate charge. The CSCM is thus a nice extension of the cell model to liquid-like dense suspensions, which should find application in the domains mentioned above. However, it fails for dilute suspensions with strong electrostatics. In this case, we show that, and explain why, perturbation methods and the rescaled mean spherical approximation are good alternatives in terms of precision, ease of implementation, computational cost, and robustness.

Adsorption of Milk Proteins (β-Casein and β-Lactoglobulin) and BSA onto Hydrophobic Surfaces

Here, we study films of proteins over planar surfaces and protein-coated microspheres obtained from the adsorption of three different proteins ( β -casein, β -lactoglobulin and bovine serum albumin (BSA)). The investigation of protein films in planar surfaces is performed by combining quartz crystal microbalance (QCM) and atomic force microscopy (AFM) measurements with all-atomic molecular dynamics (MD) simulations. We found that BSA and β -lactoglobulin form compact monolayers, almost without interstices between the proteins. However, β -casein adsorbs forming multilayers. The study of the electrokinetic mobility of protein-coated latex microspheres shows substantial condensation of ions from the buffer over the complexes, as predicted from ion condensation theories. The electrokinetic behavior of the latex-protein complexes is dominated by the charge of the proteins and the phenomenon of ion condensation, whereas the charge of the latex colloids plays only a minor role.

Drying paint: from micro-scale dynamics to mechanical instabilities

Charged colloidal dispersions make up the basis of a broad range of industrial and commercial products, from paints to coatings and additives in cosmetics. During drying, an initially liquid dispersion of such particles is slowly concentrated into a solid, displaying a range of mechanical instabilities in response to highly variable internal pressures. Here we summarize the current appreciation of this process by pairing an advection-diffusion model of particle motion with a Poisson-Boltzmann cell model of inter-particle interactions, to predict the concentration gradients in a drying colloidal film. We then test these predictions with osmotic compression experiments on colloidal silica, and small-angle X-ray scattering experiments on silica dispersions drying in Hele-Shaw cells. Finally, we use the details of the microscopic physics at play in these dispersions to explore how two macroscopic mechanical instabilities-shear-banding and fracture-can be controlled.This article is part of the themed issue 'Patterning through instabilities in complex media: theory and applications.'

Charged plate in asymmetric electrolytes: One-loop renormalization of surface charge density and Debye length due to ionic correlations

Self-consistent field theory (SCFT) is used to study the mean potential near a charged plate inside a m:-n electrolyte. A perturbation series is developed in terms of g=4πκb, where band1/κ are Bjerrum length and bare Debye length, respectively. To the zeroth order, we obtain the nonlinear Poisson-Boltzmann theory. For asymmetric electrolytes (m≠n), the first order (one-loop) correction to mean potential contains a secular term, which indicates the breakdown of the regular perturbation method. Using a renormalizaton group transformation, we remove the secular term and obtain a globally well-behaved one-loop approximation with a renormalized Debye length and a renormalized surface charge density. Furthermore, we find that if the counterions are multivalent, the surface charge density is renormalized substantially downwards and may undergo a change of sign, if the bare surface charge density is sufficiently large. Our results agrees with large MC simulation even when the density of electrolytes is relatively high.

Modeling the Electrostatics of Hollow Shell Suspensions: Ion Distribution, Pair Interactions, and Many-Body Effects

Electrostatic interactions play a key role in hollow shell suspensions as they determine their structure, stability, thermodynamics, and rheology and also the loading capacity of small charged species for nanoreservoir applications. In this work, fast, reliable modeling strategies aimed at predicting the electrostatics of hollow shells for one, two, and many colloids are proposed and validated. The electrostatic potential inside and outside a hollow shell with a finite thickness and a specific permittivity is determined analytically in the Debye-Hückel (DH) limit. An expression for the interaction potential between two such hollow shells is then derived and validated numerically. It follows a classical Yukawa form with an effective charge depending on the shell geometry, permittivity, and inner and outer surface charge densities. The predictions of the Ornstein-Zernike (OZ) equation with this pair potential to determine equations of state are then evaluated by comparison to results obtained with a Brownian dynamics algorithm coupled to the resolution of the linearized Poisson-Boltzmann and Laplace equations (PB-BD simulations). The OZ equation based on the DLVO-like potential performs very well in the dilute regime as expected, but also quite well, and more surprisingly, in the concentrated regime in which full spheres exhibit significant many-body effects. These effects are shown to vanish for shells with small thickness and high permittivity. For highly charged hollow shells, we propose and validate a charge renormalization procedure. Finally, using PB-BD simulations, we show that the cell model predicts the ion distribution inside and outside hollow shells accurately in both electrostatically dilute and concentrated suspensions. We then determine the shell loading capacity as a function of salt concentration, volume fraction, and surface charge density for nanoreservoir applications such as drug delivery, sensing, or smart coatings.

Osmotic Bulk Modulus of Charged Colloids Measured by Ensemble Optical Trapping

The optical-bottle technique is used to measure osmotic bulk moduli of colloid suspensions. The bulk modulus is determined by optically trapping an ensemble of nanoparticles and invoking a steady-state force balance between confining optical-gradient forces and repulsive osmotic-pressure forces. Osmotic bulk moduli are reported for aqueous suspensions of charged polystyrene particles in NaCl solutions as a function of particle concentration and ionic strength, and are compared to those determined by turbidity measurements under the same conditions. Effective particle charges are calculated from the bulk moduli and are found to increase as a function of ionic strength, consistent with previously reported results.

Comparison of different coupling schemes between counterions and charged nanoparticles in multiparticle collision dynamics

We applied the multiparticle collision dynamics (MPC) simulation technique to highly asymmetric electrolytes in solution, i.e., charged nanoparticles and their counterions in a solvent. These systems belong to a domain of solute size which ranges between the electrolyte and the colloidal domains, where most analytical theories are expected to fail, and efficient simulation techniques are still missing. MPC is a mesoscopic simulation method which mimics hydrodynamics properties of a fluid, includes thermal fluctuations, and can be coupled to a molecular dynamics of solutes. We took advantage of the size asymmetry between nanoparticles and counterions to treat the coupling between solutes and the solvent bath within the MPC method. Counterions were coupled to the solvent bath during the collision step and nanoparticles either through a direct interaction force or with stochastic rotation rules which mimic stick boundary conditions. Moreover, we adapted the simulation procedure to address the issue of the strong electrostatic interactions between solutes of opposite charges. We show that the short-ranged repulsion between counterions and nanoparticles can be modeled by stochastic reflection rules. This simulation scheme is very efficient from a computational point of view. We have also computed the transport coefficients for various densities. The diffusion of counterions was found in one case to increase slightly with the volume fraction of nanoparticles. The deviation of the electric conductivity from the ideal behavior (solutes at infinite dilution without any direct interactions) is found to be strong.

Nanoparticle Clusters: Assembly and Control Over Internal Order, Current Capabilities, and Future Potential

The current state of the art in the use of colloidal methods to form nanoparticle assemblies, or clusters (NPCs) is reviewed. The focus is on the two-step approach, which exploits the advantages of bottom-up wet chemical NP synthesis procedures, with subsequent colloidal destabilization to trigger assembly in a controlled manner. Recent successes in the application of functional NPCs with enhanced emergent collective properties for a wide range of applications, including in biomedical detection, surface enhanced Raman scattering (SERS) enhancement, photocatalysis, and light harvesting, are highlighted. The role of the NP-NP interactions in the formation of monodisperse ordered clusters is described and the different assembly processes from a wide range of literature sources are classified according to the nature of the perturbation from the initial equilibrium state (dispersed NPs). Finally, the future for the field and the anticipated role of computational approaches in developing next-generation functional NPCs are briefly discussed.

Hiding in Plain View: Colloidal Self-Assembly from Polydisperse Populations

We report small-angle x-ray scattering experiments on aqueous dispersions of colloidal silica with a broad monomodal size distribution (polydispersity, 14%; size, 8 nm). Over a range of volume fractions, the silica particles segregate to build first one, then two distinct sets of colloidal crystals. These dispersions thus demonstrate fractional crystallization and multiple-phase (bcc, Laves AB_{2}, liquid) coexistence. Their remarkable ability to build complex crystal structures from a polydisperse population originates from the intermediate-range nature of interparticle forces, and it suggests routes for designing self-assembling colloidal crystals from the bottom up.

Like-charged protein-polyelectrolyte complexation driven by charge patches

We study the pair complexation of a single, highly charged polyelectrolyte (PE) chain (of 25 or 50 monomers) with like-charged patchy protein models (CPPMs) by means of implicit-solvent, explicit-salt Langevin dynamics computer simulations. Our previously introduced set of CPPMs embraces well-defined zero-, one-, and two-patched spherical globules each of the same net charge and (nanometer) size with mono- and multipole moments comparable to those of globular proteins with similar size. We observe large binding affinities between the CPPM and the like-charged PE in the tens of the thermal energy, kBT, that are favored by decreasing salt concentration and increasing charge of the patch(es). Our systematic analysis shows a clear correlation between the distance-resolved potentials of mean force, the number of ions released from the PE, and CPPM orientation effects. In particular, we find a novel two-site binding behavior for PEs in the case of two-patched CPPMs, where intermediate metastable complex structures are formed. In order to describe the salt-dependence of the binding affinity for mainly dipolar (one-patched) CPPMs, we introduce a combined counterion-release/Debye-Hückel model that quantitatively captures the essential physics of electrostatic complexation in our systems.

Poisson-Boltzmann thermodynamics of counterions confined by curved hard walls

We consider a set of identical mobile pointlike charges (counterions) confined to a domain with curved hard walls carrying a uniform fixed surface charge density, the system as a whole being electroneutral. Three domain geometries are considered: a pair of parallel plates, the cylinder, and the sphere. The particle system in thermal equilibrium is assumed to be described by the nonlinear Poisson-Boltzmann theory. While the effectively one-dimensional plates and the two-dimensional cylinder have already been solved, the three-dimensional sphere problem is not integrable. It is shown that the contact density of particles at the charged surface is determined by a first-order Abel differential equation of the second kind which is a counterpart of Enig's equation in the critical theory of gravitation and combustion or explosion. This equation enables us to construct the exact series solutions of the contact density in the regions of small and large surface charge densities. The formalism provides, within the mean-field Poisson-Boltzmann framework, the complete thermodynamics of counterions inside a charged sphere (salt-free system).

General Differential Contact Identities for Macromolecules

We discuss general Maxwell identities relating a macromolecule's charge, the forces acting at its surface, and the osmotic pressure of the solution in which it sits. The identities are closely related to the contact value relations that hold for certain special geometries, but are more general. In particular, the Maxwell identities can be applied to any macromolecule geometry, and they hold both within and outside of mean-field theory. Examples illustrate that combining the identities with approximate treatments of screening can often return simple, accurate osmotic pressure estimates.

Principles for the determination of the surface potential of charged particles in mixed electrolyte solutions

The Gouy–Chapman surface potential is a key parameter for many interfacial phenomena in physical, chemical and biological systems. Existing theoretical approaches allow the determination of the surface potential at a solid–liquid interface only in single electrolyte solutions; however, mixed electrolytes are often encountered in practical applications. The development of a theoretical approach for the determination of the surface potential in mixed electrolyte solutions is therefore a desirable goal. In this study, this important issue was resolved for the first time. Based on the analytical solutions of the nonlinear Poisson–Boltzmann equation in different mixed electrolyte solutions, corresponding mathematical relationships were developed between the surface potential and the mean ionic concentration in the diffuse layer. As the mean ionic concentration in the diffuse layer can be easily determined, the surface potential could be calculated using the newly derived equations. The effects of electrolyte composition on the surface potential were theoretically quantified in the new equations, while only counterionic type was taken into account for mixed electrolyte solutions in the current studies.