Heterogeneous nucleation in the random field Ising model

We investigate the nucleation dynamics of the three-dimensional random field Ising model under an external field. We use umbrella sampling to compute the free-energy cost of a critical nucleus and use forward flux sampling for the direct estimation of nucleation rates. For moderate to strong disorder, our results indicate that the size of the nucleating cluster is not a good reaction coordinate, contrary to the pure Ising model. We rectify this problem by introducing a coordinate that also accounts for the location of the nucleus. Using the free energy barrier to predict the nucleation rate, we find reasonable agreement, although deviations become stronger as disorder increases. We attribute this effect to cluster shape fluctuations. We also discuss finite-size effects on the nucleation rate.

Anomalous Vapor and Ice Nucleation in Water at Negative Pressures: A Classical Density Functional Theory Study

In contrast to the abundance of work on the anomalous behavior of water, the relationship between the water's thermodynamic anomalies and kinetics of phase transition from metastable water is relatively unexplored. In this work, we have employed classical density functional theory to provide a unified and coherent picture of nucleation (both vapor and ice) from metastable water at negative pressure conditions. Our results suggest a peculiar nonmonotonic temperature dependence of vapor-liquid surface tension at temperatures where vapor-liquid coexistence is metastable with respect to the ice phase. The vapor nucleation barrier on isochoric cooling also shows a nonmonotonic temperature dependence. We further report that, for low density isochores, the temperature of the minimum vapor nucleation barrier (TΔΩv/min*) does not coincide with the temperature of maximum density (TMD) where metastability is maximum. The difference between the TΔΩv/min* and the TMD, however, decreases with increasing the density of the isochore. The vapor nucleation barrier along isobars shows an interesting crossover behavior in the vicinity of the Widom line on lowering the temperature. Our results on the ice nucleation suggest an anomalous retracing behavior of the nucleation barrier along isotherms at negative pressures and theoretically validate the recent findings that the reentrant ice(Ih)-liquid coexistence line can induce a drastic change in the kinetics of ice nucleation. Thus, this study establishes a direct connection between the metastable water's thermodynamic anomalies and the (vapor and ice) nucleation kinetics. In addition, this study provides deeper insights into the origin of the isothermal compressibility maximum on isochoric cooling.

Minimum in the pressure dependence of the interfacial free energy between ice Ih and water

Despite the importance of ice nucleation, this process has been barely explored at negative pressures. Here, we study homogeneous ice nucleation in stretched water by means of molecular dynamics seeding simulations using the TIP4P/Ice model. We observe that the critical nucleus size, interfacial free energy, free energy barrier, and nucleation rate barely change between isobars from -2600 to 500 bars when they are represented as a function of supercooling. This allows us to identify universal empirical expressions for homogeneous ice nucleation in the pressure range from -2600 to 500 bars. We show that this universal behavior arises from the pressure dependence of the interfacial free energy, which we compute by means of the mold integration technique, finding a shallow minimum around -2000 bars. Likewise, we show that the change in the interfacial free energy with pressure is proportional to the excess entropy and the slope of the melting line, exhibiting in the latter a reentrant behavior also at the same negative pressure. Finally, we estimate the excess internal energy and the excess entropy of the ice Ih-water interface.

Homogeneous nucleation rate of methane hydrate formation under experimental conditions from seeding simulations

In this work, we shall estimate via computer simulations the homogeneous nucleation rate for the methane hydrate at 400 bars for a supercooling of about 35 K. The TIP4P/ICE model and a Lennard-Jones center were used for water and methane, respectively. To estimate the nucleation rate, the seeding technique was employed. Clusters of the methane hydrate of different sizes were inserted into the aqueous phase of a two-phase gas-liquid equilibrium system at 260 K and 400 bars. Using these systems, we determined the size at which the cluster of the hydrate is critical (i.e., it has 50% probability of either growing or melting). Since nucleation rates estimated from the seeding technique are sensitive to the choice of the order parameter used to determine the size of the cluster of the solid, we considered several possibilities. We performed brute force simulations of an aqueous solution of methane in water in which the concentration of methane was several times higher than the equilibrium concentration (i.e., the solution was supersaturated). From brute force runs, we infer the value of the nucleation rate for this system rigorously. Subsequently, seeding runs were carried out for this system, and it was found that only two of the considered order parameters were able to reproduce the value of the nucleation rate obtained from brute force simulations. By using these two order parameters, we estimated the nucleation rate under experimental conditions (400 bars and 260 K) to be of the order of log₁₀ (J/(m³ s)) = -7(5).

Thermodynamics of diamond formation from hydrocarbon mixtures in planets

