Dynamics of electron localization in warm versus cold water clusters

From unbound to bound states: Ab initio molecular dynamics of ammonia clusters with an excess electron

Ab initio molecular dynamics simulations of negatively charged clusters of 2-48 ammonia molecules were performed to elucidate the electronic stability of the excess electron as a function of cluster size. We show that while the electronic stability of finite temperature clusters increases with cluster size, as few as 5-7 ammonia molecules can bind an excess electron, reaching a vertical binding energy slightly less than half of the bulk value for the largest system studied. These results, which are in agreement with previous studies wherever available, allowed us to analyze the excess electron binding patterns in terms of its radius of gyration and shape anisotropy and provide a qualitative interpretation based on a particle-in-a-spherical-well model.

Understanding the Temperature Dependence and Finite Size Effects in Ab Initio MD Simulations of the Hydrated Electron

The hydrated electron is of interest to both theorists and experimentalists as a paradigm solution-phase quantum system. Although the bulk of the theoretical work studying the hydrated electron is based on mixed quantum/classical (MQC) methods, recent advances in computer power have allowed several attempts to study this object using ab initio methods. The difficulty with employing ab initio methods for this system is that even with relatively inexpensive quantum chemistry methods such as density functional theory (DFT), such calculations are still limited to at most a few tens of water molecules and only a few picoseconds duration, leaving open the question as to whether the calculations are converged with respect to either system size or dynamical fluctuations. Moreover, the ab initio simulations of the hydrated electron that have been published to date have provided only limited analysis. Most works calculate the electron's vertical detachment energy, which can be compared to experiment, and occasionally the electronic absorption spectrum is also computed. Structural features, such as pair distribution functions, are rare in the literature, with the majority of the structural analysis being simple statements that the electron resides in a cavity, which are often based only on a small number of simulation snapshots. Importantly, there has been no ab initio work examining the temperature-dependent behavior of the hydrated electron, which has not been satisfactorily explained by MQC simulations. In this work, we attempt to remedy this situation by running DFT-based ab initio simulations of the hydrated electron as a function of both box size and temperature. We show that the calculated properties of the hydrated electron are not converged even with simulation sizes up to 128 water molecules and durations of several tens of picoseconds. The simulations show significant changes in the water coordination and solvation structure with box size. Our temperature-dependent simulations predict a red-shift of the absorption spectrum (computed using TD-DFT with an optimally tuned range-separated hybrid functional) with increasing temperature, but the magnitude of the predicted red-shift is larger than that observed experimentally, and the absolute position of the calculated spectra are off by over half an eV. The spectral red-shift at high temperatures is accompanied by both a partial loss of structure of the electron's central cavity and an increased radius of gyration that pushes electron density onto and beyond the first solvation shell. Overall, although ab initio simulations can provide some insights into the temperature-dependent behavior of the hydrated electron, the simulation sizes and level of quantum chemistry theory that are currently accessible are inadequate for correctly describing the experimental properties of this fascinating object.

Degradation of Per- and Polyfluoroalkyl Substances with Hydrated Electrons: A New Mechanism from First-Principles Calculations

Per- and polyfluoroalkyl substances (PFASs) are synthetic contaminants found in drinking groundwater sources and a wide variety of consumer products. Because of their adverse environmental and human health effects, remediation of these persistent compounds has attracted significant recent attention. To gain mechanistic insight into their remediation, we present the first ab initio study of PFAS degradation via hydrated electrons─a configuration that has not been correctly considered in previous computational studies up to this point. To capture these complex dynamical effects, we harness ab initio molecular dynamics (AIMD) simulations to probe the reactivities of perfluorooctanoic (PFOA) and perfluorooctane sulfonic acid (PFOS) with hydrated electrons in explicit water. We complement our AIMD calculations with advanced metadynamics sampling techniques to compute free energy profiles and detailed statistical analyses of PFOA/PFOS dynamics. Although our calculations show that the activation barrier for C-F bond dissociation in PFOS is three times larger than that in PFOA, all the computed free energy barriers are still relatively low, resulting in a diffusion-limited process. We discuss our results in the context of recent studies on PFAS degradation with hydrated electrons to give insight into the most efficient remediation strategies for these contaminants. Most importantly, we show that the degradation of PFASs with hydrated electrons is markedly different from that with excess electrons/charges, a common (but largely incomplete) approach used in several earlier computational studies.

Benzene Radical Anion Microsolvated in Ammonia Clusters: Modeling the Transition from an Unbound Resonance to a Bound Species

