Pith. sign in

REVIEW 2 major objections 4 minor 1 cited by

{\it Fermi}-LAT Stacking Analysis Technique: An Application to Extreme Blazars and Prospects for their CTA Detection

T0 review · 2 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Stacking the likelihood profiles of 172 Fermi-LAT-undetected extreme blazars produces a cumulative gamma-ray detection at $32\sigma$, with an average 10--1000 GeV flux of $6.51\times10^{-12}$ ph cm$^{-2}$ s$^{-1}$ and photon index $2.08$.

desk verdict A well-validated stacking analysis that delivers a robust cumulative detection of 172 undetected extreme blazars; the EGB lower limit carries a selection-bias caveat, but the central result holds. read the letter →

arxiv 1908.02496 v1 pith:ENAJHIUT submitted 2019-08-07 astro-ph.HE

Vaidehi S. Paliya , Alberto Dominguez , Marco Ajello , Anna Franckowiak , Dieter Hartmann This is my paper · ORCID
classification astro-ph.HE
keywords gamma-raystackingextremeblazarsFermi-LATlikelihoodprofileextragalacticbackgroundCherenkovTelescopeArrayBLLacertaeobjectsEBLabsorption
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper claims that a population of 172 extreme blazars not individually detected by Fermi-LAT nevertheless emits gamma rays as a class, and that the emission can be pulled out by stacking their individual likelihood profiles. The combined signal is detected at about $32\sigma$, which the authors take as proof that the method works and that these unresolved objects are real GeV emitters. The stacked average spectrum is hard, and after correcting for extragalactic background light absorption it shows no softening up to 1 TeV, implying the inverse-Compton peak lies above that energy. Comparing the stacked flux with CTA sensitivity, the paper concludes that the unresolved population as a whole will probably stay below CTA's threshold, while Fermi-LAT-detected extreme blazars are the better targets. Finally, the paper estimates that at 100 GeV at least about 10% of the extragalactic gamma-ray background comes from these unresolved extreme blazars.

What carries the argument

The method is likelihood-profile stacking. For each undetected source, the pipeline builds a two-dimensional grid of gamma-ray photon flux and photon index, fits the Fermi-LAT data at every grid point with all other model parameters fixed except the diffuse backgrounds, and subtracts the log-likelihood at the lowest flux to obtain a test-statistic profile. Because log-likelihoods are additive, the individual TS profiles are summed across the 172 sources and the location of the combined peak, found by spline fitting, gives the average flux and spectral index of the population. The paper validates the machinery against empty-sky positions and injected simulated sources.

What would settle it

Run the same likelihood-profile stack on the full extreme-blazar catalog, including sources with photometric redshift limits or no redshift measurement. If the recovered average flux falls significantly below $6.51\times10^{-12}$ ph cm$^{-2}$ s$^{-1}$ or the 100 GeV EGB share drops well below 10%, the known-redshift subsample is not representative.

Watch

Extended reading notes

Core claim

The central discovery is a cumulative gamma-ray signal from 172 Fermi-LAT-undetected extreme blazars: stacking their likelihood profiles yields TS = 1062, corresponding to more than $32\sigma$ for 2 degrees of freedom. The best-fit average spectral parameters are $F_{10-1000\,{\rm GeV}} = 6.51^{+0.36}_{-0.35}\times10^{-12}$ ph cm$^{-2}$ s$^{-1}$ and $\Gamma = 2.08^{+0.07}_{-0.06}$. This signal is validated by stacking empty sky positions, which reproduce the null hypothesis, and by simulations of 100 injected sources, whose input flux and index are recovered. The authors use the average spectrum to argue that the unresolved extreme-blazar population contributes at least $\sim$10% of the extragalactic gamma-ray background at 100 GeV, and that the EBL-corrected spectrum remains hard to 1 TeV, placing the inverse-Compton peak above that energy.

Load-bearing premise

The 172 unresolved sources with known redshifts are representative of all extreme blazars; if the redshift requirement systematically excludes faint or distant members of the class, the measured average flux and the derived 10% contribution to the gamma-ray background would be biased upward.

Editorial extensions

