Theoretical evidence for a dense fluid precursor to crystallization

Crystallization and arrest mechanisms of model colloids

We performed dynamic simulations of spheres with short-range attractive interactions for many values of interaction strength and range. Fast crystallization occurs in a localized region of this parameter space, but the character of crystallization pathways is not uniform within this region. Pathways range from one-step, in which a crystal nucleates directly from a gas, to two-step, in which substantial liquid-like clusters form and only subsequently become crystalline. Crystallization can fail because of slow nucleation from either gas or liquid, or because of dynamic arrest caused by strong interactions. Arrested states are characterized by the formation of networks of face-sharing tetrahedra that can be detected by a local common neighbor analysis.

Free-energy landscape and nucleation pathway of polymorphic minerals from solution in a Potts lattice-gas model

Metastable minerals commonly form during reactions between water and rock. The nucleation mechanism of polymorphic phases from solution are explored here using a two-dimensional Potts model. The model system is composed of a solvent and three polymorphic solid phases. The local state and position of the solid phase are updated by Metropolis dynamics. Below the critical temperature, a large cluster of the least stable solid phase initially forms in the solution before transitioning into more-stable phases following the Ostwald step rule. The free-energy landscape as a function of the modal abundance of each solid phase clearly reveals that before cluster formation, the least stable phase has an energetic advantage because of its low interfacial energy with the solution, and after cluster formation, phase transformation occurs along the valley of the free-energy landscape, which contains several minima for the regions of three phases. Our results indicate that the solid-solid and solid-liquid interfacial energy contribute to the formation of the complex free-energy landscape and nucleation pathways following the Ostwald step rule.

Optimization of crystal nucleation close to a metastable fluid-fluid phase transition

The presence of a metastable fluid-fluid critical point is thought to dramatically influence the crystallization pathway, increasing the nucleation rate by many orders of magnitude over the predictions of classical nucleation theory. We use molecular dynamics simulations to study the kinetics of crystallization in the vicinity of this metastable critical point and throughout the metastable fluid-fluid phase diagram. To quantitatively understand how the fluid-fluid phase separation affects the crystal nucleation, we evaluate accurately the kinetics and reconstruct the thermodynamic free-energy landscape of crystal formation. Contrary to expectations, we find no special advantage of the proximity of the metastable critical point on the crystallization rates. However, we find that the ultrafast formation of a dense liquid phase causes the crystallization to accelerate both near the metastable critical point and almost everywhere below the fluid-fluid spinodal line. These results unveil three different scenarios for crystallization that could guide the optimization of the process in experiments

A two-parameter extension of classical nucleation theory

A two-variable stochastic model for diffusion-limited nucleation is developed using a formalism derived from fluctuating hydrodynamics. The model is a direct generalization of the standard classical nucleation theory (CNT). The nucleation rate and pathway are calculated in the weak-noise approximation and are shown to be in good agreement with direct numerical simulations for the weak-solution/strong-solution transition in globular proteins. We find that CNT underestimates the time needed for the formation of a critical cluster by two orders of magnitude and that this discrepancy is due to the more complex dynamics of the two variable model and not, as often is assumed, a result of errors in the estimation of the free energy barrier.

On the question of two-step nucleation in protein crystallization

We report a real-time study on protein crystallization in the presence of multivalent salts using small angle X-ray scattering (SAXS) and optical microscopy, focusing particularly on the nucleation mechanism as well as on the role of the metastable intermediate phase (MIP). Using bovine beta-lactoglobulin as a model system in the presence of the divalent salt CdCl2, we have monitored the early stage of crystallization kinetics which demonstrates a two-step nucleation mechanism: protein aggregates form a MIP, which is followed by the nucleation of crystals within the MIP. Here we focus on characterizing and tuning the structure of the MIP using salt and the related effects on the two-step nucleation kinetics. The results suggest that increasing the salt concentration near the transition zone pseudo-c** enhances the energy barrier for both MIPs and crystal nucleation, leading to slow growth. The structural evolution of the MIP and its effect on subsequent nucleation is discussed based on the growth kinetics. The observed kinetics can be well described, using a rate-equation model based on a clear physical two-step picture. This real-time study not only provides evidence for a two-step nucleation process for protein crystallization, but also elucidates the role and the structural signature of the MIPs in the nonclassical process of protein crystallization.

