The abundance and narrow magnitude dispersion of Red Clump (RC) stars make them a popular candidate for mapping the morphology of the bulge region of the Milky Way. Using an estimate of the RC's intrinsic luminosity function, we extracted the three-dimensional density distribution of the RC from deep photometric catalogues of the VISTA Variables in the Via Lactea (VVV) survey. We used maximum entropy based deconvolution to extract the spatial distribution of the bulge from Ks-band star counts. We obtained our extrapolated non-parametric model of the bulge over the inner 40 by 40 degrees squared region of the Galactic centre. Our reconstruction also naturally matches onto a parametric fit to the bulge outside the VVV region and inpaints overcrowded and high extinction regions. We found a range of bulge properties consistent with other recent investigations based on the VVV data. In particular, we estimated the bulge mass to be in the range 13 to 17 billion solar masses, the X-component to be between 18% and 25% of the bulge mass, and the bulge angle with respect to the Sun-Galactic centre line to be between 18 and 32 degrees. Studies of the Fermi Large Area Telescope (LAT) gamma-ray Galactic centre excess suggests that the excess may be traced by Galactic bulge distributed sources. We applied our deconvolved density in a template fitting analysis of this Fermi-LAT GeV excess and found an improvement in the fit compared to previous parametric based templates.

Cosmic rays in galaxy clusters are unique probes of energetic processes operating with large-scale structures in the Universe. Precise measurements of cosmic rays in galaxy clusters are essential for improving our understanding of non-thermal components in the intracluster-medium (ICM) as well as the accuracy of cluster mass estimates in cosmological analyses. In this paper, we perform a cross-correlation analysis with the extragalactic gamma-ray background and the thermal Sunyaev-Zel'dovich (tSZ) effect in the cosmic microwave background. The expected cross-correlation signal would contain rich information about the cosmic-ray-induced gamma-ray emission in the most massive galaxy clusters at $z\sim0.1-0.2$. We analyze the gamma-ray background map with 8 years of data taken by the Large Area Telescope onboard Fermi satellite and the publicly available tSZ map by Planck. We confirm that the measured cross-correlation is consistent with a null detection, and thus it enables us to put the tightest constraint on the acceleration efficiency of cosmic ray protons at shocks in and around galaxy clusters. We find the acceleration efficiency must be below 5\% with a $2\sigma$ confidence level when the hydrostatic mass bias of clusters is assumed to be 30\%, and our result is not significantly affected by the assumed value of the hydrostatic mass bias. Our constraint implies that the non-thermal cosmic-ray pressure in the ICM can introduce only a $\le 3\%$ level of the hydrostatic mass bias, highlighting that cosmic rays alone do not account for the mass bias inferred by the Planck analyses. Finally, we discuss future detectability prospects of cosmic-ray-induced gamma rays from the Perseus cluster for the Cherenkov Telescope Array.

Axion dark matter (DM) may convert to radio-frequency electromagnetic radiation in the strong magnetic fields around neutron stars. The radio signature of such a process would be an ultra-narrow spectral peak at a frequency determined by the mass of the axion particle. We analyze data we collected from the Robert C. Byrd Green Bank Telescope in the L-band and the Effelsberg 100-m Telescope in the L-Band and S-band from a number of sources expected to produce bright signals of axion-photon conversion, including the Galactic Center of the Milky Way and the nearby isolated neutron stars RX J0720.4-3125 and RX J0806.4-4123. We find no evidence for axion DM and are able to set some of the strongest constraints to-date on the existence of axion DM in the highly-motivated mass range between ~5-11 $\mu$eV.

The extended excess towards the Galactic Center (GC) in gamma rays inferred from Fermi-LAT observations has been interpreted as being due to dark matter (DM) annihilation. Here, we perform new likelihood analyses of the GC and show that when including templates for the stellar galactic and nuclear bulges, the GC shows no significant detection of a DM annihilation template, even after generous variations in the Galactic diffuse emission models and a wide range of DM halo profiles. We include Galactic diffuse emission models with combinations of 3D inverse Compton maps, variations of interstellar gas maps, and a central source of electrons. For the DM profile, we include both spherical and ellipsoidal DM morphologies and a range of radial profiles from steep cusps to kiloparsec-sized cores, motivated in part by hydrodynamical simulations. Our derived upper limits on the dark matter annihilation flux place strong constraints on DM properties. In the case of the pure $b$-quark annihilation channel, our limits on the annihilation cross section are more stringent than those from the Milky Way dwarfs up to DM masses of $\sim$TeV, and rule out the thermal relic cross section up to $\sim$300 GeV. Better understanding of the DM profile, as well as the Fermi-LAT data at its highest energies, would further improve the sensitivity to DM properties.

