Dynamics of Excited Solvated Electrons in Aqueous Solution Monitored with Femtosecond-Time and Polarization Resolution

Simulating the Competitive Ion Pairing of Hydrated Electrons with Chaotropic Cations

Experiments show that the absorption spectrum of the hydrated electron (ehyd-) blue-shifts in electrolyte solutions compared with what is seen in pure water. This shift has been assigned to the ehyd-'s competitive ion-pairing interactions with the salt cation relative to the salt anion based on the ions' positions on the Hofmeister series. Remarkably, little work has been done investigating the ehyd-'s behavior when the salts have chaotropic cations, which should greatly change the ion-pairing interactions given that the ehyd- is a champion chaotrope. In this work, we remedy this by using mixed quantum/classical simulations to analyze the behavior of two different models of the ehyd- in aqueous RbF and RbI electrolyte solutions as a function of salt concentration. We find that the magnitude of the salt-induced spectral blue-shift is determined by a combination of the number of chaotropic Rb+ cations near the ehyd- and the number of salt anions near those cations so that the spectrum of the ehyd- directly reflects its local environment. We also find that the use of a soft-cavity ehyd- model predicts stronger competitive interactions with Rb+ relative to I- than a more traditional hard cavity model, leading to different predicted spectral shifts that should provide a way to distinguish between the two models experimentally. Our simulations predict that at the same concentration, salts with chaotropic cations should produce larger spectral blue-shifts than salts with kosmotropic cations. We also found that at high salt concentrations with chaotropic cations, the predicted blue-shift is greater when the salt anion is kosmotropic instead of chaotropic. Our goal is for this work to inspire experimentalists to make such measurements, which will help provide a spectroscopic means to distinguish between simulations models that predict different hydration structures for the ehyd-.

First-principles simulation of an ejected electron produced by monochromatic deposition energy to water at the femtosecond order

This study uses a time-dependent first-principles simulation code to investigate the transient dynamics of an ejected electron produced in the monochromatic deposition energy from 11 to 19 eV in water. The energy deposition forms a three-body single spur comprising a hydroxyl radical (OH˙), hydronium ion (H3O+), and hydrated electron (eaq-). The earliest formation involves electron thermalization and delocalization dominated by the molecular excitation of water. Our simulation results show that the transient electron dynamics primarily depends on the amount of deposition energy to water; the thermalization time varies from 200 to 500 fs, and the delocalization varies from 3 to 10 nm in this energy range. These features are crucial for determining the earliest single-spur formation and facilitating a sequential simulation from an energy deposition to a chemical reaction in water photolysis or radiolysis. The spur radius obtained from the simulation correlates reasonably with the experimental-based estimations. Our results should provide universalistic insights for analysing ultrafast phenomena dominated by the molecular excitation of water in the femtosecond order.

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

Initial yield of hydrated electron production from water radiolysis based on first-principles calculation

Many scientific insights into water radiolysis have been applied for developing life science, including radiation-induced phenomena, such as DNA damage and mutation induction or carcinogenesis. However, the generation mechanism of free radicals due to radiolysis remains to be fully understood. Consequently, we have encountered a crucial problem in that the initial yields connecting radiation physics to chemistry must be parameterized. We have been challenged in the development of a simulation tool that can unravel the initial free radical yields, from physical interaction by radiation. The presented code enables the first-principles calculation of low energy secondary electrons resulting from the ionization, in which the secondary electron dynamics are simulated while considering dominant collision and polarization effects in water. In this study, using this code, we predicted the yield ratio between ionization and electronic excitation from a delocalization distribution of secondary electrons. The simulation result presented a theoretical initial yield of hydrated electrons. In radiation physics, the initial yield predicted from parameter analysis of radiolysis experiments in radiation chemistry was successfully reproduced. Our simulation code helps realize a reasonable spatiotemporal connection from radiation physics to chemistry, which would contribute to providing new scientific insights for precise understanding of underlying mechanisms of DNA damage induction.

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.

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.

Ultrafast Processes Occurring in Radiolysis of Highly Concentrated Solutions of Nucleosides/Tides

