Pith. sign in

REVIEW 3 major objections 4 minor 16 references

On the TeV gamma-ray and X-ray correlations exhibited in high-energy peaked BL Lacs: Mrk 501 and 1ES 1959+650

T0 review · 3 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 For two more blazars, X-ray and TeV fluxes trace one straight line, with no high-flux breakdown.

desk verdict Preliminary two-source extension of the Mrk 421 correlation analysis; the Mrk 501 fit is solid, but the 'no break-down' claim for 1ES 1959+650 is untested because high-flux points were excluded. read the letter →

arxiv 1909.00897 v1 pith:DD2JREUK submitted 2019-09-03 astro-ph.HE

Mabel Osorio , Rodrigo Sacahui , M. Magdalena González , Nissim Fraija , José Andrés García-González This is my paper · ORCID
classification astro-ph.HE
keywords BLLacertaeobjectshigh-energypeakedLacsMarkarian5011ES1959+650X-ray/TeVfluxcorrelationsynchrotronself-ComptonmodelmaximumlikelihoodfittingTegamma-rayastronomy
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 tests whether the strong X-ray–TeV flux correlation seen in the blazar Mrk 421, which breaks down at high gamma-ray fluxes, is a general property of high-energy-peaked BL Lacs (HBLs). Using archival observations of Mrk 501 and 1ES 1959+650, the authors fit a straight line $F_\gamma = a F_X + b$ to each source's X-ray and TeV fluxes with a maximum-likelihood method that allows for intrinsic scatter. The fits are strong ($R^2 = 0.970$ for Mrk 501 and $R^2 = 0.889$ for 1ES 1959+650), and no high-flux point falls outside $3\sigma_d$ in either fitted sample. The authors interpret this as evidence that the high-flux break-down reported for Mrk 421 is not necessarily universal among HBLs, while noting that the 1ES 1959+650 fit excludes its strongest flare states.

What carries the argument

The central device is the flux–flux correlation plot, in which each quasi-simultaneous pair of X-ray and TeV gamma-ray flux measurements is located, and a straight line $F_\gamma = a F_X + b$ is fitted with a maximum-likelihood procedure that includes measurement errors on both axes and an intrinsic scatter $\sigma_d$ of unknown origin. The output of that fit — slope $a$, intercept $b$, scatter $\sigma_d$, and the associated $p$-value and $R^2$ — carries the argument: the claim that all points stay within $3\sigma_d$ of the line at high fluxes is the evidence for 'no break-down.' The synchrotron self-Compton (SSC) model supplies the interpretive setting, since a tight, single correlation is the signature expected when X-rays and TeV gamma-rays come from the same electron population, and it is what would later let the magnetic field be constrained under a flux-independent electron spectral index.

What would settle it

Recompute the 1ES 1959+650 correlation including the May 16–17, 2002 strong flare and the June 4 orphan TeV flare; if those points lie more than $3\sigma_d$ from the fitted line, the claimed absence of a high-flux break-down is falsified for that source.

Watch

Extended reading notes

Core claim

The paper's central claim is that the X-ray and TeV gamma-ray fluxes of Mrk 501 and 1ES 1959+650 follow a single linear relation $F_\gamma = a F_X + b$ over their observed ranges, with no break-down at high gamma-ray fluxes. For Mrk 501 the fitted parameters are $a = 194 \pm 9.20$, $b = -16.9 \pm 1.47$, $\sigma_d = 2.12 \pm 0.67$, with $R^2 = 0.970$; for 1ES 1959+650 they are $a = 71.10 \pm 6.48$, $b = 0.36 \pm 0.07$, $\sigma_d = 0.23 \pm 0.06$, with $R^2 = 0.889$. These correlations survive despite the X-ray and TeV measurements not being simultaneous, with average time offsets of about 30 minutes for Mrk 501 and 1.7 hours for 1ES 1959+650. The paper takes the absence of points beyond $3\sigma_d$ at high fluxes as evidence that the break-down seen in Mrk 421 is not a general property of HBLs, at least for the data sets analyzed.