The benzene radical anion, well-known in organic chemistry as the first intermediate in the Birch reduction of benzene in liquid ammonia, exhibits intriguing properties from the point of view of quantum chemistry. Notably, it has the character of a metastable shape resonance in the gas phase, while measurements in solution find it to be experimentally detectable and stable. In this light, our previous calculations performed in bulk liquid ammonia explicitly reveal that solvation leads to stabilization. Here, we focus on the transition of the benzene radical anion from an unstable gas-phase ion to a fully solvated bound species by explicit ionization calculations of the radical anion solvated in molecular clusters of increasing size. The computational cost of the largest systems is mitigated by combining density functional theory with auxiliary methods including effective fragment potentials or approximating the bulk by polarizable continuum models. Using this methodology, we obtain the cluster size dependence of the vertical binding energy of the benzene radical anion converging to the value of -2.3 eV at a modest computational cost.

Ab Initio Molecular Dynamics Simulations of Solvated Electrons in Ammonia Clusters

We investigated excess electron solvation dynamics in (NH₃)_n^- ammonia clusters in the n = 8-32 size range by performing finite temperature molecular dynamics simulations. In particular, we focused on three possible scenarios. The first case is designed to model electron attachment to small neutral ammonia clusters (n ≤ ∼10) that form hydrogen-bonded chains. The excess electron is bound to the clusters via dipole bound states, and persists with a VDE of ∼160 meV at 100 K for the n = 8 cluster. The coupled nuclear and electronic relaxation is fast (within ∼100 fs) and takes place predominantly by intermolecular librational motions and the intramolecular umbrella mode. The second scenario illustrates the mechanism of excess electron attachment to cold compact neutral clusters of medium size (8 ≤ n ≤ 32). The neutral clusters show increasing tendency with size to bind the excess electron on the surface of the clusters in weakly bound, diffuse, and highly delocalized states. Anionic relaxation trajectories launched from these initial states provide no indication for excess electron stabilization for sizes n < 24. Excess electrons are likely to autodetach from these clusters. The two largest investigated clusters (n = 24 and 32) can accommodate the excess electron in electronic states that are mainly localized on the surface, but they may be partly embedded in the cluster. In the last 500 fs of the simulated trajectories, the VDE of the solvated electron fluctuates around ∼200 meV for n = 24 and ∼500 meV for n = 32, consistent with the values extrapolated from the experimentally observed linear VDE-n^-1/3 trend. In the third case, we prepared neutral ammonia cluster configurations, including an n = 48 cluster, that contain possible electron localization sites within the interior of the cluster. Excess electrons added to these clusters localize in cavities with high VDEs up to 1.9 eV for n = 48. The computed VDE values for larger clusters are considerably higher than the experimentally observed photoelectric threshold energy for the solvated electron in bulk ammonia (∼1.4 eV). Molecular dynamics simulations launched from these initial cavity states strongly indicate, however, that these cavity structures exist only for ∼200 fs. During the relaxation the electron leaves the cavity and becomes delocalized, while the cluster loses its initial compactness.

Excess electron solvation in ammonia clusters

We performed a combination of quantum chemical calculations and molecular dynamics simulations to assess the stability of various size NH_{3 n} ^- ammonia cluster anions up to n = 32 monomers. In the n = 3-8 size range, cluster anions are optimized and the vertical detachment energy of the excess electron (VDE) from increasing size clusters is computed using various level methods including density functional theory, MP2, and coupled-cluster singles doubles with perturbative triples. These clusters bind the electrons in nonbranched hydrogen bonding chains in dipole bound states. The VDE increases with size from a few millielectron volt up to ∼200 meV. The electron binding energy is weaker than that in water clusters but comparable to small methanol cluster VDEs. We located the first branched hydrogen bonding cluster that binds the excess electron at n = 7. For larger (n = 8-32) clusters, we generated cold, neutral clusters by semiempirical and ab initio molecular dynamics simulations and added an extra electron to selected neutral configurations. VDE calculations on the adiabatic and the relaxed anionic structures suggest that the n = 12-32 neutral clusters weakly bind the excess electron. Electron binding energies for these clusters (∼100 meV) appear to be significantly weaker than those extrapolated from experimental data. The observed excess electron states are diffuse and localized outside the molecular frame (surface states) with minor (∼1%) penetration to the nitrogen frontier orbitals. Stable minima with excess electron states surrounded by solvent molecules (cavity states) were not found in this size regime.

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.

Thermal Equilibration Controls H-Bonding and the Vertical Detachment Energy of Water Cluster Anions

