Instantaneous normal mode analysis of hydrated electron solvation dynamics

Ab Initio Investigation of the Resonance Raman Spectrum of the Hydrated Electron

According to the conventional picture, the aqueous or "hydrated" electron, e^-(aq), occupies an excluded volume (cavity) in the structure of liquid water. However, simulations with certain one-electron models predict a more delocalized spin density for the unpaired electron, with no distinct cavity structure. It has been suggested that only the latter (non-cavity) structure can explain the hydrated electron's resonance Raman spectrum, although this suggestion is based on calculations using empirical frequency maps developed for neat liquid water, not for e^-(aq). All-electron ab initio calculations presented here demonstrate that both cavity and non-cavity models of e^-(aq) afford significant red-shifts in the O-H stretching region. This effect is nonspecific and arises due to electron penetration into frontier orbitals of the water molecules. Only the conventional cavity model, however, reproduces the splitting of the H-O-D bend (in isotopically mixed water) that is observed experimentally and arises due to the asymmetric environments of the hydroxyl moieties in the electron's first solvation shell. We conclude that the cavity model of e^-(aq) is more consistent with the measured resonance Raman spectrum than is the delocalized, non-cavity model, despite previous suggestions to the contrary. Furthermore, calculations with hybrid density functionals and with Hartree-Fock theory predict that non-cavity liquid geometries afford only unbound (continuum) states for an extra electron, whereas in reality this energy level should lie more than 3 eV below vacuum level. As such, the non-cavity model of e^-(aq) appears to be inconsistent with available vibrational spectroscopy, photoelectron spectroscopy, and quantum chemistry.

Pressure response of the THz spectrum of bulk liquid water revealed by intermolecular instantaneous normal mode analysis

The radial distribution functions of liquid water are known to change significantly their shape upon hydrostatic compression from ambient conditions deep into the kbar pressure regime. It has been shown that despite their eye-catching changes, the fundamental locally tetrahedral fourfold H-bonding pattern that characterizes ambient water is preserved up to about 10 kbar (1 GPa), which is the stability limit of liquid water at 300 K. The observed increase in coordination number comes from pushing water molecules into the first coordination sphere without establishing an H-bond, resulting in roughly two such additional interstitial molecules at 10 kbar. THz spectroscopy has been firmly established as a powerful experimental technique to analyze H-bonding in aqueous solutions given that it directly probes the far-infrared lineshape and thus the prominent H-bond network mode around 180 cm^-1. We, therefore, set out to assess pressure effects on the THz response of liquid water at 10 kbar in comparison to the 1 bar (0.1 MPa) reference, both at 300 K, with the aim to trace back the related lineshape changes to the structural level. To this end, we employ the instantaneous normal mode approximation to rigorously separate the H-bonding peak from the large background arising from the pronounced librational tail. By exactly decomposing the total molecular dynamics into hindered translations, hindered rotations, and intramolecular vibrations, we find that the H-bonding peak arises from translation-translation and translation-rotation correlations, which are successively decomposed down to the level of distinct local H-bond environments. Our utmost detailed analysis based on molecular pair classifications unveils that H-bonded double-donor water pairs contribute most to the THz response around 180 cm^-1, whereas interstitial waters are negligible. Moreover, short double-donor H-bonds have their peak maximum significantly shifted toward higher frequencies with respect to such long H-bonds. In conjunction with an increasing relative population of these short H-bonds versus the long ones (while the population of other water pair classes is essentially pressure insensitive), this explains not only the blue-shift of the H-bonding peak by about 20-30 cm^-1 in total from 1 bar to 10 kbar but also the filling of the shallow local minimum of the THz lineshape located in between the network peak and the red-wing of the librational band at 1 bar. Based on the changing populations as a function of pressure, we are also able to roughly estimate the pressure-dependence of the H-bond network mode and find that its pressure response and thus the blue-shifting are most pronounced at low kbar pressures.

Structure of the aqueous electron

