We derive an expression for the superfluid density of a uniform two-component Fermi gas through the BCS-BEC crossover in terms of the thermodynamic potential in the presence of an imposed superfluid flow. Treating the pairing fluctuations in a Gaussian approximation following the approach of Nozi\`eres and Schmitt-Rink, we use this definition of ${\ensuremath{\rho}}_{s}$ to obtain an explicit result which is valid at finite temperatures and over the full BCS-BEC crossover. It is crucial that the BCS gap $\ensuremath{\Delta}$, the chemical potential $\ensuremath{\mu}$, and ${\ensuremath{\rho}}_{s}$ all include the effect of fluctuations at the same level in a self-consistent manner. We show that the normal fluid density ${\ensuremath{\rho}}_{n}\ensuremath{\equiv}n\ensuremath{-}{\ensuremath{\rho}}_{s}$ naturally separates into a sum of contributions from Fermi BCS quasiparticles $({\ensuremath{\rho}}_{n}^{F})$ and Bose collective modes $({\ensuremath{\rho}}_{n}^{B})$. The expression for ${\ensuremath{\rho}}_{n}^{F}$ is just Landau's formula for a BCS Fermi superfluid but now calculated over the BCS-BEC crossover. The expression for the Bose contribution ${\ensuremath{\rho}}_{n}^{B}$ is more complicated and only reduces to Landau's formula for a Bose superfluid in the extreme BEC limit, where all the fermions have formed stable Bose pairs and the Bogoliubov excitations of the associated molecular Bose condensate are undamped. In a companion paper, we present numerical calculations of ${\ensuremath{\rho}}_{s}$ using an expression equivalent to the one derived in this paper, over the BCS-BEC crossover, including unitarity, and at finite temperatures.
The superfluid density is a fundamental quantity describing the response to a rotation as well as in two-fluid collisional hydrodynamics. We present extensive calculations of the superfluid density ${\ensuremath{\rho}}_{s}$ in the BCS-BEC crossover regime of a uniform superfluid Fermi gas at finite temperatures. We include strong-coupling or fluctuation effects on these quantities within a Gaussian approximation. We also incorporate the same fluctuation effects into the BCS single-particle excitations described by the superfluid order parameter $\ensuremath{\Delta}$ and Fermi chemical potential $\ensuremath{\mu}$, using the Nozi\`eres--Schmitt-Rink approximation. This treatment is shown to be necessary for consistent treatment of ${\ensuremath{\rho}}_{s}$ over the entire BCS-BEC crossover. We also calculate the condensate fraction
Nc as a function of the temperature, a quantity which is quite different from the superfluid density ${\ensuremath{\rho}}_{s}$. We show that the mean-field expression for the condensate fraction
Nc is a good approximation even in the strong-coupling BEC regime. Our numerical results show how ${\ensuremath{\rho}}_{s}$ and
Nc depend on temperature, from the weak-coupling BCS region to the BEC region of tightly bound Cooper pair molecules. In a companion paper [Phys. Rev. A 74, 063626 (2006)], we derive an equivalent expression for ${\ensuremath{\rho}}_{s}$ from the thermodynamic potential, which exhibits the role of the pairing fluctuations in a more explicit manner.