Search

Complex conjugation symmetry breaking and restoration generate two non-orthogonal 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 H$_4$. By adding spin flip projection, we can also recover time reversal symmetry in the same coupled cluster framework.

Addressing both dynamic and static correlation accurately is a primary goal in electronic structure theory. Non-orthogonal configuration interaction (NOCI) is a versatile tool for treating static correlation, offering chemical insights by combining diverse reference states. Nevertheless, achieving quantitative accuracy requires the inclusion of missing dynamic correlation. This work introduces a framework for compressing orthogonal single and double excitations into an NOCI of a much smaller dimension. This compression is repeated with each Slater determinant in a reference NOCI, resulting in another NOCI that includes all its single and double excitations (NOCISD), effectively recovering the missing dynamic correlations from the reference. This compressed NOCISD is further refined through a selection process using metric and energy tests (SNOCISD). We validate the effectiveness of SNOCISD through its application to the dissociation of the nitrogen molecule and the hole-doped two-dimensional Hubbard model at various interaction strengths.

We present an innovative cluster-based method employing linear combinations of diverse cluster mean-field (cMF) states, and apply it to describe the ground state of strongly-correlated spin systems. In cluster mean-field theory, the ground state wavefunction is expressed as a factorized tensor product of optimized cluster states. While our prior work concentrated on a single cMF tiling, this study removes that constraint by combining different tilings of cMF states. Selection criteria, including translational symmetry and spatial proximity, guide this process. We present benchmark calculations for the one- and two-dimensional $J_1-J_2$ and $XXZ$ Heisenberg models. Our findings highlight two key aspects. First, the method offers a semi-quantitative description of the $0.4 \lessapprox J_2/J_1 \lessapprox 0.6$ regime of the $J_1-J_2$ model - a particularly challenging regime for existing methods. Second, our results demonstrate the capability of our method to provide qualitative descriptions for all the models and regimes considered, establishing it as a valuable reference. However, the inclusion of additional (weak) correlations is necessary for quantitative agreement, and we explore methods to incorporate these extra correlations.

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 transformation is used to map physical spin Hamiltonians to Hamiltonians of spinless fermions, terms which contain an odd number of fermion operators may appear. The resulting fermionic Hamiltonian thus does not have number parity symmetry, and requires wave functions which 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(2$M$ + 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 Jordan-Wigner transformation of chemically relevant spin Hamiltonians., Comment: Resubmitted to J. Chem. Phys

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., Comment: 14 pages, 10 figures, revised

The Jordan-Wigner transformation is a powerful tool for converting systems of spins into systems of fermions, or vice versa. While this mapping is exact, the transformation itself depends on the labeling of the spins. One consequence of this dependence is that approximate solutions of a Jordan-Wigner--transformed Hamiltonian may depend on the (physically inconsequential) labeling of the spins. In this work, we turn to an extended Jordan-Wigner transformation which remedies this problem and which may also introduce some correlation atop the Hartree-Fock solution of a transformed spin Hamiltonian. We demonstrate that this extended Jordan-Wigner transformation can be thought of as arising from a unitary version of the Lie algebraic similarity transformation (LAST) theory. We show how these ideas, particularly in combination with the standard (non-unitary) version of LAST, can provide a potentially powerful tool for the treatment of the XXZ and $J_1-J_2$ Heisenberg Hamiltonians., Comment: The Version of Record of this manuscript has been published and is available in Molecular Physics on September 5, 2023, available at https://www.tandfonline.com/doi/full/10.1080/00268976.2023.2254857

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 paper, 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.

We introduce $S_z$ spin-projection based on cluster mean-field theory and apply it to the ground state of strongly-correlated spin systems. In cluster mean-field, 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 $S_z$ symmetry adapted. We here remove this restriction by introducing a generalized cluster mean-field (GcMF) theory, where each cluster is allowed to access all $S_z$ sectors, breaking $S_z$ symmetry. In addition, a projection scheme is used to restore global $S_z$, which gives rise to $S_z$ spin-projected generalized cluster mean-field (S$_z$GcMF). Both of these extensions contribute to accounting for inter-cluster correlations. We benchmark these methods on the 1D, quasi-2D, and 2D $J_1-J_2$ and $XXZ$ Heisenberg models. Our results indicate that the new methods (GcMF and S$_z$GcMF) 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.

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. Analogs 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 pre-factor 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., Comment: 10 pages, 5 figures, mini-review

The antisymmetrized geminal power (AGP), a wave function equivalent to number-projected Hartree--Fock--Bogoliubov (HFB), and number-projected Bardeen--Cooper--Schrieffer (BCS) when working in the paired (natural orbitals) basis, has proven to be an excellent reference for strong pairing interactions. Several correlation methods have also been applied on top of AGP. In this work, we show how AGP can also be applied to spin systems by simply basing its formulation on a spin $su(2)$ algebra. We here implement spin AGP and spin AGP-based correlation techniques and benchmark them on the XXZ and $\mathrm{J_1-J_2}$ Heisenberg models, both in 1 and 2 dimensions. Our results indicate that spin AGP is a promising starting point for modeling spin systems., Comment: 16 pages, 14 figures

The antisymmetrized geminal power (AGP) is equivalent to the number projected Bardeen-Cooper-Schrieffer (PBCS) wavefunction. 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 non-trivial entanglement properties that are hardly accessible with traditional Hartree-Fock based electronic structure methods., Comment: 11 pages, 6 figures, manuscript revised

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 J$_1$-J$_2$ Heisenberg (in 1D and 2D) indicate 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., Comment: First revision

Electronic structure methods typically benefit from symmetry breaking and restoration, specially in the strong correlation regime. The same goes for Ans\"atze on a quantum computer. We develop a unitary coupled cluster method on the antisymmetrized geminal power (AGP) -- a state formally equivalent to the number-projected Bardeen--Cooper--Schrieffer wavefunction. We demonstrate our method for the single-band Fermi--Hubbard Hamiltonian in one and two dimensions. We also explore post-selection as a state preparation step to obtain correlated AGP and prove that it scales no worse than $\mathcal{O}(\sqrt{M})$ in the number of measurements, thereby making it a less expensive alternative to gauge integration to restore particle number symmetry.

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 approaches 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.

Coupled cluster (CC) has established itself as a powerful theory to study correlated quantum many-body systems. Finite-temperature generalizations of CC theory have attracted considerable interest and have been shown to work as nicely as the ground-state theory. However, most of these recent developments address only fermionic or bosonic systems. The distinct structure of the $su(2)$ algebra requires the development of a similar thermal CC theory for spin degrees of freedom. In this paper, we provide a formulation of our thermofield-inspired thermal CC for SU(2) systems. We apply the thermal CC to the Lipkin-Meshkov-Glick system as well as the one-dimensional transverse field Ising model as benchmark applications to highlight the accuracy of thermal CC in the study of finite-temperature phase diagrams in SU(2) systems., Comment: 14 pages, 6 figures

We introduce perturbation and coupled-cluster theories based on a cluster mean-field reference for describing the ground state of strongly-correlated spin systems. In cluster mean-field, the ground state wavefunction is written as a simple tensor product of optimized cluster states. The cluster language and the mean-field nature of the ansatz allows for a straightforward improvement which uses perturbation theory and coupled-cluster to account for inter-cluster correlations. We present benchmark calculations on the 1D chain and 2D square $J_1-J_2$ Heisenberg model, using cluster mean-field, perturbation theory and coupled-cluster. We also present an extrapolation scheme that allows us to compute thermodynamic limit energies accurately. Our results indicate that, with sufficiently large clusters, the correlated methods (cPT2, cPT4 and cCCSD) can provide a relatively accurate description of the Heisenberg model in the regimes considered, which suggests that the methods presented can be used for other strongly-correlated systems. Some ways to improve upon the methods presented in this work are discussed.

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., Comment: 13 pages, 8 figures, 1 table, 3 appendices

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 ansatze in this paper. We investigate two approaches in particular. Similar to 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 Schrodinger's 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.

For variational algorithms on the near term quantum computing hardware, it is highly desirable to use very accurate ansatze with low implementation cost. Recent studies have shown that the antisymmetrized geminal power (AGP) wavefunction can be an excellent starting point for ansatze describing systems with strong pairing correlations, as those occurring in superconductors. In this work, we show how AGP can be efficiently implemented on a quantum computer with circuit depth, number of CNOTs, and number of measurements being linear in system size. Using AGP as the initial reference, we propose and implement a unitary correlator on AGP and benchmark it on the ground state of the pairing Hamiltonian. The results show highly accurate ground state energies in all correlation regimes of this model Hamiltonian.

The antisymmetrized geminal power (AGP) wavefunction has a long history and is known by different names in various chemical and physical problems. There has been recent interest in using AGP as a starting point for strongly correlated electrons. Here, we show that in a seniority-conserving regime, different AGP based correlator representations based on generators of the algebra, killing operators, and geminal replacement operators are all equivalent. We implement one representation that uses number operators as correlators and has linearly independent curvilinear metrics to distinguish the regions of Hilbert space. This correlation method called J-CI, provides excellent accuracy in energies when applied to the pairing Hamiltonian., Comment: 10 pages, 4 figures, 1 table

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 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., Comment: 8 pages, 3 figures, 1 supplemental information

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.

An intrinsic wavefunction with a broken continuous symmetry can be rotated with no energy penalty leading to an infinite set of degenerate states known as a Goldstone manifold. In this work, we show that a dual representation of such manifold exists that is sampled by an infinite set of non-degenerate states. A proof that both representations are equivalent is provided. From the work of Peierls and Yoccoz (Proc. Phys. Soc. A {\bf 70}, 381 (1957)), it is known that collective states with good symmetries can be obtained from the Goldstone manifold using a generator coordinate trial wavefunction. We show that an analogous generator coordinate can be used in the dual representation; we provide numerical evidence using an intrinsic wavefunction with particle number symmetry-breaking for the electronic structure of the Be atom and one with $\hat{S}^z$ symmetry-breaking for a H$_5$ ring. We discuss how the dual representation can be used to evaluate expectation values of symmetry-projected states when the norm $|\langle \Phi | \hat{P}^q | \Phi \rangle|$ becomes very small., Comment: 8 pages, 7 figures

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 disentangled (factorized) 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 are shown 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 derive an exact disentangled UCC parameterization that employs an infinite sequence of particle-hole or general one- and two-body substitution operators.

We present a coupled cluster and linear response theory to compute properties of many-electron systems at non-zero temperatures. For this purpose, we make use of the thermofield dynamics, which allows for a compact wavefunction representation of the thermal density matrix, and extend our recently developed framework [J. Chem. Phys. 150, 154109 (2019)] to parameterize the so-called thermal state using an exponential ansatz with cluster operators that create thermal quasiparticle excitations on a mean-field reference. As benchmark examples, we apply this method to both model (one-dimensional Hubbard and Pairing) as well as ab-initio (atomic Beryllium and molecular Hydrogen) systems, while comparing with exact results., Comment: 11 pages, 6 figures

The antisymmetrized geminal power (AGP) wave function has a long history and considerable conceptual appeal, but in many situations its accuracy is wanting. Here, we consider a form of configuration interaction (CI) based upon the AGP wave function and taking advantage of its killing operators to construct an excitation manifold. Our geminal CI reduces to standard single-determinant--based CI in the limit in which AGP reduces to a single determinant. It substantially improves upon AGP in the reduced BCS Hamiltonian, which serves as a prototype for the kinds of strong pairing correlations relevant in Bardeen-Cooper-Schrieffer--style superconductivity. Moreover, our geminal CI naturally generalizes to add correlation to more general geminal-based wave functions than AGP., Comment: Revised manuscript

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 Schrodinger 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.

The impact of applying state-of-the-art tensor factorization techniques to modern nuclear Hamiltonians derived from chiral effective field theory is investigated. Subsequently, the error induced by the tensor decomposition of the input Hamiltonian on ground-state energies of closed-shell nuclei calculated via second-order many-body perturbation theory is benchmarked. With the aid of the factorized Hamiltonian, the second-order perturbative correction to ground-state energies is decomposed and the scaling properties of the underlying tensor network are discussed. The employed tensor formats are found to lead to an efficient data compression of two-body matrix elements of the nuclear Hamiltonian. In particular, the sophisticated \emph{tensor hypercontraction} (THC) scheme yields low tensor ranks with respect to both harmonic-oscillator and Hartree-Fock single-particle bases. It is found that the tensor rank depends on the two-body total angular momentum $J$ for which one performs the decomposition, which is itself directly related to the sparsity the corresponding tensor. Furthermore, including normal-ordered two-body contributions originating from three-body interactions does not compromise the efficient data compression. Ultimately, the use of factorized matrix elements authorizes controlled approximations of the exact second-order ground-state energy corrections. In particular, a small enough error is obtained from low-rank factorizations in $^{4}$He, $^{16}$O and $^{40}$Ca., Comment: 16 pages, 13 figures, 1 table

Methods which aim at universal applicability must be able to describe both weak and strong electronic correlation with equal facility. Such methods are in short supply. The combination of symmetry projection for strong correlation and coupled cluster theory for weak correlation offers tantalizing promise to account for both on an equal footing. In order to do so, however, the coupled cluster portion of the wave function must be optimized in the presence of the symmetry projection. This paper discusses how this may be accomplished, and shows the importance of doing so for both the Hubbard model Hamiltonian and the molecular Hamiltonian, all with a computational scaling comparable to that of traditional coupled cluster theory., Comment: revised version

There have been assertions in the literature that the variational and unitary forms of coupled cluster theory lead to the same energy functional. Numerical evidence from previous authors was inconsistent with this claim, yet the small energy differences found between the two methods and the relatively large number of variational parameters precluded an unequivocal conclusion. Using the Lipkin Hamiltonian, we here present conclusive numerical evidence that the two theories yield different energies. The ambiguities arising from the size of the cluster parameter space are absent in the Lipkin model, particularly when truncating to double excitations. We show that in the symmetry adapted basis under strong correlation the differences between the variational and unitary models are large, whereas they yield quite similar energies in the weakly correlated regime previously explored. We also provide a qualitative argument rationalizing why these two models cannot be the same. Additionally, we study a generalized non-unitary and non-hermitian variant that contains excitation, de-excitation and mixed operators with different amplitudes and show that it works best when compared to the traditional, variational, unitary, and extended forms of coupled cluster doubles theories.

We study the behavior of Hartree-Fock (HF) solutions in the vicinity of conical intersections. These are here understood as regions of a molecular potential energy surface characterized by degenerate or nearly-degenerate eigenfunctions with identical quantum numbers (point group, spin, and electron number). Accidental degeneracies between states with different quantum numbers are known to induce symmetry breaking in HF. The most common closed-shell restricted HF instability is related to singlet-triplet spin degeneracies that lead to collinear unrestricted HF (UHF) solutions. Adding geometric frustration to the mix usually results in noncollinear generalized HF (GHF) solutions, identified by orbitals that are linear combinations of up and down spins. Near conical intersections, we observe the appearance of coplanar GHF solutions that break all symmetries, including complex conjugation and time-reversal, which do not carry good quantum numbers. We discuss several prototypical examples taken from the conical intersection literature. Additionally, we utilize a recently introduced a magnetization diagnostic to characterize these solutions, as well as a solution of a Jahn-Teller active geometry of H$_8^{+2}$., Comment: accepted to JCP December 2017, published online January 2018

The spin structure of wave functions is reflected in the magnetic structure of the one-particle density matrix. Indeed, for single determinants we can use either one to determine the other. In this work we discuss how one can simply examine the one-particle density matrix to faithfully determine whether the spin magnetization density vector field is collinear, coplanar, or noncoplanar. For single determinants, this test suffices to distinguish collinear determinants which are eigenfunctions of $\hat{S}_{\hat{n}}$ from noncollinear determinants which are not. We also point out the close relationship between noncoplanar magnetism on the one hand and complex conjugation symmetry breaking on the other. Finally, we use these ideas to classify the various ways single determinant wave functions break and respect symmetries of the Hamiltonian in terms of their one-particle density matrix., Comment: Submitted to J. Chem. Theory Comput

Doubly occupied configuration interaction (DOCI), the exact diagonalization of the Hamiltonian in the paired (seniority zero) sector of the Hilbert space, is a combinatorial cost wave function that can be very efficiently approximated by pair coupled cluster doubles (pCCD) at mean-field computational cost. As such, it is a very interesting candidate as a starting point for building the full configuration interaction (FCI) ground state eigenfunction belonging to all (not just paired) seniority sectors. The true seniority zero sector of FCI (referred to here as FCI${}_0$) includes the effect of coupling between all seniority sectors rather than just seniority zero, and is, in principle, different from DOCI. We here study the accuracy with which DOCI approximates FCI${}_0$. Using a set of small Hubbard lattices, where FCI is possible, we show that DOCI $\sim$ FCI${}_0$ under weak correlation. However, in the strong correlation regime, the nature of the FCI${}_0$ wavefunction can change significantly, rendering DOCI and pCCD a less than ideal starting point for approximating FCI.

We derive and implement a new way of solving coupled cluster equations with lower computational scaling. Our method is based on decomposition of both amplitudes and two electron integrals, using a combination of tensor hypercontraction and canonical polyadic decomposition. While the original theory scales as $O(N^6)$ with respect to the number of basis functions, we demonstrate numerically that we achieve sub-millihartree difference from the original theory with $O(N^4)$ scaling. This is accomplished by solving directly for the factors that decompose the cluster operator. The proposed scheme is quite general and can be easily extended to other many-body methods.

Coupled cluster theory is the method of choice for weakly correlated systems. But in the strongly correlated regime, it faces a symmetry dilemma, where it either completely fails to describe the system, or has to artificially break certain symmetries. On the other hand, projected Hartree-Fock theory captures the essential physics of many kinds of strong correlations via symmetry breaking and restoration. In this work, we combine and try to retain the merits of these two methods by applying symmetry projection to broken symmetry coupled cluster wavefunctions. The non-orthogonal nature of states resulting from the application of symmetry projection operators furnishes particle-hole excitations to all orders, thus creating an obstacle for the exact evaluation of overlaps. Here we provide a solution via a disentanglement framework theory that can be approximated rigorously and systematically. Results of projected coupled cluster theory are presented for molecules and the Hubbard model, showing that spin projection significantly improves unrestricted coupled cluster theory while restoring good quantum numbers. The energy of projected coupled cluster theory reduces to the unprojected one in the thermodynamic limit, albeit at a much slower rate than projected Hartree-Fock., Comment: Submitted to JCP. Extra figures appear in the ancillary file

The past several years have seen renewed interest in the use of symmetry-projected Hartree-Fock for the description of strong correlations. Unfortunately, these symmetry-projected mean-field methods do not adequately account for dynamic correlation. Presumably, this shortcoming could be addressed if one could combine symmetry-projected Hartree-Fock with a many-body method such as coupled cluster theory, but this is by no means straightforward because the two techniques are formulated in very different ways. However, we have recently shown that the singlet $S^2$-projected unrestricted Hartree-Fock wave function can in fact be written in a coupled cluster-like wave function: that is, the spin-projected unrestricted Hartree-Fock wave function can be written as a polynomial of a double-excitation operator acting on some closed-shell reference determinant. Here, we extend this result and show that the spin-projected generalized Hartree-Fock wave function (which has both $S^2$ and $S_z$ projection) is likewise a polynomial of low-order excitation operators acting on a closed-shell determinant, and provide a closed-form expression for the resulting polynomial coefficients. We include a few preliminary applications of the combination of this spin-projected Hartree-Fock and coupled cluster theory to the Hubbard Hamiltonian, and comment on generalizations of the methodology. Results here are not for production level, but a similarity transformed theory that combines the two offers the promise of being accurate for both weak and strong correlation, and particularly may offer significant improvements in the intermediate correlation regime where neither projected Hartree-Fock nor coupled cluster is particularly accurate., Comment: accepted by Phys. Rev. A

We present numerical results for the equation of state of an infinite chain of hydrogen atoms. A variety of modern many-body methods are employed, with exhaustive cross-checks and validation. Approaches for reaching the continuous space limit and the thermodynamic limit are investigated, proposed, and tested. The detailed comparisons provide a benchmark for assessing the current state of the art in many-body computation, and for the development of new methods. The ground-state energy per atom in the linear chain is accurately determined versus bondlength, with a confidence bound given on all uncertainties., Comment: 29 pages, 16 figures -- revised and updated version

The failures of single-reference coupled cluster for strongly correlated many-body systems is flagged at the mean-field level by the spontaneous breaking of one or more physical symmetries of the Hamiltonian. Restoring the symmetry of the mean-field determinant by projection reveals that coupled cluster fails because it factorizes high-order excitation amplitudes incorrectly. However, symmetry-projected mean-field wave functions do not account sufficiently for dynamic (or weak) correlation. Here we pursue a merger of symmetry projection and coupled cluster theory, following previous work along these lines that utilized the simple Lipkin model system as a testbed [J. Chem. Phys. 146, 054110 (2017)]. We generalize the concept of a symmetry-projected mean-field wave function to the concept of a symmetry projected state, in which the factorization of high-order excitation amplitudes in terms of low-order ones is guided by symmetry projection and is not exponential, and combine them with coupled cluster theory in order to model the ground state of the Agassi Hamiltonian. This model has two separate channels of correlation and two separate physical symmetries which are broken under strong correlation. We show how the combination of symmetry collective states and coupled cluster is effective in obtaining correlation energies and order parameters of the Agassi model throughout its phase diagram.

Projected Hartree-Fock theory provides an accurate description of many kinds of strong correlation but does not properly describe weakly correlated systems. Coupled cluster theory, in contrast, does the opposite. It therefore seems natural to combine the two so as to describe both strong and weak correlations with high accuracy in a relatively black-box manner. Combining the two approaches, however, is made more difficult by the fact that the two techniques are formulated very differently. In earlier work, we showed how to write spin-projected Hartree-Fock in a coupled-cluster-like language. Here, we fill in the gaps in that earlier work. Further, we combine projected Hartree-Fock and coupled cluster theory in a variational formulation and show how the combination performs for the description of the Hubbard Hamiltonian and for several small molecular systems., Comment: Published version

In electronic structure theory, restricted single-reference coupled cluster (CC) captures weak correlation but fails catastrophically under strong correlation. Spin-projected unrestricted Hartree-Fock (SUHF), on the other hand, misses weak correlation but captures a large portion of strong correlation. The theoretical description of many important processes, e.g. molecular dissociation, requires a method capable of accurately capturing both weak- and strong correlation simultaneously, and would likely benefit from a combined CC-SUHF approach. Based on what we have recently learned about SUHF written as particle-hole excitations out of a symmetry-adapted reference determinant, we here propose a heuristic coupled cluster doubles model to attenuate the dominant spin collective channel of the quadratic terms in the coupled cluster equations. Proof of principle results presented here are encouraging and point to several paths forward for improving the method further., Comment: 20 pages, 10 figures

Coupled cluster and symmetry projected Hartree-Fock are two central paradigms in electronic structure theory. However, they are very different. Single reference coupled cluster is highly successful for treating weakly correlated systems, but fails under strong correlation unless one sacrifices good quantum numbers and works with broken-symmetry wave functions, which is unphysical for finite systems. Symmetry projection is effective for the treatment of strong correlation at the mean-field level through multireference non-orthogonal configuration interaction wavefunctions, but unlike coupled cluster, it is neither size extensive nor ideal for treating dynamic correlation. We here examine different scenarios for merging these two dissimilar theories. We carry out this exercise over the integrable Lipkin model Hamiltonian, which despite its simplicity, encompasses non-trivial physics for degenerate systems and can be solved via diagonalization for a very large number of particles. We show how symmetry projection and coupled cluster doubles individually fail over different correlation limits, whereas models that merge these two theories are highly successful over the entire phase diagram. Despite the simplicity of the Lipkin Hamiltonian, the lessons learned in this work will be useful for building an ab initio symmetry projected coupled cluster theory that we expect to be accurate over the weakly and strongly correlated limits, as well as the recoupling regime.

We present a method incorporating biorthogonal orbital-optimization, symmetry projection, and double-occupancy screening with a non-unitary similarity transformation generated by the Gutzwiller factor $ n_{i\uparrow}n_{i\downarrow}$, and apply it to the Hubbard model. Energies are calculated with mean-field computational scaling with high-quality results comparable to coupled cluster singles and doubles. This builds on previous work performing similarity transformations with more general, two-body Jastrow-style correlators. The theory is tested on two-dimensional lattices ranging from small systems into the thermodynamic limit and is compared to available reference data.

Exchange-correlation hole is a central concept in density functional theory. It not only provides justification for an exchange-correlation energy functional, but also serves as a local ingredient in nonlocal range-separation density functional. However, due to the nonlocal nature, modelig the conventional exact exchange hole presents a great challenge to density functional theory. In this work, we propose a semilocal exchange hole underlying the Tao-Perdew-Staroverov-Scuseria (TPSS) meta-GGA functional. The present model is distinct from previous models at small separation between an electron and the hole around the electron. It is also different in the way it interpolates between the rapidly varying iso-orbital density and the slowly varying density, which is determined by the wave vector analysis based on the exactly solvable infinite barrier model for jellium surface. Our numerical tests show that the exchange hole generated from this model mimics the conventional exact exchange hole quite well for atoms. Finally, as a simple application, we apply the hole model to construct a TPSS-based range-separation functional. Our tests show that this TPSS-based range-separation functional can substantially improve TPSS band gaps and barrier heights, without losing much accuracy of molecular atomization energies., Comment: 12 pages, 9 figures

The fundamental energy gap of a periodic solid distinguishes insulators from metals and characterizes low-energy single-electron excitations. But 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 new 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 meta-generalized 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.

We compare the ability of four popular hybrid density functionals (B3LYP, B3PW91, HSE, and PBE0) for predicting band gaps of semiconductors and insulators over a large benchmark set using a consistent methodology. We observe no significant statistical dference in their overall performance although the screened hybrid HSE is more accurate for typical semiconductors. HSE can improve its accuracy for large large band gap materials --without affecting that of semiconductors-- by including a larger portion of Hartree-Fock exchange in its short range. Given that screened hybrids are computationally much less expensive than their global counterparts, we conclude that they are a better option for the black box prediction of band gaps.

Spin-projected Hartree-Fock is introduced as a particle-hole excitation ansatz over a symmetry-adapted reference determinant. Remarkably, this expansion has an analytic expression that we were able to decipher. While the form of the polynomial expansion is universal, the excitation amplitudes need to be optimized. This is equivalent to the optimization of orbitals in the conventional projected Hartree-Fock framework of non-orthogonal determinants. Using the inverse of the particle-hole expansion, we similarity transform the Hamiltonian in a coupled-cluster style theory. The left eigenvector of the non-hermitian Hamiltonian is constructed in a similar particle-hole expansion fashion, and we show that to numerically reproduce variational projected Hartree-Fock results, one needs as many pair excitations in the bra as the number of strongly correlated entangled pairs in the system. This single-excitation polynomial similarity transformation theory is an alternative to our recently presented double excitation theory, but supports projected Hartree-Fock and coupled cluster simultaneously rather than interpolating between them.

Recently, Tao and Mo developed a new nonempirical semilocal exchange-correlation density functional. The exchange part of this functional is derived from a density matrix expansion corrected to reproduce the fourth-order gradient expansion in the slowly varying limit, while the correlation part is based on the TPSS correlation model with a modification for the low-density limit. In the present work, the Tao-Mo functional is assessed by calculations on a variety of solids and jellium surfaces. This includes 22 lattice constants and bulk moduli, 7 cohesive energies, and jellium surface exchange and correlation energies for the density parameter rs in the range from 2 to 3 bohrs. Our calculations show that this meta-generalized gradient approximation can yield consistently remarkable accuracy for the properties considered here, with mean absolute errors of 0.017 {\AA} for lattice constants, 7.0 GPa for bulk moduli, 0.08 eV for cohesive energies, and 35 erg/cm2 for surface exchange-correlation energies, substantially improving upon existing nonempirical semilocal density functionals.