Resolving Nonadiabatic Dynamics of Hydrated Electrons Using Ultrafast Photoemission Anisotropy

Theoretical investigation of Aryl/Alkyl halide reduction with hydrated electrons from energy and AIMD aspects

CONTEXT

In this study, the reactions of hydrated electron (e^-(aq)) with alkyl and aryl halides were simulated with an ab initial molecular dynamics (AIMD) method to reveal the underlying mechanism. An original protocol was developed for preparing the proper initial wavefunction guess of AIMD, in which a single electron was curled in a tetrahedral cavity of four water molecules. Our results show that the stability of e^-(aq) increases with the hydrogen bond grid integrity. The organic halides prefer to react with e^-(aq) in neutral or alkaline environment, while they are more likely to react with hydrogen radical (the product of e^-(aq) and proton) under acidic conditions. The reaction between fluorobenzene/fluoromethane and hydrogen radical is considered as the least favorable reaction due to the highest reaction barriers. The bond dissociation energy (BDE) suggested that the cleavage of the carbon-halogen bond of their anion radical might be a thermodynamically favorable reaction. AIMD results indicated that the LUMO or higher orbitals were the e^-(aq) migration destination. The transplanted electron enhanced carbon-halogen bond vibration intensively, leading to bond cleavage. The solvation process of the departing halogen anions was observed in both fluorobenzene and fluoromethane AIMD simulation, indicating that it might have a significant effect on enthalpy. Side reactions and byproducts obtained during the AIMD simulation suggested the complexity of the e^-(aq) reactions and further investigation was needed to fully understand the reaction mechanisms. This study provided theoretical insight into the pollutant environmental fate and constructed a methodological foundation for AIMD simulation of analogous free radical reactions.

METHODS

The theoretical calculation was conducted on the combination of Gaussian16 and ORCA5.0.3 software packages. The initial geometries, as well as the wavefunction initial guesses, were obtained at PBE0/ma-def2-TZVP/IEFPCM-water level in Gaussian16 unless otherwise stated. AIMD simulations were performed at the same level in ORCA. Wavefunction analysis was carried out with Multiwfn. The details methods were described in the section "Computational details" section.

Distortion Correction of Low-Energy Photoelectron Spectra of Liquids Using Spectroscopic Data for Solvated Electrons

Time-resolved photoelectron spectroscopy (TRPES) enables real-time observation of ultrafast electronic dynamics in solutions. When extreme ultraviolet (EUV) probe pulses are employed, they can ionize solutes from all electronic states involved in the dynamics. However, EUV pulses also produce a strong ionization signal from a solvent that is typically 6 orders of magnitude greater than the pump-probe photoelectron signal of solutes. Alternatively, UV probe pulses enable highly sensitive and selective observation of photoexcited solutes because typical solvents such as water are transparent to UV radiation. An obstacle in such UV-TRPES measurements is spectral distortion caused by electron scattering and a yet to be identified mechanism in liquids. We have previously proposed the spectral retrieval (SR) method as an a posteriori approach to removing the distortion and overcoming this difficulty in UV-TRPES; however, its accuracy has not yet been verified by comparison with EUV-TRPES results. In the present study, we perform EUV-TRPES for charge transfer reactions in water, methanol, and ethanol, and verify SR analysis of UV-TRPES. We also estimate a previously undetermined energy-dependent intensity factor and expand the basis sets for SR analysis. The refined SR method is employed for reanalyzing the UV-TRPES data for the formation and relaxation dynamics of solvated electrons in various systems. The electron binding energy distributions for solvated electrons in liquid water, methanol, and ethanol are confirmed to be Gaussian centered at 3.78, 3.39, and 3.25 eV, respectively, in agreement with Nishitani et al. [ Sci. Adv. 2019, 5(8), eaaw6896]. An effective energy gap between the conduction band and the vacuum level at the gas-liquid interface is estimated to be 0.2 eV for liquid water and 0.1 eV for methanol and ethanol.

Photoionization of 3-Methylindole Embedded in Sodium Dodecyl Sulfate and Dodecyltrimethylammonium Chloride Micelles: Migration of Electrons Generated in Micelle Cores and Their Solvation in Outside Water

