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

Porous crystals in charged sphere suspensions by aggregate-driven phase separation

The kinetics of phase transition processes often governs the resulting material microstructure. Using optical microscopy, we here investigate the formation and stabilization of a porous crystalline microstructure forming in low-salt suspensions of charged colloidal spheres containing aggregates comprising some 5-10 of these colloids. We observe the transformation of an initially crystalline colloidal solid with homogeneously incorporated aggregates to individual, compositionally refined crystallites of perforated morphology coexisting with an aggregate-enriched fluid phase filling the holes and separating individual crystallites. A preliminary kinetic characterization suggests that the involved processes follow power laws. We show that this route to porous materials is neither restricted to nominally single component systems nor to a particular microstructure to start from. However, it necessitates an early rapid solidification stage during which the aggregates become trapped in the bulk of the host-crystals. The thermodynamic stability of the reconstructed crystalline scaffold against melting under increased salinity was found comparable to that of pure phase crystallites grown very slowly from a melt. Future implications of this novel route to porous colloidal crystals are discussed.

Polymorph Stability and Free Energy of Crystallization of Freely-Jointed Polymers of Hard Spheres

The free energy of crystallization of monomeric hard spheres as well as their thermodynamically stable polymorph have been known for several decades. In this work, we present semianalytical calculations of the free energy of crystallization of freely-jointed polymers of hard spheres as well as of the free energy difference between the hexagonal closed packed (HCP) and face-centered cubic (FCC) polymorphs. The phase transition (crystallization) is driven by an increase in translational entropy that is larger than the loss of conformational entropy of chains in the crystal with respect to chains in the initial amorphous phase. The conformational entropic advantage of the HCP polymer crystal over the FCC one is found to be ΔschHCP-FCC≈0.331×10-5k per monomer (expressed in terms of Boltzmann's constant k). This slight conformational entropic advantage of the HCP crystal of chains is by far insufficient to compensate for the larger translational entropic advantage of the FCC crystal, which is predicted to be the stable one. The calculated overall thermodynamic advantage of the FCC over the HCP polymorph is supported by a recent Monte Carlo (MC) simulation on a very large system of 54 chains of 1000 hard sphere monomers. Semianalytical calculations using results from this MC simulation yield in addition a value of the total crystallization entropy for linear, fully flexible, athermal polymers of Δs≈0.93k per monomer.

Self-Avoiding Random Walks as a Model to Study Athermal Linear Polymers under Extreme Plate Confinement

Monte Carlo (MC) simulations, built around chain-connectivity-altering moves and a wall-displacement algorithm, allow us to simulate freely-jointed chains of tangent hard spheres of uniform size under extreme confinement. The latter is realized through the presence of two impenetrable, flat, and parallel plates. Extreme conditions correspond to the case where the distance between the plates approaches the monomer size. An analysis of the local structure, based on the characteristic crystallographic element (CCE) norm, detects crystal nucleation and growth at packing densities well below the ones observed in bulk analogs. In a second step, we map the confined polymer chains into self-avoiding random walks (SAWs) on restricted lattices. We study all realizations of the cubic crystal system: simple, body centered, and face centered cubic crystals. For a given chain size (SAW length), lattice type, origin of SAW, and level of confinement, we enumerate all possible SAWs (equivalently all chain conformations) and calculate the size distribution. Results for intermediate SAW lengths are used to predict the behavior of long, fully entangled chains through growth formulas. The SAW analysis will allow us to determine the corresponding configurational entropy, as it is the driving force for the observed phase transition and the determining factor for the thermodynamic stability of the corresponding crystal morphologies.

Influence of Hydrodynamic Interactions on Colloidal Crystallization

One of the biggest unresolved problems in crystallization phenomena is the significant discrepancy in the nucleation rate between experiments and simulations even for the simplest liquid, i.e., the hard-sphere system. A popular explanation for this discrepancy is the neglect of hydrodynamic interactions (HI) in simulation studies. By comparing simulations with and without HI, we show that the long-time diffusive dynamics of the colloids is slowed down more rapidly by hydrodynamic lubrication effects with increasing volume fraction. We find that the kinetics of both nucleation and growth are controlled by this long-time diffusion and that it is possible to account for most of the effects of HI by rescaling with this timescale. Therefore, we conclude that HI is not the primary cause of the accelerated nucleation rates observed in experiments.

Dynamic Crystallization and Phase Transition in Evaporating Colloidal Droplets

Evaporating colloidal droplets are an omnipresent phenomenon in nature and engage in many scientific and commercial technologies. Despite their apparent importance, many of the fundamental aspects remain unknown, particularly the relationships between evaporation kinetics, volume fraction, crystallization, and phase transition. Here, we follow the structural evolution and drying dynamics across the liquid-to-solid transition of evaporating colloidal droplets containing polystyrene nanospheres with both spatial and temporal resolutions through the in situ small-angle X-ray scattering and ex situ electron microscopy techniques. We find the unconventional evaporation-driven heterogeneous crystallization and the sequential stacking of face-centered cubic (fcc), random hexagonal close-packed (rhcp), and random close-packed (rcp) superlattice structures. The crystallization and phase transition processes are further elucidated and coordinated with the real-time volume fraction variation, which constitutes a rich and dynamic picture of the self-assembly process. Starting with the Onsager principle, we provide quantitative analysis to the evaporation kinetics, including concentration gradient, gelation, and cavitation. Our findings impart a new mechanism of dynamic nucleation and crystallization and reveal the intimate link between structural heterogeneity and evaporation kinetics.

Formation of Laves phases in buoyancy matched hard sphere suspensions

Colloidal Laves phases (LPs) are promising precursors for photonic materials. Laves phases have not yet been observed to form in experiments on colloidal suspensions of hard spheres (HS), even though they have been reported in computer simulations. LP formation so far has been achieved only for binary mixtures of colloidal charged spheres or ligand-stabilized nano-particles after drying. Using static light scattering, we monitored LP formation and annealing in a binary mixture of buoyant hard sphere approximants (size ratio Γ = 0.77, number or molar fraction of small spheres xS = 0.76) for volume fractions in the fluid-crystal coexistence regions. All samples spontaneously formed MgZn2 type LPs on the time scale of weeks to months via bulk nucleation and growth. Irrespective of the initial suspension volume fractions, the LP volume fraction at coexistence is ΦCOEX = 0.59 which is significantly below the close packing limit ΦMAX = 0.615 and remarkably close to the expectation from simulation. At low volume fractions, crystals anneal to high quality during coarsening which is in line with recent theoretical expectations for the thermodynamic stability of different LP types. At large volume fractions, however, the diffractograms evolve towards a more MgCu2-like appearance which we attribute to the formation of randomly stacked LPs. Such structures are not known from atomic systems.

Domain size polydispersity effects on the structural and dynamical properties in lipid monolayers with phase coexistence

In lipid monolayers with phase coexistence, domains of the liquid-condensed phase always present size polydispersity. However, very few theoretical works consider size distribution effects on the monolayer properties. Because of the difference in surface densities, domains have excess dipolar density with respect to the surrounding liquid expanded phase, originating a dipolar inter-domain interaction. This interaction depends on the domain area, and hence the presence of a domain size distribution is associated with interaction polydispersity. Inter-domain interactions are fundamental to understanding the structure and dynamics of the monolayer. For this reason, it is expected that polydispersity significantly alters monolayer properties. By means of Brownian dynamics simulations, we study the radial distribution function (RDF), the average mean square displacement and the average time-dependent self-diffusion coefficient, D(t), of lipid monolayers with normally distributed size domains. For this purpose, we vary the relevant system parameters, polydispersity and interaction strength, within a range of experimental interest. We also analyze the consequences of using a monodisperse model to determine the interaction strength from an experimental RDF. We find that polydispersity strongly affects the value of the interaction strength, which is greatly underestimated if polydispersity is not considered. However, within a certain range of parameters, the RDF obtained from a polydisperse model can be well approximated by that of a monodisperse model, by suitably fitting the interaction strength, even for 40% polydispersities. For small interaction strengths or small polydispersities, the polydisperse systems obtained from fitting the experimental RDF have an average mean square displacement and D(t) in good agreement with that of the monodisperse system.

Size-Dependent Localization in Polydisperse Colloidal Glasses

We have investigated concentrated suspensions of polydisperse hard spheres and have determined the dynamics and sizes of individual particles using confocal microscopy. With increasing concentration, the dynamics of the small and large particles start to differ. The large particles exhibit slower dynamics and stronger localization. Moreover, as the particle size increases, the local volume fraction ϕ_{loc} also increases. In the glass state, the localization length significantly decreases beyond ϕ_{loc}≈0.67. This suggests a link between local crowding and dynamical heterogeneities. However dynamical arrest of subpopulations seems not directly linked to a large value of ϕ_{loc}, indicating the importance of collective effects.

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.

Colloidal crystal growth monitored by Bragg diffraction interference fringes

We monitored the crystal growth kinetics of crystallization of a shear melted crystalline colloidal array (CCA). The fcc CCA heterogeneously nucleates at the flow cell wall surface. We examined the evolution of the (1 1 1) Bragg diffraction peak, and, for the first time, quantitatively monitored growth by measuring the temporal evolution of the Bragg diffraction interference fringes. Modeling of the evolution of the fringe patterns exposes the time dependence of the increasing crystal thickness. The initial diffusion-driven linear growth is followed by ripening-driven growth. Between 80 and 90 microM NaCl concentrations the fcc crystals first linearly grow at rates between 1.9 and 4.2 microm/s until they contact homogeneously nucleated crystals in the bulk. At lower salt concentrations interference fringes are not visible because the strong electrostatic interactions between particles result in high activation barriers, preventing defect annealing and leading to a lower crystal quality. The fcc crystals melt to a liquid phase at >90 microM NaCl concentrations. Increasing NaCl concentrations slow the fcc CCA growth rate consistent with the expectation of the classical Wilson-Frenkel growth theory. The final thickness of wall-nucleated CCA, that is determined by the competition between growth of heterogeneously and homogenously nucleated CCA, increases with higher NaCl concentrations.

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.

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.