Stability of a fully magnetized ferromagnetic state in repulsively interacting ultracold Fermi gases

We construct a variational wave function to study whether a fully polarized Fermi sea of ultracold atoms is energetically stable against a single spin flip. Our variational wave function contains short-range correlations at least to the same level as Gutzwiller's projected wave function. For the Hubbard lattice model and the continuum model with pure repulsive interaction, we show that a fully polarized Fermi sea is generally unstable even for infinite repulsive strength. By contrast, for a resonance model, the ferromagnetic state is possible if the s-wave scattering length is positive and sufficiently large and the system is prepared to be orthogonal to the molecular bound state. However, we cannot rule out the possibility that more exotic correlations can destabilize the ferromagnetic state.