© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:1409.1103 · Physical Review B

Ground-state phase diagram of the half-filled bilayer Hubbard model

Michael Golor, Timo Reckling, Laura Classen, Michael M. Scherer, Stefan Weßel

Published Nov 18, 2014• 23 citations

Abstract

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.

DOI: 10.1103/physrevb.90.195131

Related papers

The phase diagram of the square lattice bilayer Hubbard model: a variational Monte Carlo study

Robert Rüger, Luca F. Tocchio, Roser Valentí +1 more

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.

arXiv:1311.6504•New Journal of Physics•Mar 7, 2014•28 citations

Phase diagram of the square lattice bilayer Hubbard model: a variational Monte Carlo study

Rüger, R., Luca F. Tocchio, Valenti, R. +1 more

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.

•Jan 1, 2014•26 citations

Band insulator to Mott insulator transition in a bilayer Hubbard model

S. S. Kancharla, Satoshi Okamoto

The ground-state phase diagram of the half-filled repulsive Hubbard model in a bilayer is investigated using cluster dynamical mean-field theory. For weak to intermediate values of Coulomb repulsion $U$, the system undergoes a transition from a Mott-insulating phase to a metallic phase and then onto a band-insulating phase as the interlayer hopping is increased. In the strong-coupling case, the model exhibits a direct crossover from a Mott-insulating phase to a band-insulating phase. These results are robust with respect to the presence or absence of magnetic order.

arXiv:cond-mat/0703728•Physical Review B•May 15, 2007•86 citations

Phase diagram of the half-filled Hubbard chain with next-nearest-neighbor hopping

S. Daul, R. M. Noack

We investigate the ground-state phase diagram of the half-filled one-dimensional Hubbard model with next-nearest-neighbor hopping using the density-matrix renormalization-group technique as well as an unrestricted Hartree-Fock approximation. We find commensurate and incommensurate disordered magnetic insulating phases and a spin-gapped metallic phase in addition to the one-dimensional Heisenberg phase. At large on-site Coulomb repulsion U, we make contact with the phase diagram of the frustrated Heisenberg chain, which has spin-gapped phases for sufficiently large frustration. For weak U, sufficiently large next-nearest-neighbor hopping ${t}_{2}$ leads to a band structure with four Fermi points rather than two, producing a spin-gapped metallic phase. As U is increased in this regime, the system undergoes a Mott-Hubbard transition to a frustrated antiferromagnetic insulator.

arXiv:cond-mat/9907256•Physical review. B, Condensed matter•Jan 15, 2000•53 citations

Bilayer Hubbard model: Analysis based on the fermionic sign problem

Yingping Mou, Rubem Mondaini, Richard T. Scalettar

The bilayer Hubbard model describes the antiferromagnet to spin singlet transition and, potentially, aspects of the physics of unconventional superconductors. Despite these important applications, significant aspects of its phase diagram in the interplane hopping ${t}_{\ensuremath{\perp}}$ versus on-site interaction $U$ parameter space, at half filling, are largely in disagreement. Here we provide an analysis making use of the average sign of weights over the course of the importance sampling in quantum Monte Carlo simulations to resolve several central open questions. Specifically, this metric of the weights clarifies the finite-sized metallic regimes at small $U$. Furthermore, at strong interactions, it points to the existence of a crossover from a correlated to uncorrelated band insulator not yet explored in a variety of existing, unbiased numerical methods. Our paper demonstrates the versatility of using properties of the weights in quantum Monte Carlo simulations to reveal important physical characteristics of the models under study.

arXiv:2205.02377•Physical review. B./Physical review. B•Sep 12, 2022•13 citations

Quantum Monte Carlo Study of an Interaction-Driven Band-Insulator–to–Metal Transition

Norman Paris, Karim Bouadim, F. Hébert +2 more

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.

arXiv:cond-mat/0607707•Physical Review Letters•Jan 25, 2007•102 citations

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.

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.

The ground-state phase diagram of the half-filled repulsive Hubbard model in a bilayer is investigated using cluster dynamical mean-field theory. For weak to intermediate values of Coulomb repulsion

U

, the system undergoes a transition from a Mott-insulating phase to a metallic phase and then onto a band-insulating phase as the interlayer hopping is increased. In the strong-coupling case, the model exhibits a direct crossover from a Mott-insulating phase to a band-insulating phase. These results are robust with respect to the presence or absence of magnetic order.

We investigate the ground-state phase diagram of the half-filled one-dimensional Hubbard model with next-nearest-neighbor hopping using the density-matrix renormalization-group technique as well as an unrestricted Hartree-Fock approximation. We find commensurate and incommensurate disordered magnetic insulating phases and a spin-gapped metallic phase in addition to the one-dimensional Heisenberg phase. At large on-site Coulomb repulsion U, we make contact with the phase diagram of the frustrated Heisenberg chain, which has spin-gapped phases for sufficiently large frustration. For weak U, sufficiently large next-nearest-neighbor hopping

{t}_{2}

leads to a band structure with four Fermi points rather than two, producing a spin-gapped metallic phase. As U is increased in this regime, the system undergoes a Mott-Hubbard transition to a frustrated antiferromagnetic insulator.

The bilayer Hubbard model describes the antiferromagnet to spin singlet transition and, potentially, aspects of the physics of unconventional superconductors. Despite these important applications, significant aspects of its phase diagram in the interplane hopping ${t}_{\ensuremath{\perp}}$ versus on-site interaction

U

parameter space, at half filling, are largely in disagreement. Here we provide an analysis making use of the average sign of weights over the course of the importance sampling in quantum Monte Carlo simulations to resolve several central open questions. Specifically, this metric of the weights clarifies the finite-sized metallic regimes at small

U

. Furthermore, at strong interactions, it points to the existence of a crossover from a correlated to uncorrelated band insulator not yet explored in a variety of existing, unbiased numerical methods. Our paper demonstrates the versatility of using properties of the weights in quantum Monte Carlo simulations to reveal important physical characteristics of the models under study.

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.