Mass classification of dark matter perturbers of stellar tidal streams
  Stellar streams formed by tidal stripping of progenitors orbiting around the
Milky Way are expected to be perturbed by encounters with dark matter subhalos.
Recent studies have shown that they are an excellent proxy to infer properties
of the perturbers, such as their mass. Here we present two different
methodologies that make use of the fully non-Gaussian density distribution of
stellar streams: a Bayesian model selection based on the probability density
function (PDF) of stellar density, and a likelihood-free gradient boosting
classifier. While the schemes do not assume a specific dark matter model, we
are mainly interested in discerning the primordial black holes cold dark matter
(PBH CDM) hypothesis form the standard particle dark matter one. Therefore, as
an application we forecast model selection strength of evidence for cold dark
matter clusters of masses $10^3$ - $10^5 M_{\odot}$ and $10^5$ - $10^9
M_{\odot}$, based on a GD-1-like stellar stream and including realistic
observational errors. Evidence for the smaller mass range, so far
under-explored, is particularly interesting for PBH CDM. We expect weak to
strong evidence for model selection based on the PDF analysis, depending on the
fiducial model. Instead, the gradient boosting model is a highly efficient
classifier (99\% accuracy) for all mass ranges here considered. As a further
test of the robustness of the method, we reach similar conclusions when
performing forecasts further dividing the largest mass range into $10^5$ -
$10^7 M_{\odot}$ and $10^7$ - $10^9 M_{\odot}$ ranges.
