Diffraction pattern and structure of amorphous solid water at 10 and 77 °K

Abiotic formation of alkylsulfonic acids in interstellar analog ices and implications for their detection on Ryugu

For the last century, the source of sulfur in Earth's very first organisms has remained a fundamental, unsolved enigma. While sulfates and their organic derivatives with sulfur in the S(+VI) oxidation state represent core nutrients in contemporary biochemistry, the limited bioavailability of sulfates during Earth's early Archean period proposed that more soluble S(+IV) compounds served as the initial source of sulfur for the first terrestrial microorganisms. Here, we reveal via laboratory simulation experiments that the three simplest alkylsulfonic acids-water soluble organic S(+IV) compounds-can be efficiently produced in interstellar, sulfur-doped ices through interaction with galactic cosmic rays. This discovery opens a previously elusive path into the synthesis of vital astrobiological significance and untangles fundamental mechanisms of a facile preparation of sulfur-containing, biorelevant organics in extraterrestrial ices; these molecules can be eventually incorporated into comets and asteroids before their delivery and detection on Earth such as in the Murchison, Tagish Lake, and Allende meteorites along with the carbonaceous asteroid Ryugu.

The Importance of Nuclear Quantum Effects on the Thermodynamic and Structural Properties of Low-Density Amorphous Ice: A Comparison with Hexagonal Ice

We study the nuclear quantum effects (NQE) on the thermodynamic properties of low-density amorphous ice (LDA) and hexagonal ice (I_h) at P = 0.1 MPa and T ≥ 25 K. Our results are based on path-integral molecular dynamics (PIMD) and classical MD simulations of H₂O and D₂O using the q-TIP4P/F water model. We show that the inclusion of NQE is necessary to reproduce the experimental properties of LDA and ice I_h. While MD simulations (no NQE) predict that the density ρ(T) of LDA and ice I_h increases monotonically upon cooling, PIMD simulations indicate the presence of a density maximum in LDA and ice I_h. MD and PIMD simulations also predict a qualitatively different T-dependence for the thermal expansion coefficient α_P(T) and bulk modulus B(T) of both LDA and ice I_h. Remarkably, the ρ(T), α_P(T), and B(T) of LDA are practically identical to those of ice I_h. The origin of the observed NQE is due to the delocalization of the H atoms, which is identical in LDA and ice I_h. H atoms delocalize considerably (over a distance ≈ 20-25% of the OH covalent-bond length) and anisotropically (preferentially perpendicular to the OH covalent bond), leading to less linear hydrogen bonds HB (larger HOO angles and longer OO separations) than observed in classical MD simulations.

Acetone-Water Interactions in Crystalline and Amorphous Ice Environments

We present research that systematically examines acetone interacting with various D₂O ices of terrestrial and astrophysical interest using time-resolved, in situ reflection absorption infrared spectroscopy (RAIRS). We examine acetone deposited on top of different D₂O ice films: high-density, nonporous amorphous (np-ASW), and crystalline (CI) films as well as porous amorphous (p-ASW) with various pore morphologies. Analysis of RAIR spectra changes after acetone exposure, and we find that more hydrogen bonding occurs between acetone and p-ASW ices as compared to acetone and np-ASW or CI ices. Hydrogen bonding quantification occurred by two independent RAIR spectral changes: a greater relative intensity of the 1703 cm^-1 feature at low acetone coverage as part of a 14 cm^-1 shift in the C═O region and an ∼30% integrated dangling bond area reduction after acetone exposure. Interestingly, when changing the water structure to be more porous (deposited at 70° compared to 30°), there is a further reduction in the amount of hydrogen bonding that occurs. This suggests that there is a lack of access to surface sites with dangling bonds in the pores as initial layers of acetone block the pores and acetone is unable to diffuse within the structure at low temperatures. In general, these results offer a clearer picture of the mechanisms that can occur when small organic hydrocarbons interact with various icy interfaces; a quantitative understanding of these interactions is essential for the accurate modeling of many astrophysical processes occurring on the surface of icy dust particles.

Structural Analysis of Molecular Materials Using the Pair Distribution Function

