We study the Gutzwiller method for the spinless fermion model in one dimension, which is one of the simplest models that incorporates the intersite Coulomb interaction. The Gutzwiller solution of this model has been studied in the literature but with differing results. We obtain the Gutzwiller solution of the problem by a careful enumeration of the many-particle configurations and explain the origin of the discrepancy in the literature to be due to the neglect of the correlation between the neighboring bond occupancies. The correct implementation of the Gutzwiller approach yields results different from the slave-boson mean-field theory, unlike for the Hubbard model with on-site interaction, where both methods are known to be equivalent. The slave-boson and the Gutzwiller results are compared to the exact solution, available for the half-filled case, and to the numerical exact diagonalization results for the case of general filling.
We examine the two existing slave-boson methods for the finite-U Hubbard model by Kotliar and Ruckenstein (KR) and by Newns et al. and Zou and Anderson (NZA) in the light of the Fermi-liquid quasiparticle picture. It is shown that to treat doubly occupied states the KR slave-boson scheme is more consistent with this picture and that the well-known erroneous gauge-symmetry breaking encountered in the mean-field treatment of the NZA double-occupancy boson is due to the presence of the hidden double-fermion operators. Then, the effective intersite antiferromagnetic interaction between quasiparticles is derived in the KR slave-boson approach by taking into account empty and double-occupancy slave-boson fluctuations to O(1/
N2) in the 1/N expansion of the functional-integral formalism. The obtained effective exchange coupling constant
Jeff(\ensuremath{\delta})=J(2
p02-\ensuremath{\delta}/2) with J\ensuremath{\equiv}
tij2/U is strongly dependent upon doping \ensuremath{\delta} and is scaled with the reduced Coulomb repulsion u\ensuremath{\equiv}U/
Uc in the two-dimensional (2D) and 3D Hubbard models, where
Uc is the critical value for the metal-to-insulator transition. The roles of this intersite correlation and the spin-dependent band-renormalization correlation included in the mean-field theory for the paramagnetic-to-antiferromagnetic instability are compared in the two Hubbard models. The intersite contribution is more significant in the 2D case than in the 3D case, but is generally small compared with the band effect except for small dopings \ensuremath{\delta} and smaller U values. The maximum doping ${\mathrm{\ensuremath{\delta}}}_{\mathrm{cr}}$ for the paramagnetic-to-ferromagnetic instability in the strong-coupling limit (U\ensuremath{\rightarrow}\ensuremath{\infty}) of the 2D and 3D Hubbard models is found to be ${\mathrm{\ensuremath{\delta}}}_{\mathrm{cr}}$\ensuremath{\sim}0.28 and ${\mathrm{\ensuremath{\delta}}}_{\mathrm{cr}}$\ensuremath{\sim}0.38, respectively.
In the first part of the paper, we study the stability of antiferromagnetic (AF), charge density wave (CDW), and superconducting (SC) states within the t-J-U-V model of strongly correlated electrons by using the statistically consistent Gutzwiller approximation (SGA). We concentrate on the role of the intersite Coulomb interaction term V in stabilizing the CDW phase. In particular, we show that the charge ordering appears only above a critical value of V in a limited hole-doping range δ. The effect of the V term on SC and AF phases is that a strong interaction suppresses SC, whereas the AF order is not significantly influenced by its presence. In the second part, separate calculations for the case of a pure SC phase have been carried out within an extended approach (the diagrammatic expansion for the Gutzwiller wave function, DE-GWF) in order to analyze the influence of the intersite Coulomb repulsion on the SC phase with the higher-order corrections included beyond the SGA method. The upper concentration for the SC disappearance decreases with increasing V, bringing the results closer to experiment. In appendices A and B we discuss the ambiguity connected with the choice of the Gutzwiller renormalization factors within the renormalized mean filed theory when either AF or CDW orders are considered. At the end, we overview briefly the possible extensions of the current models to put descriptions of the SC, AF, and CDW states on equal footing.