Precursor-mediated crystallization process in suspensions of hard spheres

Entropy-Driven Crystallization of Hard Colloidal Mixtures of Polymers and Monomers

The most trivial example of self-assembly is the entropy-driven crystallization of hard spheres. Past works have established the similarities and differences in the phase behavior of monomers and chains made of hard spheres. Inspired by the difference in the melting points of the pure components, we study, through Monte Carlo simulations, the phase behavior of athermal mixtures composed of fully flexible polymers and individual monomers of uniform size. We analyze how the relative number fraction and the packing density affect crystallization and the established ordered morphologies. As a first result, a more precise determination of the melting point for freely jointed chains of tangent hard spheres is extracted. A synergetic effect is observed in the crystallization leading to synchronous crystallization of the two species. Structural analysis of the resulting ordered morphologies shows perfect mixing and thus no phase separation. Due to the constraints imposed by chain connectivity, the local environment of the individual spheres, as quantified by the Voronoi polyhedron, is systematically more spherical and more symmetric compared to that of spheres belonging to chains. In turn, the local environment of the ordered phase is more symmetric and more spherical compared to that of the initial random packing, demonstrating the entropic origins of the phase transition. In general, increasing the polymer content reduces the degree of crystallinity and increases the melting point to higher volume fractions. According to the present findings, relative concentration is another determining factor in controlling the phase behavior of hard colloidal mixtures based on polymers.

Homogeneous crystallization in four-dimensional Lennard-Jones liquids

We observe homogeneous crystallization in simulated high-dimensional (d>3) liquids that follow physically realistic dynamics and have system sizes that are large enough to eliminate the possibility that crystallization was induced by the periodic boundary conditions. Supercooled four-dimensional (4D) Lennard-Jones (LJ) liquids maintained at zero pressure and constant temperatures 0.59 Collapse

Key Words Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Robert S Hoy Department of Physics, University of South Florida, Tampa, Florida 33620, USA

Collapse

Resolving Temperature-Dependent Hydrate Nucleation Pathway: The Role of "Transition Layer"

Understanding the nucleation of natural gas hydrate (NGH) at different conditions has important implications to NGH recovery and other industrial applications, such as gas storage and separation. Herein, vast numbers of hydrate nucleation events are traced via molecular dynamics (MD) simulations at different degrees of supercooling (or driving forces). Specifically, to precisely characterize a hydrate nucleus from an aqueous system during the MD simulation, we develop an evolutionary order parameter (OP) to recognize the nucleus size and shape. Subsequently, the free energy landscapes of hydrate during nucleation are explored by using the newly developed OP. The results suggest that at 270 K (or 0.92 Tm supercooling, where Tm is the melting point), the near-rounded nucleus prevails during the nucleation, as described from the classical nucleation theory. In contrast, at relatively strong driving forces of 0.85 and 0.88 Tm, nonclassical nucleation events arise. Specifically, the pathway toward an elongated nucleus becomes as important as the pathway toward a near-rounded nucleus. To explain the distinct nucleation phenomena at different supercoolings, a notion of a "transition layer" (or liquid-blob-like layer) is proposed. Here, the transition layer is to describe the interfacial region between the nucleus and aqueous solution, and this layer entails two functionalities: (1) it tends to retain CH4 depending on the degrees of supercooling and (2) it facilitates collision among CH4, which thus promote the incorporation of CH4 into nucleus. Our simulation indicates that compared to the near-rounded nucleus, the transition layer surrounding the elongated nucleus is more evident with the higher collision rate among CH4 molecules. As such, the transition layer tends to promote the elongated nucleus pathway, while offsetting the cost of larger surface free energy associated with the elongated nucleus. At 0.92 Tm, however, the transition layer gradually disappears, and classical nucleation events dominate. Overall, the notion of "transition layer" offers deeper insight into the NGH nucleation at different degrees of supercooling and could be extended to describe other types of hydrate nucleation.

Optothermal crystallization of hard spheres in an effective bidimensional geometry

