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

Novel stochastic approach for Gutzwiller matrix elements: Two-dimensional Hubbard model

P. G. McQueen, C. S. Wang

Published Jun 1, 1989• 6 citations
Publisher

Abstract

A new approach is presented for variational Monte Carlo simulations of Hubbard-Anderson models using Gutzwiller projection operators to restrict charge fluctuations in noninteracting wave functions. The Gutzwiller matrix elements are mapped onto a statistical model, which is evaluated using a stochastic algorithm that is normally applied to finite-temperature simulation. The method has been successfully tested for one-dimensional (1D) Hubbard and 1D Anderson models, and results are presented for the 2D Hubbard model. We find a region close to the half-filled limit where the Gutzwiller states of different magnetic ordering are very close in energy.

DOI
10.1103/physrevb.39.12414

Related papers

Variational Monte Carlo evaluations of Gutzwiller states for the Anderson-lattice model

P. G. McQueen, C. S. Wang

The electronic ground states for two- and three-dimensional Anderson-lattice clusters with large Hubbard repulsion are determined variationally within the Gutzwiller ansatz by use of a stochastic algorithm for exact evaluation of matrix elements of Gutzwiller states. The algorithm is similar to one normally used for finite-temperature simulation and is described in some detail. Variational states for a variety of Hamiltonian parameters and particle densities were considered. Of particular interest is the regime of small interband hopping, which is antiferromagnetic at densities near two particles per unit cell. We found that for some Hamiltonian-parameter regimes, doping the antiferromagnetic cluster with enough particles or holes can destroy the magnetic ordering and produce a state that resembles a heavy-Fermi-liquid state, a metal with a strongly-frequency-dependent self-energy. For a 4\ifmmode\times\else\texttimes\fi{}4 Anderson lattice, we present results of low-temperature quantum Monte Carlo simulations to evaluate the validity of the variational Gutzwiller ansatz.

Physical review. B, Condensed matter•Nov 1, 1991•6 citations

Variation studies of the 2D Hubbard model: less than half filled

Robert B. Jones, W Yeung

The authors consider the variational Monte Carlo method of studying fermion systems in condensed matter physics. This method is a powerful tool for studying the properties of interacting electrons which are responsible for phenomena like the new high-Tc superconductors. Using a parallel algorithm for updating the determinants of large complex matrices on transputer arrays they study the energy and spin correlations for a 12*12 and a 22*22 lattice. Using the Hubbard model and the Gutzwiller projected wavefunction they have evaluated these quantities as functions of band filling and coupling constant. They find that although the Gutzwiller state is stable for densities near half filling it is not at low densities. For the stable states they find antiferromagnetic correlations. Statistically the thermodynamic limit is already attained with a 12*12 lattice.

Journal of Physics Condensed Matter•Apr 2, 1990•1 citation

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

Variational Monte Carlo study of the spiral state in the Hubbard model

Yongjun Liu, Jinming Dong, Chang-De Gong +1 more

Choosing the Gutzwiller-type variational wave function and using the variational Monte Carlo (VMC) method, we have calculated the ground-state energies of two-dimensional (2D) Hubbard squares at temperature T=0 and have obtained the phase diagram in the small doping region (\ensuremath{\delta}=0--0.15). We do not find the spiral phase until at least U\ensuremath{\approxeq}20 in the half-filled and dilute-doping cases. This is different from the results obtained by C. Jayaprakash et al., who calculated the phase diagram of the 2D t-J model at temperature T=0 with the slave-fermion mean-field theory and found the spiral phase at much lower values of U/t.

Physical review. B, Condensed matter•Jul 1, 1993•1 citation

Finite-temperature Gutzwiller projection for strongly correlated electron systems

Wansheng Wang, Xiang-Mei He, Da Wang +3 more

