Combining Complex Conjugation, Time-Reversal, and Spin-Flip Symmetry Projection of Coupled Cluster Wave Functions.

Complex conjugation symmetry breaking and restoration generates two nonorthogonal configurations at the Hartree-Fock level that can capture static correlation naturally. In conjunction with broken spin-symmetry coupled cluster theory, the symmetry-projected wave function shows good agreement with full configuration interaction in beryllium hydride insertion, lithium fluoride dissociation, and symmetric stretching of tetrahedral H4. By adding spin flip projection, we can also recover time reversal symmetry in the same coupled cluster framework. We also show results including point group symmetry projection.

Correlated pair ansatz with a binary tree structure.

We develop an efficient algorithm to implement the recently introduced binary tree state (BTS) ansatz on a classical computer. BTS allows a simple approximation to permanents arising from the computationally intractable antisymmetric product of interacting geminals and respects size-consistency. We show how to compute BTS overlap and reduced density matrices efficiently. We also explore two routes for developing correlated BTS approaches: Jastrow coupled cluster on BTS and linear combinations of BT states. The resulting methods show great promise in benchmark applications to the reduced Bardeen-Cooper-Schrieffer Hamiltonian and the one-dimensional XXZ Heisenberg Hamiltonian.

Hartree-Fock-Bogoliubov theory for number-parity-violating fermionic Hamiltonians.

It is usually asserted that physical Hamiltonians for fermions must contain an even number of fermion operators. This is indeed true in electronic structure theory. However, when the Jordan-Wigner (JW) transformation is used to map physical spin Hamiltonians to Hamiltonians of spinless fermions, terms that contain an odd number of fermion operators may appear. The resulting fermionic Hamiltonian thus does not have number parity symmetry and requires wave functions that do not have this symmetry either. In this work, we discuss the extension of standard Hartree-Fock-Bogoliubov (HFB) theory to the number-parity-nonconserving case. These ideas had appeared in the literature before but, perhaps for lack of practical applications, had, to the best of our knowledge, never been employed. We here present a useful application for this more general HFB theory based on coherent states of the SO(2M + 1) Lie group, where M is the number of orbitals. We also show how using these unusual mean-field states can provide significant improvements when studying the JW transformation of chemically relevant spin Hamiltonians.

Thermofield Theory for Finite-Temperature Electronic Structure.

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. They also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

State Preparation of Antisymmetrized Geminal Power on a Quantum Computer without Number Projection.

The antisymmetrized geminal power (AGP) is equivalent to the number projected Bardeen-Cooper-Schrieffer (PBCS) wave function. It is also an elementary symmetric polynomial (ESP) state. We generalize previous research on deterministically implementing the Dicke state to a state preparation algorithm for an ESP state, or equivalently AGP, on a quantum computer. Our method is deterministic and has polynomial cost, and it does not rely on number symmetry breaking and restoration. We also show that our circuit is equivalent to a disentangled unitary paired coupled cluster operator and a layer of unitary Jastrow operator acting on a single Slater determinant. The method presented herein highlights the ability of disentangled unitary coupled cluster to capture nontrivial entanglement properties that are hardly accessible with traditional Hartree-Fock based electronic structure methods.

Robust formulation of Wick's theorem for computing matrix elements between Hartree-Fock-Bogoliubov wavefunctions.

Numerical difficulties associated with computing matrix elements of operators between Hartree-Fock-Bogoliubov (HFB) wavefunctions have plagued the development of HFB-based many-body theories for decades. The problem arises from divisions by zero in the standard formulation of the nonorthogonal Wick's theorem in the limit of vanishing HFB overlap. In this Communication, we present a robust formulation of Wick's theorem that stays well-behaved regardless of whether the HFB states are orthogonal or not. This new formulation ensures cancellation between the zeros of the overlap and the poles of the Pfaffian, which appears naturally in fermionic systems. Our formula explicitly eliminates self-interaction, which otherwise causes additional numerical challenges. A computationally efficient version of our formalism enables robust symmetry-projected HFB calculations with the same computational cost as mean-field theories. Moreover, we avoid potentially diverging normalization factors by introducing a robust normalization procedure. The resulting formalism treats even and odd number of particles on equal footing and reduces to Hartree-Fock as a natural limit. As proof of concept, we present a numerically stable and accurate solution to a Jordan-Wigner-transformed Hamiltonian, whose singularities motivated the present work. Our robust formulation of Wick's theorem is a most promising development for methods using quasiparticle vacuum states.