Effect of mixing, concentration and temperature on the formation of mesostructured solutions and their role in the nucleation of dl-valine crystals

We report investigations on the formation of mesostructured solutions in dl-valine–water–2-propanol mixtures, and the crystallization of dl-valine from these solutions. Mesostructured liquid phases, similar to those previously observed in aqueous solutions of glycine and dl-alanine, were observed using Dynamic Light Scattering and Brownian microscopy, in both undersaturated and supersaturated solutions below a certain transition temperature. Careful experimentation was used to demonstrate that the optically clear mesostructured liquid phase, comprising colloidal mesoscale clusters dispersed within bulk solution, is thermodynamically stable and present in equilibrium with the solid phase at saturation conditions. Solutions prepared by slow cooling contained mesoscale clusters with a narrow size distribution and a mean hydrodynamic diameter of around 200 nm. Solutions of identical composition prepared by rapid isothermal mixing of valine aqueous solutions with 2-propanol contained mesoscale clusters which were significantly larger than those observed in slowly cooled solutions. The presence of larger mesoscale clusters was found to correspond to faster nucleation. Observed induction times were strongly dependent on the rapid initial mixing step, although solutions were left undisturbed afterwards and the induction times observed were up to two orders of magnitude longer than the initial mixing period. We propose that mesoscale clusters above a certain critical size are likely to be the location of productive nucleation events.

Crystal templating through liquid–liquid phase separation

Fluff-like crystal growth of water in 1,2-dichloroethylene.

Observing classical nucleation theory at work by monitoring phase transitions with molecular precision

It is widely accepted that many phase transitions do not follow nucleation pathways as envisaged by the classical nucleation theory. Many substances can traverse intermediate states before arriving at the stable phase. The apparent ubiquity of multi-step nucleation has made the inverse question relevant: does multistep nucleation always dominate single-step pathways? Here we provide an explicit example of the classical nucleation mechanism for a system known to exhibit the characteristics of multi-step nucleation. Molecular resolution atomic force microscopy imaging of the two-dimensional nucleation of the protein glucose isomerase demonstrates that the interior of subcritical clusters is in the same state as the crystalline bulk phase. Our data show that despite having all the characteristics typically associated with rich phase behaviour, glucose isomerase 2D crystals are formed classically. These observations illustrate the resurfacing importance of the classical nucleation theory by re-validating some of the key assumptions that have been recently questioned.

Many nanoscale systems can form ordered microphases through non-classical multistep nucleation. Here, the authors report that glucose isomerase, which is known to exhibit the characteristics of multi-step nucleation in 3D, nucleates along the pathway predicted by classical nucleation theory in 2D.

Theoretical calculation of the melting curve of Cu-Zr binary alloys

Helmholtz free energies of the dominant binary crystalline solids found in the Cu-Zr system at high temperatures close to the melting curve are calculated. Our theoretical approach combines fundamental measure density functional theory (applied to the hard-sphere reference system) and a perturbative approach to include the attractive interactions. The studied crystalline solids are Cu(fcc), Cu_{51}Zr_{14}(β), CuZr(B2), CuZr_{2}(C11b), Zr(hcp), and Zr(bcc). The calculated Helmholtz free energies of crystalline solids are in good agreement with results from molecular-dynamics (MD) simulations. Using the same perturbation approach, the liquid phase free energies are calculated as a function of composition and temperature, from which the melting curve of the entire composition range of this system can be obtained. Phase diagrams are determined in this way for two leading embedded atom method potentials, and the results are compared with experimental data. Theoretical melting temperatures are compared both with experimental values and with values obtained directly from MD simulations at several compositions.

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.

Solidification of a colloidal hard sphere like model system approaching and crossing the glass transition

We investigated the process of vitrification and crystallization in a model system of colloidal hard spheres. The kinetics of the solidification process was measured using time resolved static light scattering, while the time evolution of the dynamic properties was determined using time resolved dynamic light scattering. By performing further analysis we confirm that solidification of hard sphere colloids is mediated by precursors. Analyzing the dynamic properties we can show that the long time dynamics and thus the shear rigidity of the metastable melt is highly correlated with the number density of solid clusters (precursors) nucleated. In crystallization these objects convert into highly ordered crystals whereas in the case of vitrification this conversion is blocked and the system is (temporarily) locked in the metastable precursor state. From the early stages of solidification one cannot clearly conclude whether the melt will crystallize or vitrify. Furthermore our data suggests that colloidal hard sphere glasses can crystallize via homogeneous nucleation.

