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.     

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.     

Collisions, caging, thermodynamics, and jamming in the barrier hopping theory of glassy hard sphere fluids

An ultralocal limit of the microscopic single particle barrier hopping theory of glassy dynamics is proposed which allows explicit analytic expressions for the characteristic length scales, energy scales, and nonequilibrium free energy to be derived. All properties are shown to be controlled by a single coupling constant determined by the fluid density and contact value of the radial distribution function. This parameter quantifies an effective mean square force exerted on a tagged particle due to collisions with its surroundings. The analysis suggests a conceptual basis for previous surprising findings of multiple inter-relationships between characteristics of the transient localized state, the early stages of cage escape, non-Gaussian or dynamic heterogeneity effects, and the barrier hopping process that defines the alpha relaxation event. The underlying physical picture is also relevant to fluids of nonspherical molecules and sticky colloidal suspensions. The possibility of a unified view of liquid dynamics is suggested spanning the range from dense gases to the zero mobility jammed state.     

A lattice model of vitrification and gelation

In this work we introduce a simple lattice model with T-shaped molecules in two dimensions that exhibits a rich range of morphological behaviors. Depending on the volume fraction and quench path, this system can adopt uniform liquid, solution, and phase-separated states, as well as inhomogeneous glass or gel-like states, as revealed by dynamic mean-field simulations. An important characteristic of this system is the existence of a large number of degenerate low-energy states with small barriers that leads to a broad, kinetically explored landscape. The mean-field stability and phase diagram of this model is constructed and provides a useful guide for understanding the complex behaviors of the system. One striking feature is that there is a cascade of instabilities that converge to mark the onset of what we identify as the glass transition. Both dynamic mean-field and Monte Carlo simulations reveal glass-like relaxation dynamics. Our results lead to a picture of gelation as a continuation of the glass transition into the two-phase region, or equivalently, as an incomplete phase separation arrested by the onset of the glass transition.     

Melting and vitrification of Lennard-Jones spheres

Melting and vitrification processes of different types of mixtures of Lennard-Jones spheres are considered; mechanisms and geometrical models of transition from the solid to liquid state are defined and established.     

On the impossibility of defining adhesive hard spheres as sticky limit of a hard-sphere-Yukawa potential

For fluids of molecules with short-ranged hard-sphere-Yukawa (HSY) interactions, it is proven that the Noro-Frenkel "extended law of corresponding states" cannot be applied down to the vanishing attraction range, since the exact HSY second virial coefficient diverges in such a limit. It is also shown that, besides Baxter's original approach, a fully correct alternative definition of "adhesive hard spheres" can be obtained by taking the vanishing-range-limit (sticky limit) not of a Yukawa tail, as is commonly done, but of a slightly different potential with a logarithmic-Yukawa attraction.     

Phase behavior of weakly polydisperse sticky hard spheres: perturbation theory for the Percus-Yevick solution

We study the effects of size polydispersity on the gas-liquid phase behavior of mixtures of sticky hard spheres. To achieve this, the system of coupled quadratic equations for the contact values of the partial cavity functions of the Percus-Yevick solution [R. J. Baxter, J. Chem. Phys. 49, 2770 (1968)] is solved within a perturbation expansion in the polydispersity, i.e., the normalized width of the size distribution. This allows us to make predictions for various thermodynamic quantities which can be tested against numerical simulations and experiments. In particular, we determine the leading order effects of size polydispersity on the cloud curve delimiting the region of two-phase coexistence and on the associated shadow curve; we also study the extent of size fractionation between the coexisting phases. Different choices for the size dependence of the adhesion strengths are examined carefully; the Asakura-Oosawa model [J. Chem. Phys. 22, 1255 (1954)] of a mixture of polydisperse colloids and small polymers is studied as a specific example.     

Scaled particle theory of mixtures of hard spheres with negatively non-additive diameters

The scaled particle theory is solved for mixtures of hard spheres with arbitrary negative non-additive diameters in three dimensions. The results obtained with the Gibbs—Duhem relation show excellent agreement with the Monte-Carlo calculations. The results from the virial relation for high densities and high non-additivity are less satisfactory.     

Phase diagram of mixtures of hard colloidal spheres and discs: a free-volume scaled-particle approach

Phase diagrams of mixtures of colloidal hard spheres with hard discs are calculated by means of the free-volume theory. The free-volume fraction available to the discs is determined from scaled-particle theory. The calculations show that depletion induced phase separation should occur at low disc concentrations in systems now experimentally available. The gas-liquid equilibrium of the spheres becomes stable at comparable size ratios as with bimodal mixtures of spheres or mixtures of rods and spheres. Introducing finite thickness of the platelets gives rise to a significant lowering of the fluid branch of the binodal.     

A new generalization of the Carnahan-Starling equation of state to additive mixtures of hard spheres

