Benchmark Calculations of the Energies for Binding Excess Electrons to Water Clusters

Single-Pole Polarization Models: Rapid Evaluation of Electron Affinities of Solvated-Electron and Superatomic Molecular Anionic States

We propose a single-parameter effective one-particle potential, termed the single-pole exchange-correlation (1p-XC), to rapidly evaluate electron affinities (EAs) of nonvalence electronic states of molecular clusters and nanoassemblies. The model combines exact-exchange and the random phase approximation (RPA) correlation potential with a single-pole approximation to model the frequency-dependent polarization function. It captures long-range static and dynamic-frequency effects in the correlation potential, with mean absolute errors of 0.06 eV for EAs of hydrated- and ammoniated-electron clusters with EA values in the range 0.24-1.77 eV. The 1p-XC approximation enables EA estimation with a computational wall-time similar to that of hybrid functionals. The model also provides a compressed-basis, which significantly reduces the rank of higher-level parameter-free one-particle Hamiltonians and further simplifies the computation of EAs. The compressed-basis approach is used to model the hybridization of superatomic molecular states of (C60)2- and (C60)3-, thereby verifying previous model Hamiltonian studies.

On the performance of second-order approximate coupled-cluster singles and doubles methods for non-valence anions

We investigate the capability of several variants of the second-order approximate coupled-cluster singles and doubles (CC2) method to describe dipole-bound, quadrupole-bound, and correlation-bound molecular anions. The binding energy of anions formed by electron attachment to closed-shell molecules is computed using the electron attachment variant of CC2 (EA-CC2), whereas anions with a closed-shell ground state are treated with the standard CC2 method that preserves the number of particles. We find that EA-CC2 captures the binding energies of dipole-bound radical anions quite well, whereas results for other types of non-valence anions are less reliable. We also test the performance of semi-empirical spin-scaling factors for all types of non-valence anions and observe that the spin-scaled CC2 variants generally do not provide more accurate binding energies for dipole-bound anions, while the binding energies of quadrupole-bound and correlation-bound anions are improved. As exemplary applications of EA-CC2, we investigate the dipole-bound anions of the steroids cortisol, progesterone, and testosterone. In addition, we characterize electron attachment to sym-tetracyanonaphthalene, a molecule that supports five anionic states, two of which can be interpreted as hitherto unobserved π-type quadrupole-bound states.

Non-Valence Anions of Pyridine and the Diazines

The dipole-bound anions of pyridine, pyridazine, and pyrimidine are characterized using equation of motion coupled cluster singles and doubles calculations. These calculations predict that the anions of pyridine, pyrimidine and pyridazine are bound in the Born-Oppenheimer approximation by 0.05, 0.8, and 19.0 meV, respectively. The binding energies of pyrimidine and pyridazine are large enough that the anions will remain bound even when allowing for corrections to the Born-Oppenheimer approximation, while that of pyridine is a borderline case. We were unable to find a stable non-valence correlation-bound anion for pyrazine, which has a zero dipole moment.

Phenanthrene: establishing lower and upper bounds to the binding energy of a very weakly bound anion

Quite a few molecules do not form stable anions that survive the time needed for their detection; their electron affinities (EA) are either very small or negative. How does one measure the EA if the anion cannot be observed? Or, at least, can one establish lower and upper bounds to their EA? We propose two approaches that provide lower and upper bounds. We choose the phenanthrene (Ph) molecule whose EA is controversial. Through competition between helium evaporation and electron detachment in HenPh- clusters, formed in helium nanodroplets, we estimate the lower bound of the vertical detachment energy (VDE) of Ph- as about -3 meV. In the second step, Ph is complexed with calcium whose electron affinity is just 24.55 meV. When CaPh- ions are collided with a thermal gas of argon, one observes Ca- product ions but no Ph-, suggesting that the EA of Ph is below that of Ca.

Dynamics of the Bulk Hydrated Electron from Many-Body Wave-Function Theory

