Diffuse Extragalactic Gamma-ray Background Research Articles

Contribution of starburst nuclei to the diffuse gamma-ray and neutrino flux

ABSTRACT In nuclei of starburst galaxies (SBGs), the combination of an enhanced rate of supernova explosions and a high gas density suggests that cosmic rays (CRs) can be efficiently produced, and that most of them lose their energy before escaping these regions, resulting in a large flux of secondary products, including neutrinos. Although the flux inferred from an individual starburst region is expected to be well below the sensitivity of current neutrino telescopes, such sources may provide a substantial contribution to the diffuse neutrino flux measured by IceCube. Here, we compute the gamma-ray and neutrino flux due to SBGs based on a physical model of CR transport in a starburst nucleus, and accounting for the redshift evolution of the number density of starburst sources as inferred from recent measurements of the star formation rate. The model accounts for gamma-ray absorption both inside the sources and in the intergalactic medium. The latter process is responsible for electromagnetic cascades, which also contribute to the diffuse gamma-ray background at lower energies. The conditions for acceleration of CR protons up to energies exceeding $\sim 10 \, \rm PeV$ in starburst regions, necessary for the production of PeV neutrinos, are investigated in a critical way. We show that starburst nuclei can account for the diffuse neutrino flux above $\sim 200 \, \rm TeV$, thereby producing $\lesssim 40 {\rm { per\, cent}}$ of the extragalactic diffuse gamma-ray background. Below $\sim 200 \, \rm TeV$, the flux from starburst appears to be somewhat lower than the observed one, where both the Galactic contribution and the flux of atmospheric neutrinos may account for the difference.

A multi-component model for observed astrophysical neutrinos

Aims. We investigated the origin of observed astrophysical neutrinos. Methods. We propose a multi-component model for the observed diffuse neutrino flux. The model includes residual atmospheric backgrounds, a Galactic contribution (e.g., from cosmic ray interactions with gas), an extragalactic contribution from pp interactions (e.g., from starburst galaxies), and a hard extragalactic contribution from photo-hadronic interactions at the highest energies (e.g., from tidal disruption events or active galactic nuclei). Results. We demonstrate that this model can address the key problems of astrophysical neutrino data, such as the different observed spectral indices in the high-energy starting and through-going muon samples, a possible anisotropy due to Galactic events, the non-observation of point sources, and the constraint from the extragalactic diffuse gamma-ray background. Furthermore, the recently observed muon track with a reconstructed muon energy of 4.5 PeV might be interpreted as evidence for the extragalactic photo-hadronic contribution. We perform the analysis based on the observed events instead of the unfolded fluxes by computing the probability distributions for the event type and reconstructed neutrino energy. As a consequence, we give the probability of each of these astrophysical components on an event-to-event basis.

Minimal model for extragalactic cosmic rays and neutrinos

We aim to explain in a unified way the experimental data on ultrahigh energy cosmic rays (UHECR) and neutrinos, using a single source class and obeying limits on the extragalactic diffuse gamma-ray background (EGRB). If UHECRs only interact hadronically with gas around their sources, the resulting diffuse CR flux can be matched well to the observed one, providing at the same time large neutrino fluxes. However, air showers in the Earth's atmosphere induced by UHECRs with energies $E>3\times 10^{18}$ eV would reach in such a case their maxima too high. Therefore additional photo-hadronic interactions of UEHCRs close to the accelerator have to be present, in order to modify the nuclear composition of CRs in a relatively narrow energy interval. We include thus both photon and gas backgrounds, and combine the resulting CR spectra with the high-energy part of the Galactic CR fluxes predicted by the escape model. As result, we find a good description of experimental data on the total CR flux, the mean shower depth $X_\max$ and its width $RMS(X_\max)$ in the whole energy range above $E\simeq 10^{17}$ eV. The predicted high-energy neutrino flux matches IceCube measurements, while the contribution to the EGRB is of order 30%.

A determination of the gamma-ray flux and photon spectral index distributions of blazars from theFermi-LAT 3LAC

