Degree of initial hole localization/delocalization in ionized water clusters

Attosecond spectroscopy of size-resolved water clusters

Electron dynamics in water are of fundamental importance for a broad range of phenomena^1-3, but their real-time study faces numerous conceptual and methodological challenges^4-6. Here, we introduce attosecond size-resolved cluster spectroscopy and build up a molecular-level understanding of the attosecond electron dynamics in water. We measure the effect that the addition of single water molecules has on the photoionization time delays^7-9 of water clusters. We find a continuous increase of the delay for clusters containing up to 4-5 molecules and little change towards larger clusters. We show that these delays are proportional to the spatial extension of the created electron hole, which first increases with cluster size and then partially localizes through the onset of structural disorder that is characteristic of large clusters and bulk liquid water. These results suggest a previously unknown sensitivity of photoionization delays to electron-hole delocalization and indicate a direct link between electronic structure and attosecond photoionization dynamics. Our results offer novel perspectives for studying electron/hole delocalization and its attosecond dynamics.

Electronic Structure and Vibrational Signatures of the Delocalized Radical in Hydrated Clusters of Copper("II") Hydroxide CuOH+(H2O)0-2

The copper hydroxide ion, CuOH⁺, serves as the catalytic core in several recently developed water-splitting catalysts, and an understanding of its chemistry is critical to determining viable catalytic mechanisms. In spite of its importance, the electronic structure of this open-shell ion has remained ambiguous in the literature. In particular, computed values for both the thermodynamics of hydration and the vibrational signatures of the mono- and dihydrates have shown prohibitively large errors compared to values from recent experimental measurements. In this work, the source of this discrepancy is demonstrated to be the propensity of this ion to exist between traditional Cu(I) and Cu(II) oxidation-state limits. The spin density of the radical is accordingly shown to delocalize between the metal center and surrounding ligands, and increasing the hydration serves to exacerbate this behavior. Equation-of-motion coupled-cluster methods demonstrated the requisite accuracy to resolve the thermodynamic discrepancies. Such methods were also needed for spectral simulations, although the latter also required a direct simulation of the role of the deuterium "tag" molecules that are used in modern predissociation spectroscopy experiments. This nominally benign tag molecule underwent direct complexation with the open-valence metal ion, thereby forming a species akin to known metal-H₂ complexes and strongly impacting the resulting spectrum. Thermal populations of this configuration and other more traditional noncovalently bound isomers led to a considerable broadening of the spectral lineshapes. Therefore, at least for the CuOH⁺(H₂O)_0-2 hydrates, these benchmark ions should be considered to be delocalized radical systems with some degree of multireference character at equilibrium. They also serve as a cautionary tale for the spectroscopy community, wherein the role of the D₂ tag is far from benign.

Probing the Electronic Structure of Bulk Water at the Molecular Length Scale with Angle-Resolved Photoelectron Spectroscopy

We report a combined experimental and theoretical study of bulk water photoionization. Angular distributions of photoelectrons produced by ionizing the valence bands of neat water using X-ray radiation (250-750 eV) show a limited (∼20%) decrease in the β anisotropy parameter compared to the gas phase, indicating that the electronic structure of the individual water molecules can be probed. We show that, in the high-energy regime, photoionization of bulk can be described using an incoherent superposition of individual molecules, in contrast to a low-energy regime where photoionization probes delocalized entangled states of molecular aggregates. The two regimes-low versus high energy-are limiting cases where the de Broglie wavelength of the photoelectron is larger or smaller than the intermolecular distance between water molecules, respectively. The comparison of measured and computed anisotropies reveals that the reduction in β at high kinetic energies is mostly due to scattering rather than rehybridization due to solvation.

Thermal Equilibration Controls H-Bonding and the Vertical Detachment Energy of Water Cluster Anions

One of the outstanding puzzles in the photoelectron spectroscopy of water anion clusters, which serve as precursors to the hydrated electron, is that the excess electron has multiple vertical detachment energies (VDEs), with different groups seeing different distributions of VDEs. We have studied the photoelectron spectroscopy of water cluster anions using simulation techniques designed to mimic the different ways that water cluster anions are produced experimentally. Our simulations take advantage of density functional theory-based Born-Oppenheimer molecular dynamics with an optimally tuned range-separated hybrid functional that is shown to give outstanding accuracy for calculating electron binding energies for this system. We find that our simulations are able to accurately reproduce the experimentally observed VDEs for cluster anions of different sizes, with different VDE distributions observed depending on how the water cluster anions are prepared. For cluster anion sizes up to 20 water molecules, we see that the excess electron always resides on the surface of the cluster and that the different discrete VDEs result from the discrete number of hydrogen bonds made to the electron by water molecules on the surface. Clusters that are less thermally equilibrated have surface waters that tend to make single H-bonds to the electron, resulting in lower VDEs, while clusters that are more thermally equilibrated have surface waters that prefer to make two H-bonds to the electron, resulting in higher VDEs.

