Search

Ice crystals in the atmosphere nucleate from supercooled liquid water and grow by vapor uptake. The structure of the ice polymorph grown has strong impact on the morphology and light scattering of the ice crystals, modulates the amount of water vapor in ice clouds, and can impact the molecular uptake and reactivity of atmospheric aerosols. Experiments and molecular simulations indicate that ice nucleated and grown from deeply supercooled liquid water is metastable stacking disordered ice. The ice polymorph grown from vapor has not yet been determined. Here we use large-scale molecular simulations to determine the structure of ice that grows as a result of uptake of water vapor in the temperature range relevant to cirrus and mixed-phase clouds, elucidate the molecular mechanism of the formation of ice at the vapor interface, and compute the free energy difference between cubic and hexagonal ice interfaces with vapor. We find that vapor deposition results in growth of stacking disordered ice only under conditions of extreme supersaturation, for which a nonequilibrium liquid layer completely wets the surface of ice. Such extreme conditions have been used to produce stacking disordered frost ice in experiments and may be plausible in the summer polar mesosphere. Growth of ice from vapor at moderate supersaturations in the temperature range relevant to cirrus and mixed-phase clouds, from 200 to 260 K, produces exclusively the stable hexagonal ice polymorph. Cubic ice is disfavored with respect to hexagonal ice not only by a small penalty in the bulk free energy (3.6 ± 1.5 J molv-1 at 260 K) but also by a large free energy penalty at the ice-vapor interface (89.7 ± 12.8 J molv-1 at 260 K). The latter originates in higher vibrational entropy of the hexagonal-terminated ice-vapor interface. We predict that the free energy penalty against the cubic ice interface should decrease strongly with temperature, resulting in some degree of stacking disorder in ice grown from vapor in the tropical tropopause layer, and in polar stratospheric and noctilucent clouds. Our findings support and explain the evolution of the morphology of ice crystals from hexagonal to trigonal symmetry with decreasing temperature, as reported by experiments and in situ measurements in clouds. We conclude that selective growth of the elusive cubic ice polymorph by manipulation of the interfacial properties can likely be achieved at the ice-liquid interface but not at the ice-vapor interface. [ABSTRACT FROM AUTHOR]

Atmospheric aerosols have a strong influence on Earth's climate. Elucidating the physical state and intemal structure of atmospheric aqueous aerosols is essential to predict their gas and water uptake, and the locus and rate of atmospherically important heterogeneous reactions. Ultrafine aerosols with sizes between 3 and 15 nm have been detected in large numbers in the troposphere and tropopause. Nanoscopic aerosols arising from bubble bursting of natural and artificial seawater have been identified in laboratory and field experiments. The internal structure and phase state of these aerosols, however, cannot yet be determined in experiments. Here we use molecular simulations to investigate the phase behavior and internal structure of liquid, vitrified, and crystallized water-salt ultrafine aerosols with radii from 2.5 to 9.5 nm and with up to 10% moles of ions. We find that both ice crystallization and vitrification of the nanodroplets lead to demixing of pure water from the solutions. Vitrification of aqueous nanodroplets yields nanodomains of pure low-density amorphous ice in coexistence with vitrified solute rich aqueous glass. The melting temperature of ice in the aerosols decreases monotonically with an increase of solute fraction and decrease of radius. The simulations reveal that nucleation of ice occurs homogeneously at the subsurface of the water-salt nanoparticles. Subsequent ice growth yields phase-segregated, internally mixed, aerosols with two phases in equilibrium: a concentrated water-salt amorphous mixture and a spherical cap-like ice nanophase. The surface of the crystallized aerosols is heterogeneous, with ice and solution exposed to the vapor. Free energy calculations indicate that as the concentration of salt in the particles, the advance of the crystallization, or the size of the particles increase, the stability of the spherical cap structure increases with respect to the alternative structure in which a core of ice is fuUy surrounded by solution. We predict that micrometer-sized particles and nanoparticles have the same equilibrium internal structure. The variation of liquid-vapor surface tension with solute concentration is a key factor in determining whether a solution-embedded ice core or vapor-exposed ice cap is the equilibrium structure of the aerosols. In agreement vnth experiments, we predict that the structure of mixed-phase HNO3-water particles, representative of polar stratospheric clouds, consists of an ice core surrounded by freeze-concentrated solution. The results of this work are important to determine the phase state and internal structure of sea spray ultrafine aerosols and other mixed-phase particles under atmospherically relevant conditions. [ABSTRACT FROM AUTHOR]

