Towards a TeV blazar sequence and its physical interpretation

C. Righi; E. Prandini; F. Bovolon; I. Viale; N. Sahakyan

REVIEW 3 major objections 5 minor 1 cited by

Towards a TeV blazar sequence and its physical interpretation

T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

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

Pith's one-line read For PKS 2155-304, choosing which X-ray data count as 'average' doubles the fitted Doppler factor and flips the magnetic field by three orders of magnitude.

desk verdict Honest progress note on a frequency-binned TeV blazar sequence, but the one quantitative result is undermined by SSC degeneracy the authors themselves flag. read the letter →

arxiv 2506.18500 v1 pith:NMS3WQGY submitted 2025-06-23 astro-ph.HE

I. Viale , E. Prandini , C. Righi , F. Bovolon , N. Sahakyan This is my paper · ORCID

classification astro-ph.HE

keywords blazarsBLLacobjectsblazarsequenceSynchrotronSelfComptonPKS2155-304TeVspectralenergydistributionX-rayvariability

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 is an early step toward explaining the blazar sequence—the anticorrelation between blazar luminosity and the frequency of their spectral peaks—in physical terms, by modeling TeV-detected BL Lac objects one by one in a synchrotron self-Compton (SSC) framework. For the representative high-energy-peaked source PKS 2155-304, the authors show that the fitted SSC parameters depend strongly on how the archival X-ray data are selected: averaging the X-ray data versus using single observations in an average state yields a Doppler factor that doubles, from 21.47 to 43.06, and a magnetic field that changes by roughly three orders of magnitude. Since both selections are reasonable and the source is highly variable in X-rays, the message is that an 'average-state' SED of a TeV blazar is not well-defined from archival data, and this ambiguity must be resolved before any physical interpretation of the TeV blazar sequence can be trusted. The wider goal is to model the remaining representative sources in five synchrotron-peak bins and look for trends in the fitted jet parameters.

What carries the argument

The machinery is the one-zone Synchrotron Self-Compton (SSC) model, in which a single relativistic blob of electrons produces both the low-energy synchrotron bump and the high-energy bump by up-scattering those same synchrotron photons through inverse Compton. The model is applied through an online fitting tool whose convolutional neural network, trained on spectra simulated with the SOPRANO code, returns best-fit parameters and corner plots over the free parameters (electron luminosity, injection index, Doppler factor, magnetic field, blob size, variability timescale). The contrast between two X-ray data selections for PKS 2155-304 is what carries the argument: the same model, the same source, and two defensible data choices produce widely separated parameter values.

What would settle it

Run both SED selections through the full SOPRANO simulation code rather than the convolutional-neural-network emulator and re-fit; if the best-fit Doppler factor and magnetic field no longer show the separation in Table 2, the claimed data-selection sensitivity would be an artifact of the surrogate model. In parallel, an SED built from strictly simultaneous X-ray and very-high-energy observations of PKS 2155-304 taken during a single quiescent epoch would test which archival selection, if either, approximates the source's true average state.

Watch

Extended reading notes

Core claim

The paper's central claim, stated on its own terms, is that data selection materially affects the modeled jet properties of PKS 2155-304 within a simple SSC model. Using the MMDC fitting tool, the authors fit the same source twice: once with averaged X-ray data and once with single X-ray observations corresponding to an average activity state. The resulting best-fit parameters are substantially different—most notably the Doppler factor (δ = 21.47 vs 43.06) and the magnetic field (log(B/G) = 1.70 vs −1.51)—while the electron luminosity and energy-break parameters stay close. The authors attribute the large Doppler-factor shift to a poorly constrained parameter distribution, and note that it is not trivial to say which selection is correct. Thus the paper establishes a caution: any physical reading of the blazar sequence from SSC fits of archival SEDs is sensitive to the choice of which observations represent the source.

Load-bearing premise

The load-bearing assumption is that the non-simultaneous archival data selected for PKS 2155-304, under either X-ray treatment, truly represents the source in an average activity state; if neither selection does, the fitted jet parameters are not physical.

Editorial extensions

If this is right

