Probing the Local Structure of Liquid Water by X-ray Absorption Spectroscopy†

Ion association in hydrothermal aqueous NaCl solutions: implications for the microscopic structure of supercritical water

Knowledge of the microscopic structure of fluids and changes thereof with pressure and temperature is important for the understanding of chemistry and geochemical processes. In this work we investigate the influence of sodium chloride on the hydrogen-bond network in aqueous solution up to supercritical conditions. A combination of in situ X-ray Raman scattering and ab initio molecular dynamics simulations is used to probe the oxygen K-edge of the alkali halide aqueous solution in order to obtain unique information about the oxygen's local coordination around the ions, e.g. solvation-shell structure and the influence of ion pairing. The measured spectra exhibit systematic temperature dependent changes, which are entirely reproduced by calculations on the basis of structural snapshots obtained via ab initio molecular dynamics simulations. Analysis of the simulated trajectories allowed us to extract detailed structural information. This combined analysis reveals a net destabilizing effect of the dissolved ions which is reduced with rising temperature. The observed increased formation of contact ion pairs and occurrence of larger polyatomic clusters at higher temperatures can be identified as a driving force behind the increasing structural similarity between the salt solution and pure water at elevated temperatures and pressures with drawback on the role of hydrogen bonding in the hot fluid. We discuss our findings in view of recent results on hot NaOH and HCl aqueous fluids and emphasize the importance of ion pairing in the interpretation of the microscopic structure of water.

PSIXAS: A Psi4 plugin for efficient simulations of X-ray absorption spectra based on the transition-potential and Δ-Kohn-Sham method

Near edge X-ray absorption fine structure (NEXAFS) spectra and their pump-probe extension (PP-NEXAFS) offer insights into valence- and core-excited states. We present PSIXAS, a recent implementation for simulating NEXAFS and PP-NEXAFS spectra by means of the transition-potential and the Δ-Kohn-Sham method. The approach is implemented in form of a software plugin for the Psi4 code, which provides access to a wide selection of basis sets as well as density functionals. We briefly outline the theoretical foundation and the key aspects of the plugin. Then, we use the plugin to simulate PP-NEXAFS spectra of thymine, a system already investigated by others and us. It is found that larger, extended basis sets are needed to obtain more accurate absolute resonance positions. We further demonstrate that, in contrast to ordinary NEXAFS simulations, where the choice of the density functional plays a minor role for the shape of the spectrum, for PP-NEXAFS simulations the choice of the density functional is important. Especially hybrid functionals (which could not be used straightforwardly before to simulate PP-NEXAFS spectra) and their amount of "Hartree-Fock like" exact exchange affects relative resonance positions in the spectrum.

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.

The simulation of X-ray absorption spectra from ground and excited electronic states using core-valence separated DFT/MRCI

We present an extension of the combined density functional theory (DFT) and multireference configuration interaction (MRCI) method (DFT/MRCI) [S. Grimme and M. Waletzke, J. Chem. Phys. 111, 5645 (1999)] for the calculation of core-excited states based on the core-valence separation (CVS) approximation. The resulting method, CVS-DFT/MRCI, is validated via the simulation of the K-edge X-ray absorption spectra of 40 organic chromophores, amino acids, and nucleobases, ranging in size from CO₂ to tryptophan. Overall, the CVS-DFT/MRCI method is found to yield accurate X-ray absorption spectra (XAS), with consistent errors in peak positions of ∼2.5-3.5 eV. Additionally, we show that the CVS-DFT/MRCI method may be employed to simulate XAS from valence excited states and compare the simulated spectra to those computed using the established wave function-based approaches [ADC(2) and ADC(2)x]. In general, each of the methods yields excited state XAS spectra in qualitative and often quantitative agreement. In the instances where the methods differ, the CVS-DFT/MRCI simulations predict intensity for transitions for which the underlying electronic states are characterized by doubly excited configurations relative to the ground state configuration. Here, we aim to demonstrate that the CVS-DFT/MRCI approach occupies a specific niche among numerous other electronic structure methods in this area, offering the ability to treat initial states of arbitrary electronic character while maintaining a low computational cost and comparatively black box usage.

Probing hydrogen bond strength in liquid water by resonant inelastic X-ray scattering

Local probes of the electronic ground state are essential for understanding hydrogen bonding in aqueous environments. When tuned to the dissociative core-excited state at the O1s pre-edge of water, resonant inelastic X-ray scattering back to the electronic ground state exhibits a long vibrational progression due to ultrafast nuclear dynamics. We show how the coherent evolution of the OH bonds around the core-excited oxygen provides access to high vibrational levels in liquid water. The OH bonds stretch into the long-range part of the potential energy curve, which makes the X-ray probe more sensitive than infra-red spectroscopy to the local environment. We exploit this property to effectively probe hydrogen bond strength via the distribution of intramolecular OH potentials derived from measurements. In contrast, the dynamical splitting in the spectral feature of the lowest valence-excited state arises from the short-range part of the OH potential curve and is rather insensitive to hydrogen bonding.

