Effect of polydispersity on the crystallization kinetics of suspensions of colloidal hard spheres when approaching the glass transition

Phase behavior of binary mixtures of hollow and regular microgels

Soft colloids are widely used to study glass transition, aging and jamming. A high size polydispersity is typically introduced in these systems to avoid crystal formation. Here, we use binary mixtures of hollow and regular microgels with comparable sizes to inhibit crystallization. The phase behavior of the mixture is probed as a function of the number fraction of hollow microgels and characterized by small-angle X-ray scattering. Molecular dynamic simulations are used to extract the particle-particle pair potential and obtain insight on their deformation. The results suggest that the high deformability of the hollow microgels offers an alternative route to maximize the entropy without crystal formation.

In Situ X-ray Scattering Reveals Coarsening Rates of Superlattices Self-Assembled from Electrostatically Stabilized Metal Nanocrystals Depend Nonmonotonically on Driving Force

Self-assembly of colloidal nanocrystals (NCs) into superlattices (SLs) is an appealing strategy to design hierarchically organized materials with promising functionalities. Mechanistic studies are still needed to uncover the design principles for SL self-assembly, but such studies have been difficult to perform due to the fast time and short length scales of NC systems. To address this challenge, we developed an apparatus to directly measure the evolving phases in situ and in real time of an electrostatically stabilized Au NC solution before, during, and after it is quenched to form SLs using small-angle X-ray scattering. By developing a quantitative model, we fit the time-dependent scattering patterns to obtain the phase diagram of the system and the kinetics of the colloidal and SL phases as a function of varying quench conditions. The extracted phase diagram is consistent with particles whose interactions are short in range relative to their diameter. We find the degree of SL order is primarily determined by fast (subsecond) initial nucleation and growth kinetics, while coarsening at later times depends nonmonotonically on the driving force for self-assembly. We validate these results by direct comparison with simulations and use them to suggest dynamic design principles to optimize the crystallinity within a finite time window. The combination of this measurement methodology, quantitative analysis, and simulation should be generalizable to elucidate and better control the microscopic self-assembly pathways of a wide range of bottom-up assembled systems and architectures.

Characterization and Optimization of Fluorescent Organosilica Colloids for 3D Confocal Microscopy Prepared Under "Zero-Flow" Conditions

We optimize and characterize the preparation of 3-trimethoxysilyl propyl methacrylate (TPM) colloidal suspensions for three-dimensional confocal microscopy. We revisit a simple synthesis of TPM microspheres by nucleation of droplets from prehydrolyzed TPM oil in a "zero-flow" regime and demonstrate how precise and reproducible control of particle size may be achieved via single-step nucleation with a focus on how the reagents are mixed. We also revamp the conventional dyeing method for TPM particles to achieve uniform transfer of a fluorophore to the organosilica droplets, improving particle identification. Finally, we illustrate how a ternary mixture of tetralin, trichloroethylene, and tetrachloroethylene may be used as a suspension medium which matches the refractive index of these particles while allowing independent control of the density mismatch between particle and solvent.

Choice of diameters in a polydisperse model glassformer: Deterministic or stochastic?

In particle-based computer simulations of polydisperse glassforming systems, the particle diameters σ=σ_{1},⋯,σ_{N} of a system with N particles are chosen with the intention to approximate a desired distribution density f with the corresponding histogram. One method to accomplish this is to draw each diameter randomly and independently from the density f. We refer to this stochastic scheme as model S. Alternatively, one can apply a deterministic method, assigning an appropriate set of N values to the diameters. We refer to this method as model D. We show that, for sample-to-sample fluctuations, especially for the glassy dynamics at low temperatures, it matters whether one chooses model S or model D. Using molecular dynamics computer simulations, we investigate a three-dimensional polydisperse nonadditive soft-sphere system with f(s)∼s^{-3}. The swap Monte Carlo method is employed to obtain equilibrated samples at very low temperatures. We show that for model S the sample-to-sample fluctuations due to the quenched disorder imposed by the diameters σ can be explained by an effective packing fraction. Dynamic susceptibilities in model S can be split into two terms: one that is of thermal nature and can be identified with the susceptibility of model D, and another one originating from the disorder in σ. At low temperatures the latter contribution is the dominating term in the dynamic susceptibility. Our study clarifies the pros and cons of the use of models S and D in practice.