Electrons were generated in the core of micelles formed by negatively charged sodium dodecyl sulfate (SDS) or positively charged dodecyltrimethylammonium chloride (DTAC) by photoionization of 3-methylindole embedded in the core. The electrons were hydrated after they moved out of the core to the outer aqueous phase. These processes were monitored with femtosecond time-resolved absorption spectroscopy. The migration of electrons from the micelle core to the outer aqueous phase was faster than the instrumental response time of 200 fs. Hot electrons in the aqueous phase were produced in ≤320 fs. There was no significant difference observed for the micellar solutions of negatively charged SDS and positively charged DTAC, or for water. The geminate recombination between the electrons and the radical cations was hindered to a large extent once the electrons hydrated at the outer aqueous phase were separated from the radical cations remaining in the micelle core.

Nonadiabatic Dynamics Studied by Liquid-Jet Time-Resolved Photoelectron Spectroscopy

The development of the liquid microjet technique by Faubel and co-workers has enabled the investigation of high vapor pressure liquids and solutions utilizing high-vacuum methods. One such method is photoelectron spectroscopy (PES), which allows one to probe the electronic properties of a sample through ionization in a state-specific manner. Liquid microjets consisting of pure solvents and solute-solvent systems have been studied with great success utilizing PES and, more recently, time-resolved PES (TRPES). Here, we discuss progress made over recent years in understanding the solvation and excited state dynamics of the solvated electron and nucleic acid constituents (NACs) using these methods, as well as the prospect for their future.The solvated electron is of particular interest in liquid microjet experiments as it represents the simplest solute system. Despite this simplicity, there were still many unresolved questions about its binding energy and excited state relaxation dynamics that are ideal problems for liquid microjet PES. In the work discussed in this Account, accurate binding energies were measured for the solvated electron in multiple high vapor pressure solvents. The advantages of liquid jet PES were further highlighted in the femtosecond excited state relaxation studies on the solvated electron in water where a 75 ± 20 fs lifetime attributable to internal conversion from the excited p-state to a hot ground state was measured, supporting a nonadiabatic relaxation mechanism.Nucleic acid constituents represent a class of important solutes with several unresolved questions that the liquid microjet PES method is uniquely suited to address. As TRPES is capable of tracking dynamics with state-specificity, it is ideal for instances where there are multiple excited states potentially involved in the dynamics. Time-resolved studies of NAC relaxation after excitation using ultraviolet light identified relaxation lifetimes from multiple excited states. The state-specific nature of the TRPES method allowed us to identify the lack of any signal attributable to the ¹nπ* state in thymine derived NACs. The femtosecond time resolution of the technique also aided in identifying differences between the excited state lifetimes of thymidine and thymidine monophosphate. These have been interpreted, aided by molecular dynamics simulations, as an influence of conformational differences leading to a longer excited state lifetime in thymidine monophosphate.Finally, we discuss advances in tabletop light sources extending into the extreme ultraviolet and soft X-ray regimes that allow expansion of liquid jet TRPES to full valence band and potentially core level studies of solutes and pure liquids in liquid microjets. As most solutes have ground state binding energies in the range of 10 eV, observation of both excited state decay and ground state recovery using ultraviolet pump-ultraviolet probe TRPES has been intractable. With high-harmonic generation light sources, it will be possible to not only observe complete relaxation pathways for valence level dynamics but to also track dynamics with element specificity by probing core levels of the solute of interest.

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.

Size-Resolved Electron Solvation in Neutral Water Clusters

Cluster-size-resolved ultrafast dynamics of the solvated electron in neutral water clusters with n = 3 to ∼200 molecules are studied with pump–probe time-of-flight mass spectrometry after below band gap excitation. For the smallest clusters, no longer-lived (>100–200 fs) hydrated electrons were detected, indicating a minimum size of n ∼ 14 for being able to sustain hydrated electrons. Larger clusters show a systematic increase of the number of hydrated electrons per molecule on the femtosecond to picosecond time scale. We propose that with increasing cluster size the underlying dynamics is governed by more effective electron formation processes combined with less effective electron loss processes, such as ultrafast hydrogen ejection and recombination. It appears unlikely that any size dependence of the solvent relaxation dynamics would be reflected in the observed time-resolved ion yields.

