Comment on "Energetics of hydrogen bond network rearrangements in liquid water"

On the existence of soliton-like collective modes in liquid water at the viscoelastic crossover

The problem of large-density variations in supercooled and ambient water has been widely discussed in the past years. Recent studies have indicated the possibility of nanometer-sized density variations on the subpicosecond and picosecond time scales. The nature of fluctuating density heterogeneities remains a highly debated issue. In the present work, we address the problem of possible association of such density variations with the dynamics of terahertz longitudinal acoustic-like modes in liquid water. Our study is based on the fact that the subpicosecond dynamics of liquid water are essentially governed by the structural relaxation. Using a mode coupling theory approach, we found that for typical values of parameters of liquid water, the dynamic mechanism coming from the combination of the structural relaxation process and the finiteness of the amplitude of terahertz longitudinal acoustic-like mode gives rise to a soliton-like collective mode on a temperature-dependent nanometer length scale. The characteristics of this mode are consistent with the estimates of the amplitudes and temperature-dependent correlation lengths of density fluctuations in liquid water obtained in experiments and simulations. Thus, the fully dynamic mechanism could contribute to the formation and dynamics of fluctuating density heterogeneities. The soliton-like collective excitations suggested by our analysis may be relevant to different phenomena connected with supercooled water and can be expected to be associated with some ultrafast biological processes.

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.

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.

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.

Towards a unified description of the hydrogen bond network of liquid water: a dynamics based approach

The definition of a hydrogen bond (H-bond) is intimately related to the topological and dynamic properties of the hydrogen bond network within liquid water. The development of a universal H-bond definition for water is an active area of research as it would remove many ambiguities in the network properties that derive from the fixed definition employed to assign whether a water dimer is hydrogen bonded. This work investigates the impact that an electronic-structure based definition, an energetic, and a geometric definition of the H-bond has upon both topological and dynamic network behavior of simulated water. In each definition, the use of a cutoff (either geometric or energetic) to assign the presence of a H-bond leads to the formation of transiently bonded or broken dimers, which have been quantified within the simulation data. The relative concentration of transient species, and their duration, results in two of the three definitions sharing similarities in either topological or dynamic features (H-bond distribution, H-bond lifetime, etc.), however no two definitions exhibit similar behavior for both classes of network properties. In fact, two networks with similar local network topology (as indicated by similar average H-bonds) can have dramatically different global network topology (as indicated by the defect state distributions) and altered H-bond lifetimes. A dynamics based correction scheme is then used to remove artificially transient H-bonds and to repair artificially broken bonds within the network such that the corrected network exhibits the same structural and dynamic properties for two H-bond definitions (the properties of the third definition being significantly improved). The algorithm described represents a significant step forward in the development of a unified hydrogen bond network whose properties are independent of the original hydrogen bond definition that is employed.

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.

A different view of structure-making and structure-breaking in alkali halide aqueous solutions through x-ray absorption spectroscopy

X-ray absorption spectroscopy measured in transmission mode was used to study the effect of alkali and halide ions on the hydrogen-bonding (H-bonding) network of water. Cl(-) and Br(-) are shown to have insignificant effect on the structure of water while I(-) locally weakens the H-bonding, as indicated by a sharp increase of the main-edge feature in the x-ray absorption spectra. All alkali cations act as structure-breakers in water, weakening the H-bonding network. The spectral changes are similar to spectra of high density ices where the 2nd shell has collapsed due to a break-down of the tetrahedral structures, although here, around the ions, the breakdown of the local tetrahedrality is rather due to non-directional H-bonding to the larger anions. In addition, results from temperature-dependent x-ray Raman scattering measurements of NaCl solution confirm the H-bond breaking effect of Na(+) and the effect on the liquid as similar to an increase in temperature.

Temperature dependent soft x-ray absorption spectroscopy of liquids

