No Arabic abstract
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.
We investigate the ionic Hubbard model (IHM) at half-filling in the limit of strong correlations and large ionic potential. The low energy effective Hamiltonian in this limit, obtained by a similarity transformation, is a modified $t-J$ model with effective second neighbour hopping terms. We explore the possibilities of d-wave pairing and extended s-wave pairing superconducting (SC) phases on a two dimensional square lattice at zero temperature within a Gutzwiller projected renormalized mean field theory. In the sector of solutions that forbid spin ordering, the system shows a finite non-zero d-wave as well as extended s-wave pairing amplitude for $Delta sim U gg t$. The width of the superconducting phase in $U-Delta$ regime shrinks with increase in $U$ and $Delta$, though the extended s-wave pairing phase is higher in energy than the d-wave pairing superconducting phase. But in a spin resolved renormalized mean field calculation, which allows for an antiferromagnetic (AF) order along with the d-wave or extended s-wave pairing, the SC phase is no longer viable and the system shows a direct transition from an AF ordered phase to a paramagnetic band insulator. Except for a thin sliver of a half-metallic AF phase close to the AF transition point, most of the AF ordered phase is a Mott insulator. We benchmarked the AF Mott insulator to band insulator transition within the Gutzwiller projected renormalized mean field theory against the dynamical mean field theory (DMFT) solved using continuous time quantum Monte-Carlo (CTQMC). Our work suggests that the ground state phase diagram of the IHM at half-filling in the limit of extreme correlations does not have any SC phase. The SC phase seen in the paramagnetic sector is a metastable phase, being higher in energy than the AF Mott insulator phase.
We demonstrate that supervised machine learning (ML) with entanglement spectrum can give useful information for constructing phase diagram in the half-filled one-dimensional extended Hubbard model. Combining ML with infinite-size density-matrix renormalization group, we confirm that bond-order-wave phase remains stable in the thermodynamic limit.
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.
We present electron and phonon spectral functions calculated from determinant quantum Monte Carlo simulations of the half-filled two-dimensional Hubbard-Holstein model on a square lattice. By tuning the relative electron-electron ($e$-$e$) and electron-phonon ($e$-$ph$) interaction strengths, we show the electron spectral function evolving between antiferromagnetic insulating, metallic, and charge density wave insulating phases. The phonon spectra concurrently gain a strong momentum dependence and soften in energy upon approaching the charge density wave phase. In particular, we study how the $e$-$e$ and $e$-$ph$ interactions renormalize the spectra, and analyze how the interplay of these interactions influence the spectral renormalizations. We find that the presence of both interactions suppresses the amount of renormalization at low energy, thus allowing the emergence of a metallic phase. These findings demonstrate the importance of considering the influence of multiple interactions in spectroscopically determining any one interaction strength in strongly correlated materials.
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.