Effective screening of hydrodynamic interactions in charged colloidal suspensions

Hydrodynamic interactions hinder transport of flow-driven colloidal particles

The flow-driven transport of interacting micron-sized particles occurs in many soft matter systems spanning from the translocation of proteins to moving emulsions in microfluidic devices. Here we combine experiments and theory to investigate the collective transport properties of colloidal particles along a rotating ring of optical traps. In the corotating reference frame, the particles are driven by a vortex flow of the surrounding fluid. When increasing the depth of the optical potential, we observe a jamming behavior that manifests itself in a strong reduction of the current with increasing particle density. We show that this jamming is caused by hydrodynamic interactions that enhance the energetic barriers between the optical traps. This leads to a transition from an over- to an under-critical tilting of the potential in the corotating frame. Based on analytical considerations, the enhancement effect is estimated to increase with increasing particle size or decreasing radius of the ring of traps. Measurements for different ring radii and Stokesian dynamics simulations for corresponding particle sizes confirm this. The enhancement of potential barriers in the flow-driven system is contrasted to the reduction of barriers in a force-driven one. This diverse behavior demonstrates that hydrodynamic interactions can have a very different impact on the collective dynamics of many-body systems. Applications to soft matter and biological systems require careful consideration of the driving mechanism and of the role of hydrodynamic interactions.

Reentrant melting of lanes of rough circular disks

We consider binary suspension of rough, circular particles in two dimensions under athermal conditions. The suspension is subject to a time-independent external drive in response to which half of the particles are pulled along the field direction, whereas the other half is pushed in the opposite direction. Simulating the system with different magnitude of external drive in steady state, we obtain oppositely moving macroscopic lanes only for a moderate range of external drive. Below as well as above the range we obtain states with no lane. Hence we find that the no-lane state reenters along the axis of the external drive in the nonequilibrium phase diagram corresponding to the laning transition, with varying roughness of individual particles and external drive. Interparticle friction (contact dissipation) due to the roughness of the individual particle is the main player behind the reentrance of the no-lane state at high external drives.

Dynamical effective field model for interacting ferrofluids: II. The proper relaxation time and effects of dynamic correlations

The recently proposed dynamical effective field model (DEFM) is quantitatively accurate for ferrofluid dynamics. It is derived in paper I within the framework of dynamical density functional theory (DDFT) along with a phenomenological description of nonadiabatic effects. However, it remains to clarify how the characteristic rotational relaxation time of a dressed particle, denoted byτ_r, is quantitatively related to that of a bare particle, denoted byτr0. By building macro-micro connections via two different routes, I reveal that under some gentle assumptionsτ_rcan be identified with the mean time characterizing long-time rotational self-diffusion. I further introduce two simple but useful integrated correlation factors, describing the effects of quasi-static (adiabatic) and dynamic (nonadiabatic) inter-particle correlations, respectively. In terms of both the dynamic magnetic susceptibility is expressed in an illuminating and elegant form. Remarkably, it shows that the macro-micro connection is established via two successive steps: a dynamical coarse-graining with nonadiabatic effects accounted for by the dynamic factor, followed by equilibrium ensemble averaging captured by the static factor. By analyzing data from Brownian dynamics simulations on monodisperse interacting ferrofluids, I findτr/τr0is, somehow unexpectedly, insensitive to changes of particle volume fraction. A physical picture is proposed to explain it. Furthermore, an empirical formula is proposed to characterize the dependence ofτr/τr0on dipole-dipole interaction strength. The DEFM supplemented with this formula leads to parameter-free predictions in good agreement with results from Brownian dynamics simulations. The theoretical developments presented in this paper may have important consequences to studies of ferrofluid dynamics in particular and other systems modeled by DDFTs in general.

Microsecond hydrodynamic interactions in dense colloidal dispersions probed at the European XFEL

