We provide a general framework for self-interacting warm dark matter (WDM) in cosmological perturbations, by deriving from first principles a Boltzmann hierarchy which retains certain independence from a particular interaction Lagrangian. We consider elastic interactions among the massive particles, and obtain a hierarchy which is more general than the ones usually obtained for non-relativistic (as for cold DM) or for ultra-relativistic (as for neutrinos) approximations. The more general momentum-dependent kernel integrals in the Boltzmann collision terms, are explicitly calculated for different field-mediator models, including examples of a scalar field (either massive or massless) or a massive vector field. As an application, we study the evolution of the interaction rate per particle under the relaxation time approximation, and assess when a given self-interaction is relevant in comparison with the Hubble expansion rate. Our framework aims to be a useful starting point to evaluate DM self-interaction effects in the linear power spectrum, necessary to then study its evolution all the way to non-linear stages of structure formation, where certain DM interactions were proven to be relevant.

We calculate the most stringent constraints up to date on the parameter space for sterile neutrino warm dark matter models possessing a radiative decay channel into X-rays. These constraints arise from the X-ray flux observations from the Galactic center (central parsec), taken by the XMM and NuSTAR missions. We compare the results obtained from using different dark matter density profiles for the Milky Way, such as NFW, Burkert or Einasto, to that produced by the Ruffini-Arg\"uelles-Rueda (RAR) fermionic model, which has the distinct feature of depending on the particle mass. We show that due to the novel core-halo morphology present in the RAR profile, the allowed particle mass window is narrowed down to $m_s\sim 10-15$ keV, when analyzed within the $\nu$MSM sterile neutrino model. We further discuss on the possible effects in the sterile neutrino parameter-space bounds due to a self-interacting nature of the dark matter candidates.