The Fermi Large Area Telescope has observed an excess of ~GeV energy gamma rays from the center of the Milky Way, which may arise from near-thermal dark matter annihilation. Firmly establishing the dark matter origin for this excess is however complicated by challenges in modeling diffuse cosmic-ray foregrounds as well as unresolved astrophysical sources, such as millisecond pulsars. Non-Poissonian Template Fitting (NPTF) is one statistical technique that has previously been used to show that at least some fraction of the GeV excess is likely due to a population of dim point sources. These results were recently called into question by Leane and Slatyer (2019), who showed that a synthetic dark matter annihilation signal injected on top of the real Fermi data is not recovered by the NPTF procedure. In this work, we perform a dedicated study of the Fermi data and explicitly show that the central result of Leane and Slatyer (2019) is likely driven by the fact that their choice of model for the Galactic foreground emission does not provide a sufficiently good description of the data. We repeat the NPTF analyses using a state-of-the-art model for diffuse gamma-ray emission in the Milky Way and introduce a novel statistical procedure, based on spherical-harmonic marginalization, to provide an improved description of the Galactic diffuse emission in a data-driven fashion. With these improvements, we find that the NPTF results continue to robustly favor the interpretation that the Galactic Center excess is due, in part, to unresolved astrophysical point sources across the analysis variations that we have explored.

Searches for dark matter annihilation signals have been carried out in a number of target regions such as the Galactic Center and Milky Way dwarf spheroidal galaxies (dSphs), among a few others. Here we propose low surface brightness galaxies (LSBGs) asnovel targets for the indirect detection of dark matter emission. In particular, LSBGs are known to have very large dark matter contents and be less contaminated by extragalactic gamma-ray sources (e.g., blazars) compared to star forming galaxies. We report on an analysis that uses eight LSBGs (detected by Subaru Hyper Suprime-Cam survey data) with known redshifts to conduct a search for gamma-ray emission at the positions of these new objects in Fermi Large Area Telescope data. We found no excesses of gamma-ray emission and set constraints on the dark matter annihilation cross-section. We exclude (at the 95% C.L.) dark matter scenarios predicting a cross-section higher than 10^-23[cm^3/s] for dark matter particles of mass 10 GeV self-annihilating in the b_b channel. Although this constraint is weaker than the ones reported in recent studies using other targets, we note that in the near future, the number of detections of new LSBGs will increase by a few orders of magnitude. We forecast that with the use of the full catalog of soon-to-be-detected LSBGs the constraint will reach cross-section sensitivities of ~ 3*10^-25 [cm^3/s] for dark matter particles with masses less than 10 GeV.

We set conservative constraints on decaying dark matter particles with masses spanning a very wide range ($10^4-10^{16}$ GeV). For this we use multimessenger observations of cosmic-ray (CR) protons/antiprotons, electrons/positrons, neutrinos/antineutrinos and gamma rays. Focusing on decays into the $\bar{b}b$ channel, we simulate the spectra of dark matter yields by using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equations and the $\texttt{Pythia}$ package. We then propagate the CRs of dark matter origin till Earth by using the state-of-the-art numerical frameworks $\texttt{CRPropa}$, $\texttt{GALPROP}$ and $\texttt{HelMod}$ for the solution of the CR transport equation in the extragalactic, Galactic region and the heliosphere, respectively. Conservative limits are obtained by requiring that the predicted dark matter spectra at Earth be less than the observed CR spectra. Overall, we exclude dark matter lifetimes of $10^{28}$ s or shorter for all the masses investigated in this work. The most stringent constraints reach $10^{30}$ s for very heavy dark matter particles with masses in the range $10^{11}-10^{14}$ GeV.

