Two-stage crystallization of charged colloids under low supersaturation conditions

A minimal colloid model of solution crystallization nucleates crystals classically

A fundamental assumption of the classical theories of crystal nucleation is that the individual molecules from the "old" phase associate to an emerging nucleus individually and sequentially. Numerous recent studies of crystal nucleation in solution have revealed nonclassical pathways, whereby crystal nuclei are hosted and fed by amorphous clusters pre-formed in the solution. A sizable knowledge gap has persisted, however, in the definition of the molecular-level parameters that direct a solute towards classical or nonclassical nucleation. Here we construct a suspension of colloid particles of hydrodynamic diameter 1.1 μm and monitor their individual motions towards a quasi-two-dimensional crystal by scanning confocal microscopy. We combine electrostatic repulsion and polymer-induced attraction to obtain a simple isotropic pair interaction potential with a single attractive minimum of tunable depth between 1.2kBT and 2.7kBT. We find that even the smallest aggregates that form in this system structure as hexagonal two-dimensional crystals and grow and maturate by the association and exchange of single particles from the solution, signature behaviors during classical nucleation. The particles in the suspension equilibrate with those in the clusters and the volume fractions of suspensions at equilibrium correspond to straightforward thermodynamic predictions based on depth of the interparticle attraction. These results demonstrate that classical nucleation is selected by particles interacting with a minimal potential and present a benchmark for future modifications of the molecular interactions that may induce nonclassical nucleation.

Template-induced crystallization of charged colloids: a molecular dynamics study

By using a large enough number of particles and implementing a parallel algorithm on the CUDA platform, we have performed brute-force molecular dynamics simulations to study the template-induced heterogeneous crystallization in charged colloids. Six kinds of templates, whose patterns include the planes of fcc(100), fcc(110), fcc(111), bcc(100), bcc(110) and bcc(111), have been implanted into the middle of the simulation box. Except the fcc(111) template, whose structure benefits not only fcc but also hcp crystals resulting in a similar behavior to homogeneous crystallization, bcc-type templates favor the formation of bcc crystals and bcc-like precursors while fcc-type templates favor the formation of fcc crystals and fcc-like precursors. Therefore, for fcc(100) and fcc(110) templates, heterogeneous crystallization will definitely result in a fcc crystallite. However, the results of heterogeneous crystallization that are induced by bcc-type templates are subtly different at different state points. At the state points where the interaction strength of charged colloids is weak and the fcc phase is thermodynamically stable, the bcc crystals formed with the promotion of bcc-type templates are not stable so as to tend to transform into fcc or hcp crystals. When the interaction strength of charged colloids is high, the predominant bcc crystals formed with the promotion of bcc-type templates can always persist within the time scale of simulation although not bcc but fcc crystals are thermodynamically stable.

Kinetic Phase Diagram for Nucleation and Growth of Competing Crystal Polymorphs in Charged Colloids

We determine the kinetic phase diagram for nucleation and growth of crystal phases in a suspension of charged colloids. Exploiting the seeding approach in extensive simulations, we calculate nucleation barrier heights for face-centered cubic (fcc) and body-centered cubic (bcc) phases for varying screening lengths and supersaturations. We determine for the entire metastable fluid region the crystal polymorph with the lowest nucleation barrier, and find a regime close to the triple point where metastable bcc can form due to a lower nucleation barrier, even though fcc is the stable phase. For higher supersaturation, we find that the difference in barrier heights decreases and we observe a mix of hexagonal close-packed, fcc, and bcc structures in the growth of crystalline seeds as well as in spontaneously formed crystals. Our kinetic phase diagram rationalizes the different crystallization mechanisms observed in previous work.

Salt-concentration-dependent nucleation rates in low-metastability colloidal charged sphere melts containing small amounts of doublets

We determined bulk crystal nucleation rates in aqueous suspensions of charged spheres at low metastability. Experiments were performed in dependence on electrolyte concentration and for two different particle number densities. The time-dependent nucleation rate shows a pronounced initial peak, while postsolidification crystal size distributions are skewed towards larger crystallite sizes. At each concentration, the nucleation rate density initially drops exponentially with increasing salt concentration. The full data set, however, shows an unexpected scaling of the nucleation rate densities with metastability times the number density of particles. Parameterization of our results in terms of classical nucleation theory reveals unusually low interfacial free energies of the nucleus surfaces and nucleation barriers well below the thermal energy. We tentatively attribute our observations to the presence of doublets introduced by the employed conditioning technique. We discuss the conditions under which such small seeds may induce nucleation.

Fast crystal growth at ultra-low temperatures

