  We present a theoretical study of the density and spin (representing the two
components) linear response of Fermi superfluids with tunable attractive
interactions and population imbalance. In both linear response theories, we
find that the fluctuations of the order parameter must be treated on equal
footing with the gauge transformations associated with the symmetries of the
Hamiltonian so that important constraints including various sum rules can be
satisfied. Both theories can be applied to the whole BCS-Bose-Einstein
condensation crossover. The spin linear responses are qualitatively different
with and without population imbalance because collective-mode effects from the
fluctuations of the order parameter survive in the presence of population
imbalance, even though the associated symmetry is not broken by the order
parameter. Since a polarized superfluid becomes unstable at low temperatures in
the weak and intermediate coupling regimes, we found that the density and spin
susceptibilities diverge as the system approaches the unstable regime, but the
emergence of phase separation preempts the divergence.
