Diffraction and IR/Raman data do not prove tetrahedral water

Temperature Dependence of Hydrogen Bond Networks of Liquid Water: Thermodynamic Properties and Structural Heterogeneity from Topological Descriptors

Topological analyses of hydrogen bond networks were performed based on the complex network and island statistics of liquid water at different temperatures. The influence of temperature on the liquid water structures and the topological properties of the hydrogen bond networks was investigated by Metropolis Monte Carlo simulations with the TIP4P/2005 potential model. The bilinear behavior of the second peak in the radial distribution function with the temperature was properly reproduced by these simulations. The average connectivity also displayed a bilinear behavior consistent with being a local descriptor. The semiglobal average path length (or geodesic distance) descriptor showed an unprecedented trimodal distribution, whose areas were dependent on the temperature. Considering equilibrium between these three sets of networks, standard enthalpy and entropy of equilibrium were determined for the first time, providing new insights into the structural heterogeneities of liquid water with interesting perspectives for modeling these quantitative properties of hydrogen bond networks.

Polarizable Multipolar Molecular Dynamics Using Distributed Point Charges

Distributed point charge models (DCM) and their minimal variants (MDCM) have been integrated with tools widely used for condensed-phase simulations, including a virial-based barostat and a slow-growth algorithm for thermodynamic integration. Minimal DCM is further developed in a systematic fashion to reduce fitting errors in the electrostatic interaction energy, and a new fragment-based approach offers considerable speedup of the MDCM fitting process for larger molecules with increased numbers of off-centered charged sites. Finally, polarizable (M)DCM is also introduced in the present work. The developments are used in condensed-phase simulations of popular force fields with commonly applied simulation conditions. (M)DCM equivalents for a range of widely used water force fields and for fluorobenzene (PhF) are developed and applied along with the original models to evaluate the impact of reformulating the electrostatic term. Comparisons of the molecular electrostatic potential (MEP), electrostatic interaction energies, and bulk properties from molecular dynamics simulations for a range of models from simple TIPnP (n = 3-5) to the polarizable, multipolar iAMOEBA models for water and an existing quadrupolar model for PhF confirm that DCMs retain the accuracy of the original models, providing a homogeneous, efficient, and generic point charge alternative to a multipolar electrostatic model for force field development and multilevel simulations.

Liquid-phase mega-electron-volt ultrafast electron diffraction

The conversion of light into usable chemical and mechanical energy is pivotal to several biological and chemical processes, many of which occur in solution. To understand the structure-function relationships mediating these processes, a technique with high spatial and temporal resolutions is required. Here, we report on the design and commissioning of a liquid-phase mega-electron-volt (MeV) ultrafast electron diffraction instrument for the study of structural dynamics in solution. Limitations posed by the shallow penetration depth of electrons and the resulting information loss due to multiple scattering and the technical challenge of delivering liquids to vacuum were overcome through the use of MeV electrons and a gas-accelerated thin liquid sheet jet. To demonstrate the capabilities of this instrument, the structure of water and its network were resolved up to the 3 rd hydration shell with a spatial resolution of 0.6 Å; preliminary time-resolved experiments demonstrated a temporal resolution of 200 fs.

Measurements of ultrafast dissociation in resonant inelastic x-ray scattering of water

There has been a discussion on the interpretation of the resonant inelastic x-ray scattering (RIXS) spectra of liquid water in terms of either different structural environments or that core hole dynamics can generate well-resolved dissociative spectral components. We have used RIXS with high resolution in the OH stretch vibration energy part, at extremely high overtones going toward the continuum of full OH bond breakage, to identify the amount of dissociative contributions in the valence band RIXS spectra at different excitation energies. We observe that at low excitation energies, corresponding to population of states with strongly antibonding character, the valence band RIXS spectra have a large contribution from a well-resolved dissociative feature. Instead, at higher excitations, this spectral component diminishes and becomes a weak structure on the high-energy side of one of the spectral peaks related to the 1b₁ state from tetrahedral configurations. This result brings both interpretations to be essential for the understanding of RIXS spectra of liquid water.