A cavity or excluded-volume structure best explains the experimental properties of the aqueous or “hydrated” electron.

Infrared photodissociation spectroscopy of ion-radical networks in cationic dimethylamine complexes

Infrared photodissociation spectroscopy was employed to establish the general trends in the stepwise growth motif of cationic dimethylamine (DMA)n+ (n = 4-13) complexes. Electronic structure calculations were performed to identify the structure of the low-lying isomers and to assign the observed spectral features. The results showed the preference of the formation of the proton-transferred (CH3)2NH2+ ion core. The (CH3)2NH2+-[(CH3)2N] ion-radical pair contact and the ion-radical separated pair could coexist at n = 4. The [(CH3)2N] radical is separated from the (CH3)2NH2+ ion core by one DMA molecule at n = 4-6 and by two or more DMA molecules in the larger clusters. This suggests that the (CH3)2NH2+-[(CH3)2N] ion-radical contact pair is not stable in the subsequent radiation-induced processes of DMA, and the [(CH3)2N] radical is released from the charged site in the cationic DMA networks.

Interior- and surface-bound excess electron states in large water cluster anions

We present the results of mixed quantum/classical simulations on relaxed thermal nanoscale water cluster anions, (H(2)O)(n)(-), with n=200, 500, 1000, and 8000. By using initial equilibration with constraints, we investigate stable/metastable negatively charged water clusters with both surface-bound and interior-bound excess electron states. Characterization of these states is performed in terms of geometrical parameters, energetics, and optical absorption spectroscopy of the clusters. The calculations provide data characterizing these states in the gap between previously published calculations and experiments on smaller clusters and the limiting cases of either an excess electron in bulk water or an excess electron at an infinite water/air interface. The present results are in general agreement with previous simulations and provide a consistent picture of the evolution of the physical properties of water cluster anions with size over the entire size range, including results for vertical detachment energies and absorption spectra that would signify their presence. In particular, the difference in size dependence between surface-bound and interior-bound state absorption spectra is dramatic, while for detachment energies the dependence is qualitatively the same.

Dynamics of electrons in ammonia cages: the discovery system of solvation

Two centuries ago solvated electrons were discovered in liquid ammonia and a century later the concept of the solvent cage was introduced. Here, we report a real time study of the dynamics of size-selected clusters, n=20 to 60, of electrons in ammonia, and, for comparison, that of electrons in water cages. Unlike the water case, the observed dynamics for ammonia indicates the formation, through a 100 fs temperature jump, of a solvent collective motion in a 500 fs relaxation process. The agreement of the experimental results-obtained for a well-defined n, gated electron kinetic energy, and time delay-with molecular dynamics theory suggests the critical and different role of the kinetic energy and the librational motions involved in solvation.

Probing the electron delocalization in liquid water and ice at attosecond time scales

We determine electron delocalization rates in liquid water and ice using core-hole decay spectroscopy. The hydrogen-bonded network delocalizes the electrons in less than 500 as. Broken or weak hydrogen bonds--in the liquid or at the surface of ice--provide states where the electron remains localized longer than 20 fs. These asymmetrically bonded water species provide electron traps, acting as a strong precursor channel to the hydrated electron.

Quantized time correlation function approach to nonadiabatic decay rates in condensed phase: Application to solvated electrons in water and methanol

A new, alternative form of the golden rule formula defining the nonadiabatic transition rate between two quantum states in condensed phase is presented. The formula involves the quantum time correlation function of the energy gap, of the nonadiabatic coupling, and their cross terms. Those quantities can be inferred from their classical counterparts, determined via molecular dynamics simulations. The formalism is applied to the problem of the nonadiabatic p-->s relaxation of an equilibrated p-electron in water and methanol. We find that, in both solvents, the relaxation is induced by the coupling to the vibrational modes and the quantum effects modify the rate by a factor of 2-10 depending on the quantization procedure applied. The resulting p-state lifetime for a hypothetical equilibrium excited state appears extremely short, in the sub-100 fs regime. Although this result is in contrast with all previous theoretical predictions, we also illustrate that the lifetimes computed here are very sensitive to the simulated electronic quantum gap and to the strongly correlated nonadiabatic coupling.