Waiting-time dependent non-equilibrium phase diagram of simple glass- and gel-forming liquids

Under numerous circumstances, many soft and hard materials are present in a puzzling wealth of non-equilibrium amorphous states, whose properties are not stationary and depend on preparation. They are often summarized in unconventional "phase diagrams" that exhibit new "phases" and/or "transitions" in which time, however, is an essential variable. This work proposes a solution to the problem of theoretically defining and predicting these non-equilibrium phases and their time-evolving phase diagrams, given the underlying molecular interactions. We demonstrate that these non-equilibrium phases and the corresponding non-stationary (i.e., aging) phase diagrams can indeed be defined and predicted using the kinetic perspective of a novel non-equilibrium statistical mechanical theory of irreversible processes. This is illustrated with the theoretical description of the transient process of dynamic arrest into non-equilibrium amorphous solid phases of an instantaneously quenched simple model fluid involving repulsive hard-sphere plus attractive square well pair interactions.

Assembly Control at a Low Péclet Number in Ultracentrifugation for Uniformly Sized Nanoparticles

The intrinsic high diffusion rate of colloids at low Péclet number results in an extremely fast crystallization process and instant formation of colloidal crystals, even at an ultracentrifugal field of extremely high intensity. By introducing a small number of clusters in sedimention, it should be possible to slow down the crystallization process, thus making the assembly order tunable in preparative ultracentrifugation experiments. Here, we used sodium dodecyl sulfate-stabilized polystyrene nanoparticles (with a size dispersity of 1.07) dispersed in a solution of high ionic strength. Sedimentation and assembly of these nanoparticles were done using preparative ultracentrifugation at various angular velocities. The sedimentation process was also analyzed in situ by analytical ultracentrifugation in real time. By creating as low as 3% of clusters into these nearly uniformly sized polystyrene nanoparticle dispersions during the sedimentation process, the superstructure order becomes easily tunable between glassy and crystalline. Theoretical calculations complemented the experiments to explain the mechanism of cluster formation in sedimentation. This work provides a novel methodology to produce superstructures with a tunable packing order for colloids at low Péclet number.

"Dense diffusion" in colloidal glasses: short-ranged long-time self-diffusion as a mechanistic model for relaxation dynamics

Despite decades of exploration of the colloidal glass transition, mechanistic explanation of glassy relaxation processes has remained murky. State-of-the-art theoretical models of the colloidal glass transition such as random first order transition theory, active barrier hopping theory, and non-equilibrium self-consistent generalized Langevin theory assert that relaxation reported at volume fractions above the ideal mode coupling theory prediction φg,MCT requires some sort of activated process, and that cooperative motion plays a central role. However, discrepancies between predicted and measured values of φg and ambiguity in the role of cooperative dynamics persist. Underlying both issues is the challenge of conducting deep concentration quenches without flow and the difficulty in accessing particle-scale dynamics. These two challenges have led to widespread use of fitting methods to identify divergence, but most a priori assume divergent behavior; and without access to detailed particle dynamics, it is challenging to produce evidence of collective dynamics. We address these limitations by conducting dynamic simulations accompanied by experiments to quench a colloidal liquid into the putative glass by triggering an increase in particle size, and thus volume fraction, at constant particle number density. Quenches are performed from the liquid to final volume fractions 0.56 ≤ φ ≤ 0.63. The glass is allowed to age for long times, and relaxation dynamics are monitored throughout the simulation. Overall, correlated motion acts to release dynamics from the glassy plateau - but only over length scales much smaller than a particle size - allowing self-diffusion to re-emerge; self-diffusion then relaxes the glass into an intransient diffusive state, which persists for φ < 0.60. We observe similar relaxation dynamics up to φ = 0.63 before achieving the intransient state. We find that this long-time self-diffusion is short-ranged: analysis of mean-square displacement reveals a glassy cage size a fraction of a particle size that shrinks with quench depth, i.e. increasing volume fraction. Thus the equivalence between cage size and particle size found in the liquid breaks down in the glass, which we confirm by examining the self-intermediate scattering function over a range of wave numbers. The colloidal glass transition can hence be viewed mechanistically as a shift in the long-time self-diffusion from long-ranged to short-ranged exploration of configurations. This shift takes place without diverging dynamics: there is a smooth transition as particle mobility decreases dramatically with concomitant emergence of a dense local configuration space that permits sampling of many configurations via local particle motion.