Self-assembly at a nonequilibrium critical point

We use analytic theory and computer simulation to study patterns formed during the growth of two-component assemblies in two and three dimensions. We show that these patterns undergo a nonequilibrium phase transition, at a particular growth rate, between mixed and demixed arrangements of component types. This finding suggests that principles of nonequilibrium statistical mechanics can be used to predict the outcome of multicomponent self-assembly, and suggests an experimental route to the self-assembly of multicomponent structures of a qualitatively defined nature.

Kinetic trapping of metastable amino acid polymorphs

Second harmonic generation (SHG) microscopy measurements indicate that inkjet-printed racemic solutions of amino acids can produce nanocrystals trapped in metastable polymorph forms upon rapid solvent evaporation. Polymorphism impacts the composition, distribution, and physico-kinetic properties of organic solids, with energetic arguments favoring the most stable polymorph. In this study, unfavored noncentrosymmetric crystal forms were observed by SHG microscopy. Polarization-dependent SHG measurement and synchrotron X-ray microdiffraction analysis of individual printed drops are consistent with formation of homochiral crystal production. Fundamentally, these results provide evidence supporting the ubiquity of Ostwald’s Rule of Stages, describing the hypothesized transitioning of crystals between metastable polymorphic forms in the early stages of crystal formation. Practically, the presence of homochiral metastable forms has implications on chiral resolution and on solid form preparations relying on rapid solvent evaporation.

Reentrant condensation, liquid–liquid phase separation and crystallization in protein solutions induced by multivalent metal ions

We briefly summarize the recent progress in tuning protein interactions as well as phase behavior in protein solutions using multivalent metal ions. We focus on the influence of control parameters and the mechanism of reentrant condensation, the metastable liquid–liquid phase separation and classical vs. non-classical pathways of protein crystallization.

Role of clusters in nonclassical nucleation and growth of protein crystals

The development of multistep nucleation theory has spurred on experimentalists to find intermediate metastable states that are relevant to the solidification pathway of the molecule under interest. A great deal of studies focused on characterizing the so-called "precritical clusters" that may arise in the precipitation process. However, in macromolecular systems, the role that these clusters might play in the nucleation process and in the second stage of the precipitation process, i.e., growth, remains to a great extent unknown. Therefore, using biological macromolecules as a model system, we have studied the mesoscopic intermediate, the solid end state, and the relationship that exists between them. We present experimental evidence that these clusters are liquid-like and stable with respect to the parent liquid and metastable compared with the emerging crystalline phase. The presence of these clusters in the bulk liquid is associated with a nonclassical mechanism of crystal growth and can trigger a self-purifying cascade of impurity-poisoned crystal surfaces. These observations demonstrate that there exists a nontrivial connection between the growth of the macroscopic crystalline phase and the mesoscopic intermediate which should not be ignored. On the other hand, our experimental data also show that clusters existing in protein solutions can significantly increase the nucleation rate and therefore play a relevant role in the nucleation process.

Diversity of transition pathways in the course of crystallization into ice VII

A new metastable ice phase is discovered in the freezing pathway into ice VII in accord with Ostwald's step rule.

Experimental modelling of single-particle dynamic processes in crystallization by controlled colloidal assembly

A comprehensive review of the experimental modeling of single particle dynamics in crystallization is presented.

How Crystals Nucleate and Grow in Aqueous NaCl Solution

Large-scale molecular dynamics simulations (64 000 particles) are used to examine the microscopic mechanism of crystal nucleation and growth in a slightly supersaturated solution of NaCl in water at 300 K and 1 atm. Early-stage nucleation is observed, and the growth of a single crystal is followed for ∼140 ns. It is shown that the nucleation and growth process is better described by Ostwald's rule of stages than by classical nucleation theory. Crystal nucleation originates in a region where the local salt concentration exceeds that of the bulk solution. The early-stage nucleus is a loosely ordered arrangement of ions that retains a significant amount of water. The residual water is slowly removed as the crystal grows and evolves toward its stable anhydrous state.

Effect of chain chirality on the self-assembly of sickle hemoglobin