The structure of the hydrated electron is a matter of debate as it evades direct experimental observation owing to the short life time and low concentrations of the species. Herein, the first molecular dynamics simulation of the bulk hydrated electron based on correlated wave‐function theory provides conclusive evidence in favor of a persistent tetrahedral cavity made up by four water molecules, and against the existence of stable non‐cavity structures. Such a cavity is formed within less than a picosecond after the addition of an excess electron to neat liquid water, with less regular cavities appearing as intermediates. The cavities are bound together by weak H−H bonds, the number of which correlates well with the number of coordinated water molecules, each type of cavity leaving a distinct spectroscopic signature. Simulations predict regions of negative spin density and a gyration radius that are both in agreement with experimental data.

Ultrafast Intermolecular Energy Transfer from Vibrations to Electronic Motion

It is discussed how vibrationally excited molecules in their electronic ground state can transfer their vibrational energy to the electronic motion of neighbors and ionize them. Based on explicit examples of vibrationally excited molecules and anionic neighbors, it is demonstrated that the transfer can be extremely efficient at intermolecular distances much beyond distances at which the molecule and its neighbor can form a bond.

Observation and ultrafast dynamics of a nonvalence correlation-bound state of an anion

Nonvalence states of molecular anions play key roles in processes, such as electron mobility, in rare-gas liquids, radiation-induced damage to DNA, and the formation of anions in the interstellar medium. Recently, a class of nonvalence bound anion state has been predicted by theory in which correlation forces are predominantly responsible for binding the excess electron. We present a direct spectroscopic observation of this nonvalence correlation-bound state (CBS) in the para-toluquinone trimer cluster anion. Time-resolved photoelectron velocity map imaging shows that photodetachment of the CBS produces a narrow and highly anisotropic photoelectron distribution, consistent with detachment from an s-like orbital. The CBS is bound by ~50 meV and decays by vibration-mediated autodetachment with a lifetime of 700 ± 100 fs. These states are likely to be common in large and/or polarizable anions and clusters and may act as doorway states in electron attachment processes.

Efficient floating diffuse functions for accurate characterization of the surface-bound excess electrons in water cluster anions

In this work, the effect of diffuse function types (atom-centered diffuse functions versus floating functions and s-type versus p-type diffuse functions) on the structures and properties of three representative water cluster anions featuring a surface-bound excess electron is studied and we find that an effective combination of such two kinds of diffuse functions can not only reduce the computational cost but also, most importantly, considerably improve the accuracy of results and even avoid incorrect predictions of spectra and the EE shape. Our results indicate that (a) simple augmentation of atom-centered diffuse functions is beneficial for the vertical detachment energy convergence, but it leads to very poor descriptions for the singly occupied molecular orbital (SOMO) and lowest unoccupied molecular orbital (LUMO) distributions of the water cluster anions featuring a surface-bound excess electron and thus a significant ultraviolet spectrum redshift; (b) the ghost-atom-based floating diffuse functions can not only contribute to accurate electronic calculations of the ground state but also avoid poor and even incorrect descriptions of the SOMO and the LUMO induced by excessive augmentation of atom-centered diffuse functions; (c) the floating functions can be realized by ghost atoms and their positions could be determined through an optimization routine along the dipole moment vector direction. In addition, both the s- and p-type floating functions are necessary to supplement in the basis set which are responsible for the ground (s-type character) and excited (p-type character) states of the surface-bound excess electron, respectively. The exponents of the diffuse functions should also be determined to make the diffuse functions cover the main region of the excess electron distribution. Note that excessive augmentation of such diffuse functions is redundant and even can lead to unreasonable LUMO characteristics.

Excess electrons in methanol clusters: Beyond the one-electron picture

