Aqueous solvation dynamics with a quantum mechanical Solute: Computer simulation studies of the photoexcited hydrated electron

Hydration structures near finite-sized nanoscopic objects reconstructed using inelastic x-ray scattering measurements

Recent work has shown that it is possible to use high resolution dynamical structure factor S(q,ω) data measured with inelastic x-ray scattering to reconstruct the Green's function of water, which describes its dynamical density response to a point charge. Here, we generalize this approach and describe a strategy for reconstructing hydration behavior near simple charge distributions with excluded volumes, with the long term goal of engaging hydration processes in complex molecular systems. We use this Green's function based imaging of dynamics method to generate hydration structures and show that they are consistent with those of well-studied model systems.

The precursors of the solvated electron in methanol studied by femtosecond pump-repump-probe spectroscopy

Using UV photoionization and delayed near-infrared reexcitation pulses, a novel time-, frequency-, and polarization-resolved pump-repump-probe spectroscopy is conducted in the probing range of 450-2400 nm with improved experimental accuracy. Both the generation process and relaxation dynamics following selective repumping of intermediate species of the solvated electron are investigated and analyzed self-consistently with the help of a kinetic model. New insight in the intermediates of the trapped electron is gained leading to a unique microscopic picture.

Electron solvation in water-ammonia mixed clusters: Structure, energetics, and the nature of localization states of the excess electron

The structure and energetics of water-ammonia mixed clusters with an excess electron, [(H2O)n(NH3)m]- with m=1, n=2-6 and m=2, n=2, and also the corresponding neutral clusters are investigated in detail by means of ab initio quantum chemical calculations. The authors focus on the localization structure of the excess electron with respect to its surface versus interiorlike states, its binding to ammonia versus water molecules, the spatial and orientational arrangement of solvent molecules around the excess electron, the changes of the overall hydrogen-bonded structure of the clusters as compared to those of the neutral ones and associated dipole moment changes, vertical detachment energies of the anionic clusters, and also the vertical attachment energies of the neutral clusters. It is found that the hydrogen-bonded structure of the anionic clusters are very different from those of the neutral clusters unlike the case of water-ammonia dimer anion, and these changes in structural arrangements lead to drastically different dipole moments of the anionic and the neutral clusters. The spatial distribution of the singly occupied molecular orbital holding the excess electron shows only surface states for the smaller clusters. However, for n=5 and 6, both surface and interiorlike binding states are found to exist for the excess electron. For the surface states, the excess electron can be bound to the dangling hydrogens of either an ammonia or a water molecule with different degrees of stability and vertical detachment energies. The interiorlike states, wherever they exist, are found to have a higher vertical detachment energy than any of the surface states of the same cluster. Also, for interiorlike states, the ammonia molecule with its dangling hydrogens is always found to stay on top or on a far side of the charge density of the excess electron without participating in the hydrogen bond network of the cluster; the intermolecular hydrogen bonds are formed by the water molecules only which add to the overall stability of these anionic clusters.

Excess electron relaxation dynamics at water/air interfaces

We have performed mixed quantum-classical molecular dynamics simulations of the relaxation of a ground state excess electron at interfaces of different phases of water with air. The investigated systems included ambient water/air, supercooled water/air, Ih ice/air, and amorphous solid water/air interfaces. The present work explores the possible connections of the examined interfacial systems to finite size cluster anions and the three-dimensional infinite, fully hydrated electron. Localization site analyses indicate that in the absence of nuclear relaxation the electron localizes in a shallow potential trap on the interface in all examined systems in a diffuse, surface-bound (SB) state. With relaxation, the weakly bound electron undergoes an ultrafast localization and stabilization on the surface with the concomitant collapse of its radius. In the case of the ambient liquid interface the electron slowly (on the 10 ps time scale) diffuses into the bulk to form an interior-bound state. In each other case, the excess electron persists on the interface in SB states. The relaxation dynamics occur through distinct SB structures which are easily distinguishable by their energetics, geometries, and interactions with the surrounding water bath. The systems exhibiting the most stable SB excess electron states (supercooled water/air and Ih ice/air interfaces) are identified by their characteristic hydrogen-bonding motifs which are found to contain double acceptor-type water molecules in the close vicinity of the electron. These surface states correlate reasonably with those extrapolated to infinite size from simulated water cluster anions.

Computational studies of aqueous-phase photochemistry and the hydrated electron in finite-size clusters

A survey of recent ab initio calculations on excited electronic states of water clusters and various chromophore-water clusters is given. Electron and proton transfer processes in these systems have been characterized by the determination of electronic wave functions, minimum-energy reaction paths and potential-energy profiles. It is pointed out that the transfer of a neutral hydrogen atom (leading to biradicals) rather than the transfer of a proton (leading to ion pairs) is the generic excited-state reaction mechanism in these systems. The hydrated hydronium radical, (H3O)(aq), plays a central role in this scenario. The electronic and vibrational spectra of H3O(H2O)(n) clusters and the decay mechanism of these metastable species have been investigated in some detail. The results suggest that (H3O)(aq) could be the carrier of the characteristic spectroscopic properties of the hydrated electron in liquid water.

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.

