The $\Lambda$CDM cosmological framework predicts the existence of thousands of subhalos in our own Galaxy not massive enough to retain baryons and become visible. Yet, some of them may shine in gamma rays provided that the dark matter (DM) is made of weakly interacting massive particles (WIMPs), that would self-annihilate and would appear as unidentified gamma-ray sources (unIDs) in gamma-ray catalogs. Indeed, unIDs have proven to be competitive targets for DM searches with gamma rays. In this work, we focus on the three high-latitude ($|b|\geq 10^\circ$) sources present in the 2HWC catalog of the High Altitude Water Cherenkov (HAWC) observatory with no associations at other wavelenghts. Indeed, only one of these sources, 2HWC J1040+308, is found to be above the HAWC detection threshold when considering 760 days of data, a factor 1.5 more exposure time than in the original 2HWC catalog. Other instruments such as Fermi-LAT or VERITAS at lower energies do not detect this source. Also, this unID is reported as spatially extended, making it even more interesting in a DM search context. While waiting for more data that may shed further light on the nature of this source, we set competitive upper limits on the annihilation cross section by comparing this HAWC unID to expectations based on state-of-the-art N-body cosmological simulations of the Galactic subhalo population. We find these constraints to be particularly competitive for heavy WIMPs, i.e., masses above $\sim 25$ (40) TeV in the case of the $b\bar{b}$ ($\tau^+\tau^-$) annihilation channel, reaching velocity-averaged cross section values of $2\cdot10^{-25}$ ($5\cdot10^{-25}$) $cm^3s^{-1}$. Although far from the thermal relic cross section value, the obtained limits are independent and nicely complementary to those from radically different DM analyses and targets, demonstrating again the high potential of this DM search approach.

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.