Electronic Levels of Excess Electrons in Liquid Water

Nuclear Quantum Effects on the Electronic Structure of Water and Ice

The electronic properties and optical response of ice and water are intricately shaped by their molecular structure, including the quantum mechanical nature of the hydrogen atoms. Despite numerous previous studies, a comprehensive understanding of the nuclear quantum effects (NQEs) on the electronic structure of water and ice at finite temperatures remains elusive. Here, we utilize molecular simulations that harness efficient machine-learning potentials and many-body perturbation theory to assess how NQEs impact the electronic bands of water and hexagonal ice. By comparing path-integral and classical simulations, we find that NQEs lead to a larger renormalization of the fundamental gap of ice, compared to that of water, ultimately yielding similar bandgaps in the two systems, consistent with experimental estimates. Our calculations suggest that the increased quantum mechanical delocalization of protons in ice, relative to water, is a key factor leading to the enhancement of NQEs on the electronic structure of ice.

Dynamics of the charge transfer to solvent process in aqueous iodide

Charge-transfer-to-solvent states in aqueous halides are ideal systems for studying the electron-transfer dynamics to the solvent involving a complex interplay between electronic excitation and solvent polarization. Despite extensive experimental investigations, a full picture of the charge-transfer-to-solvent dynamics has remained elusive. Here, we visualise the intricate interplay between the dynamics of the electron and the solvent polarization occurring in this process. Through the combined use of ab initio molecular dynamics and machine learning methods, we investigate the structure, dynamics and free energy as the excited electron evolves through the charge-transfer-to-solvent process, which we characterize as a sequence of states denoted charge-transfer-to-solvent, contact-pair, solvent-separated, and hydrated electron states, depending on the distance between the iodine and the excited electron. Our assignment of the charge-transfer-to-solvent states is supported by the good agreement between calculated and measured vertical binding energies. Our results reveal the charge transfer process in terms of the underlying atomic processes and mechanisms.

Partial Molar Solvation Volume of the Hydrated Electron Simulated Via DFT

Different simulation models of the hydrated electron produce different solvation structures, but it has been challenging to determine which simulated solvation structure, if any, is the most comparable to experiment. In a recent work, Neupane et al. [J. Phys. Chem. B 2023, 127, 5941-5947] showed using Kirkwood-Buff theory that the partial molar volume of the hydrated electron, which is known experimentally, can be readily computed from an integral over the simulated electron-water radial distribution function. This provides a sensitive way to directly compare the hydration structure of different simulation models of the hydrated electron with experiment. Here, we compute the partial molar volume of an ab-initio-simulated hydrated electron model based on density-functional theory (DFT) with a hybrid functional at different simulated system sizes. We find that the partial molar volume of the DFT-simulated hydrated electron is not converged with respect to the system size for simulations with up to 128 waters. We show that even at the largest simulation sizes, the partial molar volume of DFT-simulated hydrated electrons is underestimated by a factor of 2 with respect to experiment, and at the standard 64-water size commonly used in the literature, DFT-based simulations underestimate the experimental solvation volume by a factor of ∼3.5. An extrapolation to larger box sizes does predict the experimental partial molar volume correctly; however, larger system sizes than those explored here are currently intractable without the use of machine-learned potentials. These results bring into question what aspects of the predicted hydrated electron radial distribution function, as calculated by DFT-based simulations with the PBEh-D3 functional, deviate from the true solvation structure.

Absolute energy levels of liquid water from many-body perturbation theory with effective vertex corrections

We demonstrate the importance of addressing the Γ vertex and thus going beyond the GW approximation for achieving the energy levels of liquid water in many-body perturbation theory. In particular, we consider an effective vertex function in both the polarizability and the self-energy, which does not produce any computational overhead compared with the GW approximation. We yield the band gap, the ionization potential, and the electron affinity in good agreement with experiment and with a hybrid functional description. The achieved electronic structure and dielectric screening further lead to a good description of the optical absorption spectrum, as obtained through the solution of the Bethe-Salpeter equation. In particular, the experimental peak position of the exciton is accurately reproduced.

Tracking Cavity Formation in Electron Solvation: Insights from X-ray Spectroscopy and Theory

We present time-resolved X-ray absorption spectra of ionized liquid water and demonstrate that OH radicals, H3O+ ions, and solvated electrons all leave distinct X-ray-spectroscopic signatures. Particularly, this allows us to characterize the electron solvation process through a tool that focuses on the electronic response of oxygen atoms in the immediate vicinity of a solvated electron. Our experimental results, supported by ab initio calculations, confirm the formation of a cavity in which the solvated electron is trapped. We show that the solvation dynamics are governed by the magnitude of the random structural fluctuations present in water. As a consequence, the solvation time is highly sensitive to temperature and to the specific way the electron is injected into water.

Exploring the Unusual Reactivity of the Hydrated Electron with CO2

Many questions remain about the reactions of the hydrated electron despite decades of study. Of particular note is that they do not appear to follow the Marcus theory of electron transfer reactions, a feature that is yet to be explained. To investigate these issues, we used ab initio molecular dynamics (AIMD) simulations to investigate one of the better studied reactions, the hydrated electron reduction of CO2. The rate constant for the hydrated electron-CO2 reaction complex to react to form CO2- is for the first time estimated from AIMD simulations. Results at 298 and 373 K show the rate constant is insensitive to temperature, consistent with the low measured activation energy for the reaction, and the implications of this behavior are examined. The sampling provided by the simulations yields insight into the reaction mechanism. The reaction is found to involve both solvent reorganization and changes in the carbon dioxide structure. The latter leads to significant vibrational excitation of the bending and symmetric stretch vibrations in the CO2- product, indicating the reaction is vibrationally nonadiabatic. The former is estimated from the calculation of an approximate collective solvent coordinate and the free energy in this coordinate is determined. These results indicate that AIMD simulations can reasonably estimate hydrated electron reaction activation energies and provide new insight into the mechanism that can help illuminate the features of this unusual chemistry.

Multifunctional Roles of Clathrate Hydrate Nanoreactors for CO2 Reduction