Load-bearing premise

The analysis assumes that X-ray and TeV measurements taken up to about 1.7 hours apart, with error bars estimated from published plot symbol sizes, can stand in for truly simultaneous fluxes; for 1ES 1959+650 the strongest flare states are excluded, so the high-flux portion of the claim is untested.

Editorial extensions

If this is right

  • If the correlations are genuine, the high-flux break-down observed in Mrk 421 is not a universal feature of high-energy-peaked BL Lacs, at least among the sources studied here.
  • A single linear relation from low to high gamma-ray fluxes means that, for these two objects, one set of correlation parameters describes the full observed state, so model studies need not invoke separate spectral-index behaviors at different flux levels.
  • The correlations holding across average time offsets of 30 minutes to 1.7 hours suggest that quasi-simultaneous multiwavelength monitoring can recover the underlying X-ray/TeV relation for these HBLs.
  • For 1ES 1959+650, the low-state fit leaves the high-flux behavior unexplored, because the May 16–17 strong flare and the June 4 orphan TeV flare were deliberately excluded.

Reading between the lines

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

  • A direct extension would be to add other HBLs with archival X-ray and TeV monitoring to the same maximum-likelihood fit; a sample in which most sources keep a single line would make the Mrk 421 high-flux break-down look like the exception rather than the norm.
  • Because the two sources are measured in different energy bands and flux units, the fitted $a$ and $b$ values are not physically comparable as-is; recasting both data sets in common physical units would test whether the underlying correlation is the same mechanism.
  • Re-running the 1ES 1959+650 fit with the excluded flare points included would settle whether the no-breakdown behavior extends to high states or whether strong flares behave like the Mrk 421 high-flux outliers.
  • If truly simultaneous pairs with minute-scale binning were used, a smaller intrinsic scatter $\sigma_d$ would indicate that the 30-minute to 1.7-hour offsets add variance rather than trace the same physical state, which would strengthen or weaken the quasi-simultaneity assumption.
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

3 major / 4 minor

Summary. The paper extends a previous X-ray/TeV gamma-ray correlation analysis (originally performed for Mrk 421 by González et al. 2019) to two other high-energy peaked BL Lac objects, Mrk 501 and 1ES 1959+650. For each source the authors compile flux points from published multiwavelength campaigns (MAGIC/VERITAS/Swift for Mrk 501; Whipple/RXTE for 1ES 1959+650) and fit a linear relation F_gamma = a F_X + b with an intrinsic scatter sigma_d using the D'Agostini (2005) maximum-likelihood method. For Mrk 501 they report a = 194 +/- 9.20, b = -16.9 +/- 1.47, sigma_d = 2.12 +/- 0.67, R^2 = 0.970; for 1ES 1959+650 they report a = 71.10 +/- 6.48, b = 0.36 +/- 0.07, sigma_d = 0.23 +/- 0.06, R^2 = 0.889. The central claim in Section 5 is that, unlike Mrk 421, these two HBLs show no break-down of the correlation at high gamma-ray fluxes, and a single linear correlation holds from low to high fluxes. The paper explicitly labels these results preliminary and states that a larger sample and a study of quasi-simultaneity are planned.

Significance. If the claim of no high-flux break-down is correct for Mrk 501 and 1ES 1959+650, it would demonstrate that the Mrk 421 break-down is not a universal HBL phenomenon and would have direct implications for one-zone SSC modeling, particularly for the allowed range of magnetic field strengths when the correlation extends to high fluxes. The statistical method is clearly described and reproducible, and the authors are transparent about several data limitations, including lack of exact simultaneity and error bars estimated from symbol sizes. However, the central physical claim is currently not supported for one of the two sources, because the high-flux states of 1ES 1959+650 are excluded from the fit. The paper is nonetheless a useful preliminary contribution, and the Mrk 501 result alone is a valid data point. The strengths are the clean exposition of the likelihood procedure and the honest acknowledgment of caveats; the weakness is that the conclusion extrapolates beyond the fitted data.