A novel sample holder is introduced which allows for temperature dependent soft x-ray absorption spectroscopy of liquids in transmission mode. The setup is based on sample cells with x-ray transmissive silicon nitride windows. A cooling circuit allows for temperature regulation of the sample liquid between -10 °C and +50 °C. The setup enables to record soft x-ray absorption spectra of liquids in transmission mode with a temperature resolution of 0.5 K and better. Reliability and reproducibility of the spectra are demonstrated by investigating the characteristic temperature-induced changes in the oxygen K-edge x-ray absorption spectrum of liquid water. These are compared to the corresponding changes in the oxygen K-edge spectra from x-ray Raman scattering.

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.

Solvation structures of protons and hydroxide ions in water

X-ray Raman spectroscopy (XRS) combined with small-angle x-ray scattering (SAXS) were used to study aqueous solutions of HCl and NaOH. Hydrated structures of H(+) and OH(-) are not simple mirror images of each other. While both ions have been shown to strengthen local hydrogen bonds in the hydration shell as indicated by XRS, SAXS suggests that H(+) and OH(-) have qualitatively different long-range effects. The SAXS structure factor of HCl (aq) closely resembles that of pure water, while NaOH (aq) behaves similar to NaF (aq). We propose that protons only locally enhance hydrogen bonds while hydroxide ions induce tetrahedrality in the overall hydrogen bond network of water.

Polarization dependent resonant x-ray emission spectroscopy of D2O and H2O water: assignment of the local molecular orbital symmetry

The polarization dependence of the split two peaks in the lone-pair region in the x-ray emission spectra has been determined at several different excitation energies for both D(2)O and H(2)O water. In contrast to predictions based on a narrow range of local water structures where the two peaks would be of different molecular orbital symmetry and arise from, respectively, intact and dissociated molecules, we show that the two peaks in the lone-pair region are both of lone-pair 1b(1) orbital symmetry. The results support the interpretation that the two peaks appear due to fluctuations between two distinct different main structural environments.

A sample holder for soft x-ray absorption spectroscopy of liquids in transmission mode

A novel sample holder for soft x-ray absorption spectroscopy of liquids in transmission mode based on sample cells with x-ray transparent silicon nitride membranes is introduced. The sample holder allows for a reliable preparation of ultrathin liquid films with an adjustable thickness in the nm-μm range. This enables measurements of high quality x-ray absorption spectra of liquids in transmission mode, as will be shown for the example of liquid H(2)O, aqueous solutions of 3d-transition metal ions and alcohol-water mixtures. The fine structure of the x-ray absorption spectra is not affected by the sample thickness. No effects of the silicon nitride membranes were observed in the spectra. It is shown how an inhomogeneous thickness of the sample affects the spectra and how this can be avoided.

Increased fraction of low-density structures in aqueous solutions of fluoride

X-ray absorption spectroscopy (XAS) and small angle x-ray scattering (SAXS) were utilized to study the effect of fluoride (F(-)) anion in aqueous solutions. XAS spectra show that F(-) increases the number of strong H-bonds, likely between F(-) and water in the first hydration shell. SAXS data show a low-Q scattering intensity increase similar to the effect of a temperature decrease, suggesting an enhanced anomalous scattering behavior in F(-) solutions. Quantitative analysis revealed that fluoride solutions have larger correlation lengths than chloride solutions with the same cations but shorter compared to pure water. This is interpreted as an increased fraction of tetrahedral low-density structures in the solutions due to the presence of the F(-) ions, which act as nucleation centers replacing water in the H-bonding network and forming stronger H-bonds, but the presence of the cations restricts the extension of strong H-bonds.

Approaches to the solvation of the molecular probe N-methyl-6-quinolone in its excited state

The molecular probe N-methyl-6-quinolone (MQ) gives experimental access to its local chemical environment, e.g. inside a biomolecule. Using ab initio molecular dynamics (MD), it is possible to simulate the time evolution of the Stokes shift as a function of the actual atomistic coupling to the surrounding hydrogen bond network and thus obtain a comprehensive view of the local environment. In contrast to ground state ab initio MD simulations, the choice of a method for excited state MD is nontrivial. Here, we develop a simple and accurate model for the solvation dynamics of MQ in its first excited state.

