Arrested phase separation in a short-ranged attractive colloidal system: a numerical study

We numerically investigate the competition between phase separation and dynamical arrest in a colloidal system interacting via a short-ranged attractive potential. Equilibrium fluid configurations are quenched at two different temperatures below the critical temperature and followed during their time evolution. At the lowest studied T, the phase-separation process is interrupted by the formation of an attractive glass in the dense phase. At the higher T, no arrest is observed and the phase-separation process proceeds endlessly in the simulated time window. The final structure of the glass retains memory of the interrupted phase-separation process in the form of a frozen spinodal decomposition peak, whose location and amplitude is controlled by the average packing fraction. We also discuss the time evolution of the nonergodicity parameter, providing evidence of a progressively decreasing localization length on increasing the packing fraction. Finally, we confirm that the reported results are independent of the microscopic dynamics.     

Aging in short-ranged attractive colloids: a numerical study

We study the aging dynamics in a model for dense simple liquids, in which particles interact through a hard-core repulsion complemented by a short-ranged attractive potential, of the kind found in colloidal suspensions. In this system, at large packing fractions, kinetically arrested disordered states can be created both on cooling (attractive glass) and on heating (repulsive glass). The possibility of having two distinct glasses, at the same packing fraction, with two different dynamics offers the unique possibility of comparing-within the same model-the differences in aging dynamics. We find that, while the aging dynamics of the repulsive glass is similar to the one observed in atomic and molecular systems, the aging dynamics of the attractive glass shows novel unexpected features.     

Dynamics in the presence of attractive patchy interactions

We report extensive Monte Carlo and event-driven molecular dynamics simulations of a liquid composed of particles interacting via hard-sphere interactions complemented by four tetrahedrally coordinated short-range attractive ("sticky") spots, a model introduced several years ago by Kolafa and Nezbeda (Kolafa, J.; Nezbeda, I. Mol. Phys. 1987, 87, 161). To access the dynamic properties of the model, we introduce and implement a new event-driven molecular dynamics algorithm suited to study the evolution of hard bodies interacting, beside the repulsive hard-core, with a short-ranged interpatch square well potential. We evaluate the thermodynamic properties of the model in deep supercooled states, where the bond network is fully developed, providing evidence of density anomalies. Different from models of spherically symmetric interacting particles, the liquid can be supercooled without encountering the gas-liquid spinodal in a wide region of packing fractions phi. Around an optimal phi, a stable fully connected tetrahedral network of bonds develops. By analyzing the dynamics of the model we find evidence of anomalous behavior: around the optimal packing, dynamics accelerate on both increasing and decreasing phi. We locate the shape of the isodiffusivity lines in the (phi - T) plane and establish the shape of the dynamic arrest line in the phase diagram of the model. Results are discussed in connection with colloidal dispersions of sticky particles and gel-forming proteins and their ability to form dynamically arrested states.     

Viscoelasticity and Stokes-Einstein relation in repulsive and attractive colloidal glasses

We report a numerical investigation of the viscoelastic behavior in models for steric repulsive and short-ranged attractive colloidal suspensions, along different paths in the attraction strength vs packing fraction plane. More specifically, we study the behavior of the viscosity (and its frequency dependence) on approaching the repulsive glass, the attractive glass, and in the reentrant region where viscosity shows a nonmonotonic behavior on increasing attraction strength. On approaching the glass lines, the increase of the viscosity is consistent with a power-law divergence with the same exponent and critical packing fraction previously obtained for the divergence of the density fluctuations. Based on mode-coupling calculations, we associate the increase of the viscosity with specific contributions from different length scales. We also show that the results are independent of the microscopic dynamics by comparing Newtonian and Brownian simulations for the same model. Finally, we evaluate the Stokes-Einstein relation approaching both glass transitions, finding a clear breakdown which is particularly strong for the case of the attractive glass.     

Equation of state and structural properties of the Weeks-Chandler-Andersen fluid

