Nucleation kinetics in deionized charged colloidal model systems: a quantitative study by means of classical nucleation theory

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.

Which Surface Is More Scaling Resistant? A Closer Look at Nucleation Theories for Heterogeneous Gypsum Nucleation in Aqueous Solutions

Developing engineered surfaces with scaling resistance is an effective means to inhibit surface-mediated mineral scaling in various industries including desalination. However, contrasting results have been reported on the relationship between scaling potential and surface hydrophilicity. In this study, we combine a theoretical analysis with experimental investigation to clarify the effect of surface wetting property on heterogeneous gypsum (CaSO₄·2H₂O) formation on surfaces immersed in aqueous solutions. Theoretical prediction derived from classical nucleation theory (CNT) indicates that an increase of surface hydrophobicity reduces scaling potential, which contrasts our experimental results that more hydrophilic surfaces are less prone to gypsum scaling. We further consider the possibility of nonclassical pathway of gypsum nucleation, which proceeds by the aggregation of precursor clusters of CaSO₄. Accordingly, we investigate the affinity of CaSO₄ to substrate surfaces of varied wetting properties via calculating the total free energy of interaction, with the results perfectly predicting experimental observations of surface scaling propensity. This indicates that the interactions between precursor clusters of CaSO₄ and substrate surfaces might play an important role in regulating heterogeneous gypsum formation. Our findings provide evidence that CNT might not be applicable to describing gypsum scaling in aqueous solutions. The fundamental insights we reveal on gypsum scaling mechanisms have the potential to guide rational design of scaling-resistant engineered surfaces.

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.

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.

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.

Two-stage crystallization of charged colloids under low supersaturation conditions

We report simulations on the homogeneous liquid-fcc nucleation of charged colloids for both low and high contact energy values. As a precursor for crystal formation, we observe increased local order at the position where the crystal will form, but no correlations with the local density. Thus, the nucleation is driven by order fluctuations rather than density fluctuations. Our results also show that the transition involves two stages in both cases, first a transition of liquid → bcc, followed by a bcc → hcp/fcc transition. Both transitions have to overcome free energy barriers, so that a spherical bcc-like cluster is formed first, in which the final fcc structure is nucleated mainly at the surface of the crystallite. This means that the second stage bcc-fcc phase transition is a heterogeneous nucleation in the partially grown solid phase, even though we start from a homogeneous bulk liquid. The height of the bcc → hcp/fcc free energy barrier strongly depends on the contact energies of the colloids. For low contact energy this barrier is low, so that the bcc → hcp/fcc transition occurs spontaneously. For the higher contact energy, the second barrier is too high to be crossed spontaneously by the colloidal system. However, it was possible to ratchet the system over the second barrier and to transform the bcc nuclei into the stable hcp/fcc phase. The transitions are dominated by the first liquid-bcc transition and can be described by classical nucleation theory using an effective surface tension.

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.

Micro-structure evolution of wall based crystals after casting of model suspensions as obtained from Bragg microscopy

Growth of heterogeneously nucleated, wall based crystals plays a major role in determining the micro-structure during melt casting. This issue is here addressed using a model system of charged colloidal spheres in deionized aqueous suspension observed by Bragg microscopy which is a combination of light scattering and microscopy. We examine the evolution of the three-dimensional size, shape, and orientation of twin domains in monolithic crystals growing from two opposing planar walls into a meta-stable (shear-) melt. At each wall crystal orientation and twinning emerges during nucleation with small domains. During growth these widen and merge. From image analysis we observe the lateral coarsening velocities to follow a power law behaviour L(XY) ∝ t(1/2) as long as the vertical growth continues at constant speed. Lateral coarsening terminates upon intersection of the two solids and hardly any further ripening is seen. Initial lateral coarsening velocities show a Wilson Frenkel type dependence on the melt meta-stability.

Easy-use and low-cost fiber-based two-color dynamic light-scattering apparatus

We present a small and compact two-color cross-correlation light-scattering setup designed to study the structure and dynamics of colloidal suspensions in the regime of considerable turbidity. Using a homemade concentric four-arm goniometer, combined with fiber-optical illumination and detection devices, allows an easy and long time stable alignment-even under a temperature variation of ∼5 °C. Tests with a nearly multiple scattering free suspension of small Rayleigh scatterers show intercepts in cross correlation near 0.4 over a wide range of scattering angles from 20° to 150°. Measuring slightly turbid samples in cross-correlation mode multiple scattering is sufficiently suppressed allowing the determination of multiple scattering free dynamic structure factors as well as particle form factors.

