Pseudogap at hot spots in the two-dimensional Hubbard model at weak coupling

We analyze the interaction-induced renormalization of single-particle excitations in the two-dimensional Hubbard model at weak coupling using the Wick-ordered version of the functional renormalization group. The self energy is computed for real frequencies by integrating a flow equation with renormalized two-particle interactions. In the vicinity of hot spots, that is points where the Fermi surface intersects the umklapp surface, self energy effects beyond the usual quasi-particle renormalizations and damping occur near instabilities of the normal, metallic phase. Strongly enhanced renormalized interactions between particles at different hot spots generate a pronounced low-energy peak in the imaginary part of the self energy, leading to a pseudogap-like double-peak structure in the spectral function for single-particle excitations.Comment: 14 pages, 7 figure

Luttinger liquids with boundaries: Power-laws and energy scales

We present a study of the one-particle spectral properties for a variety of models of Luttinger liquids with open boundaries. We first consider the Tomonaga-Luttinger model using bosonization. For weak interactions the boundary exponent of the power-law suppression of the spectral weight close to the chemical potential is dominated by a term linear in the interaction. This motivates us to study the spectral properties also within the Hartree-Fock approximation. It already gives power-law behavior and qualitative agreement with the exact spectral function. For the lattice model of spinless fermions and the Hubbard model we present numerically exact results obtained using the density-matrix renormalization-group algorithm. We show that many aspects of the behavior of the spectral function close to the boundary can again be understood within the Hartree-Fock approximation. For the repulsive Hubbard model with interaction U the spectral weight is enhanced in a large energy range around the chemical potential. At smaller energies a power-law suppression, as predicted by bosonization, sets in. We present an analytical discussion of the crossover and show that for small U it occurs at energies exponentially (in -1/U) close to the chemical potential, i.e. that bosonization only holds on exponentially small energy scales. We show that such a crossover can also be found in other models.Comment: 16 pages, 9 figures included, submitted for publicatio

Renormalization-group analysis of the one-dimensional extended Hubbard model with a single impurity

We analyze the one-dimensional extended Hubbard model with a single static impurity by using a computational technique based on the functional renormalization group. This extends previous work for spinless fermions to spin-1/2 fermions. The underlying approximations are devised for weak interactions and arbitrary impurity strengths, and have been checked by comparing with density-matrix renormalization-group data. We present results for the density of states, the density profile and the linear conductance. Two-particle backscattering leads to striking effects, which are not captured if the bulk system is approximated by its low-energy fixed point, the Luttinger model. In particular, the expected decrease of spectral weight near the impurity and of the conductance at low energy scales is often preceded by a pronounced increase, and the asymptotic power laws are modified by logarithmic corrections.Comment: 36 pages, 13 figures, revised version as publishe

Effects of spin fluctuations in the t-J model

Recent experiments on the Fermi surface and the electronic structure of the cuprate-supercondutors showed the importance of short range antiferromagnetic correlations for the physics in these systems. Theoretically, features like shadow bands were predicted and calculated mainly for the Hubbard model. In our approach we calculate an approximate selfenergy of the $t$-$J$ model. Solving the $U=\infty$ Hubbard model in the Dynamical Mean Field Theory (DMFT) yields a selfenergy that contains most of the local correlations as a starting point. Effects of the nearest neighbor spin interaction $J$ are then included in a heuristical manner. Formally like in $J$-perturbation theory all ring diagrams, with the single bubble assumed to be purely local, are summed to get a correction to the DMFT-self engergy This procedure causes new bands and can furnish strong deformation of quasiparticle bands. % Our results are finally compared with %former approaches to the Hubbard model.Comment: 3 Pages, Latex, 2 Postscript-Figures submitted to Physica

Quantum phase transitions and collapse of the Mott gap in the d=1+ϵd=1+\epsilon dimensional half-filled Hubbard model

We study the low-energy asymptotics of the half-filled Hubbard model with a circular Fermi surface in $d=1+\epsilon$ continuous dimensions, based on the one-loop renormalization-group (RG) method. Peculiarity of the $d=1+\epsilon$ dimensions is incorporated through the mathematica structure of the elementary particle-partcile (PP) and particle-hole (PH) loops: infrared logarithmic singularity of the PH loop is smeared for $\epsilon>0$. The RG flows indicate that a quantum phase transition (QPT) from a metallic phase to the Mott insulator phase occurs at a finite on-site Coulomb repulsion $U$ for $\epsilon>0$. We also discuss effects of randomness.Comment: 12 pages, 10 eps figure

