Freezing lines of colloidal Yukawa spheres. II. Local structure and characteristic lengths

Quasi-two-dimensional dispersions of Brownian particles with competitive interactions: phase behavior and structural properties

Competing short-range attractive (SA) and long range repulsive (LR) particle interactions can be used to describe three-dimensional charge-stabilized colloid or protein dispersions at low added salt concentrations, as well as membrane proteins with interaction contributions mediated by lipid molecules. Using Langevin dynamics (LD) simulations, we determine the generalized phase diagram, cluster shapes and size distributions of a generic quasi-two-dimensional (Q2D) dispersion of spherical SALR particles confined to in-plane motion inside a bulk fluid. The SA and LR interaction parts are modelled by a generalized Lennard-Jones potential and a screened Coulomb potential, respectively. The microstructures of the detected equilibrium and non-equilibrium Q2D phases are distinctly different from those observed in three-dimensional (3D) SALR systems, by exhibiting different levels of hexagonal ordering. We discuss a thermodynamic perturbation theory prediction for the metastable binodal line of a reference system of particles with SA interactions only, which in the explored Q2D-SALR phase diagram region separates cluster from non-clustered phases. The transition from the high-temperature (small SA) dispersed fluid (DF) phase to the lower-temperature equilibrium cluster (EC) fluid phase is characterised by a low-wavenumber peak height of the static structure factor (corresponding to a thermal correlation length of about twice the particle diameter) featuring a distinctly smaller value (≈1.4) than in 3D SALR systems. With decreasing temperature (increasing SA), the cluster morphology changes from disk-like shapes in the equilibrium cluster phase, to double-stranded anisotropic hexagonal cluster segments formed in a cluster-percolated (CP) gel-like phase. This transition can be quantified by a hexagonal order parameter distribution function. The mean cluster size and coordination number of particles in the CP phase are insensitive to changes in the attraction strength.

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.

Short- and long-time diffusion and dynamic scaling in suspensions of charged colloidal particles

We report on a comprehensive theory-simulation-experimental study of collective and self-diffusion in concentrated suspensions of charge-stabilized colloidal spheres. In theory and simulation, the spheres are assumed to interact directly by a hard-core plus screened Coulomb effective pair potential. The intermediate scattering function, f_c(q, t), is calculated by elaborate accelerated Stokesian dynamics (ASD) simulations for Brownian systems where many-particle hydrodynamic interactions (HIs) are fully accounted for, using a novel extrapolation scheme to a macroscopically large system size valid for all correlation times. The study spans the correlation time range from the colloidal short-time to the long-time regime. Additionally, Brownian Dynamics (BD) simulation and mode-coupling theory (MCT) results of f_c(q, t) are generated where HIs are neglected. Using these results, the influence of HIs on collective and self-diffusion and the accuracy of the MCT method are quantified. It is shown that HIs enhance collective and self-diffusion at intermediate and long times. At short times self-diffusion, and for wavenumbers outside the structure factor peak region also collective diffusion, are slowed down by HIs. MCT significantly overestimates the slowing influence of dynamic particle caging. The dynamic scattering functions obtained in the ASD simulations are in overall good agreement with our dynamic light scattering (DLS) results for a concentration series of charged silica spheres in an organic solvent mixture, in the experimental time window and wavenumber range. From the simulation data for the time derivative of the width function associated with f_c(q, t), there is indication of long-time exponential decay of f_c(q, t), for wavenumbers around the location of the static structure factor principal peak. The experimental scattering functions in the probed time range are consistent with a time-wavenumber factorization scaling behavior of f_c(q, t) that was first reported by Segrè and Pusey [Phys. Rev. Lett. 77, 771 (1996)] for suspensions of hard spheres. Our BD simulation and MCT results predict a significant violation of exact factorization scaling which, however, is approximately restored according to the ASD results when HIs are accounted for, consistent with the experimental findings for f_c(q, t). Our study of collective diffusion is amended by simulation and theoretical results for the self-intermediate scattering function, f_s(q, t), and its non-Gaussian parameter α₂(t) and for the particle mean squared displacement W(t) and its time derivative. Since self-diffusion properties are not assessed in standard DLS measurements, a method to deduce W(t) approximately from f_c(q, t) is theoretically validated.

