Excess Electrons Bound to Small Ammonia 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.

Microsolvation of electrons by a handful of ammonia molecules

Microsolvation of electrons in ammonia is studied here via anionic NH_{3 n} ^- clusters with n = 2-6. Intensive samplings of the corresponding configurational spaces using second-order perturbation theory with extended basis sets uncover rich and complex energy landscapes, heavily populated by many local minima in tight energy windows as calculated from highly correlated coupled cluster methods. There is a marked energetical preference for structures that place the excess electron external to the molecular frame, effectively coordinating it with the three protons from a single ammonia molecule. Overall, as the clusters grow in size, the lowest energy dimer serves as the basic motif over which additional ammonia molecules are attached via unusually strong charge-assisted hydrogen bonds. This is a priori quite unexpected because, on electrostatic grounds, the excess electron would be expected to be in contact with as many protons as possible. Accordingly, a full quantum mechanical treatment of the bonding interactions under the tools provided by the quantum theory of atoms in molecules is carried out in order to dissect and understand the nature of intermolecular contacts. Vertical detachment energies reveal bound electrons even for n = 2.

Dynamic role of the correlation effect revealed in the exceptionally slow autodetachment rates of the vibrational Feshbach resonances in the dipole-bound state

Real-time autodetachment dynamics of the loosely bound excess electron from the vibrational Feshbach resonances of the dipole-bound states (DBS) of 4-bromophonoxide (4-BrPhO-) and 4-chlorophenoxide (4-ClPhO-) anions have been thoroughly investigated. The state-specific autodetachment rate measurements obtained by the picosecond time-resolved pump-probe method on the cryogenically cooled anions exhibit an exceptionally long lifetime (τ) of ∼823 ± 156 ps for the 11'1 vibrational mode of the 4-BrPhO- DBS. Strong mode-dependency in the wide dynamic range has also been found, giving τ ∼ 5.3 ps for the 10'1 mode, for instance. Though it is nontrivial to get the state-specific rates for the 4-ClPhO- DBS, the average autodetachment lifetime of the 19'120'1/11'1 mode has been estimated to be ∼548 ± 108 ps. Observation of these exceptionally slow autodetachment rates of vibrational Feshbach resonances strongly indicates that the correlation effect may play a significant role in the DBS photodetachment dynamics. Fermi's golden rule has been invoked so that the correlation effect is taken into account in the form of the interaction between the charge and the induced dipole where the latter is given by the polarizable counterparts of the electron-rich halogenated compound and the diffuse non-valence electron. This report suggests that one may measure, from the real-time autodetachment dynamics, the extent of the correlation effect contribution to the stabilization and/or dynamics of the excess non-valence electron among many different types of long-range interactions of the DBS.

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 electrons bound to H2S trimer and tetramer clusters

We have prepared the hydrogen sulfide trimer and tetramer anions, (H2S)3− and (H2S)4−, measured their anion photoelectron spectra, and applied high-level quantum chemical calculations to interpret the results.

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.

High-resolution photoelectron imaging and resonant photoelectron spectroscopy via noncovalently bound excited states of cryogenically cooled anions

Valence-bound anions with polar neutral cores (μ > ∼2.5 D) can support dipole-bound excited states below the detachment threshold. These dipole-bound states (DBSs) are highly diffuse and the weakly bound electron in the DBS can be readily autodetached via vibronic coupling. Excited DBSs can be observed in photodetachment spectroscopy using a tunable laser. Tuning the detachment laser to above-threshold vibrational resonances yields vibrationally enhanced resonant photoelectron spectra, which are highly non-Franck-Condon with much richer vibrational information. This perspective describes recent advances in the studies of excited DBSs of cryogenically cooled anions using high-resolution photoelectron imaging (PEI) and resonant photoelectron spectroscopy (rPES). The basic features of dipole-bound excited states and highly non-Franck-Condon resonant photoelectron spectra will be discussed. The power of rPES to yield rich vibrational information beyond conventional PES will be highlighted, especially for low-frequency and Franck-Condon-inactive vibrational modes, which are otherwise not accessible from non-resonant conventional PES. Mode-selectivity and intra-molecular rescattering have been observed during the vibrationally induced autodetachment. Conformer-specific rPES is possible due to the different dipole-bound excited states of molecular conformers with polar neutral cores. For molecules with μ ≪ 2.5 D or without dipole moments, but large quadrupole moments, excited quadrupole-bound states can exist, which can also be used to conduct rPES.

