The two-dimensional Hubbard model with a single spin-up electron interacting with a finite density of spin-down electrons is studied by use of the quantum Monte Carlo technique, the standard second-order perturbation theory and a very good variational wave function introduced recently by Edwards. We also present a new conjugate gradient scheme for the evaluation of this variational state. We performed simulations up to 242 sites at U/t=4 reaching the zero-temperature properties with no ``fermion-sign problem'' and found a surprisingly good accuracy of the Edwards state at low density or low doping. The conjugate gradient method was then applied to the system up to 3698 sites and infinite U. Fermi-liquid theory seems to remain stable in two dimensions for all cases studied with the exception of the half-filling case where a Luttinger-like behavior survives in the Hubbard model, yielding a vanishing quasiparticle weight in the thermodynamic limit.
The two-dimensional Hubbard model with a single spin-up electron interacting with a finite density of spin-down electrons is studied by use of the quantum Monte Carlo technique, the standard second-order perturbation theory and a very good variational wave function introduced recently by Edwards. We also present a new conjugate gradient scheme for the evaluation of this variational state. We performed simulations up to 242 sites at U/t=4 reaching the zero-temperature properties with no ``fermion-sign problem'' and found a surprisingly good accuracy of the Edwards state at low density or low doping. The conjugate gradient method was then applied to the system up to 3698 sites and infinite U. Fermi-liquid theory seems to remain stable in two dimensions for all cases studied with the exception of the half-filling case where a Luttinger-like behavior survives in the Hubbard model, yielding a vanishing quasiparticle weight in the thermodynamic limit.
An unbiased zero-temperature auxiliary-field quantum Monte Carlo method is employed to analyze the nature of the semimetallic phase of the two-dimensional Hubbard model on the honeycomb lattice at half filling. It is shown that the quasiparticle weight Z of the massless Dirac fermions at the Fermi level, which characterizes the coherence of zero-energy single-particle excitations, can be evaluated in terms of the long-distance equal-time single-particle Green's function. If this quantity remains finite in the thermodynamic limit, the low-energy single-particle excitations of the correlated semimetallic phase are described by a Fermi-liquid-type single-particle Green's function. Based on the unprecedentedly large-scale numerical simulations on finite-size clusters containing more than 10 000 sites, we show that the quasiparticle weight remains finite in the semimetallic phase below a critical interaction strength. This is also supported by the long-distance algebraic behavior ($\ensuremath{\sim}{r}^{\ensuremath{-}2}$, where r is distance) of the equal-time single-particle Green's function that is expected for the Fermi liquid. Our result thus provides a numerical confirmation of Fermi-liquid theory in two-dimensional correlated metals.
We study the Hubbard model in dimension d= infinity , within an approximation that simultaneously reproduces the atomic limit and the weak-coupling limit up to second order in the interaction strength U. This method, which is a slightly modified version of an approximation suggested by Edwards and Hertz, is applied to the calculation of the self-energy and of the spectral function in the paramagnetic regime for different U, band fillings n, and temperatures T. We find Fermi-liquid behaviour for small U and (for half-filling n=0.5) a metal-insulator (Mott) transition at a critical value U cr , but also a metallic non-Fermi-liquid phase for intermediate values U<U cr . We compare our results with results obtained within the second-order perturbation theory in U and (at finite temperature) with the (essentially exact) quantum Monte Carlo (QMC) results. Qualitatively, the approximation reproduces the correct tendencies obtained by QMC, but in detail there are still discrepancies, in particular for intermediate values of U.
The Fermi liquid theory is a central concept in modern condensed matter physics used to describe conventional metals. A state of this universality class consists of well-defined quasiparticles, which occupy a finite number of states within a Fermi volume defined in momentum space. According to the Luttinger’s theorem, such a state encloses a Fermi volume proportional to the electron density modulo the filled bands if no symmetry is broken. In the past decades strong deviations from Fermi liquid theory have been observed in the scaling behavior of thermodynamic observables in different materials, among which are for instance cuprates and iron pnictides. This thesis is about two prototypical systems that can not be described by Fermi liquid theory. The first part of this work investigates two-dimensional effective models, which are in a Mott insulating state at half-filling and have a finite conductivity when doped with holes. In such a Mott insulating state, the charge carriers are strongly localized due to the Coulomb repulsion. Hence, strong correlations are assumed to have a strong impact on the formation of the ground state, also in the regime of small hole concentration. In a first project we have used a so-called spinon dopon mean-field theory to represent the two-dimensional Fermi-Hubbard model with strong on-site repulsive interaction in effective degrees of freedom, in which the holes can be embedded into a quantum spin liquid. The corresponding SU(2) invariant ground state belongs to the class of fractionalized Fermi liquids. In a second project we investigate a quantum dimer model, an effective model based on a Hilbert space spanned by short range singlets and bound states of holes and spins. The focus here is on the calculation of the hole-part of the electron spectral function by using exact diagonalization and its comparison with two analytic methods, a diagrammatic computation based on the Bethe-Salpeter equation and a so-called two-mode approximation. The electron spectral function shows a similar analytic form in momentum space between nodal and antinodal point when compared to results from photoemission spectroscopy experiments on cuprates. Furthermore, in a subsequent work we calculate the exact ground state of the quantum dimer model along a certain parameter line. In order to analyze the behavior of the elctron spectral function when increasing the density of holes, we investigate the Fermi-Hubbard model in a current project using a dynamical mean-field approach. To solve the 4-site cluster impurity problem, we use a numerical renormalization group approach. However, numerical limitations force us to restrict the analysis to spin-polarized baths. According to the spectral data, the system is similar to the SU(2) invariant case at half-filling in a Mott insulating state and posseses a momentum-selective energy gap at finite doping. The topology of the Fermi surface shows a Lifshitz transition when increasing the hole concentration. Here, the curvature of the Fermi surface changes from electron- to hole-like. For comparison we apply the dynamical mean-field theory also to the quantum dimer model and observe that the electron spectral functions at finite doping are qualitatively similar to that of the two-dimensional Fermi-Hubbard model. The second half of the work is about Tomonaga-Luttinger liquid theory, which is used to describe the low-energy effective degrees of freedom of one-dimensional systems. Here, we first discuss a conceptional extension of the operator-based bosonization theory for one-dimensional systems. This extension is especially suited for inhomogeneous one-dimensional systems. First, we investigate a one-dimensional system with a local interation potential and compute an exact solution of the single particle propagator at T = 0. The critical exponent of the single particle propagator has an unconventional form as a function of the microscopic Tomonaga-Luttinger parameters, which is not covered by the original Luttinger paradigm postulated by F. Duncan M. Haldane. In a second project on one-dimensional systems, we study the impact of scattering processes among bosonic low-energy excitations on the thermalization process. Such scattering processes are irrelevant on large length scales, however strongly affect the dynamics. In our analytic analysis we focus on a experimental setup, where a one-dimensional Bose gas is instantaneously splitted in two identical, however strongly correlated, halves of one-dimensional electronic systems. In the following, the corresponding non-equilibrium state runs through multiple regimes in time, such as a metastable prethermalization regime. However, above a certain threshold in time such scattering processes cause an effective thermalization of the system. In order to demonstrate this, we compute the kinetic equation in the Keldysh field integral formalism from a diagrammatic expansion based on a self-consistent Born approximation.
A two-dimensional (2D) spin-up ideal Fermi gas interacting attractively with a spin-down impurity in the continuum undergoes, at zero temperature, a first-order phase transition from a polaron to a dimeron state. Here, we study a similar system on a square lattice, by considering the attractive 2D Fermi-Hubbard model with a single spin-down and a finite filling fraction of spin-up fermions. We study polaron and dimeron quasiparticle properties via variational ansatz up to one particle-hole excitation. Moreover, we develop a determinant diagrammatic Monte Carlo algorithm for this problem based on expansion in bare on-site coupling U . This algorithm turns out to be sign-problem-free at any filling of spin-up fermions, allowing one to sample very high diagram order (larger than 200 in our study) and to do simulations for large U / t (we go up to U / t = − 20 with t the hopping strength). Both methods give qualitatively consistent results. With variational ansatz we go to even larger on-site attraction. At very low spin-up filling fraction, we observe the polaron-to-dimeron transition, in agreement with the continuum case. Upon increasing the filling fraction, however, the transition shifts to higher values of | U | / t and the transition disappears beyond a filling fraction of about 20 % . In this region, the polaron state always gives a lower energy and has a finite quasiparticle residue. Our findings are directly relevant to cold atom experiments with 2D optical lattices: A small and finite density of spin-down impurities in the ground state will form a superfluid at strong coupling at low spin-up filling fraction. Above some critical spin-up filling fraction, the system is expected to remain a normal Fermi liquid, even in the strong-coupling limit.
Numerical results obtained from Monte Carlo simulation in the ground state and at finite temperatures as well as the exact diagonalization and the transfer matrix method are reported. Efficiency of the Monte Carlo method in the ground state is examined. Spin, charge and superconducting correlations are investigated for the Hubbard and the t - J model in one and two dimensions. The momentum distribution in the one-dimensional Hubbard model shows fermi-liquid-like behavior at least in two orders of magnitude smaller energy scale than the band width. The short-range incommensurate spin correlation is observed and analyzed in the doped Hubbard and the t - J model both in one and two dimensions. In the ground state, the superconducting correlation shows the absence of system size dependence in two dimensions at the filling smaller than 0.8 and in one dimension at any filling. The binding energy of two fermions and the spin and charge distortion around an itinerant fermion are also discussed in the t - J and the coupled spin-fermion model.
Low energy properties of the metallic state of the two-dimensional t - J model are presented for second neighbor hopping with hole-doping ( ) and electron-doping ( ), with various superexchange energy J . We use a closed set of equations for the Greens functions obtained from the extremely correlated Fermi liquid theory . These equations reproduce the known low energies features of the large U Hubbard model in infinite dimensions. The density and temperature dependent quasiparticle weight, decay rate and the peak spectral heights over the Brillouin zone are calculated. We also calculate the resistivity, Hall conductivity, Hall number and cotangent Hall angle. The spectral features display high thermal sensitivity at modest T for density , implying a suppression of the effective Fermi-liquid temperature by two orders of magnitude relative to the bare bandwidth. The cotangent Hall angle exhibits a T 2 behavior at low T , followed by an interesting kink at higher T . The Hall number exhibits strong renormalization due to correlations. Flipping the sign of changes the curvature of the resistivity versus T curves between convex and concave. Our results provide a natural route for understanding the observed difference in the temperature dependent resistivity of strongly correlated electron-doped and hole-doped matter.