Electron hydration dynamics in water clusters: A direct ab initio molecular dynamics approach

Electron attachment dynamics of excess electron in water cluster (H2O)n (n = 2 and 3) have been investigated by means of full-dimensional direct ab initio molecular dynamics (MD) method at the MP26-311++G(d,p) level. It was found that the hydrogen bond breaking due to the excess electron is an important process in the first stage of electron capture in water trimer. Time scale of electron localization and hydrogen bond breaking were determined by the direct ab initio MD simulation. The initial process of hydration in water cluster is clearly visualized in the present study. In n = 3, an excess electron is first trapped around the cyclic water trimer with a triangular form, where the excess electron is equivalently distributed on the three water molecules at time zero. After 50 fs, the excess electron is concentrated into two water molecules, while the potential energy of the system decreases by -1.5 kcal/mol from the vertical point. After 100 fs, the excess electron is localized in one of the water molecules and the potential energy decreases by -5.3 kcal/mol, but the triangular form still remained. After that, one of the hydrogen bonds in the triangular form is gradually broken by the excess electron, while the structure becomes linear at 100-300 fs after electron capture. The time scale of hydrogen bond breaking due to the excess electron is calculated to be about 300 fs. Finally, a dipole bound state is formed by the linear form of three water molecules. In the case of n = 2, the dipole bound anion is formed directly. The mechanism of electron hydration dynamics was discussed on the basis of theoretical results.

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.

Projections of quantum observables onto classical degrees of freedom in mixed quantum-classical simulations: understanding linear response failure for the photoexcited hydrated electron

We present a general analytic method for understanding how specific motions of a classical bath influence the dynamics of quantum-mechanical observables in mixed quantum-classical molecular dynamics simulations. We apply our method and develop expressions for the special case of quantum solvation, allowing us to examine how specific classical solvent motions couple to the equilibrium energy fluctuations and nonequilibrium energy relaxation of a quantum-mechanical solute. As a first application of our formalism, we investigate the motions of classical water underlying the equilibrium and nonequilibrium excited-state solvent response functions of the hydrated electron; the results allow us to explain why the linear response approximation fails for this system.

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.

Nonadiabatic Molecular Dynamics Simulations of Correlated Electrons in Solution. 1. Full Configuration Interaction (CI) Excited-State Relaxation Dynamics of Hydrated Dielectrons

The hydrated dielectron is composed of two excess electrons dissolved in liquid water that occupy a single cavity; in both its singlet and triplet spin states there is a significant exchange interaction so the two electrons cannot be considered to be independent. In this paper and the following paper,we present the results of mixed quantum/classical molecular dynamics simulations of the nonadiabatic relaxation dynamics of photoexcited hydrated dielectrons, where we use full configuration interaction (CI) to solve for the two-electron wave function at every simulation time step. To the best of our knowledge, this represents the first systematic treatment of excited-state solvation dynamics where the multiple-electron problem is solved exactly. The simulations show that the effects of exchange and correlation contribute significantly to the relaxation dynamics. For example, spin-singlet dielectrons relax to the ground state on a time scale similar to that of single electrons excited at the same energy, but spin-triplet dielectrons relax much faster. The difference in relaxation dynamics is caused by exchange and correlation: The Pauli exclusion principle imposes very different electronic structure when the electrons' spins are singlet paired than when they are triplet paired, altering the available nonadiabatic relaxation pathways. In addition, we monitor how electronic correlation changes dynamically during nonadiabatic relaxation and show that solvent dynamics cause electron correlation to evolve quite differently for singlet and triplet dielectrons. Despite such differences, our calculations show that both spin states are stable to excited-state dissociation, but that the excited-state stability has different origins for the two spin states. For singlet dielectrons, the stability depends on whether the solvent structure can rearrange to create a second cavity before the ground state is reached. For triplet dielectrons, in contrast, electronic correlation ensures that the two electrons do not dissociate, even if the dielectron is artificially kept from reaching the ground state. In addition, both singlet and triplet dielectrons change shape dramatically during relaxation, so that linear response fails to describe the solvation dynamics for either spin state. In the following paper (Larsen, R. E.; Schwartz, B. J. J. Phys. Chem. B 2006, 110, 9692), we use these simulations to calculate the pump-probe spectroscopic signal expected for photoexcited hydrated dielectrons and to predict an experiment to observe hydrated dielectrons directly.

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.

Theory of electron solvation in polar liquids: A continuum model

The solvation of electrons in polar liquids is analyzed on the basis of an extended continuum model. In addition to the long-range electron-dipole interaction two short-range interactions are introduced. Among others one accounts for interactions with groups capable of forming hydrogen bonds and the second for quadrupolar characteristics of the liquid molecules. Both are induced by the orientation of the molecular dipole. Applying the scaling method a proper reaction coordinate is introduced and the solvation dynamics are discussed for the electron in the electronic ground state and after excitation to the p-type excited state. The observed spectral evolution of the transient absorption spectra, after two photon excitations for electrons in water and in methanol, is well described by this theory. An analytic estimate for the nonradiative deactivation from the electronically excited solvated electron is found to be consistent with an observed lifetime of 50 fs for the electron in water. The theory predicts an about three times slower internal conversion in methanol as solvent in comparison with water.

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.

