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

arXiv:1104.0516 · Journal of the Physical Society of Japan

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

Hisatoshi Yokoyama, Tomoaki Miyagawa, Masao Ogata

Published Aug 15, 2011• 27 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1143/jpsj.80.084607

Related papers

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

Mott Transition and Spin Structures of Spin-1 Bosons in Two-Dimensional Optical Lattice at Unit Filling

Yuta Toga, Hiroki Tsuchiura, Makoto Yamashita +2 more

We study the ground state properties of spin-1 bosons in a two-dimensional optical lattice, by applying a variational Monte Carlo method to the S=1 Bose-Hubbard model on a square lattice at unit filling. A doublon-holon binding factor introduced in the trial state provides a noticeable improvement in the variational energy over the conventional Gutzwiller wave function and allows us to deal effectively with the inter-site correlations of particle densities and spins. We systematically show how spin-dependent interactions modify the superfluid-Mott insulator transitions in the S=1 Bose-Hubbard model due to the interplay between the density and spin fluctuations of bosons. Furthermore, regarding the magnetic phases in the Mott region, the calculated spin structure factor elucidates the emergence of nematic and ferromagnetic spin orders for antiferromagnetic ($U_2>0$) and ferromagnetic ($U_2<0$) couplings, respectively.

arXiv:1204.1175•Journal of the Physical Society of Japan•May 2, 2012•13 citations

Variational Monte Carlo Studies of the Hubbard Model in One- and Two-Dimensions –Off-Diagonal Intersite Correlation Effects–

Hiromi Otsuka

We investigate the ground state properties of the one- and the two-dimensional Hubbard models using a new variational Monte Carlo method. By generalizing the Gutzwiller variational wave function (GF) to take into account an off-diagonal intersite correlation, we examine the metal-insulator transition mechanism at the half filling. As Kaplan et al. first pointed out, the binding mechanism between doubly occupied sites and empty sites plays a crucial role for this transition; we confirm this by examining the real space correlation function between these two kinds of sites. Results for the spin and the charge structure factors are presented as well as the momentum distribution function and are discussed in connection with the pairing effect. We also find that in less-than half filling, the improvements of the structure factors at 2 k F are realized as compared with GF results.

Journal of the Physical Society of Japan•May 15, 1992•51 citations

Mott Transitions and d-Wave Superconductivity in Half-Filled-Band Hubbard Model on Square Lattice with Geometric Frustration

Hisatoshi Yokoyama, Masao Ogata, Yukio Tanaka

Mechanisms of Mott transitions and d x 2 Ày 2 -wave superconductivity (SC) are studied in the half-filledband Hubbard model on square lattices with a diagonal hopping term (t 0 ), using an optimization (or correlated) variational Monte Carlo method.In the trial wave functions, a doublon-holon binding effect is introduced in addition to the onsite Gutzwiller projection.We mainly treat a d-wave singlet state and a projected Fermi sea.In both wave functions, first-order Mott transitions without direct relevance to magnetic orders occur at U ¼ U c , which is approximately the bandwidth, for arbitrary t 0 =t.These transitions originate in the binding or unbinding of a doublon to a holon.d-wave SC appears in a narrow range immediately below U c .The robust d-wave superconducting correlation is necessarily accompanied by enhanced antiferromagnetic correlation; the strength of SC decreases, as t 0 =t increases.

arXiv:cond-mat/0607470•Journal of the Physical Society of Japan•Nov 15, 2006•118 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

Strongly correlated nature of d-wave superconductivity in cuprates

H. Yokoyama, Shun Tamura, T. Miyagawa +2 more

