Muonium diffusion in ice

Muonium (Mu, μ⁺e⁻) is generally regarded as a light isotope of hydrogen. The procession signals of muonium in single crystals of H₂O and D₂O ice have been studied from 8 to 263 K using the muon spin rotation (μSR) technique. Transverse spin relaxation rates have been extracted and interpreted in terms of modulation of the dipolar interaction between muonium and the protons/deuterons in the lattice by translational diffusion of muonium. In contrast to the situation for H and a previous claim for Mu, muonium is found to be diffusing at temperatures as low as 8 K. An activation energy of 40 meV is obtained by fitting the highest temperature data to an Arrhenius expression. At low temperature muonium is thought to diffuse by quantum tunnelling.     

Partial spin depolarization of muonium in ice

Muon spin precession signals arising from both muonium and a diamagnetic muon species have been studied in single crystal of ice over the temperature range 90–263 K. Conversion of initial signal amplitudes to fractions of muon polarization reveals that for temperatures above 200 K part of the original muon polarization is unaccounted for. Such a missing fraction is well known for liquid water. but was not found in earlier work on polycrystalline ice. Muonium signals in Polycrystalline ice were reanalyzed using a non-exponential decay function appropriate to the powder pattern spectrum. Smaller muonium fractions were found. consistent with those determined for single crystall, thus confirming the existence of the missing fraction. The origin of the missing fraction is discussed, and it is proposed that non-reactive spin exchange encounters between muonium and hydrogen atoms may be the cause.     

Quantum path integral simulation of isotope effects in the melting temperature of ice Ih

The isotope effect in the melting temperature of ice Ih has been studied by free energy calculations within the path integral formulation of statistical mechanics. Free energy differences between isotopes are related to the dependence of their kinetic energy on the isotope mass. The water simulations were performed by using the q-TIP4P/F model, a point charge empirical potential that includes molecular flexibility and anharmonicity in the OH stretch of the water molecule. The reported melting temperature at ambient pressure of this model (T=251?K) increases by 6.5±0.5 and 8.2±0.5?K upon isotopic substitution of hydrogen by deuterium and tritium, respectively. These temperature shifts are larger than the experimental ones (3.8 and 4.5 K, respectively). In the classical limit, the melting temperature is nearly the same as that for tritiated ice. This unexpected behavior is rationalized by the coupling between intermolecular interactions and molecular flexibility. This coupling makes the kinetic energy of the OH stretching modes larger in the liquid than in the solid phase. However, the opposite behavior is found for intramolecular modes, which display larger kinetic energy in ice than in liquid water.     

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.     

Consequences of chain networks on thermodynamic, dielectric and structural properties for liquid water

A vast array of experimental data on water provides a global view of the liquid that implicates its tetrahedral hydrogen-bonding network as the unifying molecular connection to its observed structural, thermodynamic, and dielectric property trends with temperature. Recently the classification of water as a tetrahedral liquid has been challenged based on X-ray absorption (XAS) experiments on liquid water (Ph. Wernet et al., Science, 2004, 304, 995), which have been interpreted to show a hydrogen-bonding network that replaces tetrahedral structure with chains or large rings of water molecules. We examine the consequences of tetrahedral vs. chain networks using three different modified water models that exhibit a local hydrogen bonding environment of two hydrogen bonds (2HB) and therefore networks of chains. Using these very differently parameterized models we evaluate their bulk densities, enthalpies of vaporization, heat capacities, isothermal compressibilities, thermal expansion coefficients, and dielectric constants, over the temperature range of 235-323 K. We also evaluate the entropy of the 2HB models at room temperature and whether such models support an ice I(h) structure. All show poor agreement with experimentally measured thermodynamic and dielectric properties over the same temperature range, and behave similarly in most respects to normal liquids.     

Electrofreezing of confined water

We report results from molecular dynamics simulations of the freezing transition of TIP5P water molecules confined between two parallel plates under the influence of a homogeneous external electric field, with magnitude of 5 V/nm, along the lateral direction. For water confined to a thickness of a trilayer we find two different phases of ice at a temperature of T=280 K. The transformation between the two, proton-ordered, ice phases is found to be a strong first-order transition. The low-density ice phase is built from hexagonal rings parallel to the confining walls and corresponds to the structure of cubic ice. The high-density ice phase has an in-plane rhombic symmetry of the oxygen atoms and larger distortion of hydrogen bond angles. The short-range order of the two ice phases is the same as the local structure of the two bilayer phases of liquid water found recently in the absence of an electric field [J. Chem. Phys. 119, 1694 (2003)]. These high- and low-density phases of water differ in local ordering at the level of the second shell of nearest neighbors. The results reported in this paper, show a close similarity between the local structure of the liquid phase and the short-range order of the corresponding solid phase. This similarity might be enhanced in water due to the deep attractive well characterizing hydrogen bond interactions. We also investigate the low-density ice phase confined to a thickness of 4, 5, and 8 molecular layers under the influence of an electric field at T=300 K. In general, we find that the degree of ordering decreases as the distance between the two confining walls increases.     

