We present the continued fraction method (CFM) as a microscopic approximation to the spectral density of the Hubbard model in the correlated metal phase away from half filling. The quantity expanded as a continued fraction is the single-particle Green's function. Leading spectral moments are taken into account through a set of real expansion coefficients, as known from the projection technique. Further stages are added to the continued fraction, with complex coefficients, thus defining a terminator function. This enables us to treat the entire spectral range of the Green's function on equal footing and determine the energy scale of the Fermi liquid quasiparticles by minimizing the total energy. The solution is free of phenomenological parameters and remains well defined in the strong-coupling limit, near the doping-controlled metal-insulator transition. Our results for the density of states agree reasonably with those from several variants of the dynamical mean-field theory. The CFM requires minimal numerical effort and can be generalized in several ways that are interesting for applications to real materials.
We present the continued fraction method (CFM) as a microscopic approximation to the spectral density of the Hubbard model in the correlated metal phase away from half filling. The quantity expanded as a continued fraction is the single-particle Green's function. Leading spectral moments are taken into account through a set of real expansion coefficients, as known from the projection technique. Further stages are added to the continued fraction, with complex coefficients, thus defining a terminator function. This enables us to treat the entire spectral range of the Green's function on equal footing and determine the energy scale of the Fermi liquid quasiparticles by minimizing the total energy. The solution is free of phenomenological parameters and remains well defined in the strong-coupling limit, near the doping-controlled metal-insulator transition. Our results for the density of states agree reasonably with those from several variants of the dynamical mean-field theory. The CFM requires minimal numerical effort and can be generalized in several ways that are interesting for applications to real materials.
A new non-perturbative framework for many-body correlated systems is formulated by extending the operator projection method (OPM). This method offers a systematic expansion which enables us to project into the low-energy structure after extracting the higher-energy hierarchy. This method also opens a way to systematically take into account the effects of collective excitations. The Mott-Hubbard metal-insulator transition in the Hubbard model is studied by means of this projection beyond the second order by taking into account magnetic and charge fluctuations in the presence of the high-energy Mott-Hubbard structure. At half filling, the Mott-Hubbard gap is correctly eproduced between the separated two bands. Near half filling, a strongly renormalized low-energy single-particle excitations coexisting with the Mott-Hubbard bands are shown to appear. Signifcance of momentum-dependent self-energy in the results is stressed.
We implement a highly efficient strong-coupling expansion for the Green's function of the Hubbard model. In the limit of extreme correlations, where the onsite interaction is infinite, the evaluation of diagrams simplifies dramatically enabling us to carry out the expansion to the eighth order in powers of the hopping amplitude. We compute the finite-temperature Green's function analytically in the momentum and Matsubara frequency space as a function of the electron density. Employing Padé approximations, we study the equation of state, Kelvin thermopower, momentum distribution function, quasiparticle fraction, and quasiparticle lifetime of the system at temperatures lower than, or of the order of, the hopping amplitude. We also discuss several different approaches for obtaining the spectral functions through analytic continuation of the imaginary frequency Green's function, and show results for the system near half filling. We benchmark our results for the equation of state against those obtained from a numerical linked-cluster expansion carried out to the eleventh order.
We introduce an alternative approach for obtaining smooth finite-temperature spectral functions of quantum impurity models using the numerical renormalization group (NRG) technique. It is based on calculating first the Green's function on the imaginary-frequency axis, followed by an analytic continuation to the real-frequency axis using Pad\'e approximants. The arbitrariness in choosing a suitable kernel in the conventional broadening approach is thereby removed and, furthermore, we find that the Pad\'e method is able to resolve fine details in spectral functions with less artifacts on the scale of $\ensuremath{\omega}\ensuremath{\sim}T$. We discuss the convergence properties with respect to the NRG calculation parameters (discretization $\ensuremath{\Lambda}$, z-averaging, truncation cutoff) and the number of Matsubara points taken into account in the analytic continuation. We test the technique on the single-impurity Anderson model and the Hubbard model (within the dynamical mean-field theory). For the Anderson impurity model, we discuss the shape of the Kondo resonance and its temperature dependence. For the Hubbard model, we discuss the inner structure of the Hubbard bands in metallic and insulating solutions at half-filling, as well as in the doped Mott insulator. Based on these test cases, we conclude that the Pad\'e approximant approach provides improved results for spectral functions at low-frequency scales of $\ensuremath{\omega}\ensuremath{\lesssim}T$ and that it is capable of resolving sharp spectral features also at high frequencies.
Implementing an improved method for analytic continuation and working with imaginary-time correlation functions computed using quantum Monte Carlo simulations, we resolve the single-particle dispersion relation and the density of states (DOS) of the two-dimensional Hubbard model at half filling. At intermediate interactions of U/t=4,6, we find quadratic dispersion around the gap minimum at wave vectors $\mathbf{k}=(\ifmmode\pm\else\textpm\fi{}\ensuremath{\pi}/2,\ifmmode\pm\else\textpm\fi{}\ensuremath{\pi}/2)$ (the $\mathrm{\ensuremath{\Sigma}}$ points). We find saddle points at $\mathbf{k}=(\ifmmode\pm\else\textpm\fi{}\ensuremath{\pi},0),(0,\ifmmode\pm\else\textpm\fi{}\ensuremath{\pi})$ (the X points), where the dispersion is approximately quartic, leading to a sharp DOS maximum above the almost flat ledge arising from the states close to $\mathrm{\ensuremath{\Sigma}}$. The fraction of quasiparticle states within the ledge is ${n}_{\mathrm{ledge}}\ensuremath{\approx}0.15$. Upon doping away from half filling, within the rigid-band approximation, these results support Fermi pockets around the $\mathrm{\ensuremath{\Sigma}}$ points, with states around the X points becoming filled only at doping fractions $x\ensuremath{\ge}{n}_{\mathrm{ledge}}$. The high density of states away from the $\mathrm{\ensuremath{\Sigma}}$ gap edge may be an important clue for a finite minimum doping level for superconductivity and other instabilities of doped Mott insulators.
A dynamical mean-field theory analysis of the attractive Hubbard model in the normal phase is carried out upon restricting to solutions where superconducting order is not allowed. A clear first-order pairing transition as a function of the coupling takes place at all the electron densities out of half filling between a Fermi liquid, stable for U<U(c), and an insulating bound pairs phase for U>U(c), and it is accompanied by phase separation. The spectral function in the metallic phase is constituted by a low-energy structure around the Fermi level, which disappears discontinuously at U = U(c), and two high-energy features (Hubbard bands), which persist in the insulating phase.
The local moment approach (LMA) has presented itself as a powerful semianalytical quantum impurity solver (QIS) in the context of the dynamical mean-field theory (DMFT) for the periodic Anderson model and it correctly captures the low-energy Kondo scale for the single impurity model, having excellent agreement with the Bethe ansatz and numerical renormalization group (NRG) results. However, the most common correlated lattice model, the Hubbard model, has not been explored well within the LMA+DMFT framework beyond the insulating phase. Here in our work, within the framework we complete the filling-interaction phase diagram of the single band Hubbard model at zero temperature. Our formalism is generic to any particle filling and can be extended to finite temperature. We contrast our results with another QIS, namely the iterated perturbation theory (IPT) and show that the second spectral moment sum rule improves better as the Hubbard interaction strength grows stronger in LMA, whereas it severely breaks down after the Mott transition in IPT. For the metallic case, the Fermi liquid (FL) scaling agreement with the NRG spectral density supports the fact that the FL scale emerges from the inherent Kondo physics of the impurity model. We also show that, in the metallic phase, the FL scaling of the spectral density leads to universality which extends to infinite frequency range at infinite correlation strength (strong coupling). At large interaction strength, the off half-filling spectral density forms a pseudogap near the Fermi level and filling-controlled Mott transition occurs as one approaches the half-filling. As a response property, we finally study the zero temperature optical conductivity and find universal features such as absorption peak position governed by the FL scale and a doping independent crossing point, often dubbed the isosbestic point in experiments.