Effect of N Atom Substitution on Electronic Resonances: A 2D Photoelectron Spectroscopic and Computational Study of Anthracene, Acridine, and Phenazine Anions.

The accommodation of an excess electron by polycyclic aromatic hydrocarbons (PAHs) has important chemical and technological implications ranging from molecular electronics to charge balance in interstellar molecular clouds. Here, we use two-dimensional photoelectron spectroscopy and equation-of-motion coupled-cluster calculations of the radical anions of acridine (C13H9N-) and phenazine (C12H8N2-) and compare our results for these species to those for the anthracene anion (C14H10-). The calculations predict the observed resonances and additionally find low-energy two-particle-one-hole states, which are not immediately apparent in the spectra, and offer a slightly revised interpretation of the resonances in anthracene. The study of acridine and phenazine allows us to understand how N atom substitution affects electron accommodation. While the electron affinity associated with the ground state anion undergoes a sizable increase with the successive substitution of N atoms, the two lowest energy excited anion states are not affected significantly by the substitution. The net result is that there is an increase in the energy gap between the two lowest energy resonances and the bound ground electronic state of the radical anion from anthracene to acridine to phenazine. Based on an energy gap law for the rate of internal conversion, this increased gap makes ground state formation progressively less likely, as evidenced by the photoelectron spectra.

Use of bound state methods to calculate partial and total widths of shape resonances.

In this work we study the 2Π resonances of a two-site model system designed to mimic a smooth transition from the 2Πg temporary anion of N2 to the 2Π temporary anion of CO. The model system possesses the advantage that scattering and bound state (L2) methods can be directly compared without obfuscating electron-correlation effects. Specifically, we compare resonance parameters obtained with the complex Kohn variational (CKV) method with those from stabilization, complex absorbing potential, and regularized analytical continuation calculations. The CKV calculations provide p-wave and d-wave widths, the sum of which provides a good approximation of the total width. Then we demonstrate that the width obtained with modified bound state methods depends on the basis set employed: It can be the total width, a partial width, or an ill-defined sum of partial widths. Provided the basis set is chosen appropriately, widths from bound state methods agree well with the CKV results.

Mapping the temperature-dependent and network site-specific onset of spectral diffusion at the surface of a water cluster cage.

We explore the kinetic processes that sustain equilibrium in a microscopic, finite system. This is accomplished by monitoring the spontaneous, time-dependent frequency evolution (the frequency autocorrelation) of a single OH oscillator, embedded in a water cluster held in a temperature-controlled ion trap. The measurements are carried out by applying two-color, infrared-infrared photodissociation mass spectrometry to the D3O+·(HDO)(D2O)19 isotopologue of the "magic number" protonated water cluster, H+·(H2O)21 The OH group can occupy any one of the five spectroscopically distinct sites in the distorted pentagonal dodecahedron cage structure. The OH frequency is observed to evolve over tens of milliseconds in the temperature range (90 to 120 K). Starting at 100 K, large "jumps" are observed between two OH frequencies separated by ∼300 cm-1, indicating migration of the OH group from the bound OH site at 3,350 cm-1 to the free position at 3,686 cm-1 Increasing the temperature to 110 K leads to partial interconversion among many sites. All sites are observed to interconvert at 120 K such that the distribution of the unique OH group among them adopts the form one would expect for a canonical ensemble. The spectral dynamics displayed by the clusters thus offer an unprecedented view into the molecular-level processes that drive spectral diffusion in an extended network of water molecules.

Crystallographic Mapping and Tuning of Water Adsorption in Metal-Organic Frameworks Featuring Distinct Open Metal Sites.

Crucial steps toward designing water sorption materials and fine-tuning their properties for specific applications include precise identification of adsorption sites and establishment of rigorous molecular-level insight into the water adsorption process. We report stepwise crystallographic mapping and density functional theory computations of adsorbed water molecules in ALP-MOF-1, a metal-organic framework decorated with distinct open metal sites and carbonyl functional groups that serve as water anchoring sites for seeding the nucleation of a complex water network. Identification of an unusual water adsorption step in ALP-MOF-1 motivated the tuning of metal ion composition to adjust water uptake. These studies provide direct evidence that the identity of the open metal sites in MOFs can dramatically affect water adsorption behavior between 0 and ∼20% RH and that multiple proximal water anchoring sites along the MOF skeleton facilitate water uptake which could be potentially useful for applications requiring rapid and energetically facile water sorption.

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.

