Planned intervention: On Wednesday April 3rd 05:30 UTC Zenodo will be unavailable for up to 2-10 minutes to perform a storage cluster upgrade.
Published June 7, 2023 | Version 1.0
Dataset Open

Data and simulation files for "Constraints on the intergalactic magnetic field using Fermi-LAT and H.E.S.S. blazar observations"

  • 1. University of Southern Denmark

Description

In this repository, we provide data files in connection to our paper “Constraints on the intergalactic magnetic field using Fermi-LAT and H.E.S.S. blazar observations” accepted for publication in the Astrophysical Journals and soon available on Arxiv.

In the publication, we perform a joint analysis of observations of five blazars with the Fermi Large Area Telescope (LAT) and the High Energy Stereoscopic System (H.E.S.S.) in order to search for signatures of a gamma-ray halo around these sources. The non-detection of such extended emission allows us to place lower limits on the intergalactic magnetic field (IGMF).

In this repository, we provide our data analysis products of both H.E.S.S. and LAT data for the case when a template for the halo flux is not included in the data. Furthermore, we provide files that contain the log likelihood profiles as functions of the IGMF in case the halo emission is included. Lastly, we also provide our template files for the halo, generated with CRPropa 3.

Below, we provide minimal code examples to demonstrate how to read in the specific files.

H.E.S.S. observational results

We provide the best-fit spectral parameters as well as the flux points (spectral energy distribution; SED) for the H.E.S.S. observations of the five blazars under consideration. The corresponding files are:

  • hess_fit_result_*.fits which contain the best-fit parameters,
  • hess_sed_file_*.fits which contain the flux points.

In the file names above, the '*' should be replaced with a the corresponding source name, e.g. 1ES0229+200. The files can be read in using astropy:

from astropy.table import Table

src = "1ES0229+200"

best_fit_pars = Table.read("hess_fit_result_1ES0229+200.fits")
sed = Table.read("hess_sed_file_1ES0229+200.fits")

Fermi observational results

For Fermi-LAT, we provide the SED files as well as the best-fit models for the region of interests. These files are called:

  • fermi_avg_file_*.npy provides the best-fit ROI model
  • fermi_sed_file_*.npy provides the SED.

Both of these files are generated with fermipy and can be read-in the following way:

import numpy as np

# first a little helper function since the 
# fermipy analysis was run under python 2.7
def convert(data):
    if isinstance(data, bytes):  return data.decode('ascii')
    if isinstance(data, dict):   return dict(map(convert, data.items()))
    if isinstance(data, tuple):  return map(convert, data)
    return data

# Load the ROI fit 
roi_fit_file = "fermi_avg_file_1ES0229+200.npy"
roi_fit = np.load(avg_file, allow_pickle=True, encoding="latin1").flat[0]

# if you want to inspect the dictionaries in python 3, you need to run the convert function.
# For example, to inspect the central source of the ROI
# you would first get the source name 
src_fgl_name = roi_fit['config']['selection']['target']

# and then you can get the dictionary for the central source
src_dict = convert(roi_fit['sources'])[src_fgl_name]

# Load the SED
sed_file = "fermi_sed_file_1ES0229+200.npy"
sed = np.load(sed_file, allow_pickle=True, encoding='latin1').flat[0]


# to plot the SED, you can use the SEDPlotter class from fermipy
from fermipy.plotting import SEDPlotter

SEDPlotter.plot_sed(sed)

Likelihood profiles

The likelihood profiles as function of the IGMF strengths are provided in the files logl_profile_*_*yr.npz. Their are provided for all five sources and all tested blazar activity times of 10, 104, and 107 years. They can be read in with the following code snippet:

import numpy as np

logl = dict(np.load("logl_profile_1ES0229+200_1.0e+07yr.npz"))
b_fields = np.array([1.00000e-16, 3.16228e-16, 1.00000e-15,
                     3.16228e-15, 1.00000e-14,
                     3.16228e-14, 1.00000e-13])

for k, v in logl.items():
    print(k,v)

As the print command shows, the python dictionary contains 3 entries: "fermi_only" are the likelihood values for the Fermi data as a function of magnetic field, "combined" are the likelihood values from Fermi and H.E.S.S. combined, and "ps" is the likelihood value of the Fit without halo to the H.E.S.S. data only.

Halo simulations

