Hydrated electrons as a probe of local anisotropy: Simulations of ultrafast polarization-dependent spectral hole burning

Simulating ultrafast transient absorption spectra from first principles using a time-dependent configuration interaction probe

Transient absorption spectroscopy (TAS) is among the most common ultrafast photochemical experiments, but its interpretation remains challenging. In this work, we present an efficient and robust method for simulating TAS signals from first principles. Excited-state absorption and stimulated emission (SE) signals are computed using time-dependent complete active space configuration interaction (TD-CASCI) simulations, leveraging the robustness of time-domain simulation to minimize electronic structure failure. We demonstrate our approach by simulating the TAS signal of 1'-hydroxy-2'-acetonapthone (HAN) from ab initio multiple spawning nonadiabatic molecular dynamics simulations. Our results are compared to gas-phase TAS data recorded from both jet-cooled (T ∼ 40 K) and hot (∼403 K) molecules via cavity-enhanced TAS (CE-TAS). Decomposition of the computed spectrum allows us to assign a rise in the SE signal to excited-state proton transfer and the ultimate decay of the signal to relaxation through a twisted conical intersection. The total cost of computing the observable signal (∼1700 graphics processing unit hours for ∼4 ns of electron dynamics) was markedly less than that of performing the ab initio multiple spawning calculations used to compute the underlying nonadiabatic dynamics.

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

Failure of molecular dynamics to provide appropriate structures for quantum mechanical description of the aqueous chloride ion charge-transfer-to-solvent ultraviolet spectrum

The lowest band in the charge-transfer-to-solvent ultraviolet absorption spectrum of aqueous chloride ion is studied by experiment and computation. Interestingly, the experiments indicate that at concentrations up to at least 0.25 M, where calculations indicate ion pairing to be significant, there is no notable effect of ionic strength on the spectrum. The experimental spectra are fitted to aid comparison with computations. Classical molecular dynamic simulations are carried out on dilute aqueous Cl-, Na+, and NaCl, producing radial distribution functions in reasonable agreement with experiment and, for NaCl, clear evidence of ion pairing. Clusters are extracted from the simulations for quantum mechanical excited state calculations. Accurate ab initio coupled-cluster benchmark calculations on a small number of representative clusters are carried out and used to identify and validate an efficient protocol based on time-dependent density functional theory. The latter is used to carry out quantum mechanical calculations on thousands of clusters. The resulting computed spectrum is in excellent agreement with experiment for the peak position, with little influence from ion pairing, but is in qualitative disagreement on the width, being only about half as wide. It is concluded that simulation by classical molecular dynamics fails to provide an adequate variety of structures to explain the experimental CTTS spectrum of aqueous Cl-.

Structure of the aqueous electron

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

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.

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.

The ultrafast charge-transfer-to-solvent dynamics of iodide in tetrahydrofuran. 1. Exploring the roles of solvent and solute electronic structure in condensed-phase charge-transfer reactions

Although they represent the simplest possible charge-transfer reactions, the charge-transfer-to-solvent (CTTS) dynamics of atomic anions exhibit considerable complexity. For example, the CTTS dynamics of iodide in water are very different from those of sodide (Na-) in tetrahydrofuran (THF), leading to the question of the relative importance of the solvent and solute electronic structures in controlling charge-transfer dynamics. In this work, we address this issue by investigating the CTTS spectroscopy and dynamics of I- in THF, allowing us to make detailed comparisons to the previously studied I-/H2O and Na-/THF CTTS systems. Since THF is weakly polar, ion pairing with the counterion can have a substantial impact on the CTTS spectroscopy and dynamics of I- in this solvent. In this study, we have isolated "counterion-free" I- in THF by complexing the Na+ counterion with 18-crown-6 ether. Ultrafast pump-probe experiments reveal that THF-solvated electrons (e-THF) appear 380 +/- 60 fs following the CTTS excitation of "free" I- in THF. The absorption kinetics are identical at all probe wavelengths, indicating that the ejected electrons appear with no significant dynamic solvation but rather with their equilibrium absorption spectrum. After their initial appearance, ejected electrons do not exhibit any additional dynamics on time scales up to approximately 1 ns, indicating that geminate recombination of e-THF with its iodine atom partner does not occur. Competitive electron scavenging measurements demonstrate that the CTTS excited state of I- in THF is quite large and has contact with scavengers that are several nanometers away from the iodide ion. The ejection time and lack of electron solvation observed for I- in THF are similar to what is observed following CTTS excitation of Na- in THF. However, the relatively slow ejection time, the complete lack of dynamic solvation, and the large ejection distance/lack of recombination dynamics are in marked contrast to the CTTS dynamics observed for I- in water, in which fast electron ejection, substantial solvation, and appreciable recombination have been observed. These differences in dynamical behavior can be understood in terms of the presence of preexisting, electropositive cavities in liquid THF that are a natural part of its liquid structure; these cavities provide a mechanism for excited electrons to relocate to places in the liquid that can be nanometers away, explaining the large ejection distance and lack of recombination following the CTTS excitation of I- in THF. We argue that the lack of dynamic solvation observed following CTTS excitation of both I- and Na- in THF is a direct consequence of the fact that little additional relaxation is required once an excited electron nonadiabatically relaxes into one of the preexisting cavities. In contrast, liquid water contains no such cavities, and CTTS excitation of I- in water leads to local electron ejection that involves substantial solvent reorganization.

