## DWARF GALAXY ANNIHILATION AND DECAY EMISSION PROFILES FOR DARK MATTER EXPERIMENTS

ALEX GERINGER-SAMETH 1 , SAVVAS M. KOUSHIAPPAS 2 , AND MATTHEW WALKER 3

Draft version March 12, 2015

## ABSTRACT

Gamma-ray searches for dark matter annihilation and decay in dwarf galaxies rely on an understanding of the dark matter density profiles of these systems. Conversely, uncertainties in these density profiles propagate into the derived particle physics limits as systematic errors. In this paper we quantify the expected dark matter signal from 20 Milky Way dwarfs using a uniform analysis of the most recent stellar-kinematic data available. Assuming that the observed stellar populations are equilibrium tracers of spherically-symmetric gravitational potentials that are dominated by dark matter, we find that current stellar-kinematic data can predict the amplitudes of annihilation signals to within a factor of a few for the ultra-faint dwarfs of greatest interest. On the other hand, the expected signal from several classical dwarfs (with high-quality observations of large numbers of member stars) can be localized to the ∼ 20% level. These results are important for designing maximally sensitive searches in current and future experiments using space and ground-based instruments.

Keywords: dark matter - galaxies: dwarf - galaxies: fundamental parameters - galaxies: kinematics and dynamics

## 1. INTRODUCTION

The search for cosmological dark matter annihilation or decay is a major effort in contemporary astrophysics. Educing the dark matter particle physics from observations requires a detailed understanding of the dark matter distribution in the systems under study. A productive avenue of approach has been to search for gamma-rays generated by dark matter annihilation in Milky Way dwarf spheroidal galaxies (e.g. Scott et al. 2010; Essig et al. 2010; Aleksi´ c et al. 2011; Geringer-Sameth &amp; Koushiappas 2011; Ackermann et al. 2011; Geringer-Sameth &amp; Koushiappas 2012; Aliu et al. 2012; Ackermann et al. 2014; Geringer-Sameth et al. 2014). Such systems are nearby, dark matter-dominated, and contain no conventional sources of astrophysical backgrounds (e.g. cosmic ray generation and propagation through interstellar gas). Many such dwarf galaxies have been discovered in recent years (Willman et al. 2005; Zucker et al. 2006b,a; Walsh et al. 2007; Belokurov et al. 2007, 2008; Belokurov et al. 2009, 2010) with the prospect of more discoveries from ongoing and future sky surveys like Pan-Starrs (Kaiser et al. 2002), the Vista Hemisphere Survey (Ashby et al. 2013, 2014), the Dark Energy Survey (Flaugher 2005), and eventually the Large Synoptic Survey Telescope (Tyson et al. 2003).

Previous studies of dwarf galaxies have begun to constrain the physical properties of dark matter (Geringer-Sameth &amp; Koushiappas 2011; Ackermann et al. 2011, 2014; GeringerSameth et al. 2014). The lack of any significant gamma-ray excess lead to the exclusion of generic dark matter candidates with annihilation cross sections on the order of the benchmark value for a thermal relic ( ∼ 3 × 10 -26 cm 3 s -1 ) and with masses less than a few tens of GeV. Despite current non-detections, dwarf galaxies-and their lack of astrophysical contaminating sources-offer the cleanest possible signature of dark matter

- 1 Department of Physics, Brown University, Providence, RI 02912 and McWilliams Center for Cosmology, Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213; alexgs@cmu.edu
- 2 Department of Physics, Brown University, Providence, RI 02912; koushiappas@brown.edu
- 3 McWilliams Center for Cosmology, Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213; mgwalker@andrew.cmu.edu

annihilation or decay compared with other targets. This is especially interesting in the context of recent claims of a Galactic center gamma-ray excess and associated dark matter interpretation (e.g. Hooper &amp; Goodenough 2011; Boyarsky et al. 2011; Abazajian &amp; Kaplinghat 2012, 2013; Abazajian et al. 2014; Daylan et al. 2014). Observations of dwarf galaxies have the potential to either confirm or rule out such an interpretation.

The dark matter distribution within a target system is a necessary ingredient for placing constraints on any particle theory that predicts dark matter annihilation or decay. Knowledge of the relative signal strengths amongst different targets as well as the spatial distribution of the emission is required for designing maximally sensitive searches in current and future experiments. The overall emission rate from annihilation is described by the ' J value', the integral along the line-ofsight and over an aperture of the square of the dark matter density. The amplitude of J helps to identify which dwarfs are the most promising for searches (i.e. are the 'brightest').

Different groups have devised various methods for estimating dark matter distributions (and uncertainties) using observations of line-of-sight velocities of dwarf galaxy member stars. Some authors use the kinematic data to fit for the mass and/or concentration of dark matter density profiles that are assumed to follow an analytic form typically used to describe low-mass 'subhalos' (virial mass ∼ 10 9 -10 M glyph[circledot] ) that form around Milky-Way-like galaxies in dissipationless cosmological simulations based on cold dark matter (e.g. Strigari et al. 2007; Martinez et al. 2009; Martinez 2013). Some studies make an explicit assumption of a cored profile (Cholis &amp; Salucci 2012; Salucci et al. 2012), while others take a more agnostic approach, fitting relatively flexible density profiles that are not restricted to the form used to describe simulated halos (e.g. Charbonnier et al. 2011).

In addition, different groups use different techniques for propagating the uncertainties in those dark matter distributions when incorporating gamma-ray non-detections to derive limits on the annihilation cross section as a function of particle mass. For example, in their joint analysis of stellarkinematic and gamma-ray data for several dwarfs, Ackermann

et al. (2011) take the uncertainty in J to be described by a lognormal distribution that is subsequently folded into a likelihood for a given particle physics model. Working with similar data, Geringer-Sameth &amp; Koushiappas (2011) separated that same log-normal distribution for J from statistical uncertainties in the gamma-ray data, allowing this systematic J uncertainty to be compared with the derived cross section limits. The results shown by Geringer-Sameth &amp; Koushiappas (2011) are alarming: the uncertainty in J can affect limits on the cross section by a factor of 10; nevertheless, the conservative edge of the 95% systematic band still rules out light dark matter with the benchmark thermal cross section (see Fig. 2 in Geringer-Sameth &amp; Koushiappas (2011)).

Building on these previous works, we have devised a novel statistical approach that operates on stellar-kinematic and gamma-ray data in order to maximize sensitivity to annihilation/decay signals. We present results in a pair of papers. Here (Paper I), we use published stellar-kinematic data to estimate the dark matter density profiles of 20 Milky Way dwarfs. These systems lie at heliocentric distances 23 glyph[lessorsimilar] D / kpc glyph[lessorsimilar] 250 and span a range in luminosity 3 glyph[lessorsimilar] log 10 [ L L / glyph[circledot] ] glyph[lessorsimilar] 7 (McConnachie 2012, and references therein). The primary goal of this work is to quantify uncertainties in the dark matter density profiles-and hence in the J values-that reflect statistical errors due to finite data as well as uncertainties regarding the shapes of dark matter density profiles. We adopt the relatively agnostic modeling approach of Charbonnier et al. (2011) and extend it to the least luminous of the Milky Way's dwarf satellites, where statistical uncertainties can become dominant over systematics; for a thorough examination of the behaviors of random and systematic errors that are inherent to this analysis, we refer the reader to the recent study by Bonnivard et al. (2015). In Paper II, we present a new analysis of gamma-ray data from the Fermi Large Area Telescope (LAT) that, in combination with the estimates of J values presented here, places new limits on dark matter particle interactions.

This paper is structured as follows. In Sec. 2 we introduce the quantities we seek to constrain by briefly reviewing the physics of dark matter annihilation and decay. Sec. 3 describes the parameterization of the density profile and the equations relating this profile to astronomical observables. Sec. 4 introduces the observational datasets used in this analysis and Sec. 5 describes the fitting procedure. We discuss some physical considerations and the importance of truncating the dark matter potentials in Sec. 6. The results of the analysis are presented in Sec. 7. We discuss their relevance in the context of previous work in Sec. 8 and conclude in Sec. 9. Appendix A presents constraints on integrated emission profiles and containment fractions for each dwarf, which are also provided in machine-readable form.

## 2. DARK MATTER SIGNAL

The distribution of dark matter within a system determines the flux of photons and other products of dark matter annihilation and decay. The annihilation rate per volume per time is given by

$$R _ { \text{ann} } = \frac { 1 } { 2 } n ^ { 2 } \langle \sigma v \rangle, \quad \quad \quad ( 1 ) \quad \begin{array} { c } \sin \\ e \vdots \\ \sin \end{array}$$

where n is the number density of dark matter particles and 〈 σ v 〉 is the velocity-averaged annihilation cross section (e.g. Jungman et al. 1996; Gondolo &amp; Gelmini 1991). (Note that n represents the total number density of dark matter irrespective of particle vs. antiparticle; an additional factor of 1 2 appears if particles and antiparticles each constitute half the total abundance.) As dynamical observations offer a handle on the mass density of dark matter but not its number density, it is convenient to write n = ρ/ M χ , where ρ is the mass density and M χ is the mass of a single dark matter particle.

If dark matter decays, the decay rate per volume is given by

$$R _ { \text{dec} } = n \Gamma,$$

where Γ is the decay rate of an individual dark matter particle.

Each annihilation (or decay) gives rise to photons described by the spectrum dN γ ( E ) / dE : the number of photons per energy produced per annihilation (or decay). Photons produced in annihilation or decay travel along straight lines and so the expected flux of photons is determined by dark matter annihilation (or decay) taking place along the entire line of sight in a certain direction. The number of photons per solid angle per energy, area, and time coming from sky-direction ˆ n , in the case of annihilation, is given by

$$\frac { d F ( \hat { \mathbf n }, E ) } { d \Omega d E } = \frac { \langle \sigma v \rangle } { 8 \pi M _ { \chi } ^ { 2 } } \frac { d N _ { \gamma } ( E ) } { d E } \int _ { \ell = 0 } ^ { \infty } d \ell \left [ \rho ( \ell \hat { \mathbf n } ) \right ] ^ { 2 }, \quad ( 3 )$$

while for decay it is given by

$$\frac { d F ( \hat { \mathbf n }, E ) } { d \Omega d E } = \frac { \Gamma } { 4 \pi M _ { \chi } } \frac { d N _ { \gamma } ( E ) } { d E } \int _ { \ell = 0 } ^ { \infty } d \ell \, \rho ( \ell \hat { \mathbf n } ). \quad \ ( 4 )$$

Here, F has the units of photons per area per time, glyph[lscript] is the line of sight distance from Earth and ρ glyph[lscript] ( ˆ n ) is the dark matter mass density at location glyph[lscript] ˆ n . The integrals over the line of sight should be thought of as including the units of [solid angle] -1 . For the case of dark matter annihilation, we define the J -profile to be

$$\frac { d J ( \hat { \mathbf n } ) } { d \Omega } = \int _ { \ell = 0 } ^ { \infty } d \ell \left [ \rho ( \ell \hat { \mathbf n } ) \right ] ^ { 2 }.$$

The corresponding quantity in the case of dark matter decay is the projected mass density along the line of sight:

$$\frac { d J _ { \text{decay} } ( \hat { \mathbf n } ) } { d \Omega } = \int _ { \ell = 0 } ^ { \infty } d \ell \, \rho ( \ell \hat { \mathbf n } ). \quad \quad ( 6 )$$

The terms before the integrals in Eqs. (3) and (4) describe the microscopic physics of dark matter while the J -profile reflects its distribution on large scales. The goal of indirect detection is to use knowledge of dJ / d Ω (or dJ decay / d Ω ) along with observations of photons ( dF / d Ω dE ) to learn something about dark matter particle physics { M χ , 〈 σ v 〉 , Γ } . In the following we will present results for both annihilation and decay.