We generalized the Gutzwiller projectional variational method for the ground state of strongly correlated electron systems to the case of finite temperature. Under the Gutzwiller approximation, we show that this maps to a finite temperature renormalized mean-field theory. As one of the key ingredients in the theory, we obtained an explicit expression of the projection entropy or the entropy change due to the projection. We illustrate the application of the theory to the Anderson impurity problem and the half-filled Hubbard model and compare the theory to more elaborate techniques. We find qualitative agreement. The theory can be applied to a wide variety of Hubbard, $t\text{\ensuremath{-}}J$, and Anderson impurity models.

Physical Review B•Sep 7, 2010•26 citations

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
A new approach is presented for variational Monte Carlo simulations of Hubbard-Anderson models using Gutzwiller projection operators to restrict charge fluctuations in noninteracting wave functions. The Gutzwiller matrix elements are mapped onto a statistical model, which is evaluated using a stochastic algorithm that is normally applied to finite-temperature simulation. The method has been successfully tested for one-dimensional (1D) Hubbard and 1D Anderson models, and results are presented for the 2D Hubbard model. We find a region close to the half-filled limit where the Gutzwiller states of different magnetic ordering are very close in energy.
The electronic ground states for two- and three-dimensional Anderson-lattice clusters with large Hubbard repulsion are determined variationally within the Gutzwiller ansatz by use of a stochastic algorithm for exact evaluation of matrix elements of Gutzwiller states. The algorithm is similar to one normally used for finite-temperature simulation and is described in some detail. Variational states for a variety of Hamiltonian parameters and particle densities were considered. Of particular interest is the regime of small interband hopping, which is antiferromagnetic at densities near two particles per unit cell. We found that for some Hamiltonian-parameter regimes, doping the antiferromagnetic cluster with enough particles or holes can destroy the magnetic ordering and produce a state that resembles a heavy-Fermi-liquid state, a metal with a strongly-frequency-dependent self-energy. For a 4\ifmmode\times\else\texttimes\fi{}4 Anderson lattice, we present results of low-temperature quantum Monte Carlo simulations to evaluate the validity of the variational Gutzwiller ansatz.
The authors consider the variational Monte Carlo method of studying fermion systems in condensed matter physics. This method is a powerful tool for studying the properties of interacting electrons which are responsible for phenomena like the new high-Tc superconductors. Using a parallel algorithm for updating the determinants of large complex matrices on transputer arrays they study the energy and spin correlations for a 12*12 and a 22*22 lattice. Using the Hubbard model and the Gutzwiller projected wavefunction they have evaluated these quantities as functions of band filling and coupling constant. They find that although the Gutzwiller state is stable for densities near half filling it is not at low densities. For the stable states they find antiferromagnetic correlations. Statistically the thermodynamic limit is already attained with a 12*12 lattice.
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.
Choosing the Gutzwiller-type variational wave function and using the variational Monte Carlo (VMC) method, we have calculated the ground-state energies of two-dimensional (2D) Hubbard squares at temperature T=0 and have obtained the phase diagram in the small doping region (\ensuremath{\delta}=0--0.15). We do not find the spiral phase until at least U\ensuremath{\approxeq}20 in the half-filled and dilute-doping cases. This is different from the results obtained by C. Jayaprakash et al., who calculated the phase diagram of the 2D t-J model at temperature T=0 with the slave-fermion mean-field theory and found the spiral phase at much lower values of U/t.
We generalized the Gutzwiller projectional variational method for the ground state of strongly correlated electron systems to the case of finite temperature. Under the Gutzwiller approximation, we show that this maps to a finite temperature renormalized mean-field theory. As one of the key ingredients in the theory, we obtained an explicit expression of the projection entropy or the entropy change due to the projection. We illustrate the application of the theory to the Anderson impurity problem and the half-filled Hubbard model and compare the theory to more elaborate techniques. We find qualitative agreement. The theory can be applied to a wide variety of Hubbard, $t\text{\ensuremath{-}}J$, and Anderson impurity models.
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.