Nanosecond freezing of water under multiple shock wave compression: optical transmission and imaging measurements

Water samples were subjected to multiple shock wave compressions, generating peak pressures of 1-5 GPa on nanosecond time scales. This loading process approximates isentropic compression and leads to temperatures where the ice VII phase is more stable than the liquid phase above 2 GPa. Time resolved optical transmission and imaging measurements were performed to determine the solidification rate under such conditions. Freezing occurred faster at higher pressures as water was compressed further into the ice VII phase, in agreement with classical micleation theory. Water consistently froze when in contact with a silica window, whereas no solidification occurred in the presence of sapphire windows. The transition was determined to be a surface initiated process--freezing began via heterogeneous nucleation at the water/window interface and propagated over thicknesses greater than 0.01 mm. The first optical images of freezing on nanosecond time scales were obtained. These images demonstrate heterogeneous nucleation and irregular solid growth over 0.01-0.10 mm lateral length scales and are consistent with latent heat emission during the transformation. The combination of optical transmission and imaging measurements presented here provide the first consistent evidence for freezing on short time scales.     

Nanosecond Rapid Crystallization of Water Induced by Quartz Glass under Dynamic Compression

Optical transmission characteristics of water between quartz glass under shock compression are in situ observed by using the technique of missile-borne light source. Through these transmission properties, the phase transition of liquid water is studied. The experimental results show that liquid water exhibits transparency decline phenomenon when the pressure is lower than 2 GPa under shock compression process, and the transparency variation is related to the existence of quartz glass. So, the transparency decline is attributed to a quartz-induced freezing phenomenon of water.     

Heat of freezing for supercooled water: measurements at atmospheric pressure

Unlike reversible phase transitions, the amount of heat released upon freezing of a metastable supercooled liquid depends on the degree of supercooling. Although terrestrial supercooled water is ubiquitous and has implications for cloud dynamics and nucleation, measurements of its heat of freezing are scarce. We have performed calorimetric measurements of the heat released by freezing water at atmospheric pressure as a function of supercooling. Our measurements show that the heat of freezing can be considerably below one predicted from a reversible hydrostatic process. Our measurements also indicate that the state of the resulting ice is not fully specified by the final pressure and temperature; the ice is likely to be strained on a variety of scales, implying a higher vapor pressure. This would reduce the vapor gradient between supercooled water and ice in mixed phase atmospheric clouds.     

Crystallization, melting, and structure of water nanoparticles at atmospherically relevant temperatures

Water nanoparticles play an important role in atmospheric processes, yet their equilibrium and nonequilibrium liquid-ice phase transitions and the structures they form on freezing are not yet fully elucidated. Here we use molecular dynamics simulations with the mW water model to investigate the nonequilibrium freezing and equilibrium melting of water nanoparticles with radii R between 1 and 4.7 nm and the structure of the ice formed by crystallization at temperatures between 150 and 200 K. The ice crystallized in the particles is a hybrid form of ice I with stacked layers of the cubic and hexagonal ice polymorphs in a ratio approximately 2:1. The ratio of cubic ice to hexagonal ice is insensitive to the radius of the water particle and is comparable to that found in simulations of bulk water around the same temperature. Heating frozen particles that contain multiple crystallites leads to Ostwald ripening and annealing of the ice structures, accompanied by an increase in the amount of ice at the expense of the liquid water, before the particles finally melt from the hybrid ice I to liquid, without a transition to hexagonal ice. The melting temperatures T(m) of the nanoparticles are not affected by the ratio of cubic to hexagonal layers in the crystal. T(m) of the ice particles decreases from 255 to 170 K with the particle size and is well described by the Gibbs-Thomson equation, T(m)(R) = T(m)(bulk) - K(GT)/(R - d), with constant K(GT) = 82 ± 5 K·nm and a premelted liquid of width d = 0.26 ± 0.05 nm, about one monolayer. The freezing temperatures also decrease with the particles' radii. These results are important for understanding the composition, freezing, and melting properties of ice and liquid water particles under atmospheric conditions.     

Initiation of the ice phase by marine biogenic surfaces in supersaturated gas and supercooled aqueous phases

