Electronic structure effects in liquid water studied by photoelectron spectroscopy and density functional theory

Attosecond-pump attosecond-probe x-ray spectroscopy of liquid water

Attosecond-pump/attosecond-probe experiments have long been sought as the most straightforward method for observing electron dynamics in real time. Although there has been much success with overlapped near-infrared femtosecond and extreme ultraviolet attosecond pulses combined with theory, true attosecond-pump/attosecond-probe experiments have been limited. We used a synchronized attosecond x-ray pulse pair from an x-ray free-electron laser to study the electronic response to valence ionization in liquid water through all x-ray attosecond transient absorption spectroscopy (AX-ATAS). Our analysis showed that the AX-ATAS response is confined to the subfemtosecond timescale, eliminating any hydrogen atom motion and demonstrating experimentally that the 1b1 splitting in the x-ray emission spectrum is related to dynamics and is not evidence of two structural motifs in ambient liquid water.

Core and Valence Level Photoelectron Spectroscopy of Nanosolvated KCl

The solvation of alkali and halide ions in the aqueous environment has been a subject of intense experimental and theoretical research with multidisciplinary interests; yet, a comprehensive molecular-level understanding has still not been obtained. In recent years, electron spectroscopy has been increasingly applied to study the electronic and structural properties of aqueous ions with implications, especially in atmospheric chemistry. In this work, we report core and valence level (Cl 2p, Cl 3p, and K 3p) photoelectron spectra of the common alkali halide, KCl, doped in gas-phase water clusters in the size range of a few hundred water molecules. The results indicate that the electronic structure of these nanosolutions shows a distinct character from that observed at the liquid-vapor interface in liquid microjets and ambient pressure setups. Insights are provided into the unique solvation properties of ions in a nanoaqueous environment, emerging properties of bulk electrolyte solutions with growing cluster size, and sensitivity of the electronic structure to varying solvation configurations.

Ionization Energy of Liquid Water Revisited

The ionization energy of liquid water is one of its most fundamental properties, an important benchmark for first-principles electronic-structure calculations and a crucial reference in the growing field of liquid-phase photoelectron spectroscopy. Despite this significance, a consensus on its value appears to be missing in the literature. Therefore, we use a monochromatized high-harmonic light source to perform detailed measurements of the ionization energy of liquid water in the presence of a tunable bias voltage applied to the liquid jet. Our results suggest that this simple method is sufficient to simultaneously compensate the effects of the streaming potential and that of the vacuum-level offset between the liquid and the photoelectron spectrometer. Our measurements yield corrected values of the vertical and adiabatic ionization energies of the 1b₁ band of bulk liquid water of 11.67(15) and 10.12(15) eV, respectively. Our method is broadly applicable and is likely to result in corrections to the measured ionization energies of solvated species as well.

Interdependent Electronic Structure, Protonation, and Solvatization of Aqueous 2-Thiopyridone

2-Thiopyridone (2-TP), a common model system for excited-state proton transfer, has been simulated in aqueous solution with ab initio molecular dynamics. The interplay of electronic structure, protonation, and solvatization is investigated by comparison of three differently protonated molecular forms and between the lowest singlet and triplet electronic states. An interdependence clearly manifests in the mixed-character T₁ state for the 2-TP form, systematic structural distortions of the 2-mercaptopyridine (2-MP) form, and photobase protolysis of the 2-TP^- form, in the aqueous phase. In comparison, simplified continuum models for the solvatization are found to be significantly inaccurate for several of the species. To facilitate future computational studies, we therefore present a minimal representative solvatization complex for each stable form and electronic state. Our findings demonstrate the importance of explicit solvatization of the compound and sets the stage for including it also in future studies.

Do water's electrons care about electrolytes?