Using colloids effectively confined in two dimensions by a cell with a thickness comparable to the particle size, we investigate the nucleation and growth of crystallites induced by locally heating the solvent with a near-infrared laser beam. The particles, which are "thermophilic," move towards the laser spot solely because of thermophoresis with no convection effects, forming dense clusters whose structure is monitored using two order parameters that gauge the local density and the orientational ordering. We find that ordering takes place when the cluster reaches an average surface density that is still below the upper equilibrium limit for the fluid phase of hard disks, meaning that we do not detect any sign of a proper "two-stage" nucleation from a glass or a polymorphic crystal structure. The crystal obtained at late growth stage displays a remarkable uniformity with a negligible amount of defects, arguably because the incoming particles diffuse, bounce, and displace other particles before settling at the crystal interface. This "fluidization" of the outer crystal edge may resemble the surface enhanced mobility giving rise to ultra-stable glasses by physical vapor deposition.

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.

Controlling crystallization: what liquid structure and dynamics reveal about crystal nucleation mechanisms

Over recent years, molecular simulations have provided invaluable insights into the microscopic processes governing the initial stages of crystal nucleation and growth. A key aspect that has been observed in many different systems is the formation of precursors in the supercooled liquid that precedes the emergence of crystalline nuclei. The structural and dynamical properties of these precursors determine to a large extent the nucleation probability as well as the formation of specific polymorphs. This novel microscopic view on nucleation mechanisms has further implications for our understanding of the nucleating ability and polymorph selectivity of nucleating agents, as these appear to be strongly linked to their ability in modifying structural and dynamical characteristics of the supercooled liquid, namely liquid heterogeneity. In this perspective, we highlight recent progress in exploring the connection between liquid heterogeneity and crystallization, including the effects of templates, and the potential impact for controlling crystallization processes. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

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.

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.

In Situ Kinetic Observations on Crystal Nucleation and Growth

Nucleation and growth are critical steps in crystallization, which plays an important role in determining crystal structure, size, morphology, and purity. Therefore, understanding the mechanisms of nucleation and growth is crucial to realize the controllable fabrication of crystalline products with desired and reproducible properties. Based on classical models, the initial crystal nucleus is formed by the spontaneous aggregation of ions, atoms, or molecules, and crystal growth is dependent on the monomer's diffusion and the surface reaction. Recently, numerous in situ investigations on crystallization dynamics have uncovered the existence of nonclassical mechanisms. This review provides a summary and highlights the in situ studies of crystal nucleation and growth, with a particular emphasis on the state-of-the-art research progress since the year 2016, and includes technological advances, atomic-scale observations, substrate- and temperature-dependent nucleation and growth, and the progress achieved in the various materials: metals, alloys, metallic compounds, colloids, and proteins. Finally, the forthcoming opportunities and challenges in this fascinating field are discussed.

Planar defect-free pure red perovskite light-emitting diodes via metastable phase crystallization

Solution-processable all-inorganic CsPbI_3-xBr_x perovskite holds great potential for pure red light-emitting diodes. However, the widely existing defects in this mixed halide perovskite markedly limit the efficiency and stability of present light-emitting diode devices. We here identify that intragrain Ruddlesden-Popper planar defects are primary forms of such defects in the CsPbI_3-xBr_x thin film owing to the lattice strain caused by inhomogeneous halogen ion distribution. To eliminate these defects, we develop a stepwise metastable phase crystallization strategy to minimize the CsPbI_3-xBr_x perovskite lattice strain, which brings planar defect-free CsPbI_3-xBr_x thin film with improved radiative recombination, narrowed emission band, and enhanced spectral stability. Using these high-quality thin films, we fabricate spectrally stable pure red perovskite light-emitting diodes, showing 17.8% external quantum efficiency and 9000 candela meter^-2 brightness with color coordinates required by Rec. 2020.

Template-Induced Precursor Formation in Heterogeneous Nucleation: Controlling Polymorph Selection and Nucleation Efficiency

We present an atomistic study of heterogeneous nucleation in Ni employing transition path sampling, which reveals a template precursor-mediated mechanism of crystallization. Most notably, we find that the ability of tiny templates to modify the structural features of the liquid and promote the formation of precursor regions with enhanced bond-orientational order is key to determining their nucleation efficiency and the polymorphs that crystallize. Our results reveal an intrinsic link between structural liquid heterogeneity and the nucleating ability of templates, which significantly advances our understanding toward the control of nucleation efficiency and polymorph selection.

Kinetics of Polymer Crystallization with Aggregating Small Crystallites