Vibrational Signatures of HNO3 Acidity When Complexed with Microhydrated Alkali Metal Ions, M+·(HNO3)(H2O)n=5 (M = Li, K, Na, Rb, Cs), at 20 K.

The speciation of strong acids like HNO3 under conditions of restricted hydration is an important factor in the rates of chemical reactions at the air-water interface. Here, we explore the trade-offs at play when HNO3 is attached to alkali ions (Li+-Cs+) with four water molecules in their primary hydration shells. This is achieved by analyzing the vibrational spectra of the M+·(HNO3)(H2O)5 clusters cooled to about 20 K in a cryogenic photofragmentation mass spectrometer. The local acidity of the acidic OH group is estimated by the extent of the red shift in its stretching frequency when attached to a single water molecule. The persistence of this local structural motif (HNO3-H2O) in all of these alkali metal clusters enables us to determine the competition between the effect of the direct complexation of the acid with the cation, which acts to enhance acidity, and the role of the water network in the first hydration shell around the ions, which acts to counter (screen) the intrinsic effect of the ion. Analysis of the vibrational features associated with the acid molecule, as well as those of the water network, reveals how cooperative interactions in the microhydration regime conspire to effectively offset the intrinsic enhancement of HNO3 acidity afforded by attachment to the smaller cations.

Water Network Shape-Dependence of Local Interactions with the Microhydrated -NO2- and -CO2- Anionic Head Groups by Cold Ion Vibrational Spectroscopy.

We report the structural evolutions of water networks and solvatochromic response of the CH3NO2- radical anion in the OH and CH stretching regions by analysis of the vibrational spectra displayed by cryogenically cooled CH3NO2-·(H2O)n=1-6 clusters. The OH stretching bands evolve with a surprisingly large discontinuity at n = 6, which features the emergence of an intense, strongly red-shifted band along with a weaker feature that appears in the region assigned to a free OH fundamental. Very similar behavior is displayed by the perdeuterated carboxylate clusters, RCO2-·(H2O)n=5-7 (R = CD3CD2), indicating that this behavior is a general feature in the microhydration of the triatomic anionic domain and not associated with CH oscillators. Electronic structure calculations trace this behavior to the formation of a "book" isomer of the water hexamer that adopts a configuration in which one of the water molecules resides in an acceptor-acceptor-donor (AAD) (A = acceptor, D = donor) H-bonding site. Excitation of the bound OH in the AAD site explores the local network topology best suited to stabilize an incipient -XO2H-OH-(H2O)2 intracluster proton-transfer reaction. These systems thus provide particularly clear examples where the network shape controls the potential energy landscape that governs water network-mediated, intracluster proton transfer. The CH stretching bands of the CH3NO2-·(H2O)n=1-6 clusters also exhibit strong solvatochromic shifts, but in this case, they smoothly blue-shift with increasing hydration with no discontinuity at n = 6. This behavior is analyzed in the context of the solute-ion polarizability response and partial charge transfer to the water networks.

The binding of atomic hydrogen on graphene from density functional theory and diffusion Monte Carlo calculations.

In this work, density functional theory (DFT) and diffusion Monte Carlo (DMC) methods are used to calculate the binding energy of a H atom chemisorbed on the graphene surface. The DMC value of the binding energy is about 16% smaller in magnitude than the Perdew-Burke-Ernzerhof (PBE) result. The inclusion of exact exchange through the use of the Heyd-Scuseria-Ernzerhof functional brings the DFT value of the binding energy closer in line with the DMC result. It is also found that there are significant differences in the charge distributions determined using PBE and DMC approaches.

Molecular-level origin of the carboxylate head group response to divalent metal ion complexation at the air-water interface.

We exploit gas-phase cluster ion techniques to provide insight into the local interactions underlying divalent metal ion-driven changes in the spectra of carboxylic acids at the air-water interface. This information clarifies the experimental findings that the CO stretching bands of long-chain acids appear at very similar energies when the head group is deprotonated by high subphase pH or exposed to relatively high concentrations of Ca2+ metal ions. To this end, we report the evolution of the vibrational spectra of size-selected [Ca2+·RCO2-]+·(H2O) n=0to12 and RCO2-·(H2O) n=0to14 cluster ions toward the features observed at the air-water interface. Surprisingly, not only does stepwise hydration of the RCO2- anion and the [Ca2+·RCO2-]+ contact ion pair yield solvatochromic responses in opposite directions, but in both cases, the responses of the 2 (symmetric and asymmetric stretching) CO bands to hydration are opposite to each other. The result is that both CO bands evolve toward their interfacial asymptotes from opposite directions. Simulations of the [Ca2+·RCO2-]+·(H2O) n clusters indicate that the metal ion remains directly bound to the head group in a contact ion pair motif as the asymmetric CO stretch converges at the interfacial value by n = 12. This establishes that direct metal complexation or deprotonation can account for the interfacial behavior. We discuss these effects in the context of a model that invokes the water network-dependent local electric field along the C-C bond that connects the head group to the hydrocarbon tail as the key microscopic parameter that is correlated with the observed trends.