Design and characterization of a magnetic bottle electron spectrometer for time-resolved extreme UV and X-ray photoemission spectroscopy of liquid microjets

We describe a magnetic bottle time-of-flight electron spectrometer designed for time-resolved photoemission spectroscopy of a liquid microjet using extreme UV and X-ray radiation. The spectrometer can be easily reconfigured depending on experimental requirements and the energy range of interest. To improve the energy resolution at high electron kinetic energy, a retarding potential can be applied either via a stack of electrodes or retarding mesh grids, and a flight-tube extension can be attached to increase the flight time. A gated electron detector was developed to reject intense parasitic signal from light scattered off the surface of the cylindrically shaped liquid microjet. This detector features a two-stage multiplication with a microchannel plate plus a fast-response scintillator followed by an image-intensified photon detector. The performance of the spectrometer was tested at SPring-8 and SACLA, and time-resolved photoelectron spectra were measured for an ultrafast charge transfer to solvent reaction in an aqueous NaI solution with a 200 nm UV pump pulses from a table-top ultrafast laser and the 5.5 keV hard X-ray probe pulses from SACLA.

Below Band Gap Formation of Solvated Electrons in Neutral Water Clusters?

Below band gap formation of solvated electrons in neutral water clusters using pump-probe photoelectron imaging is compared with recent data for liquid water and with above band gap excitation studies in liquid and clusters. Similar relaxation times on the order of 200 fs and 1-2 ps are retrieved for below and above band gap excitation, in both clusters and liquid. The independence of the relaxation times from the generation process indicates that these times are dominated by the solvent response, which is significantly slower than the various solvated electron formation processes. The analysis of the temporal evolution of the vertical electron binding energy and the electron binding energy at half-maximum suggests a dependence of the solvation time on the binding energy.

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.

Real-time observation of water radiolysis and hydrated electron formation induced by extreme-ultraviolet pulses

The dominant pathway of radiation damage begins with the ionization of water. Thus far, however, the underlying primary processes could not be conclusively elucidated. Here, we directly study the earliest steps of extreme ultraviolet (XUV)-induced water radiolysis through one-photon excitation of large water clusters using time-resolved photoelectron imaging. Results are presented for H2O and D2O clusters using femtosecond pump pulses centered at 133 or 80 nm. In both excitation schemes, hydrogen or proton transfer is observed to yield a prehydrated electron within 30 to 60 fs, followed by its solvation in 0.3 to 1.0 ps and its decay through geminate recombination on a ∼10-ps time scale. These results are interpreted by comparison with detailed multiconfigurational non-adiabatic ab-initio molecular dynamics calculations. Our results provide the first comprehensive picture of the primary steps of radiation chemistry and radiation damage and demonstrate new approaches for their study with unprecedented time resolution.

Probing the Structural Evolution of the Hydrated Electron in Water Cluster Anions (H2O)n-, n ≤ 200, by Electronic Absorption Spectroscopy

Electronic absorption spectra of water cluster anions (H₂O)_n^–, n ≤ 200, at T = 80 K are obtained by photodissociation spectroscopy and compared with simulations from literature and experimental data for bulk hydrated electrons. Two almost isoenergetic electron binding motifs are seen for cluster sizes 20 ≤ n ≤ 40, which are assigned to surface and partially embedded isomers. With increasing cluster size, the surface isomer becomes less populated, and for n ≥ 50, the partially embedded isomer prevails. The absorption shifts to the blue, reaching a plateau at n ≈ 100. In this size range, the absorption spectrum is similar to that of the bulk hydrated electron but is slightly red-shifted; spectral moment analysis indicates that these clusters are reasonable model systems for hydrated electrons near the liquid–vacuum interface.

Relaxation Dynamics and Genuine Properties of the Solvated Electron in Neutral Water Clusters

We have investigated the solvation dynamics and the genuine binding energy and photoemission anisotropy of the solvated electron in neutral water clusters with a combination of time-resolved photoelectron velocity map imaging and electron scattering simulations. The dynamics was probed with a UV probe pulse following above-band-gap excitation by an EUV pump pulse. The solvation dynamics is completed within about 2 ps. Only a single band is observed in the spectra, with no indication for isomers with distinct binding energies. Data analysis with an electron scattering model reveals a genuine binding energy in the range of 3.55-3.85 eV and a genuine anisotropy parameter in the range of 0.51-0.66 for the ground-state hydrated electron. All of these observations coincide with those for liquid bulk, which is rather unexpected for an average cluster size of 300 molecules.