We present a determination of the distributions of gamma-ray photon flux – the so-called LogN–LogS relation – and photon spectral index for blazars, based on the third extragalactic source catalogue of the FermiGamma-ray Space Telescope's Large Area Telescope, and considering the photon energy range from 100 MeV to 100 GeV. The data set consists of the 774 blazars in the so-called Clean sample detected with a greater than approximately 7σ detection threshold and located above ±20° Galactic latitude. We use non-parametric methods verified in previous works to reconstruct the intrinsic distributions from the observed ones which account for the data truncations introduced by observational bias and includes the effects of the possible correlation between the flux and photon index. The intrinsic flux distribution can be represented by a broken power law with a high-flux power-law index of −2.43 ± 0.08 and a low-flux power-law index of −1.87 ± 0.10. The intrinsic photon index distribution can be represented by a Gaussian with mean of 2.62 ± 0.05 and width of 0.17 ± 0.02. We also report the intrinsic distributions for the subpopulations of BL Lac and FSRQ (Flat Spectrum Radio Quasar)-type blazars separately and these differ substantially. We then estimate the contribution of FSRQs and BL Lacs to the diffuse extragalactic gamma-ray background radiation. Under the simplistic assumption that the flux distributions probed in this analysis continue to arbitrary low flux, we calculate that the best-fitting contribution of FSRQs is 35 per cent and BL Lacs 17 per cent of the total gamma-ray output of the Universe in this energy range.

Gamma-ray luminosity function of BL Lac objects

The gamma-ray luminosity function (GLF) of BL Lac objects is constructed by using a sample of BL Lac objects with redshifts selected from the Second LAT AGN catalog. The GLFs of BL Lacs in the frame of the pure density evolution (PDE), the pure luminosity evolution (PLE), and the luminosity-dependent density (LDDE) models are determined by using the Markov Chain Monte Carlo (MCMC) technique, respectively. Our results suggest that the PDE model can give best description for BL Lac GLF based on the combination of constraints of model parameters and good fits to the observed data of Fermi-Large Area Telescope (LAT) BL Lacs, but other two models (PLE and LDDE) cannot be excluded. Based on our constructed GLFs, the contribution to the extragalactic diffuse gamma-ray background (EGRB) from BL Lacs is estimated, and ∼1–5 per cent of the EGRB in the 0.1–100 GeV band is found to come from unresolved BL Lacs (including the cascade emission). In addition, it is found that the BL Lac GLF is very different from flat spectrum radio quasar GLF and then the contribution of blazars to the EGRB should be estimated separately.

Constraints on diffuse gamma-ray emission from structure formation processes in the Coma cluster

We analyze 5-year (63 months) data of the Large Area Telescope on board Fermi satellite from the Coma galaxy cluster in the energy range between 100 MeV and 100 GeV. The likelihood analyses are performed with several templates motivated by models predicting gamma-ray emission due to structure formation processes. We find no excess emission and derive the most stringent constraints to date on the Coma cluster above 100 MeV, and on the tested scenarios in general. The upper limits on the integral flux range from 10^-10 to 10^-9 cm^-2s^-1, and are stringent enough to challenge different scenarios. We find that the acceleration efficiency of cosmic ray protons and electrons at shocks must be below approximately 15% and 1%, respectively. Additionally, we argue that the proton acceleration efficiency should be lower than 5% in order to be consistent with radio data. This, however, relays on magnetic field estimates in the cluster. In particular, this implies that the contribution to the diffuse extragalactic gamma-ray background due to gamma-rays from structure formation processes in clusters of galaxies is negligible, below 1%. Finally, we discuss future detectability prospects for Astro-H, Fermi after 10-yr of operation, and the Cherenkov Telescope Array.

The sensitivity of Cherenkov telescopes to dark matter and astrophysical anisotropies in the diffuse gamma-ray background

In this article, the capability of present (H.E.S.S., MAGIC, VERITAS) and planned (CTA) ground-based Cherenkov telescope systems for detecting angular anisotropies in the diffuse gamma-ray background is investigated. Following up on a study of the impact of instrumental characteristics (effective area, field of view, angular resolution, and background rejection efficiency), the first part examines the influence of different observational strategies, i.e. whether a single deep observation or a splitting over multiple shallow fields is preferred. In the second part, the sensitivity to anisotropies generated by self-annihilating dark matter is studied for different common dark matter models. We find that a relative contribution of ∼ 10% from dark matter annihilation to the extra-galactic diffuse gamma-ray background can be detected with planned configurations of CTA. In terms of the thermally-averaged self-annihilation cross section, the sensitivity of CTA corresponds to values below the thermal freeze-out expectation ⟨σv⟩ = 3 × 10−26 cm3 s−1 for dark matter particles lighter than ∼ 200 GeV. We stress the importance of constraining anisotropies from unresolved astrophysical sources with currently operating instruments already, as a novel and complementary method for investigating the properties of TeV sources.

