-
Effective interactions, structure, and pressure in charge-stabilized colloidal suspensions: Critical assessment of charge renormalization methods
Authors:
Mariano E. Brito,
Gerhard Nägele,
Alan R. Denton
Abstract:
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 c…
▽ More
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.
△ Less
Submitted 15 December, 2023;
originally announced December 2023.
-
Clustering in quasi-two-dimensional dispersions of Brownian particles with competitive interactions: Phase diagram and structural properties
Authors:
Zihan Tan,
Vania Calandrini,
Jan K. G. Dhont,
Gerhard Nägele
Abstract:
Competing short-range attractive (SA) and long range repulsive (LR) interactions have been invoked to describe colloid or protein solutions, as well as membrane proteins interactions mediated by lipid molecules. Using Langevin dynamics simulations, we determine the generalized phase diagram, the cluster shapes and size distributions of a generic Q2D dispersion of spherical SALR particles confined…
▽ More
Competing short-range attractive (SA) and long range repulsive (LR) interactions have been invoked to describe colloid or protein solutions, as well as membrane proteins interactions mediated by lipid molecules. Using Langevin dynamics simulations, we determine the generalized phase diagram, the cluster shapes and size distributions of a generic Q2D dispersion of spherical SALR particles confined to in-plane motion. SA and LR interactions are modelled by a generalized Lenard-Jones potential and a screened Coulomb potential, respectively. The microstructures of the various equilibrium and non-equilibrium phases turn out to be distinctly different from the ones observed in three-dimensional (3D) SALR systems. We discuss perturbation theory predictions for the metastable binodal line of a reference system of particles with SA interactions only, which in the Q2D-SALR phase diagram separates cluster from non-cluster phases. The transition from the high-temperature (low SA) dispersed fluid phase to the lower-temperature equilibrium cluster phase is characterised by a low-wavenumber peak of the static structure factor (corresponding to a thermal correlation length of about twice the particle diameter) featuring a distinctly smaller height ($\approx1.4$) than in 3D SALR systems. By further decreasing the temperature (increasing SA), the cluster morphology changes from disk-like shapes in the equilibrium cluster phase, to double-stranded anisotropic hexagonal cluster forms in the cluster-percolated gel phase. This transition is quantified by the hexagonal order parameter distribution function. The mean cluster size and coordination number of particles in the gel phase are insensitive to changes in the attraction strength.
△ Less
Submitted 8 October, 2023;
originally announced October 2023.
-
Hydrodynamics of Immiscible Binary Fluids with Viscosity Contrast: A multiparticle collision dynamics approach
Authors:
Zihan Tan,
Vania Calandrini,
Jan K. G. Dhont,
Gerhard Nägele,
Roland G. Winkler
Abstract:
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 trans…
▽ More
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.
△ Less
Submitted 5 September, 2021; v1 submitted 4 May, 2021;
originally announced May 2021.
-
Modeling deswelling, thermodynamics, structure, and dynamics in ionic microgel suspensions
Authors:
Mariano E. Brito,
Alan R. Denton,
Gerhard Nägele
Abstract:
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 comprehensiv…
▽ More
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 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- and collective diffusion and the osmotic pressure, lowers the suspension viscosity, and significantly shifts the suspension crystallization point to higher concentrations. The paper presents a bottom-up approach to efficiently computing suspension properties of crowded ionic microgels using single-particle characteristics.
△ Less
Submitted 2 October, 2019;
originally announced October 2019.
-
Short- and long-time diffusion, and dynamic scaling in suspensions of charged colloidal particles
Authors:
Adolfo J. Banchio,
Marco Heinen,
Peter Holmqvist,
Gerhard Nägele
Abstract:
We report on a comprehensive theory-simulation-experimental study of collective and self-diffusion in suspensions of charge-stabilized colloidal spheres. In simulation and theory, the spheres interact by a hard-core plus screened Coulomb pair potential. Intermediate and self-intermediate scattering functions are calculated by accelerated Stokesian Dynamics simulations where hydrodynamic interactio…
▽ More
We report on a comprehensive theory-simulation-experimental study of collective and self-diffusion in suspensions of charge-stabilized colloidal spheres. In simulation and theory, the spheres interact by a hard-core plus screened Coulomb pair potential. Intermediate and self-intermediate scattering functions are calculated by accelerated Stokesian Dynamics simulations where hydrodynamic interactions (HIs) are fully accounted for. The study spans the range from the short-time to the colloidal long-time regime. Additionally, Brownian Dynamics simulation and mode-coupling theory (MCT) results are generated where HIs are neglected. It is shown that HIs enhance collective and self-diffusion at intermediate and long times, whereas at short times self-diffusion, and for certain wavenumbers also collective diffusion, are slowed down. MCT significantly overestimate the slowing influence of dynamic particle caging. The simulated scattering functions are in decent agreement with our dynamic light scattering (DLS) results for suspensions of charged silica spheres. Simulation and theoretical results are indicative of a long-time exponential decay of the intermediate scattering function. The approximate validity of a far-reaching time-wavenumber factorization of the scattering function is shown to be a consequence of HIs. Our study of collective diffusion is amended by simulation and theoretical results for the self-intermediate scattering function and the particle mean squared displacement (MSD). Since self-diffusion is not assessed in DLS measurements, a method to deduce the MSD approximately in DLS is theoretically validated.
△ Less
Submitted 27 November, 2017;
originally announced November 2017.
-
Ultrafiltration of charge-stabilized dispersions at low salinity
Authors:
Rafael Roa,
Daniel Menne,
Jonas Riest,
Pompilia Buzatu,
Emiliy K. Zholkovskiy,
Jan K. G. Dhont,
Matthias Wessling,
Gerhard Nägele
Abstract:
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 comp…
▽ More
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 $D_c$. A strong electro-hydrodynamic enhancement of $D_c$ 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 $D_c$ and $η$. According to our filtration model, $D_c$ 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.
△ Less
Submitted 18 May, 2016; v1 submitted 18 March, 2016;
originally announced March 2016.
-
Ultrafiltration modeling of non-ionic microgels
Authors:
Rafael Roa,
Emiliy K. Zholkovskiy,
Gerhard Nägele
Abstract:
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 contribution…
▽ More
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.
△ Less
Submitted 17 May, 2015; v1 submitted 23 March, 2015;
originally announced March 2015.
-
Dynamics of suspensions of hydrodynamically structured particles: Analytic theory and experiment
Authors:
Jonas Riest,
Thomas Eckert,
Walter Richtering,
Gerhard Nägele
Abstract:
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…
▽ More
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 solvent permeability of the microgel particles has a significant dynamic effect at larger concentrations.
△ Less
Submitted 12 January, 2015;
originally announced January 2015.
-
Rotational self-diffusion in suspensions of charged particles: Revised Beenakker-Mazur and Pairwise Additivity methods versus numerical simulations
Authors:
Karol Makuch,
Marco Heinen,
Gustavo C. Abade,
Gerhard Nägele
Abstract:
To the present day, the Beenakker-Mazur (BM) method is the most comprehensive statistical physics approach to the calculation of short-time transport properties of colloidal suspensions. A revised version of the BM method with an improved treatment of hydrodynamic interactions is presented and evaluated regarding the rotational short-time self-diffusion coefficient, $D^r$ , of suspensions of charg…
▽ More
To the present day, the Beenakker-Mazur (BM) method is the most comprehensive statistical physics approach to the calculation of short-time transport properties of colloidal suspensions. A revised version of the BM method with an improved treatment of hydrodynamic interactions is presented and evaluated regarding the rotational short-time self-diffusion coefficient, $D^r$ , of suspensions of charged particles interacting by a hard-sphere plus screened Coulomb (Yukawa) pair potential. To assess the accuracy of the method, elaborate simulations of $D^r$ have been performed, covering a broad range of interaction parameters and particle concentrations. The revised BM method is compared in addition with results by a simplifying pairwise additivity (PA) method in which the hydrodynamic interactions are treated on a two-body level. The static pair correlation functions re- quired as input to both theoretical methods are calculated using the Rogers-Young integral equation scheme. While the revised BM method reproduces the general trends of the simulation results, it systematically and significantly underestimates the rotational diffusion coefficient. The PA method agrees well with the simulation data at lower volume fractions, but at higher concentrations $D^r$ is likewise underestimated. For a fixed value of the pair potential at mean particle distance comparable to the thermal energy, $D^r$ increases strongly with increasing Yukawa potential screening parameter.
△ Less
Submitted 7 January, 2015;
originally announced January 2015.
-
Electrokinetic and hydrodynamic properties of charged-particles systems: From small electrolyte ions to large colloids
Authors:
G. Nägele,
M. Heinen,
A. J. Banchio,
C. Contreras-Aburto
Abstract:
Dynamic processes in dispersions of charged spherical particles are of importance both in fundamental science, and in technical and bio-medical applications. There exists a large variety of charged-particles systems, ranging from nanometer-sized electrolyte ions to micron-sized charge-stabilized colloids. We review recent advances in theoretical methods for the calculation of linear transport coef…
▽ More
Dynamic processes in dispersions of charged spherical particles are of importance both in fundamental science, and in technical and bio-medical applications. There exists a large variety of charged-particles systems, ranging from nanometer-sized electrolyte ions to micron-sized charge-stabilized colloids. We review recent advances in theoretical methods for the calculation of linear transport coefficients in concentrated particulate systems, with the focus on hydrodynamic interactions and electrokinetic effects. Considered transport properties are the dispersion viscosity, self- and collective diffusion coefficients, sedimentation coefficients, and electrophoretic mobilities and conductivities of ionic particle species in an external electric field. Advances by our group are also discussed, including a novel mode-coupling-theory method for conduction-diffusion and viscoelastic properties of strong electrolyte solutions. Furthermore, results are presented for dispersions of solvent-permeable particles, and particles with non-zero hydrodynamic surface slip. The concentration-dependent swelling of ionic microgels is discussed, as well as a far-reaching dynamic scaling behavior relating colloidal long- to short-time dynamics.
△ Less
Submitted 13 September, 2013;
originally announced September 2013.
-
Effect of glycerol and dimethyl sulfoxide on the phase behavior of lysozyme: Theory and experiments
Authors:
Christoph Gögelein,
Dana Wagner,
Frédéric Cardinaux,
Gerhard Nägele,
Stefan U. Egelhaaf
Abstract:
Salt, glycerol and dimethyl sulfoxide (DMSO) are used to modify the properties of protein solutions. We experimentally determined the effect of these additives on the phase behavior of lysozyme solutions. Upon the addition of glycerol and DMSO, the fluid-solid transition and the gas-liquid coexistence curve (binodal) shift to lower temperatures and the gap between them increases. The experimentall…
▽ More
Salt, glycerol and dimethyl sulfoxide (DMSO) are used to modify the properties of protein solutions. We experimentally determined the effect of these additives on the phase behavior of lysozyme solutions. Upon the addition of glycerol and DMSO, the fluid-solid transition and the gas-liquid coexistence curve (binodal) shift to lower temperatures and the gap between them increases. The experimentally observed trends are consistent with our theoretical predictions based on the thermodynamic perturbation theory (TPT) and the Derjaguin-Landau-Verwey-Overbeek (DLVO) model for the lysozyme-lysozyme pair interactions. The values of the parameters describing the interactions, namely the refractive indices, dielectric constants, Hamaker constant and cut-off length, are extracted from literature or are experimentally determined by independent experiments, including static light scattering to determine the second virial coefficient. We observe that both, glycerol and DMSO, render the potential more repulsive, while sodium chloride reduces the repulsion.
△ Less
Submitted 23 December, 2011; v1 submitted 12 December, 2011;
originally announced December 2011.
-
Dynamics of charged gibbsite platelets in the isotropic phase
Authors:
Dzina Kleshchanok,
Marco Heinen,
Gerhard Nägele,
Peter Holmqvist
Abstract:
We report on depolarized and non-depolarized dynamic light scattering, static light scattering, and static viscosity measurements on interacting charged gibbsite platelets suspended in dimethyl sulfoxide (DMSO). The average collective and (long-time) translational self-diffusion coefficients, and the rotational diffusion coefficient, have been measured as functions of the platelet volume fraction…
▽ More
We report on depolarized and non-depolarized dynamic light scattering, static light scattering, and static viscosity measurements on interacting charged gibbsite platelets suspended in dimethyl sulfoxide (DMSO). The average collective and (long-time) translational self-diffusion coefficients, and the rotational diffusion coefficient, have been measured as functions of the platelet volume fraction φ, up to the isotropic-liquid crystal (I/LC) transition. The non-depolarized intensity autocorrelation function, measured at low scattering wavenumbers, consists of a fast and a slowly decaying mode which we interpret as the orientationally averaged collective and translational self-diffusion coefficients, respectively. Both the rotational and the long-time self-diffusion coefficients decrease very strongly, by more than two orders of magnitude, in going from the very dilute limit to the I/LC transition concentration. A similarly strong decrease, with increasing φ, is observed for the inverse zero-strain limiting static shear viscosity. With increasing φ, increasingly strong shear-thinning is observed, accompanied by a shrinking of the low shear-rate Newtonian plateau. The measured diffusion coefficients are interpreted theoretically in terms of a simple model of effective charged spheres interacting by a screened Coulomb potential, with hydrodynamic interactions included. The disk-like particle shape, and the measured particle radius and thickness polydispersities, enter into the model calculations via the scattering amplitudes. The interaction-induced enhancement of the collective diffusion coefficient by more than a factor of 20 at larger volume fractions is well captured in the effective sphere model, whereas the strong declines both of the experimental translational and rotational self-diffusion coefficients are underestimated.
△ Less
Submitted 15 September, 2011;
originally announced September 2011.
-
Viscosity and Diffusion: Crowding and Salt Effects in Protein Solutions
Authors:
Marco Heinen,
Fabio Zanini,
Felix Roosen-Runge,
Diana Fedunová,
Fajun Zhang,
Marcus Hennig,
Tilo Seydel,
Ralf Schweins,
Michael Sztucki,
Marián Antalík,
Frank Schreiber,
Gerhard Nägele
Abstract:
We report on a joint experimental-theoretical study of collective diffusion in, and static shear viscosity of solutions of bovine serum albumin (BSA) proteins, focusing on the dependence on protein and salt concentration. Data obtained from dynamic light scattering and rheometric measurements are compared to theoretical calculations based on an analytically treatable spheroid model of BSA with iso…
▽ More
We report on a joint experimental-theoretical study of collective diffusion in, and static shear viscosity of solutions of bovine serum albumin (BSA) proteins, focusing on the dependence on protein and salt concentration. Data obtained from dynamic light scattering and rheometric measurements are compared to theoretical calculations based on an analytically treatable spheroid model of BSA with isotropic screened Coulomb plus hard-sphere interactions. The only input to the dynamics calculations is the static structure factor obtained from a consistent theoretical fit to a concentration series of small-angle X-ray scattering (SAXS) data. This fit is based on an integral equation scheme that combines high accuracy with low computational cost. All experimentally probed dynamic and static properties are reproduced theoretically with an at least semi-quantitative accuracy. For lower protein concentration and low salinity, both theory and experiment show a maximum in the reduced viscosity, caused by the electrostatic repulsion of proteins. The validity range of a generalized Stokes-Einstein (GSE) relation connecting viscosity, collective diffusion coefficient, and osmotic compressibility, proposed by Kholodenko and Douglas [PRE 51, 1081 (1995)] is examined. Significant violation of the GSE relation is found, both in experimental data and in theoretical models, in semi-dilute systems at physiological salinity, and under low-salt conditions for arbitrary protein concentrations.
△ Less
Submitted 14 September, 2011;
originally announced September 2011.
-
Short-time Rheology and Diffusion in Suspensions of Yukawa-type Colloidal Particles
Authors:
Marco Heinen,
Adolfo J. Banchio,
Gerhard Nägele
Abstract:
A comprehensive study is presented on the short-time dynamics in suspensions of charged colloidal spheres. The explored parameter space covers the major part of the fluid-state regime, with colloid concentrations extending up to the freezing transition. The particles are assumed to interact directly by a hard-core plus screened Coulomb potential, and indirectly by solvent-mediated hydrodynamic int…
▽ More
A comprehensive study is presented on the short-time dynamics in suspensions of charged colloidal spheres. The explored parameter space covers the major part of the fluid-state regime, with colloid concentrations extending up to the freezing transition. The particles are assumed to interact directly by a hard-core plus screened Coulomb potential, and indirectly by solvent-mediated hydrodynamic interactions (HIs). By comparison with accurate accelerated Stokesian Dynamics (ASD) simulations of the hydrodynamic function H(q), and the high-frequency viscosity, we investigate the accuracy of two fast and easy-to-implement analytical schemes. The first scheme, referred to as the pairwise additive (PA) scheme, uses exact two-body hydrodynamic mobility tensors. It is in good agreement with the ASD simulations of H(q) and the high-frequency viscosity, for smaller volume fractions up to about 10% and 20%, respectively. The second scheme is a hybrid method combining the virtues of the δγ-scheme by Beenakker and Mazur with those of the PA scheme. It leads to predictions in good agreement with the simulation data, for all considered concentrations, combining thus precision with computational efficiency. The hybrid method is used to test the accuracy of a generalized Stokes-Einstein (GSE) relation proposed by Kholodenko and Douglas, showing its severe violation in low salinity systems. For hard spheres, however, this GSE relation applies decently well.
△ Less
Submitted 6 September, 2011;
originally announced September 2011.
-
Rotational and translational self-diffusion in concentrated suspensions of permeable particles
Authors:
Gustavo C. Abade,
Bogdan Cichocki,
Maria L. Ekiel-Jezewska,
Gerhard Naegele,
Eligiusz Wajnryb
Abstract:
In our recent work on concentrated suspensions of uniformly porous colloidal spheres with excluded volume interactions, a variety of short-time dynamic properties were calculated, except for the rotational self-diffusion coefficient. This missing quantity is included in the present paper. Using a precise hydrodynamic force multipole simulation method, the rotational self-diffusion coefficient is e…
▽ More
In our recent work on concentrated suspensions of uniformly porous colloidal spheres with excluded volume interactions, a variety of short-time dynamic properties were calculated, except for the rotational self-diffusion coefficient. This missing quantity is included in the present paper. Using a precise hydrodynamic force multipole simulation method, the rotational self-diffusion coefficient is evaluated for concentrated suspensions of permeable particles. Results are presented for particle volume fractions up to 45%, and for a wide range of permeability values. From the simulation results and earlier results for the first-order virial coefficient, we find that the rotational self-diffusion coefficient of permeable spheres can be scaled to the corresponding coefficient of impermeable particles of the same size. We also show that a similar scaling applies to the translational self-diffusion coefficient considered earlier. From the scaling relations, accurate analytic approximations for the rotational and translational self-diffusion coefficients in concentrated systems are obtained, useful to the experimental analysis of permeable-particle diffusion. The simulation results for rotational diffusion of permeable particles are used to show that a generalized Stokes-Einstein-Debye relation between rotational self-diffusion coefficient and high-frequency viscosity is not satisfied.
△ Less
Submitted 19 March, 2011;
originally announced March 2011.
-
First-order virial expansion of short-time diffusion and sedimentation coefficients of permeable particles suspensions
Authors:
Bogdan Cichocki,
Maria L. Ekiel-Jezewska,
G. Naegele,
E. Wajnryb
Abstract:
For suspensions of permeable particles, the short-time translational and rotational self-diffusion coefficients, and collective diffusion and sedimentation coefficients are evaluated theoretically. An individual particle is modeled as a uniformly permeable sphere of a given permeability, with the internal solvent flow described by the Debye-Bueche-Brinkman equation. The particles are assumed to in…
▽ More
For suspensions of permeable particles, the short-time translational and rotational self-diffusion coefficients, and collective diffusion and sedimentation coefficients are evaluated theoretically. An individual particle is modeled as a uniformly permeable sphere of a given permeability, with the internal solvent flow described by the Debye-Bueche-Brinkman equation. The particles are assumed to interact non-hydrodynamically by their excluded volumes. The virial expansion of the transport properties in powers of the volume fraction is performed up to the two-particle level. The first-order virial coefficients corresponding to two-body hydrodynamic interactions are evaluated with very high accuracy by the series expansion in inverse powers of the inter-particle distance. Results are obtained and discussed for a wide range of the ratio, x, of the particle radius to the hydrodynamic screening length inside a permeable sphere. It is shown that for x >= 10, the virial coefficients of the transport properties are well-approximated by the hydrodynamic radius (annulus) model developed by us earlier for the effective viscosity of porous-particle suspensions.
△ Less
Submitted 23 January, 2011;
originally announced January 2011.
-
A simple patchy colloid model for the phase behavior of lysozyme dispersions
Authors:
Christoph Gogelein,
Gerhard Nagele,
Remco Tuinier,
Thomas Gibaud,
Anna Stradner,
Peter Schurtenberger
Abstract:
We propose a minimal model for spherical proteins with aeolotopic pair interactions to describe the equilibrium phase behavior of lysozyme. The repulsive screened Coulomb interactions between the particles are taken into account assuming that the net charges are smeared out homogeneously over the spherical protein surfaces. We incorporate attractive surface patches, with the interactions between…
▽ More
We propose a minimal model for spherical proteins with aeolotopic pair interactions to describe the equilibrium phase behavior of lysozyme. The repulsive screened Coulomb interactions between the particles are taken into account assuming that the net charges are smeared out homogeneously over the spherical protein surfaces. We incorporate attractive surface patches, with the interactions between patches on different spheres modeled by an attractive Yukawa potential. The parameters entering the attractive Yukawa potential part are determined using information on the experimentally accessed gas-liquid-like critical point. The Helmholtz free energy of the fluid and solid phases is calculated using second-order thermodynamic perturbation theory. Our predictions for the solubility curve are in fair agreement with experimental data. In addition, we present new experimental data for the gas-liquid coexistence curves at various salt concentrations and compare these with our model calculations. In agreement with earlier findings, we observe that the strength and the range of the attractive potential part only weakly depend on the salt content.
△ Less
Submitted 30 January, 2009;
originally announced January 2009.
-
Non-monotonic density dependence of the diffusion of DNA fragments in low-salt suspensions
Authors:
M. G. McPhie,
G. Naegele
Abstract:
The high linear charge density of 20-base-pair oligomers of DNA is shown to lead to a striking non-monotonic dependence of the long-time self-diffusion on the concentration of the DNA in low-salt conditions. This generic non-monotonic behavior results from both the strong coupling between the electrostatic and solvent-mediated hydrodynamic interactions, and from the renormalization of these elec…
▽ More
The high linear charge density of 20-base-pair oligomers of DNA is shown to lead to a striking non-monotonic dependence of the long-time self-diffusion on the concentration of the DNA in low-salt conditions. This generic non-monotonic behavior results from both the strong coupling between the electrostatic and solvent-mediated hydrodynamic interactions, and from the renormalization of these electrostatic interactions at large separations, and specifically from the dominance of the far-field hydrodynamic interactions caused by the strong repulsion between the DNA fragments.
△ Less
Submitted 26 November, 2008;
originally announced November 2008.
-
Super-heterodyne light scattering on interacting colloidal suspensions: theory and experiment
Authors:
Thomas Palberg,
Holger Reiber,
Tetyana Köller,
Martin Medebach,
Gerhard Nägele
Abstract:
In soft matter structure couples to flow and vice versa. Complementary to structural investigations, we here are interested in the determination of particle velocities of charged colloidal suspensions of different structure under flow. In a combined effort of theory and experiment we determine the Fourier transform of the super-heterodyne field auto-correlation function (power spectrum) which in…
▽ More
In soft matter structure couples to flow and vice versa. Complementary to structural investigations, we here are interested in the determination of particle velocities of charged colloidal suspensions of different structure under flow. In a combined effort of theory and experiment we determine the Fourier transform of the super-heterodyne field auto-correlation function (power spectrum) which in frequency space is found to be well separated from homodyne contributions and low frequency noise. Under certain conditions the power spectrum is dominated by incoherently scattered light, originating from the unavoidable size polydispersity of colloidal particles. A simple approximate form for the low-wavenumber self-intermediate scattering function is proposed, reminiscent to the case of non-interacting particles. We experimentally scrutinize the range of applicability of these simplified calculations on employing a parabolic electro-osmotic flow profile. Both for non-interacting and strongly interacting fluid particle systems, the spectra are well described as diffusion-broadened velocity distributions comprising an osmotic flow-averaged superposition of Lorentzians at distinct locations. We discuss the performance and scope of this approach with particular focus on moderately strong interactions and on multiphase flow. In addition, we point to some remaining theoretical challenges in connection to the observed linear increase of the effective diffusion constant and the integrated spectral power with increasing electric field strength.
△ Less
Submitted 14 November, 2008;
originally announced November 2008.
-
Structure and Short-time Dynamics in Suspensions of Charged Silica Spheres in the entire Fluid Regime
Authors:
J. Gapinski,
A. Patkowski,
A. J. Banchio,
J. Buitenhuis,
P. Holmqvist,
M. P. Lettinga,
G. Meier,
G. Naegele
Abstract:
We present an experimental study of short-time diffusion properties in fluid-like suspensions of monodisperse charge-stabilized silica spheres suspended in DMF. The static structure factor S(q), the short-time diffusion function, D(q), and the hydrodynamic function, H(q), in these systems have been probed by combining X-ray photon correlation spectroscopy experiments with static small-angle X-ra…
▽ More
We present an experimental study of short-time diffusion properties in fluid-like suspensions of monodisperse charge-stabilized silica spheres suspended in DMF. The static structure factor S(q), the short-time diffusion function, D(q), and the hydrodynamic function, H(q), in these systems have been probed by combining X-ray photon correlation spectroscopy experiments with static small-angle X-ray scattering. Our experiments cover the full liquid-state part of the phase diagram, including deionized systems right at the liquid-solid phase boundary. We show that the dynamic data can be consistently described by the renormalized density fluctuation expansion theory of Beenakker and Mazur over a wide range of concentrations and ionic strengths. In accord with this theory and Stokesian dynamics computer simulations, the measured short-time properties cross over monotonically, with increasing salt content, from the bounding values of salt-free suspensions to those of neutral hard spheres. Moreover, we discuss an upper bound for the hydrodynamic function peak height of fluid systems based on the Hansen-Verlet freezing criterion.
△ Less
Submitted 7 November, 2008;
originally announced November 2008.
-
Sedimentation of strongly and weakly charged colloidal particles: Prediction of fractional density dependence
Authors:
M. Watzlawek,
G. Naegele
Abstract:
We report on calculations of the reduced sedimentation velocity $U/U_{0}$ in homogenous suspensions of strongly and weakly charged colloidal spheres as a function of particle volume fraction $φ$. For dilute suspensions of strongly charged spheres at low salinity, $U/U_{0}$ is well represented by the parametric form $1-pφ^α$ with a fractional exponent $α=1/3$ and a parameter $p\simeq 1.8$, which…
▽ More
We report on calculations of the reduced sedimentation velocity $U/U_{0}$ in homogenous suspensions of strongly and weakly charged colloidal spheres as a function of particle volume fraction $φ$. For dilute suspensions of strongly charged spheres at low salinity, $U/U_{0}$ is well represented by the parametric form $1-pφ^α$ with a fractional exponent $α=1/3$ and a parameter $p\simeq 1.8$, which is essentially independent from the macroion charge $Z$. This non-linear volume fraction dependence can be quantitatively understood in terms of a model of effective hard spheres with $φ$-dependent diameter. For weakly charged spheres in a deionized solvent, we show that the exponent $α$ can be equal to 1/2, if an expression for $U/U_0$ given by Petsev and Denkov [J. Colloid Interface Sci. 149, 329 (1992)] is employed. We further show that the range of validity of this expression is limited to very small values of $φ$ and $Z$, which are probably not accessible in sedimentation experiments. The presented results might also hold for other systems like spherical proteins or ionic micelles.
△ Less
Submitted 9 July, 1999; v1 submitted 5 March, 1999;
originally announced March 1999.
-
Diffusion of colloids at short times
Authors:
M. Watzlawek,
G. Naegele
Abstract:
We study the combined effects of electrostatic and hydrodynamic interactions (HI) on the short-time dynamics of charge-stabilized colloidal spheres. For this purpose, we calculate the translational and the rotational self-diffusion coefficients, $D^t_s$ and $D^r_s$, as function of volume fraction $φ$ for various values of the effective particle charge $Z$ and various concentrations $n_s$ of adde…
▽ More
We study the combined effects of electrostatic and hydrodynamic interactions (HI) on the short-time dynamics of charge-stabilized colloidal spheres. For this purpose, we calculate the translational and the rotational self-diffusion coefficients, $D^t_s$ and $D^r_s$, as function of volume fraction $φ$ for various values of the effective particle charge $Z$ and various concentrations $n_s$ of added 1--1 electrolyte. Our results show that the self-diffusion coefficients in deionized suspensions are less affected by HI than in suspensions with added electrolyte. For very large $n_s$, we recover the well-known results for hard spheres, i.e. a linear $φ$-dependence of $D^t_s$ and $D^r_s$ at small $φ$. In contrast, for deionized charged suspensions at small $φ$, we observe the interesting non-linear scaling properties $D^t_s\propto1-a_tφ^{4/3}$ and $D^r_s\propto 1-a_rφ^2$. The coefficients $a_t$ and $a_r$ are found to be nearly independent of $Z$. The qualitative differences between the dynamics of charged and uncharged particles can be well explained in terms of an effective hard sphere (EHS) model.
△ Less
Submitted 28 July, 1997;
originally announced July 1997.
-
Self-diffusion coefficients of charged particles: Prediction of Nonlinear volume fraction dependence
Authors:
M. Watzlawek,
G. Naegele
Abstract:
We report on calculations of the translational and rotational short-time self-diffusion coefficients $D^t_s$ and $D^r_s$ for suspensions of charge-stabilized colloidal spheres. These diffusion coefficients are affected by electrostatic forces and many-body hydrodynamic interactions (HI). Our computations account for both two-body and three-body HI. For strongly charged particles, we predict inte…
▽ More
We report on calculations of the translational and rotational short-time self-diffusion coefficients $D^t_s$ and $D^r_s$ for suspensions of charge-stabilized colloidal spheres. These diffusion coefficients are affected by electrostatic forces and many-body hydrodynamic interactions (HI). Our computations account for both two-body and three-body HI. For strongly charged particles, we predict interesting nonlinear scaling relations $D^t_s\propto 1-a_tφ^{4/3}$ and $D^r_s\propto 1-a_rφ^2$ depending on volume fraction $φ$, with essentially charge-independent parameters $a_t$ and $a_r$. These scaling relations are strikingly different from the corresponding results for hard spheres. Our numerical results can be explained using a model of effective hard spheres. Moreover, we perceptibly improve the known result for $D^t_s$ of hard sphere suspensions.
△ Less
Submitted 28 July, 1997;
originally announced July 1997.
-
Short-time rotational diffusion in monodisperse charge-stabilized colloidal suspensions
Authors:
M. Watzlawek,
G. Naegele
Abstract:
We investigate the combined effects of electrostatic interactions and hydrodynamic interactions on the short-time rotational self-diffusion coefficient in charge-stabilized suspensions. We calculate this coefficient as a function of volume fraction for various effective particle charges and amounts of added electrolyte. The influence of the hydrodynamic interactions on the rotational diffusion c…
▽ More
We investigate the combined effects of electrostatic interactions and hydrodynamic interactions on the short-time rotational self-diffusion coefficient in charge-stabilized suspensions. We calculate this coefficient as a function of volume fraction for various effective particle charges and amounts of added electrolyte. The influence of the hydrodynamic interactions on the rotational diffusion coefficient is less pronounced for charged particles than for uncharged ones. Salt-free suspensions are weakly influenced by hydrodynamic interactions. For these strongly correlated systems we obtain a quadratic volume fraction-dependence of the diffusion coefficient, which is well explained in terms of an effective hard sphere model.
△ Less
Submitted 16 April, 1996;
originally announced April 1996.