Hydrocarbon mixtures are extremely abundant in the Universe, and diamond formation from them can play a crucial role in shaping the interior structure and evolution of planets. With first-principles accuracy, we first estimate the melting line of diamond, and then reveal the nature of chemical bonding in hydrocarbons at extreme conditions. We finally establish the pressure-temperature phase boundary where it is thermodynamically possible for diamond to form from hydrocarbon mixtures with different atomic fractions of carbon. Notably, here we show a depletion zone at pressures above 200 GPa and temperatures below 3000 K-3500 K where diamond formation is thermodynamically favorable regardless of the carbon atomic fraction, due to a phase separation mechanism. The cooler condition of the interior of Neptune compared to Uranus means that the former is much more likely to contain the depletion zone. Our findings can help explain the dichotomy of the two ice giants manifested by the low luminosity of Uranus, and lead to a better understanding of (exo-)planetary formation and evolution.

Phase diagrams-Why they matter and how to predict them

Understanding the thermodynamic stability and metastability of materials can help us to, for example, gauge whether crystalline polymorphs in pharmaceutical formulations are likely to be durable. It can also help us to design experimental routes to novel phases with potentially interesting properties. In this Perspective, we provide an overview of how thermodynamic phase behavior can be quantified both in computer simulations and machine-learning approaches to determine phase diagrams, as well as combinations of the two. We review the basic workflow of free-energy computations for condensed phases, including some practical implementation advice, ranging from the Frenkel-Ladd approach to thermodynamic integration and to direct-coexistence simulations. We illustrate the applications of such methods on a range of systems from materials chemistry to biological phase separation. Finally, we outline some challenges, questions, and practical applications of phase-diagram determination which we believe are likely to be possible to address in the near future using such state-of-the-art free-energy calculations, which may provide fundamental insight into separation processes using multicomponent solvents.

Extending and validating bubble nucleation rate predictions in a Lennard-Jones fluid with enhanced sampling methods and transition state theory

We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach—based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction—allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Phase Transition between Crystalline Variants of Ordinary Ice

Water is one of the most abundant molecules on Earth. However, this common and "simple" material has more than 18 different phases, which poses a great challenge to theoretically study the nature of water and ice. We designed two reaction coordinates that can distinguish between water and various ice states and used them to efficiently sample all possible states of the system in all-atom molecular dynamics simulation at ambient temperature and pressure. Various structural and thermodynamics properties, including the water-ice phase diagrams, can thus be calculated. We also present a simple model that successfully explains the thermodynamic stability of different ice states. Our work provides effective methods and data for theoretical studies of different phases of water and ice.

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

Homogeneous ice nucleation rates for mW and TIP4P/ICE models through Lattice Mold calculations

Water freezing is the most common liquid-to-crystal phase transition on Earth, however, despite its critical implications on climate change and cryopreservation among other disciplines, its characterization through experimental and computational techniques remains elusive. In this work, we make use of computer simulations to measure the nucleation rate (J) of water at normal pressure under different supercooling conditions, ranging from 215 to 240K. We employ two different water models, mW, a coarse-grained potential for water, and TIP4P/ICE, an atomistic non-polarizable water model that provides one of the most accurate representations of the different ice phases. To evaluate J, we apply the Lattice Mold technique, a computational method based on the use of molds to induce the nucleus formation from the metastable liquid under conditions at which observing spontaneous nucleation would be unfeasible. With this method, we obtain estimates of the nucleation rate for ice Ih, Ic and a stacking mixture of ice Ih/Ic; reaching consensus with most of the previously reported rates, although differing with some others. Furthermore, we confirm that the predicted nucleation rates by the TIP4P/ICE model are in better agreement with experimental data than those obtained through the mW potential. Taken together, our study provides a reliable methodology to measure nucleation rates in a simple and computationally efficient manner which contributes to benchmarking the freezing behaviour of two popular water models.

Nucleation rates from small scale atomistic simulations and transition state theory

The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.

Is Ice Nucleation by Organic Crystals Nonclassical? An Assessment of the Monolayer Hypothesis of Ice Nucleation

Potent ice nucleating organic crystals display an increase in nucleation efficiency with pressure and memory effect after pressurization that set them apart from inorganic nucleants. These characteristics were proposed to arise from an ordered water monolayer at the organic-water interface. It was interpreted that ordering of the monolayer is the limiting step for ice nucleation on organic crystals, rendering their mechanism of nucleation nonclassical. Despite the importance of organics in atmospheric ice nucleation, that explanation has never been investigated. Here we elucidate the structure of interfacial water and its role in ice nucleation at ambient pressure on phloroglucinol dihydrate, the paradigmatic example of outstanding ice nucleating organic crystal, using molecular simulations. The simulations confirm the existence of an interfacial monolayer that orders on cooling and becomes fully ordered upon ice formation. The monolayer does not resemble any ice face but seamlessly connects the distinct hydrogen-bonding orders of ice and the organic surface. Although large ordered patches develop in the monolayer before ice nucleates, we find that the critical step is the formation of the ice crystallite, indicating that the mechanism is classical. We predict that the fully ordered, crystalline monolayer nucleates ice above -2 °C and could be responsible for the exceptional ice nucleation by the organic crystal at high pressures. The lifetime of the fully ordered monolayer around 0 °C, however, is too short to account for the memory effect reported in the experiments. The latter could arise from an increase in the melting temperature of ice confined by strongly ice-binding surfaces.

