Radial distribution functions of water: Models vs experiments.

We study the temperature behavior of the first four peaks of the oxygen-oxygen radial distribution function of water, simulated by the TIP4P/2005, MB-pol, TIP5P, and SPC/E models and compare to experimental X-ray diffraction data, including a new measurement which extends down to 235 K [H. Pathak et al., J. Chem. Phys. 150, 224506 (2019)]. We find the overall best agreement using the MB-pol and TIP4P/2005 models. We observe, upon cooling, a minimum in the position of the second shell simulated with TIP4P/2005 and SPC/E potentials, located close to the temperature of maximum density. We also calculated the two-body entropy and the contributions coming from the first, second, and outer shells to this quantity. We show that, even if the main contribution comes from the first shell, the contribution of the second shell can become important at low temperature. While real water appears to be less ordered at short distance than obtained by any of the potentials, the different water potentials show more or less order compared to the experiments depending on the considered length-scale.

Translational and rotational dynamics of high and low density TIP4P/2005 water.

We use molecular dynamics simulations using TIP4P/2005 to investigate the self- and distinct-van Hove functions for different local environments of water, classified using the local structure index as an order parameter. The orientational dynamics were studied through the calculation of the time-correlation functions of different-order Legendre polynomials in the OH-bond unit vector. We found that the translational and orientational dynamics are slower for molecules in a low-density local environment and correspondingly the mobility is enhanced upon increasing the local density, consistent with some previous works, but opposite to a recent study on the van Hove function. From the analysis of the distinct dynamics, we find that the second and fourth peaks of the radial distribution function, previously identified as low density-like arrangements, show long persistence in time. The analysis of the time-dependent interparticle distance between the central molecule and the first coordination shell shows that particle identity persists longer than distinct van Hove correlations. The motion of two first-nearest-neighbor molecules thus remains coupled even when this correlation function has been completely decayed. With respect to the orientational dynamics, we show that correlation functions of molecules in a low-density environment decay exponentially, while molecules in a local high-density environment exhibit bi-exponential decay, indicating that dynamic heterogeneity of water is associated with the heterogeneity among high-density and between high-density and low-density species. This bi-exponential behavior is associated with the existence of interstitial waters and the collapse of the second coordination sphere in high-density arrangements, but not with H-bond strength.

Coherent X-rays reveal the influence of cage effects on ultrafast water dynamics.

The dynamics of liquid water feature a variety of time scales, ranging from extremely fast ballistic-like thermal motion, to slower molecular diffusion and hydrogen-bond rearrangements. Here, we utilize coherent X-ray pulses to investigate the sub-100 fs equilibrium dynamics of water from ambient conditions down to supercooled temperatures. This novel approach utilizes the inherent capability of X-ray speckle visibility spectroscopy to measure equilibrium intermolecular dynamics with lengthscale selectivity, by measuring oxygen motion in momentum space. The observed decay of the speckle contrast at the first diffraction peak, which reflects tetrahedral coordination, is attributed to motion on a molecular scale within the first 120 fs. Through comparison with molecular dynamics simulations, we conclude that the slowing down upon cooling from 328 K down to 253 K is not due to simple thermal ballistic-like motion, but that cage effects play an important role even on timescales over 25 fs due to hydrogen-bonding.

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

Despite very significant developments in scattering experiments like X-ray and neutron diffraction, it has been challenging to elucidate the nature of tetrahedral molecular configurations in liquid water. A key question is whether the pair correlation functions, which can be obtained from scattering experiments, are sufficient to describe the tetrahedral ordering of water molecules. In our previous study (Dhabal et al., J. Chem. Phys., 2014, 141, 174504), using data-sets generated from reverse Monte Carlo and molecular dynamics simulations, we showed that the triplet correlation functions contain important information on the tetrahedrality of water in the liquid state. In the present study, X-ray scattering experiments and molecular dynamics (MD) simulations are used to link the isothermal pressure derivative of the structure factor with the triplet correlation functions for water. Triplet functions are determined for water up to 3.3 kbar at 298 K to display the effect of pressure on the water structure. The results suggest that triplet functions (H[combining tilde](q)) obtained using a rigid-body TIP4P/2005 water model are consistent with the experimental results. The triplet functions obtained in experiment as well as in simulations evince that in the case of tetrahedral liquids, exertion of higher pressure leads to a better agreement with the Kirkwood superposition approximation (KSA). We further validate this observation using the triplet correlation functions (g(3)(r,s,t)) calculated directly from simulation trajectory, revealing that both H[combining tilde](q) in q-space and g(3)(r,s,t) in real-space contain similar information on the tetrahedrality of liquids. This study demonstrates that the structure factor, even though it has only pair correlation information of the liquid structure, can shed light on three-body correlations in liquid water through its isothermal pressure derivative term.

Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions.

