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.     

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.     

Dielectron attachment and hydrogen evolution reaction in water clusters

Binding of excess electrons to nanosize water droplets, with a focus on the hitherto largely unexplored properties of doubly-charged clusters, were investigated experimentally using mass spectrometry and theoretically with large-scale first-principles simulations based on spin-density-functional theory, with all the valence electrons (that is, 8e per water molecule) and excess electrons treated quantum mechanically. Singly-charged clusters (H(2)O)(n)(-1) were detected for n = 6-250, and our calculated vertical detachment energies agree with previously measured values in the entire range 15 ≤ n ≤ 105, giving a consistent interpretation in terms of internal, surface and diffuse states of the excess electron. Doubly-charged clusters were measured in the range of 83 ≤ n ≤ 123, with (H(2)O)(n)(-2) clusters found for 83 ≤ n < 105, and mass-shifted peaks corresponding to (H(2)O)(n-2)(OH(-))(2) detected for n ≥ 105. The simulations revealed surface and internal dielectron, e(-)(2), localization modes and elucidated the mechanism of the reaction (H(2)O)(n)(-2) → (H(2)O)(n-2) (OH(-))(2) + H(2) (for n ≥ 105), which was found to occur via concerted approach of a pair of protons belonging to two water molecules located in the first shell of the dielectron internal hydration cavity, culminating in formation of a hydrogen molecule 2H(+) + e(-)(2) → H(2). Instability of the dielectron internal localization impedes the reaction for smaller (n < 105) doubly-charged clusters.     

Hydrated Electron Transfer to Nucleobases in Aqueous Solutions Revealed by Ab Initio Molecular Dynamics Simulations

We present an ab initio molecular dynamics (AIMD) simulation study into the transfer dynamics of an excess electron from its cavity‐shaped hydrated electron state to a hydrated nucleobase (NB)‐bound state. In contrast to the traditional view that electron localization at NBs (G/A/C/T), which is the first step for electron‐induced DNA damage, is related only to dry or prehydrated electrons, and a fully hydrated electron no longer transfers to NBs, our AIMD simulations indicate that a fully hydrated electron can still transfer to NBs. We monitored the transfer dynamics of fully hydrated electrons towards hydrated NBs in aqueous solutions by using AIMD simulations and found that due to solution‐structure fluctuation and attraction of NBs, a fully hydrated electron can transfer to a NB gradually over time. Concurrently, the hydrated electron cavity gradually reorganizes, distorts, and even breaks. The transfer could be completed in about 120–200 fs in four aqueous NB solutions, depending on the electron‐binding ability of hydrated NBs and the structural fluctuation of the solution. The transferring electron resides in the π*‐type lowest unoccupied molecular orbital of the NB, which leads to a hydrated NB anion. Clearly, the observed transfer of hydrated electrons can be attributed to the strong electron‐binding ability of hydrated NBs over the hydrated electron cavity, which is the driving force, and the transfer dynamics is structure‐fluctuation controlled. This work provides new insights into the evolution dynamics of hydrated electrons and provides some helpful information for understanding the DNA‐damage mechanism in solution.     

The structure of the hydrated electron. Part 2. A mixed quantum/classical molecular dynamics embedded cluster density functional theory: single-excitation configuration interaction study

Adiabatic mixed quantum/classical (MQC) molecular dynamics (MD) simulations were used to generate snapshots of the hydrated electron in liquid water at 300 K. Water cluster anions that include two complete solvation shells centered on the hydrated electron were extracted from the MQC MD simulations and embedded in a roughly 18 Ax18 Ax18 A matrix of fractional point charges designed to represent the rest of the solvent. Density functional theory (DFT) with the Becke-Lee-Yang-Parr functional and single-excitation configuration interaction (CIS) methods were then applied to these embedded clusters. The salient feature of these hybrid DFT(CIS)/MQC MD calculations is significant transfer (approximately 18%) of the excess electron's charge density into the 2p orbitals of oxygen atoms in OH groups forming the solvation cavity. We used the results of these calculations to examine the structure of the singly occupied and the lower unoccupied molecular orbitals, the density of states, the absorption spectra in the visible and ultraviolet, the hyperfine coupling (hfcc) tensors, and the infrared (IR) and Raman spectra of these embedded water cluster anions. The calculated hfcc tensors were used to compute electron paramagnetic resonance (EPR) and electron spin echo envelope modulation (ESEEM) spectra for the hydrated electron that compared favorably to the experimental spectra of trapped electrons in alkaline ice. The calculated vibrational spectra of the hydrated electron are consistent with the red-shifted bending and stretching frequencies observed in resonance Raman experiments. In addition to reproducing the visible/near IR absorption spectrum, the hybrid DFT model also accounts for the hydrated electron's 190-nm absorption band in the ultraviolet. Thus, our study suggests that to explain several important experimentally observed properties of the hydrated electron, many-electron effects must be accounted for: one-electron models that do not allow for mixing of the excess electron density with the frontier orbitals of the first-shell solvent molecules cannot explain the observed magnetic, vibrational, and electronic properties of this species. Despite the need for multielectron effects to explain these important properties, the ensemble-averaged radial wavefunctions and energetics of the highest occupied and three lowest unoccupied orbitals of the hydrated electrons in our hybrid model are close to the s- and p-like states obtained in one-electron models. Thus, one-electron models can provide a remarkably good approximation to the multielectron picture of the hydrated electron for many applications; indeed, the two approaches appear to be complementary.     

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.     

