Pith. sign in

REVIEW 4 major objections 4 minor 42 references

Do multifrequency polarimetric observations of BL Lac rule out a hadronic origin for its X-ray emission?

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

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

Pith's one-line read A hadronic jet model with proton-synchrotron X-rays can reproduce BL Lac's high optical and low X-ray polarization, so the IXPE measurements do not rule out a hadronic origin.

desk verdict A serious but over-claimed attempt to revive a hadronic explanation for BL Lac's X-ray polarization; the physical mechanism is real, but the quantitative match depends on a cooling-volume inconsistency and post-hoc parameters. read the letter →

arxiv 2505.19784 v2 pith:OGJTFAYG submitted 2025-05-26 astro-ph.HE

F. Tavecchio , F. Bolis , E. Sobacchi , S. Boula , A. Sciaccaluga This is my paper · ORCID
classification astro-ph.HE
keywords BLLacblazarjetshadronicmodelsprotonsynchrotronX-raypolarizationIXPEmagneticallydominateddegree
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

The paper argues that the multifrequency polarization data for BL Lac do not force a leptonic interpretation of its X-ray emission. It constructs a hadronic model in which relativistic electrons produce the optical synchrotron component while protons, radiating through proton synchrotron, produce the X-rays and gamma-rays. In this model the X-ray polarization is naturally small, mostly because the proton spectrum is hard (slope $p_p = 2.3$) and because slow-cooling protons fill a much larger jet volume than fast-cooling electrons. The authors conclude that a proton-synchrotron origin for BL Lac's X-rays remains viable, with jet powers below the Eddington luminosity.

What carries the argument

The central mechanism is the polarization of synchrotron radiation computed from the MHD model of magnetically dominated, stationary, axisymmetric outflows of Lyubarsky (2009), in which the jet is parabolic ($R_0 \propto z_0^q$) and carries a helical magnetic field. For such fields, the polarization degree depends strongly on the slope of the emitting particle energy distribution, more so than in the textbook uniform-field case, and also on the size and shape of the emission region. The optical versus X-ray contrast comes from two different cooling regimes: electrons emitting in the optical band are fast-cooling and radiate from an extremely thin shell ($\Delta z \to 0$), while protons are slow-cooling and are advected outward until adiabatic losses set a scale $\Delta z \simeq 3z_0/(2q)$, a few times $z_0$. The hard proton distribution and the larger proton-emitting volume jointly push $\Pi_X$ below the IXPE limit, while the steeper electron distribution keeps $\Pi_O$ high.

What would settle it

Measure the X-ray polarization of BL Lac during a similar state to the IXPE campaign and see whether it stays below about $7\%$; a value above roughly $15\%$ would contradict the model. Independently, determine the proton injection slope from the gamma-ray SED and the shape of the high-energy hump; a slope steeper than $p_p \approx 2.3$ would remove the parameter choice that keeps $\Pi_X$ low.

Watch

Extended reading notes

Core claim

On the paper's own terms, the low X-ray polarization upper limit of about $7\%$ measured by IXPE does not rule out hadronic X-ray production. Using a magnetically dominated, stationary, axisymmetric jet with a helical magnetic field, the authors compute synchrotron polarization from two particle populations: electrons with spectral slope $p_e = 4.6$ radiating from a very thin shell at the acceleration site, and protons with slope $p_p = 2.3$ radiating from a region extending to $\Delta z \simeq 3\text{--}4\,z_0$ set by adiabatic cooling. The calculation yields optical polarization $\Pi_O \sim 25\text{--}45\%$ and X-ray polarization $\Pi_X \lesssim 7\%$, matching the observed averages, with the electric vector position angle aligned with the jet axis in both bands as observed. The earlier conclusion that hadronic models are ruled out is attributed to single-zone models with uniform fields and radiative-only cooling, both of which this paper drops.

Load-bearing premise

The match to the IXPE upper limit depends on adopting a hard proton spectrum with slope 2.3 and letting the protons radiate from a region roughly three to four times the distance from the black hole; if the real proton spectrum were steeper or the proton emission region smaller, the predicted X-ray polarization would rise above 7%.

Editorial extensions