Temporary Anion Resonances of Pyrene: A 2D Photoelectron Imaging and Computational Study.

The low-energy electron-scattering resonances of pyrene were characterized using experimental and computational methods. Experimentally, a two-dimensional photoelectron imaging of the pyrene anion was used to probe the dynamics of resonances over the first 4 eV of the continuum. Computationally, the energies and character of the anion states were determined using equation-of-motion coupled cluster calculations, while taking specific care to avoid the collapse onto discretized continuum levels, and an application of the pairing theorem. Our results are in good agreement with the predictions of electron-scattering calculations that include an offset and with the pyrene anion absorption spectrum in a glass matrix. Taken together, we offer an assignment of the first five electronic resonances of pyrene. Some of the population in the lowest-energy 2B1u resonance was observed to decay to the ground electronic state of the anion, while all other resonances decay by a direct autodetachment. The astronomical relevance of a ground-state electron capture proceeding via a low-energy resonance in pyrene is discussed.

Role of Overlap between the Discrete State and Pseudocontinuum States in Stabilization Calculations of Metastable States.

In a diabatic picture metastable states subject to decay by electron detachment can be viewed as arising from the coupling between a discrete state and a continuum. In treating such states with bound-state quantum chemical methods, the continuum is discretized. In this study, we elucidate the role of overlap in this interaction in the application of the stabilization method to temporary anion states. This is accomplished by use of a minimalist stabilization calculation on the lowest energy l=2 (D) resonance of the finite spherical well potential using two basis functions, one describing the diabatic discrete state and the other a diabatic discretized continuum state. We show that even such a simple treatment predicts a complex resonance energy in good agreement with the exact result. If the energy of the discrete state is assumed to be constant, which is tantamount to orthogonalizing the discretized continuum state to the discrete state, it is demonstrated that the square of the off-diagonal coupling has a maximum close to the crossing point of the orthogonalized diabatic curves and that the curvature in the coupling is responsible for the complex stationary point associated with the resonance. Moreover, this curvature is a consequence of the overlap between the two diabatic states.

Frontiers of stochastic electronic structure calculations.

In recent years there has been a rapid growth in the development and application of new stochastic methods in electronic structure. These methods are quite diverse, from many-body wave function techniques in real space or determinant space to being used to sum perturbative expansions. This growth has been spurred by the more favorable scaling with the number of electrons and often better parallelization over large numbers of central processing unit (CPU) cores or graphical processing units (GPUs) than for high-end non-stochastic wave function based methods. This special issue of the Journal of Chemical Physics includes 33 papers that describe recent developments and applications in this area. As seen from the articles in the issue, stochastic electronic structure methods are applicable to both molecules and solids and can accurately describe systems with strong electron correlation. This issue was motivated, in part, by the 2019 Telluride Science Research Center workshop on Stochastic Electronic Structure Methods that we organized. Below we briefly describe each of the papers in the special issue, dividing the papers into six subtopics.

Analysis of the Contributions to the Kinetic and Potential Energies of an H Atom in the Presence of a Point Charge: The Molecular Virial Theorem Revisited.

The molecular virial theorem states that for a diatomic molecule or for an atom in the presence of a point charge, the changes in the average kinetic energy and average potential energy are equal to [Formula: see text] and [Formula: see text], respectively, where U is the interaction energy and R is the internuclear separation or the atom-point charge separation. In this paper we directly evaluate the ⟨T⟩ and ⟨V⟩ expectation values of an H atom in the presence of a distant point charge, obtaining exact analytical expressions by use of Dalgarno-Lewis perturbation theory.

The role of high-order electron correlation effects in a model system for non-valence correlation-bound anions.

