Physical review. B, Condensed matter
Finite-Uslave-boson methods and the intersite magnetic correlations in the Hubbard models
Hiroyuki Kaga
Abstract
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/${\mathit{N}}^{2}$) in the 1/N expansion of the functional-integral formalism. The obtained effective exchange coupling constant ${\mathit{J}}_{\mathrm{eff}}$(\ensuremath{\delta})=J(2${\mathit{p}}_{0}^{2}$-\ensuremath{\delta}/2) with J\ensuremath{\equiv}${\mathit{t}}_{\mathit{i}\mathit{j}}^{2}$/U is strongly dependent upon doping \ensuremath{\delta} and is scaled with the reduced Coulomb repulsion u\ensuremath{\equiv}U/${\mathit{U}}_{\mathit{c}}$ in the two-dimensional (2D) and 3D Hubbard models, where ${\mathit{U}}_{\mathit{c}}$ 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.