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

A Variational Method to the Hubbard Model

Ichiya Sadakata

Published Jun 1, 1977• 1 citation
Publisher

Abstract

An approximate variational method of calculating the ground state energy of the Hubbard model is presented. The method is based on the trial wave function of Wigner. The result for the one-dimensional model is good compared with other approximate results. The extension of the method to finite temperatures is straightforward.
DOI
10.1143/jpsj.42.1810

Related papers

Correlation between perturbation and variation methods in the study of strongly correlated electron systems

Robinson O. Okanigbuan, JO Idiodi

The ground–state wave function and energy are calculated for two electrons subject to a one-band Hubbard Hamiltonian on a one dimensional lattice containing N electronic sites, N = 2,3,4,5,6, and a 3 x 3 cluster of the square lattice, using perturbation and variational methods. The results from these two approximation methods are then compared with the result from exact calculational method. Journal of the Nigerian Association of Mathematical Physics Vol. 10 2006: pp. 583-594

Journal of the Nigerian Association of Mathematical Physics•May 22, 2008•0 citations

High-order perturbation expansion for the two-dimensional Hubbard model using the Gutzwiller wave function

Zsolt Gulácsi, M. Gulácsi, Boldizsár Jankó

A variational approach of the two-dimensional Hubbard model is presented, based on the Gutzwiller wave function, making no use of the Gutzwiller approximation. We calculated exactly the coefficients of the (${\mathit{g}}^{2}$-1${)}^{\mathit{m}}$ expansion of the ground-state energy (g is the Gutzwiller variational parameter) up to the m=8 order. The m>8 terms were deduced with an acuracy O[(n/2${)}^{18}$] as a function of the electron concentration n. Double occupancy, expectation value of the interaction energy, average kinetic energy, and the minimized ground-state energy are determined and their electron concentration and dimensionality dependence are discussed. An evaluation of the nonanalyticities for any dimension is also presented.

Physical review. B, Condensed matter•Feb 15, 1993•10 citations

Pair theory of the Hubbard Hamiltonian: variational approach

H. Barentzen, Christo I. Ivanov, Marvin D. Girardeau

A novel variational approach is outlined for treating electron correlation in many-particle systems and is applied to the one-band Hubbard model. The authors' treatment is based on the concept of fermion pairs and their representation in terms of ideal-boson operators. In this representation approximate eigen-states of the Hamiltonian are obtained in the form of an antisymmetrised product of pair orbitals. For the ground state of the Hubbard model these orbitals are then determined by the use of the variational principle. Because of the complexity of the resulting equation a simplified ansatz is used for the orbitals, which leads to the same energy expression as is obtained by the traditional Bogolubov method. The trial states used in both approaches are, however, quite different. They find reasonable qualitative agreement with well established results for the Hubbard model in one dimension even in those cases where the Bogolubov or related approaches fail. They are also able to show that in the atomic limit the Bogolubov energy bound is the best upper bound that can be obtained by using a single determinant of pair orbitals. Therefore, the only way to improve the ground-state energy in the atomic limit is by means of configuration interaction.

Journal of Physics C Solid State Physics•Oct 20, 1987•3 citations

Ground-state energy of the one- and two-dimensional Hubbard model calculated by the method of Singwi, Tosi, Land, and Sjölander

M. R. Hedayati, Giovanni Vignale

By using the classic approximation method of Singwi, Tosi, Land, and Sj\"olander (STLS), we calculate the ground-state energy of the one- and two-dimensional Hubbard model in the half-filled antiferromagnetic state. In the one-dimensional case, by comparison with the exact result of Lieb and Wu, we find that the STLS correlation energy is generally better than that obtained from the random-phase approximation (RPA) and competes in quality with the best values obtained from variational calculations based on the antiferromagnetic Gutzwiller wave function (AFGWF). In particular, the STLS approximation yields a maximum in the absolute value of the correlation energy as a function of U, which is also found in the exact treatment, but not in the RPA. However, the numerical values of the STLS correlation energy are not as accurate as in the case of the uniform electron gas. In two dimensions, where the exact energy is not known, we find a similar relationship between the results of STLS, RPA, and AFGWF calculations, but the absolute value of the correlation energy is smaller. For large values of U, in both dimensions, the STLS method (like the RPA) runs into difficulties, due to the insufficient treatment of transverse spin fluctuations.