If this is right

  • The unresolved population of extreme blazars is established as a genuine collective GeV source, so it must be counted in models of the gamma-ray sky and of the extragalactic gamma-ray background.
  • Taking the stacked spectrum at face value, CTA's 50-hour sensitivity is unlikely to detect the average unresolved extreme blazar; any CTA strategy should prioritise the 165 Fermi-LAT-detected extreme blazars.
  • The absence of EBL-corrected softening up to 1 TeV locates the inverse-Compton peak of extreme blazars above 1 TeV, implying intrinsically hard TeV spectra.
  • At 100 GeV, unresolved extreme blazars make up at least about 10% of the extragalactic gamma-ray background, so they cannot be ignored in background accounting.
  • The stacking pipeline can extract signals roughly an order of magnitude fainter than the Fermi-LAT point-source threshold, making it suitable for other unresolved astrophysical populations.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the known-redshift requirement biases the 172-source sample toward brighter or lower-redshift objects, the true population average could be fainter than $6.5\times10^{-12}$ ph cm$^{-2}$ s$^{-1}$, and the 10% EGB share could be an overestimate rather than a floor; the paper does not correct for this selection.
  • The same profile-stacking approach could be applied to other position-catalogued but gamma-ray-undetected classes, such as star-forming galaxies or galaxy clusters; the simulation test the paper performs for extreme blazars would need repeating for each class.
  • A testable extension is to stack the known extreme blazars in two redshift slices; if the average spectrum is redshift-independent after EBL correction, the claim of a hard intrinsic TeV spectrum is strengthened.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The paper presents a stacking technique for Fermi-LAT data that combines individual source likelihood profiles to measure the average gamma-ray properties of a population of unresolved sources. The method is applied to 172 gamma-ray-undetected extreme blazars selected from the 3HSP catalog with known redshift, yielding a cumulative detection with TS = 1062 (quoted as ~32σ for 2 dof), an average 10–1000 GeV flux of 6.51e-12 ph/cm2/s, and a photon index of 2.08. The authors validate the method with 172 empty-sky positions that show null stacked TS and with simulated sources whose input flux and index are recovered. They further compare the stacked spectra with CTA sensitivity limits, concluding that the unresolved population lies below CTA's detection threshold whereas the Fermi-LAT-detected extreme blazars are promising CTA targets. Finally, they use the average spectrum to estimate that unresolved extreme blazars contribute at least ~10% of the extragalactic gamma-ray background at 100 GeV.

Significance. If correct, the paper provides a computationally efficient and flexible stacking tool that extends Fermi-LAT reach to populations roughly an order of magnitude fainter than individual-source detection thresholds. The validation strategy is a strength: 172 empty-sky positions yield a stacked TS consistent with the null hypothesis, and simulated sources are recovered with input parameters, giving confidence that the method is not stacking diffuse background fluctuations. The result that extreme blazars form a genuinely unresolved GeV-emitting population and potentially contribute ~10% of the EGB at 100 GeV is astrophysically interesting. The work is also timely for CTA planning, as it identifies Fermi-LAT-detected extreme blazars as better targets than the currently unresolved population.

major comments (2)
  1. [Section 3, sample selection (also Section 5/footnote 7)] The stacked sample is built from 337 extreme blazars with known redshift drawn from the 2011-source 3HSP catalog. Redshift determination requires optical spectroscopy, so the known-redshift subset is plausibly biased toward optically brighter, lower-redshift, and possibly gamma-ray-brighter sources. The paper does not quantify the redshift completeness of 3HSP, nor does it compare the redshift/flux/synchrotron-peak distributions of the known-redshift subset with those of the full 3HSP extreme-blazar population. This is directly relevant to the headline claim that unresolved extreme blazars contribute 'at least ~10%' of the EGB at 100 GeV, because that estimate uses the average flux of the 172 stacked sources as if it represented the whole population. If the known-redshift subset is systematically brighter, the EGB lower limit is inflated. I request that the authors either (a) test representativeness by comparing the known-redshift sample with the full 3HSP sample in observable properties (e.g., radio flux, X-ray flux, synchrotron peak frequency, redshift distribution) and propagate any offset, or (b) explicitly restrict the EGB claim to the known-redshift subset and soften the population-level wording in the abstract and Section 5.
  2. [Section 2.2 and Section 4.2] The stacking pipeline assumes that the population is described by a single power law with a common flux and index, and the simulation validation in Section 4.2 draws simulated sources from exactly that model (Gaussian index around 2.1, log-normal flux). This confirms the method recovers input parameters under the assumed model, but it does not test how the technique behaves when the true population has a dispersion in spectral indices or spectral curvature (e.g., log-parabolic or EBL-absorbed spectra). Since the reported average index 2.08 is used to infer the position of the inverse-Compton peak and to compute the EGB contribution, it would be valuable to run a simulation with a realistic range of indices and with EBL attenuation to show that the stacked likelihood-profile method still recovers the flux-averaged index without significant bias. Such a test is within the scope of the paper and would strengthen the applicability of the method beyond the single-template case.