We study the prospects for the Large Synoptic Survey Telescope (LSST) to find new dwarf spheroidal galaxies in the Milky Way. Adopting models of Milky-Way halo substructure and phenomenological prescriptions connecting subhalos and satellite galaxies, we obtain surface brightness distributions of $V$-band magnitude that lead us to predict that LSST will discover tens to hundreds of dwarf spheroidal galaxies above its sensitivity. The soon-to-be-discovered dwarfs will be interesting targets for indirect searches of dark matter annihilation yields. We forecast the distribution function of gamma-ray emission from dark matter annihilation in these objects, and discuss the detectability of these signals at both Fermi Large Area Telescope (LAT) and Cherenkov Telescope Array (CTA). By combining information from the predicted dwarf galaxies, we obtain an expected sensitivity to the annihilation cross section $\langle \sigma v \rangle$ of $10^{-26}$ cm$^3$ s$^{-1}$ (for dark matter particles of mass 10 GeV with Fermi-LAT) and $5\times 10^{-24}$ cm$^3$ s$^{-1}$ (for dark matter particles of mass 500 GeV with CTA). We find that the current uncertainties in the mass measurement of the Milky-Way halo are relatively minor compared with the Poisson errors associated to drawing the most promising dwarfs from the underlying flux distribution.

There is growing evidence that the Galactic Center Excess identified in the $\textit{Fermi}$-LAT gamma-ray data arises from a population of faint astrophysical sources. We provide compelling supporting evidence by showing that the morphology of the excess traces the stellar over-density of the Galactic bulge. By adopting a template of the bulge stars obtained from a triaxial 3D fit to the diffuse near-infrared emission, we show that it is detected at high significance. The significance deteriorates when either the position or the orientation of the template is artificially shifted, supporting the correlation of the gamma-ray data with the Galactic bulge. In deriving these results, we have used more sophisticated templates at low-latitudes for the $\textit{Fermi}$ bubbles compared to previous work and the three-dimensional Inverse Compton (IC) maps recently released by the ${\tt GALPROP}$ team. Our results provide strong constraints on Millisecond Pulsar (MSP) formation scenarios proposed to explain the excess. We find that an $\textit{admixture formation}$ scenario, in which some of the relevant binaries are $\textit{primordial}$ and the rest are formed $\textit{dynamically}$, is preferred over a primordial-only formation scenario at $7.6\sigma$ confidence level. Our detailed morphological analysis also disfavors models of the disrupted globular clusters scenario that predict a spherically symmetric distribution of MSPs in the Galactic bulge. For the first time, we report evidence of a high energy tail in the nuclear bulge spectrum that could be the result of IC emission from electrons and positrons injected by a population of MSPs and star formation activity from the same site.

Analyses of Fermi Gamma-Ray Space Telescope data have revealed a source of excess diffuse gamma rays towards the Galactic center that extends up to roughly $\pm$20 degrees in latitude. The leading theory postulates that this GeV excess is the aggregate emission from a large number of faint millisecond pulsars (MSPs). The electrons and positrons ($e^\pm$) injected by this population could produce detectable inverse-Compton (IC) emissions by up-scattering ambient photons to gamma-ray energies. In this work, we calculate such IC emissions using GALPROP. A triaxial three-dimensional model of the bulge stars obtained from a fit to infrared data is used as a tracer of the putative MSP population. This model is compared against one in which the MSPs are spatially distributed as a Navarro-Frenk-White squared profile. We show that the resulting spectra for both models are indistinguishable, but that their spatial morphologies have salient recognizable features. The IC component above $\sim$TeV energies carries information on the spatial morphology of the injected $e^\pm$. Such differences could potentially be used by future high-energy gamma-ray detectors such as the Cherenkov Telescope Array to provide a viable multiwavelength handle for the MSP origin of the GeV excess.

The cross-correlation study of the unresolved $\gamma$-ray background (UGRB) with galaxy clusters has a potential to reveal the nature of the UGRB. In this paper, we perform a cross-correlation analysis between $\gamma$-ray data by the Fermi Large Area Telescope (Fermi-LAT) and a galaxy cluster catalogue from the Subaru Hyper Suprime-Cam (HSC) survey. The Subaru HSC cluster catalogue provides a wide and homogeneous large-scale structure distribution out to the high redshift at $z=1.1$, which has not been accessible in previous cross-correlation studies. We conduct the cross-correlation analysis not only for clusters in the all redshift range ($0.1 < z < 1.1$) of the survey, but also for subsamples of clusters divided into redshift bins, the low redshift bin ($0.1 < z < 0.6$) and the high redshift bin ($0.6 < z < 1.1$), to utilize the wide redshift coverage of the cluster catalogue. We find the evidence of the cross-correlation signals with the significance of 2.0-2.3$\sigma$ for all redshift and low-redshift cluster samples. On the other hand, for high-redshift clusters, we find the signal with weaker significance level (1.6-1.9$\sigma$). We also compare the observed cross-correlation functions with predictions of a theoretical model in which the UGRB originates from $\gamma$-ray emitters such as blazars, star-forming galaxies and radio galaxies. We find that the detected signal is consistent with the model prediction.