Dramatic slowing of compositional relaxations in the approach to the glass transition for a bimodal colloidal suspension

Molecular dynamics simulation was used to study a model colloidal suspension with two species of slightly different sized colloidal particles in an explicit solvent. In this work we calculated the four interdiffusion coefficients for the ternary system, which were then used to calculate the decay coefficients D_{±} of the two independent diffusive modes. We found that the slower D_{-} decay mode, which is associated with the system's ability to undergo compositional changes, was responsible for the long-time decay in the intermediate scattering function. We also found that a decrease in D_{-} to negligible values at a packing fraction of Φ_{g}=0.592 resulted in an extreme slow-down in the long-time decay of the intermediate scattering function often associated with the glass transition. Above Φ_{g}, the system formed a long-lived metastable state that did not relax to its equilibrium crystal state within the simulation time window. We concluded that the inhibition of crystallization was caused by the inability of the quenched fluid to undergo the compositional changes needed for the formation of the equilibrium crystal.

Introduction to the Colloidal Glass Transition

Colloids are suspensions of small solid particles in a liquid and exhibit glassy behavior when the particle concentration is high. In these samples, the particles are roughly analogous to individual molecules in a traditional glass. This model system has been used to study the glass transition since the 1980s. In this Viewpoint I summarize some of the intriguing behaviors of the glass transition in colloids and discuss open questions.

Dynamics of highly polydisperse colloidal suspensions as a model system for bacterial cytoplasm

There are various kinds of macromolecules in bacterial cell cytoplasm. The size polydispersity of the macromolecules is so significant that the crystallization and the phase separation could be suppressed, thus stabilizing the liquid state of bacterial cytoplasm. On the other hand, recent experiments suggested that the macromolecules in bacterial cytoplasm should exhibit glassy dynamics, which should be also affected significantly by the size polydispersity of the macromolecules. In this work, we investigate the anomalous and slow dynamics of highly polydisperse colloidal suspensions, of which size distribution is chosen to mimic Escherichia coli cytoplasm. We find from our Langevin dynamics simulations that the diffusion coefficient (D_{tot}) and the displacement distribution functions (P(r,t)) averaged over all colloids of different sizes do not show anomalous and glassy dynamic behaviors until the system volume fraction ϕ is increased up to 0.82. This indicates that the intrinsic polydispersity of bacterial cytoplasm should suppress the glass transition and help maintain the liquid state of the cytoplasm. On the other hand, colloids of each kind show totally different dynamic behaviors depending on their size. The dynamics of colloids of different size becomes non-Gaussian at a different range of ϕ, which suggests that a multistep glass transition should occur. The largest colloids undergo the glass transition at ϕ=0.65, while the glass transition does not occur for smaller colloids in our simulations even at the highest value of ϕ. We also investigate the distribution (P(θ,t)) of the relative angles of displacement for macromolecules and find that macromolecules undergo directionally correlated motions in a sufficiently dense system.

Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations

The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insights into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that, in the past few decades, have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state-of-the-art computational methods by reviewing simulations of such processes as ice nucleation and the crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insights into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that, as a result, the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that, by improving (i) existing interatomic potentials and (ii) currently available enhanced sampling methods, the community can move toward accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.

Dynamic decoupling and local atomic order of a model multicomponent metallic glass-former

The dynamics of multicomponent metallic alloys is spatially heterogeneous near glass transition. The diffusion coefficient of one component of the metallic alloys may also decouple from those of other components, i.e., the diffusion coefficient of each component depends differently on the viscosity of metallic alloys. In this work we investigate the dynamic heterogeneity and decoupling of a model system for multicomponent Pd43Cu27Ni10P20 melts by using a hard sphere model that considers the size disparity of alloys but does not take chemical effects into account. We also study how such dynamic behaviors would relate to the local atomic structure of metallic alloys. We find, from molecular dynamics simulations, that the smallest component P of multicomponent Pd43Cu27Ni10P20 melts becomes dynamically heterogeneous at a translational relaxation time scale and that the largest major component Pd forms a slow subsystem, which has been considered mainly responsible for the stabilization of amorphous state of alloys. The heterogeneous dynamics of P atoms accounts for the breakdown of Stokes-Einstein relation and also leads to the dynamic decoupling of P and Pd atoms. The dynamically heterogeneous P atoms decrease the lifetime of the local short-range atomic orders of both icosahedral and close-packed structures by orders of magnitude.

