BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Thomas M. Henderson

Thomas M. Henderson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.chem-phcond-mat.str-elcond-mat.mtrl-scicond-mat.othercond-mat.quant-gasnucl-thquant-ph

Catalog footprint

What is connected

28works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2024arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Thermal coupled cluster theory for SU(2) systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Geminal replacement models based on AGP

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On a dual representation of the Goldstone manifold

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 Φ| \hat{P}^q | Φ\rangle|$ becomes very small.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Geminal-Based Configuration Interaction

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Thermofield theory for finite-temperature coupled cluster

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Projected Coupled Cluster Theory: Optimization of cluster amplitudes in the presence of symmetry projection

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.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Influence of broken-pair excitations on the exact pair wavefunction

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.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

On the Magnetic Structure of Density Matrices

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.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Projected Coupled Cluster Theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Tensor-Structured Coupled Cluster Theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Polynomial Similarity Transformation Theory: A smooth interpolation between coupled cluster doubles and projected BCS applied to the reduced BCS Hamiltonian

We present a similarity transformation theory based on a polynomial form of a particle-hole pair excitation operator. In the weakly correlated limit, this polynomial becomes an exponential, leading to coupled cluster doubles. In the opposite strongly correlated limit, the polynomial becomes an extended Bessel expansion and yields the projected BCS wavefunction. In between, we interpolate using a single parameter. The effective Hamiltonian is non-hermitian and this Polynomial Similarity Transformation Theory follows the philosophy of traditional coupled cluster, left projecting the transformed Hamiltonian onto subspaces of the Hilbert space in which the wave function variance is forced to be zero. Similarly, the interpolation parameter is obtained through minimizing the next residual in the projective hierarchy. We rationalize and demonstrate how and why coupled cluster doubles is ill suited to the strongly correlated limit whereas the Bessel expansion remains well behaved. The model provides accurate wave functions with energy errors that in its best variant are smaller than 1\% across all interaction stengths. The numerical cost is polynomial in system size and the theory can be straightforwardly applied to any realistic Hamiltonian.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Projected Hartree Fock Theory as a Polynomial Similarity Transformation Theory of Single Excitations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Using full configuration interaction quantum Monte Carlo in a seniority zero space to investigate the correlation energy equivalence of pair coupled cluster doubles and doubly occupied configuration interaction

Over the past few years pair coupled cluster doubles (pCCD) has shown promise for the description of strong correlation. This promise is related to its apparent ability to match results from doubly occupied configuration interaction (DOCI), even though the latter method has exponential computational cost. Here, by modifying the full configuration interaction quantum Monte Carlo (FCIQMC) algorithm to sample only the seniority zero sector of Hilbert space, we show that the DOCI and pCCD energies are in agreement for a variety of 2D Hubbard models, including for systems well out of reach for conventional configuration interaction algorithms. Our calculations are aided by the sign problem being much reduced in the seniority zero space compared with the full space. We present evidence for this, and then discuss the sign problem in terms of the wave function of the system which appears to have a simplified sign structure.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Can single-reference coupled cluster theory describe static correlation?

While restricted single-reference coupled cluster theory truncated to singles and doubles (CCSD) provides very accurate results for weakly correlated systems, it usually fails in the presence of static or strong correlation. This failure is generally attributed to the qualitative breakdown of the reference, and can accordingly be corrected by using a multi-determinant reference, including higher-body cluster operators in the ansatz, or allowing symmetry breaking in the reference. None of these solutions are ideal; multi-reference coupled cluster is not black box, including higher-body cluster operators is computationally demanding, and allowing symmetry breaking leads to the loss of good quantum numbers. It has long been recognized that quasidegeneracies can instead be treated by modifying the coupled cluster ansatz. The recently introduced pair coupled cluster doubles (pCCD) approach is one such example which avoids catastrophic failures and accurately models strong correlations in a symmetry-adapted framework. Here we generalize pCCD to a singlet-paired coupled cluster model (CCD0) intermediate between coupled cluster doubles and pCCD, yielding a method that possesses the invariances of the former and much of the stability of the latter. Moreover, CCD0 retains the full structure of coupled cluster theory, including a fermionic wave function, antisymmetric cluster amplitudes, and well-defined response equations and density matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pair extended coupled cluster doubles

The accurate and efficient description of strongly correlated systems remains an important challenge for computational methods. Doubly occupied configuration interaction (DOCI), in which all electrons are paired and no correlations which break these pairs are permitted, can in many cases provide an accurate account of strong correlations, albeit at combinatorial computational cost. Recently, there has been significant interest in a method we refer to as pair coupled cluster doubles (pCCD), a variant of coupled cluster doubles in which the electrons are paired. This is simply because pCCD provides energies nearly identical to those of DOCI, but at mean-field computational cost (disregarding the cost of the two-electron integral transformation). Here, we introduce the more complete pair extended coupled cluster doubles (pECCD) approach which, like pCCD, has mean-field cost and reproduces DOCI energetically. We show that unlike pCCD, pECCD also reproduces the DOCI wave function with high accuracy. Moreoever, pECCD yields sensible albeit inexact results even for attractive interactions where pCCD breaks down.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