If this is right

  • If this model is right, the IXPE polarization measurement no longer uniquely favors leptonic X-rays; proton-synchrotron X-rays remain observationally viable for BL Lac.
  • The standard objections to hadronic models, namely long variability timescales and too-high predicted X-ray polarization, are resolved once adiabatic cooling and an extended proton emission region are included.
  • Because the proton cooling length is essentially energy-independent below the SED peak, the MeV band should show a similarly low degree of polarization, while above the roughly 100 MeV peak the polarization should rise, reaching about $15\text{--}20\%$ in the constant-density case.
  • The required proton energy flux, about $3\times 10^{45}\,\mathrm{erg\,s^{-1}}$, keeps the total jet power below the Eddington luminosity of the estimated black hole mass, so the hadronic scenario is not energetically extreme.
  • The model predicts a very low neutrino flux, well below current detector sensitivities, because photomeson production is negligible in the adopted parameter regime.

Reading between the lines

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

  • The parameter combination that reproduces the observations ($p_p = 2.3$ and $\Delta z \simeq 3\text{--}4\,z_0$) is selected after the IXPE results were known; an independent test would be to derive the proton injection spectrum from the gamma-ray SED and check whether it is actually as hard as 2.3.
  • A testable extension of the jet geometry is that an optical flaring episode produced by a small moving region on a helical path should show a large increase in optical polarization with little change in X-ray polarization; simultaneous IXPE and optical monitoring during a flare would distinguish this picture from single-zone models.
  • If a future MeV polarimeter measures the band below the high-energy peak and finds polarization significantly above the X-ray value, that would challenge the model's prediction that adiabatic cooling keeps the polarization uniformly low across the whole high-energy hump.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 4 minor

Summary. The paper revisits IXPE multifrequency polarimetry of BL Lac and argues that the observed low X-ray polarization (ΠX ≲ 7%) with high optical polarization (ΠO ~ 25–30%) does not rule out a hadronic, proton-synchrotron origin of the X-ray emission. The authors construct a stationary, axisymmetric, magnetically dominated MHD jet with a helical magnetic field, assume that electrons emit optical synchrotron radiation from a thin shell while protons emit X-rays from an extended region, and compute Stokes parameters. They claim that for jet shape parameter q = 0.6–0.7 and proton emission length Δz ≳ 3 z0, the model reproduces the observed polarization in both a constant-density and a constant-magnetization prescription. The paper concludes that a hadronic model remains viable.

Significance. If the calculation were self-consistent, this paper would be an important counterpoint to the widely cited conclusion that IXPE observations of BL Lac favor a leptonic origin of its X-ray emission. The MHD field setup, the treatment of synchrotron polarization, and the resulting Table 1 are transparent and allow the reader to check the parameter dependence; the appendices provide useful derivations of the adiabatic cooling length, jet power, and Stokes integrals. However, the central claim is presently carried by an internally inconsistent choice of the proton emission volume and by a misstatement of the results in Table 1. After correcting these issues, only a narrow parameter subset survives, so the paper is a useful contribution in need of substantial revision rather than a definitive demonstration that hadronic models are not ruled out.

major comments (4)
  1. [Sect. 2.1, Appendix B, Table 1] The proton emission lengths Δz = 3–4 z0 adopted to obtain the low X-ray polarization are longer than the adiabatic cooling length derived in the paper. Appendix B gives l_adiab = 3z/(2q), which is 2.5 z0 for q = 0.6 and 2.1 z0 for q = 0.7, and Sect. 2.1 identifies the proton cooling length with this scale; Fig. 1's caption even states Δz ∼ z0. At the physically allowed Δz = 2 z0, the constant-density entries in Table 1 give ΠX = 0.164 (q = 0.6) and 0.101 (q = 0.7), both above the IXPE upper limit of 0.07; among the q = 0.6–0.7 cases highlighted in the text, only the constant-magnetization q = 0.7 entry passes (ΠX = 0.032). The claim that the low X-ray polarization is naturally produced by the larger proton cooling volume is therefore not supported for the parameter range emphasized in the text.
  2. [Sect. 2.3, Table 1] The statement that the observed polarization can be reproduced 'in both cases ... for q = 0.6−0.7 and Δz ≳ 3' is contradicted by Table 1. For q = 0.6 with constant density, Δz = 4 z0 gives ΠX = 0.118, which exceeds the IXPE upper limit; for q = 0.6 with constant magnetization, Δz = 2 z0 gives ΠX = 0.112, also above 0.07. Only a subset of the quoted parameter combinations actually satisfies ΠX < 0.07. The text, Table 1, and Fig. 2 must be brought into agreement, and the claimed parameter range must be corrected.
  3. [Appendix D] The Stokes integrals in Appendix D integrate a fixed power-law proton distribution over the entire volume z0 < z < z0 + Δz, with no depletion or spectral evolution of the protons due to adiabatic losses. Since the adiabatic cooling length is only about 2 z0, the outer part of a 4 z0 emission region should have a lower normalization and a high-energy cutoff; its contribution to the computed low polarization is therefore spurious. The calculation should either restrict Δz to the cooling length or include a distance-dependent proton distribution that accounts for adiabatic cooling.
  4. [Sect. 2.3] The proton spectral slope pp = 2.3 is adopted from the observed X-ray spectrum (Agudo et al. 2025), and the emission region size Δz is chosen after the fact to reach ΠX < 7%. Consequently, the low X-ray polarization is a consistency check rather than an independent prediction of the hadronic model. The paper would be strengthened by a sensitivity study showing how ΠX varies with pp and by an explicit statement that pp and Δz are tuned to the observations.
