The formation of sea ice in polar regions is possible because a salinity gradient or halocline keeps the water column stable despite intense cooling. Here, we demonstrate that a unique water property is central to the maintenance of the polar halocline, namely, that the thermal expansion coefficient (TEC) of seawater increases by one order of magnitude between polar and tropical regions. Using a fully coupled climate model, it is shown that, even with excess precipitations, sea ice would not form at all if the near-freezing temperature TEC was not well below its ocean average value. The leading order dependence of the TEC on temperature is essential to the coexistence of the mid/low-latitude thermally stratified and the high-latitude sea ice–covered oceans that characterize our planet. A key implication is that nonlinearities of water properties have a first-order impact on the global climate of Earth and possibly exoplanets. 

Understanding the mechanism responsible for the protein low-temperature crossover observed at T approximate to 220 K can help us improve current cryopreservation technologies. This crossover is associated with changes in the dynamics of the system, such as in the mean-squared displacement, whereas experimental evidence of structural changes is sparse. Here we investigate hydrated lysozyme proteins by using a combination of wide-angle X-ray scattering and molecular dynamics (MD) simulations. Experimentally we suppress crystallization by accurate control of the protein hydration level, which allows access to temperatures down to T = 175 K. The experimental data indicate that the scattering intensity peak at Q = 1.54 angstrom(-1), attributed to interatomic distances, exhibits temperature-dependent changes upon cooling. In the MD simulations it is possible to decompose the water and protein contributions and we observe that, while the protein component is nearly temperature independent, the hydration water peak shifts in a fashion similar to that of bulk water. The observed trends are analysed by using the water-water and water-protein radial distribution functions, which indicate changes in the local probability density of hydration water.

We study the structural dynamics of liquid water by time-resolved anisotropic x-ray scattering under the optical Kerr effect condition. In this way, we can separate the anisotropic scattering decay of 160 fs from the delayed temperature increase of similar to 0.1 K occurring at 1 ps and quantify transient changes in the O-O pair distribution function. Polarizable molecular dynamics simulations reproduce well the experiment, indicating transient alignment of molecules along the electric field, which shortens the nearest-neighbor distances. In addition, analysis of the simulated water local structure provides evidence that two hypothesized fluctuating water configurations exhibit different polarizability.

Studying the freezing of saltwater on a molecular level is of fundamental importance for improving freeze desalination techniques. In this study, we investigate the freezing process of NaCl solutions using a combination of X-ray diffraction and molecular dynamics simulations (MD) for different salt-water concentrations, ranging from seawater conditions to saturation. A linear superposition model reproduces well the brine rejection due to hexagonal ice Ih formation and allows us to quantify the fraction of ice and brine. Furthermore, upon cooling at T = 233 K, we observe the formation of NaCl center dot 2H(2)O hydrates (hydrohalites), which coexist with ice Ih. MD simulations are utilized to model the formation of NaCl crystal hydrates. From the simulations, we estimate that the salinity of the newly produced ice is 0.5% mass percent (m/m) due to ion inclusions, which is within the salinity limits of fresh water. In addition, we show the effect of ions on the local ice structure using the tetrahedrality parameter and follow the crystallite formation using the ion coordination parameter and cluster analysis.

Accommodation of vapor-phase water molecules into ice crystal surfaces is a fundamental process controlling atmospheric ice crystal growth. Experimental studies investigating the accommodation process with various techniques report widely spread values of the water accommodation coefficient on ice, αice, and the results on its potential temperature dependence are inconclusive. We run molecular dynamics simulations of molecules condensing onto the basal plane of ice Ih using the TIP4P/Ice empirical force field and characterize the accommodated state from this molecular perspective, utilizing the interaction energy, the tetrahedrality order parameter, and the distance below the instantaneous interface as criteria. Changes of the order parameter turn out to be a suitable measure to distinguish between the surface and bulk states of a molecule condensing onto the disordered interface. In light of the findings from the molecular dynamics, we discuss and re-analyze a recent experimental data set on αice obtained with an environmental molecular beam (EMB) setup [Kong, X.; J. Phys. Chem. A 2014, 118 (22), 3973−3979] using kinetic molecular flux modeling, aiming at a more comprehensive picture of the accommodation process from a molecular perspective. These results indicate that the experimental observations indeed cannot be explained by evaporation alone. At the same time, our results raise the issue of rapidly growing relaxation times upon decreasing temperature, challenging future experimental efforts to cover relevant time scales. Finally, we discuss the relevance of the water accommodation coefficient on ice in the context of atmospheric cloud particle growth processes. 

Based on recent experimental data that can be interpreted as indicating the presence of specific structures in liquid water, we build and optimize two structural models which we compare with the available experimental data. To represent the proposed high-density liquid structures, we use a model consisting of chains of water molecules, and for low-density liquid, we investigate fused dodecahedra as templates for tetrahedral fluctuations. The computed infrared spectra of the models are in very good agreement with the extracted experimental spectra for the two components, while the extracted structures from molecular dynamics (MD) simulations give spectra that are intermediate between the experimentally derived spectra. Computed x-ray absorption and emission spectra as well as the O-O radial distribution functions of the proposed structures are not contradicted by experiment. The stability of the proposed dodecahedral template structures is investigated in MD simulations by seeding the starting structure, and remnants found to persist on an similar to 30 ps time scale. We discuss the possible significance of such seeds in simulations and whether they can be viable candidates as templates for structural fluctuations below the compressibility minimum of liquid water.