Reactivity and energy level of a localized hole in liquid water

Reaction and redox level of hole capture in liquid water from first principles.

Quantum molecular motion in the mixed ion-radical complex, [(H2O)(H2S)]. Phys Chem Chem Phys 2016;18:27450-27459. [PMID: 27711703 DOI: 10.1039/c6cp05299a]

The cation dimer of water and hydrogen sulfide, [(H₂O)(H₂S)]⁺, serves as a fundamental model for the oxidation chemistry of H₂S. The known oxidative metabolism of H₂S by biological species in sulfur-rich environments has motivated the study of the inherent properties of this benchmark complex, with possible mechanistic implications for modern water oxidation chemistry. The low-energy isomer of this open-shell ion is a proton-transferred (PT) structure, H₃O⁺SH˙. An alternative PT structure, H₃S⁺OH˙, and a hemibonded (HB) isomer, [H₂O·SH₂]⁺, are also stable isomers, placing this complex intermediate to known (H₂O)₂⁺ (PT) and (H₂S)₂⁺ (HB) limiting regimes. This intermediate character suggested the possibility of unique molecular motion, even in the vibrational ground state. Path integral molecular dynamics and anharmonic vibrational spectroscopy simulations have been performed in this study, in order to understand the inherent quantum molecular motion of this complex. The resulting structural distributions were found to deviate significantly from both classical and harmonic analyses, including the observation of large-amplitude anharmonic motion of the central proton and nearly free rotation of the terminal hydrogens. The predicted vibrational spectra are particularly unique and suggest characteristic signatures of the strong electronic interactions and anharmonic vibrational mode couplings in this radical cation.

Nonadiabatic dynamics of floppy hydrogen bonded complexes: the case of the ionized ammonia dimer

Non-adiabatic dynamics of a floppy hydrogen bonded ammonia dimer was studied by ab initio molecular dynamics simulations.

Ionization-Induced Structural Changes in Uracil Dimers and Their Spectroscopic Signatures

The electronic structure of the three representative isomers of the ionized uracil dimers is characterized by high-level electronic structure calculations. Noncovalent interactions between the fragments lower the vertical ionization energies by 0.13-0.35 eV, the largest drop being observed for the stacked and the T-shaped isomers. The initial hole is delocalized in the stacked and the H-bonded isomers and is localized in the T-shaped one. The ionization induces significant structural relaxation and increases the binding energies. The stacked dimer cation relaxes to the symmetric structure bound by 22.7 kcal/mol. The T-shaped dimer cation has a binding energy of 25.1 kcal/mol. Thus, the relative order of the stacked and T-shaped isomers is reversed upon ionization. Finally, the H-bonded isomer, which relaxes to the proton-transferred structure, is bound by 37.0 kcal/mol. The electronic spectra of all three isomers characterized at the vertical and the relaxed geometries show different patterns, which may be exploited in spectroscopic probing of ionization-induced dynamics in these species.

Structural evolution and solvation of the OH radical in ionized water radical cations (H2O)n+, n = 5–8

Structural evolution of ionized water radical cations (H₂O)_n⁺, n = 5–8, is studied by ab intio methods.

Four Bases Score a Run: Ab Initio Calculations Quantify a Cooperative Effect of H-Bonding and π-Stacking on the Ionization Energy of Adenine in the AATT Tetramer

Benchmark calculations of the lowest ionized state of the (A:T)2 (mixed adenine-thymine) cluster at the geometry taken from the DNA X-ray structure are presented. Vertical ionization energies (IEs) computed by the equation-of-motion coupled-cluster method with single and double substitutions are reported and analyzed. The shift in IE relative to the monomer (A) is -0.7 eV. The performance of the widely used B3LYP, ωB97X-D, and M06-2X functionals with respect to their ability to describe energetics and the character (localization versus delocalization) of the ionized states is also investigated. The shifts in IEs caused by H-bonding and stacking interactions are analyzed in terms of additive versus cooperative effects. It is found that the cooperative effect accounts for more than 20% of the shift in IE relative to the monomer. The cooperative effect and, consequently, the magnitude of the shift are well reproduced by the hybrid quantum mechanics/molecular mechanics scheme in which neutral thymine bases are represented by point charges.

First-principle protocol for calculating ionization energies and redox potentials of solvated molecules and ions: theory and application to aqueous phenol and phenolate

