Critical behaviour near the metal-insulator transition of a doped Mott insulator

We have studied the critical behaviour of a doped Mott insulator near the metal-insulator transition for the infinite-dimensional Hubbard model using a linearized form of dynamical mean-field theory. The discontinuity in the chemical potential in the change from hole to electron doping, for U larger than a critical value U _c, has been calculated analytically and is found to be in good agreement with the results of numerical methods. We have also derived analytic expressions for the compressibility, the quasiparticle weight, the double occupancy and the local spin susceptibility near half-filling as functions of the on-site Coulomb interaction and the doping. Received 15 March 2001 and Received in final form 22 May 2001     

Anisotropy of the energy gap in the insulating phase of the U-t-t' Hubbard model

We apply a diagrammatic expansion method around the atomic limit () for the U-t-t ' Hubbard model at half filling and finite temperature by means of a continued fraction representation of the one-particle Green's function. From the analysis of the spectral function we find an energy dispersion relation with a modulation of the energy gap in the insulating phase. This anisotropy is compared with experimental ARPES results on insulating cuprates. Received 18 May 2000 and Received in final form 9 August 2000     

Unusually flat hole dispersion relation in the two-dimensional Hubbard model and restoration of coherence by addition of pair-hopping processes

The dispersion relation of a doped hole in the half-filled 2D Hubbard model is shown to follow a law around the and points in the Brillouin zone. Upon addition of pair-hopping processes this dispersion relation is unstable towards a law. The above follows from T=0 Quantum Monte-Carlo calculations of the single particle spectral function on lattices. We discuss finite dopings and argue that the added term restores coherence to charge dynamics and drives the system towards a d _{x2 - y2} superconductor. Received 22 March 1999     

Phase diagram of the three-dimensional Hubbard model at half filling

We investigate the phase diagram of the three-dimensional Hubbard model at half filling using quantum Monte Carlo (QMC) simulations. The antiferromagnetic Néel temperature is determined from the specific heat maximum in combination with finite-size scaling of the magnetic structure factor. Our results interpolate smoothly between the asymptotic solutions for weak and strong coupling, respectively, in contrast to previous QMC simulations. The location of the metal-insulator transition in the paramagnetic phase above is determined using the electronic compressibility as criterion. Received 11 April 2000 and Revised in final form 29 June 2000     

Remnant Fermi surface in a pseudogap regime of the two-dimensional Hubbard model at finite temperature

A precursor effect on the Fermi surface in the two-dimensional Hubbard model at finite temperatures near the antiferromagnetic instability is studied using three different itinerant approaches: the second order perturbation theory, the paramagnon theory (PT), and the two-particle self-consistent (TPSC) approach. In general, at finite temperature, the Fermi surface of the interacting electron systems is not sharply defined due to the broadening effects of the self-energy. In order to take account of those effects we consider the single-particle spectral function A(, 0) at the Fermi level, to describe the counterpart of the Fermi surface at T = 0. We find that the Fermi surface is destroyed close to the pseudogap regime due to the spin-fluctuation effects in both PT and TPSC approaches. Moreover, the top of the effective valence band is located around = (π/2,π/2) in agreement with earlier investigations on the single-hole motion in the antiferromagnetic background. A crossover behavior from the Fermi-liquid regime to the pseudogap regime is observed in the electron concentration dependence of the spectral function and the self-energy. Received 8 September 2000 and Received in final form 20 December 2000     

“Linearized” dynamical mean-field theory for the Mott-Hubbard transition

The Mott-Hubbard metal-insulator transition is studied within a simplified version of the Dynamical Mean-Field Theory (DMFT) in which the coupling between the impurity level and the conduction band is approximated by a single pole at the Fermi energy. In this approach, the DMFT equations are linearized, and the value for the critical Coulomb repulsion can be calculated analytically. For the symmetric single-band Hubbard model at zero temperature, the critical value is found to be given by 6 times the square root of the second moment of the free (U=0) density of states. This result is in good agreement with the numerical value obtained from the Projective Selfconsistent Method and recent Numerical Renormalization Group calculations for the Bethe and the hypercubic lattice in infinite dimensions. The generalization to more complicated lattices is discussed. The “linearized DMFT” yields plausible results for the complete geometry dependence of the critical interaction. Received 6 May 1999 and Received in final form 2 July 1999     

Analytical approach to the quantum-phase transition in the one-dimensional spinless Holstein model

We study the one-dimensional Holstein model of spinless fermions interacting with dispersion-less phonons by using a recently developed projector-based renormalization method (PRM). At half-filling the system shows a metal-insulator transition to a Peierls distorted state at a critical electron-phonon coupling where both phases are described within the same theoretical framework. The transition is accompanied by a phonon softening at the Brillouin zone boundary and a gap in the electronic spectrum. For different filling, the phonon softening appears away from the Brillouin zone boundary and thus reflects a different type of broken symmetry state.     

