Revealing roles of competing local structural orderings in crystallization of polymorphic systems

In search of a precursor for crystal nucleation of hard and charged colloids

The interplay between crystal nucleation and the structure of the metastable fluid has been a topic of significant debate over recent years. In particular, it has been suggested that even in simple model systems such as hard or charged colloids, crystal nucleation might be foreshadowed by significant fluctuations in local structure around the location where the nucleus first arises. We investigate this using computer simulations of spontaneous nucleation events in both hard and charged colloidal systems. To detect local structural variations, we use both standard and unsupervised machine learning methods capable of finding hidden structures in the metastable fluid phase. We track numerous nucleation events for the face-centered cubic and body-centered cubic crystals on a local level and demonstrate that all signs of crystallinity emerge simultaneously from the very start of the nucleation process. We thus conclude that we observe no precursor for the crystal nucleation of hard and charged colloids.

Controlling polymorph selection during nucleation by tuning the structure of metallic melts

Controlling the polymorphism of crystals is crucial to the design of novel metallic materials with specific properties; however, the atomistic mechanism underlying polymorph selection during crystallization remains unclear. In this work, molecular dynamics simulations combined with well-tempered metadynamics simulations are employed to explore the atomic mechanisms of polymorph selection during the nucleation process of FCC aluminum and copper. Simulation results suggest that the distinct nucleation pathways of both FCC metals originate from different free-energy surfaces of nucleation processes and diverse symmetries of nucleation precursors. The initially forming phase from undercooled melts is most likely to be the one that has the symmetry closest to the precursors. Besides, tiny seeds with diverse crystal symmetries could induce the formation of preordered precursors for nucleation around the seed, leading to the reduction of free-energy barrier and thus the promotion of nucleation. Controlling polymorph selection with tiny seeds is realized by tuning the symmetry of precursors. Our findings not only shed significant light on understanding polymorph selection, but also provide theoretical guidance for better controlling the nucleation pathway in practice.

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.

Multiple Scenarios of Low-Temperature Nucleation in Simple Liquids

Usually, sufficient supercooling of a liquid is employed to bypass the free energy barrier and speed up crystallization. However, lowering the temperature T induces new issues competing with the crystallization, e.g., slow particle motion, geometric frustration, and the glass formation, which complicates our understanding of crystal growth. Here we systematically study the low-temperature nucleation kinetics discriminated by the maximum nucleation rate temperature T_{d} and the glass transition temperature T_{g}. At T_{d}, the ratio of the precursor and geometrically frustrated particles reaches the maximum. When T_{g}<T<T_{d}, nucleation kinetics is characterized by the subdiffusive slow particle motion, the high degrees of geometric frustration, and the saturation of precursors. In this regime, nucleation can proceed through the diffusionless-like ordering of precursors. Near T_{g}, there is a crossover regime, where geometrically frustrated particles percolate and the glass formation strongly slows down the nucleation. When T<T_{g}, diffusionless nucleation is obstructed due to the weak vibrational motion and the mechanical stability of the glassy state.

New Interpretation of Glass Formation in Isomeric Substances: Shifting from Melting-Point to Melting-Entropy

Revealing the critical thermodynamic parameters determining the glass formation of substances is of great significance for understanding the glass transition and guiding the composition design of glass-forming materials. Nevertheless, the direct access to glass-forming ability (GFA) by thermodynamics for various substances remains to be substantiated. The strategy to seek the fundamental properties of glass formation is explored several decades ago, as pioneered by Angell, arguing that the GFA in isomeric xylenes depends on the low lattice energy manifested by the low melting point. Here, an in-depth study is advanced using two more isomeric systems. Surprisingly, the results do not constantly support the reported relationship between the melting point and glass formation among isomeric molecules. Instead, molecules with enhanced glass formability are featured by the properties of low melting entropy without exception. Comprehensive studies of isomeric molecules find that the low melting entropy is roughly accompanied by the low melting point, explaining the apparent link between melting point and glass formation. Progressively, the viscosity measurements of the isomers uncover a strong dependence of the melting viscosity on melting entropy. These results emphasize the significance of the melting entropy in governing the glass formability of substances.

Revealing the role of liquid preordering in crystallisation of supercooled liquids

The recent discovery of non-classical crystal nucleation pathways has revealed the role of fluctuations in the liquid structural order, not considered in classical nucleation theory. On the other hand, classical crystal growth theory states that crystal growth is independent of interfacial energy, but this is questionable. Here we elucidate the role of liquid structural ordering in crystal nucleation and growth using computer simulations of supercooled liquids. We find that suppressing the crystal-like structural order in the supercooled liquid through a new order-killing strategy can reduce the crystallisation rate by several orders of magnitude. This indicates that crystal-like liquid preordering and the associated interfacial energy reduction play an essential role in nucleation and growth processes, forcing critical modifications of the classical crystal growth theory. Furthermore, we evaluate the importance of this additional factor for different types of liquids. These findings shed new light on the fundamental understanding of crystal growth kinetics.