Seed- and wall-induced heterogeneous nucleation in charged colloidal model systems under microgravity

Understanding the process that drives an undercooled fluid to the crystal state is still a challenging issue for condensed matter physics and plays a key role in designing new materials. The crystallization kinetics and the resulting polycrystalline morphology are given by a complex interplay of crystal nucleation, growth, and ripening. A great deal of progress has been made in recent years using colloidal suspensions as model systems in the study of crystallization. Close analogies to atomic systems are observed which can be exploited to address questions not accessible in atomic solidification. Here we present systematic measurements of the crystallization kinetics of a charged colloidal model system adding small amounts of seeds using time resolved scattering techniques. Large seeds show strong sedimentation under gravity even on the time scale of the crystallization process. To avoid this problem we performed our measurements under microgravity during parabolic flights. We report how the seed size and structure affect crystal nucleation and growth as functions of metastability giving the possibility to modify the crystallization process and the resulting microstructure of the polycrystal.

Colloids as model systems for metals and alloys: a case study of crystallization

Metallic systems are widely used as materials in daily human life. Their properties depend very much on the production route. In order to improve the production process and even develop novel materials a detailed knowledge of all physical processes involved in crystallization is mandatory. Atomic systems like metals are characterized by very high relaxation rates, which make direct investigations of crystallization very difficult and in some cases impossible. In contrast, phase transitions in colloidal systems are very sluggish and colloidal suspensions are optically transparent. Therefore, colloidal systems are often discussed as model systems for metals. In the present work, we study the process of crystallization of charged colloidal systems from the very beginning. Charged colloids offer the advantage that the interaction potential can be systematically tuned by a variation of the particle number density and the salt concentration. We use light scattering and ultra-small angle x-ray scattering to investigate the formation of short-range order in the liquid state even far from equilibrium, crystal nucleation and crystal growth. The results are compared with those of equivalent studies on metallic systems. They are critically assessed as regards similarities and differences.

Phase behaviour of deionized binary mixtures of charged colloidal spheres

We review recent work on the phase behaviour of binary charged sphere mixtures as a function of particle concentration and composition. Both size ratios Γ and charge ratios Λ are varied over a wide range. Unlike the case for hard spheres, the long-ranged Coulomb interaction stabilizes the crystal phase at low particle concentrations and shifts the occurrence of amorphous solids to particle concentrations considerably larger than the freezing concentration. Depending on Γ and Λ, we observe upper azeotrope, spindle, lower azeotrope and eutectic types of phase diagrams, all known well from metal systems. Most solids are of body centred cubic structure. Occasionally stoichiometric compounds are formed at large particle concentrations. For very low Γ, entropic effects dominate and induce a fluid-fluid phase separation. Since for charged spheres the charge ratio Λ is also decisive for the type of phase diagram, future experiments with charge variable silica spheres are suggested.

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.

Enhanced crystal stability in a binary mixture of charged colloidal spheres

We report on the phase behavior of a binary charged sphere mixture of size ratio Gamma=0.68 and charge ratio Lambda approximately 1 as a function of composition p and number density n . For p=0.1-0.3 we observe freezing at densities well below the freezing densities of the pure components. At all compositions our data indicate the formation of substitutional alloy crystals of body centered cubic structure. No indications for compound formation were observed. Rather, our findings point at the first observation of an upper azeotrope. Measurements of the crystallization kinetics reveal a combined density and composition dependence of growth velocities and nucleation rate densities, with small but significant anomalies at p*=0.2 . These correlated deviations can be rationalized within classical theories of solidification and suggest an increased similarity between melt and solid for this particular composition.

Surface tension of liquid ternary Fe-Cu-Mo alloys measured by electromagnetic levitation oscillating drop method

For the liquid Fe-Cu-Mo ternary system, the surface tensions of three selected alloys, i.e., three typical monotectic alloys, were measured by the electromagnetic levitation oscillating drop method over a broad temperature range, including both superheated and undercooled states. The maximum undercooling attained is up to 173 K. The experimental results show a good linear correlation between the surface tension and the temperature. By applying on the Butler equation, the surface tensions were also calculated and they are in good agreement with the measured ones, except that in the undercooled state, the calculated value is slightly larger than the measured results. Interestingly, both the measured and calculated results indicate that the enriched element on the droplet surface is much more conspicuous than other elements in influencing the surface tension. Besides, the viscosity and the density of the liquid Fe-Cu-Mo ternary alloys are also derived on the grounds of the experimentally measured surface tensions.