Ions have a profound effect on the geometrical structure of liquid water and an aqueous environment is known to change the electronic structure of ions. Here we combine photoelectron spectroscopy measurements from liquid microjets with molecular dynamical and quantum chemical calculations to address the reverse question, to what extent do ions affect the electronic structure of liquid water? We study aqueous solutions of sodium iodide (NaI) over a wide concentration range, from nearly pure water to 8 M solutions, recording spectra in the 5 to 60 eV binding energy range to include all water valence and the solute Na+ 2p, I- 4d, and I- 5p orbital ionization peaks. We observe that the electron binding energies of the solute ions change only slightly as a function of electrolyte concentration, less than 150 ± 60 meV over an ∼8 M range. Furthermore, the photoelectron spectrum of liquid water is surprisingly mildly affected as we transform the sample from a dilute aqueous salt solution to a viscous, crystalline-like phase. The most noticeable spectral changes are a negative binding energy shift of the water 1b2 ionizing transition (up to -370 ± 60 meV) and a narrowing of the flat-top shape water 3a1 ionization feature (up to 450 ± 90 meV). A novel computationally efficient technique is introduced to calculate liquid-state photoemission spectra using small clusters from molecular dynamics (MD) simulations embedded in dielectric continuum. This theoretical treatment captured the characteristic positions and structures of the aqueous photoemission peaks, reproducing the experimentally observed narrowing of the water 3a1 feature and weak sensitivity of the water binding energies to electrolyte concentration. The calculations allowed us to attribute the small binding energy shifts to ion-induced disruptions of intermolecular electronic interactions. Furthermore, they demonstrate the importance of considering concentration-dependent screening lengths for a correct description of the electronic structure of solvated systems. Accounting for electronic screening, the calculations highlight the minimal effect of electrolyte concentration on the 1b1 binding energy reference, in accord with the experiments. This leads us to a key finding that the isolated, lowest-binding-energy, 1b1, photoemission feature of liquid water is a robust energetic reference for aqueous liquid microjet photoemission studies.

Femtosecond Extreme Ultraviolet Photoelectron Spectroscopy of Organic Molecules in Aqueous Solution

Time-resolved valence photoelectron spectroscopy is an established tool for studies of ultrafast molecular dynamics in the gas phase. Here we demonstrate time-resolved XUV photoelectron spectroscopy from dilute aqueous solutions of organic molecules, paving the way to application of this method to photodynamics studies of organic molecules in natural environments, which so far have only been accessible to all-optical transient spectroscopies. We record static and time-resolved photoelectron spectra of a sample molecule, quinoline yellow WS, analyze its electronic structure, and follow the relaxation dynamics upon excitation with 400 nm pulses. The dynamics exhibit three time scales, of which a 250 ± 70 fs time scale is attributed to solvent rearrangement. The two longer time scales of 1.3 ± 0.4 and 90 ± 20 ps can be correlated to the recently proposed ultrafast excited-state intramolecular proton transfer in a closely related molecule, quinophthalone.

Aqueous Solvation of Ammonia and Ammonium: Probing Hydrogen Bond Motifs with FT-IR and Soft X-ray Spectroscopy

In a multifaceted investigation combining local soft X-ray and vibrational spectroscopic probes with ab initio molecular dynamics simulations, hydrogen-bonding interactions of two key principal amine compounds in aqueous solution, ammonia (NH₃) and ammonium ion (NH₄⁺), are quantitatively assessed in terms of electronic structure, solvation structure, and dynamics. From the X-ray measurements and complementary determination of the IR-active hydrogen stretching and bending modes of NH₃ and NH₄⁺ in aqueous solution, the picture emerges of a comparatively strongly hydrogen-bonded NH₄⁺ ion via N-H donating interactions, whereas NH₃ has a strongly accepting hydrogen bond with one water molecule at the nitrogen lone pair but only weakly N-H donating hydrogen bonds. In contrast to the case of hydrogen bonding among solvent water molecules, we find that energy mismatch between occupied orbitals of both the solutes NH₃ and NH₄⁺ and the surrounding water prevents strong mixing between orbitals upon hydrogen bonding and, thus, inhibits substantial charge transfer between solute and solvent. A close inspection of the calculated unoccupied molecular orbitals, in conjunction with experimentally measured N K-edge absorption spectra, reveals the different nature of the electronic structural effects of these two key principal amine compounds imposed by hydrogen bonding to water, where a pH-dependent excitation energy appears to be an intrinsic property. These results provide a benchmark for hydrogen bonding of other nitrogen-containing acids and bases.

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.