minor comments (4)
  1. [Abstract and Section 3] The conversion from TS = 1062 to 'more than 32σ' is only valid in the Gaussian limit with 2 dof approximately as sqrt(TS); the authors should state this explicitly or quote the chi-square tail probability instead of implying exact 32σ Gaussian significance.
  2. [Figure 1, right panel] The color-bar label 'TS (Extreme blazars)' is redundant and slightly confusing because the panel already shows TS contours; consider labeling it as 'TS' or removing it.
  3. [Section 5 and Figure 2, right panel] The text says the EGB contribution is derived from the average spectrum but does not specify whether the EBL-corrected or observed spectrum is used in that calculation; since the EGB is an observed quantity, the observed (intrinsic minus EBL absorption) spectrum should be used, and this should be clarified.
  4. [Section 2.1] The threshold of TS > 25 for adding unmodeled sources is mentioned, but the paper does not state the resulting number of such excesses in the preprocessing of the 172 undetected sources; a brief note on how many tentative sources were added and then excluded would improve reproducibility.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the stacked flux/index are measured from Fermi-LAT likelihood profiles and validated against empty-sky and simulated-source checks; EGB and CTA statements are derived inferences, not fitted inputs.

full rationale

The paper's central quantities—the stacked 10–1000 GeV flux of 6.51e-12 ph cm−2 s−1 and photon index 2.08—are obtained directly from a combined Fermi-LAT likelihood profile over the 172 undetected extreme blazars (Section 2.2, Section 3). There is no parameter fitted to an external target; the EGB fraction and CTA detectability are computed afterward from the measured average spectrum using independent inputs (Ackermann et al. 2015 for the total EGB, CTA sensitivity curves, and the Dominguez et al. 2011 EBL model). The pipeline is not defined in terms of any of these derived claims. The assumed power-law spectral shape is a modeling choice, not a circular reduction, and the paper independently validates the method: stacking 172 empty sky positions gives a TS distribution compatible with the null hypothesis, and stacking 100 simulated sources recovers input flux 1.04e-11 vs 1.0e-11 ph cm−2 s−1 and index 2.11 vs 2.1. The known-redshift sample selection is a possible bias in extrapolating to the full population, but it does not make the derivation circular. Self-citations (e.g., Dominguez & Ajello 2015; Paliya et al. 2019) are contextual and not load-bearing for the stacking result, and no uniqueness theorem or ansatz is smuggled in through author-only citations. The paper is therefore self-contained against its own data and external benchmarks.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The paper's central claims rest on standard gamma-ray sky modeling and two explicit simplifications: a single power-law average spectrum and a uniform sky distribution for the EGB estimate. No new physical entities are introduced. The fitted flux and index are the measured output of the stacked likelihood, not ad hoc inputs.

free parameters (2)
  • Average photon flux (10-1000 GeV) = 6.51e-12 ph cm^-2 s^-1
    Best-fit peak of the stacked TS profile for 172 undetected extreme blazars; this is the measured average flux, not an ad hoc input.
  • Average photon index (10-1000 GeV) = 2.08
    Best-fit photon index from the same stacked likelihood profile; this is the measured average spectral slope.
assumptions (5)
  • domain assumption The 4FGL catalog, interstellar emission model, and isotropic template accurately describe the gamma-ray sky in each ROI.
    Used throughout Section 2.1 to model sources and diffuse emission; residual mismodeling could bias the stacked signal, though the background-stacking test with empty positions yields null results.
  • domain assumption The average spectral shape of the source population is a single power law over 10-1000 GeV.
    Stated in Section 2.2 step 1; if the true average spectrum is curved, the fitted flux and index would be biased, affecting derived quantities.
  • domain assumption The EBL model of Dominguez et al. (2011) correctly describes gamma-ray absorption.
    Used in Section 5 to derive EBL-corrected spectra; an incorrect EBL would alter the conclusion that the inverse Compton peak lies above 1 TeV.
  • domain assumption Undetected extreme blazars are uniformly distributed on the sky outside the Galactic plane.
    Stated in Section 5; the authors note this yields a lower limit and a proper treatment requires a luminosity function.
  • domain assumption The 172 undetected objects with known redshift are representative of all gamma-ray undetected extreme blazars.
    Sample selection in Section 3; the known-redshift requirement may bias against faint, high-redshift sources, limiting the representativeness of the average spectrum and EGB contribution.

