Source author record

Gustavo E. Scuseria

Gustavo E. Scuseria appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.str-el physics.chem-ph cond-mat.mtrl-sci cond-mat.other physics.comp-ph cond-mat.quant-gas quant-ph cond-mat.mes-hall nucl-th cond-mat.supr-con

Catalog footprint

What is connected

70works

10topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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.

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.

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.

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.

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.

preprint2020arXiv

Correlating AGP on a quantum computer

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.

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.

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.

preprint2019arXiv

Exact Parameterization of Fermionic Wave Functions via Unitary Coupled Cluster Theory

A formal analysis is conducted on the exactness of various forms of unitary coupled cluster (UCC) theory based on particle-hole excitation and de-excitation operators. Both the conventional single exponential UCC parameterization and a 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.

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.

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.

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.

preprint2018arXiv

On the difference between variational and unitary coupled cluster theories

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.

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.

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.

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.

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.

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.

preprint2016arXiv

Biorthogonal projected energies of a Gutzwiller similarity transformed Hamiltonian

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.

preprint2016arXiv

Performance of a Nonempirical Density Functional on Molecules and Hydrogen-Bonded Complexes

Recently, Tao and Mo (TM) derived a meta-generalized gradient approximation functional based on a model exchange-correlation hole. In this work, the performance of this functional is assessed on standard test sets, using the 6-311++G(3df,3pd) basis set. These test sets include 223 G3/99 enthalpies of formation, 99 atomization energies, 76 barrier heights, 58 electron affinities, 8 proton affinities, 96 bond lengths, 82 harmonic vibrational frequencies, 10 hydrogen-bonded molecular complexes, and 22 atomic excitation energies. Our calculations show that the TM functional can achieve high accuracy for most properties considered, relative to the LSDA, PBE, and TPSS functionals. In particular, it yields the best accuracy for proton affinities, harmonic vibrational frequencies, hydrogen-bonded dissociation energies and bond lengths, and atomic excitation energies.

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.

preprint2016arXiv

Predicting Band Gaps with Hybrid Density Functionals

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.

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.

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.

preprint2015arXiv

A cluster-based mean-field and perturbative description of strongly correlated fermion systems. Application to the 1D and 2D Hubbard model

We introduce a mean-field and perturbative approach, based on clusters, to describe the ground state of fermionic strongly-correlated systems. In cluster mean-field, the ground state wavefunction is written as a simple tensor product over optimized cluster states. The optimization of the single-particle basis where the cluster mean-field is expressed is crucial in order to obtain high-quality results. The mean-field nature of the ansatz allows us to formulate a perturbative approach to account for inter-cluster correlations; other traditional many-body strategies can be easily devised in terms of the cluster states. We present benchmark calculations on the half-filled 1D and (square) 2D Hubbard model, as well as the lightly-doped regime in 2D, using cluster mean-field and second-order perturbation theory. Our results indicate that, with sufficiently large clusters or to second-order in perturbation theory, a cluster-based approach can provide an accurate description of the Hubbard model in the considered regimes. Several avenues to improve upon the results presented in this work are discussed.

preprint2015arXiv

Actinide chemistry using singlet-paired coupled cluster and its combinations with density functionals

Singlet-paired coupled cluster doubles (CCD0) is a simplification of CCD that relinquishes a fraction of dynamic correlation in order to be able to describe static correlation. Combinations of CCD0 with density functionals that recover specifically the dynamic correlation missing in the former have also been developed recently. Here, we assess the accuracy of CCD0 and CCD0+DFT (and variants of these using Brueckner orbitals) as compared to well-established quantum chemical methods for describing ground-state properties of singlet actinide molecules. The $f^0$ actinyl series (UO$_2^{2+}$, NpO$_2^{2+}$, PuO$_2^{2+}$), the isoelectronic NUN, and Thorium (ThO, ThO$^{2+}$) and Nobelium (NoO, NoO$_2$) oxides are studied.

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.

preprint2015arXiv

Combinations of coupled cluster, density functionals, and the random phase approximation for describing static and dynamic correlation, and van der Waals interactions

Contrary to standard coupled cluster doubles (CCD) and Brueckner doubles (BD), singlet-paired analogues of CCD and BD (denoted here as CCD0 and BD0) do not break down when static correlation is present, but neglect substantial amounts of dynamic correlation. In fact, CCD0 and BD0 do not account for any contributions from multielectron excitations involving only same-spin electrons at all. We exploit this feature to add---without introducing double counting, self-interaction, or increase in cost---the missing correlation to these methods via meta-GGA density functionals (TPSS and SCAN). Furthermore, we improve upon these CCD0+DFT blends by invoking range separation: the short- and long-range correlations absent in CCD0/BD0 are evaluated with DFT and the direct random phase approximation (dRPA), respectively. This corrects the description of long-range van der Waals forces. Comprehensive benchmarking shows that the combinations presented here are very accurate for weakly correlated systems, while also providing a reasonable description of strongly correlated problems without resorting to symmetry breaking.

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.

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.

preprint2014arXiv

Composite fermion-boson mapping for fermionic lattice models