Dual analyzer system for surface analysis dedicated for angle-resolved photoelectron spectroscopy at liquid surfaces and interfaces

The investigation of liquid surfaces and interfaces with the powerful toolbox of ultra-high vacuum (UHV)-based surface science techniques generally has to overcome the issue of liquid evaporation within the vacuum system. In the last decade, however, new classes of liquids with negligible vapor pressure at room temperature-in particular, ionic liquids (ILs)-have emerged for surface science studies. It has been demonstrated that particularly angle-resolved X-ray Photoelectron Spectroscopy (ARXPS) allows for investigating phenomena that occur at gas-liquid and liquid-solid interfaces on the molecular level. The results are not only relevant for IL systems but also for liquids in general. In all of these previous ARXPS studies, the sample holder had to be tilted in order to change the polar detection angle of emitted photoelectrons, which restricted the liquid systems to very thin viscous IL films coating a flat solid support. We now report on the concept and realization of a new and unique laboratory "Dual Analyzer System for Surface Analysis (DASSA)" which enables fast ARXPS, UV photoelectron spectroscopy, imaging XPS, and low-energy ion scattering at the horizontal surface plane of macroscopically thick non-volatile liquid samples. It comprises a UHV chamber equipped with two electron analyzers mounted for simultaneous measurements in 0° and 80° emission relative to the surface normal. The performance of DASSA on a first macroscopic liquid system will be demonstrated.

Valence Electronic Structure of Aqueous Solutions: Insights from Photoelectron Spectroscopy

The valence orbital electron binding energies of water and of embedded solutes are crucial quantities for understanding chemical reactions taking place in aqueous solution, including oxidation/reduction, transition-metal coordination, and radiation chemistry. Their experimental determination based on liquid-photoelectron spectroscopy using soft X-rays is described, and we provide an overview of valence photoelectron spectroscopy studies reported to date. We discuss principal experimental aspects and several theoretical approaches to compute the measured binding energies of the least tightly bound molecular orbitals. Solutes studied are presented chronologically, from simple electrolytes, via transition-metal ion solutions and several organic and inorganic molecules, to biologically relevant molecules, including aqueous nucleotides and their components. In addition to the lowest vertical ionization energies, the measured valence photoelectron spectra also provide information on adiabatic ionization energies and reorganization energies for the oxidation (ionization) half-reaction. For solutes with low solubility, resonantly enhanced ionization provides a promising alternative pathway.

Application of Koopmans' theorem for density functional theory to full valence-band photoemission spectroscopy modeling

In this work, Koopmans' theorem for Kohn-Sham density functional theory (KS-DFT) is applied to the photoemission spectra (PES) modeling over the entire valence-band. To examine the validity of this application, a PES modeling scheme is developed to facilitate a full valence-band comparison of theoretical PES spectra with experiments. The PES model incorporates the variations of electron ionization cross-sections over atomic orbitals and a linear dispersion of spectral broadening widths. KS-DFT simulations of pristine rubrene (5,6,11,12-tetraphenyltetracene) and potassium-rubrene complex are performed, and the simulation results are used as the input to the PES models. Two conclusions are reached. First, decompositions of the theoretical total spectra show that the dissociated electron of the potassium mainly remains on the backbone and has little effect on the electronic structures of phenyl side groups. This and other electronic-structure results deduced from the spectral decompositions have been qualitatively obtained with the anionic approximation to potassium-rubrene complexes. The qualitative validity of the anionic approximation is thus verified. Second, comparison of the theoretical PES with the experiments shows that the full-scale simulations combined with the PES modeling methods greatly enhance the agreement on spectral shapes over the anionic approximation. This agreement of the theoretical PES spectra with the experiments over the full valence-band can be regarded, to some extent, as a collective validation of the application of Koopmans' theorem for KS-DFT to valence-band PES, at least, for this hydrocarbon and its alkali-adsorbed complex.

The accurate calculation of the band gap of liquid water by means of GW corrections applied to plane-wave density functional theory molecular dynamics simulations