Evidence of Electron Capture of an Outgoing Photoelectron Wave by a Nonvalence State in (C6F6) n. J Phys Chem Lett 2018;9:2504-2509. [PMID: 29694047 DOI: 10.1021/acs.jpclett.8b00739]

Frequency-resolved photoelectron spectra are presented for (C₆F₆) _n^- with n = 1-5 that show that C₆F₆^- is solvated by neutral C₆F₆ molecules. Direct photodetachment channels of C₆F₆^- are observed for all n, leaving the neutral in the S₀ ground state or triplet states, T₁ and T₂. For n ≥ 2, an additional indirect electron loss channel is observed when the triplet-state channels open. This indirect emission appears to arise from the electron capture of the outgoing photoelectron s-wave by a neutral solvent molecule through an anion nonvalence state. The same process is not observed for the S₀ detachment channel because the outgoing electron wave is predominantly a p-wave. Our results show that anion nonvalence states can act as electron-accepting states in cluster environments and can be viewed as precursor states for diffuse states of liquid C₆F₆.

Ultrafast dynamics of low-energy electron attachment via a non-valence correlation-bound state

The primary electron-attachment process in electron-driven chemistry represents one of the most fundamental chemical transformations with wide-ranging importance in science and technology. However, the mechanistic detail of the seemingly simple reaction of an electron and a neutral molecule to form an anion remains poorly understood, particularly at very low electron energies. Here, time-resolved photoelectron imaging was used to probe the electron-attachment process to a non-polar molecule using time-resolved methods. An initially populated diffuse non-valence state of the anion that is bound by correlation forces evolves coherently in ∼30 fs into a valence state of the anion. The extreme efficiency with which the correlation-bound state serves as a doorway state for low-energy electron attachment explains a number of electron-driven processes, such as anion formation in the interstellar medium and electron attachment to fullerenes.

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.

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.

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.

A first principles molecular dynamics study of excess electron and lithium atom solvation in water-ammonia mixed clusters: structural, spectral, and dynamical behaviors of [(H2O)5NH3]- and Li(H2O)5NH3 at finite temperature

First principles molecular dynamics simulations are carried out to investigate the solvation of an excess electron and a lithium atom in mixed water-ammonia cluster (H(2)O)(5)NH(3) at a finite temperature of 150 K. Both [(H(2)O)(5)NH(3)](-) and Li(H(2)O)(5)NH(3) clusters are seen to display substantial hydrogen bond dynamics due to thermal motion leading to many different isomeric structures. Also, the structures of these two clusters are found to be very different from each other and also very different from the corresponding neutral cluster without any excess electron or the metal atom. Spontaneous ionization of Li atom occurs in the case of Li(H(2)O)(5)NH(3). The spatial distribution of the singly occupied molecular orbital shows where and how the excess (or free) electron is primarily localized in these clusters. The populations of single acceptor (A), double acceptor (AA), and free (NIL) type water and ammonia molecules are found to be significantly high. The dangling hydrogens of these type of water or ammonia molecules are found to primarily capture the free electron. It is also found that the free electron binding motifs evolve with time due to thermal fluctuations and the vertical detachment energy of [(H(2)O)(5)NH(3)](-) and vertical ionization energy of Li(H(2)O)(5)NH(3) also change with time along the simulation trajectories. Assignments of the observed peaks in the vibrational power spectra are done and we found a one to one correlation between the time-averaged populations of water and ammonia molecules at different H-bonding sites with the various peaks of power spectra. The frequency-time correlation functions of OH stretch vibrational frequencies of these clusters are also calculated and their decay profiles are analyzed in terms of the dynamics of hydrogen bonded and dangling OH modes. It is found that the hydrogen bond lifetimes in these clusters are almost five to six times longer than that of pure liquid water at room temperature.