Ross E. Larsen

Does the Hydrated Electron Occupy a Cavity

Filling a Cavity Unlike liquid ammonia, water cannot sustain a steady concentration of isolated electrons. Nonetheless, high-energy irradiation can introduce a small number of free charges that engage in potent reductive chemistry and have clear spectroscopic signatures. The manner in which water solubilizes these hydrated electrons has remained uncertain, but the general consensus has been that repulsive interactions drive the nearest water molecules away, leaving the electron in a nearly spherical empty cavity. Larsen et al. (p. 65; see the Perspective by Jordan and Johnson) upend this consensus with simulations based on a more thorough potential function for modeling the competing attractions and repulsions between the electron and surrounding water. The calculations suggest that the hydrated electron actually draws water in, occupying a region denser than the pure bulk liquid. The model reproduces experimental spectral and dynamic observations as effectively as, and in some cases better than, the cavity framework. A long-standing model of the solvent geometry surrounding a free charge in water is questioned by new numerical simulations. Since the discovery of the hydrated electron more than 40 years ago, a general consensus has emerged that the hydrated electron occupies a quasispherical cavity in liquid water. We simulated the electronic structure and dynamics of the hydrated electron using a rigorously derived pseudopotential to treat the electron-water interaction, which incorporates attractive oxygen and repulsive hydrogen features that have not been included in previous pseudopotentials. What emerged was a hydrated electron that did not reside in a cavity but instead occupied a ~1-nanometer-diameter region of enhanced water density. Both the calculated ground-state absorption spectrum and the excited-state spectral dynamics after simulated photoexcitation of this noncavity hydrated electron showed excellent agreement with experiment. The relaxation pathway involves a rapid internal conversion followed by slow ground-state cooling, the opposite of the mechanism implicated by simulations in which the hydrated electron occupies a cavity.

Instantaneous Pair Theory for High-Frequency Vibrational Energy Relaxation in Fluids

Notwithstanding the long and distinguished history of studies of vibrational energy relaxation, exactly how it is that high frequency vibrations manage to relax in a liquid remains somewhat of a mystery. Both experimental and theoretical approaches seem to say that there is a natural frequency range associated with intermolecular motion in liquids, typically spanning no more than a few hundred cm−1. Landau–Teller-type theories explain rather easily how a solvent can absorb any vibrational energy within this “band,” but how is it that molecules can rid themselves of superfluous vibrational energies significantly in excess of these values? In this paper we develop a theory for such processes based on the idea that the crucial liquid motions are those that most rapidly modulate the force on the vibrating coordinate — and that by far the most important of these motions are those involving what we have called the mutual nearest neighbors of the vibrating solute. Specifically, we suggest that whenever there is ...

A computationally efficient exact pseudopotential method. I. Analytic reformulation of the Phillips-Kleinman theory

Even with modern computers, it is still not possible to solve the Schrodinger equation exactly for systems with more than a handful of electrons. For many systems, the deeply bound core electrons serve merely as placeholders and only a few valence electrons participate in the chemical process of interest. Pseudopotential theory takes advantage of this fact to reduce the dimensionality of a multielectron chemical problem: the Schrodinger equation is solved only for the valence electrons, and the effects of the core electrons are included implicitly via an extra term in the Hamiltonian known as the pseudopotential. Phillips and Kleinman (PK) [Phys. Rev. 116, 287 (1959)]. demonstrated that it is possible to derive a pseudopotential that guarantees that the valence electron wave function is orthogonal to the (implicitly included) core electron wave functions. The PK theory, however, is expensive to implement since the pseudopotential is nonlocal and its computation involves iterative evaluation of the full Hamiltonian. In this paper, we present an analytically exact reformulation of the PK pseudopotential theory. Our reformulation has the advantage that it greatly simplifies the expressions that need to be evaluated during the iterative determination of the pseudopotential, greatly increasing the computational efficiency. We demonstrate our new formalism by calculating the pseudopotential for the 3s valence electron of the Na atom, and in the subsequent paper, we show that pseudopotentials for molecules as complex as tetrahydrofuran can be calculated with our formalism in only a few seconds. Our reformulation also provides a clear geometric interpretation of how the constraint equations in the PK theory, which are required to obtain a unique solution, are themselves sufficient to calculate the pseudopotential.

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

Resonance Raman and temperature-dependent electronic absorption spectra of cavity and noncavity models of the hydrated electron

