Gamma-ray Luminosity Function Research Articles

High-energy Emission Component, Population, and Contribution to the Extragalactic Gamma-Ray Background of Gamma-Ray-emitting Radio Galaxies

In this study, we systematically studied the X-ray to GeV gamma-ray spectra of 61 Fermi Large Area Telescope detected radio galaxies. We found an anticorrelation between peak frequency and peak luminosity in the high-energy spectral component of radio galaxies, similar to blazars. With this sample, we also constructed a gamma-ray luminosity function (GLF) of gamma-ray-loud radio galaxies. We found that blazar-like GLF shapes can reproduce their redshift and luminosity distribution, but the log N–log S relation prefers models with more low-z radio galaxies. Utilizing our latest GLF, the contribution of radio galaxies to the extragalactic gamma-ray background is found to be 1%–10%. We further investigated the nature of gamma-ray-loud radio galaxies. Compared to radio or X-ray flux-limited radio galaxy samples, the gamma-ray-selected sample tends to lack high radio power galaxies like FR II radio galaxies. We also found that only ∼10% of radio galaxies are GeV gamma-ray loud. Radio galaxies may contribute to the cosmic MeV gamma-ray background comparable to blazars if gamma-ray-quiet radio galaxies have X-ray to gamma-ray spectra like Cen A, with a small gamma-ray-to-X-ray flux ratio.

Testing gamma-ray models of blazars in the extragalactic sky

The global contribution of unresolved gamma-ray point sources to the extragalactic gamma-ray background has been recently measured down to gamma-ray fluxes lower than those reached with standard source detection techniques, and by employing the statistical properties of the observed gamma-ray counts. We investigate and exploit the complementarity of the information brought by the one-point statistics of photon counts (using more than 10 years of Fermi-LAT data) and by the recent measurement of the angular power spectrum of the unresolved gamma-ray background (based on 8 years of Fermi-LAT data). We determine, under the assumption that the source-count distribution of the brightest unresolved objects is dominated by blazars, their gamma-ray luminosity function and spectral energy distribution down to fluxes almost two orders of magnitude smaller than the threshold for detecting resolved sources. The different approaches provide consistent predictions for the gamma-ray luminosity function of blazars, and they show a significant complementarity.

Complementarity of Stacking and Multiplet Constraints on the Blazar Contribution to the Cumulative High-energy Neutrino Intensity

Abstract We investigate the blazar contribution to the cumulative neutrino intensity assuming a generic relationship between neutrino and gamma-ray luminosities, L ν ∝ ( L ph ) γ lw . Using the gamma-ray luminosity functions for blazars, including flat-spectrum radio quasars (FSRQs) and BL Lac objects, as well as the Fermi-LAT detection efficiency, we estimate contributions from blazars resolved by Fermi-LAT as well as the unresolved counterpart. Combining the existing upper limits from stacking analyses, the cumulative neutrino flux from all blazars (including Fermi-LAT resolved and unresolved ones) are constrained in the range 0 ≲ γ lw ≲ 2.5. We also evaluate the effects of the redshift evolution and the effective local number densities for each class of FSRQs, BL Lacs, and all blazars, on which we place another type of constraints on the blazar contribution using the nondetection of high-energy neutrino multiplets. We demonstrate that these two upper limits are complementary, and that the joint consideration of the stacking and multiplet analyses not only supports the argument that blazars are disfavored as the dominant sources of the 100 TeV neutrino background, but also extends it by including Fermi-LAT-unresolved blazars as well, for a more generic luminosity correlation L ν ∝ ( L ph ) γ lw .

Gamma-ray luminosity function of BL Lac objects and contribution to the extragalactic gamma-ray background

ABSTRACT Using a significantly enlarged Fermi-LAT BL Lac objects (BL Lacs) sample, we construct the gamma-ray luminosity function (GLF) of BL Lacs, by the joint use of the space density distribution and source counts distribution. We use three well-studied forms of the GLF, i.e. the forms of pure density evolution (PDE), pure luminosity evolution (PLE), and luminosity-dependent density evolution (LDDE). The Markov Chain Monte Carlo (MCMC) technique is used to constrain model parameters. Our results suggest that LDDE model can give the best description for the BL Lac GLF. And the model shows that the BL Lacs with a harder GeV spectrum and a less luminosity evolve as strongly as flat spectrum radio quasars, and the evolution decreases as increasing luminosity. We also model the average photon spectra of BL Lacs with a double power-laws model. Using this modelled spectra, BL Lacs contribute $\sim \!20{{\ \rm per\ cent}}$ of the total extragalactic gamma-ray background (EGB) at E > 100 MeV, $\sim \!100{{\ \rm per\ cent}}$ of the EGB at E > 50 GeV, and the unresolved BL Lacs contribute $\sim\! 20 {{\ \rm per\ cent}}$ of the isotropic diffuse gamma-ray background at E > 100 MeV. A prediction of the TeV EGB spectra are given, which may be tested by the future detectors.