We present simulation results on the self-assembly behavior of sickle hemoglobin (HbS). A coarse-grained HbS model, which contains hydrophilic and hydrophobic particles explicitly, is constructed to match the structural properties and physical description of HbS. The hydrophobic interactions are shown to be necessary with chirality being the main driver for the formation of HbS fibers. In the absence of chain chirality, only small self-assembled aggregates are observed whereas self-assembled elongated steplike bundle microstructures appear when we include chain chirality. We also investigate the effect of confinement on self-assembly, and find that elongated fibers-similar to open-space ones-can be obtained in hard confinement domains but cannot be formed within compliant red blood cell (RBC) domains under the same assumptions. We show, however, that by placing explicitly HbS fibers inside the RBCs and subjecting them to linear elongation and bending, we obtain different types of sickle-shaped RBCs as observed in sickle cell anemia.

Kinetically driven ordered phase formation in binary colloidal crystals

The aggregation of binary colloids of the same size and balanced charges is studied by Brownian dynamics simulations for dilute suspensions. It is shown that, under appropriate conditions, the formation of colloidal crystals is dominated by kinetic effects leading to the growth of well-ordered crystallites of the sodium-chloride (NaCl) bulk phase. These crystallites form with very high probability even when the cesium-chloride (CsCl) phase is more stable thermodynamically. Global optimization searches show that this result is not related to the most favorable structures of small clusters, which are either amorphous or of the CsCl structure. The formation of the NaCl phase is related to the specific kinetics of the crystallization process, which takes place by a two-step mechanism. In this mechanism, dense fluid aggregates form at first and then crystallization follows. It is shown that the type of short-range order in these dense fluid aggregates determines which phase is finally formed in the crystallites. The role of hydrodynamic effects in the aggregation process is analyzed by stochastic rotation dynamics - molecular dynamics simulations, and we find that these effects do not play a major role in the formation of the crystallites.

Self-assembly of multicomponent structures in and out of equilibrium

Theories of phase change and self-assembly often invoke the idea of a "quasiequilibrium," a regime in which the nonequilibrium association of building blocks results nonetheless in a structure whose properties are determined solely by an underlying free energy landscape. Here we study a prototypical example of multicomponent self-assembly, a one-dimensional fiber grown from red and blue blocks. We find that if the equilibrium structure possesses compositional correlations different from those characteristic of random mixing, then it cannot be generated without error at any finite growth rate: there is no quasiequilibrium regime. However, by exploiting dynamic scaling, structures characteristic of equilibrium at one point in phase space can be generated, without error, arbitrarily far from equilibrium. Our results, supported by mean-field theory in higher dimensions, thus suggest a "nonperturbative" strategy for multicomponent self-assembly in which the target structure is, by design, not the equilibrium one.

Heterogeneous nucleation of organic crystals mediated by single-molecule templates

Fundamental understanding of how crystals of organic molecules nucleate on a surface remains limited because of the difficulty of probing rare events at the molecular scale. Here we show that single-molecule templates on the surface of carbon nanohorns can nucleate the crystallization of two organic compounds from a supersaturated solution by mediating the formation of disordered and mobile molecular nanoclusters on the templates. Single-molecule real-time transmission electron microscopy indicates that each nanocluster consists of a maximum of approximately 15 molecules, that there are fewer nanoclusters than crystals in solution, and that in the absence of templates physisorption, but not crystal formation, occurs. Our findings suggest that template-induced heterogeneous nucleation mechanistically resembles two-step homogeneous nucleation.

Bond orientational order in liquids: Towards a unified description of water-like anomalies, liquid-liquid transition, glass transition, and crystallization: Bond orientational order in liquids

There are at least three fundamental states of matter, depending upon temperature and pressure: gas, liquid, and solid (crystal). These states are separated by first-order phase transitions between them. In both gas and liquid phases a complete translational and rotational symmetry exist, whereas in a solid phase both symmetries are broken. In intermediate phases between liquid and solid, which include liquid crystal and plastic crystal phases, only one of the two symmetries is preserved. Among the fundamental states of matter, the liquid state is the most poorly understood. We argue that it is crucial for a better understanding of liquids to recognize that a liquid generally has the tendency to have a local structural order and its presence is intrinsic and universal to any liquid. Such structural ordering is a consequence of many-body correlations, more specifically, bond angle correlations, which we believe are crucial for the description of the liquid state. We show that this physical picture may naturally explain difficult unsolved problems associated with the liquid state, such as anomalies of water-type liquids (water, Si, Ge, ...), liquid-liquid transition, liquid-glass transition, crystallization and quasicrystal formation, in a unified manner. In other words, we need a new order parameter representing a low local free-energy configuration, which is a bond orientational order parameter in many cases, in addition to a density order parameter for the physical description of these phenomena. Here we review our two-order-parameter model of liquid and consider how transient local structural ordering is linked to all of the above-mentioned phenomena. The relationship between these phenomena is also discussed.

