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

Phase diagram of the Hubbard model: A variational wave-function approach

S. N. Coppersmith, Clare C. Yu

Published Jun 1, 1989• 18 citations
Publisher

Abstract

We explore the phase diagram of the Hubbard model by generalizing the Gutzwiller variational wave function to take into account nearest-neighbor spin correlations. The energies of various variational states are calculated and compared. At half filling the lowest energy state has long-range antiferromagnetic correlations, in agreement with previous results. Away from half filling, when the antiferromagnet no longer has the lowest energy, the Gutzwiller approximation yields a lower energy than that found for small spin polarons.

DOI
10.1103/physrevb.39.11464

Related papers

Normal and antiferromagnetic states of an extended Hubbard model

S. N. Coppersmith

Variational Monte Carlo calculations using extensions of the Gutzwiller variational wave function are used to study the phase diagram of an extended Hubbard model on a copper-oxygen lattice in two dimensions. Antiferromagnetism appears to occur at half-filling only when the oxygen p levels and copper d levels have substantially different energies. It is also found that strong nearest-neighbor repulsions can enhance the tendency for local-moment antiferromagnetism. The dependence of the variational parameters on the Hamiltonian parameters, particularly nearest-neighbor repulsions, is examined. The implications for models of high-temperature superconductivity are discussed.

Physical review. B, Condensed matter•May 1, 1990•14 citations

Ground State of Strongly Correlated Fermions: Short-Range Order

Yu. B. Kudasov

A new variational method is developed to calculate the ground state energy of Fermi systems with strong short-range correlations. A trial wave function of Gutzwiller's type contains additional variational parameters corresponding to configurations of pairs of nearest-neighbor sites. To evaluate the ground state energy, generalized Kikuchi's pseudo-ensemble method is used. The Hubbard model at half band-filling is investigated. The ground state energy of the paramagnetic phase is calculated for a chain, square and simple cubic lattices. It is shown that the short-range order lowers drastically the ground state energy of the Hubbard model at intermediate interaction strength. The paramagnetic phase of the Kondo-Hubbard model ($S=1/2$ and $S=5/2$) at half band-filling is investigated. The ground state energy, correlation functions and effective mass are calculated for chain, square and simple cubic lattices. A quantum phase transition was found for simple cubic lattice.

arXiv:cond-mat/0012018•Dec 1, 2000•0 citations

GENERALIZED GUTZWILLER ANSATZ FOR THE HALF-FILLED HUBBARD CHAIN

P. Fazekas, Karlo Penc

The well-known Gutzwiller wave function is generalized by including new variational parameters to control nearest-neighbour charge-charge, charge-spin, and spin-spin correlations. The non-magnetic state of the one-dimensional, half-filled Hubbard model is studied. Within the Gutzwiller approximation, the expression for the ground state energy can be worked out analytically. The correlation between empty and doubly occupied sites is found to play the most essential role. Minimization in the large-U limit shows that the Brinkman-Rice transition has been pushed to U → ∞, and the leading term of the ground state energy density is of order −t2/ U . In contrast to results obtained with the Gutzwiller wave function, we find that the band occupation number nk is monotonically decreasing both above and below kF. The dominant k–dependence is given by ~(t/U) cos k, in agreement with t/U–expansion results. nk has also a weak step-function component, with the discontinuity at kF vanishing as (t/U)2 in the limit U/t ≫ 1.

International Journal of Modern Physics B•Oct 1, 1988•27 citations

Variational ground states of the two-dimensional Hubbard model

D. Baeriswyl, David Eichenberger, M. Z. Menteshashvili

Recent refinements of analytical and numerical methods have improved our understanding of the ground-state phase diagram of the two-dimensional (2D) Hubbard model. Here, we focus on variational approaches, but comparisons with both quantum cluster and Gaussian Monte Carlo methods are also made. Our own ansatz leads to an antiferromagnetic ground state at half filling with a slightly reduced staggered order parameter (as compared to simple mean-field theory). Away from half filling, we find d-wave superconductivity, but confined to densities where the Fermi surface passes through the antiferromagnetic zone boundary (if hopping between both nearest-neighbour and next-nearest-neighbour sites is considered). Our results agree surprisingly well with recent numerical studies using the quantum cluster method. An interesting trend is found by comparing gap parameters Δ (antiferromagnetic or superconducting) obtained with different variational wave functions. Δ varies by an order of magnitude and thus cannot be taken as a characteristic energy scale. In contrast, the order parameter is much less sensitive to the degree of sophistication of the variational schemes, at least at and near half filling.

arXiv:0907.1593•New Journal of Physics•Jul 23, 2009•45 citations

Approximate solution of variational wave functions for strongly correlated systems: Description of bound excitons in metals and insulators

Balázs Hetényi