Cherenkov telescope array extragalactic survey discovery potential and the impact of axion-like particles and secondary gamma rays

The Cherenkov Telescope Array (CTA) is about to enter construction phase and one of its main key science projects is to perform an unbiased survey in search of extragalactic sources. We make use of both the latest blazar gamma–ray luminosity function and spectral energy distribution to derive the expected number of detectable sources for both the planned Northern and Southern arrays of the CTA observatory. We find that a shallow, wide survey of about 0.5 hour per field of view would lead to the highest number of blazar detections. Furthermore, we investigate the effect of axion-like particles and secondary gamma rays from propagating cosmic rays on the source count distribution, since these processes predict different spectral shape from standard extragalactic background light attenuation. We can generally expect more distant objects in the secondary gamma-ray scenario, while axion-like particles do not significantly alter the expected distribution. Yet, we find that, these results strongly depend on the assumed magnetic field strength during the propagation. We also provide source count predictions for the High Altitude Water Cherenkov observatory (HAWC), the Large High Altitude Air Shower Observatory (LHAASO) and a novel proposal of a hybrid detector.

The gamma-ray pulsar population of globular clusters: implications for the GeV excess

It has been suggested that the GeV excess, observed from the region surrounding the Galactic Center, might originate from a population of millisecond pulsars that formed in globular clusters. With this in mind, we employ the publicly available Fermi data to study the gamma-ray emission from 157 globular clusters, identifying a statistically significant signal from 25 of these sources (ten of which are not found in existing gamma-ray catalogs). We combine these observations with the predicted pulsar formation rate based on the stellar encounter rate of each globular cluster to constrain the gamma-ray luminosity function of millisecond pulsars in the Milky Way's globular cluster system. We find that this pulsar population exhibits a luminosity function that is quite similar to those millisecond pulsars observed in the field of the Milky Way (i.e. the thick disk). After pulsars are expelled from a globular cluster, however, they continue to lose rotational kinetic energy and become less luminous, causing their luminosity function to depart from the steady-state distribution. Using this luminosity function and a model for the globular cluster disruption rate, we show that millisecond pulsars born in globular clusters can account for only a few percent or less of the observed GeV excess. Among other challenges, scenarios in which the entire GeV excess is generated from such pulsars are in conflict with the observed mass of the Milky Way's Central Stellar Cluster.

STATISTICAL MEASUREMENT OF THE GAMMA-RAY SOURCE-COUNT DISTRIBUTION AS A FUNCTION OF ENERGY

ABSTRACT Statistical properties of photon count maps have recently been proven as a new tool to study the composition of the gamma-ray sky with high precision. We employ the 1-point probability distribution function of six years of Fermi-LAT data to measure the source-count distribution dN/dS and the diffuse components of the high-latitude gamma-ray sky as a function of energy. To that aim, we analyze the gamma-ray emission in five adjacent energy bands between 1 and 171 GeV. It is demonstrated that the source-count distribution as a function of flux is compatible with a broken power law up to energies of ∼50 GeV. The index below the break is between 1.95 and 2.0. For higher energies, a simple power-law fits the data, with an index of 2.2 − 0.3 + 0.7 in the energy band between 50 and 171 GeV. Upper limits on further possible breaks as well as the angular power of unresolved sources are derived. We find that point-source populations probed by this method can explain 83 − 13 + 7 % ( 81 − 19 + 52 %) of the extragalactic gamma-ray background between 1.04 and 1.99 GeV (50 and 171 GeV). The method has excellent capabilities for constraining the gamma-ray luminosity function and the spectra of unresolved blazars.

Cosmological tests with the FSRQ gamma-ray luminosity function

The extensive catalog of $\gamma$-ray selected flat-spectrum radio quasars (FSRQs) produced by \emph{Fermi} during a four-year survey has generated considerable interest in determining their $\gamma$-ray luminosity function (GLF) and its evolution with cosmic time. In this paper, we introduce the novel idea of using this extensive database to test the differential volume expansion rate predicted by two specific models, the concordance $\Lambda$CDM and $R_{\rm h}=ct$ cosmologies. For this purpose, we use two well-studied formulations of the GLF, one based on pure luminosity evolution (PLE) and the other on a luminosity-dependent density evolution (LDDE). Using a Kolmogorov-Smirnov test on one-parameter cumulative distributions (in luminosity, redshift, photon index and source count), we confirm the results of earlier works showing that these data somewhat favour LDDE over PLE; we show that this is the case for both $\Lambda$CDM and $R_{\rm h}=ct$. Regardless of which GLF one chooses, however, we also show that model selection tools very strongly favour $R_{\rm h}=ct$ over $\Lambda$CDM. We suggest that such population studies, though featuring a strong evolution in redshift, may nonetheless be used as a valuable independent check of other model comparisons based solely on geometric considerations.