minor comments (4)
  1. [Fig. 1 caption] The caption says protons are advected 'up to a distance Δz ∼ z0', which is inconsistent with the values Δz = 3–4 z0 used in the calculations; please reconcile the sketch and caption with the quantitative model.
  2. [Table 1] The header is confusing: the first numerical column (labeled Δz → 0) contains ΠO, while the following columns contain ΠX, but this is not stated explicitly. Please label the columns as ΠO and ΠX to avoid ambiguity.
  3. [Appendix B] The function A(q) in Eq. (B.5) is described as 'a factor of order unity' but no explicit expression is given; please provide the expression or a precise citation to the derivation in Bolis et al. (2024b).
  4. [Sect. 2.1] The sentence 'lp,cool = Δz ≃ 3z0/2q' is at odds with both Fig. 1's caption (Δz ∼ z0) and the later use of Δz = 3–4 z0; this should be clarified in the text.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the paper explicitly presents a post-hoc consistency model rather than a prediction, and the polarization calculation is an independent Stokes integration with slopes taken from the observed SED.

full rationale

The paper's central claim is explicitly framed as a consistency demonstration: 'these observations can also be explained by a hadronic model' (Sect. 3), not as an a priori prediction. The X-ray and optical polarization degrees are computed by integrating the Stokes parameters over an MHD jet model (Appendix D), with a nontrivial dependence on the particle slope, geometry, and emission-region width; this is not equivalent to the input parameters by construction. The proton slope pp=2.3 and electron slope pe=4.6 are adopted from the observed multifrequency SED (Sect. 2.3, citing Agudo et al. 2025), and q and Δz are scanned over a physically motivated range, so this is parameter fitting rather than definitional circularity. Self-citations to Bolis et al. (2024a,b) are load-bearing for the polarization formalism, but the formalism is stated with explicit assumptions and reproduced in Appendix D, so under the rules it counts as independent support rather than circular self-citation. The principal caveat is internal consistency rather than circularity: the adiabatic cooling length derived in Appendix B, l_adiab=3z/(2q), is about 2.1-2.5 z0 for q=0.6-0.7, while the successful entries in Table 1 require Δz≥3-4 z0; at Δz=2z0, the constant-density cases give ΠX=0.164 (q=0.6) and 0.101 (q=0.7), both above the IXPE upper limit of 0.07. This weakens the claim that the low X-ray polarization is 'naturally' produced, but it is a physical-consistency problem, not a reduction of the output to the input by definition.

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

The model introduces no new particle or field. The burden is carried by adopted field geometry, isotropic power-law distributions, and the untested choices of q, Delta z, B', delta, and spectral slopes. The polarization match is obtained by varying q and Delta z after the IXPE limit was known.

free parameters (6)
  • jet shape parameter q = 0.6-0.8
    This parameter determines the opening angle and field structure; the match to the observed optical and X-ray polarization requires q roughly 0.6-0.7. It is not independently measured for BL Lac.
  • proton emission region length Delta z/z0 = approximately 3-4
    The low X-ray polarization is reached only for Delta z at least about 3 z0; this is motivated by the adiabatic cooling length but is not measured. Smaller values overproduce X-ray polarization.
  • magnetic field B' = 10 G
    Fiducial value chosen so that protons producing the 100 MeV peak have plausible energies and SSC is suppressed; not fitted to the polarization but adopted as an input.
  • Doppler factor delta and Lorentz factor Gamma = delta=10, Gamma=10, theta_obs=0.1 rad
    Fiducial blazar values; Gamma=10 is set through the condition Gamma(R0)=10 in the MHD model and the viewing angle is assumed.
  • proton spectral slope pp = 2.3
    Adopted from the observed X-ray spectral slope (Agudo et al. 2025); the computed X-ray polarization is very sensitive to this slope, so the input is load-bearing.
  • electron spectral slope pe = 4.6
    Adopted from the observed optical spectrum; the high optical polarization depends on this steep slope.