Characterization of excess electrons in water-cluster anions by quantum simulations

Water-cluster anions can serve as a bridge to understand the transition from gaseous species to the bulk hydrated electron. However, debate continues regarding how the excess electron is bound in (H2O)-n, as an interior, bulklike, or surface electronic state. To address the uncertainty, the properties of (H2O)-n clusters with 20 to 200 water molecules have been evaluated by mixed quantum-classical simulations. The theory reproduces every observed energetic, spectral, and structural trend with cluster size that is seen in experimental photoelectron and optical absorption spectra. More important, surface states and interior states each manifest a characteristic signature in the simulation data. The results strongly support assignment of surface-bound electronic states to the water-cluster anions in published experimental studies thus far.

The role of solvent structure in the absorption spectrum of solvated electrons: Mixed quantum/classical simulations in tetrahydrofuran

In polar fluids such as water and methanol, the peak of the solvated electron's absorption spectrum in the red has been assigned as a sum of transitions between an s-like ground state and three nearly degenerate p-like excited states bound in a quasispherical cavity. In contrast, in weakly polar solvents such as tetrahydrofuran (THF), the solvated electron has an absorption spectrum that peaks in the mid-infrared, but no definitive assignment has been offered about the origins of the spectrum or the underlying structure. In this paper, we present the results of adiabatic mixed quantum/classical molecular dynamic simulations of the solvated electron in THF, and provide a detailed explanation of the THF-solvated electron's absorption spectrum and electronic structure. Using a classical solvent model and a fully quantum mechanical excess electron, our simulations show that although the ground and first excited states are bound in a quasispherical cavity, a multitude of other, nearby solvent cavities support numerous, nearly degenerate, bound excited states that have little Franck-Condon overlap with the ground state. We show that these solvent cavities, which are partially polarized so that they act as electron trapping sites, are an inherent property of the way THF molecules pack in the liquid. The absorption spectrum is thus assigned to a sum of bound-to-bound transitions between a localized ground state and multiple disjoint excited states scattered throughout the fluid. Furthermore, we find that the usual spherical harmonic labels (e.g., s-like, p-like) are not good descriptors of the excited-state wave functions of the solvated electron in THF. Our observation of multiple disjoint excited states is consistent with femtosecond pump-probe experiments in the literature that suggest that photoexcitation of solvated electrons in THF causes them to relocalize into solvent cavities far from where they originated.

Hydrated Electron Dynamics: From Clusters to Bulk

The electronic relaxation dynamics of size-selected (H2O)n-/(D2O)n[25 </= n </= 50] clusters have been studied with time-resolved photoelectron imaging. The excess electron (ec-) was excited through the ec-(p)<--ec-(s) transition with an ultrafast laser pulse, with subsequent evolution of the excited state monitored with photodetachment and photoelectron imaging. All clusters exhibited p-state population decay with concomitant s-state repopulation (internal conversion) on time scales ranging from 180 to 130 femtoseconds for (H2O)n- and 400 to 225 femtoseconds for (D2O)n-; the lifetimes decrease with increasing cluster sizes. Our results support the "nonadiabatic relaxation" mechanism for the bulk hydrated electron (eaq-), which invokes a 50-femtosecond eaq-(p)-->eaq-(s(dagger)) internal conversion lifetime.

Critical evaluation of approximate quantum decoherence rates for an electronic transition in methanol solution

We present a quantum molecular dynamics calculation of a semiclassical decoherence function to evaluate the accuracy of alternative short-time approximations for coherence loss in the dynamics of condensed phase electronically nonadiabatic processes. The semiclassical function from mixed quantum-classical molecular dynamics simulations and frozen Gaussian wave packets is computed for the electronic transition of an excited state excess electron to the ground state in liquid methanol. The decoherence function decays on a 10 fs time scale that is qualitatively similar to the aqueous case. We demonstrate that it is the motion of the hydrogen atom, and, in particular, the hydrogen rotation around the oxygen-methyl bond which is predominantly responsible for destroying the quantum correlations between alternative states. Multiple time scales due to the slower diffusive nuclear modes, which dominate the solvation response of methanol, do not contribute to the coherence loss. The choice of the coordinate representation is investigated in detail and concluded to be irrelevant to the decay. Changes in both nuclear momenta and positions on the two alternative potential surfaces are found to contribute to decoherence, the former dominating at short times (t < 5 fs), the latter controlling the decay at longer times. Various short-time approximations to the full dynamics for the decoherence function are tested for the first time. The present treatment rigorously develops the short-time description and establishes its range of validity. Whereas the lowest-order short-time approximation proves to be a very good approximation up to about 5 fs, we also find that it bounds the decay of the decoherence function. After 5 fs, the coherence decay in fact becomes faster than the single Gaussian predicted in the lowest-order short-time limit. This decay is well reflected by an enhanced low-order approximation, which is also easily computed from equilibrium classical forces.