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

arXiv:1611.05174 · Journal of the Physical Society of Japan

Variational Wavefunction for the Periodic Anderson Model with Onsite Correlation Factors

Katsunori Kubo, Hiroaki Onishi

Published Nov 30, 2016• 4 citations
arXivHTMLPDF

Abstract

We propose a variational wavefunction containing parameters to tune the probabilities of all the possible onsite configurations for the periodic Anderson model. We call it the full onsite-correlation wavefunction (FOWF). This is a simple extension of the Gutzwiller wavefunction (GWF), in which one parameter is included to tune the double occupancy of the f electrons at the same site. We compare the energy of the GWF and the FOWF evaluated by the variational Monte Carlo method and that obtained with the density-matrix renormalization group method. We find that the energy is considerably improved in the FOWF. On the other hand, the physical quantities do not change significantly between these two wavefunctions as long as they describe the same phase, such as the paramagnetic phase. From these results, we not only demonstrate the improvement by the FOWF, but we also gain insights on the applicability and limitation of the GWF to the periodic Anderson model.

DOI
10.7566/jpsj.86.013701

Related papers

Wave Functions of Correlated-Electron State in the Periodic Anderson Model

Takashi Yanagisawa

We investigate the periodic Anderson model by a variational Monte Carlo method. A new wave function is proposed based on the well-known Gutzwiller wave function. We calculate the physical properties using the Monte Carlo algorithm to consider the strong correlation effects. The ground state energy is improved considerably by taking into account electron-hole excitation effects.

Journal of the Physical Society of Japan•Mar 15, 1999•1 citation

Periodic Anderson model with correlated conduction electrons: Variational and exact diagonalization study

Imre Hagymási, K. Itai, Jenő Sólyom

We investigate an extended version of the periodic Anderson model (the so-called periodic Anderson-Hubbard model) with the aim to understand the role of interaction between conduction electrons in the formation of the heavy-fermion and mixed-valence states. Two methods are used: (i) variational calculation with the Gutzwiller wave function optimizing numerically the ground-state energy and (ii) exact diagonalization of the Hamiltonian for short chains. The $f$-level occupancy and the renormalization factor of the quasiparticles are calculated as a function of the energy of the $f$ orbital for a wide range of the interaction parameters. The results obtained by the two methods are in reasonably good agreement for the periodic Anderson model. The agreement is maintained even when the interaction between band electrons, ${U}_{d}$, is taken into account, except for the half-filled case. This discrepancy can be explained by the difference between the physics of the one- and higher-dimensional models. We find that this interaction shifts and widens the energy range of the bare $f$ level, where heavy-fermion behavior can be observed. For large-enough ${U}_{d}$ this range may lie even above the bare conduction band. The Gutzwiller method indicates a robust transition from Kondo insulator to Mott insulator in the half-filled model, while ${U}_{d}$ enhances the quasiparticle mass when the filling is close to half filling.

arXiv:1106.4999•Physical Review B•Jun 11, 2012•15 citations

Hubbard physics in the symmetric half-filled periodic anderson-hubbard model

I. Hagymási, K. Itai, J. Sólyom

Two very different methods — exact diagonalization on finite chains and a variational method — are used to study the possibility of a metal-insulator transition in the symmetric half-filled periodic Anderson-Hubbard model. With this aim we calculate the density of doubly occupied d sites (gn d ) as a function of various parameters. In the absence of on-site Coulomb interaction (U f ) between f electrons, the two methods yield similar results. The double occupancy of d levels remains always finite just as in the one-dimensional Hubbard model. Exact diagonalization on finite chains gives the same result for finite U f , while the Gutzwiller method leads to a Brinkman-Rice transition at a critical value (U ), which depends on U f and V.

arXiv:1207.0381•Journal of the Korean Physical Society•May 1, 2013•0 citations

Hubbard physics in the symmetric half-filled periodic anderson-hubbard model

Imre Hagymási, Itai, Kazumasa, Sólyom, Jenő

Two very different methods – exact diagonalization on finite chains and a variational method – are used to study the possibility of a metal-insulator transition in the symmetric half-filled periodic Anderson-Hubbard model. With this aim we calculate the density of doubly occupied d sites (νd) as a function of various parameters. In the absence of on-site Coulomb interaction (Uf ) between f electrons, the two methods yield similar results. The double occupancy of d levels remains always finite just as in the one-dimensional Hubbard model. Exact diagonalization on finite chains gives the same result for finite Uf , while the Gutzwiller method leads to a Brinkman-Rice transition at a critical value (U c d ), which depends on Uf and V .

•Jan 1, 2013•0 citations

An efficient random-sampling method for calculating double occupancy of Gutzwiller wave function in single-band 1D and 2D lattices

Feng Zhang, Zhuo Ye, Yongxin Yao +2 more

