No Arabic abstract
Colloidal model systems allow studying crystallization kinetics under fairly ideal conditions with rather well characterized pair interactions and minimized external influences. In complementary approaches therefore experiment, analytic theory and simulation have been employed to study colloidal solidification in great detail. These studies were based on advanced optical methods, careful system characterization and sophisticated numerical methods. Both the effects of the type, strength and range of the pair-interaction between the colloidal particles and those of the colloid-specific polydispersity were addressed in a quantitative way. Key parameters of crystallization were derived and compared to those of metal systems. These systematic investigations significantly contributed to an enhanced understanding of the crystallization processes in general. Further, new fundamental questions have arisen and (partially) been solved over the last decade including e.g. a two step nucleation mechanism in homogeneous nucleation, choice of the crystallization pathway or the subtle interplay of boundary conditions in heterogeneous nucleation. On the other side, via the application of both gradients and external fields the competition between different nucleation and growth modes can be controlled and the resulting micro-structure be influenced. The present review attempts an account of the interesting developments occurred since the turn of the millennium and an identification of important novel trends with particular focus on experimental aspects.
Microscopic relaxation timescales are estimated from the autocorrelation functions obtained by dynamic light scattering experiments for Laponite suspensions with different concentrations ($C_{L}$), added salt concentrations ($C_{S}$) and temperatures ($T$). It has been shown in an earlier work [Soft Matter, 10, 3292-3300 (2014)] that the evolutions of relaxation timescales of colloidal glasses can be compared with molecular glass formers by mapping the waiting time ($t_{w}$) of the former with the inverse of thermodynamic temperature ($1/T$) of the latter. In this work, the fragility parameter $D$, which signifies the deviation from Arrhenius behavior, is obtained from fits to the time evolutions of the structural relaxation timescales. For the Laponite suspensions studied in this work, $D$ is seen to be independent of $C_{L}$ and $C_{S}$, but is weakly dependent on $T$. Interestingly, the behavior of $D$ corroborates the behavior of fragility in molecular glass formers with respect to equivalent variables. Furthermore, the stretching exponent $beta$, which quantifies the width $w$ of the spectrum of structural relaxation timescales is seen to depend on $t_{w}$. A hypothetical Kauzmann time $t_{k}$, analogous to the Kauzmann temperature for molecular glasses, is defined as the timescale at which $w$ diverges. Corresponding to the Vogel temperature defined for molecular glasses, a hypothetical Vogel time $t^{infty}_{alpha}$ is also defined as the time at which the structural relaxation time diverges. Interestingly, a correlation is observed between $t_{k}$ and $t^{infty}_{alpha}$, which is remarkably similar to that known for fragile molecular glass formers. A coupling model that accounts for the $t_{w}$-dependence of the stretching exponent is used to analyse and explain the observed correlation between $t_{k}$ and $t^{infty}_{alpha}$.
We report on a large scale computer simulation study of crystal nucleation in hard spheres. Through a combined analysis of real and reciprocal space data, a picture of a two-step crystallization process is supported: First dense, amorphous clusters form which then act as precursors for the nucleation of well-ordered crystallites. This kind of crystallization process has been previously observed in systems that interact via potentials that have an attractive as well as a repulsive part, most prominently in protein solutions. In this context the effect has been attributed to the presence of metastable fluid-fluid demixing. Our simulations, however, show that a purely repulsive system (that has no metastable fluid-fluid coexistence) crystallizes via the same mechanism.
We study the collective dynamics of colloidal suspensions in the presence of a time-dependent potential, by means of dynamical density functional theory. We consider a non-linear diffusion equation for the density and show that spatial patterns emerge from a sinusoidal external potential with a time-dependent wavelength. These patterns are characterized by a sinusoidal density with the average wavelength and a Bessel-function envelope with an induced wavelength that depends only on the amplitude of the temporal oscillations. As a generalization of this result, we propose a design strategy to obtain a family of spatial patterns using time-dependent potentials of practically arbitrary shape.
The Dynamic Monte Carlo (DMC) method is an established molecular simulation technique for the analysis of the dynamics in colloidal suspensions. An excellent alternative to Brownian Dynamics or Molecular Dynamics simulation, DMC is applicable to systems of spherical and/or anisotropic particles and to equilibrium or out-of-equilibrium processes. In this work, we present a theoretical and methodological framework to extend DMC to the study of heterogeneous systems, where the presence of an interface between coexisting phases introduces an additional element of complexity in determining the dynamic properties. In particular, we simulate a Lennard-Jones fluid at the liquid-vapor equilibrium and determine the diffusion coefficients in the bulk of each phase and across the interface. To test the validity of our DMC results, we also perform Brownian Dynamics simulations and unveil an excellent quantitative agreement between the two simulation techniques.
We present a comprehensive study of the slip and flow of concentrated colloidal suspensions using cone-plate rheometry and simultaneous confocal imaging. In the colloidal glass regime, for smooth, non-stick walls, the solid nature of the suspension causes a transition in the rheology from Herschel-Bulkley (HB) bulk flow behavior at large stress to a Bingham-like slip behavior at low stress, which is suppressed for sufficient colloid-wall attraction or colloid-scale wall roughness. Visualization shows how the slip-shear transition depends on gap size and the boundary conditions at both walls and that partial slip persist well above the yield stress. A phenomenological model, incorporating the Bingham slip law and HB bulk flow, fully accounts for the behavior. Microscopically, the Bingham law is related to a thin (sub-colloidal) lubrication layer at the wall, giving rise to a characteristic dependence of slip parameters on particle size and concentration. We relate this to the suspensions osmotic pressure and yield stress and also analyze the influence of van der Waals interaction. For the largest concentrations, we observe non-uniform flow around the yield stress, in line with recent work on bulk shear-banding of concentrated pastes. We also describe residual slip in concentrated liquid suspensions, where the vanishing yield stress causes coexistence of (weak) slip and bulk shear flow for all measured rates.