Jorge R. Espinosa

Homogeneous ice nucleation evaluated for several water models

In this work, we evaluate by means of computer simulations the rate for ice homogeneous nucleation for several water models such as TIP4P, TIP4P/2005,TIP4P/ICE, and mW (following the same procedure as in Sanz et al. [J. Am. Chem. Soc. 135, 15008 (2013)]) in a broad temperature range. We estimate the ice-liquid interfacial free-energy, and conclude that for all water models γ decreases as the temperature decreases. Extrapolating our results to the melting temperature, we obtain a value of the interfacial free-energy between 25 and 32 mN/m in reasonable agreement with the reported experimental values. Moreover, we observe that the values of γ depend on the chosen water model and this is a key factor when numerically evaluating nucleation rates, given that the kinetic prefactor is quite similar for all water models with the exception of the mW (due to the absence of hydrogens). Somewhat surprisingly the estimates of the nucleation rates found in this work for TIP4P/2005 are slightly higher than those of the mW model, even though the former has explicit hydrogens. Our results suggest that it may be possible to observe in computer simulations spontaneous crystallization of TIP4P/2005 at about 60 K below the melting point.

Seeding approach to crystal nucleation

We present a study of homogeneous crystal nucleation from metastable fluids via the seeding technique for four different systems: mW water, Tosi-Fumi NaCl, Lennard-Jones, and Hard Spheres. Combining simulations of spherical crystal seeds embedded in the metastable fluid with classical nucleation theory, we are able to successfully describe the nucleation rate for all systems in a wide range of metastability. The crystal-fluid interfacial free energy extrapolated to coexistence conditions is also in good agreement with direct calculations of such parameter. Our results show that seeding is a powerful technique to investigate crystal nucleation.

On fluid-solid direct coexistence simulations: the pseudo-hard sphere model.

We investigate methodological issues concerning the direct coexistence method, an increasingly popular approach to evaluate the solid-fluid coexistence by means of computer simulations. The first issue is the impact of the simulation ensemble on the results. We compare the NpT ensemble (easy to use but approximate) with the NpzT ensemble (rigorous but more difficult to handle). Our work shows that both ensembles yield similar results for large systems (>5000 particles). Another issue, which is usually disregarded, is the stochastic character of a direct coexistence simulation. Here, we assess the impact of stochasticity in the determination of the coexistence point. We demonstrate that the error generated by stochasticity is much larger than that caused by the use of the NpT ensemble, and can be minimized by simply increasing the system size. To perform this study we use the pseudo hard-sphere model recently proposed by Jover et al. [J. Chem. Phys. 137, 144505 (2012)], and obtain a coexistence pressure of p∗ = 11.65(1), quite similar to that of hard spheres (only about 0.6% higher). Therefore, we conclude that this model can be reliably used to investigate the physics of hard spheres in phenomena like crystal nucleation.

Competition between ices Ih and Ic in homogeneous water freezing

The role of cubic ice, ice Ic, in the nucleation of ice from supercooled water has been widely debated in the past decade. Computer simulations can provide insightful information about the mechanism of ice nucleation at a molecular scale. In this work, we use molecular dynamics to study the competition between ice Ic and hexagonal ice, ice Ih, in the process of ice nucleation. Using a seeding approach, in which classical nucleation theory is combined with simulations of ice clusters embedded in supercooled water, we estimate the nucleation rate of ice for a pathway in which the critical nucleus has an Ic structure. Comparing our results with those previously obtained for ice Ih [Sanz et al., J. Am. Chem. Soc. 135, 15008 (2013)], we conclude that within the accuracy of our calculations both nucleation pathways have the same rate for the studied water models (TIP4P/Ice and TIP4P/2005). We examine in detail the factors that contribute to the nucleation rate and find that the chemical potential difference with the fluid, the attachment rate of particles to the cluster, and the ice-water interfacial free energy are the same within the estimated margin of error for both ice polymorphs. Furthermore, we study the morphology of the ice clusters and conclude that they have a spherical shape.

On the calculation of solubilities via direct coexistence simulations: Investigation of NaCl aqueous solutions and Lennard-Jones binary mixtures

Direct coexistence molecular dynamics simulations of NaCl solutions and Lennard-Jones binary mixtures were performed to explore the origin of reported discrepancies between solubilities obtained by direct interfacial simulations and values obtained from the chemical potentials of the crystal and solution phases. We find that the key cause of these discrepancies is the use of crystal slabs of insufficient width to eliminate finite-size effects. We observe that for NaCl crystal slabs thicker than 4 nm (in the direction perpendicular to the interface), the same solubility values are obtained from the direct coexistence and chemical potential routes, namely, 3.7 ± 0.2 molal at T = 298.15 K and p = 1 bar for the JC-SPC/E model. Such finite-size effects are absent in the Lennard-Jones system and are likely caused by surface dipoles present in the salt crystals. We confirmed that μs-long molecular dynamics runs are required to obtain reliable solubility values from direct coexistence calculations, provided that the initial solution conditions are near the equilibrium solubility values; even longer runs are needed for equilibration of significantly different concentrations. We do not observe any effects of the exposed crystal face on the solubility values or equilibration times. For both the NaCl and Lennard-Jones systems, the use of a spherical crystallite embedded in the solution leads to significantly higher apparent solubility values relative to the flat-interface direct coexistence calculations and the chemical potential values. Our results have broad implications for the determination of solubilities of molecular models of ionic systems.

The crystal-fluid interfacial free energy and nucleation rate of NaCl from different simulation methods