Many soft-matter systems are composed of macromolecules or nanoparticles suspended in water. The characteristic times at intrinsic length scales of a few nanometres fall therefore in the microsecond and sub-microsecond time regimes. With the development of free-electron lasers (FELs) and fourth-generation synchrotron light-sources, time-resolved experiments in such time and length ranges will become routinely accessible in the near future. In the present work we report our findings on prototypical soft-matter systems, composed of charge-stabilized silica nanoparticles dispersed in water, with radii between 12 and 15 nm and volume fractions between 0.005 and 0.2. The sample dynamics were probed by means of X-ray photon correlation spectroscopy, employing the megahertz pulse repetition rate of the European XFEL and the Adaptive Gain Integrating Pixel Detector. We show that it is possible to correctly identify the dynamical properties that determine the diffusion constant, both for stationary samples and for systems driven by XFEL pulses. Remarkably, despite the high photon density the only observable induced effect is the heating of the scattering volume, meaning that all other X-ray induced effects do not influence the structure and the dynamics on the probed timescales. This work also illustrates the potential to control such induced heating and it can be predicted with thermodynamic models.

From Femtoseconds to Hours—Measuring Dynamics over 18 Orders of Magnitude with Coherent X-rays

X-ray photon correlation spectroscopy (XPCS) enables the study of sample dynamics between micrometer and atomic length scales. As a coherent scattering technique, it benefits from the increased brilliance of the next-generation synchrotron radiation and Free-Electron Laser (FEL) sources. In this article, we will introduce the XPCS concepts and review the latest developments of XPCS with special attention on the extension of accessible time scales to sub-μs and the application of XPCS at FELs. Furthermore, we will discuss future opportunities of XPCS and the related technique X-ray speckle visibility spectroscopy (XSVS) at new X-ray sources. Due to its particular signal-to-noise ratio, the time scales accessible by XPCS scale with the square of the coherent flux, allowing to dramatically extend its applications. This will soon enable studies over more than 18 orders of magnitude in time by XPCS and XSVS.

Structure and Dynamics of Ribonuclease A during Thermal Unfolding: The Failure of the Zimm Model

Disordered regions as found in intrinsically disordered proteins (IDP) or during protein folding define response time to stimuli and protein folding times. Neutron spin-echo spectroscopy is a powerful tool to directly access the collective motions of the unfolded chain to enlighten the physical origin of basic conformational relaxation. During the thermal unfolding of native ribonuclease A, we examine the structure and dynamics of the disordered state within a two-state transition model using polymer models, including internal friction, to describe the chain dynamics. The presence of four disulfide bonds alters the disordered configuration to a more compact configuration compared to a Gaussian chain that is defined by the additional links, as demonstrated by coarse-grained simulation. The dynamics of the disordered chain is described by Zimm dynamics with internal friction (ZIF) between neighboring amino acids. Relaxation times are dominated by mode-independent internal friction. Internal friction relaxation times show an Arrhenius-like behavior with an activation energy of 33 kJ/mol. The Zimm dynamics is dominated by internal friction and suggest that the characteristic motions correspond to overdamped elastic modes similar to the motions observed for folded proteins but within a pool of disordered configurations spanning the configurational space. For IDP, internal friction dominates while solvent friction and hydrodynamic interactions are smaller corrections.

Yield stress fluids and fundamental particle statistics

Yield stress in complex fluids is described by resorting to fundamental statistical mechanics for clusters with different particle occupancy numbers. Probability distribution functions are determined for canonical ensembles of volumes displaced at the incipient motion in three representative states (single, double, and multiple occupancies). The statistical average points out an effective solid fraction by which the yield stress behavior is satisfactorily described in a number of aqueous (Si₃N₄, Ca₃(PO₄)₂, ZrO₂, and TiO₂) and non-aqueous (Al₂O₃/decalin and MWCNT/PC) disperse systems. Interestingly, the only two model coefficients (maximum packing fraction and stiffness parameter) turn out to be correlated with the relevant suspension quantities. The latter relates linearly with (Young’s and bulk) mechanical moduli, whereas the former, once represented versus the Hamaker constant of two particles in a medium, returns a good linear extrapolation of the packing fraction for the simple cubic cell, here recovered within a relative error ≈ 1.3%.

