{"id":"c125d4a5-1ad8-4022-9df6-ca770e4f248a","arxiv_id":"1908.06826","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"H.E.S.S. non-detections of the z=1.388 blazar PKS 2023-07 during flares, together with Fermi-LAT spectra, imply an intrinsic cutoff below about 35-128 GeV or absorption in the broad line region.","lead":"This paper reports that H.E.S.S. did not detect gamma rays above 100 GeV from the distant quasar PKS 2023-07 during three flares in 2016 and 2017, and derives upper limits on its very high energy emission. The non-detection, combined with Fermi-LAT spectra, suggests the quasar's emission either cuts off at tens of GeV or is absorbed by gas near the supermassive black hole, giving new constraints on this distant blazar.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"H.E.S.S. upper limits are computed for a Γ=-3 power law but then used to constrain cutoff spectra without forward-folding; the derived E_cut limits may shift.","rationale":"The reader's weakest_assumption identifies the same soft spot, and it is genuinely load-bearing: every number in Section 4 derives from the H.E.S.S. upper limits, and those limits are spectral-shape dependent. My stress-test sharpens the reader's point into a forward-folding consistency test, which is the specific check needed before the E_cut constraints can be taken as reliable. The existing CONDITIONAL verdict is appropriate; I would not change it on the basis of this pass. I did not find a stronger concern that would overturn the argument: the Fermi-LAT power-law assumption is explicitly tested against a log-parabolic shape, the EBL scenario is explored by varying the EBL normalization, and the BLR interpretation is clearly labeled as depending on assumed L_disk and black-hole mass. One smaller internal inconsistency is that r = 1.8×10^17 cm corresponds to roughly 1200 r_g, not 120 r_g, for M = 10^9 M_sun, but this affects only the secondary BLR branch, not the intrinsic-cutoff constraint.","tokens_in":8379,"tokens_out":22040,"duration_ms":240986,"concrete_test":"Recompute the H.E.S.S. constraints by forward-folding each tested model (Fermi-LAT power law + EBL + exponential cutoff, with the photon index fixed to the LAT best fit and also varied over its 1σ range) through the H.E.S.S. mono-mode response, deriving the 95% C.L. exclusion on E_cut directly from the binned likelihood instead of using the Γ=-3 differential upper limits and a Gaussian translation. If the resulting E_cut upper limits for the April 2016 and September 2017 flares move by more than about 20%, or if the 95% incompatibility with a no-cutoff power law reverses for either flare, the central quantitative claim is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"All quantitative constraints in Section 4 (E_cut < 128 GeV for April 2016, E_cut < 35 GeV for September 2017, and the EBL normalization factors) rest on the H.E.S.S. 95% C.L. upper limits of Section 2.1. Those limits are explicitly derived assuming a source spectrum that is a power law with photon index -3. In Section 4, the authors convert these fixed-shape limits into a Gaussian likelihood and compare them with predicted fluxes from Fermi-LAT power laws (generally with different photon indices) modified by EBL absorption and an exponential cutoff. For an IACT, the mapping from a count upper limit to a differential flux at a quoted energy depends on the assumed spectral shape through the energy threshold, effective area, and energy-bin migration; the effect is largest for the lowest bin near threshold (~0.11 TeV), which carries much of the constraining power. The analysis is therefore not self-consistent: the H.E.S.S. data are not re-analyzed for each trial spectral shape. If the true VHE spectrum is steeper than -3, or contains a cutoff within the H.E.S.S. band, the differential upper limits at the quoted energies change, and the derived E_cut upper limits and the claimed incompatibility with a no-cutoff power-law extrapolation can shift.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":8642,"tokens_out":4793,"duration_ms":49392,"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":[{"comment":"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.","section":"Sec. 2.1 and Sec. 4"},{"comment":"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.","section":"Sec. 4"},{"comment":"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.","section":"Table 2"},{"comment":"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.","section":"Sec. 4, BLR absorption paragraph"}],"minor_comments":[{"comment":"The abbreviation 'FRSQ' should be 'FSRQ' (flat-spectrum radio quasar).","section":"Introduction"},{"comment":"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.'","section":"Sec. 2.1"},{"comment":"The text refers to the 'Thompson regime'; the correct spelling is 'Thomson regime'.","section":"Sec. 4"},{"comment":"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.","section":"Sec. 4"}],"recommendation":"major_revision","confidential_remarks":"This is a conference proceedings paper (ICRC2019) rather than a full journal article. The central observational result—non-detection and the tension with power-law extrapolation—is interesting and likely correct, but the spectral-shape dependence of the H.E.S.S. upper limits and the incompletely specified likelihood conversion are load-bearing for the quantitative E_cut and r constraints. The authors should be given the opportunity to clarify or correct these points. The overlap in authorship with Ref. [22] (Meyer et al. 2019) is not problematic because the BLR model is used as an external input and the analysis does not reduce to that model."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: this paper gives the first VHE upper limits for PKS 2023-07, the highest-redshift FSRQ, and the derived cutoff constraints are interesting but not as precise as the numbers suggest. The H.E.S.S. limits are computed for a -3 spectral index and then used to constrain other spectral shapes without refitting.\n\nThe new thing is the first VHE non-detections for this z=1.388 FSRQ during three flaring epochs, plus constraints from a profile-likelihood comparison with Fermi-LAT spectra: an intrinsic cutoff below 128 GeV (April 2016) and below 35 GeV (September 2017), or BLR absorption if the emission region sits within about 120-650 gravitational radii. That is a legitimate single-source result, and the paper is honest about its assumptions.\n\nWhat it does well: the H.E.S.S. analysis is run with two independent chains (Model and ImPACT) with consistent results; the Fermi-LAT analysis is standard and checks whether a log-parabola is needed; the multi-wavelength light curves provide useful context. The non-detection is not oversold.\n\nThe soft spots are real but not fatal. The H.E.S.S. differential upper limits are computed assuming a photon index of -3 (Section 2.1, Table 1). In Section 4 those fixed-shape limits are translated into a Gaussian likelihood and compared with predicted fluxes from Fermi-LAT power laws with different indices, modified by EBL absorption and an exponential cutoff. For an IACT, the conversion from counts to flux at a given energy depends on the assumed spectrum through the threshold and effective area; the lowest energy bin around 0.11 TeV carries much of the constraining power, and that is where the dependence is strongest. So the E_cut and EBL-normalization results are not self-consistent: the H.E.S.S. data are not re-fit for each trial spectral shape. The direction and size of the shift are not quantified. I would trust the qualitative conclusion (simple power-law extrapolation is ruled out for April and September) but treat the specific E_cut numbers as indicative, not precise.\n\nThe BLR absorption limits are conditional on assumed values for L_disk, L_Hbeta, and black hole mass, which the authors state they could not find in the literature. That is fine as an exploratory statement but should not be over-read. The citation pattern is okay; the BLR model from Meyer et al. 2019 has overlapping authorship, but it is a cited model, not a hidden assumption.\n\nWho is this for? Observational gamma-ray astronomers working on FSRQs, EBL, and the location of the VHE emission region. It is a modest single-source paper, but it is a legitimate data point for the highest-redshift FSRQ targeted by an IACT. I would send it to a referee if it were submitted as a journal paper; the proceedings format already gives some vetting, but the spectral-index issue deserves a comment.\n\nBottom line: accept with minor revisions, but ask the authors to re-derive the limits for a range of intrinsic indices or at least discuss the systematic. As it stands, the key numbers are provisional.","headline":"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.","tokens_in":9264,"tokens_out":4137,"would_cite":false,"duration_ms":38968,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["PKS 2023-07","flat spectrum radio quasar","very high energy gamma rays","H.E.S.S.","Fermi-LAT","intrinsic spectral cutoff","extragalactic background light","broad line region absorption"],"falsifier":"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.","tokens_in":8154,"feed_emoji":"🔭","tokens_out":15295,"duration_ms":140700,"temperature":0.7,"pith_summary":"PKS 2023-07 is a flat-spectrum radio quasar, a blazar whose relativistic jet points toward Earth, at redshift 1.388 — more distant than any blazar yet detected at very high energies. The paper combines H.E.S.S. upper limits from three flares (April 2016, September 2017, October 2017) with near-simultaneous Fermi-LAT spectra to test whether the GeV flare spectrum can simply extend into the TeV band. For the two brightest flares, that simple power-law extension, after standard EBL absorption, predicts more very-high-energy flux than the H.E.S.S. 95% confidence upper limits allow. The authors conclude that the source must either have an intrinsic exponential cutoff, with $E_\\mathrm{cut}<128\\ \\mathrm{GeV}$ in April 2016 and $E_\\mathrm{cut}<35\\ \\mathrm{GeV}$ in September 2017, or absorb its TeV photons in the broad-line region, which requires the emission region to lie within roughly 120--650 gravitational radii of the black hole. The result matters because it places a very distant blazar under direct constraints on particle acceleration and on where gamma-ray flares are produced.","feed_headline":"Missing TeV light from quasar flare points to cutoff below 35 GeV","feed_subtitle":"H.E.S.S. limits clash with Fermi-LAT extrapolations, forcing an intrinsic cutoff or absorption near the black hole.","key_machinery":"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.","core_discovery":"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).","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the extragalactic background light model used for the EBL attenuation and normalisation test.","marker":"[17]"},{"why":"Supplies the broad-line-region photon field and scaling relations used for BLR absorption and external-Compton electron energy estimates.","marker":"[21]"},{"why":"Provides the treatment of gamma-ray absorption by broad-line-region photons that the paper follows to constrain the emission-region distance.","marker":"[22]"},{"why":"Provides the profile-likelihood-ratio method used to convert H.E.S.S. upper limits into confidence-level exclusions.","marker":"[18]"},{"why":"Provides the external-Compton framework used to translate the intrinsic cutoff bound into electron Lorentz factor bounds.","marker":"[20]"},{"why":"Describes the H.E.S.S. monoscopic reconstruction chain used to derive the differential upper limits.","marker":"[8]"}],"fun_headline_variants":["No TeV detection in quasar flare hints at intrinsic cutoff below 35 GeV","Fermi-LAT extrapolation overshoots H.E.S.S. limits, forcing cutoff","Quasar PKS 2023-07 flare: H.E.S.S. silence sets cutoff at 35 GeV","TeV darkness from quasar flare points to absorption or cutoff","Lack of TeV light in flare constrains quasar emission zone"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["No TeV detection in quasar flare hints at intrinsic cutoff below 35 GeV","Fermi-LAT extrapolation overshoots H.E.S.S. limits, forcing cutoff","Quasar PKS 2023-07 flare: H.E.S.S. silence sets cutoff at 35 GeV","TeV darkness from quasar flare points to absorption or cutoff","Lack of TeV light in flare constrains quasar emission zone"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0003,"raw_usage":{"total_tokens":1827,"prompt_tokens":1134,"completion_tokens":693,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":750,"completion_tokens_details":{"reasoning_tokens":583}},"tokens_in":750,"tokens_out":693,"duration_ms":7080,"temperature":1.0,"reasoning_tokens":583,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:33:06.974451+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"External Compton Scattering in Blazar Jets and the Location of the Gamma-Ray Emitting Region","cited_arxiv_id":"1607.03907","evidence_quote":"Supplies the broad-line-region photon field and scaling relations used for BLR absorption and external-Compton electron energy estimates."},{"cited_title":"Madejski and M","cited_arxiv_id":null,"evidence_quote":"Provides the external-Compton framework used to translate the intrinsic cutoff bound into electron Lorentz factor bounds."}],"review_version":1}