Graphene Oxide Restricts Growth and Recrystallization of Ice Crystals

We show graphene oxide (GO) greatly suppresses the growth and recrystallization of ice crystals, and ice crystals display a hexagonal shape in the GO dispersion. Preferred adsorption of GO on the ice crystal surface in liquid water leads to curved ice crystal surface. Therefore, the growth of ice crystal is suppressed owing to the Gibbs–Thompson effect, that is, the curved surface lowers the freezing temperature. Molecular dynamics simulation analysis reveals that oxidized groups on the basal plane of GO form more hydrogen bonds with ice in comparison with liquid water because of the honeycomb hexagonal scaffold of graphene, giving a molecular‐level mechanism for controlling ice formation. Application of GO for cryopreservation shows that addition of only 0.01 wt % of GO to a culture medium greatly increases the motility (from 24.3 % to 71.3 %) of horse sperms. This work reports the control of growth of ice with GO, and opens a new avenue for the application of 2D materials.     

Interaction of atomic and molecular deuterium with a nonporous amorphous water ice surface between 8 and 30 K

Molecular and atomic interactions of hydrogen on dust grains covered with ice at low temperatures are key mechanisms for star formation and chemistry in dark interstellar clouds. We have experimentally studied the interaction of atomic and molecular deuterium on nonporous amorphous water ice surfaces between 8 and 30 K, in conditions compatible with an extrapolation to an astrophysical context. The adsorption energy of D(2) presents a wide distribution, as already observed on porous water ice surfaces. At low coverage, the sticking coefficient of D(2) increases linearly with the number of deuterium molecules already adsorbed on the surface. Recombination of atomic D occurs via a prompt reaction that releases molecules into the gas phase. Part of the newly formed molecules are in vibrationally excited states (v=1-7). The atomic recombination efficiency increases with the presence of D(2) molecules already adsorbed on the water ice, probably because these increase the sticking coefficient of the atoms, as in the case of incident D(2). We have measured the atomic recombination efficiency in the presence of already absorbed D(2), as it is expected to occur in the interstellar medium. The recombination efficiency decreases rapidly with increasing temperature and is zero at 13 K. This allows us to estimate an upper limit to the value of the atom adsorption energy E(a) approximately 29 meV, in agreement with previous calculations.     

Phase diagrams and water activities of aqueous ammonium salts of malonic acid

Malonic acid has been observed in the free troposphere and as a component of tropospheric aerosol, among other dicarboxylic acids. These aerosols can uptake ammonia, which partially or completely neutralizes the acids. Therefore, the impact of ammoniated dicarboxylic acids on the phases that can exist in aerosols at atmospheric temperatures needs investigation. To that end, the low temperature, solid/liquid phase diagrams of ammonium hydrogen malonate/water, ammonium malonate/water, and triammonium hydrogen malonate/water have been investigated with differential scanning calorimetry and infrared spectroscopy of thin films. Results show that the order of increasing solubility is triammonium hydrogen malonate, ammonium hydrogen malonate, malonic acid, and ammonium malonate. We have also determined a hydrate may form in the ammonium malonate system and decompose below 240 K. We report water activities at the ice melting points for each system up to the respective eutectic concentrations, and find for a given mole fraction of water, increasing ammonium content leads to decreasing water activity coefficients.     

Neutron diffraction study of water freezing on aircraft engine combustor soot