High energy neutrinos from gamma-ray burst fireballs

The diffuse high-energy neutrino emission from long and short gamma-ray bursts (GRBs) is studied within the fireball emission model. By requiring that the GRB high-energy neutrino counterparts follow up-to-date gamma-ray luminosity functions and redshift evolutions, we find that GRBs could contribute up to a few percents to the observed IceCube high-energy neutrino flux for sub-PeV energies, if the latter has a diffuse origin. Our findings suggest that larger exposure is mandatory to detect neutrinos from GRBs in future stacking searches.

The gamma-ray luminosity function of millisecond pulsars and implications for the GeV excess

It has been proposed that a large population of unresolved millisecond pulsars (MSPs) could potentially account for the excess of GeV-scale gamma-rays observed from the region surrounding the Galactic Center. The viability of this scenario depends critically on the gamma-ray luminosity function of this source population, which determines how many MSPs Fermi should have already detected as resolved point sources. In this paper, we revisit the gamma-ray luminosity function of MSPs, without relying on uncertain distance measurements. Our determination, based on a comparison of models with the observed characteristics of the MSP population, suggests that Fermi should have already detected a significant number of sources associated with such a hypothesized Inner Galaxy population. We cannot rule out a scenario in which the MSPs residing near the Galactic Center are systematically less luminous than those present in the Galactic Plane or within globular clusters.

Diffuse emission of high-energy neutrinos from gamma-ray burst fireballs

Gamma-ray bursts (GRBs) have been suggested as possible sources of the high-energy neutrino flux recently detected by the IceCube telescope. We revisit the fireball emission model and elaborate an analytical prescription to estimate the high-energy neutrino prompt emission from pion and kaon decays, assuming that the leading mechanism for the neutrino production is lepto-hadronic. To this purpose, we include hadronic, radiative and adiabatic cooling effects and discuss their relevance for long- (including high- and low-luminosity) and short-duration GRBs. The expected diffuse neutrino background is derived, by requiring that the GRB high-energy neutrino counterparts follow up-to-date gamma-ray luminosity functions and redshift evolutions of the long and short GRBs. Although dedicated stacking searches have been unsuccessful up to now, we find that GRBs could contribute up to a few % to the observed IceCube high-energy neutrino flux for sub-PeV energies, assuming that the latter has a diffuse origin. Gamma-ray bursts, especially low-luminosity ones, could however be the main sources of the IceCube high-energy neutrino flux in the PeV range. While high-luminosity and low-luminosity GRBs have comparable intensities, the contribution from the short-duration component is significantly smaller. Our findings confirm the most-recent IceCube results on the GRB searches and suggest that larger exposure is mandatory to detect high-energy neutrinos from high-luminosity GRBs in the near future.

Diffuse emission of high-energy neutrinos from gamma-ray burst fireballs

Gamma-ray bursts (GRBs) have been suggested as possible sources of the high-energy neutrino flux recently detected by the IceCube telescope. We revisit the fireball emission model and elaborate an analytical prescription to estimate the high-energy neutrino prompt emission from pion and kaon decays, assuming that the leading mechanism for the neutrino production is lepto-hadronic. To this purpose, we include hadronic, radiative and adiabatic cooling effects and discuss their relevance for long- (including high- and low-luminosity) and short-duration GRBs. The expected diffuse neutrino background is derived, by requiring that the GRB high-energy neutrino counterparts follow up-to-date gamma-ray luminosity functions and redshift evolutions of the long and short GRBs. Although dedicated stacking searches have been unsuccessful up to now, we find that GRBs could contribute up to a few % to the observed IceCube high-energy neutrino flux for sub-PeV energies, assuming that the latter has a diffuse origin. Gamma-ray bursts, especially low-luminosity ones, could however be the main sources of the IceCube high-energy neutrino flux in the PeV range. While high-luminosity and low-luminosity GRBs have comparable intensities, the contribution from the short-duration component is significantly smaller. Our findings confirm the most-recent IceCube results on the GRB searches and suggest that larger exposure is mandatory to detect high-energy neutrinos from high-luminosity GRBs in the near future.

