Surface structure of thin ice films

Imaging surface structure and premelting of ice Ih with atomic resolution

Ice surfaces are closely relevant to many physical and chemical properties, such as melting, freezing, friction, gas uptake and atmospheric reaction1-8. Despite extensive experimental and theoretical investigations9-17, the exact atomic structures of ice interfaces remain elusive owing to the vulnerable hydrogen-bonding network and the complicated premelting process. Here we realize atomic-resolution imaging of the basal (0001) surface structure of hexagonal water ice (ice Ih) by using qPlus-based cryogenic atomic force microscopy with a carbon monoxide-functionalized tip. We find that the crystalline ice-Ih surface consists of mixed Ih- and cubic (Ic)-stacking nanodomains, forming 19 × 19 periodic superstructures. Density functional theory reveals that this reconstructed surface is stabilized over the ideal ice surface mainly by minimizing the electrostatic repulsion between dangling OH bonds. Moreover, we observe that the ice surface gradually becomes disordered with increasing temperature (above 120 Kelvin), indicating the onset of the premelting process. The surface premelting occurs from the defective boundaries between the Ih and Ic domains and can be promoted by the formation of a planar local structure. These results put an end to the longstanding debate on ice surface structures and shed light on the molecular origin of ice premelting, which may lead to a paradigm shift in the understanding of ice physics and chemistry.

Intrinsic reconstruction of ice-I surfaces

Understanding the precise atomic structure of ice surfaces is critical for revealing the mechanisms of physical and chemical phenomena at the surfaces, such as ice growth, melting, and chemical reactions. Nevertheless, no conclusive structure has been established. In this study, noncontact atomic force microscopy was used to address the characterization of the atomic structures of ice Ih(0001) and Ic(111) surfaces. The topmost hydrogen atoms are arranged with a short-range (2 × 2) order, independent of the ice thickness and growth substrates used. The electrostatic repulsion between non-hydrogen-bonded water molecules at the surface causes a reduction in the number of the topmost hydrogen atoms together with a distortion of the ideal honeycomb arrangement of water molecules, leading to a short-range-ordered surface reconstruction.

Molecular Structure and Modeling of Water-Air and Ice-Air Interfaces Monitored by Sum-Frequency Generation

From a glass of water to glaciers in Antarctica, water-air and ice-air interfaces are abundant on Earth. Molecular-level structure and dynamics at these interfaces are key for understanding many chemical/physical/atmospheric processes including the slipperiness of ice surfaces, the surface tension of water, and evaporation/sublimation of water. Sum-frequency generation (SFG) spectroscopy is a powerful tool to probe the molecular-level structure of these interfaces because SFG can specifically probe the topmost interfacial water molecules separately from the bulk and is sensitive to molecular conformation. Nevertheless, experimental SFG has several limitations. For example, SFG cannot provide information on the depth of the interface and how the orientation of the molecules varies with distance from the surface. By combining the SFG spectroscopy with simulation techniques, one can directly compare the experimental data with the simulated SFG spectra, allowing us to unveil the molecular-level structure of water-air and ice-air interfaces. Here, we present an overview of the different simulation protocols available for SFG spectra calculations. We systematically compare the SFG spectra computed with different approaches, revealing the advantages and disadvantages of the different methods. Furthermore, we account for the findings through combined SFG experiments and simulations and provide future challenges for SFG experiments and simulations at different aqueous interfaces.

Desorption of neutrals, cations, and anions from core-excited amorphous solid water

Core-excitation of water ice releases many different molecules and ions in the gas phase. Studying these desorbed species and the underlying desorption mechanisms can provide useful information on the effects of x-ray irradiation in ice. We report a detailed study of the x-ray induced desorption of a number of neutral, cationic, and anionic species from amorphous solid water. We discuss the desorption mechanisms and the relative contributions of Auger and secondary electrons (x-ray induced electron stimulated desorption) and initial excitation (direct desorption) as well as the role of photochemistry. Anions are shown to desorb not just through processes linked with secondary electrons but also through direct dissociation of the core-excited molecule. The desorption spectra of oxygen ions (O⁺, OH⁺, H₂O⁺, O^-, and OH^-) give a new perspective on their previously reported very low desorption yields for most types of irradiations of water, showing that they mostly originate from the dissociation of photoproducts such as H₂O₂.

Ionization of Nitric Acid on Crystalline Ice: The Role of Defects and Collective Proton Movement