assumptions (6)
  • domain assumption Magnetically dominated stationary axisymmetric outflow with external pressure profile P_ext proportional to z^-kappa
    The magnetic field structure and jet shape are taken from Lyubarsky 2009; invoked in Sect. 2.2.
  • domain assumption The jet bulk flow velocity equals the drift velocity E x B / B^2
    Stated in Sect. 2.2 as appropriate for Poynting-dominated outflows; not independently confirmed for BL Lac.
  • domain assumption Particle distributions are isotropic power laws in the fluid frame
    Assumed in Sect. 2.2 for both electrons and protons; this is the standard synchrotron modeling assumption.
  • ad hoc to paper Electrons emit in a thin shell while protons fill a region up to Delta z about 4 z0
    Cooling length arguments place electrons in a thin slice and protons in a larger volume, but the size of the proton region is tuned to match the observed polarization.
  • domain assumption SSC and proton-photon losses are negligible relative to proton synchrotron
    Appendix A and footnote 2 justify neglecting these channels; this assumption makes proton synchrotron the sole high-energy emission channel.
  • standard math Standard synchrotron Stokes parameter formulas for a stationary jet apply
    The integrals in Appendix D follow Lyutikov et al. 2003 and Del Zanna et al. 2006; used without modification.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Do multifrequency polarimetric observations of BL Lac rule out a hadronic origin for its X-ray emission?." pith.science (2026). https://pith.science/paper/OGJTFAYG

@misc{pith2026250519784,
  author       = {Pith},
  title        = {Pith review of: Do multifrequency polarimetric observations of BL Lac rule out a hadronic origin for its X-ray emission?},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OGJTFAYG}},
  note         = {Machine review of arXiv:2505.19784}
}
abstract

Recent multifrequency polarimetric observations of the eponymous blazar BL Lac reveal an extremely large degree of polarization in the optical band (average of $25\%$, reaching $45\%$), together with a small ($\lesssim 7\%$) degree of polarization in the X-ray band. This has been interpreted as evidence that the X-rays are produced through inverse Compton emission by relativistic electrons, thus ruling out alternative models based on hadronic processes. Here we revisit the observational evidence, interpreting it in a framework where the observed radiation is entirely produced through synchrotron emission. Electrons produce the radio-to-optical component and protons produce the X-rays and the gamma-rays. We determine the jet magnetic fields from an MHD model of magnetically dominated stationary axisymmetric outflows, and show that the X-ray emission from the protons is naturally less polarized than the optical emission from the electrons. The model parameters required to reproduce the multifrequency polarimetric observations are fully compatible with blazar jets.

Figures

Figures reproduced from arXiv: 2505.19784 by the authors.

Figure 1
Figure 1. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. X-ray polarization degree, ΠX (top panels) and chromaticity, ΠO/ΠX (bottom panels) as a function of the width of the proton emission region, ∆z. In the left panels, the proper number density is assumed to be constant, whereas in the right panels, the magnetization is assumed to be constant. The X-ray polarization degree is ΠX = Πpp=2.3, and the optical polarization degree is ΠO = Πpe=4.6. The viewing angle is assume… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

42 extracted references · 33 canonical work pages

  1. [1]

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

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

    work page 2010

  2. [2]

    2025, arXiv e-prints, arXiv:2505.01832

    Agudo , I., Liodakis , I., Otero-Santos , J., et al. 2025, arXiv e-prints, arXiv:2505.01832

    arXiv 2025

  3. [3]

    Aharonian , F. A. 2000, , 5, 377

    work page 2000

  4. [4]

    & Petruk , O

    Bandiera , R. & Petruk , O. 2016, , 459, 178

    work page 2016

  5. [5]

    & Petruk , O

    Bandiera , R. & Petruk , O. 2024, , 689, A137

    work page 2024

  6. [6]

    2019, , 57, 467

    Blandford , R., Meier , D., & Readhead , A. 2019, , 57, 467

    2019

  7. [7]

    2024 a , , 690, A14

    Bolis , F., Sobacchi , E., & Tavecchio , F. 2024 a , , 690, A14

    work page 2024

  8. [8]

    2024 b , , 110, 123032

    Bolis , F., Sobacchi , E., & Tavecchio , F. 2024 b , , 110, 123032

    work page 2024