The isothermal crystallization near the glass transition temperature from the melt state of poly(trimethylene terephthalate) has been studied by wide-angle x-ray diffraction (WAXD), small-angle x-ray scattering (SAXS), and optical microscopy. The SAXS and WAXD results show the crystallization mechanism in which the crystalline nodules cover the entire sample with the formation of aggregation regions. The analysis of the SAXS results using Kolmogorov-Johnson-Mehl-Avrami theory indicates that the formation kinetics of the aggregation regions is of three-dimensional homogeneous nucleation type. The analysis of the SAXS profiles using Sekimoto's theory provides the growth velocity and the nucleation rate of the aggregation region. The temperature dependence of the growth velocity of the aggregation region is a natural extrapolation of that of spherulite to the high supercooling region. The temperature dependence of the nucleation rate of the aggregation region is also represented by the parameters of the spherulitic growth rate. The result of the growth velocities of the aggregation region and the spherulite suggests the existence of precursors at the front of the crystal growth.

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.

Crystallisation and polymorph selection in active Brownian particles

We explore crystallisation and polymorph selection in active Brownian particles with numerical simulation. In agreement with previous work (Wysocki et al. in Europhys Lett 105:48004, 2014), we find that crystallisation is suppressed by activity and occurs at higher densities with increasing Péclet number ([Formula: see text]). While the nucleation rate decreases with increasing activity, the crystal growth rate increases due to the accelerated dynamics in the melt. As a result of this competition, we observe the transition from a nucleation and growth regime at high [Formula: see text] to "spinodal nucleation" at low [Formula: see text]. Unlike the case of passive hard spheres, where preference for FCC over HCP polymorphs is weak, activity causes the annealing of HCP stacking faults, thus strongly favouring the FCC symmetry at high [Formula: see text]. When freezing occurs more slowly, in the nucleation and growth regime, this tendency is much reduced and we see a trend towards the passive case of little preference for either polymorph.

Fcc vs. hcp competition in colloidal hard-sphere nucleation: on their relative stability, interfacial free energy and nucleation rate

Hard-sphere crystallization has been widely investigated over the last six decades by means of colloidal suspensions and numerical methods. However, some aspects of its nucleation behaviour are still under debate. Here, we provide a detailed computational characterisation of the polymorphic nucleation competition between the face-centered cubic (fcc) and the hexagonal-close packed (hcp) hard-sphere crystal phases. By means of several state-of-the-art simulation techniques, we evaluate the melting pressure, chemical potential difference, interfacial free energy and nucleation rate of these two polymorphs, as well as of a random stacking mixture of both crystals. Our results highlight that, despite the fact that both polymorphs have very similar stability, the interfacial free energy of the hcp phase could be marginally higher than that of the fcc solid, which in consequence, mildly decreases its propensity to nucleate from the liquid compared to the fcc phase. Moreover, we analyse the abundance of each polymorph in grown crystals from different types of inserted nuclei: fcc, hcp and stacking disordered fcc/hcp seeds, as well as from those spontaneously emerged from brute force simulations. We find that post-critical crystals fundamentally grow maintaining the polymorphic structure of the critical nucleus, at least until moderately large sizes, since the only crystallographic orientation that allows stacking close-packed disorder is the fcc (111) plane, or equivalently the hcp (0001) one. Taken together, our results contribute with one more piece to the intricate puzzle of colloidal hard-sphere crystallization.

Temperature-dependent kinetic pathways of heterogeneous ice nucleation competing between classical and non-classical nucleation

Ice nucleation on the surface plays a vital role in diverse areas, ranging from physics and cryobiology to atmospheric science. Compared to ice nucleation in the bulk, the water-surface interactions present in heterogeneous ice nucleation complicate the nucleation process, making heterogeneous ice nucleation less comprehended, especially the relationship between the kinetics and the structures of the critical ice nucleus. Here we combine Markov State Models and transition path theory to elucidate the ensemble pathways of heterogeneous ice nucleation. Our Markov State Models reveal that the classical one-step and non-classical two-step nucleation pathways can surprisingly co-exist with comparable fluxes at T = 230 K. Interestingly, we find that the disordered mixing of rhombic and hexagonal ice leads to a favorable configurational entropy that stabilizes the critical nucleus, facilitating the non-classical pathway. In contrast, the favorable energetics promotes the formation of hexagonal ice, resulting in the classical pathway. Furthermore, we discover that, at elevated temperatures, the nucleation process prefers to proceed via the classical pathway, as opposed to the non-classical pathway, since the potential energy contributions override the configurational entropy compensation. This study provides insights into the mechanisms of heterogeneous ice nucleation and sheds light on the rational designs to control crystallization processes.

