skip to main content

Title: Ab Initio Simulations of Poorly- and Well-Equilibrated (CH3CN)n¯ Cluster Anions: Assigning Experimental Photoelectron Peaks to Surface-Bound Electrons and Solvated Monomer and Dimer Anions
Excess electrons in liquid acetonitrile are of particular interest because they exist in two different forms in equilibrium: they can be present as traditional solvated electrons in a cavity, and they can form some type of solvated molecular anion. Studies of small acetonitrile cluster anions in the gas phase show two isomers with distinct vertical detachment energies, and it is tempting to presume that the two gas-phase cluster anion isomers are precursors of the two excess electron species present in bulk solution. In this paper, we perform DFT-based ab initio molecular dynamics simulations of acetonitrile cluster anions to understand the electronic species that are present and why they have different binding energies. Using a long-range-corrected density functional that was optimally tuned to describe acetonitrile cluster anion structures, we have theoretically explored the chemistry of (CH3CN)n¯ cluster anions with sizes n=5,7 and 10. Since the temperature of the experimental cluster anions is not known, we performed two sets of simulations that investigated how the way in which the cluster anions are prepared affects the excess electron binding motif: one set of simulations simply attached excess electrons to neutral (CH3CN)n clusters, providing little opportunity for the clusters to relax in the presence more » of the excess electron, while the other set allowed the cluster anions to thermally equilibrate near room temperature. We find that both sets of simulations show three distinct electron binding motifs: electrons can attach to the surface of the cluster (dipole-bound) or be present as either solvated monomer anions, CH3CN¯, or as solvated molecular dimer anions, (CH3CN)2¯. All three species have higher binding energies at larger cluster sizes. Thermal equilibration strongly favors the formation of the valence-bound molecular anions relative to surface-bound excess electrons, and the dimer anion becomes more stable than the monomer anion and surface-bound species as the cluster size increases. The calculated photoelectron spectra from our simulations in which there was poor thermal equilibration are in good agreement with experiment, suggesting assignment of the two experimental cluster anion isomers as the surface-bound electron and the solvated molecular dimer anion. The simulations also suggest that the shoulder seen experimentally on the low-energy isomer's detachment peak is not part of a vibronic progression but instead results from molecular monomer anions. Nowhere in the size range that we explore do we see evidence for a non-valence, cavity-bound interior-solvated electron, indicating that this species is likely only accessible at larger sizes with good thermal equilibration. « less
Award ID(s):
Publication Date:
Journal Name:
The journal of physical chemistry
Sponsoring Org:
National Science Foundation
More Like this
  1. Quadrupole-bound anions are negative ions in which their excess electrons are loosely bound by long-range electron-quadrupole attractions. Experimental evidence for quadrupole-bound anions has been scarce; until now, only trans -succinonitrile had been experimentally confirmed to form a quadrupole-bound anion. In this study, we present experimental evidence for a new quadrupole-bound anion. Our combined Rydberg electron transfer/anion photoelectron spectroscopy study demonstrates that the ee conformer of 1,4-dicyanocyclohexane (DCCH) supports a quadrupole-bound anion state, and that the cis -DCCH conformer forms a dipole-bound anion state. The electron binding energies of the quadrupole- and dipole-bound anions are measured as 18 and 115 meV,more »respectively, both of which are in excellent agreement with theoretical calculations by Sommerfeld.« less
  2. Despite its importance in electron transfer reactions and radiation chemistry, there has been disagreement over the fundamental nature of the hydrated electron, such as whether or not it resides in a cavity. Mixed quantum/classical simulations of the hydrated electron give different structures depending on the pseudopotential employed, and ab initio models of computational necessity use small numbers of water molecules and/or provide insufficient statistics to compare to experimental observables. A few years ago, Kumar et al. (J. Phys. Chem. A 2015, 119, 9148) proposed a minimalist ab initio model of the hydrated electron with only a small number of explicitlymore »treated water molecules plus a polarizable continuum model (PCM). They found that the optimized geometry had four waters arranged tetrahedrally around a central cavity, and that the calculated vertical detachment energy and radius of gyration agreed well with experiment, results that were largely independent of the level of theory employed. The model, however, is based on a fixed structure at 0 K and does not explicitly incorporate entropic contributions or the thermal fluctuations that should be associated with the room-temperature hydrated electron. Thus, in this paper, we extend the model of Kumar et al. by running Born−Oppenheimer molecular dynamics (BOMD) of a small number of water molecules with an excess electron plus PCM at room temperature. We find that when thermal fluctuations are introduced, the level of theory chosen becomes critical enough when only four waters are used that one of the waters dissociates from the cluster with certain density functionals. Moreover, even with an optimally tuned range-separated hybrid functional, at room temperature the tetrahedral orientation of the 0 K first-shell waters is entirely lost and the central cavity collapses, a process driven by the fact that the explicit water molecules prefer to make H-bonds with each other more than with the excess electron. The resulting average structure is quite similar to that produced by a noncavity mixed quantum/classical model, so that the minimalist 4-water BOMD models suffer from problems similar to those of noncavity models, such as predicting the wrong sign of the hydrated electron’s molar solvation volume. We also performed BOMD with 16 explicit water molecules plus an extra electron and PCM. We find that the inclusion of an entire second solvation shell of explicit water leads to little change in the outcome from when only four waters were used. In fact, the 16-water simulations behave much like those of water cluster anions, in which the electron localizes at the cluster surface, showing that PCM is not acceptable for use in minimalist models to describe the behavior of the bulk hydrated electron. For both the 4- and 16-water models, we investigate how the introduction of thermal motions alters the predicted absorption spectrum, vertical detachment energy, and resonance Raman spectrum of the simulated hydrated electron. We also present a set of structural criteria that can be used to numerically determine how cavity-like (or not) a particular hydrated electron model is. All of the results emphasize that the hydrated electron is a statistical object whose properties are inadequately captured using only a small number of explicit waters, and that a proper treatment of thermal fluctuations is critical to understanding the hydrated electron’s chemical and physical behavior.« less
  3. We have prepared the hydrogen sulfide trimer and tetramer anions, (H 2 S) 3 − and (H 2 S) 4 − , measured their anion photoelectron spectra, and applied high-level quantum chemical calculations to interpret the results. The sharp peaks at low electron binding energies in their photoelectron spectra and their diffuse Dyson orbitals are evidence for them both being dipole-bound anions. While the dipole moments of the neutral (H 2 S) 3 and (H 2 S) 4 clusters are small, the excess electron induces structural distortions that enhance the charge-dipolar attraction and facilitate the binding of diffuse electrons.
  4. Few anions exhibit electronically excited states, and, if they do, the one or two possible excitations typically transpire beyond the visible spectrum into the near-infrared. These few, red-shifted electronic absorption features make anions tantalizing candidates as carriers of the diffuse interstellar bands (DIBs), a series of mostly unknown, astronomically ubiquitous absorption features documented for over a century. The recent interstellar detection of benzonitrile implies that cyano-functionalized polycyclic aromatic hydrocarbon (PAH) anions may be present in space. The presently reported quantum chemical work explores the electronic properties of deprotonated benzene, naphthalene, and anthracene anions functionalized with a single cyano group. Bothmore »the absorption and emission properties of the electronically excited states are explored. The findings show that the larger anions absorption and emission energies possess both valence and dipole bound excitations in the 450–900 nm range with oscillator strengths for both types of >1×10−4. The valence and dipole bound excited state transitions will produce slightly altered substructure from one another making them appear to originate with different molecules. The known interstellar presence of related molecules, the two differing natures of the excited states for each, and the wavelength range of peaks for these cyano-functionalized PAH anions are coincident with DIB properties. Finally, the methods utilized appear to be able to predict the presence of dipole-bound excited states to within a 1.0 meV window relative to the electron binding energy.« less
  5. Anionic states of benzonitrile are investigated by high-level electronic structure methods. The calculations using equation-of-motion coupled-cluster theory for electron-attached states confirm earlier conclusions drawn from the photodetachment experiments wherein the ground state of the anion is the valence 2 B 1 state, while the dipole bound state lies adiabatically ∼0.1 eV above. Inclusion of triple excitations and zero-point vibrational energies is important for recovering relative state correct ordering. The computed Franck–Condon factors and photodetachment cross-sections further confirm that the observed photodetachment spectrum originates from the valence anion. The valence anion is electronically bound at its equilibrium geometry, but it ismore »metastable at the equilibrium geometry of the neutral. The dipole-bound state, which is the only bound anionic state at the neutral equilibrium geometry, may serve as a gateway state for capturing the electron. Thus, the emerging mechanistic picture entails electron capture via a dipole bound state, followed by non-adiabatic relaxation forming valence anions.« less