Yield stress in complex fluids is described by resorting to fundamental statistical mechanics for clusters with different particle occupancy numbers.

Static structure and dynamical behavior of colloidal liquid crystals consisting of hydroxyapatite-based nanorod hybrids

Biominerals such as bones and teeth have elaborate nanostructures composed of aligned anisotropic hydroxyapatite (HAp) nanocrystals, which results in excellent mechanical properties. Construction of such ordered structures of HAp nanocrystals in synthetic materials is challenging. Recently, we reported that HAp-nanorod-based colloidal liquid crystals could be obtained. In the present study, the static structure and dynamics of liquid-crystalline (LC) colloidal dispersions of HAp nanorods are investigated by using small-angle X-ray scattering (SAXS) and X-ray photon correlation spectroscopy (XPCS). The SAXS results reveal that the interparticle distance decreases with increasing HAp concentration, φHAp, and the decrease of the interparticle distance for the short-axis direction is significantly smaller in the LC phase than the interparticle distance in the isotropic phase. In the dynamical studies of the LC phase using XPCS, we observe the diffusive motion of the HAp colloids, with the diffusion coefficient being dependent on the wave number. The diffusive motion slows down with increasing φHAp. We observe anisotropic dynamics after long-term storage (160 days after sealing), whereas only isotropic dynamics are observed in the initial XPCS measurements after short-term storage (14 days after sealing). Moreover, we have found that the dynamics slows down with increasing storage time.

Preparation and characterization of a biodegradable polyurethane hydrogel and the hybrid gel with soy protein for 3D cell-laden bioprinting

Thermo-responsive hydrogels of a polyurethane–soy protein hybrid provide unique rheological properties for 3D bioprinting and a biomimetic environment for neural repair.

Correlated two-particle diffusion in dense colloidal suspensions at early times: Theory and comparison to experiment

The spatially resolved diffusive dynamic cross correlations of a pair of colloids in dense quasi-two-dimensional monolayers of identical particles are studied experimentally and theoretically at early times where motion is Fickian. In very dense systems where strong oscillatory equilibrium packing correlations are present, we find an exponential decay of the dynamic cross correlations on small and intermediate length scales. At large separations where structure becomes random, an apparent power law decay with an exponent of approximately -2.2 is observed. For a moderately dense suspension where local structural correlations are essentially absent, this same apparent power law decay is observed over all probed interparticle separations. A microscopic nonhydrodynamic theory is constructed for the dynamic cross correlations which is based on interparticle frictional effects and effective structural forces. Hydrodynamics enters only via setting the very short-time single-particle self-diffusion constant. No-adjustable-parameter quantitative predictions of the theory for the dynamic cross correlations are in good agreement with experiment over all length scales. The origin of the long-range apparent power law is the influence of the constraint of fixed interparticle separation on the amplitude of the mean square force exerted on the two tagged particles by the surrounding fluid. The theory is extended to study high-packing-fraction 3D hard sphere fluids. The same pattern of an oscillatory exponential form of the dynamic cross correlation function is predicted in the structural regime, but the long-range tail decays faster than in monolayers with an exponent of -3.

Length-scale dependent transport properties of colloidal and protein solutions for prediction of crystal nucleation rates

We propose a scaling equation describing transport properties (diffusion and viscosity) in the solutions of colloidal particles. We apply the equation to 23 different systems including colloids and proteins differing in size (range of diameters: 4 nm to 1 μm), and volume fractions (10(-3)-0.56). In solutions under study colloids/proteins interact via steric, hydrodynamic, van der Waals and/or electrostatic interactions. We implement contribution of those interactions into the scaling law. Finally we use our scaling law together with the literature values of the barrier for nucleation to predict crystal nucleation rates of hard-sphere like colloids. The resulting crystal nucleation rates agree with existing experimental data.

