We provide a semi-analytic study of the small scale aspects of the power
spectra of warm dark matter (WDM) candidates that decoupled while relativistic
with arbitrary distribution functions. These are characterized by two widely
different scales keqββΌ0.01(Mpc)β1 and k_{fs}=
\sqrt{3}\,k_{eq}/2\,^{1/2} with 1/2βͺ1 the
velocity dispersion at matter radiation equality. Density perturbations evolve
through three stages: radiation domination when the particle is relativistic
and non-relativistic and matter domination. An early ISW effect during the
first stage leads to an enhancement of density perturbations and a plateau in
the transfer function for kβ²kfsβ. An effective fluid description
emerges at small scales which includes the effects of free streaming in initial
conditions and inhomogeneities. The transfer function features
\emph{WDM-acoustic oscillations} at scales kβ³2kfsβ. We study the
power spectra for two models of sterile neutrinos with mβΌkeV
produced non-resonantly, at the QCD and EW scales respectively. The latter case
yields acoustic oscillations on mass scales βΌ108Mββ. Our
results reveal a \emph{quasi-degeneracy} between the mass, distribution
function and decoupling temperature suggesting caveats on the constraints on
the mass of a sterile neutrino from current WDM N-body simulations and
Lyman-Ξ± forest data. A simple analytic interpolation of the power
spectra between large and small scales and its numerical implementation is
given.Comment: 47 pages, 17 figures, section with comparison with Boltzmann code