No Arabic abstract
We investigate the ground-state phase diagram of the Hubbard model for the AB$_{N-1}$ chain with filling 1/N, where $N$ is the number of atoms per unit cell. In the strong-coupling limit, a charge transition takes place from a band insulator (BI) to a correlated insulator (CI) for increasing on-site repulsion $U$ and positive on-site energy difference $Delta$ (energy at A sites lower than at B sites). In the weak-coupling limit, a bosonization analysis suggests that for $N > 2$ the physics is qualitatively similar to the case $N = 2$ which has already been studied: an intermediate phase emerges, which corresponds to a bond-ordered ferroelectric insulator (FI) with spontaneously broken inversion symmetry. We have determined the quantum phase diagram for the cases $N = 3$ and $N = 4$ from the crossings of energy levels of appropriate excited states, which correspond to jumps in the charge and spin Berry phases, and from the change of sign of the localization parameter $z_{L}^{c}$. From these techniques we find that, quantitatively, the BI and FI phases are broader for $N > 2$ than when $N = 2$, in agreement with the bosonization analysis. Calculations of the Drude weight and $z_{L}^{c}$ indicate that the system is insulating for all parameters, with the possible exception of the boundary between the BI and FI phases.
We study the phase diagram of the ionic Hubbard model (IHM) at half-filling using dynamical mean field theory (DMFT), with two impurity solvers, namely, iterated perturbation theory (IPT) and continuous time quantum Monte Carlo (CTQMC). The physics of the IHM is governed by the competition between the staggered potential $Delta$ and the on-site Hubbard U. In both the methods we find that for a finite $Delta$ and at zero temperature, anti-ferromagnetic (AFM) order sets in beyond a threshold $U=U_{AF}$ via a first order phase transition below which the system is a paramagnetic band insulator. Both the methods show a clear evidence for a transition to a half-metal phase just after the AFM order is turned on, followed by the formation of an AFM insulator on further increasing U. We show that the results obtained within both the methods have good qualitative and quantitative consistency in the intermediate to strong coupling regime. On increasing the temperature, the AFM order is lost via a first order phase transition at a transition temperature $T_{AF}(U, Delta)$ within both the methods, for weak to intermediate values of U/t. But in the strongly correlated regime, where the effective low energy Hamiltonian is the Heisenberg model, IPT is unable to capture the thermal (Neel) transition from the AFM phase to the paramagnetic phase, but the CTQMC does. As a result, at any finite temperature T, DMFT+CTQMC shows a second phase transition (not seen within DMFT+IPT) on increasing U beyond $U_{AF}$. At $U_N > U_{AF}$, when the Neel temperature $T_N$ for the effective Heisenberg model becomes lower than T, the AFM order is lost via a second order transition. In the 3-dimensonal parameter space of $(U/t,T/t,Delta/t)$, there is a line of tricritical points that separates the surfaces of first and second order phase transitions.
The finite-temperature phase diagram of the Hubbard model in $d=3$ is obtained from renormalization-group analysis. It exhibits, around half filling, an antiferromagnetic phase and, between 30%--40% electron or hole doping from half filling, a new $tau $ phase in which the electron hopping strength $t$ asymptotically becomes infinite under repeated rescalings. Next to the $tau $ phase, a first-order phase boundary with very narrow phase separation (less than 2% jump in electron density) occurs. At temperatures above the $tau $ phase, an incommensurate spin modulation phase is indicated. In $d=2$, we find that the Hubbard model has no phase transition at finite temperature.
We investigate the phases of the ionic Hubbard model in a two-dimensional square lattice using determinant quantum Monte Carlo (DQMC). At half-filling, when the interaction strength or the staggered potential dominate we find Mott and band insulators, respectively. When these two energies are of the same order we find a metallic region. Charge and magnetic structure factors demonstrate the presence of antiferromagnetism only in the Mott region, although the externally imposed density modulation is present everywhere in the phase diagram. Away from half-filling, other insulating phases are found. Kinetic energy correlations do not give clear signals for the existence of a bond-ordered phase.
We investigate paramagnetic metal-insulator transitions in the infinite-dimensional ionic Hubbard model at finite temperatures. By means of the dynamical mean-field theory with an impurity solver of the continuous-time quantum Monte Carlo method, we show that an increase in the interaction strength brings about a crossover from a band insulating phase to a metallic one, followed by a first-order transition to a Mott insulating phase. The first-order transition turns into a crossover above a certain critical temperature, which becomes higher as the staggered lattice potential is increased. Further, analysis of the temperature dependence of the energy density discloses that the intermediate metallic phase is a Fermi liquid. It is also found that the metallic phase is stable against strong staggered potentials even at very low temperatures.
We investigate the phase diagram of the half-filled SU(N) Hubbard-Heisenberg model with hopping t, exchange J and Hubbard U, on a square lattice. In the large-N limit, and as a function of decreasing values of t/J, the model shows a transition from a d-density wave state to a spin dimerized insulator. A similar behavior is observed at N=6 whereas at N=2 a spin density wave insulating ground state is stabilized. The N=4 model, has a d-density wave ground state at large values of t/J which as a function of decreasing values of t/J becomes unstable to an insulating state with no apparent lattice and spin broken symmetries. In this state, the staggered spin-spin correlations decay as a power-law,resulting in gapless spin excitations at q = (pi,pi). Furthermore, low lying spin modes with small spectral weight are apparent around the wave vectors q = (0,pi) and q = (pi,0). This gapless spin liquid state is equally found in the SU(4) Heisenberg model in the self-adjoint antisymmetric representation. An interpretation of this state in terms of a pi-flux phase is offered. Our results stem from projective (T=0) quantum Monte-Carlo simulations on lattice sizes ranging up to 24 X 24.