In this study, we explore a possible platform for the CO2 reduction (CO2 R) in one of water's solid phases, namely clathrate hydrates (CHs), by ab initio molecular dynamics and well-tempered metadynamics simulations with periodic boundary conditions. We found that the stacked H2 O nanocages in CHs help to initialize CO2 R by increasing the electron-binding ability of CO2 . The substantial CO2 R processes are further influenced by the hydrogen bond networks in CHs. The first intermediate CO2 - in this process can be stabilized through cage structure reorganization into the H-bonded [CO2 - ⋅⋅⋅H-OHcage ] complex. Further cooperative structural dynamics enables the complex to convert into a vital transient [CO2 2- ⋅⋅⋅H-OHcage ] intermediate in a low-barrier disproportionation-like process. Such a highly reactive intermediate spontaneously triggers subsequent double proton transfer along its tethering H-bonds, finally converting it into HCOOH. These hydrogen-bonded nanoreactors feature multiple functions in facilitating CO2 R such as confining, tethering, H-bond catalyzing and proton pumping. Our findings have a general interest and extend the knowledge of CO2 R into porous aqueous systems.

Empirically Optimized One-Electron Pseudopotential for the Hydrated Electron: A Proof-of-Concept Study

Mixed quantum-classical molecular dynamics simulations have been important tools for studying the hydrated electron. They generally use a one-electron pseudopotential to describe the interactions of an electron with the water molecules. This approximation shows both the strength and weakness of the approach. On the one hand, it enables extensive statistical sampling and large system sizes that are not possible with more accurate ab initio molecular dynamics methods. On the other hand, there has (justifiably) been much debate about the ability of pseudopotentials to accurately and quantitatively describe the hydrated electron properties. These pseudopotentials have largely been derived by fitting them to ab initio calculations of an electron interacting with a single water molecule. In this paper, we present a proof-of-concept demonstration of an alternative approach in which the pseudopotential parameters are determined by optimizing them to reproduce key experimental properties. Specifically, we develop a new pseudopotential, using the existing TBOpt model as a starting point, which correctly describes the hydrated electron vertical detachment energy and radius of gyration. In addition to these properties, this empirically optimized model displays a significantly modified solvation structure, which improves, for example, the prediction of the partial molar volume.

Influence of H-bonds on the photoionization of aromatic chromophores in water: The aniline molecule

We have conducted time-resolved experiments (pump-probe and pump-repump-probe) on a model aromatic chromophore, aniline, after excitation in water at 267 nm. In the initial spectra recorded, in addition to the absorption corresponding to the bright ππ* excitation, the fingerprint of a transient state with the electron located on the solvent molecule is identified. We postulate that the latter corresponds to the πσ* state along the N-H bond, whose complete relaxation with a ∼500 ps lifetime results in the formation of the fully solvated electron and cation. This ionization process occurs in parallel with the ππ* photophysical channel that yields the characteristic ∼1 ns fluorescence lifetime. The observed branched pathway is rationalized in terms of the different H-bonds that the water establishes with the amino group. The proposed mechanism could be common for aromatics in water containing N-H or O-H bonds and would allow the formation of separated charges after excitation at the threshold of their electronic absorptions.

Many-body screening effects in liquid water

The screening arising from many-body excitations is a crucial quantity for describing absorption and inelastic X-ray scattering (IXS) of materials. Similarly, the electron screening plays a critical role in state-of-the-art approaches for determining the fundamental band gap. However, ab initio studies of the screening in liquid water have remained limited. Here, we use a combined analysis based on the Bethe-Salpeter equation and time-dependent density functional theory. We first show that absorption spectra at near-edge energies are insufficient to assess the accuracy by which the screening is described. Next, when the energy range under scrutiny is extended, we instead find that the IXS spectra are highly sensitive and allow for the selection of the optimal theoretical scheme. This leads to good agreement with experiment over a large range of transferred energies and momenta, and enables establishing the elusive fundamental band gap of liquid water at 9.3 eV.

How Ions Break Local Symmetry: Simulations of Polarized Transient Hole Burning for Different Models of the Hydrated Electron in Contact Pairs with Na. J Phys Chem Lett 2023;14:3014-3022. [PMID: 36943261 DOI: 10.1021/acs.jpclett.3c00220]

The hydrated electron (e_aq^-) is known via polarized transient hole-burning (pTHB) experiments to have a homogeneously broadened absorption spectrum. Here, we explore via quantum simulation how the pTHB spectroscopy of different e_aq^- models changes in the presence of electrolytes. The idea is that cation-e_aq^- pairing can break the local symmetry and, thus, induce persistent inhomogeneity. We find that a "hard" cavity model shows a modest increase in the pTHB recovery time in the presence of salt, while a "soft" cavity model remains homogeneously broadened independent of the salt concentration. We also explore the orientational anisotropy of a fully ab initio density functional theory-based model of the e_aq^-, which is strongly inhomogeneously broadened without salt and which becomes significantly more inhomogeneously broadened in the presence of salt. The results provide a direct prediction for experiments that can distinguish between different models and, thus, help pin down the hydration structure and dynamics of the e_aq^-.

Hydrated electrons as nodes in porous clathrate hydrates

We investigate the structures of hydrated electrons (e^- _aq) in one of water's solid phases, namely, clathrate hydrates (CHs). Using density functional theory (DFT) calculations, DFT-based ab initio molecular dynamics (AIMD), and path-integral AIMD simulations with periodic boundary conditions, we find that the structure of the e^- _aq@node model is in good agreement with the experiment, suggesting that an e^- _aq could form a node in CHs. The node is a H₂O defect in CHs that is supposed to be composed of four unsaturated hydrogen bonds. Since CHs are porous crystals that possess cavities that can accommodate small guest molecules, we expect that these guest molecules can be used to tailor the electronic structure of the e^- _aq@node, and it leads to experimentally observed optical absorption spectra of CHs. Our findings have a general interest and extend the knowledge of e^- _aq into porous aqueous systems.

Birth of the Hydrated Electron via Charge-Transfer-to-Solvent Excitation of Aqueous Iodide

A primary means to generate hydrated electrons in laboratory experiments is excitation to the charge-transfer-to-solvent (CTTS) state of a solute such as I^-(aq), but this initial step in the genesis of e^-(aq) has never been simulated directly using ab initio molecular dynamics. We report the first such simulations, combining ground- and excited-state simulations of I^-(aq) with a detailed analysis of fluctuations in the Coulomb potential experienced by the nascent solvated electron. What emerges is a two-step picture of the evolution of e^-(aq) starting from the CTTS state: I^-(aq) + hν → I^-*(aq) → I^•(aq) + e^-(aq). Notably, the equilibrated ground state of e^-(aq) evolves from I^-*(aq) without any nonadiabatic transitions, simply as a result of solvent reorganization. The methodology used here should be applicable to other photochemical electron transfer processes in solution, an important class of problems directly relevant to photocatalysis and energy transfer.