major comments (3)
  1. [Sections 3 and 5; Figure 4 caption] The statement in Section 5 that 'we notice the lack of a break-down at high gamma-ray fluxes' and that a linear correlation holds 'from low to high gamma-ray fluxes' is not supported for 1ES 1959+650. Section 3 explicitly states that 'the correlation was calculated only in the low state' because the May 16-17, 2002 strong flare and the June 4, 2002 orphan flare were excluded, and the Figure 4 caption repeats that the data are 'for only a low state.' Consequently, the highest-flux points for this source are absent from the fit, so the absence of a break-down is an artifact of the selection, not a measured property. This is the load-bearing point for the paper's central conclusion, since the claim is about behavior at high fluxes. The authors should either re-fit 1ES 1959+650 including the excluded flare points (provided the X-ray counterpart is available and quasi-simultaneous pairing is possible), or explicitly restrict the no-break-down conclusion to Mrk 501 and revise the abstract and Section 5 accordingly.
  2. [Sections 2 and 3; Figure captions 1-4] The analysis pairs non-simultaneous X-ray and TeV gamma-ray observations, with average time offsets of about 30 minutes for Mrk 501 and 1.7 hours for 1ES 1959+650 as stated in the figure captions. The paper claims that the correlations 'exist even though there is not an exact time simultaneity,' but no test or quantitative assessment of the effect of these offsets on the fitted parameters is presented. If the intrinsic variability timescale is comparable to or shorter than the offset, the derived slope, intercept, and intrinsic scatter could be biased. The authors should either pair data with a stricter simultaneity criterion, model the effect of the offsets, or at minimum soften the claim that simultaneity is not important to the result.
  3. [Figure 1 and 2 captions; Section 4] The error bars for both sources are not taken from the original published errors but 'estimated by the size of the symbol in the plot reported by A. Furniss et al. 2015' and by Krawczynski et al. 2004, and the captions state that they are 'must likely overestimated.' The maximum-likelihood fit uses exactly these uncertainties to weight the data and to estimate the parameter errors, p-values, and the intrinsic scatter sigma_d. If the errors are overestimated, the reported parameter uncertainties are correspondingly unreliable, and a fit with more realistic errors could change the conclusions, particularly the significance of any apparent deviations or lack thereof. The authors should obtain the actual published errors or perform a sensitivity analysis to show that the conclusions are robust to the assumed error bars.
minor comments (4)
  1. [Throughout] There are several typographical issues: 'must likely overestimated' should be 'most likely overestimated'; 'out layers' should be 'outliers'; 'el al.' in the reference to M.M. González et al. should be 'et al.'; and 'the data is' should be 'the data are' or 'the data set is.' These should be corrected in a revision.
  2. [Reference list] Reference [9] contains a stray overleaf URL inserted into the citation: 'astro-ph/9707179]https://www.overleaf.com/project/...'. This appears to be an unintended artifact and should be removed.
  3. [Figure 3 and 4] The correlation plots are central to the paper, but the captions do not state whether the axes are the X-ray flux and gamma-ray flux in the same units as the text, nor do they show the fitted line with its 3-sigma_d band. Adding the best-fit lines and the 3-sigma_d envelopes would make the 'no break-down' claim directly visible and easier to assess.
  4. [Section 5] The sentence 'the correlations do not show out layers at high gamma-ray fluxes' is stated as a result, but the preceding sentence acknowledges that the data set is reduced. I recommend rephrasing to make explicit that the absence of outliers is an upper limit based on the current sample, not a definitive conclusion.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the paper reports descriptive maximum-likelihood fits to external multiwavelength data and does not recycle fitted parameters as inputs to a derived result.

full rationale