With the high-Tc cuprates in mind, properties of correlated d-wave superconducting (SC) states are studied for a Hubbard model on square lattices with a diagonal transfer (t-t’-U model), using a variational Monte Carlo method. We employ a simple wave function, which includes crucial parameters, in particular, a doublon-holon (D-H) binding factor important for correlated SC and normal states as doped Mott insulators. We first check that the range of dominant superconductivity is limited to a strongly correlated regime (U>W, U: onsite correlation strength, W: band width), and coincides with the effective range of the D-H binding factors. In this range of U/t and δ (doping rate), holons (in hole-doped cases) are classified into two types: doped holons and ones created as D-H pairs. Only the former holons participate in current, whereas the latter contribute to singlet-pair formation. Next, we show that the SC properties undergo a crossover at U=Uco∼W. For U<Uco, the behavior is appropriately interpreted with the conventional BCS mechanism, whereas for U>Uco (the regime of cuprates), a new idea is needed to understand the peculiar SC behavior.

Physics Procedia•Jan 1, 2012•1 citation
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.
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.
We study the ground state properties of spin-1 bosons in a two-dimensional optical lattice, by applying a variational Monte Carlo method to the S=1 Bose-Hubbard model on a square lattice at unit filling. A doublon-holon binding factor introduced in the trial state provides a noticeable improvement in the variational energy over the conventional Gutzwiller wave function and allows us to deal effectively with the inter-site correlations of particle densities and spins. We systematically show how spin-dependent interactions modify the superfluid-Mott insulator transitions in the S=1 Bose-Hubbard model due to the interplay between the density and spin fluctuations of bosons. Furthermore, regarding the magnetic phases in the Mott region, the calculated spin structure factor elucidates the emergence of nematic and ferromagnetic spin orders for antiferromagnetic (U2>0U_2>0U2​>0) and ferromagnetic (U2<0U_2<0U2​<0) couplings, respectively.
We investigate the ground state properties of the one- and the two-dimensional Hubbard models using a new variational Monte Carlo method. By generalizing the Gutzwiller variational wave function (GF) to take into account an off-diagonal intersite correlation, we examine the metal-insulator transition mechanism at the half filling. As Kaplan et al. first pointed out, the binding mechanism between doubly occupied sites and empty sites plays a crucial role for this transition; we confirm this by examining the real space correlation function between these two kinds of sites. Results for the spin and the charge structure factors are presented as well as the momentum distribution function and are discussed in connection with the pairing effect. We also find that in less-than half filling, the improvements of the structure factors at 2 k F are realized as compared with GF results.
Mechanisms of Mott transitions and d x 2 Ày 2 -wave superconductivity (SC) are studied in the half-filledband Hubbard model on square lattices with a diagonal hopping term (t 0 ), using an optimization (or correlated) variational Monte Carlo method.In the trial wave functions, a doublon-holon binding effect is introduced in addition to the onsite Gutzwiller projection.We mainly treat a d-wave singlet state and a projected Fermi sea.In both wave functions, first-order Mott transitions without direct relevance to magnetic orders occur at U ¼ U c , which is approximately the bandwidth, for arbitrary t 0 =t.These transitions originate in the binding or unbinding of a doublon to a holon.d-wave SC appears in a narrow range immediately below U c .The robust d-wave superconducting correlation is necessarily accompanied by enhanced antiferromagnetic correlation; the strength of SC decreases, as t 0 =t increases.
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.
With the high-Tc cuprates in mind, properties of correlated d-wave superconducting (SC) states are studied for a Hubbard model on square lattices with a diagonal transfer (t-t’-U model), using a variational Monte Carlo method. We employ a simple wave function, which includes crucial parameters, in particular, a doublon-holon (D-H) binding factor important for correlated SC and normal states as doped Mott insulators. We first check that the range of dominant superconductivity is limited to a strongly correlated regime (U>W, U: onsite correlation strength, W: band width), and coincides with the effective range of the D-H binding factors. In this range of U/t and δ (doping rate), holons (in hole-doped cases) are classified into two types: doped holons and ones created as D-H pairs. Only the former holons participate in current, whereas the latter contribute to singlet-pair formation. Next, we show that the SC properties undergo a crossover at U=Uco∼W. For U<Uco, the behavior is appropriately interpreted with the conventional BCS mechanism, whereas for U>Uco (the regime of cuprates), a new idea is needed to understand the peculiar SC behavior.