Among the radicals (hydroxyl radical (^•OH), hydrogen atom (H^•), and solvated electron (e_sol⁻)) that are generated via water radiolysis, ^•OH has been shown to be the main transient species responsible for radiation damage to DNA via the indirect effect. Reactions of these radicals with DNA-model systems (bases, nucleosides, nucleotides, polynucleotides of defined sequences, single stranded (ss) and double stranded (ds) highly polymeric DNA, nucleohistones) were extensively investigated. The timescale of the reactions of these radicals with DNA-models range from nanoseconds (ns) to microseconds (µs) at ambient temperature and are controlled by diffusion or activation. However, those studies carried out in dilute solutions that model radiation damage to DNA via indirect action do not turn out to be valid in dense biological medium, where solute and water molecules are in close contact (e.g., in cellular environment). In that case, the initial species formed from water radiolysis are two radicals that are ultrashort-lived and charged: the water cation radical (H₂O^•+) and prethermalized electron. These species are captured by target biomolecules (e.g., DNA, proteins, etc.) in competition with their inherent pathways of proton transfer and relaxation occurring in less than 1 picosecond. In addition, the direct-type effects of radiation, i.e., ionization of macromolecule plus excitations proximate to ionizations, become important. The holes (i.e., unpaired spin or cation radical sites) created by ionization undergo fast spin transfer across DNA subunits. The exploration of the above-mentioned ultrafast processes is crucial to elucidate our understanding of the mechanisms that are involved in causing DNA damage via direct-type effects of radiation. Only recently, investigations of these ultrafast processes have been attempted by studying concentrated solutions of nucleosides/tides under ambient conditions. Recent advancements of laser-driven picosecond electron accelerators have provided an opportunity to address some long-term puzzling questions in the context of direct-type and indirect effects of DNA damage. In this review, we have presented key findings that are important to elucidate mechanisms of complex processes including excess electron-mediated bond breakage and hole transfer, occurring at the single nucleoside/tide level.

Electronic Excited State Lifetimes of Anionic Water Clusters: Dependence on Charge Solvation Motif

An ongoing controversy about water cluster anions concerns the electron-binding motif, whether the charge center is localized at the surface or within the cluster interior. Here, mixed quantum-classical dynamics simulations have been carried out for a wide range of cluster sizes (n ≤ 1000) for (H₂O)_n^- and (D₂O)_n^-, based on a nonequilibrium first-order rate constant approach. The computed data are in good general agreement with time-resolved photoelectron imaging results (n ≤ 200). The analysis reveals that, for surface state electrons, the cluster size dependence of the excited state electronic energy gap and the magnitude of the nonadiabatic couplings have compensating influences on the excited state lifetimes: the excited state lifetime for surface states reaches a minimum for n ∼ 150 and then increases for larger clusters. It is concluded that the electron resides in a surface-localized motif in all of these measured clusters, dominating at least up to n = 200.

Resolving Nonadiabatic Dynamics of Hydrated Electrons Using Ultrafast Photoemission Anisotropy

We have studied ultrafast nonadiabatic dynamics of excess electrons trapped in the band gap of liquid water using time- and angle-resolved photoemission spectroscopy. Anisotropic photoemission from the first excited state was discovered, which enabled unambiguous identification of nonadiabatic transition to the ground state in 60 fs in H_{2}O and 100 fs in D_{2}O. The photoelectron kinetic energy distribution exhibited a rapid spectral shift in ca. 20 fs, which is ascribed to the librational response of a hydration shell to electronic excitation. Photoemission anisotropy indicates that the electron orbital in the excited state is depolarized in less than 40 fs.

To be or not to be in a cavity: the hydrated electron dilemma