Molecular dynamics simulations have been carried out for the equation of state and percolation properties of the Weeks-Chandler-Andersen (WCA) system in its fluid phase as functions of density and temperature. The compressibility factor Z collapses well for the various isotherms, using an effective particle diameter for the WCA particle which is (in the usual WCA reduced units) sigma(e)=2(16)(1+T)(16), where T is the temperature. A corresponding "effective" packing fraction is zeta(e)=pisigma(e) (3)N6V, for N particles in volume V, which therefore scales out the effects of temperature. Using zeta(e) the simulation derived Z can be fitted to a simple analytic form which is similar to the Carnahan-Starling hard sphere equation of state and which is valid at all temperatures and densities where the WCA fluid is thermodynamically stable. The data, however, are not scalable onto the hard sphere equation of state for the complete packing fraction range. We explored the continuum percolation behavior of the WCA fluids. The percolation distance sigma(p) for the various states collapses well onto a single curve when plotted as sigma(p)sigma(e) against zeta(e). The ratio sigma(p)sigma(e) exhibits a monotonic decrease with increasing zeta(e) between the percolation line for permeable spheres and the glass transition limit, where sigma(p)sigma(e) approximately 1. The percolation packing fraction was calculated as a function of effective packing fraction and fitted to an empirical expression. The local coordination number at the percolation threshold showed a transition between the soft core and hard core limits from ca. 2:74 to 1:5, as previously demonstrated in the literature for true hard spheres. A number of simple analytic expressions that represent quite well the percolation characteristics of the WCA system are proposed.     

Dynamical arrest, percolation, gelation, and glass formation in model nanoparticle dispersions with thermoreversible adhesive interactions

We report an experimental study of the dynamical arrest transition for a model system consisting of octadecyl coated silica suspended in n-tetradecane from dilute to concentrated conditions spanning the state diagram. The dispersion's interparticle potential is tuned by temperature affecting the brush conformation leading to a thermoreversible model system. The critical temperature for dynamical arrest, T*, is determined as a function of dispersion volume fraction by small-amplitude dynamic oscillatory shear rheology. We corroborate this transition temperature by measuring a power-law decay of the autocorrelation function and a loss of ergodicity via fiber-optic quasi-elastic light scattering. The structure at T* is measured using small-angle neutron scattering. The scattering intensity is fit to extract the interparticle pair-potential using the Ornstein-Zernike equation with the Percus-Yevick closure approximation, assuming a square-well interaction potential with a short-range interaction (1% of particle diameter). (1) The strength of attraction is characterized using the Baxter temperature (2) and mapped onto the adhesive hard sphere state diagram. The experiments show a continuous dynamical arrest transition line that follows the predicted dynamical percolation line until ? ≈ 0.41 where it subtends the predictions toward the mode coupling theory attractive-driven glass line. An alternative analysis of the phase transition through the reduced second virial coefficient B(2)* shows a change in the functional dependence of B(2)* on particle concentration around ? ≈ 0.36. We propose this signifies the location of a gel-to-glass transition. The results presented herein differ from those observed for depletion flocculated dispersion of micrometer-sized particles in polymer solutions, where dynamical arrest is a consequence of multicomponent phase separation, suggesting dynamical arrest is sensitive to the physical mechanism of attraction.     

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.     

The ideal glass transition of hard spheres

We use the replica method to study the ideal glass transition of a liquid of identical hard spheres. We obtain estimates of the configurational entropy in the liquid phase, of the Kauzmann packing fraction phi(K), in the range of 0.58-0.62, and of the random close packing density phi(c), in the range of 0.64-0.67, depending on the approximation we use for the equation of state of the liquid. We also compute the pair-correlation function in the glassy states (i.e., dense amorphous packings) and we find that the mean coordination number at phi(c) is equal to 6. All these results compare well with numerical simulations and with other existing theories.     

Slow dynamics in a primitive tetrahedral network model

We report extensive Monte Carlo and event-driven molecular dynamics simulations of the fluid and liquid phase of a primitive model for silica recently introduced by Ford et al. [J. Chem. Phys. 121, 8415 (2004)]. We evaluate the isodiffusivity lines in the temperature-density plane to provide an indication of the shape of the glass transition line. Except for large densities, arrest is driven by the onset of the tetrahedral bonding pattern and the resulting dynamics is strong in Angell's classification scheme [J. Non-Cryst. Solids 131-133, 13 (1991)]. We compare structural and dynamic properties with corresponding results of two recently studied primitive models of network forming liquids-a primitive model for water and an angular-constraint-free model of four-coordinated particles-to pin down the role of the geometric constraints associated with bonding. Eventually we discuss the similarities between "glass" formation in network forming liquids and "gel" formation in colloidal dispersions of patchy particles.     