A combined 3D and 2D light scattering study on aqueous colloidal model systems with tunable interactions

In this article we report on the synthesis and characterization of a system of colloidal spheres suspended in an aqueous solvent which can be refractive index-matched, thus allowing for investigations of the particle near-wall dynamics by evanescent wave dynamic light scattering at concentrations up to the isotropic to ordered transition and beyond. The particles are synthesized by copolymerization of a fluorinated acrylic ester monomer with a polyethylene-glycol (PEG) oligomer by surfactant free emulsion polymerization. Static and dynamic light scattering experiments in combination with cryo transmission electron microscopy reveal that the particles have a core shell structure with a significant enrichment of the PEG chains on the particles surface. In index-matching DMSO/water suspensions the particles arrange in an ordered phase at volume fraction above 7%, if no additional electrolyte is present. The near-wall dynamics at low volume fraction are quantitatively described by the combination of electrostatic repulsion and hydrodynamic interaction between the particles and the wall. At volume fractions close to the isotropic to ordered transition, the near-wall dynamics are more complex and qualitatively reminiscent of the behaviour which was observed in hard sphere suspensions at high concentrations.

Ultrafiltration of charge-stabilized dispersions at low salinity

We present a comprehensive study of cross-flow ultrafiltration (UF) of charge-stabilized suspensions, under low-salinity conditions of electrostatically strongly repelling colloidal particles. The axially varying permeate flux, near-membrane concentration-polarization (CP) layer and osmotic pressure profiles are calculated using a macroscopic diffusion-advection boundary layer method, and are compared with filtration experiments on aqueous suspensions of charge-stabilized silica particles. The theoretical description based on the one-component macroion fluid model (OCM) accounts for the strong influence of surface-released counterions on the renormalized colloid charge and suspension osmotic compressibility, and for the influence of the colloidal hydrodynamic interactions and electric double layer repulsion on the concentration-dependent suspension viscosity η, and collective diffusion coefficient Dc. A strong electro-hydrodynamic enhancement of Dc and η, and likewise of the osmotic pressure, is predicted theoretically, as compared with their values for a hard-sphere suspension. We also point to the failure of generalized Stokes-Einstein relations describing reciprocal relations between Dc and η. According to our filtration model, Dc is of dominant influence, giving rise to an only weakly developed CP layer having practically no effect on the permeate flux. This prediction is quantitatively confirmed by our UF measurements of the permeate flux using an aqueous suspension of charged silica spheres as the feed system. The experimentally detected fouling for the largest considered transmembrane pressure values is shown not to be due to filter cake formation by crystallization or vitrification.

Rotational self-diffusion in suspensions of charged particles: simulations and revised Beenakker-Mazur and pairwise additivity methods

We present a comprehensive joint theory-simulation study of rotational self-diffusion in suspensions of charged particles whose interactions are modeled by the generic hard-sphere plus repulsive Yukawa (HSY) pair potential. Elaborate, high-precision simulation results for the short-time rotational self-diffusion coefficient, D(r), are discussed covering a broad range of fluid-phase state points in the HSY model phase diagram. The salient trends in the behavior of D(r) as a function of reduced potential strength and range, and particle concentration, are systematically explored and physically explained. The simulation results are further used to assess the performance of two semi-analytic theoretical methods for calculating D(r). The first theoretical method is a revised version of the classical Beenakker-Mazur method (BM) adapted to rotational diffusion which includes a highly improved treatment of the salient many-particle hydrodynamic interactions. The second method is an easy-to-implement pairwise additivity (PA) method in which the hydrodynamic interactions are treated on a full two-body level with lubrication corrections included. The static pair correlation functions required as the only input to both theoretical methods are calculated using the accurate Rogers-Young integral equation scheme. While the revised BM method reproduces the general trends of the simulation results, it significantly underestimates D(r). In contrast, the PA method agrees well with the simulation results for D(r) even for intermediately concentrated systems. A simple improvement of the PA method is presented which is applicable for large concentrations.