Diffusion of spheres in isotropic and nematic suspensions of rods

Diffusion of a small tracer sphere (apoferritin) in isotropic and nematic networks [of fd virus] is discussed. For a tracer sphere that is smaller than the mesh size of the network, screened hydrodynamic interactions between the sphere and the network determine its diffusion coefficient. A theory is developed for such interactions as well as their relation to the long-time self-diffusion coefficient. Fluorescence correlation spectroscopy measurements on mixtures of apoferritin and fd virus are presented. The long-time self-diffusion coefficient of apoferritin is measured as a function of the fd-virus concentration, both in the isotropic and nematic state, in directions parallel and perpendicular to the nematic director. The hydrodynamic screening length of the fd-virus network as a function of fd concentration is obtained by combining these experimental data with the theory. Surprisingly, the screening length increases with increasing concentration in nematic networks. This is due to the increase in the degree of alignment, which apparently leads to a strong increase of the screening length. Hydrodynamic screening is thus strongly diminished by alignment. A self-consistent calculation of the screening length does not work at higher concentrations, probably due to the strong variation of the typical incident flow fields over the contour of a rod.     

Diffusion of spheres in isotropic and nematic networks of rods: electrostatic interactions and hydrodynamic screening

Translational diffusion of a small charged tracer sphere in isotropic and nematic suspensions of long and thin charged rods is investigated as a function of ionic strength and rod concentration. A theory for the diffusive properties of a small sphere is developed, where both (screened) hydrodynamic interactions and charge interactions between the tracer sphere and the rod network are analyzed. Hydrodynamic interactions are formulated in terms of the hydrodynamic screening length. As yet, there are no independent theoretical predictions for the hydrodynamic screening length for rod networks. Experimental tracer-diffusion data are presented for various ionic strengths as a function of the rod concentration, both in the isotropic and nematic states. Orientational order parameters are measured for the same ionic strengths as a function of the rod concentration. The hydrodynamic screening length is determined from these experimental data and scaling relations obtained from the above mentioned theory. For the isotropic networks, a master curve is found for the hydrodynamic screening length as a function of the rod concentration. For the nematic networks the screening length turns out to be a very sensitive function of the orientational order parameter.     

Percolation of randomly centered rods and spheres

The percolation properties of randomly centered rods and spheres are studied. The approach is based on the detailed study of frequencies of cluster occurrences. For random rods, the analytic expressions are derived for all cluster frequencies. It is then shown that one-dimensional systems of random rods exhibit critical behaviour with? _c = ∞,γ = 1. For randomly centered spheres, we designed a numerical method for calculating the cluster frequencies. The approach is based on the principles of the Monte Carlo method. It can cope with clusters containing up to seven particles, which should suffice for the evaluation of accurate values of critical density and critical exponents.     

Rheology of contcentrated suspensions of slender rods

The rheology og suspensions of slender rods is discussed for various concentration regimes (i.e. dilute, semi-dilute, concentrated and liquid crystalline regimes) on the basis of the general kinetic theory of Brownian particles.     

Depletion potentials in colloidal mixtures of hard spheres and rods

The depletion potential between a hard sphere and a planar hard wall, or two hard spheres, imposed by suspended rigid spherocylindrical rods is computed by the acceptance ratio method through the application of Monte Carlo simulation. The accurate results and ideal-gas approximation results of the depletion potential are determined with the acceptance ratio method in our simulations. For comparison, the depletion potentials are also studied by using both the density functional theory and Derjaguin approximations. The density profile as a function of positions and orientations of rods, used in the density functional theory, is calculated by Monte Carlo simulation. The potential obtained by the acceptance ratio method is in good agreement with that of density functional theory under the ideal-gas approximation. The comparison between our results and those of other theories suggests that the acceptance ratio method is the only efficient method used to compute the depletion potential induced by nonspherical colloids with the volume fraction beyond the ideal-gas approximation.     

Tracer diffusion in colloidal suspensions under dilute and crowded conditions with hydrodynamic interactions