Models of the contribution of blazars to the anisotropy of the extragalactic diffuse gamma-ray background

We study the relation between the measured anisotropies in the extragalactic diffuse gamma-ray background (DGRB) and the DGRB spectral intensity, and their potential origin from the unresolved blazar population. Using a physical-evolution model for blazars with a luminosity dependent density evolution (LDDE) and an observationally-determined luminosity-dependent blazar spectral energy distribution (SED), we find that blazars can account for the observed anisotropy of the DGRB consistent with their observed source-count distribution, but are in turn constrained in contributing significantly to the observed DGRB intensity. For the best-fit LDDE model accounting for the DGRB anisotropy and source-count distribution, blazars only contribute 5.9+2.1−1.0% (68% CL) of the DGRB intensity above 1 GeV. Requiring a higher fraction of the DGRB intensity contribution by blazars overproduces the DGRB anisotropy, and therefore blazars in the LDDE+SED-sequence model cannot simultaneously account for the DGRB intensity as well as anisotropy.We discuss the limitations of LDDE models. However, these models do not require the many unjustified and observationally-inconsistent simplifying assumptions — including a single power law for all blazar spectra and a simple broken power-law model for their source-count distribution — that are present in much previous work.

EXTRAGALACTIC VERY HIGH ENERGY GAMMA-RAY BACKGROUND

We study the origin of the extragalactic diffuse gamma-ray background using the data from the Fermi telescope. To estimate the background level, we count photons at high Galactic latitudes |b|>60 degrees. Subtracting photons associated to known sources and the residual cosmic ray and Galactic diffuse backgrounds, we estimate the Extragalactic Gamma-ray Background (EGB) flux. We find that the spectrum of EGB in the very-high-energy (VHE) band above 30 GeV follows the stacked spectrum of BL Lacs. LAT data reveal the positive (1+z)^k, 1<k<4 cosmological evolution of the BL Lac source population consistent with that of their parent population, FR I radio galaxies. We show that EGB at E>30 GeV could be completely explained by emission from unresolved BL Lacs if k~3.

FLUX AND PHOTON SPECTRAL INDEX DISTRIBUTIONS OFFERMI-LAT BLAZARS AND CONTRIBUTION TO THE EXTRAGALACTIC GAMMA-RAY BACKGROUND

We present a determination of the distributions of photon spectral index and gamma-ray flux - the so called LogN-LogS relation - for the 352 blazars detected with a greater than approximately seven sigma detection threshold and located above +/- 20 degrees Galactic latitude by the Large Area Telescope of the Fermi Gamma-ray Space Telescope in its first year catalog. Because the flux detection threshold depends on the photon index, the observed raw distributions do not provide the true LogN-LogS counts or the true distribution of the photon index. We use the non-parametric methods developed by Efron and Petrosian to reconstruct the intrinsic distributions from the observed ones which account for the data truncations introduced by observational bias and includes the effects of the possible correlation between the two variables. We demonstrate the robustness of our procedures using a simulated data set of blazars and then apply these to the real data and find that for the population as a whole the intrinsic flux distribution can be represented by a broken power law with high and low indexes of -2.37 +/- 0.13 and -1.70 +/- 0.26, respectively, and the intrinsic photon index distribution can be represented by a Gaussian with mean of 2.41 +/- 0.13 and width of 0.25 +/- 0.03. We also find the intrinsic distributions for the sub-populations of BL Lac and FSRQs type blazars separately. We then calculate the contribution of Fermi blazars to the diffuse extragalactic gamma-ray background radiation. Under the assumption that the flux distribution of blazars continues to arbitrarily low fluxes, we calculate the best fit contribution of all blazars to the total extragalactic gamma-ray output to be 60%, with a large uncertainty.

