The authors present here a comprehensive and numerically unbiased study for finite-temperature properties of the 3D fermionic Hubbard model. At half-filling, a complete temperature-interaction strength phase diagram of the model is established, including the highly accurate results of the N\'eel transition temperature and a multimessenger characterization of the metal-insulator crossover in the paramagnetic phase. The temperature dependence of various thermodynamic quantities and their underlying physics are thoroughly explored. The persistence of N\'eel AFM ordered phase to finite doping is further revealed.
Employing a combination of functional renormalization group calculations and projective determinantal quantum Monte Carlo simulations, we examine the Hubbard model on the square lattice bilayer at half filling. From this combined analysis, we obtain a comprehensive account of the ground-state phase diagram with respect to the extent of the system's metallic and (antiferromagnetically ordered) Mott-insulating as well as band-insulating regions. By means of an unbiased functional renormalization group approach, we exhibit the antiferromagnetic Mott-insulating state as the relevant instability of the free metallic state, induced by any weak finite on-site repulsion. Upon performing a careful analysis of the quantum Monte Carlo data, we resolve the difficulty of identifying this antiferromagnetic ground state for finite interlayer hopping in the weak-coupling regime, where nonmonotonic finite-size corrections are shown to relate to the two-sheeted Fermi surface structure of the metallic phase. On the other hand, quantum Monte Carlo simulations are well suited to identify the transition between the Mott-insulating phase and the band insulator in the intermediate-to-strong coupling regime. Here, we compare our numerical findings to indications for the transition region obtained from the functional renormalization group procedure.
The three-dimensional repulsive Hubbard model is investigated using the strong-coupling diagram technique. For half-filling, the boundary between paramagnetic and antiferromagnetic states is determined for the range of the Hubbard repulsion
4t≤U≤12t, where
t is the hopping integral between neighboring sites. Along this boundary, the density of states is calculated, and it demonstrates the Mott transition at
U≈9t. For
U≥6t and half-filling, the density of states has the shape inherent in the strong electron repulsion, while for
U=4t, its shape points to weak coupling. The dependence of the Néel temperature
TN on the electron concentration
nˉ is investigated for the cases
U=4t and
12t. In the former case,
TN decreases monotonously with
nˉ, while in the latter case, there is a plateau in the dependence near
nˉ=0.87. The plateau is connected to a reconstruction of the density of states caused by an effective weakening of electron coupling due to electron depopulation. For
U=12t and half-filling, the magnetic critical exponent
γ≈1.4, which is close to the value in the Heisenberg model. Some features resembling the first-order phase transition, revealing themselves in a finite crystal, are discussed.
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 thermodynamics of the three-dimensional Hubbard model at half filling on approach to the N\'eel transition by means of large-scale unbiased diagrammatic determinant Monte Carlo simulations. We obtain the transition temperature in the strongly correlated regime, as well as the temperature dependence of the energy, entropy, double occupancy, and nearest-neighbor spin correlation function. Our results improve the accuracy of previous unbiased studies and present accurate benchmarks in the ongoing effort to realize the antiferromagnetic state of matter with ultracold atoms in optical lattices.
The Hubbard model at temperatures above the Néel transition, despite being a paramagnet, can exhibit rich physics due to the interplay of Fermi surface, on-site interaction U and thermal fluctuations. Nevertheless, the understanding of the crossover physics remains only at a qualitative level, because of the intrinsically smooth behavior. Employing an improved variant of the numerically exact auxiliary-field quantum Monte Carlo algorithm equipped with numerical analytic continuation, we obtain a broad variety of thermodynamic and dynamical properties of the three-dimensional Hubbard model at half filling, quantitatively determine the crossover boundaries, and observe that the metal-insulator crossover regime, in which antiferromagnetic spin correlations appear strongest, exists over an extended regime in between the Fermi liquid for small U and the Mott insulator for large U. In particular, the location of the most rapid suppression of double occupancy as U increases, is found to fully reside in the metallic Fermi liquid regime, in contrast to the conventional intuition that it is a representative feature for entering the Mott insulator. Beside providing a reliable methodology for numerical study of crossover physics, our work can also serve as a timely and important guideline for the most recent optical lattice experiments.