Homogeneous ice nucleation from supercooled water

Homogeneous ice nucleation from supercooled water was studied in the temperature range of 220-240 K through combining the forward flux sampling method (Allen et al., J. Chem. Phys., 2006, 124, 024102) with molecular dynamics simulations (FFS/MD), based on a recently developed coarse-grained water model (mW) (Molinero et al., J. Phys. Chem. B, 2009, 113, 4008). The calculated ice nucleation rates display a strong temperature dependence, ranging from 2.148 ± 0.635 × 10(25) m(-3) s(-1) at 220 K to 1.672 ± 0.970 × 10(-7) m(-3) s(-1) at 240 K. These rates can be fitted according to the classical nucleation theory, yielding an estimate of the effective ice-water interface energy γ(ls) of 31.01 ± 0.21 mJ m(-2) for the mW water model. Compared to experiments, our calculation underestimates the homogeneous ice nucleation rate by a few orders of magnitude. Possible reasons for the discrepancy are discussed. The nucleating ice embryo contains both cubic ice Ic and hexagonal ice Ih, with the fraction of each structure being roughly 50% when the critical size is reached. In particular, a novel defect structure containing nearly five-fold twin boundaries is identified in the ice clusters formed during nucleation. The way such defect structure is formed is found to be different from mechanisms proposed for the formation of the same defect in metallic nanoparticles and thin film. The quasi five-fold twin boundary structure found here is expected to occur in the crystallization of a wide range of materials with the diamond cubic structure, including ice.     

Two-dimensional water and ice layers: neutron diffraction studies at 278, 263, and 20 k

Neutron diffraction elucidates the structures of two-dimensional (2D) water layers (278 K) or 2D ice layers confined in an organic slit-shaped nanospace. The two-dimensional ice phases reported here consist of individual eight-membered rings or folded-chain segments (263 K) and condensed twelve-membered irregular rings (20 K). This is quite different from bulk or other 2D ice structures; the latter usually form hexagonal honeycomb lattices. Both low-temperature structures typically feature water molecules which are surrounded by two or three other water molecules. Neutron diffraction and thermochemical studies indicate a liquid-solid-phase transition around 277 K for two-dimensional D2O layers. A further solid-solid-phase transition occurs between 263 and 20 K.     

Dynamics of ice nucleation on water repellent surfaces

Prevention of ice accretion and adhesion on surfaces is relevant to many applications, leading to improved operation safety, increased energy efficiency, and cost reduction. Development of passive nonicing coatings is highly desirable, since current antiicing strategies are energy and cost intensive. Superhydrophobicity has been proposed as a lead passive nonicing strategy, yet the exact mechanism of delayed icing on these surfaces is not clearly understood. In this work, we present an in-depth analysis of ice formation dynamics upon water droplet impact on surfaces with different wettabilities. We experimentally demonstrate that ice nucleation under low-humidity conditions can be delayed through control of surface chemistry and texture. Combining infrared (IR) thermometry and high-speed photography, we observe that the reduction of water-surface contact area on superhydrophobic surfaces plays a dual role in delaying nucleation: first by reducing heat transfer and second by reducing the probability of heterogeneous nucleation at the water-substrate interface. This work also includes an analysis (based on classical nucleation theory) to estimate various homogeneous and heterogeneous nucleation rates in icing situations. The key finding is that ice nucleation delay on superhydrophobic surfaces is more prominent at moderate degrees of supercooling, while closer to the homogeneous nucleation temperature, bulk and air-water interface nucleation effects become equally important. The study presented here offers a comprehensive perspective on the efficacy of textured surfaces for nonicing applications.     

Heterogeneous ice nucleation in aqueous solutions: the role of water activity