THE DIFFUSE GAMMA-RAY BACKGROUND FROM TYPE Ia SUPERNOVAE

The origin of the diffuse extragalactic gamma-ray background (EGB) has been intensively studied but remains unsettled. Current popular source candidates include unresolved star-forming galaxies, starburst galaxies, and blazars. In this paper we calculate the EGB contribution from the interactions of cosmic rays accelerated by Type Ia supernovae (SNe), extending earlier work which only included core-collapse SNe. We consider Type Ia events in star-forming galaxies, but also in quiescent galaxies that lack star formation. For star-forming galaxies, consistently including Type Ia events makes little change to the star-forming EGB prediction, so long as both SN types have the same cosmic-ray acceleration efficiencies in star-forming galaxies. Thus, our updated EGB estimate continues to show that star-forming galaxies can represent a substantial portion of the signal measured by Fermi. For quiescent galaxies, conversely, we find a wide range of possibilities for the EGB contribution. The dominant uncertainty we investigated comes from the mass in hot gas, which provides targets for cosmic rays; total gas masses are as yet poorly known, particularly at larger radii. Additionally, the EGB estimation is very sensitive to the cosmic-ray acceleration efficiency and confinement, especially in quiescent galaxies. In the most optimistic allowed scenarios, quiescent galaxies can be an important source of the EGB. In this case, star-forming galaxies and quiescent galaxies together will dominate the EGB and leave little room for other contributions. If other sources, such as blazars, are found to have important contributions to the EGB, then either the gas mass or cosmic-ray content of quiescent galaxies must be significantly lower than in their star-forming counterparts. In any case, improved Fermi EGB measurements will provide important constraints on hot gas and cosmic rays in quiescent galaxies.

ELMAG: A Monte Carlo simulation of electromagnetic cascades on the extragalactic background light and in magnetic fields

A Monte Carlo program for the simulation of electromagnetic cascades initiated by high-energy photons and electrons interacting with extragalactic background light (EBL) is presented. Pair production and inverse Compton scattering on EBL photons as well as synchrotron losses and deflections of the charged component in extragalactic magnetic fields (EGMF) are included in the simulation. Weighted sampling of the cascade development is applied to reduce the number of secondary particles and to speed up computations. As final result, the simulation procedure provides the energy, the observation angle, and the time delay of secondary cascade particles at the present epoch. Possible applications are the study of TeV blazars and the influence of the EGMF on their spectra or the calculation of the contribution from ultrahigh energy cosmic rays or dark matter to the diffuse extragalactic gamma-ray background. As an illustration, we present results for deflections and time delays relevant for the derivation of limits on the EGMF. Program summaryProgram title: ELMAG 1.01Catalogue identifier: AELK_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AELK_v1_0.htmlProgram obtainable from: CPC Program Library, Queenʼs University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 4290No. of bytes in distributed program, including test data, etc.: 459 305Distribution format: tar.gzProgramming language: Fortran 95Computer: Any computer with Fortran 95 compilerOperating system: Any system with Fortran 95 compilerHas the code been vectorised or parallelised?: Yes, using the MPI versionRAM: 4 MbytesClassification: 11.3, 11.4Nature of problem: Calculation of secondaries produced by electromagnetic cascades on the extragalactic background light (EBL).Solution method: Monte Carlo simulation of pair production and inverse Compton scattering on EBL photons; two parametrisations from Ref. [1] can be chosen as EBL; weighted sampling of the cascading secondaries; recording of energy, observation angle and time delay of secondary particles at the present epoch.Restrictions: Deflections and time delays are calculated in the small-angle approximation.Running time: 400 seconds for 103 photons injected at redshift z=0.2 with energy E=100 TeV using one Intel(R) Core(TM) i7 CPU with 2.8 GHz.

Contribution of blazars to the extragalactic diffuse gamma-ray background and their future spatial resolution