Ab Initio Studies of Hydrated Electron/Cation Contact Pairs: Hydrated Electrons Simulated with Density Functional Theory Are Too Kosmotropic

We have performed the first DFT-based ab initio MD simulations of a hydrated electron (e_aq^-) in the presence of Na⁺, a system chosen because ion-pairing behavior in water depends sensitively on the local hydration structure. Experiments show that e_aq^-'s interact weakly with Na⁺; the e_aq^-'s spectrum blue shifts by only a few tens of meV upon ion pairing without changing shape. We find that the spectrum of the DFT-simulated e_aq^- red shifts and changes shape upon interaction with Na⁺, in contrast with experiment. We show that this is because the hydration structure of the DFT-simulated e_aq^- is too ordered or kosmotropic. Conversely, simulations that produce e_aq^-'s with a less ordered or chaotropic hydration structure form weaker ion pairs with Na⁺, yielding predicted spectral blue shifts in better agreement with experiment. Thus, ab initio simulations based on hybrid GGA DFT functionals fail to produce the correct solvation structure for the hydrated electron.

Observation of a transient intermediate in the ultrafast relaxation dynamics of the excess electron in strong-field-ionized liquid water

A unified picture of the electronic relaxation dynamics of ionized liquid water has remained elusive despite decades of study. Here, we employ sub-two-cycle visible to short-wave infrared pump-probe spectroscopy and ab initio nonadiabatic molecular dynamics simulations to reveal that the excess electron injected into the conduction band (CB) of ionized liquid water undergoes sequential relaxation to the hydrated electron s ground state via an intermediate state, identified as the elusive p excited state. The measured CB and p-electron lifetimes are 0.26 ± 0.02 ps and 62 ± 10 fs, respectively. Ab initio quantum dynamics yield similar lifetimes and furthermore reveal vibrational modes that participate in the different stages of electronic relaxation, with initial relaxation within the dense CB manifold coupled to hindered translational motions whereas subsequent p-to-s relaxation facilitated by librational and even intramolecular bending modes of water. Finally, energetic considerations suggest that a hitherto unobserved trap state resides ~0.3-eV below the CB edge of liquid water. Our results provide a detailed atomistic picture of the electronic relaxation dynamics of ionized liquid water with unprecedented time resolution.

First-principles redox energy estimates under the condition of satisfying the general form of Koopmans’ theorem: An atomistic study of aqueous iron

In this work, we explore the relative accuracy to which a hybrid functional, in the context of density functional theory, may predict redox properties under the constraint of satisfying the general form of Koopmans’ theorem. Taking aqueous iron as our model system within the framework of first-principles molecular dynamics, direct comparison between computed single-particle energies and experimental ionization data is assessed by both (1) tuning the degree of hybrid exchange, to satisfy the general form of Koopmans’ theorem, and (2) ensuring the application of finite-size corrections. These finite-size corrections are benchmarked through classical molecular dynamics calculations, extended to large atomic ensembles, for which good convergence is obtained in the large supercell limit. Our first-principles findings indicate that while precise quantitative agreement with experimental ionization data cannot always be attained for solvated systems, when satisfying the general form of Koopmans’ theorem via hybrid functionals, theoretically robust estimates of single-particle redox energies are most often arrived at by employing a total energy difference approach. That is, when seeking to employ a value of exact exchange that does not satisfy the general form of Koopmans’ theorem, but some other physical metric, the single-particle energy estimate that would most closely align with the general form of Koopmans’ theorem is obtained from a total energy difference approach. In this respect, these findings provide important guidance for the more general comparison of redox energies computed via hybrid functionals with experimental data.

Investigation of the Failure of Marcus Theory for Hydrated Electron Reactions

Reactions of the hydrated electron with a wide variety of substrates have been found to exhibit unusually similar activation energies in a manner incompatible with Marcus electron transfer theory. Given the fundamental linear response assumption of Marcus theory, one possible explanation for this apparent failure is that the underlying free energy surfaces governing the reactions are not harmonic; i.e., hydrated electron structural fluctuations exhibit non-Gaussian behavior. In this work, we test this hypothesis by using simulations to calculate the hydrated electron vertical detachment energy distribution. We consider both cavity and noncavity models for the hydrated electron, between which the actual hydrated electron behavior is expected to lie. Our results identify a possible origin for non-Gaussian behavior of the hydrated electron but show that it is not of sufficient magnitude to explain the failure of Marcus theory to describe its reactions. Thus, other explanations must be sought.

Competitive Ion Pairing and the Role of Anions in the Behavior of Hydrated Electrons in Electrolytes

Experiments have shown that in the presence of electrolytes, the hydrated electron's absorption spectrum experiences a blue shift whose magnitude depends on both the salt concentration and chemical identity. Previous computer simulations have suggested that the spectral blue shift results from the formation of (cation, electron) contact pairs and that the concentration dependence arises because the number of cations simultaneously paired with the electron increases with increasing concentration. In this work, we perform new simulations to build an atomistic picture that explains the effect of salt identity on the observed hydrated electron spectral shifts. We simulate hydrated electrons in the presence of both monovalent (Na⁺) and divalent (Ca²⁺) cations paired with both Cl^- and a spherical species representing ClO₄^- anions. Our simulations reproduce the experimental observations that divalent ions produce larger blue shifts of the hydrated electron's spectrum than monovalent ions with the same anion and that perchlorate salts show enhanced blue shifts compared to chloride salts with the same cation. We find that these observations can be explained by competitive ion pairing. With small kosmotropic cations such as Na⁺ and Ca²⁺, aqueous chloride salts tend to form (cation, anion) contact pairs, whereas there is little ion pairing between these cations and chaotropic perchlorate anions. Hydrated electrons also strongly interact with these cations, but if the cations are also paired with anions, this affects the free energy of the electron-cation interaction. With chloride salts, hydrated electrons end up in complexes containing multiple cations plus a few anions as well as the electron. Repulsive interactions between the electron and the nearby Cl^- anions reduce the cation-induced spectral blue shift of the hydrated electron. With perchlorate salts, hydrated electrons pair with multiple cations without any associated anions, leading to the largest possible cation-induced spectral blue shift. We also see that the reason multivalent cations produce larger spectral blue shifts than monovalent cations is because hydrated electrons are able to simultaneously pair with a larger number of multivalent cations due to a larger free energy of interaction. Overall, the interaction of hydrated electrons with electrolytes fits well with the Hofmeister series, where the electron behaves as an anion that is slightly more able to break water's H-bond structure than chloride.