It is believed that the slow liquid diffusion and geometric frustration brought by a rapid, deep quench inhibit fast crystallization and promote vitrification. Here we report fast crystal growth in charged colloidal systems under deep supercooling, where liquid diffusion is extremely low. By combining experiments and simulations, we show that this process occurs via wall-induced barrierless ordering consisting of two coupled steps: the step-like advancement of the rough interface that disintegrates frustration, followed by defect repairing inside the newly formed solid phase. The former is a diffusionless collective process, whereas the latter controls crystal quality. We further show that the intrinsic mechanical instability of a disordered glassy state subject to the crystal growth front allows for domino-like fast crystal growth even at ultra-low temperatures. These findings contribute to a deeper understanding of fast crystal growth and may be useful for applications related to vitrification prevention and crystal-quality control.

Thermodynamics of two-dimensional Yukawa systems across coupling regimes

Thermodynamics of two-dimensional Yukawa (screened Coulomb or Debye-Hückel) systems is studied systematically using molecular dynamics (MD) simulations. Simulations cover very broad parameter range spanning from weakly coupled gaseous states to strongly coupled fluid and crystalline states. Important thermodynamic quantities, such as internal energy and pressure, are obtained and accurate physically motivated fits are proposed. This allows us to put forward simple practical expressions to describe thermodynamic properties of two-dimensional Yukawa systems. For crystals, in addition to numerical simulations, the recently developed shortest-graph interpolation method is applied to describe pair correlations and hence thermodynamic properties. It is shown that the finite-temperature effects can be accounted for by using simple correction of peaks in the pair correlation function. The corresponding correction coefficients are evaluated using MD simulation. The relevance of the obtained results in the context of colloidal systems, complex (dusty) plasmas, and ions absorbed to interfaces in electrolytes is pointed out.

Effect of particle-size dynamics on properties of dense spongy-particle systems: Approach towards equilibrium

Open-porous deformable particles, often envisaged as sponges, are ubiquitous in biological and industrial systems (e.g., casein micelles in dairy products and microgels in cosmetics). The rich behavior of these suspensions is owing to the elasticity of the supporting network of the particle, and the viscosity of permeating solvent. Therefore, the rate-dependent size change of these particles depends on their structure, i.e., the permeability. This work aims at investigating the effect of the particle-size dynamics and the underlying particle structure, i.e., the particle permeability, on the transient and long-time behavior of suspensions of spongy particles in the absence of applied deformation, using the dynamic two-scale model developed by Hütter et al. [Farad. Discuss. 158, 407 (2012)1359-664010.1039/c2fd20025b]. In the high-density limit, the transient behavior is found to be accelerated by the particle-size dynamics, even at average size changes as small as 1%. The accelerated dynamics is evidenced by (i) the higher short-time diffusion coefficient as compared to elastic-particle systems and (ii) the accelerated formation of the stable fcc crystal structure. Furthermore, after long times, the particle-size dynamics of spongy particles is shown to result in lower stationary values of the energy and normal stresses as compared to elastic-particle systems. This dependence of the long-time behavior of these systems on the permeability, that essentially is a transport coefficient and hence must not affect the equilibrium properties, confirms that full equilibration has not been reached.

Formation of a transient amorphous solid in low density aqueous charged sphere suspensions

Colloidal glasses formed from hard spheres, nearly hard spheres, ellipsoids and platelets or their attractive variants, have been studied in great detail. Complementing and constraining theoretical approaches and simulations, the many different types of model systems have significantly advanced our understanding of the glass transition in general. Despite their early prediction, however, no experimental charged sphere glasses have been found at low density, where the competing process of crystallization prevails. We here report the formation of a transient amorphous solid formed from charged polymer spheres suspended in thoroughly deionized water at volume fractions of 0.0002-0.01. From optical experiments, we observe the presence of short-range order and an enhanced shear rigidity as compared to the stable polycrystalline solid of body centred cubic structure. On a density dependent time scale of hours to days, the amorphous solid transforms into this stable structure. We further present preliminary dynamic light scattering data showing the evolution of a second slow relaxation process possibly pointing to a dynamic heterogeneity known from other colloidal glasses and gels. We compare our findings to the predicted phase behaviour of charged sphere suspensions and discuss possible mechanisms for the formation of this peculiar type of colloidal glass.

Interpolation method for pair correlations in classical crystals