The planned fits of representative sources in the other synchrotron-peak bins will only support a physical interpretation of the blazar sequence if the fitted parameters are stable against equivalent data-selection choices; the PKS 2155-304 result puts the onus on showing that stability.
Reported SSC parameters for PKS 2155-304 should carry a note about the X-ray selection: the Doppler factor and magnetic field cannot be quoted as single numbers without specifying whether X-ray data were averaged or taken from single observations.
Because the Doppler factor distribution is poorly constrained, conclusions that use the fitted δ for this source—for example, estimates of jet beaming or jet power—should be treated as provisional.
Future continuous X-ray monitoring of TeV BL Lacs, or tools that filter observations by epoch, would replace the two-selection ambiguity with a well-defined average state.

Reading between the lines

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

If the sensitivity seen in PKS 2155-304 holds for other representatives, then the trend that a future 'TeV blazar sequence' reveals could be partly an artifact of how each source's archival X-ray data are selected; binning by synchrotron peak frequency does not by itself control for this.
A natural stress test would be to fit multiple alternative realizations of the same source—say X-ray data split by flux state—and check whether the inferred SSC parameters trace the activity state continuously; if they do, the fits could be inverted to estimate the state of the source at the time of each observation.
The large swing in magnetic field from log(B/G) = 1.70 to −1.51, if real, has consequences beyond the sequence: it changes estimates of the jet's particle energy budget and of the γ-ray emission region's location, so the data-selection ambiguity could affect TeV BL Lac population studies generally.

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.

Desk Editor's Note

read the letter

This is a short proceedings paper, and it reads like one: a clear description of a planned program, one worked example, and a heavily caveated result. The frequency-binning of the blazar sequence is a sensible alternative to luminosity binning, and the paper is refreshingly honest that the two X-ray selections give quite different best-fit parameters—and that the Doppler factor, in particular, is poorly constrained.

What is genuinely new: binning by synchrotron peak frequency rather than luminosity, and choosing representative sources rather than averaging SEDs. That is a reasonable methodological twist and worth keeping in mind. The comparison of averaged versus single-observation X-ray data for PKS 2155-304 is also a useful caution for anyone building SEDs from archival data, and the corner plot lets the reader see the degeneracies directly.

The soft spots are real, though. The paper's central claim—that data selection materially changes the fitted jet parameters—is not supported by the numbers. Table 2 gives point estimates with no error bars. The Doppler factor shifts from 21.47 to 43.06, but the authors themselves say it is poorly constrained and the results are less reliable. The magnetic field swings from log B = 1.70 to -1.51, a 3.2 dex jump, which is exactly what you expect if the fit is just sliding along the degenerate B–δ ridge in a one-zone SSC model. So the two ``best-fit'' points may correspond to the same physical solution, not two distinct states. The stress-test note is right about this, and the manuscript's own admission about δ strengthens the point.

There is also a mismatch between the title and the content: no TeV blazar sequence is delivered here, just one source and a plan. That is fine for a proceedings, but the framing oversells. The MMDC circularity is a minor issue—the tool is published, but it was built by the same group, so independent cross-checks with a traditional fitting code would add confidence. The use of non-simultaneous archival data for a highly variable source also makes the ``average state'' concept fragile, as the authors concede.

Bottom line: as a progress note, this is acceptable and should not be desk-rejected. But it should be treated as a description of work in progress, not as evidence that X-ray selection changes the physics. The referee should ask for credible intervals or an explicit statement that the two selections cannot be distinguished given the degeneracies. A light referee is enough; a heavy one is overkill.

Referee Report

3 major / 5 minor

Summary. This proceedings paper proposes a new way to revisit the blazar sequence by binning TeV-detected BL Lacs into one-decade synchrotron-peak-frequency bins and selecting one representative source per bin, then modeling each source with a one-zone synchrotron self-Compton (SSC) model. The paper focuses on PKS 2155-304, the representative source for the bin 10^15 < nu_syn < 10^16 Hz, and compares SSC fits obtained with two different X-ray data selections: an average of archival X-ray data versus a single X-ray observation. The resulting best-fit parameters (Table 2) differ substantially, most notably in the Doppler factor delta (21.47 versus 43.06) and the magnetic field log B (1.70 versus -1.51), and the authors interpret this as evidence that data selection materially affects the modeled jet properties, with implications for the physical interpretation of the TeV blazar sequence. The paper explicitly acknowledges that delta is poorly constrained and that the source is highly variable, and it frames the work as a preliminary step toward modeling a larger sample.

