We present a new estimate of the mass of the Milky Way, inferred via a Bayesian approach by making use of tracers of the circular velocity in the disk plane and stars in the stellar halo, as from the publicly available {\tt galkin} compilation. We use the rotation curve method to determine the dark matter distribution and total mass under different assumptions for the dark matter profile, while the total stellar mass is constrained by surface stellar density and microlensing measurements. We also include uncertainties on the baryonic morphology via Bayesian model averaging, thus converting a potential source of systematic error into a more manageable statistical uncertainty. We evaluate the robustness of our result against various possible systematics, including rotation curve data selection, uncertainty on the Sun's velocity $V_0$, dependence on the dark matter profile assumptions, and choice of priors. We find the Milky Way's dark matter virial mass to be $\log_{10}M_{200}^{\rm DM}/ {\rm M_\odot} = 11.92^{+0.06}_{-0.05}{\rm(stat)}\pm{0.28}\pm0.27{\rm(syst)}$ ($M_{200}^{\rm DM}=8.3^{+1.2}_{-0.9}{\rm(stat)}\times10^{11}\,{\rm M_\odot}$). We also apply our framework to Gaia DR2 rotation curve data and find good statistical agreement with the above results.

Dwarf spheroidal galaxies that form in halo substructures provide stringent constraints on dark matter annihilation. Many ultrafaint dwarfs discovered with modern surveys contribute significantly to these constraints. At present, because of the lack of abundant stellar kinematic data for the ultrafaints, non-informative prior assumptions are usually made for the parameters of the density profiles. Based on semi-analytic models of dark matter subhalos and their connection to satellite galaxies, we present more informative and realistic satellite priors. We show that our satellite priors lead to constraints on the annihilation rate that are between a factor of 2 and a factor of 7 weaker than under non-informative priors. As a result, the thermal relic cross section can at best only be excluded (with 95% probability) for dark matter masses of $\lesssim 40$ GeV from dwarf spheroidal data, assuming annihilation into $b\bar{b}$.

We search for a dark matter signal in 11 years of Fermi-LAT gamma-ray data from 27 Milky Way dwarf spheroidal galaxies with spectroscopically measured $J$-factors. Our analysis includes uncertainties in $J$-factors and background normalisations and compares results from a Bayesian and a frequentist perspective. We revisit the dwarf spheroidal galaxy Reticulum II, confirming that the purported gamma-ray excess seen in Pass 7 data is much weaker in Pass 8, independently of the statistical approach adopted. We introduce for the first time posterior predictive distributions to quantify the probability of a dark matter detection from another dwarf galaxy given a tentative excess. A global analysis including all 27 dwarfs shows no indication for a signal in nine annihilation channels. We present stringent new Bayesian and frequentist upper limits on the dark matter cross section as a function of dark matter mass. The best-fit dark matter parameters associated with the Galactic Centre excess are excluded by at least 95% confidence level/posterior probability in the frequentist/Bayesian framework in all cases. However, from a Bayesian model comparison perspective, dark matter annihilation within the dwarfs is not strongly disfavoured compared to a background-only model. These results constitute the highest exposure analysis on the most complete sample of dwarfs to date. Posterior samples and likelihood maps from this study are publicly available.

We develop a novel Bayesian methodology aimed at reliably and precisely inferring the distribution of dark matter within the Milky Way using rotation curve data. We identify a subset of the available rotation curve tracers that are mutually consistent with each other, thus eliminating data sets that might suffer from systematic bias. We investigate three different models for the mass distribution of the luminous (baryonic) component that bracket the range of likely morphologies. We demonstrate the statistical performance of our method on simulated data in terms of coverage, fractional distance, and mean squared error. Applying it to Milky Way data we measure the local dark matter density at the solar circle $\rho_0$ to be $\rho_0 = 0.43\pm 0.02(\rm{stat})\pm0.01(\rm{sys})$ GeV/cm$^3$, with an accuracy $\sim$ 6%. This result is robust to the assumed baryonic morphology. The scale radius and inner slope of the dark matter profile are degenerate and cannot be individually determined with high accuracy. We show that these results are robust to several possible residual systematic errors in the rotation curve data.

