Fast Emulation of Fermionic Circuits with Matrix Product States.

We describe a matrix product state (MPS) extension for the Fermionic Quantum Emulator (FQE) software library. We discuss the theory behind symmetry-adapted MPSs for approximating many-body wave functions of spin-1/2 Fermions, and we present an open-source, MPS-enabled implementation of the FQE interface (MPS-FQE). The software uses the open-source pyblock3 and block2 libraries for most elementary tensor operations, and it can largely be used as a drop-in replacement for FQE that allows for more efficient but approximate emulation of larger Fermionic circuits.

Quantum harmonic free energies for biomolecules and nanomaterials.

Obtaining the free energy of large molecules from quantum mechanical energy functions is a long-standing challenge. We describe a method that allows us to estimate, at the quantum mechanical level, the harmonic contributions to the thermodynamics of molecular systems of large size, with modest cost. Using this approach, we compute the vibrational thermodynamics of a series of diamond nanocrystals, and show that the error per atom decreases with system size in the limit of large systems.

Fast Exchange with Gaussian Basis Set Using Robust Pseudospectral Method.

In this article, we present an algorithm to efficiently evaluate the exchange matrix in periodic systems when a Gaussian basis set with pseudopotentials is used. The usual algorithm for evaluating exchange matrix scales cubically with the system size because one has to perform () fast Fourier transform (FFT). Here, we introduce an algorithm that retains the cubic scaling but reduces the prefactor significantly by eliminating the need to do FFTs during each exchange build.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package.

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods.

Conservation laws in coupled cluster dynamics at finite temperature.

