We present methods to assess whether gamma-ray excesses towards Milky Way dwarf galaxies can be attributed to astrophysical sources rather than to dark matter annihilation. As a case study we focus on Reticulum II, the dwarf which shows the strongest evidence for a gamma-ray signal in Fermi data. Dark matter models and those with curved energy spectra provide good fits to the data, while a simple power law is ruled out at 97.5% confidence. We compare RetII's spectrum to known classes of gamma-ray sources and find a useful representation in terms of spectral curvature and the energy at which the spectral energy distribution peaks. In this space the blazar classes appear segregated from the confidence region occupied by RetII. Pulsars have similar gamma-ray spectra to RetII but we show that RetII is unlikely to host a pulsar population detectable in gamma rays. Tensions with astrophysical explanations are stronger when analyzing 6.5 years of Pass 7 than with the same amount of Pass 8 data, where the excess is less significant. These methods are applicable to any dwarf galaxy which is a promising dark matter target and shows signs of gamma-ray emission along its line of sight.

The expected gamma-ray flux coming from dark matter annihilation in dwarf spheroidal (dSph) galaxies depends on the so-called `J-factor', the integral of the squared dark matter density along the line-of-sight. We examine the degree to which estimates of J are sensitive to contamination (by foreground Milky Way stars and stellar streams) of the stellar-kinematic samples that are used to infer dark matter densities in `ultrafaint' dSphs. Applying standard kinematic analyses to hundreds of mock data sets that include varying levels of contamination, we find that mis-classified contaminants can cause J-factors to be overestimated by orders of magnitude. Stellar-kinematic data sets for which we obtain such biased estimates tend 1) to include relatively large fractions of stars with ambiguous membership status, and 2) to give estimates for J that are sensitive to specific choices about how to weight and/or to exclude stars with ambiguous status. Comparing publicly-available stellar-kinematic samples for the nearby dSphs Reticulum~II and Segue~I, we find that only the latter displays both of these characteristics. Estimates of Segue~I's J-factor should therefore be regarded with a larger degree of caution when planning and interpreting gamma-ray observations. Moreover, robust interpretations regarding dark matter annihilation in dSph galaxies in general will require explicit examination of how interlopers might affect the inferred dark matter density profile.

