Crystal nucleation in binary hard sphere mixtures: A Monte Carlo simulation study

Nanocrystal Assemblies: Current Advances and Open Problems

We explore the potential of nanocrystals (a term used equivalently to nanoparticles) as building blocks for nanomaterials, and the current advances and open challenges for fundamental science developments and applications. Nanocrystal assemblies are inherently multiscale, and the generation of revolutionary material properties requires a precise understanding of the relationship between structure and function, the former being determined by classical effects and the latter often by quantum effects. With an emphasis on theory and computation, we discuss challenges that hamper current assembly strategies and to what extent nanocrystal assemblies represent thermodynamic equilibrium or kinetically trapped metastable states. We also examine dynamic effects and optimization of assembly protocols. Finally, we discuss promising material functions and examples of their realization with nanocrystal assemblies.

Mechanism of Nucleation Pathway Selection in Binary Lennard-Jones Solution: A Combined Study of Molecular Dynamics Simulation and Free Energy Analysis

The nucleation process, which is the initial step in particle synthesis, determines the properties of the resultant particles. Although recent studies have observed various nucleation pathways, the physical factors that determine these pathways have not been fully elucidated. Herein, we conducted molecular dynamics simulations in a binary Lennard-Jones system as a model solution and found that the nucleation pathway can be classified into four types depending on microscopic interactions. The key parameters are (1) the strength of the solute-solute interaction and (2) the difference between the strengths of the like-pair and unlike-pair interactions. The increment of the former alters the nucleation mechanism from a two-step to a one-step pathway, whereas that of the latter causes quick assembly of solutes. Moreover, we developed a thermodynamic model based on the formation of core-shell nuclei to calculate the free energy landscapes. Our model successfully described the pathway observed in the simulations and demonstrated that the two parameters, (1) and (2), define the degree of supercooling and supersaturation, respectively. Thus, our model interpreted the microscopic insights from a macroscopic point of view. Because the only inputs required for our model are the interaction parameters, our model can a priori predict the nucleation pathway.

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.

Polymorphism and Perfection in Crystallization of Hard Sphere Polymers

We present results on polymorphism and perfection, as observed in the spontaneous crystallization of freely jointed polymers of hard spheres, obtained in an unprecedentedly long Monte Carlo (MC) simulation on a system of 54 chains of 1000 monomers. Starting from a purely amorphous configuration, after an initial dominance of the hexagonal closed packed (HCP) polymorph and a transitory random hexagonal close packed (rHCP) morphology, the system crystallizes in a final, stable, face centered cubic (FCC) crystal of very high perfection. An analysis of chain conformational characteristics, of the spatial distribution of monomers and of the volume accessible to them shows that the phase transition is caused by an increase in translational entropy that is larger than the loss of conformational entropy of the chains in the crystal, compared to the amorphous state. In spite of the significant local re-arrangements, as reflected in the bending and torsion angle distributions, the average chain size remains unaltered during crystallization. Polymers in the crystal adopt ideal random walk statistics as their great length renders local conformational details, imposed by the geometry of the FCC crystal, irrelevant.

Efficient event-driven simulations of hard spheres

 Hard spheres are arguably one of the most fundamental model systems in soft matter physics, and hence a common topic of simulation studies. Event-driven simulation methods provide an efficient method for studying the phase behavior and dynamics of hard spheres under a wide range of different conditions. Here, we examine the impact of several optimization strategies for speeding up event-driven molecular dynamics of hard spheres and present a light-weight simulation code that outperforms existing simulation codes over a large range of system sizes and packing fractions. The presented differences in simulation speed, typically a factor of five to ten, save significantly on both CPU time and energy consumption and may be a crucial factor for studying slow processes such as crystal nucleation and glassy dynamics.

Machine-Learned Free Energy Surfaces for Capillary Condensation and Evaporation in Mesopores

