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.

Geometrical Influence on Particle Transport in Cross-Flow Ultrafiltration: Cylindrical and Flat Sheet Membranes.

Cross-flow membrane ultrafiltration (UF) is used for the enrichment and purification of small colloidal particles and proteins. We explore the influence of different membrane geometries on the particle transport in, and the efficiency of, inside-out cross-flow UF. For this purpose, we generalize the accurate and numerically efficient modified boundary layer approximation (mBLA) method, developed in recent work by us for a hollow cylindrical membrane, to parallel flat sheet geometries with one or two solvent-permeable membrane sheets. Considering a reference dispersion of Brownian hard spheres where accurate expressions for its transport properties are available, the generalized mBLA method is used to analyze how particle transport and global UF process indicators are affected by varying operating parameters and the membrane geometry. We show that global process indicators including the mean permeate flux, the solvent recovery indicator, and the concentration factor are strongly dependent on the membrane geometry. A key finding is that irrespective of the many input parameters characterizing an UF experiment and its membrane geometry, the process indicators are determined by three independent dimensionless variables only. This finding can be very useful in the design, optimization, and scale-up of UF processes.

Hydrodynamics of immiscible binary fluids with viscosity contrast: a multiparticle collision dynamics approach.

We present a multiparticle collision dynamics (MPC) implementation of layered immiscible fluids A and B of different shear viscosities separated by planar interfaces. The simulated flow profile for imposed steady shear motion and the time-dependent shear stress functions are in excellent agreement with our continuum hydrodynamics results for the composite fluid. The wave-vector dependent transverse velocity auto-correlation functions (TVAF) in the bulk-fluid regions of the layers decay exponentially, and agree with those of single-phase isotropic MPC fluids. In addition, we determine the hydrodynamic mobilities of an embedded colloidal sphere moving steadily parallel or transverse to a fluid-fluid interface, as functions of the distance from the interface. The obtained mobilities are in good agreement with hydrodynamic force multipoles calculations, for a no-slip sphere moving under creeping flow conditions near a clean, ideally flat interface. The proposed MPC fluid-layer model can be straightforwardly implemented, and it is computationally very efficient. Yet, owing to the spatial discretization inherent to the MPC method, the model can not reproduce all hydrodynamic features of an ideally flat interface between immiscible fluids.

Modeling cross-flow ultrafiltration of permeable particle dispersions.

Cross-flow ultrafiltration is a pressure-driven separation and enrichment process of small colloidal particles where a colloidal feed dispersion is continuously pumped through a membrane pipe permeable to the solvent only. We present a semi-analytic modified boundary layer approximation (mBLA) method for calculating the inhomogeneous concentration-polarization (CP) layer of particles near the membrane and the dispersion flow in a cross-flow filtration setup with a hollow fiber membrane. Conditions are established for which unwarranted axial flow and permeate flow reversal are excluded, and non-monotonic CP profiles are observed. The permeate flux is linked to the particle concentration on the membrane wall using the Darcy-Starling expression invoking axially varying osmotic and trans-membrane pressures. Results are discussed for dispersions of hard spheres serving as a reference system and for solvent-permeable particles mimicking non-ionic microgels. Accurate analytic expressions are employed for the concentration and solvent permeability dependent dispersion viscosity and gradient diffusion coefficient entering into the effective Stokes flow and advection-diffusion equations. We show that the mBLA concentration and flow profiles are in quantitative agreement with results by a finite element method. The mBLA results are compared with predictions by an earlier CP layer similarity solution, showing the higher precision of the former method.

Modeling deswelling, thermodynamics, structure, and dynamics in ionic microgel suspensions.

Ionic microgel particles in a good solvent swell to an equilibrium size determined by a balance of electrostatic and elastic forces. When crowded, ionic microgels deswell owing to a redistribution of microions inside and outside the particles. The concentration-dependent deswelling affects the interactions between the microgels and, consequently, the suspension properties. We present a comprehensive theoretical study of crowding effects on thermodynamic, structural, and dynamic properties of weakly cross-linked ionic microgels in a good solvent. The microgels are modeled as microion- and solvent-permeable colloidal spheres with fixed charge uniformly distributed over the polymer gel backbone, whose elastic and solvent-interaction free energies are described using the Flory-Rehner theory. Two mean-field methods for calculating the crowding-dependent microgel radius are investigated and combined with calculations of the net microgel charge characterizing the electrostatic part of an effective microgel pair potential, with charge renormalization accounted for. Using this effective pair potential, thermodynamic and static suspension properties are calculated, including the osmotic pressure and microgel pair distribution function. The latter is used in our calculations of dynamic suspension properties, where we account for hydrodynamic interactions. Results for diffusion and rheological properties are presented over ranges of microgel concentration and charge. We show that deswelling mildly enhances self- diffusion and collective diffusion and the osmotic pressure, lowers the suspension viscosity, and significantly shifts the suspension crystallization point to higher concentrations. This paper presents a bottom-up approach to efficiently computing suspension properties of crowded ionic microgels using single-particle characteristics.

