REVIEW 4 major objections 4 minor 22 references
Observations of blazar PKS 2023-07 in flaring state with HESS and Fermi-LAT in 2016-2017 and constraints on an intrinsic cut-off
T0 review · 4 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read During gamma-ray flares of the distant quasar PKS 2023-07, H.E.S.S. non-detection forces an intrinsic cutoff below 35 GeV or absorption within a few hundred gravitational radii.
desk verdict First VHE upper limits for the highest-redshift FSRQ, with plausible cutoff constraints that carry an unquantified systematic from the fixed -3 spectral index used to derive the H.E.S.S. limits. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing comparison is between the Fermi-LAT power-law spectrum extrapolated to very high energies (attenuated by the nominal EBL model) and the H.E.S.S. 95% confidence differential upper limits, which are derived assuming a photon index of $-3$. Each H.E.S.S. upper limit is converted into a Gaussian likelihood, and a profile-likelihood-ratio test is scanned over three physical parameters: the intrinsic cutoff energy $E_\mathrm{cut}$, the EBL normalisation factor $n_\mathrm{EBL}$, and the distance $r$ of the emission region from the central black hole using a broad-line-region photon model. This scan is what turns non-detections into quantitative constraints on the source spectrum and emission geometry.
What would settle it
Re-analyse the April 2016 and September 2017 H.E.S.S. data assuming a steeper photon index (for example $-4$) and re-run the profile likelihood; if the standard EBL-attenuated Fermi-LAT power law then falls below the revised 95% upper limits, the claimed incompatibility and the $E_\mathrm{cut}$ bounds are not robust. Alternatively, a future detection of very-high-energy photons from PKS 2023-07 during a comparable flare at a flux at or above the extrapolated power-law level would disprove the need for an intrinsic cutoff or broad-line-region absorption in those states.
Extended reading notes
Core claim
The paper's central claim is that the absence of a H.E.S.S. detection during the April 2016 and September 2017 flares is an informative measurement. A power-law extrapolation of the contemporaneous Fermi-LAT spectrum, after applying a standard EBL attenuation model, overshoots the 95% confidence H.E.S.S. differential upper limits; a profile-likelihood test shows the mismatch is significant. The authors therefore infer either an intrinsic exponential cutoff in the emitted spectrum, with $E_\mathrm{cut}<128\ \mathrm{GeV}$ (April 2016) and $E_\mathrm{cut}<35\ \mathrm{GeV}$ (September 2017), or gamma-gamma absorption by broad-line-region photons, which requires the emission region to be closer than $r=1.8\times10^{17}\ \mathrm{cm}\approx120\,r_g$ (April 2016) and $r=9.5\times10^{17}\ \mathrm{cm}\approx650\,r_g$ (September 2017). Raising the EBL density instead would require factors above 1.39 and 3.07, which the authors reject as incompatible with independent EBL constraints. The October 2017 flare is consistent with an unbroken power law. Under the intrinsic-cutoff interpretation, the bound $E_\mathrm{cut}<35\ \mathrm{GeV}$ translates, for a Doppler factor of 20, to parent electron Lorentz factors $\gamma'\lesssim 3.2\times10^3$ (10 eV BLR seed photons), $\lesssim1.4\times10^3$ (50 eV disk photons), and $\lesssim3.2\times10^4$ (1000 K dust-torus photons).
Load-bearing premise
The H.E.S.S. upper limits are computed assuming the source spectrum is a power law with photon index $-3$; if the true very-high-energy spectrum is steeper, the quoted limits and the inferred cutoff bounds could change.
Editorial extensions
If this is right
- During the two brightest flares, the very-high-energy spectrum of PKS 2023-07 cannot be a simple extension of the Fermi-LAT power law with standard EBL absorption; additional suppression is required.
- If the suppression is intrinsic, it corresponds to an exponential cutoff below about 35 GeV in September 2017, which in external-Compton models bounds the emitting electron Lorentz factor to roughly $10^3$--$10^4$.
- If the suppression is environmental, the gamma-ray emission region must lie within a few hundred gravitational radii of the black hole, placing the flare inside the broad-line region.
- The October 2017 data remain compatible with an unbroken power law, showing the suppression is not present in every flaring state.
- EBL densities high enough to explain the non-detection are excluded, so the simple power-law model cannot be saved by appealing to a stronger extragalactic background light.
Reading between the lines
- If the broad-line-region absorption explanation is correct, the same physical picture predicts that the high-energy gamma-ray flare should be accompanied by reprocessed or absorbed signatures at other wavelengths; a joint optical-UV to gamma-ray variability campaign could test whether the emitting region really sits inside the BLR.
- The two flare epochs give different cutoff bounds (128 and 35 GeV); if a single particle population is responsible, the tighter September bound may indicate that the cutoff itself varies with flare state, which could be checked by correlating $E_\mathrm{cut}$ with simultaneous synchrotron or seed-photon fluxes.
- Applying the same upper-limit-plus-profile-likelihood technique to other high-redshift flat-spectrum radio quasars could map intrinsic cutoffs as a function of redshift and separate source-intrinsic acceleration limits from EBL absorption.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper reports H.E.S.S. and Fermi-LAT observations of the flat-spectrum radio quasar PKS 2023−07 (z = 1.388) during three flaring episodes in April 2016 and September/October 2017. H.E.S.S. did not detect the source in any of the campaigns, and 95% confidence-level upper limits were derived assuming a power-law spectrum with photon index −3. The authors compare these limits with near-simultaneous Fermi-LAT power-law spectra extrapolated to very high energies with nominal EBL absorption, and find that a simple power-law extrapolation is incompatible with the upper limits for the April 2016 and September 2017 flares. Using a profile likelihood approach, they derive constraints on an intrinsic exponential cutoff (E_cut < 128 GeV for April 2016 and E_cut < 35 GeV for September 2017), on the EBL normalization (needing an increase by factors 3.07 and 1.39), and on the distance of the gamma-ray emission region from the black hole if absorption occurs in the broad line region (r < 1.8×10^17 cm and r < 9.5×10^17 cm). They conclude that the non-detection is best explained by an intrinsic cutoff or BLR absorption rather than an EBL correction.
Significance. If the constraints are correct, this is the highest-redshift FSRQ for which VHE upper limits have been used to constrain the intrinsic spectrum and the location of the emission region, providing a valuable albeit indirect probe of jet physics and the EBL at z~1.4. The use of two independent H.E.S.S. analysis chains (Model and ImPACT) strengthens the reliability of the non-detection, and the multi-wavelength context is well documented. The results are, however, upper-limit-driven and rely on a series of stated modeling assumptions; the paper does not claim a detection. The main value is in demonstrating the feasibility of such constraints and in quantifying the tension between LAT extrapolations and H.E.S.S. non-detections for a distant FSRQ.
major comments (4)
- [Sec. 2.1 and Sec. 4] The H.E.S.S. differential upper limits in Fig. 2 and Table 1 are explicitly derived assuming a power-law spectrum with photon index −3, yet in Section 4 they are used to constrain models with exponential cutoffs and to compare with Fermi-LAT power laws of different spectral indices. The mapping from IACT count upper limits to differential flux upper limits depends on the assumed spectral shape through the energy threshold, effective area, and bin-to-bin migration, with the lowest-energy bin (~0.11 TeV) being most sensitive. Since the derived E_cut and n_EBL constraints rely on these upper limits, the analysis is not self-consistent: the H.E.S.S. data are not re-analyzed for each trial spectral shape. Please either re-derive the upper limits for the tested spectral models (forward-folding) or quantify the expected shift in the resulting E_cut and n_EBL constraints from the index dependence.
- [Sec. 4] The conversion of H.E.S.S. upper limits into a Gaussian likelihood is not specified. The text states that the upper limits 'can be translated into a Gaussian likelihood profile' but does not state how the 95% C.L. upper limit defines the Gaussian width, whether the likelihood is truncated at zero flux, or how the profile likelihood ratio test statistic is calibrated. Without this detail, the reported confidence intervals on E_cut, n_EBL, and r are not reproducible. Please provide the full likelihood construction, including the relation between the one-sided 95% upper limit and the Gaussian sigma, and the treatment of bins where the predicted flux is below zero.
- [Table 2] The Fermi-LAT spectral photon indices for the analyzed periods are not reported in Table 2, which only lists the model flux integrated from 0.1 to 500 GeV. The extrapolation to VHE, and thus all quantitative constraints in Section 4, depends critically on the photon index of the power-law fit. Please include the fitted index and its uncertainty for each period used in the analysis. In addition, the September 2017 intervals labelled 'sept17-a' and 'sept17-c' have no model flux or significance values in the table; clarify whether these intervals are included in the spectral analysis or only in the light curve.
- [Sec. 4, BLR absorption paragraph] The constraints on the emission-region distance r from BLR absorption are derived assuming fixed values of L_disk = 10^46 erg/s, L_Hβ = 10^43 erg/s, and M_BH = 10^9 M_sun, for which no measurements for PKS 2023−07 are cited. These quantities are degenerate with r, and the ranges of plausible values for this source are likely broad. The derived radii (1.8×10^17 cm and 9.5×10^17 cm) should be accompanied by a systematic exploration over the plausible ranges of these parameters, or at least a statement of how r scales with L_disk, L_Hβ, and M_BH, so that the reader can gauge the robustness of the BLR-absorption hypothesis.
minor comments (4)
- [Introduction] The abbreviation 'FRSQ' should be 'FSRQ' (flat-spectrum radio quasar).
- [Sec. 2.1] There is a duplicated word in the sentence 'H.E.S.S. can also operate in monoscopic mode using using only events detected with the CT5 telescope.'
- [Sec. 4] The text refers to the 'Thompson regime'; the correct spelling is 'Thomson regime'.
- [Sec. 4] The reference to the profile likelihood ratio test [18] is a paper on confidence intervals with nuisance parameters; please clarify the test statistic and the distribution used to derive the reported one-sided 95% limits.
Circularity Check
No significant circularity: the E_cut and BLR constraints are upper limits derived from the same observations rather than independent predictions, and the only self-citation is minor and not load-bearing.
full rationale
The derivation chain is self-contained. The Fermi-LAT power-law spectra (Section 2.2) and the H.E.S.S. differential upper limits (Section 2.1) are independent measurements from different instruments; the profile-likelihood analysis in Section 4 tests trial spectra against those limits and derives one-sided 95% upper limits on E_cut, the EBL normalization factor, and the BLR emission-region distance. These are constraints on parameters from the data, not predictions of independent quantities, so the derivation does not reduce to its inputs by construction. The assumed photon index of -3 used to compute the H.E.S.S. upper limits is explicitly stated; using those fixed-shape upper limits to evaluate cutoff trial spectra without full forward-folding is a methodological caveat that could shift the quoted E_cut values, but it is not a circular definitional reduction. The only self-citation, Ref. [22] (Meyer et al. 2019, with overlapping authorship), is invoked to adopt the BLR model of the independent Ref. [21] (Finke 2016); the BLR constraints are driven by the H.E.S.S./LAT likelihood, so this citation is not load-bearing. No circularity score above 2 is warranted.
Assumptions & free parameters
free parameters (7)
- E_cut (intrinsic cutoff energy) =
April 2016: < 128 GeV; September 2017: < 35 GeV
- n_EBL (EBL normalization) =
Minimum 1.39 (April 2016), 3.07 (September 2017)
- r (gamma-ray emission region distance) =
April 2016: < 1.8e17 cm (~120 r_g); September 2017: < 9.5e17 cm (~650 r_g)
- L_disk (accretion disk luminosity) =
1e46 erg/s (assumed)
- L_Hbeta (H-beta line luminosity) =
1e43 erg/s (assumed)
- M_bh (black hole mass) =
1e9 solar masses (assumed)
- Doppler factor delta_D =
20 (assumed)
assumptions (6)
- domain assumption The Dominguez et al. (2011) EBL model provides the correct energy-dependent EBL density, and deviations can be described by a single multiplicative normalization n_EBL.
- domain assumption The intrinsic gamma-ray spectrum of PKS 2023-07 is a power law with an exponential cutoff, with no other spectral features.
- domain assumption H.E.S.S. upper limits can be translated into Gaussian likelihood profiles.
- domain assumption The Fermi-LAT power-law spectrum can be extrapolated beyond the LAT energy range up to the H.E.S.S. band with no additional spectral features.
- domain assumption The H.E.S.S. upper limits are evaluated under the assumption of a power-law spectrum with photon index -3.
- domain assumption The broad line region model of Finke (2016) with ring geometry, and the scaling relation of total BLR luminosity = 0.03 L_disk, apply to PKS 2023-07.
Cite this review
Pith. "Pith review of Observations of blazar PKS 2023-07 in flaring state with HESS and Fermi-LAT in 2016-2017 and constraints on an intrinsic cut-off." pith.science (2026). https://pith.science/paper/K6QDIMOA
@misc{pith2026190806826,
author = {Pith},
title = {Pith review of: Observations of blazar PKS 2023-07 in flaring state with HESS and Fermi-LAT in 2016-2017 and constraints on an intrinsic cut-off},
year = {2026},
howpublished = {\url{https://pith.science/paper/K6QDIMOA}},
note = {Machine review of arXiv:1908.06826}
}
abstract
PKS 2023-07 is a flat spectrum radio quasar located at a redshift $z=1.388$, farther than any source currently detected at very high energies ($E>100$ GeV). At such energies, absorption by the extragalactic background light (EBL) renders the detection of distant sources particularly challenging. The High Energy Stereoscopic System (H.E.S.S.) observed the source following reports from AGILE (April 2016) and Fermi-LAT (April 2016, September and October 2017) on high-flux states in gamma rays. During each of the three flaring periods, near-simultaneous observations were obtained with H.E.S.S., Fermi-LAT and multiple telescopes at other wavelengths. Though the source was not significantly detected by H.E.S.S., upper limits were derived for each observation period. Through constraints given by Fermi-LAT in the MeV--GeV domain and differential upper limits by H.E.S.S., we searched for an intrinsic cutoff in the EBL-corrected gamma ray spectrum of PKS 2023-07.
Figures
Reference graph
Works this paper leans on
-
[1]
M. G. Hauser and E. Dwek, The Cosmic Infrared Background: Measurements and Implications, ARA&A 39 (2001) 249 [astro-ph/0105539]
arXiv 2001
- [2]
- [3]
-
[4]
P. Padovani, D. M. Alexander, R. J. Assef, B. De Marco, P. Giommi, R. C. Hickox et al., Active galactic nuclei: what’s in a name?, A&A~Rev. 25 (2017) 2 [1707.07134]
arXiv 2017
- [5]
-
[6]
S. Ciprini and Fermi Large Area Telescope Collaboration, Fermi-LAT detection of a GeV gamma-ray flare from the blazar PKS 2023-07, The Astronomer’s Telegram8932 (2016) . 6 PKS 2023−07 in flaring state with H.E.S.S. and Fermi-LAT in 2016–2017 G. Emery
work page 2016
-
[7]
J.-P. Lenain, FLaapLUC: A pipeline for the generation of prompt alerts on transient Fermi-LATγ-ray sources, Astronomy and Computing 22 (2018) 9 [1709.04065]
arXiv 2018
-
[8]
M. de Naurois and L. Rolland, A high performance likelihood reconstruction of γ-rays for imaging atmospheric Cherenkov telescopes, Astroparticle Physics 32 (2009) 231 [0907.2610]
arXiv 2009
Show all 22 references
-
[9]
Murach, M
T. Murach, M. Gajdus and R. D. Parsons, A Neural Network-Based Monoscopic Reconstruction Algorithm for H.E.S.S. II, arXiv e-prints (2015) arXiv:1509.00794 [1509.00794]
2015 arXiv
-
[10]
W. B. Atwood, A. A. Abdo, M. Ackermann, W. Althouse, B. Anderson, M. Axelsson et al., The Large Area Telescope on the Fermi Gamma-Ray Space Telescope Mission, ApJ 697 (2009) 1071 [0902.1089]
2009 arXiv
-
[11]
Acero, M
F. Acero, M. Ackermann, M. Ajello, A. Albert, W. B. Atwood, M. Axelsson et al., Fermi Large Area Telescope Third Source Catalog, ApJS 218 (2015) 23 [1501.02003]
2015 arXiv
-
[12]
The Fermi-LAT collaboration, Fermi Large Area Telescope Fourth Source Catalog, arXiv e-prints (2019) arXiv:1902.10045 [1902.10045]
2019 arXiv
-
[13]
Gehrels, G
N. Gehrels, G. Chincarini, P. Giommi, K. O. Mason, J. A. Nousek, A. A. Wells et al., The Swift Gamma-Ray Burst Mission, ApJ 611 (2004) 1005 [astro-ph/0405233]
2004 arXiv
-
[14]
D. N. Burrows, J. E. Hill, J. A. Nousek, J. A. Kennea, A. Wells, J. P. Osborne et al., The Swift X-Ray Telescope, Space~Sci.~Rev. 120 (2005) 165 [astro-ph/0508071]
2005 arXiv
-
[15]
P. M. W. Kalberla, W. B. Burton, D. Hartmann, E. M. Arnal, E. Bajaja, R. Morras et al., The Leiden/Argentine/Bonn (LAB) Survey of Galactic HI. Final data release of the combined LDS and IAR surveys with improved stray-radiation corrections, A&A 440 (2005) 775 [astro-ph/0504140]
2005 arXiv
-
[16]
K. C. Chambers, E. A. Magnier, N. Metcalfe, H. A. Flewelling, M. E. Huber, C. Z. Waters et al., The Pan-STARRS1 Surveys, arXiv e-prints (2016) arXiv:1612.05560 [1612.05560]
2016 arXiv
-
[17]
Domínguez, J
A. Domínguez, J. R. Primack, D. J. Rosario, F. Prada, R. C. Gilmore, S. M. Faber et al., Extragalactic background light inferred from AEGIS galaxy-SED-type fractions, MNRAS 410 (2011) 2556 [1007.1459]
2011 arXiv
-
[18]
W. A. Rolke, A. M. López and J. Conrad, Limits and confidence intervals in the presence of nuisance parameters, Nuclear Instruments and Methods in Physics Research A 551 (2005) 493 [physics/0403059]
2005 arXiv
-
[19]
Abdollahi, M
Fermi-LAT Collaboration, S. Abdollahi, M. Ackermann, M. Ajello, W. B. Atwood, L. Baldini et al., A gamma-ray determination of the Universe’s star formation history, Science 362 (2018) 1031 [1812.01031]
2018 arXiv
-
[20]
Madejski and M
G. . Madejski and M. Sikora, Gamma-Ray Observations of Active Galactic Nuclei, ARA&A 54 (2016) 725
2016
-
[21]
J. D. Finke, External Compton Scattering in Blazar Jets and the Location of the Gamma-Ray Emitting Region, ApJ 830 (2016) 94 [1607.03907]
2016 arXiv
-
[22]
Meyer, J
M. Meyer, J. D. Scargle and R. D. Blandford, Characterizing the Gamma-Ray Variability of the Brightest Flat Spectrum Radio Quasars Observed with the Fermi LAT, ApJ 877 (2019) 1 [1902.02291]. 7
2019 arXiv
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.