Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

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

Rüger, R., Luca F. Tocchio, Valenti, R., Claudius Gros

Published Jan 1, 2014• 26 citations
OpenAlex

Abstract

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.

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

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

Assessing the orbital selective Mott transition with variational wave functions

Luca F. Tocchio, Federico Arrigoni, Sandro Sorella +1 more

We study the Mott metal-insulator transition in the two-band Hubbard model with different hopping amplitudes t1 and t2 for the two orbitals on the two-dimensional square lattice by using non-magnetic variational wave functions, similarly to what has been considered in the limit of infinite dimensions by dynamical mean-field theory. We work out the phase diagram at half filling (i.e. two electrons per site) as a function of R = t2/t1 and the on-site Coulomb repulsion U, for two values of the Hund's coupling J = 0 and J/U = 0.1. Our results are in good agreement with previous dynamical mean-field theory calculations, demonstrating that the non-magnetic phase diagram is only slightly modified from infinite to two spatial dimensions. Three phases are present: a metallic one, for small values of U, where both orbitals are itinerant; a Mott insulator, for large values of U, where both orbitals are localized because of the Coulomb repulsion; and the so-called orbital-selective Mott insulator (OSMI), for small values of R and intermediate Us, where one orbital is localized while the other one is still itinerant. The effect of the Hund's coupling is two-fold: on one side, it favors the full Mott phase over the OSMI; on the other side, it stabilizes the OSMI at larger values of R.

arXiv:1505.07006•Journal of Physics Condensed Matter•Feb 16, 2016•16 citations

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

Michael Golor, Timo Reckling, Laura Classen +2 more

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.

arXiv:1409.1103•Physical Review B•Nov 18, 2014•23 citations

Quantum Monte Carlo Simulation of the two-dimensional ionic Hubbard model

B. Martinie

The Quantum Monte Carlo simulations of the ionic Hubbard model on a two-dimensional square lattice at half filling were performed. The method based on the direct-space, proposed by Suzuki and al., Hirsch and al., was used. Cycles of increasing and decreasing values of the Coulomb interaction $ U $ were performed for fixed temperature ($ kT=0.01 $). Results indicate that, at low temperature, the two insulator phases are separated by a metallic phase for weak to intermediate values of the staggered potential $ Δ$. For large Coulomb repulsion the system is in a Mott insulator with an antiferromagnetism order. On increasing and decreasing the Coulomb interaction $ U $ the metal-Mott insulator transition shows an hysteresis phenomenon while the metal-band insulator transition is continue. For large $ Δ$ it seems that the metallic region shrinks to a single metallic point. However, the band insulator to the Mott insulator transition is not direct for the studied model. A phase diagram is drawn for the temperature $ kT=0.01 $. For $ Δ=0.5 $ cycles of increasing and decreasing temperature were programmed for different values of the Coulomb interaction $ U $ . A behaviour change appears for $ U\simeq 1.75 $. This suggests that a crossover line divides the metallic region of the phase diagram.

arXiv:1002.0325•Feb 1, 2010•0 citations

Variational Monte Carlo Study of Two-Dimensional Multi-Orbital Hubbard Model on Square Lattice

Y Takenaka, N. Kawakami

We investigate the Mott transition in a two-dimensional two-orbital Hubbard model on a square lattice at half filling by using the variational Monte Carlo method. As a variational trial wave function, we consider Gutzwiller's on-site correlations and nearest-neighbor doublon-holon correlations. We calculate the ground-state energy and the other physical quantities, all of which clearly signal a first-order Mott transition. We thus obtain the phase diagram at zero temperature, where we find the coexistence region of the metal and the Mott insulator. We discuss how Hund's coupling affects the nature of the Mott transition in this system.

Journal of Physics Conference Series•Dec 17, 2012•5 citations
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.
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.
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 UUU, 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 study the Mott metal-insulator transition in the two-band Hubbard model with different hopping amplitudes t1 and t2 for the two orbitals on the two-dimensional square lattice by using non-magnetic variational wave functions, similarly to what has been considered in the limit of infinite dimensions by dynamical mean-field theory. We work out the phase diagram at half filling (i.e. two electrons per site) as a function of R = t2/t1 and the on-site Coulomb repulsion U, for two values of the Hund's coupling J = 0 and J/U = 0.1. Our results are in good agreement with previous dynamical mean-field theory calculations, demonstrating that the non-magnetic phase diagram is only slightly modified from infinite to two spatial dimensions. Three phases are present: a metallic one, for small values of U, where both orbitals are itinerant; a Mott insulator, for large values of U, where both orbitals are localized because of the Coulomb repulsion; and the so-called orbital-selective Mott insulator (OSMI), for small values of R and intermediate Us, where one orbital is localized while the other one is still itinerant. The effect of the Hund's coupling is two-fold: on one side, it favors the full Mott phase over the OSMI; on the other side, it stabilizes the OSMI at larger values of R.
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 Quantum Monte Carlo simulations of the ionic Hubbard model on a two-dimensional square lattice at half filling were performed. The method based on the direct-space, proposed by Suzuki and al., Hirsch and al., was used. Cycles of increasing and decreasing values of the Coulomb interaction $ U $ were performed for fixed temperature ($ kT=0.01 $). Results indicate that, at low temperature, the two insulator phases are separated by a metallic phase for weak to intermediate values of the staggered potential $ Δ$. For large Coulomb repulsion the system is in a Mott insulator with an antiferromagnetism order. On increasing and decreasing the Coulomb interaction $ U $ the metal-Mott insulator transition shows an hysteresis phenomenon while the metal-band insulator transition is continue. For large $ Δ$ it seems that the metallic region shrinks to a single metallic point. However, the band insulator to the Mott insulator transition is not direct for the studied model. A phase diagram is drawn for the temperature $ kT=0.01 $. For $ Δ=0.5 $ cycles of increasing and decreasing temperature were programmed for different values of the Coulomb interaction $ U $ . A behaviour change appears for $ U\simeq 1.75 $. This suggests that a crossover line divides the metallic region of the phase diagram.
We investigate the Mott transition in a two-dimensional two-orbital Hubbard model on a square lattice at half filling by using the variational Monte Carlo method. As a variational trial wave function, we consider Gutzwiller's on-site correlations and nearest-neighbor doublon-holon correlations. We calculate the ground-state energy and the other physical quantities, all of which clearly signal a first-order Mott transition. We thus obtain the phase diagram at zero temperature, where we find the coexistence region of the metal and the Mott insulator. We discuss how Hund's coupling affects the nature of the Mott transition in this system.