One of the outstanding puzzles in the photoelectron spectroscopy of water anion clusters, which serve as precursors to the hydrated electron, is that the excess electron has multiple vertical detachment energies (VDEs), with different groups seeing different distributions of VDEs. We have studied the photoelectron spectroscopy of water cluster anions using simulation techniques designed to mimic the different ways that water cluster anions are produced experimentally. Our simulations take advantage of density functional theory-based Born-Oppenheimer molecular dynamics with an optimally tuned range-separated hybrid functional that is shown to give outstanding accuracy for calculating electron binding energies for this system. We find that our simulations are able to accurately reproduce the experimentally observed VDEs for cluster anions of different sizes, with different VDE distributions observed depending on how the water cluster anions are prepared. For cluster anion sizes up to 20 water molecules, we see that the excess electron always resides on the surface of the cluster and that the different discrete VDEs result from the discrete number of hydrogen bonds made to the electron by water molecules on the surface. Clusters that are less thermally equilibrated have surface waters that tend to make single H-bonds to the electron, resulting in lower VDEs, while clusters that are more thermally equilibrated have surface waters that prefer to make two H-bonds to the electron, resulting in higher VDEs.

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.

Ab initio molecular dynamics study of solvated electrons in methanol clusters

Stable surface excess electronic states in small methanol cluster anions were identified and characterized in ab initio molecular dynamics simulations.

The mechanism for the formation of OH radicals in condensed-phase water under ultraviolet irradiation

OH radicals can be produced via direct water photolysis through concerted proton and electron transfer.

Hydrated Electrons in Water Clusters: Inside or Outside, Cavity or Noncavity?

In this work, we compare the applicability of three electron–water molecule pseudopotentials in modeling the physical properties of hydrated electrons. Quantum model calculations illustrate that the recently suggested Larsen–Glover–Schwartz (LGS) model and its modified m-LGS version have a too-attractive potential in the vicinity of the oxygen. As a result, LGS models predict a noncavity hydrated electron structure in clusters at room temperature, as seen from mixed one-electron quantum–classical molecular dynamics simulations of water cluster anions, with the electron localizing exclusively in the interior of the clusters. Comparative calculations using the cavity-preferring Turi–Borgis (TB) model predict interior-state and surface-state cluster isomers. The computed associated physical properties are also analyzed and compared to available experimental data. We find that the LGS and m-LGS potentials provide results that appear to be inconsistent with the size dependence of the experimental data. The simulated TB tendencies are qualitatively correct. Furthermore, ab initio calculations on static LGS noncavity structures indicate weak stabilization of the excess electron in regions where the LGS potential preferably and strongly binds the electron. TB calculations give stabilization energies that are in line with the ab initio results. In conclusion, we observe that the cavity-preferring pseudopotential model predicts cluster physical properties in better agreement with experimental data and ab initio calculations than the models predicting noncavity structures for the hydrated electron.

On the nature of the solvated electron in ice Ih

The water-solvated excess electron (EE) is a key chemical agent whose hallmark signature, its asymmetric optical absorption spectrum, continues to be a topic of debate.

Potential Paths for the Hydrogen-Bond Relaxing With (H2O)N Cluster Size

Relaxation of the inter- and intra-molecular interactions for the hydrogen bond (O:H-O) between undercoordinated molecules determines the unusual behavior of water nanodroplets and nanobubbles. However, probing such potentials remains unreality. Here we show that the Lagrangian solution [Huang et al., J. Phys. Chem. B, 2013. 117: 13639] transforms the observed H-O bond (x = H) and O:H nonbond (x = L) lengths and phonon frequencies (dx, x) [Sun et al., J. Phys. Chem. Lett., 2013. 4: 2565] into the respective force constants and bond energies (kx, Ex) and hence enables the mapping of the potential paths for the O:H-O bond relaxing with water cluster size. Results show that molecular undercoordination not only reduces the molecular size (dH) with enhanced H-O energy from the bulk value of 3.97 to 5.10 eV for a H2O monomer, but also enlarges the molecular separation (dL) with reduced O:H energy from 95 to 35 meV for a dimer. The H-O energy gain raises the melting point from bulk value 273 to 310 K for the skin and the O:H energy loss lowers the freezing temperature from bulk value 258 to 202 K for 1.4 nm sized droplet, by dispersing the quasisolid phase boundaries.

Time- and angle-resolved photoemission spectroscopy of hydrated electrons near a liquid water surface

We present time- and angle-resolved photoemission spectroscopy of trapped electrons near liquid surfaces. Photoemission from the ground state of a hydrated electron at 260 nm is found to be isotropic, while anisotropic photoemission is observed for the excited states of 1,4-diazabicyclo[2,2,2]octane and I- in aqueous solutions. Our results indicate that surface and subsurface species create hydrated electrons in the bulk side. No signature of a surface-bound electron has been observed.