Crystallization seeds favour crystallization only during initial growth

Crystallization represents the prime example of a disorder-order transition. In realistic situations, however, container walls and impurities are frequently present and hence crystallization is heterogeneously seeded. Rarely the seeds are perfectly compatible with the thermodynamically favoured crystal structure and thus induce elastic distortions, which impede further crystal growth. Here we use a colloidal model system, which not only allows us to quantitatively control the induced distortions but also to visualize and follow heterogeneous crystallization with single-particle resolution. We determine the sequence of intermediate structures by confocal microscopy and computer simulations, and develop a theoretical model that describes our findings. The crystallite first grows on the seed but then, on reaching a critical size, detaches from the seed. The detached and relaxed crystallite continues to grow, except close to the seed, which now prevents crystallization. Hence, crystallization seeds facilitate crystallization only during initial growth and then act as impurities.

Probing planar defects in nanoparticle superlattices by 3D small-angle electron diffraction tomography and real space imaging

We demonstrate how the acquisition and processing of 3D electron diffraction data can be extended to characterize structural features on the mesoscale, and show how lattice distortions in superlattices of self-assembled spherical Pd nanoparticles can be quantified by three-dimensional small-angle electron diffraction tomography (3D SA-EDT). Transmission electron microscopy real space imaging and 3D SA-EDT reveal a high density of stacking faults that was related to a competition between fcc and hcp arrangements during assembly. Information on the orientation of the stacking faults was used to make analogies between planar defects in the superlattices and Shockley partial dislocations in metallic systems.

Crystallization kinetics of colloidal model suspensions: recent achievements and new perspectives

Colloidal model systems allow studying crystallization kinetics under fairly ideal conditions, with rather well-characterized pair interactions and minimized external influences. In complementary approaches 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. Over the last decade, both the effects of the type, strength and range of the pair-interaction between the colloidal particles and those of the colloid-specific polydispersity have been addressed in a quantitative way. Key parameters of crystallization have been 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, for example, 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 hand, via the application of both gradients and external fields the competition between different nucleation and growth modes can be controlled and the resulting microstructure be influenced. The present review attempts to cover the interesting developments that have occurred since the turn of the millennium and to identify important novel trends, with particular focus on experimental aspects.

Solidification of a colloidal hard sphere like model system approaching and crossing the glass transition

We investigated the process of vitrification and crystallization in a model system of colloidal hard spheres. The kinetics of the solidification process was measured using time resolved static light scattering, while the time evolution of the dynamic properties was determined using time resolved dynamic light scattering. By performing further analysis we confirm that solidification of hard sphere colloids is mediated by precursors. Analyzing the dynamic properties we can show that the long time dynamics and thus the shear rigidity of the metastable melt is highly correlated with the number density of solid clusters (precursors) nucleated. In crystallization these objects convert into highly ordered crystals whereas in the case of vitrification this conversion is blocked and the system is (temporarily) locked in the metastable precursor state. From the early stages of solidification one cannot clearly conclude whether the melt will crystallize or vitrify. Furthermore our data suggests that colloidal hard sphere glasses can crystallize via homogeneous nucleation.

Phonons in two-dimensional colloidal crystals with bond-strength disorder

We study phonon modes in two-dimensional colloidal crystals composed of soft microgel particles with hard polystyrene particle dopants distributed randomly on the triangular lattice. This experimental approach produces close-packed lattices of spheres with random bond strength disorder, i.e., the effective springs coupling nearest neighbors are very stiff, very soft, or of intermediate stiffness. Particle tracking video microscopy and covariance matrix techniques are then employed to derive the phonon modes of the corresponding "shadow" crystals with bond strength disorder as a function of increasing dopant concentration. At low frequencies, hard and soft particles participate equally in the phonon modes, and the samples exhibit Debye-like density of states behavior characteristic of crystals. For mid- and high-frequency phonons, the relative participation of hard versus soft particles in each mode is found to vary systematically with dopant concentration. Additionally, a few localized modes, primarily associated with hard particle motions, are found at the highest frequencies.