We performed a series of comparative quantum chemical calculations on various size negatively charged methanol clusters, CH₃OH_n^-. The clusters are examined in their optimized geometries (n = 2-4), and in geometries taken from mixed quantum-classical molecular dynamics simulations at finite temperature (n = 2-128). These latter structures model potential electron binding sites in methanol clusters and in bulk methanol. In particular, we compute the vertical detachment energy (VDE) of an excess electron from increasing size methanol cluster anions using quantum chemical computations at various levels of theory including a one-electron pseudopotential model, several density functional theory (DFT) based methods, MP2 and coupled-cluster CCSD(T) calculations. The results suggest that at least four methanol molecules are needed to bind an excess electron on a hydrogen bonded methanol chain in a dipole bound state. Larger methanol clusters are able to form stronger interactions with an excess electron. The two simulated excess electron binding motifs in methanol clusters, interior and surface states, correlate well with distinct, experimentally found VDE tendencies with size. Interior states in a solvent cavity are stabilized significantly stronger than electron states on cluster surfaces. Although we find that all the examined quantum chemistry methods more or less overestimate the strength of the experimental excess electron stabilization, MP2, LC-BLYP, and BHandHLYP methods with diffuse basis sets provide a significantly better estimate of the VDE than traditional DFT methods (BLYP, B3LYP, X3LYP, PBE0). A comparison to the better performing many electron methods indicates that the examined one-electron pseudopotential can be reasonably used in simulations for systems of larger size.

Benchmark calculations of excess electrons in water cluster cavities: balancing the addition of atom-centered diffuse functions versus floating diffuse functions

Diffuse functions have been proved to be especially crucial for the accurate characterization of excess electrons which are usually bound weakly in intermolecular zones far away from the nuclei. To examine the effects of diffuse functions on the nature of the cavity-shaped excess electrons in water cluster surroundings, both the HOMO and LUMO distributions, vertical detachment energies (VDEs) and visible absorption spectra of two selected (H2O)24(-) isomers are investigated in the present work. Two main types of diffuse functions are considered in calculations including the Pople-style atom-centered diffuse functions and the ghost-atom-based floating diffuse functions. It is found that augmentation of atom-centered diffuse functions contributes to a better description of the HOMO (corresponding to the VDE convergence), in agreement with previous studies, but also leads to unreasonable diffuse characters of the LUMO with significant red-shifts in the visible spectra, which is against the conventional point of view that the more the diffuse functions, the better the results. The issue of designing extra floating functions for excess electrons has also been systematically discussed, which indicates that the floating diffuse functions are necessary not only for reducing the computational cost but also for improving both the HOMO and LUMO accuracy. Thus, the basis sets with a combination of partial atom-centered diffuse functions and floating diffuse functions are recommended for a reliable description of the weakly bound electrons. This work presents an efficient way for characterizing the electronic properties of weakly bound electrons accurately by balancing the addition of atom-centered diffuse functions and floating diffuse functions and also by balancing the computational cost and accuracy of the calculated results, and thus is very useful in the relevant calculations of various solvated electron systems and weakly bound anionic systems.

Nonvalence Correlation-Bound Anion States of Polycyclic Aromatic Hydrocarbons

In this work, we characterize the nonvalence correlation-bound anion states of several polycyclic aromatic hydrocarbon (PAH) molecules. Unlike the analogous image potential states of graphene that localize the charge density of the excess electron above and below the plane of the sheet, we find that for PAHs, much of the charge distribution of the excess electron is localized around the periphery of the molecule. This is a consequence of the electrostatic interaction of the electron with the polar CH groups. By replacing the H atoms by F atoms or the CH groups by N atoms, the charge density of the excess electron shifts from the periphery to above and below the plane of the ring systems.

Barrierless Single-Electron-Induced cis-trans Isomerization

Lowering the activation energy of a chemical reaction is an essential part in controlling chemical reactions. By attaching a single electron, a barrierless path for the cis-trans isomerization of maleonitrile on the anionic surface is formed. The anionic activation can be applied in both reaction directions, yielding the desired isomer. We identify the microscopic mechanism that leads to the formation of the barrierless route for the electron-induced isomerization. The generalization to other chemical reactions is discussed.