We introduce an expansion of the equation of state for additive hard-sphere mixtures in powers of the total packing fraction with coefficients which depend on a set of weighted densities used in scaled particle theory and fundamental measure theory. We demand that the mixture equation of state recovers the quasiexact Carnahan-Starling [J. Chem. Phys. 51, 635 (1969)] result in the case of a one-component fluid and show from thermodynamic considerations and consistency with an exact scaled particle relation that the first and second orders of the expansion lead unambiguously to the Boublik-Mansoori-Carnahan-Starling-Leland [J. Chem. Phys. 53, 471 (1970); J. Chem. Phys. 54, 1523 (1971)] equation and the extended Carnahan-Starling equation introduced by Santos et al. [Mol. Phys. 96, 1 (1999)]. In the third order of the expansion, our approach allows us to define a new equation of state for hard-sphere mixtures which we find to be more accurate than the former equations when compared to available computer simulation data for binary and ternary mixtures. Using the new mixture equation of state, we calculate expressions for the surface tension and excess adsorption of the one-component fluid at a planar hard wall and compare its predictions to available simulation data.     

Stacking in sediments of colloidal hard spheres

We use computer simulations to investigate the crystallization dynamics of sedimenting hard spheres in large systems (hundreds of thousands of particles). We show that slow sedimentation results primarily in face-centered cubic (fcc) stacked crystals, instead of random hexagonal close packed or hexagonal close packed (hcp) crystals. We also find slanted stacking faults, in the fcc regions. However, we attribute the formation of fcc to the free energy difference between fcc and hcp and not to the presence of these slanted stacking faults. Although the free energy difference between hcp and fcc per particle is small (only 10(-3) times the thermal energy), it can become considerable, when multiplied by the number of particles in each domain. The ratio of fcc to hcp obtained from dynamic simulations is in excellent agreement with well-equilibrated Monte Carlo simulations, in which no slanted stacking faults were found. Our results explain a range of experiments on colloids, in which the amount of fcc increases upon lowering the sedimentation rate or decreasing the initial volume fraction.     

Enthalpy of mixing for polymer solutions based on Gibbs excess function limit of hard spheres mixtures

Enthalpy of mixing of polymer with solvent has been evaluated by the new polymer/solvent theory proposed by the authors in a previous article. The new theory was based on the excess Gibbs function limit of hard sphere mixtures with infinite size difference. The calculated enthalpy of mixing for polymer/solvent mixtures by the new theory, agreed with experimental data with good accuracy and indicated that the theory is capable to produce the enthalpy of mixing in the whole concentration range of polymer compared with Flory–Huggins theory. Also the calculations provided information on the studied polymer chains and the molecular interaction effects which were consistent with the properties of polymers and solvents used in the mixtures.     

Metastable states in the system of hard spheres

Approximate equations of the theory of liquids were applied to analyze the special features of the structural behavior of the system of hard spheres in the region of densities at which crystals are stable. Correlation functions for the liquid metastable state were obtained. The correlation functions of crystal-like metastable states that can exist under the same conditions as liquid metastable states are also considered. The reliability of the results is substantiated by a thermodynamic comparison with the data obtained using analytic equations of state.     

Freezing of hard spheres confined in narrow cylindrical pores

Monte Carlo simulations for the equation of state and phase behavior of hard spheres confined inside very narrow hard tubes are presented. For pores whose radii are greater than 1.1 hard sphere diameters, a sudden change in the density and the microscopic structure of the fluid is neatly observed, indicating the onset of freezing. In the high-density structure the particles rearrange in such a way that groups of three particles fit in sections across the pore.     

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.     

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.     

On the relation between virial coefficients and the close-packing of hard disks and hard spheres

The question of whether the known virial coefficients are enough to determine the packing fraction η(∞) at which the fluid equation of state of a hard-sphere fluid diverges is addressed. It is found that the information derived from the direct Pade? approximants to the compressibility factor constructed with the virial coefficients is inconclusive. An alternative approach is proposed which makes use of the same virial coefficients and of the equation of state in a form where the packing fraction is explicitly given as a function of the pressure. The results of this approach both for hard-disk and hard-sphere fluids, which can straightforwardly accommodate higher virial coefficients when available, lends support to the conjecture that η(∞) is equal to the maximum packing fraction corresponding to an ordered crystalline structure.     

On the caging number of two- and three-dimensional hard spheres

Local structural arrest in random packings of colloidal or granular spheres is quantified by a caging number, defined as the average minimum number of randomly placed spheres on a single sphere that immobilize all its translations. We present an analytic solution for the caging number for two-dimensional hard disks immobilized by neighbor disks which are placed at random positions under the constraint of a nonoverlap condition. Immobilization of a disk with radius r = 1 by arbitrary larger neighbor disks with radius r > or = 1 is solved analytically, whereas for contacting neighbors with radius 0 < r < 1, the caging number can be evaluated accurately with an approximate excluded volume model that also applies to spheres in higher Euclidean dimension. Comparison of our exact two-dimensional caging number with studies on random disk packing indicates that it relates to the average coordination number of random loose packing, whereas the parking number is more indicative for coordination in random dense packing of disks.