Solid particles in protoplanetary discs can grow by direct vapour deposition outside of ice lines. The presence of microscopic silicate particles may nevertheless hinder growth into large pebbles, since the available vapour is deposited predominantly on the small grains that dominate the total surface area. Experiments on heterogeneous ice nucleation, performed to understand ice clouds in the Martian atmosphere, show that the formation of a new ice layer on a silicate surface requires a substantially higher water vapour pressure than the deposition of water vapour on an existing ice surface. In this paper, we investigate how the difference in partial vapour pressure needed for deposition of vapour on water ice versus heterogeneous ice nucleation on silicate grains influences particle growth close to the water ice line. We developed and tested a dynamical 1D deposition and sublimation model, where we include radial drift, sedimentation, and diffusion in a turbulent protoplanetary disc. We find that vapour is deposited predominantly on already ice-covered particles, since the vapour pressure exterior of the ice line is too low for heterogeneous nucleation on bare silicate grains. Icy particles can thus grow to centimetre-sized pebbles in a narrow region around the ice line, whereas silicate particles stay dust-sized and diffuse out over the disc. The inhibition of heterogeneous ice nucleation results in a preferential region for growth into planetesimals close to the ice line where we find large icy pebbles. The suppression of heterogeneous ice nucleation on silicate grains may also be the mechanism behind some of the observed dark rings around ice lines in protoplanetary discs, as the presence of large ice pebbles outside ice lines leads to a decrease in the opacity there.

Atmospheric black carbon (BC) particles are a concern due to their impact on air quality and climate. Their net climate effect 15 is, however, still uncertain. This uncertainty is partly related to the contribution of coated BC-particles to the global CCN budgets. In this study, laboratory measurements were performed to investigate cloud condensation nuclei (CCN) activity of BC (Regal black) particles, in pure state or coated through evaporating and subsequent condensation of glutaric acid, levoglucosan (both water-soluble organics) or oleic acid (an organic compound with low solubility). A combination of Soot Particle Aerosol Mass Spectrometer (SP-AMS) measurements and size distribution measurements with Scanning Mobility 20 Particle Sizer (SMPS) showed that the studied BC particles were nearly spherical agglomerates with a fractal dimension of 2.79 and that they were coated evenly by the organic species. The CCN activity of BC particles increased after coating with all the studied compounds and was governed by the fraction of organic material. The CCN activation of the BC particles coated by glutaric acid and levoglucosan were in good agreement with the theoretical calculations using shell-and-core model, which is based on a combination of the CCN activities of the pure compounds. The oleic acid coating enhanced the CCN 25 activity of the BC particles, even though the pure oleic acid particles were CCN inactive. The surprising effect of oleic acid might be related to the arrangement of the oleic acid molecules on the surface of the BC cores or other surface phenomena facilitating water condensation onto the coated particles. Our results show potential in accurately predicting the CCN activity of atmospheric BC coated with organic species by present theories, given that the identities and amount of the coating species are known. Furthermore, our results suggest that even relatively thin soluble coatings (around 2 nm for the compounds studied here) are enough to make the insoluble BC particles CCN active at typical atmospheric supersaturations and thus be efficiently taken up by cloud droplets. This highlights the need of an accurate description of the composition of atmospheric particles containing BC to unravel their net impact on climate.

The structure factor and oxygen-oxygen pair distribution functions of amorphous ices at liquid nitrogen temperature (T = 77 K) have been derived from wide-angle X-ray scattering (WAXS) up to interatomic distances of r = 23 angstrom, where local structure differences between the amorphous ices can be seen for the entire range. The distances to the first coordination shell for low-, high-, and very-high-density amorphous ice (LDA, HDA, VHDA) were determined to be 2.75, 2.78, and 2.80 angstrom, respectively, with high accuracy due to measurements up to a large momentum transfer of 23 angstrom(-1). Similarities in pair-distribution functions between LDA and supercooled water at 254.1 K, HDA and liquid water at 365.9 K, and VHDA and high-pressure liquid water were found up to around 8 angstrom, but beyond that at longer distances, the similarities were lost. In addition, the structure of the high-density amorphous ices was compared to high-pressure crystalline ices IV, IX, and XII, and conclusions were drawn about the local ordering.

Water exists in high- and low-density amorphous ice forms (HDA and LDA), which could correspond to the glassy states of high(HDL) and low-density liquid (LDL) in the metastable part of the phase diagram. However, the nature of both the glass transition and the high-to-low-density transition are debated and new experimental evidence is needed. Here we combine wide-angle X-ray scattering (WAXS) with X-ray photon-correlation spectroscopy (XPCS) in the small-angle X-ray scattering (SAXS) geometry to probe both the structural and dynamical properties during the high-to-low-density transition in amorphous ice at 1 bar. By analyzing the structure factor and the radial distribution function, the coexistence of two structurally distinct domains is observed at T = 125 K. XPCS probes the dynamics in momentum space, which in the SAXS geometry reflects structural relaxation on the nanometer length scale. The dynamics of HDA are characterized by a slow component with a large time constant, arising from viscoelastic relaxation and stress release from nanometer-sized heterogeneities. Above 110 K a faster, strongly temperature-dependent component appears, with momentum transfer dependence pointing toward nanoscale diffusion. This dynamical component slows down after transition into the low-density form at 130 K, but remains diffusive. The diffusive character of both the high- and low-density forms is discussed among different interpretations and the results are most consistent with the hypothesis of a liquid-liquid transition in the ultraviscous regime.