Data for the paper "Breakdown of Hawking Evaporation opens new Mass Window for Primordial Black Holes as Dark Matter Candidate" (Thoss et. al. 2024)

The files provided in this respository are provided to supplement the research article "Light Primordial Black Holes as Dark Matter due to the breakdown of Hawking Evaporation", Valentin Thoss et. al. 2024.They are sufficient to reproduce the figures and results presented in the paper.

For any questions please feel free to contact me: [e-mail]

All files are saved using the "numpy.save" command and can be imported by using "numpy.load" in Python

General:

k.npy gives the range of values of k for which the data has been computed

m_init.npy gives the range of values of M_0 for which the data has been computed

Values of "Inf" for f_PBH imply that there is no constraint

Values of "0.0" for f_PBH imply that PBH can not exist today because they would have fully evaporated

Values of "Inf" for k_PBH imply that PBH are excluded from being 100% of the dark matter independent of k for the given mass

Values of "0.0" for k_PBH imply that PBH can be 100% of the dark matter even for k=0.0

Galactic gamma ray emission:

Files:

f_galactic_primary.npy - Constraint on f_PBH from primary galactic gamma ray emission

f_galactic_secondary.npy - Constraint on f_PBH from secondary galactic gamma ray emission

k_galactic_primary.npy - Value of k for which f_PBH=1, obtained from primary galactic gamma ray emission

k_galactic_secondary.npy - Value of k for which f_PBH=1, obtained from secondary galactic gamma ray emissio

The arrays for f_PBH have the shape 5x300x300 = (number of observational data)x(number of masses)x(number of k values)

The first index iterates the observational data used in the following order:

INTEGRAL

COMPTEL

EGRET

FERMI-LAT

LHAASO

The second index iterates over the black hole masses, given in m_init.npy

The third index iterates over the values of k, given in k.npy

The arrays for k have the shape 5x300 = (number of observational data)x(number of masses)

The first index iterates the observational data in the order given above

The second iterates over the black hole masses, given in m_init.npy

Extragalactic gamma ray emission:

Files:

f_extragalactic_primary.npy - Constraint on f_PBH from primary extragalactic gamma ray emission

f_extragalactic_secondary.npy - Constraint on f_PBH from secondary extragalactic gamma ray emission

k_extragalactic_primary.npy - Value of k for which f_PBH=1, obtained from primary extragalactic gamma ray emission

k_extragalactic_secondary.npy - Value of k for which f_PBH=1, obtained from secondary extragalactic gamma ray emissio

The arrays for f_PBH have the shape 5x300x300 = (number of observational data)x(number of masses)x(number of k values)

The first index iterates the observational data used in the following order:

HEAO

COMPTEL

EGRET

FERMI-LAT

LHAASO

The second index iterates over the black hole masses, given in m_init.npy

The third index iterates over the values of k, given in k.npy

The arrays for k have the shape 5x300 = (number of observational data)x(number of masses)

The first index iterates the observational data in the order given above

The second iterates over the black hole masses, given in m_init.npy

CMB anisotropies:

Files:

f_cmb.npy - Constraint on f_PBH from CMB anisotropies

k_cmb.npy - Value of k for which f_PBH=1, obtained from CMB anisotropies

The array for f_PBH has the shape 300x300 = (number of masses)x(number of k values)

The first index iterates over the black hole masses, given in m_init.npy

The second index iterates over the values of k, given in k.npy

The array for k has the length 300 = (number of masses)

The index iterates over the black hole masses, given in m_init.npy

BBN:

Files:

f_bbn.npy - Constraint on f_PBH from BBN

k_bbn.npy - Value of k for which f_PBH=1, obtained from CMB anisotropies

The array for f_PBH has the shape 300x300 = (number of masses)x(number of k values)

The first index iterates over the black hole masses, given in m_init.npy

The second index iterates over the values of k, given in k.npy

The array for k has the length 300 = (number of masses)

The index iterates over the black hole masses, given in m_init.npy