Biogenic particles have the potential to affect the formation of ice crystals in the atmosphere with subsequent consequences for the hydrological cycle and climate. We present laboratory observations of heterogeneous ice nucleation in immersion and deposition modes under atmospherically relevant conditions initiated by Nannochloris atomus and Emiliania huxleyi, marine phytoplankton with structurally and chemically distinct cell walls. Temperatures at which freezing, melting, and water uptake occur are observed using optical microscopy. The intact and fragmented unarmoured cells of N. atomus in aqueous NaCl droplets enhance ice nucleation by 10-20 K over the homogeneous freezing limit and can be described by a modified water activity based ice nucleation approach. E. huxleyi cells covered by calcite plates do not enhance droplet freezing temperatures. Both species nucleate ice in the deposition mode at an ice saturation ratio, S(ice), as low as ~1.2 and below 240 K, however, for each, different nucleation modes occur at warmer temperatures. These observations show that markedly different biogenic surfaces have both comparable and contrasting effects on ice nucleation behaviour depending on the presence of the aqueous phase and the extent of supercooling and water vapour supersaturation. We derive heterogeneous ice nucleation rate coefficients, J(het), and cumulative ice nuclei spectra, K, for quantification and analysis using time-dependent and time-independent approaches, respectively. Contact angles, α, derived from J(het)via immersion freezing depend on T, a(w), and S(ice). For deposition freezing, α can be described as a function of S(ice) only. The different approaches yield different predictions of atmospheric ice crystal numbers primarily due to the time evolution allowed for the time-dependent approach with implications for the evolution of mixed-phase and ice clouds.     

Atomistic molecular dynamics simulations of shock compressed quartz

Atomistic non-equilibrium molecular dynamics simulations of shock wave compression of quartz have been performed using the so-called BKS semi-empirical potential of van Beest, Kramer, and van Santen [Phys. Rev. B 43, 5068 (1991)] to construct the Hugoniot of quartz. Our scheme mimics the real world experimental set up by using a flyer-plate impactor to initiate the shock wave and is the first shock wave simulation that uses a geometry optimised system of a polar slab in a three-dimensional system employing periodic boundary conditions. Our scheme also includes the relaxation of the surface dipole in the polar quartz slab which is an essential pre-requisite to a stable simulation. The original BKS potential is unsuited to shock wave calculations and so we propose a simple modification. With this modification, we find that our calculated Hugoniot is in good agreement with experimental shock wave data up to 25 GPa, but significantly diverges beyond this point. We conclude that our modified BKS potential is suitable for quartz under representative pressure conditions of the Earth core, but unsuitable for high-pressure shock wave simulations. We also find that the BKS potential incorrectly prefers the β-quartz phase over the α-quartz phase at zero-temperature, and that there is a β → α phase-transition at 6 GPa.     

Time-resolving analysis of cryotropic gelation of water/poly(vinyl alcohol) solutions via small-angle neutron scattering

The structural transformations occurring in initially homogeneous aqueous solutions of poly(vinyl alcohol) (PVA) through application of freezing (-13 degrees C) and thawing (20 degrees C) cycles is investigated by time resolving small-angle neutron scattering (SANS). These measurements indicate that formation of gels of complex hierarchical structure arises from occurrence of different elementary processes, involving different length and time scales. The fastest process that could be detected by our measurements during the first cryotropic treatment consists of the crystallization of the solvent. However, solvent crystallization is incomplete, and an unfrozen liquid microphase more concentrated in PVA than the initial solution is also formed. Crystallization of PVA takes place inside the unfrozen liquid microphase and is slowed down because of formation of a microgel fraction. Water crystallization takes place in the early 10 min of the treatment of the solution at subzero temperatures, and although below 0 degrees C the PVA solutions used for preparation of cryogels should be below the spinodal curve, occurrence of liquid-liquid phase separation could not be detected in our experiments. Upon thawing, ice crystals melt, and transparent gels are obtained that become opaque in approximately 200 min, due to a slow and progressive increase of the size of microheterogeneities (dilute and dense regions) imprinted during the fast freezing by the crystallization of water. During the permanence of these gels at room temperature (for hours), the presence of a high content of water (higher than 85% by mass) prevents further crystallization of PVA. Crystallization of PVA, in turn, is resumed by freezing the gels at subzero temperatures, after water crystallization and consequent formation of an unfrozen microphase. The kinetic parameters of PVA crystallization during the permanence of these gels at subzero temperatures are the same shown by PVA during the first freezing step of the solutions.     

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.     