Spinodal fractionation in a polydisperse square-well fluid

Using kinetic Monte Carlo simulation, we model gas-liquid spinodal decomposition in a size-polydisperse square well fluid, representing a "near-monodisperse" colloidal dispersion. We find that fractionation (demixing) of particle sizes between the phases begins asserting itself shortly after the onset of phase ordering. Strikingly, the direction of size fractionation can be reversed by a seemingly trivial choice between two interparticle potentials which, in the monodisperse case, are identical--we rationalize this in terms of a perturbative, equilibrium theory of polydispersity. Furthermore, our quantitative results show that kinetic Monte Carlo simulation can provide detailed insight into the role of fractionation in real colloidal systems.

The physics of the colloidal glass transition

As one increases the concentration of a colloidal suspension, the system exhibits a dramatic increase in viscosity. Beyond a certain concentration, the system is said to be a colloidal glass; structurally, the system resembles a liquid, yet motions within the suspension are slow enough that it can be considered essentially frozen. For several decades, colloids have served as a valuable model system for understanding the glass transition in molecular systems. The spatial and temporal scales involved allow these systems to be studied by a wide variety of experimental techniques. The focus of this review is the current state of understanding of the colloidal glass transition, with an emphasis on experimental observations. A brief introduction is given to important experimental techniques used to study the glass transition in colloids. We describe features of colloidal systems near and in glassy states, including increases in viscosity and relaxation times, dynamical heterogeneity and ageing, among others. We also compare and contrast the glass transition in colloids to that in molecular liquids. Other glassy systems are briefly discussed, as well as recently developed synthesis techniques that will keep these systems rich with interesting physics for years to come.

Kinetics study of crystallization with the disorder-bcc-fcc phase transition of charged colloidal dispersions

Structure transformation (disorder-bcc-fcc) in charged colloidal dispersions, as a manifestation of the Ostwald's step rule, was confirmed by means of reflection spectrum (RS) measurements in our previous study. By taking advantage of a reflection spectrum containing plenty of information about the crystallization behaviors, time-dependent changes of parameters associated with the crystal structure and composition during the disorder-bcc-fcc transition are reported by treating the data from RS in this article. In addition, Avrami's model is adopted to analyze the transition process and investigate the transition rate. On the basis of the above investigations, associated kinetic features of crystallization with the disorder-bcc-fcc transition are described.

Crystallization in suspensions of hard spheres: a Monte Carlo and molecular dynamics simulation study

The crystallization of a metastable melt is one of the most important non-equilibrium phenomena in condensed matter physics, and hard sphere colloidal model systems have been used for several decades to investigate this process by experimental observation and computer simulation. Nevertheless, there is still an unexplained discrepancy between the simulation data and experimental nucleation rate densities. In this paper we examine the nucleation process in hard spheres using molecular dynamics and Monte Carlo simulation. We show that the crystallization process is mediated by precursors of low orientational bond-order and that our simulation data fairly match the experimental data sets.

Densest local sphere-packing diversity. II. Application to three dimensions