Watching Na Atoms Solvate into (Na+,e-) Contact Pairs: Untangling the Ultrafast Charge-Transfer-to-Solvent Dynamics of Na- in Tetrahydrofuran (THF)

With the large dye molecules employed in typical studies of solvation dynamics, it is often difficult to separate the intramolecular relaxation of the dye from the relaxation associated with dynamic solvation. One way to avoid this difficulty is to study solvation dynamics using an atom as the solvation probe; because atoms have only electronic degrees of freedom, all of the observed spectroscopic dynamics must result from motions of the solvent. In this paper, we use ultrafast transient absorption spectroscopy to investigate the solvation dynamics of newly created sodium atoms that are formed following the charge transfer to solvent (CTTS) ejection of an electron from sodium anions (sodide) in liquid tetrahydrofuran (THF). Because the absorption spectra of the sodide reactant, the sodium atom, and the solvated electron products overlap, we first examined the dynamics of the ejected CTTS electron in the infrared to build a detailed model of the CTTS process that allowed us to subtract the spectroscopic contributions of the sodide bleach and the solvated electron and cleanly reveal the spectroscopy of the solvated atom. We find that the neutral sodium species created following CTTS excitation of sodide initially absorbs near 590 nm, the position of the gas-phase sodium D-line, suggesting that it only weakly interacts with the surrounding solvent. We then see a fast solvation process that causes a red-shift of the sodium atom's spectrum in approximately 230 fs, a time scale that matches well with the results of MD simulations of solvation dynamics in liquid THF. After the fast solvation is complete, the neutral sodium atoms undergo a chemical reaction that takes place in approximately 740 fs, as indicated by the observation of an isosbestic point and the creation of a species with a new spectrum. The spectrum of the species created after the reaction then red-shifts on a approximately 10-ps time scale to become the equilibrium spectrum of the THF-solvated sodium atom, which is known from radiation chemistry experiments to absorb near approximately 900 nm. There has been considerable debate as to whether this 900-nm absorbing species is better thought of as a solvated atom or a sodium cation:solvated electron contact pair, (Na+,e-). The fact that we observe the initially created neutral Na atom undergoing a chemical reaction to ultimately become the 900-nm absorbing species suggests that it is better assigned as (Na+,e-). The approximately 10-ps solvation time we observe for this species is an order of magnitude slower than any other solvation process previously observed in liquid THF, suggesting that this species interacts differently with the solvent than the large molecules that are typically used as solvation probes. Together, all of the results allow us to build the most detailed picture to date of the CTTS process of Na- in THF as well as to directly observe the solvation dynamics associated with single sodium atoms in solution.

Moving solvated electrons with light: Nonadiabatic mixed quantum/classical molecular dynamics simulations of the relocalization of photoexcited solvated electrons in tetrahydrofuran (THF)

Motivated by recent ultrafast spectroscopic experiments [Martini et al., Science 293, 462 (2001)], which suggest that photoexcited solvated electrons in tetrahydrofuran (THF) can relocalize (that is, return to equilibrium in solvent cavities far from where they started), we performed a series of nonequilibrium, nonadiabatic, mixed quantum/classical molecular dynamics simulations that mimic one-photon excitation of the THF-solvated electron. We find that as photoexcited THF-solvated electrons relax to their ground states either by continuous mixing from the excited state or via nonadiabatic transitions, approximately 30% of them relocalize into cavities that can be over 1 nm away from where they originated, in close agreement with the experiments. A detailed investigation shows that the ability of excited THF-solvated electrons to undergo photoinduced relocalization stems from the existence of preexisting cavity traps that are an intrinsic part of the structure of liquid THF. This explains why solvated electrons can undergo photoinduced relocalization in solvents like THF but not in solvents like water, which lack the preexisting traps necessary to stabilize the excited electron in other places in the fluid. We also find that even when they do not ultimately relocalize, photoexcited solvated electrons in THF temporarily visit other sites in the fluid, explaining why the photoexcitation of THF-solvated electrons is so efficient at promoting recombination with nearby scavengers. Overall, our study shows that the defining characteristic of a liquid that permits the photoassisted relocalization of solvated electrons is the existence of nascent cavities that are attractive to an excess electron; we propose that other such liquids can be found from classical computer simulations or neutron diffraction experiments.