Quantum theoretical study of electron solvation dynamics in ice layers on a Cu(111) surface

Recent experiments using time- and angle-resolved two-photon photoemission (2PPE) spectroscopy at metal/polar adsorbate interfaces succeeded in time-dependent analysis of the process of electron solvation. A fully quantum mechanical, two-dimensional simulation of this process, which explicitly includes laser excitation, is presented here, confirming the origin of characteristic features, such as the experimental observation of an apparently negative dispersion. The inference of the spatial extent of the localized electron states from the angular dependence of the 2PPE spectra has been found to be non-trivial and system-dependent.

Electrons in Finite-Sized Water Cavities: Hydration Dynamics Observed in Real Time

We directly observed the hydration dynamics of an excess electron in the finite-sized water clusters of (H2O)n- with n = 15, 20, 25, 30, and 35. We initiated the solvent motion by exciting the hydrated electron in the cluster. By resolving the binding energy of the excess electron in real time with femtosecond resolution, we captured the ultrafast dynamics of the electron in the presolvated ("wet") and hydrated states and obtained, as a function of cluster size, the subsequent relaxation times. The solvation time (300 femtoseconds) after the internal conversion [140 femtoseconds for (H2O)35-] was similar to that of bulk water, indicating the dominant role of the local water structure in the dynamics of hydration. In contrast, the relaxation in other nuclear coordinates was on a much longer time scale (2 to 10 picoseconds) and depended critically on cluster size.

Critical evaluation of approximate quantum decoherence rates for an electronic transition in methanol solution

We present a quantum molecular dynamics calculation of a semiclassical decoherence function to evaluate the accuracy of alternative short-time approximations for coherence loss in the dynamics of condensed phase electronically nonadiabatic processes. The semiclassical function from mixed quantum-classical molecular dynamics simulations and frozen Gaussian wave packets is computed for the electronic transition of an excited state excess electron to the ground state in liquid methanol. The decoherence function decays on a 10 fs time scale that is qualitatively similar to the aqueous case. We demonstrate that it is the motion of the hydrogen atom, and, in particular, the hydrogen rotation around the oxygen-methyl bond which is predominantly responsible for destroying the quantum correlations between alternative states. Multiple time scales due to the slower diffusive nuclear modes, which dominate the solvation response of methanol, do not contribute to the coherence loss. The choice of the coordinate representation is investigated in detail and concluded to be irrelevant to the decay. Changes in both nuclear momenta and positions on the two alternative potential surfaces are found to contribute to decoherence, the former dominating at short times (t < 5 fs), the latter controlling the decay at longer times. Various short-time approximations to the full dynamics for the decoherence function are tested for the first time. The present treatment rigorously develops the short-time description and establishes its range of validity. Whereas the lowest-order short-time approximation proves to be a very good approximation up to about 5 fs, we also find that it bounds the decay of the decoherence function. After 5 fs, the coherence decay in fact becomes faster than the single Gaussian predicted in the lowest-order short-time limit. This decay is well reflected by an enhanced low-order approximation, which is also easily computed from equilibrium classical forces.

Ultrafast dynamics of electron localization and solvation in ice layers on Cu(111)

The femtosecond dynamics of localization and solvation of photoinjected electrons in ultrathin layers of amorphous solid H2O and D2O have been studied by time- and angle-resolved two-photon-photoelectron spectroscopy. After electron transfer from the metal substrate into the conduction band of ice, the excess electron localizes within the first 100 fs in a state at 2.9 eV above E(F), which is further stabilized by 300 meV on a time scale of 0.5-1 ps due to molecular rearrangements in the adlayer. A pronounced change of the solvation dynamics at a coverage of approximately 2 bilayers is attributed to different rigidity of the solvation shell in the bulk and near the surface of ice.