Molecular dynamics simulations of electron-alkali cation pairs in bulk water

The structural, dynamic, and thermodynamic properties of an excess electron interacting with an alkali cation (Na+, K+, Li+) in bulk water were investigated by means of a mixed quantum-classical molecular dynamics simulation technique. This study includes a reparametrization of the electron-cation pseudopotentials. The free energy calculations for all three systems show that a contact electron-cation pair can be observed, which is either as stable as the dissociated pair (Li+) or more stable by only a few kT (Na+, K+). Given that the dissociation barrier is also quite small, we suggest that the average cation-electron distance in the experiments at room temperature will not depend on this free energy profile but rather on the minimization of the Coulombic repulsive interaction between like charges in the solvent medium. This enables us to compare the present molecular dynamics simulations with the spectroscopic data obtained for different ionic strengths. The overall trend of the UV-vis hydrated absorption spectra, namely, the shift toward shorter wavelengths at high ionic strengths, is fairly well reproduced. This confirms our hypothesis of statistical distribution of the cations and solvated electrons.     

Solvation of excess electrons in LiF ionic pair matrix: evidence for a solvated dielectron from ab initio molecular dynamics simulations and calculations

Ab initio molecular dynamics simulations and first-principles calculations reveal the existence of a solvated dielectron species, (2e)s, in an LiF ionic matrix. The nature of the solvation mechanism and the stability of the species was explored. In addition to electrostatic interactions, a hole-orbital coupling among solvent molecules may significantly enhance the stability of the solvated electrons and govern the extent of electron solvation. This hole-orbital coupling is different from either an electrostatic coupling or conventional chemical bonding, and it may be described as a transition between them.     

Model potential approaches for describing the interaction of excess electrons with water clusters: incorporation of long-range correlation effects

In this work we focus on the binding of excess electrons to water clusters, a problem for which dispersion interactions, which originate from long-range correlation effects, are especially important. Two different model potential approaches, one using quantum Drude oscillators and the other using polarization potentials, are investigated for describing the long-range correlation effects between the weakly bound excess electron and the more tightly bound electrons of the monomers. We show that these two approaches are related in that the polarization potential models can be derived from the quantum Drude model approach by use of an adiabatic separation between the excess electron and the Drude oscillators. The model potential approaches are applied to clusters containing up to 45 water monomers. Where possible, comparison is made with the results of ab initio electronic structure calculations. Overall, the polarization potential approach is found to give electron binding energies in good agreement with those from the Drude model and ab initio calculations, with the greatest discrepancies being found for "cavity-bound" anion states.     

Second-order post-Hartree-Fock perturbation theory for the electron current

Based on the super-fermion representation of quantum kinetic equations we develop nonequilibrium, post-Hartree-Fock many-body perturbation theory for the current through a region of interacting electrons. We apply the theory to out of equilibrium Anderson model and discuss practical implementation of the approach. Our calculations show that nonequilibrium electronic correlations may produce significant quantitative and qualitative corrections to mean-field electronic transport properties.     

Electron localization and radiation chemistry of amides