We consider spherically symmetric density profiles and so dJ / d Ω is a function only of θ , the angular separation between the line of sight ˆ n and the direction towards the center of the dwarf galaxy. The integration in Eqs. (5) and (6) along the line of sight is carried out numerically for every value of θ . We let the variable x denote the distance, along the line of sight, from the point where the line of sight makes its closest approach to the center of the dwarf. That is, the line of sight corresponding to an angular separation θ has an impact parameter b = D sin( θ ), where D is the distance from Earth to the center of the dwarf. The limits of integration in Eqs. (5) and (6) are x = -∞ and x = + ∞ √ , and dx = d glyph[lscript] . The dark matter density is a function of r = b 2 + x 2 , the distance from the center of the dwarf, so that ρ glyph[lscript] ( ˆ n ) is given by ρ ( √ b 2 + x 2 ).

## 3. RECONSTRUCTING THE DARK MATTER POTENTIAL WITH STELLAR KINEMATICS

## 3.1. Dark matter density

In order to accurately quantify uncertainties in the spatial distribution of dark matter it is necessary to use a suitably flexible functional form for the density profile (Bonnivard et al. 2015). Following Charbonnier et al. (2011), we adopt the functional form introduced by Zhao (1996) to generalize the Hernquist (1990) profile. In this spherically symmetric model, the density of dark matter at halo-centric radius r is

$$\rho ( r ) = \frac { \rho _ { s } } { \left ( r / r _ { s } \right ) ^ { \gamma } \left [ 1 + \left ( r / r _ { s } \right ) ^ { \alpha } \right ] ^ { ( \beta - \gamma ) / \alpha } }. \quad \quad ( 7 )$$

This five-parameter profile, normalized by the scale density ρ s , describes a split power law with inner logarithmic slope d log ρ/ d log r | r glyph[lessmuch] rs = -γ and outer logarithmic slope d log ρ/ d log r | r glyph[greatermuch] rs = -β . The transition happens near the scale radius rs , with α specifying its sharpness. For ( α,β,γ ) = (1 3 1) , , one recovers the two-parameter NFW profile that characterizes cold dark matter (CDM) halos formed in dissipationless numerical simulations (Navarro et al. 1997). However, the profile can also describe halos with even steeper central 'cusps' ( γ &gt; 1), or halos with 'cores' of uniform central density ( γ ∼ 0), as are usually inferred from observations of real galaxies (de Blok 2010, and references therein; Walker et al. 2011; Donato et al. 2009). This flexibility lets us explore a wide range of physically plausible dark matter profiles.

## 3.2. Estimation of dark matter profile parameters

From Eq. (3), the flux of annihilation by-products depends on the density of dark matter particles within the source, and thus on the source's gravitational potential. For collisionless stellar systems like dwarf galaxies, the gravitational potential is related fundamentally to the phase-space density of stars f ( r u , ), defined such that f ( r u , ) d 3 r d 3 u gives the expected number of stars lying within the phase-space volume d 3 r d 3 u centered on ( r u , ). However, dwarf galaxies are sufficiently far away that current instrumentation resolves only the projection of their internal phase-space distributions, effectively providing information in just three dimensions: position as projected onto the plane perpendicular to the line of sight, and velocity along the line of sight (from Doppler redshift). Given these limitations, it is common to infer the gravitational potential Φ by considering its relation to moments of the phase-space distribution: the stellar density profile,

$$\nu ( r ) \equiv \int f ( { \mathbf r }, { \mathbf u } ) \, d ^ { 3 } { \mathbf u },$$

and the stellar velocity dispersion profile,

$$\overline { u ^ { 2 } } ( r ) = \overline { u _ { r } ^ { 2 } } ( r ) + \overline { u _ { \theta } ^ { 2 } } ( r ) + \overline { u _ { \phi } ^ { 2 } } ( r ) \quad \quad ( 9 ) \quad R \epsilon$$

$$= \frac { 1 } { \nu ( r ) } \int u ^ { 2 } f ( { \mathbf r }, { \mathbf u } ) \, d ^ { 3 } { \mathbf u }. \quad \quad ( 1 0 ) \quad \frac { \mathbf a } { \mathbf w }$$

Assuming dynamic equilibrium and spherical symmetry, these quantities are related according to the spherical Jeans equation (Binney &amp; Tremaine 2008),

$$\frac { 1 } { \nu ( r ) } \frac { d } { d r } [ \nu ( r ) \overline { u _ { r } ^ { 2 } } ( r ) ] + 2 \frac { \beta _ { a } ( r ) \overline { u _ { r } ^ { 2 } } ( r ) } { r } = - \frac { d \Phi } { d r } = - \frac { G M ( r ) } { r ^ { 2 } }, \ \ ( 1 1 ) \quad \text{recer}$$

where

$$\beta _ { a } ( r ) \equiv 1 - \frac { \overline { 2 u _ { \theta } ^ { 2 } } ( r ) } { \overline { u _ { r } ^ { 2 } } ( r ) }$$

characterizes the orbital anisotropy and the enclosed mass profile

$$M ( r ) = 4 \pi \int _ { 0 } ^ { r } s ^ { 2 } \rho ( s ) d s$$

includes contributions from the dark matter halo.

Equation 11 has the general solution (van der Marel 1994; Mamon &amp; Łokas 2005)

$$\nu ( r ) \overline { u _ { r } ^ { 2 } } ( r ) = \frac { 1 } { f ( r ) } \int _ { r } ^ { \infty } f ( s ) \, \nu ( s ) \, \frac { G M ( s ) } { s ^ { 2 } } \, d s, \quad ( 1 4 )$$

where

$$f ( r ) = 2 \, f ( r _ { 1 } ) \, \exp \left [ \int _ { r _ { 1 } } ^ { r } \beta _ { a } ( s ) s ^ { - 1 } \, d s \right ] \,. \quad \ \ ( 1 5 )$$

Projecting along the line of sight, the mass profile relates to observable profiles, the projected stellar density Σ ( R ), and line-of-sight velocity dispersion σ ( R ), according to (Binney &amp;Tremaine 2008)

$$\int _ { 0 } ^ { \shortmid } \int _ { 0 } ^ { \longmid } \left ( 1 - \beta _ { a } ( r ) \frac { R ^ { 2 } } { r ^ { 2 } } \right ) \frac { \nu ( r ) \overline { u _ { r } ^ { 2 } } ( r ) \ r } { \sqrt { r ^ { 2 } - R ^ { 2 } } } \, d r. \quad ( 1 6 )$$

We use Eq. (16) to fit models for ρ ( r ) and β a ( r ) to observed velocity dispersion and surface brightness profiles under the following assumptions:

- · dynamic equilibrium and spherical symmetry, both implicit in the use of Eq. (11);
- · the stars are distributed according to a Plummer (1911) profile,

$$\nu ( r ) = \frac { 3 L } { 4 \pi R _ { e } ^ { 3 } } \frac { 1 } { ( 1 + R ^ { 2 } / R _ { e } ^ { 2 } ) ^ { 5 / 2 } }, \quad \quad ( 1 7 )$$

implying surface brightness profiles of the form

$$\Sigma ( R ) = \frac { L } { \pi R _ { e } ^ { 2 } } \, \frac { 1 } { ( 1 + R ^ { 2 } / R _ { e } ^ { 2 } ) ^ { 2 } },$$

where L is the total luminosity and Re is the projected halflight radius;

- · the stars contribute negligibly to the gravitational potential, such that Re is the only meaningful parameter in ν ( r ) and Σ ( R );
- · β a = constant;
- · the distribution of stellar velocities is not significantly influenced by the presence of binary stars.

Real galaxies violate all of these assumptions at some level and it is important to consider that the error distributions that we derive for J values will not include the resulting systematic errors. For the present work, we are concerned primarily with quantifying statistical uncertainties that arise from finite sizes of stellar-kinematic samples. For a thorough study of systematic errors that can arise due to different stellar density profiles, non-spherical symmetry, and more complicated behaviors of the velocity anisotropy, we refer the reader to the recent study by Bonnivard et al. (2015).

## 4. OBSERVATIONS

## 4.1. Classical dwarfs

For the Milky Way's eight most luminous 'classical' dwarf galaxies, we adopt projected halflight radii listed in Table 1 of Walker, Mateo &amp; Olszewski (2009, original source is Irwin &amp; Hatzidimitriou 1995). We use the stellar-kinematic data published by Mateo et al. (2008) for Leo I, and by Walker, Mateo &amp; Olszewski (2009) for Carina, Fornax, Sculptor and Sextans.

For Draco, Leo II, and Ursa Minor we use stellar-kinematic data acquired with the Hectochelle spectrograph at the 6.5-m MMT. These data have previously been analyzed by Walker et al. (2009b) and Charbonnier et al. (2011), and will soon be made public (Walker, Olszewski &amp; Mateo, in preparation).

## 4.2. Ultra-faint satellites

For the Milky Way's less luminous 'ultra-faint' satellites discovered over the past seven years (e.g., Willman et al. 2005; Zucker et al. 2006b; Belokurov et al. 2007, 2008; Belokurov et al. 2009), we make use of published data from a variety of sources. We adopt projected halflight radii from the review of McConnachie (2012). Original sources are Martin et al. (2008) and/or discovery papers for satellites discovered after that work (Belokurov et al. 2008; Belokurov et al. 2009).

For Coma Berenices, Canes Venatici I, Canes Venatici II, Leo IV, Leo T, Ursa Major I and Ursa Major II, we utilize the Keck/Deimos velocity data of Simon &amp; Geha (2007), generously provided by Marla Geha (private communication).

For Hercules we adopt the stellar-kinematic data published by Adén et al. (2009b). After using Stromgren photometry to identify and remove foreground contamination independently of velocity, Adén et al. (2009a) use these data to measure a global velocity dispersion of 3 7 . ± 0 9 km s . -1 , slightly smaller than the earlier estimate of 5 1 . ± 0 9 km s . -1 by Simon &amp; Geha (2007). In the interest of deriving conservative estimates on the dark matter density (and ultimately the exclusion limits for particle models), we adopt the Hercules data of Adén et al. (2009b).

For Boötes I, we adopt the stellar-kinematic data published by Koposov et al. (2011). The observing strategy of Koposov et al. (2011) provided ∼ 10 independent velocity measurements of each star over the course of one month, thereby enabling a direct examination of intrinsic velocity variability (e.g., due to unresolved binary stars) that can potentially inflate observed velocity dispersions (and hence the inferred dark matter density) above the values attributable to the galaxy's gravitational potential (McConnachie &amp; Côté 2010). Koposov et al. (2011) directly resolve velocity variability-including near-constant accelerations of ∼ 10 km s -1 month -1 -for ∼ 10% of the stars in their sample. In order to obtain conservative estimates of the dark matter density, we adopt the most restrictive sample of Boötes I members identified by Koposov et al. (2011, marked with a 'B' in their Table 1), which includes only the 37 stars that show no evidence for velocity variability and have small ( ≤ 2 5 km s . -1 ) velocity measurement errors.

For Leo V, we adopt the stellar-kinematic data published by Walker et al. (2009a). These include seven stars identified as likely members. However, while five of these stars lie within three times the projected half-light radius of Leo V ( rh ∼ 40 pc; Belokurov et al. 2008), the other two lie glyph[greaterorsimilar] 10 rh ∼ 400pc away from Leo V's center and in the direction toward Leo IV, which itself lies only ∼ 20 kpc from Leo V. This configuration is highly improbable for a dynamically relaxed system-even given the small number of stars-fueling speculation that Leo IV and Leo V are interacting gravitationally. In that case, the two outermost stars in the Leo V sample may trace a stellar 'bridge' of low surface brightness. However, deep photometric studies have not yielded unambiguous evidence for such a structure (de Jong et al. 2010; Jin et al. 2012; Sand et al. 2010, 2012). In any case, once again in the interest of conservatism, we consider only the five innermost members in the analysis of Leo V.

For Segue 1, we adopt the stellar-kinematic data published by Simon et al. (2011). These data include repeat measurements for tens of stars with velocities originally measured by Geha et al. (2009), enabling an analysis of velocity variability. Martinez et al. (2011) perform a Bayesian analysis and conclude that the presence of unresolved binary stars is likely to have only mild (a ∼ 10% effect) influence on estimates of Segue 1's intrinsic velocity dispersion. Along with velocity measurements, we adopt the 'Bayesian' membership probabilities listed for individual stars in Table 3 of Simon et al. (2011).