The microscopic pathway to crystallization in supercooled liquids

Despite its fundamental and technological importance, a microscopic understanding of the crystallization process is still elusive. By computer simulations of the hard-sphere model we reveal the mechanism by which thermal fluctuations drive the transition from the supercooled liquid state to the crystal state. In particular we show that fluctuations in bond orientational order trigger the nucleation process, contrary to the common belief that the transition is initiated by density fluctuations. Moreover, the analysis of bond orientational fluctuations shows that these not only act as seeds of the nucleation process, but also i) determine the particular polymorph which is to be nucleated from them and ii) at high density favour the formation of fivefold structures which can frustrate the formation of crystals. These results can shed new light on our understanding of the relationship between crystallization and vitrification.

Stochastic resonance in multistable systems: the role of dimensionality

The theory of stochastic resonance in multistable systems is extended to account for both direct transitions between all stable states present and indirect ones involving intermediate states. It is shown that to satisfy these requirements the dynamics needs to be embedded in phase spaces of dimension equal to at least two. Under well defined conditions, the conjunction of the presence of intermediate states and the multidimensional character of the process leads to an enhancement of the response of the system to an external periodic forcing.

Amorphous nucleation precursor in highly nonequilibrium fluids

Dynamical density-functional simulations reveal structural aspects of crystal nucleation in undercooled liquids: The first appearing solid is amorphous, which promotes the nucleation of bcc crystals but suppresses the appearance of the fcc and hcp phases. These findings are associated with features of the effective interaction potential deduced from the amorphous structure.

Protein crystallization facilitated by molecularly imprinted polymers

We present a previously undescribed initiative and its application, namely the design of molecularly imprinted polymers (MIPs) for producing protein crystals that are essential for determining high-resolution 3D structures of proteins. MIPs, also referred to as "smart materials," are made to contain cavities capable of rebinding protein; thus the fingerprint of the protein created on the polymer allows it to serve as an ideal template for crystal formation. We have shown that six different MIPs induced crystallization of nine proteins, yielding crystals in conditions that do not give crystals otherwise. The incorporation of MIPs in screening experiments gave rise to crystalline hits in 8-10% of the trials for three target proteins. These hits would have been missed using other known nucleants. MIPs also facilitated the formation of large single crystals at metastable conditions for seven proteins. Moreover, the presence of MIPs has led to faster formation of crystals in all cases where crystals would appear eventually and to major improvement in diffraction in some cases. The MIPs were effective for their cognate proteins and also for other proteins, with size compatibility being a likely criterion for efficacy. Atomic force microscopy (AFM) measurements demonstrated specific affinity between the MIP cavities and a protein-functionalized AFM tip, corroborating our hypothesis that due to the recognition of proteins by the cavities, MIPs can act as nucleation-inducing substrates (nucleants) by harnessing the proteins themselves as templates.

Novel approach to controlled protein crystallization through ligandation of yttrium cations

Crystal structure determination of macromolecules is often hampered by the lack of crystals suitable for diffraction experiments. This article describes a protocol to crystallize the acidic protein bovine β-lactoglobulin in the presence of yttrium to yield high-quality crystals that belong to a new space group. The yttrium ions not only are used to engineer the crystallization, but are an integral part of the crystal lattice and can therefore be used to solve the phase problem using anomalous dispersion methods. Protein crystallization conditions were first optimized using an experimental phase diagram in the protein and salt concentration plane. Crystal growth strongly depends on the position in the phase diagram, and the best crystals grow near the phase transition boundaries. The structure analysis demonstrates the specific binding of yttrium ions to surface-exposed glutamate and aspartate side chains contributed by different molecules in the crystal lattice. By bridging molecules in this manner, contacts between molecules are formed that enable the formation of a stable crystal lattice. The potential application of this strategy to the crystallization of other acidic proteins is discussed on the basis of the universal features of the phase behavior of these proteins and the interactions induced by multivalent ions.