Ionization of nitric acid (HNO3) on a model ice surface is studied using ab initio molecular dynamics at temperatures of 200 and 40 K with a surface slab model that consists of the ideal ice basal plane with locally optimized and annealed defects. Pico- and subpicosecond ionization of nitric acid can be achieved in the defect sites. Key features of the rapid ionization are (a) the efficient solvation of the polyatomic nitrate anion, by stealing hydrogen bonds from the weakened hydrogen bonds at defect sites, (b) formation of contact ion pairs to stable "presolvated" molecular species that are present at the defects,

Some fundamental properties and reactions of ice surfaces at low temperatures

Ice surfaces offer a unique chemical environment in which reactions occur quite differently from those in liquid water or gas phases. In this article, we examine the basic properties of ice surfaces below the surface premelting temperature and discuss some of the recent investigations carried out on reactions at the ice surfaces. The static and dynamic properties of an ice surface as a reaction medium, such as its structure, molecule diffusion and proton transfer dynamics, and the surface preference of hydronium and hydroxide ions, are discussed in relation to the reactivity of the surface.

SpecSwap-RMC: a novel reverse Monte Carlo approach using a discrete set of local configurations and pre-computed properties

We present a novel approach to reverse Monte Carlo (RMC) modeling, SpecSwap-RMC, specifically applicable to structure modeling based on properties that require significant computer time to evaluate. In this approach pre-computed property data from a discrete set of local configurations are used and the configuration space is expressed in this basis. Atomistic moves are replaced with swap moves of contributions to a sample set representing the state of the simulated system. We demonstrate the approach by fitting jointly and separately the EXAFS signal and x-ray absorption spectrum (XAS) of ice Ih using a SpecSwap sample set of 80 configurations from a library of 1382 local structures with associated pre-computed spectra. As an additional demonstration we compare SpecSwap and FEFFIT fits of EXAFS data on crystalline copper, finding excellent agreement. SpecSwap-RMC thus extends RMC structure modeling to any property that can be computed from a structure irrespective of computational expense, but at the cost of a reduced configuration space. The method is general enough that it can be applied to any sets of computed properties, not necessarily limited to structure determination.

The structure of mixed H2O-OH monolayer films on Ru(0001)

Scanning tunneling microscopy (STM) and x-ray absorption spectroscopy (XAS) have been used to study the structures produced by water on Ru(0001) at temperatures above 140 K. It was found that while undissociated water layers are metastable below 140 K, heating above this temperature produces drastic transformations, whereby a fraction of the water molecules partially dissociate and form mixed H(2)O-OH structures. X-ray photoelectron spectroscopy and XAS revealed the presence of hydroxyl groups with their O-H bond essentially parallel to the surface. STM images show that the mixed H(2)O-OH structures consist of long narrow stripes aligned with the three crystallographic directions perpendicular to the close-packed atomic rows of the Ru(0001) substrate. The internal structure of the stripes is a honeycomb network of H-bonded water and hydroxyl species. We found that the metastable low temperature molecular phase can also be converted to a mixed H(2)O-OH phase through excitation by the tunneling electrons when their energy is 0.5 eV or higher above the Fermi level. Structural models based on the STM images were used for density functional theory optimizations of the stripe geometry. The optimized geometry was then utilized to calculate STM images for comparison with the experiment.

Electron damage to supported ice investigated by scanning tunneling microscopy and spectroscopy

We study the interaction of low-energy electrons with crystalline ice (D2O) on Cu(111) by low-temperature scanning tunneling microscopy and spectroscopy. Electrons induce dissociation of the molecules with a threshold of approximately 3 eV. The large dissociation yield of the order of 10(-8)/electron and the extended area of dissociation are attributed to a shift in conduction band during the dissociation. Voltage dependent differences in imaging of ice and dissociated ice are reflected in the spectroscopic signature.

The local electronic structure of alpha-Li3N

New theoretical and experimental investigations of the occupied and unoccupied local electronic densities of states (DOS) are reported for alpha-Li(3)N. Band-structure and density-functional theory calculations confirm the absence of covalent bonding character. However, real-space full-multiple-scattering (RSFMS) calculations of the occupied local DOS find less extreme nominal valences than have previously been proposed. Nonresonant inelastic x-ray scattering, RSFMS calculations, and calculations based on the Bethe-Salpeter equation are used to characterize the unoccupied electronic final states local to both the Li and N sites. There is a good agreement between experiment and theory. Throughout the Li 1s near-edge region, both experiment and theory find strong similarities in the s-and p-type components of the unoccupied local final DOS projected onto an orbital angular momentum basis (l-DOS). An unexpected, significant correspondence exists between the near-edge spectra for the Li 1s and N 1s initial states. We argue that both spectra are sampling essentially the same final DOS due to the combination of long core-hole lifetimes, long photoelectron lifetimes, and the fact that orbital angular momentum is the same for all relevant initial states. Such considerations may be generally applicable for low atomic number compounds.