Seeing real-space dynamics of liquid water through inelastic x-ray scattering

Water is ubiquitous on earth, but we know little about the real-space motion of molecules in liquid water. We demonstrate that high-resolution inelastic x-ray scattering measurement over a wide range of momentum and energy transfer makes it possible to probe real-space, real-time dynamics of water molecules through the so-called Van Hove function. Water molecules are found to be strongly correlated in space and time with coupling between the first and second nearest-neighbor molecules. The local dynamic correlation of molecules observed here is crucial to a fundamental understanding of the origin of the physical properties of water, including viscosity. The results also suggest that the quantum-mechanical nature of hydrogen bonds could influence its dynamics. The approach used here offers a powerful experimental method for investigating real-space dynamics of liquids.

Probing the OH Stretch in Different Local Environments in Liquid Water

We use resonant inelastic X-ray scattering (RIXS) to resolve vibrational losses corresponding to the OH stretch where the X-ray absorption process allows us to selectively probe different structural subensembles in liquid water. The results point to a unified interpretation of X-ray and vibrational spectroscopic data in line with a picture of two classes of structural environments in the liquid at ambient conditions with predominantly close-packed high-density liquid (HDL) and occasional local fluctuations into strongly tetrahedral low-density liquid (LDL).

Probing the triplet correlation function in liquid water by experiments and molecular simulations

Three-body information of liquid water is extracted using X-ray diffraction experiment as well as in molecular simulations via isothermal pressure derivative of structure factor term.

X-ray and Electron Spectroscopy of Water

Here we present an overview of recent developments of X-ray and electron spectroscopy to probe water at different temperatures. Photon-induced ionization followed by detection of electrons from either the O 1s level or the valence band is the basis of photoelectron spectroscopy. Excitation between the O 1s and the unoccupied states or occupied states is utilized in X-ray absorption and X-ray emission spectroscopies. These techniques probe the electronic structure of the liquid phase and show sensitivity to the local hydrogen-bonding structure. Both experimental aspects related to the measurements and theoretical simulations to assist in the interpretation are discussed in detail. Different model systems are presented such as the different bulk phases of ice and various adsorbed monolayer structures on metal surfaces.

The structural origin of anomalous properties of liquid water

Water is unique in its number of unusual, often called anomalous, properties. When hot it is a normal simple liquid; however, close to ambient temperatures properties, such as the compressibility, begin to deviate and do so increasingly on further cooling. Clearly, these emerging properties are connected to its ability to form up to four well-defined hydrogen bonds allowing for different local structural arrangements. A wealth of new data from various experiments and simulations has recently become available. When taken together they point to a heterogeneous picture with fluctuations between two classes of local structural environments developing on temperature-dependent length scales.

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.

Interaction between bound water molecules and local protein structures: A statistical analysis of the hydrogen bond structures around bound water molecules

Water molecules play an important role in protein folding and protein interactions through their structural association with proteins. Examples of such structural association can be found in protein crystal structures, and can often explain protein functionality in the context of structure. We herein report the systematic analysis of the local structures of proteins interacting with water molecules, and the characterization of their geometric features. We first examined the interaction of water molecules with a large local interaction environment by comparing the preference of water molecules in three regions, namely, the protein-protein interaction (PPI) interfaces, the crystal contact (CC) interfaces, and the non-interfacial regions. High preference of water molecules to the PPI and CC interfaces was found. In addition, the bound water on the PPI interface was more favorably associated with the complex interaction structure, implying that such water-mediated structures may participate in the shaping of the PPI interface. The pairwise water-mediated interaction was then investigated, and the water-mediated residue-residue interaction potential was derived. Subsequently, the types of polar atoms surrounding the water molecules were analyzed, and the preference of the hydrogen bond acceptor was observed. Furthermore, the geometries of the structures interacting with water were analyzed, and it was found that the major structure on the protein surface exhibited planar geometry rather than tetrahedral geometry. Several previously undiscovered characteristics of water-protein interactions were unfolded in this study, and are expected to lead to a better understanding of protein structure and function.