Random dispersion approximation for the Hubbard model

We use the Random Dispersion Approximation (RDA) to study the Mott-Hubbard transition in the Hubbard model at half band filling. The RDA becomes exact for the Hubbard model in infinite dimensions. We implement the RDA on finite chains and employ the Lanczos exact diagonalization method in real space to calculate the ground-state energy, the average double occupancy, the charge gap, the momentum distribution, and the quasi-particle weight. We find a satisfactory agreement with perturbative results in the weak- and strong-coupling limits. A straightforward extrapolation of the RDA data for L ≤ 14 lattice results in a continuous Mott-Hubbard transition at U_c≈W. We discuss the significance of a possible signature of a coexistence region between insulating and metallic ground states in the RDA that would correspond to the scenario of a discontinuous Mott-Hubbard transition as found in numerical investigations of the Dynamical Mean-Field Theory for the Hubbard model.     

Renormalized quasiparticles in antiferromagnetic states of the Hubbard model

We analyze the properties of the quasiparticle excitations of metallic antiferromagnetic states in a strongly correlated electron system. The study is based on dynamical mean field theory (DMFT) for the infinite dimensional Hubbard model with antiferromagnetic symmetry breaking. Self-consistent solutions of the DMFT equations are calculated using the numerical renormalization group (NRG). The low energy behavior in these results is then analyzed in terms of renormalized quasiparticles. The parameters for these quasiparticles are calculated directly from the NRG derived self-energy, and also from the low energy fixed point of the effective impurity model. From these the quasiparticle weight and the effective mass are deduced. We show that the main low energy features of the k-resolved spectral density can be understood in terms of the quasiparticle picture. We also find that Luttinger's theorem is satisfied for the total electron number in the doped antiferromagnetic state.     

Localization properties of two interacting particles in a quasi-periodic potential with a metal-insulator transition

We study the influence of many-particle interactions on a metal-insulator transition. We consider the two-interacting-particle problem for onsite interacting particles on a one-dimensional quasiperiodic chain, the so-called Aubry-André model. We show numerically by the decimation method and finite-size scaling that the interaction does not modify the critical parameters such as the transition point and the localization-length exponent. We compare our results to the case of finite density systems studied by means of the density-matrix renormalization scheme. Received 28 June 2001     

Extended Hubbard model with renormalized Wannier wave functions in the correlated state: beyond the parametrized models

The method used earlier for analysis of correlated nanoscopic systems is extended to infinite (periodic) s-band-like systems described by the Hubbard model. The optimized single-particle Wannier wave functions contained in the parameters of the extended Hubbard model (in the nearest-neghbor hopping (-t), in the magnitude of the intraatomic interaction U, and in other parameters) are determined explicitly in the correlated state for the electronic systems of various symmetries and dimensions: Hubbard chain, square and triangular planar lattices, and the three cubic lattices (SC, BCC, FCC). In effect, the evolution of the electronic properties as a function of interatomic distance R is obtained. The model parameters in most cases do not scale linearly with the lattice spacing and hence, their solution as a function of microscopic parameters reflects only qualitatively the system evolution. Also, the atomic energy changes with R and therefore should be included in the model analysis. The solutions in one dimension (D = 1) can be analyzed both rigorously (by making use of the Lieb–Wu solution) and compared with the approximate Gutzwiller treatments. In higher dimensions (D = 2 and 3) only the latter approach is possible to implement within the scheme. The renormalized single particle wave functions are almost independent of the choice of the scheme selected to diagonalize the Hamiltonian in the Fock space in D = 1 case. For dimensions D > 1 the qualitative behavior is independent of the structure considered. The wave-function size increases above the Mott-Hubbard localization threshold and gradually reaches the atomic limit value. The method can be extended to other approximation schemes, as stressed at the end.     

Dynamical mean-field study of the Mott transition in thin films

The correlation-driven transition from a paramagnetic metal to a paramagnetic Mott-Hubbard insulator is studied within the half-filled Hubbard model for a thin-film geometry. We consider simple-cubic films with different low-index surfaces and film thickness d ranging from d=1 (two-dimensional) up to d=8. Using the dynamical mean-field theory, the lattice (film) problem is self-consistently mapped onto a set of d single-impurity Anderson models which are indirectly coupled via the respective baths of conduction electrons. The impurity models are solved at zero temperature using the exact-diagonalization algorithm. We investigate the layer and thickness dependence of the electronic structure in the low-energy regime. Effects due to the finite film thickness are found to be the more pronounced the lower is the film-surface coordination number. For the comparatively open sc(111) geometry we find a strong layer dependence of the quasi-particle weight while it is much less pronounced for the sc(110) and the sc(100) film geometries. For a given geometry and thickness d there is a unique critical interaction strength U _c2 (d) at which all effective masses diverge and there is a unique strength U _c1 (d) where the insulating solution disappears. U _c2 (d) and U _c1 (d) gradually increase with increasing thickness eventually approaching their bulk values. A simple analytical argument explains the complete geometry and thickness dependence of U_c2. U_c1 is found to scale linearly with U_c2. Received 19 August 1998     