Using molecular simulations, we study the processes of capillary condensation and capillary evaporation in model mesopores. To determine the phase transition pathway, as well as the corresponding free energy profile, we carry out enhanced sampling molecular simulations using entropy as a reaction coordinate to map the onset of order during the condensation process and of disorder during the evaporation process. The structural analysis shows the role played by intermediate states, characterized by the onset of capillary liquid bridges and bubbles. We also analyze the dependence of the free energy barrier on the pore width. Furthermore, we propose a method to build a machine learning model for the prediction of the free energy surfaces underlying capillary phase transition processes in mesopores.

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.

The seven deadly sins: When computing crystal nucleation rates, the devil is in the details

The formation of crystals has proven to be one of the most challenging phase transformations to quantitatively model-let alone to actually understand-be it by means of the latest experimental technique or the full arsenal of enhanced sampling approaches at our disposal. One of the most crucial quantities involved with the crystallization process is the nucleation rate, a single elusive number that is supposed to quantify the average probability for a nucleus of critical size to occur within a certain volume and time span. A substantial amount of effort has been devoted to attempt a connection between the crystal nucleation rates computed by means of atomistic simulations and their experimentally measured counterparts. Sadly, this endeavor almost invariably fails to some extent, with the venerable classical nucleation theory typically blamed as the main culprit. Here, we review some of the recent advances in the field, focusing on a number of perhaps more subtle details that are sometimes overlooked when computing nucleation rates. We believe it is important for the community to be aware of the full impact of aspects, such as finite size effects and slow dynamics, that often introduce inconspicuous and yet non-negligible sources of uncertainty into our simulations. In fact, it is key to obtain robust and reproducible trends to be leveraged so as to shed new light on the kinetics of a process, that of crystal nucleation, which is involved into countless practical applications, from the formulation of pharmaceutical drugs to the manufacturing of nano-electronic devices.

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.

Self-assembly and entropic effects in pear-shaped colloid systems. II. Depletion attraction of pear-shaped particles in a hard-sphere solvent

We consider depletion effects of a pear-shaped colloidal particle in a hard-sphere solvent for two different model realizations of the pear-shaped colloidal particle. The two models are the pear hard Gaussian overlap (PHGO) particles and the hard pears of revolution (HPR). The motivation for this study is to provide a microscopic understanding for the substantially different mesoscopic self-assembly properties of these pear-shaped colloids, in dense suspensions, that have been reported in the previous studies. This is done by determining their differing depletion attractions via Monte Carlo simulations of PHGO and HPR particles in a pool of hard spheres and comparing them with excluded volume calculations of numerically obtained ideal configurations on the microscopic level. While the HPR model behaves as predicted by the analysis of excluded volumes, the PHGO model showcases a preference for splay between neighboring particles, which can be attributed to the special non-additive characteristics of the PHGO contact function. Lastly, we propose a potentially experimentally realizable pear-shaped particle model, the non-additive hard pear of revolution model, which is based on the HPR model but also features non-additive traits similar to those of PHGO particles to mimic their depletion behavior.

Spontaneous Crystallization in Systems of Binary Hard Sphere Colloids

Computer simulations of the fluid-to-solid phase transition in the hard sphere system were instrumental for our understanding of crystallization processes. But while colloid experiments and theory have been predicting the stability of several binary hard sphere crystals for many years, simulations were not successful to confirm this phenomenon. Here, we report the growth of binary hard sphere crystals isostructural to Laves phases, AlB_{2}, and NaZn_{13} in simulation directly from the fluid. We analyze particle kinetics during Laves phase growth using event-driven molecular dynamics simulations with and without swap moves that speed up diffusion. The crystallization process transitions from nucleation and growth to spinodal decomposition already deep within the fluid-solid coexistence regime. Finally, we present packing fraction-size ratio state diagrams in the vicinity of the stability regions of three binary crystals.

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.

Nucleation of Capillary Bridges and Bubbles in Nanoconfined CO2