The hydrated electron-the species that results from the addition of a single excess electron to liquid water-has been the focus of much interest both because of its role in radiation chemistry and other chemical reactions, and because it provides for a deceptively simple system that can serve as a means to confront the predictions of quantum molecular dynamics simulations with experiment. Despite all this interest, there is still considerable debate over the molecular structure of the hydrated electron: does it occupy a cavity, have a significant number of interior water molecules, or have a structure somewhere in between? The reason for all this debate is that different computer simulations have produced each of these different structures, yet the predicted properties for these different structures are still in reasonable agreement with experiment. In this Feature Article, we explore the reasons underlying why different structures are produced when different pseudopotentials are used in quantum simulations of the hydrated electron. We also show that essentially all the different models for the hydrated electron, including those from fully ab initio calculations, have relatively little direct overlap of the electron's wave function with the nearby water molecules. Thus, a non-cavity hydrated electron is better thought of as an "inverse plum pudding" model, with interior waters that locally expel the surrounding electron's charge density. Finally, we also explore the agreement between different hydrated electron models and certain key experiments, such as resonance Raman spectroscopy and the temperature dependence and degree of homogeneous broadening of the optical absorption spectrum, in order to distinguish between the different simulated structures. Taken together, we conclude that the hydrated electron likely has a significant number of interior water molecules.

Femtosecond relaxation dynamics of solvated electrons in liquid ammonia

The ultrafast relaxation dynamics of the well-known solvated electron in liquid ammonia solutions are investigated with femtosecond near-infrared pump-probe absorption spectroscopy. Immediately after photoexcitation, the dynamic absorption spectrum of the electron is substantially red-shifted with respect to its stationary spectrum. A subsequent dynamic blue shift of the pump-probe spectrum occurs on a timescale of 150 fs. The data are understood in terms of ground-state "cooling" and can be quantitatively simulated by an intuitive temperature-jump model employing a dynamically evolving Kubo line shape for the electronic resonance. A simple estimate implies that, on average, the electron in the liquid is coordinated to six nearest-neighbor ammonia molecules. An equivalent analysis of the data based on a bubble-formation/cavity-contraction mechanism is briefly outlined.

Temperature effect on the absorption spectrum of the hydrated electron paired with a lithium cation in deuterated water

The absorption spectra of the hydrated electron in 1.0 to 4.0 M LiCl or LiClO4 deuterated water solutions were measured by pulse radiolysis techniques from room temperature to 300 degrees C at a constant pressure of 25 MPa. The results show that when the temperature is increased and the density is decreased, the absorption spectrum of the electron in the presence of a lithium cation is shifted to lower energies. Quantum classical molecular dynamics (QCMD) simulations of an excess electron in bulk water and in the presence of a lithium cation have been performed to compare with the experimental results. According to the QCMD simulations, the change in the shape of the spectrum is due to one of the three p-like excited states of the solvated electron destabilized by core repulsion. The study of s --> p transition energies for the three p-excited states reveals that for temperatures higher than room temperature, there is a broadening of each individual s --> p absorption band due to a less structured water solvation shell.

Probing chemical dynamics with negative ions

Experiments are reviewed in which key problems in chemical dynamics are probed by experiments based on photodetachment and/or photoexcitation of negative ions. Examples include transition state spectroscopy of biomolecular reactions, spectroscopy of open shell van der Waals complexes, photodissociation of free radicals, and time-resolved dynamics in clusters. The experimental methods used in these investigations are described along with representative systems that have been studied.

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

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

Photodetachment of Ferrocyanide in Reverse Micelles

Solvated electrons have been generated in reverse micelles (RMs) through photodetachment of ferrocyanide (Fe(CN)(6)(4-)) in sodium bis(2-ethylhexyl) sulfosuccinate (AOT) RMs. We have measured both bleach recovery of the parent ferrocyanide CN stretch in the infrared and the decay of the solvated electron absorption at 800 nm. The bleach recovery has been fit to a diffusion model for the geminate recombination process. The fit parameters suggest a narrowing of the spatial distribution of ejected electrons due to confinement in the RMs when compared to bulk water. The diffusion coefficient of the solvated electron does not appear to be significantly affected by RM confinement. The decay of the solvated electron absorption exhibits an additional decay component that is not observed in bulk water and is smaller for larger RMs. No corresponding additional component is seen in the parent ferrocyanide IR bleach recovery, which supports our interpretation that the confinement-induced new decay process in RMs is due to electrons reacting with AOT headgroups.

Spectroscopic Study of Laser-Induced Phase Transition of Gold Nanoparticles on Nanosecond Time Scales and Longer