Ultrafast Internal Conversion and Solvation of Electrons in Water, Methanol, and Ethanol

Ultrafast internal conversion from the first excited state of a solvated electron in water, methanol, and ethanol is investigated using time-resolved photoelectron spectroscopy of liquid microjets and a spectral retrieval method. Photoelectron spectra corrected for inelastic scattering clearly reveal well-separated signals from the excited and ground states, and the latter enables us to analyze the solvation dynamics in the ground state after internal conversion. Measurements with 25 fs time resolution identify a rapid increase in the vertical electron binding energy of the solvated electron owing to nuclear wave packet motions in the excited state and allow us to precisely determine the internal conversion time.

Binding energy of solvated electrons and retrieval of true UV photoelectron spectra of liquids

The electronic energy and dynamics of solvated electrons, the simplest yet elusive chemical species, is of interest in chemistry, physics, and biology. Here, we present the electron binding energy distributions of solvated electrons in liquid water, methanol, and ethanol accurately measured using extreme ultraviolet (EUV) photoelectron spectroscopy of liquids with a single-order high harmonic. The distributions are Gaussian in all cases. Using the EUV and UV photoelectron spectra of solvated electrons, we succeeded in retrieving sharp electron kinetic energy distributions from the spectra broadened and energy shifted by inelastic scattering in liquids, overcoming an obstacle in ultrafast UV photoelectron spectroscopy of liquids. The method is demonstrated for the benchmark systems of charge transfer to solvent reaction and ultrafast internal conversion of hydrated electron from the first excited state.

Design and characterization of a recirculating liquid-microjet photoelectron spectrometer for multiphoton ultraviolet photoelectron spectroscopy

A new recirculating liquid-microjet photoelectron spectrometer for multiphoton ultraviolet photoelectron spectroscopy is described. A recirculating system is essential for studying samples that are only available in relatively small quantities. The reduction in background pressure when using the recirculating system compared to a liquid-nitrogen cold-trap results in a significant improvement in the quality of the photoelectron spectra. Moreover, the recirculating system results in a negligible streaming potential. The instrument design, operation, and characterization are described in detail, and its performance is illustrated by comparing a photoelectron spectrum of aqueous phenol recorded using the recirculating system with one recorded using a liquid nitrogen cold-trap.

Magic Numbers for the Photoelectron Anisotropy in Li-Doped Dimethyl Ether Clusters

Photoelectron velocity map imaging of Li(CH₃OCH₃)_n clusters (1 ≤ n ≤ 175) is used to search for magic numbers related to the photoelectron anisotropy. Comparison with density functional calculations reveals magic numbers at n = 4, 5, and 6, resulting from the symmetric charge distribution with high s-character of the highest occupied molecular orbital. Since each of these three cluster sizes correspond to the completion of a first coordination shell, they can be considered as “isomeric motifs of the first coordination shell”. Differences in the photoelectron anisotropy, the vertical ionization energies and the enthalpies of vaporization between Li(CH₃OCH₃)_n and Na(CH₃OCH₃)_n can be rationalized in terms of differences in their solvation shells, atomic ionization energies, polarizabilities, metal–oxygen bonds, ligand–ligand interactions and by cooperative effects.

Observation of dissociative quasi-free electron attachment to nucleoside via excited anion radical in solution

Damage to DNA via dissociative electron attachment has been well-studied in both the gas and condensed phases; however, understanding this process in bulk solution at a fundamental level is still a challenge. Here, we use a picosecond pulse of a high energy electron beam to generate electrons in liquid diethylene glycol and observe the electron attachment dynamics to ribothymidine at different stages of electron relaxation. Our transient spectroscopic results reveal that the quasi-free electron with energy near the conduction band effectively attaches to ribothymidine leading to a new absorbing species that is characterized in the UV-visible region. This species exhibits a nearly concentration-independent decay with a time constant of ~350 ps. From time-resolved studies under different conditions, combined with data analysis and theoretical calculations, we assign this intermediate to an excited anion radical that undergoes N1-C1′ glycosidic bond dissociation rather than relaxation to its ground state.

