Breakdown of the Nagaoka phase at finite doping

The Nagaoka $(U=\ensuremath{\infty})$ limit of the Hubbard model on a square lattice is mapped onto the itinerant-localized Kondo model at infinitely strong coupling. Such a model is well suited to perform quantum Monte Carlo (QMC) simulations to compute spin correlation functions. For periodic boundary conditions, this model is shown to exhibit no short-range ferromagnetic (FM) spin correlations at any doping $\ensuremath{\delta}\ensuremath{\ge}0.01$ and at finite temperature $T=0.1t.$ Our simulations give no indication that there is a tendency towards ferromagnetic ordering in the ground state, with more than one hole. Employing on the other hand the open boundary conditions (or mixed boundary conditions) may result in the qualitatively different results for the thermodynamic limit depending on the way one chooses to approach this limit. These observations imply that the relevant thermodynamic limit remains unclear.