We present results from spectroscopic observations with the Michigan/Magellan Fiber System (M2FS) of $147$ stellar targets along the line of sight to the newly-discovered `ultrafaint' stellar systems Tucana 2 (Tuc 2) and Grus 1 (Gru 1). Based on simultaneous estimates of line-of-sight velocity and stellar-atmospheric parameters, we identify 8 and 7 stars as probable members of Tuc 2 and and Gru 1, respectively. Our sample for Tuc 2 is sufficient to resolve an internal velocity dispersion of $8.6_{-2.7}^{+4.4}$ km s$^{-1}$ about a mean of $-129.1_{-3.5}^{+3.5}$ km s$^{-1}$ (solar rest frame), and to estimate a mean metallicity of [Fe/H]= $-2.23_{-0.12}^{+0.18}$. These results place Tuc 2 on chemodynamical scaling relations followed by dwarf galaxies, suggesting a dominant dark matter component with dynamical mass $2.7_{-1.3}^{+3.1}\times 10^6$ $\mathrm{M}_{\odot}$ enclosed within the central $\sim 160$ pc, and dynamical mass-to-light ratio $1900_{-900}^{+2200}$ $\mathrm{M}_{\odot}/L_{V,\odot}$. For Gru 1 we estimate a mean velocity of $-140.5_{-1.6}^{+2.4}$ km s$^{-1}$ and a mean metallicity of [Fe/H]=$-1.42_{-0.42}^{+0.55}$, but our sample does not resolve Gru 1's velocity dispersion. The radial coordinates of Tuc 2 and Gru 1 in Galactic phase space suggest that their orbits are among the most energetic within distance $\leq 300$ kpc. Moreover, their proximity to each other in this space arises naturally if both objects are trailing the Large Magellanic Cloud.

We present an update of the CLUMPY code for the calculation of the astrophysical J-factors (from dark matter annihilation/decay) for any Galactic or extragalactic dark matter halo including substructures: halo-to-halo concentration scatter may now be enabled, boost factors can include several levels of substructures, and triaxiality is a new option for dark matter haloes. This new version takes advantage of the CFITSIO and HEALPix libraries to propose FITS output maps using the HEALPix pixelisation scheme. Skymaps for $\gamma$-ray and $\nu$ signals from generic annihilation/decay spectra are now direct outputs of CLUMPY. Making use of HEALPix routines, smoothing by a user-defined instrumental Gaussian beam and computing the angular power spectrum of the maps is now possible. In addition to these improvements, the main novelty is the implementation of a Jeans analysis module, to obtain dark matter density profiles from kinematic data in relaxed spherical systems (e.g., dwarf spheroidal galaxies). The code is also interfaced with the GreAT toolkit designed for Markov Chain Monte Carlo analyses, from which probability density functions and credible intervals can be obtained for velocity dispersions, dark matter profiles, and J-factors.

The dwarf spheroidal galaxies (dSph) of the Milky Way are among the most attractive targets for indirect searches of dark matter. In this work, we reconstruct the dark matter annihilation (J-factor) and decay profiles for the newly discovered dSph Reticulum II. Using an optimized spherical Jeans analysis of kinematic data obtained from the Michigan/Magellan Fiber System (M2FS), we find Reticulum II's J-factor to be among the largest of any Milky Way dSph. We have checked the robustness of this result against several ingredients of the analysis. Unless it suffers from tidal disruption or significant inflation of its velocity dispersion from binary stars, Reticulum II may provide a unique window on dark matter particle properties.

Dwarf spheroidal (dSph) galaxies are prime targets for present and future gamma-ray telescopes hunting for indirect signals of particle dark matter. The interpretation of the data requires careful assessment of their dark matter content in order to derive robust constraints on candidate relic particles. Here, we use an optimised spherical Jeans analysis to reconstruct the `astrophysical factor' for both annihilating and decaying dark matter in 21 known dSphs. Improvements with respect to previous works are: (i) the use of more flexible luminosity and anisotropy profiles to minimise biases, (ii) the use of weak priors tailored on extensive sets of contamination-free mock data to improve the confidence intervals, (iii) systematic cross-checks of binned and unbinned analyses on mock and real data, and (iv) the use of mock data including stellar contamination to test the impact on reconstructed signals. Our analysis provides updated values for the dark matter content of 8 `classical' and 13 `ultrafaint' dSphs, with the quoted uncertainties directly linked to the sample size; the more flexible parametrisation we use results in changes compared to previous calculations. This translates into our ranking of potentially-brightest and most robust targets---viz., Ursa Minor, Draco, Sculptor---, and of the more promising, but uncertain targets---viz., Ursa Major 2, Coma---for annihilating dark matter. Our analysis of Segue 1 is extremely sensitive to whether we include or exclude a few marginal member stars, making this target one of the most uncertain. Our analysis illustrates challenges that will need to be addressed when inferring the dark matter content of new `ultrafaint' satellites that are beginning to be discovered in southern sky surveys.

Dwarf spheroidal (dSph) galaxies are among the most promising targets for the indirect detection of dark matter (DM) from annihilation and/or decay products. Empirical estimates of their DM content - and hence the magnitudes of expected signals - rely on inferences from stellar-kinematic data. However, various kinematic analyses can give different results and it is not obvious which are most reliable. Using extensive sets of mock data of various sizes (mimicking 'ultra-faint' and 'classical' dSphs) and an MCMC engine, here we investigate biases, uncertainties, and limitations of analyses based on parametric solutions to the spherical Jeans equation. For a variety of functional forms for the tracer and DM density profiles, as well as the orbital anisotropy profile, we examine reliability of estimates for the astrophysical J- and D-factors for annihilation and decay, respectively. For large (N > 1000) stellar-kinematic samples typical of 'classical' dSphs, errors tend to be dominated by systematics, which can be reduced through the use of sufficiently general and flexible functional forms. For small (N < 100) samples typical of 'ultrafaints', statistical uncertainties tend to dominate systematic errors and flexible models are less necessary. We define an optimal strategy that would mitigate sensitivity to priors and other aspects of analyses based on the spherical Jeans equation. We also find that the assumption of spherical symmetry can bias estimates of J (with the 95% credibility intervals not encompassing the true J-factor) when the object is mildly triaxial (axis ratios b/a = 0.8, c/a = 0.6). A concluding table summarises the typical error budget and biases for the different sample sizes considered.