We explore two possible scenarios to explain the observed gamma-ray emission associated with the atypical globular cluster Omega-Centauri: emission from millisecond pulsars (MSP) and dark matter (DM) annihilation. In the first case the total number of MSPs needed to produce the gamma-ray flux is compatible with the known (but not confirmed) MSP candidates observed in X-rays. A DM interpretation is motivated by the possibility of Omega-Centauri being the remnant core of an ancient dwarf galaxy hosting a surviving DM component. At least two annihilation channels, light quarks and muons, can plausibly produce the observed gama-ray spectrum. We outline constraints on the parameter space of DM mass versus the product of the pair-annihilation cross section and integrated squared DM density (the so-called J-factor). We translate upper limits on the dark matter content of Omega-Centauri into lower limits on the annihilation cross section. This shows s-wave annihilation into muons to be inconsistent with CMB observations, while a small window for annihilation into light quarks is allowed. Further analysis of Omega-Centauri's internal kinematics, and/or additional information on the resident MSP population will yield much stronger constraints and shed light about the origin of this otherwise mysterious gamma-ray source.

We report the discovery of a Milky Way satellite in the constellation of Antlia. The Antlia 2 dwarf galaxy is located behind the Galactic disc at a latitude of $b\sim 11^{\circ}$ and spans $1.26$ degrees, which corresponds to $\sim2.9$ kpc at its distance of 130 kpc. While similar in spatial extent to the Large Magellanic Cloud, Antlia 2 is orders of magnitude fainter at $M_V=-9$ mag, making it by far the lowest surface brightness system known (at $\sim31.9$ mag/arcsec$^2$), $\sim100$ times more diffuse than the so-called ultra diffuse galaxies. The satellite was identified using a combination of astrometry, photometry and variability data from \textit{Gaia} Data Release 2, and its nature confirmed with deep archival DECam imaging, which revealed a conspicuous BHB signal. We have also obtained follow-up spectroscopy using AAOmega on the AAT, identifying 159 member stars, and we used them to measure the dwarf's systemic velocity, $290.9\pm0.5$km/s, its velocity dispersion, $5.7\pm1.1$ km/s, and mean metallicity, [Fe/H]$=-1.4$. From these properties we conclude that Antlia 2 inhabits one of the least dense Dark Matter (DM) halos probed to date. Dynamical modelling and tidal-disruption simulations suggest that a combination of a cored DM profile and strong tidal stripping may explain the observed properties of this satellite. The origin of this core may be consistent with aggressive feedback, or may even require alternatives to cold dark matter (such as ultra-light bosons).

We develop and implement a model to analyze the internal kinematics of galaxy clusters that may contain subpopulations of galaxies that do not independently trace the cluster potential. The model allows for substructures within the cluster environment, disentangles cluster members from contaminating foreground and background galaxies, and includes an overall cluster rotation term as part of the cluster kinematics. We estimate the cluster velocity dispersion and/or mass while marginalizing over uncertainties in all of the above complexities. In a first application to our published data for Abell 267 (A267), we find no evidence for cluster rotation but we identify up to five distinct galaxy subpopulations. We use these results to explore the sensitivity of inferred cluster properties to the treatment of substructure. Compared to a model that assumes no substructure, our substructure model reduces the dynamical mass of A267 by $\sim 20\%$ and shifts the cluster mean velocity by $\sim 100$ km s$^{-1}$, approximately doubling the offset with respect to the velocity of A267's brightest cluster galaxy. Embedding the spherical Jeans equation within this framework, we infer for A267 a dark matter halo of mass $M_{200}=6.77\pm1.06\times10^{14}M_\odot/h$, concentration $\log_{10}c_{200}=0.61\pm0.39$, consistent with the mass-concentration relation found in cosmological simulations.

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

