Effective Hamiltonians offer a realistic and practical approach to the representation of low-lying electronic excitations in strongly correlated charge transfer insulators, such as the high-Tc parent La2CuO4. Not only do such Hamiltonians minimize the required number of orbital degrees of freedom, but it is suggested that the fundamental problem of charged excitations in these materials may be folded into renormalized Hamiltonian parameters which already reflect solid state screening, and that such parameters may easily be obtained from constrained local density functional calculations. Beneficial consequences of the renormalization include reasonable insulating gaps already from the Hartree-Fock one-electron spectra, and a restriction of configuration interaction to those essential correlations associated with the spin degrees of freedom. In large periodic clusters where states have well characterized crystal momentum, such spin correlations in the vicinity of particles doped into the parent insulator are still reminiscent of the total spin eigenstates found in impurity model solutions.
Effective Hamiltonians offer a realistic and practical approach to the representation of low-lying electronic excitations in strongly correlated charge transfer insulators, such as the high-Tc parent La2CuO4. Not only do such Hamiltonians minimize the required number of orbital degrees of freedom, but it is suggested that the fundamental problem of charged excitations in these materials may be folded into renormalized Hamiltonian parameters which already reflect solid state screening, and that such parameters may easily be obtained from constrained local density functional calculations. Beneficial consequences of the renormalization include reasonable insulating gaps already from the Hartree-Fock one-electron spectra, and a restriction of configuration interaction to those essential correlations associated with the spin degrees of freedom. In large periodic clusters where states have well characterized crystal momentum, such spin correlations in the vicinity of particles doped into the parent insulator are still reminiscent of the total spin eigenstates found in impurity model solutions.
Special solutions of the Hartree-Fock (HF) problem for Coulomb interacting electrons, being described by a simple model of the Cu-O planes in La2CuO4, are presented. One of the mean field states obtained, is able to predict some of the basic properties of this material, such as its insulator character and the antiferromagnetic order. The natural appearance of pseudogaps in some states of this compound is also indicated by another of the HF states obtained. These surprising results follow after eliminating spin and crystal symmetry restrictions which are usually imposed on the single particle HF orbitals, by means of employing a rotational invariant formulation of the HF scheme which was originally introduced by Dirac. Therefore, it is exemplified how, up to now being considered strong correlation effects, can be described by improving the HF solution of the physical systems. In other words, defining the correlation effects as such ones shown by the physical system and which are not predicted by the best HF (lowest energy) solution, allows to explain currently assumed as strong correlation properties, as simple mean field ones. The discussion also helps to clarify the role of the antiferromagnetism and pseudogaps in the physics of the HTSC materials and indicates a promising way to start conciliating the Mott and Slater pictures for the description of the transition metal oxides and other strongly correlated electron systems.
The principal purpose of this topical review is to single out some of the universal features of high-temperature (high-Tc) and related strongly-correlated systems, which can be compared with experiment in a quantitative manner. The description starts with the concept of exchange-interaction-mediated (real-space) pairing, combined with strong correlations among narrow band electrons, for which the reference state is that of the Mott–Hubbard insulator. The physical discussion of concrete properties relies on variational approach, starting from generalized renormalized mean-field theory (RMFT) in the form of statistically-consistent Gutzwiller approximation (SGA), and its subsequent generalization, i.e., a systematic Diagrammatic Expansion of the Variational (Gutzwiller-type) Wave Function (DE-GWF). The solution leads to the two energy scales, one involving unconventional quasiparticles close to the Fermi energy, and the other reflecting the fully correlated state, involving electrons deeper below the Fermi surface. Those two regimes are separated by a kink in the dispersion relation, which is observed in photoemission. As a result, one obtains both the doping dependent properties in the correlated state, and effective renormalized quasiparticles. The reviewed ground-state characteristics for high-Tc systems encompass high-Tc superconductivity, nematicity, charge- (and pair-) density-wave effects, as well as non-BCS kinetic energy gain in the paired state. The calculated dynamic properties are: the universal Fermi velocity, the Fermi wave-vector, the effective mass enhancement, the pseudogap; and the d-wave gap magnitude. We discuss that, within the variational approach, the minimal realistic model is represented by the so-called t-J-U Hamiltonian. Inadequacy of the t-J and Hubbard models, particularly of their RMFT versions, is discussed explicitly. For heavy fermion systems, modeled by the Anderson lattice model and with the DE-GWF approach, we discuss the phase diagram, encompassing superconducting, Kondo insulating, ferro- and anti-ferromagnetic states. The superconducting state is then a two-d-wave gap system. If the orbital degeneracy of f-electrons is included, the coexistent ferromagnetic- (spin–triplet) superconducting phases appear and match those observed for UGe2 in a semiquantitative manner. Finally, in the second part, we generalize our approach to the collective spin and charge fluctuations in high-Tc systems, starting from variational approach, defining the saddle-point state, and combining it with 1/N expansion. The present scheme differs essentially from that starting from the saddle-point Hartree–Fock approximation, and incorporating the fluctuations in the random phase approximation (RPA). The spectrum of collective spin and charge excitations is determined for the Hubbard and t-J-U models, and subsequently compared quantitatively with recent experiments. The Appendices provide formal details to make this review self-contained.
Electronic structure of materials from first principles is essential for one-to-one comparisons between theories and experiments. However, the many-body effects resulting from electron correlation and their interplay with other degrees of freedom pose a great challenge to design an ab initio approximate method with good accuracy and efficiency at the same time. In this thesis, we present our recent developments of ab initio many-body theory for solids based on the finite-temperature Green’s function formalism. We begin by reviewing the basic concepts and definitions of realistic electronic Hamiltonians (Ch. 2) as well as Green’s functions at finite temperature (Ch. 3). Overviews of self-consistent many-body perturbation theories and quantum embedding theories are then presented in Chs. 4 and 5 as complementary materials for the following chapters. In Ch. 6, we present algorithms and implementation details for the finite-temperature self-consistent GW (scGW) method in Gaussian-type orbitals for solids. Our implementation is based on the finite-temperature Green’s function formalism in which all equations are solved on the imaginary axis, without resorting to analytical continuation during the self-consistency. No quasiparticle approximation is employed and all matrix elements of self-energies are explicitly evaluated. We show agreement with other, differently formulated finite-temperature scGW implementations. By migrating computationally intensive calculations to GPUs, we obtain scalable results on large supercomputers with nearly optimal performance. Inclusion of relativistic effects in scGW is discussed in Ch. 7. Our formulation of relativistic scGW is based on the exact two-component formalism with one-electron approximation (X2C1e) and non-relativistic Coulomb interactions. Our theory allows us to study scalar relativistic effects, spin-orbit coupling, and their interplay with electron correlation without adjustable parameters. We examine the effect of the X2C1e by comparison to the established four-component formalism and reach excellent agreement. The simplicity of X2C1e enables construction of higher order theories, such as embedding theories, on top of perturbative calculations. In Ch. 8, we use self-energy embedding theory (SEET) to study the spectra of the prototypical undistorted cubic perovskites SrVO3 and SrMnO3. In the strongly correlated metallic SrVO3, we find that the usual attribution of the satellite peaks at -1.8 eV to Hund or Hubbard physics in the t2g orbitals is inconsistent with our calculations. In the strongly correlated insulator SrMnO3 we recover insulating behavior due to a feedback effect between the strongly correlated orbitals and the weakly correlated environment. Our calculation shows a systematic convergence of spectral features as the space of strongly correlated orbitals is enlarged, paving the way to a systematic parameter-free study of correlated perovskites. Lastly, in Ch. 9, we apply the Green’s function coupled cluster singles and doubles (GFCCSD) impurity solver to realistic impurity problems in strongly correlated solids within the framework of SEET. We describe the implementation details of our Green’s function coupled cluster (GFCC) solver, investigate its performance, and highlight potential advantages and problems on examples of impurities created for antiferromagnetic MnO and paramagnetic SrMnO3. GFCCSD provides satisfactory descriptions for weakly and moderately correlated impurities with sizes that are intractable by existing accurate impurity solvers such as exact diagonalization. However, when correlations become strong, the singles and doubles approximation used in GFCC could lead to instabilities in searching for the particle number present in impurity problems.
Strong correlations arising from interactions between particles are responsible for many phenomena in condensed matter physics, from superconductivity to heavy fermions. An accurate theoretical description of such complex systems requires to solve so-called many-body problems, which are inherently difficult to tackle due to the exponential growth of the underlying Fock space with the size of the system. Usually, single-particle descriptions or perturbation theory at first orders fail to provide satisfactory results in correlated regimes, and pertinent Ansätze or cutting edge numerical methods are required to solvethese problems.Among the wide variety of strongly correlated electronic systems, quantum impurity systems are of central interest. These systems, where a few localized degrees of freedom that undergo strong interactions hybridize with a larger system of free particles, have been widely studied for decades. Beyond their relevance to study dilute dynamic impurities in metals, quantum dots, or more complex nanostructures, they are exploited in a broader class of problems through dynamical mean field theory (DMFT), that is designed for fully interacting lattice problems. Despite the effort devoted to studying of these problems,some questions remain open and new approaches are required to address them.In this thesis, quantum impurity problems are tackled with the use of natural orbitals, a representation rooted in the framework of quantum chemistry. The numerical renormalization group, the standard non-perturbative method to treat these systems, is employed to demonstrate the hierarchical structure of correlations unveiled by natural orbitals. An ansatz for the wavefunction of the ground state is then proposed, separating explicitely a small set of correlated orbitals to a large uncorrelated one, leaving the problem to a few-body complexity. This simplification brought new insights for simulating quantumimpurity problems, and a new algorithm is developed to generate an optimized set of correlated orbitals independently of existing methods, going beyond their usual limitations.In the first place, the algorithm is used to compute the screening cloud in large one-dimensional chains and in a realistic 2D square lattice coupled to a localized adatom. This observable is of great interest to understand the screening mechanism of the magnetic moment of the impurity that plays a central role in transport properties, for unidimensional systems or in bulk metals experiencing RKKY interactions. Spatial correlations were found to be highly non-symmetric, breaking the s-wave picture. In one direction, corre-lations vanish way faster than the usual isotropic correlation length, while on the other one they spread further than this same expected length. This anisotropy suggests that impurities in bulk systems may interact differently depending on their relative positions with respect to the lattice.Taking further advantage of the efficiency of this new algorithm allowed us to perform a statistical study on large systems to investigate correlations around a quantum impurity in a disordered host, averaging over a large amount of realizations, which was not possible with previous approaches. The question motivating this study was the fate of the correlation cloud spreading into the disordered electronic bath. Although spatial correlations are expected to vanish when the disorder increases, long range correlationssurvive to some configurations of the potential, inducing a bimodal structure of probability distribution functions of relevant scales.Finally, these disordered problems are investigated in the non-interacting situation with random matrix theory, which reproduced qualitatively the results obtained before. The joint distribution of one-particle orbitals energies and amplitudes is calculated, which allows to compute any disordered averaged local correlation functions, and calculations in the large-N limit are presented.
Mott insulators, in which strong Coulomb interactions fully localize electrons on single atomic sites, play host to an incredibly rich and exciting array of strongly correlated physics. One can naturally extend this concept to cluster Mott insulators, wherein electrons localize not on single atoms but across clusters of atoms, forming "molecules in solids". The resulting localized degrees of freedom incorporate the full spectrum of electronic degrees of freedom, spin, orbital, and charge. These serve as the building blocks for cluster Mott insulators, and understanding them is an important first step toward understanding the many-body physics that emerges in candidate cluster Mott insulators. Here, we focus on elementary building blocks, neglecting some of the complexity present in real materials which can often obfuscate the underlying principles at play. Through an extensive set of exact theoretical calculations on clusters of varying geometry, number of orbitals, and number of electrons, we uncover some of the basic organizing principles of cluster Mott phases, particularly when interactions dominate and negate a simple single-particle picture. These include, for example, the identification of an additional "cluster Hund's rule", of cluster ground states that are best understood from a purely interacting perspective, and of several localized degrees of freedom which are protected by an unusual combination of discrete spatial or orbital symmetries. Finally, we discuss the impact of adding additional terms, relevant to material candidates, on the phase diagrams presented throughout, as well as the potential next steps in the journey to building a more complete picture of cluster Mott insulators.
In the previous chapters DMFT has been developed as a methodology for calculation of dynamical and thermodynamic properties of many-body systems, with applications to model systems. In order for DMFT to be a general method for quantitative calculations, there must be systematic procedures to derive all the ingredients from first principles. In this chapter we consider the most prominent examples, materials with partially filled 3 d , 4 f , or 5 f localized, atomic-like states, and set up the characteristic form of the hamiltonian. This is the basis for DFT+DMFT methods discussed here and used in the examples in Ch. 20, and for the combination of DMFT with many-body perturbation theory in Ch. 21. Simpler techniques, including self-interaction correction, DFT+U, and DFT+Gutzwiller, are also useful and are summarized in the last sections of this chapter. Up to this point, dynamical mean-field theory has been developed for model systems with Hubbard-type hamiltonians, where the hopping matrix elements and interactions are considered as parameters. For problems with strong interactions, DMFT leads to renormalized bands, satellites in the one-particle spectra, magnetic phase transitions and local moments, metal–insulator transitions, and many other phenomena observed in real materials, as brought out in Ch. 2. The purpose of this chapter is to provide the framework for a quantitative theory of these materials. These are difficult problems, experimentally, theoretically, and computationally, and the first step is to identify the aspects of the materials that must be treated explicitly as the local strongly interacting system, and the “rest” that are treated by a different method. This division of the problem is useful because a real material involves many degrees of freedom, i.e., the electronic states must be described in a basis sufficiently large to describe all the relevant states. However, the part that can be treated by non-perturbative solvers (Ch. 18) is limited by the computational cost that scales rapidly with the number of states included in the embedded site or cluster. Thus one must choose a small number of localized orbitals that inevitably are non-unique. This may seem like a step backward for quantitative methods; however, it opens a window for methods to treat these interesting problems and it presents a challenge to overcome the limitations.