For Segue 2, we adopt the stellar-kinematic data published by Kirby et al. (2013) for 25 member stars. These measurements do not resolve Segue 2's internal velocity dispersion, instead placing a 95% upper limit of σ &lt; 2 6 km s . -1 . Aprevious study by Belokurov et al. (2009) reported a velocity dispersion of σ ∼ 3 5 km s . -1 , based on a sample of ∼ 5 member stars that they cautioned might be contaminated by members of a dynamically hotter stream in the same vicinity. Again in the interest of placing conservative limits on the expected dark matter signal, we adopt the sample of Kirby et al. (2013)-not only because it implies a smaller velocity dispersion, but also because it provides a larger number of member stars. While the small velocity dispersion estimated by Kirby et al. (2013) does not, by itself, require that Segue 2 is embedded within a dark matter halo, Kirby et al. (2013) argue that the variance in metallicity among Segue 2's stars constitutes indirect evidence for a dark matter halo (whose deeper potential would help to retain chemically-enriched gas despite strong winds generated by star formation).

Finally, we note that while high-quality stellar-kinematic data sets are available for the Milky-Way satellites Sagittarius and Willman 1, these objects show strong evidence for tidal disruption and/or non-equilibrium kinematics (Ibata et al. 1997; Willman et al. 2011). Since any mass inference based on the Jeans equation relies fundamentally on the assumption of dynamic equilibrium, we do not consider these objects here.

Table 1 lists central coordinates, distances (from the Sun), absolute V-band magnitudes, and projected halflight radii for the Milky Way satellites we consider here. The last column gives the number of member stars with velocity measurements available for the kinematic analysis. Figures 1 and 2 display projected velocity dispersion profiles for each dwarf. These binned profiles are included only for the purpose of display; the fitting that we describe below uses the unbinned data directly.

## 5. FITTING PROCEDURE

Given the available kinematic data and adopted estimates of Re (which fixes Σ ( R ) under the assumption of Plummer profiles), we fit models for ρ ( r ) and β a ( r ) (see Sec. 3) following the procedure of Strigari et al. (2007). Specifically, we assume that the velocity data sample a line-of-sight velocity

Figure 1. Line-of-sight stellar velocity dispersion profiles observed for the Milky Way's eight classical dwarf spheroidal satellites, adopted from Walker et al. (2009b). Solid curves indicate, at each projected radius, the median velocity dispersion of models sampled in the Markov-Chain Monte Carlo analysis. Dashed and dotted curves enclose the central 68% and 95% of velocity dispersion values from the sampled models. The model profiles are fit to the unbinned kinematic data, but clearly show good agreement with the binned data plotted here.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0002v2-figure-9e37b37f-6c92-4354-baff-79c47bde8a12.jpg)

distribution that is Gaussian 4 . Thus we adopt the likelihood function in Charbonnier et al. (2011)) over the following ranges:

$$L = \prod _ { i = 1 } ^ { N } \frac { 1 } { ( 2 \pi ) ^ { 1 / 2 } \left [ \delta _ { u, i } ^ { 2 } + \sigma ^ { 2 } ( R _ { i } ) \right ] ^ { 1 / 2 } } \exp \left [ - \frac { 1 } { 2 } \frac { ( u _ { i } - \langle u \rangle ) ^ { 2 } } { \delta _ { u, i } ^ { 2 } + \sigma ^ { 2 } ( R _ { i } ) } \right ], \ ( 1 9 ) \quad \quad \quad$$

