The ground-state energy of the two-sublattice two-dimensional Hubbard model is calculated in the static-fluctuation approximation with allowance for electron transfer from sites to next-to-nearest neighbor sites. In a specific case, the energy of the one-dimensional Hubbard model is calculated and compared with an exact solution.
We have calculated the ground-state energy and the fluctuation contribution to the free-energy for the Hubbard model at half filling. We have used a gauge that leads to the Hartree-Fock solution in the saddle-point approximation. The mean-field ground-state energy calculated here agrees with the results obtained by Dagotto et al. [Phys. Rev. B 45, 10 741 (1992)] with use of the exact-diagonalization techniques on a 4\ifmmode\times\else\texttimes\fi{}4 site cluster. In the large-U limit, the ground-state energy per site is obtained as Eg/N=-U/4-J+O(J2/U). The fluctuation contribution to the free-energy was found to have a stronger temperature and U dependence than the mean-field contribution, and to play a dominant role for strongly correlated systems.
The problem of the ground state of the electronic system in the Hubbard model for U=∞ is discussed. The author investigates the normal (singlet or nonmagnetic) N state of the electronic system over the entire range of electron densities n⩽1. It is shown that the energy of the N state ɛ 0 (1) (n) in a one-particle approximation, such as (e.g.) the extended Hartree-Fock approximation, is lower than the energy of the saturated ferromagnetic FM state ɛ FM(n) for all n. The dynamic magnetic susceptibility is calculated in the random phase approximation, and it is shown that the N state is stable over the entire range of electron densities: The static susceptibility (ω=0) does not have a band singularity in the zero-wave vector limit q→0. A formally exact representation is obtained for the mass operator of the one-particle Green’s function, and an approximation of this operator is proposed: M k(E)⋍λF(E), where λ=n(1−n)/(1−n/2)z is the kinematic interaction parameter, z is the number of nearest neighbors, and F(E) is the total single-site Green’s function. For an elliptical density of states the integral equation for F(E) is solved exactly, ad it is shown that the spectral intensity rigorously satisfies the sum rule. The calculated energy of the strongly correlated N state ɛ 0(n)<ɛ FM(n) for all n, and in light of this relationship the author discusses the hypothesis that the ground state of the system is the normal (singlet) state in the thermodynamic limit. The electron distribution function at T=0 differs significantly from the Fermi step; it is “smeared” along the entire energy spectrum, and discontinuities do not occur in the region of the chemical potential m.
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.