Abstract The mass content of the universe is predominantly composed of Dark Matter (DM), whose unknown nature has to be determined via extensions of the Standard Model of particle physics (SM). A class of candidates frequently taken into account are Weakly Interactive Massive Particles (WIMPs): stable and massive particles coupled with SM via weak interaction. According to various models, WIMPs may annihilate producing SM particles, resulting in high-energy radiation observable in regions with high DM density, such as the Galactic core or dwarf galaxies. The main challenge lies in defining the spectral and spatial distribution of this diffuse source. This work introduces a new computational method to define the spatial distribution of the extensive gamma-ray flux resulting from WIMP pair annihilation, commonly referred to as the J-factor. Implemented in Python, this method is compatible with Gammapy, the official analysis tool for the Cherenkov Telescope Array Observatory (CTAO) Collaboration. We also present a detectability study of this signal with the CTAO in the Large Magellanic Cloud, considering an astrophysical background of various known sources. For WIMP masses around TeV, the CTAO could detect this radiation, shedding light on DM nature.
Ambrosino et al. (Tue,) studied this question.