where ui and Ri are the line-of-sight velocity and magnitude of the projected position vector (with respect to the center of the dwarf) of the i th star in the kinematic data set, δ u i , is the observational error in the velocity, and σ ( R ) is the velocity dispersion at projected position R , as specified by model parameters and calculated from Eq. (16). We consider only stars for which published probabilities of membership are greater than 0.95. The bulk velocity of the system 〈 u 〉 is a nuisance parameter that we marginalize over with a flat prior. Besides 〈 u 〉 , the model has six free parameters and we adopt uniform priors (as

4 Given that we allow models with anisotropic and inherently nonGuassian velocity dispersions, this assumption of Gaussianity introduces an internal inconsistency. However, by enabling the simple likelihood function given by Eq. (19), it avoids problems (e.g., arbitrariness of bin boundaries, unresolved dispersions) associated with analyses of binned profiles. A more rigorous treatment would generate the likelihood function directly from a 6-D phase-space distribution function (M. Wilkinson, in preparation).

- · -1 ≤ -log 10 [1 -β a ] ≤ + 1;
- · -4 ≤ log 10 [ ρ / s ( M glyph[circledot] pc -3 )] ≤ + 4;
- · 0 ≤ log 10 [ rs / pc] ≤ + 5;
- · 0 5 . ≤ α ≤ 3;
- · 3 ≤ β ≤ 10;
- · 0 ≤ γ ≤ 1 2. .

In order to sample the parameter space efficiently, we use the nested-sampling Monte Carlo algorithm introduced by Skilling (2004) and implemented in the software package MultiNest (Feroz &amp; Hobson 2008; Feroz et al. 2009), which outputs samples from the model's posterior probability distribution function (PDF).

## 6. PHYSICAL CONSIDERATIONS AND TRUNCATION OF HALO PROFILES

Figure 3 displays samples from the posterior PDFs returned by MultiNest for Fornax and Segue 1-the most luminous

Figure 2. Same as Figure 1, but for the Milky Way's ultra-faint satellites. In many bins the estimated velocity dispersion is zero because the actual dispersion is unresolved by the available data. As in Fig. 1 the points with error bars are for illustration; binned velocity dispersion estimates are not used in the fitting procedure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0002v2-figure-6cd90fab-26f4-48fd-8b90-ed8600b2dfc7.jpg)

classical dwarf and one of the least luminous ultra-faints, respectively. As the model that we adopt for the halo density profile is free (and unconstrained by, e.g., N-body considerations) the kinematic data of each dwarf is compatible with a wide range of profiles. Therefore, we apply three additional filters to the kinematically-allowed dark matter density profiles. The first two involve identifying an outer boundary for a given halo, while the third is a requirement that the halo formed in a cosmologically plausible way.

## 6.1. Halo truncation

tant consequences for the expected dark matter signal from a dwarf galaxy. The data for nearly all dwarf galaxies are consistent with density profiles described by single power laws with logarithmic slopes d log ρ/ d log r &gt; -3 - indeed, despite its unphysically infinite mass, the 'isothermal sphere', characterized by ρ ( r ) ∝ r -2 , has long been used to model kinematics of spheroidal galaxies (Binney &amp; Tremaine 2008). It is therefore important to define some means for preventing the outer parts of a halo - i.e., regions outside the orbits of the observed stellar populations - from dominating the integral used to calculate the J -profile (Eq. (5)).

Given the form of Eq. (7), the annihilation rate will drop rapidly at galactocentric distances r glyph[greatermuch] rs where the density profile is steeply falling ( β &gt; 3). However, within rs the radial distribution of the emission is determined by the slope of the inner density profile γ . For sufficiently cuspy profiles ( γ &gt; 1) the emission is dominated by annihilation near the halo center. For γ ∼ 1 the annihilation rate receives approximately equal contributions from all radii and for γ &lt; 1 the emission comes primarily from the largest radii within rs .

Unfortunately, the current data sets - even for the classical dwarfs - do not place strong upper bounds on rs , thereby allowing emission that extends to an arbitrarily large radius 5 . Therefore, the question of where the halo ends has impor-

5 For spherically symmetric halos the mass exterior to a star's orbit exerts zero net force on that star (Newton 1687); thus stellar kinematics in general carry no information about the mass distribution beyond the orbits of the stars.

## 6.2. Truncating at the outermost observed star

An obvious choice for a conservative truncation radius is that of the outermost member star used to estimate the velocity dispersion profile. For stars well beyond the luminous scale radius Re , the projected distances that we observe are likely to be similar to the de-projected distances r . However, if we observe enough stars close to the center it becomes likely that some of these stars lie at large galactocentric distances. Therefore, we use the entire distribution of projected radii of the kinematic sample to estimate the maximum galactocentric distance r max among those stars.

We estimate r max in the following way. Given spherical symmetry, it is straightforward to find the probability distribution for the unprojected distance to the outermost observed star given the projected distances to the observed stars. We start by considering an individual star. Given its projected

Table 1 Properties of Milky Way Satellites * and Stellar-Kinematic Samples

| object            | RA (J2000) [hh:mm:ss]   | Dec. (J2000) [dd:mm:ss]   | Distance [kpc]   | M V [mag]        | R half [pc]   |   N sample | r max [pc]         |
|-------------------|-------------------------|---------------------------|------------------|------------------|---------------|------------|--------------------|
| Carina            | 06:41:36.7              | - 50:57:58                | 105 ± 6          | - 9 . 1 ± 0 . 5  | 250 ± 39      |        774 | 2224 + 885 - 441   |
| Draco             | 17:20:12.4              | +57:54:55                 | 76 ± 6           | - 8 . 8 ± 0 . 3  | 221 ± 19      |        292 | 1866 + 715 - 317   |
| Fornax            | 02:39:59.3              | - 34:26:57                | 147 ± 12         | - 13 . 4 ± 0 . 3 | 710 ± 77      |       2483 | 6272 + 2616 - 1366 |
| Leo I             | 10:08:28.1              | +12:18:23                 | 254 ± 15         | - 12 . 0 ± 0 . 3 | 251 ± 27      |        267 | 1948 + 794 - 407   |
| Leo II            | 11:13:28.8              | +22:09:06                 | 233 ± 14         | - 9 . 8 ± 0 . 3  | 176 ± 42      |        126 | 824 + 345 - 178    |
| Sculptor          | 01:00:09.4              | - 33:42:33                | 86 ± 6           | - 11 . 1 ± 0 . 5 | 283 ± 45      |       1365 | 2673 + 1099 - 569  |
| Sextans           | 10:13:03.0              | - 01:36:53                | 86 ± 4           | - 9 . 3 ± 0 . 5  | 695 ± 44      |        441 | 2544 + 1109 - 587  |
| Ursa Minor        | 15:09:08.5              | +67:13:21                 | 76 ± 3           | - 8 . 8 ± 0 . 5  | 181 ± 27      |        313 | 1580 + 626 - 312   |
| Bootes I          | 14:00:06.0              | +14:30:00                 | 66 ± 2           | - 6 . 3 ± 0 . 2  | 242 ± 21      |         37 | 544 + 252 - 135    |
| Canes Venatici I  | 13:28:03.5              | +33:33:21                 | 218 ± 10         | - 8 . 6 ± 0 . 2  | 564 ± 36      |        214 | 2030 + 884 - 468   |
| Canes Venatici II | 12:57:10.0              | +34:19:15                 | 160 ± 4          | - 4 . 9 ± 0 . 5  | 74 ± 14       |         25 | 352 + 105 - 28     |
| Coma Berenices    | 12:26:59.0              | +23:54:15                 | 44 ± 4           | - 4 . 1 ± 0 . 5  | 77 ± 10       |         59 | 238 + 103 - 53     |
| Hercules          | 16:31:02.0              | +12:47:30                 | 132 ± 12         | - 6 . 6 ± 0 . 4  | 330 + 75 - 52 |         30 | 638 + 295 - 147    |
| Leo IV            | 11:32:57.0              | - 00:32:00                | 154 ± 6          | - 5 . 8 ± 0 . 4  | 206 ± 37      |         18 | 443 + 197 - 95     |
| Leo V             | 11:31:09.6              | +02:13:12                 | 178 ± 10         | - 5 . 2 ± 0 . 4  | 135 ± 32      |          5 | 201 + 95 - 43      |
| Leo T             | 09:34:53.4              | +17:03:05                 | 417 ± 19         | - 8 . 0 ± 0 . 5  | 120 ± 9       |         19 | 534 + 183 - 60     |
| Segue 1           | 10:07:04.0              | +16:04:55                 | 23 ± 2           | - 1 . 5 ± 0 . 8  | 29 + 8 - 5    |         70 | 139 + 56 - 28      |
| Segue 2           | 02:19:16.0              | +20:10:31                 | 35 ± 2           | - 2 . 5 ± 0 . 3  | 35 ± 3        |         25 | 119 + 45 - 18      |
| Ursa Major I      | 10:34:52.8              | +51:55:12                 | 97 ± 4           | - 5 . 5 ± 0 . 3  | 319 ± 50      |         39 | 732 + 338 - 181    |
| Ursa Major II     | 08:51:30.0              | +63:07:48                 | 32 ± 4           | - 4 . 2 ± 0 . 6  | 149 ± 21      |         20 | 294 + 139 - 74     |

* Central coordinates, distances, absolute magnitudes and projected half-light radii are adopted from the review of McConnachie (2012, see references to original sources therein).

distance R we take the probability of its line of sight distance z (relative to the halo center) to be proportional to the (deprojected) Plummer density profile:

$$P ( z | R ) \infty \left ( 1 + \frac { z ^ { 2 } + R ^ { 2 } } { R _ { e } ^ { 2 } } \right ) ^ { - 5 / 2 }, \quad \quad ( 2 0 ) \quad \text{ an}$$

where Re is the projected halflight radius (Sec. 3). Once the above probability has been normalized (by integrating over z ) we can construct the probability distribution for the unprojected distance r given the projected distance R :

$$P ( r | R ) = \underset { z } { \int } P ( r | z, R ) \, P ( z | R ) \, d z. \quad \quad ( 2 1 ) \quad \text{cai} \\ \text{sta} _ { | }$$

In the above P( r z | , R ) is simply the Dirac delta function δ ( r -√ z 2 + R 2 ). Note that the above integral is zero unless r &gt; R , in which case the delta function picks out two values of z. The result we will need is the cumulative distribution function (CDF) of r given R :

$$C D F ( r | R ) = \int _ { 0 } ^ { r } P ( r ^ { \prime } | R ) d r ^ { \prime } = \frac { \left ( r ^ { 2 } - R ^ { 2 } \right ) ^ { 1 / 2 } \left ( r ^ { 2 } + \frac { 1 } { 2 } ( 3 R _ { h } ^ { 2 } + R ^ { 2 } ) \right ) } { \left ( r ^ { 2 } + R _ { h } ^ { 2 } \right ) ^ { 3 / 2 } }, \quad \text{sensi} \\ \text{for } r > R \text{ and CDF} ( r | R ) = 0 \text{ for } r < R \text{ } \text{ } \text{$$

for r &gt; R and CDF( r R | ) = 0 for r &lt; R .

To find the CDF for the distance to the outermost of n observed stars we simply multiply the CDFs for each of the n stars:

$$\mu _ { 0 } \text{ CDF} _ { \max } ( r | R _ { 1 }, \dots, R _ { n } ) = \text{ CDF} ( r | R _ { 1 } ) \cdots \text{ CDF} ( r | R _ { n } ), \quad ( 2 3 )$$

where each term on the right-hand side is given by Eq. (22) and Ri is the measured projected distance to the i th member star.

For each dwarf, Eq. (23) can be used to estimate the distance to the outermost star used in the Jeans analysis. The median estimate and ± 1 σ confidence intervals for the distance to the outermost member star for each dwarf are shown in the last column of Table 1. When computing J -profiles we truncate all halo profiles obtained from the Jeans/MultiNest sampling analysis at the median estimate to the outermost member star.

Note that this truncation is not imposed on the mass profile in Eq. (14) when calculating the integral in Eq. (16) during the Jeans/MultiNest sampling. However, the integral in Eq. (16) is dominated by the contribution from radii r &lt; Re , such that as long as the truncation radius is larger than the luminous effective radius (as it is for every dwarf galaxy we consider), the result from the Jeans/MultiNest analysis is insensitive to whether or not we truncate the halo density profile at the outermost star in the kinematic sample. We have verified this argument by a controlled experiment and found that the significant effect of truncation is on the subsequent integration over the density profile that enters the calculation of the annihilation signal (see Eq. (5)). For the purpose of this work (quantifying the expected dark matter flux from a

Figure 3. Samples from the posterior probability distribution functions of dark matter halo density profile and velocity anisotropy parameters for Segue 1 (black) and Fornax (red).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0002v2-figure-3b8db11f-83f3-4cf1-9201-cc5f7a0c4729.jpg)

dwarf) this particular choice of truncation is a conservative one. A spherical Jeans analysis cannot, in principle, constrain the mass distribution far beyond the outermost member stars and we therefore set the density to zero at these distances.

## 6.3. Tidal radius

nal and external potentials (i.e., those of the dwarf and Milky Way, respectively), the orbit of the dwarf, and on the orbital configuration within the dwarf (orbits that are prograde with respect to the dwarf's orbit are more easily stripped than those that are retrograde; Read et al. 2006).

For some allowed models, however, the halo density at the galactocentric radius of the outermost star is far smaller than that expected for the Milky Way halo at the same location. This situation would be inconsistent, as the outermost star (and dark matter particles) would likely be lost due to tides. Therefore we impose an additional, physically-motivated filter by requiring that the tidal radius of any acceptable halo be larger than the distance to the outermost star.

The magnitude of the tidal radius rt depends on the inter-

For each kinematically allowed halo profile we follow von Hoerner (1957) and King (1962) and estimate a tidal radius rt by solving

$$\text{ar} \\ \text{is.} \quad \ r _ { t } ^ { 3 } - D ^ { 3 } \frac { M ( r _ { t } ) } { M _ { \text{MW} } ( D ) } \left [ 2 + \frac { \omega ^ { 2 } D ^ { 3 } } { G M _ { \text{MW} } ( D ) } - \frac { d \ln M _ { \text{MW} } } { d \ln r } \Big | _ { r = D } \right ] ^ { - 1 } = 0. \\ \text{before} \ \ D \text{is the distance between the Millvarpoonright with your center and the}$$

Here, D is the distance between the Milky Way center and the dwarf, M rt ( ) is the mass within the tidal radius of the dwarf

galaxy, M MW( ) is the mass enclosed within radius r r in the Milky Way, and ω is the angular speed of the dwarf about the Galactic center (which we take to be a linear speed of 200 km/s divided by the distance from the dwarf to the Galactic center). The Milky Way mass model is taken to be an NFW profile with virial mass M MW = 10 12 M glyph[circledot] , a scale radius rs , MW = 21 5kpc, and a concentration . c MW = 12 (Klypin et al. 2002) (for a more recent treatment of the subject see, e.g., Nesti &amp; Salucci 2013; the results are not sensitive to the uncertainties in the Milky Way mass model).

The expression in Eq. (24) should be taken as a crude approximation for several reasons (see, for example, discussions in Binney &amp; Tremaine (2008); Mo et al. (2010)). First, it makes the assumption that the dwarf galaxy is on a circular orbit (tidal radii for systems with eccentric orbits are not well defined). Second, the three dimensional tidal surface is not of constant radius and thus does not correspond to one unique value for rt . Third, Eq. (24) does not include the effect of orbital dynamics of the particles within the dwarf galaxy itself (manifested as some variance about the angular velocity ω ). Nevertheless, the utility of this prescription has been thoroughly explored in studies of globular clusters, dark matter substructure, and dwarf galaxies (Johnston 1998; Johnston et al. 2001; Taylor &amp; Babul 2001, 2004; Zentner &amp; Bullock 2002, 2003), and gives a reasonably intuitive definition of 'tidal radius'.

Therefore, we reject any halo profiles for which this estimate of the tidal radius is smaller than the radius we estimate for the outermost member in the kinematic samples 6 . This consistency condition turns out to affect only two dwarfs, the ultra-faints Segue 2 and Leo IV, removing profiles with small values of both ρ s and rs (i.e. profiles that do not lie on the typical ρ s vs. rs constraint seen in Fig. 3).

## 6.4. Cosmological considerations

As a final filter, we reject halo profiles that would be cosmologically 'implausible' in the sense that their formation would have required extremely rare peaks in the primordial matter density field. We estimate the rareness of a candidate halo in the following way. Cosmological simulations show that the density profile near the center of a dark matter halo is more or less set at the time of its formation (modulo feedback effects from baryon-physical processes). The outer profile then evolves as mass is accumulated by accretion in a hierarchical scenario (see e.g., Navarro et al. (1997); Wechsler et al. (2006); Gao et al. (2008)). To a first approximation, then, the primordial mass of a halo (the virial mass of the collapsing overdensity) is greater (or on the order of) the mass M within rs today. In addition, the density of the inner parts of a halo is also roughly set at the time of formation, following a distribution in concentration, c ≡ R vir / rs , where R vir is the radius of the virialized region (Navarro et al. 1996).

For example, the scale density of an NFW profile is ρ s = δ ρ c crit f ( c v), where δ c is the characteristic overdensity of a collapsed object and ρ crit (scaling as (1 + z ) 3 in the matterdominated era) is the critical density of the Universe at the time of collapse (e.g. Gao et al. 2008). The function f ( c v) = ln(1 + c v) + c v / (1 + c v) is a function of order unity. There is thus an approximately one-to-one mapping between the collapse redshift and the scale density ρ s of a halo (ignoring the in-

6 We use the minimum possible distance to the outermost star: the largest projected distance to any of the members. This has a conservative impact in the resulting J values.

trinsic variation in concentration Wechsler et al. (2006); Gao et al. (2008)). For simplicity, and due to a lack of knowledge of the initial conditions that give rise to the intrinsic properties of the dwarf galaxies at their formation time, we use the approximation (Navarro et al. 1996)

$$\frac { 1 + z _ { \text{col} } } { 2 0 } & \gtrsim \left ( \frac { \rho _ { s } } { M _ { \odot } p c ^ { - 3 } } \right ) ^ { 1 / 3 }. \\ \text{$r$ limit on the value of the characteristic distance and} \text{$of$}$$

as an upper limit on the value of the characteristic density of each halo in the chains.

The collapse redshift, along with an estimate of a halo's mass at collapse M , can be used to quantify the rarity of the density perturbation that collapsed to form the halo. In order to estimate the mass of the collapsed perturbation we use the conventional theory of spherical collapse. A region of space has collapsed and virialized when its average overdensity (with respect to the background as computed in linear perturbation theory) is δ collapse ∼ 1 686. . The overdensity at redshift z , averaged over a region of mass M , is a Gaussian random variable with mean 0 and standard deviation σ ( M ) D ( z ), where σ ( M ) = 〈 δ 2 M 〉 1 2 / is the standard deviation of the overdensity today when smoothed over a region of mass M and D ( z ), the growth function, quantifies the linear growth of perturbations such that D ( z = 0) = 1. Therefore, the rarity (in standard deviations) of a region of mass M collapsing at redshift z is given by (see e.g. Bond et al. (1991); Bower (1991); Lacey &amp; Cole (1993)),

$$\nu ( M, z ) = \frac { \delta _ { \text{collapse} } } { \sigma ( M ) \mathcal { D } ( z ) }.$$

The probability that such a region has collapsed by redshift z is the tail probability of a standard normal distribution to the right of ν . We can make a rough estimate of the probability that the Milky Way contains any objects of mass M which formed at (or before) a redshift z by incorporating a trials factor which counts the number of 'independent' regions of the Universe of mass M that eventually make up the Milky Way today. We estimate this quantity simply as M MW / M , where M MW = 10 12 M glyph[circledot] is the mass of the Milky Way. We take the mass M of a halo to be the lesser of the mass within rs and the mass within the outermost star. For the halos generated by the Monte Carlo sampling the masses range from about 10 5 to 10 7 M glyph[circledot] . This leads to trials factors between 10 5 and 10 7 . The probability that the Milky Way contains a halo that collapsed with mass M by redshift z is then

$$\mathsf P ( M, z ) = 1 - \Phi \left [ \nu ( M, z ) \right ] ^ { M _ { M W } / M }, \text{ \quad \ \ } ( 2 7 )$$

where Φ ( x ) is the cumulative distribution function of a standard normal distribution and ν ( M z , ) is given by Eq. (26). We apply the constraint by demanding that the Milky Way is 'typical' at the 3 σ level, i.e. we reject a halo if P( M z , ) &lt; . 0 003. For a trials factor of 10 6 this constraint is equivalent to demanding that ν glyph[lessorsimilar] 6. The results are insensitive to the threshold value for P( M z , ) - the probability of the Milky Way containing a halo goes to zero extremely rapidly as the halo's ν value increases beyond about ν ≈ 6.

In applying this cosmological argument a posteriori on the chains we find that in the case of classical dwarfs (for which we have a large number of gravitational stellar tracers) stellar kinematics alone does not allow for rare peaks and this filter essentially has no effect. However for the ultra-faint dwarfs the stellar kinematic data allow dark matter halos that appear

to be extremely rare, suggesting that the origin of such rare peaks is due to the small size of kinematic samples available compared with the classical dwarfs. In practice, the halos which are eliminated are those with ρ s glyph[greaterorsimilar] 1 M glyph[circledot] / pc 3 .

## 7. RESULTS

## 7.1. J-profiles for annihilation and decay

The gamma-ray flux from dark matter annihilation is fully described by the function dJ ( θ / ) d Ω and any dark matter search is properly conducted using this function (e.g. Geringer-Sameth et al. (2014)). However, it is useful to explore the dark matter distributions in dwarfs using summary quantities based on the J -profiles.

Concerning indirect detection, the most important properties of the dwarf galaxies are the amplitude and spatial extent of the J -profile. The amplitude is often given as the integral of dJ ( θ / ) d Ω over some solid angle; i.e. we define

$$J ( \theta ) \equiv \int _ { 0 } ^ { \theta } \frac { d J ( \theta ^ { \prime } ) } { d \Omega } 2 \pi \sin ( \theta ^ { \prime } ) d \theta ^ { \prime }, \quad \ \ ( 2 8 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

with an analogous definition for J decay ( θ ). A detectorindependent amplitude is the J -profile integrated out to the truncation radius of the halo. This corresponds to integrating dJ ( θ / ) d Ω out to an angle θ max = arcsin( r max / D ), where r max is the distance from the center of the dwarf to the outermost member star (Sec. 6.2) and D is the distance from Earth to the dwarf. Integrating the J -profile within θ max gives the total dark matter flux expected from a halo . We use J or J decay to denote this scalar quantity when there is no confusion.

Using the kinematic data discussed in Sec. 4, and applying the formalism described in Secs. 3, 5, and 6, we apply equations (5) and (6) to estimate the kinematically allowed values of J and J decay . Figures 4 and 5 along with Table 2 show the main results of this work.

Figure 4 shows the derived J values for all of the dwarf galaxies considered in this analysis. Error bars represent the 1 σ uncertainty in the value of J (e.g. the 16 th and 84 th percentiles of the posterior distriution). There are many interesting features of this distribution of J values among the dwarf galaxies. Concerning the overall amplitude of an annihilation signal, we find that among the classical dwarfs Draco and Ursa Minor have the largest expected flux. However, Sculptor and, to a lesser extent, Fornax and Carina have the most constraining kinematic sample, giving a very small range of allowed values for J and thus the smallest uncertainties (in agreement with the analysis of Charbonnier et al. 2011). Among the ultra-faint dwarf galaxies, Segue 1 and Ursa Major II have central J values that are higher then those of any of the classical dwarfs. However, the uncertainties in these values are also larger than those of all the classical dwarfs. Of all the dwarf galaxies in this sample, Sculptor has the smallest uncertainty in its total J value (about 15%). On the other hand, all ultra-faint dwarf galaxies exhibit 1 σ errors that are about an order of magnitude (in some cases several orders of magnitude - e.g., Leo IV). This is an outcome of the limited size of the kinematic samples available to date. From the entire ultra-faint sample, Segue 1, Coma Berenices, and Ursa Major II exhibit well-constrained (less than an order of magnitude) J values, and due to their overall high amplitude they should be considered for current and future annihilation searches.

Figure 5 shows the corresponding J decay values for the twenty dwarf galaxies analyzed. As in Fig. 4, error bars represent the 1 σ range in the value of J decay . In the case of dark mat- ter decay scenarios, the emission profile is set by the amount of mass along the line of sight (Eq. (4)) and is less sensitive to the inner slope of the density profile. As a result, the preferential ordering of dwarfs according to their emission amplitude is different than when searching for annihilation products. In the case of decay, the kinematic analysis shows that Draco, and to a lesser extent, Sextans and Sculptor are the sources with the largest expected emission. As in the annihilation case, Fornax and Sculptor have the least systematic uncertainty in their emission amplitudes compared with the rest of the dwarfs.

Many of the ultra-faint dwarfs do not appear to be as promising targets as several of the classical dwarfs. This is due to the conservative way we chose to truncate the halos when computing J and J decay . The outermost observed member stars of the ultra-faints are much closer in than those in the classical dwarfs (see the θ max values along the bottom of Figs. 4 and 5 and the last column of Table 1). Therefore, we simply do not allow the ultra-faint halos to be as extended as those of the classical dwarfs. Since the emission profile from decay is more sensitive to the total mass of the halo, whereas the annihilation profile is more sensitive to the inner slope, the effect of the truncation has different effects when considering annihilation and decay. Whether this truncation is conservative or if it truly reflects the physical sizes of the ultra-faints' dark matter halos is unknown.

It is important to note that the ranking of the various J values implies a consistency check for any detection claim. First, it is likely that if a signal is seen, it will be seen in multiple dwarfs. Consider the pairs Draco/Ursa Minor (classical) and Segue 1/Ursa Major II (ultra-faint). The similar J values among the members of a pair imply that if an annihilation signal is detected in any of these dwarfs there should be a detection of similar amplitude in the other member of the pair (modulo differences in the diffuse gamma-ray background between the dwarfs).

In the farther future, this argument may be turned around to provide an example of 'dark matter particle astronomy': the relative annihilation fluxes measured in multiple dwarfs can be used to constrain the dark matter distributions in their halos. For example, the detection of a signal in a highly constrained dwarf like Sculptor, Draco, or Ursa Minor would immediately tell us something about the dark matter distribution in Segue 1, an object less luminous by 3 orders of magnitude.

## 7.2. Spatial extent

The question of whether a dark matter halo will appear as an extended source for a gamma-ray instrument is important as the additional information available from the angular distribution of emission can be used to increase the signal-to-noise ratio. More point-like sources of emission are more straightforward to detect. However, the detection of a spatially varying annihilation signal immediately reveals properties of the dark matter halo.

Figure 6 shows the constraints on the emission profiles for two benchmark dwarf galaxies: the classical dwarf Draco and the ultra-faint dwarf Segue 1. The envelopes show the ± 1 σ and median values of the emission profiles as a function of the angular separation from the center of the dwarf. It is important to emphasize that none of the curves in the figure correspond to an individual dark matter profile. Rather, the envelope should be thought of as constraining the value of dJ / d Ω or dJ decay / d Ω at a given angular separation. In the lower panels we show the emission profiles integrated over solid an-

Figure 4. Annihilation J -profiles integrated out to θ max for all dwarf galaxies. The error bars show the 1 σ allowed range in the value of J based on the kinematic analysis described in the text. The most prominent dwarf galaxies for annihilation studies are Draco and Ursa Minor (classical) and Coma Berenices, Segue 1, and Ursa Major II (ultra-faint). Numerical values are provided in Table 2 and in machine-readable form as described in Appendix A.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0002v2-figure-b404ef4b-a0bf-4f54-81b5-297a712e69df.jpg)

Figure 5. Same as Fig. 4 for but for dark matter decay ( J decay ). The most prominent dwarf galaxy for decay studies is Draco, with Sculptor and Sextans both a factor of about 3 fainter. Numerical values are provided in Table 2 and in machine-readable form as described in Appendix A.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0002v2-figure-2bebfb06-8f4e-4f9e-b7cf-2913ae49e3d1.jpg)

Table 2 J values for annihilation and decay in dwarf galaxies *

| Dwarf         | θ max [deg]   | θ 0 . 5 [deg]            | θ 0 . 5decay [deg]       | log 10 J ( θ max ) [GeV 2 cm - 5 ]   | log 10 J (0 . 5 ◦ ) [GeV 2 cm - 5 ]   | log 10 J decay ( θ max ) [GeVcm - 2 ]   | log 10 J decay (0 . 5 ◦ ) [GeVcm - 2 ]   |
|---------------|---------------|--------------------------|--------------------------|--------------------------------------|---------------------------------------|-----------------------------------------|------------------------------------------|
| Carina        | 1 . 26        | 0 . 15 + 0 . 15 - 0 . 07 | 0 . 46 + 0 . 16 - 0 . 12 | 17 . 92 + 0 . 19 - 0 . 11            | 17 . 87 + 0 . 10 - 0 . 09             | 18 . 15 + 0 . 34 - 0 . 25               | 17 . 90 + 0 . 17 - 0 . 16                |
| Draco         | 1 . 30        | 0 . 40 + 0 . 16 - 0 . 15 | 0 . 64 + 0 . 06 - 0 . 14 | 19 . 05 + 0 . 22 - 0 . 21            | 18 . 84 + 0 . 12 - 0 . 13             | 18 . 97 + 0 . 17 - 0 . 24               | 18 . 53 + 0 . 10 - 0 . 12                |
| Fornax        | 2 . 61        | 0 . 13 + 0 . 04 - 0 . 05 | 0 . 31 + 0 . 08 - 0 . 05 | 17 . 84 + 0 . 11 - 0 . 06            | 17 . 83 + 0 . 12 - 0 . 06             | 17 . 99 + 0 . 11 - 0 . 08               | 17 . 86 + 0 . 04 - 0 . 05                |
| Leo I         | 0 . 45        | 0 . 13 + 0 . 05 - 0 . 05 | 0 . 22 + 0 . 02 - 0 . 04 | 17 . 84 + 0 . 20 - 0 . 16            | 17 . 84 + 0 . 20 - 0 . 16             | 17 . 91 + 0 . 15 - 0 . 20               | 17 . 91 + 0 . 15 - 0 . 20                |
| Leo II        | 0 . 23        | 0 . 04 + 0 . 05 - 0 . 02 | 0 . 09 + 0 . 03 - 0 . 05 | 17 . 97 + 0 . 20 - 0 . 18            | 17 . 97 + 0 . 20 - 0 . 18             | 17 . 24 + 0 . 35 - 0 . 48               | 17 . 24 + 0 . 35 - 0 . 48                |
| Sculptor      | 1 . 94        | 0 . 15 + 0 . 05 - 0 . 05 | 0 . 48 + 0 . 14 - 0 . 11 | 18 . 57 + 0 . 07 - 0 . 05            | 18 . 54 + 0 . 06 - 0 . 05             | 18 . 47 + 0 . 16 - 0 . 14               | 18 . 19 + 0 . 07 - 0 . 06                |
| Sextans       | 1 . 70        | 0 . 58 + 0 . 32 - 0 . 47 | 0 . 87 + 0 . 10 - 0 . 53 | 17 . 92 + 0 . 35 - 0 . 29            | 17 . 52 + 0 . 28 - 0 . 18             | 18 . 56 + 0 . 25 - 0 . 73               | 17 . 89 + 0 . 13 - 0 . 23                |
| Ursa Minor    | 1 . 37        | 0 . 06 + 0 . 07 - 0 . 03 | 0 . 25 + 0 . 14 - 0 . 09 | 18 . 95 + 0 . 26 - 0 . 18            | 18 . 93 + 0 . 27 - 0 . 19             | 18 . 13 + 0 . 26 - 0 . 18               | 18 . 03 + 0 . 16 - 0 . 13                |
| Boötes I      | 0 . 47        | 0 . 22 + 0 . 05 - 0 . 10 | 0 . 26 + 0 . 02 - 0 . 04 | 18 . 24 + 0 . 40 - 0 . 37            | 18 . 24 + 0 . 40 - 0 . 37             | 17 . 90 + 0 . 23 - 0 . 26               | 17 . 90 + 0 . 23 - 0 . 26                |
| Coma          | 0 . 31        | 0 . 16 + 0 . 02 - 0 . 05 | 0 . 17 + 0 . 01 - 0 . 02 | 19 . 02 + 0 . 37 - 0 . 41            | 19 . 02 + 0 . 37 - 0 . 41             | 17 . 96 + 0 . 20 - 0 . 25               | 17 . 96 + 0 . 20 - 0 . 25                |
| CVnI          | 0 . 53        | 0 . 11 + 0 . 15 - 0 . 09 | 0 . 23 + 0 . 07 - 0 . 17 | 17 . 44 + 0 . 37 - 0 . 28            | 17 . 43 + 0 . 37 - 0 . 28             | 17 . 57 + 0 . 37 - 0 . 73               | 17 . 57 + 0 . 36 - 0 . 72                |
| CVnII         | 0 . 13        | 0 . 07 + 0 . 01 - 0 . 02 | 0 . 07 + 0 . 00 - 0 . 01 | 17 . 65 + 0 . 45 - 0 . 43            | 17 . 65 + 0 . 45 - 0 . 43             | 16 . 97 + 0 . 24 - 0 . 23               | 16 . 97 + 0 . 24 - 0 . 23                |
| Hercules      | 0 . 28        | 0 . 07 + 0 . 08 - 0 . 06 | 0 . 12 + 0 . 03 - 0 . 09 | 16 . 86 + 0 . 74 - 0 . 68            | 16 . 86 + 0 . 74 - 0 . 68             | 16 . 66 + 0 . 42 - 0 . 40               | 16 . 66 + 0 . 42 - 0 . 40                |
| Leo IV        | 0 . 16        | 0 . 05 + 0 . 03 - 0 . 04 | 0 . 08 + 0 . 01 - 0 . 06 | 16 . 32 + 1 . 06 - 1 . 69            | 16 . 32 + 1 . 06 - 1 . 69             | 16 . 12 + 0 . 71 - 1 . 14               | 16 . 12 + 0 . 71 - 1 . 14                |
| Leo V         | 0 . 07        | 0 . 03 + 0 . 01 - 0 . 02 | 0 . 04 + 0 . 00 - 0 . 01 | 16 . 37 + 0 . 94 - 0 . 87            | 16 . 37 + 0 . 94 - 0 . 87             | 15 . 86 + 0 . 46 - 0 . 47               | 15 . 86 + 0 . 46 - 0 . 47                |
| Leo T         | 0 . 08        | 0 . 03 + 0 . 01 - 0 . 02 | 0 . 04 + 0 . 00 - 0 . 01 | 17 . 11 + 0 . 44 - 0 . 39            | 17 . 11 + 0 . 44 - 0 . 39             | 16 . 48 + 0 . 22 - 0 . 25               | 16 . 48 + 0 . 22 - 0 . 25                |
| Segue 1       | 0 . 35        | 0 . 13 + 0 . 05 - 0 . 07 | 0 . 18 + 0 . 01 - 0 . 05 | 19 . 36 + 0 . 32 - 0 . 35            | 19 . 36 + 0 . 32 - 0 . 35             | 17 . 99 + 0 . 20 - 0 . 31               | 17 . 99 + 0 . 20 - 0 . 31                |
| Segue 2       | 0 . 19        | 0 . 07 + 0 . 03 - 0 . 05 | 0 . 10 + 0 . 01 - 0 . 05 | 16 . 21 + 1 . 06 - 0 . 98            | 16 . 21 + 1 . 06 - 0 . 98             | 15 . 89 + 0 . 56 - 0 . 37               | 15 . 89 + 0 . 56 - 0 . 37                |
| Ursa Major I  | 0 . 43        | 0 . 15 + 0 . 08 - 0 . 12 | 0 . 22 + 0 . 03 - 0 . 14 | 17 . 87 + 0 . 56 - 0 . 33            | 17 . 87 + 0 . 56 - 0 . 33             | 17 . 61 + 0 . 20 - 0 . 38               | 17 . 61 + 0 . 20 - 0 . 38                |
| Ursa Major II | 0 . 53        | 0 . 24 + 0 . 06 - 0 . 11 | 0 . 29 + 0 . 02 - 0 . 04 | 19 . 42 + 0 . 44 - 0 . 42            | 19 . 42 + 0 . 44 - 0 . 42             | 18 . 39 + 0 . 25 - 0 . 27               | 18 . 38 + 0 . 25 - 0 . 27                |

* Errors correspond to the 1 σ range of the kinematic analysis (i.e. the 16 th and 84 th percentiles). J ( θ ) is the J -profile (Eq. (5)) integrated over a cone with radius θ (Eq. (28)); θ 0 5 . is the 'half-light radius' for dark matter emission (i.e. the angle containing 50% of the total emission: J ( θ 0 5 . ) = 0 5 . × J ( θ max), see Sec. 7.2). Analogous definitions apply for dark matter decay. Halos are truncated at a radius corresponding to θ max (Sec. 6.2).

Figure 6. Expected emission profiles for annihilation (purple) and decay (green) for Draco and Segue 1. At each angle the solid and dashed lines show the median profiles and the shaded band corresponds to the ± 1 σ distribution as derived in the kinematic analysis. The top panels show log 10 dJ ( θ ) / d Ω (purple), and log 10 J decay ( θ ) / d Ω (green) (see Eqs. (5) and (6)) in units of GeV 2 cm -5 and GeVcm -2 respectively. The lower panels show these quantities integrated over a solid angle of radius θ (Eq. (28)). These envelopes should be thought of as giving the uncertainty in the J -profile and integrated J -profile at each value of θ . (Integrated J vs. θ and J decay vs. θ constraints for all the dwarfs are available in machine-readable form as described in Appendix A.)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0002v2-figure-25fb80da-2cc1-4b2e-8f74-4fe0d8e154b9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0002v2-figure-6d80d326-efbc-4b65-943f-ad9220a9429d.jpg)

