Indirect detection strategies of particle Dark Matter (DM) in Dwarf spheroidal satellite galaxies (dSphs) typically entail searching for annihilation signals above the astrophysical background. To robustly compare model predictions with the observed fluxes of product particles, most analyses of astrophysical data -- which are generally frequentist -- rely on estimating the abundance of DM by calculating the so-called $\textit{J-factor}$. This quantity is usually inferred from the kinematic properties of the stellar population of a dSph using Jeans equation, commonly by means of Bayesian techniques which entail the presence (and additional systematic uncertainty) of prior choice. Here, extending earlier work, we develop a scheme to derive the profile likelihood for $J$-factors of dwarf spheroidals for models with five or more free parameters. We validate our method on a publicly available simulation suite, released by the Gaia Challenge, finding satisfactory statistical properties for coverage and bias. We present the profile likelihood function and maximum likelihood estimates for the $J$-factor of ten dSphs . As an illustration, we apply these profiles likelihood to recently published analyses of gamma-ray data with the Fermi Large Area Telescope to derive new, consistent upper limits on the DM annihilation cross-section. We do this for a subset of systems, generally referred to as $\textit{classical dwarfs}$. The implications of these findings for DM searches are discussed, together with future improvements and extensions of this technique.

Dwarf spheroidal galaxies are among the most promising targets for indirect dark matter (DM) searches in $\gamma$-rays. The $\gamma$-ray flux from DM annihilation in a dwarf spheroidal galaxy is proportional to the $J$-factor of the source. The $J$-factor of a dwarf spheroidal galaxy is the line-of-sight integral of the DM mass density squared times $\langle \sigma_{\rm ann} v_{\rm rel} \rangle/(\sigma_{\rm ann} v_{\rm rel})_0$, where $\sigma_{\rm ann} v_{\rm rel}$ is the DM annihilation cross-section times relative velocity $v_{\rm rel}=|{\bf v}_{\rm rel}|$, angle brackets denote average over ${\bf v}_{\rm rel}$, and $(\sigma_{\rm ann} v_{\rm rel})_0$ is the $v_{\rm rel}$-independent part of $\sigma_{\rm ann} v_{\rm rel}$. If $\sigma_{\rm ann} v_{\rm rel}$ is constant in $v_{\rm rel}$, $J$-factors only depend on the DM space distribution in the source. However, if $\sigma_{\rm ann} v_{\rm rel}$ varies with $v_{\rm rel}$, as in the presence of DM self-interactions, $J$-factors also depend on the DM velocity distribution, and on the strength and range of the DM self-interaction. Models for self-interacting DM are increasingly important in the study of the small scale clustering of DM, and are compatible with current cosmological observations. Here we derive the $J$-factor of 20 dwarf spheroidal galaxies from stellar kinematic data under the assumption of Yukawa DM self-interactions. $J$-factors are derived through a profile Likelihood approach, assuming either NFW or cored DM profiles. We also compare our results with $J$-factors derived assuming the same velocity for all DM particles in the target galaxy. We find that this common approximation overestimates the $J$-factors by up to one order of magnitude. $J$-factors for a sample of DM particle masses, self-interaction coupling constants and density profiles are provided electronically, ready to be used in other projects.

Line-of-sight integrals of the squared density, commonly called the J-factor, are essential for inferring dark matter annihilation signals. The J-factors of dark matter-dominated dwarf spheroidal satellite galaxies (dSphs) have typically been derived using Bayesian techniques, which for small data samples implies that a choice of priors constitutes a non-negligible systematic uncertainty. Here we report the development of a new fully frequentist approach to construct the profile likelihood of the J-factor. Using stellar kinematic data from several classical and ultra-faint dSphs, we derive the maximum likelihood value for the J-factor and its confidence intervals. We validate this method, in particular its bias and coverage, using simulated data from the Gaia Challenge. We find that the method possesses good statistical properties. The J-factors and their uncertainties are generally in good agreement with the Bayesian-derived values, with the largest deviations restricted to the systems with the smallest kinematic datasets. We discuss improvements, extensions, and future applications of this technique.