We consider tracer diffusion in colloidal suspensions under solid loading conditions, where hydrodynamic interactions play an important role. To this end, we carry out computer simulations based on the hybrid stochastic rotation dynamics-molecular dynamics (SRD-MD) technique. Many details of the simulation method are discussed in detail. In particular, our choices for the SRD-MD parameters and for the different scales are adapted to simulating colloidal suspensions under realistic conditions. Our simulation data are compared with published theoretical, experimental and numerical results and compared to Brownian dynamics simulation data. We demonstrate that our SRD-MD simulations reproduce many features of the hydrodynamics in colloidal fluids under finite loading. In particular, finite-size effects and the diffusive behavior of colloids for a range of volume fractions of the suspension show that hydrodynamic interactions are correctly included within the SRD-MD technique.     

Isotropic-nematic phase transition of nonaqueous suspensions of natural clay rods

A novel model system for studying the behavior of hard colloidal rods is presented, consisting of sterically stabilized particles of natural sepiolite clay. Electron microscopy and scattering results confirmed that the organophilic clay particles were individual, rigid rods when dispersed in organic solvents. With a length-to-diameter ratio of approximately 27, the particles showed nematic ordering for volume fractions phi > 0.06. Polarizing microscopy revealed that the phase separation process involved nucleation, growth, and coalescence of nematic domains. The phase volumes and particle concentrations in the coexisting phases were determined. The dependence of these quantities on the total concentration of the suspension agrees well with Onsager's [Ann. N. Y. Acad. Sci. 51, 627 (1949)] isotropic-nematic phase transition theory extended to bidisperse and polydisperse rod systems, and with previous experimental results for rigid rodlike particles. Particle size distributions were obtained by analyzing transmission electron microscopy images. A significant fractionation with respect to rod length (but not diameter) was observed in the coexisting isotropic and nematic phases. The relative polydispersity of both daughter phases was distinctly smaller than that of the parent suspension. The phase behavior of these daughter fractions agrees well with the predictions for hard spherocylinders of corresponding aspect ratios. An isotropic-nematic-nematic phase equilibrium was seen to develop in phase separated samples after 1 month standing and is ascribed to the effect of polydispersity and possibly gravity. The second nematic phase appearing is dominated by very long rods.     

Sedimentation velocity and potential in concentrated suspensions of charged porous spheres

The body-force-driven migration in a homogeneous suspension of polyelectrolyte molecules or charged flocs in an electrolyte solution is analyzed. The model used for the particle is a porous sphere in which the density of the hydrodynamic frictional segments, and therefore also that of the fixed charges, is constant. The effects of particle interactions are taken into account by employing a unit cell model. The overlap of the electric double layers of adjacent particles is allowed and the relaxation effect in the double layer surrounding each particle is considered. The electrokinetic equations which govern the electrostatic potential profile, the ionic concentration (or electrochemical potential energy) distributions, and the fluid velocity field inside and outside the porous particle in a unit cell are linearized by assuming that the system is only slightly distorted from equilibrium. Using a regular perturbation method, these linearized equations are solved for a symmetrically charged electrolyte with the density of the fixed charges as the small perturbation parameter. An analytical expression for the settling velocity of the charged porous sphere is obtained from a balance among its gravitational, electrostatic, and hydrodynamic forces. A closed-form formula for the sedimentation potential in a suspension of identical charged porous spheres is also derived by using the requirement of zero net electric current. The dependence of the sedimentation velocity and potential of the suspension on the particle volume fraction and other properties of the particle-solution system is found to be quite complicated.     

Diffusion and viscosity in a crowded environment: from nano- to macroscale

Although water is the chief component of living cells, food, and personal care products, the supramolecular components make their viscosity larger than that of water by several orders of magnitude. Using fluorescence correlation spectroscopy (FCS), photon correlation spectroscopy (PCS), NMR, and rheology data, we show how the viscosity changes from the value for water at the molecular scale to the large macroviscosity. We determined the viscosity experienced by nanoprobes (of sizes from 0.28 to 190 nm) in aqueous micellar solution of hexaethylene-glycol-monododecyl-ether (in a range of concentration from 0.1% w/w to 35% w/w) and identified a clear crossover at the length scale of 17 +/- 2 nm (slightly larger than persistence length of micelles) at which viscosity acquires its macroscopic value. The sharp dependence of the viscosity coefficients on the size of the probe in the nanoregime has important consequences for diffusion-limited reactions in crowded environments (e.g., living cells).     

Nanonecklaces assembled from gold rods, spheres, and bipyramids

Gold nanorods, nanospheres, and bipyramids have been assembled using glutathione and cysteine into three types of necklace structures, which might be useful for the fabrication of nanoscale photonic, electronic, and optoelectronic devices.     