The effect of hydration on the lowest vertical ionization energy (VIE) of phenol and phenolate solvated in bulk water was characterized using the equation-of-motion ionization potential coupled-cluster (EOM-IP-CCSD) and effective fragment potential (EFP) methods (referred to as EOM/EFP) and determined experimentally by valence photoemission measurements using microjets and synchrotron radiation. The computed solvent-induced shifts in VIEs (ΔVIEs) are -0.66 and +5.72 eV for phenol and phenolate, respectively. Our best estimates of the absolute values of VIEs (7.9 and 7.7 eV for phenol and phenolate) agree reasonably well with the respective experimental values (7.8 ± 0.1 and 7.1 ± 0.1 eV). The EOM/EFP scheme was benchmarked against full EOM-IP-CCSD using microsolvated phenol and phenolate clusters. A protocol for calculating redox potentials with EOM/EFP was developed based on linear response approximation (LRA) of free energy determination. The oxidation potentials of phenol and phenolate calculated using LRA and EOM/EFP are 1.32 and 0.89 V, respectively; they agree well with experimental values.

Chasing charge localization and chemical reactivity following photoionization in liquid water

The ultrafast dynamics of the cationic hole formed in bulk liquid water following ionization is investigated by ab initio molecular dynamics simulations and an experimentally accessible signature is suggested that might be tracked by femtosecond pump-probe spectroscopy. This is one of the fastest fundamental processes occurring in radiation-induced chemistry in aqueous systems and biological tissue. However, unlike the excess electron formed in the same process, the nature and time evolution of the cationic hole has been hitherto little studied. Simulations show that an initially partially delocalized cationic hole localizes within ~30 fs after which proton transfer to a neighboring water molecule proceeds practically immediately, leading to the formation of the OH radical and the hydronium cation in a reaction which can be formally written as H(2)O(+) + H(2)O → OH + H(3)O(+). The exact amount of initial spin delocalization is, however, somewhat method dependent, being realistically described by approximate density functional theory methods corrected for the self-interaction error. Localization, and then the evolving separation of spin and charge, changes the electronic structure of the radical center. This is manifested in the spectrum of electronic excitations which is calculated for the ensemble of ab initio molecular dynamics trajectories using a quantum mechanics/molecular mechanics (QM∕MM) formalism applying the equation of motion coupled-clusters method to the radical core. A clear spectroscopic signature is predicted by the theoretical model: as the hole transforms into a hydroxyl radical, a transient electronic absorption in the visible shifts to the blue, growing toward the near ultraviolet. Experimental evidence for this primary radiation-induced process is sought using femtosecond photoionization of liquid water excited with two photons at 11 eV. Transient absorption measurements carried out with ~40 fs time resolution and broadband spectral probing across the near-UV and visible are presented and direct comparisons with the theoretical simulations are made. Within the sensitivity and time resolution of the current measurement, a matching spectral signature is not detected. This result is used to place an upper limit on the absorption strength and/or lifetime of the localized H(2)O(+) ((aq)) species.

Ionization of Water Clusters Mediated by Exciton Energy Transfer from Argon Clusters

The exciton energy deposited in an argon cluster (Arn, ⟨n = 20⟩) using VUV radiation is transferred to softly ionize doped water clusters ((H2O)n, n = 1-9), leading to the formation of nonfragmented clusters. Following the initial excitation, electronic energy is channeled to ionize the doped water cluster while evaporating the Ar shell, allowing identification of fragmented and complete water cluster ions. Examination of the photoionization efficiency curve shows that cluster evaporation from excitons located above 12.6 eV is not enough to cool the energized water cluster ion and leads to their dissociation to (H2O)n-2H(+) (protonated) clusters.

Effect of solvation on the vertical ionization energy of thymine: from microhydration to bulk

The effect of hydration on the vertical ionization energy (VIE) of thymine was characterized using equation-of-motion ionization potential coupled-cluster (EOM-IP-CCSD) and effective fragment potential (EFP) methods. We considered several microsolvated clusters as well as thymine solvated in bulk water. The VIE in bulk water was computed by averaging over solvent-solute configurations obtained from equilibrium molecular dynamics trajectories at 300 K. The effect of microsolvation was analyzed and contrasted against the combined effect of the first solvation shell in bulk water. Microsolvation reduces the ionization energy (IE) by about 0.1 eV per water molecule, while the first solvation shell increases the IE by 0.1 eV. The subsequent solvation lowers the IE, and the bulk value of the solvent-induced shift of thymine's VIE is approximately -0.9 eV. The combined effect of the first solvation shell was explained in terms of specific solute-solvent interactions, which were investigated using model structures. The convergence of IE to the bulk value requires the hydration sphere of approximately 13.5 Å radius. The performance of the EOM-IP-CCSD/EFP scheme was benchmarked against full EOM-IP-CCSD using microhydrated structures. The errors were found to be less than 0.01-0.02 eV. The relative importance of the polarization and higher multipole moments in EFP model was also investigated.