We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [A. F. White and G.

A coupled cluster framework for electrons and phonons.

We describe a coupled cluster framework for coupled systems of electrons and harmonic phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard-Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations, which generally perform well for weak to moderate coupling.

Finite-temperature coupled cluster: Efficient implementation and application to prototypical systems.

We discuss the theory and implementation of the finite temperature coupled cluster singles and doubles (FT-CCSD) method including the equations necessary for an efficient implementation of response properties. Numerical aspects of the method including the truncation of the orbital space and integration of the amplitude equations are tested on some simple systems, and we provide some guidelines for applying the method in practice. The method is then applied to the 1D Hubbard model, the uniform electron gas (UEG) at warm, dense conditions, and some simple materials.

Generalized single excitation configuration interaction: an investigation into the impact of the inclusion of non-orthogonality on the calculation of core-excited states.

In this paper, we investigate different non-orthogonal generalizations of the configuration interaction with single substitutions (CIS) method for the calculation of core-excited states. Fully non-orthogonal CIS (NOCIS) has been described previously for species with singlet and doublet ground states, and this paper reports the extension to molecules in their triplet ground state. In addition to NOCIS, we present a novel method, one-center NOCIS (1C-NOCIS), for open-shell molecules which is intermediate between NOCIS and the computationally less demanding static exchange approximation (STEX).

Time-Dependent Coupled Cluster Theory on the Keldysh Contour for Nonequilibrium Systems.

We leverage the Keldysh formalism to extend our implementation of finite temperature coupled cluster theory [ 5690-5700] to thermal systems that have been driven out of equilibrium. The resulting Keldysh coupled cluster theory is discussed in detail. We describe the implementation of the equations necessary to perform Keldysh coupled cluster singles and doubles calculations of finite-temperature dynamics, and we apply the method to some simple systems including a Hubbard model with a Peierls phase and an ab initio model of warm-dense silicon subject to an ultrafast XUV pulse.

Non-Orthogonal Configuration Interaction with Single Substitutions for Core-Excited States: An Extension to Doublet Radicals.

In this paper, we present an open-shell extension of the non-orthogonal configuration interaction singles (NOCIS) method for the calculation of core-excited states, intended for peak assignment in XAS spectra of doublet radicals. This extension requires the consideration of additional configurations due to the singly occupied open-shell orbital, and the addition of essential orbital relaxation effects is found to provide a significant improvement on standard CIS, while maintaining the desirable properties of spin purity, variationality, and size consistency. We apply this method to the calculation of core excitations for several open-shell molecules and demonstrate that it performs competitively with other available methods, despite a lack of dynamic correlation.

A Time-Dependent Formulation of Coupled-Cluster Theory for Many-Fermion Systems at Finite Temperature.

We present a time-dependent formulation of coupled cluster theory. This theory allows for direct computation of the free energy of quantum systems at finite temperature by imaginary time integration and is closely related to the thermal cluster cumulant theory of Mukherjee and co-workers [ Chem. Phys.

Non-orthogonal configuration interaction with single substitutions for the calculation of core-excited states.

In this paper, we present the non-orthogonal configuration interaction singles (NOCIS) method for calculating core-excited states of closed-shell molecules. NOCIS is a black-box variant of NOCI, which uses A different core-ionized determinants for a molecule with A atoms of a given element to form single substitutions. NOCIS is a variational, spin-pure, size-consistent ab initio method that dramatically improves on standard CIS by capturing essential orbital relaxation effects, in addition to essential configuration interaction.

Second order Møller-Plesset and coupled cluster singles and doubles methods with complex basis functions for resonances in electron-molecule scattering.

The method of complex basis functions is applied to molecular resonances at correlated levels of theory. Møller-Plesset perturbation theory at second order and equation-of-motion electron attachment coupled-cluster singles and doubles (EOM-EA-CCSD) methods based on a non-Hermitian self-consistent-field reference are used to compute accurate Siegert energies for shape resonances in small molecules including N, CO, CO, and CHO. Analytic continuation of complex 𝜃-trajectories is used to compute Siegert energies, and the 𝜃-trajectories of energy differences are found to yield more consistent results than those of total energies.

Simulating the absorption spectra of helium clusters (N = 70, 150, 231, 300) using a charge transfer correction to superposition of fragment single excitations.

Simulations of the n = 2 absorption spectra of He (N = 70, 150, 231, 300) clusters are reported, with nuclear configurations sampled by path integral molecular dynamics. The electronic structure is treated by a new approach, ALMO-CIS+CT, which is a formulation of configuration interaction singles (CIS) based on absolutely localized molecular orbitals (ALMOs). The method generalizes the previously reported ALMO-CIS model [K.

Stabilizing potentials in bound state analytic continuation methods for electronic resonances in polyatomic molecules.

The computation of Siegert energies by analytic continuation of bound state energies has recently been applied to shape resonances in polyatomic molecules by several authors. We critically evaluate a recently proposed analytic continuation method based on low order (type III) Padé approximants as well as an analytic continuation method based on high order (type II) Padé approximants. We compare three classes of stabilizing potentials: Coulomb potentials, Gaussian potentials, and attenuated Coulomb potentials.

Restricted and unrestricted non-Hermitian Hartree-Fock: Theory, practical considerations, and applications to metastable molecular anions.

This work describes the implementation and applications of non-Hermitian self-consistent field (NH-SCF) theory with complex basis functions for the ab initio computation of positions and widths of shape resonances in molecules. We utilize both the restricted open-shell and the previously unexplored spin-unrestricted variants to compute Siegert energies of several anionic shape resonances in small diatomic and polyatomic molecules including carbon tetrafluoride which has been the subject of several recent experimental studies. The computation of general molecular properties from a non-Hermitian wavefunction is discussed, and a density-based analysis is applied to the (2)B1 shape resonance in formaldehyde.

Complex basis functions revisited: implementation with applications to carbon tetrafluoride and aromatic N-containing heterocycles within the static-exchange approximation.

The method of complex basis functions for computing positions and widths of molecular resonances is revisited. An open-ended and efficient implementation is described. The basis set requirements of the complex basis are investigated within the computationally inexpensive static-exchange approximation, and the results of this investigation lead to a hierarchy of basis sets for complex basis function calculations on small molecules.