Prospect for future MeV gamma-ray active galactic nuclei population studies

Abstract While the X-ray, GeV gamma-ray, and TeV gamma-ray skies have been extensively studied, the MeV gamma-ray sky has not yet been well investigated since the Imaging Compton Telescope (COMPTEL) scanned the sky about two decades ago. In this paper, we investigate the prospects for active galactic nuclei population studies with future MeV gamma-ray missions using recent spectral models and luminosity functions of Seyferts and flat-spectrum radio quasars (FSRQs). Both of them are plausible candidates for the origins of the cosmic MeV gamma-ray background. If the cosmic MeV gamma-ray background radiation is dominated by non-thermal emission from Seyferts, a sensitivity of 10−12 erg cm−2 s−1 is required to detect several hundred Seyferts in the entire sky. If FSRQs make up the cosmic MeV gamma-ray background, a sensitivity of ∼ 4 × 10−12 erg cm−2 s−1 is required to detect several hundred FSRQs following the recent FSRQ X-ray luminosity function. However, based on the latest FSRQ gamma-ray luminosity function, with which FSRQs can explain up to ∼ 30% of the MeV background, we can expect several hundred FSRQs even with a sensitivity of 10−11 erg cm−2 s−1, which is almost the same as the sensitivity goal of the next-generation MeV telescopes.

EVIDENCE OF CROSS-CORRELATION BETWEEN THE CMB LENSING AND THE γ -RAY SKY

We report the measurement of the angular power spectrum of cross-correlation between the unresolved component of the Fermi-LAT gamma-ray sky-maps and the CMB lensing potential map reconstructed by the Planck satellite. The matter distribution in the Universe determines the bending of light coming from the last scattering surface. At the same time, the matter density drives the growth history of astrophysical objects, including their capability at generating non-thermal phenomena, which in turn give rise to gamma-ray emissions. The Planck lensing map provides information on the integrated distribution of matter, while the integrated history of gamma-ray emitters is imprinted in the Fermi-LAT sky maps. We report here the first evidence of their correlation. We find that the multipole dependence of the cross-correlation measurement is in agreement with current models of the gamma-ray luminosity function for AGN and star forming galaxies, with a statistical evidence of 3.0$\sigma$. Moreover, its amplitude can in general be matched only assuming that these extra-galactic emitters are also the bulk contribution of the measured isotopic gamma-ray background (IGRB) intensity. This leaves little room for a big contribution from galactic sources to the IGRB measured by Fermi-LAT, pointing toward a direct evidence of the extragalactic origin of the IGRB.

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.

Astrophysical ZeV acceleration in the relativistic jet from an accreting supermassive blackhole

An accreting supermassive blackhole, the central engine of active galactic nucleus (AGN), is capable of exciting extreme amplitude Alfven waves whose wavelength (wave packet) size is characterized by its clumpiness. The pondermotive force and wakefield are driven by these Alfven waves propagating in the AGN (blazar) jet, and accelerate protons/nuclei to extreme energies beyond Zetta-electron volt (ZeV=1021eV). Such acceleration is prompt, localized, and does not suffer from the multiple scattering/bending enveloped in the Fermi acceleration that causes excessive synchrotron radiation loss beyond 1019eV. The production rate of ZeV cosmic rays is found to be consistent with the observed gamma-ray luminosity function of blazars and their time variabilities.

Contribution from blazar cascade emission to the extragalactic gamma-ray background: what role does the extragalactic magnetic field play?

We estimate the contribution to the extragalactic gamma-ray background (EGRB) from both intrinsic and cascade emissions produced by blazars using a simple semi- analysis method for two models of the blazar gamma-ray luminosity function (GLF). For the cascade emission, we consider two possible contributions: one is due to that the flux of the cascade emission is lower than the Fermi LAT sensitivity (case I), which is independent on the extragalactic magnetic field (EGMF), another is due to the fact that the flux of the cascade emission is larger than the Fermi LAT sensitivity but the emission angle is larger than LAT point-spread-function (PSF) angle (case II), which depends on the EGMF. Our results indicate that (1) blazar contribution to the EGRB is dominant although it depends on the GLF model and the EGMF; (2) the EGMF plays an important role in estimating the contribution from the cascade emission produced by blazars, the contribution from the cascade emission will significantly alter the EGRB spectrum when the strength of the EGMF is large enough (say BEGMF > 10-12 G); and (3) since the cascade emission in case II reaches a saturation when the strength of the EGMF is ? 10-11 G, it is very possible that the contribution from the cascade emission produced by blazars can be considered as another method to probe the upper limit of the strength of the EGMF.