Impact of atomistic or crystallographic descriptors for classification of gold nanoparticles

Machine learning models are known to be sensitive to the features used to train them, but there is currently no way to predict the impact of using different features prior to feature extraction. This is particularly important to fields such as nanotechnology that are highly multi-disciplinary, and samples can be characterised many different ways depending on the preferences of individual researchers. Does it matter if nanomaterials are described using the interatomic coordinations or more complex order parameters? In this study we compare results of supervised and unsupervised learning on a single set of gold nanoparticles that has been characterised by two different descriptors, each with a unique feature space. We find that there are some consistencies, and model selection is descriptor-agnostic, but the level of detail and the type of information that can be extracted from the results is sensitive to the way the particles are described. Unsupervised clustering revealed that an atomistic descriptor provides a finer-grained interpretation and clusters that are sub-clusters of a more sophisticated crystallographic descriptor, which is consistent with both how the features were calculated, and how they are interpreted in the domain. A supervised classifier revealed that the types of features responsible for the separation are related to the bulk structure, regardless of the descriptor, but capture different types of information. For both the atomistic and crystallographic descriptor the gradient boosting decision tree classifier gave superior results of F1-scores of 0.96 and 0.98, respectively, with excellent precision and recall, even though the clustering presented a challenging multi-classification problem.

Entropy-Driven Heterogeneous Crystallization of Hard-Sphere Chains under Unidimensional Confinement

We investigate, through Monte Carlo simulations, the heterogeneous crystallization of linear chains of tangent hard spheres under confinement in one dimension. Confinement is realized through flat, impenetrable, and parallel walls. A wide range of systems is studied with respect to their average chain lengths (N = 12 to 100) and packing densities (ϕ = 0.50 to 0.61). The local structure is quantified through the Characteristic Crystallographic Element (CCE) norm descriptor. Here, we split the phenomenon into the bulk crystallization, far from the walls, and the projected surface crystallization in layers adjacent to the confining surfaces. Once a critical volume fraction is met, the chains show a phase transition, starting from regions near the hard walls. The established crystal morphologies consist of alternating hexagonal close-packed or face-centered cubic layers with a stacking direction perpendicular to the confining walls. Crystal layer perfection is observed with an increasing concentration. As in the case of the unconstrained phase transition of athermal polymers at high densities, crystal nucleation and growth compete with the formation of sites of a fivefold local symmetry. While surface crystallites show perfection with a predominantly triangular character, the morphologies of square crystals or of a mixed type are also formed. The simulation results show that the rate of perfection of the surface crystallization is not significantly faster than that of the bulk crystallization.

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.

Parasitic crystallization of colloidal electrolytes: growing a metastable crystal from the nucleus of a stable phase

Colloidal particles have been extensively used to comprehend the main principles governing liquid-crystal nucleation. Multiple mechanisms and frameworks have been proposed, through either experiments or computational approaches, to rationalise the ubiquitous formation of colloidal crystals. In this work, we elucidate the nucleation scenario behind the crystallization of oppositely charged colloids. By performing molecular dynamics simulations of colloidal electrolytes in combination with the Seeding technique, we evaluate the fundamental factors, such as the nucleation rate, free energy barrier, surface tension and kinetic pre-factor, that determine the liquid-to-solid transition of several crystalline polymorphs. Our results show that at a high packing fraction, there is a cross-over between the nucleation of the CsCl structure and that of a substitutionally disordered fcc phase, despite the CuAu crystal being the most stable phase. We demonstrate that the crucial factor in determining which phase nucleates the fastest is the free energy cost of the cluster formation rather than their kinetic ability to grow from the liquid. While at a low packing fraction, the stable phase, CsCl, is the one that nucleates and subsequently grows, we show how at moderate and high packing fractions, a disordered fcc phase subsequently grows regardless of the nature of the nucleating phase, termed parasitic crystallization. Taken together, our results provide a panoramic perspective of the complex nucleation scenario of oppositely charged colloids at moderate temperature and rationalise the different thermodynamic and kinetic aspects behind it.

Interplay of structural and dynamical heterogeneity in the nucleation mechanism in Ni

Gaining fundamental understanding of crystal nucleation processes in metal alloys is crucial for the development and design of high-performance materials with targeted properties. Yet, crystallizationis a complex non-equilibrium process and,...

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

