Short-Range Electron Correlation Stabilizes Noncavity Solvation of the Hydrated Electron

2-in-1 Phase Space Sampling for Calculating the Absorption Spectrum of the Hydrated Electron

The investigation of vibrational effects on absorption spectrum calculations often employs Wigner sampling or thermal sampling. While Wigner sampling incorporates zero-point energy, it may not be suitable for flexible systems. Thermal sampling is applicable to anharmonic systems yet treats nuclei classically. The application of generalized smoothed trajectory analysis (GSTA) as a postprocessing method allows for the incorporation of nuclear quantum effects (NQEs), combining the advantages of both sampling methods. We demonstrate this approach in computing the absorption spectrum of a hydrated electron. Theoretical exploration of the hydrated electron and its embryonic forms, such as water cluster anions, poses a significant challenge due to the diffusivity of the excess electron and the continuous motion of water molecules. In many previous studies, the wave nature of atomic nuclei is often neglected, despite the substantial impact of NQEs on thermodynamic and spectroscopic properties, particularly for hydrogen atoms. In our studies, we examine these NQEs for the excess electrons in various water systems. We obtained structures from mixed classical-quantum simulations for water cluster anions and the hydrated electron by incorporating the quantum effects of atomic nuclei with the filtration of the classical trajectories. Absorption spectra were determined at different theoretical levels. Our results indicate significant NQEs, red shift, and broadening of the spectra for hydrated electron systems. This study demonstrates the applicability of GSTA to complex systems, providing insights into NQEs on energetic and structural properties.

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.

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.

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^-.

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.

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.

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.

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.

Flexible boundary layer using exchange for embedding theories. I. Theory and implementation

Embedding theory is a powerful computational chemistry approach to exploring the electronic structure and dynamics of complex systems, with Quantum Mechanical/Molecular Mechanics (QM/MM) being the prime example. A challenge arises when trying to apply embedding methodology to systems with diffusible particles, e.g., solvents, if some of them must be included in the QM region, for example, in the description of solvent-supported electronic states or reactions involving proton transfer or charge-transfer-to-solvent: without a special treatment, inter-diffusion of QM and MM particles will eventually lead to a loss of QM/MM separation. We have developed a new method called Flexible Boundary Layer using Exchange (FlexiBLE) that solves the problem by adding a biasing potential to the system that closely maintains QM/MM separation. The method rigorously preserves ensemble averages by leveraging their invariance to an exchange of identical particles. With a careful choice of the biasing potential and the use of a tree algorithm to include only important QM and MM exchanges, we find that the method has an MM-forcefield-like computational cost and thus adds negligible overhead to a QM/MM simulation. Furthermore, we show that molecular dynamics with the FlexiBLE bias conserves total energy, and remarkably, sub-diffusional dynamical quantities in the inner QM region are unaffected by the applied bias. FlexiBLE thus widens the range of chemistry that can be studied with embedding theory.

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.

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.

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.

Unprecedented O:⇔:O compression and H↔H fragilization in Lewis solutions

Charge injection in terms of protons, lone pairs, cations and anions by acid and base solvation mediates the HB network and properties of Lewis solutions through H↔H fragilization, O:⇔:O compression and polarization, ionic polarization and hydrating H₂O dipolar screen shielding, anion–anion repulsion, compressed solvent H–O bond elongation and undercoordinated solute H–O bond contraction.

Structure of the aqueous electron

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

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.

Electronic Levels of Excess Electrons in Liquid Water

We provide a consistent description of the electronic levels associated with localized and delocalized excess electrons in liquid water by combining hybrid-functional molecular dynamics simulations with a grand canonical formulation of solutes in aqueous solution. The excess electron localizes in a cavity with an average radius of 1.8 Å and a majority coordination of five water molecules. The vertical binding energy, the optical s-p transitions, and the adiabatic redox level are found to agree closely with their experimental counterparts. The energy level associated with electron delocalization V₀ is inferred to lie at -0.97 eV with respect to the vacuum level.