No Arabic abstract
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 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 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 report a density functional theory plus dynamical mean field theory (DFT+DMFT) study of an oxide heterostructure of LaNiO$_3$ (LNO) bilayers in [111] direction interleaved with four atomic monolayers of LaAlO$_3$. DFT+$U$ optimizations yield two stable solutions: a uniform structure with equivalent NiO$_6$ octahedra, as well as a bond-disproportionated (BD) structure featuring a breathing distortion. For both structures, we construct the low-energy models describing the Ni $e_g$ states by means of Wannier projections supplemented by the Kanamori interaction, and solve them by DMFT. Using the continuous-time quantum Monte Carlo algorithm in the hybridization expansion, we study the temperature range between 145 and 450 K. For the uniform and the BD structure, we find similar phase diagrams that comprise four phases: a ferromagnetic metal (FM), a paramagnetic metal (PM), an antiferro-orbitally-ordered insulator (AOI), as well as a paramagnetic insulator (PI). By calculating momentum-resolved spectral functions on a torus and a cylinder, we demonstrate that the FM phase is not a Dirac metal, while both insulating phases are topologically trivial. By a comparison with available experimental data and model DMFT studies for the two-orbital Hubbard model, we suggest that LNO bilayers are in the AOI phase at room temperature.
We study finite-temperature transport properties of the one-dimensional Hubbard model using the density matrix renormalization group. Our aim is two-fold: First, we compute both the charge and the spin current correlation function of the integrable model at half filling. The former decays rapidly, implying that the corresponding Drude weight is either zero or very small. Second, we calculate the optical charge conductivity sigma(omega) in presence of small integrability-breaking next-nearest neighbor interactions (the extended Hubbard model). The DC conductivity is finite and diverges as the temperature is decreased below the gap. Our results thus suggest that the half-filled, gapped Hubbard model is a normal charge conductor at finite temperatures. As a testbed for our numerics, we compute sigma(omega) for the integrable XXZ spin chain in its gapped phase.
We study the charge conductivity of the one-dimensional repulsive Hubbard model at finite temperature using the method of dynamical quantum typicality, focusing at half filling. This numerical approach allows us to obtain current autocorrelation functions from systems with as many as 18 sites, way beyond the range of standard exact diagonalization. Our data clearly suggest that the charge Drude weight vanishes with a power law as a function of system size. The low-frequency dependence of the conductivity is consistent with a finite dc value and thus with diffusion, despite large finite-size effects. Furthermore, we consider the mass-imbalanced Hubbard model for which the charge Drude weight decays exponentially with system size, as expected for a non-integrable model. We analyze the conductivity and diffusion constant as a function of the mass imbalance and we observe that the conductivity of the lighter component decreases exponentially fast with the mass-imbalance ratio. While in the extreme limit of immobile heavy particles, the Falicov-Kimball model, there is an effective Anderson-localization mechanism leading to a vanishing conductivity of the lighter species, we resolve finite conductivities for an inverse mass ratio of $eta gtrsim 0.25$.