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

Journal of Physics Conference Series

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

Y Takenaka, N. Kawakami

Published Dec 17, 2012• 5 citations
Publisher

Abstract

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.

DOI
10.1088/1742-6596/400/3/032099

Related papers

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

Akihisa Koga, N. Kawakami, H. Yokoyama +1 more

We investigate electron correlations in the two‐dimensional degenerate Hubbard model on the square lattice. By means of variational Monte Carlo simulations including nearest‐neighbor doublon‐holon correlations as well as on‐site correlations, we discuss the Mott transitions in the two‐orbital system at zero temperature.

AIP conference proceedings•Jan 1, 2006•6 citations

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

Effect of Doublon–Holon Binding on Mott Transition–Variational Monte Carlo Study of Two-Dimensional Bose Hubbard Models

Hisatoshi Yokoyama, Tomoaki Miyagawa, Masao Ogata

To understand the mechanism of Mott transitions in case of no magnetic influence, superfluid-insulator (Mott) transitions in the S=0 Bose Hubbard model at unit filling are studied on the square and triangular lattices, using a variational Monte Carlo method. In trial many-body wave functions, we introduce various types of attractive correlation factors between a doubly-occupied site (doublon, D) and an empty site (holon, H), which play a central role for Mott transitions, in addition to the onsite repulsive (Gutzwiller) factor. By optimizing distance-dependent parameters, we study various properties of this type of wave functions. With a hint from the Mott transition arising in a completely D-H bound state, we propose an improved picture of Mott transitions, by introducing two characteristic length scales, the D-H binding length $ξ_{\rm dh}$ and the minimum D-D exclusion length $ξ_{\rm dd}$. Generally, a Mott transition occurs when $ξ_{\rm dh}$ becomes comparable to $ξ_{\rm dd}$. In the conductive (superfluid) state, domains of D-H pairs overlap with each other ($ξ_{\rm dh}>ξ_{\rm dd}$); thereby D and H can propagate independently as density carriers by successively exchanging the partners. In contrast, intersite repulsive Jastrow (D-D and H-H) factors have little importance for the Mott transition.

arXiv:1104.0516•Journal of the Physical Society of Japan•Aug 15, 2011•27 citations

Order parameter for the Mott transition in two dimensions: Small-world network

Chung-Pin Chou

We show how to define the appropriate network measure that characterizes the nonmagnetic Mott transition starting from the Hubbard model on a square lattice. To simulate the Mott transition, we employ the variational Monte Carlo method based on the Gutzwiller wave function with nearest-neighbor doublon-holon binding correlations. In order to search the Mott transition point, we use the doublon-holon correlation function as the weighted network representation for analyzing the half-filled electronic system. We then show that a network quantity, small-worldness, can be applied to defining the order parameter describing a hidden symmetry broken in the network representation, which allows us to capture the Mott transition point. As a result, we suggest that the picture of the Mott transition can be obtained using a mass of data entirely contained in topologies of the proper weighted networks.

arXiv:1402.2057•Feb 10, 2014•0 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

Effects of Long-Range Correlations on Nonmagnetic Mott Transitions in Hubbard Model on Square Lattice

智章 宮川, Hisatoshi Yokoyama

The mechanism of Mott transition in the Hubbard model on the square lattice is studied without explicit introduction of magnetic and superconducting correlations, using a variational Monte Carlo method. In the trial wave functions, we consider various types of binding factors between a doubly-occupied site (doublon, D) and an empty site (holon, H), like a long-range type as well as a conventional nearest-neighbor type, and add independent long-range D-D (H-H) factors. It is found that a wide choice of D-H binding factor leads to Mott transitions at critical values near the band width. We renew the D-H binding picture of Mott transitions by introducing two characteristic length scales, the D-H binding length l_{DH} and the minimum D-D distance l_{DD}, which we appropriately estimate. A Mott transition takes place at l_{DH}=l_{DD}. In the metallic regime (l_{DH}>l_{DD}), the domains of D-H pairs overlap with one another, thereby doublons and holons can move independently by exchanging the partners one after another. In contrast, the D-D factors give only a minor contribution to the Mott transition.

arXiv:1104.0109•Journal of the Physical Society of Japan•Jul 25, 2011•20 citations
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.
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.
To understand the mechanism of Mott transitions in case of no magnetic influence, superfluid-insulator (Mott) transitions in the S=0 Bose Hubbard model at unit filling are studied on the square and triangular lattices, using a variational Monte Carlo method. In trial many-body wave functions, we introduce various types of attractive correlation factors between a doubly-occupied site (doublon, D) and an empty site (holon, H), which play a central role for Mott transitions, in addition to the onsite repulsive (Gutzwiller) factor. By optimizing distance-dependent parameters, we study various properties of this type of wave functions. With a hint from the Mott transition arising in a completely D-H bound state, we propose an improved picture of Mott transitions, by introducing two characteristic length scales, the D-H binding length ξdhξ_{\rm dh}ξdh​ and the minimum D-D exclusion length ξddξ_{\rm dd}ξdd​. Generally, a Mott transition occurs when ξdhξ_{\rm dh}ξdh​ becomes comparable to ξddξ_{\rm dd}ξdd​. In the conductive (superfluid) state, domains of D-H pairs overlap with each other (ξdh>ξddξ_{\rm dh}>ξ_{\rm dd}ξdh​>ξdd​); thereby D and H can propagate independently as density carriers by successively exchanging the partners. In contrast, intersite repulsive Jastrow (D-D and H-H) factors have little importance for the Mott transition.
We show how to define the appropriate network measure that characterizes the nonmagnetic Mott transition starting from the Hubbard model on a square lattice. To simulate the Mott transition, we employ the variational Monte Carlo method based on the Gutzwiller wave function with nearest-neighbor doublon-holon binding correlations. In order to search the Mott transition point, we use the doublon-holon correlation function as the weighted network representation for analyzing the half-filled electronic system. We then show that a network quantity, small-worldness, can be applied to defining the order parameter describing a hidden symmetry broken in the network representation, which allows us to capture the Mott transition point. As a result, we suggest that the picture of the Mott transition can be obtained using a mass of data entirely contained in topologies of the proper weighted networks.
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.
The mechanism of Mott transition in the Hubbard model on the square lattice is studied without explicit introduction of magnetic and superconducting correlations, using a variational Monte Carlo method. In the trial wave functions, we consider various types of binding factors between a doubly-occupied site (doublon, D) and an empty site (holon, H), like a long-range type as well as a conventional nearest-neighbor type, and add independent long-range D-D (H-H) factors. It is found that a wide choice of D-H binding factor leads to Mott transitions at critical values near the band width. We renew the D-H binding picture of Mott transitions by introducing two characteristic length scales, the D-H binding length l_{DH} and the minimum D-D distance l_{DD}, which we appropriately estimate. A Mott transition takes place at l_{DH}=l_{DD}. In the metallic regime (l_{DH}>l_{DD}), the domains of D-H pairs overlap with one another, thereby doublons and holons can move independently by exchanging the partners one after another. In contrast, the D-D factors give only a minor contribution to the Mott transition.