This is a review of atomic pair distribution function (PDF) analysis as applied to the study of molecular materials. The PDF method is a powerful approach to study short- and intermediate-range order in materials on the nanoscale. It may be obtained from total scattering measurements using X-rays, neutrons, or electrons, and it provides structural details when defects, disorder, or structural ambiguities obscure their elucidation directly in reciprocal space. While its uses in the study of inorganic crystals, glasses, and nanomaterials have been recently highlighted, significant progress has also been made in its application to molecular materials such as carbons, pharmaceuticals, polymers, liquids, coordination compounds, composites, and more. Here, an overview of applications toward a wide variety of molecular compounds (organic and inorganic) and systems with molecular components is presented. We then present pedagogical descriptions and tips for further implementation. Successful utilization of the method requires an interdisciplinary consolidation of material preparation, high quality scattering experimentation, data processing, model formulation, and attentive scrutiny of the results. It is hoped that this article will provide a useful reference to practitioners for PDF applications in a wide realm of molecular sciences, and help new practitioners to get started with this technique.

On the spectral features of dangling bonds in CH4/H2O amorphous ice mixtures

Dangling bond (DB) bands in IR spectra, above 3600 cm-1, are a source of information on structural properties of amorphous water ice, and especially on ice mixtures of water and other frozen gases. We deal in this paper with the spectroscopic behavior of DB bands of CH4/H2O mixtures. We use ab initio methodology to predict theoretical results which are compared with experimental results. Our model mixtures are created by inserting a variable number of molecules of either species into a cell of appropriate size to reach an initial density of 1 g cm-3, which can be modified by including an empty space at the top, to simulate pores. The cell is taken as a unit cell for a solid state calculation The structure of the mixture is optimized and the IR spectrum is calculated for the converged geometry. We find two different kinds of dangling bonds, in which the O-H stretching responsible for this mode is directed either to an empty space of a pore or towards a nearby CH4 molecule, with which some interaction takes place. The spectral characteristics of these two DB types are clearly different, and follow satisfactorily the pattern observed in experimental spectra. Estimated band strengths for these DB bands are given for the first time.

Gaseous "nanoprobes" for detecting gas-trapping environments in macroscopic films of vapor-deposited amorphous ice

Vapor-deposited amorphous ice, traditionally called amorphous solid water (ASW), is one of the most abundant materials in the universe and a prototypical material for studying physical vapor-deposition processes. Its complex nature arises from a strong tendency to form porous structures combined with complicated glass transition, relaxation, and desorption behavior. To gain further insights into the various gas-trapping environments that exist in ASW and hence its morphology, films in the 25-100 μm thickness range were codeposited with small amounts of gaseous "nanoprobes" including argon, methane, helium, and carbon dioxide. Upon heating in the 95-185 K temperature range, three distinct desorption processes are observed which we attribute to the gas desorption out of open cracks above 100 K, from internal voids that collapse due to the glass transition at ∼125 K and finally from fully matrix-isolated gas induced by the irreversible crystallization to stacking disordered ice (ice Isd) at ∼155 K. Nanoscale films of ASW have only displayed the latter desorption process which means that the first two desorption processes arise from the macroscopic dimensions of our ASW films. Baffling the flow of water vapor toward the deposition plate greatly reduces the first desorption feature, and hence the formation of cracks, but it significantly increases the amount of matrix-isolated gas. The complex nature in which ASW can trap gaseous species is thought to be relevant for a range of cosmological processes.

Resolving solution conformations of the model semi-flexible polyelectrolyte homogalacturonan using molecular dynamics simulations and small-angle x-ray scattering

The conformation of polyelectrolytes in the solution state has long been of interest in polymer science. Herein we utilize all atom molecular dynamics simulations (MD) and small-angle x-ray scattering experiments (SAXS) to elucidate the molecular structure of the model polyelectrolyte homogalacturonan. Several degrees of polymerization were studied and in addition partial methylesterification of the otherwise charge-carrying carboxyl groups was used in order to generate samples with varying intra-chain charge distributions. It is shown that at length scales above around 1nm the conformation of isolated chains has surprisingly little dependence on the charge distribution or the concentration of attendant monovalent salts, reflective of the intrinsic stiffness of the saccharide rings and the dynamical constraints of the glycosidic linkage. Indeed the conformation of isolated chains over all accessible length scales is well described by the atomic coordinates available from fibre diffraction studies. Furthermore, in more concentrated systems it is shown that, after careful analysis of the SAXS data, the form of the inter-particle effects heralded by the emergence of a so-called polyelectrolyte peak, can be extracted, and that this phenomena can be reproduced by multiple chain MD simulations.

Machine learning coarse grained models for water