Phase separation in the strongly correlated Falicov-Kimball model in infinite dimensions

Phase separation in the strongly correlated Falicov-Kimball model in infinite dimensions is examined. We show that the phase separation can occur for any values of the interaction constant J^* when the site energy of the localized electrons is equal to zero. Electron-poor regions always have homogeneous state and electron-rich regions have chessboard state for , chessboard state or homogeneous state in dependence upon temperature for 0<J ^* <0.03 and homogeneous state for J ^* =0. For J ^* =0 and T=0, phase separation (segregation) occurs at .The obtained results are exact for the Bethe lattice with infinite number of the nearest neighbours. Received 1 December 1998 and Received in final form 12 April 1999     

Ground states of an extended Falicov-Kimball model

We present a non-perturbative study of an extended Falicov-Kimball model in one dimension. Working within the binary alloy interpretation, we include the spin of the itinerant electrons and a Hubbard interaction to model the inter-electron correlations. We derive an effective Ising model for the atomic configuration in order to show how the Hubbard term affects the stability of the phase separated states. Furthermore, we investigate the competition between the Mott insulator state of the itinerant electrons and the checkerboard phase of the spinless Falicov-Kimball model.     

Band effects on the conductivity for a two-band Hubbard model

The band effects on the conductivity of a one-dimensional two-band Hubbard model is studied based on the ground state energy analysis. It is found that the system with filling factor one is a metal at zero temperature if the on-site interaction U is smaller than a critical value U_c, and is an insulator if U is larger than U_c. The value of metal-insulator transition point U_c is obtained. This result is different from that of 1D single-band Hubbard model where the quantum phase transition point U_c=0. Therefore, the orbital degree of freedom plays an essential role in the states of matter.     

Phase diagram of the extended hubbard model with pair hopping interaction

A one-dimensional model of interacting electrons with on-site U, nearest-neighbor V, and pair-hopping interaction W is studied at half-filling using the continuum limit field theory approach. The ground state phase diagram is obtained for a wide range of coupling constants. In addition to the insulating spin-density wave (SDW) and charge-density wave (CDW) phases for large U and V, respectively, we identify a bond-charge-density-wave (BCDW) phase W < 0, | U - 2V| < | 2W| and a bond-spin-density-wave (BSDW) for W > 0, | U - 2V| < W. The possibility of bond-located ordering results from the site-off-diagonal nature of the pair-hopping term and is a special feature of the half-filled band case. The BCDW phase corresponding to an enhanced Peierls instability in the system. The BdSDW is an unconventional insulating magnetic phase, characterized by a gapless spin excitation spectrum and a staggered magnetization located on bonds between sites. The general ground state phase diagram including insulating, metallic, and superconducting phases is discussed. A transition to the η-superconducting phase at | U - 2V| ≪ 2t?W is briefly discussed. Received 20 February 2002 / Received in final form 11 April 2002 Published online 19 July 2002     

Magnetism and phase separation in the ground state of the Hubbard model

We discuss the ground state magnetic phase diagram of the Hubbard model off half filling within the dynamical mean-field theory. The effective single-impurity Anderson model is solved by Wilson's numerical renormalization group calculations, adapted to symmetry broken phases. We find a phase separated, antiferromagnetic state up to a critical doping for small and intermediate values of U, but could not stabilize a Néel state for large U and finite doping. At very large U, the phase diagram exhibits an island with a ferromagnetic ground state. Spectral properties in the ordered phases are discussed. Received 9 January 2002 Published online 25 June 2002     

Breakdown of Luttinger's theorem in two-orbital Mott insulators

An analysis of Luttinger's theorem shows that – contrary to recent claims – it is not valid for a generic Mott insulator. For a two-orbital Hubbard model with two electrons per site the crossover from a non-magnetic correlated insulating phase (Mott or Kondo insulator) to a band insulator is investigated. Mott insulating phases are characterized by poles of the self-energy and corresponding zeros in the Greens functions defining a “Luttinger surface” which is absent for band insulators. Nevertheless, the ground states of such insulators with two electrons per unit cell are adiabatically connected.     

The semiconductor-to-ferromagnetic-metal transition in FeSb2

We propose FeSb₂ to be a nearly ferromagnetic small gap semiconductor, hence a direct analog of FeSi. We find that despite different compositions and crystal structures, in the local density approximation with on-site Coulomb repulsion correction (LDA+U) method magnetic and semiconducting solutions for U=2.6 eV are energetically degenerate similar to the case of FeSi. For both FeSb₂ and FeSi (FeSi_1-xGe_x alloys) the underlying transition mechanism allows one to switch from a small gap semiconductor to a ferromagnetic metal with magnetic moment ≈1 μ_B per Fe ion with external magnetic field.