Structure and dynamics in suspensions of soft core-shell colloids in the fluid regime.

We report on a detailed experimental study of the structure and short-time dynamics in fluid-regime suspensions of soft core-shell spherical particles with different molecular weights of the chains forming the soft outer shell, and therefore different degrees of particle softness, using 3D dynamic light scattering (3D-DLS). Owing to the particle softness, the liquid-crystal coexistence regime is found to be broader than that of hard-sphere (HS) suspensions. Static light scattering in the dilute regime yields form factors that can be described using a spherical core-shell model and second virial coefficients A2 > 0 indicative of purely repulsive interactions. The particle-particle interactions are longer ranged for all considered systems except those of the smaller molecular weight chain grafted particles which show a HS-like behavior. 3D-DLS experiments in the concentrated regime up to the liquid-crystal transition provide the short-time diffusion function, D(q), in a broad range of scattering wavenumbers, q, from which the structural (cage) and short-time self-diffusion coefficients D(qm) and DS = D(q ≫ qm), respectively, are deduced as functions of the effective particle volume fraction, ϕ = c/c*, where c* is the overlap concentration, calculated using the hydrodynamic particle radius, RH. The size of the nearest-neighbor cage of particles is characterized by 2π/qm, with D(q) and the static structure factor S(q) attaining at qm the smallest and largest values, respectively. Experimental data of D(qm) and DS are contrasted with analytic theoretical predictions based on a simplifying hydrodynamic radius model where the internal hydrodynamic structure of the core-shell particles is mapped on a single hydrodynamic radius parameter γ = RH/Reff, for constant direct interactions characterized by an (effective) hard-core radius Reff. The particle softness is reflected, in particular, in the corresponding shape of the static structure factor, while the mean solvent (Darcy) permeability of the particles related to γ is reflected in the dynamic properties only. For grafted particles with longer polymer chains, D(qm) and DS are indicative of larger permeability values while particles with shorter chains are practically nonpermeable. The particle softness is also evident in the effective random close packing fraction estimated from the extrapolated zero-value limit of the cage diffusion coefficient D(qm).

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, fc(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 fc(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 fc(q, t), there is indication of long-time exponential decay of fc(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 fc(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 fc(q, t). Our study of collective diffusion is amended by simulation and theoretical results for the self-intermediate scattering function, fs(q, t), and its non-Gaussian parameter α2(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 fc(q, t) is theoretically validated.

Short-time dynamics of lysozyme solutions with competing short-range attraction and long-range repulsion: Experiment and theory.

Recently, atypical static features of microstructural ordering in low-salinity lysozyme protein solutions have been extensively explored experimentally and explained theoretically based on a short-range attractive plus long-range repulsive (SALR) interaction potential. However, the protein dynamics and the relationship to the atypical SALR structure remain to be demonstrated. Here, the applicability of semi-analytic theoretical methods predicting diffusion properties and viscosity in isotropic particle suspensions to low-salinity lysozyme protein solutions is tested. Using the interaction potential parameters previously obtained from static structure factor measurements, our results of Monte Carlo simulations representing seven experimental lysoyzme samples indicate that they exist either in dispersed fluid or random percolated states. The self-consistent Zerah-Hansen scheme is used to describe the static structure factor, S(q), which is the input to our calculation schemes for the short-time hydrodynamic function, H(q), and the zero-frequency viscosity η. The schemes account for hydrodynamic interactions included on an approximate level. Theoretical predictions for H(q) as a function of the wavenumber q quantitatively agree with experimental results at small protein concentrations obtained using neutron spin echo measurements. At higher concentrations, qualitative agreement is preserved although the calculated hydrodynamic functions are overestimated. We attribute the differences for higher concentrations and lower temperatures to translational-rotational diffusion coupling induced by the shape and interaction anisotropy of particles and clusters, patchiness of the lysozyme particle surfaces, and the intra-cluster dynamics, features not included in our simple globular particle model. The theoretical results for the solution viscosity, η, are in qualitative agreement with our experimental data even at higher concentrations. We demonstrate that semi-quantitative predictions of diffusion properties and viscosity of solutions of globular proteins are possible given only the equilibrium structure factor of proteins. Furthermore, we explore the effects of changing the attraction strength on H(q) and η.

Clustering and dynamics of particles in dispersions with competing interactions: theory and simulation.

Dispersions of particles with short-range attractive and long-range repulsive interactions exhibit rich equilibrium microstructures and a complex phase behavior. We present theoretical and simulation results for structural and, in particular, short-time diffusion properties of a colloidal model system with such interactions, both in the dispersed-fluid and equilibrium-cluster phase regions. The particle interactions are described by a generalized Lennard-Jones-Yukawa pair potential. For the theoretical-analytical description, we apply the hybrid Beenakker-Mazur pairwise additivity (BM-PA) scheme. The static structure factor input to this scheme is calculated self-consistently using the Zerah-Hansen integral equation theory approach. In the simulations, a hybrid simulation method is adopted, combing molecular dynamics simulations of colloids with the multiparticle collision dynamics approach for the fluid, which fully captures hydrodynamic interactions. The comparison of our theoretical and simulation results confirms the high accuracy of the BM-PA scheme for dispersed-fluid phase systems. For particle attraction strengths exceeding a critical value, our simulations yield an equilibrium cluster phase. Calculations of the mean lifetime of the appearing clusters and the comparison with the analytical prediction of the dissociation time of an isolated particle pair reveal quantitative differences pointing to the importance of many-particle hydrodynamic interactions for the cluster dynamics. The cluster lifetime in the equilibrium-cluster phase increases far stronger with increasing attraction strength than that in the dispersed-fluid phase. Moreover, significant changes in the cluster shapes are observed in the course of time. Hence, an equilibrium-cluster dispersion cannot be treated dynamically as a system of permanent rigid bodies.

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.

Short-time dynamics in dispersions with competing short-range attraction and long-range repulsion.

Dynamic clustering of globular Brownian particles in dispersions exhibiting competing short-range attraction and long-range repulsion (SALR) such as low-salinity protein solutions has gained a lot of interest over the past few years. While the structure of the various cluster phases has been intensely explored, little is known about the dynamics of SALR systems. We present the first systematic theoretical study of short-time diffusion and rheological transport properties of two-Yukawa potential SALR systems in the single-particle dominated dispersed-fluid phase, using semi-analytic methods where the salient hydrodynamic interactions are accounted for. We show that the dynamics has unusual features compared to reference systems with pure repulsion or attraction. Results are discussed for the hydrodynamic function characterizing short-time diffusion that reveals an intermediate-range-order (cluster) peak, self-diffusion and sedimentation coefficients, and high-frequency viscosity. As important applications, we discuss the applicability of two generalized Stokes-Einstein relations, and assess the wavenumber range required for the determination of self-diffusion in a dynamic scattering experiment.

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.

Ultrafiltration modeling of non-ionic microgels.

Membrane ultrafiltration (UF) is a pressure driven process allowing for the separation and enrichment of protein solutions and dispersions of nanosized microgel particles. The permeate flux and the near-membrane concentration-polarization (CP) layer in this process is determined by advective-diffusive dispersion transport and the interplay of applied and osmotic transmembrane pressure contributions. The UF performance is thus strongly dependent on the membrane properties, the hydrodynamic structure of the Brownian particles, their direct and hydrodynamic interactions, and the boundary conditions. We present a macroscopic description of cross-flow UF of non-ionic microgels modeled as solvent-permeable spheres. Our filtration model involves recently derived semi-analytic expressions for the concentration-dependent collective diffusion coefficient and viscosity of permeable particle dispersions [Riest et al., Soft Matter, 2015, 11, 2821]. These expressions have been well tested against computer simulation and experimental results. We analyze the CP layer properties and the permeate flux at different operating conditions and discuss various filtration process efficiency and cost indicators. Our results show that the proper specification of the concentration-dependent transport coefficients is important for reliable filtration process predictions. We also show that the solvent permeability of microgels is an essential ingredient to the UF modeling. The particle permeability lowers the particle concentration at the membrane surface, thus increasing the permeate flux.

Dynamics of suspensions of hydrodynamically structured particles: analytic theory and applications to experiments.

We present an easy-to-use analytic toolbox for the calculation of short-time transport properties of concentrated suspensions of spherical colloidal particles with internal hydrodynamic structure, and direct interactions described by a hard-core or soft Hertz pair potential. The considered dynamic properties include self-diffusion and sedimentation coefficients, the wavenumber-dependent diffusion function determined in dynamic scattering experiments, and the high-frequency shear viscosity. The toolbox is based on the hydrodynamic radius model (HRM) wherein the internal particle structure is mapped on a hydrodynamic radius parameter for unchanged direct interactions, and on an existing simulation data base for solvent-permeable and spherical annulus particles. Useful scaling relations for the diffusion function and self-diffusion coefficient, known to be valid for hard-core interaction, are shown to apply also for soft pair potentials. We further discuss extensions of the toolbox to long-time transport properties including the low-shear zero-frequency viscosity and the long-time self-diffusion coefficient. The versatility of the toolbox is demonstrated by the analysis of a previous light scattering study of suspensions of non-ionic PNiPAM microgels [Eckert et al., J. Chem. Phys., 2008, 129, 124902] in which a detailed theoretical analysis of the dynamic data was left as an open task. By the comparison with Hertz potential based calculations, we show that the experimental data are consistently and accurately described using the Verlet-Weis corrected Percus-Yevick structure factor as input, and for a solvent penetration length equal to three percent of the excluded volume radius. This small amount of solvent permeability of the microgel particles has a significant dynamic effect at larger concentrations.

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

Using the Rogers-Young (RY) integral equation scheme for the static pair correlation functions combined with the liquid-phase Hansen-Verlet freezing rule, we study the generic behavior of the radial distribution function and static structure factor of monodisperse charge-stabilized suspensions with Yukawa-type repulsive particle interactions at freezing. In a related article, labeled Paper I [J. Gapinski, G. Nägele, and A. Patkowski, J. Chem. Phys. 136, 024507 (2012)], this hybrid method was used to determine two-parameter freezing lines for experimentally controllable parameters, characteristic of suspensions of charged silica spheres in dimethylformamide. A universal scaling of the RY radial distribution function maximum is shown to apply to the liquid-bcc and liquid-fcc segments of the universal freezing line. A thorough analysis is made of the behavior of characteristic distances and wavenumbers, next-neighbor particle coordination numbers, osmotic compressibility factor, and the Ravaché-Mountain-Streett minimum-maximum radial distribution function ratio.

A unifying mode-coupling theory for transport properties of electrolyte solutions. I. General scheme and limiting laws.

We develop a general method for calculating conduction-diffusion transport properties of strong electrolyte mixtures, including specific conductivities, steady-state electrophoretic mobilities, and self-diffusion coefficients. The ions are described as charged Brownian spheres, and the solvent-mediated hydrodynamic interactions (HIs) are also accounted for in the non-instantaneous ion atmosphere relaxation effect. A linear response expression relating long-time partial mobilities to associated dynamic structure factors is employed in our derivation of a general mode coupling theory (MCT) method for the conduction-diffusion properties. A simplified solution scheme for the MCT method is discussed. Analytic results are obtained for transport coefficients of pointlike ions which, for very low ion concentrations, reduce to the Deby-Falkenhagen-Onsager-Fuoss limiting law expressions. As an application, an unusual non-monotonic concentration dependence of the polyion electrophoretic mobility in a mixture of two binary electrolytes is discussed. In addition, leading-order extensions of the limiting law results are derived with HIs included. The present method complements a related MCT method by the authors for the electrolyte viscosity and shear relaxation function [C. Contreras-Aburto and G. Nägele, J. Phys.: Condens. Matter 24, 464108 (2012)], so that a unifying scheme for conduction-diffusion and viscoelastic properties is obtained. We present here the general framework of the method, illustrating its versatility for conditions where fully analytic results are obtainable. Numerical results for conduction-diffusion properties and the viscosity of concentrated electrolytes are presented in Paper II [C. Contreras Aburto and G. Nägele, J. Chem. Phys. 139, 134110 (2013)].

A unifying mode-coupling theory for transport properties of electrolyte solutions. II. Results for equal-sized ions electrolytes.

On the basis of a versatile mode-coupling theory (MCT) method developed in Paper I [C. Contreras Aburto and G. Nägele, J. Chem. Phys. 139, 134109 (2013)], we investigate the concentration dependence of conduction-diffusion linear transport properties for a symmetric binary electrolyte solution. The ions are treated in this method as charged Brownian spheres, and the solvent-mediated ion-ion hydrodynamic interactions are accounted for also in the ion atmosphere relaxation effect. By means of a simplified solution scheme, convenient semi-analytic MCT expressions are derived for the electrophoretic mobilities, and the molar conductivity, of an electrolyte mixture with equal-sized ions. These expressions reduce to the classical Debye-Falkenhagen-Onsager-Fuoss results in the limit of very low ion concentration. The MCT expressions are numerically evaluated for a binary electrolyte, and compared to experimental data and results by another theoretical method. Our analysis encloses, in addition, the electrolyte viscosity. To analyze the dynamic influence of the hydration shell, the significance of mixed slip-stick hydrodynamic surface boundary conditions, and the effect of solvent permeability are explored. For the stick boundary condition employed in the hydrodynamic diffusivity tensors, our theoretical results for the molar conductivity and viscosity of an aqueous 1:1 electrolyte are in good overall agreement with reported experimental data for aqueous NaCl solutions, for concentrations extending even up to two molar.

Sterically stabilized colloids with tunable repulsions.

When studying tunable electrostatic repulsions in aqueous suspensions of charged colloids, irreversible colloid aggregation or gelation may occur at high salt concentrations. For many commonly used synthetic colloids, such as polystyrene and silica particles, the reason for coagulation is the presence of unbalanced, strongly attractive, and short-ranged van der Waals (VDW) forces. Here, we present an aqueous polystyrene model colloid that is sterically stabilized against VDW attractions. We show that the synthesis procedure, based on a neutral initiator couple and a nonionic surfactant, introduces surface charges that can be further increased by the addition of charged comonomer methacrylic acid. Thus, the interactions between the polystyrene spheres can be conveniently tuned from hard-sphere-like to charge-stabilized with long-ranged electrostatic repulsions described by a Yukawa-type pair potential. The particle size, grafting density, core-shell structure, and surface charge are characterized by light and neutron scattering. Using X-ray and neutron scattering in combination with an accurate analytic integral equation scheme for the colloidal static structure factor, we deduce effective particle charges for colloid volume fractions ≥0.1 and salt concentrations in the range of 1.5 to 50 mM.

Viscosity of electrolyte solutions: a mode-coupling theory.

We present a versatile theoretical method for calculating the steady-state viscosity and shear relaxation function of strong electrolyte solutions. In this method, the ions are described on a primitive model level as charged Brownian spheres, and the essential ion-ion hydrodynamic interactions (HIs) are accounted for in the shear relaxation effect of the ionic atmosphere. The method combines a many-component mode-coupling theory (MCT) approach by Nägele et al (1998 J. Chem. Phys. 108 9893) with a simplified solution scheme, leading to an analytic expression for the shear relaxation contribution to the viscosity. This expression accounts for both the excluded volumes of the ions and their HIs. We show that the limiting law results for the viscosity of electrolyte mixtures by Falkenhagen and by Onsager and Fuoss are recovered at very low concentrations, and we discuss HIs corrections appearing at higher concentrations. Our numerical results for a 1:1 electrolyte reveal a strong enlargement of the viscosity caused by the HIs. The high-frequency viscosity gives the largest contribution to the total viscosity at higher concentrations.

Electro-kinetics of charged-sphere suspensions explored by integral low-angle super-heterodyne laser Doppler velocimetry.

We investigated the flow behaviour of colloidal charged-sphere suspensions using a newly designed integral low-angle super-heterodyne laser Doppler velocimetry instrument, which combines the advantages of several previous approaches. Sample conditions ranged from strong electrostatic interactions with pronounced short-range order to individual particles with no spatial correlations. The obtained power spectra correspond to diffusion broadened velocity distributions across the complete sample cross section. The excellent performance of the instrument is highlighted in detail by the example of electro-kinetic flow of suspensions in a closed cell of a rectangular cross section. We demonstrate the excellent performance of our approach with the example of electro-phoretic-electro-osmotic experiments, showing that a comprehensive flow characterization becomes possible by analysing the measured electro-kinetic mobilities, the flow-profile, an effective diffusion coefficient and the integrated scattering density. We briefly discuss present limitations, possible extensions and interesting applications in other fields.