THE EFFECT OF BLAZAR SPECTRAL BREAKS ON THE BLAZAR CONTRIBUTION TO THE EXTRAGALACTIC GAMMA-RAY BACKGROUND

The spectral shapes of the contributions of different classes of unresolved gamma-ray emitters can provide insight into their relative contributions to the extragalactic gamma-ray background (EGB) and the natures of their spectra at GeV energies. We calculate the spectral shapes of the contributions to the EGB arising from BL Lacertae type objects (BL Lacs) and flat-spectrum radio quasars (FSRQs) assuming blazar spectra can be described as broken power laws. We fit the resulting total blazar spectral shape to the Fermi Large Area Telescope measurements of the EGB, finding that the best-fit shape reproduces well the shape of the Fermi EGB for various break scenarios. We conclude that a scenario in which the contribution of blazars is dominant cannot be excluded on spectral grounds alone, even if spectral breaks are shown to be common among Fermi blazars. We also find that while the observation of a featureless (within uncertainties) power-law EGB spectrum by Fermi does not necessarily imply a single class of contributing unresolved sources with featureless individual spectra, such an observation and the collective spectra of the separate contributing populations determine the ratios of their contributions. As such, a comparison with studies including blazar gamma-ray luminosity functions could have profound implications for the blazar contribution to the EGB, blazar evolution, and blazar gamma-ray spectra and emission.

CONTRIBUTION OF GAMMA-RAY-LOUD RADIO GALAXIES’ CORE EMISSIONS TO THE COSMIC MeV AND GeV GAMMA-RAY BACKGROUND RADIATION

The Fermi gamma-ray satellite has recently detected gamma-ray emissions from radio galaxy cores. From these samples, we first examine the correlation between the luminosities at 5 GHz, L_{5GHz}, and at 0.1-10 GeV, L_{gamma}, of these gamma-ray loud radio galaxies. We find that the correlation is significant with L_{gamma} \propto L_{5GHz}^{1.16} based on a partial correlation analysis. Using this correlation and the radio luminosity function (RLF) of radio galaxies, we further explore the contribution of gamma-ray loud radio galaxies to the unresolved extragalactic gamma-ray background (EGRB). The gamma-ray luminosity function is obtained by normalizing the RLF to reproduce the source count distribution of the Fermi gamma-ray loud radio galaxies. We find that gamma-ray loud radio galaxies will explain ~25% of the unresolved Fermi EGRB flux above 100 MeV and will also make a significant contribution to the EGRB in the 1-30 MeV energy band. Since blazars explain 22% of the EGRB above 100 MeV, radio loud active galactic nuclei (AGNs) population explains ~47% of the unresolved EGRB. We further make an interpretation on the origin of the EGRB. The observed EGRB spectrum at 0.2-100 GeV does not show an absorption signature by the extragalactic background light. Thus, the dominant population of the origin of EGRB at very high energy (>30 GeV) might be nearby gamma-ray emitting sources or sources with very hard gamma-ray spectrum.

ERRATUM: “THE BLAZAR SEQUENCE AND THE COSMIC GAMMA-RAY BACKGROUND RADIATION IN THEFERMIERA” (2009, ApJ, 702, 523)

(Abridged) We present a new model of the blazar gamma-ray luminosity function (GLF) and the spectrum of the extragalactic gamma-ray background (EGRB), which is consistent with the observed distributions of EGRET blazars. The unified sequence of blazar spectral energy distribution (SED) is taken into account to make a non-trivial prediction for the EGRB spectrum and more realistic comparison with the data than previous studies. We then try to explain the EGRB data by the two AGN populations: one is blazars, and the other is non-blazar AGNs that are responsible for the EGRB in the MeV band. We find that ~80% of the EGRB photon flux at > 100 MeV can be explained by the sum of the two populations, while ~45% can be accounted for only by blazars. The predicted EGRB spectrum is in agreement with a wide range of the observed data from X-ray to GeV, within the systematic uncertainties in the EGRB determination by EGRET. These results indicate that AGNs including blazars are the primary source of EGRB. Blazars are dominant in EGRB at higher energy bands of >100 MeV, while non-blazar AGNs dominate at <100 MeV. Almost all of the EGRB flux from blazars will be resolved into discrete sources by the Fermi Gamma-ray Space Telescope, while that from non-blazar AGNs will largely remain unresolved. Therefore, comparison between the integrated source counts and diffuse EGRB flux as a function of photon energy will give a simple and clear test of our model. Various quantitative predictions for Fermi observations are also made. Especially, our model predicts 600-1200 blazars in all sky down to 2e-9 photons/cm^-2/s (>100 MeV), which is considerably smaller than most of previous studies.