Potential Paths for the Hydrogen-Bond Relaxing With (H2O)N Cluster Size

Relaxation of the inter- and intra-molecular interactions for the hydrogen bond (O:H-O) between undercoordinated molecules determines the unusual behavior of water nanodroplets and nanobubbles. However, probing such potentials remains unreality. Here we show that the Lagrangian solution [Huang et al., J. Phys. Chem. B, 2013. 117: 13639] transforms the observed H-O bond (x = H) and O:H nonbond (x = L) lengths and phonon frequencies (dx, x) [Sun et al., J. Phys. Chem. Lett., 2013. 4: 2565] into the respective force constants and bond energies (kx, Ex) and hence enables the mapping of the potential paths for the O:H-O bond relaxing with water cluster size. Results show that molecular undercoordination not only reduces the molecular size (dH) with enhanced H-O energy from the bulk value of 3.97 to 5.10 eV for a H2O monomer, but also enlarges the molecular separation (dL) with reduced O:H energy from 95 to 35 meV for a dimer. The H-O energy gain raises the melting point from bulk value 273 to 310 K for the skin and the O:H energy loss lowers the freezing temperature from bulk value 258 to 202 K for 1.4 nm sized droplet, by dispersing the quasisolid phase boundaries.

Triplet correlation functions in liquid water

Triplet correlations have been shown to play a crucial role in the transformation of simple liquids to anomalous tetrahedral fluids [M. Singh, D. Dhabal, A. H. Nguyen, V. Molinero, and C. Chakravarty, Phys. Rev. Lett. 112, 147801 (2014)]. Here we examine triplet correlation functions for water, arguably the most important tetrahedral liquid, under ambient conditions, using configurational ensembles derived from molecular dynamics (MD) simulations and reverse Monte Carlo (RMC) datasets fitted to experimental scattering data. Four different RMC data sets with widely varying hydrogen-bond topologies fitted to neutron and x-ray scattering data are considered [K. T. Wikfeldt, M. Leetmaa, M. P. Ljungberg, A. Nilsson, and L. G. M. Pettersson, J. Phys. Chem. B 113, 6246 (2009)]. Molecular dynamics simulations are performed for two rigid-body effective pair potentials (SPC/E and TIP4P/2005) and the monatomic water (mW) model. Triplet correlation functions are compared with other structural measures for tetrahedrality, such as the O-O-O angular distribution function and the local tetrahedral order distributions. In contrast to the pair correlation functions, which are identical for all the RMC ensembles, the O-O-O triplet correlation function can discriminate between ensembles with different degrees of tetrahedral network formation with the maximally symmetric, tetrahedral SYM dataset displaying distinct signatures of tetrahedrality similar to those obtained from atomistic simulations of the SPC/E model. Triplet correlations from the RMC datasets conform closely to the Kirkwood superposition approximation, while those from MD simulations show deviations within the first two neighbour shells. The possibilities for experimental estimation of triplet correlations of water and other tetrahedral liquids are discussed.

Comparison of x-ray absorption spectra between water and ice: new ice data with low pre-edge absorption cross-section

The effect of crystal growth conditions on the O K-edge x-ray absorption spectra of ice is investigated through detailed analysis of the spectral features. The amount of ice defects is found to be minimized on hydrophobic surfaces, such as BaF2(111), with low concentration of nucleation centers. This is manifested through a reduction of the absorption cross-section at 535 eV, which is associated with distorted hydrogen bonds. Furthermore, a connection is made between the observed increase in spectral intensity between 544 and 548 eV and high-symmetry points in the electronic band structure, suggesting a more extended hydrogen-bond network as compared to ices prepared differently. The spectral differences for various ice preparations are compared to the temperature dependence of spectra of liquid water upon supercooling. A double-peak feature in the absorption cross-section between 540 and 543 eV is identified as a characteristic of the crystalline phase. The connection to the interpretation of the liquid phase O K-edge x-ray absorption spectrum is extensively discussed.