Using molecular simulation, we examine the capillary condensation and the capillary evaporation of CO₂ in cylindrical nanopores. More specifically, we employ the recently developed μV T-S method to determine the microscopic mechanism associated with these processes and the corresponding free energy profiles. We calculate the free energy barrier for capillary condensation and identify that the key step consists in the nucleation of a liquid bridge of a critical size. Similarly, the free energy maximum found for the capillary evaporation process is found to correspond to the nucleation of a vapor bubble of a critical size. In addition, we assess the impact of the strength of the wall-fluid on the height of the free energy barrier and on the critical size of liquid bridges (condensation process) and vapor bubbles (evaporation process). We observe that the height of the free energy barrier increases with the strength of the wall-fluid interactions. Finally, we build a theoretical model, based on capillary theory, to rationalize our findings. In particular, the simulation results reveal a linear scaling of the free energy barrier with the critical size, in excellent agreement with the theoretical predictions.

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.

Correction: Polymer-enforced crystallization of a eutectic binary hard sphere mixture

Correction for 'Polymer-enforced crystallization of a eutectic binary hard sphere mixture' by Anna Kozina et al., Soft Matter, 2012, 8, 627-630.

Composition dependence of the glass forming ability in binary mixtures: The role of demixing entropy

We present a comparative study of the glass forming ability of binary systems with varying composition, where the systems have similar global crystalline structure (CsCl+fcc). Biased Monte Carlo simulations using umbrella sampling technique show that the free energy cost to create a CsCl nucleus increases as the composition of the smaller particles is decreased. We find that systems with comparatively lower free energy cost to form CsCl nucleus exhibit more pronounced pre-crystalline demixing near the liquid/crystal interface. The structural frustration between the CsCl and fcc crystal demands this demixing. We show that closer to the equimolar mixture, the entropic penalty for demixing is lower and a glass forming system may crystallize when seeded with a nucleus. This entropic penalty as a function of composition shows a non-monotonic behaviour with a maximum at a composition similar to the well known Kob-Anderson (KA) model. Although the KA model shows the maximum entropic penalty and thus maximum frustration against CsCl formation, it also shows a strong tendency towards crystallization into fcc lattice of the larger "A" particles which can be explained from the study of the energetics. Thus for systems closer to the equimolar mixture although it is the requirement of demixing which provides their stability against crystallization, for KA model it is not demixing but slow dynamics and the presence of the "B" particles make it a good glass former. The locally favoured structure around "B" particles is quite similar to the CsCl structure and the incompatibility of CsCl and fcc hinders the fcc structure growth in the KA model. Although the glass forming binary systems studied here are quite similar, differing only in composition, we find that their glass forming ability cannot be attributed to a single phenomenon.

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.

Enhancement of nucleation of protein crystals on nano-wrinkled surfaces

The synthesis of high quality protein crystals is essential for determining their structure. Hence the development of strategies to facilitate the nucleation of protein crystals is of prime importance. Recently, Ghatak and Ghatak [Langmuir 2013, 29, 4373] reported heterogeneous nucleation of protein crystals on nano-wrinkled surfaces. Through a series of experiments on different proteins, they were able to obtain high quality protein crystals even at low protein concentrations and sometimes without the addition of a precipitant. In this study, the mechanism of protein crystal nucleation on nano-wrinkled surfaces is studied through Monte Carlo simulations. The wrinkled surface is modeled by a sinusoidal surface. Free-energy barriers for heterogeneous crystal nucleation on flat and wrinkled surfaces are computed and compared. The study reveals that the enhancement of nucleation is closely related to the two step nucleation process seen during protein crystallization. There is an enhancement of protein concentration near the trough of the sinusoidal surface which aids in nucleation. However, the high curvature at the trough acts as a deterrent to crystal nucleus formation. Hence, significant lowering of the free-energy barrier is seen only if the increase in the protein concentration at the trough is very high.

Extension of the BMCSL equation of state for hard spheres to the metastable disordered region: Application to the SAFT approach

