Nature's most squishy ion: The important role of solvent polarization in the description of the hydrated electron

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.

Evidence for Water Antibonding Orbital Mixing in the Hydrated Electron from Its Oxygen 1s X-ray Absorption Spectrum

The X-ray absorption spectrum (XAS) of the hydrated electron (e_(aq)^-) has been simulated using time-dependent density functional theory with a quantum mechanics/molecular mechanics description. A unique XAS peak at 533 eV is observed with an energy and intensity in quantitative agreement with recent time-resolved experiments, allowing its assignment as arising from water O1s transitions to the singly occupied molecular orbital (SOMO) in which the excess electron resides. The transitions acquire oscillator strength due to the SOMO comprising an admixture of a cavity-localized orbital and water 4a₁ and 2b₂ antibonding orbitals. The mixing of antibonding orbitals has implications for the strength of couplings between e_(aq)^- and intramolecular modes of water.

Coherent polaron dynamics of electrons solvated in polar liquids

An electron solvated in a polar liquid is an elementary quantum system with properties governed by electric interactions with a fluctuating molecular environment. In the prevailing single particle picture, the quantum ground and excited states are determined by a self-consistent potential, as defined by the particular local configuration of the solvation shell. This description neglects collective many-body excitations, which arise from the coupling of electronic degrees of freedom and nuclear motions of the environment. While recent experiments have demonstrated collective nonequilbrium electronic-nuclear motion, i.e. polaron excitations in liquid water, their relevance in the broader context of polar liquids has remained unexplored. Here, we study the nonequilibrium dielectric response of the, compared to water, less polar alcohols isopropanol, and ethylene glycol, that also display a different hydrogen bond pattern. We demonstrate that ultrafast relaxation of photogenerated electrons impulsively induces coherent charge oscillations, which persist for some 10 ps. They emit electric waves in a frequency range from 0.1 to 2 THz, depending on electron concentration. Oscillation frequencies and line shapes are reproduced by a unified polaron picture for alcohols and water, which is based on a Clausius-Mossotti local field approach for the THz dielectric function. The analysis suggests a longitudinal character of many-body polaron excitations and a weak coupling to transverse excitations, supported by the underdamped character of charge oscillations. Polaron dynamics are governed by the long-range Coulomb interaction between an excess electron and several thousands of polar solvent molecules, while local electron solvation geometries play a minor role.

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

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.

Partial Molar Volume of the Hydrated Electron

The partial molar volume of the hydrated electron was investigated with pulse radiolysis and transient absorption by measuring the pressure dependence of the equilibrium constant for e^-_aq + NH₄⁺ ⇔ H + NH₃. At 2 kbar pressure, the equilibrium constant decreases relative to 1 bar by only 6%. Using tabulated molar volumes for ammonia and ammonium, we have the result V̅(e^-_aq) - V̅(H) = 11.3 cm³/mol at 25 °C, confirming that V̅(e^-_aq) is positive and even larger than the hydrophobic H atom. Assuming on the basis of recent molecular dynamics simulations that the molar volume of the H atom is somewhat less than that of H₂, we estimate V̅(e^-_aq) = 26 ± 6 cm³/mol. The positive molar volume is consistent with an electron that exists largely in a small solvent void (cavity), ruling out a recent model ( Larsen , R. E. ; Glover , W. J. ; Schwartz , B. J. Science 2010 , 329 , 65 - 69 ) that suggests a noncavity structure with negative molar volume.

Selectivity in Electron Attachment to Water Clusters

Electron attachment onto water clusters to form water cluster anions is studied by varying the point of electron attachment along a molecular beam axis and probing the produced cluster anions using photoelectron spectroscopy. The results show that the point of electron attachment has a clear effect on the final distribution of isomers for a cluster containing 78 water molecules, with isomer I formed preferentially near the start of the expansion and isomer II formed preferentially once the molecular beam has progressed for several millimeters. These changes can be accounted for by the cluster growth rate along the beam. Near the start of the expansion, cluster growth is proceeding rapidly with condensing water molecules solvating the electron, while further along the expansion, the growth has terminated and electrons are attached to large and cold preformed clusters, leading to the isomer associated with a loosely bound surface state.

Structure of the aqueous electron

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

Prediction of rate constants of important chemical reactions in water radiation chemistry in sub and supercritical water – non-equilibrium reactions

