In the absence of a clear hint of dark matter (DM) signals in the GeV regime so far, heavy, $\mathcal{O}$(TeV) DM candidates are gradually earning more and more attention within the community. Among others, extra-dimensional \textit{brane-world} models may produce thermal DM candidates with masses up to 100 TeV. These heavy DM candidates could be detected with the next generation of very-high-energy gamma-ray observatories such as the Cherenkov Telescope Array (CTA). In this work, we study the sensitivity of CTA to branon DM via the observation of representative astrophysical DM targets, namely dwarf spheroidal galaxies. In particular, we focus on Draco and Sculptor, two well-known dwarfs visible from the Northern and Southern Hemisphere, respectively. For each of these targets, we simulated 300\,h of CTA observations and studied the sensitivity of both CTA-North and CTA-South to branon annihilations using the latest publicly available instrument response functions and most recent analysis tools. We computed annihilation cross section values needed to reach a $5\sigma$ detection as a function of the branon mass. Additionally, in the absence of a predicted DM signal, we obtained $2\sigma$ upper limits on the annihilation cross section. These limits lie $1.5-2$ orders of magnitude above the thermal relic cross section value, depending on the considered branon mass. Yet, CTA will allow to exclude a significant portion of the brane tension-mass parameter space in the $0.1-60$ TeV branon mass range, and up to tensions of $\sim 10$ TeV. More importantly, CTA will significantly enlarge the region already excluded by AMS and CMS, and will provide valuable complementary information to future SKA radio observations. We conclude that CTA will possess potential to constrain brane-world models and, more in general, TeV DM candidates.

Fermi-LAT unidentified sources (unIDs) have proven to be compelling targets for performing indirect dark matter (DM) searches. In a previous work, we found that among the 1235 unIDs in Fermi-LAT catalogs (3FGL, 2FHL and 3FHL) only 44 of those are DM subhalos candidates. We now implement a spectral analysis to test whether these remaining sources are compatible or not with DM origin. This analysis is executed using almost 10 years of Pass 8 Fermi-LAT data. None of the unIDs are found to significantly prefer DM-induced emission compared to other, more conventional, astrophysical sources. In order to discriminate between pulsar and DM sources, we developed a new method which is based on the source's spectral curvature, peak energy, and its detection significance. We also look for spatial extension, which may be a hint for a DM origin according to our N-body simulation studies of the subhalo population. In addition, we used Gaia DR2 data to search for a potential stellar counterpart to our best DM subhalo candidates and, although no firm associations could be found, one of them coincides with the Sagittarius stream. Finally, previous constraints on the DM annihilation cross section are updated with the new number of remaining DM subhalo candidates among unIDs. Our limits now rule out canonical thermal WIMPs up to masses of 10 GeV for $b\bar{b}$ and 20 GeV for $\tau^+\tau^-$ annihilation channels, in this way being as sensitive and complementary to those obtained from other targets and probes.

One possible and natural derivation from the collisionless cold dark matter (CDM) standard cosmological framework is the assumption of the existence of interactions between dark matter (DM) and photons or neutrinos. Such possible interacting dark matter (IDM) model would imply a suppression of small-scale structures due to a large collisional damping effect, even though the weakly interacting massive particle (WIMP) can still be the DM candidate. Because of this, IDM models can help alleviate alleged tensions between standard CDM predictions and observations at small mass scales. In this work, we investigate the properties of DM halo substructure or subhalos formed in a high-resolution cosmological N-body simulation specifically run within these alternative models. We also run its CDM counterpart, which allowed us to compare subhalo properties in both cosmologies. We show that, in the lower mass range covered by our simulation runs, both subhalo concentrations and abundances are systematically lower in IDM compared to the CDM scenario. Yet, as in CDM, we find that median IDM subhalo concentration values increase towards the innermost regions of their hosts for same mass subhalos. Also similarly to CDM, we find IDM subhalos to be more concentrated than field halos of the same mass. Our work has a direct application on studies aimed at the indirect detection of DM where subhalos are expected to boost the DM signal of their host halos significantly. From our results, we conclude that the role of halo substructure in DM searches will be less important in interacting scenarios than in CDM, but is nevertheless far from being negligible.

One of the predictions of the $\Lambda$CDM cosmological framework is the hierarchical formation of structure, giving rise to dark matter (DM) halos and subhalos. When the latter are massive enough they retain gas (i.e., baryons) and become visible. This is the case of the dwarf satellite galaxies in the Milky Way (MW). Below a certain mass, halos may not accumulate significant amounts of baryons and remain completely dark. However, if DM particles are Weakly Interacting Massive Particles (WIMPs), we expect them to annihilate in subhalos, producing gamma rays which can be detected with the Fermi satellite. Using the three most recent point-source Fermi Large Area Telescope (LAT) catalogs (3FGL, 2FHL and 3FHL), we search for DM subhalo candidates among the unidentified sources, i.e., sources with no firm association to a known astrophysical object. We apply several selection criteria based on the expected properties of the DM-induced emission from subhalos, which allow us to significantly reduce the list of potential candidates. Then, by characterizing the minimum detection flux of the instrument and comparing our sample to predictions from the Via Lactea II (VL-II) N-body cosmological simulation, we place conservative and robust constraints on the $\langle\sigma v\rangle-m_{DM}$ parameter space. For annihilation via the $\tau^+\tau^-$ channel, we put an upper limit of $4\times 10^{-26}~(5\times 10^{-25})~cm^3~s^{-1}$ for a mass of 10 (100) GeV. A critical improvement over previous treatments is the repopulation we made to include low-mass subhalos below the VL-II mass resolution. With more advanced subhalo candidate filtering the sensitivity reach of our method can potentially improve these constraints by a factor 3 (2) for $\tau^+\tau^-$ ($b \bar{b}$) channel.