We present a mapping of elementary fermion operators onto a quadratic form of composite fermionic and bosonic operators. The mapping is an exact isomorphism as long as the physical constraint of one composite particle per cluster is satisfied. This condition is treated on average in a composite particle mean-field approach, which consists of an ansatz that decouples the composite fermionic and bosonic sectors. The theory is tested on the one- and two-dimensional Hubbard models. Using a Bogoliubov determinant for the composite fermions and either a coherent or Bogoliubov state for the bosons, we obtain a simple and accurate procedure for treating the Mott insulating phase of the Hubbard model with mean-field computational cost.

preprint2014arXiv

Coupled Cluster Channels in the Homogeneous Electron Gas

We discuss diagrammatic modifications to the coupled cluster doubles (CCD) equations, wherein different groups of terms out of rings, ladders, crossed-rings and mosaics can be removed to form approximations to the coupled cluster method, of interest due to their similarity with various types of random phase approximations. The finite uniform electron gas is benchmarked for 14- and 54-electron systems at the complete basis set limit over a wide density range and performance of different flavours of CCD are determined. These results confirm that rings generally overcorrelate and ladders generally undercorrelate; mosaics-only CCD yields a result surprisingly close to CCD. We use a recently developed numerical analysis [J. J. Shepherd and A. Grüneis, Phys. Rev. Lett. 110, 226401 (2013)] to study the behaviours of these methods in the thermodynamic limit. We determine that the mosaics, on forming the Brueckner Hamltonian, open a gap in the effective one-particle eigenvalues at the Fermi energy. Numerical evidence is presented which shows that methods based on this renormalisation have convergent energies in the thermodynamic limit including mosaic-only CCD, which is just a renormalised MP2. All other methods including only a single channel, namely ladder-only CCD, ring-only CCD and crossed-ring-only CCD, appear to yield divergent energies; incorporation of mosaic terms prevents this from happening.

preprint2014arXiv

Electron correlation in solids via density embedding theory