gle out to an angular separation θ (Eq. (28)). For Draco we see a familiar result: there is a particular radius at which the differential flux profile is most tightly constrained, and another (slightly larger) angle within which the total annihilation flux is best constrained (Walker et al. 2011; Charbonnier et al. 2011; Bonnivard et al. 2015). The uncertainty in the flux within 0 01 . ◦ is about a factor of 5 and decreases to about 20% when integrating within about 0 3 . ◦ , an angle corresponding to twice the projected half-light radius. For Segue 1, however, the situation is somewhat different. While the integrated J value within 0 01 . ◦ can be inferred to within a factor of 6, similar to the case of Draco, and the minimum uncertainty again occurs when integrating within about twice the halflight radius ( θ ≈ 0 15 . ◦ ), even there J can only be determined to within a factor of 3.5. We do not see the drastic decrease in the uncertainty of Segue 1's expected emission that we see with most of the classical dwarfs. The larger uncertainty for Segue 1 is a direct consequence of the relatively small size of its available kinematic sample.

We can quantify the extent to which halos can be spatially resolved in gamma-ray telescopes by comparing the derived emission profiles for either annihilation or decay with the point spread function (PSF) of specific instruments.

Figures 7 and 8 show the angular distribution of dark matter annihilation and decay. The bands show constraints on the 'containment fraction' curves for the different dwarfs. The containment fraction, at angle θ , is defined simply as J ( θ / ) J ( θ max), where J ( θ ) is given by Eq. (28). Each halo profile gives rise to a containment fraction curve and the dotted line corresponds to the median value of the containment fraction among all the allowed halos, computed at each θ . The shaded band corresponds to the 16 th and 84 th percentiles. For example, the constraint on the 'half-light radius' of the dark matter emission profile is the intersection of the horizontal line y = 0 5 with the shaded band. We use . θ 0 5 . and θ 0 5decay . to denote the half-light radii for J - and J decay -profiles and tabulate them in Table 2.