Diffusion and residence time of hydrogen peroxide and water in crowded protein environments

Reactive oxygen species (ROS) have important functions in cell signaling and, when present at overly high levels, may cause oxidation of important biological molecules. Kinetic models to study diffusion of ROS inside of mitochondria often assume dynamics similar to that in solution. However, it is well-known that separation of proteins in the cytosol or inside of mitochondria, where ROS are most predominant, can be smaller than 1 nm. Diffusion of small molecules can be better regarded as a percolation process. In this article, we report results of diffusivity and residence of water and hydrogen peroxide in the proximity of proteins. In carrying out this study, we found some issues with the conventional way of computing residence times by means of survival time correlation functions. The main problem is that particles remaining on the surface of a protein for long times and for which one has very poor statistics contribute significantly to the short time behavior of the survival time correlation function. We mathematically describe this problem and propose methodology to overcome it.     

Glassy dynamics and kinetic vitrification of isotropic suspensions of hard rods

A nonlinear Langevin equation (NLE) theory for the translational center-of-mass dynamics of hard nonspherical objects has been applied to isotropic fluids of rigid rods. The ideal kinetic glass transition volume fraction is predicted to be a monotonically decreasing function beyond an aspect ratio of two. The functional form of the decrease is weaker than the inverse aspect ratio. Vitrification occurs at lower volume fractions for corrugated tangent bead rods compared to their smooth spherocylinder analogs. The ideal glass transition signals a crossover to activated dynamics, which is estimated to be observable before the nematic phase boundary is encountered if the aspect ratio is less than roughly 25. Calculations of the glassy elastic shear modulus and absolute yield stress reveal a roughly exponential growth with volume fraction. The dependence of entropic barriers and mean barrier hopping times on concentration for rods of variable aspect ratios can be collapsed quite well based on a difference volume fraction variable that quantifies the distance from the ideal glass boundary. Full numerical solution of the NLE theory via stochastic trajectory simulation was performed for tangent bead rods, and the results were compared to their hard sphere analogs. With increasing shape anisotropy the characteristic length scales of the nonequilibrium free energy increase and the magnitude of the localization well and entropic barrier curvatures decreases. These changes result in a significant aspect ratio dependence of dynamical properties and time correlation functions including weaker intermediate time subdiffusive transport, stronger two-step decay of the incoherent dynamic structure factor, longer mean alpha relaxation time, and stronger wavevector-dependent decoupling of relaxation times and the self-diffusion constant. The theoretical results are potentially testable via computer simulation, confocal microscopy, and dynamic light scattering.     

Entropic forces and directed alignment of hard squares in suspensions of rods and disks

We use Monte Carlo simulations in two dimensions to study the depletion forces between two hard squares in a suspension of hard rods or disks. We determine the effects of size and concentration of rods and disks on the potential of mean force between the squares. Both rods and disks produce a short-range depletion attraction between the two squares. The depletion interaction can be strong enough to outweigh the (rotational) entropic repulsion between the squares at certain sizes and concentrations of the rods and disks. We also probe the relative orientation that two squares adopt as they approach each other and we observe rich behavior, in which the relative orientation depends on the size, concentration, and shape of the depletion agent. Simple models based on the ideas of Asakura and Oosawa [J. Chem. Phys. 22, 1255 (1954)] can explain trends in the potentials of mean force obtained from the simulations.     

On the influence of a patterned substrate on crystallization in suspensions of hard spheres

We present a computer simulation study on crystal nucleation and growth in supersaturated suspensions of mono-disperse hard spheres induced by a triangular lattice substrate. The main result is that compressed substrates are wet by the crystalline phase (the crystalline phase directly appears without any induction time), while for stretched substrates we observe heterogeneous nucleation. The shapes of the nucleated crystallites fluctuate strongly. In the case of homogeneous nucleation amorphous precursors have been observed [T. Schilling et al., Phys. Rev. Lett. 105(2), 025701 (2010)]. For heterogeneous nucleation we do not find such precursors. The fluid is directly transformed into highly ordered crystallites.     

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.     

Percolation, phase separation, and gelation in fluids and mixtures of spheres and rods