We examine the constraints on the luminosity-dependent density evolution model for the evolution of blazars given the observed spectrum of the diffuse gamma-ray background (DGRB), blazar source-count distribution, and the blazar spectral energy distribution sequence model, which relates the observed the blazar spectrum to its luminosity. We show that the DGRB observed by the Large Area Telescope (LAT) aboard the Fermi Gamma Ray Space Telescope can be produced entirely by gamma-ray emission from blazars and nonblazar active galactic nuclei, and that our blazar evolution model is consistent with and constrained by the spectrum of the DGRB and flux source-count distribution function of blazars observed by Fermi-LAT. Our results are consistent with previous work that used EGRET spectral data to forecast the Fermi-LAT DGRB. The model includes only three free parameters, and forecasts that >~ 95% of the flux from blazars will be resolved into point sources by Fermi-LAT with 5 years of observation, with a corresponding reduction of the flux in the DGRB by a factor of ~2 to 3 (95% confidence level), which has implications for the Fermi-LAT's sensitivity to dark matter annihilation photons.

A cross-correlation study of the Fermi-LAT γ-ray diffuse extragalactic signal

In this work, starting from 21 months of data from the Fermi-Large Area Telescope, we derive maps of the residual isotropic gamma-ray emission, a relevant fraction of which is expected to be contributed by the extragalactic diffuse gamma-ray background. We compute the angular two-point auto-correlation function of the residual Fermi-LAT maps at energies E>1GeV, E>3GeV and E>30GeV well above the Galactic plane and find no significant correlation signal. This is, indeed, what is expected if the EGB were contributed by BL Lacertae, Flat Spectrum Radio Quasars or star-forming galaxies, since, in this case, the predicted signal is very weak. Then, we search for the Integrated Sachs-Wolfe signature by cross-correlating the Fermi-LAT maps with the WMAP7-Cosmic Microwave Background map. We find a cross-correlation consistent with zero, even though the expected signal is larger than that of the EGB auto-correlation. Finally, in an attempt to constrain the nature of the gamma-ray background we cross-correlate the Fermi-LAT maps with the angular distributions of objects that may contribute to the EGB: QSOs in the SDSS-DR6 catalog, NVSS galaxies, 2MASS galaxies and LRG in the SDSS catalog. The cross-correlation is always consistent with zero, in agreement with theoretical expectations, but we find (with low statistical significance) some interesting features that may indicate that some specific classes of objects contribute to the EGB. A chi2 analysis confirms that the correlation properties of the 21-month data do not provide strong constraints of the EGB origin. However, the results suggest that the situation will significantly improve with the 5- and 10-year Fermi-LAT data. The future EGB analysis will then allow placing significant constraints on the nature of the EGB and might provide in addition a detection of the ISW signal. (Abridged)

CONSTRAINING THE EMISSIVITY OF ULTRAHIGH ENERGY COSMIC RAYS IN THE DISTANT UNIVERSE WITH THE DIFFUSE GAMMA-RAY EMISSION

Ultra-high cosmic rays (UHECRs) with energies >10^19 eV emitted at cosmological distances will be attenuated by cosmic microwave and infrared background radiation through photohadronic processes. Lower energy extra-galactic cosmic rays (~10^18-10^19 eV) can only travel a linear distance smaller than ~Gpc in a Hubble time due to the diffusion if the extra-galactic magnetic fields are as strong as nano Gauss. These prevent us from directly observing most of the UHECRs in the universe, and thus the observed UHECR intensity reflects only the emissivity in the nearby universe within hundreds of Mpc. However, UHECRs in the distant universe, through interactions with the cosmic background photons, produce UHE electrons and gamma-rays that in turn initiate electromagnetic cascades on cosmic background photons. This secondary cascade radiation forms part of the extragalactic diffuse GeV-TeV gamma-ray radiation and, unlike the original UHECRs, is observable. Motivated by new measurements of extragalactic diffuse gamma-ray background radiation by Fermi/LAT, we obtained upper limits placed on the UHECR emissivity in the distant universe by requiring that the cascade radiation they produce not exceed the observed levels. By comparison with the gamma-ray emissivity of candidate UHECR sources (such as GRBs and AGNs) at high-redshifts, we find that the obtained upper limit for a flat proton spectrum is ~10^1.5 times larger than the gamma-ray emissivity in GRBs and ~10 times smaller than the gamma-ray emissivity in BL Lac objects. In the case of iron nuclei composition, the derived upper limit of the UHECR emissivity is a factor of 3-5 times higher. Robust upper limit on the cosmogenic neutrino flux is further obtained, which is marginally reachable by the Icecube detector and the next-generation detector JEM-EUSO.