The rate constants for reactions involved in the radiolysis of water under relevant thermodynamic conditions in supercritical water-cooled reactors are estimated for inputs in simulations of the radiation chemistry in Generation IV nuclear reactors. We have discussed the mechanism of each chemical reaction with a focus on non-equilibrium reactions. We found most of the reactions are activation controlled above the critical point and that the rate constants are not significantly pressure dependent below 300 °C. This work will aid industry with developing chemical control strategies to suppress the concentration of eroding species.

A Simple ab Initio Model for the Solvated Electron in Methanol

The solvation structure of a solvated electron in methanol is investigated with ab initio calculations of small anion methanol clusters in a polarized dielectric continuum. We find that the lowest-energy structure in best agreement with experiment, calculated with CCSD, MP2, and B3LYP methods with aug-cc-pvdz basis set, is a tetrahedral arrangement of four methanol molecules with OH bonds oriented toward the center. The optimum distance from the tetrahedron center to the hydroxyl protons is ∼1.8 Å, significantly smaller than previous estimates. We are able to reproduce experimental radius of gyration Rg (deduced from optical absorption), vertical detachment energy, and resonance Raman frequencies. The electron paramagnetic resonance g-factor shift is qualitatively reproduced using density functional theory.

Hydrated Electrons in Water Clusters: Inside or Outside, Cavity or Noncavity?

In this work, we compare the applicability of three electron–water molecule pseudopotentials in modeling the physical properties of hydrated electrons. Quantum model calculations illustrate that the recently suggested Larsen–Glover–Schwartz (LGS) model and its modified m-LGS version have a too-attractive potential in the vicinity of the oxygen. As a result, LGS models predict a noncavity hydrated electron structure in clusters at room temperature, as seen from mixed one-electron quantum–classical molecular dynamics simulations of water cluster anions, with the electron localizing exclusively in the interior of the clusters. Comparative calculations using the cavity-preferring Turi–Borgis (TB) model predict interior-state and surface-state cluster isomers. The computed associated physical properties are also analyzed and compared to available experimental data. We find that the LGS and m-LGS potentials provide results that appear to be inconsistent with the size dependence of the experimental data. The simulated TB tendencies are qualitatively correct. Furthermore, ab initio calculations on static LGS noncavity structures indicate weak stabilization of the excess electron in regions where the LGS potential preferably and strongly binds the electron. TB calculations give stabilization energies that are in line with the ab initio results. In conclusion, we observe that the cavity-preferring pseudopotential model predicts cluster physical properties in better agreement with experimental data and ab initio calculations than the models predicting noncavity structures for the hydrated electron.

On the nature of the solvated electron in ice Ih

The water-solvated excess electron (EE) is a key chemical agent whose hallmark signature, its asymmetric optical absorption spectrum, continues to be a topic of debate.

Reactivity of Hydrated Electron in Finite Size System: Sodium Pickup on Mixed N2O-Water Nanoparticles

We investigate the reactivity of hydrated electron generated by alkali metal deposition on small water particles with nitrous oxide dopant by means of mass spectrometry and ab initio molecular dynamics simulations. The mixed nitrous oxide/water clusters were generated in a molecular beam and doped with Na atoms in a pickup experiment, and investigated by mass spectrometry using two different ionization schemes: an electron ionization (EI), and UV photoionization after the Na doping (NaPI). The NaPI is a soft-ionization nondestructive method, especially for water clusters provided that a hydrated electron es– is formed in the cluster. The missing signal for the doped clusters indicates that the hydrated electron is not present in the N2O containing clusters. The simulations reveal that the hydrated electron is formed, but it immediately reacts with N2O, forming first N2O– radical anion, later O–, and finally an OH• and OH– pair.

Quantifying solvated electrons' delocalization

The electron delocalization range EDR(r;u_av) (left) captures the spin density (right) of an electron delocalized over u_av = 5.77 Å on the surface of an (H₂O)₂₀⁻ cluster.

Glucose-Promoted Localization Dynamics of Excess Electrons in Aqueous Glucose Solution Revealed by Ab Initio Molecular Dynamics Simulation

Ab initio molecular dynamics simulations reveal that an excess electron (EE) can be more efficiently localized as a cavity-shaped state in aqueous glucose solution (AGS) than in water. Compared with that (∼1.5 ps) in water, the localization time is shortened by ∼0.7-1.2 ps in three AGSs (0.56, 1.12, and 2.87 M). Although the radii of gyration of the solvated EEs are all close to 2.6 Å in the four solutions, the solvated EE cavities in the AGSs become more compact and can localize ∼80% of an EE, which is considerably larger than that (∼40-60% and occasionally ∼80%) in water. These observations are attributed to a modification of the hydrogen-bonded network by the introduction of glucose molecules into water. The water acts as a promoter and stabilizer, by forming voids around glucose molecules and, in this fashion, favoring the localization of an EE with high efficiency. This study provides important information about EEs in physiological AGSs and suggests a new strategy to efficiently localize an EE in a stable cavity for further exploration of biological function.