The curves in Figs. 7 and 8 illustrate the point spread functions (PSFs) of two gamma-ray experiments. The containment fraction of a PSF is simply the probability that a gammaray will be reconstructed within an angle θ of its true origin. The solid blue, magenta, red, and green lines correspond to the PSF of the Fermi-LAT at photon energies of 0.5, 1, 2, and 10 GeV (computed using gtpsf - see software and documentation at the Fermi Science Support Center 7 ). The dashed orange line corresponds to a 2-dimensional Gaussian PSF with a 68% containment angle of 0 1 . ◦ (e.g. a Rayleigh distribution with a mean of 0 083 . ◦ ). This corresponds to the benchmark PSF of current-generation Atmospheric Cerenkov ˇ Telescopes (ACTs). Figure 8 is identical to Fig. 7 but shows the containment fractions for J decay .

We find that for many of the classical dwarfs (Carina, Draco, Fornax, Leo I, Sculptor, Sextans) ACTs should be able to detect extended emission from dark matter annihilation (if the emission can be detected at all) and similarly for some of the ultra-faint dwarfs (Boötes I, Coma Berenices, and Ursa Major II). Regarding Fermi-LAT, at the highest energies ( &gt; 10 GeV) only Draco and, perhaps, Ursa Major II appear to be extended enough to be detected, and therefore any limits derived using Fermi-LAT data will not be affected significantly by the assumption of point sources when it comes to dwarf galaxies (in agreement with Ackermann et al. (2014)).

## 8. COMPARISON WITH OTHER WORK

In order to compare the expected signals derived in this paper with the predictions from other work, Figure 9 shows the distributions of the J -profile integrated within a cone of radius 0 5 . ◦ for all the dwarf galaxies in the sample. In this figure, the green diamonds are the median values of J from the sampled halos in this work with ± 1 σ error bars. The red and blue points show the J values integrated within 0 5 . ◦ reported by Ackermann et al. (2011) and Ackermann et al. (2014, NFW profiles) respectively. The J values in the latter study come from Martinez (2013). The error bars on these points correspond to the 1 σ errors quoted in those studies.

We find that to within an order of magnitude the constraints on J values are consistent with those derived by Ackermann

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0002v2-figure-ca7a1767-a1bf-4a89-aed4-1f43ba5e1b02.jpg)

θ

[deg]

Figure 7. The containment fraction for annihilation as a function of angular distance from the center of each dwarf. The x -axis is in degrees up to θ max for each dwarf. The containment fraction is defined as J ( θ ) / J ( θ max). The dotted line shows the median value of the containment fraction while the shaded band corresponds to its 16 th and 84 th percentiles at each angle. The solid blue, magenta, red, and green lines show the containment fraction of the Fermi-LAT PSF at 0.5, 1, 2 &amp; 10 GeV respectively. The dashed orange line corresponds to the PSF of a typical ACT (68% containment of 0 1 . ◦ ). This figure together with Fig. 4 can be used to estimate the proper normalization of the expected emission within any aperture for each dwarf. The data used to construct this figure and Fig. 8 are available in machine-readable form as described in Appendix A.

et al. (2011) and Ackermann et al. (2014); Martinez (2013). The differences, however, appear to be systematic and not random. In particular, for the ultra-faint dwarfs the central J values from Ackermann et al. (2014); Martinez (2013) are almost always larger than those we find. Of the eight classical dwarfs, which are the least dependent on priors in the analysis presented here, the results are inconsistent by more than 1 σ for Fornax, Leo II, and Sextans when compared to Ackermann et al. (2014). For Fornax and Ursa Minor the results of Ackermann et al. (2011) also disagree at similar significance.

