Effective static and high-frequency viscosities of concentrated suspensions of soft particles

We obtain an analytic expression that allows to determine the static η and high-frequency η(∞) viscosities as function of the volume fraction φ of a concentrated suspension of soft spherical particles in a liquid of viscosity η(0). The particles consist of a hard core of radius a covered by a porous layer of thickness d. Suspensions of hard spheres and homogeneous porous particles are limiting cases of the model. The proposed expression incorporates the results for the intrinsic viscosity obtained on the basis of a cell model [H. Ohshima, Langmuir 26, 6287 (2010)] into a recently obtained relation for the effective viscosity of concentrated colloidal suspensions [C. I. Mendoza and I. Santamari?a-Holek, J. Chem. Phys. 130, 044904 (2009); J. Colloid. Interface Sci. 346, 118 (2010)]. In this model, the correlations between the particles due to crowding effects are introduced through an effective volume fraction φ(eff) which is then used as integration variable in a differential effective medium procedure. The final expression is simple, accurate, and allows to collapse all the data in a universal master curve that is independent of the parameters characterizing the system. The only difference between the static and high-frequency cases is that in the later case φ(eff) also incorporates hydrodynamic interactions arising from the so-called relaxation term. We have tested the accuracy of our model comparing with experimental results for spherical polymeric brushes and simulations for the high-frequency viscosity of homogeneous porous particles. In all cases the agreement with the data is extremely good.     

Short-time transport properties in dense suspensions: from neutral to charge-stabilized colloidal spheres

We present a detailed study of short-time dynamic properties in concentrated suspensions of charge-stabilized and of neutral colloidal spheres. The particles in many of these systems are subject to significant many-body hydrodynamic interactions. A recently developed accelerated Stokesian dynamics (ASD) simulation method is used to calculate hydrodynamic functions, wave-number-dependent collective diffusion coefficients, self-diffusion and sedimentation coefficients, and high-frequency limiting viscosities. The dynamic properties are discussed in dependence on the particle concentration and salt content. Our ASD simulation results are compared with existing theoretical predictions, notably those of the renormalized density fluctuation expansion method of Beenakker and Mazur [Physica A 126, 349 (1984)], and earlier simulation data on hard spheres. The range of applicability and the accuracy of various theoretical expressions for short-time properties are explored through comparison with the simulation data. We analyze, in particular, the validity of generalized Stokes-Einstein relations relating short-time diffusion properties to the high-frequency limiting viscosity, and we point to the distinctly different behavior of de-ionized charge-stabilized systems in comparison to hard spheres.     

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.     

Pair structure of the hard-sphere Yukawa fluid: an improved analytic method versus simulations, Rogers-Young scheme, and experiment

We present a comprehensive study of the equilibrium pair structure in fluids of nonoverlapping spheres interacting by a repulsive Yukawa-like pair potential, with special focus on suspensions of charged colloidal particles. The accuracy of several integral equation schemes for the static structure factor, S(q), and radial distribution function, g(r), is investigated in comparison to computer simulation results and static light scattering data on charge-stabilized silica spheres. In particular, we show that an improved version of the so-called penetrating-background corrected rescaled mean spherical approximation (PB-RMSA) by Snook and Hayter [Langmuir 8, 2880 (1992)], referred to as the modified PB-RMSA (MPB-RMSA), gives pair structure functions which are in general in very good agreement with Monte Carlo simulations and results from the accurate but nonanalytical and therefore computationally more expensive Rogers-Young integral equation scheme. The MPB-RMSA preserves the analytic simplicity of the standard rescaled mean spherical (RMSA) solution. The combination of high accuracy and fast evaluation makes the MPB-RMSA ideally suited for extensive parameter scans and experimental data evaluation, and for providing the static input to dynamic theories. We discuss the results of extensive parameter scans probing the concentration scaling of the pair structure of strongly correlated Yukawa particles, and we determine the liquid-solid coexistence line using the Hansen-Verlet freezing rule.     

Thermodynamically consistent closure approximation for hard spheres systems

We present a new closure relation that is an extension of the Percus-Yevick approximation. In the proposed closure, we introduce an additional term and a mixing coefficient that can be determined by imposing a condition of thermodynamic self-consistency. Moreover, the mixing coefficient is calculated analytically within a linear approximation. In the case of a monodisperse system of hard spheres, we compare the results of our model to well-established thermodynamic expressions and also to the structural properties of fairly known closure approximations. In the second case, and using an equivalent scheme, the new closure relation is extended to the depletion potential between two large hard spheres immersed in a liquid of small hard spheres. In both cases, the results of our model are in good agreement with numerical simulations performed at intermediate concentrations.     

Diffusion of spheres in crowded suspensions of rods

