The Quantum Monte-Carlo simulations of the two-dimensional Hubbard model are presented for the half filling. The method based on the direct-space proposed by Suzuki and al., and Hirsch and al. was used. The states generated by this method are basis states in occupation number representation built with Wannier states localised on each site of the square array. The configurations of fermions can be observed on the real 2D array. An antiferromagnetic factor is defined and calculated for each temperature. The curves of energy, specific heat, conducivity and antiferromagnetic factor are presented for different values of the repulsive coulombian on site interaction U. There is a metal-insulator transition at low temperature for small values of U. This transition corresponds with a paramagnetic-ferromagnetic first order transition. Indeed, for these interaction values, the energy curves show a gap which is a characteristic of a first order transition. An hysteresis phenomenon appears on the conductivity curves. There is a behaviour change for U/t=3.5. For the values U>3.5 there is ferromagnetic-paramagnetic change without observable effect on the energy and the specific heat. The metal-insulator transition does not exist any more, the conductivity stays very small. Isotherms of the physical quantities versus U/t show a transition which seems to be the metal-insulator Mott transition. These results allow to draw a phase diagram with two first order transition lines.
We study the phase diagram of the ionic Hubbard model (IHM) at half filling on a Bethe lattice of infinite connectivity 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 ionic potential $\mathrm{\ensuremath{\Delta}}$ and the on-site Hubbard
U. We find that for a finite $\mathrm{\ensuremath{\Delta}}$ and at zero temperature, long-range antiferromagnetic (AFM) order sets in beyond a threshold
U=UAF via a first-order phase transition. For
U smaller than
UAF the system is a correlated band insulator. Both methods show a clear evidence for a quantum transition to a half-metal (HM) 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 methods have good qualitative and quantitative consistency in the intermediate-to-strong-coupling regime at zero temperature as well as at finite temperature. On increasing the temperature, the AFM order is lost via a first-order phase transition at a transition temperature ${T}_{AF}(U,\mathrm{\ensuremath{\Delta}})$ [or, equivalently, on decreasing
U below ${U}_{AF}(T,\mathrm{\ensuremath{\Delta}})$], within both methods, for weak to intermediate values of
U/t. 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. At a finite temperature
T,DMFT+CTQMC shows a second phase transition (not seen within
DMFT+IPT) on increasing
U beyond
UAF. At
UN>UAF, when the Neel temperature
TN for the effective Heisenberg model becomes lower than
T, the AFM order is lost via a second-order transition. For $U\ensuremath{\gg}\mathrm{\ensuremath{\Delta}}, {T}_{N}\ensuremath{\sim}{t}^{2}/U(1\ensuremath{-}{x}^{2})$, where $x=2\mathrm{\ensuremath{\Delta}}/U$ and thus
TN increases with increase in $\mathrm{\ensuremath{\Delta}}/U$. In the three-dimensional parameter space of $(U/t,\phantom{\rule{0.28em}{0ex}}T/t,\phantom{\rule{4.pt}{0ex}}\text{and}\phantom{\rule{4.pt}{0ex}}\mathrm{\ensuremath{\Delta}}/t)$, as
T increases, the surface of first-order transition at ${U}_{AF}(T,\mathrm{\ensuremath{\Delta}})$ and that of the second-order transition at ${U}_{N}(T,\mathrm{\ensuremath{\Delta}})$ approach each other, shrinking the range over which the AFM order is stable. There is a line of tricritical points that separates the surfaces of first- and second-order phase transitions.
We investigate the phases of the ionic Hubbard model in a two-dimensional square lattice using determinant quantum Monte Carlo. At half-filling, when the interaction strength or the staggered potential dominates, 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 study the transitions from band insulator to metal to Mott insulator in the ionic Hubbard model on a two-dimensional square lattice using determinant quantum Monte Carlo. Evaluation of the temperature dependence of the conductivity demonstrates that the metallic region extends for a finite range of interaction values. The Mott phase at strong coupling is accompanied by antiferromagnetic order. Inclusion of these intersite correlations changes the phase diagram qualitatively compared to dynamical mean field theory.
We investigate the phase diagram of the square lattice bilayer Hubbard model at half-filling with the variational Monte Carlo method for both the magnetic and the paramagnetic case as a function of the interlayer hopping and on-site Coulomb repulsion U . With this study we resolve some discrepancies in previous calculations based on the dynamical mean-field theory, and we are able to determine the nature of the phase transitions between metal, Mott insulator and band insulator. In the magnetic case we find only two phases: an antiferromagnetic Mott insulator at small for any value of U and a band insulator at large . At large U values we approach the Heisenberg limit. The paramagnetic phase diagram shows at small a metal to Mott insulator transition at moderate U values and a Mott to band insulator transition at larger U values. We also observe a re-entrant Mott insulator to metal transition and metal to band insulator transition for increasing in the range of . Finally, we discuss the phase diagrams obtained in relation to findings from previous studies based on different many-body approaches.
Abstract. We investigate the phase diagram of the square lattice bilayer Hubbard model at half filling with the variational Monte Carlo method for both the magnetic and the paramagnetic case as a function of the interlayer hopping t ⊥ and on-site Coulomb repulsion U. With this study we resolve some discrepancies in previous calculations based on the dynamical mean field theory, and we are able to determine the nature of the phase transitions between metal, Mott insulator and band insulator. In the magnetic case we find only two phases: An antiferromagnetic Mott insulator at small t ⊥ for any value of U and a band insulator at large t⊥. At large U values we approach the Heisenberg limit. The paramagnetic phase diagram shows at small t ⊥ a metal to Mott insulator transition at moderate U values and a Mott to band insulator transition at larger U values. We also observe a re-entrant Mott insulator to metal transition and metal to band insulator transition for increasing t ⊥ in the range of 5.5t < U < 7.5t. Finally, we discuss the obtained phase diagrams in relation to previous studies based on different many-body approaches. PACS numbers: 71.10.Fd, 71.27.+a, 71.30.+h 1.