Symmetry-projected cluster mean-field theory applied to spin systems.

We introduce Sz spin-projection based on cluster mean-field theory and apply it to the ground state of strongly correlated spin systems. In cluster mean-fields, the ground state wavefunction is written as a factorized tensor product of optimized cluster states. In previous work, we have focused on unrestricted cluster mean-field, where each cluster is Sz symmetry adapted. We here remove this restriction by introducing a generalized cluster mean-field (GcMF) theory, where each cluster is allowed to access all Sz sectors, breaking Sz symmetry. In addition, a projection scheme is used to restore global Sz, which gives rise to the Sz spin-projected generalized cluster mean-field (SzGcMF). Both of these extensions contribute to accounting for inter-cluster correlations. We benchmark these methods on the 1D, quasi-2D, and 2D J1 - J2 and XXZ Heisenberg models. Our results indicate that the new methods (GcMF and SzGcMF) provide a qualitative and semi-quantitative description of the Heisenberg lattices in the regimes considered, suggesting them as useful references for further inter-cluster correlations, which are discussed in this work.

A power series approximation in symmetry projected coupled cluster theory.

Projected Hartree-Fock theory provides an accurate description of many kinds of strong correlations but does not properly describe weakly correlated systems. On the other hand, single-reference methods, such as configuration interaction or coupled cluster theory, can handle weakly correlated problems but cannot properly account for strong correlations. Ideally, we would like to combine these techniques in a symmetry-projected coupled cluster approach, but this is far from straightforward. In this work, we provide an alternative formulation to identify the so-called disentangled cluster operators, which arise when we combine these two methodological strands. Our formulation shows promising results for model systems and small molecules.

Strong-weak duality via Jordan-Wigner transformation: Using fermionic methods for strongly correlated su(2) spin systems.

The Jordan-Wigner transformation establishes a duality between su(2) and fermionic algebras. We present qualitative arguments and numerical evidence that when mapping spins to fermions, the transformation makes strong correlation weaker, as demonstrated by the Hartree-Fock approximation to the transformed Hamiltonian. This result can be rationalized in terms of rank reduction of spin shift terms when transformed to fermions. Conversely, the mapping of fermions to qubits makes strong correlation stronger, complicating its solution when one uses qubit-based correlators. The presence of string operators poses challenges to the implementation of quantum chemistry methods on classical computers, but these can be dealt with using established techniques of low computational cost. Our proof of principle results for XXZ and J1-J2 Heisenberg (in 1D and 2D) indicates that the JW transformed fermionic Hamiltonian has reduced complexity in key regions of their phase diagrams and provides a better starting point for addressing challenging spin problems.

Construction of linearly independent non-orthogonal AGP states.

We show how to construct a linearly independent set of antisymmetrized geminal power (AGP) states, which allows us to rewrite our recently introduced geminal replacement models as linear combinations of non-orthogonal AGPs. This greatly simplifies the evaluation of matrix elements and permits us to introduce an AGP-based selective configuration interaction method, which can reach arbitrary excitation levels relative to a reference AGP, balancing accuracy and cost as we see fit.

Exploring non-linear correlators on AGP.

Single-reference methods such as Hartree-Fock-based coupled cluster theory are well known for their accuracy and efficiency for weakly correlated systems. For strongly correlated systems, more sophisticated methods are needed. Recent studies have revealed the potential of the antisymmetrized geminal power (AGP) as an excellent initial reference for the strong correlation problem. While these studies improved on AGP by linear correlators, we explore some non-linear exponential Ansätze in this paper. We investigate two approaches in particular. Similar to Wahlen-Strothman et al. [Phys. Rev. B 91, 041114(R) (2015)], we show that the similarity transformed Hamiltonian with a Hilbert-space Jastrow operator is summable to all orders and can be solved over AGP by projecting the Schrödinger equation. The second approach is based on approximating the unitary pair-hopper Ansatz recently proposed for application on a quantum computer. We report benchmark numerical calculations against the ground state of the pairing Hamiltonian for both of these approaches.

Correlating the antisymmetrized geminal power wave function.

Strong pairing correlations are responsible for superconductivity and off-diagonal long-range order in the two-particle density matrix. The antisymmetrized geminal power wave function was championed many years ago as the simplest model that can provide a reasonable qualitative description for these correlations without breaking number symmetry. The fact remains, however, that the antisymmetrized geminal power is not generally quantitatively accurate in all correlation regimes. In this work, we discuss how we might use this wave function as a reference state for a more sophisticated correlation technique such as configuration interaction, coupled cluster theory, or the random phase approximation.

Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems.