Ion-Solvation-Induced Molecular Reorganization in Liquid Water Probed by Resonant Inelastic Soft X-ray Scattering

The molecular structure of liquid water is susceptible to changes upon admixture of salts due to ionic solvation, which provides the basis of many chemical and biochemical processes. Here we demonstrate how the local electronic structure of aqueous potassium chloride (KCl) solutions can be studied by resonant inelastic soft X-ray scattering (RIXS) to monitor the effects of the ion solvation on the hydrogen-bond (HB) network of liquid water. Significant changes in the oxygen K-edge emission spectra are observed with increasing KCl concentration. These changes can be attributed to modifications in the proton dynamics, caused by a specific coordination structure around the salt ions. Analysis of the spectator decay spectra reveals a spectral signature that could be characteristic of this structure.

A continuous mixture of two different dimers in liquid water

Liquid water is formed by a continuous mixture of two different dimers (cis and trans) with distinct energies related to different relative water molecule orientations.

The Radial Distribution Functions of Water as Derived from Radiation Total Scattering Experiments: Is There Anything We Can Say for Sure?

The present paper reviews the investigation of ambient water structure and focusses in particular on the determination of the radial distribution functions of water from total experimental radiation scattering experiments. A novel method for removing the inelastic scattering from neutron data is introduced, and the effect of Compton scattering on X-ray data is discussed. In addition the extent to which quantum effects can be discerned between heavy and light water is analysed against these more recent data. It is concluded that, with the help of modern data analysis and computer simulation tools to interrogate the scattering data, a considerable degree of consistency can be obtained between recent and past scattering experiments on water. That consistency also gives a realistic estimate of the likely uncertainties in the extracted radial distribution functions, as well as offering a benchmark against which future experiments can be judged.

Wide-angle X-ray diffraction and molecular dynamics study of medium-range order in ambient and hot water

We have developed wide-angle X-ray diffraction measurements with high energy-resolution and accuracy to study water structure at three different temperatures (7, 25 and 66 °C) under normal pressure. Using a spherically curved Ge crystal an energy resolution better than 15 eV has been achieved which eliminates influence from Compton scattering. The high quality of the data allows for a reliable Fourier transform of the experimental data resolving shell structure out to ~12 Å, i.e. 5 hydration shells. Large-scale molecular dynamics (MD) simulations using the TIP4P/2005 force-field reproduce excellently the experimental shell-structure in the range 4-12 Å although less agreement is seen for the first peak in the intermolecular pair-correlation function (PCF). The Shiratani-Sasai Local Structure Index [J. Chem. Phys. 104, 7671 (1996)] identifies a tetrahedral minority giving the intermediate-range oscillations in the O-O PCF and a disordered majority providing a more featureless background in this range. The current study supports the proposal that the structure of liquid water, even at high temperatures, can be described in terms of a two-state fluctuation model involving local structures related to the high-density and low-density forms of liquid water postulated in the liquid-liquid phase transition hypothesis.

Hydrogen bond dynamics in heavy water studied with quantum dynamical simulations

The structure and dynamics of the hydrogen-bond network in heavy water (D(2)O) is studied as a function of the temperature using quantum dynamical simulations. Our approach combines an ab initio-based representation of the water interactions with an explicit quantum treatment of the molecular motion. A direct connection between the calculated linear and nonlinear vibrational spectra and the underlying molecular dynamics is made, which provides new insights into the rearrangement of the hydrogen-bond network in heavy water. A comparison with previous calculations on liquid H(2)O suggests that tunneling does not effectively contribute to the dynamics of the water hydrogen-bond network above the melting point. However, the effects of nuclear quantization are not negligible at all temperatures and become increasingly important near the melting point, which is in agreement with recent experimental analysis of the structural properties of liquid water as well as of the proton momentum distribution in supercooled water.