Introducing Ionic-Current Detection for X-ray Absorption Spectroscopy in Liquid Cells

Photons and electrons are two common relaxation products upon X-ray absorption, enabling fluorescence yield and electron yield detections for X-ray absorption spectroscopy (XAS). The ions that are created during the electron yield process are relaxation products too, which are exploited in this study to produce ion yield for XA detection. The ionic currents measured in a liquid cell filled with water or iron(III) nitrate aqueous solutions exhibit characteristic O K-edge and Fe L-edge absorption profiles as a function of excitation energy. Application of two electrodes installed in the cell is crucial for obtaining the XA spectra of the liquids behind membranes. Using a single electrode can only probe the species adsorbed on the membrane surface. The ionic-current detection, termed as total ion yield (TIY) in this study, also produces an undistorted Fe L-edge XA spectrum, indicating its promising role as a novel detection method for XAS studies in liquid cells.

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.

Essential dynamics for the study of microstructures in liquids

Essential Dynamics (ED) is a powerful tool for analyzing molecular dynamics (MD) simulations and it is widely adopted for conformational analysis of large molecular systems such as, for example, proteins and nucleic acids. In this study, we extend the use of ED to the study of clusters of arbitrary size constituted by weakly interacting particles, for example, atomic clusters and supramolecular systems. The key feature of the method we present is the identification of the relevant atomic-molecular clusters to be analyzed by ED for extracting the information of interest. The application of this computational approach allows a straightforward and unbiased conformational study of the local microstructures in liquids, as emerged from semiclassical MD simulations. The good performance of the method is demonstrated by calculating typical observables of liquid water, that is, NMR, NEXAFS O1s, and IR spectra, known to be rather sensitive both to the presence and to the conformational features of hydrogen-bonded clusters.

Mass accommodation of water: bridging the gap between molecular dynamics simulations and kinetic condensation models

The condensational growth of submicrometer aerosol particles to climate relevant sizes is sensitive to their ability to accommodate vapor molecules, which is described by the mass accommodation coefficient. However, the underlying processes are not yet fully understood. We have simulated the mass accommodation and evaporation processes of water using molecular dynamics, and the results are compared to the condensation equations derived from the kinetic gas theory to shed light on the compatibility of the two. Molecular dynamics simulations were performed for a planar TIP4P-Ew water surface at four temperatures in the range 268-300 K as well as two droplets, with radii of 1.92 and 4.14 nm at T = 273.15 K. The evaporation flux from molecular dynamics was found to be in good qualitative agreement with that predicted by the simple kinetic condensation equations. Water droplet growth was also modeled with the kinetic multilayer model KM-GAP of Shiraiwa et al. [Atmos. Chem. Phys. 2012, 12, 2777]. It was found that, due to the fast transport across the interface, the growth of a pure water droplet is controlled by gas phase diffusion. These facts indicate that the simple kinetic treatment is sufficient in describing pure water condensation and evaporation. The droplet size was found to have minimal effect on the value of the mass accommodation coefficient. The mass accommodation coefficient was found to be unity (within 0.004) for all studied surfaces, which is in agreement with previous simulation work. Additionally, the simulated evaporation fluxes imply that the evaporation coefficient is also unity. Comparing the evaporation rates of the mass accommodation and evaporation simulations indicated that the high collision flux, corresponding to high supersaturation, present in typical molecular dynamics mass accommodation simulations can under certain conditions lead to an increase in the evaporation rate. Consequently, in such situations the mass accommodation coefficient can be overestimated, but in the present cases the corrected values were still close to unity with the lowest value at ≈0.99.

Structure and dynamics of water dangling OH bonds in hydrophobic hydration shells. Comparison of simulation and experiment

Molecular dynamics and electric field strength simulations are performed in order to quantify the structural, dynamic, and vibrational properties of non-H-bonded (dangling) OH groups in the hydration shell of neopentane, as well as in bulk water. The results are found to be in good agreement with the experimentally observed high-frequency (∼3660 cm(-1)) OH band arising from the hydration shell of neopentanol dissolved in HOD/D(2)O, obtained by analyzing variable concentration Raman spectra using multivariate curve resolution (Raman-MCR). The simulation results further indicate that hydration shell dangling OH groups preferentially point toward the central carbon atom of neopentane to a degree that increases with the lifetime of the dangling OH.

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.

High-resolution soft X-ray photoelectron spectroscopy of liquid water

High-resolution soft X-ray photoelectron spectra of liquid water (H(2)O and D(2)O) were measured using a liquid beam photoelectron spectrometer. The 1a(1) (O1s) band and the lowest valence 1b(1) band had single peaks, which is not consistent with the split 1b(1)→ 1a(1) of the X-ray emission band of liquid water if the splitting is assumed to originate from level shifts in two different hydrogen bonding structures. The second valence 3a(1) band of liquid water exhibited a flat top implying that two bands exist underneath a broad feature, which is similar to the case of the 3a(1) band of amorphous ice. The energy splitting between the two 3a(1) bands is estimated to be 1.38 eV (H(2)O) and 1.39 eV (D(2)O). Ab initio calculations suggest that the large splitting of the 3a(1) band is characteristic of water molecules that function as both proton donor and acceptor. The overall result is consistent with the conventional model of a tetrahedral hydrogen-bonding network in liquid water.