Trap-Seeking or Trap-Digging? Photoinjection of Hydrated Electrons into Aqueous NaCl Solutions

It is well-known that when excess electrons are injected into an aqueous solution, they localize and solvate in ∼1 ps. Still debated is whether localization occurs via "trap-digging", in which the electron carves out a suitable localization site, or by "trap-seeking", where the electron prefers to localize at pre-existing low-energy trap sites in solution. To distinguish between these two possible mechanisms, we study the localization dynamics of excess electrons in aqueous NaCl solutions using both ultrafast spectroscopy and mixed quantum-classical molecular dynamics simulations. By introducing pre-existing traps in the form of Na⁺ ions, we can use the cation-induced blue-shift of the hydrated electron's absorption spectrum to directly monitor the site of electron localization. Our experimental and computational results show that the electron prefers to localize directly at the sites of Na⁺ traps; the presence of concentrated electrolytes otherwise has little impact on the way trap-seeking hydrated electrons relax following injection.

Temperature Dependent Properties of the Aqueous Electron

The temperature‐dependent properties of the aqueous electron have been extensively studied using mixed quantum‐classical simulations in a wide range of thermodynamic conditions based on one‐electron pseudopotentials. While the cavity model appears to explain most of the physical properties of the aqueous electron, only a non‐cavity model has so far been successful in accounting for the temperature dependence of the absorption spectrum. Here, we present an accurate and efficient description of the aqueous electron under various thermodynamic conditions by combining hybrid functional‐based molecular dynamics, machine learning techniques, and multiple time‐step methods. Our advanced simulations accurately describe the temperature dependence of the absorption maximum in the presence of cavity formation. Specifically, our work reveals that the red shift of the absorption maximum results from an increasing gyration radius with temperature, rather than from global density variations as previously suggested.

Understanding the Temperature Dependence and Finite Size Effects in Ab Initio MD Simulations of the Hydrated Electron

The hydrated electron is of interest to both theorists and experimentalists as a paradigm solution-phase quantum system. Although the bulk of the theoretical work studying the hydrated electron is based on mixed quantum/classical (MQC) methods, recent advances in computer power have allowed several attempts to study this object using ab initio methods. The difficulty with employing ab initio methods for this system is that even with relatively inexpensive quantum chemistry methods such as density functional theory (DFT), such calculations are still limited to at most a few tens of water molecules and only a few picoseconds duration, leaving open the question as to whether the calculations are converged with respect to either system size or dynamical fluctuations. Moreover, the ab initio simulations of the hydrated electron that have been published to date have provided only limited analysis. Most works calculate the electron's vertical detachment energy, which can be compared to experiment, and occasionally the electronic absorption spectrum is also computed. Structural features, such as pair distribution functions, are rare in the literature, with the majority of the structural analysis being simple statements that the electron resides in a cavity, which are often based only on a small number of simulation snapshots. Importantly, there has been no ab initio work examining the temperature-dependent behavior of the hydrated electron, which has not been satisfactorily explained by MQC simulations. In this work, we attempt to remedy this situation by running DFT-based ab initio simulations of the hydrated electron as a function of both box size and temperature. We show that the calculated properties of the hydrated electron are not converged even with simulation sizes up to 128 water molecules and durations of several tens of picoseconds. The simulations show significant changes in the water coordination and solvation structure with box size. Our temperature-dependent simulations predict a red-shift of the absorption spectrum (computed using TD-DFT with an optimally tuned range-separated hybrid functional) with increasing temperature, but the magnitude of the predicted red-shift is larger than that observed experimentally, and the absolute position of the calculated spectra are off by over half an eV. The spectral red-shift at high temperatures is accompanied by both a partial loss of structure of the electron's central cavity and an increased radius of gyration that pushes electron density onto and beyond the first solvation shell. Overall, although ab initio simulations can provide some insights into the temperature-dependent behavior of the hydrated electron, the simulation sizes and level of quantum chemistry theory that are currently accessible are inadequate for correctly describing the experimental properties of this fascinating object.

Ultrafast proton transfer of the aqueous phenol radical cation

Proton transfer (PT) reactions are fundamental to numerous chemical and biological processes. While sub-picosecond PT involving electronically excited states has been extensively studied, little is known about ultrafast PT triggered by photoionization. Here, we employ femtosecond optical pump-probe spectroscopy and quantum dynamics calculations to investigate the ultrafast proton transfer dynamics of the aqueous phenol radical cation (PhOH˙⁺). Analysis of the vibrational wave packet dynamics reveals unusually short dephasing times of 0.18 ± 0.02 ps and 0.16 ± 0.02 ps for the PhOH˙⁺ O-H wag and bend frequencies, respectively, suggestive of ultrafast PT occurring on the ∼0.1 ps timescale. The reduced potential energy surface obtained from ab initio calculations shows that PT is barrierless when it is coupled to the intermolecular hindered translation between PhOH˙⁺ and the proton-acceptor water molecule. Quantum dynamics calculations yield a lifetime of 193 fs for PhOH˙⁺, in good agreement with the experimental results and consistent with the PT reaction being mediated by the intermolecular O⋯O stretch. These results suggest that photoionization can be harnessed to produce photoacids that undergo ultrafast PT. In addition, they also show that PT can serve as an ultrafast deactivation channel for limiting the oxidative damage potential of radical cations.

Hydrated Electrons in High-Concentration Electrolytes Interact with Multiple Cations: A Simulation Study

