Data for the paper "High-energy gamma-ray emission from memory-burdened primordial black holes" (Chianese 2025)

The files in this repository are provided to supplement the research article "High-energy gamma-ray emission from memory-burdened primordial black holes", Chianese 2025. They contain the constraints at 95% CL on the today abundance f_PBH of memory-burdended Primordial Black Holes (PBHs) placed by high-energy gamma-ray (Chianese 2025) and neutrino (Chianese et al. 2024) measurements, as a function of the PBHs mass at formation M_PBH and the memory-burden parameter k (the parameter q is fixed to 0.5). For any questions please feel free to contact me: m.chianese[AT]ssmeridionale.it All files are saved using the "numpy.save" command and can be imported by using "numpy.load" in Python. General The file k_values.npy gives the range of values of k for which the constraints have been computed, specificaly from k = 1.0 to k = 5.0 (number of values = 21). For each value of k, the lowest value for M_PBH corresponds to memory-burdened PBHs with a lifetime equal to the age of the Universe. Values of 1e2 for f_PBH mean that there is no constraint. Gamma-ray constraints Files: gamma_constraints_mPBH.npy - Values of M_PBH for different values of k.Shape 21 x 101 = (number of k values) x (number of masses) gamma_constraints_fPBH.npy - Constraint on f_PBH.Shape 21 x 101 x 9 = (number of k values) x (number of masses) x (different data sample) gamma_constraints_fPBH=1.npy - PBH masses correspoding to f_PBH = 1 as a function of k, i.e. the function M_PBH = M_PBH(k). The values for k are given in the file k_values.npy. Values NaN mean that there is no constraint.Shape 21 x 9 = (number of k values) x (different data sample) The data samples analyzed are: index 0 - Isotropic Gamma-Ray Background (IGRB) from Fermi-LAT index 1 - Galactic emission from LHAASO (inner region) index 2 - Galactic emission from LHAASO (outer region) index 3 - CASA-MIA index 4 - KASCADE index 5 - KASCADE-Grande index 6 - Pierre Auger Observatory (2015) index 7 - Pierre Auger Observatory hybrid (2016) index 8 - Telescope Array Surface Detector (2018) Neutrino constraints Files: neutrino_constraints_mPBH.npy - Values of M_PBH for different values of k.Shape 21 x 101 = (number of k values) x (number of masses) neutrino_constraints_fPBH.npy - Constraint on f_PBH.Shape 21 x 101 x 4 = (number of k values) x (number of masses) x (different data sample) neutrino_constraints_fPBH=1.npy - PBH masses correspoding to f_PBH = 1 as a function of k, i.e. the function M_PBH = M_PBH(k). The values for k are given in the file k_values.npy. Values NaN mean that there is no constraint.Shape 21 x 4 = (number of k values) x (different data sample) The data samples analyzed are: index 0 - IceCube HESE data (7.5 years) index 1 - IceCube EHE data (7 years) index 2 - IceCube-gen2 forecast (3 years) index 3 - GRAND200k forecast (3 years)