Homogeneous Ice Nucleation Rate in Water Droplets

Abstract

To predict the radiative forcing of clouds it is necessary to know the rate with which ice homogeneously nucleates in supercooled water. Such rate is often measured in drops to avoid the presence of impurities. At large supercooling small (nanoscopic) drops must be used to prevent simultaneous nucleation events. The pressure inside such drops is larger than the atmospheric one by virtue of the Laplace equation. In this work, we take into account such pressure raise in order to predict the nucleation rate in droplets using the TIP4P/Ice water model. We start from a recent estimate of the maximum drop size that can be used at each supercooling avoiding simultaneous nucleation events [Espinosa et al. J. Chem. Phys., 2016]. We then evaluate the pressure inside the drops with the Laplace equation. Finally, we obtain the rate as a function of the supercooling by interpolating our previous results for 1 and 2000 bar [Espinosa et al. Phys. Rev. Lett. 2016] using the Classical Nucleation Theory expression for the rate.This requires, in turn, interpolating the ice-water interfacial free energy and chemical potential difference. The TIP4P/Ice rate curve thus obtained is in good agreement with most droplet-based experiments. In particular, we find a good agreement with measurements performed using nanoscopic drops, that are currently under debate. The

Crystallization of hard spheres revisited. II. Thermodynamic modeling, nucleation work, and the surface of tension

“…parametrize the solid equation of state by shifting P s (ρ) so that P s (ρ m ) = P l (ρ f ). This yields a coexistence pressure P coex = 11.639(3) in good agreement with previous estimates [73,74]. We have then performed a thermodynamic integration to evaluate µ s , taking the coexistence as new reference state point where µ s (ρ m ) = µ l (ρ f ).…”

“…These simulations lead to the determination of the nucleation work ∆F c , which can then be used to quantify possible deviations from the classical picture on CNT. The second approach is based on the "seeding" of droplets [6,[37][38][39][40], whereby an initial droplet is placed into the metastable liquid at given supersaturation. This initial seed either grows or dissolves, with the crossing determining the size of the critical droplet n c .…”

Shear viscosity of pseudo hard-spheres

“…We limit ourselves to this temperature because according to Refs. [15,16] the equilibrium properties of the PHS and HS models are similar at this temperature. In what follows, all quantities are given in reduced units:…”

“…Using this nearly hard, but smooth potential, Jover et al have been able to reproduce structural and thermodynamic properties accurately compared with available simulation data for the original HS system. It has also been shown to provide reliable results for studying liquid-solid coexistence [16].…”

Effect of samarium doping on the nucleation of fcc-aluminum in undercooled liquids

“…Once the critical nucleus is obtained, we are able to compute the attachment rate at the critical nucleus. According to Auer and Frenkel [27,43,44], the attachment rate + can be measured in MD simulations as the diffusion coefficient for the size change at the critical point by (4)…”