Classical nucleation theory predicts the shape of the nucleus in homogeneous solidification

Macroscopic models of nucleation provide powerful tools for understanding activated phase transition processes. These models do not provide atomistic insights and can thus sometimes lack material-specific descriptions. Here, we provide a comprehensive framework for constructing a continuum picture from an atomistic simulation of homogeneous nucleation. We use this framework to determine the equilibrium shape of the solid nucleus that forms inside bulk liquid for a Lennard-Jones potential. From this shape, we then extract the anisotropy of the solid-liquid interfacial free energy, by performing a reverse Wulff construction in the space of spherical harmonic expansions. We find that the shape of the nucleus is nearly spherical and that its anisotropy can be perfectly described using classical models.

Selective nucleation of ice crystals depending on the inclination angle of nanostructures

Heterogeneous nucleation is decided by many factors, and surface morphology is one of the most important elements. This paper reports the selective ice nucleation and growth process on a series of nanorods with different inclinations, which were rarely mentioned in previous research studies. It is found that the nanorods with special inclinations can cause the selective nucleation of ice crystals because of the spatial geometry matching. On this basis, we can regulate the ice crystal types (mainly including cubic ice and hexagonal ice) accordingly and even improve the freezing efficiency via controlling the inclinations of surface nanorods. In particular, cubic ice occupies the dominant role in the ice crystal on the surface of 45°-inclination nanorods, yet 90°-inclination nanorods are more beneficial for the formation of hexagonal ice. The shape of the nanorods not only controls the type of ice crystal, but also changes the freezing efficiency because different ice crystals have an unequal nucleation energy barrier. There are no apparent differences in the freezing efficiency on nanostructures with 45°, 75° and 90° inclination nanorods, and 60°-inclination nanorods are more favorable for ice nucleation. Our studies can promote the understanding on the selective nucleation of ice crystals and provide a theoretical basis for achieving the regulation of freezing efficiency.

Temperature Dependence of Homogeneous Nucleation in Ice

Ice nucleation is a process of great relevance in physics, chemistry, technology, and environmental sciences; much theoretical effort has been devoted to its understanding, but it still remains a topic of intense research. We shed light on this phenomenon by performing atomistic based simulations. Using metadynamics and a carefully designed set of collective variables, reversible transitions between water and ice are able to be simulated. We find that water freezes into a stacking disordered structure with the all-atom transferable intermolecular potential with 4 points/ice (TIP4P/ice) model, and the features of the critical nucleus of nucleation at the microscopic level are revealed. We have also estimated the ice nucleation rates along with other nucleation parameters at different undercoolings. Our results are in agreement with recent experimental and other theoretical works, and they confirm that nucleation is preceded by a large increase in tetrahedrally coordinated water molecules.

Ice Nucleation of Confined Monolayer Water Conforms to Classical Nucleation Theory

We confirmed that monolayer water confined by parallel graphene sheets spontaneously crystallizes from a structurally and dynamically heterogeneous liquid phase under moderate supercooling via direct molecular dynamics simulation. Square-lattice-like geometric order is observed at the early stage of nucleation and is preserved during the entire nucleus growth process. The diffusion coefficient and free energy profile in the cluster space extracted from a Bayesian trajectory analysis agree well with the classical nucleation theory (CNT) prediction and yield thermodynamic quantities exhibiting linear temperature dependence. The effectiveness of maximum cluster size as the descriptor of ice nucleation dynamics in the CNT framework can be attributed to the dynamical time scale decoupling and strong structural pattern dependence of density fluctuation in the liquid phase.

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.

Ab initio thermodynamics of liquid and solid water

Thermodynamic properties of liquid water as well as hexagonal (Ih) and cubic (Ic) ice are predicted based on density functional theory at the hybrid-functional level, rigorously taking into account quantum nuclear motion, anharmonic fluctuations, and proton disorder. This is made possible by combining advanced free-energy methods and state-of-the-art machine-learning techniques. The ab initio description leads to structural properties in excellent agreement with experiments and reliable estimates of the melting points of light and heavy water. We observe that nuclear-quantum effects contribute a crucial [Formula: see text] to the stability of ice Ih, making it more stable than ice Ic. Our computational approach is general and transferable, providing a comprehensive framework for quantitative predictions of ab initio thermodynamic properties using machine-learning potentials as an intermediate step.