An approximate solution scheme, similar to the Gutzwiller approximation, is presented for the Baeriswyl and the Baeriswyl-Gutzwiller variational wave functions. The phase diagram of the one-dimensional Hubbard model as a function of interaction strength and particle density is determined. For the Baeriswyl wave function a metal-insulator transition is found at half filling, where the metallic phase $(U<{U}_{c})$ corresponds to the Hartree-Fock solution, the insulating phase is one with finite double occupations arising from bound excitons. This transition can be viewed as the ``inverse'' of the Brinkman-Rice transition. Close to but away from half filling, the $U>{U}_{c}$ phase displays a finite Fermi step, as well as double occupations originating from bound excitons. As the filling is changed away from half-filling bound excitons are suppressed. For the Baeriswyl-Gutzwiller wave function at half filling a metal-insulator transition between the correlated metallic and excitonic insulating state is found. Away from half-filling bound excitons are suppressed quicker than for the Baeriswyl wave function.

arXiv:1008.1272•Physical Review B•Sep 3, 2010•21 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 explore the phase diagram of the Hubbard model by generalizing the Gutzwiller variational wave function to take into account nearest-neighbor spin correlations. The energies of various variational states are calculated and compared. At half filling the lowest energy state has long-range antiferromagnetic correlations, in agreement with previous results. Away from half filling, when the antiferromagnet no longer has the lowest energy, the Gutzwiller approximation yields a lower energy than that found for small spin polarons.
Variational Monte Carlo calculations using extensions of the Gutzwiller variational wave function are used to study the phase diagram of an extended Hubbard model on a copper-oxygen lattice in two dimensions. Antiferromagnetism appears to occur at half-filling only when the oxygen p levels and copper d levels have substantially different energies. It is also found that strong nearest-neighbor repulsions can enhance the tendency for local-moment antiferromagnetism. The dependence of the variational parameters on the Hamiltonian parameters, particularly nearest-neighbor repulsions, is examined. The implications for models of high-temperature superconductivity are discussed.
A new variational method is developed to calculate the ground state energy of Fermi systems with strong short-range correlations. A trial wave function of Gutzwiller's type contains additional variational parameters corresponding to configurations of pairs of nearest-neighbor sites. To evaluate the ground state energy, generalized Kikuchi's pseudo-ensemble method is used. The Hubbard model at half band-filling is investigated. The ground state energy of the paramagnetic phase is calculated for a chain, square and simple cubic lattices. It is shown that the short-range order lowers drastically the ground state energy of the Hubbard model at intermediate interaction strength. The paramagnetic phase of the Kondo-Hubbard model (S=1/2S=1/2S=1/2 and S=5/2S=5/2S=5/2) at half band-filling is investigated. The ground state energy, correlation functions and effective mass are calculated for chain, square and simple cubic lattices. A quantum phase transition was found for simple cubic lattice.
The well-known Gutzwiller wave function is generalized by including new variational parameters to control nearest-neighbour charge-charge, charge-spin, and spin-spin correlations. The non-magnetic state of the one-dimensional, half-filled Hubbard model is studied. Within the Gutzwiller approximation, the expression for the ground state energy can be worked out analytically. The correlation between empty and doubly occupied sites is found to play the most essential role. Minimization in the large-U limit shows that the Brinkman-Rice transition has been pushed to U → ∞, and the leading term of the ground state energy density is of order −t2/ U . In contrast to results obtained with the Gutzwiller wave function, we find that the band occupation number nk is monotonically decreasing both above and below kF. The dominant k–dependence is given by ~(t/U) cos k, in agreement with t/U–expansion results. nk has also a weak step-function component, with the discontinuity at kF vanishing as (t/U)2 in the limit U/t ≫ 1.
Recent refinements of analytical and numerical methods have improved our understanding of the ground-state phase diagram of the two-dimensional (2D) Hubbard model. Here, we focus on variational approaches, but comparisons with both quantum cluster and Gaussian Monte Carlo methods are also made. Our own ansatz leads to an antiferromagnetic ground state at half filling with a slightly reduced staggered order parameter (as compared to simple mean-field theory). Away from half filling, we find d-wave superconductivity, but confined to densities where the Fermi surface passes through the antiferromagnetic zone boundary (if hopping between both nearest-neighbour and next-nearest-neighbour sites is considered). Our results agree surprisingly well with recent numerical studies using the quantum cluster method. An interesting trend is found by comparing gap parameters Δ (antiferromagnetic or superconducting) obtained with different variational wave functions. Δ varies by an order of magnitude and thus cannot be taken as a characteristic energy scale. In contrast, the order parameter is much less sensitive to the degree of sophistication of the variational schemes, at least at and near half filling.
An approximate solution scheme, similar to the Gutzwiller approximation, is presented for the Baeriswyl and the Baeriswyl-Gutzwiller variational wave functions. The phase diagram of the one-dimensional Hubbard model as a function of interaction strength and particle density is determined. For the Baeriswyl wave function a metal-insulator transition is found at half filling, where the metallic phase (U<Uc)(U<{U}_{c})(U<Uc​) corresponds to the Hartree-Fock solution, the insulating phase is one with finite double occupations arising from bound excitons. This transition can be viewed as the ``inverse'' of the Brinkman-Rice transition. Close to but away from half filling, the U>UcU>{U}_{c}U>Uc​ phase displays a finite Fermi step, as well as double occupations originating from bound excitons. As the filling is changed away from half-filling bound excitons are suppressed. For the Baeriswyl-Gutzwiller wave function at half filling a metal-insulator transition between the correlated metallic and excitonic insulating state is found. Away from half-filling bound excitons are suppressed quicker than for the Baeriswyl wave function.
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.