A view of hydrogen/hydroxide flux across lipid membranes

A topic emerging roughly 30 years ago and engendering an incompletely resolved controversy is reviewed in this article: the relatively high permeability and pH independence associated with H(+)/OH(-) passive movements across lipid membranes. We summarize the expected characteristics of simple H(+)/OH(-) diffusion and those of a reaction between H(+) and OH(-) being attracted from opposite surfaces and condensing in an interfacial zone of the membrane. An interfacial H(+)/OH(-) reaction mechanism gives the experimentally observed behavior of an H(+)/OH(-) flux that is independent of the pH measurement range. This mechanism assumes that H(+) and OH(-) within the interfacial zone become electrostatically aligned on opposite sides of the hydrophobic membrane core. Electrostatic attraction and charge delocalization among a small cluster of water molecules surrounding the ions reduce the Born energy for H(+)/OH(-) insertion into lipid. This transmembrane condensation model predicts the magnitude of the experimentally determined H(+)/OH(-) flux, which is significantly greater than that of other monovalent ions. The consequences of an elevated H(+)/OH(-) permeability compared to other ions and the relative pH independence of this flux have consequences for understanding the chemical evolution of life.

Temperature-dependent structure of liquid water investigated by means of energy-dispersive x-ray diffraction and molecular dynamics simulations

X-ray diffraction is investigated in the supercooled region (263-273 K) and, for comparison, at ambient conditions. By means of an energy-dispersive set-up coherent x-ray scattering cross sections within a wide range of momentum transfer values are monitored simultaneously using in-house laboratory equipment. The data are corrected for the energy-dependent detector response, geometric broadening effects, sample absorption and Compton scattering. The corrected diffraction curves are in good agreement with results from the Advanced Light Source (ALS). Additional comparisons with available literature data and with computer simulation results of different rigid water models are presented, relating the scattering intensities to the microscopic H-bond structure and dynamics.

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.

Coefficients of evaporation and gas phase diffusion of low-volatility organic solvents in nitrogen from interferometric study of evaporating droplets

Evaporation of motionless, levitating droplets of pure, low-volatility liquids was studied with interferometric methods. Experiments were conducted on charged droplets in the electrodynamic trap in nitrogen at atmospheric pressure at 298 K. Mono-, di-, tri-, and tetra(ethylene glycols) and 1,3-dimethyl-2-imidazolidinone were studied. The influence of minute impurities (<0.1%) upon the process of droplet evaporation was observed and discussed. The gas phase diffusion and evaporation coefficients were found from droplet radii evolution under the assumption of known vapor pressure. Diffusion coefficients were compared with independent measurements and calculations (in air). Good agreement was found for mono- and di(ethylene glycols), and for 1,3-dimethyl-2-imidazolidinone, which confirmed the used vapor pressure values. The value of equilibrium vapor pressure for tri(ethylene glycol) was proposed to be 0.044 +/- 0.008 Pa. The evaporation coefficient was found to increase from 0.035 to 0.16 versus the molecular mass of the compound.

Role of non-hydrogen-bonded molecules in the oxygen K-edge spectrum of ice

We report the oxygen K-edge spectra of ices Ih, VI, VII, and VIII measured with X-ray Raman scattering. The pre-edge and main-edge contributions increase strongly with density, even though the hydrogen bond arrangements are very similar in these phases. While the near-edge spectral features in water and ice have often been linked to hydrogen bonding, we show that the spectral changes in the phases studied here can be quantitatively related to structural changes in the second coordination shell. Density-functional theory calculations reproduce the experimental results and support the conclusion. Our results suggest that non-hydrogen-bonded neighbors can have a significant effect also in the liquid water spectrum. We discuss the implications of the results for the actively debated interpretation of the liquid water spectrum in terms of local structure.

Probing the hydrogen-bond network of water via time-resolved soft X-ray spectroscopy

We report time-resolved studies of hydrogen bonding in liquid H(2)O, in response to direct excitation of the O-H stretch mode at 3 mum, probed via soft X-ray absorption spectroscopy at the oxygen K-edge. This approach employs a newly developed nanofluidic cell for transient soft X-ray spectroscopy in the liquid phase. Distinct changes in the near-edge spectral region (XANES) are observed, and are indicative of a transient temperature rise of 10 K following transient laser excitation and rapid thermalization of vibrational energy. The rapid heating occurs at constant volume and the associated increase in internal pressure, estimated to be 8 MPa, is manifested by distinct spectral changes that differ from those induced by temperature alone. We conclude that the near-edge spectral shape of the oxygen K-edge is a sensitive probe of internal pressure, opening new possibilities for testing the validity of water models and providing new insight into the nature of hydrogen bonding in water.