We report the discovery of a nearby dwarf galaxy in the constellation of Hydrus, between the Large and the Small Magellanic Clouds. Hydrus 1 is a mildy elliptical ultra-faint system with luminosity $M_V\sim$ -4.7 and size $\sim$ 50 pc, located 28 kpc from the Sun and 24 kpc from the LMC. From spectroscopy of $\sim$ 30 member stars, we measure a velocity dispersion of 2.7 km/s and find tentative evidence for a radial velocity gradient consistent with 3 km/s rotation. Hydrus 1's velocity dispersion indicates that the system is dark matter dominated, but its dynamical mass-to-light ratio M/L $\sim$ 66 is significantly smaller than typical for ultra-faint dwarfs at similar luminosity. The kinematics and spatial position of Hydrus~1 make it a very plausible member of the family of satellites brought into the Milky Way by the Magellanic Clouds. While Hydrus 1's proximity and well-measured kinematics make it a promising target for dark matter annihilation searches, we find no evidence for significant gamma-ray emission from Hydrus 1. The new dwarf is a metal-poor galaxy with a mean metallicity [Fe/H]=-2.5 and [Fe/H] spread of 0.4 dex, similar to other systems of similar luminosity. Alpha-abundances of Hyi 1 members indicate that star-formation was extended, lasting between 0.1 and 1 Gyr, with self-enrichment dominated by SN Ia. The dwarf also hosts a highly carbon-enhanced extremely metal-poor star with [Fe/H] $\sim$ -3.2 and [C/Fe] $\sim$ +3.0.