Influence of individual ionic components on the agglomeration kinetics of silver nanoparticles

The precise characteristic of the agglomeration behavior of colloidal suspensions is of paramount interest to many current studies in nanoscience. This work seeks to elucidate the influence that differently charged salts have on the agglomeration state of a Lee-Meisel-type silver colloid. Moreover, we investigate the influence of the chemical nature of individual ions on their potential to induce agglomeration. Raman spectroscopy and surface-enhanced Raman spectroscopy are used to give insights into mechanistic aspects of the agglomeration process and to assess the differences in the influence of different salts on the agglomeration behavior. Finally, we demonstrate the potential of the measurement procedure used in this work to determine the elementary charge on colloidal NPs.

NMR chemical shifts as a tool to analyze first principles molecular dynamics simulations in condensed phases: the case of liquid water

We present (1)H NMR chemical shift calculations of liquid water based on first principles molecular dynamics simulations under periodic boundary conditions. We focus on the impact of computational parameters on the structural and spectroscopic data, which is an important question for understanding how sensitive the computed (1)H NMR resonances are upon variation of the simulation setup. In particular, we discuss the influence of the exchange-correlation functional and the size of the basis set, the choice for the fictitious electronic mass and the use of pseudopotentials for the nuclear magnetic resonance (NMR) calculation on one hand and the underlying Car-Parrinello-type molecular dynamics simulations on the other hand. Our findings show that the direct effect of these parameters on (1)H shifts is not big, whereas the indirect dependence via the structural data is more important. The (1)H NMR chemical shifts clearly reflect the induced structural changes, illustrating once again the sensitivity of (1)H NMR observables on small changes in the local chemical structure of complex hydrogen-bonded liquids.

Properties of liquid water from a systematic refinement of a high-rank multipolar electrostatic potential

We build on previous work [S. Y. Liem and P. L. A. Popelier, J. Chem. Theory Comput. 4, 353 (2008)], where for the first time, a high-rank multipolar electrostatic potential was used in molecular dynamics simulations of liquid water at a wide range of pressures and temperatures, and using a multipolar Ewald summation. Water is represented as a rigid body, with atomic multipole moments defined by quantum chemical topology partitioning its gas phase electron density. The effect of the level of theory on the local structure of liquid water is systematically addressed. Values for Lennard-Jones (LJ) parameters are optimized, for both oxygen and hydrogen atoms, against bulk properties. The best LJ parameters were then used in a set of simulations at 30 different temperatures (1 atm) and another set at 11 different pressures (at 298 K). Inclusion of the hydrogen LJ parameters significantly increases the self-diffusion coefficient. The behavior of bulk properties was studied and the local water structure analyzed by both radial and spatial distribution functions. Comparisons with familiar point-charge potentials, such as TIP3P, TIP4P, TIP5P, and simple point charge, show the benefits of multipole moments.

Probing the Equilibrium Size and Hydrogen Bonding Structure in Aqueous Aerosol Droplets

Comparative measurements on two or more optically trapped aerosol droplets can allow the hygroscopic properties of aerosol to be investigated in exquisite detail. In this paper we consider two methods for assessing the dependence on relative humidity of the water activity and solute concentration in the condensed phase. We demonstrate in this paper that using a control droplet as a microprobe of relative humidity in the vicinity of a second droplet of interest can allow the RH to be determined with an accuracy of between 0.01 to 0.1 %. Not only is such a probe highly responsive, but the accuracy of the RH determination improves with increase in GF or RH, counter to the performance of macroscopic RH probes. We also demonstrate that spontaneous Raman scattering recorded from excitation of the OH stretching vibration in supersaturated solution droplets of sodium chloride shows the considerable perturbation induced by the chloride ion on the hydrogen bonding network. Chloride is recognised as a hydrogen bonding structure breaker and these observations are shown to be consistent with measurements on more dilute solutions.