Atmospheric aerosols can promote the heterogeneous nucleation of ice, impacting the radiative properties of clouds and Earth's climate. The experimental investigation of heterogeneous freezing of water droplets by carbonaceous particles reveals widespread ice freezing temperatures. It is not known which structural and chemical characteristics of soot account for the variability in ice nucleation efficiency. Here we use molecular dynamics simulations to investigate the nucleation of ice from liquid water in contact with graphitic surfaces. We find that atomically flat carbon surfaces promote heterogeneous nucleation of ice, while molecularly rough surfaces with the same hydrophobicity do not. Graphitic surfaces and other surfaces that promote ice nucleation induce layering in the interfacial water, suggesting that the order imposed by the surface on liquid water may play an important role in the heterogeneous nucleation mechanism. We investigate a large set of graphitic surfaces of various dimensions and radii of curvature and find that variations in nanostructures alone could account for the spread in the freezing temperatures of ice on soot in experiments. We conclude that a characterization of the nanostructure of soot is needed to predict its ice nucleation efficiency. [ABSTRACT FROM AUTHOR]

Cold-adapted organisms produce antifreeze and ice-nucleating proteins to prevent and promote ice formation. The crystal structure of hyperactive bacterial antifreeze protein (AFP) MpAFP suggests that this protein binds ice through an anchored clathrate motif. It is not known whether other hyperactive AFPs and ice-nucleating proteins (INPs) use the same motif to recognize or nucleate ice. Here we use molecular simulations to elucidate the ice-binding motifs of hyperactive insect AFPs and a model INP of Pseudomonas syringae. We find that insect AFPs recognize ice through anchored clathrate motifs distinct from that of MpAFP. By performing simulations of ice nucleation by PsINP, we identify two distinct ice-binding sites on opposite sides of the β-helix. The ice-nucleating sequences identified in the simulations agree with those previously proposed for the closely related INP of Pseudomonas borealis based on the structure of the protein. The simulations indicate that these sites have comparable ice-nucleating efficiency, but distinct binding motifs, controlled by the amino acid sequence: one is an anchored clathrate and the other ice-like. We conclude that anchored clathrate and ice-like motifs can be equally effective for binding proteins to ice and promoting ice nucleation.

The ice-air interface is an important locus of environmental chemical reactions. The structure and dynamics of the ice surface impact the uptake of trace gases and kinetics of reactions in the atmosphere and snowpack. At tropospheric temperatures, the ice surface is partially premelted. Experiments indicate that ions increase the liquidity of the ice surface but hydrophilic organics do not. However, it is not yet known the extent of the perturbation solutes induce at the ice surface and what is the role of the disordered liquid-like layer in modulating the interaction between solutes and their mobility and aggregation at the ice surface. Here we use large-scale molecular simulations to investigate the effect of ions and glyoxal, one of the most abundant oxygenated volatile organic compounds in the atmosphere, on the structure, dynamics, and solvation properties of the ice surface. We find that the premelted surface of ice has unique solvation properties, different from those of liquid water. The increase in surface liquidity resulting from the hydration of ions leads to a water-mediated attraction of ions at the ice surface. Glyoxal molecules, on the other hand, perturb only slightly the surface of ice and do not experience water-driven attraction. They nonetheless accumulate as dry agglomerates at the ice surface, driven by direct interactions between the organic molecules. The enhanced attraction and clustering of ions and organics at the ice surface may play a significant role in modulating the mechanism and rate of heterogeneous chemical reactions occurring at the surface of atmospheric ice particles.

Heterogeneous nucleation of ice induced by organic materials is of fundamental importance for climate, biology, and industry. Among organic ice-nucleating surfaces, monolayers of long chain alcohols are particularly effective, while monolayers of fatty acids are significantly less so. As these monolayers expose to water hydroxyl groups with an order that resembles the one in the basal plane of ice, it was proposed that lattice matching between ice and the surface controls their ice-nucleating efficiency. Organic monolayers are soft materials and display significant fluctuations. It has been conjectured that these fluctuations assist in the nucleation of ice. Here we use molecular dynamic simulations and laboratory experiments to investigate the relationship between the structure and fluctuations of hydroxylated organic surfaces and the temperature at which they nucleate ice. We find that these surfaces order interfacial water to form domains with ice-like order that are the birthplace of ice. Both mismatch and fluctuations decrease the size of the preordered domains and monotonously decrease the ice freezing temperature. The simulations indicate that fluctuations depress the freezing efficiency of monolayers of alcohols or acids to half the value predicted from lattice mismatch alone. The model captures the experimental trend in freezing efficiencies as a function of chain length and predicts that alcohols have higher freezing efficiency than acids of the same chain length. These trends are mostly controlled by the modulation of the structural mismatch to ice. We use classical nucleation theory to show that the freezing efficiencies of the monolayers are directly related to their free energy of binding to ice. This study provides a general framework to relate the equilibrium thermodynamics of ice binding to a surface and the nonequilibrium ice freezing temperature and suggests that these could be predicted from the structure of interfacial water.