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 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.
The two-dimensional Hubbard model on the square lattice is studied in the presence of lattice distortions in the adiabatic approximation. The self energy is computed within perturbation theory up to second order, which provides a way for studying the quasiparticle dispersion. We compute numerically the second order contribution to the self-energy using a standard Fast Fourier Transform Algorithm for finite system sizes. The stability of the lattice distortions is investigated and a schematic phase diagram is drawn. The Fermi surface is analyzed for densities near to half filling, the presence of lattice distortion change some spectral properties of the model and gives an anisotropic interacting FS. The spectral function is calculated along several lines in momentum space and the renormalized quasiparticle dispersion is obtained. The behavior of the density of states is shown for different values of the intrasite repulsion U in the different phases.
The ground-state wave function and energy are calculated for the one-band Hubbard Hamiltonian on a 4\ifmmode\times\else\texttimes\fi{}4 cluster of the square lattice up to first and second order of the on-site Coulomb repulsion U, respectively. Changes induced by a deviation of the number of electrons from half-filling (one electron per lattice site) are closely followed. Charge-charge and spin-spin correlation functions are directly obtained from the ground-state wave function and show the existence of antiferromagnetic correlations that extend over the whole cluster at half-filling. Correlations among sites separated by more than a lattice site diminish very rapidly as holes are introduced into the cluster. The analysis is extended to the negative-U case. Ground-state properties change qualitatively.
Though most fermionic Mott insulators order at low temperatures, ordering is ancillary to their insulating behavior. Our emphasis here is on disentangling ordering from the intrinsic strongly correlated physics of a doped half-filled band. To this end, we focus on the two-dimensional Hubbard model. Because the charge gap arises from on-site correlations, we have been refining the nonperturbative approach of Matsumoto and Mancini [Phys. Rev. B 55, 2095 (1997)] which incorporates local physics. Crucial to this method is a self-consistent two-site dynamical cluster expansion which builds in the nearest-neighbor energy scale J. At half-filling, we find that the spectral function possesses a gap of order U and is devoid of any coherent quasiparticle peaks although ordering or charge fractionalization are absent. At low temperatures, local antiferromagnetic correlations emerge. In the doped case, we find that the Fermi surface exceeds the Luttinger volume. The breakdown of Luttinger's theorem in the underdoped regime is traced both to the dynamically generated Mott gap as well to a nonvanishing of the imaginary part of the self-energy at the Fermi level. Spectral weight transfer across the Mott gap also emerges as a ubiquitous feature of a doped Mott insulator and suggests that high- and low-energy scales are inseparable. Additionally in the underdoped regime, we find that a pseudogap exists in the single-particle density of states as well as in the heat capacity. The pseudogap (which is set by the energy scale
t2/U) is argued to be a ubiquitous feature of a lightly doped Mott state and simply represents the fact that hole transport involves double occupancy. In analogy with the Mott gap and antiferromagnetism, we propose that ordering may also accompany the formation of a pseudogap. We suggest a current pattern within a one-band model that preserves translational but breaks time-reversal symmetry along the canonical x and y axes but not along $x=\ifmmode\pm\else\textpm\fi{}y$ that is consistent with the experimental observations. Finally, we show that the Hall coefficient in a doped Mott insulator must change sign at a doping level
x<1/3. The sign change is tied to a termination of strong correlation physics in the doped Mott state.
Using exact diagonalization, we calculate the density of states of the two-dimensional Hubbard model on a 4\ifmmode\times\else\texttimes\fi{}4 square lattice at U/t=0.5, 4, and 10, and even number of electrons with filling factors n ranging from a quarter filling up to half filling. We compare the ground-state energy and density of states at U/t=0.5 and 4 with second-order perturbation theory in U/t in the paramagnetic phase, and find that while the agreement is reasonable at U/t=0.5, it becomes worse as the perturbatively determined (i.e., using Stoner's criterion) boundary of the paramagnetic to spin-density-wave instability is approached. In the strong coupling regime (U/t=10), we find reasonable agreement between the density of states of the Hubbard and the t-J model especially for low doping fractions. In general, we find that at half filling the filled states are separated from the empty states by a gap. At U/t=10, the density of states shows two bands clearly separated by a Mott-Hubbard gap of order U.
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.