Numerical results for ground state and excited state properties (energies, double occupancies, and Matsubara-axis self energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed node approximation, unrestricted coupled cluster theory, and multi-reference projected Hartree-Fock. Comparison of results obtained by different methods allows for the identification of uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic limit values is emphasized. Cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Lie Algebraic Similarity Transformed Hamiltonians for Lattice Model Systems

We present a class of Lie algebraic similarity transformations generated by exponentials of two-body on-site hermitian operators whose Hausdorff series can be summed exactly without truncation. The correlators are defined over the entire lattice and include the Gutzwiller factor $n_{i\uparrow}n_{i\downarrow}$, and two-site products of density $(n_{i\uparrow} + n_{i\downarrow})$ and spin $(n_{i\uparrow}-n_{i\downarrow})$ operators. The resulting non-hermitian many-body Hamiltonian can be solved in a biorthogonal mean-field approach with polynomial computational cost. The proposed similarity transformation generates locally weighted orbital transformations of the reference determinant. Although the energy of the model is unbound, projective equations in the spirit of coupled cluster theory lead to well-defined solutions. The theory is tested on the 1D and 2D repulsive Hubbard model where we find accurate results across all interaction strengths.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Seniority-based coupled cluster theory

Doubly occupied configuration interaction (DOCI) with optimized orbitals often accurately describes strong correlations while working in a Hilbert space much smaller than that needed for full configuration interaction. However, the scaling of such calculations remains combinatorial with system size. Pair coupled cluster doubles (pCCD) is very successful in reproducing DOCI energetically, but can do so with low polynomial scaling ($N^3$, disregarding the two-electron integral transformation from atomic to molecular orbitals). We show here several examples illustrating the success of pCCD in reproducing both the DOCI energy and wave function, and show how this success frequently comes about. What DOCI and pCCD lack are an effective treatment of dynamic correlations, which we here add by including higher-seniority cluster amplitudes which are excluded from pCCD. This frozen pair coupled cluster approach is comparable in cost to traditional closed-shell coupled cluster methods with results that are competitive for weakly correlated systems and often superior for the description of strongly correlated systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Particle-particle and quasiparticle random phase approximations: Connections to coupled cluster theory

We establish a formal connection between the particle-particle (pp) random phase approximation (RPA) and the ladder channel of the coupled cluster doubles (CCD) equations. The relationship between RPA and CCD is best understood within a Bogoliubov quasiparticle (qp) RPA formalism. This work is a follow-up to our previous formal proof on the connection between particle-hole (ph) RPA and ring-CCD. Whereas RPA is a quasibosonic approximation, CC theory is a correct bosonization in the sense that the wavefunction and Hilbert space are exactly fermionic. Coupled cluster theory achieves this goal by interacting the ph (ring) and pp (ladder) diagrams via a third channel that we here call "crossed-ring" whose presence allows for full fermionic antisymmetry. Additionally, coupled cluster incorporates what we call "mosaic" terms which can be absorbed into defining a new effective one-body Hamiltonian. The inclusion of these mosaic terms seems to be quite important. The pp-RPA an d qp-RPA equations are textbook material in nuclear structure physics but are largely unknown in quantum chemistry, where particle number fluctuations and Bogoliubov determinants are rarely used. We believe that the ideas and connections discussed in this paper may help design improved ways of incorporating RPA correlation into density functionals based on a CC perspective.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Projected Hartree-Fock Theory

Projected Hartree-Fock theory (PHF) has a long history in quantum chemistry. PHF is here understood as the variational determination of an N-electron broken symmetry Slater determinant that minimizes the energy of a projected state with the correct quantum numbers. The method was actively pursued for several decades but seems to have been abandoned. We here derive and implement a "variation after projection" PHF theory using techniques different from those previously employed in quantum chemistry. Our PHF methodology has modest mean-field computational cost, yields relatively simple expressions, can be applied to both collinear and non-collinear spin cases, and can be used in conjunction with deliberate symmetry breaking and restoration of other molecular symmetries like complex conjugation and point group. We present several benchmark applications to dissociation curves and singlet-triplet energy splittings, showing that the resulting PHF wavefunctions are of high quality multireference character. We also provide numerical evidence that in the thermodynamic limit, the energy in PHF is not lower than that of broken-symmetry HF, a simple consequence of the lack of size consistency and extensivity of PHF.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Improved Semiconductor Lattice Parameters and Band Gaps from a Middle-Range Screened Hybrid 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 commercially. 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 geometries predicted by semilocal functionals. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Projected Quasiparticle Theory for Molecular Electronic Structure

We derive and implement symmetry-projected Hartree-Fock-Bogoliubov (HFB) equations and apply them to the molecular electronic structure problem. All symmetries (particle number, spin, spatial, and complex conjugation) are deliberately broken and restored in a self-consistent variation-after-projection approach. We show that the resulting method yields a comprehensive black-box treatment of strong correlations with effective one-electron (mean-field) computational cost. The ensuing wave function is of multireference character and permeates the entire Hilbert space of the problem. The energy expression is different from regular HFB theory but remains a functional of an independent quasiparticle density matrix. All reduced density matrices are expressible as an integration of transition density matrices over a gauge grid. We present several proof-of-principle examples demonstrating the compelling power of projected quasiparticle theory for electronic structure theory.

0 citations0 reviewsNo code linkedOpen access