In this work, we calculate the crystal-fluid interfacial free energy, γ(cf), for the Tosi-Fumi model of NaCl using three different simulation techniques: seeding, umbrella sampling, and mold integration. The three techniques give an orientationaly averaged γ(cf) of about 100 mJ/m(2). Moreover, we observe that the shape of crystalline clusters embedded in the supercooled fluid is spherical. Using the mold integration technique, we compute γ(cf) for four different crystal orientations. The obtained interfacial free energies range from 100 to 114 mJ/m(2), being (100) and (111) the crystal planes with the lowest and highest γ(cf), respectively. Within the accuracy of our calculations, the interfacial free energy either does not depend on temperature or changes very smoothly with it. Combining the seeding technique with classical nucleation theory, we also estimate nucleation free energy barriers and nucleation rates for a wide temperature range (800-1040 K). The obtained results compare quite well with brute force calculations and with previous results obtained with umbrella sampling [Valeriani et al., J. Chem. Phys, 122, 194501 (2005)].

The mold integration method for the calculation of the crystal-fluid interfacial free energy from simulations

The interfacial free energy between a crystal and a fluid, γcf, is a highly relevant parameter in phenomena such as wetting or crystal nucleation and growth. Due to the difficulty of measuring γcf experimentally, computer simulations are often used to study the crystal-fluid interface. Here, we present a novel simulation methodology for the calculation of γcf. The methodology consists in using a mold composed of potential energy wells to induce the formation of a crystal slab in the fluid at coexistence conditions. This induction is done along a reversible pathway along which the free energy difference between the initial and the final states is obtained by means of thermodynamic integration. The structure of the mold is given by that of the crystal lattice planes, which allows to easily obtain the free energy for different crystal orientations. The method is validated by calculating γcf for previously studied systems, namely, the hard spheres and the Lennard-Jones systems. Our results for the latter show that the method is accurate enough to deal with the anisotropy of γcf with respect to the crystal orientation. We also calculate γcf for a recently proposed continuous version of the hard sphere potential and obtain the same γcf as for the pure hard sphere system. The method can be implemented both in Monte Carlo and Molecular Dynamics. In fact, we show that it can be easily used in combination with the popular Molecular Dynamics package GROMACS.

On the time required to freeze water

By using the seeding technique the nucleation rate for the formation of ice at room pressure will be estimated for the TIP4P/ICE model using longer runs and a smaller grid of temperatures than in the previous work. The growth rate of ice will be determined for TIP4P/ICE and for the mW model of water. Although TIP4P/ICE and mW have a similar melting point and melting enthalpy, they differ significantly in the dynamics of freezing. The nucleation rate of mW is lower than that of TIP4P/ICE due to its higher interfacial free energy. Experimental results for the nucleation rate of ice are between the predictions of these two models when obtained from the seeding technique, although closer to the predictions of TIP4P/ICE. The growth rate of ice for the mW model is four orders of magnitude larger than for TIP4P/ICE. Avramis expression is used to estimate the crystallization time from the values of the nucleation and growth rates. For mW the minimum in the crystallization time is found at approximately 85 K below the melting point and its value is of about a few ns, in agreement with the results obtained from brute force simulations by Moore and Molinero. For the TIP4P/ICE the minimum is found at about 55 K below the melting point, but its value is about ten microseconds. This value is compatible with the minimum cooling rate required to avoid the formation of ice and obtaining a glass phase. The crossover from the nucleation controlled crystallization to the growth controlled crystallization will be discussed for systems of finite size. This crossover could explain the apparent discrepancy between the values of J obtained by different experimental groups for temperatures below 230 K and should be considered as an alternative hypothesis to the two previously suggested: internal pressure and/or surface freezing effects. A maximum in the compressibility was found for the TIP4P/ICE model in supercooled water. The relaxation time is much smaller than the crystallization time at the temperature at which this maximum occurs, so this maximum is a real thermodynamic feature of the model. At the temperature of minimum crystallization time, the crystallization time is larger than the relaxation time by just two orders of magnitude.

NaCl nucleation from brine in seeded simulations: Sources of uncertainty in rate estimates

This work reexamines seeded simulation results for NaCl nucleation from a supersaturated aqueous solution at 298.15 K and 1 bar pressure. We present a linear regression approach for analyzing seeded simulation data that provides both nucleation rates and uncertainty estimates. Our results show that rates obtained from seeded simulations rely critically on a precise driving force for the model system. The driving force vs. solute concentration curve need not exactly reproduce that of the real system, but it should accurately describe the thermodynamic properties of the model system. We also show that rate estimates depend strongly on the nucleus size metric. We show that the rate estimates systematically increase as more stringent local order parameters are used to count members of a cluster and provide tentative suggestions for appropriate clustering criteria.

Role of Salt, Pressure, and Water Activity on Homogeneous Ice Nucleation

Pure water can be substantially supercooled below the melting temperature without transforming into ice. The achievable supercooling can be enhanced by adding solutes or by applying hydrostatic pressure. Avoiding ice formation is of great importance in the cryopreservation of food or biological samples. In this Letter, we investigate the similarity between the effects of pressure and salt on ice formation using a combination of state-of-the-art simulation techniques. We find that both hinder ice formation by increasing the energetic cost of creating the ice-fluid interface. Moreover, we examine the widely accepted proposal that the ice nucleation rate for different pressures and solute concentrations can be mapped through the activity of water [ Koop , L. ; Tsias , P. Nature , 2000 , 406 , 611 ]. We show that such a proposal is not consistent with the nucleation rates predicted in our simulations because it does not include all parameters affecting ice nucleation. Therefore, even though salt and pressure have a qualitatively similar effect on ice formation, they cannot be quantitatively mapped onto one another.