Relaxation mechanism of the hydrated electron

The relaxation dynamics of the photoexcited hydrated electron have been subject to conflicting interpretations. Here, we report time-resolved photoelectron spectra of hydrated electrons in a liquid microjet with the aim of clarifying ambiguities from previous experiments. A sequence of three ultrashort laser pulses (~100 femtosecond duration) successively created hydrated electrons by charge-transfer-to-solvent excitation of dissolved anions, electronically excited these electrons via the s→p transition, and then ejected them into vacuum. Two distinct transient signals were observed. One was assigned to the initially excited p-state with a lifetime of ~75 femtoseconds, and the other, with a lifetime of ~400 femtoseconds, was attributed to s-state electrons just after internal conversion in a nonequilibrated solvent environment. These assignments support the nonadiabatic relaxation model.

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.

Electron at the Surface of Water: Dehydrated or Not?

The hydrated electron is a crucial species in radiative processes, and it has been speculated that its behavior at the water surface could lead to specific interfacial chemical properties. Here, we address fundamental questions concerning the structure and energetics of an electron at the surface of water. We use the method of ab initio molecular dynamics, which was shown to provide a faithful description of solvated electrons in large water clusters and in bulk water. The present results clearly demonstrate that the surface electron is mostly buried in the interfacial water layer, with only about 10 % of its density protruding into the vapor phase. Consequently, it has a structure that is very similar to that of an electron solvated in the aqueous bulk. This points to a general feature of charges at the surface of water, namely, that they do not behave as half-dehydrated but rather as almost fully hydrated species.

Unraveling the Complex Nature of the Hydrated Electron

The structure of the hydrated electron, which is a key species in radiative processes in water, has remained elusive. The traditional cavity model has been questioned recently, but the newly suggested picture of an electron delocalized over a region of enhanced water density is controversial. Here, we present results from ab initio molecular dynamics simulations, where not only the excess electron but also the valence electrons of the surrounding water molecules are described quantum mechanically. Unlike in previous one-electron pseudopotential calculations, many-electron interactions are explicitly accounted for. The present approach allows for partitioning of the electron solvated in liquid water into contributions from an inner cavity, neighboring water molecules, and a diffuse tail. We demonstrate that all three of these contributions are sizable and, consequently, important, which underlines the complex nature of the hydrated electron and warns against oversimplified interpretations based on pseudopotential models.

Solvation and evolution dynamics of an excess electron in supercritical CO2

We present an ab initio molecular dynamics simulation of the dynamics of an excess electron solvated in supercritical CO2. The excess electron can exist in three types of states: CO2-core localized, dual-core localized, and diffuse states. All these states undergo continuous state conversions via a combination of long lasting breathing oscillations and core switching, as also characterized by highly cooperative oscillations of the excess electron volume and vertical detachment energy. All of these oscillations exhibit a strong correlation with the electron-impacted bending vibration of the core CO2, and the core-switching is controlled by thermal fluctuations.

Solvated electrons at the water-air interface: surface versus bulk signal in low kinetic energy photoelectron spectroscopy

Time-resolved photoelectron spectroscopy at low kinetic energies (≲5 eV) is applied to dilute iodide solutions with different surface and bulk contributions. The results indicate a pronounced surface sensitivity. Signals assigned to solvated electrons near the liquid surface decay rapidly on a sub-ps timescale. In contrast to the literature, a long-lived surface solvated electron at 1.6 eV binding energy is not observed.

Dynamics and reactivity of trapped electrons on supported ice crystallites

The solvation dynamics and reactivity of localized excess electrons in aqueous environments have attracted great attention in many areas of physics, chemistry, and biology. This manifold attraction results from the importance of water as a solvent in nature as well as from the key role of low-energy electrons in many chemical reactions. One prominent example is the electron-induced dissociation of chlorofluorocarbons (CFCs). Low-energy electrons are also critical in the radiation chemistry that occurs in nuclear reactors. Excess electrons in an aqueous environment are localized and stabilized by the local rearrangement of the surrounding water dipoles. Such solvated or hydrated electrons are known to play an important role in systems such as biochemical reactions and atmospheric chemistry. Despite numerous studies over many years, little is known about the microscopic details of these electron-induced chemical processes, and interest in the fundamental processes involved in the reactivity of trapped electrons continues. In this Account, we present a surface science study of the dynamics and reactivity of such localized low-energy electrons at D(2)O crystallites that are supported by a Ru(001) single crystal metal surface. This approach enables us to investigate the generation and relaxation dynamics as well as dissociative electron attachment (DEA) reaction of excess electrons under well-defined conditions. They are generated by photoexcitation in the metal template and transferred to trapping sites at the vacuum interface of crystalline D(2)O islands. In these traps, the electrons are effectively decoupled from the electronic states of the metal template, leading to extraordinarily long excited state lifetimes on the order of minutes. Using these long-lived, low-energy electrons, we study the DEA to CFCl(3) that is coadsorbed at very low concentrations (∼10(12) cm(-2)). Using rate equations and direct measurement of the change of surface dipole moment, we estimated the electron surface density for DEA, yielding cross sections that are orders of magnitude higher than the electron density measured in the gas phase.