Martinez (2013) introduced a new analysis that effectively models the entire population of Milky Way dwarfs simultaneously. As with the study of Ackermann et al. (2011), all halos are assumed to follow the NFW form, but the (assumed linear) relationship between log v max and log r max is modeled simultaneously with an (assumed linear) relationship between log L and log v max, where L is optical luminosity. Empirical information about kinematics enters only in the form of published estimates of masses enclosed within dwarf half-light radii, which are determined given v max and r max. With respect to the J values estimated by Martinez (2013), the analysis presented here yields tighter constraints for the classical dwarfs and looser constraints for ultra-faints. The former makes intuitive sense as the analysis presented here uses the greater amount of information that is available in the larger, unbinned kinematic data sets that are available for classical dwarfs (e.g. using velocity dispersion as a function of radius). We suspect that the discrepency for the ultra-faints results from the assumption by Martinez (2013) that the scatter about the linear relation between log L and log v max is independent of luminosity and that between log v max and log r max is independent of log v max. Given the hierarchical nature of the model, this assumption enables 'sharing' of information among classical and ultra-faint dwarfs. That is, the model is constrained primarily by the better-sampled classical dwarfs, but because scatter about the modeled relations is assumed to be independent of luminosity, it is impossible for the inferred constraints for ultra-faints to be looser than those inferred at the luminous end (notice the uniformity of error bars among nearly all the blue points in Figure 9, especially for the ultra-faints). Thus, it may be possible that the error bars obtained by Martinez (2013) for the ultra-faint dwarfs are suppressed by the unduly restrictive hierarchical model priors.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0002v2-figure-457542f9-8328-4c14-aab3-868c58ec5eec.jpg)

θ

[deg]

Figure 8. Same as Fig. 7 but for dark matter decay.

## 9. CONCLUSION

Dwarf galaxies represent the cleanest laboratory in which to search for a signal from dark matter annihilation and decay. Understanding the dark matter distribution in these systems is therefore of paramount importance.

Using the latest kinematic observations from 20 dwarf galaxies, we have presented a thorough study of the dark matter content and distribution in these systems. We have shown that, owing to the quality and size of available data sets, the classical dwarf galaxies are better constrained than ultra-faint dwarfs. Relevant to dark matter annihilation searches, Draco and Ursa Minor are the most promising dwarfs, with Sculptor being only a factor 3 fainter but with the smallest uncertainty in its emission. From the ultra-faint dwarf sample, Segue 1, Ursa Major II, and Coma Berenices are by far the brightest dwarf galaxies, albeit with significantly larger systematic uncertainties. Relevant to dark matter decay studies, Draco is the brightest dwarf galaxy in the sample, with Sculptor, Sextans and Ursa Major II a factor of roughly 3 fainter.

In addition, we have explored the angular distribution of annihilation and decay emission from these dwarfs, and show that Atmospheric Cerenkov Telescopes have the necessary anˇ gular resolution to detect extended emission from many of the classical and ultra-faint dwarf galaxies. With the future of gamma-ray astronomy efforts focused on the construction of the Cherenkov Telescope Array (CTA) (Actis et al. 2011), this work provides the necessary ingredients that can be used to derive the expected reach of CTA in the context of dark matter searches.

In summary, the results presented here provide one of the key pieces of information required to test hypotheses about dark matter annihilation and decay using dwarf galaxies. Such searches are currently even more important in light of the recent claims of dark matter annihilation evidence from the Galactic center region. Dwarf galaxies are the most reliable sources where these claims can be tested. This shall be the focus of Paper II.

We acknowledge useful conversations with Gordon Blackadder, Vincent Bonnivard, Celine Combet, David Maurin, Justin Read, Louie Strigari, Mark Wilkinson, and Andrew Zentner. SMK is supported by DOE DE-SC0010010, NSF PHYS-1417505 and NASA NNX13AO94G. MGW acknowledges support from NSF grant AST-1313045. SMK and MGW thank the Aspen Center for Physics for hospitality where part of this work was completed.

Figure 9. J -profiles integrated within 0 5 . ◦ for all the dwarf galaxies. The green diamonds with error bars show the range of J values derived in this work with 1 σ error bars. The red and blue points are the integrated J values for NFW profiles reported in Ackermann et al. (2011) and Ackermann et al. (2014), respectively, with error bars corresponding to the 1 σ uncertainties.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0002v2-figure-126d6ff3-231d-4828-813f-0e3b1e738211.jpg)

## APPENDIX

## A. TABULATED CONSTRAINTS FOR INDIVIDUAL DWARFS

In this appendix we present an additional table with integrated J -profiles and containment fractions for each dwarf as a function of θ . The data in this tables can be used to reconstruct Figs. 4, 5, 7, 8, 9, the lower panels of Fig. 6, and Table 2.

The first few lines of the table are shown in Table 3. The full table, containing all the dwarfs, is available in machine-readable format as an ancillary file. For each value of θ there are four constrained quantities: the integrated J -profiles J ( θ ) and J decay ( θ ), and the containment fractions J ( θ / ) J ( θ max) and J decay ( θ / ) J decay ( θ max). The columns labeled -2 σ, -1 σ,..., + 2 σ correspond to quantiles of the distribution of a quantity among the Monte Carlo realizations of halos. For example, for each halo realization we find the value of J (0 5 . ◦ ). The 2 3 16 50 84, and 97 7 percentiles of this collection of . , , , . J (0 5 . ◦ )'s are listed in columns 3-7 of the row corresponding to θ = 0 5 . ◦ . The table has 50 rows per dwarf, corresponding to log-spaced angles between 0 01 . ◦ and 2 6 . ◦ (the largest θ max of any of the dwarfs).

## B. HALO PROFILE PARAMETER CONSTRAINTS

In order to facilitate the reproducibility of these results we also provide constraints on the individual halo profile parameters for each dwarf. For each of the parameters ρ , s rs , α, β, γ, and β a (see Eqs. 7 and 12), we provide the -2 σ, -1 σ , median, + 1 σ , and + 2 σ quantiles of the posterior distribution (marginalized over the other parameters). The consistency conditions described in Secs. 6.3 and 6.4 have been applied to the posterior samples. Table 4 shows the constraints for the first few dwarfs. The full table is available in machine-readable format as an ancillary file.

Table 3 Integrated J-profiles and Containment Fractions *

| Name         | θ [deg]     | log 10 J ( θ ) [GeV 2 cm - 5 ]   | log 10 J ( θ ) [GeV 2 cm - 5 ]   | log 10 J ( θ ) [GeV 2 cm - 5 ]   |       |       |                 |       |         |         | J ( θ ) / J ( θ max ) -   | J ( θ ) / J ( θ max ) -   | J ( θ ) / J ( θ max ) -   | J ( θ ) / J ( θ max ) -   | J decay ( θ ) / J decay   | J decay ( θ ) / J decay   | J decay ( θ ) / J decay   | J decay ( θ ) / J decay   | J decay ( θ ) / J decay   |
|--------------|-------------|----------------------------------|----------------------------------|----------------------------------|-------|-------|-----------------|-------|---------|---------|---------------------------|---------------------------|---------------------------|---------------------------|---------------------------|---------------------------|---------------------------|---------------------------|---------------------------|
| Name         | -           | 2 σ - 1 σ                        | median                           | + 1 σ + 2 σ                      | - 2 σ | - 1 σ | [GeVcm - median | ] + 1 | σ + 2 σ | - 2 σ   | - 1 σ                     | median                    | + 1 σ                     | + 2 σ                     | - 2 σ                     | - 1 σ                     | - median + 1 σ            | + 2 σ                     |                           |
| Carina 0.010 | 15.38       | 15.88                            | 16.52                            | 16.93 17.33                      | 14.98 | 15.09 | 15.23           | 15.32 | 15.44   | 2.3e-03 | 8.8e-03                   | 3.8e-02                   | 1.0e-01                   | 2.3e-01                   | 3.2e-04                   | 5.4e-04                   | 1.1e-03 2.3e-03           | 4.6e-03                   |                           |
| Carina       | 0.011 15.47 | 15.96                            | 16.57                            | 16.97 17.37                      | 15.08 | 15.19 | 15.32           | 15.41 | 15.52   | 2.8e-03 | 1.0e-02                   | 4.2e-02                   | 1.1e-01                   | 2.5e-01                   | 4.0e-04                   | 6.7e-04                   | 1.3e-03 2.8e-03           | 5.6e-03                   |                           |
| Carina 0.013 | 15.57       | 16.03                            | 16.62                            | 17.01 17.40                      | 15.17 | 15.28 | 15.41           | 15.49 | 15.60   | 3.5e-03 | 1.2e-02                   | 4.7e-02                   | 1.3e-01                   | 2.7e-01                   | 5.0e-04                   | 8.2e-04                   | 1.6e-03 3.4e-03           | 6.7e-03                   |                           |
| Carina 0.014 | 15.66       | 16.11                            | 16.67                            | 17.04 17.44                      | 15.27 | 15.38 | 15.49           | 15.57 | 15.68   | 4.3e-03 | 1.5e-02                   | 5.3e-02                   | 1.4e-01                   | 2.9e-01                   | 6.2e-04                   | 1.0e-03                   | 2.0e-03 4.2e-03           | 8.1e-03                   |                           |
| Carina 0.016 | 15.75       | 16.19                            | 16.72                            | 17.08 17.47                      | 15.37 | 15.47 | 15.58           | 15.66 | 15.76   | 5.3e-03 | 1.7e-02                   | 6.0e-02                   | 1.5e-01                   | 3.1e-01                   | 7.6e-04                   | 1.2e-03                   | 2.5e-03 5.1e-03           | 9.8e-03                   |                           |
| Carina 0.018 |             | 15.84 16.26                      | 16.77                            | 17.11 17.50                      | 15.46 | 15.56 | 15.67           | 15.74 | 15.84   | 6.4e-03 | 2.1e-02                   | 6.7e-02                   | 1.6e-01                   | 3.3e-01                   | 9.5e-04                   | 1.5e-03                   | 3.0e-03 6.2e-03           | 1.2e-02                   |                           |
| Carina       | 0.020       | 15.93 16.34                      | 16.82                            | 17.15 17.53                      | 15.56 | 15.66 | 15.76           | 15.83 | 15.92   | 7.7e-03 | 2.4e-02                   | 7.5e-02                   | 1.8e-01                   | 3.6e-01                   | 1.2e-03                   | 1.9e-03                   | 3.7e-03 7.5e-03           | 1.4e-02                   |                           |

*

This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content.

Table 4 *

## Halo Profile Parameter Constraints

| Name -                                 | log 10 ρ s [M glyph[circledot] pc - 3 ] - 2 σ - 1 σ median + 1 σ + 2 σ - 2 σ                                                                                                  | log 10 r s [pc] - 1 σ median + 1 σ + 2 σ          | α - - 2 σ                                                                  | - 1 σ                    | median + 1 σ + 2 σ - 2 σ   | median + 1 σ + 2 σ - 2 σ                          | β - - 1 σ                | median                   | + 1 σ + 2 σ - 2 σ - 1 σ                                                                             | γ - median + 1 σ + 2 σ - 2 σ   | γ - median + 1 σ + 2 σ - 2 σ                                                       | - log 10 (1 - - - 1 σ median   | β a ) + 1 σ + 2 σ                                 |    |
|----------------------------------------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|---------------------------------------------------|----------------------------------------------------------------------------|--------------------------|----------------------------|---------------------------------------------------|--------------------------|--------------------------|-----------------------------------------------------------------------------------------------------|--------------------------------|------------------------------------------------------------------------------------|--------------------------------|---------------------------------------------------|----|
| Carina Draco Fornax Leo I Leo II . . . | -3.74 -2.86 -1.96 -1.28 -0.73 2.61 -3.36 -2.66 -1.74 -1.09 -0.78 2.85 -2.11 -1.83 -1.49 -1.20 -0.76 2.73 -3.69 -3.08 -2.18 -1.38 -0.92 2.91 -3.24 -2.31 -0.92 -0.03 0.47 1.95 | 2.90 3.31 3.11 3.57 2.93 3.09 3.23 3.80 2.29 2.89 | 3.96 4.63 0.62 4.34 4.83 0.75 3.25 3.51 0.86 4.55 4.91 0.71 4.03 4.77 0.64 | 0.87 1.18 1.39 1.12 1.01 | 1.47 2.01 2.13 1.93 1.76   | 2.89 3.12 2.95 3.19 2.96 3.30 2.94 3.16 2.90 3.16 | 3.75 4.16 4.54 3.99 3.89 | 5.60 6.34 6.97 6.15 5.95 | 8.36 9.71 0.13 0.54 8.69 9.74 0.06 0.29 9.02 9.84 0.03 0.19 8.70 9.76 0.09 0.41 8.56 9.73 0.08 0.35 | 0.95 0.71 0.61 0.84 0.82       | 1.19 -0.41 -0.23 1.16 0.01 0.25 1.17 -0.25 -0.15 1.18 -0.20 -0.03 1.18 -0.88 -0.51 | -0.07 0.54 -0.06 0.14 -0.01    | 0.08 0.25 0.81 0.97 0.02 0.09 0.37 0.72 0.51 0.88 |    |