Translational tracer diffusion of spherical macromolecules in crowded suspensions of rodlike colloids is investigated. Experiments are done using several kinds of spherical tracers in fd-virus suspensions. A wide range of size ratios L/2a of the length L of the rods and the diameter 2a of the tracer sphere is covered by combining several experimental methods: fluorescence correlation spectroscopy for small tracer spheres, dynamic light scattering for intermediate sized spheres, and video microscopy for large spheres. Fluorescence correlation spectroscopy is shown to measure long-time diffusion only for relatively small tracer spheres. Scaling of diffusion coefficients with a/xi, predicted for static networks, is not found for our dynamical network of rods (with xi the mesh size of the network). Self-diffusion of tracer spheres in the dynamical network of freely suspended rods is thus fundamentally different as compared to cross-linked networks. A theory is developed for the rod-concentration dependence of the translational diffusion coefficient at low rod concentrations for freely suspended rods. The proposed theory is based on a variational solution of the appropriate Smoluchowski equation without hydrodynamic interactions. The theory can, in principle, be further developed to describe diffusion through dynamical networks at higher rod concentrations with the inclusion of hydrodynamic interactions. Quantitative agreement with the experiments is found for large tracer spheres, and qualitative agreement for smaller spheres. This is probably due to the increasing importance of hydrodynamic interactions as compared to direct interactions as the size of the tracer sphere decreases.     

Diffusion, sedimentation, and rheology of concentrated suspensions of core-shell particles

Short-time dynamic properties of concentrated suspensions of colloidal core-shell particles are studied using a precise force multipole method which accounts for many-particle hydrodynamic interactions. A core-shell particle is composed of a rigid, spherical dry core of radius a surrounded by a uniformly permeable shell of outer radius b and hydrodynamic penetration depth κ(-1). The solvent flow inside the permeable shell is described by the Brinkman-Debye-Bueche equation, and outside the particles by the Stokes equation. The particles are assumed to interact non-hydrodynamically by a hard-sphere no-overlap potential of radius b. Numerical results are presented for the high-frequency shear viscosity, η(∞), sedimentation coefficient, K, and the short-time translational and rotational self-diffusion coefficients, D(t) and D(r). The simulation results cover the full three-parametric fluid-phase space of the composite particle model, with the volume fraction extending up to 0.45, and the whole range of values for κb, and a/b. Many-particle hydrodynamic interaction effects on the transport properties are explored, and the hydrodynamic influence of the core in concentrated systems is discussed. Our simulation results show that for thin or hardly permeable shells, the core-shell systems can be approximated neither by no-shell nor by no-core models. However, one of our findings is that for κ(b - a) ? 5, the core is practically not sensed any more by the weakly penetrating fluid. This result is explained using an asymptotic analysis of the scattering coefficients entering into the multipole method of solving the Stokes equations. We show that in most cases, the influence of the core grows only weakly with increasing concentration.     

Langevin dynamics for rigid bodies of arbitrary shape

We present an algorithm for carrying out Langevin dynamics simulations on complex rigid bodies by incorporating the hydrodynamic resistance tensors for arbitrary shapes into an advanced rotational integration scheme. The integrator gives quantitative agreement with both analytic and approximate hydrodynamic theories for a number of model rigid bodies and works well at reproducing the solute dynamical properties (diffusion constants and orientational relaxation times) obtained from explicitly solvated simulations.     

Monte Carlo simulation and self-consistent integral equation theory for polymers in quenched random media

The conformational properties and static structure of freely jointed hard-sphere chains in matrices composed of stationary hard spheres are studied using Monte Carlo simulations and integral equation theory. The simulations show that the chain size is a nonmonotonic function of the matrix density when the matrix spheres are the same size as the monomers. When the matrix spheres are of the order of the chain size the chain size decreases monotonically with increasing matrix volume fraction. The simulations are used to test the replica-symmetric polymer reference interaction site model (RSP) integral equation theory. When the simulation results for the intramolecular correlation functions are input into the theory, the agreement between theoretical predictions and simulation results for the pair-correlation functions is quantitative only at the highest fluid volume fractions and for small matrix sphere sizes. The RSP theory is also implemented in a self-consistent fashion, i.e., the intramolecular and intermolecular correlation functions are calculated self-consistently by combining a field theory with the integral equations. The theory captures qualitative trends observed in the simulations, such as the nonmonotonic dependence of the chain size on media fraction.     

Mesoscopic model for diffusion-influenced reaction dynamics

A hybrid mesoscopic multiparticle collision model is used to study diffusion-influenced reaction kinetics. The mesoscopic particle dynamics conserves mass, momentum, and energy so that hydrodynamic effects are fully taken into account. Reactive and nonreactive interactions with catalytic solute particles are described by full molecular dynamics. Results are presented for large-scale, three-dimensional simulations to study the influence of diffusion on the rate constants of the A + C <==> B + C reaction. In the limit of a dilute solution of catalytic C particles, the simulation results are compared with diffusion equation approaches for both the irreversible and reversible reaction cases. Simulation results for systems where the volume fraction phi of catalytic spheres is high are also presented, and collective interactions among reactions on catalytic spheres that introduce volume fraction dependence in the rate constants are studied.     