Nucleation in short-range attractive colloids: ordering and symmetry of clusters

Results from extensive Brownian dynamics simulations are presented for nucleation in a system of colloidal particles interacting via a short-range attractive potential. Our analysis shows that, even though the system is not in equilibrium, structures of small size clusters compare well with the theoretically predicted and experimentally observed ground state structures for short-range colloidal systems. In addition, the distribution of the symmetric structures in nucleation is comparable to the distribution seen in equilibrium. We also investigate how the shape and structure of fluctuating clusters in the prenucleation regime affect the formation of a stable nucleating cluster.

High contrast x-ray speckle from atomic-scale order in liquids and glasses

The availability of ultrafast pulses of coherent hard x rays from the Linac Coherent Light Source opens new opportunities for studies of atomic-scale dynamics in amorphous materials. Here, we show that single ultrafast coherent x-ray pulses can be used to observe the speckle contrast in the high-angle diffraction from liquid Ga and glassy Ni(2)Pd(2)P and B(2)O(3). We determine the thresholds above which the x-ray pulses disturb the atomic arrangements. Furthermore, high contrast speckle is observed in scattering patterns from the glasses integrated over many pulses, demonstrating that the source and optics are sufficiently stable for x-ray photon correlation spectroscopy studies of dynamics over a wide range of time scales.

Hydrodynamic interactions in active colloidal crystal microrheology

In dense colloids it is commonly assumed that hydrodynamic interactions do not play a role. However, a found theoretical quantification is often missing. We present computer simulations that are motivated by experiments where a large colloidal particle is dragged through a colloidal crystal. To qualify the influence of long-ranged hydrodynamics, we model the setup by conventional Langevin dynamics simulations and by an improved scheme with limited hydrodynamic interactions. This scheme significantly improves our results and allows to show that hydrodynamics strongly impacts the development of defects, the crystal regeneration, as well as the jamming behavior.

Shear history independence in colloidal aggregation

Stimulated by experiments, we have carried out detailed simulations of aggregation in the presence of shear in a model colloidal system with a short-range attractive potential. For weak shear rates, we find that the shear enhanced the aggregation and that the long-time state of the system is independent of the shear history. For strong shear rates, precipitous fragmentation occurred after the shear was turned on and, after an induction period, the aggregation quickly rebounded in a stochastic manner similar to classical nucleation phenomena. However, the long-time state of the system is, once again, independent of the shear history. Thus, for both weak and strong shear cases, the shear rate acts as a state variable of the aggregating system. Shear rates employed in the simulations can be attained in laboratory experiments, as confirmed by computing the dimensionless Péclet numbers.

Spinodal fractionation in a polydisperse square-well fluid

Using kinetic Monte Carlo simulation, we model gas-liquid spinodal decomposition in a size-polydisperse square well fluid, representing a "near-monodisperse" colloidal dispersion. We find that fractionation (demixing) of particle sizes between the phases begins asserting itself shortly after the onset of phase ordering. Strikingly, the direction of size fractionation can be reversed by a seemingly trivial choice between two interparticle potentials which, in the monodisperse case, are identical--we rationalize this in terms of a perturbative, equilibrium theory of polydispersity. Furthermore, our quantitative results show that kinetic Monte Carlo simulation can provide detailed insight into the role of fractionation in real colloidal systems.

Short-time rheology and diffusion in suspensions of Yukawa-type colloidal particles

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. 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 η(∞), 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.

Dynamics in dense hard-sphere colloidal suspensions