Probing the electron delocalization in liquid water and ice at attosecond time scales

We determine electron delocalization rates in liquid water and ice using core-hole decay spectroscopy. The hydrogen-bonded network delocalizes the electrons in less than 500 as. Broken or weak hydrogen bonds--in the liquid or at the surface of ice--provide states where the electron remains localized longer than 20 fs. These asymmetrically bonded water species provide electron traps, acting as a strong precursor channel to the hydrated electron.

Oxygen K-edge X-ray absorption fine structure studies of vacuum-deposited ice films

Multilayer ice films deposited on polycrystalline Au(111) substrates at 90 and 153 K under ultra-high-vacuum conditions are studied using O K-edge X-ray absorption fine structure spectroscopy in the surface-sensitive partial electron yield mode. Both near-edge and extended oscillatory fine structures are analyzed in combination with the theoretical real-space full-multiple-scattering simulations based on the FEFF8 code. The experimental data consistently indicate that the local structure of the near-surface regions in ice films at both substrate temperatures resembles that of high-density crystalline modifications of ice (e.g., ice II, ice III, or high-density amorphous ice). In addition, the ice films deposited at 153 K most probably contain a minor fraction of low-density cubic ice (I(c) phase).

Surface trapped excess electrons on ice

Local trapping of excess electrons at the surface of solid water systems has recently been observed in large water clusters and at the ice/vacuum interface. The existence of stable surface-bound states for the excess electron may have important implications in atmospheric chemistry, electrochemistry, and radiation physics. By means of first-principles molecular dynamics we find that excess electrons induce a structural reconstruction of the ice surface on a time scale of a fraction of a picosecond. The surface molecular rearrangement leads to an increase of the number of dangling OH bonds pointing towards the vacuum and to the appearance of an electrostatic barrier preventing the penetration of the electron in the bulk. Both factors imply a remarkable stability for the surface-bound excess electron, with respect to its decay into the bulk solvated state.

Crystalline ice growth on PT(111): observation of a hydrophobic water monolayer

The growth of crystalline water films on Pt(111) is investigated using rare gas physisorption. The water monolayer wets Pt(111) at all temperatures investigated (20-155 K). At low temperatures (T< or =120 K), additional water layers kinetically wet the monolayer surface. However, crystalline ice films grown at higher temperatures (T > 135 K) do not wet the water monolayer. These results are consistent with recent theory and experiments suggesting that the molecules in the water monolayer form a surface with no dangling OH bonds or lone pair electrons, giving rise to a hydrophobic water monolayer on Pt(111).

The electronic structure of liquid water within density-functional theory

In the last decade, computational studies of liquid water have mostly concentrated on ground-state properties. However, recent spectroscopic measurements have been used to infer the structure of water, and the interpretation of optical and x-ray spectra requires accurate theoretical models of excited electronic states, not only of the ground state. To this end, we investigate the electronic properties of water at ambient conditions using ab initio density-functional theory within the generalized gradient approximation (DFT/GGA), focusing on the unoccupied subspace of Kohn-Sham eigenstates. We generate long (250 ps) classical trajectories for large supercells, up to 256 molecules, from which uncorrelated configurations of water molecules are extracted for use in DFT/GGA calculations of the electronic structure. We find that the density of occupied states of this molecular liquid is well described with 32-molecule supercells using a single k point (k=0) to approximate integration over the first Brillouin zone. However, the description of the unoccupied electronic density of states (u-EDOS) is sensitive to finite size effects. Small, 32-molecule supercell calculations, using the Gamma-point approximation, yield a spuriously isolated state above the Fermi level. Nevertheless, the more accurate u-EDOS of large, 256-molecule supercells may be reproduced using smaller supercells and increased k-point sampling. This indicates that the electronic structure of molecular liquids such as water is relatively insensitive to the long-range disorder in the molecular structure. These results have important implications for efficiently increasing the accuracy of spectral calculations for water and other molecular liquids.

X-ray Absorption Spectroscopy Measurements of Liquid Water