A simple modification of the Boublík-Mansoori-Carnahan-Starling-Leland equation of state is proposed for an application to the metastable disordered region. The new model has a positive pole at the jamming limit and can accurately describe the molecular simulation data of pure hard in the stable fluid region and along the metastable branch. The new model has also been applied to binary mixtures hard spheres, and an excellent description of the fluid and metastable branches can be obtained by adjusting the jamming packing fraction. The new model for hard sphere mixtures can be used as the repulsive term of equations of state for real fluids. In this case, the modified equations of state give very similar predictions of thermodynamic properties as the original models, and one can remove the multiple liquid density roots observed for some versions of the Statistical Associating Fluid Theory (SAFT) at low temperature without any modification of the dispersion term.

The role of bond tangency and bond gap in hard sphere crystallization of chains

We report results from Monte Carlo simulations on dense packings of linear, freely-jointed chains of hard spheres of uniform size. In contrast to our past studies where bonded spheres along the chain backbone were tangent, in the present work a finite tolerance in the bond is allowed. Bond lengths are allowed to fluctuate in the interval [σ, σ + dl], where σ is the sphere diameter. We find that bond tolerance affects the phase behaviour of hard-sphere chains, especially in the close vicinity of the melting transition. First, a critical dl(crit) exists marking the threshold for crystallization, whose value decreases with increasing volume fraction. Second, bond gaps enhance the onset of phase transition by accelerating crystal nucleation and growth. Finally, bond tolerance has an effect on crystal morphologies: in the tangent limit the majority of structures correspond to stack-faulted random hexagonal close packing (rhcp). However, as bond tolerance increases a wealth of diverse structures can be observed: from single fcc (or hcp) crystallites to random hcp/fcc stackings with multiple directions. By extending the simulations over trillions of MC steps (10(12)) we are able to observe crystal-crystal transitions and perfection even for entangled polymer chains in accordance to the Ostwald's rule of stages in crystal polymorphism. Through simple geometric arguments we explain how the presence of rigid or flexible constraints affects crystallization in general atomic and particulate systems. Based on the present results, it can be concluded that proper tuning of bond gaps and of the connectivity network can be a controlling factor for the phase behaviour of model, polymer-based colloidal and granular systems.

A design equation for low dosage additives that accelerate nucleation

Additives are used to control nucleation in many natural and industrial environments. However, the mechanisms by which additives inhibit or accelerate solute precipitate nucleation are not well understood. We propose an equation that predicts changes in nucleation barriers based on the adsorption properties and concentrations of trace additives. The equation shows that nucleant efficacy depends on the product of an adsorption equilibrium constant and the reduction in interfacial tension. Moreover, the two factors that determine the potency of additives are related to each other, suggesting that assays of just one property might facilitate additive design. We test the design equation for a Potts lattice gas model with surfactant-like additives in addition to solutes and solvents.

Crystallization kinetics of colloidal binary mixtures with depletion attraction

In this work the crystallization kinetics of colloidal binary mixtures with attractive interaction potential (Asakura-Oosawa) has been addressed. Parameters such as fraction of crystals, linear crystal dimension and crystal packing have been quantified in order to understand how the crystal formation is driven in terms of the depth of the attractive potential and the composition of the binary mixture (described by the number ratio). It was found that inside the eutectic triangle, crystallization is mainly governed by nucleation and the crystal packing is close to the close-packing of hard spheres. Moving out from the eutectic triangle towards small component results in the crystallization of small spheres. Enrichment of the eutectic mixture with large component results in the crystallization of both large and small spheres, however, the kinetics are completely different from those of the eutectic composition. Crosslinked polystyrene microgels with nearly hard sphere interactions were used as model systems. Attraction was introduced by addition of linear polystyrene. The time evolution of crystallization has been followed by static light scattering.

Reaction coordinate of incipient methane clathrate hydrate nucleation