The pulsed laser induced phase transition of gold nanoparticles in aqueous solution was observed via a transient absorption on nanosecond time scales and longer. Gold nanoparticles were excited with an intense picosecond laser pulse (355 nm, 30 ps), and the subsequent changes were monitored using two continuous wave laser wavelengths (488 and 635 nm). On the nanosecond time scale, below 6.3 mJ cm(-2), no change was observed; however, in the low fluence region between 6.3 and 17 mJ cm(-2), gold nanoparticles produced a bleach signal (488 nm) attributed to the melting of the gold nanoparticles, which decreased linearly with increasing laser fluence. Laser fluences above 17 mJ cm(-2) resulted in a strong absorption at both wavelengths, which is ascribed to vaporization of gold nanoparticles rather than solvated electrons (ejected from gold nanoparticles) or light scattering. The decay of both signals was faster than the 5 ns time resolution used in our experimental system. On the microsecond time scale, increase in absorbance at 635 nm was observed with a time constant of 1.0 micros, while no change was observed at 488 nm. It is considered that this increase is attributed to the formation of smaller gold nanoparticles resulting from pulsed laser induced size reduction of initial gold nanoparticles.

Electronic Relaxation Dynamics of Water Cluster Anions

The electronic relaxation dynamics of water cluster anions, (H(2)O)(n)(-), have been studied with time-resolved photoelectron imaging. In this investigation, the excess electron was excited through the p<--s transition with an ultrafast laser pulse, with subsequent electronic evolution monitored by photodetachment. All excited-state lifetimes exhibit a significant isotope effect (tau(D)2(O)/tau(H)2(O) approximately 2). Additionally, marked dynamical differences are found for two classes of water cluster anions, isomers I and II, previously assigned as clusters with internally solvated and surface-bound electrons, respectively. Isomer I clusters with n > or = 25 decay exclusively by internal conversion, with relaxation times that extrapolate linearly with 1/n toward an internal conversion lifetime of 50 fs in bulk water. Smaller isomer I clusters (13 < or = n < or = 25) decay through a combination of excited-state autodetachment and internal conversion. The relaxation of isomer II clusters shows no significant size dependence over the range of n = 60-100, with autodetachment an important decay channel following excitation of these clusters. Photoelectron angular distributions (PADs) were measured for isomer I and isomer II clusters. The large differences in dynamical trends, relaxation mechanisms, and PADs between large isomer I and isomer II clusters are consistent with their assignment to very different electron binding motifs.

Solvation and Thermalization of Electrons Generated by above-the-Gap (12.4 eV) Two-Photon Ionization of Liquid H2O and D2O

Temporal evolution of transient absorption (TA) spectra of electrons generated by above-the-gap (12.4 eV total energy) two-photon ionization of liquid H2O and D2O has been studied on femto- and picosecond time scales. The spectra were obtained at intervals of 50 nm between 0.5 and 1.7 mum. Two distinct regimes of the spectral evolution were observed: t < 1 ps and t > 1 ps. In both of these regimes, the spectral profile changes considerably with the delay time of the probe pulse. The "continuous blue shift" and the "temperature jump" models, in which the spectral profile does not change as it progressively shifts, as a whole, to the blue, are not supported by our data. Furthermore, no p-state electron, postulated by several authors to be a short-lived intermediate of the photoionization process, was observed by the end of the 300 fs, 200 nm pump pulse. For t < 1 ps, two new TA features (the 1.15 microm peak and 1.4 mum shoulder) were observed for the electron in the spectral region where O-H overtones appear in the spectra of light water. These two features were not observed for the electron in D2O. The 1.4 mum peak observed in D2O may be the isotope-shift analogue of the 1.15 microm feature in H2O. Vibronic coupling to the modes of water molecules lining the solvation cavity is a possible origin of these features. On the sub-picosecond time scale, the absorption band of solvated electron progressively shifts to the blue. At later delay times (t > 1 ps), the position of the band maximum is "locked", but the spectral profile continues to change by narrowing on the red side and broadening on the blue side; the oscillator strength is constant within 10%. The time constant of this narrowing is ca. 0.56 ps for H2O and 0.64 ps for D2O. Vibrational relaxation and time-dependent decrease in the size and sphericity of the solvation cavity are suggested as possible causes for the observed spectral transformations in both of these regimes.