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

Physical review. B, Condensed matter

Study of the accuracy of the Gutzwiller wave function for the two-dimensional Hubbard model

X. Q. Hong, J. E. Hirsch

Published Mar 1, 1990• 4 citations
Publisher

Abstract

We use a variational Monte Carlo method (VMC) and a quantum Monte Carlo method to study the two-dimensional half-filled Hubbard model. We use the Gutzwiller wave function as a trial wave function in VMC and measure energy, correlation functions, and momentum distribution for lattices up to 12\ifmmode\times\else\texttimes\fi{}12 sites. We then use a Monte Carlo algorithm to study the Hubbard model with the same boundary conditions as in VMC at very low temperatures. By comparing the results from these two methods, we assess the accuracy of the Gutzwiller wave function as a representation of the ground state of the Hubbard model. Major differences are found, implying that the Gutzwiller wave function is not a useful approximation to the ground state of the Hubbard model in two dimensions.

DOI
10.1103/physrevb.41.4410

Related papers

Gutzwiller wave function in the three-band Hubbard model: A variational Monte Carlo study

Akira Oguri, Tetsuro Asahata, Sadamichi Maekawa

The Gutzwiller wave function for the three-band Hubbard model on a two-dimensional ${\mathrm{CuO}}_{2}$ plane is studied by using the variational Monte Carlo (VMC) method in the limit ${\mathit{U}}_{\mathit{d}}$\ensuremath{\rightarrow}\ensuremath{\infty}, where ${\mathit{U}}_{\mathit{d}}$ is the Coulomb repulsion between holes on a Cu site. The VMC results for the energy and the fraction of d holes are compared with those of the Rice-Ueda type Gutzwiller approximation (GA). The difference between the VMC and the GA results are most pronounced at half-filling, and away from half-filling the two results agree well in both the hole-doped and the electron-doped cases. The doping dependence of the momentum distribution function is also studied.

arXiv:cond-mat/9307032•Physical review. B, Condensed matter•Mar 1, 1994•18 citations

Improved multi-variable variational Monte Carlo method examined by high-precision calculations of one-dimensional Hubbard model

Ryui Kaneko, Satoshi Morita, Masatoshi Imada

We revisit the accuracy of the variational Monte Carlo (VMC) method by taking an example of ground state properties for the one-dimensional Hubbard model. We start from the variational wave functions with the Gutzwiller and long-range Jastrow factor introduced by Capello et al. [Phys. Rev. B 72 , 085121 (2005)] and further improve it by considering several quantum-number projections and a generalized one-body wave function. We find that the quantum spin projection and total momentum projection greatly improve the accuracy of the ground state energy within 0.5% error, for both small and large systems at half filling. Besides, the momentum distribution function n ( k ) at quarter filling calculated up to 196 sites allows us direct estimate of the critical exponents of the charge correlations from the power-law behavior of n ( k ) near the Fermi wave vector. Estimated critical exponents well reproduce those predicted by the Tomonaga-Luttinger theory.

arXiv:1303.6041•Journal of Physics Conference Series•Aug 12, 2013•7 citations

Off-Diagonal Wave Function Monte Carlo Studies of Hubbard Model I

Takashi Yanagisawa, Soh Koike, Yamaji K.

We propose a Monte Carlo method, which is a hybrid method of the quantum Monte Carlo method and variational Monte Carlo theory, to study the Hubbard model. The theory is based on the off-diagonal and the Gutzwiller type correlation factors which are taken into account by a Monte Carlo algorithm. In the 4x4 system our method is able to reproduce the exact results obtained by the diagonalization. An application is given to investigate the half-filled band case of two-dimensional square lattice. The energy is favorably compared with quantum Monte Carlo data.

arXiv:cond-mat/0206033•Journal of the Physical Society of Japan•Nov 1, 1998•47 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

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

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
We use a variational Monte Carlo method (VMC) and a quantum Monte Carlo method to study the two-dimensional half-filled Hubbard model. We use the Gutzwiller wave function as a trial wave function in VMC and measure energy, correlation functions, and momentum distribution for lattices up to 12\ifmmode\times\else\texttimes\fi{}12 sites. We then use a Monte Carlo algorithm to study the Hubbard model with the same boundary conditions as in VMC at very low temperatures. By comparing the results from these two methods, we assess the accuracy of the Gutzwiller wave function as a representation of the ground state of the Hubbard model. Major differences are found, implying that the Gutzwiller wave function is not a useful approximation to the ground state of the Hubbard model in two dimensions.
The Gutzwiller wave function for the three-band Hubbard model on a two-dimensional CuO2{\mathrm{CuO}}_{2}CuO2​ plane is studied by using the variational Monte Carlo (VMC) method in the limit Ud{\mathit{U}}_{\mathit{d}}Ud​\ensuremath{\rightarrow}\ensuremath{\infty}, where Ud{\mathit{U}}_{\mathit{d}}Ud​ is the Coulomb repulsion between holes on a Cu site. The VMC results for the energy and the fraction of d holes are compared with those of the Rice-Ueda type Gutzwiller approximation (GA). The difference between the VMC and the GA results are most pronounced at half-filling, and away from half-filling the two results agree well in both the hole-doped and the electron-doped cases. The doping dependence of the momentum distribution function is also studied.
We revisit the accuracy of the variational Monte Carlo (VMC) method by taking an example of ground state properties for the one-dimensional Hubbard model. We start from the variational wave functions with the Gutzwiller and long-range Jastrow factor introduced by Capello et al. [Phys. Rev. B 72 , 085121 (2005)] and further improve it by considering several quantum-number projections and a generalized one-body wave function. We find that the quantum spin projection and total momentum projection greatly improve the accuracy of the ground state energy within 0.5% error, for both small and large systems at half filling. Besides, the momentum distribution function n ( k ) at quarter filling calculated up to 196 sites allows us direct estimate of the critical exponents of the charge correlations from the power-law behavior of n ( k ) near the Fermi wave vector. Estimated critical exponents well reproduce those predicted by the Tomonaga-Luttinger theory.
We propose a Monte Carlo method, which is a hybrid method of the quantum Monte Carlo method and variational Monte Carlo theory, to study the Hubbard model. The theory is based on the off-diagonal and the Gutzwiller type correlation factors which are taken into account by a Monte Carlo algorithm. In the 4x4 system our method is able to reproduce the exact results obtained by the diagonalization. An application is given to investigate the half-filled band case of two-dimensional square lattice. The energy is favorably compared with quantum Monte Carlo data.
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.
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.
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.