Recent studies, based on X-ray absorption spectroscopy (XAS) and X-ray Raman scattering (XRS), have shown that the hydrogen bond network in liquid water consists mainly of water molecules with only two strong hydrogen bonds. Since this result is controversial, it is important to demonstrate the reliability of the experimental data, which is the purpose of this paper. Here we compare X-ray absorption spectra of liquid water recorded with five very different techniques sensitive to the local environment of the absorbing molecule. Overall, the spectra obtained with photon detection show a very close similarity and even the observable minor differences can be understood. The comparison demonstrates that XAS and XRS can indeed be applied reliably to study the local bonding of the water molecule and thus to reveal the hydrogen bond situation in bulk water.

X-ray Absorption Spectroscopy of Liquid Methanol Microjets: Bulk Electronic Structure and Hydrogen Bonding Network

We have measured the X-ray absorption (XA) spectrum of liquid (298 K) methanol at the oxygen and carbon K edges. The 4a(1) orbital at the O K edge exhibits a pronounced sensitivity to the formation of intermolecular hydrogen bonds, with significant differences observed between the vapor and bulk spectra, whereas the C K edge reveals only subtle corresponding spectral changes. Comparison with DFT computed spectra of model methanol clusters indicates that the bulk liquid comprises long chains (n > 6) and rings of hydrogen-bonded monomers.

Water Adsorption, Desorption, and Clustering on FeO(111)

The adsorption of water on FeO(111) is investigated using temperature programmed desorption (TPD) and infrared reflection absorption spectroscopy (IRAS). Well-ordered 2 ML thick FeO(111) films are grown epitaxially on a Pt(111) substrate. Water adsorbs molecularly on FeO(111) and desorbs with a well resolved monolayer peak. IRAS measurements as a function of coverage are performed for water deposited at 30 and 135 K. For all coverages (0.2 ML and greater), the adsorbed water exhibits significant hydrogen bonding. Differences in IRAS spectra for water adsorbed at 30 and 135 K are subtle but suggest that water adsorbed at 135 K is well ordered. Monolayer nitrogen TPD spectra from water covered FeO(111) surfaces are used to investigate the clustering of the water as a function of deposition or annealing temperature. Temperature dependent water overlayer structures result from differences in water diffusion rates on bare FeO(111) and on water adsorbed on FeO(111). Features in the nitrogen TPD spectra allow the monolayer wetting and 2-dimensional (2D) ordering of water on FeO(111) to be followed. Voids in a partially disordered first water layer exist for water deposited below 120 K and ordered 2D islands are found when depositing water above 120 K.

The hydrogen bond in ice probed by soft x-ray spectroscopy and density functional theory

We combine photoelectron and x-ray absorption spectroscopy with density functional theory to derive a molecular orbital picture of the hydrogen bond in ice. We find that the hydrogen bond involves donation and back-donation of charge between the oxygen lone pair and the O-H antibonding orbitals on neighboring molecules. Together with internal s-p rehybridization this minimizes the repulsive charge overlap of the connecting oxygen and hydrogen atoms, which is essential for a strong attractive electrostatic interaction. Our joint experimental and theoretical results demonstrate that an electrostatic model based on only charge induction from the surrounding medium fails to properly describe the internal charge redistributions upon hydrogen bonding.

Half or full core hole in density functional theory X-ray absorption spectrum calculations of water?

We analyze the performance of two different core-hole potentials in the theoretical modeling of XAS of ice, liquid and gas phase water; the use of a full core-hole (FCH) in the calculations, as suggested by Hetenyi et al. [B. Hetenyi, F. De Angelis, P. Giamozzi and R. Car, J. Chem. Phys., 2004, 120(18), 8632], gives poor agreement with experiment in terms of intensity distribution as well as transition energies, while the half core hole (HCH) potential, in the case of water, provides a better compromise between initial and final state effects, leading to good agreement with the experimental data.

Ultrafast molecular dissociation of water in ice

Using x-ray emission and photoemission spectroscopies to measure the occupied valence levels in a thin crystalline ice film, we resolve the ionization-induced dissociation of water in ice on a femtosecond time scale. Isotope substitution confirms proton transfer during the core-hole lifetime in spite of the nonresonant excitation. Through ab initio molecular dynamics on the core-ionized state, the dissociation and spectrum evolution are followed at femtosecond intervals. The theoretical simulations confirm the experimental analysis and allow for a detailed study of the dissociative reaction path.

Water? What's so special about it?

What is so special about water? Why does it have the properties it has, and how might these reasons be relevant to its apparent biological importance? By exploring the structure and dynamics of water, from the isolated molecule and its interactions, through its many crystalline phases and to its so-called anomalous liquid phase, some of its apparently unusual behaviour is rationalized. The way in which it interacts with some relatively simple interfaces is also discussed. As a result of this exploration, a checklist of possible molecular-level reasons for its biological importance is devised.