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

Capturing Strong Correlation in Molecules with Phaseless Auxiliary-Field Quantum Monte Carlo Using Generalized Hartree-Fock Trial Wave Functions

Generalized Hartree-Fock (GHF) is a long-established electronic structure method that can lower the energy (compared to spin-restricted variants) by breaking physical wave function symmetries, namely S ^ 2 and S ^ z . After an exposition of GHF theory, we assess the use of GHF trial wave functions in phaseless auxiliary field quantum Monte Carlo (ph-AFQMC-G) calculations of strongly correlated molecular systems including symmetrically stretched hydrogen rings, carbon dioxide, and dioxygen. Imaginary time propagation is able to restore S ^ 2 symmetry and yields energies of comparable or better accuracy than CCSD(T) with unrestricted HF and GHF references, and consistently smooth dissociation curves─a remarkable result given the relative scalability of ph-AFQMC-G to larger system sizes. The present exploration of model strongly correlated systems marks a promising starting point for future studies of more chemically relevant molecules, and demonstrates that ph-AFQMC-G provides a highly accurate (and, in contrast to active-space-based trials, relatively black box and always size-consistent) description of challenging systems exhibiting, e.g., antiferromagnetic coupling and/or geometric spin frustration.

Capturing Correlation Effects in Positron Binding to Atoms and Molecules

A major challenge in contemporary electronic structure theory involves the development of methods to describe in a balanced manner the contribution of correlation effects to energy differences. This challenge can be even greater for multicomponent systems containing more than one type of quantum particle. In the present work, we describe a flexible code for carrying out self-consistent field and configuration interaction (CI) calculations on multicomponent systems and use it to generate trial wave functions for use in diffusion Monte Carlo (DMC) calculations of the positron affinity of Be, Be2, Be4, Mg, CS2, and benzene. The resulting positron affinities (PAs) are in good agreement with the best values from the literature.

Demonstrating the Connection between the Nonvalence Correlation-Bound Anions of Polyaromatic Hydrocarbons and the Image Potential States of Graphene Using a One-Electron Model Hamiltonian

The ground and excited state nonvalence correlation-bound (NVCB) anion states of the C 6 n 2 H 6 n hexagonal polycylic aromatic hydrocarbons and of hexagonal C 6 n 2 graphene nanoflakes are characterized using a one-electron model Hamiltonian which incorporates atomic electrostatic moments up to the quadrupole, coupled inducible charges and dipoles, and atom-centered Gaussians to describe the short-range repulsive interactions. Extrapolation of the calculated electron binding energies of the lowest energy symmetric and antisymmetric (with respect to the molecular plane) NVCB anions of both the polycylic aromatic hydrocarbons and the carbon nanoflakes to the n → ∞ limit yields binding energies that are in good agreement with those of the most stable symmetric and antisymmetric image potential states of freestanding graphene as determined from two-photon photoemission spectroscopy (2PPE) experiments.

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.

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.