An investigation of the kinetics and thermodynamics of NaCl nucleation through composite clusters

Having a good understanding of nucleation is critical for the control of many important processes, such as polymorph selection during crystallization. However, a complete picture of the molecular-level mechanisms of nucleation remains elusive. In this work, we take an in-depth look at the NaCl homogeneous nucleation mechanism through thermodynamics. Distinguished from the classical nucleation theory, we calculate the free energy of nucleation as a function of two nucleus size coordinates: crystalline and amorphous cluster sizes. The free energy surface reveals a thermodynamic preference for a nonclassical mechanism of nucleation through a composite cluster, where the crystalline nucleus is surrounded by an amorphous layer. The thickness of the amorphous layer increases with an increase in supersaturation. The computed free energy landscape agrees well with the composite cluster-free energy model, through which phase specific thermodynamic properties are evaluated. As the supersaturation increases, there is a change in stability of the amorphous phase relative to the solution phase, resulting in a change from one-step to two-step mechanism, seen clearly from the free energy profile along the minimum free energy path crossing the transition curve. By obtaining phase-specific diffusion coefficients, we construct the full mesoscopic model and present a clear roadmap for NaCl nucleation.

Revealing thermally-activated nucleation pathways of diffusionless solid-to-solid transition

Solid-to-solid transitions usually occur via athermal nucleation pathways on pre-existing defects due to immense strain energy. However, the extent to which athermal nucleation persists under low strain energy comparable to the interface energy, and whether thermally-activated nucleation is still possible are mostly unknown. To address these questions, the microscopic observation of the transformation dynamics is a prerequisite. Using a charged colloidal system that allows the triggering of an fcc-to-bcc transition while enabling in-situ single-particle-level observation, we experimentally find both athermal and thermally-activated pathways controlled by the softness of the parent crystal. In particular, we reveal three new transition pathways: ingrain homogeneous nucleation driven by spontaneous dislocation generation, heterogeneous nucleation assisted by premelting grain boundaries, and wall-assisted growth. Our findings reveal the physical principles behind the system-dependent pathway selection and shed light on the control of solid-to-solid transitions through the parent phase's softness and defect landscape.

An Artificial Neural Network Reveals the Nucleation Mechanism of a Binary Colloidal AB13 Crystal

Colloidal suspensions of two species have the ability to form binary crystals under certain conditions. The hunt for these functional materials and the countless investigations on their formation process are justified by the plethora of synergetic and collective properties these binary superlattices show. Among the many crystal structures observed over the past decades, the highly exotic colloidal icosahedral AB₁₃ crystal was predicted to be stable in binary hard-sphere mixtures nearly 30 years ago, yet the kinetic pathway of how homogeneous nucleation occurs in this system is still unknown. Here we investigate binary nucleation of the AB₁₃ crystal from a binary fluid phase of nearly hard spheres. We calculate the nucleation barrier and nucleation rate as a function of supersaturation and draw a comparison with nucleation of single-component and other binary crystals. To follow the nucleation process, we employ a neural network to identify the AB₁₃ phase from the binary fluid phase and the competing fcc crystal with single-particle resolution and significant accuracy in the case of bulk phases. We show that AB₁₃ crystal nucleation proceeds via a coassembly process where large spheres and icosahedral small-sphere clusters simultaneously attach to the nucleus. Our results lend strong support for a classical pathway that is well-described by classical nucleation theory, even though the binary fluid phase is highly structured and exhibits local regions of high bond orientational order.

Quantifying patterns in optical micrographs of one- and two-dimensional ellipsoidal particle assemblies

Current developments in colloidal science include the assembly of anisotropic colloids with broad geometric diversity. As the complexity of particle assemblies increases, the need for ubiquitous algorithms that quantitatively analyze images of the assemblies to deliver key information such as quantification of crystal structures becomes more urgent. This contribution describes algorithms capable of image analysis for classifying colloidal structures based on abstracted interparticle relationship information and quantitatively analyzing the abundance of each structure in mixed pattern assemblies. The algorithm parameters can be adjusted, allowing for the algorithms to be adapted for different image analyses. Three different ellipsoidal particle assembly images are presented to demonstrate the effectiveness of the algorithms: a one-dimensional (1D) particle chain assembly and two two-dimensional (2D) polymorphic crystals each consisting of assemblies of two distinct plane symmetry groups. Angle relationships between neighbouring particles are calculated and neighbour counts of each particle are determined. Combining these two parameters as rules for classification criteria allows for the labeling and quantification of each particle into a defined symmetry class within an assembly. The algorithms provide a labelled image comprising classification results and particle counts of each defined class. For multiple images or individual frames from a video, the script can be looped to achieve automatic processing. The yielded classification data allow for more in-depth image analysis of mixed pattern particle assemblies. We envision that these algorithms will have utility in quantitative analysis of images comprising ellipsoidal colloidal materials, nanoparticles, or biological matter.