Significance. If validated, the proposed frequency-binned selection protocol would provide a systematic framework for modeling TeV BL Lacs in the context of the blazar sequence, and the emphasis on average-state data would address a known weakness of archival SED modeling. The paper's strengths include its transparent presentation of the source-selection criteria, its use of publicly available databases (4LAC-DR2, MMDC, STeVECat), and its explicit caveats about X-ray variability and the poor constraint on delta. The manuscript also honestly reports the tension between two reasonable data selections rather than hiding it. However, the central quantitative claim--that the X-ray selection changes the physical jet parameters--currently rests on point estimates without uncertainty quantification, and the paper itself flags the degeneracy-related unreliability of delta. Because the paper aims at a physical interpretation of the blazar sequence, the absence of credible intervals or a model-comparison test limits the evidentiary weight of the reported parameter shifts.

major comments (3)

[Sect. 5 and Table 2] The central claim that the two X-ray selections yield materially different SSC parameters is not robustly established because Table 2 reports only point estimates without uncertainties or credible intervals. The text itself states that the difference in delta 'can be ascribed to its poorly constrained distribution in the modeling (see Fig. 4b), which makes these results less reliable,' and the change in log B from 1.70 to -1.51 (3.2 dex) is presented without demonstrating that log B is independently constrained. Since one-zone SSC models are known to have strong degeneracies among B, delta, and source size, the two best-fit points could lie on the same degenerate ridge; a quantitative comparison of the posterior distributions (e.g., credible intervals, profile likelihoods) or fits with delta fixed to a grid of values is needed to attribute the parameter shifts to the data selection rather than to model degeneracy.
[Sect. 3] The claim that the assembled non-simultaneous SEDs represent an 'average activity state' is not demonstrated quantitatively. The paper acknowledges that simultaneous multi-wavelength data were not always available and that X-ray observations were gathered from different periods, but no criterion (e.g., X-ray flux relative to the long-term distribution, or an activity-state indicator) is given for either the averaged selection or the selected single observations. Given the high variability of PKS 2155-304 shown in Fig. 3, if neither selection actually corresponds to the average state, the fitted jet parameters and the implication for the TeV blazar sequence are not physically meaningful. The authors should either validate the selected states against the source's long-term X-ray light curve or frame the results explicitly as bracketing states rather than average-state properties.
[Sect. 4] The modeling tool MMDC is trained by the same group (including a co-author) on SOPRANO simulations, and no independent cross-check against other SSC fitting codes or Monte Carlo methods is presented. Because the paper aims at a physical interpretation, the absolute parameter values in Table 2 inherit any systematic bias of the training set; a validation run with an independent code (e.g., a standard Markov chain Monte Carlo implementation) on at least PKS 2155-304 would demonstrate that the reported differences are not artifacts of the specific CNN architecture or training data.

minor comments (5)

[General] The in-text references to table numbers are inconsistent: Sect. 4 refers to 'Table 4' for the modeling results while the displayed table is labeled Table 2, and Sect. 2 refers to 'Table 2' for representative sources while the displayed table is labeled Table 1.
[Sect. 3] There is a typo in the sentence 'VHEγ-ray data were taken from from the Spectral TeV Extragalactic Catalog'--the word 'from' is duplicated.
[Sect. 4] The word 'Netowrks' in the description of MMDC is misspelled and should be 'Networks'.
[Sect. 5] The phrase 'continuous of sources in X-rays' is unclear and should be reworded, for example as 'continuous monitoring of sources in X-rays'.
[Fig. 4] The corner plot in Fig. 4 would be easier to interpret if the parameter ranges and contours were labeled with units (e.g., delta is dimensionless, B in Gauss), and if the posterior distributions were overlaid for both X-ray selections.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the reported parameter differences are empirical fit outcomes, not predictions forced by construction or by self-citation.

full rationale

The paper makes no derived prediction; it reports best-fit SSC parameters obtained by fitting two archival data selections for PKS 2155-304. The qualitative result that averaged versus single-observation X-ray data yield different best-fit parameters is an empirical outcome of the fitting procedure and does not reduce to the inputs by construction. The same tool (MMDC) and the same model were applied to both selections, so the comparison is not an artifact of the fitting code. The authors explicitly flag that delta is poorly constrained and that these results are less reliable, which is a degeneracy/uncertainty caveat rather than circularity. The only self-citation is to MMDC (Sahakyan et al. 2024) and to the SOPRANO simulation code; this is transparent methodology, not a load-bearing uniqueness theorem, and the paper's central claim does not depend on the tool's internal validity. No equation is equivalent to another by definition, and no fitted quantity is renamed as a prediction. Lack of independent validation of MMDC is a validation concern, not circularity.