The dynamic behavior of a hard-sphere colloidal suspension was studied by x-ray photon correlation spectroscopy and small-angle x-ray scattering over a wide range of particle volume fractions. The short-time mobility of the particles was found to be smaller than that of free particles even at relatively low concentrations, showing the importance of indirect hydrodynamic interactions. Hydrodynamic functions were derived from the data, and for moderate particle volume fractions (Φ≤ 0.40) there is good agreement with earlier many-body theory calculations by Beenakker and Mazur [Physica A 120, 349 (1984)]. Important discrepancies appear at higher concentrations, above Φ≈ 0.40, where the hydrodynamic effects are overestimated by the Beenakker-Mazur theory, but predicted accurately by an accelerated Stokesian dynamics algorithm developed by Banchio and Brady [J. Chem. Phys. 118, 10323 (2003)]. For the relaxation rates, good agreement was also found between the experimental data and a scaling form predicted by the mode coupling theory. In the high concentration range, with the fluid suspensions approaching the glass transition, the long-time diffusion coefficient was compared with the short-time collective diffusion coefficient to verify a scaling relation previously proposed by Segrè and Pusey [Phys. Rev. Lett. 77, 771 (1996)]. We discuss our results in view of previous experimental attempts to validate this scaling law [L. Lurio et al., Phys. Rev. Lett. 84, 785 (2000)].

Short-time diffusion of charge-stabilized colloidal particles: generic features

Analytical theory and Stokesian dynamics simulations are used in conjunction with dynamic light scattering to investigate the role of hydrodynamic interactions in short-time diffusion in suspensions of charge-stabilized colloidal particles. The particles are modeled as solvent-impermeable charged spheres, repelling each otherviaa screened Coulomb potential. Numerical results for self-diffusion and sedimentation coefficients, as well as hydrodynamic and short-time diffusion functions, are compared with experimental data for a wide range of volume fractions. The theoretical predictions for the generic behavior of short-time properties obtained from this model are shown to be in full accord with experimental data. In addition, the effects of microion kinetics, nonzero particle porosity and residual attractive forces on the form of the hydrodynamic function are estimated. This serves to rule out possible causes for the strikingly small hydrodynamic function values determined in certain synchrotron radiation experiments.

Exploring soft matter with x-rays: from the discovery of the DNA structure to the challenges of free electron lasers

X-rays have long been a precious tool for the study of the structure of matter. While the short wavelength makes them ideal for investigating materials down to the atomic scale, their high penetration power allows for the exploration of opaque samples at a multitude of length scales. We give an overview of the x-ray techniques suited for the characterization of soft matter and of their application to systems of current interest. We describe the advantages and limitations of existing x-ray methods and outline the possible developments following the introduction of a new kind of coherent source: the x-ray free electron laser.

Generic behavior of the hydrodynamic function of charged colloidal suspensions

We discuss the generic behavior of the hydrodynamic function H(q) and diffusion function D(q) characterizing the short-time diffusion in suspensions of charge-stabilized colloidal spheres, by covering the whole fluid regime. Special focus is given to the behavior of these functions at the freezing transition specified by the Hansen-Verlet freezing rule. Results are presented in dependence on scattering wavenumber q, effective particle charge, volume fraction, salt concentration, and particle size, by considering both the low-charge and high-charge branch solutions of static structure factors. The existence of two charge branches leads to the prediction of a re-entrant melting-freezing-melting transition for increasing particle concentration at very low salinity. A universal limiting contour line is derived for the principal peak height value of H(q), independent of particle charge and diameter, and concentration and salinity, which separates the fluid from the fluid-solid coexistence region. This line is only weakly dependent on the value of the structure factor peak height entering the Hansen-Verlet rule. A dynamic freezing criterion is derived in terms of the short-time cage diffusion coefficient, a quantity easily measurable in a scattering experiment. The higher-dimensional parameter scans underlying this study make use of the fast and highly efficient deltagamma-scheme in conjunction with the analytic rescaled mean spherical approximation input for the static structure factor. Our results constitute a comprehensive database useful to researchers performing dynamic scattering experiments on charge-stabilized dispersions.

Self-assembly of ligated gold nanoparticles: phenomenological modeling and computer simulations