Most systems have more than two stable crystalline states in the phase diagram, which is known as polymorphism. Crystallization in such a system is often under strong influence of competing orderings linked to those crystals. However, how such competition affects crystal nucleation and ordering toward the final crystalline state is largely unknown. This is primarily because the competition takes place locally and thus is masked by large positional fluctuations. We develop a unique method to correctly identify local symmetries by removing their distortions due to positional fluctuations. This allows us to experimentally access the spatiotemporal fluctuations of local symmetries at a single-particle level in crystallization of a charged colloidal system near the body-centered cubic-face-centered cubic border. Thus, we successfully reveal the crucial roles of competing ordering in the initial selection of polymorphs and the final grain boundary motion toward the most stable state from a microscopic perspective.

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.

Elucidating optical field directed hierarchical self-assembly of homogenous versus heterogeneous nanoclusters with femtosecond optical tweezers

Insights into the morphology of nanoclusters would facilitate the design of nano-devices with improved optical, electrical, and magnetic responses. We have utilized optical gradient forces for the directed self-assembly of colloidal clusters using high-repetition-rate femtosecond laser pulses to delineate their structure and dynamics. We have ratified our experiments with theoretical models derived from the Langevin equation and defined the valid ranges of applicability. Our femtosecond optical tweezer-based technique characterizes the in-situ formation of hierarchical self-assembled clusters of homomers as well as heteromers by analyzing the back focal plane displacement signal. This technique is able to efficiently distinguish between nano-particles in heterogeneous clusters and is in accordance with our theory. Herein, we report results from our technique, and also develop a model to describe the mechanism of such processes where corner frequency changes. We show how the corner frequency changes enables us to recognize the structure and dynamics of the coagulation of colloidal homogeneous and heterogeneous clusters in condensed media over a broad range of nanoparticle sizes. The methods described here are advantageous, as the backscatter position-sensitive detection probes the in-situ self-assembly process while other light scattering approaches are leveraged for the characterization of isolated clusters.

Can the pathway of stepwise nucleation be predicted and controlled?

Predicting the critical nucleus size and the nucleation barrier is of central importance in controlling the dynamics of nucleation. However, as the nucleation of a crystal involves intermediate states, the prediction becomes inaccessible with currently available models. Here, we show that based on single-particle level observations, the properties of crystal nuclei in a microscopic stepwise nucleation (MSN) can be well-quantified by incorporating the size and structure order parameter into the formula of free energy without prior knowledge of interfacial tension. The quantified free energy reveals that the intermediate structures arise from thermodynamics rather than kinetics. Precritical and postcritical nuclei are distinct not only in structure but also in the mechanism of crystalline ordering. The relative stability of intermediate structures and the pathway of nucleation can be well-controlled by supercooling. Our studies offer a successful approach to quantify MSN and shed new light on resolving the long-standing discrepancies between simulations and experiments.

Molecular dynamics investigation of the local structure in iron melts and its role in crystal nucleation during rapid solidification

Nucleation process of a bcc crystal after the formation of an MRO cluster.

Nucleation of pseudo hard-spheres and dumbbells at moderate metastability: appearance of A15 Frank–Kasper phase at intermediate elongations

Crystal nucleation of repulsive hard-dumbbells from the sphere to the two tangent spheres limit is investigated at moderately high metastability by brute-force molecular dynamics simulations.

Correlating anisotropy and disorder with the surface structure of platinum nanoparticles

Due to the competition between numerous physicochemical variables during formation and processing, platinum nanocatalysts typically contain a mixture of shapes, distributions of sizes, and a considerable degree of surface imperfection. Structural imperfection and sample polydispersivity are inevitable at scale, but accepting bulk and surface diversity as legitimate design features provides new opportunities for nanoparticle design. In recent years disorder and anisotropy have been embraced as useful design parameters but predicting the impact of uncontrollable imperfection a priori is challenging. In the present work we have created an ensemble of uniquely imperfect nanoparticles extracted from classical molecular dynamics trajectories and applied statistical filters to restrict the ensemble in ways that reflect common industrial design principles. We find that targeting different sizes and size distributions may be an effective way of promoting or suppressing internal disorder or crystallinity (as required), but the degree of anisotropy of the particle as a whole has a greater impact on the population of different types of surface ordering and active sites. These results indicate that tuning of disordered and anisotropic Pt nanoparticles is possible, but it is not as straightforward as geometrically ideal nanoparticles with a high degree of crystallinity. It is unlikely that a synthesis strategy could eliminate this diversity entirely, or ensure this type of structural complexity does not develop post-synthesis under operational conditions, but it may be possible to bias the formation of specific bulk structures and the surface anisotropy.

