Using a self-consistent Hartree-Fock approximation we investigate the relative stability of various stripe phases in the extended $t$-$t$-$U$ Hubbard model. One finds that a negative ratio of next- to nearest-neighbor hopping $t/t<0$ expells holes from antiferromagnetic domains and reinforces the stripe order. Therefore the half-filled stripes not only accommodate holes but also redistribute them so that the kinetic energy is gained, and these stripes take over in the regime of $t/tsimeq -0.3$ appropriate for YBa$_2$Cu$_3$O$_{6+delta}$.
We investigate melting of stripe phases in the overdoped regime x>0.3 of the two-dimensional t-t-U Hubbard model, using a spin rotation invariant form of the slave boson representation. We show that the spin and charge order disappear simultaneously, and discuss a mechanism stabilizing bond-centered and site-centered stripe structures.
We study the half filled Hubbard chain including next-nearest-neighbor hopping $t$. The model has three phases: one insulating phase with dominant spin-density-wave correlations at large distances (SDWI), another phase with dominant spin-dimer correlations or dimerized insulator (DI), and a third one in which long distance correlations indicate singlet superconductivity (SS). The boundaries of the SDWI are accurately determined numerically through a crossing of excited energy levels equivalent to the jump in the spin Berry phase. The DI-SS boundary is studied using several indicators like correlation exponent $K_{rho}$, Drude weight $D_{c}$, localization parameter $z_{L}$ and charge gap $Delta_{c}$.
Two very different methods -- exact diagonalization on finite chains and a variational method -- are used to study the possibility of a metal-insulator transition in the symmetric half-filled periodic Anderson-Hubbard model. With this aim we calculate the density of doubly occupied $d$ sites as a function of various parameters. In the absence of on-site Coulomb interaction ($U_f$) between $f$ electrons, the two methods yield similar results. The double occupancy of $d$ levels remains always finite just as in the one-dimensional Hubbard model. Exact diagonalization on finite chains gives the same result for finite $U_f$, while the Gutzwiller method leads to a Brinkman-Rice transition at a critical value ($U_d^c$), which depends on $U_f$ and $V$.
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.
By means of the dynamical vertex approximation (D$Gamma$A) we include spatial correlations on all length scales beyond the dynamical mean field theory (DMFT) for the half-filled Hubbard model in three dimensions. The most relevant changes due to non-local fluctuations are: (i) a deviation from the mean-field critical behavior with the same critical exponents as for the three dimensional Heisenberg (anti)-ferromagnet and (ii) a sizable reduction of the Neel temperature ($T_N$) by $sim 30%$ for the onset of antiferromagnetic order. Finally, we give a quantitative estimate of the deviation of the spectra between D$Gamma$A and DMFT in different regions of the phase-diagram.