We study the assembly of ligated gold nanoparticles by both phenomenological modeling and computer simulations for various ligand chain lengths. First, we develop an effective nanoparticle-nanoparticle pair potential by treating the ligands as flexible polymer chains. Besides van der Waals interactions, we incorporate both the free energy of mixing and elastic contributions from compression of the ligands in our effective pair potentials. The separation of the nanoparticles at the potential minimum compares well with experimental results of gold nanoparticle superlattice constants for various ligand lengths. Next, we use the calculated pair potentials as input to Brownian dynamics simulations for studying the formation of nanoparticle assembly in three dimensions. For dodecanethiol ligated nanoparticles in toluene, our model gives a relatively shallower well depth and the clusters formed after a temperature quench are compact in morphology. Simulation results for the kinetics of cluster growth in this case are compared with phase separations in binary mixtures. For decanethiol ligated nanoparticles, the model well depth is found to be deeper, and simulations show hybrid, fractal-like morphology for the clusters. Cluster morphology in this case shows a compact structure at short length scales and a fractal structure at large length scales. Growth kinetics for this deeper potential depth is compared with the diffusion-limited cluster-cluster aggregation (DLCA) model.

Dynamics in dense suspensions of charge-stabilized colloidal particles

The dynamic behavior of charge-stabilized colloidal particles in suspension was studied by photon correlation spectroscopy with coherent X-rays (XPCS). The short-time diffusion coefficient, D(Q) , was measured for volume concentrations phi < or = 0.18 and compared to the free particle diffusion constant D(0) and the static structure factor S(Q) . The data show that indirect, hydrodynamic interactions are relevant for the system and hydrodynamic functions were derived. The results are in striking contrast to the predictions of the PA (pairwise-additive approximation) model, but show features typical for a hard-sphere system. The observed mobility is however considerably smaller than the one of a respective hard-sphere system. The hydrodynamic functions can be modelled quantitatively if one allows for an increased effective viscosity relative to the hard-sphere case.

Tunable colloids: control of colloidal phase transitions with tunable interactions

Systems of spherical colloidal particles mimic the thermodynamics of atomic crystals. Control of interparticle interactions in colloids, which has recently begun to be extensively exploited, gives rise to rich phase behaviours as well as crystal structures with nanoscale and micron-scale lattice spacings. This provides model systems in which to study fundamental problems in condensed matter physics, such as the dynamics of crystal nucleation and melting, and the nature of the glass transition, at experimentally accessible lengthscales and timescales. Tunable control of these interactions provides reversible control. This will enable quantitative studies of phase transition kinetics as well as the creation of advanced materials with switchability of function and properties.

Collective diffusion in charge-stabilized suspensions: Concentration and salt effects

The authors present a joint experimental-theoretical study of collective diffusion properties in aqueous suspensions of charge-stabilized fluorinated latex spheres. Small-angle x-ray scattering and x-ray photon correlation spectroscopy have been used to explore the concentration and ionic-strength dependence of the static and short-time dynamic properties including the hydrodynamic function H(q), the wave-number-dependent collective diffusion coefficient D(q), and the intermediate scattering function over the entire accessible range. They show that all experimental data can be quantitatively described and explained by means of a recently developed accelerated Stokesian dynamics simulation method, in combination with a modified hydrodynamic many-body theory. In particular, the behavior of H(q) for de-ionized and dense suspensions can be attributed to the influence of many-body hydrodynamics, without any need for postulating hydrodynamic screening to be present, as it was done in earlier work. Upper and lower boundaries are provided for the peak height of the hydrodynamic function and for the short-time self-diffusion coefficient over the entire range of added salt concentrations.

Diffraction with a coherent X-ray beam: dynamics and imaging