Nonadiabatic Molecular Dynamics Simulations of Correlated Electrons in Solution. 2. A Prediction for the Observation of Hydrated Dielectrons with Pump−Probe Spectroscopy

The hydrated dielectron is a highly correlated, two-electron, solvent-supported state consisting of two spin-paired electrons confined to a single cavity in liquid water. Although dielectrons have been predicted to exist theoretically and have been used to explain the lack of ionic strength effect in the bimolecular reaction kinetics of hydrated electrons, they have not yet been observed directly. In this paper, we use the extensive nonadiabatic mixed quantum/classical excited-state molecular dynamics simulations from the previous paper to calculate the transient spectroscopy of hydrated dielectrons. Because our simulations use full configuration interaction (CI) to determine the ground and excited state two-electron wave functions at every instant, our nonequilibrium simulations allow us to compute the absorption, stimulated emission (SE), and bleach spectroscopic signals of both singlet and triplet dielectrons following excitation by ultraviolet light. Excited singlet dielectrons are predicted to display strong SE in the mid infrared and a transient absorption in the near-infrared. The near-infrared transient absorption of the singlet dielectron, which occurs near the peak of the (single) hydrated electron's equilibrium absorption, arises because the two electrons tend to separate in the excited state. In contrast, excitation of the hydrated electron gives a bleach signal in this wavelength region. Thus, our calculations suggest a clear pump-probe spectroscopic signature that may be used in the laboratory to distinguish hydrated singlet dielectrons from hydrated electrons: By choosing an excitation energy that is to the blue of the peak of the hydrated electron's absorption spectrum and probing near the maximum of the single electron's absorption, the single electron's transient bleach signal should shrink or even turn into a net absorption as sample conditions are varied to produce more dielectrons.

Mapping CTTS dynamics of Na−in tetrahydrofurane with ultrafast multichannel pump–probe spectroscopy

Using multichannel femtosecond spectroscopy we have followed Na- charge transfer to solvent (CTTS) dynamics in THF solution. Absorption of the primary photoproducts in the visible, resolved here for the first time, consists of an asymmetric triplet centered at 595 nm, which we assign to a metastable incompletely solvated neutral atomic sodium species. Decay of this feature within approximately 1 ps to a broad and structureless solvated neutral is accompanied by broadening and loss of spectral detail. Kinetic analysis shows that both the spectral structure and the decay of this band are independent of the excitation photon frequency in the range 400-800 nm. With different pump-probe polarizations the anisotropy in transient transmission has been charted and its variation with excitation wavelength surveyed. The anisotropies are assigned to the reactant bleach, indicating that due to solvent-induced symmetry breaking, the CTTS absorption band of Na- is made up of discreet orthogonally polarized sub bands. None of the anisotropy in transient absorption could be associated with the photoproduct triplet band even at the earliest measurable time delays. Along with the documented differences in the spatial distribution of ejected electrons across the tested excitation wavelength range, these results lead us to conclude that photoejection is extremely rapid, and that loss of correlations between the departing electron and its neutral core is faster than our time resolution of approximately 60 fs.

ULTRAFAST SOLVATION DYNAMICS EXPLORED BY FEMTOSECOND PHOTON ECHO SPECTROSCOPIES

Chemical reaction and optical dynamics in the liquid phase are strongly affected by specific solute-solvent interactions. The dynamical part of this coupling leads to energy fluctuations. The associated energy gap dynamics can be probed by using various nonlinear optical spectroscopies. We discuss various forms of photon echo--time-integrated, time-gated, and heterodyne-detected photon echo--as well as Fourier transform spectral interferometry. It is shown that for solutions of the dye molecule DTTCI, a system-bath correlation function can be acquired that provides a quantitative description of all (non)linear spectroscopic experiments. The deduced correlation function is projected onto the multimode Brownian oscillator model, which allows for a physical interpretation of the multiple-time correlation function and a determination of the spectral density relevant to the solvation process. The following applications of photon echo to condensed phase dynamics are discussed: enhanced vibrational mode suppression, Liouville pathways interference, and dynamical Stokes shift. Recent results of echo-peak shift experiments on the hydrated electron are also presented. The review concludes that photon echo should be useful as a novel tool to explore transition state dynamics.