Experimental studies have demonstrated that the hydrated electron's absorption spectrum undergoes a concentration-dependent blue-shift in the presence of electrolytes such as NaCl. The blue-shift increases roughly linearly at low salt concentration but saturates as the solubility limit of the salt is approached. Previous attempts to understand the origin of the concentration-dependent spectral shift using molecular simulation have only examined the interaction between the hydrated electron and a single sodium cation, and these simulations predicted a spectral blue-shift that was an order of magnitude larger than that seen experimentally. Thus, in this paper, we first explore the reasons for the exaggerated spectral blue-shift when a simulated hydrated electron interacts with a single Na⁺. We find that the issue arises from nonpairwise additivity of the Na⁺-e^- and H₂O-e^- pseudopotentials used in the simulation. This effect arises because the solvating water molecules donate charge into the empty orbitals of Na⁺, lowering the effective charge of the cation and thus reducing the excess electron-cation interaction. Careful analysis shows, however, that although this nonpairwise additivity changes the energetics of the electron-Na⁺ interaction, the forces between the electron, Na⁺, and water are unaffected, so that mixed quantum/classical (MQC) simulations produce the correct structure and dynamics. With this in hand, we then use MQC simulations to explore the behavior of the hydrated electron as an explicit function of NaCl salt concentration. We find that the simulations correctly reproduce the observed experimental spectral shifting behavior. The reason for the spectral shift is that as the electrolyte concentration increases, the average number of cations simultaneously interacting in contact pairs with the hydrated electron increases from 1.0 at low concentrations to ∼2.5 near the saturation limit. As the number of cations that interact with the electron increases, the cation/electron interactions becomes slightly weaker, so that the corresponding Na⁺-e^- distance increases with increasing salt concentration. We also find that the dielectric constant of the solution plays little role in the observed spectroscopy, so that the electrolyte-dependent spectral shifts of the hydrated electron are directly related to the concentration-dependent number of closely interacting cations.

Field-Induced Electron Generation in Water: Solvation Dynamics and Many-Body Interactions

The solvated electron represents an elementary quantum system in a liquid environment. Electrons solvated in water have raised strong interest because of their prototypical properties, their role in radiation chemistry, and their relevance for charge separation and transport. Nonequilibrium dynamics of photogenerated electrons in water occur on ultrafast time scales and include charge transfer, localization, and energy dissipation processes. We present new insight into the role of fluctuating electric fields of the liquid for generating electrons in the presence of an external terahertz field and address polaronic many-body properties of solvated electrons. This Perspective combines a review of recent results from experiment and theory with a discussion of basic electric interactions of electrons in water.

Spectroscopic observation and ultrafast coherent vibrational dynamics of the aqueous phenylalanine radical

The phenylalanine radical (Phe˙) has been proposed to mediate biological electron transport (ET) and exhibit long-lived electronic coherences following attosecond photoionization. However, the coupling of ultrafast structural reorganization to the oxidation/ionization of biomolecules such as phenylalanine remains unexplored. Moreover, studies of ET involving Phe˙ are hindered by its hitherto unobserved electronic spectrum. Here, we report the spectroscopic observation and coherent vibrational dynamics of aqueous Phe˙, prepared by sub-6 fs photodetachment of phenylalaninate anions. Sub-picosecond transient absorption spectroscopy reveals the ultraviolet absorption signature of Phe˙. Ultrafast structural reorganization drives coherent vibrational motion involving nine fundamental frequencies and one overtone. DFT calculations rationalize the absence of the decarboxylation reaction, a photodegradation pathway previously identified for Phe˙. Our findings guide the interpretation of future attosecond experiments aimed at elucidating coherent electron motion in photoionized aqueous biomolecules and pave way for the spectroscopic identification of Phe˙ in studies of biological ET.

Shallow and deep trap states of solvated electrons in methanol and their formation, electronic excitation, and relaxation dynamics

We present condensed-phase first-principles molecular dynamics simulations to elucidate the presence of different electron trapping sites in liquid methanol and their roles in the formation, electronic transitions, and relaxation of solvated electrons (e_met⁻) in methanol. Excess electrons injected into liquid methanol are most likely trapped by methyl groups, but rapidly diffuse to more stable trapping sites with dangling OH bonds. After localization at the sites with one free OH bond (1OH trapping sites), reorientation of other methanol molecules increases the OH coordination number and the trap depth, and ultimately four OH bonds become coordinated with the excess electrons under thermal conditions. The simulation identified four distinct trapping states with different OH coordination numbers. The simulation results also revealed that electronic transitions of e_met⁻ are primarily due to charge transfer between electron trapping sites (cavities) formed by OH and methyl groups, and that these transitions differ from hydrogenic electronic transitions involving aqueous solvated electrons (e_aq⁻). Such charge transfer also explains the alkyl-chain-length dependence of the photoabsorption peak wavelength and the excited-state lifetime of solvated electrons in primary alcohols.

Condensed-phase first-principles molecular dynamics simulations elucidate the presence of different electron trapping sites in liquid methanol and their roles in the formation, electronic transitions, and relaxation of solvated electrons.

Flexible boundary layer using exchange for embedding theories. II. QM/MM dynamics of the hydrated electron

The FlexiBLE embedding method introduced in Paper I [Z. Shen and W. J. Glover, J. Chem. Phys. 155, 224112 (2021)] is applied to explore the structure and dynamics of the aqueous solvated electron at an all-electron density functional theory Quantum Mechanics/Molecular Mechanics level. Compared to a one-electron mixed quantum/classical description, we find the dynamics of the many-electron model of the hydrated electron exhibits enhanced coupling to water OH stretch modes. Natural bond orbital analysis reveals this coupling is due to significant population of water OH σ^* orbitals, reaching 20%. Based on this, we develop a minimal frontier orbital picture of the hydrated electron involving a cavity orbital and important coupling to 4-5 coordinating OH σ^* orbitals. Implications for the interpretation of the spectroscopy of this interesting species are discussed.

How Water-Ion Interactions Control the Formation of Hydrated Electron:Sodium Cation Contact Pairs

Although solvated electrons are a perennial subject of interest, relatively little attention has been paid to the way they behave in aqueous electrolytes. Experimentally, it is known that the hydrated electron's (e_aq^-) absorption spectrum shifts to the blue in the presence of salts, and the magnitude of the shift depends on the ion concentration and the identities of both the cation and anion. Does the blue-shift result from some type of dielectric effect from the bulk electrolyte, or are there specific interactions between the hydrated electron and ions in solution? Previous work has suggested that e_aq^- forms contact pairs with aqueous ions such as Na⁺, leading to the question of what controls the stability of such contact pairs and their possible connection to the observed spectroscopy. In this work, we use mixed quantum/classical simulations to examine the nature of Na⁺:e^- contact pairs in water, using a novel method for quantum umbrella sampling to construct e_aq^--ion potentials of mean force (PMF). We find that the nature of the contact pair PMF depends sensitively on the choice of the classical interactions used to describe the Na⁺-water interactions. When the ion-water interactions are slightly stronger, the corresponding cation:e^- contact pairs form at longer distances and become free energetically less stable. We show that this is because there is a delicate balance between solvation of the cation, solvation of e_aq^- and the direct electronic interaction between the cation and the electron, so that small changes in this balance lead to large changes in the formation and stability of e^--ion contact pairs. In particular, strengthening the ion-water interactions helps to maintain a favorable local solvation environment around Na⁺, which in turn forces water molecules in the first solvation shell of the cation to be unfavorably oriented toward the electron in a contact pair; stronger solvation of the cation also reduces the electronic overlap of e_aq^- with Na⁺. We also find that the calculated spectra of different models of Na⁺:e^- contact pairs do not shift monotonically with cation-electron distance, and that the calculated spectral shifts are about an order of magnitude larger than experiment, suggesting that isolated contact pairs are not the sole explanation for the blue-shift of the hydrated electron's spectrum in the presence of electrolytes.