Physical review. B, Condensed matter•Nov 1, 1989•22 citations

Evaluation of the BCS approximation for the attractive Hubbard model in one dimension

F. Marsiglio

The ground-state energy and energy gap to the first excited state are calculated for the attractive Hubbard model in one dimension using both the Bethe ansatz equations and the variational BCS wave function. Comparisons are provided as a function of coupling strength and electron density. While the ground-state energies are always in very good agreement, the BCS energy gap is sometimes incorrect by an order of magnitude, particularly at half-filling. Finite-size effects are also briefly discussed for cases where an exact solution in the thermodynamic limit is not possible. In general, the BCS result for the energy gap is poor compared to the exact result.

arXiv:cond-mat/9609275•Physical review. B, Condensed matter•Jan 1, 1997•39 citations

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

X. Q. Hong, J. E. Hirsch

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.

Physical review. B, Condensed matter•Mar 1, 1990•4 citations
A variational approach of the two-dimensional Hubbard model is presented, based on the Gutzwiller wave function, making no use of the Gutzwiller approximation. We calculated exactly the coefficients of the (g2{\mathit{g}}^{2}g2-1)m{)}^{\mathit{m}})m expansion of the ground-state energy (g is the Gutzwiller variational parameter) up to the m=8 order. The m>8 terms were deduced with an acuracy O[(n/2)18{)}^{18})18] as a function of the electron concentration n. Double occupancy, expectation value of the interaction energy, average kinetic energy, and the minimized ground-state energy are determined and their electron concentration and dimensionality dependence are discussed. An evaluation of the nonanalyticities for any dimension is also presented.
A novel variational approach is outlined for treating electron correlation in many-particle systems and is applied to the one-band Hubbard model. The authors' treatment is based on the concept of fermion pairs and their representation in terms of ideal-boson operators. In this representation approximate eigen-states of the Hamiltonian are obtained in the form of an antisymmetrised product of pair orbitals. For the ground state of the Hubbard model these orbitals are then determined by the use of the variational principle. Because of the complexity of the resulting equation a simplified ansatz is used for the orbitals, which leads to the same energy expression as is obtained by the traditional Bogolubov method. The trial states used in both approaches are, however, quite different. They find reasonable qualitative agreement with well established results for the Hubbard model in one dimension even in those cases where the Bogolubov or related approaches fail. They are also able to show that in the atomic limit the Bogolubov energy bound is the best upper bound that can be obtained by using a single determinant of pair orbitals. Therefore, the only way to improve the ground-state energy in the atomic limit is by means of configuration interaction.
By using the classic approximation method of Singwi, Tosi, Land, and Sj\"olander (STLS), we calculate the ground-state energy of the one- and two-dimensional Hubbard model in the half-filled antiferromagnetic state. In the one-dimensional case, by comparison with the exact result of Lieb and Wu, we find that the STLS correlation energy is generally better than that obtained from the random-phase approximation (RPA) and competes in quality with the best values obtained from variational calculations based on the antiferromagnetic Gutzwiller wave function (AFGWF). In particular, the STLS approximation yields a maximum in the absolute value of the correlation energy as a function of U, which is also found in the exact treatment, but not in the RPA. However, the numerical values of the STLS correlation energy are not as accurate as in the case of the uniform electron gas. In two dimensions, where the exact energy is not known, we find a similar relationship between the results of STLS, RPA, and AFGWF calculations, but the absolute value of the correlation energy is smaller. For large values of U, in both dimensions, the STLS method (like the RPA) runs into difficulties, due to the insufficient treatment of transverse spin fluctuations.
The ground-state energy and energy gap to the first excited state are calculated for the attractive Hubbard model in one dimension using both the Bethe ansatz equations and the variational BCS wave function. Comparisons are provided as a function of coupling strength and electron density. While the ground-state energies are always in very good agreement, the BCS energy gap is sometimes incorrect by an order of magnitude, particularly at half-filling. Finite-size effects are also briefly discussed for cases where an exact solution in the thermodynamic limit is not possible. In general, the BCS result for the energy gap is poor compared to the exact result.
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.