The paper's derivation chain is a maximum-likelihood linear correlation fit (F_gamma = a F_X + b with intrinsic scatter sigma_d) applied to previously published X-ray and TeV data for Mrk 501 and 1ES 1959+650. The fitted parameters a, b, and sigma_d are outputs of the statistical procedure, not inputs; they are not subsequently used to derive a theoretical quantity or to justify the same correlation. The statistical method is explicitly attributed to D'Agostini 2005, a standard external reference, and the data come from Furniss et al. 2015 and Krawczynski et al. 2004, also external. The paper cites the authors' prior work on Mrk 421 (Gonzalez et al. 2019) only as motivation and as a comparison point, not as evidence for the new fits, so this self-citation is not load-bearing. The one substantive weakness is that for 1ES 1959+650 the correlation was computed only on the low-state subset because the May 16-17 flare and the June 4 orphan flare were excluded, making the later claim of 'no break-down at high gamma-ray fluxes' an overgeneralization from limited coverage. That is a data-selection and external-validity concern, not a circularity: the excluded high-flux points are not used in the fit, and the fit is not defined in terms of the conclusion. No equation in the paper reduces to another equation by construction, and no fitted parameter is renamed as a prediction. Therefore the circularity score is 0.

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

The central claim is an empirical correlation, so the free parameters are the fitted slope, intercept, and intrinsic scatter for each source. The main axioms are the applicability of the linear-with-intrinsic-scatter model and the adequacy of quasi-simultaneous data. No new entities are introduced.

free parameters (6)
  • Mrk 501 correlation slope a = 194 ± 9.20
    Maximum-likelihood slope of the linear correlation F_gamma = a F_X + b, fitted to MAGIC/VERITAS/Swift data.
  • Mrk 501 correlation intercept b = -16.9 ± 1.47
    Intercept of the linear correlation fit for Mrk 501, fitted to the same data.
  • Mrk 501 intrinsic scatter sigma_d = 2.12 ± 0.67
    Fitted intrinsic scatter of unknown nature in the D'Agostini maximum-likelihood method.
  • 1ES 1959+650 correlation slope a = 71.10 ± 6.48
    Maximum-likelihood slope of the linear correlation fit for 1ES 1959+650, using Whipple/RXTE data.
  • 1ES 1959+650 correlation intercept b = 0.36 ± 0.07
    Intercept of the linear correlation fit for 1ES 1959+650.
  • 1ES 1959+650 intrinsic scatter sigma_d = 0.23 ± 0.06
    Fitted intrinsic scatter of unknown nature for the 1ES 1959+650 correlation.
assumptions (3)
  • domain assumption The D'Agostini maximum-likelihood method assumes a linear relation with Gaussian intrinsic scatter and known measurement errors on both axes.
    Used in Section 4 to fit the correlation and to quantify significance.
  • domain assumption Quasi-simultaneous pairing of X-ray and TeV observations is adequate to recover the intrinsic correlation, despite average time offsets of 30 minutes for Mrk 501 and 1.7 hours for 1ES 1959+650.
    The paper states the data lack simultaneity but still uses them as correlated pairs (Sections 2, 3, 5).
  • domain assumption The original data published by Furniss et al. (2015) and Krawczynski et al. (2004) are reliable, and error bars estimated from symbol sizes are conservative.
    The authors digitized the light curves from the original papers and state error bars are likely overestimated.

how reviews work

0 comments
Cite this review

Pith. "Pith review of On the TeV gamma-ray and X-ray correlations exhibited in high-energy peaked BL Lacs: Mrk 501 and 1ES 1959+650." pith.science (2026). https://pith.science/paper/DD2JREUK

@misc{pith2026190900897,
  author       = {Pith},
  title        = {Pith review of: On the TeV gamma-ray and X-ray correlations exhibited in high-energy peaked BL Lacs: Mrk 501 and 1ES 1959+650},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/DD2JREUK}},
  note         = {Machine review of arXiv:1909.00897}
}
read the original abstract