An anomalous gamma-ray excess emission has been found in Fermi Large Area Telescope data covering the centre of the Galaxy. Several theories have been proposed for this `Galactic Centre Excess'. They include self-annihilation of dark matter particles, an unresolved population of millisecond pulsars, an unresolved population of young pulsars, or a series of burst events. Here we report on a new analysis that exploits hydrodynamical modelling to register the position of interstellar gas associated with diffuse Galactic gamma-ray emission. We find evidence that the Galactic Centre Excess gamma rays are statistically better described by the stellar over-density in the Galactic bulge and the nuclear stellar bulge, rather than a spherical excess. Given its non-spherical nature, we argue that the Galactic Centre Excess is not a dark matter phenomenon but rather associated with the stellar population of the Galactic bulge and nuclear bulge.

The extragalactic $\gamma$-ray background (EGB) arises from the accumulation of $\gamma$-ray emissions from resolved and unresolved extragalactic sources as well as diffuse processes. It is important to study the statistical properties of the EGB in the context of cosmological structure formation. Known astrophysical $\gamma$-ray sources such as blazars, star-forming galaxies, and radio galaxies are expected to trace the underlying cosmic matter density distribution. We explore the correlation of the EGB from Fermi-LAT data with the large-scale matter density distribution from the Subaru Hyper Suprime-Cam (HSC) SSP survey. We reconstruct an unbiased surface matter density distribution $\kappa$ at $z<1$ by applying weak-gravitational lensing analysis to the first-year HSC data. We then calculate the $\gamma - \kappa$ cross-correlation. Our measurements are consistent with a null detection, but a weak correlation is found at angular scales of 30-60 arcmin, especially when distant source galaxies at $z > 1$ are used for the lensing $\kappa$ reconstruction. The large-scale correlation suggests strong clustering of high-redshift $\gamma$-ray sources such as blazars. However, the inferred bias factor of $4-5$ is larger by about a factor of two than results from other clustering analyses. The final HSC data covering 1,400 squared degrees will play an essential role to determine accurately the blazar bias at $z > 0.5$.

Fermi Large Area Telescope data reveal an excess of GeV gamma rays from the direction of the Galactic Center and bulge. Several explanations have been proposed for this excess including an unresolved population of millisecond pulsars (MSPs) and self-annihilating dark matter. It has been claimed that a key discriminant for or against the MSP explanation can be extracted from the properties of the luminosity function describing this source population. Specifically, is the luminosity function of the putative MSPs in the Galactic Center consistent with that characterizing the resolved MSPs in the Galactic disk? To investigate this we have used a Bayesian Markov Chain Monte Carlo to evaluate the posterior distribution of the parameters of the MSP luminosity function describing both resolved MSPs and the Galactic Center excess. At variance with some other claims, our analysis reveals that, within current uncertainties, both data sets can be well fit with the same luminosity function.

We derive constraints on dark matter (DM) annihilation cross section and decay lifetime from cross-correlation analyses of the data from Fermi-LAT and weak lensing surveys that cover a wide area of $\sim660$ squared degrees in total. We improve upon our previous analyses by using an updated extragalactic $\gamma$-ray background data reprocessed with the Fermi Pass 8 pipeline, and by using well-calibrated shape measurements of about twelve million galaxies in the Canada-France-Hawaii Lensing Survey (CFHTLenS) and Red-Cluster-Sequence Lensing Survey (RCSLenS). We generate a large set of full-sky mock catalogs from cosmological $N$-body simulations and use them to estimate statistical errors accurately. The measured cross correlation is consistent with null detection, which is then used to place strong cosmological constraints on annihilating and decaying DM. For leptophilic DM, the constraints are improved by a factor of $\sim100$ in the mass range of O(1) TeV when including contributions from secondary $\gamma$ rays due to the inverse-Compton upscattering of background photons. Annihilation cross-sections of $\langle \sigma v \rangle \sim 10^{-23}\, {\rm cm}^3/{\rm s}$ are excluded for TeV-scale DM depending on channel. Lifetimes of $\sim 10^{25}$ sec are also excluded for the decaying TeV-scale DM. Finally, we apply this analysis to wino DM and exclude the wino mass around 200 GeV. These constraints will be further tightened, and all the interesting wino DM parameter region can be tested, by using data from future wide-field cosmology surveys.

The DUNE/LBNF program aims to address key questions in neutrino physics and astroparticle physics. Realizing DUNE's potential to reconstruct low-energy particles in the 10-100 MeV energy range will bring significant benefits for all DUNE's science goals. In neutrino physics, low-energy sensitivity will improve neutrino energy reconstruction in the GeV range relevant for the kinematics of DUNE's long-baseline oscillation program. In astroparticle physics, low-energy capabilities will make DUNE's far detectors the world's best apparatus for studying the electron-neutrino flux from a supernova. This will open a new window to unrivaled studies of the dynamics and neutronization of a star's central core in real time, the potential discovery of the neutrino mass hierarchy, provide new sensitivity to physics beyond the Standard Model, and evidence of neutrino quantum-coherence effects. The same capabilities will also provide new sensitivity to `boosted dark matter' models that are not observable in traditional direct dark matter detectors.

