The nature of dark matter remains uncertain despite several decades of dedicated experimental searches. The lack of tangible evidence for its non-gravitational interactions with ordinary matter gives good motivation for exploring new avenues of inferring its properties through purely gravitational probes. In particular, addressing its small-scale distribution could provide valuable new insights into its particle nature, either confirming the predictions of cold dark matter hypothesis or favouring models with suppressed small-scale matter power spectrum. In this work a machine learning technique for constraining the abundance of DM subhalos through the analysis of galactic stellar field is proposed. While the study is the first of its kind and hence applied only to a simplified synthetic datasets, the obtained results show promising potential for addressing the amount of DM substructure present within Milky Way. Using accurate astrometric observations, which became available only recently and are expected to rapidly improve in the near future, there is good hope to reach sensitivity needed for detecting DM subhalos with masses down to $10^7 M_\odot$.

Milky Way dwarf spheroidal satellites are a prime target for Dark Matter (DM) indirect searches. There have been recent reassessments of the expected DM gamma-ray signals in case of long-range interactions, commonly known as Sommerfeld enhancement. Since details of the underlying DM phase-space distribution function become critical, there are potentially large uncertainties in the final result. We provide here a first attempt towards a comprehensive investigation of these systematics, addressing the impact on the expected DM flux from Milky Way dwarfs via Bayesian inference on the available stellar kinematic datasets. After reconsidering the study case of ergodic systems, we investigate for the first time scenarios where DM particle orbits may have a radial or tangential bias. We consider both cuspy and cored parametric DM density profiles, together with the case of a non-parametric halo modelling directly connected to line-of-sight observable quantities. The main findings of our work highlight the relevance of the assumed phase-space distribution: Referring to a generalized J-factor, namely the line-of-sight convolution of the spatial part in case of velocity-dependent annihilation rate, an enhancement (suppression) with respect to the limit of isotropic phase-space distributions is obtained for the case of tangentially (radially) biased DM particle orbits. We provide new estimates for J-factors for the eight brightest Milky Way dwarfs also in the limit of velocity-independent DM annihilation, in good agreement with previous results in literature, and derive data-driven lower-bounds based on the non-parametric modelling of the halo density. The outcome of our broad study stands out as a representative of the state-of-the-art in the field, and falls within the interest of current and future experimental collaborations involved in DM indirect detection programs.

We address the problem of reconstructing the phase-space distribution function for an extended collisionless system, with known density profile and in equilibrium within an axisymmetric gravitational potential. Assuming that it depends on only two integrals of motion, namely the energy and the component of the angular momentum along the axis of symmetry $L_z$, there is a one-to-one correspondence between the density profile and the component of the distribution function that is even in $L_z$, as well as between the weighted azimuthal velocity profile and the odd component. This inversion procedure was originally proposed by Lynden-Bell and later refined in its numerical implementation by Hunter & Qian; after overcoming a technical difficulty, we apply it here for the first time in presence of a strongly flattened component, as a novel approach of extracting the phase-space distribution function for dark matter particles in the halo of spiral galaxies. We compare results obtained for realistic axisymmetric models to those in the spherical symmetric limit as assumed in previous analyses, showing the rather severe shortcomings in the latter. We then apply the scheme to the Milky Way and discuss the implications for the direct dark matter searches. In particular, we reinterpret the null results of the Xenon1T experiment for spin-(in)dependent interactions and make predictions for the annual modulation of the signal for a set of axisymmetric models, including a self-consistently defined co-rotating halo.