We investigate the normal selfenergy and the mass enhancement factor in the Hubbard model on the two-dimensional square lattice. Our purpose in this paper is to evaluate the mass enhancement factor more quantitatively than the conventional third order perturbation theory. We calculate it by expanding perturbatively up to the fourth order with respect to the on-site repulsion
U. We consider the cases that the system is near the half-filling, which are similar situations to high-
Tc cuprates. As results of the calculations, we obtain the large mass enhancement on the Fermi surface by introducing the fourth order terms. This is mainly originated from the fourth order particle-hole and particle-particle diagrams. Although the other fourth order terms have effect of reducing the effective mass, this effect does not cancel out the former mass enhancement completely and there remains still a large mass enhancement effect. In addition, we find that the mass enhancement factor becomes large with increasing the on-site repulsion
U and the density of state (DOS) at the Fermi energy
ρ(0). According to many current reseaches, such large
U and
ρ(0) enhance the effective interaction between quasiparticles, therefore the superconducting transition temperature
Tc increases. On the other hand, the large mass enhancement leads the reduction of the energy scale of quasiparticles, as a result,
Tc is reduced. When we discuss
Tc, we have to estimate these two competitive effects.
We investigate the Hubbard model on a two-dimensional square lattice by the perturbation expansion to the fourth order in the on-site Coulomb repulsion U. Numerically calculating all diagrams up to the fourth order in self-energy, we examine the convergence of perturbation series in the lattice system. We indicate that the coefficient of each order term rapidly decreases as in the impurity Anderson model for T > 0.1t in the half-filled case, but it holds in the doped case even at lower temperatures. Thus, we can expect that the convergence of perturbation expansion in U is very good in a wide parameter region also in the lattice system, except for T < 0.1t in the half-filled case. We next calculate the density of states in the fourth-order perturbation. In the half-filled case, the shape in a moderate correlation regime is quite different from the three peak structure in the second-order perturbation. Remarkable upper and lower Hubbard bands locate at w = +(-)U/2, and a pseudogap appears at the Fermi level w=0. This is considered as the precursor of the Mott-Hubbard antiferromagnetic structure. In the doped case, quasiparticles with very heavy mass are formed at the Fermi level. Thus, we conclude that the fourth-order perturbation theory overall well explain the asymptotic behaviors in a strong correlation regime.
We investigate the two-dimensional Hubbard model with the long-range Coulomb interaction. We assume the form of the effective Coulomb interaction, which is treated as a perturbation. For the half-filling case, the self-energy is numerically evaluated up to the second-order terms of the perturbation expansion with the full momentum dependence included. The k -mass is calculated from this self-energy. It is found that the first-order term makes a dominant contribution to the reduction of the k -mass.
The electron self-energy Sigma(k,\ensuremath{\omega}) corresponding to the one-band Hubbard Hamiltonian on clusters of the square lattice has been obtained as a function of the band filling up to second order in the interaction parameter U. The k dependence of the self-energy is completely taken into account. As a consequence, the imaginary part of the self-energy shows a linear behavior for fillings close to one electron per site whereas for large doping rates the quadratic \ensuremath{\omega} dependence characterizing Fermi liquids is recovered. The origin of this behavior has been investigated analytically: linear terms associated with nesting are shown to exist in any finite-dimensional lattice although its numerical relevance decreases very rapidly for space dimensionality larger than two. Implications of these results on the renormalization factor Z have been analyzed. From a practical point of view, we conclude that a standard second-order perturbative treatment is quantitatively precise for U values of the on-site Coulomb interaction smaller than the bandwidth of the noninteracting spectrum.
Single-particle excitation spectra of the two-dimensional Hubbard model on the square lattice near half filling and at zero temperature are investigated on the basis of the self-consistent projection operator method. The method guarantees a high accuracy of the spectra with high energy and high momentum resolutions. It takes into account long-range intersite correlations as well as the strong on-site correlations. Effects of nonlocal excitations are clarified by comparing the results with those of the single-site approximation. The calculated spectra verify the quantum Monte-Carlo results for finite temperatures. The spectra at the Fermi level yield a hole-like Fermi surface in the underdoped region and an electron-like Fermi surface in the overdoped region. From a numerical analysis of the momentum dependent effective mass and self-energy, it is concluded that a marginal Fermi-liquid like state persists even at finite doping concentrations in the strongly correlated region because a van Hove singularity is pinned to the Fermi surface. It is also found that a kink structure appears in the quasiparticle energy band in the same region. The kink is shown to be caused by a mixing between the quasiparticle band and an excitation band with strong short-range antiferromagnetic correlations. These results suggest an explanation for some of the unusual properties of the normal state in high-
Tc cuprates.
The density of states (DOS) of the two-dimensional Hubbard model is calculated in second-order perturbation with respect to the on-site Coulomb interaction U . The direct analytic continuation formula is presented to calculate the self-energy on the real-frequency axis. The momentum as well as the energy integral in the self-energy is seriously carried out with use of the Fast-Fourier-Transformation algorithm. Near the half-filling, the imaginary part of the self-energy becomes proportional to the energy in the vicinity of the Fermi level in our weak-coupling calculation. In the interacting case, the van Hove singularity (VHS) in DOS is smeared by the self-energy effect, but there still exists a sharp peak structure due to VHS.