Knowledge about the intrinsic electronic properties of water is imperative for understanding the behaviour of aqueous solutions that are used throughout biology, chemistry, physics, and industry.

Insights into the ultraviolet spectrum of liquid water from model calculations: the different roles of donor and acceptor hydrogen bonds in water pentamers

With a view toward a better understanding of changes in the peak position and shape of the first absorption band of water with condensation or temperature, results from electronic structure calculations using high level wavefunction based and time-dependent density functional methods are reported for water pentamers. Excitation energies, oscillator strengths, and redistributions of electron density are determined for the quasitetrahedral water pentamer in its C(2v) equilibrium geometry and for many pentamer configurations sampled from molecular simulation of liquid water. Excitations associated with surface molecules are removed in order to focus on those states associated with the central molecule, which are the most representative of the liquid environment. The effect of hydrogen bonding on the lowest excited state associated with the central molecule is studied by adding acceptor or donor hydrogen bonds to tetramer and trimer substructures of the C(2v) pentamer, and by sampling liquid-like configurations having increasing number of acceptor or donor hydrogen bonds of the central molecule. Our results provide clear evidence that the blueshift of excitation energies upon condensation is essentially determined by acceptor hydrogen bonds, and the magnitudes of these shifts are determined by the number of such, whereas donor hydrogen bonds do not induce significant shifts in excitation energies. This qualitatively different role of donor and acceptor hydrogen bonds is understood in terms of the different roles of the 1b(1) monomer molecular orbitals, which establishes an intimate connection between the valence hole and excitation energy shifts. Since the valence hole of the lowest excitation associated with the central molecule is found to be well localized in all liquid-like hydrogen bonding environments, with an average radius of gyration of ~1.6 Å that is much lower than the nearest neighbor O-O distance, a clear and unambiguous connection between hydrogen bonding environments and excitation energy shifts can be established. Based on these results, it is concluded that peak position of the first absorption band is mainly determined by the relative distribution of single and double acceptor hydrogen bonding environments, whereas the shape of the first absorption band is mainly determined by the relative distribution of acceptor and broken acceptor hydrogen bonding environments. The temperature dependence of the peak position and shape of the first absorption band can be readily understood in terms of changes to these relative populations.

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.

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.

Insights into the ultraviolet spectrum of liquid water from model calculations

With a view toward a better molecular level understanding of the effects of hydrogen bonding on the ultraviolet absorption spectrum of liquid water, benchmark electronic structure calculations using high level wave function based methods and systematically enlarged basis sets are reported for excitation energies and oscillator strengths of valence excited states in the equilibrium water monomer and dimer and in a selection of liquid-like dimer structures. Analysis of the electron density redistribution associated with the two lowest valence excitations of the water dimer shows that these are usually localized on one or the other monomer, although valence hole delocalization can occur for certain relative orientations of the water molecules. The lowest excited state is mostly associated with the hydrogen bond donor and the significantly higher energy second excited state mostly with the acceptor. The magnitude of the lowest excitation energies is strongly dependent on where the valence hole is created, and only to a lesser degree on the perturbation of the excited electron density distribution by the neighboring water molecule. These results suggest that the lowest excitation energies in clusters and liquid water can be associated with broken acceptor hydrogen bonds, which provide energetically favorable locations for the formation of a valence hole. Higher valence excited states of the dimer typically involve delocalization of the valence hole and/or delocalization of the excited electron and/or charge transfer. Two of the higher valence excited states that involve delocalized valence holes always have particularly large oscillator strengths. Due to the pervasive delocalization and charge transfer, it is suggested that most condensed phase water valence excitations intimately involve more than one water molecule and, as a consequence, will not be adequately described by models based on perturbation of free water monomer states. The benchmark calculations are further used to evaluate a series of representative semilocal, global hybrid, and range separated hybrid functionals used in efficient time-dependent density functional methods. It is shown that such an evaluation is only meaningful when comparison is made at or near the complete basis set limit of the wave function based reference method. A functional is found that quantitatively describes the two lowest excitations of water dimer and also provides a semiquantitative description of the higher energy valence excited states. This functional is recommended for use in further studies on the absorption spectrum of large water clusters and of condensed phase water.

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.

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.