The study of the formation of condensation trails and cirrus clouds on aircraft emitted soot particles is important because of its possible effects on climate. In the present work we studied the freezing of water on aircraft engine combustor (AEC) soot particles under conditions of pressure and temperature similar to the upper troposphere. The microstructure of the AEC soot was found to be heterogeneous containing both primary particles of soot and metallic impurities (Fe, Cu, and Al). We also observed various surface functional groups such as oxygen-containing groups, including sulfate ions, that can act as active sites for water adsorption. Here we studied the formation of ice on the AEC soot particles by using neutron diffraction. We found that for low amount of adsorbed water, cooling even up to 215 K did not lead to the formation of hexagonal ice. Whereas, larger amount of adsorbed water led to the coexistence of liquid water (or amorphous ice) and hexagonal ice (I(h)); 60% of the adsorbed water was in the form of ice I(h) at 255 K. Annealing of the system led to the improvement of the crystal quality of hexagonal ice crystals as demonstrated from neutron diffraction.     

Melting points and thermal expansivities of proton-disordered hexagonal ice with several model potentials

A method of free energy calculation is proposed, which enables to cover a wide range of pressure and temperature. The free energies of proton-disordered hexagonal ice (ice Ih) and liquid water are calculated for the TIP4P [J. Chem. Phys. 79, 926 (1983)] model and the TIP5P [J. Chem. Phys. 112, 8910 (2000)] model. From the calculated free energy curves, we determine the melting point of the proton-disordered hexagonal ice at 0.1 MPa (atmospheric pressure), 50 MPa, 100 MPa, and 200 MPa. The melting temperatures at atmospheric pressure for the TIP4P ice and the TIP5P ice are found to be about T(m)=229 K and T(m)=268 K, respectively. The melting temperatures decrease as the pressure is increased, a feature consistent with the pressure dependence of the melting point for realistic proton-disordered hexagonal ice. We also calculate the thermal expansivity of the model ices. Negative thermal expansivity is observed at the low temperature region for the TIP4P ice, but not for the TIP5P ice at the ambient pressure.     

Comparative Adsorption of Acetone on Water and Ice Surfaces

Small organic molecules on ice and water surfaces are ubiquitous in nature and play a crucial role in many environmentally relevant processes. Herein, we combine surface‐specific vibrational spectroscopy and a controllable flow cell apparatus to investigate the molecular adsorption of acetone onto the basal plane of single‐crystalline hexagonal ice with a large surface area. By comparing the adsorption of acetone on the ice/air and the water/air interface, we observed two different types of acetone adsorption, as apparent from the different responses of both the free O?H and the hydrogen‐bonded network vibrations for ice and liquid water. Adsorption on ice occurs preferentially through interactions with the free OH group, while the interaction of acetone with the surface of liquid water appears less specific.     

Is it cubic? Ice crystallization from deeply supercooled water

Ice crystallized below 200 K has the diffraction pattern of a faulty cubic ice, and not of the most stable hexagonal ice polymorph. The origin and structure of this faulty cubic ice, presumed to form in the atmosphere, has long been a puzzle. Here we use large-scale molecular dynamics simulations with the mW water model to investigate the crystallization of water at 180 K and elucidate the development of cubic and hexagonal features in ice as it nucleates, grows and consolidates into crystallites with characteristic dimensions of a few nanometres. The simulations indicate that the ice crystallized at 180 K contains layers of cubic ice and hexagonal ice in a ratio of approximately 2 to 1. The stacks of hexagonal ice are very short, mostly one and two layers, and their frequency does not seem to follow a regular pattern. In spite of the high fraction of hexagonal layers, the diffraction pattern of the crystals is, as in the experiments, almost identical to that of cubic ice. Stacking of cubic and hexagonal layers is observed for ice nuclei with as little as 200 water molecules, but a preference for cubic ice is already well developed in ice nuclei one order of magnitude smaller: the critical ice nuclei at 180 K contain approximately ten water molecules in their core and are already rich in cubic ice. The energies of the cubic-rich and hexagonal-rich nuclei are indistinguishable, suggesting that the enrichment in cubic ice does not have a thermodynamic origin.     

Onsager heat of transport for water vapour at the surface of water and ice: thermal accommodation coefficients for water vapour on a stainless-steel surface

The Onsager heat of transport Q* has been measured for water vapour at the surface of water, supercooled water, and ice, over the temperature range -8 to +10 degrees C. For liquid water, Q* is constant at -24.7 +/- 3.6 kJ mol(-1) (two standard deviations) over the pressure range 4-9.5 Torr. Provided the ice is suitably aged, the |Q*| values are very similar for water and ice, a result which is consistent with the presence of a liquid-like layer at the surface of ice. The values are slightly larger for ice, in proportion to the ratio of the heat of sublimation of ice to the heat of vaporization of the liquid. Departures from linearity of plots of P against DeltaT are attributed to temperature jumps at the surface of the dry upper plate. Hence jump coefficients and thermal accommodation coefficients have been derived as a function of temperature for collisions of water molecules with type-304 stainless steel.     