We report a random-sampling method for computing the expectation value of physical quantities based on the Gutzwiler variational wave function. As the first application, we calculated the double occupancy, which is a critical quantity for understanding the correlation effects in many-body systems, for single-band 1D and 2D lattices. We demonstrated that the random-sampling scheme is more efficient than an existing Metropolis Monte-Carlo algorithm. For the 1D Hubbard model with only nearest-neighbour hopping, our results are almost identical to the exact analytic solution. We have also studied systems to which analytic solutions are not available, including the 1D lattices with next-nearest-neighbour hopping and 2D lattices. In addition, constraints on real-space configurations can be easily implemented in the current scheme to further improve the Gutzwiller wave function. As an example, we calculated the double occupancy for 1D Hubbard model by applying the constraint that all double-occupied sites are paired with an empty site. With enhanced correlation between double-occupied and empty sites, the constraint results in much improved ground-state energy for 1D Hubbard model with strong on-site repulsion.

Molecular Physics•Aug 28, 2020•3 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 propose a variational wavefunction containing parameters to tune the probabilities of all the possible onsite configurations for the periodic Anderson model. We call it the full onsite-correlation wavefunction (FOWF). This is a simple extension of the Gutzwiller wavefunction (GWF), in which one parameter is included to tune the double occupancy of the f electrons at the same site. We compare the energy of the GWF and the FOWF evaluated by the variational Monte Carlo method and that obtained with the density-matrix renormalization group method. We find that the energy is considerably improved in the FOWF. On the other hand, the physical quantities do not change significantly between these two wavefunctions as long as they describe the same phase, such as the paramagnetic phase. From these results, we not only demonstrate the improvement by the FOWF, but we also gain insights on the applicability and limitation of the GWF to the periodic Anderson model.
We investigate an extended version of the periodic Anderson model (the so-called periodic Anderson-Hubbard model) with the aim to understand the role of interaction between conduction electrons in the formation of the heavy-fermion and mixed-valence states. Two methods are used: (i) variational calculation with the Gutzwiller wave function optimizing numerically the ground-state energy and (ii) exact diagonalization of the Hamiltonian for short chains. The fff-level occupancy and the renormalization factor of the quasiparticles are calculated as a function of the energy of the fff orbital for a wide range of the interaction parameters. The results obtained by the two methods are in reasonably good agreement for the periodic Anderson model. The agreement is maintained even when the interaction between band electrons, Ud{U}_{d}Ud​, is taken into account, except for the half-filled case. This discrepancy can be explained by the difference between the physics of the one- and higher-dimensional models. We find that this interaction shifts and widens the energy range of the bare fff level, where heavy-fermion behavior can be observed. For large-enough Ud{U}_{d}Ud​ this range may lie even above the bare conduction band. The Gutzwiller method indicates a robust transition from Kondo insulator to Mott insulator in the half-filled model, while Ud{U}_{d}Ud​ enhances the quasiparticle mass when the filling is close to half filling.
Two very different methods — exact diagonalization on finite chains and a variational method — are used to study the possibility of a metal-insulator transition in the symmetric half-filled periodic Anderson-Hubbard model. With this aim we calculate the density of doubly occupied d sites (gn d ) as a function of various parameters. In the absence of on-site Coulomb interaction (U f ) between f electrons, the two methods yield similar results. The double occupancy of d levels remains always finite just as in the one-dimensional Hubbard model. Exact diagonalization on finite chains gives the same result for finite U f , while the Gutzwiller method leads to a Brinkman-Rice transition at a critical value (U ), which depends on U f and V.
Two very different methods – exact diagonalization on finite chains and a variational method – are used to study the possibility of a metal-insulator transition in the symmetric half-filled periodic Anderson-Hubbard model. With this aim we calculate the density of doubly occupied d sites (νd) as a function of various parameters. In the absence of on-site Coulomb interaction (Uf ) between f electrons, the two methods yield similar results. The double occupancy of d levels remains always finite just as in the one-dimensional Hubbard model. Exact diagonalization on finite chains gives the same result for finite Uf , while the Gutzwiller method leads to a Brinkman-Rice transition at a critical value (U c d ), which depends on Uf and V .
We report a random-sampling method for computing the expectation value of physical quantities based on the Gutzwiler variational wave function. As the first application, we calculated the double occupancy, which is a critical quantity for understanding the correlation effects in many-body systems, for single-band 1D and 2D lattices. We demonstrated that the random-sampling scheme is more efficient than an existing Metropolis Monte-Carlo algorithm. For the 1D Hubbard model with only nearest-neighbour hopping, our results are almost identical to the exact analytic solution. We have also studied systems to which analytic solutions are not available, including the 1D lattices with next-nearest-neighbour hopping and 2D lattices. In addition, constraints on real-space configurations can be easily implemented in the current scheme to further improve the Gutzwiller wave function. As an example, we calculated the double occupancy for 1D Hubbard model by applying the constraint that all double-occupied sites are paired with an empty site. With enhanced correlation between double-occupied and empty sites, the constraint results in much improved ground-state energy for 1D Hubbard model with strong on-site repulsion.
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.