Structure, dynamics, and reactivity of hydrated electrons by ab initio molecular dynamics

Understanding the properties of hydrated electrons, which were first observed using pulse radiolysis of water in 1962, is crucial because they are key species in many radiation chemistry processes. Although time-resolved spectroscopic studies and molecular simulations have shown that an electron in water (prepared, for example, by water photoionization) relaxes quickly to a localized, cavity-like structure ∼2.5 Å in radius, this picture has recently been questioned. In another experimental approach, negatively charged water clusters of increasing size were studied with photoelectron and IR spectroscopies. Although small water clusters can bind an excess electron, their character is very different from bulk hydrated species. As data on electron binding in liquid water have become directly accessible experimentally, the cluster-to-bulk extrapolations have become a topic of lively debate. Quantum electronic structure calculations addressing experimental measurables have, until recently, been largely limited to small clusters; extended systems were approached mainly with pseudopotential calculations combining a classical description of water with a quantum mechanical treatment of the excess electron. In this Account, we discuss our investigations of electrons solvated in water by means of ab initio molecular dynamics simulations. This approach, applied to a model system of a negatively charged cluster of 32 water molecules, allows us to characterize structural, dynamical, and reactive aspects of the hydrated electron using all of the system's valence electrons. We show that under ambient conditions, the electron localizes into a cavity close to the surface of the liquid cluster. This cavity is, however, more flexible and accessible to water molecules than an analogous area around negatively charged ions. The dynamical process of electron attachment to a neutral water cluster is strongly temperature dependent. Under ambient conditions, the electron relaxes in the liquid cluster and becomes indistinguishable from an equilibrated, solvated electron on a picosecond time scale. In contrast, for solid, cryogenic systems, the electron only partially localizes outside of the cluster, being trapped in a metastable, weakly bound "cushion-like" state. Strongly bound states under cryogenic conditions could only be prepared by cooling equilibrated, liquid, negatively charged clusters. These calculations allow us to rationalize how different isomers of electrons in cryogenic clusters can be observed experimentally. Our results also bring into question the direct extrapolation of properties of cryogenic, negatively charged water clusters to those of electrons in the bulk liquid. Ab initio molecular dynamics represents a unique computational tool for investigating the reactivity of the solvated electron in water. As a prototype, the electron-proton reaction was followed in the 32-water cluster. In accord with experiment, the molecular mechanism is a proton transfer process that is not diffusion limited, but rather controlled by a proton-induced deformation of the excess electron's solvent shell. We demonstrate the necessary ingredients of a successful density functional methodology for the hydrated electron that avoids potential pitfalls, such as self-interaction error, insufficient basis set, or lack of dispersion interactions. We also benchmark the density functional theory methods and outline the path to faithful ab initio simulations of dynamics and reactivity of electrons solvated in extended aqueous systems.

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.

Excess electron localization in solvated DNA bases

We present a first-principles molecular dynamics study of an excess electron in condensed phase models of solvated DNA bases. Calculations on increasingly large microsolvated clusters taken from liquid phase simulations show that adiabatic electron affinities increase systematically upon solvation, as for optimized gas-phase geometries. Dynamical simulations after vertical attachment indicate that the excess electron, which is initially found delocalized, localizes around the nucleobases within a 15 fs time scale. This transition requires small rearrangements in the geometry of the bases.

Unveiling Electron Promiscuity

Although the wave-like proclivity of electrons for delocalization is familiar to every student of chemistry, it seems that electrons may have less respect for atomic and molecular boundaries than one might have considered proper. The boundaries in question include those between H-bonded dimers and within hydrated clusters, as well as those of aqueous cavities, colloidal suspensions, and macroscopic air-water and oil-water interfaces. Unveiling the promiscuous behavior of electrons at such frontiers may both raise eyebrows and demand acknowledgment.