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.

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.

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.

Coupled Cluster and Perturbation Theories Based on a Cluster Mean-Field Reference Applied to Strongly Correlated Spin Systems

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 wave function is written as a simple tensor product of optimized cluster states. The cluster language and the mean-field nature of the ansatz allow for a straightforward improvement which uses perturbation theory and coupled-cluster to account for intercluster correlations. We present benchmark calculations on the 1D chain and 2D square J₁-J₂ 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.

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.

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.

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.

Singlet-paired coupled cluster theory for open shells

Restricted single-reference coupled cluster theory truncated to single and double excitations accurately describes weakly correlated systems, but often breaks down in the presence of static or strong correlation. Good coupled cluster energies in the presence of degeneracies can be obtained by using a symmetry-broken reference, such as unrestricted Hartree-Fock, but at the cost of good quantum numbers. A large body of work has shown that modifying the coupled cluster ansatz allows for the treatment of strong correlation within a single-reference, symmetry-adapted framework. The recently introduced singlet-paired coupled cluster doubles (CCD0) method is one such model, which recovers correct behavior for strong correlation without requiring symmetry breaking in the reference. Here, we extend singlet-paired coupled cluster for application to open shells via restricted open-shell singlet-paired coupled cluster singles and doubles (ROCCSD0). The ROCCSD0 approach retains the benefits of standard coupled cluster theory and recovers correct behavior for strongly correlated, open-shell systems using a spin-preserving ROHF reference.

Assessment of a Middle-Range Hybrid Functional

While hybrid functionals are largely responsible for the utility of modern Kohn-Sham density functional theory, they are not without their weaknesses. In the solid state, the slow decay of their nonlocal Hartree-Fock-type exchange makes hybrids computationally demanding and can introduce unphysical effects. Both problems can be remedied by a screened hybrid which uses exact exchange only at short-range. Many molecular properties, in contrast, benefit from the inclusion of long-range exact exchange. Recently, the authors reconciled these two seemingly contradictory requirements by introducing the HISS functional [ J. Chem. Phys. 2007 , 127 , 221103 ], which uses exact exchange only in the middle range. In this paper, we expand upon our previous work, benchmarking the performance of the HISS functional for several simple properties and applying it to the dissociation of homonuclear diatomic cations and to the polarizability of linear H2 chains to determine the importance of middle-range exact exchange for these systems, which are expected to be sensitive to the asymptotic exchange potential.

Range separated hybrids of pair coupled cluster doubles and density functionals

Pair coupled cluster doubles (pCCD) is a size-consistent, size-extensive, low-cost simplification of CCD that has been shown to be able to describe static correlation without breaking symmetry. We combine pCCD with Kohn-Sham functionals of the density and the local pair density in order to incorporate dynamic correlation in pCCD while maintaining its low cost. Double counting is eliminated by splitting the (interelectron) Coulomb operator into complementary short- and long-range parts, and evaluating the two-body energy with pCCD in the long-range and with density functionals in the short-range. This simultaneously suppresses self-interaction in the Hartree-exchange term of the functionals. Generalizations including a fraction of wavefunction two-body energy in the short-range are also derived and studied. The improvement of our pCCD+DFT hybrids over pCCD is demonstrated in calculations on benchmarks where both types of correlation are important.

Range-separated Brueckner coupled cluster doubles theory

We introduce a range-separation approximation to coupled cluster doubles (CCD) theory that successfully overcomes limitations of regular CCD when applied to the uniform electron gas. We combine the short-range ladder channel with the long-range ring channel in the presence of a Bruckner renormalized one-body interaction and obtain ground-state energies with an accuracy of 0.001 a.u./electron across a wide range of density regimes. Our scheme is particularly useful in the low-density and strongly correlated regimes, where regular CCD has serious drawbacks. Moreover, we cure the infamous overcorrelation of approaches based on ring diagrams (i.e., the particle-hole random phase approximation). Our energies are further shown to have appropriate basis set and thermodynamic limit convergence, and overall this scheme promises energetic properties for realistic periodic and extended systems which existing methods do not possess.

Improved semiconductor lattice parameters and band gaps from a middle-range screened hybrid exchange functional

We show that the middle-range exchange-correlation hybrid of Henderson, Izmaylov, Scuseria and Savin (HISS) performs extremely well for elemental and binary semiconductors with narrow or visible spectrum band gaps, as well as some wider gap or more ionic systems used in devices. The lattice parameters are superior to those predicted by the screened hybrid functional of Heyd, Scuseria and Ernzerhof (HSE), and provide a significant improvement over the geometries predicted by typical semilocal functionals, yielding results competitive with PBEsol, which was specially tuned for solids. HISS also yields band gaps superior to those produced by functionals developed specifically for the solid state. Timings indicate that HISS is more computationally efficient than HSE, implying that the high quality lattice constants coupled with improved optical band gap predictions render HISS a useful adjunct to HSE in the modeling of geometry-sensitive semiconductors.