Single-filament dynamics and long-range ordering of semiflexible biopolymers under flow and confinement

We report the collective and single-filament dynamics of long semiflexible actin filaments flowing in an evaporating droplet adhering on glass and accumulating along the physical barrier constituted by the droplet triple line. The observation of fluorescent reporter filaments embedded in the entangled network enables us to relate the final collective organization of the accumulated filaments to the individual filament dynamics. Three areas corresponding to distinct filament organizations are observed in the region of the initial triple line pinning, after complete evaporation of the droplet. A nematic liquid-crystal-like alignment of the filaments is observed at the edge of the droplet because of the dynamic filament alignment, whereas a less-ordered packing is generated because of the bending and folding of most of the filaments. The latter unconventional dynamics is analyzed in terms of the amplification of undulation modes typical of semiflexible polymers. The receding regime of the droplet triple line leads finally to a remaining film of actin filaments showing random organization.     

Self-diffusion coefficients and shear viscosity of inverse power fluids: from hard- to soft-spheres

Molecular dynamics computer simulation has been used to compute the self-diffusion coefficient, D, and shear viscosity, eta(s), of soft-sphere fluids, in which the particles interact through the soft-sphere or inverse power pair potential, phi(r) = epsilon(sigma/r)(n), where n measures the steepness or stiffness of the potential, and epsilon and sigma are a characteristic energy and distance, respectively. The simulations were carried out on monodisperse systems for a range of n values from the hard-sphere (n --> infinity) limit down to n = 4, and up to densities in excess of the fluid-solid co-existence value. A new analytical procedure is proposed which reproduces the transport coefficients at high densities, and can be used to extrapolate the data to densities higher than accurately accessible by simulation or experiment, and tending to the glass transition. This formula, DX(c-1) proportional, variant A/X + B, where c is an adjustable parameter, and X is either the packing fraction or the pressure, is a development of one proposed by Dymond. In the expression, -A/B is the value of X at the ideal glass transition (i.e., where D and eta(s)(-1) --> 0). Estimated values are presented for the packing fraction and the pressure at the glass transition for n values between the hard and soft particle limits. The above expression is also shown to reproduce the high density viscosity data of supercritical argon, krypton and nitrogen. Fits to the soft-sphere simulation transport coefficients close to solid-fluid co-existence are also made using the analytic form, ln(D) = alpha(X)X, and n-dependence of the alpha(X) is presented (X is either the packing fraction or the pressure).     

Flow properties of freshly prepared ettringite suspensions in water at 25 degrees C

The rheology of a complex, heterogeneous mineral colloid was rationalised using models devised for model rod systems. Mixing a calcium hydroxide slurry with an aluminium sulphate solution produces a suspension of rod-shaped ettringite particles. Ettringite rod suspensions exhibit non-Newtonian flow behaviour, which depends on the shape of the particles, their size distribution, concentration and surface properties as well as the suspension medium characteristics. We have measured the shear viscosity of suspensions of ettringite rods with a median aspect ratio, r(i) approximately 8, at 25 degrees C as a function of particle volume fraction, phi, in the range 0.0001-0.08. It was found that the viscosity of the suspensions increased with phi, and showed a marked change of slope at phi approximately 0.01, which we identified as the minimum overlap concentration phi(*). Above phi(*), the system is in the semi-dilute regime. At phi>phi(*), when Pe(rot)>1, hydrodynamic interactions between rods become increasingly significant, and we observe shear-thinning behaviour. The high effective hydrodynamic volume of rotating rods, resulting in much lower values of the maximum packing fraction, phi(c), than for spheres, dominates the rheological behaviour of ettringite suspensions.     

Phase and glass transitions in short-range central potential model systems: the case of C60