X-ray absorption signatures of the molecular environment in water and ice

The x-ray absorption spectra of water and ice are calculated with a many-body approach for electron-hole excitations. The experimental features, including the effects of temperature change in the liquid, are reproduced from configurations generated by ab initio molecular dynamics. The spectral difference between the solid and the liquid is due to two major short-range order effects. One, due to breaking of hydrogen bonds, enhances the pre-edge intensity in the liquid. The other, due to a nonbonded molecular fraction in the first coordination shell, affects the main spectral edge in the conversion of ice to water. This effect may not involve hydrogen bond breaking as shown by experiment in high-density amorphous ice.

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.

Polarizable interaction potential for water from coupled cluster calculations. II. Applications to dimer spectra, virial coefficients, and simulations of liquid water

The six-dimensional CC-pol interaction potential for the water dimer was used to predict properties of the dimer and of liquid water, in the latter case after being supplemented by a nonadditive potential. All the results were obtained purely from first principles, i.e., without any fitting to experimental data. Calculations of the vibration-rotation-tunneling levels of (H(2)O)(2) and (D(2)O)(2), a very sensitive test of the potential surface, gave results in good agreement with experimental high-resolution spectra. Also the virial coefficients and properties of liquid water agree well with measured values. The present model performs better than published force fields for water in a simultaneous reproduction of experimental data for dimer spectra, virials, and properties of the liquid.

Vacuum Ultraviolet (VUV) Photoionization of Small Water Clusters

Tunable vacuum ultraviolet (VUV) photoionization studies of water clusters are performed using 10-14 eV synchrotron radiation and analyzed by reflectron time-of-flight (TOF) mass spectrometry. Photoionization efficiency (PIE) curves for protonated water clusters (H2O)(n)H+ are measured with 50 meV energy resolution. The appearance energies of a series of protonated water clusters are determined from the photoionization threshold for clusters composed of up to 79 molecules. These appearance energies represent an upper limit of the adiabatic ionization energy of the corresponding parent neutral water cluster in the supersonic molecular beam. The experimental results show a sharp drop in the appearance energy for the small neutral water clusters (from 12.62 +/- 0.05 to 10.94 +/- 0.06 eV, for H2O and (H2O)4, respectively), followed by a gradual decrease for clusters up to (H2O)23 converging to a value of 10.6 eV (+/-0.2 eV). The dissociation energy to remove a water molecule from the corresponding neutral water cluster is derived through thermodynamic cycles utilizing the dissociation energies of protonated water clusters reported previously in the literature. The experimental results show a gradual decrease of the dissociation energy for removal of one water molecule for small neutral water clusters (3 <or= n <or= 9). This dissociation energy is discussed within the context of hydrogen bond breaking in a neutral water cluster.

The importance of tetrahedrally coordinated molecules for the explanation of liquid water properties

Ab initio calculations on molecular clusters and a quantum statistical model are used to probe the structure of liquid water and its anomalies. Characteristic temperature dependent mixtures of ring and three-dimensional, voluminous water clusters provide the famous density maximum. The mixture model also reproduces the shift of the density maximum as a function of pressure and isotopic substitution. This finding is consistent with femtosecond spectroscopy data suggesting that two distinct molecular species exist in liquid water. The given structures also reproduce the oxygen-oxygen pair correlation function and the vibrational IR spectrum of liquid water. The results underline the importance of three-dimensional, tetrahedrally coordinated structures for the understanding of water anomalies and the existence of two liquid phases in the supercooled region.

Hydrogen bonds in liquid water studied by photoelectron spectroscopy