Charge gaps and quasiparticle bands of the ionic Hubbard model

The ionic Hubbard model on a cubic lattice is investigated using analytical approximations and Wilson's renormalization group for the charge excitation spectrum. Near the Mott insulating regime, where the Hubbard repulsion starts to dominate all energies, the formation of correlated bands is described. The corresponding partial spectral weights and local densities of states show characteristic features, which compare well with a hybridized-band picture appropriate for the regime at small $U$, which at half-filling is known as a band insulator. In particular, a narrow charge gap is obtained at half-filling, and the distribution of spectral quasi-particle weight reflects the fundamental hybridization mechanism of the model

Equation of motion approach to the Hubbard model in infinite dimensions

We consider the Hubbard model on the infinite-dimensional Bethe lattice and construct a systematic series of self-consistent approximations to the one-particle Green's function, $G^{(n)}(\omega),\ n=2,3,\dots\$ . The first $n-1$ equations of motion are exactly fullfilled by $G^{(n)}(\omega)$ and the $n$'th equation of motion is decoupled following a simple set of decoupling rules. $G^{(2)}(\omega)$ corresponds to the Hubbard-III approximation. We present analytic and numerical results for the Mott-Hubbard transition at half filling for $n=2,3,4$.Comment: 10pager, REVTEX, 8-figures not available in postscript, manuscript may be understood without figure

Crossover from Luttinger- to Fermi-liquid behavior in strongly anisotropic systems in large dimensions

We consider the low-energy region of an array of Luttinger liquids coupled by a weak interchain hopping. The leading logarithmic divergences can be re-summed to all orders within a self-consistent perturbative expansion in the hopping, in the large-dimension limit. The anomalous exponent scales to zero below the one-particle crossover temperature. As a consequence, coherent quasiparticles with finite weight appear along the whole Fermi surface. Extending the expansion self-consistently to all orders turns out to be crucial in order to restore the correct Fermi-liquid behavior.Comment: Shortened version to appear in Physical Review Letter

Hole dynamics in generalized spin backgrounds in infinite dimensions

We calculate the dynamical behaviour of a hole in various spin backgrounds in infinite dimensions, where it can be determined exactly. We consider hypercubic lattices with two different types of spin backgrounds. On one hand we study an ensemble of spin configurations with an arbitrary spin probability on each sublattice. This model corresponds to a thermal average over all spin configurations in the presence of staggered or uniform magnetic fields. On the other hand we consider a definite spin state characterized by the angle between the spins on different sublattices, i.e a classical spin system in an external magnetic field. When spin fluctuations are considered, this model describes the physics of unpaired particles in strong coupling superconductors.Comment: Accepted in Phys. Rev. B. 18 pages of text (1 fig. included) in Latex + 2 figures in uuencoded form containing the 2 postscripts (mailed separately

Application of the Density Matrix Renormalization Group in momentum space

We investigate the application of the Density Matrix Renormalization Group (DMRG) to the Hubbard model in momentum-space. We treat the one-dimensional models with dispersion relations corresponding to nearest-neighbor hopping and $1/r$ hopping and the two-dimensional model with isotropic nearest-neighbor hopping. By comparing with the exact solutions for both one-dimensional models and with exact diagonalization in two dimensions, we first investigate the convergence of the ground-state energy. We find variational convergence of the energy with the number of states kept for all models and parameter sets. In contrast to the real-space algorithm, the accuracy becomes rapidly worse with increasing interaction and is not significantly better at half filling. We compare the results for different dispersion relations at fixed interaction strength over bandwidth and find that extending the range of the hopping in one dimension has little effect, but that changing the dimensionality from one to two leads to lower accuracy at weak to moderate interaction strength. In the one-dimensional models at half-filling, we also investigate the behavior of the single-particle gap, the dispersion of spinon excitations, and the momentum distribution function. For the single-particle gap, we find that proper extrapolation in the number of states kept is important. For the spinon dispersion, we find that good agreement with the exact forms can be achieved at weak coupling if the large momentum-dependent finite-size effects are taken into account for nearest-neighbor hopping. For the momentum distribution, we compare with various weak-coupling and strong-coupling approximations and discuss the importance of finite-size effects as well as the accuracy of the DMRG.Comment: 15 pages, 11 eps figures, revtex