Holes in a Mott insulator are represented by spinless fermions in the fermion-boson model introduced by Edwards. Although the physically interesting regime is for low to moderate fermion density, the model has interesting properties over the whole density range. It has previously been studied at half-filling in the one-dimensional (1D) case by numerical methods, in particular using exact diagonalization and the density matrix renormalization group (DMRG). In the present study the one-particle Green's function is calculated analytically by means of a decoupling scheme for the equations of motion, valid for arbitrary density in 1D, 2D and 3D with fairly large boson energy and zero boson relaxation parameter. The Green's function is used to compute some ground state properties, and the one-fermion spectral function, for fermion densities n = 0.1, 0.5 and 0.9 in the 1D case. The results are generally in good agreement with numerical results obtained using the DMRG and dynamical DMRG, and new light is shed on the nature of the ground state at different fillings. The Green's function approximation is sufficiently successful in 1D to justify future application to the 2D and 3D cases.
Holes in a Mott insulator are represented by spinless fermions in the fermion-boson model introduced by Edwards. Although the physically interesting regime is for low to moderate fermion density, the model has interesting properties over the whole density range. It has previously been studied at half-filling in the one-dimensional (1D) case by numerical methods, in particular using exact diagonalization and the density matrix renormalization group (DMRG). In the present study the one-particle Green's function is calculated analytically by means of a decoupling scheme for the equations of motion, valid for arbitrary density in 1D, 2D and 3D with fairly large boson energy and zero boson relaxation parameter. The Green's function is used to compute some ground state properties, and the one-fermion spectral function, for fermion densities n = 0.1, 0.5 and 0.9 in the 1D case. The results are generally in good agreement with numerical results obtained using the DMRG and dynamical DMRG, and new light is shed on the nature of the ground state at different fillings. The Green's function approximation is sufficiently successful in 1D to justify future application to the 2D and 3D cases.
We study the spectral function of interacting one-dimensional fermions for an integrable lattice model away from half-filling. The divergent power-law singularity of the spectral function near the single-particle or single-hole energy is described by an effective x-ray edge type model. At low densities and for momentum near the zone boundary, we find a second divergent singularity at higher energies which is associated with a two-particle bound state. We use the Bethe ansatz solution of the model to calculate the exact singularity exponents for any momentum and for arbitrary values of chemical potential and interaction strength in the critical regime. We relate the singularities of the spectral function to the long-time decay of the fermion Green's function and compare our predictions with numerical results from the time-dependent density-matrix renormalization group (tDMRG). Our results show that the tDMRG method is able to provide accurate time decay exponents in the cases of power-law decay of the Green's function. Some implications for the line shape of the dynamical structure factor away from half-filling are also discussed. In addition, we show that 85% of the total spectral weight in the zero field Heisenberg model near wave-vector $q=\ensuremath{\pi}$ is given by the two-spinon contribution.
Starting from the known Lieb and Wu solution of the one-dimension-Hubbard model in the U to infinity limit, the authors show how the spin-charge decoupling of the elementary excitations is responsible for several peculiar features in one-particle properties, such as momentum distribution, quasiparticle weight and the Green function. In particular they analyse in detail the structure of the one-hole Green function at half-filling, which has not been previously calculated by field theory methods due to the breakdown of conformal invariance. A rich structure is found with branch cut singularities at omega =+or-2 sin k but no simple poles. The non-trivial dependence on the momentum of the hole allows for hole propagation although the analytic structure of G(k, omega ) is quite different from that usually characterizing band insulators. These results provide a precise characterization of one-dimensional Mott insulators. The relationship between the branch cuts of the Green function and the finite-size scaling of the quasiparticle weight is also discussed together with its implications for the analysis of numerical data.
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.
We introduce a slave-fermion formulation in which to study the charge dynamics of the half-filled Hubbard model on the square lattice. In this description, the charge degrees of freedom are represented by fermionic holons and doublons and the Mott-insulating characteristics of the ground state are the consequence of holon–doublon bound-state formation. The bosonic spin degrees of freedom are described by the antiferromagnetic Heisenberg model, yielding long-ranged (Néel) magnetic order at zero temperature. Within this framework and in the self-consistent Born approximation, we perform systematic calculations of the average double occupancy, the electronic density of states, the spectral function and the optical conductivity. Qualitatively, our method reproduces the lower and upper Hubbard bands, the spectral-weight transfer into a coherent quasiparticle band at their lower edges and the renormalisation of the Mott gap, which is associated with holon–doublon binding, due to the interactions of both quasiparticle species with the magnons. The zeros of the Green function at the chemical potential give the Luttinger volume, the poles of the self-energy reflect the underlying quasiparticle dispersion with a spin-renormalised hopping parameter and the optical gap is directly related to the Mott gap. Quantitatively, the square-lattice Hubbard model is one of the best-characterised problems in correlated condensed matter and many numerical calculations, all with different strengths and weaknesses, exist with which to benchmark our approach. From the semi-quantitative accuracy of our results for all but the weakest interaction strengths, we conclude that a self-consistent treatment of the spin-fluctuation effects on the charge degrees of freedom captures all the essential physics of the antiferromagnetic Mott–Hubbard insulator. We remark in addition that an analytical approximation with these properties serves a vital function in developing a full understanding of the fundamental physics of the Mott state, both in the antiferromagnetic insulator and at finite temperatures and dopings.
We calculate the single-particle excitation spectrum of holes in the Emery model thereby extending and improving previous calculations. The system is considered at half filling (nh=1, one hole per CuO2 unit) and for hole doping, where the on-site hole-hole repulsions are kept finite. A paramagnetic form of the ground state is used. For the determination of the retarded Green's functions of copper and oxygen holes, the projection technique is applied solving the resulting equations of motions self-consistently. At half filling, the excitation spectrum exhibits a charge-transfer gap bounded by Zhang-Rice singlet states and the upper Hubbard band. Upon hole doping the flat singlet band crosses the Fermi level giving rise to a large Fermi surface at a hole concentration of nh=1.25. Moreover, spectral weight is shifted from the upper Hubbard band to the states near the Fermi energy. The calculated spectral densities, the singlet dispersion for the doped system, and the transfer of spectral weight are in good quantitative agreement with exact diagonalization results for 2\ifmmode\times\else\texttimes\fi{}2 CuO2 cluster.
We study the dynamical response of the half-filled one-dimensional Hubbard model for a range of interaction strengths U and temperatures T by a combination of numerical and analytical techniques. Using time-dependent density matrix renormalization group computations we find that the single-particle spectral function undergoes a crossover to a spin-incoherent Luttinger liquid regime at temperatures $T\ensuremath{\sim}J=4{t}^{2}/U$ for sufficiently large U>4t. At smaller values of U and elevated temperatures the spectral function is found to exhibit two thermally broadened bands of excitations, reminiscent of what is found in the Hubbard-I approximation. The dynamical density-density response function is shown to exhibit a finite-temperature resonance at low frequencies inside the Mott gap, with a physical origin similar to the Villain mode in gapped quantum spin chains. We complement our numerical computations by developing an analytic strong-coupling approach to the low-temperature dynamics in the spin-incoherent regime.