Robust three-body water simulation model

The most common potentials used in classical simulations of liquid water assume a pairwise additive form. Although these models have been very successful in reproducing many properties of liquid water at ambient conditions, none is able to describe accurately water throughout its complicated phase diagram. The primary reason for this is the neglect of many-body interactions. To this end, a simulation model with explicit three-body interactions was introduced recently [R. Kumar and J. L. Skinner, J. Phys. Chem. B 112, 8311 (2008)]. This model was parameterized to fit the experimental O-O radial distribution function and diffusion constant. Herein we reparameterize the model, fitting to a wider range of experimental properties (diffusion constant, rotational correlation time, density for the liquid, liquid/vapor surface tension, melting point, and the ice Ih density). The robustness of the model is then verified by comparing simulation to experiment for a number of other quantities (enthalpy of vaporization, dielectric constant, Debye relaxation time, temperature of maximum density, and the temperature-dependent second and third virial coefficients), with good agreement.     

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.     

Isotope effect in the diffusion of hydrogen and deuterium in polymers

Temperature dependences of diffusion and permeation coefficients of hydrogen and deuterium in glassy and rubbery polymer films have been measured. The size of the free volume element in rubbery polymers has been calculated according to the theory of Frisch and Rogers for the quantum isotope effect, but the free volume is too large for precise calculation below the glass-transition temperature. The cooperative movement of segments is also discussed using the ratio of preexponential factors for diffusion mechanisms above and below the glass-transition temperature.     

Positive muons in condensed phase ammonia

The chemical behavior of positive muons in condensed phase ammonia has been investigated in order to elucidate the phase and temperature effects on the chemical and physical behavior of the muon and muonium formation in a simple binary compound. Diamagnetic muon yield (P_D) was constant at 0.67±0.01 in both solid and liquid above 125 K. Muonium formation in solids were observed above 100 K with slow muonium spin relaxation. In liquids, the muonium yield and its spin relaxation rate showed temperature dependence. Addition of metallic sodium increased P_D in liquids.     

Carbodiimide production from cyanamide by UV irradiation and thermal reaction on amorphous water ice

Cyanamide (NH(2)CN), an interstellar molecule, is a relevant molecule in prebiotic chemistry, because it can be converted into urea in liquid water. Carbodiimide (HNCNH), the most stable cyanamide isomer, is able to assemble amino acids into peptides. In this work, using FTIR spectroscopy, we show that carbodiimide can be formed from cyanamide at low temperature (10 K), by a photochemical process in argon matrix, in water matrix, or in solid film. We also report experimental evidence about the carbodiimide formation when cyanamide is condensed at low temperature (50-140 K) on an amorphous water ice surface, or when it is trapped in the water ice. The water ice acts as a catalyst. This isomerization reaction occurs at low temperature (T < 100 K), which agrees with those expected in the interstellar clouds composed of dust grains in which water is the most predominant compound. Finally, the hydrolysis reaction of cyanamide or carbodiimide leading to urea or isourea formation is not observed under our experimental conditions.     

Temperature dependence of solvation dynamics of probe molecules in methanol-doped ice and in liquid ethanol

We have studied the solvation statics and dynamics of coumarin 343 and a strong photoacid (pK* approximately 0.7) 2-naphthol-6, 8-disulfonate (2N68DS) in methanol-doped ice (1% molar concentration of methanol) and in cold liquid ethanol in the temperature range of 160-270 K. Both probe molecules show a relatively fast solvation dynamics in ice, ranging from a few tens of picoseconds at about 240 K to nanoseconds at about 160 K. At about 160 K in doped ice, we observe a sharp decrease of the dynamic Stokes shift of both coumarin 343 and 2N68DS. Its value is approximately only 200 cm-1 at approximately 160 K compared to about 1100 cm-1 at T >/= 200 K (at times longer than t > 10 ps). We find a good correlation between the inefficient and slow excited-state proton-transfer rate at low-temperature ice, T < 180 K, and the dramatic decrease of the solvation energy, as measured by the dynamic band shift, at these low temperatures. We find that the average solvation rate in ice is similar to its value in liquid ethanol at all given temperatures in the range of 200-250 K. The surprisingly fast solvation rate in ice is explained by the relatively large freedom of the water hydrogen rotation in ice Ih.