Experimental investigation of the homogeneous freezing of aqueous ammonium sulfate droplets

We have measured the light scattering intensity and homogeneous ice nucleation temperatures from water droplets containing 0-33 wt % ammonium sulfate. In these laboratory experiments, we used a free-fall freezing tube technique to determine the fraction of frozen droplets at a particular droplet temperature by measuring the depolarized light scattering intensity from the droplets in free-fall. Previously reported freezing temperatures for solution concentrations greater than 5 wt % display a larger spread than can be accounted for by the reported experimental errors. We find freezing temperatures in good agreement with the lowest temperature freezing results reported by previous experiments. Our ammonium sulfate freezing temperature data set with water activity less than 0.98 is consistent with a curve that deviates in activity shift by about 5% from the best-fit ice nucleation temperature versus water activity curve found by Koop et al. in 2000, but the significance of this deviation will only be known with further high-precision ice nucleation temperature measurements for other aqueous solutions.     

Magnetic freezing of confined water

We report results from molecular dynamic simulations of the freezing transition of liquid water in the nanoscale hydrophobic confinement under the influence of a homogeneous external magnetic field of 10 T along the direction perpendicular to the parallel plates. A new phase of bilayer crystalline ice is obtained at an anomalously high freezing temperature of 340 K. The water-to-ice translation is found to be first order. The bilayer ice is built from alternating rows of hexagonal rings and rhombic rings parallel to the confining plates, with a large distortion of the hydrogen bonds. We also investigate the temperature shifts of the freezing transition due to the magnetic field. The freezing temperature, below which the freezing of confined water occurs, shifts to a higher value as the magnetic field enhances. Furthermore, the temperature of the freezing transition of confined water is proportional to the denary logarithm of the external magnetic field.     

Phase diagram of water between hydrophobic surfaces

Molecular dynamics simulations demonstrate that there are at least two classes of quasi-two-dimensional solid water into which liquid water confined between hydrophobic surfaces freezes spontaneously and whose hydrogen-bond networks are as fully connected as those of bulk ice. One of them is the monolayer ice and the other is the bilayer solid which takes either a crystalline or an amorphous form. Here we present the phase transformations among liquid, bilayer amorphous (or crystalline) ice, and monolayer ice phases at various thermodynamic conditions, then determine curves of melting, freezing, and solid-solid structural change on the isostress planes where temperature and intersurface distance are variable, and finally we propose a phase diagram of the confined water in the temperature-pressure-distance space.     

Formation of high density amorphous ice by decompression of ice VII and ice VIII at 135 K

Monte Carlo computer simulations of ice VII and ice VIII phases have been undertaken using the four-point transferable intermolecular potential model of water. By following thermodynamic paths similar to those used experimentally, ice is decompressed resulting in an amorphous phase. These phases are compared to the high density amorphous phase formed upon compression of ice Ih and are found to have very similar structures. By cooling liquid water along the water/Ih melting line a high density amorphous phase was also generated.     

Particle tracking microrheology of lyotropic liquid crystals

We present comprehensive results on the microrheological study of lyotropic liquid crystalline phases of various space groups constituted by water-monoglyceride (Dimodan) mixtures. In order to explore the viscoelastic properties of these systems, we use particle tracking of probe colloidal particles suitably dispersed in the liquid crystals and monitored by diffusing wave spectroscopy. The identification of the various liquid crystalline phases was separately carried out by small-angle X-ray scattering. The restricted motion of the particles was monitored and identified by the decay time of intensity autocorrelation function and the corresponding time-dependent mean square displacement (MSD), which revealed space group-dependent behavior. The characteristic time extracted by the intersection of the slopes of the MSD at short and long time scales, provided a characteristic time which could be directly compared with the relaxation time obtained by microrheology. Further direct comparison of microrheology and bulk rheology measurements was gained via the Laplace transform of the generalized time-dependent MSD, yielding the microrheology storage and loss moduli, G'(ω) and G'(ω), in the frequency domain ω. The general picture emerging from the microrheology data is that all liquid crystals exhibit viscoelastic properties in line with results from bulk rheology and the transition regime (elastic to viscous) differs according to the specific liquid crystal considered. In the case of the lamellar phase, a plastic fluid is measured by bulk rheology, while microrheology indicates viscoelastic behavior. Although we generally find good qualitative agreement between the two techniques, all liquid crystalline systems are found to relax faster when studied with microrheology. The most plausible explanation for this difference is due to the different length scales probed by the two techniques: that is, microscopical relaxation on these structured fluids, is likely to occur at shorter time scales which are more suitably probed by microrheology, whereas bulk, macroscopic relaxations occurring at longer time scales can only be probed by bulk rheology.     