Excess GeV gamma rays from the Galactic Center (GC) have been measured with the Fermi Large Area Telescope (LAT). The presence of the GC excess (GCE) appears to be robust with respect to changes in the diffuse galactic background modeling. The three main proposals for the GCE are an unresolved population of millisecond pulsars (MSPs), outbursts of cosmic rays from the GC region, and self-annihilating dark matter (DM). The injection of secondary electrons and positrons into the interstellar medium (ISM) by an unresolved population of MSPs or DM annihilations can lead to observable gamma-ray emission via inverse Compton scattering or bremsstrahlung. Here we investigate how to determine whether secondaries are important in a model for the GCE. We develop a method of testing model fit which accounts for the different spatial morphologies of the secondary emission. We examine several models which give secondary emission and illustrate a case where a broadband analysis is not sufficient to determine the need for secondary emission.

The singlet-doublet fermion dark matter model (SDFDM) provides a good DM candidate as well as the possibility of generating neutrino masses radiatively. The search and identification of DM requires the combined effort of both indirect and direct DM detection experiments in addition to the LHC. Remarkably, an excess of GeV gamma rays from the Galactic Center (GCE) has been measured with the \textit{Fermi} Large Area Telescope (LAT) which appears to be robust with respect to changes in the diffuse galactic background modeling. Although several astrophysical explanations have been proposed, DM remains a simple and well motivated alternative. In this work, we examine the sensitivities of dark matter searches in the SDFDM scenario using $\textit{Fermi}$-LAT, CTA, IceCube/DeepCore, LUX, PICO and LHC with an emphasis on exploring the regions of the parameter space that can account for the GCE. We find that DM particles present in this model with masses close to $\sim 99$ GeV and $\sim (173-190)$ GeV annihilating predominantly into the $W^+W^-$ channel and $t\bar{t}$ channel respectively, provide an acceptable fit to the GCE while being consistent with different current experimental bounds. We also find that much of the obtained parameter space can be ruled out by future direct search experiments like LZ and XENON-1T, in case of null results by these detectors. Interestingly, we show that the most recent data by LUX is starting to probe the best fit region in the SDFDM model.

The High-Energy Stereoscopic System (HESS) has detected diffuse TeV emission correlated with the distribution of molecular gas along the Ridge at the Galactic Center. Diffuse, non-thermal emission is also seen by the Fermi large area telescope (Fermi-LAT) in the GeV range and by radio telescopes in the GHz range. Additionally, there is a distinct, spherically symmetric excess of gamma rays seen by Fermi-LAT in the GeV range. A cosmic ray flare, occurring in the Galactic Center, $10^4$ years ago has been proposed to explain the TeV Ridge. An alternative, steady-state model explaining all three data sets (TeV, GeV, and radio) invokes purely leptonic processes. We show that the flare model from the Galactic Center also provides an acceptable fit to the GeV and radio data, provided the diffusion coefficient is energy independent. However, if Kolmogorov-type turbulence is assumed for the diffusion coefficient, we find that two flares are needed, one for the TeV data (occurring approximately $10^4 $ years ago) and an older one for the GeV data (approximately $10^5$ years old). We find that the flare models we investigate do not fit the spherically symmetric GeV excess as well as the usual generalized Navarro-Frenk-White spatial profile, but are better suited to explaining the Ridge. We also show that a range of single-zone, steady-state models are able to explain all three spectral data sets. Large gas densities equal to the volumetric average in the region can be accommodated by an energy independent diffusion or streaming based steady-state model. Additionally, we investigate how the flare and steady-state models may be distinguished with future gamma-ray data looking for a spatial dependence of the gamma-ray spectral index.