Observation of intra- and intermolecular vibrational coherences of the aqueous tryptophan radical induced by photodetachment

The study of the photodetachment of amino acids in aqueous solution is pertinent to the understanding of elementary processes that follow the interaction of ionizing radiation with biological matter. In the case of tryptophan, the tryptophan radical that is produced by electron ejection also plays an important role in numerous redox reactions in biology, although studies of its ultrafast molecular dynamics are limited. Here, we employ femtosecond optical pump-probe spectroscopy to elucidate the ultrafast structural rearrangement dynamics that accompany the photodetachment of the aqueous tryptophan anion by intense, ∼5-fs laser pulses. The observed vibrational wave packet dynamics, in conjunction with density functional theory calculations, identify the vibrational modes of the tryptophan radical, which participate in structural rearrangement upon photodetachment. Aside from intramolecular vibrational modes, our results also point to the involvement of intermolecular modes that drive solvent reorganization about the N-H moiety of the indole sidechain. Our study offers new insight into the ultrafast molecular dynamics of ionized biomolecules and suggests that the present experimental approach can be extended to investigate the photoionization- or photodetachment-induced structural dynamics of larger biomolecules.

Ultrafast vibrational wave packet dynamics of the aqueous tyrosyl radical anion induced by photodetachment

The ultrafast dynamics triggered by the photodetachment of the tyrosinate dianion in aqueous environment shed light on the elementary processes that accompany the interaction of ionizing radiation with biological matter. Photodetachment of the tryosinate dianion yields the tyrosyl radical anion, an important intermediate in biological redox reactions, although the study of its ultrafast dynamics is limited. Here, we utilize femtosecond optical pump-probe spectroscopy to investigate the ultrafast structural reorganization dynamics that follow the photodetachment of the tyrosinate dianion in aqueous solution. Photodetachment of the tyrosinate dianion leads to vibrational wave packet motion along seven vibrational modes that are coupled to the photodetachment process. The vibrational modes are assigned with the aid of density functional theory (DFT) calculations. Our results offer a glimpse of the elementary dynamics of ionized biomolecules and suggest the possibility of extending this approach to investigate the ionization-induced structural rearrangement of other aromatic amino acids and larger biomolecules.

Electronic Structure of Water from Koopmans-Compliant Functionals

Obtaining a precise theoretical description of the spectral properties of liquid water poses challenges for both molecular dynamics (MD) and electronic structure methods. The lower computational cost of the Koopmans-compliant functionals with respect to Green's function methods allows the simulations of many MD trajectories, with a description close to the state-of-art quasi-particle self-consistent GW plus vertex corrections method (QSGW + f_xc). Thus, we explore water spectral properties when different MD approaches are used, ranging from classical MD to first-principles MD, and including nuclear quantum effects. We have observed that different MD approaches lead to up to 1 eV change in the average band gap; thus, we focused on the band gap dependence with the geometrical properties of a system to explain such spread. We have evaluated the changes in the band gap due to variations in the intramolecular O-H bond distance and HOH angle, as well as the intermolecular hydrogen bond O···O distance and the OHO angles. We have observed that the dominant contribution comes from the O-H bond length; the O···O distance plays a secondary role, and the other geometrical properties do not significantly influence the gap. Furthermore, we analyze the electronic density of states (DOS), where the KIPZ functional shows good agreement with the DOS obtained with state-of-art approaches employing quasi-particle self-consistent GW plus vertex corrections. The O-H bond length also significantly influences the DOS. When nuclear quantum effects are considered, broadening of the peaks driven by the broader distribution of the O-H bond lengths is observed, leading to a closer agreement with the experimental photoemission spectra.

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low-energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. Intensity borrowing by these spurious states can affect intensities of the valence excitations, altering electronic bandshapes. At higher excitation energies, it is difficult to distinguish spurious charge-transfer states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces "natural charge-transfer orbitals" that provide a means to isolate orbitals that are most likely to participate in a CTTS excitation. Projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation and be used as an automated subspace selection algorithm for projection-based embedding of a high-level description of the CTTS state in a lower-level description of its environment. We apply this method to an ab initio molecular dynamics trajectory of I^-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with the experiment, and only one-third of the water molecules in the I^-(H₂O)₉₆ simulation cell need to be described with LR-TDDFT to obtain excitation energies that are converged to <0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments.

Simulating the ghost: quantum dynamics of the solvated electron

The nature of the bulk hydrated electron has been a challenge for both experiment and theory due to its short lifetime and high reactivity, and the need for a high-level of electronic structure theory to achieve predictive accuracy. The lack of a classical atomistic structural formula makes it exceedingly difficult to model the solvated electron using conventional empirical force fields, which describe the system in terms of interactions between point particles associated with atomic nuclei. Here we overcome this problem using a machine-learning model, that is sufficiently flexible to describe the effect of the excess electron on the structure of the surrounding water, without including the electron in the model explicitly. The resulting potential is not only able to reproduce the stable cavity structure but also recovers the correct localization dynamics that follow the injection of an electron in neat water. The machine learning model achieves the accuracy of the state-of-the-art correlated wave function method it is trained on. It is sufficiently inexpensive to afford a full quantum statistical and dynamical description and allows us to achieve accurate determination of the structure, diffusion mechanisms, and vibrational spectroscopy of the solvated electron.

Molecular Dynamics Characterization of Dielectron Hydration in Liquid Water with Unique Double Proton Transfers