Methods for carrying out coherent X-ray scattering experiments are reviewed. The brilliance of the available synchrotron sources, the characteristics of the existing optics, the various ways of obtaining a beam of controlled coherence properties and the detectors used are summarized. Applications in the study of the dynamics of speckle patterns are described. In the case of soft condensed matter, the movement of inclusions like fillers in polymers or colloidal particles can be observed and these can reflect polymer or liquid-crystal fluctuations. In hard condensed-matter problems, like phase transitions, charge-density waves or phasons in quasicrystals, the study of speckle fluctuations provides new time-resolved methods. In the domain of lensless imaging, the coherent beam gives the modulus of the sample Fourier transform. If oversampling conditions are fulfilled, the phase can be obtained and the image in the direct space can be reconstructed. The forthcoming improvements of all these techniques are discussed.

Confinement-induced liquid ordering investigated by x-ray phase retrieval

Using synchrotron x-ray diffraction, we have determined the ensemble-averaged density profile of colloidal fluids within confining channels of different widths. We observe an oscillatory ordering-disordering behavior of the colloidal particles as a function of the channel width, while the colloidal solution remains in the liquid state. This phenomenon has been suggested by surface force studies of hard-sphere fluids and also theoretically predicted, but here we see it by direct measurements of the structure for comparable systems.

Many-body hydrodynamic interactions in charge-stabilized suspensions

In this joint experimental-theoretical work we study hydrodynamic interaction effects in dense suspensions of charged colloidal spheres. Using x-ray photon correlation spectroscopy we have determined the hydrodynamic function H(q), for a varying range of electrosteric repulsion. We show that H(q) can be quantitatively described by means of a novel Stokesian dynamics simulation method for charged Brownian spheres, and by a modification of a many-body theory developed originally by Beenakker and Mazur. Very importantly, we can explain the behavior of H(q) for strongly correlated particles without resorting to the controversial concept of hydrodynamic screening, as was attempted in earlier work by Riese [Phys. Rev. Lett. 85, 5460 (2000)].

Crossover from intermittent to continuum dynamics for locally driven colloids

We simulate a colloid with charge q(d) driven through a disordered assembly of interacting colloids with charge q and show that, for q(d) approximately q, the velocity-force relation is nonlinear and the velocity fluctuations of the driven particle are highly intermittent with a 1/f characteristic. When g(d) >>q , the average velocity drops, the velocity-force relation becomes linear, and the velocity fluctuations are Gaussian. We discuss the results in terms of a crossover from strongly intermittent heterogeneous dynamics to continuum dynamics. We also make several predictions for the transient response in the different regimes.

Diffusion and microstructural properties of solutions of charged nanosized proteins: Experiment versus theory

We have reanalyzed our former static small-angle x-ray scattering and photon correlation spectroscopy results on dense solutions of charged spherical apoferritin proteins using theories recently developed for studies of colloids. The static structure factors S(q), and the small-wave-number collective diffusion coefficient D(c) determined from those experiments are interpreted now in terms of a theoretical scheme based on a Derjaguin-Landau-Verwey-Overbeek-type continuum model of charged colloidal spheres. This scheme accounts, in an approximate way, for many-body hydrodynamic interactions. Stokesian dynamics computer simulations of the hydrodynamic function have been performed for the first time for dense charge-stabilized dispersions to assess the accuracy of the theoretical scheme. We show that the continuum model allows for a consistent description of all experimental results, and that the effective particle charge is dependent upon the protein concentration relative to the added salt concentration. In addition, we discuss the consequences of small ions dynamics for the collective protein diffusion within the framework of the coupled-mode theory.

Spontaneously formed monodisperse biomimetic unilamellar vesicles: the effect of charge, dilution, and time