The Fermi Large Area Telescope data appear to have an excess of gamma rays from the inner 150 pc of the Galactic Center (GC). The main explanations proposed for this are: an unresolved population of millisecond pulsars (MSPs), dark matter (DM) annihilation, and cosmic rays interacting with molecular clouds at the GC. In this conference proceeding article we highlight some of the cosmic-ray results of [1]. The MSPs and DM explanations were modeled with spatial templates well fitted by the square of a generalized Navarro-Frenk-White (NFW) profile with inner slope $\gamma=1.2$. The cosmic-ray option was modeled with a 20-cm continuum emission Galactic Ridge template. A template based on the HESS residuals were shown to give similar results. The gamma-ray excess was found to be best fit by a combination of the generalized NFW squared template and a Galactic Ridge template. We also found the spectra of each template was not significantly affected in the combined fit and is consistent with previous single template fits. Due to the ridge like morphology of the molecular gas at the GC, the cosmic-ray option was not able to explain all the excess gamma rays, indicating that DM annihilation and/or an unresolved population of MSPs is also needed.

The Fermi-LAT data appear to have an excess of gamma rays from the inner 150 pc of the Galactic Center. The main explanations proposed for this are: an unresolved population of millisecond pulsars (MSPs), dark matter (DM) annihilation, and nonthermal bremsstrahlung produced by a population of electrons interacting with neutral gas in molecular clouds. The first two options have spatial templates well fitted by the square of a generalized Navarro-Frenk-White (NFW) profile with inner slope $\gamma=1.2$. We model the third option with a 20-cm continuum emission Galactic Ridge template. A template based on the HESS residuals is shown to give similar results. The gamma-ray excess is found to be best fit by a combination of the generalized NFW squared template and a Galactic Ridge template. We also find the spectra of each template is not significantly affected in the combined fit and is consistent with previous single template fits. That is the generalized NFW squared spectrum can be fit by either of order 1000 unresolved MSPs or DM with mass around 30 GeV, a thermal cross section, and mainly annihilating to $b \bar{b}$ quarks. While the Galactic Ridge continues to have a spectrum consistent with a population of nonthermal electrons whose spectrum also provides a good fit to synchrotron emission measurements. We also show that the current DM fit may be hard to test, even with 10 years of Fermi-LAT data, especially if there is a mixture of DM and MSPs contributing to the signal, in which case the implied DM cross section will be suppressed.

Employing Fermi-LAT gamma ray observations, several independent groups have found excess extended gamma ray emission at the Galactic center (GC). Both, annihilating dark matter (DM) or a population of $\sim 10^3$ unresolved millisecond pulsars (MSPs) are regarded as well motivated possible explanations. However, there is significant uncertainties in the diffuse galactic background at the GC. We have performed a revaluation of these two models for the extended gamma ray source at the GC by accounting for the systematic uncertainties of the Galactic diffuse emission model. We also marginalize over point source and diffuse background parameters in the region of interest. We show that the excess emission is significantly more extended than a point source. We find that the DM (or pulsars population) signal is larger than the systematic errors and therefore proceed to determine the sectors of parameter space that provide an acceptable fit to the data. We found that a population of order a 1000 MSPs with parameters consistent with the average spectral shape of Fermi-LAT measured MSPs was able to fit the GC excess emission. For DM, we found that a pure $\tau^+\tau^-$ annihilation channel is not a good fit to the data. But a mixture of $\tau^+\tau^-$ and $b\bar{b}$ with a $\left<\sigma v\right>$ of order the thermal relic value and a DM mass of around 20 to 60 GeV provides an adequate fit.