Most of what is known about the structure of the hydrated electron comes from mixed quantum/classical simulations, which depend on the pseudopotential that couples the quantum electron to the classical water molecules. These potentials usually are highly repulsive, producing cavity-bound hydrated electrons that break the local water H-bonding structure. However, we recently developed a more attractive potential, which produces a hydrated electron that encompasses a region of enhanced water density. Both our noncavity and the various cavity models predict similar experimental observables. In this paper, we work to distinguish between these models by studying both the temperature dependence of the optical absorption spectrum, which provides insight into the balance of the attractive and repulsive terms in the potential, and the resonance Raman spectrum, which provides a direct measure of the local H-bonding environment near the electron. We find that only our noncavity model can capture the experimental red shift of the hydrated electron’s absorption spectrum with increasing temperature at constant density. Cavity models of the hydrated electron predict a solvation structure similar to that of the larger aqueous halides, leading to a Raman O–H stretching band that is blue-shifted and narrower than that of bulk water. In contrast, experiments show the hydrated electron has a broader and red-shifted O–H stretching band compared with bulk water, a feature recovered by our noncavity model. We conclude that although our noncavity model does not provide perfect quantitative agreement with experiment, the hydrated electron must have a significant degree of noncavity character.

Resolution and polarization in apertureless near-field microscopy

We present numerical calculations of the electromagnetic fields produced around a sharp metal tip held above a dielectric surface and illuminated from above by visible light. This geometry is that of an apertureless near-field scanning optical microscope (A-NSOM). The electric fields produced by irradiation are calculated using a real-space Green’s function technique, also known as the frequency-domain method of moments. We investigate numerically the convergence properties of this method with respect to the grid size used to discretize the equations and find that the usual criterion invoked to ensure convergence does not apply inside small metallic scatterers; for this reason, the fields near the border of the tip keep changing as we make the grid smaller and smaller. Although the calculations do not show convergence (for the intensity of the fields near the borders of the probe), the intensity and the polarization of the field in other regions of space are converged; they do not change as the grid spaci...

The role of electronic symmetry in charge-transfer-to-solvent reactions: Quantum nonadiabatic computer simulation of photoexcited sodium anions

Since charge-transfer-to-solvent (CTTS) reactions represent the simplest class of solvent-driven electron transfer reactions, there has been considerable interest in understanding the solvent motions responsible for electron ejection. The major question that we explore in this paper is what role the symmetry of the electronic states plays in determining the solvent motions that account for CTTS. To this end, we have performed a series of one-electron mixed quantum/classical nonadiabatic molecular dynamics simulations of the CTTS dynamics of sodide, Na−, which has its ground-state electron in an s orbital and solvent-supported CTTS excited states of p-like symmetry. We compare our simulations to previous theoretical work on the CTTS dynamics of the aqueous halides, in which the ground state has the electron in a p orbital and the CTTS excited state has s-like symmetry. We find that the key motions for Na− relaxation involve translations of solvent molecules into the node of the p-like CTTS excited state. T...

Close Packing of Nitroxide Radicals in Stable Organic Radical Polymeric Materials

The relationship between the polymer network and electronic transport properties for stable radical polymeric materials has come under investigation owing to their potential application in electronic devices. For the radical polymer poly(2,2,6,6-tetramethylpiperidine-4-yl-1-oxyl methacrylate), it is unclear whether the radical packing is optimal for charge transport partially because the relationship between radical packing and molecular structure is not well-understood. Using the paramagnetic nitroxide radical as a probe of the polymer and synthetic techniques to control the radical concentration on the methyl methacrylate backbone, we investigate the dependence of radical concentration on molecular structure. The electron paramagnetic resonance data indicate that radicals in the PTMA assume a closest approach distance to each other when more than 60% of the backbone is populated with radical pendant groups. Below 60% coverage, the polymer rearranges to accommodate larger radical-radical spacing. These findings are consistent with theoretical calculations and help explain some experimentally determined electron-transport properties.

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.

Quenching of the Perylene Fluorophore by Stable Nitroxide Radical-Containing Macromolecules

Stable nitroxide radical bearing organic polymer materials are attracting much attention for their application as next generation energy storage materials. A greater understanding of the inherent charge transfer mechanisms in such systems will ultimately be paramount to further advancements in the understanding of both intrafilm and interfacial ion- and electron-transfer reactions. This work is focused on advancing the fundamental understanding of these dynamic charge transfer properties by exploiting the fact that these species are efficient fluorescence quenchers. We systematically incorporated fluorescent perylene dyes into solutions containing the 2,2,6,6-tetramethylpiperidine-N-oxyl (TEMPO) radical and controlled their interaction by binding the TEMPO moiety into macromolecules with varying morphologies (e.g., chain length, density of radical pendant groups). In the case of the model compound, 4-oxo-TEMPO, quenching of the perylene excited state was found to be dominated by a dynamic (collisional) process, with a contribution from an apparent static process that is described by an ∼2 nm quenching sphere of action. When we incorporated the TEMPO unit into a macromolecule, the quenching behavior was altered significantly. The results can be described by using two models: (A) a collisional quenching process that becomes less efficient, presumably due to a reduction in the diffusion constant of the quenching entity, with a quenching sphere of action similar to 4-oxo-TEMPO or (B) a collisional quenching process that becomes more efficient as the radius of interaction grows larger with increasing oligomer length. This is the first study that definitively illustrates that fluorophore quenching by a polymer system cannot be explained using merely a classical Stern-Volmer approach but rather necessitates a more complex model.