Show all 42 references
  1. [9]

    & Saggion , A

    Bonometto , S. & Saggion , A. 1973, , 23, 9

    1973

  2. [10]

    2019, Galaxies, 7, 20

    B \"o ttcher , M. 2019, Galaxies, 7, 20

    2019

  3. [11]

    M., & Buttiglione , S

    Capetti , A., Raiteri , C. M., & Buttiglione , S. 2010, , 516, A59

    2010

  4. [12]

    2006, , 453, 621

    Del Zanna , L., Volpi , D., Amato , E., & Bucciantini , N. 2006, , 453, 621

    2006

  5. [13]

    L., Ehlert , S

    Di Gesu , L., Marshall , H. L., Ehlert , S. R., et al. 2023, Nature Astronomy, 7, 1245

    2023

  6. [14]

    1998, , 299, 433

    Fossati , G., Maraschi , L., Celotti , A., Comastri , A., & Ghisellini , G. 1998, , 299, 433

    1998

  7. [15]

    1998, , 301, 451

    Ghisellini , G., Celotti , A., Fossati , G., Maraschi , L., & Comastri , A. 1998, , 301, 451

    1998

  8. [16]

    Hillas , A. M. 1984, , 22, 425

    1984

  9. [17]

    2012, , 744, 30

    Krawczynski , H. 2012, , 744, 30

    2012

  10. [18]

    & Petropoulou , M

    Liodakis , I. & Petropoulou , M. 2020, , 893, L20

    2020

  11. [19]

    2025, arXiv e-prints, arXiv:2505.13603

    Liodakis , I., Zhang , H., Boula , S., et al. 2025, arXiv e-prints, arXiv:2505.13603

    2025 arXiv

  12. [20]

    2009, , 698, 1570

    Lyubarsky , Y. 2009, , 698, 1570

    2009

  13. [21]

    I., & Blandford , R

    Lyutikov , M., Pariev , V. I., & Blandford , R. D. 2003, , 597, 998

    2003

  14. [22]

    L., Liodakis , I., Marscher , A

    Marshall , H. L., Liodakis , I., Marscher , A. P., et al. 2024, , 972, 74

    2024

  15. [23]

    2023, , 942, L10

    Middei , R., Liodakis , I., Perri , M., et al. 2023, , 942, L10

    2023

  16. [24]

    D., Takami , H., & Migliori , G

    Murase , K., Dermer , C. D., Takami , H., & Migliori , G. 2012, , 749, 63

    2012

  17. [25]

    E., Middei , R., et al

    Pacciani , L., Kim , D. E., Middei , R., et al. 2025, , 983, 78

    2025

  18. [26]

    L., Liodakis , I., & Romani , R

    Peirson , A. L., Liodakis , I., & Romani , R. W. 2022, , 931, 59

    2022

  19. [27]

    L., Negro , M., Liodakis , I., et al

    Peirson , A. L., Negro , M., Liodakis , I., et al. 2023, , 948, L25

    2023

  20. [28]

    1994, , 92, 607

    Poutanen , J. 1994, , 92, 607

    1994

  21. [29]

    E., Boettcher , M., Markoff , S., & Tavecchio , F

    Romero , G. E., Boettcher , M., Markoff , S., & Tavecchio , F. 2017, , 207, 5

    2017

  22. [30]

    C., & Rees , M

    Sikora , M., Begelman , M. C., & Rees , M. J. 1994, , 421, 153

    1994

  23. [31]

    & Zech , A

    Sol , H. & Zech , A. 2022, Galaxies, 10, 105

    2022

  24. [32]

    2003, , 403, 83

    Tavecchio , F., Ghisellini , G., & Celotti , A. 2003, , 403, 83

    2003

  25. [33]

    & COSI Collaboration

    Tomsick , J. & COSI Collaboration . 2022, in 37th International Cosmic Ray Conference, 652

    2022

  26. [34]

    2004, , 600, 324

    Vlahakis , N. 2004, , 600, 324

    2004

  27. [35]

    Webb , G. M. 1989, , 340, 1112

    1989

  28. [36]

    C., Soffitta , P., Baldini , L., et al

    Weisskopf , M. C., Soffitta , P., Baldini , L., et al. 2022, Journal of Astronomical Telescopes, Instruments, and Systems, 8, 026002

    2022

  29. [37]

    Zdziarski , A. A. & Bottcher , M. 2015, , 450, L21

    2015

  30. [38]

    2017, , 602, A25

    Zech , A., Cerruti , M., & Mazin , D. 2017, , 602, A25

    2017

  31. [39]

    2017, Galaxies, 5, 32

    Zhang , H. 2017, Galaxies, 5, 32

    2017

  32. [40]

    & B \"o ttcher , M

    Zhang , H. & B \"o ttcher , M. 2013, , 774, 18

    2013

  33. [41]

    , " * 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.sent...

  34. [42]

    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 gl...

Pith tools

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