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.     

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.     

Glassy dynamics and mechanical response in dense fluids of soft repulsive spheres. I. Activated relaxation, kinetic vitrification, and fragility

The microscopic nonlinear Langevin equation theory of activated glassy dynamics is applied to dense fluids of spherical particles that interact via a finite range Hertzian contact soft repulsion. The activation barrier and mean alpha relaxation time are predicted to be rich functions of volume fraction and particle stiffness, exhibiting a non-monotonic variation with concentration at high volume fractions. The latter is due to a structural "soft jamming" crossover where the real space local cage order weakens when soft particles significantly overlap. The highly variable dependences of the relaxation time on temperature and volume fraction are reasonably well collapsed onto two distinct master curves that are qualitatively consistent with a recent scaling ansatz and computer simulation study. A kinetic vitrification diagram is constructed and compared to its dynamic crossover analog. Intersection of the dynamic crossover and soft jamming threshold boundaries occurs for particles that are sufficiently soft, implying the nonexistence of a clear activated dynamics regime or kinetic arrest transition for such particles. The isothermal dynamic fragility is predicted to vary over a wide range as a function of particle stiffness, and soft particles behave as strong glasses. Qualitative comparisons with simulations and microgel experiments reveal good agreement.     

Heterogeneous dynamics in columnar liquid crystals of parallel hard rods

In the wake of previous studies on the rattling-and-jumping diffusion in smectic liquid crystal phases of colloidal rods, we analyze here for the first time the heterogeneous dynamics in columnar phases. More specifically, we perform computer simulations to investigate the relaxation dynamics of a binary mixture of perfectly aligned hard spherocylinders. We detect that the columnar arrangement of the system produces free-energy barriers that the particles should overcome to jump from one column to another, thus determining a hopping-type diffusion. This phenomenon accounts for the non-Gaussian intercolumn diffusion and shows a two-step structural relaxation that is remarkably analogous to that of out-of-equilibrium glass-forming systems and gels. Surprisingly enough, slight deviations from the behavior of simple liquids due to transient cages is also observed in the direction perpendicular to this plane, where the system is usually referred to as liquidlike.     

Viscosity of hard sphere suspensions

A simple functional representation of the concentration dependence of the low-shear viscosity eta of hard sphere suspensions is proposed. The representation, which agrees with published literature at all volume fractions phi, has a hitherto-unremarked transition in its functional form at phi approximately 0.42 identical with phi(t). phi(t) is definitely less than the volume fraction 0.49 of the hard sphere melting transition.     

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.     

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.     

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.     

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.     

Structure and dynamics of polymer-grafted clay suspensions

The structure and dynamics of polymer-grafted two-dimensional silicate layers in solution were investigated. The geometry of the individual silicate layers was examined by looking at both polarized and depolarized light scattering from dilute solutions, while higher-concentration systems were used to study the interaction and dynamics of polymer-grafted silicate layers in suspension. The form factor for an oblate ellipsoid was used to fit the polarized intensity profile, and values of a approximately 80 nm and b approximately 380 nm for the semi-axes were obtained. The 80 nm value compares reasonably with the dimensions of the polymer brushes grafted on the surface of the silicate layers. The modulus of the grafted silicate in solution, as determined by Brillouin scattering, is of the order of 10 GPa. The cooperative diffusion mechanism, typical of interacting polymer chains, is suppressed due to the high polymer osmotic pressure. The osmotic pressure is also responsible for the weak interpenetration of the densely grafted polymer chains on the surface of the silicate layers. The scattering data indicates that the polymer-grafted nanoparticles move via collective diffusion and experience significant decrease in mobility above their overlap concentration.     

Isotropic-nematic transition of hard rods immersed in random sphere matrices

Using replica density functional theory and Monte Carlo computer simulations we investigate a system of annealed hard spherocylinders adsorbed in a matrix of quenched hard spheres. Theoretical predictions for the partition coefficient, defined as the ratio of density of rods in the matrix and that in a reservoir, agree well with simulation results. Theory predicts the isotropic-nematic transition to remain first order upon increasing sphere packing fraction, and to shift towards lower rod densities. This scenario is consistent with our simulation results that clearly show a jump in the nematic order parameter upon increasing the rod density at constant matrix packing fraction, corresponding to the isotropic-nematic transition, even for sphere matrix packing fractions < or approximately equal to 0.3.     

Glassy protein dynamics and gigantic solvent reorganization energy of plastocyanin