The authors report on photoelectron emission spectroscopy measurements of the oxygen 1s orbital of liquid water, using a liquid microjet in ultrahigh vacuum. By suitably changing the soft x-ray photon energy, within 600-1200 eV, the electron probing depth can be considerably altered as to either predominantly access the surface or predominantly bulk water molecules. The absolute probing depth in liquid water was inferred from the evolution of the O1s signal and from comparison with aqueous salt solution. The presence of two distinctive components in the core-level photoelectron spectrum, with significantly different binding energies, is revealed. The dominant contribution, at a vertical binding energy of 538.1 eV, was found in bulk and surface sensitive spectra. A weaker component at 536.6 eV binding energy appears to be present only in bulk water. Hartree-Fock calculations of O1s binding energies in different geometric arrangements of the water network are presented to rationalize the experimental distribution of O1s electron binding energies.

Are recent water models obtained by fitting diffraction data consistent with infrared/Raman and x-ray absorption spectra?

X-ray absorption (XA) spectra have been computed based on water structures obtained from a recent fit to x-ray and neutron diffraction data using models ranging from symmetrical to asymmetrical local coordination of the water molecules [A. K. Soper, J. Phys.: Condens. Matter 17, S3273 (2005)]. It is found that both the obtained symmetric and asymmetric structural models of water give similar looking XA spectra, which do not match the experiment. The fitted models both contain unphysical structures that are allowed by the diffraction data, where, e.g., hydrogen-hydrogen interactions may occur. A modification to the asymmetric model, in which the non-hydrogen-bonded OH intramolecular distance is allowed to become shorter while the bonded OH distance becomes longer, improves the situation somewhat, but the overall agreement is still unsatisfactory. The electric field (E-field) distributions and infrared (IR) spectra are also calculated using two established theoretical approaches, which, however, show significant discrepancies in their predictions for the asymmetric structural models. Both approaches predict the Raman spectrum of the symmetric model fitted to the diffraction data to be significantly blueshifted compared to experiment. At the moment no water model exists that can equally well describe IR/Raman, x-ray absorption spectroscopy, and diffraction data.

Predictions of the Properties of Water from First Principles

A force field for water has been developed entirely from first principles, without any fitting to experimental data. It contains both pairwise and many-body interactions. This force field predicts the properties of the water dimer and of liquid water in excellent agreement with experiments, a previously elusive objective. Precise knowledge of the intermolecular interactions in water will facilitate a better understanding of this ubiquitous substance.

Nuclear Magnetic Resonance Chemical Shifts and Quadrupole Couplings for Different Hydrogen-Bonding Cases Occurring in Liquid Water: A Computational Study

Nuclear magnetic resonance (NMR) parameters are determined theoretically for the oxygen and hydrogen/deuterium nuclei of differently hydrogen-bonded water molecules in liquid water at 300 K. The parameters are the chemical shift, the shielding anisotropy, the asymmetry parameter of shielding, the nuclear quadrupole coupling constant, and the asymmetry parameter of the nuclear quadrupole coupling. We sample instantaneous configurations from a Car-Parrinello molecular dynamics simulation and feed nuclear coordinates into a quantum chemical program for the calculation of NMR parameters using density-functional theory with the three-parameter hybrid exchange-correlation (B3LYP) functional. In the subsequent analysis, molecules are divided into groups according to the number of hydrogen bonds they possess, and the full average NMR tensors are calculated separately for each group. The classification of the hydrogen-bonding cases is performed using a simple distance-based criterion. The analysis reveals in detail how the NMR tensors evolve as the environment changes gradually from gas to liquid upon increasing the number of hydrogen bonds to the molecule of interest. Liquid-state distributions of the instantaneous values of the NMR properties show a wide range of values for each hydrogen-bonding species with significant overlap between the different cases. Our study shows how local changes in the environment, along with classical thermal averaging, affect the NMR parameters in liquid water. For example, a broken or alternatively extra hydrogen bond induces major changes in the NMR tensors, and the effect is more pronounced for hydrogen or deuterium than for oxygen. The data sheds light on the usefulness of NMR experiments in investigating the local coordination of liquid water.