The relationship between kinetic arrest, connectivity percolation, structure and phase separation in protein, nanoparticle, and colloidal suspensions is a rich and complex problem. Using a combination of integral equation theory, connectivity percolation methods, nai?ve mode coupling theory, and the activated dynamics nonlinear Langevin equation approach, we study this problem for isotropic one-component fluids of spheres and variable aspect ratio rigid rods, and also percolation in rod-sphere mixtures. The key control parameters are interparticle attraction strength and its (short) spatial range, total packing fraction, and mixture composition. For spherical particles, formation of a homogeneous one-phase kinetically stable and percolated physical gel is predicted to be possible, but depends on non-universal factors. On the other hand, the dynamic crossover to activated dynamics and physical bond formation, which signals discrete cluster formation below the percolation threshold, almost always occurs in the one phase region. Rods more easily gel in the homogeneous isotropic regime, but whether a percolation or kinetic arrest boundary is reached first upon increasing interparticle attraction depends sensitively on packing fraction, rod aspect ratio and attraction range. Overall, the connectivity percolation threshold is much more sensitive to attraction range than either the kinetic arrest or phase separation boundaries. Our results appear to be qualitatively consistent with recent experiments on polymer-colloid depletion systems and brush mediated attractive nanoparticle suspensions.     

Diffusion in macromolecular crowded media: Monte Carlo simulation of obstructed diffusion vs. FRAP experiments

The diffusion of tracer particles in 3D macromolecular crowded media has been studied using two methodologies, simulation and experimental, with the aim of comparing their results. First, the diffusion of a tracer in an obstructed 3D lattice with mobile and big size obstacles has been analyzed through a Monte Carlo (MC) simulation procedure. Secondly, fluorescence recovery after photobleaching (FRAP) experiments have been carried out to study the diffusion of a model protein (alpha-chymotrypsin) in in vitro crowded solution where two type of Dextran molecules are used as crowder agents. To facilitate the comparison, the relative size between the tracer and the crowder is the same in both studies. The results indicate a qualitative agreement between the diffusional behaviors observed in the two studies. The dependence of the anomalous diffusion exponent and the limiting diffusion coefficient on the obstacle size and excluded volume shows, in both cases, a similar tendency. The introduction of a reduced mobility parameter in the simulation model accounting for the short-range tracer?Cobstacle interactions allows obtaining a quantitative agreement between the limiting diffusion coefficient values yielded by both procedures. The simulation?Cexperiment quantitative agreement for the anomalous diffusion exponent requires further improvements. As far as we know, this is the first reported work where both techniques are used in parallel to study the diffusion in macromolecular crowded media.     

Theory of dielectric relaxations due to the interfacial polarization for two-component suspensions of spheres

A theoretical formula of dielectric relaxation in a form of complex relative permittivity is derived for dilute suspensions of spherical particles of two kinds on the basis of the Maxwell-Wagner theory of interfacial polarization. Another theoretical formula is derived further for concentrated suspensions of spheres of two kinds on condition that the formula derived above holds for the infinitesimally increasing process in concentration of the dispersed spheres. Furthermore a theoretical formula is derived for concentrated suspensions of shelled spheres of two kinds as the extension of the formula for concentrated suspensions. By use of the theoretical formulas proposed, values of the permittivities and the conductivities of the two-component suspensions were calculated for some examples with different sets of phase parameters. Results of the numerical calculation demonstrates dielectric relaxation profiles full of variety and characteristic of the suspensions containing two kinds of spheres covered with or without shells.     

Phase diagram of alloy crystal in the exhaustively deionized suspensions of binary mixtures of colloidal spheres

Phase diagrams of liquidlike, alloy crystal-like and amorphous solid-like(AS) structures have been obtained for the exhaustively deionized aqueous suspensions of the binary mixtures of polystyrene or silica spheres. Diameter, polydispersity index (standard deviation of diameter divided by the mean diameter) and size ratio of the binary spheres (diameter of small sphere divided by that of large one) range from 85 to 136 nm, 0.07 to 0.26 and 0.76 to 0.93, respectively. Close-up color photographs of the alloy crystals are taken and the crystal structure has been analysed from reflection spectroscopy. Most of the alloy crystals aresubstitutional solid-solution (sss) type and body-contered cubic lattice structure. Formation of the alloy crystals is attributed to the important role of the expanded electrical double layers in the deionized condition and increase toward unity in the effective size ratio, which is the effective diameter of small sphere including double layer divided by that of large sphere AS structure is formed at the rather high concentrations of two spheres, where the thickness of the electrical double layer is thin and the effective size rado is comparatively small.     

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.