Assumptions & free parameters 6 free parameters · 3 assumptions · 0 invented entities

The analysis rests on three categories of assumptions: the validity of one-zone SSC for this source, the representativeness of non-simultaneous data as an average state, and the as-yet-unvalidated choice of one representative source per frequency bin. No new physical entities are introduced.

free parameters (6)

t_var (variability timescale) = 4.88e5 s (X-ray average), 4.86e4 s (single obs)
Best fit from MMDC SSC model; strongly dependent on X-ray selection.
p1 (electron index) = 2.24 (average), 2.34 (single)
Best fit from MMDC SSC model; dependent on X-ray selection.
log(Le) (electron luminosity) = 44.90 (average), 44.77 (single) erg/s
Best fit from MMDC SSC model; dependent on X-ray selection.
log(gamma_max) = 5.45 (average), 5.28 (single)
Best fit from MMDC SSC model; dependent on X-ray selection.
delta (Doppler factor) = 21.47 (average), 43.06 (single)
Best fit from MMDC SSC model; poorly constrained per Fig. 4b.
log(B) (magnetic field) = 1.70 (average), -1.51 (single) [G]
Best fit from MMDC SSC model; strongly dependent on X-ray selection.

assumptions (3)

domain assumption One-zone SSC with homogeneous emitting region and electron distribution described by a smooth broken power law describes PKS 2155-304 emission.
Invoked in Sect. 4 via MMDC; external photon fields and more complex geometries are neglected.
ad hoc to paper The non-simultaneous multi-wavelength data assembled in Sect. 3 represent the average activity state of the source.
The paper gathers data from different epochs and either averages X-rays or selects single observations; this is not independently justified.
ad hoc to paper Splitting TeV BL Lacs into one-decade synchrotron peak frequency bins with one representative source per bin preserves the blazar sequence information.
Proposed in Sect. 2 but not yet validated; the current paper only models one of the five representatives.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Towards a TeV blazar sequence and its physical interpretation." pith.science (2026). https://pith.science/paper/NMS3WQGY

@misc{pith2026250618500,
  author       = {Pith},
  title        = {Pith review of: Towards a TeV blazar sequence and its physical interpretation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NMS3WQGY}},
  note         = {Machine review of arXiv:2506.18500}
}

read the original abstract

Blazars, a highly energetic subclass of jetted active galactic nuclei, show a broad band spectral energy distribution (SED) with two bumps, resulting from non-thermal jet emission. In 1998, an anticorrelation between the SED luminosity and the peaks frequency was found, later confirmed in 2017, called the blazar sequence. Its origin is still unclear. This work is part of a broader effort aiming at giving a physical interpretation to the blazar sequence, by modeling the sources emission in a Synchrotron Self Compton framework, concentrating only on TeV-detected blazars of BL Lac type. Unlike the original sequence, sources were binned by synchrotron peak frequency. The SED of one representative source is modeled for each bin, using data from average activity state to ensure consistency. Here, we focus on the study of one of the selected representative sources, PKS 2155-304, a high-energy-peaked BL Lac, showing the performed data selection and preliminary modeling results.

Figures

Figures reproduced from arXiv: 2506.18500 by the authors.

**Figure 1.** Synchrotron peak frequency distribution of 4LAC-DR2 blazars. All sources are shown in grey, FSRQs in orange, BL Lacs in blue, TeV-detected BL Lacs in cyan. modeling of the sources emission, getting a final sample of 55 objects. In previous works on the blazar sequence, sources were divided into bins based on their luminosity in a pre-defined band. Then, the single-blazar SEDs in each bin were averaged in order to g… view at source ↗

**Figure 2.** Synchrotron peak frequency-based distribution and bins defined for our study. Only TeV BL Lac objects are shown. The vertical green lines indicate the νsyn of the selected representative sources. Bin [Hz] Source νsyn [Hz] Bin 1 νsyn < 1014 BL Lac 3.86 × 1013 Bin 2 1014 < νsyn < 1015 TXS 0506+056 3.55 × 1014 Bin 3 1015 < νsyn < 1016 PKS 2155-304 5.69 × 1015 Bin 4 1016 < νsyn < 1017 PG 1218+304 1.85 × 1016 Bin 5 νsyn … view at source ↗