Using small-angle neutron scattering and dynamic light scattering, we have constructed partial structural phase diagrams of lipid mixtures composed of the phosphatidylcholines dimyristoyl and dihexanoyl doped with calcium ions (Ca2+) and/or the negatively charged lipid, dimyristoyl phosphatidylglycerol (DMPG). For dilute solutions (lipid concentration < or =1 wt %), spontaneously forming unilamellar vesicles (ULVs) were found, and their polydispersity was determined to be approximately 20%. The stability of the Ca2+- or DMPG-doped ULVs was monitored over a period of 4 days and their structural parameters (e.g., average outer radius, <Ro>) were found to be insensitive to the lipid concentration (Clp). However, doping the dimyristoyl/dihexanoyl system with both Ca2+ and DMPG resulted in ULVs whose <Ro> was found to be Clp dependent. The <Ro> of DMPG-doped ULVs remained unchanged over an extended period of time (at least 4 days), a good indication of their stability.

Structure factor scaling in colloidal phase separation

The dynamical scaling hypothesis for the structure factor, S (q) , in depletion-driven colloidal phase separation is studied by carrying out Brownian dynamics simulations. A true dynamical scaling is observed for shallow quenches into the two-phase coexistence region. In such a quench, compact clusters nucleate and grow with time and there is only one characteristic length scale in the system after an initial transient period. Scaling is satisfied beyond this initial period. In contrast, deep quenches lead to fractal cluster growth, and the system is controlled by two characteristic lengths that evolve differently in time [Huang, Oh, and Sorensen (HOS), Phys. Rev. E 57, 875 (1998)]. True dynamical scaling thus cannot be expected to hold. However, an apparent scaling for the structure factor is observed over some period of time when these two characteristic length scales become comparable to each other. We compare our simulation results for the total structure factor to theoretical predictions by HOS by writing it as a product of cluster-cluster and the averaged single-cluster structure factors, each with its own characteristic length.

Kinetics of phase transformations in depletion-driven colloids

We present results from a detailed numerical study of the kinetics of phase transformations in a model two-dimensional depletion-driven colloidal system. Transition from a single, dispersed phase to a two-phase coexistence of monomers and clusters is obtained as the depth of the interaction potential among the colloidal particles is changed. Increasing the well depth further, fractal clusters are observed in the simulation. These fractal clusters have a hybrid structure in the sense that they show hexagonal closed-packed crystalline ordering at short length scales and a ramified fractal nature at larger length scales. For sufficiently deep potential wells, the diffusion-limited cluster-cluster aggregation model is recovered in terms of the large-scale fractal dimension Df of the clusters, the kinetic exponent z, and the scaling form of the cluster size distribution. For shallower well depths inside the two-phase coexistence region, simulation results for the kinetics of cluster growth are compared with intermediate-stage phase separation in binary mixtures. In the single-phase region, growth kinetics agree well with a mean-field aggregation-fragmentation model of Sorensen, Zhang, and Taylor.

Molecular dynamics simulation of the transition from dispersed to solid phase

Molecular dynamics simulations in two dimensions of particles interacting with finite potentials comparable to k(B)T yield aggregates which cross over from a fractal to a compact crystalline morphology. Growth kinetics and aggregate size distributions evolve from nonequilibrium to equilibrium limits.

Average hydrodynamic correction for the Brownian dynamics calculation of flocculation rates in concentrated dispersions

In order to account for the hydrodynamic interaction (HI) between suspended particles in an average way, Honig et al. [J. Colloid Interface Sci. 36, 97 (1971)] and more recently Heyes [Mol. Phys. 87, 287 (1996)] proposed different analytical forms for the diffusion constant. While the formalism of Honig et al. strictly applies to a binary collision, the one from Heyes accounts for the dependence of the diffusion constant on the local concentration of particles. However, the analytical expression of the latter approach is more complex and depends on the particular characteristics of each system. Here we report a combined methodology, which incorporates the formula of Honig et al. at very short distances and a simple local volume-fraction correction at longer separations. As will be shown, the flocculation behavior calculated from Brownian dynamics simulations employing the present technique, is found to be similar to that of Batchelor's tensor [J. Fluid. Mech. 74, 1 (1976); 119, 379 (1982)]. However, it corrects the anomalous coalescence found in concentrated systems as a result of the overestimation of many-body HI.