Radiation chemistry of water and aqueous solutions has always been an interesting scientific issue owing to involving electronic excitations, ionization of solvated species, and formation of radiolytic species and many elementary reactions, but the underlying mechanisms are still poorly understood. Here, we for the first time molecular dynamics characterize the hydration dynamics of two correlated electrons and their triggered unique phenomena in liquid water associated with radiolysis of water using the combined hybrid functional and nonlocal dispersion functional. Hydration of two electrons may experience two distinctly different mechanisms, one forming a spin-paired closed-shell unicaged dielectron hydrate (e₂^2-_aq) and the other forming a spin-paired metastable open-shell bicaged hydrated electron pair (e^-_aq···e^-_aq) which exhibits intriguing antiferromagnetic spin coupling dynamics (in a range of -40 cm^-1 to -500 cm^-1). e^-_aq···e^-_aq can recombine to e₂^2-_aq through a unique solvent fluctuation-controlled gradual-flowing mechanism, and enlarging fluctuation can promote the conversion. Interestingly, we directly observe that e₂^2-_aq as the precursor can trigger hydrogen evolution via unique continuous spontaneous double proton transfer to the dielectron with a short-lived H^-_aq intermediate, but e^-_aq···e^-_aq does not directly. This is the first direct observation for the connection between e₂^2-_aq and spontaneous hydrogen evolution including participation of H^-_aq in aqueous solution, bridging relevant experimental phenomena. This work also evidences an unnoticed process, the double proton transfer mediated charge separation, and presents the first detailed analysis regarding the evolution dynamics of e₂^2-_aq for the understanding of the radiolysis reactions in aqueous solutions.

Mechanism of Aqueous Carbon Dioxide Reduction by the Solvated Electron

Aqueous solvated electron (e_aq^-), a key species in radiation and plasma chemistry, can efficiently reduce CO₂ in a potential green chemistry application. Here, the mechanism of this reaction is unravelled by condensed-phase molecular dynamics based on the correlated wave function and an accurate density functional theory (DFT) approximation. Here, we design and apply the holistic protocol for solvated electron's reactions encompassing all relevant reaction stages starting from diffusion. The carbon dioxide reduction proceeds via a cavity intermediate, which is separated from the product (CO₂^-) by an energy barrier due to the bending of CO₂ and the corresponding solvent reorganization energy. The formation of the intermediate is caused by solvated electron's diffusion, whereas the intermediate transformation to CO₂^- is triggered by hydrogen bond breaking in the second solvation shell of the solvated electron. This picture of an activation-controlled e_aq^- reaction is very different from both rapid barrierless electron transfer and proton-coupled electron transfer, where key transformations are caused by proton migration.

Evaluating Simple Ab Initio Models of the Hydrated Electron: The Role of Dynamical Fluctuations

Despite its importance in electron transfer reactions and radiation chemistry, there has been disagreement over the fundamental nature of the hydrated electron, such as whether or not it resides in a cavity. Mixed quantum/classical simulations of the hydrated electron give different structures depending on the pseudopotential employed, and ab initio models of computational necessity use small numbers of water molecules and/or provide insufficient statistics to compare to experimental observables. A few years ago, Kumar et al. (J. Phys. Chem. A 2015, 119, 9148) proposed a minimalist ab initio model of the hydrated electron with only a small number of explicitly treated water molecules plus a polarizable continuum model (PCM). They found that the optimized geometry had four waters arranged tetrahedrally around a central cavity, and that the calculated vertical detachment energy and radius of gyration agreed well with experiment, results that were largely independent of the level of theory employed. The model, however, is based on a fixed structure at 0 K and does not explicitly incorporate entropic contributions or the thermal fluctuations that should be associated with the room-temperature hydrated electron. Thus, in this paper, we extend the model of Kumar et al. by running Born-Oppenheimer molecular dynamics (BOMD) of a small number of water molecules with an excess electron plus PCM at room temperature. We find that when thermal fluctuations are introduced, the level of theory chosen becomes critical enough when only four waters are used that one of the waters dissociates from the cluster with certain density functionals. Moreover, even with an optimally tuned range-separated hybrid functional, at room temperature the tetrahedral orientation of the 0 K first-shell waters is entirely lost and the central cavity collapses, a process driven by the fact that the explicit water molecules prefer to make H-bonds with each other more than with the excess electron. The resulting average structure is quite similar to that produced by a noncavity mixed quantum/classical model, so that the minimalist 4-water BOMD models suffer from problems similar to those of noncavity models, such as predicting the wrong sign of the hydrated electron's molar solvation volume. We also performed BOMD with 16 explicit water molecules plus an extra electron and PCM. We find that the inclusion of an entire second solvation shell of explicit water leads to little change in the outcome from when only four waters were used. In fact, the 16-water simulations behave much like those of water cluster anions, in which the electron localizes at the cluster surface, showing that PCM is not acceptable for use in minimalist models to describe the behavior of the bulk hydrated electron. For both the 4- and 16-water models, we investigate how the introduction of thermal motions alters the predicted absorption spectrum, vertical detachment energy, and resonance Raman spectrum of the simulated hydrated electron. We also present a set of structural criteria that can be used to numerically determine how cavity-like (or not) a particular hydrated electron model is. All of the results emphasize that the hydrated electron is a statistical object whose properties are inadequately captured using only a small number of explicit waters, and that a proper treatment of thermal fluctuations is critical to understanding the hydrated electron's chemical and physical behavior.

Probing the Birth and Ultrafast Dynamics of Hydrated Electrons at the Gold/Liquid Water Interface via an Optoelectronic Approach

The hydrated electron has fundamental and practical significance in radiation and radical chemistry, catalysis, and radiobiology. While its bulk properties have been extensively studied, its behavior at solid/liquid interfaces is still unclear due to the lack of effective tools to characterize this short-lived species in between two condensed matter layers. In this study, we develop a novel optoelectronic technique for the characterization of the birth and structural evolution of solvated electrons at the metal/liquid interface with a femtosecond time resolution. Using this tool, we record for the first time the transient spectra (in a photon energy range from 0.31 to 1.85 eV) in situ with a time resolution of 50 fs revealing several novel aspects of their properties at the interface. Especially the transient species show state-dependent optical transition behaviors from being isotropic in the hot state to perpendicular to the surface in the trapped and solvated states. The technique will enable a better understanding of hot electron driven reactions at electrochemical interfaces.