Nucleation from solution is a ubiquitous phenomenon with relevance to myriad scientific disciplines, including pharmaceuticals, biomineralization, and disease. One prominent example is the nucleation of clathrate hydrates, multicomponent crystalline inclusion compounds relevant to the energy industry where they block pipelines and also constitute a potential vast energy resource. Despite their importance, the molecular mechanism of incipient hydrate formation remains unknown. Herein, we employ advanced molecular simulation tools (pB histogram, equilibrium path sampling) to provide a statistical-mechanical basis for extracting physical insight into the molecular steps by which clathrates form. Through testing the Mutually Coordinated Guest (MCG) order parameter, we demonstrate that both guest (methane) and host (water) structuring are crucial to accurately describe the nucleation of hydrates and determine a critical nucleus size of MCG-1 = 16 at 255 K and 500 bar. Equipped with a validated (and novel) reaction coordinate, subsequent equilibrium path sampling simulations yield the free energy barrier and nucleation rate. The resulting quantitative nucleation process is described by the MCG clustering mechanism. This constitutes a significant advance in the field of hydrates research, as the fitness of a molecular descriptor has never been statistically verified. More broadly, this work has significance to a wide range of multicomponent nucleation contexts wherein the formation mechanism depends on contributions from both solute and solvent.

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.

Numerical study of vapor condensation on patterned hydrophobic surfaces using the string method

Vapor condensation on solid surfaces plays a crucial role across a wide range of industrial applications. Recent advances of nanotechnology have made possible the manipulation of the condensation process through the control of surface structures. In this work, we study vapor condensation on hydrophobic surfaces patterned with microscale pillars. The critical nuclei, the activation barriers, and the minimum energy paths are computed using the climbing string method. The effects of pillar height, interpillar spacing, the level of supersaturation, and the intrinsic wettability of the solid surface on the nucleation process are investigated. Two nucleation scenarios are obtained from the computation. In the case of high pillar, narrow interpillar spacing, low supersaturation, and/or low surface wettability, the critical nucleus prefers the suspended Cassie state; otherwise, it prefers the impaled Wenzel state. A comparison of the nucleation barrier with that on a flat surface of the same material reveals that vapor condensation is inhibited by the microstructures in the former case, while enhanced in the latter case. The critical values of the pillar height, the interpillar spacing, and the supersaturation at which the critical nucleus changes from the Cassie state to the Wenzel state are identified from the phase diagram of the critical nucleus. It is found that the dependence of the critical interpillar spacing on the supersaturation follows closely the curve of the critical radii in a homogeneous nucleation. The relaxation dynamics of the condensate after the critical nucleus is formed is computed by solving the steepest descent equation. It is observed that when the pillar is low and/or the interpillar spacing is wide, a condensate initially in the Cassie state may evolve into the Wenzel state during the relaxation.

Condensation of droplets on nanopillared hydrophobic substrates

Using the constrained lattice density functional theory, we investigated the mechanism of droplet condensation, including droplet nucleation and growth, on nanopillared substrates. We find that similar to a macroscopic droplet on such a substrate, the critical nucleus also exhibits either the Wenzel or Cassie wetting state, depending on both the pillar height and the interpillar spacing. Our calculations show that there exists a critical value of the interpillar spacing, above which the critical nucleus is always in the Wenzel state and the pillared substrate always promotes the nucleation as compared to the smooth substrate, regardless of the pillar height. Below the critical interpillar spacing, however, the pillars always inhibit the nucleation, and the wetting state of the critical nucleus depends on the pillar height. Furthermore, our results demonstrate that the wetting state of the critical nuclei is not necessarily the wetting state of the formed microdroplets: droplets originated from the critical nuclei in the Wenzel state may change into the Cassie state spontaneously during the droplet growth process if the pillar height is greater than a critical value.

Calculation of the interfacial free energy of a binary hard-sphere fluid at a planar hard wall

Using molecular-dynamics simulation and Gibbs-Cahn Integration, we calculate the interfacial free energy γ of a binary hard-sphere fluid mixture at a structureless, planar hard wall. The calculation is performed as a function of packing fraction (density) for several values of the diameter ratio α = σ2∕σ1, where σ1 and σ2 are the diameters of the two particle types in the mixture. Our results are compared to those obtained from the bulk version of the White Bear Mark II (WBII) classical density-functional theory, which is a modification of the Fundamental-Measure Theory of Rosenfeld. The WBII bulk theory is shown to be in very good agreement with the simulation results, with significant deviation only at the very highest packing fractions.