The diffusion Monte Carlo (DMC), auxiliary field quantum Monte Carlo (AFQMC), and equation-of-motion coupled cluster (EOM-CC) methods are used to calculate the electron binding energy (EBE) of the non-valence anion state of a model (H2O)4 cluster. Two geometries are considered, one at which the anion is unbound and the other at which it is bound in the Hartree-Fock (HF) approximation. It is demonstrated that DMC calculations can recover from the use of a HF trial wave function that has collapsed onto a discretized continuum solution, although larger EBEs are obtained when using a trial wave function for the anion that provides a more realistic description of the charge distribution and, hence, of the nodal surface. For the geometry at which the cluster has a non-valence correlation-bound anion, both the inclusion of triples in the EOM-CC method and the inclusion of supplemental diffuse d functions in the basis set are important. DMC calculations with suitable trial wave functions give EBE values in good agreement with our best estimate EOM-CC result. AFQMC using a trial wave function for the anion with a realistic electron density gives a value of the EBE nearly identical to the EOM-CC result when using the same basis set. For the geometry at which the anion is bound in the HF approximation, the inclusion of triple excitations in the EOM-CC calculations is much less important. The best estimate EOM-CC EBE value is in good agreement with the results of DMC calculations with appropriate trial wave functions.

Toward a systematic improvement of the fixed-node approximation in diffusion Monte Carlo for solids-A case study in diamond.

While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for ab initio electronic structure calculations, in practice, the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros). This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. This approach is illustrated in the case of the carbon diamond using Slater-Jastrow trial wavefunctions including up to one million Slater determinants. Fixed-node DMC energies obtained with such large expansions are much improved, and the fixed-node error is found to decrease monotonically and smoothly as a function of the number of determinants in the trial wavefunction, a property opening the way to a better control of this error. The cohesive energy extrapolated to the thermodynamic limit is in close agreement with the estimated experimental value. Interestingly, this is also the case at the single-determinant level, thus, indicating a very good error cancellation in carbon diamond between the bulk and atomic total fixed-node energies when using single-determinant nodes.

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo.

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

Model potential study of non-valence correlation-bound anions of (C60)n clusters: the role of electric field-induced charge transfer.

A polarization model which accounts for electric field-induced charge transfer between fullerene molecules is introduced. Application of this model to the C60 dimer and trimer shows that intermolecular charge transfer makes a significant contribution to the polarizabilities of these clusters. This polarization model is incorporated into a one-electron Hamiltonian for describing non-valence correlation-bound anions, allowing us to further demonstrate that intermolecular charge transfer also results in increased stability of these anion states.

One-Dimensional Adiabatic Model Approach for Calculating Progressions in Vibrational Spectra of Ion-Water Complexes.

Many water-anion complexes of the form X-·(H2O), where X- is a polyatomic anion, display a peak progression in the OH stretch region of the vibrational spectra with spacings of 65-85 cm-1. These progressions result from strong anharmonic coupling between the OH stretch and a low-frequency intermolecular rock vibration. In this study, we calculate these progressions in HCO2-·(H2O), NO3-·(H2O), and CS2-·(H2O) by use of a one-dimensional adiabatic model with rock potentials generated from ab initio energies and frequencies. The importance of using a geometry-dependent reduced mass in calculating the peak spacings is demonstrated. We find that the one-dimensional adiabatic model is more successful in predicting peak spacings in the spectrum of HCO2-·(H2O) than for NO3-·(H2O), for which the rock vibration is highly anharmonic.

Prediction of a Nonvalence Temporary Anion Shape Resonance for a Model (H2O)4 System.

Ab initio calculations are used to demonstrate the existence of a nonvalence temporary anion shape resonance for a model (H2O)4 cluster system with no net dipole moment. The cluster is composed of two water dimers, the distance between which is varied. Each dimer possesses a weakly bound nonvalence anion state. For large separations of the dimer subunits, there are two bound nonvalence anion states (of Ag and B2u symmetry) corresponding to the symmetric and asymmetric combinations of the nonvalence anion states of the two dimer subunits. As the separation between the dimer subunits is decreased, the B2u anion increases in energy and becomes a temporary anion shape resonance. The real part of the resonance energy is determined as a function of the distance between the dimers and is found to increase monotonically from just above threshold to 28 meV for the range of geometries considered. Over this same range of geometries, the resonance half-width varies from 0 to 21 meV. The B2u anion, both when bound and when temporary, has a very diffuse charge distribution. The effective radial potential for the interaction of the excess electron with the cluster has a barrier at large distance arising from the electron-quadrupole interaction in combination with the repulsive angular momentum ( l = 1) contribution. This barrier impacts both the resonance energy and its lifetime.