Heterogeneous ice nucleation experiments have been performed with four different ice nuclei (IN), namely nonadecanol, silica, silver iodide and Arizona test dust. All IN are either immersed in the droplets or located at the droplets surface. The IN were exposed to various aqueous solutions, which consist of (NH4)2SO4, H2SO4, MgCl2, NaCl, LiCl, Ca(NO3)2, K2CO3, CH3COONa, ethylene glycol, glycerol, malonic acid, PEG300 or a NaCl/malonic acid mixture. Freezing was studied using a differential scanning calorimeter and a cold finger cell. The results show that the heterogeneous ice freezing temperatures decrease with increasing solute concentration; however, the magnitude of this effect is solute dependent. In contrast, when the results are analyzed in terms of the solution water activity a very consistent behavior emerges: heterogeneous ice nucleation temperatures for all four IN converge each onto a single line, irrespective of the nature of the solute. We find that a constant offset with respect to the ice melting point curve, Deltaaw,het, can describe the observed freezing temperatures for each IN. Such a behavior is well-known for homogeneous ice nucleation from supercooled liquid droplets and has led to the development of water-activity-based ice nucleation theory. The large variety of investigated solutes together with different general types of ice nuclei studied (monolayers, ionic crystals, covalently bound network-forming compounds, and a mixture of chemically different crystallites) underlines the general applicability of water-activity-based ice nucleation theory also for heterogeneous ice nucleation in the immersion mode. Finally, the ice nucleation efficiencies of the various IN, as well as the atmospheric implication of the developed parametrization are discussed.     

The nucleation rate of crystalline ice in amorphous solid water

The kinetics of crystalline ice nucleation and growth in nonporous, molecular beam deposited amorphous solid water (ASW) films are investigated at temperatures near 140 K. We implement an experimental methodology and corresponding model of crystallization kinetics to decouple growth from nucleation and quantify the temperature dependence and absolute rates of both processes. Nucleation rates are found to increase from approximately 3x10(13) m(-3) s(-1) at 134 K to approximately 2x10(17) m(-3) s(-1) at 142 K, corresponding to an Arrhenius activation energy of 168 kJ/mol. Over the same temperature range, the growth velocity increases from approximately 0.4 to approximately 4 A s(-1), also exhibiting Arrhenius behavior with an activation energy of 47 kJ/mol. These nucleation rates are up to ten orders of magnitude larger than in liquid water near 235 K, while growth velocities are approximately 10(9) times smaller. Crystalline ice nucleation kinetics determined in this study differ significantly from those reported previously for porous, background vapor deposited ASW, suggesting the nucleation mechanism is dependent upon film morphology.     

Three dimensional water and ice quantification using neutron imaging

Neutron radiography (NR) and computed tomography (NCT) are important non-destructive testing tools which determine information about an object’s interior structure and material properties. Although water quantification using 2-D neutron imaging has been applied to fuel cells for over a decade, the development of an accurate 3-D method has only recently been demonstrated. The 3-D water quantification technique developed at the Pennsylvania State University’s Radiation Science and Engineering Center has been applied to the quantification of both liquid and ice phase water. Quantification results of water and ice inside a known small channel test object were accurate to within 2%. This capability allows the quantification of ice within a fuel cell flow field under cold-start conditions.     

No evidence for large-scale proton ordering in Antarctic ice from powder neutron diffraction

We have examined a sample of 3000 year old Antarctic ice, collected at the Kohnen Station, by time-of-flight powder neutron diffraction to test the hypothesis of Fukazawa et al. [e.g., Ann. Glaciol. 31, 247 (2000)] that such ice may be partially proton ordered. Great care was taken to keep our sample below the proposed ordering temperature (237 K) at all times, but we did not observe any evidence of proton ordering.     

Molecular dynamics simulation of quasi-two-dimensional water clusters on ice nucleation protein

We performed molecular dynamics simulations of systems that consisted of the ice nucleation protein and the quasi-two-dimensional water cluster on it. The angle distributions, percolation probabilities, mean cluster sizes, cluster size distributions, and hydrogen bond relaxation times were analyzed. We concluded that the behavior of the water clusters on the ice nucleation protein was elaborately intertwined by the interaction between the ice nucleation protein and water, the interaction between the water molecules and the effect of temperature. The percolation probability and mean cluster size depended on the interactions and temperatures.     

Metadynamics simulations of ice nucleation and growth

The metadynamics method for accelerating rate events in molecular simulations is applied to the problem of ice freezing. We demonstrate homogeneous nucleation and growth of ice at 180 K in the isothermal-isobaric ensemble without the presence of external fields or surfaces. This result represents the first report of continuous and dynamic ice nucleation in a system of freely evolving density. Simulations are conducted using a variety of periodic simulation domains. In all cases the cubic polymorph ice I(c) is grown. The influence of boundary effects on estimates of the nucleation free energy barrier are discussed in relation to differences between this and earlier work.     