how reviews work

0 comments
Cite this review

Pith. "Pith review of {\it Fermi}-LAT Stacking Analysis Technique: An Application to Extreme Blazars and Prospects for their CTA Detection." pith.science (2026). https://pith.science/paper/ENAJHIUT

@misc{pith2026190802496,
  author       = {Pith},
  title        = {Pith review of: \it Fermi-LAT Stacking Analysis Technique: An Application to Extreme Blazars and Prospects for their CTA Detection},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ENAJHIUT}},
  note         = {Machine review of arXiv:1908.02496}
}
abstract

We present a likelihood profile stacking technique based on the {\it Fermi}-Large Area Telescope (LAT) data to explore the $\gamma$-ray characteristics of {\it Fermi}-LAT undetected astrophysical populations. The pipeline is applied to a sample of $\gamma$-ray unresolved extreme blazars, i.e., sources with the highest synchrotron peak frequencies ($\nu_{\rm Syn}^{\rm peak}\geqslant 10^{17}$ Hz), and we report a cumulative $\gamma$-ray detection with more than 32$\sigma$ confidence for 2 degrees of freedom. Comparing the generated stacked $\gamma$-ray spectrum with the sensitivity limits of the upcoming Cherenkov Telescope Array (CTA), we find that the {\it Fermi}-LAT undetected population of such extreme blazars, on average, may remain well below the CTA detection threshold due to their faintness and extragalactic background light (EBL) absorption. However, $\gamma$-ray detected blazars belonging to the same class are promising candidates for CTA observations. The EBL corrected stacked spectra of these sources do not show any softening up to 1 TeV. This finding suggests the inverse Compton peak of extreme blazars to lie above 1 TeV, thus indicating a hard intrinsic TeV spectrum. Our analysis also predicts that at 100 GeV, at least $\sim$10\% of the diffuse extragalactic $\gamma$-ray background originates from the $\gamma$-ray undetected extreme blazars. These results highlight the effectiveness of the developed stacking technique to explore the uncharted territory of $\gamma$-ray undetected astrophysical objects.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. An 18-year Fermi-LAT stacking limit on GeV $\gamma$-ray emission from particle-accelerating colliding-wind binaries

    astro-ph.HE 2026-07 accept novelty 6.0 of 10

    Clean stack of six high-latitude PACWBs is null (p=0.83), yielding F(>1 GeV)≲1.1e-11 ph cm^{-2} s^{-1} and η≲4e-6 (d/kpc)^2, two orders below η Car.

Reference graph

Works this paper leans on

30 extracted references · 21 canonical work pages · cited by 1 Pith paper

  1. [1]

    A., Ackermann, M., Agudo, I., et al

    Abdo, A. A., Ackermann, M., Agudo, I., et al. 2010, ApJ, 716, 30

    2010

  2. [2]

    2018, Science, 362, 1031

    Abdollahi, S., Ackermann, M., Ajello, M., et al. 2018, Science, 362, 1031

    work page 2018

  3. [3]

    2016, ApJS, 223, 26

    Acero, F., Ackermann, M., Ajello, M., et al. 2016, ApJS, 223, 26

    2016

  4. [4]

    2011, Physical Review Letters, 107, 241302

    Ackermann, M., Ajello, M., Albert, A., et al. 2011, Physical Review Letters, 107, 241302

    work page 2011

  5. [5]

    2012, ApJ, 755, 164

    Ackermann, M., Ajello, M., Allafort, A., et al. 2012, ApJ, 755, 164

    2012

  6. [6]

    2014, PhRvD, 89, 042001

    Ackermann, M., Albert, A., Anderson, B., et al. 2014, PhRvD, 89, 042001

    work page 2014

  7. [7]

    2015, ApJ, 799, 86 —

    Ackermann, M., Ajello, M., Albert, A., et al. 2015, ApJ, 799, 86 —. 2016, PhRvL, 116, 151105

    work page 2015

  8. [8]

    2011, Experimental Astronomy, 32, 193

    Actis, M., Agnetta, G., Aharonian, F., et al. 2011, Experimental Astronomy, 32, 193

    2011