Reversible Self-Assembly (fcc-bct) Crystallization of Confined Granular Spheres via a Shear Dimensionality Mechanism

By combining vibrational annealing and shear dimensionality, we experimentally show (1) a fast reversible crystallization fcc-bct (face-centered cubic-body-centered tetragonal) in a granular system that is composed of dissipative millimeter-sized dry spheres, (2) a two-dimensional (planar) shear promotes self-assembly into an fcc crystal, while one-dimensional shear produces a bct crystal, and (3) in analogy with heterogeneous nucleation, a granular temperature gradient leads to the formation of crystal domains showing competition of polymorphic phases in the cold regions. Our findings suggest that controlling the directionality of the interactions steers to reversible crystallization of hard spheres, adds clues for theoretical studies, and provides a novel mechanism for the technological development of the applications of self-assembling phononic crystals.

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.

Crystallization of hard spheres revisited. I. Extracting kinetics and free energy landscape from forward flux sampling

We investigate the kinetics and the free energy landscape of the crystallization of hard spheres from a supersaturated metastable liquid though direct simulations and forward flux sampling. In this first paper, we describe and test two different ways to reconstruct the free energy barriers from the sampled steady state probability distribution of cluster sizes without sampling the equilibrium distribution. The first method is based on mean first passage times, and the second method is based on splitting probabilities. We verify both methods for a single particle moving in a double-well potential. For the nucleation of hard spheres, these methods allow us to probe a wide range of supersaturations and to reconstruct the kinetics and the free energy landscape from the same simulation. Results are consistent with the scaling predicted by classical nucleation theory although a quantitative fit requires a rather large effective interfacial tension.

Structural covariance in the hard sphere fluid

We study the joint variability of structural information in a hard sphere fluid biased to avoid crystallisation and form five-fold symmetric geometric motifs. We show that the structural covariance matrix approach, originally proposed for on-lattice liquids [P. Ronceray and P. Harrowell, J. Stat. Mech.: Theory Exp. 2016(8), 084002], can be meaningfully employed to understand structural relationships between different motifs and can predict, within the linear-response regime, structural changes related to motifs distinct from that used to bias the system.

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.

Nonclassical nucleation pathways in protein crystallization

Classical nucleation theory (CNT), which was established about 90 years ago, has been very successful in many research fields, and continues to be the most commonly used theory in describing the nucleation process. For a fluid-to-solid phase transition, CNT states that the solute molecules in a supersaturated solution reversibly form small clusters. Once the cluster size reaches a critical value, it becomes thermodynamically stable and favored for further growth. One of the most important assumptions of CNT is that the nucleation process is described by one reaction coordinate and all order parameters proceed simultaneously. Recent studies in experiments, computer simulations and theory have revealed nonclassical features in the early stage of nucleation. In particular, the decoupling of order parameters involved during a fluid-to-solid transition leads to the so-called two-step nucleation mechanism, in which a metastable intermediate phase (MIP) exists between the initial supersaturated solution and the final crystals. Depending on the exact free energy landscapes, the MIPs can be a high density liquid phase, mesoscopic clusters, or a pre-ordered state. In this review, we focus on the studies of nonclassical pathways in protein crystallization and discuss the applications of the various scenarios of two-step nucleation theory. In particular, we focus on protein solutions in the presence of multivalent salts, which serve as a model protein system to study the nucleation pathways. We wish to point out the unique features of proteins as model systems for further studies.

Triple Junction at the Triple Point Resolved on the Individual Particle Level

At the triple point of a repulsive screened Coulomb system, a fcc crystal, a bcc crystal, and a fluid phase coexist. At their intersection, these three phases form a liquid groove, the triple junction. Using confocal microscopy, we resolve the triple junction on a single-particle level in a model system of charged PMMA colloids in a nonpolar solvent. The groove is found to be extremely deep and the incommensurate solid-solid interface to be very broad. Thermal fluctuations hence appear to dominate the solid-solid interface. This indicates a very low interfacial energy. The fcc-bcc interfacial energy is quantitatively determined based on Young's equation and, indeed, it is only about 1.3 times higher than the fcc-fluid interfacial energy close to the triple point.