Extensive molecular dynamics simulations show that a short-range central potential, suited to model C60, undergoes a high temperature transition to a glassy phase characterized by the positional disorder of the constituent particles. Crystallization, melting, and sublimation, which also take place during the simulation runs, are illustrated in detail. It turns out that vitrification and the mentioned phase transitions occur when the packing fraction of the system-defined in terms of an effective hard-core diameter-equals that of hard spheres at their own glass and melting transition, respectively. A close analogy also emerges between our findings and recent mode coupling theory calculations of structural arrest lines in a similar model of protein solutions. We argue that the conclusions of the present study might hold for a wide class of potentials currently employed to mimic interactions in complex fluids (some of which are of biological interest), suggesting how to achieve at least qualitative predictions of vitrification and crystallization in those systems.     

Colloid Vibration Potential in a Concentrated Suspension of Spherical Colloidal Particles

A relation between the dynamic electrophoretic mobility of spherical colloidal particles in a concentrated suspension and the colloid vibration potential (CVP) generated in the suspension by a sound wave is obtained from the analogy with the corresponding Onsager relation between electrophoretic mobility and sedimentation potential in concentrated suspensions previously derived on the basis of Kuwabara's cell model. The obtained expression for CVP is applicable to the case where the particle zeta potential is low, the particle relative permittivity is very small, and the overlapping of the electrical double layers of adjacent particles is negligible. It is found that CVP shows much stronger dependence on the particle volume fraction φ than predicted from the φ dependence of the dynamic electrophoretic mobility. It is also suggested that the same relation holds between the electrokinetic sonic amplitude of a concentrated suspension of spherical colloidal particles and the dynamic electrophoretic mobility. Copyright 1999 Academic Press.     

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.     

One-dimensional cluster growth and branching gels in colloidal systems with short-range depletion attraction and screened electrostatic repulsion

We report extensive numerical simulations of a simple model for charged colloidal particles in suspension with small nonadsorbing polymers. The chosen effective one-component interaction potential is composed of a short-range attractive part complemented by a Yukawa repulsive tail. We focus on the case where the screening length is comparable to the particle radius. Under these conditions, at low temperature, particles locally cluster into quasi one-dimensional aggregates which, via a branching mechanism, form a macroscopic percolating gel structure. We discuss gel formation and contrast it with the case of longer screening lengths, for which previous studies have shown that arrest is driven by the approach to a Yukawa glass of spherical clusters. We compare our results with recent experimental work on charged colloidal suspensions (Phys. Rev. Lett. 2005, 94, 208301).     

Particle-stabilized emulsions comprised of solid droplets

We kinetically stabilize oil-in-water emulsions comprising paraffin crystals by adsorbing solid particles (silica) of colloidal size at the oil/water interface. We obtain a set of emulsions that are quiescently stable for a long period of time (months), while the same emulsions are destabilized after only a few hours in the presence of surfactant molecules alone. The emulsions are submitted to a shear stress in order to probe their stability under flow conditions. Partial coalescence and gelation occur when the shear is applied for a sufficiently long period of time. The experiments reveal the existence of a critical droplet mass fraction, phi*, that defines a sharp transition between slow and fast gelation. The process of gelation is rather slow for phi < phi*, occurring at the scale of hours, and becomes almost instantaneous above phi*.     

Binary Mixtures in Linear Convection Arrays

We numerically investigated the dynamics of a mixture of finite-size active and passive disks in a linear array of two-dimensional convection rolls. The interplay of advection and steric interactions produces a number of interesting effects, like the stirring of a passive colloidal fluid by a small fraction of slow active particles, or the separation of the mixture active and passive colloidal fractions by increasing the motility of the active one, which eventually clusters in stagnation areas along the array walls. These mechanisms are quantitatively characterized by studying the dependence of the diffusion constants of the active and passive particles on the parameters of the active mixture fraction.     

Thermodynamic interaction parameters for the system water/NMMO hydrate

Vapor pressures of water were measured for aqueous solutions of N-methyl-morpholine N-oxide (NMMO) at 80, 90 and 100 degrees C. The Flory-Huggins interaction parameters, chi, calculated from these data as a function of phi, the volume fraction of NMMO, are negative at all concentrations; at low phi, they decrease by more than a factor of 2 as T is raised, whereas they remain almost unchanged as phi approaches unity. Accordingly, the heat of mixing is pronouncedly endothermal at low NMMO concentrations but close to athermal at low water content. The composition dependence of chi can be equally well described by the Redlich-Kister equation and by an approach subdividing the mixing process into two separate steps. The opportunities of the latter modeling for a better molecular understanding of the mixing processes are discussed.