An accurate and computationally efficient molecular level description of mesoscopic behavior of ice-water systems remains a major challenge. Here, we introduce a set of machine-learned coarse-grained (CG) models (ML-BOP, ML-BOP_dih, and ML-mW) that accurately describe the structure and thermodynamic anomalies of both water and ice at mesoscopic scales, all at two orders of magnitude cheaper computational cost than existing atomistic models. In a significant departure from conventional force-field fitting, we use a multilevel evolutionary strategy that trains CG models against not just energetics from first-principles and experiments but also temperature-dependent properties inferred from on-the-fly molecular dynamics (~ 10's of milliseconds of overall trajectories). Our ML BOP models predict both the correct experimental melting point of ice and the temperature of maximum density of liquid water that remained elusive to-date. Our ML workflow navigates efficiently through the high-dimensional parameter space to even improve upon existing high-quality CG models (e.g. mW model).

Structure of sticky-hard-sphere random aggregates: The viewpoint of contact coordination and tetrahedra

We study more than 10^{4} random aggregates of 10^{6} monodisperse sticky hard spheres each, generated by various static algorithms. Their packing fraction varies from 0.370 up to 0.593. These aggregates are shown to be based on two types of disordered structures: random regular polytetrahedra and random aggregates, the former giving rise to δ peaks on pair distribution functions. Distortion of structural (Delaunay) tetrahedra is studied by two parameters, which show some similarities and some differences in terms of overall tendencies. Isotropy of aggregates is characterized by the nematic order parameter. The overall structure is then studied by distinguishing spheres in function of their contact coordination number (CCN). Distributions of average CCN around spheres of a given CCN value show trends that depend on packing fraction and building algorithms. The radial dependence of the average CCN turns out to be dependent upon the CCN of the central sphere and shows discontinuities that resemble those of the pair distribution function. Moreover, it is shown that structural details appear when the CCN is used as pseudochemical parameter, such as various angular distribution of bond angles, partial pair distribution functions, Ashcroft-Langreth and Bhatia-Thornton partial structure factors. These allow distinguishing aggregates with the same values of packing fraction or average tetrahedral distortion or even similar global pair distribution function, indicative of the great interest of paying attention to contact coordination numbers to study more precisely the structure of random aggregates.

Probing RbBr solvation in freestanding sub-2 nm water clusters

Concentration dependent solvation of RbBr in freestanding sub-2 nm water clusters was studied using core level photoelectron spectroscopy with synchrotron radiation. Spectral features recorded from dilute to saturated clusters indicate that either solvent shared or contact ion pairs are present in increasing amount when the concentration exceeds 2 mol kg^-1. For comparison, spectra from anhydrous RbBr clusters are also presented.

Kinetic Monte Carlo simulations of water ice porosity: extrapolations of deposition parameters from the laboratory to interstellar space

Using an off-lattice kinetic Monte Carlo model we reproduce experimental laboratory trends in the density of amorphous solid water (ASW) for varied deposition angle, rate and surface temperature. Extrapolation of the model to conditions appropriate to protoplanetary disks and interstellar dark clouds indicate that these ices may be less porous than laboratory ices.

Adsorption of PAHs on interstellar ice viewed by classical molecular dynamics

The present work represents a complete description of PAH–ice interaction in the ground electronic state and at low temperature, providing the binding energies and barrier heights necessary to the ongoing improvement of astrochemical models.

Peeling the astronomical onion

This work presents a study of water mobility on interstellar dust grain analogues at temperatures as low as 18 K. The work indicates that water forms pure domains rather than covering the entire grain, thereby leaving bare dust grain surfaces available on which other molecules can adsorb as well as themselves providing surfaces for further adsorption from the interstellar gas.

Predicting the melting temperature of ice-Ih with only electronic structure information as input

The melting temperature of ice-Ih was calculated with only electronic structure information as input by creating a problem-specific force field. The force field, Water model by AFM for Ice and Liquid (WAIL), was developed with the adaptive force matching (AFM) method by fitting to post-Hartree-Fock quality forces obtained in quantum mechanics∕molecular mechanics calculations. WAIL predicts the ice-Ih melting temperature to be 270 K. The model also predicts the densities of ice and water, the temperature of maximum density of water, the heat of vaporizations, and the radial distribution functions for both ice and water in good agreement with experimental measurements. The non-dissociative WAIL model is very similar to a flexible version of the popular TIP4P potential and has comparable computational cost. By customizing to problem-specific configurations with the AFM approach, the resulting model is remarkably more accurate than any variants of TIP4P for simulating ice-Ih and water in the temperature range from 253 K and 293 K under ambient pressure.

Molecular arrangement in water: random but not quite

Water defines life on Earth from the cellular to the terrestrial level. Yet the molecular level arrangement in water is not well understood, posing problems in comprehending its very special chemical, physical and biological properties. Here we present high-resolution x-ray diffraction data for water clearly showing that its molecular arrangement exhibits specific correlations that are consistent with the presence of rings of H(2)O molecules linked together by hydrogen bonds into tetrahedral-like units from a continuous network. This level of molecular arrangement complexity is beyond what a simple 'two-state' model of water (Bernal and Fowler 1933 J. Chem. Phys.1 515-48) could explain. It may not be explained by the recently put forward 'chains-clusters of completely uncorrelated molecules' model (Wernet et al 2004 Science 304 995-9) either. Rather it indicates that water is homogeneous down to the molecular level where different water molecules form tetrahedral units of different perfection and/or participate in rings of different sizes, thus experiencing different local environments. The local diversity of this tetrahedral network coupled to the flexibility of the hydrogen bonds that hold it together may explain well the rich phase diagram of water and why it responds non-uniformly to external stimuli such as, for example, temperature and pressure.

Dielectric relaxation in liquid water: two fractions or two dynamics?

Dielectric relaxation in liquid water is studied using molecular dynamics (MD) simulations in the temperature range of 240 to 340 K at atmospheric pressure. The main dielectric and fast relaxation mode are identified in the spectra of dipole moment autocorrelation functions. The microscopic origin of the fast dielectric relaxation process, which takes place on a time scale of subpicoseconds at room temperature, is discussed. A new hypothesis for the fast dielectric mode is presented. It is based on the assumption of the intrawell rotational relaxation taking place during the waiting period between thermally activated large angle jumps occurring in the course of changing H-bond partners.

Absence of molecular mobility on nano-second time scales in amorphous ice phases

High-resolution neutron backscattering techniques are exploited to study the elastic and quasi-elastic response of the high-density amorphous (HDA), the low-density amorphous (LDA) and the crystalline ice Ic upon temperature changes. Within the temperature ranges of their structural stability (HDA at T < or = 80 K, LDA at T < or = 135 K, ice Ic at T < or = 200 K) the Debye-Waller factors and mean-square displacements characterise all states as harmonic solids. During the transformations HDA --> LDA (T approximately 100 K), LDA --> Ic (T approximately 150 K) and the supposed glass transition with Tg approximately 135 K no relaxation processes can be detected on a time scale t < 4 ns. It can be concluded from coherent scattering measurements (D2O) that LDA starts to recrystallise into ice Ic at T approximately 135 K, i.e. at the supposed Tg. In the framework of the Debye model of harmonic solids HDA reveals the highest Debye temperature among the studied ice phases, which is in full agreement with the lowest Debye level in the generalised density of states derived from time-of-flight neutron scattering experiments. The elastic results at low T indicate the presence of an excess of modes in HDA, which do not obey the Bose statistics.

Liquid water and ices: understanding the structure and physical properties

A review of the structure and some properties of condensed phases of water is given. Since the discovery of the polymorphism of crystalline ice (beginning of the twentieth century), 15 ice modifications have been found and their structures have been determined. If we do not take into consideration proton ordering or disordering, nine distinct crystalline ice modifications in which water molecules retain their individuality are known. In the tenth, ice X, there are no H(2)O molecules. It contains ions (or atoms) of oxygen and hydrogen. The structure of all these modifications is described and information about their fields of stability and about the transition between them is given. It is emphasized that there are ice modifications which are metastable at any temperature and pressure (ices Ic, IV and XII), and many modifications can exist as metastable phases beyond their fields of stability. The ability of water to exist in metastable states is one of its remarkable properties. Several amorphous ice modifications (all of them are metastable) are known. Brief information about their properties and transitions between them is given. At the end of the 1960s the conception of the water structure as a three-dimensional hydrogen-bonded network was conclusively formed. Discovery of the polymorphism of amorphous ices awakened interest in the heterogeneity of the water network. Structural and dynamical heterogeneity of liquid water is discussed in detail. Computer simulation showed that the diffusion coefficient of water molecules in dense regions of the network is lower than in the loose regions, while an increase of density of the entire network gives rise to an increase of diffusion coefficient. This finding contradicts the conceptions associated with the primitive two-state models and can be explained from pressure dependences of melting temperature and of homogeneous nucleation temperature. A brief discussion of the picture of molecular motions in liquid water based on experiment and on computer simulation is given. This picture is still very incomplete. The most fascinating idea that was put forward during the last 20 years was the second critical point conjecture. It is still not clear whether this conjecture corresponds to reality.

Characterization of porosity in vapor-deposited amorphous solid water from methane adsorption

We have characterized the porosity of vapor-deposited amorphous solid water (ice) films deposited at 30-40 K using several complementary techniques such as quartz crystal microgravimetry, UV-visible interferometry, and infrared reflectance spectrometry in tandem with methane adsorption. The results, inferred from the gas adsorption isotherms, reveal the existence of microporosity in all vapor-deposited films condensed from both diffuse and collimated water vapor sources. Films deposited from a diffuse source show a step in the isotherms and much less adsorption at low pressures than films deposited from a collimated source with the difference increasing with film thickness. Ice films deposited from a collimated vapor source at 77 degrees incidence are mesoporous, in addition to having micropores. Remarkably, mesoporosity is retained upon warming to temperatures as high as 140 K where the ice crystallized. The binding energy distribution for methane adsorption in the micropores of ice films deposited from a collimated source peaks at approximately 0.083 eV for deposition at normal incidence and at approximately 0.077 eV for deposition at >45 degrees incidence. For microporous ice, the intensity of the infrared bands due to methane molecules on dangling OH bonds on pore surfaces increases linearly with methane uptake, up to saturation adsorption. This shows that the multilayer condensation of methane does not occur inside the micropores. Rather, filling of the core volume results from coating the pore walls with the first layer of methane, indicating pore widths below a few molecular diameters. For ice deposited at 77 degrees incidence, the increase in intensity of the dangling bond absorptions modified by methane adsorption departs from linearity at large uptakes.

Layer-by-layer growth of thin amorphous solid water films on Pt(111) and Pd(111)

The growth of amorphous solid water (ASW) films on Pt(111) is investigated using rare gas (e.g., Kr) physisorption. Temperature programmed desorption of Kr is sensitive to the structure of thin water films and can be used to assess the growth modes of these films. At all temperatures that are experimentally accessible (20-155 K), the first layer of water wets Pt(111). Over a wide temperature range (20-120 K), ASW films wet the substrate and grow approximately layer by layer for at least the first three layers. In contrast to the ASW films, crystalline ice films do not wet the water monolayer on Pt(111). Virtually identical results were obtained for ASW films on epitaxial Pd(111) films grown on Pt(111). The desorption rates of thin ASW and crystalline ice films suggest that the relative free energies of the films are responsible for the different growth modes. However, at low temperatures, surface relaxation or "transient mobility" is primarily responsible for the relative smoothness of the films. A simple model of the surface relaxation semiquantitatively accounts for the observations.

Structure and hydrogen bonding in liquid and supercritical aqueous NaCl solutions at a pressure of 1000 bar and temperatures up to 500 degrees C: A comprehensive experimental and computational study

The behavior of aqueous 1.1 M NaCl solution at a constant pressure of 1000 bar in the temperature range 25-500 degrees C has been studied with the use of IR absorption, Raman scattering, X-ray diffraction, and molecular dynamics (MD) simulations. The results are compared with the data for pure water under identical external conditions. The main purpose of the experimental and theoretical studies was to understand in what way an electrolyte dissolved in water influences the hydrogen bonding and structural features of water. As was found, the vibrational spectra show no essential difference between the properties of solution and pure water. However, the experimental pair correlation functions and the results of MD simulations present an evidence for very different nature of these substances. A characteristic feature of the structure of NaCl solution is a considerable contribution of strong O-H...Cl- bonds. As the temperature increases, the number of such bonds decreases partially due to a phenomenon of ion pairing, so that at high temperatures the properties of the solution become closer to the properties of water.

Compaction of microporous amorphous solid water by ion irradiation

We have studied the compaction of vapor-deposited amorphous solid water by energetic ions at 40 K. The porosity was characterized by ultraviolet-visible spectroscopy, infrared spectroscopy, and methane adsorption/desorption. These three techniques provide different and complementary views of the structural changes in ice resulting from irradiation. We find that the decrease in internal surface area of the pores, signaled by infrared absorption by dangling bonds, precedes the decrease in the pore volume during irradiation. Our results imply that impacts from cosmic rays can cause compaction in the icy mantles of the interstellar grains, which can explain the absence of dangling bond features in the infrared spectrum of molecular clouds.