**Figure 3.** Data selection for the source PKS 2155- 304. representative of the third bin (1015 < νsyn < 1016). ternal photon fields in their environment, making the effect of the External Compton process more prominent and necessary to well interpret their emission. Therefore, in the following, we will model only the representative BL Lacs from the three higher-frequency bins and compare their results. To perform the source … view at source ↗

Figures from the paper (1 more)

**Figure 4.** Figure 4: SSC model (left) and related corner plot (right) of the source PKS 2155-304. In this SED the data selection comprises single X-ray observations. References Ajello, M., Angioni, R., Axelsson, M., et al. 2020, ApJ, 892, 105 Atwood, W. B., Abdo, A. A., Ackermann, M., et a…

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

Radio-Gamma-Ray Properties and High-Energy Implications for Fermi Blazars
astro-ph.HE 2026-08 conditional novelty 6.0 of 10

Using 1,687 Fermi blazars, the paper reports a plateau in gamma-ray loudness at high synchrotron peak frequencies, interpreted as Klein-Nishina suppression yielding magnetic field constraints for HBLs.

Reference graph

Works this paper leans on

13 extracted references · 8 canonical work pages · cited by 1 Pith paper

[1]

, " * write output.state after.block = add.period write newline

ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...
[2]

write newline

" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....
[3]

2020, , 892, 105

Ajello , M., Angioni , R., Axelsson , M., et al. 2020, , 892, 105

work page 2020
[4]

B., Abdo , A

Atwood , W. B., Abdo , A. A., Ackermann , M., et al. 2009, , 697, 1071

2009
[5]

2024, Astrophys

B e \' e gu e \' e , D., Sahakyan, N., Dereli-B e \' e gu e \' e , H., et al. 2024, Astrophys. J., 963, 71

work page 2024
[6]

1998, Monthly Notices of the Royal Astronomical Society, 299, 433

Fossati, G., Maraschi, L., Celotti, A., Comastri, A., & Ghisellini, G. 1998, Monthly Notices of the Royal Astronomical Society, 299, 433

work page 1998
[7]

2022, Mon

Gasparyan, S., B e \' e gu e \' e , D., & Sahakyan, N. 2022, Mon. Not. R. Astron. Soc., 509, 2102

work page 2022
[8]

2017, , 469, 255

Ghisellini , G., Righi , C., Costamante , L., & Tavecchio , F. 2017, , 469, 255

work page 2017

Show all 13 references

[9]

2023, arXiv e-prints, arXiv:2304.00835

Gr e \' e aux, L., Biteau, J., Hassan, T., et al. 2023, arXiv e-prints, arXiv:2304.00835

2023 arXiv
[10]

2020, arXiv e-prints, arXiv:2010.08406

Lott , B., Gasparrini , D., & Ciprini , S. 2020, arXiv e-prints, arXiv:2010.08406

2020 arXiv
[11]

1992, , 397, L5

Maraschi , L., Ghisellini , G., & Celotti , A. 1992, , 397, L5

1992
[12]

& Ghisellini, G

Prandini, E. & Ghisellini, G. 2022, Galaxies, 10, 35

2022
[13]

2024, Astron

Sahakyan, N., Vardanyan, V., Giommi, P., et al. 2024, Astron. J., 168, 289

2024

Pith tools

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

[1] [1]

, " * write output.state after.block = add.period write newline

ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...

[2] [2]

write newline

" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....

[3] [3]

2020, , 892, 105

Ajello , M., Angioni , R., Axelsson , M., et al. 2020, , 892, 105

work page 2020

[4] [4]

B., Abdo , A

Atwood , W. B., Abdo , A. A., Ackermann , M., et al. 2009, , 697, 1071

2009

[5] [5]

2024, Astrophys

B e \' e gu e \' e , D., Sahakyan, N., Dereli-B e \' e gu e \' e , H., et al. 2024, Astrophys. J., 963, 71

work page 2024

[6] [6]

1998, Monthly Notices of the Royal Astronomical Society, 299, 433

Fossati, G., Maraschi, L., Celotti, A., Comastri, A., & Ghisellini, G. 1998, Monthly Notices of the Royal Astronomical Society, 299, 433

work page 1998

[7] [7]

2022, Mon

Gasparyan, S., B e \' e gu e \' e , D., & Sahakyan, N. 2022, Mon. Not. R. Astron. Soc., 509, 2102

work page 2022

[8] [8]

2017, , 469, 255

Ghisellini , G., Righi , C., Costamante , L., & Tavecchio , F. 2017, , 469, 255

work page 2017