We present a wave function representation for the canonical ensemble thermal density matrix by projecting the thermofield double state against the desired number of particles. The resulting canonical thermal state obeys an imaginary-time evolution equation. Starting with the mean-field approximation, where the canonical thermal state becomes an antisymmetrized geminal power (AGP) wave function, we explore two different schemes to add correlation: by number-projecting a correlated grand-canonical thermal state and by adding correlation to the number-projected mean-field state. As benchmark examples, we use number-projected configuration interaction and an AGP-based perturbation theory to study the hydrogen molecule in a minimal basis and the six-site Hubbard model.

Understanding band gaps of solids in generalized Kohn-Sham theory.

The fundamental energy gap of a periodic solid distinguishes insulators from metals and characterizes low-energy single-electron excitations. However, the gap in the band structure of the exact multiplicative Kohn-Sham (KS) potential substantially underestimates the fundamental gap, a major limitation of KS density-functional theory. Here, we give a simple proof of a theorem: In generalized KS theory (GKS), the band gap of an extended system equals the fundamental gap for the approximate functional if the GKS potential operator is continuous and the density change is delocalized when an electron or hole is added. Our theorem explains how GKS band gaps from metageneralized gradient approximations (meta-GGAs) and hybrid functionals can be more realistic than those from GGAs or even from the exact KS potential. The theorem also follows from earlier work. The band edges in the GKS one-electron spectrum are also related to measurable energies. A linear chain of hydrogen molecules, solid aluminum arsenide, and solid argon provide numerical illustrations.

Thermofield theory for finite-temperature quantum chemistry.

Thermofield dynamics has proven to be a very useful theory in high-energy physics, particularly since it permits the treatment of both time- and temperature-dependence on an equal footing. We here show that it also has an excellent potential for studying thermal properties of electronic systems in physics and chemistry. We describe a general framework for constructing finite temperature correlated wave function methods typical of ground state methods. We then introduce two distinct approaches to the resulting imaginary time Schrödinger equation, which we refer to as fixed-reference and covariant methods. As an example, we derive the two corresponding versions of thermal configuration interaction theory and apply them to the Hubbard model, while comparing with exact benchmark results.

Exact parameterization of fermionic wave functions via unitary coupled cluster theory.

A formal analysis is conducted on the exactness of various forms of unitary coupled cluster (UCC) theory based on particle-hole excitation and de-excitation operators. Both the conventional single exponential UCC parameterization and a factorized (referred to here as "disentangled") version are considered. We formulate a differential cluster analysis to determine the UCC amplitudes corresponding to a general quantum state. The exactness of conventional UCC (ability to represent any state) is explored numerically, and it is formally shown to be determined by the structure of the critical points of the UCC exponential mapping. A family of disentangled UCC wave functions is proven to exactly parameterize any state, thus showing how to construct Trotter-error-free parameterizations of UCC for applications in quantum computing. From these results, we construct an exact disentangled UCC parameterization that employs an infinite sequence of particle-hole or general one- and two-body substitution operators.

Polynomial-product states: A symmetry-projection-based factorization of the full coupled cluster wavefunction in terms of polynomials of double excitations.

Our goal is to remedy the failure of symmetry-adapted coupled-cluster theory in the presence of strong correlation. Previous work along these lines has taken us from a diagram-level analysis of the coupled-cluster equations to an understanding of the collective modes which can occur in various channels of the coupled-cluster equations to the exploration of non-exponential wavefunctions in efforts to combine coupled-cluster theory with symmetry projection. In this manuscript, we extend these efforts by introducing a new, polynomial product wavefunction ansatz that incorporates information from symmetry projection into standard coupled-cluster theory in a way that attempts to mitigate the effects of the lack of size extensivity and size consistency characteristic of symmetry-projected methods. We describe the new approach in detail within the context of our previous efforts, explore some illustrative calculations, and consider one route for reducing the computational cost of the new method.

Efficient evaluation of AGP reduced density matrices.

We propose and implement an algorithm to calculate the norm and reduced density matrices (RDMs) of the antisymmetrized geminal power of any rank with polynomial cost. Our method scales quadratically per element of the RDMs. Numerical tests indicate that our method is very fast and capable of treating systems with a few thousand orbitals and hundreds of electrons reliably in double-precision. In addition, we present reconstruction formulas that allow one to decompose higher order RDMs in terms of linear combinations of lower order ones and geminal coefficients, thereby reducing the computational cost significantly.