CLUMPY: Jeans analysis, $\gamma$-ray and neutrino fluxes from dark matter (sub-)structures

TL;DR

CLUMPY is an updated code that calculates astrophysical J-factors for dark matter structures, now including substructure boosts, triaxiality, and a new Jeans analysis module for deriving density profiles from kinematic data.

Contribution

The new CLUMPY version introduces advanced features like substructure boost factors, halo triaxiality, and a Jeans analysis module, enhancing dark matter flux predictions and profile derivations.

Findings

Enabled halo-to-halo concentration scatter modeling.

Produced gamma-ray and neutrino skymaps with HEALPix.

Integrated Jeans analysis for dark matter density profiles.

Abstract

We present an update of the CLUMPY code for the calculation of the astrophysical J-factors (from dark matter annihilation/decay) for any Galactic or extragalactic dark matter halo including substructures: halo-to-halo concentration scatter may now be enabled, boost factors can include several levels of substructures, and triaxiality is a new option for dark matter haloes. This new version takes advantage of the CFITSIO and HEALPix libraries to propose FITS output maps using the HEALPix pixelisation scheme. Skymaps for $\gamma$-ray and $\nu$ signals from generic annihilation/decay spectra are now direct outputs of CLUMPY. Making use of HEALPix routines, smoothing by a user-defined instrumental Gaussian beam and computing the angular power spectrum of the maps is now possible. In addition to these improvements, the main novelty is the implementation of a Jeans analysis module, to obtain dark matter density profiles from kinematic data in relaxed spherical systems (e.g., dwarf spheroidal galaxies). The code is also interfaced with the GreAT toolkit designed for Markov Chain Monte Carlo analyses, from which probability density functions and credible intervals can be obtained for velocity dispersions, dark matter profiles, and J-factors.