Lastly, we also provide the output simulations files from CRPropa. For details how the simulations were run, please consult the accompanying paper, in particular Section 3.1 and Appendix C. For each source redshift, a tar file is provided, which in itself contains 7 hdf5 files with the simulation outputs for each tested magnetic field strength. The name of the files is casc_file_z*.tar.gz. After unpacking the files, they can be read in with your favorite hdf5 library; in python you would need to install h5py. We recommend that you check out this github repository which provides an advanced python wrapper for CRPropa and functions to read in the files. In particular, you can use this function to read in the files. It also writes a new hdf5 file with parallel transport applied. The written data is also returned together with the configuration dictionary.

from simCRpropa.cascmaps import stack_results_lso

data, config = stack_results_lso("casc_file_z0.140_B1.00e-16.hdf5", 
                                 "casc_file_z0.140_B1.00e-16_theta_obs0.0.hdf5"
                                 )

You can provide arbitrary angles between the observer and the jet angles using the theta_obs keyword. Note, however, that the simulations used a jet opening angle of 3 degrees and going beyond that value will return zero halo photons.

Files

Files (667.0 MB)

Name Size Download all
md5:ef152eafc2aa5d443389edc0cf0a7cd3
128.4 MB Download
md5:70f4f50ebbb401dd8eaa7be8521d5c30
133.5 MB Download
md5:733430908ef51c3201682b95c77a9aa6
134.3 MB Download
md5:fd85e20b4c43bb0e636974e06d57fc3b
135.0 MB Download
md5:de83cc4ba8b2c00054748b3b7a6112d1
134.8 MB Download
md5:aff48502c36aab9bc2c1a89d840ab499
96.8 kB Download
md5:c199df647b16d45442938be622a9e734
63.2 kB Download
md5:4ddc5962f6f065fa68fa80aaeececaa8
83.3 kB Download
md5:b082d3b1db4cf827b6e269656447e145
96.4 kB Download
md5:67c53971007cd304422413581340a35e
147.1 kB Download
md5:3ab86b6d9bb32348e4d75ea3f0fc99b0
127.8 kB Download
md5:1c87954db667cbb229641653a463b20c
16.4 kB Download
md5:b49ab44d39ea9715ac36b06244908dd6
32.8 kB Download
md5:a3069caa0a3c6349e16204420fdca906
32.8 kB Download
md5:1c66514630ecdd7b6bbd567343eb6252
32.8 kB Download
md5:452015c59acdeadf16f875948e0337a3
32.8 kB Download
md5:57f2efb6be75859f01aa46cb5fe12d98
32.8 kB Download
md5:bd07d931b44a71a8f8eef5a032d3e3a2
8.6 kB Download
md5:57e570649441c05dbb5656342628564e
8.6 kB Download
md5:c6300ffa53d49ada5571f8113b9aec68
8.6 kB Download
md5:711d0d4d7d72ca9849ec2aca3e7d343c
8.6 kB Download
md5:49318238b73371d1cb06470904780872
8.6 kB Download
md5:e774d041dc7fd43b26d3d6748d7f984a
17.3 kB Download
md5:8473f63df23593a31b3a42d0093cb3da
17.3 kB Download
md5:50b65759d74158bdd516d6b54f4b32da
17.3 kB Download
md5:0995916869e199d00cbc247b2db6f968
17.3 kB Download
md5:b57f30ee40a6c1ccf1b8799ee03ddc3e
17.3 kB Download
md5:067c6d3e1850d7d172cc88e4ba774485
878 Bytes Download
md5:d567892639b99a0c658f4198ce30558a
878 Bytes Download
md5:ca3426ca5e1bc60060a08132dd85ea7f
878 Bytes Download
md5:7e4a364f8b35b7d34d3775e2cd1e69a4
878 Bytes Download
md5:598d54e45c9f87e962757d45015bfece
878 Bytes Download
md5:94ded16efef8b659fbc647964394f002
878 Bytes Download
md5:9a8aa76b2b34866d81d7058475875766
878 Bytes Download
md5:5dcf96109ec63c6c6615372d8cea647d
878 Bytes Download
md5:10a8355f08d34a6e15111788b5217b2e
878 Bytes Download
md5:47cfbe739a9ec631f3f14cf836ebfa83
878 Bytes Download
md5:6279fbd6fbeabba64ce0e7ab420aa10d
878 Bytes Download
md5:62be66490e5eb177302724d6b5d46016
878 Bytes Download
md5:a5c6ccf2a8ea7eb1ea95509b3eb21c83
878 Bytes Download
md5:8f954aa85fa00e2771ed070b794b19b0
878 Bytes Download
md5:bc3589d669a07150a945027cf5b81c32
878 Bytes Download

Additional details

Funding

AxionDM – Searching for axion and axion-like-particle dark matter in the laboratory and with high-energy astrophysical observations 948689
European Commission