Abstract We discuss a variational ground state wave function for the symmetric periodic Anderson model with commensurate spin density order. The energy of this ansatz is evaluated in closed form. Our approach generalizes a variational treatment proposed recently by Strack and Vollhardt and results in significantly lower energy. Contrary to the Gutzwiller ansatz the wave function recovers the Schrieffer‐Wolff limit for large Coulomb repulsion. We clarify the relation of our approach to unrestricted Hartree‐Fock and present a comparison with existing quantum Monte Carlo calculations for one dimension.
The ground-state energy, hybridization energy, momentum distributions, local magnetic moments, and magnetic phase diagram of the symmetric periodic Anderson model are calculated using a correlated wave function in d=1,2, and 3 dimensions. The wave function is obtained by applying a Gutzwiller projector on a recently proposed single-particle product state with spin-density wave order. Evaluations are performed with an approximate scheme that becomes exact in high dimensions. Comparison with Monte Carlo data by Blankenbecler et al. for d=1 shows that the results are remarkably accurate for all values of the coupling. Hence the results for d=2,3 can be expected to be even more accurate and may serve as benchmarks for future numerical work in d>1.
We investigate Gutzwiller-correlated wave functions as variational ground states for the periodic Anderson model. We extend a recently developed technique for the Hubbard model to the case of two hybridized bands with degeneracy N. We calculate expectation values in the following limits: dilute system (single-impurity case), large degeneracy N, and large dimension d. The mean-field results using the slave-boson approach of Coleman and of Kotliar and Ruckenstein are reproduced in the large-N and large-d limit, respectively. These techniques, therefore, obey variational principles at T=0 in the respective limits and their applicability for T>0 is naturally limited to the Fermi-liquid regime. Evaluations of Gutzwiller-correlated wave functions using the ``Gutzwiller approximation'' are seen to become exact in the limit of large dimensions. A systematic expansion in 1/d can be performed for N=2. In general, corrections are small, of order 1/(dN).
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,
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
f level, where heavy-fermion behavior can be observed. For large-enough
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 enhances the quasiparticle mass when the filling is close to half filling.