While calculations and measurements of single-particle spectral properties often offer the most direct route to study correlated electron systems, the underlying physics may remain quite elusive, if information at higher particle levels is not explicitly included. Here, we present a comprehensive overview of the different approaches which have been recently developed and applied to identify the dominant two-particle scattering processes controlling the shape of the one-particle spectral functions and, in some cases, of the physical response of the system. In particular, we will discuss the underlying general idea, the common threads and the specific peculiarities of all the proposed approaches. While all of them rely on a selective analysis of the Schwinger-Dyson (or the Bethe-Salpeter) equation, the methodological differences originate from the specific two-particle vertex functions to be computed and decomposed. Finally, we illustrate the potential strength of these methodologies by means of their applications the two-dimensional Hubbard model, and we provide an outlook over the future perspective and developments of this route for understanding the physics of correlated electrons.
While calculations and measurements of single-particle spectral properties often offer the most direct route to study correlated electron systems, the underlying physics may remain quite elusive, if information at higher particle levels is not explicitly included. Here, we present a comprehensive overview of the different approaches which have been recently developed and applied to identify the dominant two-particle scattering processes controlling the shape of the one-particle spectral functions and, in some cases, of the physical response of the system. In particular, we will discuss the underlying general idea, the common threads and the specific peculiarities of all the proposed approaches. While all of them rely on a selective analysis of the Schwinger-Dyson (or the Bethe-Salpeter) equation, the methodological differences originate from the specific two-particle vertex functions to be computed and decomposed. Finally, we illustrate the potential strength of these methodologies by means of their applications the two-dimensional Hubbard model, and we provide an outlook over the future perspective and developments of this route for understanding the physics of correlated electrons.
In this thesis, we use the dynamical cluster approximation to study strongly correlated electron systems, especially from the angle of two-particle quantities, such as dynamical susceptibilities and vertex functions. The thesis starts with an introduction to the strongly correlated systems, including their definitions, prominent features, applications, difficulties in explaining them theoretically and some numerical approaches developed. The following section is an introduction to one family of strongly correlated systems we focus on in this thesis, the high temperature cuprates. The salient features in their general phase diagrams are described and discussed. Then an overview of the model we used to study high Tc cuprates is provided, with its limitations and extensions. To solve this model, the numerical method we employ for our study is the dynamical mean-field theory, the dynamical cluster approximation and the continuous time auxiliary field impurity solver. the last part of Chap.1 contains brief derivations for these algorithms. In Chap. 2, we apply dynamical cluster approximation to solve the one-band 2D Hub- bard model. The physical quantities of interest are two-particle quantitites, such as the dynamical susceptibility, irreducible vertex functions and full vertex functions. In this chapter, we describe how to obtain these susceptibilities via linear response theory and write down a detailed example for superconducting susceptibility. Finally we show how to calculate two particle quantities within DCA and obtain phase boundary with them. Chap. 3 is based on one of our publications. In this work, we specifically address the problem of optimizing the superconducting transition temperature in the 2D Hubbard model by analyzing wide regions of parameter space in density, interaction, and second- nearest-neighbor hopping strength.We mainly focus on d-wave superconductivity but show results of other symmetries in the last section. Chap. 4 follows another publication of ours. We study the temperature and doping evolution of NMR response in the normal state of the 2D Hubbard model using cluster dynamical mean-field theory. We simulate the Knight shift, the spin-echo decay rate and the spin-lattice relaxation time, and compare them to the cuprates experimental results. The last chapter extends the calculation of the NMR response to the superconducting state with the Nambu formalism. We show the detailed formulas and diagrams to calculate two-particle quantities, including the Dyson-Schwinger equation of motion.
We demonstrate how to identify which physical processes dominate the low-energy spectral functions of correlated electron systems. We obtain an unambiguous classification through an analysis of the equation of motion for the electron self-energy in its charge, spin, and particle-particle representations. Our procedure is then employed to clarify the controversial physics responsible for the appearance of the pseudogap in correlated systems. We illustrate our method by examining the attractive and repulsive Hubbard model in two dimensions. In the latter, spin fluctuations are identified as the origin of the pseudogap, and we also explain why d-wave pairing fluctuations play a marginal role in suppressing the low-energy spectral weight, independent of their actual strength.
Identifying the fingerprints of the Mott-Hubbard metal-insulator transition may be quite elusive in correlated metallic systems if the analysis is limited to the single particle level. However, our dynamical mean-field calculations demonstrate that the situation changes completely if the frequency dependence of the two-particle vertex functions is considered: The first nonperturbative precursors of the Mott physics are unambiguously identified well inside the metallic regime by the divergence of the local Bethe-Salpeter equation in the charge channel. In the low-temperature limit this occurs for interaction values where incoherent high-energy features emerge in the spectral function, while at high temperatures it is traceable up to the atomic limit.
This thesis contains a series of numerical studies of strongly correlated electron systems. In these systems, interesting emergent properties are brought about by strong electron correlations, leading to unusual fluctuations and phase transitions. We start with a review of basic concepts in many-body physics from a field theory point of view, including second quantization and path integral formalisms. We then introduce extensions of the single-particle and two-particle Green’s function formalisms in the singlet superconducting state, which enable us to extract physical information of the systems in the symmetry-broken phase. Next, we proceed to the models and numerical methods we use to study the strongly correlated systems. We introduce a low energy effective model -- the Hubbard model, which contains a subset of the electron degrees of freedom and can be solved numerically using advanced many-body methods. With the dynamical cluster approximation and continuous-time quantum Monte Carlo impurity solver, we are able to study the competing fluctuations in the paramagnetic state and analyze the fluctuations behind superconductivity in the singlet superconducting state. We introduce self-consistent diagrammatic methods used in realistic material calculations in the last part of this thesis, with an outline of group theory concepts that can be used to optimize the simulations. Numerical representations and methods for effectively solving equations in realistic material calculations are discussed afterwards. We briefly review some of the developments in this field, and then introduce spectral methods that are based on mathematical properties of Legendre polynomials for solving both the imaginary- and real-time Dyson equations.
In this work we investigate some aspects of the physics of strongly correlated systems by taking into account both electron-electron and electron-phonon interactions as basic mechanisms for reproducing electronic correlations in real materials. The relevance of the electron-electron interactions is discussed in the first part of this thesis in the framework of a self-consistent theoretical approach, named Composite Operator Method (COM), which accounts for the relevant quasi-particle excitations in terms of a set of composite operators that appear as a result of the modification imposed by the interactions on the canonical electronic fields. We show that the COM allows the calculation of all the relevant Green’s and correlation functions in terms of a number of unknown internal parameters to be determined self-consistently. Therefore, depending on the balance between unknown parameters and self-consistent equations, exact and approximate solutions can be obtained. By way of example, we discuss the application of the COM to the extended t-U- J-h model in the atomic limit, and to the two-dimensional single-band Hubbard model. In the former case, we show that the COM provides the exact solution of the model in one dimension. We study the effects of electronic correlations as responsible for the formation of a plethora of different charge and/or spin orderings. We report the phase diagram of the model, as well as a detailed analysis of both zero and finite temperature single-particle and thermodynamic properties. As far as the single-band Hubbard model is concerned, we illustrate an approximated selfconsistent scheme based on the choice of a two-field basis. We report a detailed analysis of many unconventional features that arise in single-particle properties, thermodynamics and system’s response functions. We emphasize that the accuracy of the COM in describing the effects of electronic correlations strongly relies on the choice of the basis, paving the way for possible multi-pole extensions to the twofield theory. To this purpose, we also study a three-field approach to the single-band Hubbard model, showing a significant step forward in the agreements with numerical data with respect to the two-pole results. The role of the electron-phonon interaction in the physics of strongly correlated systems is discussed in the second part of this thesis. We show that in highly polarizable lattices the competition between unscreened Coulomb and Fröhlich interactions results in a short-range polaronic exchange term Jp that favours the formation of local and light pairs of bosonic nature, named bipolarons, which condense with a critical temperature well in excess of hundred kelvins. These findings, discussed in the framework of the so-called polaronic t-Jp model, are further investigated in the presence of a finite on-site potential ~U , coming from the competition between on-site Coulomb and Fröhlich interactions. We discuss the role of ~U as the driving parameter for a small-to-large bipolaron transition, providing a possible explanation of the BEC-BCS crossover in terms of the properties of the bipolaronic ground state. Finally, we show that a hard-core bipolarons gas, studied as a charged Bose-Fermi mixture, allows for the description of many non Fermi liquid behaviours, allowing also for a microscopic explanation of pseudogap features in terms of a thermal-induced recombination of polarons and bipolarons, without any assumption on preexisting order or broken symmetries. [edited by Author]
In this thesis, we tackle various problems in strongly correlated electron systems, which can be addressed properly via the non-perturbative dynamical mean field theory (DMFT) approach using the continuous time quantum Monte Carlo method as an impurity solver. First, we revisit the old Nagaoka ferromagnetism problem in the U=oo Hubbard model and study the stability of the ferromagnetic state as a function of the temperature, the doping level, and the next-nearest-neighbor hopping t'. We then address the nature of the Mott transition in the two-dimensional Hubbard model at half-filling using cluster DMFT. Cluster DMFT can incorporate the short-range correlations beyond DMFT by extending the spatial range in which correlations are treated exactly to a finite cluster size. The non-local correlations reduce substantially the critical interaction U and modify the shape of the transition lines in the phase diagram. We then concentrate on the calculation of two-particle response functions from the ab initio perspective by means of computing the one-particle excitation spectrum using the combination of the density functional theory (DFT) and DMFT and extracting the two-particle irreducible vertex function from a local two-particle Green's function computed within DMFT. In particular, we derive the equations for calculating the magnetic/charge susceptibility and the pairing susceptibility in superconductivity. This approach is applied to the Hubbard model and the periodic Anderson model and we determine the phase diagram of magnetism and superconductivity in these models. We show that the superconducting phase is indeed stable near the magnetic phase where the pairing interaction mediated by spin fluctuations is dominantly enhanced. The non-local correlation effect to superconductivity is also discussed using the dual fermion approach and the dynamical vertex approximation. We finally apply the vertex function approach within DFT+DMFT to a Fe-based superconductor, BaFe2As2, and compute the dynamical magnetic susceptibility in this material. Our calculation results show a good agreement with the magnetic excitation spectra observed in a neutron scattering experiment. The response function calculation method derived in this thesis can capture both a localized and an itinerant nature of collective excitations in strongly correlated electron systems.