Lattice-Boltzmann simulation of the sedimentation of charged disks

We report a series of lattice-Boltzmann simulations of the sedimentation velocity of charged disks. In these simulations, we explicitly account for the hydrodynamic and electrostatic forces on disks and on their electrical double layer. By comparing our results with those for spheres with equal surface and charge, we can clarify the effect of the particle shape on the sedimentation process. We find that disks and spheres exhibit a different dependence of the sedimentation velocity on the Debye screening length. An analysis of the behavior of highly charged disks (beyond the scope of the linearized Poisson-Boltzmann equation) shows that, in that regime, the charge dependence of the sedimentation velocity of disks and spheres is similar. This suggests that, at high charge, the effective hydrodynamic shape of the disks becomes more spherical.     

Gaussian split Ewald: A fast Ewald mesh method for molecular simulation

Gaussian split Ewald (GSE) is a versatile Ewald mesh method that is fast and accurate when used with both real-space and k-space Poisson solvers. While real-space methods are known to be asymptotically superior to k-space methods in terms of both computational cost and parallelization efficiency, k-space methods such as smooth particle-mesh Ewald (SPME) have thus far remained dominant because they have been more efficient than existing real-space methods for simulations of typical systems in the size range of current practical interest. Real-space GSE, however, is approximately a factor of 2 faster than previously described real-space Ewald methods for the level of force accuracy typically required in biomolecular simulations, and is competitive with leading k-space methods even for systems of moderate size. Alternatively, GSE may be combined with a k-space Poisson solver, providing a conveniently tunable k-space method that performs comparably to SPME. The GSE method follows naturally from a uniform framework that we introduce to concisely describe the differences between existing Ewald mesh methods.     

Drag and Torque on Clusters of N Arbitrary Spheres at Low Reynolds Number

Hydrodynamics of particle clusters suspended in viscous fluids is a subject of considerable theoretical and practical importance. Using a multipole expansion of the flow velocity in a series of spherical harmonics, Lamb's fundamental solution of the Stokes flow outside a single sphere is generalized in this work to the case of N nonoverlapping spheres of arbitrary size with slip boundary conditions. The expansion coefficients are found by transforming the boundary conditions to the Lamb form and by transforming the spherical coordinates and solid spherical harmonics centered at different spheres. The problem is reduced to the solution of the linear system of equations for the expansion coefficients, which is carried out numerically. Based on the developed theory, the relation between the hydrodynamic and gyration radius of fractal-like aggregates with different structure is established. In another application, an asymptotic slip-regime dependence of the aggregate hydrodynamic radius on the Knudsen number and the number of particles is found by performing calculations of drag forces acting on the gas-borne fractal-like and straight chain aggregates. A good agreement is shown in comparing predictions of the described theory with available experimental and theoretical results on motion of various small sphere clusters in viscous fluid. Copyright 2000 Academic Press.     

Electro-osmosis at inhomogeneous charged surfaces: hydrodynamic versus electric friction

Electrokinetic methods are efficient in probing the electrostatic surface properties of charged systems. However, anomalies observed in experiments indicate that the classical electrokinetic theory should be reconsidered. Using Green's function methods and hydrodynamic simulations, we investigate electro-osmosis driven by electric-field-induced ion motion near a charged planar substrate with smooth or rough boundary. First, a reformulation of electro-osmotic theory for planar charged surfaces employing Green's functions shows that the Helmholtz-Smoluchowski (HS) relation between electrostatic potential and solvent velocity is exact for smooth surfaces, even in the presence of ion correlations. Deviations from HS theory are caused by combined hydrodynamic and electric surface friction, as our hydrodynamic simulations of ions at smooth and corrugated charged surfaces in lateral electric fields demonstrate. Within the simulations, hydrodynamic interactions are treated in the continuum limit and the presence of a no-slip boundary condition at the surface is taken into account. While electrofriction is relevant in highly charged system and/or for multivalent ions, hydrodynamic friction is dominant in systems with moderate surface charge density and/or low ionic valency. We also derive the effective electrokinetic surface charge from the electro-osmotic solvent profiles, which is substantially reduced when compared with the bare value and shows qualitative agreement with the experimental tendency.     

Dynamic Mobility of Two Spherical Particles with Thick Double Layers