The spectral energy distribution (SED) of blazars is generally understood through the standard one-zone synchrotron self-Compton (SSC) model, where a strong correlation between X-ray and TeV gamma-ray fluxes is expected. Recently, for Mrk 421 the correlation was confirmed and interpreted within the SSC model arguing that a description of the correlation up to large flux values considerably restricts the value of the magnetic field under the assumption of an electron spectral index independent of the flux state. In this work, we extend the analysis done on Mrk 421 for other two high-energy peaked BL Lacs: Mrk 501 and 1ES 1959+650.

Figures

Figures reproduced from arXiv: 1909.00897 by the authors.

Figure 1
Figure 1. Gamma-ray light curve (above) and X-ray light curve (below) of Mrk 501. The data taken was from MAGIC, VERITAS for the gamma-rays and Swift XRT (0.3-3keV) for X-rays. The error bars were estimated by the size of the symbol in the plot reported by A. Furniss, et al.,2015 (see [13]), so they are must likely overestimated. 3. 1ES 1959+650 1ES 1959+650 is a BL Lac type blazar with a reported redshift of z=0.047 [14]. It… view at source ↗
Figure 2
Figure 2. The plot above shows the gamma-ray light curve observed by Whipple at 600 GeV and below the X-ray light curve observed by RXTE in 10 keV of HBL 1ES 1959+650. The error bars were estimated by the size of the symbol in the plot reported by H. Krawczynski, et al., 2004 [14], so they are must likely overestimated. 4. Analysis The maximum likelihood approach is used to estimate the robustness of the correlations and to 3… view at source ↗
Figure 3
Figure 3. Correlation of gamma-rays and X-rays for Mrk 501 with data of MAGIC and VERITAS for the gamma-rays and the Swift telescope for the X-rays. The data lack time simultaneity in an average of 30 minutes. With these data there is no evidence of a break-down at high fluxes, all of the points are within 3σd. Acknowledgements The authors would like to thank DGAPA-UNAM for financial support by grants AG100317, IA102019 and D… view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Correlation of gamma-rays and X-rays for the HBL 1ES 1959+650 with data of Whipple for gamma-rays and of RXTE for X-rays for only a low state. The data is not simultaneous in an average time of 1.7 hours, however the correlation seems not to be affected by it. [7] V.A.…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

16 extracted references · 8 canonical work pages

  1. [1]

    M. H. Ulrich et al., Variability of Active Galactic Nuclei, (1997) ARA&A 35 (1997) 455-502

    work page 1997

  2. [2]

    A. A. Abdo et al., Fermi Large Area Telescope Observations of Markarian 421: The Missing Piece of its Spectral Energy Distribution, (2011) The Astrophysical Journal 736 (2011) 2 [astro-ph/1106.1348]

    arXiv 2011

  3. [3]

    Could a plasma in quasi-thermal equilibrium be associated to the "orphan" TeV flares ?

    N. Fraija, Could a plasma in quasi-thermal equilibrium be associated to the “orphan” TeV flares?, (2015) Astroparticle Physics 71 (2015) 1-20 [astro-ph/1501.04165]

    work page Pith review arXiv 2015

  4. [4]

    N. Fraija et al., Long-term Optical Polarization Variability and Multiwavelength Analysis of Blazar Mrk 421, (2017) The Astrophysical Journal Supplement Series 232 (2017) 7 [astro-ph/1707.05624]

    arXiv 2017

  5. [5]

    M.M. González, el al., Reconcilement of VHE γ-ray/X-ray correlation studies in Mrk 421 and break-down at high fluxes, Monthly Notices of the Royal Astronomical Society 484 (2019) 2944–2953 [arXiv:1901.02454]

    work page Pith review arXiv 2019

  6. [6]

    F. Aharonian, et al., TeV γ-ray light curve and energy spectrum of Mrk 421 during its 2001 flare as measured with HEGRA CT1, Astronomy & Astrophysics 410 (2003) 813-821 5 Correlation between γ-ray and X-ray emissions in HBLs. Rodrigo Sacahui Figure 4: Correlation of gamma-rays and X-rays for the HBL 1ES 1959+650 with data of Whipple for gamma-rays and of R...

    work page 2003

  7. [7]

    Observation of Markarian 421 in TeV gamma rays over a 14-year time span

    V .A. Acciari, et al., Observation of Markarian 421 in TeV gamma-rays over a 14-year time span, Astroparticle physics 54 (2014) 1–10 [arXiv:1310.8150]

    work page Pith review arXiv 2014

  8. [8]

    Quinn, et al., Detection of gamma-rays with E> 300 GeV from Markarian 501, The Astrophysical Journal Letters 456 (1996) L83

    J. Quinn, et al., Detection of gamma-rays with E> 300 GeV from Markarian 501, The Astrophysical Journal Letters 456 (1996) L83

    work page 1996

