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 the Physical Society of Japan

Variational Monte-Carlo Studies of Hubbard Model. III. Intersite Correlation Effects

Hisatoshi Yokoyama, Hiroyuki Shiba

Published Oct 15, 1990• 133 citations
Publisher

Abstract

New variational wave functions are introduced to describe the Hubbard model in strong correlation. Intersite correlation (or Jastrow) factors are taken into account beyond the Gutzwiller wave function (GF). A study is made on the one-dimensional lattice and the two-dimensional square lattice at half-filling as well as less-than-half-filling. Much improvement can be achieved thereby in the description of physical quantities such as the momentum distribution and spin and charge density correlation functions. In addition, the ground state energy is lowered appreciably compared with the GF. In this connection, the Mott-Hubbard-type metal-insulator transition is also re-examined.

DOI
10.1143/jpsj.59.3669

Related papers

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

Variational Monte-Carlo Studies of Hubbard Model. II

Hisatoshi Yokoyama, Hiroyuki Shiba

As a continuation of a previous paper [J. Phys. Soc. Jpn. 56 (1987) 1490], the variational Monte-Carlo method is extended to include the antiferromagnetic long-range order. The theory is based on the Gutzwiller-type correlation factor and its effect is exactly taken into account by the Monte-Carlo procedure. An application is made to the half-filled-band case of one-dimensional lattice (50 sites), two-dimensional square lattice (up to 20×20 sites) and three-dimensional simple cubic lattice (6×6×6 sites). The result is qualitatively different from previous studies relying on the random-phase-type "Gutzwiller approximation." The variational energy for two and three dimensions is favorably compared with Hirsch's quantum Monte-Carlo data.

Journal of the Physical Society of Japan•Oct 15, 1987•153 citations

Variational Description of Mott Insulators

Manuela Capello, Federico Becca, Michele Fabrizio +2 more

The Gutzwiller wave function for a strongly correlated model can, if supplemented with a long-range Jastrow factor, provide a proper variational description of Mott insulators, so far unavailable. We demonstrate this concept in the prototypical one-dimensional t-t(') Hubbard model, where at half-filling we reproduce all known phases, namely, the ordinary Mott undimerized insulator with power-law spin correlations at small t(')/t, the spin-gapped metal above a critical t(')/t and small U, and the dimerized Mott insulator at large repulsion.

arXiv:cond-mat/0403430•Physical Review Letters•Jan 20, 2005•196 citations

Variational Monte Carlo Studies of Attractive Hubbard Model. I: Normal-State Properties

H. Yokoyama

Normal states of the attractive Hubbard model, mainly in two dimensions, are studied in the context of a transition from a Fermi liquid to an insulating or gapped state. A series of variational Monte Carlo calculations with better statistics was carried out to estimate accurately the expectation values of several many-body wave functions. Although a relatively clear crossover is observed, even for the plain Gutzwiller wave function, the states in both regimes are metallic. Meanwhile, a substantial metal-insulator transition takes place at ∣U∣∼W (band width) for an improved wave function, into which an intersite correlation is introduced by taking account of virtual states in the second-order perturbation in the infinite ∣U∣ limit. The critical value compares favorably with recent results of the dynamical-mean-field approximation. In contrast, a conventional Jastrow-type wave function scarcely improves the normal state. In addition, the issue of the Brinkman-Rice metal-insulator transition is reconsidered with much larger systems.

arXiv:cond-mat/0205281•Progress of Theoretical Physics•Jul 1, 2002•45 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

Hubbard Model in Strong Correlation Regime –Variational Monte-Carlo Studies on Singlet Liquid and Néel State–

Hisatoshi Yokoyama, Hiroyuki Shiba

The ground state of the Hubbard model on one-dimensional and two-dimensional square lattices is studied with the effective Hamiltonian derived for the strong correlation regime and Gutzwiller-type wave function. The non-half-filled band as well as the half-filled band is discussed. The variational Monte-Carlo method enables us to evaluate the energy for the singlet liquid (resonating-valence-bond state) and for Néel state. The result shows that the stability of the singlet liquid depends on the dimensionality of the system.

Journal of the Physical Society of Japan•Oct 15, 1987•155 citations
New variational wave functions are introduced to describe the Hubbard model in strong correlation. Intersite correlation (or Jastrow) factors are taken into account beyond the Gutzwiller wave function (GF). A study is made on the one-dimensional lattice and the two-dimensional square lattice at half-filling as well as less-than-half-filling. Much improvement can be achieved thereby in the description of physical quantities such as the momentum distribution and spin and charge density correlation functions. In addition, the ground state energy is lowered appreciably compared with the GF. In this connection, the Mott-Hubbard-type metal-insulator transition is also re-examined.
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.
As a continuation of a previous paper [J. Phys. Soc. Jpn. 56 (1987) 1490], the variational Monte-Carlo method is extended to include the antiferromagnetic long-range order. The theory is based on the Gutzwiller-type correlation factor and its effect is exactly taken into account by the Monte-Carlo procedure. An application is made to the half-filled-band case of one-dimensional lattice (50 sites), two-dimensional square lattice (up to 20×20 sites) and three-dimensional simple cubic lattice (6×6×6 sites). The result is qualitatively different from previous studies relying on the random-phase-type "Gutzwiller approximation." The variational energy for two and three dimensions is favorably compared with Hirsch's quantum Monte-Carlo data.
The Gutzwiller wave function for a strongly correlated model can, if supplemented with a long-range Jastrow factor, provide a proper variational description of Mott insulators, so far unavailable. We demonstrate this concept in the prototypical one-dimensional t-t(') Hubbard model, where at half-filling we reproduce all known phases, namely, the ordinary Mott undimerized insulator with power-law spin correlations at small t(')/t, the spin-gapped metal above a critical t(')/t and small U, and the dimerized Mott insulator at large repulsion.
Normal states of the attractive Hubbard model, mainly in two dimensions, are studied in the context of a transition from a Fermi liquid to an insulating or gapped state. A series of variational Monte Carlo calculations with better statistics was carried out to estimate accurately the expectation values of several many-body wave functions. Although a relatively clear crossover is observed, even for the plain Gutzwiller wave function, the states in both regimes are metallic. Meanwhile, a substantial metal-insulator transition takes place at ∣U∣∼W (band width) for an improved wave function, into which an intersite correlation is introduced by taking account of virtual states in the second-order perturbation in the infinite ∣U∣ limit. The critical value compares favorably with recent results of the dynamical-mean-field approximation. In contrast, a conventional Jastrow-type wave function scarcely improves the normal state. In addition, the issue of the Brinkman-Rice metal-insulator transition is reconsidered with much larger systems.
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 ground state of the Hubbard model on one-dimensional and two-dimensional square lattices is studied with the effective Hamiltonian derived for the strong correlation regime and Gutzwiller-type wave function. The non-half-filled band as well as the half-filled band is discussed. The variational Monte-Carlo method enables us to evaluate the energy for the singlet liquid (resonating-valence-bond state) and for Néel state. The result shows that the stability of the singlet liquid depends on the dimensionality of the system.