We present constraints on the annihilation cross section of WIMP dark matter based on the joint statistical analysis of four dwarf galaxies with VERITAS. These results are derived from an optimized photon weighting statistical technique that improves on standard imaging atmospheric Cherenkov telescope (IACT) analyses by utilizing the spectral and spatial properties of individual photon events. We report on the results of $\sim$230 hours of observations of five dwarf galaxies and the joint statistical analysis of four of the dwarf galaxies. We find no evidence of gamma-ray emission from any individual dwarf nor in the joint analysis. The derived upper limit on the dark matter annihilation cross section from the joint analysis is $1.35\times 10^{-23} {\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the bottom quark ($b\bar{b}$) final state, $2.85\times 10^{-24}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the tau lepton ($\tau^{+}\tau^{-}$) final state and $1.32\times 10^{-25}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the gauge boson ($\gamma\gamma$) final state.

We present results from MMT/Hectochelle spectroscopy of red giant candidate stars along the line of sight to the recently-discovered Galactic satellite Crater 2. Modelling the joint distribution of stellar positions, velocities and metallicities as a mixture of Crater 2 and Galactic foreground populations, we identify 62 members of Crater 2, for which we resolve line-of-sight velocity dispersion 2.7 +/- 0.3 km/s about mean velocity of 87.5 +/- 0.4 km/s. We also resolve a metallicity dispersion 0.22 about a mean of [Fe/H]=-1.98 +/- 0.1 that is 0.28 +/- 0.14 poorer than is estimated from photometry. Despite Crater 2's relatively large size (projected halflight radius R(h)=1 kpc) and intermediate luminosity (M_V =-8), its velocity dispersion is the coldest that has been resolved for any dwarf galaxy. These properties make Crater 2 the most extreme low-density outlier in dynamical as well as structural scaling relations among the Milky Way's dwarf spheroidals. Even so, under assumptions of dynamical equilibrium and negligible contamination by unresolved binary stars, the observed velocity distribution implies a gravitationally dominant dark matter halo, with dynamical mass 4.4 x 10^6 Msun and mass-to-light ratio M/L=53 enclosed within a radius of 1 kpc, where the equivalent circular velocity is 4.3 km/s.

Gamma-ray experiments seeking to detect evidence of dark matter annihilation in dwarf spheroidal galaxies require knowledge of the distribution of dark matter within these systems. We analyze the effects of flattening on the annihilation (J) and decay (D) factors of dwarf spheroidal galaxies with both analytic and numerical methods. Flattening has two consequences: first, there is a geometric effect as the squeezing (or stretching) of the dark matter distribution enhances (or diminishes) the J-factor; second, the line of sight velocity dispersion of stars must hold up the flattened baryonic component in the flattened dark matter halo. We provide analytic formulae and a simple numerical approach to estimate the correction to the J- and D-factors required over simple spherical modeling. The formulae are validated with a series of equilibrium models of flattened stellar distributions embedded in flattened dark-matter distributions. We compute corrections to the J- and D-factors for the Milky Way dwarf spheroidal galaxies under the assumption that they are all prolate or all oblate and find that the hierarchy of J-factors for the dwarf spheroidals is slightly altered (typical correction factors for an ellipticity of $0.4$ are $0.75$ for the oblate case and $1.6$ for the prolate case). We demonstrate that spherical estimates of the D-factors are very insensitive to the flattening and introduce uncertainties significantly less than the uncertainties in the D-factors from the other observables for all the dwarf spheroidals (for example, ${}^{+10\mathrm{\,percent}}_{-3\mathrm{\,percent}}$ for a typical ellipticity of $0.4$). We conclude by investigating the spread in correction factors produced by triaxial figures and provide uncertainties in the J-factors for the dwarf spheroidals using different physically-motivated assumptions for their intrinsic shape and axis alignments. (abridged)

J-factors (or D-factors) describe the distribution of dark matter in an astrophysical system and determine the strength of the signal provided by annihilating (or decaying) dark matter respectively. We provide simple analytic formulae to calculate the J-factors for spherical cusps obeying the empirical relationship between enclosed mass, velocity dispersion and half-light radius. We extend the calculation to the spherical Navarro-Frenk-White (NFW) model, and demonstrate that our new formulae give accurate results in comparison to more elaborate Jeans models driven by Markov Chain Monte Carlo methods. Of the known ultrafaint dwarf spheroidals, we show that Ursa Major II, Reticulum II, Tucana II and Horologium I have the largest J-factors and so provide the most promising candidates for indirect dark matter detection experiments. Amongst the classical dwarfs, Draco, Sculptor and Ursa Minor have the highest J-factors. We show that the behaviour of the J-factor as a function of integration angle can be inferred for general dark halo models with inner slope $\gamma$ and outer slope $\beta$. The central and asymptotic behaviour of the J-factor curves are derived as a function of the dark halo properties. Finally, we show that models obeying the empirical relation on enclosed mass and velocity dispersion have J-factors that are most robust at the integration angle equal to the projected half-light radius of the dSph divided by heliocentric distance. For most of our results, we give the extension to the D-factor which is appropriate for the decaying dark matter picture.

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

We present a search for gamma-ray emission from the direction of the newly discovered dwarf galaxy Reticulum II. Using Fermi-LAT data, we detect a signal that exceeds expected backgrounds between ~2-10 GeV and is consistent with annihilation of dark matter for particle masses less than a few x 10^2 GeV. Modeling the background as a Poisson process based on Fermi-LAT diffuse models, and taking into account trials factors, we detect emission with p-value less than 9.8 x 10^-5 (>3.7 sigma). An alternative, model-independent treatment of background reduces the significance, raising the p-value to 9.7 x 10^-3 (2.3 sigma). Even in this case, however, Reticulum II has the most significant gamma-ray signal of any known dwarf galaxy. If Reticulum II has a dark matter halo that is similar to those inferred for other nearby dwarfs, the signal is consistent with the s-wave relic abundance cross section for annihilation.

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

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.

We present a new formalism designed to discover dark matter annihilation occurring in the Milky Way's dwarf galaxies. The statistical framework extracts all available information in the data by simultaneously combining observations of all the dwarf galaxies and incorporating the impact of particle physics properties, the distribution of dark matter in the dwarfs, and the detector response. The method performs maximally powerful frequentist searches and produces confidence limits on particle physics parameters. Probability distributions of test statistics under various hypotheses are constructed exactly, without relying on large sample approximations. The derived limits have proper coverage by construction and claims of detection are not biased by imperfect background modeling. We implement this formalism using data from the Fermi Gamma-ray Space Telescope to search for an annihilation signal in the complete sample of Milky Way dwarfs whose dark matter distributions can be reliably determined. We find that the observed data is consistent with background for each of the dwarf galaxies individually as well as in a joint analysis. The strongest constraints are at small dark matter particle masses. Taking the median of the systematic uncertainty in dwarf density profiles, the cross section upper limits are below the pure s-wave weak scale relic abundance value (2.2 x 10^-26 cm^3/s) for dark matter masses below 26 GeV (for annihilation into b quarks), 29 GeV (tau leptons), 35 GeV (up, down, strange, and charm quarks and gluons), 6 GeV (electrons/positrons), and 114 GeV (two-photon final state). For dark matter particle masses less than 1 TeV, these represent the strongest limits obtained to date using dwarf galaxies.

In this white paper, submitted as a part of Snowmass 2013 (subgroup CF2), we examine the current status and future prospects of the VERITAS indirect dark matter detection program. The VERITAS array of imaging atmospheric Cherenkov telescopes (IACTs), sensitive in the 0.1-50 TeV regime, is in the process of completing a multi-year program aimed at detecting signatures of neutralino dark matter. This program is spread out over a range of astrophysical targets which can potentially yield definitive signatures of neutralino self-annihilation such as dwarf spheroidal galaxies (dSphs) and the center of the Milky Way galaxy. While the program is still in progress, initial results on dSphs have produced very competitive upper limits on the thermally averaged cross-section of neutralino self-annihilation as well as strongly constraining leptophillic dark matter models (such as those inferred from the PAMELA and AMS positron fraction results). As the program continues over the next 5 years, VERITAS observations will yield some of the strongest constraints available from IACTs, pushing down the limit on the thermally averaged cross-section of neutralino self-annihilation to the expected natural thermal relic scale.

The VERITAS array of Cherenkov telescopes, designed for the detection of gamma-rays in the 100 GeV-10 TeV energy range, performs dark matter searches over a wide variety of targets. VERITAS continues to carry out focused observations of dwarf spheroidal galaxies in the Local Group, of the Milky Way galactic center, and of Fermi-LAT unidentified sources. This report presents our extensive observations of these targets, new statistical techniques, and current constraints on dark matter particle physics derived from these observations.

We perform a joint analysis of dwarf galaxy data from the Fermi Gamma-ray Space Telescope in search of dark matter annihilation into a gamma-ray line. We employ a novel statistical method that takes into account the spatial and spectral information of individual photon events from a sample of seven dwarf galaxies. Dwarf galaxies show no evidence of a gamma-ray line between 10 GeV and 1 TeV. The subsequent upper limit on the annihilation cross section to a two-photon final state is 3.9(+7.1)(-3.7) x 10^-26 cm^3/s at 130 GeV, where the errors reflect the systematic uncertainty in the distribution of dark matter within the dwarf galaxies.

We present a statistical technique which can be used to detect the presence and properties of moving sources contributing to a diffuse background. The method is a generalization of the 2-point correlation function to include temporal as well as spatial information. We develop a formalism which allows for a derivation of the spacetime 2-point function in terms of the properties of the contributing sources. We test this technique in simulated sky maps, and demonstrate its robustness in identifying the presence of moving and stationary sources. Applications of this formalism to the diffuse gamma-ray background include searches for solar system bodies, fast moving primordial black holes, and dense cores of dark matter proto-halos in the solar neighborhood. Other applications include detecting the contribution of energetic neutrinos originating in the solar system, as well as probing compact objects in long timeline lensing experiments.

Dwarf spheroidal galaxies are known to be excellent targets for the detection of annihilating dark matter. We present new limits on the annihilation cross section of Weakly Interacting Massive Particles (WIMPs) based on the joint analysis of seven Milky Way dwarfs using a frequentist Neyman construction and Pass 7 data from the Fermi Gamma-ray Space Telescope. We exclude generic WIMP candidates annihilating into b-bbar with mass less than 40 GeV that reproduce the observed relic abundance. To within 95% systematic errors on the dark matter distribution within the dwarfs, the mass lower limit can be as low as 19 GeV or as high as 240 GeV. For annihilation into tau+tau- these limits become 19 GeV, 13 GeV, and 80 GeV respectively.

We present a statistical framework which can be used to determine the contribution of an unresolved population of pulsars to the gamma-ray background. This formalism is based on the joint analysis of photon time series over extended regions of the sky. We demonstrate the robustness of this technique in controlled simulations of pulsar populations, and show that the Fermi Gamma-ray Space Telescope can be used to detect a pulsar contribution as small as 0.1% of the gamma-ray background. This technique is sensitive to pulsar populations with photon fluxes greater than ~1E-10 cm^-2 s^-1. The framework is extensible to arbitrarily complex searches for periodicity and can therefore be tailored to specific applications such as all-sky surveys and studies of the Galactic center and globular clusters.