Accurate first principles model potentials for intermolecular interactions

The general effective fragment potential (EFP) method provides model potentials for any molecule that is derived from first principles, with no empirically fitted parameters. The EFP method has been interfaced with most currently used ab initio single-reference and multireference quantum mechanics (QM) methods, ranging from Hartree-Fock and coupled cluster theory to multireference perturbation theory. The most recent innovations in the EFP model have been to make the computationally expensive charge transfer term much more efficient and to interface the general EFP dispersion and exchange repulsion interactions with QM methods. Following a summary of the method and its implementation in generally available computer programs, these most recent new developments are discussed.

Visualization of chemical reaction dynamics: toward understanding complex polyatomic reactions

Polyatomic molecules have several electronic states that have similar energies. Consequently, their chemical dynamics often involve nonadiabatic transitions between multiple potential energy surfaces. Elucidating the complex reactions of polyatomic molecules is one of the most important tasks of theoretical and experimental studies of chemical dynamics. This paper describes our recent experimental studies of the multidimensional multisurface dynamics of polyatomic molecules based on two-dimensional ion/electron imaging. It also discusses ultrafast photoelectron spectroscopy of liquids for elucidating nonadiabatic electronic dynamics in aqueous solutions. (Communicated by Hiroo INOKUCHI, M.J.A.)

Unraveling the Complex Nature of the Hydrated Electron

The structure of the hydrated electron, which is a key species in radiative processes in water, has remained elusive. The traditional cavity model has been questioned recently, but the newly suggested picture of an electron delocalized over a region of enhanced water density is controversial. Here, we present results from ab initio molecular dynamics simulations, where not only the excess electron but also the valence electrons of the surrounding water molecules are described quantum mechanically. Unlike in previous one-electron pseudopotential calculations, many-electron interactions are explicitly accounted for. The present approach allows for partitioning of the electron solvated in liquid water into contributions from an inner cavity, neighboring water molecules, and a diffuse tail. We demonstrate that all three of these contributions are sizable and, consequently, important, which underlines the complex nature of the hydrated electron and warns against oversimplified interpretations based on pseudopotential models.

Structure of the aqueous electron: assessment of one-electron pseudopotential models in comparison to experimental data and time-dependent density functional theory

The prevailing structural paradigm for the aqueous electron is that of an s-like ground-state wave function that inhabits a quasi-spherical solvent cavity, a viewpoint that is supported by numerous atomistic simulations using various one-electron pseudopotential models. This conceptual picture has recently been challenged, however, on the basis of results obtained from a new electron-water pseudopotential model that predicts a more delocalized wave function and no well-defined solvent cavity. Here, we examine this new model in comparison to two alternative, cavity-forming pseudopotential models. We find that the cavity-forming models are far more consistent with the experimental data for the electron's radius of gyration, optical absorption spectrum, and vertical electron binding energy. Calculations of the absorption spectrum using time-dependent density functional theory are in quantitative or semiquantitative agreement with experiment when the solvent geometries are obtained from the cavity-forming pseudopotential models, but differ markedly from experiment when geometries that do not form a cavity are used.

Effects of microviscosity, dry electron scavenging, and protein mobility on the radiolysis of albumen hydrogel

Nanosecond pulse radiolysis experiments performed on the oxygen and nitrogen saturated thick fraction of egg white (which has an immeasurably high macroviscosity) produced a rate constant for hydrated electron reaction with oxygen equal to 1.7 × 10(10) M(-1) s(-1) at 18 °C, indicating that, due to the very low microviscosity, hydrated electrons are as mobile and reactive in the albumen hydrogel as in neat water. Also, the radiolytic yield for the hydrated electron (G-value) in the thick fraction of egg white (measured at the end of a 14 ns electron pulse) was found to be 86% of that determined in neat water, which can be attributed to the reaction of dry electrons with the protein constituents. Steady-state γ radiolysis studies on air-saturated thick and thin fraction (true solution) revealed that the immobility of egg white proteins disfavors reactions that result in molecular mass change.