The densest local packings of N three-dimensional identical nonoverlapping spheres within a radius R(min)(N) of a fixed central sphere of the same size are obtained for selected values of N up to N=1054. In the predecessor to this paper [A. B. Hopkins, F. H. Stillinger, and S. Torquato, Phys. Rev. E 81, 041305 (2010)], we described our method for finding the putative densest packings of N spheres in d-dimensional Euclidean space R(d) and presented those packings in R(2) for values of N up to N=348. Here we analyze the properties and characteristics of the densest local packings in R(3) and employ knowledge of the R(min)(N), using methods applicable in any d, to construct both a realizability condition for pair correlation functions of sphere packings and an upper bound on the maximal density of infinite sphere packings. In R(3), we find wide variability in the densest local packings, including a multitude of packing symmetries such as perfect tetrahedral and imperfect icosahedral symmetry. We compare the densest local packings of N spheres near a central sphere to minimal-energy configurations of N+1 points interacting with short-range repulsive and long-range attractive pair potentials, e.g., 12-6 Lennard-Jones, and find that they are in general completely different, a result that has possible implications for nucleation theory. We also compare the densest local packings to finite subsets of stacking variants of the densest infinite packings in R(3) (the Barlow packings) and find that the densest local packings are almost always most similar as measured by a similarity metric, to the subsets of Barlow packings with the smallest number of coordination shells measured about a single central sphere, e.g., a subset of the fcc Barlow packing. Additionally, we observe that the densest local packings are dominated by the dense arrangement of spheres with centers at distance R(min)(N). In particular, we find two "maracas" packings at N=77 and N=93, each consisting of a few unjammed spheres free to rattle within a "husk" composed of the maximal number of spheres that can be packed with centers at respective R(min)(N).

Polymorphism, crystal nucleation and growth in the phase-field crystal model in 2D and 3D

We apply a simple dynamical density functional theory, the phase-field crystal (PFC) model of overdamped conservative dynamics, to address polymorphism, crystal nucleation, and crystal growth in the diffusion-controlled limit. We refine the phase diagram for 3D, and determine the line free energy in 2D and the height of the nucleation barrier in 2D and 3D for homogeneous and heterogeneous nucleation by solving the respective Euler-Lagrange (EL) equations. We demonstrate that, in the PFC model, the body-centered cubic (bcc), the face-centered cubic (fcc), and the hexagonal close-packed structures (hcp) compete, while the simple cubic structure is unstable, and that phase preference can be tuned by changing the model parameters: close to the critical point the bcc structure is stable, while far from the critical point the fcc prevails, with an hcp stability domain in between. We note that with increasing distance from the critical point the equilibrium shapes vary from the sphere to specific faceted shapes: rhombic dodecahedron (bcc), truncated octahedron (fcc), and hexagonal prism (hcp). Solving the equation of motion of the PFC model supplied with conserved noise, solidification starts with the nucleation of an amorphous precursor phase, into which the stable crystalline phase nucleates. The growth rate is found to be time dependent and anisotropic; this anisotropy depends on the driving force. We show that due to the diffusion-controlled growth mechanism, which is especially relevant for crystal aggregation in colloidal systems, dendritic growth structures evolve in large-scale isothermal single-component PFC simulations. An oscillatory effective pair potential resembling those for model glass formers has been evaluated from structural data of the amorphous phase obtained by instantaneous quenching. Finally, we present results for eutectic solidification in a binary PFC model.

Precursor-mediated crystallization process in suspensions of hard spheres

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.

Experimental study of random-close-packed colloidal particles

A collection of spherical particles can be packed tightly together into an amorphous packing known as "random close packing" (RCP). This structure is of interest as a model for the arrangement of molecules in simple liquids and glasses, as well as the arrangement of particles in sand piles. We use confocal microscopy to study the arrangement of colloidal particles in an experimentally realized RCP state. We image a large volume containing more than 450,000 particles with a resolution of each particle position to better than 0.02 particle diameters. While the arrangement of the particles satisfies multiple criteria for being random, we also observe a small fraction (less than 3%) of tiny crystallites (4 particles or fewer). These regions pack slightly better and are thus associated with locally higher densities. The structure factor of our sample at long length scales is nonzero, S(0)=0.049±0.008, suggesting that there are long wavelength density fluctuations in our sample. These may be due to polydispersity or tiny crystallites. Our results suggest that experimentally realizable RCP systems may be different from simulated RCP systems, in particular, with the presence of these long wavelength density fluctuations.

Optical experiments on a crystallizing hard-sphere-polymer mixture at coexistence