The dynamic electrophoretic mobility of a pair of nearby spherical particles is analyzed in the case when the thickness of the electrical double layer around each particle is comparable to the particle radius. By means of an integral reciprocal relation, a formal expression is obtained for the force and torque on N spheres subject to an oscillating electric field which may be spatially varying. Upon linearizing in the surface potential, this expression is shown to depend upon a set of purely hydrodynamic problems involving N neutral spheres, the calculation of the electric field around N neutral spheres, and the equilibrium charge distribution around N charged spheres. In the case of a single particle, the known analytic formula for the dynamic mobility is recovered. For a pair of identical particles, the dynamic mobility is calculated numerically, using known solutions to the required subproblems. An analytical expression for the mobility of a pair of widely separated spheres is also obtained by a method of reflections, and this is in excellent agreement with the numerical results outside the range of double layer overlap. Copyright 2000 Academic Press.     

Freezing lines of colloidal Yukawa spheres. I. A Rogers-Young integral equation study

Using the Rogers-Young (RY) integral equation scheme for the static structure factor combined with the one-phase Hansen-Verlet (HV) freezing rule, we study the equilibrium structure and two-parameter freezing lines of colloidal particles with Yukawa-type pair interactions representing charge-stabilized silica spheres suspended in dimethylformamide (DMF). Results are presented for a vast range of concentrations, salinities and effective charges covering particles with masked excluded-volume interactions. The freezing lines were obtained for the low-charge and high-charge solutions of the static structure factor, for various two-parameter sets of experimentally accessible system parameters. All RY-HV based freezing lines can be mapped on a universal fluid-solid coexistence line in good agreement with computer simulation predictions. The RY-HV calculations extend the freezing lines obtained in earlier simulations to a broader parameter range. The experimentally observed fluid-bcc-fluid reentrant transition of charged silica spheres in DMF can be explained using the freezing lines obtained in this work.     

Prediction of a structural transition in the hard disk fluid

Starting from the second equilibrium equation in the BBGKY hierarchy under the Kirkwood superposition closure, we implement a new method for studying the asymptotic decay of correlations in the hard disk fluid in the high density regime. From our analysis and complementary numerical studies, we find that exponentially damped oscillations can occur only up to a packing fraction η(?)～0.718, a value that is in substantial agreement with the packing fraction, η～0.723, believed to characterize the transition from the ordered solid phase to a dense fluid phase, as inferred from Mak's Monte Carlo simulations [Phys. Rev. E 73, 065104 (2006)]. Next, we show that the same method of analysis predicts that the exponential damping of oscillations in the hard sphere fluid becomes impossible when λ=4nπσ(3)[1+H(1)]≥34.81, where H(1) is the contact value of the correlation function, n is the number density, and σ is the sphere diameter in exact agreement with the condition, λ≥34.8, which is first reported in a numerical study of the Kirkwood equation by Kirkwood et al. [J. Chem. Phys. 18, 1040 (1950)]. Finally, we show that our method confirms the absence of any structural transition in hard rods for the entire range of densities below close packing.     

Rate equation of gelation of chromium(III)-polyacrylamide sol

During gelation, the time dependence of the apparent viscosity has been used to analyze the kinetics of the gelation of polyacrylamide (PA) sol with chromium (III) ions. The investigations have stablished the following facts: Before gelation the relation between the viscosity, η, of Cr(III)-PA sol and the PA concentration is η = f([PA]) = 2.36×10¹⁵ [PA]^3.15 and the rate equation of gelation in the steady-state stage during gelation is expressed as V_n = k_s f([PA]) [PA]² [Cr(III)]².     

Effective medium theory expressions for the effective diffusion in chromatographic beds filled with porous, non-porous and porous-shell particles and cylinders. Part II: Numerical verification and quantitative effect of solid core on expected B-term band broadening

The results of a numerical simulation study of the diffusion and retention in fully porous spheres and cylinders are compared with some of the high order accuracy analytical solutions for the effective diffusion coefficient that have been derived from the effective medium theory (EMT) theory in part I of the present study. A variety of different ordered (spheres and cylinders) and disordered (cylinders) packings arrangements has been considered. The agreement between simulations and theory was always excellent, lying within the (very tight) accuracy limits of the simulations over the full range of retention factor and diffusion constant values that is practically relevant for most LC applications. Subsequently filling up the spheres and cylinders with a central solid core, while keeping the same packing geometry and the same mobile phase (same thermodynamic retention equilibrium), it was found that the core induces an additional obstruction which reduces the effective intra-particle diffusion coefficient exactly with a factor γ(part)=2/(2+ρ3) for spherical particles and γ(part)=1/(1+ρ2) for cylinders (ρ is the ratio of the core to the particle diameter, ρ=d(core)/d(part)). These expressions hold independently of the packing geometry, the value of the diffusion coefficients and the equilibrium constant or the size of the core. The expressions also imply that, if considering equal mobile phase conditions, the presence of the solid core will never reduce the particle contribution to the B-term band broadening with more than 33% (50% in case of cylindrical pillars).