Can Current Experimental Data Exclude Non-Gaussian Genuine Band Shapes in Ultraviolet Photoelectron Spectra of the Hydrated Electron?

Two recent articles present results that allegedly exclude a possible multimodal distribution of the hydrated electron in ultraviolet photoelectron spectra. The first article bases its conclusion on the assumption that the non-Gaussian genuine band shape previously retrieved for the solvated electron in liquid water is an artifact arising from insufficient electron scattering cross sections used in the retrieval. The second article excludes a multimodal band shape based on a photoelectron spectrum of the solvated electron in water clusters recorded at a single ultraviolet photon energy, and it further assumes that cluster results are transferable to the liquid without further justification. Here, we show that based on current data multimodal distributions cannot be unambiguously excluded. Furthermore, the transferability of cluster results to the liquid can be neither justified nor refuted on the basis of currently available experimental ultraviolet photoelectron spectra.

The Fluxional Nature of the Hydrated Electron: Energy and Entropy Contributions to Aqueous Electron Free Energies

There has been a great deal of recent controversy over the structure of the hydrated electron and whether it occupies a cavity or contains a significant number of interior waters (noncavity). The questions we address in this work are, from a free energy perspective, how different are these proposed structures? Do the different structures all lie along a single continuum, or are there significant differences (i.e., free energy barriers) between them? To address these questions, we have performed a series of one-electron calculations using umbrella sampling with quantum biased molecular dynamics along a coordinate that directly reflects the number of water molecules in the hydrated electron's interior. We verify that a standard cavity model of the hydrated electron behaves essentially as a hard sphere: the model is dominated by repulsion at short range such that water is expelled from a local volume around the electron, leading to a water solvation shell like that of a pseudohalide ion. The repulsion is much larger than thermal energies near room temperature, explaining why such models exhibit properties with little temperature dependence. On the other hand, our calculations reveal that a noncavity model is highly fluxional, meaning that thermal motions cause the number of interior waters to fluctuate from effectively zero (i.e., a cavity-type electron) to potentially above the bulk water density. The energetic contributions in the noncavity model are still repulsive in the sense that they favor cavity formation, so the fluctuations in structure are driven largely by entropy: the entropic cost for expelling water from a region of space is large enough that some water is still driven into the electron's interior. As the temperature is lowered and entropy becomes less important, the noncavity electron's structure is predicted to become more cavity-like, consistent with the observed temperature dependence of the hydrated electron's properties. Thus, we argue that although the specific noncavity model we study overestimates the preponderance of fluctuations involving interior water molecules, with appropriate refinements to correctly capture the true average number of interior waters and molar solvation volume, a fluxional model likely makes the most sense for understanding the various experimental properties of the hydrated electron.

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.

Reaction of Electrons with DNA: Radiation Damage to Radiosensitization

This review article provides a concise overview of electron involvement in DNA radiation damage. The review begins with the various states of radiation-produced electrons: Secondary electrons (SE), low energy electrons (LEE), electrons at near zero kinetic energy in water (quasi-free electrons, (e-qf)) electrons in the process of solvation in water (presolvated electrons, e-pre), and fully solvated electrons (e-aq). A current summary of the structure of e-aq, and its reactions with DNA-model systems is presented. Theoretical works on reduction potentials of DNA-bases were found to be in agreement with experiments. This review points out the proposed role of LEE-induced frank DNA-strand breaks in ion-beam irradiated DNA. The final section presents radiation-produced electron-mediated site-specific formation of oxidative neutral aminyl radicals from azidonucleosides and the evidence of radiosensitization provided by these aminyl radicals in azidonucleoside-incorporated breast cancer cells.

Ultrafast structural rearrangement dynamics induced by the photodetachment of phenoxide in aqueous solution

The elementary processes that accompany the interaction of ionizing radiation with biologically relevant molecules are of fundamental importance. However, the ultrafast structural rearrangement dynamics induced by the ionization of biomolecules in aqueous solution remain hitherto unknown. Here, we employ femtosecond optical pump-probe spectroscopy to elucidate the vibrational wave packet dynamics that follow the photodetachment of phenoxide, a structural mimic of tyrosine, in aqueous solution. Photodetachment of phenoxide leads to wave packet dynamics of the phenoxyl radical along 12 different vibrational modes. Eight of the modes are totally symmetric and support structural rearrangement upon electron ejection. Comparison to a previous photodetachment study of phenoxide in the gas phase reveals the important role played by the solvent environment in driving ultrafast structural reorganization induced by ionizing radiation. This work provides insight into the ultrafast molecular dynamics that follow the interaction of ionizing radiation with molecules in aqueous solution.

The interaction of biomolecules with ionizing radiation induces structural changes which are still largely unknown. The authors use femtosecond wave packet spectroscopy to observe ultrafast structural dynamics that follow the photodetachment of phenoxide in aqueous solution.

Picture of the wet electron: a localized transient state in liquid water

A transient state of the excess electron in liquid water preceding the development of the solvation shell, the so-called wet electron, has been invoked to explain spectroscopic observations, but its binding energy and atomic structure have remained highly elusive. Here, we carry out hybrid functional molecular dynamics to unveil the ultrafast solvation mechanism leading to the hydrated electron. In the pre-hydrated regime, the electron is found to repeatedly switch between a quasi-free electron state in the conduction band and a localized state with a binding energy of 0.26 eV, which we assign to the wet electron. This transient state self-traps in a region of the liquid which extends up to ∼4.5 Å and involves a severe disruption of the hydrogen-bond network. Our picture provides an unprecedented view on the nature of the wet electron, which is instrumental to understanding the properties of this fundamental species in liquid water.

Dynamics of the Bulk Hydrated Electron from Many-Body Wave-Function Theory

The structure of the hydrated electron is a matter of debate as it evades direct experimental observation owing to the short life time and low concentrations of the species. Herein, the first molecular dynamics simulation of the bulk hydrated electron based on correlated wave‐function theory provides conclusive evidence in favor of a persistent tetrahedral cavity made up by four water molecules, and against the existence of stable non‐cavity structures. Such a cavity is formed within less than a picosecond after the addition of an excess electron to neat liquid water, with less regular cavities appearing as intermediates. The cavities are bound together by weak H−H bonds, the number of which correlates well with the number of coordinated water molecules, each type of cavity leaving a distinct spectroscopic signature. Simulations predict regions of negative spin density and a gyration radius that are both in agreement with experimental data.