Effects of anharmonicity on the pair correlation function of classical crystals are studied. The recently proposed shortest-graph approach using the Gaussian representation of the individual correlation peaks (the peak width is determined by the length of the shortest graph connecting a given pair of particles) is further improved, to account for anharmonic corrections due to finite temperatures and hard-sphere-like interactions. Two major effects are identified, leading to a modification of the correlation peaks at large or short distances: (i) the peaks at large distances, well described by Gaussians, should be calculated from the finite-temperature phonon spectra; (ii) at short distances, the correlation peaks deviate significantly from the Gaussian form due to the lattice discreteness. We propose the analytical interpolation method, based on the shortest-graph approach, which includes both effects. By employing the molecular dynamics simulations, the accuracy of the method is verified for three- and two-dimensional crystals with the Yukawa, inverse-power-law, and pseudo-hard-sphere pair interactions. The capabilities of the method are demonstrated by calculating the phase diagram of a three-dimensional Yukawa system.

The relationship between local density and bond-orientational order during crystallization of the Gaussian core model

Whether nucleation is triggered by density or by bond-orientational order is one of the most hotly debated issues in recent investigations of the crystallization process. Here, we present a numerical study of the relationship between them for soft particles within the isothermal-isobaric ensemble. We compress the system and thus obtain the fluid-solid transition. By investigating locally dense-packed particles and particles with a relatively high bond-orientational order in the compressing process, we find a sharp increase of the spatial correlations for both densely packed particles and highly bond-orientational ordered particles at the phase transition point, which provide new characterization methods for the liquid-crystal transition. We also find that it is the bond-orientational order rather than density that triggers the nucleation process. The relationship between the local density and the bond-orientational order parameter is strongly affected by the characterization methods used. The local bond order parameter (q6) shows clear correlation with the local density (ρ) in the fluid stage, while the coarse-grained form (q[combining macron]6) does not correlate with ρ at all, owing to the comparable spatial scales of q6 and ρ. Nevertheless, q[combining macron]6 shows an obvious advantage in distinguishing between solid and liquid particles in our work. These results may elevate our understanding of the mechanism of the crystallization process.

Equilibrium fluid-crystal interfacial free energy of bcc-crystallizing aqueous suspensions of polydisperse charged spheres

The interfacial free energy is a central quantity in crystallization from the metastable melt. In suspensions of charged colloidal spheres, nucleation and growth kinetics can be accurately measured from optical experiments. In previous work, from these data effective nonequilibrium values for the interfacial free energy between the emerging bcc nuclei and the adjacent melt in dependence on the chemical potential difference between melt phase and crystal phase were derived using classical nucleation theory (CNT). A strictly linear increase of the interfacial free energy was observed as a function of increased metastability. Here, we further analyze these data for five aqueous suspensions of charged spheres and one binary mixture. We utilize a simple extrapolation scheme and interpret our findings in view of Turnbull's empirical rule. This enables us to present the first systematic experimental estimates for a reduced interfacial free energy, σ(0,bcc), between the bcc-crystal phase and the coexisting equilibrium fluid. Values obtained for σ(0,bcc) are on the order of a few k(B)T. Their values are not correlated to any of the electrostatic interaction parameters but rather show a systematic decrease with increasing size polydispersity and a lower value for the mixture as compared to the pure components. At the same time, σ(0) also shows an approximately linear correlation to the entropy of freezing. The equilibrium interfacial free energy of strictly monodisperse charged spheres may therefore be still greater.

Crystallization pathways of liquid-bcc transition for a model iron by fast quenching

We report simulations on the local structural evolution in the liquid-bcc transition of a model iron. Fourteen main Voronoi polyhedra are chosen as the representatives of short-range orders (SROs) and their transformations during crystallization are also investigated. Thus, the crystallization pathways for the main SROs are drawn. Our results also show that the transformations between two SROs in the crystallization pathways can be classified into two categories, first the enlargement of coordination number, second the transformation of local symmetry from five-fold to four-fold. The former reduces the potential energy while the latter increases it. It is found that the potential energy cannot decease monotonously whatever crystallization pathway is chosen to transform the icosahedral SRO to bcc SRO. Therefore, the latter transformation might provide the energy barrier of crystallization. We propose two transformation styles among SROs. All the transformations in the crystallization pathways can be achieved according to the styles. Moreover, the two transformation styles indicates that the bcc structure is more similar to liquid than other crystals. That might be the reason why the first phase nucleated during a rapid cooling process should be bcc crystal.

Crystallization Kinetics of Concurrent Liquid-Metastable and Metastable-Stable Transitions, and Ostwald's Step Rule

Experimental measurements of colloidal crystallization in a wide range of volume fractions of charged particles were performed to investigate the liquid-metastable-stable transition process. To fit the obtained experimental data, we developed a theoretical model to formulate the kinetics of the concurrent liquid-metastable and metastable-stable transitions. This model is well-supported by our observations. We found that when the ratio of the metastable-stable transition rate to the liquid-metastable rate is very large, the metastable state can become undetectable, although it still exists, offering a possible explanation for very few exceptions to Ostwald's step rule.