Almost 50 years have passed from the first computer simulations of water, and a large number of molecular models have been proposed since then to elucidate the unique behavior of water across different phases. In this article, we review the recent progress in the development of analytical potential energy functions that aim at correctly representing many-body effects. Starting from the many-body expansion of the interaction energy, specific focus is on different classes of potential energy functions built upon a hierarchy of approximations and on their ability to accurately reproduce reference data obtained from state-of-the-art electronic structure calculations and experimental measurements. We show that most recent potential energy functions, which include explicit short-range representations of two-body and three-body effects along with a physically correct description of many-body effects at all distances, predict the properties of water from the gas to the condensed phase with unprecedented accuracy, thus opening the door to the long-sought "universal model" capable of describing the behavior of water under different conditions and in different environments.

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.

Selective probing of the OH or OD stretch vibration in liquid water using resonant inelastic soft-x-ray scattering.

High-resolution O 1s resonant inelastic x-ray scattering spectra of liquid H2O, D2O, and HDO, obtained by excitation near the preedge resonance show, in the elastic line region, well-separated multiple vibrational structures corresponding to the internal OH stretch vibration in the ground state of water. The energy of the first-order vibrational excitation is strongly blueshifted with respect to the main band in the infrared or Raman spectra of water, indicating that water molecules with a highly weakened or broken donating hydrogen bond are correlated with the preedge structure in the x-ray absorption spectrum. The vibrational profile of preedge excited HDO water is well fitted with 50%±20% greater OH-stretch contribution compared to OD, which strongly supports a preference for OH being the weakened or broken H-bond in agreement with the well-known picture that D2O makes stronger H-bonds than H2O. Accompanying path-integral molecular dynamics simulations show that this is particularly the case for strongly asymmetrically H-bonded molecules, i.e., those that are selected by preedge excitation.

In situ X-ray probing reveals fingerprints of surface platinum oxide.

In situ X-ray absorption spectroscopy (XAS) at the Pt L(3) edge is a useful probe for Pt-O interactions at polymer electrolyte membrane fuel cell (PEMFC) cathodes. We show that XAS using the high energy resolution fluorescence detection (HERFD) mode, applied to a well-defined monolayer Pt/Rh(111) sample where the bulk penetrating hard X-rays probe only surface Pt atoms, provides a unique sensitivity to structure and chemical bonding at the Pt-electrolyte interface. Ab initio multiple-scattering calculations using the FEFF code and complementary extended X-ray absorption fine structure (EXAFS) results indicate that the commonly observed large increase of the white-line at high electrochemical potentials on PEMFC cathodes originates from platinum oxide formation, whereas previously proposed chemisorbed oxygen-containing species merely give rise to subtle spectral changes.

Oxygen-oxygen correlations in liquid water: addressing the discrepancy between diffraction and extended x-ray absorption fine-structure using a novel multiple-data set fitting technique.

The first peak of the oxygen-oxygen pair-correlation function (O-O PCF) is a critical measure of the first coordination-shell distances in liquid water. Recently, a discrepancy has been uncovered between diffraction and extended x-ray absorption fine-structure (EXAFS) regarding the height and position of this peak, where EXAFS gives a considerably more well-defined peak at a shorter distance compared to the diffraction results. This discrepancy is here investigated through a novel multiple-data set structure modeling technique, SpecSwap-RMC, based on the reverse Monte Carlo (RMC) method. Fitting simultaneously to both EXAFS and a diffraction-based O-O PCF shows that even though the reported EXAFS results disagree with diffraction, the two techniques can be reconciled by taking into account a strong contribution from the photoelectron scattering focusing effect in EXAFS originating from nearly linear hydrogen bonds. This many-body contribution, which is usually neglected in RMC modeling of EXAFS data, is included in the fits by precomputing and storing EXAFS signals from real-space multiple-scattering calculations on a large number of unique water clusters. On the other hand, fitting also the O-O PCF from diffraction is seen to enhance the amount of structural disorder in the joint fit. Thus, both structures containing nearly linear hydrogen bonds and local structural disorder are important to reproduce diffraction and EXAFS simultaneously. This work also illustrates a few of many possible uses of the SpecSwap-RMC method in modeling disordered materials, particularly for fitting computationally demanding techniques and combining multiple data sets.

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.

Diffraction and IR/Raman data do not prove tetrahedral water.

We use the reverse Monte Carlo modeling technique to fit two extreme structure models for water to available x-ray and neutron diffraction data in q space as well as to the electric field distribution as a representation of the OH stretch Raman spectrum of dilue HOD in D(2)O; the internal geometries were fitted to a quantum distribution. Forcing the fit to maximize the number of hydrogen (H) bonds results in a tetrahedral model with 74% double H-bond donors (DD) and 21% single donors (SD). Maximizing instead the number of SD species gives 81% SD and 18% DD, while still reproducing the experimental data and losing only 0.7-1.8 kJ/mole interaction energy. By decomposing the simulated Raman spectrum we can relate the models to the observed ultrafast frequency shifts in recent pump-probe measurements. Within the tetrahedral DD structure model the assumed connection between spectrum position and H-bonding indicates ultrafast dynamics in terms of breaking and reforming H bonds while in the strongly distorted model the observed frequency shifts do not necessarily imply H-bond changes. Both pictures are equally valid based on present diffraction and vibrational experimental data. There is thus no strict proof of tetrahedral water based on these data. We also note that the tetrahedral structure model must, to fit diffraction data, be less structured than most models obtained from molecular dynamics simulations.