* This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content.

## REFERENCES

| Abazajian, K. N., Canac, N., Horiuchi, S., &Kaplinghat, M. 2014, Phys. Rev. D, 90, 023526 Abazajian, K. N., &Kaplinghat, M. 2012, Phys. Rev. D, 86,                                                                 |
|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| 083511 -. 2013, Phys. Rev. D, 87, 129902 Ackermann, M., Ajello, M., Albert, A., et al. 2011, Physical Review                                                                                                        |
| Letters, 107, 241302 M., Albert, A., Anderson, B., et al. 2014, Phys. Rev. D, 89,                                                                                                                                   |
| Ackermann, 042001                                                                                                                                                                                                   |
| Actis, M., Agnetta, G., Aharonian, F., et al. 2011, Experimental Astronomy, 32, 193 D., Wilkinson, M. I., Read, J. I., et al. 2009a, ApJ, 706, L150                                                                 |
| Adén, Adén, D., Feltzing, S., Koch, A., et al. 2009b, A&A, 506, 1147                                                                                                                                                |
| Aleksi´, c J., Alvarez, E. A., Antonelli, L. A., et al. 2011, J. Cosmology Astropart. Phys., 6, 35 Aliu, E., Archambault, S., Arlen, T., et al. 2012, Phys. Rev. D, 85, 062001                                      |
| Ashby, M. L. N., Stanford, S. A., Brodwin, M., et al. 2013, ApJS, 209, 22 -. 2014, ApJS, 212, 16 Belokurov, V., Walker, M. G., Evans, N. W., et al. 2009, MNRAS, 397,                                               |
| 1748 -. 2010, ApJ, 712, L103                                                                                                                                                                                        |
| Belokurov et al. 2007, ApJ, 654, 897 -. 2008, ApJ, 686, L83                                                                                                                                                         |
| Binney, J., &Tremaine, S. 2008, Galactic Dynamics: Second Edition (Princeton University Press) Bond, J. R., Cole, S., Efstathiou, G., &Kaiser, N. 1991, ApJ, 379, 440                                               |
| Bonnivard, V., Combet, C., Maurin, D., &Walker, M. G. 2015, MNRAS, 446, 3002                                                                                                                                        |
| Bower, R. G. 1991, MNRAS, 248, 332 Boyarsky, A., Malyshev, D., &Ruchayskiy, O. 2011, Physics Letters 705, 165                                                                                                       |
| B, Charbonnier, A., Combet, C., Daniel, M., et al. 2011, MNRAS, 418, 1526                                                                                                                                           |
| Cholis, I., &Salucci, P. 2012, Phys. Rev. D, 86, 023528 Daylan, T., Finkbeiner, D. P., Hooper, D., et al. 2014, ArXiv e-prints, arXiv:1402.6703                                                                     |
| de Blok, W. J. G. 2010, Advances in Astronomy, 2010, 789293                                                                                                                                                         |
| de Jong, J. T. A., Martin, N. F., Rix, H.-W., et al. 2010, ApJ, 710, 1664 Donato, F., Gentile, G., Salucci, P., et al. 2009, MNRAS, 397, 1169 Essig, R., Sehgal, N., Strigari, L. E., Geha, M., &Simon, J. D. 2010, |
| Phys. Rev. D, 82, 123503                                                                                                                                                                                            |
| Feroz, F., &Hobson, M. P. 2008, MNRAS, 384, 449                                                                                                                                                                     |
| Feroz, F., Hobson, M. P., &Bridges, M. 2009, MNRAS, 398,                                                                                                                                                            |
| 1601 Flaugher, B. 2005, International Journal of Modern Physics A, 20, 3121                                                                                                                                         |
| Gao, L., Navarro, J. F., Cole, S., et al. 2008, MNRAS, 387, 536                                                                                                                                                     |
| Geha, M., Willman, B., Simon, J. D., et al. 2009, ApJ, 692,                                                                                                                                                         |
| 1464 Geringer-Sameth, A., &Koushiappas, S. M. 2011, Physical 107, 241303                                                                                                                                            |
| Review Letters, -. 2012, Phys. Rev. D, 86, 021302                                                                                                                                                                   |
| Geringer-Sameth, A., Koushiappas, S. M., &Walker, M. G. 2014, ArXiv e-prints, arXiv:1410.2242                                                                                                                       |
| Gondolo, P., &Gelmini, G. 1991, Nuclear Physics B, 360, 145                                                                                                                                                         |
| Hernquist, L. 1990, ApJ, 356, 359 Hooper, D., &Goodenough, L. 2011, Physics Letters B, 697, 412 Ibata, R. A., Wyse, R. F. G., Gilmore, G., Irwin, M. J., &Suntzeff,                                                 |
| N. 1997, AJ, 113, 634                                                                                                                                                                                               |
| B. Irwin, M., &Hatzidimitriou, D. 1995, MNRAS, 277, 1354                                                                                                                                                            |
| Jin, S., Martin, N., de Jong, J., et al. 2012, in Astronomical Society of Pacific Conference Series, Vol. 458, Galactic Archaeology:                                                                                |
| Cosmology and the Formation of the Milky Way, ed. W. Aoki,                                                                                                                                                          |
| the Near-Field M. Ishigaki, T. Suda, T. Tsujimoto, &N. Arimoto, 153                                                                                                                                                 |
| K. V. 1998, ApJ, 495, 297                                                                                                                                                                                           |
| Johnston, Johnston, K. V., Sackett, P. D., &Bullock, J. S. 2001, ApJ, 557, 137                                                                                                                                      |
| G., Kamionkowski, M., &Griest, K. 1996, Phys. Rep., 267, 195                                                                                                                                                        |
| Jungman,                                                                                                                                                                                                            |

Kaiser, N., Aussel, H., Burke, B. E., et al. 2002, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 4836, Survey and Other Telescope Technologies and Discoveries, ed. J. A. Tyson &amp; S. Wolff, 154-164

King, I. 1962, AJ, 67, 471

Kirby, E. N., Boylan-Kolchin, M., Cohen, J. G., et al. 2013, ApJ, 770, 16

Klypin, A., Zhao, H., &amp; Somerville, R. S. 2002, ApJ, 573, 597

Koposov, S. E., Gilmore, G., Walker, M. G., et al. 2011, ApJ, 736, 146

Lacey, C., &amp; Cole, S. 1993, MNRAS, 262, 627

Mamon, G. A., &amp; Łokas, E. L. 2005, MNRAS, 363, 705

Martin, N. F., de Jong, J. T. A., &amp; Rix, H.-W. 2008, ApJ, 684, 1075

Martinez, G. D. 2013, ArXiv e-prints, arXiv:1309.2641

Martinez, G. D., Bullock, J. S., Kaplinghat, M., Strigari, L. E., &amp; Trotta, R.

2009, J. Cosmology Astropart. Phys., 6, 14

Martinez, G. D., Minor, Q. E., Bullock, J., et al. 2011, ApJ, 738, 55

Mateo, M., Olszewski, E. W., &amp; Walker, M. G. 2008, ApJ, 675, 201

McConnachie, A. W. 2012, AJ, 144, 4

- McConnachie, A. W., &amp; Côté, P. 2010, ApJ, 722, L209
- Mo, H., van den Bosch, F. C., &amp; White, S. 2010, Galaxy Formation and Evolution (Cambridge: Cambridge University Press)
- Navarro, J. F., Frenk, C. S., &amp; White, S. D. M. 1996, ApJ, 462, 563 -. 1997, ApJ, 490, 493
- Nesti, F., &amp; Salucci, P. 2013, J. Cosmology Astropart. Phys., 7, 16
- Newton, I. 1687, Philosophiae naturalis principia mathematica (London: Royal Society Press)

Plummer, H. C. 1911, MNRAS, 71, 460

- Read, J. I., Wilkinson, M. I., Evans, N. W., Gilmore, G., &amp; Kleyna, J. T. 2006, MNRAS, 366, 429

Salucci, P., Wilkinson, M. I., Walker, M. G., et al. 2012, MNRAS, 420, 2034

Sand, D. J., Seth, A., Olszewski, E. W., et al. 2010, ApJ, 718, 530

Sand, D. J., Strader, J., Willman, B., et al. 2012, ApJ, 756, 79

- Scott, P., Conrad, J., Edsjö, J., et al. 2010, J. Cosmology Astropart. Phys., 1, 31
- Simon, J. D., &amp; Geha, M. 2007, ApJ, 670, 313
- Simon, J. D., Geha, M., Minor, Q. E., et al. 2011, ApJ, 733, 46
- Skilling, J. 2004, in American Institute of Physics Conference Series, Vol. 735, Bayesian Inference and Maximum Entropy Methods in Science and Engineering, ed. R. Fischer, R. Preuss, &amp; U. V. Toussaint, 395-405 Strigari, L. E., Koushiappas, S. M., Bullock, J. S., &amp; Kaplinghat, M. 2007, Phys. Rev. D, 75, 083526
- Taylor, J. E., &amp; Babul, A. 2001, ApJ, 559, 716
- -. 2004, MNRAS, 348, 811

Tyson, J. A., Wittman, D. M., Hennawi, J. F., &amp; Spergel, D. N. 2003,

Nuclear Physics B Proceedings Supplements, 124, 21

van der Marel, R. P. 1994, MNRAS, 270, 271

- von Hoerner, S. 1957, ApJ, 125, 451
- Walker, M. G., Belokurov, V., Evans, N. W., et al. 2009a, ApJ, 694, L144
- Walker, M. G., Combet, C., Hinton, J. A., Maurin, D., &amp; Wilkinson, M. I. 2011, ApJ, 733, L46+
- Walker, M. G., Mateo, M., Olszewski, E. W., et al. 2009b, ApJ, 704, 1274 Walker, Mateo &amp; Olszewski. 2009, AJ, 137, 3100
- Walsh, S. M., Jerjen, H., &amp; Willman, B. 2007, ApJ, 662, L83
- Wechsler, R. H., Zentner, A. R., Bullock, J. S., Kravtsov, A. V., &amp; Allgood, B. 2006, ApJ, 652, 71

Willman, B., Geha, M., Strader, J., et al. 2011, AJ, 142, 128

Willman et al. 2005, ApJ, 626, L85

Zentner, A. R., &amp; Bullock, J. S. 2002, Phys. Rev. D, 66, 043003

- -. 2003, ApJ, 598, 49

Zhao, H. 1996, MNRAS, 278, 488

Zucker et al. 2006a, ApJ, 650, L41

- -. 2006b, ApJ, 643, L103
Alex Geringer-Sameth Savvas M. Koushiappas Matthew Walker
Dwarf galaxy annihilation and decay emission profiles for dark matter experiments
Gamma-ray searches for dark matter annihilation and decay in dwarf galaxies
rely on an understanding of the dark matter density profiles of these systems.
Conversely, uncertainties in these density profiles propagate into the derived
particle physics limits as systematic errors. In this paper we quantify the
expected dark matter signal from 20 Milky Way dwarfs using a uniform analysis
of the most recent stellar-kinematic data available. Assuming that the observed
stellar populations are equilibrium tracers of spherically-symmetric
gravitational potentials that are dominated by dark matter, we find that
current stellar-kinematic data can predict the amplitudes of annihilation
signals to within a factor of a few for the ultra-faint dwarfs of greatest
interest. On the other hand, the expected signal from several classical dwarfs
(with high-quality observations of large numbers of member stars) can be
localized to the ~20% level. These results are important for designing
maximally sensitive searches in current and future experiments using space and
ground-based instruments.