Show all 16 references
  1. [9]

    M. Catanese, et al., Multiwavelength observations of a flare from Markarian 501, The Astrophysical Journal Letters 487 (1997) L143 [astro-ph/9707179]https://www.overleaf.com/project/5cf9527b60107d05d0fb2c36

    1997 arXiv

  2. [10]

    Kataoka, et al., High-energy emission from the TeV blazar Markarian 501 during multiwavelength observations in 1996, The Astrophysical Journal 514 (1999) 138 [astro-ph/9811014]

    J. Kataoka, et al., High-energy emission from the TeV blazar Markarian 501 during multiwavelength observations in 1996, The Astrophysical Journal 514 (1999) 138 [astro-ph/9811014]

    1999 arXiv

  3. [11]

    Petry, et al., Multiwavelength observations of Markarian 501 during the 1997 high state, The Astrophysical Journal 536 (2000) 742 [astro-ph/0002255]

    D. Petry, et al., Multiwavelength observations of Markarian 501 during the 1997 high state, The Astrophysical Journal 536 (2000) 742 [astro-ph/0002255]

    2000 arXiv

  4. [12]

    A. A. Abdo, et al., Insights into the High-energy γ-ray Emission of Markarian 501 from Extensive Multifrequency Observations in the Fermi Era, The Astrophysical Journal 727 (2011) 129 [arXiv:1011.5260]

    2011 arXiv

  5. [13]

    Furniss, et al., First NuSTAR observations of mrk 501 within a radio to tev multi-instrument campaign, The Astrophysical Journal 812 (2015) 65 [arXiv:1509.04936]

    A. Furniss, et al., First NuSTAR observations of mrk 501 within a radio to tev multi-instrument campaign, The Astrophysical Journal 812 (2015) 65 [arXiv:1509.04936]

    2015 arXiv

  6. [14]

    Krawczynski, et al., Multiwavelength Observations of Strong Flares from the TeV Blazar 1ES 1959+650, The Astrophysical Journal 601 (2004) 151 [astro-ph/0310158]

    H. Krawczynski, et al., Multiwavelength Observations of Strong Flares from the TeV Blazar 1ES 1959+650, The Astrophysical Journal 601 (2004) 151 [astro-ph/0310158]

    2004 arXiv

  7. [15]

    Jahoda, et al., Calibration of the Rossi X-Ray Timing Explorer Proportional Counter Array, The Astrophysical Journal Supplement Series 163 (2006) 401-423 [astro-ph/0511531]

    K. Jahoda, et al., Calibration of the Rossi X-Ray Timing Explorer Proportional Counter Array, The Astrophysical Journal Supplement Series 163 (2006) 401-423 [astro-ph/0511531]. 6 Correlation between γ-ray and X-ray emissions in HBLs. Rodrigo Sacahui

    2006 arXiv

  8. [16]

    D’Agostini, Fits, and especially linear fits, with errors on both axes, extra variance of the data points and other complications, (2005) [arXiv preprint:physics/0511182 ]

    G. D’Agostini, Fits, and especially linear fits, with errors on both axes, extra variance of the data points and other complications, (2005) [arXiv preprint:physics/0511182 ]. 7

    2005 arXiv

Pith tools

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