Radiation-induced low-energy electrons in solution are implicated in DNA damage, but their relaxation dynamics are not well understood. Here the authors observe how quasi-free electrons dissociate glycosidic bonds via an excited nucleoside anion radical, whereas solvated electrons reside on the nucleoside as a relatively stable anion radical.

Structure of the aqueous electron

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

Tuning Solvated Electrons by Polar-Nonpolar Oxide Heterostructure

Solvated electron states at the oxide/aqueous interface represent the lowest energy charge-transfer pathways, thereby playing an important role in photocatalysis and electronic device applications. However, their energies are usually higher than the conduction band minimum (CBM), which makes the solvated electrons difficult to utilize in charge-transfer processes. Thus it is essential to stabilize the energy of the solvated electron states. Taking LaAlO₃/SrTiO₃ (LAO/STO) oxide heterostructure with H₂O-adsorbed monolayer as a prototypical system, we show using DFT and ab initio time-dependent nonadiabatic molecular dynamics simulation that the energy and dynamics of solvated electrons can be tuned by the electric field in the polar-nonpolar oxide heterostructure. In particular, for LAO/STO with p-type interface, the CBM is contributed by the solvated electron state when LAO is thicker than four unit cells. Furthermore, the solvated electron band minimum can be partially occupied when LAO is thicker than eight unit cells. We propose that the tunability of solvated electron states can be achieved on polar-nonpolar oxide heterostructure surfaces as well as on ferroelectric oxides, which is important for charge and proton transfer at oxide/aqueous interfaces.

Angle-resolved photoemission spectroscopy of liquid water at 29.5 eV

Angle-resolved photoemission spectroscopy of liquid water was performed using extreme ultraviolet radiation at 29.5 eV and a time-of-flight photoelectron spectrometer. SiC/Mg coated mirrors were employed to select the single-order 19th harmonic from laser high harmonics, which provided a constant photon flux for different laser polarizations. The instrument was tested by measuring photoemission anisotropy for rare gases and water molecules and applied to a microjet of an aqueous NaI solution. The solute concentration was adjusted to eliminate an electric field gradient around the microjet. The observed photoelectron spectra were analyzed considering contributions from liquid water, water vapor, and an isotropic background. The anisotropy parameters of the valence bands (1b1, 3a1, and 1b2) of liquid water are considerably smaller than those of gaseous water, which is primarily attributed to electron scattering in liquid water.

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.

Genuine binding energy of the hydrated electron

The unknown influence of inelastic and elastic scattering of slow electrons in water has made it difficult to clarify the role of the solvated electron in radiation chemistry and biology. We combine accurate scattering simulations with experimental photoemission spectroscopy of the hydrated electron in a liquid water microjet, with the aim of resolving ambiguities regarding the influence of electron scattering on binding energy spectra, photoelectron angular distributions, and probing depths. The scattering parameters used in the simulations are retrieved from independent photoemission experiments of water droplets. For the ground-state hydrated electron, we report genuine values devoid of scattering contributions for the vertical binding energy and the anisotropy parameter of 3.7 ± 0.1 eV and 0.6 ± 0.2, respectively. Our probing depths suggest that even vacuum ultraviolet probing is not particularly surface-selective. Our work demonstrates the importance of quantitative scattering simulations for a detailed analysis of key properties of the hydrated electron.

Size-Resolved Photoelectron Anisotropy of Gas Phase Water Clusters and Predictions for Liquid Water

We report the first measurements of size-resolved photoelectron angular distributions for the valence orbitals of neutral water clusters with up to 20 molecules. A systematic decrease of the photoelectron anisotropy is found for clusters with up to 5-6 molecules, and most remarkably, convergence of the anisotropy for larger clusters. We suggest the latter to be the result of a local short-range scattering potential that is fully described by a unit of 5-6 molecules. The cluster data and a detailed electron scattering model are used to predict the anisotropy of slow photoelectrons in liquid water. Reasonable agreement with experimental liquid jet data is found.