Neutron resonance capture and neutron diffraction analysis of Roman bronze water taps

Parts of four Roman bronze water taps, excavated in The Netherlands consisting of two cylinders and two cylinder housings, were non-destructively investigated by neutron resonance capture analysis and by time-of-flight neutron diffraction. The two complementary neutron methods have been, for the first time, combined for providing a comprehensive alloy characterization of intact, undisturbed archaeological bronze objects. Bulk chemical analyses and phase analyses of the copper alloys were carried out in order to establish relations and differences among the taps, and to address questions of standardization and fabrication.     

Free energy landscapes for homogeneous nucleation of ice for a monatomic water model

We simulate the homogeneous nucleation of ice from supercooled liquid water at 220 K in the isobaric-isothermal ensemble using the MW monatomic water potential. Monte Carlo simulations using umbrella sampling are performed in order to determine the nucleation free energy barrier. We find the Gibbs energy profile to be relatively consistent with that predicted by classical nucleation theory; the free energy barrier to nucleation was determined to be ~18 k(B)T and the critical nucleus comprised ~85 ice particles. Growth from the supercooled liquid gives clusters that are predominantly cubic, whilst starting with a pre-formed subcritical nucleus of cubic or hexagonal ice results in the growth of predominantly that phase of ice only.     

Adsorption structures of water in NaX studied by DRIFT spectroscopy and neutron powder diffraction

Diffuse reflectance infrared Fourier transform spectroscopic (DRIFTS) measurements (4000-1500 cm(-1)) and the results of neutron powder diffraction have been combined to study the structure of adsorption complexes of water in a NaX zeolite at different water loadings (25, 48, 72, and 120 water molecules per unit cell, respectively). Sharp bands corresponding to non-hydrogen-bonded OH groups of water molecules and broad associate bands due to hydrogen-bonded molecules are observed in the DRIFT spectra. We observe a remarkable downshift of the high-frequency associate band in a narrow temperature interval when the water amount decreases from 120 to 72 molecules per unit cell, which could signify some kind of "phase transition" for the water inside the zeolite cavities. Neutron powder diffraction results show that water molecules are predominantly localized in or near the 12-ring windows. Water molecules with hydrogen-bonded and non-hydrogen-bonded OH groups were found, in agreement with the observation of sharp and broad bands in the DRIFT spectra. We find strong evidence for the formation of cyclic hexamers of water molecules localized in the 12-ring windows, which are further stabilized by hydrogen bonds to framework oxygen atoms.     

Homogeneous ice nucleation from aqueous inorganic/organic particles representative of biomass burning: water activity, freezing temperatures, nucleation rates

Homogeneous ice nucleation plays an important role in the formation of cirrus clouds with subsequent effects on the global radiative budget. Here we report on homogeneous ice nucleation temperatures and corresponding nucleation rate coefficients of aqueous droplets serving as surrogates of biomass burning aerosol. Micrometer-sized (NH(4))(2)SO(4)/levoglucosan droplets with mass ratios of 10:1, 1:1, 1:5, and 1:10 and aqueous multicomponent organic droplets with and without (NH(4))(2)SO(4) under typical tropospheric temperatures and relative humidities are investigated experimentally using a droplet conditioning and ice nucleation apparatus coupled to an optical microscope with image analysis. Homogeneous freezing was determined as a function of temperature and water activity, a(w), which was set at droplet preparation conditions. The ice nucleation data indicate that minor addition of (NH(4))(2)SO(4) to the aqueous organic droplets renders the temperature dependency of water activity negligible in contrast to the case of aqueous organic solution droplets. The mean homogeneous ice nucleation rate coefficient derived from 8 different aqueous droplet compositions with average diameters of ～60 μm for temperatures as low as 195 K and a(w) of 0.82-1 is 2.18 × 10(6) cm(-3) s(-1). The experimentally derived freezing temperatures and homogeneous ice nucleation rate coefficients are in agreement with predictions of the water activity-based homogeneous ice nucleation theory when taking predictive uncertainties into account. However, the presented ice nucleation data indicate that the water activity-based homogeneous ice nucleation theory overpredicts the freezing temperatures by up to 3 K and corresponding ice nucleation rate coefficients by up to ～2 orders of magnitude. A shift of 0.01 in a(w), which is well within the uncertainty of typical field and laboratory relative humidity measurements, brings experimental and predicted freezing temperatures and homogeneous ice nucleation rate coefficients into agreement. The experimentally derived ice nucleation data are applied to constrain the water activity-based homogeneous ice nucleation theory to smaller than ±1 order of magnitude compared to the predictive uncertainty of larger than ±6 orders of magnitude. The atmospheric implications of these findings are discussed.     