Density matrix embedding theory (Phys. Rev. Lett. 109, 186404 (2012)) and density embedding theory ((Phys. Rev. B 89, 035140 (2014)) have recently been introduced for model lattice Hamiltonians and molecular systems. In the present work, the formalism is extended to the ab initio description of infinite systems. An appropriate definition of the impurity Hamiltonian for such systems is presented and demonstrated in cases of 1, 2 and 3 dimensions, using coupled cluster theory as the impurity solver. Additionally, we discuss the challenges related to disentanglement of fragment and bath states. The current approach yields results comparable to coupled cluster calculations of infinite systems even when using a single unit cell as the fragment. The theory is formulated in the basis of Wannier functions but it does not require separate localization of unoccupied bands. The embedding scheme presented here is a promising way of employing highly accurate electronic structure methods for extended systems at a fraction of their original computational cost.

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.

preprint2014arXiv

Multi-reference symmetry-projected variational approximation for the ground state of the doped one-dimensional Hubbard model

A multi-reference configuration mixing scheme is used to describe the ground state, characterized by well defined spin and space group symmetry quantum numbers as well as doping fractions $N_{e}/N_{sites}$, of one dimensional Hubbard lattices with nearest-neighbor hopping and periodic boundary conditions. Within this scheme, each ground state is expanded in a given number of nonorthogonal and variationally determined symmetry-projected configurations. The results obtained for the ground state and correlation energies of half-filled and doped lattices with 30, 34 and 50 sites, compare well with the exact Lieb-Wu solutions as well as with the ones obtained with other state-of-the-art approximations. The structure of the intrinsic symmetry-broken determinants resulting from the variational procedure is interpreted in terms of solitons whose translational and breathing motions can be regarded as basic units of quantum fluctuations. It is also shown that in the case of doped 1D lattices, a part of such fluctuations can also be interpreted in terms of polarons. In addition to momentum distributions, both spin-spin and density-density correlation functions are studied as functions of doping. The spectral functions and density of states, computed with an ansatz whose quality can be well-controlled by the number of symmetry-projected configurations used to approximate the $N_{e} \pm 1$ electron systems, display features beyond a simple quasiparticle distribution, as well as spin-charge separation trends.

preprint2014arXiv

On the equivalence of LIST and DIIS methods for convergence acceleration

Self-consistent field extrapolation methods play a pivotal role in quantum chemistry and electronic structure theory. We here demonstrate the mathematical equivalence between the recently proposed family of LIST methods [J. Chem. Phys. 134, 241103 (2011); J. Chem. Theory Comput. 7, 3045 (2011)] with Pulay's DIIS [Chem. Phys. Lett. 73, 393 (1980)]. Our results also explain the differences in performance among the various LIST methods.

preprint2014arXiv

Polyradical character and spin frustration in fullerene molecules: An ab initio non-collinear Hartree--Fock study

Most {\em ab initio} calculations on fullerene molecules have been carried out based on the paradigm of the Hückel model. This is consistent with the restricted nature of the independent-particle model underlying such calculations, even in single-reference-based correlated approaches. On the other hand, previous works on some of these molecules using model Hamiltonians have clearly indicated the importance of short-range inter-atomic spin-spin correlations. In this work, we consider {\em ab initio} non-collinear Hartree--Fock (HF) solutions for representative fullerene systems: the bowl, cage, ring, and pentagon isomers of C$_{20}$, and the larger C$_{30}$, C$_{36}$, C$_{60}$, C$_{70}$, and C$_{84}$ fullerene cages. In all cases but the ring we find that the HF minimum corresponds to a truly non-collinear solution with a torsional spin density wave. Optimized geometries at the generalized HF (GHF) level lead to fully symmetric structures, even in those cases where Jahn-Teller distortions have been previously considered. The nature of the GHF solutions is consistent with the $π$-electron space becoming polyradical in nature: each $p$-orbital remains effectively singly occupied. The spin frustration, induced by the pentagon rings in an otherwise anti-ferromagnetic background, is minimized at the HF level by aligning the spins in non-collinear arrangements. The long-range magnetic ordering observed is reminiscent of the character of broken symmetry HF solutions in polyacene systems.

preprint2014arXiv

Quasiparticle Coupled Cluster Theory for Pairing Interactions

We present an extension of the pair coupled cluster doubles (p-CCD) method to quasiparticles and apply it to the attractive pairing Hamiltonian. Near the transition point where number symmetry gets spontaneously broken, the proposed BCS-based p-CCD method yields significantly better energies than existing methods when compared to exact results obtained via solution of the Richardson equations. The quasiparticle p-CCD method has a low computational cost of $\mathcal{O}(N^3)$ as a function of system size. This together with the high quality of results here demonstrated, points to considerable promise for the accurate description of strongly correlated systems with more realistic pairing interactions.

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.

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.

preprint2014arXiv

Symmetry-projected Wave Functions in Quantum Monte Carlo Calculations

We consider symmetry-projected Hartree--Fock trial wave functions in constrained-path Monte Carlo (CPMC) calculations. Previous CPMC calculations have mostly employed Hartree--Fock (HF) trial wave functions, restricted or unrestricted. The symmetry-projected HF approach results in a hierarchy of wave functions with increasing quality: the more symmetries that are broken and restored in a self-consistent manner, the higher the quality of the trial wave function. This hierarchy is approximately maintained in CPMC calculations: the accuracy in the energy increases and the statistical variance decreases when further symmetries are broken and restored. Significant improvement is achieved in CPMC with the best symmetry-projected trial wave functions over those from simple HF. We analyze and quantify the behavior using the two-dimensional repulsive Hubbard model as an example. In the sign-problem-free region, where CPMC can be made exact but a constraint is deliberately imposed here, spin-projected wave functions remove the constraint bias. Away from half-filling, spatial symmetry restoration in addition to that of the spin lead to highly accurate results from CPMC. Since the computational cost of symmetry-projected HF trial wave functions in CPMC can be made to scale algebraically with system size, this provides a potentially general approach for accurate calculations in many-fermion systems.

preprint2014arXiv

The sign problem in full configuration interaction quantum Monte Carlo: Linear and sub-linear representation regimes for the exact wave function

We investigate the sign problem for full configuration interaction quantum Monte Carlo (FCIQMC), a stochastic algorithm for finding the ground state solution of the Schrödinger equation with substantially reduced computational cost compared with exact diagonalisation. We find $k$-space Hubbard models for which the solution is yielded with storage that grows sub-linearly in the size of the many-body Hilbert space, in spite of using a wave function that is simply linear combination of states. The FCIQMC algorithm is able to find this sub-linear scaling regime without bias and with only a choice of Hamiltonian basis. By means of a demonstration we solve for the energy of a 70-site half-filled system (with a space of $10^{38}$ determinants) in 250 core hours, substantially quicker than the $\sim$10$^{36}$ core hours that would be required by exact diagonalisation. This is the largest space that has been sampled in an unbiased fashion. The challenge for the recently-developed FCIQMC method is made clear: expand the sub-linear scaling regime whilst retaining exact on average accuracy. This result rationalizes the success of the initiator adaptation (i-FCIQMC) and offers clues to improve it. We argue that our results changes the landscape for development of FCIQMC and related methods.

preprint2014arXiv

Variational description of the ground state of the repulsive two-dimensional Hubbard model in terms of nonorthogonal symmetry-projected Slater determinants

The few determinant (FED) methodology, introduced in our previous works for 1D lattices, is here adapted for the repulsive two-dimensional Hubbard model at half-filling and with finite doping fractions. Within this configuration mixing scheme, a given ground state with well defined spin and space group quantum numbers, is expanded in terms of a nonorthogonal symmetry-projected basis determined through chains of variation-after projection calculations. The results obtained for the ground state and correlation energies of half-filled and doped 4$\times$4, 6$\times$6, 8$\times$8, and 10$\times$10 lattices, as well as momentum distributions and spin-spin correlation functions in small lattices, compare well with those obtained using other state-of-the-art approximations. The structure of the intrinsic determinants resulting from the variational strategy is interpreted in terms of defects that encode information on the basic units of quantum fluctuations in the considered 2D systems. The varying nature of the underlying quantum fluctuations, reflected in a transition to a stripe regime for increasing onsite repulsions, is discussed using the intrinsic determinants belonging to a 16$\times$4 lattice with 56 electrons. Such a transition is further illustrated by computing spin-spin and charge-charge correlation functions with the corresponding multireference FED wave functions. In good agreement with previous studies, the analysis of the pairing correlation functions reveals a weak enhancement of the extended $s$-wave and d$_{x^2-y^2}$ pairing modes. Given the quality of results here reported together with those previously obtained for 1D lattices and the parallelization properties of the FED scheme, we believe that symmetry projection techniques are very well suited for building ground state wave functions of correlated electronic systems, regardless of their dimensionality.

preprint2013arXiv

Composite Boson Mapping for Lattice Boson Systems

We present a canonical mapping transforming physical boson operators into quadratic products of cluster composite bosons that preserves matrix elements of operators when a physical constraint is enforced. We map the 2D lattice Bose-Hubbard Hamiltonian into $2\times 2$ composite bosons and solve it at mean field. The resulting Mott insulator-superfluid phase diagram reproduces well Quantum Monte Carlo results. The Higgs boson behavior along the particle-hole symmetry line is unraveled and in remarkable agreement with experiment. Results for the properties of the ground and excited states are competitive with other state-of-the-art approaches, but at a fraction of their computational cost. The composite boson mapping here introduced can be readily applied to frustrated many-body systems where most methodologies face significant hurdles.

preprint2013arXiv

Density Functionals that Recognize Covalent, Metallic, and Weak Bonds

Computationally-efficient semilocal approximations of density functional theory at the level of the local spin density approximation (LSDA) or generalized gradient approximation (GGA) poorly describe weak interactions. We show improved descriptions for weak bonds (without loss of accuracy for strong ones) from a newly-developed semilocal meta-GGA (MGGA), by applying it to molecules, surfaces, and solids. We argue that this improvement comes from using the right MGGA dimensionless ingredient to recognize all types of orbital overlap.

preprint2013arXiv

Density Matrix Embedding from Broken Symmetry Lattice Mean-Fields

Several variants of the recently proposed Density Matrix Embedding Theory (DMET) [G. Knizia and G. K-L. Chan, Phys. Rev. Lett. 109, 186404 (2012)] are formulated and tested. We show that spin symmetry breaking of the lattice mean-field allows precise control of the lattice and fragment filling while providing very good agreement between predicted properties and exact results. We present a rigorous proof that at convergence this method is guaranteed to preserve lattice and fragment filling. Differences arising from fitting the fragment one-particle density matrix alone versus fitting fragment plus bath are scrutinized. We argue that it is important to restrict the density matrix fitting to solely the fragment. Furthermore, in the proposed broken symmetry formalism, it is possible to substantially simplify the embedding procedure without sacrificing its accuracy by resorting to density instead of density matrix fitting. This simplified Density Embedding Theory (DET) greatly improves the convergence properties of the algorithm.

preprint2013arXiv

Excited electronic states from a variational approach based on symmetry-projected Hartree--Fock configurations

Recent work from our research group has demonstrated that symmetry-projected Hartree--Fock (HF) methods provide a compact representation of molecular ground state wavefunctions based on a superposition of non-orthogonal Slater determinants. The symmetry-projected ansatz can account for static correlations in a computationally efficient way. Here we present a variational extension of this methodology applicable to excited states of the same symmetry as the ground state. Benchmark calculations on the C$_2$ dimer with a modest basis set, which allows comparison with full configuration interaction results, indicate that this extension provides a high quality description of the low-lying spectrum for the entire dissociation profile. We apply the same methodology to obtain the full low-lying vertical excitation spectrum of formaldehyde, in good agreement with available theoretical and experimental data, as well as to a challenging model $C_{2v}$ insertion pathway for BeH$_2$. The variational excited state methodology developed in this work has two remarkable traits: it is fully black-box and will be applicable to fairly large systems thanks to its mean-field computational cost.

preprint2013arXiv

Multi-component symmetry-projected approach for molecular ground state correlations

The symmetry-projected Hartree--Fock ansatz for the electronic structure problem can efficiently account for static correlation in molecules, yet it is often unable to describe dynamic correlation in a balanced manner. Here, we consider a multi-component, systematically-improvable approach, that accounts for all ground state correlations. Our approach is based on linear combinations of symmetry-projected configurations built out of a set of non-orthogonal, variationally optimized determinants. The resulting wavefunction preserves the symmetries of the original Hamiltonian even though it is written as a superposition of deformed (broken-symmetry) determinants. We show how short expansions of this kind can provide a very accurate description of the electronic structure of simple chemical systems such as the nitrogen and the water molecules, along the entire dissociation profile. In addition, we apply this multi-component symmetry-projected approach to provide an accurate interconversion profile among the peroxo and bis($μ$-oxo) forms of [Cu$_2$O$_2$]$^{2+}$, comparable to other state-of-the-art quantum chemical methods.

preprint2013arXiv

Multi-reference symmetry-projected variational approaches for ground and excited states of the one-dimensional Hubbard model

We present a multi-reference configuration mixing scheme for describing ground and excited states, with well defined spin and space group symmetry quantum numbers, of the one-dimensional Hubbard model with nearest-neighbor hopping and periodic boundary conditions. Within this scheme, each state is expanded in terms of non-orthogonal and variationally determined symmetry-projected configurations. The results for lattices up to 30 and 50 sites compare well with the exact Lieb-Wu solutions as well as with results from other state-of-the-art approximations. In addition to spin-spin correlation functions in real space and magnetic structure factors, we present results for spectral functions and density of states computed with an ansatz whose quality can be well-controlled by the number of symmetry-projected configurations used to approximate the systems with $N_{e}$ and $N_{e} \pm 1$ electrons. The intrinsic symmetry-broken determinants resulting from the variational calculations have rich structures in terms of defects that can be regarded as basic units of quantum fluctuations. Given the quality of the results here reported, as well as the parallelization properties of the considered scheme, we believe that symmetry-projection techniques, which have found ample applications in nuclear structure physics, deserve further attention in the study of low-dimensional correlated many-electron systems.

preprint2013arXiv

On pair functions for strong correlations

The UHF wave function may be written as a spin-contaminated \textit{pair} wave function of the APSG form, and the overlap of the alpha and beta corresponding orbitals of the UHF solution can be taken as a proxy for the strength of the correlation captured by breaking symmetry. We demonstrate this with calculations on one- and two-dimensional hydrogen clusters and make contact with the well studied Hubbard model. The UHF corresponding orbitals pair in a manner that allows a smooth evolution from doubly occupied orbitals at small distance to one in which wave function breaks symmetry, segregating the $α$ and $β$ electrons onto distinct sublattices at large distances. By performing spin projection on these UHF solutions, we address strong correlations that are difficult to capture at intermediate distances using a single determinant. Approved for public release: LA-UR-13-22691.

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.

preprint2013arXiv

Proper and improper zero energy modes in Hartree-Fock theory and their relevance for symmetry breaking and restoration

We study the spectra of the molecular orbital Hessian (stability matrix) and random-phase approximation Hamiltonian of broken-symmetry Hartree-Fock solutions, focusing on zero eigenvalue modes. After all negative eigenvalues are removed from the Hessian by following their eigenvectors downhill, one is left with only positive and zero eigenvalues. Zero modes correspond to orbital rotations with no restoring force. These rotations determine states in the Goldstone manifold, which originates from a spontaneously broken continuous symmetry in the wave function. Zero modes can be classified as improper or proper according to their different mathematical and physical properties. Improper modes arise from symmetry breaking and their restoration always lowers the energy. Proper modes, on the other hand, correspond to degeneracies of the wave function, and their symmetry restoration does not necessarily lower the energy. We discuss how the RPA Hamiltonian distinguishes between proper and improper modes by doubling the number of zero eigenvalues associated with the latter. Proper modes in the Hessian always appear in pairs which do not double in RPA. We present several pedagogical cases exemplifying the above statements. The relevance of these results for projected Hartree-Fock methods is also addressed.

preprint2013arXiv

Renormalized Second-order Perturbation Theory for The Electron Correlation Energy: Concept, Implementation, and Benchmarks

We present a renormalized second-order perturbation theory (rPT2), based on a Kohn-Sham (KS) reference state, for the electron correlation energy that includes the random-phase approximation (RPA), second-order screened exchange (SOSEX), and renormalized single excitations (rSE). These three terms all involve a summation of certain types of diagrams to infinite order, and can be viewed as "renormalization" of the 2nd-order direct, exchange, and single excitation (SE) terms of Rayleigh-Schrödinger perturbation theory based on an KS reference. In this work we establish the concept of rPT2 and present the numerical details of our SOSEX and rSE implementations. A preliminary version of rPT2, in which the renormalized SE (rSE) contribution was treated approximately, has already been benchmarked for molecular atomization energies and chemical reaction barrier heights and shows a well balanced performance [Paier et al, New J. Phys. 14, 043002 (2012)]. In this work, we present a refined version of rPT2, in which we evaluate the rSE series of diagrams rigorously. We then extend the benchmark studies to non-covalent interactions, including the rare-gas dimers, and the S22 and S66 test sets. Despite some remaining shortcomings, we conclude that rPT2 gives an overall satisfactory performance across different chemical environments, and is a promising step towards a generally applicable electronic structure approach.

preprint2013arXiv

Semilocal and Hybrid Meta-Generalized Gradient Approximations Based on the Understanding of the Kinetic-Energy-Density Dependence

We present a global hybrid meta-generalized gradient approximation (meta-GGA) with three empirical parameters, as well as its underlying semilocal meta-GGA and a meta-GGA with only one empirical parameter. All of them are based on the new meta-GGA resulting from the understanding of kinetic-energy-density dependence [J. Chem. Phys. 137, 051101 (2012)]. The obtained functionals show robust performances on the considered molecular systems for the properties of heats of formation, barrier heights, and noncovalent interactions. The pair-wise additive dispersion corrections to the functionals are also presented.

preprint2013arXiv

The role of screened exact exchange in accurately describing properties of transition metal oxides: Modeling defects in LaAlO3

The properties of many intrinsic defects in the wide band gap semiconductor LaAlO3 are studied using the screened hybrid functional of Heyd, Scuseria, and Ernzerhof (HSE). As in pristine structures, exact exchange included in the screened hybrid functional alleviates the band gap underestimation problem, which is common to semilocal functionals; this allows accurate prediction of defect properties. We propose correction-free defect energy levels for bulk LaAlO3 computed using HSE that might serve as guide in the interpretation of photoluminescence experiments.

preprint2012arXiv

Assessment of correlation energies based on the random-phase approximation

The random-phase approximation to the ground state correlation energy (RPA) in combination with exact exchange (EX) has brought Kohn-Sham (KS) density functional theory one step closer towards a universal, "general purpose first principles method". In an effort to systematically assess the influence of several correlation energy contributions beyond RPA, this work presents dissociation energies of small molecules and solids, activation energies for hydrogen transfer and non-hydrogen transfer reactions, as well as reaction energies for a number of common test sets. We benchmark EX+RPA and several flavors of energy functionals going beyond it: second-order screened exchange (SOSEX), single excitation (SE) corrections, renormalized single excitation (rSE) corrections, as well as their combinations. Both the single excitation correction as well as the SOSEX contribution to the correlation energy significantly improve upon the notorious tendency of EX+RPA to underbind. Surprisingly, activation energies obtained using EX+RPA based on a KS reference alone are remarkably accurate. RPA+SOSEX+rSE provides an equal level of accuracy for reaction as well as activation energies and overall gives the most balanced performance, which makes it applicable to a wide range of systems and chemical reactions.

preprint2012arXiv

N-electron Slater determinants from non-unitary canonical transformations of fermion operators

Mean-field methods such as Hartree-Fock (HF) or Hartree-Fock-Bogoliubov (HFB) constitute the building blocks upon which more elaborate many-body theories are based on. The HF and HFB wavefunctions are built out of independent quasi-particles resulting from a unitary linear canonical transformation of the elementary fermion operators. Here, we discuss the possibility of allowing the HF transformation to become non-unitary. The properties of such HF vacua are discussed, as well as the evaluation of matrix elements among such states. We use a simple ansatz to demonstrate that a non-unitary transformation brings additional flexibility that can be exploited in variational approximations to many-fermion wavefunctions. The action of projection operators on non-unitary based HF states is also discussed and applied to the one-dimensional Hubbard model with periodic boundary conditions.

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.

preprint2012arXiv

Structural Phase Transitions of the Metal Oxide Perovskites SrTiO3, LaAlO3 and LaTiO3 Studied with a Screened Hybrid Functional

We have investigated the structural phase transitions of the transition metal oxide perovskites SrTiO$_{3}$, LaAlO$_{3}$ and LaTiO$_{3}$ using the screened hybrid density functional of Heyd, Scuseria and Ernzerhof (HSE06). We show that HSE06-computed lattice parameters, octahedral tilts and rotations, as well as electronic properties, are significantly improved over semilocal functionals. We predict the crystal field splitting ($Δ_{CF}$) resulting from the structural phase transition in SrTiO$_{3}$ and LaAlO$_{3}$ to be 3 meV and 10 meV, respectively, in excellent agreement with experimental results. HSE06 identifies correctly LaTiO$_{3}$ in the magnetic sates as a Mott insulator. Also, it predicts that the GdFeO$_{3}$-type distortion in non-magnetic LaTiO$_{3}$ will induce a large $Δ_{CF}$ of 410 meV. This large crystal-field splitting associated with the large magnetic moment found in the G-type antiferromagnetic state suggest that LaTiO$_{3}$ has an induced orbital order, which is confirmed by the visualisation of the highest occupied orbitals. These results strongly indicate that HSE06 is capable of efficiently and accurately modeling perovskite oxides, and promises to efficiently capture the physics at their heterointerfaces.

preprint2012arXiv

Symmetry-projected variational approach for ground and excited states of the two-dimensional Hubbard model

We present a symmetry-projected configuration mixing scheme to describe ground and excited states, with well defined quantum numbers, of the two-dimensional Hubbard model with nearestneighbor hopping and periodic boundary conditions. Results for the half-filled 2{\times}4, 4{\times}4, and 6{\times}6 lattices, as well as doped 4 {\times} 4 systems, compare well with available results, both exact and from other state-of-the-art approximations. We report spectral functions and density of states obtained from a well-controlled ansatz for the (Ne {\pm} 1)-electron system. Symmetry projected methods have been widely used for the many-body nuclear physics problem but have received little attention in the solid state community. Given their relatively low (mean-field) computational cost and the high quality of results here reported, we believe that they deserve further scrutiny.

preprint2011arXiv

Accurate Modeling of the Cubic and Antiferrodistortive Phases of SrTiO3 with Screened Hybrid Density Functional Theory

We have calculated the properties of SrTiO3 (STO) using a wide array of density functionals ranging from standard semi-local functionals to modern range-separated hybrids, combined with several basis sets of varying size/quality. We show how these combination's predictive ability varies significantly, both for STO's cubic and antiferrodistortive (AFD) phases, with the greatest variation in functional/basis set efficacy seen in modeling the AFD phase. The screened hybrid functionals we utilized predict the structural properties of both phases in very good agreement with experiment, especially if used with large (but still computationally tractable) basis sets. The most accurate results presented in this study, namely those from HSE06/modified-def2-TZVP, stand as the most accurate modeling of STO to date when compared to the literature; these results agree well with experimental structural and electronic properties as well as providing insight into the band structure alteration during the phase transition.

preprint2011arXiv

Computing the energy of a water molecule using MultiDeterminants: A simple, efficient algorithm

Quantum Monte Carlo (QMC) methods such as variational Monte Carlo and fixed node diffusion Monte Carlo depend heavily on the quality of the trial wave function. Although Slater-Jastrow wave functions are the most commonly used variational ansatz in electronic structure, more sophisticated wave-functions are critical to ascertaining new physics. One such wave function is the multiSlater-Jastrow wave function which consists of a Jastrow function multiplied by the sum of Slater determinants. In this paper we describe a method for working with these wavefunctions in QMC codes that is easy to implement, efficient both in computational speed as well as memory, and easily parallelized. The computational cost scales quadratically with particle number making this scaling no worse than the single determinant case and linear with the total number of excitations. Additionally we implement this method and use it to compute the ground state energy of a water molecule.

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.

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.

preprint2010arXiv

ROHF Theory Made Simple

Restricted open-shell Hartree-Fock (ROHF) theory is formulated as a projected self-consistent unrestricted HF (UHF) model by mathematically constraining spin density eigenvalues. The resulting constrained UHF (CUHF) wave function is identical to that obtained from Roothaan's effective Fock operator. Our $α$ and $β$ CUHF Fock operators are parameter-free and have canonical orbitals and orbital energies that are physically meaningful as in UHF, except for eliminating spin contamination. The present approach removes ambiguities in ROHF orbital energies and the non-uniqueness of methods that build upon them. We present benchmarks to demonstrate CUHF physical correctness and good agreement with experimental results.

preprint2009arXiv

Accurate Hartree-Fock energy of extended systems using large Gaussian basis sets

Calculating highly accurate thermochemical properties of condensed matter via wave function-based approaches (such as e.g. Hartree-Fock or hybrid functionals) has recently attracted much interest. We here present two strategies providing accurate Hartree-Fock energies for solid LiH in a large Gaussian basis set and applying periodic boundary conditions. The total energies were obtained using two different approaches, namely a supercell evaluation of Hartree-Fock exchange using a truncated Coulomb operator and an extrapolation toward the full-range Hartree-Fock limit of a Padé fit to a series of short-range screened Hartree-Fock calculations. These two techniques agreed to significant precision. We also present the Hartree-Fock cohesive energy of LiH (converged to within sub-meV) at the experimental equilibrium volume as well as the Hartree-Fock equilibrium lattice constant and bulk modulus.

preprint2009arXiv

Electromechanical properties of suspended Graphene Nanoribbons

Graphene nanoribbons present diverse electronic properties ranging from semiconducting to half-metallic, depending on their geometry, dimensions and chemical composition. Here we present a route to control these properties via externally applied mechanical deformations. Using state-of-the-art density functional theory calculations combined with classical elasticity theory considerations, we find a remarkable Young's modulus value of ~7 TPa for ultra-narrow graphene strips and a pronounced electromechanical response towards bending and torsional deformations. Given the current advances in the synthesis of nanoscale graphene derivatives, our predictions can be experimentally verified opening the way to the design and fabrication of miniature electromechanical sensors and devices based on ultra-narrow graphene nanoribbons.

preprint2009arXiv

Strong correlations via constrained-pairing mean-field theory

We present a mean-field approach for accurately describing strong correlations via electron number fluctuations and pairings constrained to an active space. Electron number conservation is broken and correct only on average but both spin and spatial symmetries are preserved. Optimized natural orbitals and occupations are determined by diagonalization of a mean-field Hamiltonian. This constrained-pairing mean-field theory (CPMFT) yields a two-particle density matrix ansatz that exclusively describe strong correlations. We demonstrate CPMFT accuracy with applications to the metal-insulator transition of large hydrogen clusters and molecular dissociation curves.

preprint2008arXiv

Analytical infrared intensities for periodic systems with local basis sets

We report a method for the efficient evaluation of analytic infrared (IR) intensities within generalized Kohn-Sham density functional theory using Gaussian orbitals and periodic boundary conditions. A discretized form of the Berry phase is used to evaluate a periodic dipole moment and its derivatives with respect to in-phase nuclear coordinate displacements. Benchmark calculations are presented for one-dimensional chains of water molecules and poly(paraphenylenevinylene).

preprint2007arXiv

Enhanced Half-Metallicity in Edge-Oxidized Zigzag Graphene Nanoribbons

We present a novel comprehensive first-principles theoretical study of the electronic properties and relative stabilities of edge-oxidized zigzag graphene nanoribbons. The oxidation schemes considered include hydroxyl, carboxyl, ether, and ketone groups. Using screened exchange density functional theory, we show that these oxidized ribbons are more stable than hydrogen-terminated nanoribbons except for the case of the etheric groups. The stable oxidized configurations maintain a spin-polarized ground state with antiferromagnetic ordering localized at the edges, similar to the fully hydrogenated counterparts. More important, edge oxidation is found to lower the onset electric field required to induce half-metallic behavior and extend the overall field range at which the systems remain half-metallic. Once the half-metallic state is reached, further increase of the external electric field intensity produces a rapid decrease in the spin magnetization up to a point where the magnetization is quenched completely. Finally, we find that oxygen containing edge groups have a minor effect on the energy difference between the antiferromagnetic ground state and the above-lying ferromagnetic state.

Gustavo E. Scuseria

What is connected

Connect this record

See the researcher in context

Building this map preview

70 published item(s)

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

A Power Series Approximation in Symmetry Projected Coupled Cluster Theory

Coupled cluster and perturbation theories based on a cluster mean-field reference applied to strongly correlated spin systems

Thermal coupled cluster theory for SU(2) systems

Exploring non-linear correlators on AGP

Correlating AGP on a quantum computer

Geminal replacement models based on AGP

On a dual representation of the Goldstone manifold

Exact Parameterization of Fermionic Wave Functions via Unitary Coupled Cluster Theory

Geminal-Based Configuration Interaction

Thermofield theory for finite-temperature coupled cluster

Thermofield Theory for Finite-Temperature Quantum Chemistry

On the difference between variational and unitary coupled cluster theories

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

Influence of broken-pair excitations on the exact pair wavefunction

On the Magnetic Structure of Density Matrices

Projected Coupled Cluster Theory

Tensor-Structured Coupled Cluster Theory

Biorthogonal projected energies of a Gutzwiller similarity transformed Hamiltonian

Performance of a Nonempirical Density Functional on Molecules and Hydrogen-Bonded Complexes

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

Predicting Band Gaps with Hybrid Density Functionals

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

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

A cluster-based mean-field and perturbative description of strongly correlated fermion systems. Application to the 1D and 2D Hubbard model

Actinide chemistry using singlet-paired coupled cluster and its combinations with density functionals

Can single-reference coupled cluster theory describe static correlation?

Combinations of coupled cluster, density functionals, and the random phase approximation for describing static and dynamic correlation, and van der Waals interactions

Pair extended coupled cluster doubles

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

Composite fermion-boson mapping for fermionic lattice models

Coupled Cluster Channels in the Homogeneous Electron Gas

Electron correlation in solids via density embedding theory

Lie Algebraic Similarity Transformed Hamiltonians for Lattice Model Systems

Multi-reference symmetry-projected variational approximation for the ground state of the doped one-dimensional Hubbard model

On the equivalence of LIST and DIIS methods for convergence acceleration

Polyradical character and spin frustration in fullerene molecules: An ab initio non-collinear Hartree--Fock study

Quasiparticle Coupled Cluster Theory for Pairing Interactions

Range Separated Brueckner Coupled Cluster Doubles Theory

Seniority-based coupled cluster theory

Symmetry-projected Wave Functions in Quantum Monte Carlo Calculations

The sign problem in full configuration interaction quantum Monte Carlo: Linear and sub-linear representation regimes for the exact wave function

Variational description of the ground state of the repulsive two-dimensional Hubbard model in terms of nonorthogonal symmetry-projected Slater determinants

Composite Boson Mapping for Lattice Boson Systems

Density Functionals that Recognize Covalent, Metallic, and Weak Bonds

Density Matrix Embedding from Broken Symmetry Lattice Mean-Fields

Excited electronic states from a variational approach based on symmetry-projected Hartree--Fock configurations

Multi-component symmetry-projected approach for molecular ground state correlations

Multi-reference symmetry-projected variational approaches for ground and excited states of the one-dimensional Hubbard model

On pair functions for strong correlations

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

Proper and improper zero energy modes in Hartree-Fock theory and their relevance for symmetry breaking and restoration

Renormalized Second-order Perturbation Theory for The Electron Correlation Energy: Concept, Implementation, and Benchmarks

Semilocal and Hybrid Meta-Generalized Gradient Approximations Based on the Understanding of the Kinetic-Energy-Density Dependence

The role of screened exact exchange in accurately describing properties of transition metal oxides: Modeling defects in LaAlO3

Assessment of correlation energies based on the random-phase approximation

N-electron Slater determinants from non-unitary canonical transformations of fermion operators

Projected Hartree-Fock Theory

Structural Phase Transitions of the Metal Oxide Perovskites SrTiO3, LaAlO3 and LaTiO3 Studied with a Screened Hybrid Functional

Symmetry-projected variational approach for ground and excited states of the two-dimensional Hubbard model

Accurate Modeling of the Cubic and Antiferrodistortive Phases of SrTiO3 with Screened Hybrid Density Functional Theory

Computing the energy of a water molecule using MultiDeterminants: A simple, efficient algorithm

Improved Semiconductor Lattice Parameters and Band Gaps from a Middle-Range Screened Hybrid Functional

Projected Quasiparticle Theory for Molecular Electronic Structure

ROHF Theory Made Simple

Accurate Hartree-Fock energy of extended systems using large Gaussian basis sets

Electromechanical properties of suspended Graphene Nanoribbons

Strong correlations via constrained-pairing mean-field theory

Analytical infrared intensities for periodic systems with local basis sets

Enhanced Half-Metallicity in Edge-Oxidized Zigzag Graphene Nanoribbons