Effective Hamiltonians derived from equation-of-motion coupled-cluster wave functions: Theory and application to the Hubbard and Heisenberg Hamiltonians
Effective Hamiltonians, which are commonly used for fitting experimental observables, provide a coarse-grained representation of exact many-electron states obtained in quantum chemistry calculations; however, the mapping between the two is not trivial. In this contribution, we apply Bloch's formalism to equation-of-motion coupled-cluster wave functions to rigorously derive effective Hamiltonians in Bloch's and des Cloizeaux's forms. We report the key equations and illustrate the theory by application to systems with two or three unpaired electrons, which give rise to electronic states of covalent and ionic characters. We show that Hubbard's and Heisenberg's Hamiltonians can be extracted directly from the so-obtained effective Hamiltonians. By establishing a quantitative connection between many-body states and simple models, the approach facilitates the analysis of the correlated wave functions. We propose a simple diagnostic for assessing the validity of the model space choice based on the overlaps between the target- and model-space states. Artifacts affecting the quality of electronic structure calculations such as spin contamination are also discussed.
Effective Hamiltonians, which are commonly used for fitting experimental observables, provide a coarse-grained representation of exact many-electron states obtained in quantum chemistry calculations; however, the mapping between the two is not trivial. In this contribution, we apply Bloch's formalism to equation-of-motion coupled-cluster wave functions to rigorously derive effective Hamiltonians in Bloch's and des Cloizeaux's forms. We report the key equations and illustrate the theory by application to systems with two or three unpaired electrons, which give rise to electronic states of covalent and ionic characters. We show that Hubbard's and Heisenberg's Hamiltonians can be extracted directly from the so-obtained effective Hamiltonians. By establishing a quantitative connection between many-body states and simple models, the approach facilitates the analysis of the correlated wave functions. We propose a simple diagnostic for assessing the validity of the model space choice based on the overlaps between the target- and model-space states. Artifacts affecting the quality of electronic structure calculations such as spin contamination are also discussed.
The problem of finding the spectrum of the Fock-space Hamiltonian \ifmmode \hat{H}\else \^{H}\fi{} for a system of many fermions is analyzed. The quasiparticle formalism is employed, with ``holes'' and ``particles'' defined with respect to some N-particle determinantal wave function (the model vacuum). The basic idea behind the proposed approach is to perform a similarity transformation of Hamiltonian \ifmmode \hat{H}\else \^{H}\fi{}, such that the resulting effective Hamiltonian \ifmmode \hat{G}\else \^{G}\fi{} is, unlike \ifmmode \hat{H}\else \^{H}\fi{}, a quasiparticle-number--conserving operator. It is shown that eigenvalues of \ifmmode \hat{G}\else \^{G}\fi{}, corresponding to small numbers of quasiparticles (0,1,2) can be easily calculated. This is equivalent to finding eigenvalues of \ifmmode \hat{H}\else \^{H}\fi{} for certain states of N, N\ifmmode\pm\else\textpm\fi{}1, and N\ifmmode\pm\else\textpm\fi{}2 particles. The construction of the operator transforming \ifmmode \hat{H}\else \^{H}\fi{} into \ifmmode \hat{G}\else \^{G}\fi{} (the wave operator) stems from an analysis of the structure of the algebra of linear operators acting in a (finite-dimensional) Fock space. The exponential Ansatz for the wave operator is used, resulting in a generalization of the coupled-cluster (CC) method of Coester [Nucl. Phys. 7, 421 (1958)]. The generalized CC equations determining the wave operator, and equations determining the effective Hamiltonian \ifmmode \hat{G}\else \^{G}\fi{}, are presented in a diagrammatic form. An effort has been made to obtain a concise notation for expressing these equations in an algebraic form. Approximation schemes, necessary for practical applications of the proposed method, are also studied.
In this Thesis, we develop and apply an accurate method to determine the ground-state properties of strongly-correlated electrons in the framework of lattice model Hamiltonians. In lattice density-functional theory (LDFT) the basic variable is the one-particle density matrix γ. From a generalized HK theorem, the ground state energy Egs[γgs] = minγ E[γ] is obtained by minimizing the energy-functional E[γ] over all the physical or representable γ. The energy functional can be divided into two contributions: the kinetic-energy functional T [γ], which linear dependence on γ is exactly known, and the correlation-energy functional W [γ], which dependence on γ is not explicitly known. Finding accurate approximations for W [γ] constitutes the actual challenge of this thesis. Part of this work is based on previous studies which derive an approximation of W [γ] for the Hubbard Hamiltonian based on scaling hypothesis and an exact analytical results for the Hubbard dimer. However, this approach is limited to spin independent and homogeneous systems. In order to extend the scope of LDFT, we have developed three different approaches to derive W [γ] in order to study broken symmetry problems. First we extend the previous scaling functional to system having charge transfer. By a systematic study of the dependence of W [γ] on the charge distribution we find similar scaling properties as for the homogeneous case. An extension to the Hubbard model on bipartite lattices is then derived and applied to finite and infinite systems with repulsive and attractive interactions. The high accuracy of this functional is demonstrated. Nevertheless, this approach is difficult to transfer to more complex systems, since the calculation of W [γ] applies to the system as a whole. To overcome this problem, we derive another approximation based on local scaling properties. This functional is formulated on-site so that it can be applied to any ordered or disordered Hamiltonian with on-site interaction. As applications we study the metal-insulator transition in the ionic Hubbard model in oneand two-dimensions, as well as the one-dimensional Hubbard chain with first and second nearest neighbor hoppings. Finally, we develop a numerical approach of W [γ], on the basis on exact diagonalizations of an effective many-body Hamiltonian corresponding to a cluster surrounded by an effective medium. This effective Hamiltonian depends on the density matrix γ and allows to derive approximation to W [γ] which improve systematically with increasing cluster size. The formulation is spin-dependent and allows a straightforward generalization to multi-orbital correlated systems such spd-Hamiltonians. In addition it take into account the effects of short-range charge and spin fluctuations in the functional. The accuracy of the method is demonstrated for the Hubbard model by comparison with the Bethe-Ansatz solution (1D) and quantum Monte Carlo simulations (2D). Finally, an outlook or relevant future developments of the theory is provided.
Model Hamiltonians (e.g., the Hubbard or Heisenberg models) provide a simple yet physically meaningful description of complex electronic phenomena through a small number of parameters, such as electron hopping integrals or magnetic exchange couplings. Their construction typically exploits the local nature of electron correlation effects and therefore relies on localized orbitals. Although many localization schemes are available, the resulting localized orbitals depend not only on the chosen localization functional but also on the orbital space to which the localization procedure is applied. These choices are often treated as technical details and are rarely discussed explicitly. Here, we demonstrate that different, chemically reasonable choices of localized orbitals can lead to substantially different model Hamiltonian parameters, even when they produce nearly identical electronic energies. Using Density Matrix Downfolding, we systematically investigate how the localization functional and, more importantly, the choice of model orbital space affect effective Hamiltonians derived from ab initio calculations for representative π-conjugated systems and transition metal complexes.
The standard and renormalized coupled cluster methods with singles, doubles, and noniterative triples and their generalizations to excited states, based on the equation of motion coupled cluster approach, are applied to the 4He and 16O nuclei. A comparison of coupled cluster results with the results of the exact diagonalization of the Hamiltonian in the same model space shows that the quantum chemistry inspired coupled cluster approximations provide an excellent description of ground and excited states of nuclei. The bulk of the correlation effects is obtained at the coupled cluster singles and doubles level. Triples, treated noniteratively, provide the virtually exact description.
A new approach to the diagrammatic formulation of many-body perturbation theory for open-shell systems is presented. The formalism is based on a generalized form of the Bloch equation that also generates the Rayleigh--Schroedinger perturbation expansion for a system with several open shells. Second quantization in the particle-hole formulation is used together with Wick's theorem in order to derive graphical rules in the usual way. The linked-diagram property of the wave operator and of the effective interaction is shown by expanding the wave operator in terms of normal products of connected diagram clusters, in analogy with the exp(S) formalism of Coester and Kuemmel for closed-shell systems. Self-consistent coupled-cluster equations are derived for the open-shell case from the generalized Bloch equation by a straightforward extension of the procedure of Cizek and Paldus. The application of such equations for investigating different properties of open-shell atoms is discussed. 56 references.
Coupled cluster (CC) theory is the “gold standard” of molecular electronic structure theory, yet the full arc of the subject—from its variational origins in Hartree–Fock theory, through the exponential ansatz and its diagrammatic machinery, to its modern reformulation as a system of polynomial equations on an algebraic variety, its multireference generalizations, and its emerging reincarnation as a quantum circuit—is rarely presented as a single connected narrative. This review builds that narrative completely from first principles and carries it through to today’s research frontier. Starting from the electronic Schrödinger equation, we derive Hartree–Fock theory, second quantization, and Wick’s theorem; construct the coupled cluster singles-and-doubles (CCSD) equations via the similarity-transformed Hamiltonian; and show how these same equations reappear, verbatim, as a truncated eigenvalue problem on a projective variety generalizing the Grassmannian, with a computable CC degree controlling the number of solutions. We then develop the perturbative triples correction CCSD(T), the linear-response equationof-motion (EOM-CC) framework for excited, ionized, and electron-attached states, the effective-Hamiltonian formalism underlying multireference coupled cluster theory and its notorious intruder-state problem, and the unitary coupled cluster (UCC) ansatz that restores the variational principle at the cost of an infinite operator expansion—tamed, on a quantum computer, by Trotterization and the variational quantum eigensolver. We close with a synthesis of what is rigorously known about the resulting nonlinear equations, a map of directions deliberately left out of scope, and an assessment of the open question of quantum advantage in this domain. Every claim is stated formally, with a full statement and proof or proof sketch, and cited to its origin; the computational cost of every method is compared explicitly throughout.