We report on the crystallization kinetics in an entropically attractive colloidal system using a combination of time resolved scattering methods and microscopy. Hard sphere particles are polystyrene microgels swollen in a good solvent (radius a=380 nm, starting volume fraction 0.534) with the short ranged attractions induced by the presence of short polymer chains (radius of gyration r g=3 nm, starting volume fraction 0.0224). After crystallization, stacking faulted face centered cubic crystals coexist with about 5% of melt remaining in the grain boundaries. From the Bragg scattering signal we infer the amount of crystalline material, the average crystallite size and the number density of crystals as a function of time. This allows to discriminate an early stage of conversion, followed by an extended coarsening stage. The small angle scattering (SALS) appears only long after completed conversion and exhibits Furukawa scaling for all times. Additional microscopic experiments reveal that the grain boundaries have a reduced Bragg scattering power but possess an increased refractive index. Fits of the Furukawa function indicate that the dimensionality of the scatterers decreases from 2.25 at short times to 1.65 at late times and the characteristic length scale is slightly larger than the average crystallite size. Together this suggests the SALS signal is due scattering from a foam like grain boundary network as a whole.

Observation of the disorder-induced crystal-to-glass transition

The role of frustration and quenched disorder in driving the transformation of a crystal into a glass is investigated in quasi-two-dimensional binary colloidal suspensions. Frustration is induced by added smaller particles. The crystal-glass transition is measured to differ from the liquid-glass transition in quantitative and qualitative ways. The crystal-glass transition bears structural signatures similar to those of the crystal-fluid transition: at the transition point, the persistence of orientational order decreases sharply from quasilong range to short range, and the orientational order susceptibility exhibits a maximum. The crystal-glass transition also features a sharp variation in particle dynamics: at the transition point, dynamic heterogeneity grows rapidly, and a dynamic correlation length scale increases abruptly.

Microscopic structural relaxation in a sheared supercooled colloidal liquid

The rheology of dense amorphous materials under large shear strain is not fully understood, partly due to the difficulty of directly viewing the microscopic details of such materials. We use a colloidal suspension to simulate amorphous materials and study the shear-induced structural relaxation with fast confocal microscopy. We quantify the plastic rearrangements of the particles in several ways. Each of these measures of plasticity reveals spatially heterogeneous dynamics, with localized regions where many particles are strongly rearranging by these measures. We examine the shapes of these regions and find them to be essentially isotropic, with no alignment in any particular direction. Furthermore, individual particles are equally likely to move in any direction other than the overall bias imposed by the strain.

Hard spheres: crystallization and glass formation

Motivated by old experiments on colloidal suspensions, we report molecular dynamics simulations of assemblies of hard spheres, addressing crystallization and glass formation. The simulations cover wide ranges of polydispersity s (standard deviation of the particle size distribution divided by its mean) and particle concentration. No crystallization is observed for s>0.07. For 0.02<s<0.07, we find that increasing the polydispersity at a given concentration slows down crystal nucleation. The main effect here is that polydispersity reduces the supersaturation since it tends to stabilize the fluid but to destabilize the crystal. At a given polydispersity (<0.07), we find three regimes of nucleation: standard nucleation and growth at concentrations in and slightly above the coexistence region; 'spinodal nucleation', where the free-energy barrier to nucleation appears to be negligible, at intermediate concentrations; and, at the highest concentrations, a new mechanism, still to be fully understood, which only requires small rearrangement of the particle positions. The cross-over between the second and third regimes occurs at a concentration, approximately 58 per cent by volume, where the colloid experiments show a marked change in the nature of the crystals formed and the particle dynamics indicate an 'ideal' glass transition.

Competition between heterogeneous and homogeneous nucleation near a flat wall

We studied the competition between heterogeneous and homogeneous nucleation of an aqueous suspension of charged colloidal spheres close to the container walls. Samples of equilibrium crystalline structure were shear-melted and the metastable melt left to solidify after the cessation of shear. The crystallization kinetics was monitored using time-resolved scattering techniques: at low particle number densities n we applied an improved static light scattering method while at large particle concentrations ultra-small-angle x-ray scattering was applied for the first time. Our results show some unexpected behavior: the heterogeneous nucleation at the container walls is delayed in comparison to the homogeneous bulk nucleation and its rate density appears surprisingly slightly smaller, demonstrating the complexity of the observed crystallization process.

Crystallization of hard-sphere glasses