Freezing of water and aqueous NaCl droplets coated by organic monolayers as a function of surfactant properties and water activity

This study presents heterogeneous ice nucleation from water and aqueous NaCl droplets coated by 1-nonadecanol and 1-nonadecanoic acid monolayers as a function of water activity (a(w)) from 0.8 to 1 accompanied by measurements of the corresponding pressure-area isotherms and equilibrium spreading pressures. For water and aqueous NaCl solutions of ~0-20 wt % in concentration, 1-nonadecanol exhibits a condensed phase, whereas the phase of 1-nonadecanoic acid changes from an expanded to a condensed state with increasing NaCl content of the aqueous subphase. 1-Nonadecanol-coated aqueous droplets exhibit the highest median freezing temperatures that can be described by a shift in a(w) of the ice melting curve by 0.098 according to the a(w)-based ice nucleation approach. This freezing curve represents a heterogeneous ice nucleation rate coefficient (J(het)) of 0.85 ± 0.30 cm(-2) s(-1). The median freezing temperatures of 1-nonadecanoic acid-coated aqueous droplets decrease less with increasing NaCl content compared to the homogeneous freezing temperatures. This trend in freezing temperature is best described by a linear function in a(w) and not by the a(w)-based ice nucleation approach most likely due to an increased ice nucleation efficiency of 1-nonadecanoic acid governed by the monolayer state. This freezing curve represents J(het) = 0.46 ± 0.16 cm(-2) s(-1). Contact angles (α) for 1-nonadecanol- and 1-nonadecanoic acid-coated aqueous droplets increase as temperature decreases for each droplet composition, but absolute values depend on employed water diffusivity and the interfacial energies of the ice embryo. A parametrization of log[J(het)(Δa(w))] is presented which allows prediction of freezing temperatures and heterogeneous ice nucleation rate coefficients for water and aqueous NaCl droplets coated by 1-nonadecanol without knowledge of the droplet's composition and α.     

Studies of cavitation and ice nucleation in 'doubly-metastable' water: time-lapse photography and neutron diffraction

High-speed photographic studies and neutron diffraction measurements have been made of water under tension in a Berthelot tube. Liquid water was cooled below the normal ice-nucleation temperature and was in a doubly-metastable state prior to a collapse of the liquid state. This transition was accompanied by an exothermic heat release corresponding with the rapid production of a solid phase nucleated by cavitation. Photographic techniques have been used to observe the phase transition over short time scales in which a solidification front is observed to propagate through the sample. Significantly, other images at a shorter time interval reveal the prior formation of cavitation bubbles at the beginning of the process. The ice-nucleation process is explained in terms of a mechanism involving hydrodynamically-induced changes in tension in supercooled water in the near vicinity of an expanding cavitation bubble. Previous explanations have attributed the nucleation of the solid phase to the production of high positive pressures. Corresponding results are presented which show the initial neutron diffraction pattern after ice-nucleation. The observed pattern does not exhibit the usual crystalline pattern of hexagonal ice [I(h)] that is formed under ambient conditions, but indicates the presence of other ice forms. The composite features can be attributed to a mixture of amorphous ice, ice-I(h)/I(c) and the high-pressure form, ice-III, and the diffraction pattern continues to evolve over a time period of about an hour.     

Solid/Liquid phase diagram of the ammonium sulfate/maleic acid/water system

We have studied the low temperature phase diagram and water activities of the ammonium sulfate/maleic acid/water system using differential scanning calorimetry and infrared spectroscopy of thin films. Using the results from our experiments, we have mapped the solid/liquid ternary phase diagram, determined the water activities based on the freezing point depression, and determined the ice/maleic acid phase boundary as well as the ternary eutectic composition and temperature. We also compare our results to the predictions of the extended AIM aerosol thermodynamics model and find good agreement for the ice melting points in the ice primary phase field of this system; however significant differences were found with respect to phase boundaries, maleic acid dissolution, and ammonium sulfate dissolution.     