Dissociation behavior of C2H6 hydrate at temperatures below the ice point: melting to liquid water followed by ice nucleation

The dissociation of C(2)H(6) hydrate particles by slow depressurization at temperatures slightly below the ice melting point was studied using optical microscopy and Raman spectroscopy. Visual observations and Raman measurements revealed that ethane hydrates can be present as a metastable state at pressures lower than the dissociation pressures of the three components: ice, hydrate, and free gas. However, they decompose into liquid water and gas phases once the system pressure drops to the equilibrium boundary for supercooled water, hydrate, and free gas. Structural analyses of obtained Raman spectra indicate that structures of the metastable hydrates and liquid water from the hydrate decay are fundamentally identical to those of the stable hydrates and supercooled water without experience of the hydration. These results imply a considerably high energy barrier for the direct hydrate-to-ice transition. Water solidification, probably induced by dynamic nucleation, was also observed during melting.     

Structural studies of water and other hydrogen—bonded liquids by neutron diffraction

The structural characteristics of hydrogen-bonded liquids may be studied by neutron diffraction. A brief review of the experimental techniques is given with particular emphasis on H/D isotope substitution and temperature difference techniques. These methods are illustrated with reference to recent work on water in its normal and supercooled liquid phases. The temperature variation is linked to the behaviour of the hydrogen bonds and the evolution towards the continuous random network of amorphous ice. A brief report is made for some other hydrogen-bonded liquids, contrasting the strong correlations in hydrogen fluoride which are not well understood with the study of methanol which gives good agreement with computer simulation results. Formic acid has also been investigated and novel techniques used in the analysis of the H/D datasets. The review ends with an overview of current issues and a consideration of future developments.     

Structural studies of water in hydrophilic and hydrophobic mesoporous silicas: an x-ray and neutron diffraction study at 297 K

Water confined in a sol-gel network has been characterized by x-ray and neutron diffraction for two samples of mesoporous silica: one with a hydrophilic character (a nonmodified one) and another with a hydrophobic character (a modified one with a methylated internal pore surface). The pore size has been previously characterized [J. Jelassi et al., Phys. Chem. Chem. Phys. 134, 1039 (2010)] to have a mean pore diameter of approximately 55 A?. The diffraction measurements presented in this paper have been made at room temperature [293 K] for a filling factor of 0.45, giving a mean thickness of 8-9 A? for the water layer. The results show that the local order of the confined water molecules in the intermediate region of 3-6 A? is significantly different from that of the bulk water and also for the two different environments. For the hydrophilic sample, the siloxyl groups at the surface modify the water structure through the effects of interfacial hydrogen-bonding, which influences the orientational configuration of local water molecules and creates a modified spatial arrangement in the pore. In the case of the hydrophobic sample, there is no specific interaction with the pore wall, which is primarily van der Waals type, and the water molecules at the interface are differently oriented to create a hydrogen-bonded network linked more directly to the rest of the water volume. In the present circumstances, the thickness of the water layer has a relatively small dimension so that the interpretation of the measured diffraction pattern is not as straightforward as for the bulk liquids, and it is necessary to consider the effects of diffraction-broadening from a distributed sample volume and also the contribution from cross-terms that remain after conducting a "wet-minus-dry" analysis procedure. These analytic difficulties are discussed in the context of the present measurements and compared with the work of other groups engaged in the study of water confined in different environments. The present results, again, emphasize the complexity influencing the properties of water in a confined geometry and the strong influence of surface interactions on its behavior.