We study by molecular dynamics the interplay between arrest and crystallization in hard spheres. For state points in the plane of volume fraction (0.54 <or= varphi <or= 0.63) and polydispersity (0 <or= s <or= 0.085), we delineate states that spontaneously crystallize from those that do not. For noncrystallizing (or precrystallization) samples we find isodiffusivity lines consistent with an ideal glass transition at varphi_{g} approximately 0.585, independent of s. Despite this, for s < 0.05, crystallization occurs at varphi > varphi_{g}. This happens on time scales for which the system is aging, and a diffusive regime in the mean square displacement is not reached; by those criteria, the system is a glass. Hence, contrary to a widespread assumption in the colloid literature, the occurrence of spontaneous crystallization within a bulk amorphous state does not prove that this state was an ergodic fluid rather than a glass.

Arrest of flow and emergence of activated processes at the glass transition of a suspension of particles with hard spherelike interactions

By combining aspects of the coherent and self-intermediate scattering functions, we show that the arrest of particle number density fluctuations spreads from the position of the main structure factor peak. We propose that this arrest impairs the system's ability to respond to diffusing momentum currents, leading to an enhanced resistance to flow. From the stretching of the coherent intermediate scattering functions in the glass, we read a manifestation of the undissipated thermal energy-the source of the ergodicity restoring processes that short-circuit the sharp transition to a perfect glass.

Unusual crystallization kinetics in a hard sphere colloid-polymer mixture

We investigated the crystallization kinetics of a hard sphere colloid-polymer mixture at conditions where about 95% of solid coexists with about 5% of fluid. From time resolved Bragg and small angle light scattering, we find that the crystallite size increases with a power law of exponent alpha approximately 1/3 during both the conversion and the coarsening stage. This observation points to a single conserved order parameter for both stages and cannot be explained if the mixture is regarded as an effective one-component system. We alternatively suggest that--based on local geometric demixing--the polymer density takes the role of the conserved order parameter.

Ripening-dominated crystallization in polydisperse hard-sphere-like colloids

We report on the crystal growth scenario in gravity-matched, polydisperse hard-sphere-like colloids at increasing particle concentration. In the fluid-crystal coexistence region, the crystal size as a function of time shows two separate regimes corresponding to crystal growth and crystal ripening. At higher supersaturation the crystal size grows according to the same power law through the whole experimental window of a few days: crystal growth and ripening merge together. We show that our observations cannot be explained by considering the slowing down of single-particle dynamics due to increasing volume fraction. We suggest that size fractionation occurring at the crystal-fluid interface is the dominant mechanism.

Dynamics and structure of an aging binary colloidal glass

We study aging in a colloidal suspension consisting of micron-sized particles in a liquid. This system is made glassy by increasing the particle concentration. We observe samples composed of particles of two sizes, with a size ratio of 1:2.1 and a volume fraction ratio 1:6, using fast laser scanning confocal microscopy. This technique yields real-time, three-dimensional movies deep inside the colloidal glass. Specifically, we look at how the size, motion, and structural organization of the particles relate to the overall aging of the glass. Particles move in spatially heterogeneous cooperative groups. These mobile regions tend to be richer in small particles, and these small particles facilitate the motion of nearby particles of both sizes.

Crystallization of dense binary hard-sphere mixtures with marginal size ratio

Molecular dynamics simulations are performed for binary hard-sphere mixtures with a size ratio of gamma =0.9 and a volume fraction of phi=0.58 over a range of compositions. We show how, at this high volume fraction, crystallization depends sensitively on the composition. Evidence is presented that crystallization in these mixtures does not proceed by the standard nucleation and growth paradigm. Rather, some crystallite forms almost immediately and then an interplay between compositional fluctuations and crystal growth is able to dramatically extend the time scale on which further crystallization occurs. This can be seen as a form of geometric frustration.

Reversible phase transition of colloids in a binary liquid solvent

We report fluid-fluid and fluid-solid phase transitions of charge-stabilized polystyrene particles suspended in a binary liquid mixture of 3-methylpyridine and water. These thermally reversible phase transitions occur in the homogeneous phase of the binary liquid mixture below the coexistence temperature of the two liquids. Close density matching of the particles and the solvent allows us to follow the phase behavior until complete coexistence of macroscopic phases with temperature as the control parameter. We use small angle x-ray scattering to characterize these phases as colloidal gas, liquid, fcc crystal, and glass.