Observation of water transport in the micro-porous layer of a polymer electrolyte fuel cell with a freezing method and cryo-scanning electron microscope

Micro-porous layers (MPLs) play an important role in the water management of polymer electrolyte fuel cells (PEFCs), however, the detailed mechanism of how the produced water is drained from these layers is not well understood. This paper observed the cross-sectional distribution of liquid water inside the cathode MPL to elucidate details of the phase state of the water transported through the MPL. The freezing method and cryo-scanning electron microscope (cryo-SEM) are used for the observations; the freezing method enables immobilization of the liquid water in the cell as ice forms by the freezing, and the cryo-SEM can visualize the water distribution in the vicinity of the MPL at high resolution without the ice melting. It was shown that no liquid water accumulates inside the MPL in operation at 35 °C, while the pores of the MPL are filled with liquid water under very low cell temperature operation, at 5 °C. These results indicate that the produced water passes through the MPL not as a liquid but in the vapor state in usual PEFC operation. Additionally, liquid water at the interface between the MPL and a catalyst layer (CL) was identified, and the effect of the interfacial contact on the water distribution was examined.     

Frequency response to liquid density of a piezoelectric quartz crystal sensor with longitudinal wave

Elimination of longitudinal wave effect is an important aspect in the detection cell design, although such consideration is ignored in most of references. Three detection cells were designed to investigate the influence of longitudinal wave effect on the frequency response of a piezoelectric quartz crystal (PQC) to liquid density. In the cell with horizontally mounted quartz crystal, the air/liquid interface acts the reflection surface for the longitudinal wave. The variation in liquid height by regent addition or solvent evaporation can result in fluctuation in the oscillating frequency of the PQC. The influence of the longitudinal wave is more obvious in a test liquid of lower density. In the cell with perpendicularly mounted quartz crystal, the longitudinal wave is mainly reflected back by the inner wall body. The fine structure of plotting of frequency shift (Δf) versus (ρη)^1/2 shows a wave shape, which is different from the well-known linear relationship between of Δf and (ρη)^1/2, where ρ and η are the density and viscosity of the liquid, respectively. And wave-shaped frequency-temperature curves were observed. The longitudinal wave was a kind of potential error source in the PQC measurements. The longitudinal wave effect can be efficiently eliminated by using a rough reflection surface. After eliminating the influence of reflected longitudinal wave, the stability of the sensor PQC was much improved.     

Freezing and Melting of Water Confined in Silica Nanopores

In nanosized pores, liquid water can be thermodynamically stable down to temperatures well below the limit of homogeneous nucleation of bulk water (～235 K). Studies of water in such pores therefore offer an opportunity to reveal the anomalous behavior of deeply supercooled water. Herein we focus on recent studies of the limits of freezing and melting of water in the cylindrical pores of ordered mesoporous silicas with pore diameters in the range of 2–10 nm, based on vapor sorption measurements, calorimetric studies, NMR spectroscopy and cryoporometry, and neutron diffraction studies.     

Phase transition induced by a shock wave in hard-sphere and hard-disk systems

Dynamic phase transition induced by a shock wave in hard-sphere and hard-disk systems is studied on the basis of the system of Euler equations with caloric and thermal equations of state. First, Rankine-Hugoniot conditions are analyzed. The quantitative classification of Hugoniot types in terms of the thermodynamic quantities of the unperturbed state (the state before a shock wave) and the shock strength is made. Especially Hugoniot in typical two possible cases (P-1 and P-2) of the phase transition is analyzed in detail. In the case P-1 the phase transition occurs between a metastable liquid state and a stable solid state, and in the case P-2 the phase transition occurs through coexistence states, when the shock strength changes. Second, the admissibility of the two cases is discussed from a viewpoint of the recent mathematical theory of shock waves, and a rule with the use of the maximum entropy production rate is proposed as the rule for selecting the most probable one among the possible cases, that is, the most suitable constitutive equation that predicts the most probable shock wave. According to the rule, the constitutive equation in the case P-2 is the most promising one in the dynamic phase transition. It is emphasized that hard-sphere and hard-disk systems are suitable reference systems for studying shock wave phenomena including the shock-induced phase transition in more realistic condensed matters.