X-ray absorption spectrum simulations of hexagonal ice.

We calibrate basis sets and performance of two theoretical approaches to compute X-ray absorption spectra (XAS) of condensed water by comparison to experiments on hexagonal ice Ih. We apply both the transition-potential half-core-hole approach and the complex polarization propagator using four different models of the crystal with increasing oxygen and proton disorder but find poor agreement with experiments. We note that there are large variations in experimental spectra depending on detection mode and how the ice samples were prepared, which leads us to critically investigate what structures were actually prepared and measured in each case. This is done by using a Monte Carlo-based fitting technique which fits the spectra based on a library of precomputed spectra and assigns weights to contributions from different model structures. These are then used to generate O-O and O-H radial distribution functions and tetrahedrality parameters associated with each of the measured spectra. We find that all spectra are associated with sharp peaks at the oxygen positions in the perfect lattice, but with significant disorder around these positions. We suggest that presently available XAS of hexagonal ice are not fully representative of the perfect crystalline lattice, but contain varying amounts of defects and possible contributions from low-density amorphous ice.

A proposal for the structure of high- and low-density fluctuations in liquid water.

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 ∼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.

Relationship between x-ray emission and absorption spectroscopy and the local H-bond environment in water.

The connection between specific features in the water X-ray absorption spectrum and X-ray emission spectrum (XES) and the local H-bond coordination is studied based on structures obtained from path-integral molecular dynamics simulations using either the opt-PBE-vdW density functional or the MB-pol force field. Computing the XES spectrum using all molecules in a snapshot results in only one peak in the lone-pair (1b1) region, while the experiment shows two peaks separated by 0.8-0.9 eV. Different H-bond configurations were classified based on the local structure index (LSI) and a geometrical H-bond cone criterion. We find that tetrahedrally coordinated molecules characterized by high LSI values and two strong donated and two strong accepted H-bonds contribute to the low energy 1b1 emission peak and to the post-edge region in absorption. Molecules with the asymmetric H-bond environment with one strong accepted H-bond and one strong donated H-bond and low LSI values give rise to the high energy 1b1 peak in the emission spectrum and mainly contribute to the pre-edge and main-edge in the absorption spectrum. The 1b1 peak splitting can be increased to 0.62 eV by imposing constraints on the H-bond length, i.e., for very tetrahedral structures short H-bonds (less than 2.68 Å) and for very asymmetric structures elongated H-bonds (longer than 2.8 Å). Such structures are present, but underrepresented, in the simulations which give more of an average of the two extremes.

Requirements of first-principles calculations of X-ray absorption spectra of liquid water.

A computational benchmark study on X-ray absorption spectra of water has been performed by means of transition-potential density functional theory (TP-DFT), damped time-dependent density functional theory (TDDFT), and damped coupled cluster (CC) linear response theory. For liquid water, using TDDFT with a tailored CAM-B3LYP functional and a polarizable embedding, we find that an embedding with over 2000 water molecules is required to fully converge spectral features for individual molecules, but a substantially smaller embedding can be used within averaging schemes. TP-DFT and TDDFT calculations on 100 MD structures demonstrate that TDDFT produces a spectrum with spectral features in good agreement with experiment, while it is more difficult to fully resolve the spectral features in the TP-DFT spectrum. Similar trends were also observed for calculations of bulk ice. In order to further establish the performance of these methods, small water clusters have been considered also at the CC2 and CCSD levels of theory. Issues regarding the basis set requirements for spectrum simulations of liquid water and the determination of gas-phase ionization potentials are also discussed.

X-ray emission spectroscopy of bulk liquid water in "no-man's land".

The structure of bulk liquid water was recently probed by x-ray scattering below the temperature limit of homogeneous nucleation (TH) of ∼232 K [J. A. Sellberg et al., Nature 510, 381-384 (2014)]. Here, we utilize a similar approach to study the structure of bulk liquid water below TH using oxygen K-edge x-ray emission spectroscopy (XES). Based on previous XES experiments [T. Tokushima et al., Chem. Phys. Lett. 460, 387-400 (2008)] at higher temperatures, we expected the ratio of the 1b1' and 1b1″ peaks associated with the lone-pair orbital in water to change strongly upon deep supercooling as the coordination of the hydrogen (H-) bonds becomes tetrahedral. In contrast, we observed only minor changes in the lone-pair spectral region, challenging an interpretation in terms of two interconverting species. A number of alternative hypotheses to explain the results are put forward and discussed. Although the spectra can be explained by various contributions from these hypotheses, we here emphasize the interpretation that the line shape of each component changes dramatically when approaching lower temperatures, where, in particular, the peak assigned to the proposed disordered component would become more symmetrical as vibrational interference becomes more important.