We report the results of molecular dynamics simulations of electron-transfer activation parameters of plastocyanin metalloprotein involved as an electron carrier in natural photosynthesis. We have discovered that slow, non-ergodic conformational fluctuations of the protein, coupled to hydrating water, result in a very broad distribution of donor-acceptor energy gaps far exceeding those observed for commonly studied inorganic and organic donor-acceptor complexes. The Stokes shift is not affected by these fluctuations and can be calculated from solvation models in terms of the linear response of the solvent dipolar polarization. The non-ergodic character of large-amplitude protein/water mobility breaks the strong link between the Stokes shift and the reorganization energy characteristic of equilibrium (ergodic) theories of electron transfer. This mechanism might be responsible for fast electronic transitions in natural electron-transfer proteins characterized by low reaction free energy.     

Dynamics of hard sphere suspensions using dynamic light scattering and X-ray photon correlation spectroscopy: dynamics and scaling of the intermediate scattering function

Intermediate scattering functions are measured for colloidal hard sphere systems using both dynamic light scattering and x-ray photon correlation spectroscopy. We compare the techniques, and discuss the advantages and disadvantages of each. Both techniques agree in the overlapping range of scattering vectors. We investigate the scaling behavior found by Segré and Pusey [Phys. Rev. Lett. 77, 771 (1996)] but challenged by Lurio et al. [Phys. Rev. Lett. 84, 785 (2000)]. We observe a scaling behavior over several decades in time but not in the long-time regime. Moreover, we do not observe long-time diffusive regimes at scattering vectors away from the peak of the structure factor and so question the existence of long-time diffusion coefficients at these scattering vectors.     

Activated hopping and dynamical fluctuation effects in hard sphere suspensions and fluids

Single particle Brownian dynamics simulation methods are employed to establish the full trajectory level predictions of our nonlinear stochastic Langevin equation theory of activated hopping dynamics in glassy hard sphere suspensions and fluids. The consequences of thermal noise driven mobility fluctuations associated with the barrier hopping process are determined for various ensemble-averaged properties and their distributions. The predicted mean square displacements show classic signatures of transient trapping and anomalous diffusion on intermediate time and length scales. A crossover to a stronger volume fraction dependence of the apparent nondiffusive exponent occurs when the entropic barrier is of order the thermal energy. The volume fraction dependences of various mean relaxation times and rates can be fitted by empirical critical power laws with parameters consistent with ideal mode-coupling theory. However, the results of our divergence-free theory are largely a consequence of activated dynamics. The experimentally measurable alpha relaxation time is found to be very similar to the theoretically defined mean reaction time for escape from the barrier-dominated regime. Various measures of decoupling have been studied. For fluid states with small or nonexistent barriers, relaxation times obey a simple log-normal distribution, while for high volume fractions the relaxation time distributions become Poissonian. The product of the self-diffusion constant and mean alpha relaxation time increases roughly as a logarithmic function of the alpha relaxation time. The cage scale incoherent dynamic structure factor exhibits nonexponential decay with a modest degree of stretching. A nearly universal collapse of the different volume fraction results occurs if time is scaled by the mean alpha relaxation time. Hence, time-volume fraction superposition holds quite well, despite the presence of stretching and volume fraction dependent decoupling associated with the stochastic barrier hopping process. The relevance of other origins of dynamic heterogeneity (e.g., mesoscopic domains), and comparison of our results with experiments, simulations, and alternative theories, is discussed.     

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.     

On melting dynamics and the glass transition. I. Glassy aspects of melting dynamics

The following properties are in the present literature associated with the behavior of supercooled glass-forming liquids: faster than exponential growth of the relaxation time, dynamical heterogeneities, growing point-to-set correlation length, crossover from mean-field behavior to activated dynamics. In this paper we argue that these properties are also present in a much simpler situation, namely the melting of the bulk of an ordered phase beyond a first order phase transition point. This is a promising path toward a better theoretical, numerical and experimental understanding of the above phenomena and of the physics of supercooled liquids. We discuss in detail the analogies and the differences between the glass and the bulk melting transitions.     

Crystalline ordering in lattice models of hard rods with nearest neighbour attraction

Several lattice models for hard rods with attraction between neighbouring parallel rodes are considered. It is proved that the models for sufficiently low temperature and high fugacity exhibit a phase transition to an ordered crystalline structure. Nothing is proved about the existence of a liquid crystal state.     

Phonons in suspensions of hard sphere colloids: volume fraction dependence

The propagation of sound waves in suspensions of hard sphere colloids is studied as a function of their volume fraction up to random close packing using Brillouin light scattering. The rich experimental phonon spectra of up to five phonon modes are successfully described by theoretical calculations based on the multiple scattering method. Two main types of phonon modes are revealed: Type A modes are acoustic excitations which set up deformations in both the solid (particles) and the liquid (solvent) phases; for type B modes the stress and strain are predominantly localized near the interface between the solid particles and the surrounding liquid (interface waves). While the former become harder (increase their effective sound velocity) as the particle volume fraction increases the latter become softer (the corresponding sound velocity decreases).