FUELING LOBES OF RADIO GALAXIES: STATISTICAL PARTICLE ACCELERATION AND THE EXTRAGALACTIC γ-RAY BACKGROUND

The recent discovery of the gamma-ray emission from the lobes of the closest radio galaxy Centaurus A by Fermi implies the presence of high-energy electrons at least up to gamma ~ 10^5 - 10^6. These high-energy electrons are required to interpret the observed gamma-ray radiation in terms of inverse Compton emission off the cosmic microwave background (IC/CMB); the widely accepted scenario to describe the X-ray emission of radio galaxy lobes. In this Letter, we consider the giant radio lobes of FR II radio galaxies showing that it is possible to maintain electrons at energies gamma ~ 10^5 - 10^6, assuming an acceleration scenario (driven by turbulent magnetic fields) that compensates the radiative losses. In addition, we consider the contribution to the diffuse extragalactic gamma-ray background due to the IC/CMB emission of FR IIs' lobes showing its relevance in the keV to MeV energy range.

Extragalactic diffuse gamma-ray emission at high energies

New measurements with the Fermi Large Area Telescope extend our knowledge of the extragalactic diffuse gamma-ray background and may help resolve the question of its origins.

Contribution of inverse Compton scattering to the diffuse extragalactic gamma-ray background from annihilating dark matter

In addition to gamma-rays, dark matter annihilation products can include energetic electrons which inverse Compton scatter with the cosmic microwave background to produce a diffuse extragalactic background of gamma-rays and X-rays. In models in which the dark matter particles annihilate primarily to electrons or muons, the measurements of EGRET and COMPTEL can provide significant constraints on the annihilation cross section. The Fermi Gamma-Ray Space Telescope will likely provide an even more stringent test of such scenarios.

Detecting gamma-ray anisotropies from decaying dark matter: Prospects for Fermi LAT

Decaying dark matter particles could be indirectly detected as an excess over a simple power law in the energy spectrum of the diffuse extragalactic gamma-ray background. Furthermore, since the Earth is not located at the center of the Galactic dark matter halo, the exotic contribution from dark matter decay to the diffuse gamma-ray flux is expected to be anisotropic, offering a complementary method for the indirect search for decaying dark matter particles. In this paper we discuss in detail the expected dipole-like anisotropies in the dark matter signal, taking also into account the radiation from inverse Compton scattering of electrons and positrons from dark matter decay. A different source for anisotropies in the gamma-ray flux are the dark matter density fluctuations on cosmic scales. We calculate the corresponding angular power spectrum of the gamma-ray flux and comment on observational prospects. Finally, we calculate the expected anisotropies for the decaying dark matter scenarios that can reproduce the electron/positron excesses reported by PAMELA and the Fermi LAT, and we estimate the prospects for detecting the predicted gamma-ray anisotropy in the near future.

Constraints on leptonically annihilating dark matter from reionization and extragalactic gamma background

The PAMELA, Fermi and HESS experiments (PFH) have shown anomalous excesses in the cosmic positron and electron fluxes. A very exciting possibility is that those excesses are due to annihilating dark matter (DM). In this paper we calculate constraints on leptonically annihilating DM using observational data on diffuse extragalactic gamma-ray background and measurements of the optical depth to the last-scattering surface, and compare those with the PFH favored region in the m_{DM} - <\sigma_A v> plane. Having specified the detailed form of the energy input with PYTHIA Monte Carlo tools we solve the radiative transfer equation which allows us to determine the amount of energy being absorbed by the cosmic medium and also the amount left over for the diffuse gamma background. We find that the constraints from the optical depth measurements are able to rule out the PFH favored region fully for the \tau^{-}+\tau^{+} annihilation channel and almost fully for the \mu^{-}+\mu^{+} annihilation channel. It turns out that those constraints are quite robust with almost no dependence on low redshift clustering boost. The constraints from the gamma-ray background are sensitive to the assumed halo concentration model and, for the power law model, rule out the PFH favored region for all leptonic annihilation channels. We also find that it is possible to have models that fully ionize the Universe at low redshifts. However, those models produce too large free electron fractions at z > ~100 and are in conflict with the optical depth measurements. Also, the magnitude of the annihilation cross-section in those cases is larger than suggested by the PFH data.