Show all 30 references
  1. [9]

    S., Romani, R

    Ajello, M., Shaw, M. S., Romani, R. W., et al. 2012, ApJ, 751, 108

    2012

  2. [10]

    W., Gasparrini, D., et al

    Ajello, M., Romani, R. W., Gasparrini, D., et al. 2014, ApJ, 780, 73

    2014

  3. [11]

    2015, ApJL, 800, L27

    Ajello, M., Gasparrini, D., S´anchez-Conde, M., et al. 2015, ApJL, 800, L27

    2015

  4. [12]

    2018, MNRAS, 480, 2165

    Arsioli, B., Barres de Almeida, U., Prandini, E., Fraga, B., & Foffano, L. 2018, MNRAS, 480, 2165

    2018

  5. [13]

    2013, ArXiv e-prints, arXiv:1303.3514 B¨ottcher, M

    Atwood, W., Albert, A., Baldini, L., et al. 2013, ArXiv e-prints, arXiv:1303.3514 B¨ottcher, M. 2019, Galaxies, 7, 20

    2013 arXiv

  6. [14]

    L., Arsioli, B., Giommi, P., & Padovani, P

    Chang, Y . L., Arsioli, B., Giommi, P., & Padovani, P. 2017, A&A, 598, A17

    2017

  7. [15]

    Chang, Y .-L., Arsioli, B., Giommi, P., Padovani, P., & Brandt, C. H. 2019, A&A, submitted

    2019

  8. [16]

    N., Torres, D

    Cillis, A. N., Torres, D. F., & Reimer, O. 2005, ApJ, 621, 139

    2005

  9. [17]

    2001, A&A, 371, 512 Dom´ınguez, A., & Ajello, M

    Costamante, L., , G., Giommi, P., et al. 2001, A&A, 371, 512 Dom´ınguez, A., & Ajello, M. 2015, ApJ, 813, L34 Dom´ınguez, A., Primack, J. R., Rosario, D. J., et al. 2011, MNRAS, 410, 2556

    2001

  10. [18]

    2013, A&A Rv, 21, 64

    Dubus, G. 2013, A&A Rv, 21, 64

    2013

  11. [19]

    2019, arXiv e-prints, arXiv:1907.03666

    Franceschini, A., Foffano, L., Prandini, E., & Tavecchio, F. 2019, arXiv e-prints, arXiv:1907.03666

    2019 arXiv

  12. [20]

    2012, A&A, 547, A102

    Huber, B., Farnier, C., Manalaysay, A., Straumann, U., & Walter, R. 2012, A&A, 547, A102

    2012

  13. [21]

    R., Bertsch, D

    Mattox, J. R., Bertsch, D. L., Chiang, J., et al. 1996, ApJ, 461, 396

    1996

  14. [22]

    2017, International Cosmic Ray Conference, 301, 798

    Moiseev, A., & Amego Team. 2017, International Cosmic Ray Conference, 301, 798

    2017

  15. [23]

    C., & Celotti, A

    Padovani, P., Ghisellini, G., Fabian, A. C., & Celotti, A. 1993, MNRAS, 260, L21

    1993

  16. [24]

    S., Marcotulli, L., Ajello, M., et al

    Paliya, V . S., Marcotulli, L., Ajello, M., et al. 2017, ApJ, 851, 33

    2017

  17. [25]

    S., Ajello, M., Marcotulli, L., et al

    Paliya, V . S., Ajello, M., Marcotulli, L., et al. 2019, arXiv e-prints, arXiv:1903.06106

    2019 arXiv

  18. [26]

    Petropoulou, M., & Dermer, C. D. 2016, ApJL, 825, L11

    2016

  19. [27]

    Reimer, O., Pohl, M., Sreekumar, P., & Mattox, J. R. 2003, ApJ, 588, 155

    2003

  20. [28]

    M., Fried, J

    Stickel, M., Padovani, P., Urry, C. M., Fried, J. W., & Kuehr, H. 1991, ApJ, 374, 431

    1991

  21. [29]

    2011, MNRAS, 414, 3566 The Fermi-LAT collaboration

    Tavecchio, F., Ghisellini, G., Bonnoli, G., & Foschini, L. 2011, MNRAS, 414, 3566 The Fermi-LAT collaboration. 2019a, arXiv e-prints, arXiv:1902.10045 —. 2019b, arXiv e-prints, arXiv:1905.10771

    2011 arXiv

  22. [30]

    2017, ArXiv e-prints, arXiv:1707.09551

    Wood, M., Caputo, R., Charles, E., et al. 2017, ArXiv e-prints, arXiv:1707.09551

    2017 arXiv

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.