Alkylamides (such as N,N'-dimethylformamide, N,N'-diethylformamide, and N,N'-dimethylacetamide) are aprotic solvents that are widely used in organic synthesis. These polar molecules have no electron affinity, and it is believed that irradiated liquid and solid amides stabilize excess electrons as cavity-type species analogous to hydrated and ammoniated electrons. In this study, we use isotope substitution and EPR spectroscopy to demonstrate that, in frozen amides, the suspected "cavity electron" is, in fact, a solvent-stabilized monomer anion. Our observations call into question other attributions of such features in the literature, both in low temperature solids and room temperature liquids. We also provide a general scheme describing amide radiolysis, as the related amides are used as metal ion extracting agents in nuclear separations.     

Independent pairs and Monte-Carlo simulations of the geminate recombination of solvated electrons in liquid-to-supercritical water

Independent pairs (IP) and Monte Carlo (MC) simulations are employed to model experimental femtosecond time-resolved pump-probe spectroscopic data on the geminate recombination dynamics of solvated electrons in liquid-to-supercritical water. The hydrated electron was created by two-photon ionization of the neat fluid with a total ionization energy of 9.3 eV. In both numerical approaches, the ejection length, , (i.e. the distance from the ionization core, at which the electron is thermally and spatially localized) is used as the primary adjustable fitting parameter that can bring both model simulations into quantitative agreement with the ultrafast kinetic experiment. The influence of the thermodynamic conditions on the ejection length and on the recombination mechanism is discussed. Whereas in the compressed liquid associated with a high dielectric constant (ε ≥ 20), the electron recombines predominantly with the OH radical, the dissociative recombination via charge neutralization with the hydronium cation takes over at small dielectric constants (ε < 20). The importance of charge-dipole interactions for Monte-Carlo simulations of the recombination reactions of the hydrated electrons in the low-permittivity region is stressed.     

Ammoniated electron as a solvent stabilized multimer radical anion

The excess electron in liquid ammonia ("ammoniated electron") is commonly viewed as a cavity electron in which the s-type wave function fills the interstitial void between 6 and 9 ammonia molecules. Here we examine an alternative model in which the ammoniated electron is regarded as a solvent stabilized multimer radical anion in which most of the excess electron density resides in the frontier orbitals of N atoms in the ammonia molecules forming the solvation cavity. The cavity is formed due to the repulsion between negatively charged solvent molecules. Using density functional theory calculations, we demonstrate that such core anions would semiquantitatively account for the observed pattern of Knight shifts for 1H and 14N nuclei observed by NMR spectroscopy and the downshifted stretching and bending modes observed by infrared spectroscopy. We speculate that the excess electrons in other aprotic solvents might be, in this respect, analogous to the ammoniated electron, with substantial transfer of the spin density into the frontier N and C orbitals of methyl, amino, and amide groups.     

Nuclear quantum effects on the nonadiabatic decay mechanism of an excited hydrated electron

We present a kinetic analysis of the nonadiabatic decay mechanism of an excited state hydrated electron to the ground state. The theoretical treatment is based on a quantized, gap dependent golden rule rate constant formula which describes the nonadiabatic transition rate between two quantum states. The rate formula is expressed in terms of quantum time correlation functions of the energy gap and of the nonadiabatic coupling. These gap dependent quantities are evaluated from three different sets of mixed quantum-classical molecular dynamics simulations of a hydrated electron equilibrated (a) in its ground state, (b) in its first excited state, and (c) on a hypothetical mixed potential energy surface which is the average of the ground and the first excited electronic states. The quantized, gap dependent rate results are applied in a phenomenological kinetic equation which provides the survival probability function of the excited state electron. Although the lifetime of the equilibrated excited state electron is computed to be very short (well under 100 fs), the survival probability function for the nonequilibrium process in pump-probe experiments yields an effective excited state lifetime of around 300 fs, a value that is consistent with the findings of several experimental groups and previous theoretical estimates.     

Analysis of localization sites for an excess electron in neutral methanol clusters using approximate pseudopotential quantum-mechanical calculations

We have used a recently developed electron-methanol molecule pseudopotential in approximate quantum mechanical calculations to evaluate and statistically analyze the physical properties of an excess electron in the field of equilibrated neutral methanol clusters ((CH(3)OH)(n), n=50-500). The methanol clusters were generated in classical molecular dynamics simulations at nominal 100 and 200 K temperatures. Topological analysis of the neutral clusters indicates that methyl groups cover the surface of the clusters almost exclusively, while the associated hydroxyl groups point inside. Since the initial neutral clusters are lacking polarity on the surface and compact inside, the excess electron can barely attach to these structures. Nevertheless, most of the investigated cluster configurations do support weakly stabilized cluster anion states. We find that similarly to water clusters, the pre-existing instantaneous dipole moment of the neutral clusters binds the electron. The localizing electrons occupy diffuse, weakly bound surface states that largely engulf the cluster although their centers are located outside the cluster molecular frame. The initial localization of the excess electron is reflected in its larger radius compared to water due to the lack of free OH hydrogens on the cluster surface. The stabilization of the excess electron increases, while the radius decreases monotonically as the clusters grow in size. Stable, interior bound states of the excess electron are not observed to form neither in finite size methanol clusters nor in the equilibrium bulk.     

Molecular Dynamics Insights into Electron-Catalyzed Dissociation Repair of Cyclobutane Pyrimidine Dimer

Excess electrons are not only an important source of radiation damage, but also participate in the repair process of radiation damage such as cyclobutane pyrimidine dimer (CPD). Using ab initio molecular dynamics (AIMD) simulations, we reproduce the single excess electron stepwise catalytic CPD dissociation process in detail with an emphasis on the energy levels and molecular structure details associated with excess electrons. On the basis of the AIMD simulations on the CPD aqueous solution with two vertically added excess electrons, we exclude the early-proposed [2+2]-like concerted synchronous dissociation mechanism, and analyze the difference between the symmetry of the actual reaction and the symmetry of the frontier molecular orbitals which deeply impact the mechanism. Importantly, we propose a new model of the stepwise electron-catalyzed dissociation mechanism that conforms to the reality. This work not only provides dynamics insights into the excess electron catalyzed dissociation mechanism, but also reveals different roles of two excess electrons in two bond-cleavage steps (promoting versus inhibiting).     

The solvation of two electrons in the gaseous clusters of Na- (NH3)n and Li- (NH3)n

Alkali metal ammonia clusters, in their cationic, neutral, and anionic form, are molecular models for the alkali-ammonia solutions, which have rich variation of phases with the solvated electrons playing an important role. With two s electrons, the Na(-)(NH(3))(n) and Li(-)(NH(3))(n) clusters are unique in that they capture the important aspect of the coupling between two solvated electrons. By first principles calculations, we demonstrate that the two electrons are detached from the metal by n = 10, which produces a cluster with a solvated electron pair in the vicinity of a solvated alkali cation. The coupling of the two electrons leads to either the singlet or triplet state, both of which are stable. They are also quite distinct from the hydrated anionic clusters Na(-)(H(2)O)(n) and Li(-)(H(2)O)(n), in that the solvated electrons are delocalized and widely distributed among the solvent ammonia molecules. The Na(-)(NH(3))(n) and Li(-)(NH(3))(n) series, therefore, provide another interesting type of molecular model for the investigation of solvated electron pairs.     

Evolution of lone pair of excess electrons inside molecular cages with the deformation of the cage in e2@C60F60 systems

For unusual e(2)@C(60)F(60)(I(h), D(6h), and D(5d)) cage structures with two excess electrons, it is reported that not only the lone pair in singlet state but also two single excess electrons in triplet state can be encapsulated inside the C(60)F(60) cages to form single molecular solvated dielectrons. The interesting relationship between the shape of the cage and the spin state of the system has revealed that ground states are singlet state for spherical shaped e(2)@C(60)F(60)(I(h)) and triplet states for short capsular shaped e(2)@C(60)F(60)(D(6h)) and long capsular shaped e(2)@C(60)F(60)(D(5d)), which shows a spin evolution from the singlet to triplet state with the deformation of the cage from spherical to capsular shape. For these excess electron systems, the three ground state structures have large vertical electron detachment energies (VDEs (I) of 1.720-2.283 eV and VDEs (II) of 3.959-5.288 eV), which shows their stabilities and suggests that the large C(60)F(60) cage is the efficient container of excess electrons.     

Effect of temperature and water content on the shear viscosity of the ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide as studied by atomistic simulations

Atomistic simulations are conducted to examine the dependence of the viscosity of 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide on temperature and water content. A nonequilibrium molecular dynamics procedure is utilized along with an established fixed charge force field. It is found that the simulations quantitatively capture the temperature dependence of the viscosity as well as the drop in viscosity that occurs with increasing water content. Using mixture viscosity models, we show that the relative drop in viscosity with water content is actually less than that that would be predicted for an ideal system. This finding is at odds with the popular notion that small amounts of water cause an unusually large drop in the viscosity of ionic liquids. The simulations suggest that, due to preferential association of water with anions and the formation of water clusters, the excess molar volume is negative. This means that dissolved water is actually less effective at lowering the viscosity of these mixtures when compared to a solute obeying ideal mixing behavior. The use of a nonequilibrium simulation technique enables diffusive behavior to be observed on the time scale of the simulations, and standard equilibrium molecular dynamics resulted in sub-diffusive behavior even over 2 ns of simulation time.