Infrared spectra of small anionic water clusters from density functional theory and wavefunction theory calculations

We performed systematic theoretical studies on small anionic water/deuterated water clusters W/D(-)(N=2-6) at both density functional theory (B3LYP) and wavefunction theory (MP2) levels. The focus of the study is to examine the convergence of calculated infrared (IR) spectra with respect to the increasing number of diffuse functions. It is found that at the MP2 level for larger clusters (n = 4-6), only one extra diffuse function is needed to obtain the converged relative IR intensities, while two or three more sets of extra diffuse functions are needed for smaller clusters. Such behaviour is strongly associated with the convergence of the electronic structure of corresponding clusters at the MP2 level. It is striking to observe that at the B3LYP level, the calculated relative IR intensities for all the clusters under investigations are diverse and show no trend of convergence upon increasing the number of diffuse functions. Moreover, the increasing contribution from the extra diffuse functions to the dynamic IR dipole moment indicates that the B3LYP electronic structure also fails to converge. These results manifest that MP2 is a preferential theoretical method, as compared to the widely used B3LYP, for the IR intensity of dipole bounded electron systems.

Quantifying solvated electrons' delocalization

The electron delocalization range EDR(r;u_av) (left) captures the spin density (right) of an electron delocalized over u_av = 5.77 Å on the surface of an (H₂O)₂₀⁻ cluster.

The HAlF4 superacid fragmentation induced by an excess electron attachment

The HAlF₄ superacid binds an excess electron and undergoes a spontaneous fragmentation that leads to a H atom and an AlF₄⁻ superhalogen anion.

Nonvalence correlation-bound anion states of spherical fullerenes

We present a one-electron model Hamiltonian for characterizing nonvalence correlation-bound anion states of fullerene molecules. These states are the finite system analogs of image potential states of metallic surfaces. The model potential accounts for both atomic and charge-flow polarization and is used to characterize the nonvalence correlation-bound anion states of the C60, (C60)2, C240, and C60@C240 fullerene systems. Although C60 is found to have a single (s-type) nonvalence correlation-bound anion state, the larger fullerenes are demonstrated to have multiple nonvalence correlation-bound anion states.

Excess electrons bound to molecular systems with a vanishing dipole but large molecular quadrupole

Electron attachment properties of covalent molecules and ion clusters with vanishing dipole moments but large quadrupoles are studied with coupled cluster ab initio methods. Selection of the molecules studied is driven by two goals, finding a paradigm quadrupole-bound anion and investigating whether there is a correlation between the magnitude of the molecular quadrupole and the vertical attachment energy. Out of all examined species, only the ion clusters and four of the covalent molecules are found to support bound anions. The shapes and spatial extents of the associated excess electron distributions are qualitatively and quantitatively characterized, respectively. Two of the four covalent systems are especially promising as paradigm systems because of advantageous trade-offs regarding the number of isomers and conformers as well as synthetic closeness to commercial sources. No correlation was found between the vertical attachment energy and molecular quadrupole in an analysis that included the newly identified bound anions, those molecules, which were found not to support bound anions, and succinonitrile, which had been studied before. Moreover, there is clearly no such thing as a "critical quadrupole moment". There are, however, very strong electron correlation effects involved in the binding of the excess electrons, and similar to succinonitrile, for five out of six anions identified here, the molecular quadrupole of the neutral itself is too weak to bind an excess electron, and electron correlation in the form of dynamic polarization is required to do so.

Existence of a Correlation Bound s-Type Anion State of C60

It is established using high-level electronic structure calculations that C60 has an s-type correlation-bound anion state with an electron binding energy of about 120 meV. Examination of the "singly occupied" natural orbital of the anion reveals that about 9% of the charge density of the excess electron is localized inside, and about 91% is localized outside the C60 cage. Calculations were also carried out for the He@C60, Ne@C60, and H2O@C60 endohedral complexes. For each of these species, the s-type anion is predicted to be less strongly bound than for C60 itself.