Discriminating blazar emission models with high-energy polarimetry: Multi-band predictions and detectability
Pith reviewed 2026-06-29 01:38 UTC · model grok-4.3
The pith
Simultaneous multi-wavelength polarimetry distinguishes between leptonic, hadronic, and hybrid emission models in blazars by differences in polarization degree.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper establishes that for a statistically complete sample of gamma-ray blazars, the expected polarization degrees under leptonic, hadronic, and hybrid scenarios differ enough that simultaneous observations across soft X-ray to gamma-ray energies can discriminate between the models, after accounting for instrument minimum detectable polarization and detection probabilities in blind surveys. It derives specific sensitivity requirements for future missions to increase the number of detectable sources.
What carries the argument
Polarimetric signatures modeled for leptonic, hadronic, and hybrid emission scenarios using Bjet_MCMC spectral energy distribution fitting on a RoboPol sample of blazars.
Load-bearing premise
The polarization signatures predicted for the three emission scenarios remain distinguishable after propagation effects, instrument responses, and uncertainties from spectral energy distribution fits are taken into account.
What would settle it
A set of simultaneous multi-band polarization measurements on several blazars that show polarization degrees inconsistent with the differences predicted by the three models, or where all models predict overlapping values within measurement errors.
Figures
read the original abstract
Polarimetric properties of blazars provide key constraints on the acceleration mechanisms powering their relativistic jets, the high-energy emission processes involved, and the composition of the jet itself. We present a multi-band polarimetric study spanning from soft X-rays to gamma-ray energies, considering several bands for current and future missions (0.275 keV, 2-8 keV, 0.5-10 keV, 6-35 keV, 0.2-5 MeV, and 1-100 GeV). Our sample is drawn from the RoboPol monitoring programme, a statistically complete gamma-ray sample including low-, intermediate-, and high-synchrotron-peaked blazars. Using spectral energy distribution fitting performed with the Bjet_MCMC tool, we give predictions on the flux expected for each source in the selected energy bands. We model the polarimetric signatures under three competing emission scenarios: leptonic, hadronic, and hybrid. The detectability of each source was assessed by accounting for the instruments' minimum detectable polarisation (MDP) and by computing the probability of the detection in a blind survey. We show that simultaneous multi-wavelength observations can effectively discriminate between competing emission models due to the difference in the expected polarisation degree. Finally, we derive sensitivity requirements for future gamma-ray polarimetric missions aimed at increasing the number of detectable sources.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper uses Bjet_MCMC SED fits to a RoboPol gamma-ray blazar sample to predict multi-band fluxes (0.275 keV to 1-100 GeV) and then models the expected polarization degree under leptonic, hadronic, and hybrid emission scenarios. It claims that the resulting differences in polarization degree are large enough, relative to instrument MDP values, that simultaneous multi-wavelength polarimetry can discriminate among the three models, and it derives sensitivity requirements for future gamma-ray polarimeters to increase the number of detectable sources.
Significance. If the polarization differences survive propagation of SED-fit uncertainties and instrument effects, the work supplies concrete, observationally actionable predictions that link existing multi-wavelength monitoring to high-energy polarimetry as a model discriminator for blazar jets. The use of a statistically complete sample and explicit MDP-based detectability calculations are strengths that make the forecasts directly usable for proposal planning.
major comments (2)
- [polarization modeling and detectability sections] The central discrimination claim (abstract and final section) requires that polarization-degree separations between the three scenarios exceed both MDP and the spread induced by Bjet_MCMC posterior uncertainties. The manuscript computes polarization only at the best-fit SED parameters rather than drawing from the MCMC chains or marginalizing over magnetic-field and particle-distribution parameters; this leaves the claimed distinguishability untested against realistic uncertainties.
- [results table of polarization predictions] Table of predicted polarization degrees and detection probabilities (the table that lists per-source values for the three models) reports point estimates without error bars propagated from the SED fits. Because the separation between leptonic/hadronic/hybrid values is the load-bearing quantity for the discrimination argument, the absence of these uncertainties directly affects whether the differences remain statistically significant once fit errors are folded in.
minor comments (2)
- Notation for the three emission scenarios is introduced without a compact reference table; adding one would help readers track which model corresponds to which polarization curve in the figures.
- The MDP values adopted for each instrument band are stated without citing the exact references or formulas used to compute them; a short appendix or footnote with the MDP expressions would improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for their careful and constructive review. The major comments correctly identify that our polarization predictions rely on best-fit SED parameters and that the table lacks propagated uncertainties. We address each point below and will incorporate revisions to strengthen the discrimination analysis.
read point-by-point responses
-
Referee: [polarization modeling and detectability sections] The central discrimination claim (abstract and final section) requires that polarization-degree separations between the three scenarios exceed both MDP and the spread induced by Bjet_MCMC posterior uncertainties. The manuscript computes polarization only at the best-fit SED parameters rather than drawing from the MCMC chains or marginalizing over magnetic-field and particle-distribution parameters; this leaves the claimed distinguishability untested against realistic uncertainties.
Authors: We agree that the current analysis uses best-fit parameters and does not fully propagate posterior uncertainties from Bjet_MCMC. This choice was made for computational practicality across the full sample. To address the concern, we will add a new subsection that samples from the MCMC posteriors for a representative subset of sources (e.g., 10-15 objects spanning the synchrotron-peak classes). Polarization degrees will be recomputed for each draw under the three emission models, yielding uncertainty ranges. These ranges will be compared directly to the inter-model separations and MDP values to test whether the discrimination remains viable. The results and any necessary qualifications will be reported in the revised text. revision: yes
-
Referee: [results table of polarization predictions] Table of predicted polarization degrees and detection probabilities (the table that lists per-source values for the three models) reports point estimates without error bars propagated from the SED fits. Because the separation between leptonic/hadronic/hybrid values is the load-bearing quantity for the discrimination argument, the absence of these uncertainties directly affects whether the differences remain statistically significant once fit errors are folded in.
Authors: We concur that the absence of uncertainties in the table limits the strength of the discrimination argument. In the revision we will augment the table with uncertainty estimates on the polarization degrees. For the subset of sources where full posterior sampling is performed, we will report the standard deviation across samples. For the remaining sources we will derive representative uncertainty ranges by varying key parameters (magnetic field strength, particle indices) within their posterior widths and recomputing polarization. These additions will allow direct evaluation of whether model separations exceed the combined uncertainties and MDP. revision: yes
Circularity Check
No significant circularity; polarization modeling is independent of SED fits
full rationale
The derivation uses Bjet_MCMC SED fits solely to predict fluxes in selected bands, then applies separate modeling of polarimetric signatures under leptonic, hadronic, and hybrid scenarios. No quoted equations or steps in the abstract or description show polarization degrees reducing by construction to fitted parameters, self-citations, or ansatzes imported from prior work. The central claim of distinguishability rests on differences in modeled polarization degrees after accounting for MDP and survey probabilities, which are presented as independent outputs. This is the most common honest finding for papers whose core computation chain remains externally falsifiable.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
2025, ApJ, 985, L15
Agudo, I., Liodakis, I., Otero-Santos, J., et al. 2025, ApJ, 985, L15
2025
-
[2]
2020, ApJ, 892, 105
Ajello, M., Angioni, R., Axelsson, M., et al. 2020, ApJ, 892, 105
2020
-
[3]
2022, ApJS, 263, 24
Ajello, M., Baldini, L., Ballet, J., et al. 2022, ApJS, 263, 24
2022
-
[4]
Angel, J. R. P. & Stockman, H. S. 1980, ARA&A, 18, 321
1980
-
[5]
Fermi Large Area Telescope Fourth Source Catalog Data Release 4 (4FGL-DR4)
Ballet, J., Bruel, P., Burnett, T., et al. 2023, arXiv preprint arXiv:2307.12546
work page internal anchor Pith review Pith/arXiv arXiv 2023
-
[6]
2019, ARA&A, 57, 467
Blandford, R., Meier, D., & Readhead, A. 2019, ARA&A, 57, 467
2019
-
[7]
& Pavlidou, V
Blinov, D. & Pavlidou, V . 2019, Galaxies, 7, 46
2019
-
[8]
2018, MNRAS, 474, 1296
Blinov, D., Pavlidou, V ., Papadakis, I., et al. 2018, MNRAS, 474, 1296
2018
-
[9]
E., et al
Blinov, D., Pavlidou, V ., Papadakis, I. E., et al. 2016, MNRAS, 457, 2252
2016
-
[10]
Modeling the Spectral Energy Distributions and Variability of Blazars
Boettcher, M. 2012, arXiv e-prints, arXiv:1205.0539
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[11]
& Saggion, A
Bonometto, S. & Saggion, A. 1973, A&A, 23, 9 Böttcher, M., Reimer, A., Sweeney, K., & Prakash, A. 2013, ApJ, 768, 54
1973
-
[12]
2025, A&A, 703, A19
Capecchiacci, S., Liodakis, I., Middei, R., et al. 2025, A&A, 703, A19
2025
-
[13]
A., et al
Caputo, R., Ajello, M., Kierans, C. A., et al. 2022, Journal of Astronomical Tele- scopes, Instruments, and Systems, 8, 044003
2022
-
[14]
J., Liodakis, I., Middei, R., et al
Chen, C.-T. J., Liodakis, I., Middei, R., et al. 2024, ApJ, 974, 50 de Angelis, A., Tatischeff, V ., Grenier, I. A., et al. 2018, Journal of High Energy Astrophysics, 19, 1 de Jaeger, T., Shappee, B. J., Kochanek, C. S., et al. 2023, MNRAS, 519, 6349 Di Gesu, L., Donnarumma, I., Tavecchio, F., et al. 2022, ApJ, 938, L7 Di Gesu, L., Ferrazzoli, R., Donnaru...
2024
-
[15]
R., Ferrazzoli, R., Marinucci, A., et al
Ehlert, S. R., Ferrazzoli, R., Marinucci, A., et al. 2022, ApJ, 935, 116
2022
-
[16]
P., et al
Errando, M., Liodakis, I., Marscher, A. P., et al. 2024, ApJ, 963, 5
2024
-
[17]
A., & Youngquist, A
Hervet, O., Johnson, C. A., & Youngquist, A. 2024, ApJ, 962, 140
2024
-
[18]
& Lindfors, E
Hovatta, T. & Lindfors, E. 2019, New A Rev., 87, 101541 IceCube Collaboration, Aartsen, M. G., Ackermann, M., et al. 2018, Science, 361, eaat1378
2019
-
[19]
E., Di Gesu, L., Liodakis, I., et al
Kim, D. E., Di Gesu, L., Liodakis, I., et al. 2024, A&A, 681, A12
2024
-
[20]
M., Hovatta, T., Lindfors, E., et al
Kouch, P. M., Hovatta, T., Lindfors, E., et al. 2026, A&A, 708, A383
2026
-
[21]
M., Liodakis, I., Middei, R., et al
Kouch, P. M., Liodakis, I., Middei, R., et al. 2024, A&A, 689, A119
2024
-
[22]
2012, ApJ, 744, 30
Krawczynski, H. 2012, ApJ, 744, 30
2012
-
[23]
A., Otero-Santos, J., Negro, M., et al
Latiolais, G. A., Otero-Santos, J., Negro, M., et al. 2026, The Astrophysical Jour- nal, 1004, 80
2026
-
[24]
P., Agudo, I., et al
Liodakis, I., Marscher, A. P., Agudo, I., et al. 2022, Nature, 611, 677
2022
-
[25]
W., Filippenko, A
Liodakis, I., Romani, R. W., Filippenko, A. V ., et al. 2018, MNRAS, 480, 5517
2018
-
[26]
W., Filippenko, A
Liodakis, I., Romani, R. W., Filippenko, A. V ., Kocevski, D., & Zheng, W. 2019, ApJ, 880, 32
2019
-
[27]
P., Liodakis, I., Saade, M
Maksym, W. P., Liodakis, I., Saade, M. L., et al. 2025, ApJ, 986, 230
2025
-
[28]
& Biermann, P
Mannheim, K. & Biermann, P. L. 1992, A&A, 253, L21
1992
-
[29]
2022, in 44th COSPAR Scientific As- sembly
Marshall, H., Heine, S., Garner, A., et al. 2022, in 44th COSPAR Scientific As- sembly. Held 16-24 July, V ol. 44, 1871
2022
-
[30]
L., Liodakis, I., Marscher, A
Marshall, H. L., Liodakis, I., Marscher, A. P., et al. 2024, ApJ, 972, 74
2024
-
[31]
2023b, ApJ, 953, L28 Mücke, A
Middei, R., Perri, M., Puccetti, S., et al. 2023b, ApJ, 953, L28 Mücke, A. & Protheroe, R. J. 2001, Astroparticle Physics, 15, 121
2001
-
[32]
L., Abdo, A
Nolan, P. L., Abdo, A. A., Ackermann, M., et al. 2012, ApJS, 199, 31
2012
-
[33]
E., Middei, R., et al
Pacciani, L., Kim, D. E., Middei, R., et al. 2025, ApJ, 983, 78
2025
-
[34]
2014, MNRAS, 442, 1693
Pavlidou, V ., Angelakis, E., Myserlis, I., et al. 2014, MNRAS, 442, 1693
2014
-
[35]
L., Negro, M., Liodakis, I., et al
Peirson, A. L., Negro, M., Liodakis, I., et al. 2023, ApJ, 948, L25
2023
-
[36]
Raiteri, C. M. 2025, A&A Rev., 33, 8
2025
-
[37]
N., Rajarshi, C
Ramaprakash, A. N., Rajarshi, C. V ., Das, H. K., et al. 2019, MNRAS, 485, 2355
2019
-
[38]
P., Fuhrmann, L., et al
Rani, B., Krichbaum, T. P., Fuhrmann, L., et al. 2013, A&A, 552, A11
2013
-
[39]
2023, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, V ol
Soffitta, P., Baldini, L., Baumgartner, W., et al. 2023, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, V ol. 12678, UV , X- Ray, and Gamma-Ray Space Instrumentation for Astronomy XXIII, ed. O. H. Siegmund & K. Hoadley, 1267803
2023
-
[40]
2025, A&A, 700, A185
Tavecchio, F., Bolis, F., Sobacchi, E., Boula, S., & Sciaccaluga, A. 2025, A&A, 700, A185
2025
-
[41]
2021, in American Astronomical So- ciety Meeting Abstracts, V ol
Tomsick, J., Boggs, S., Zoglauer, A., et al. 2021, in American Astronomical So- ciety Meeting Abstracts, V ol. 237, American Astronomical Society Meeting Abstracts #237, 315.01
2021
-
[42]
2014, in 40th COSPAR Scientific As- sembly, V ol
Tomsick, J., Jean, P., Chang, H.-K., et al. 2014, in 40th COSPAR Scientific As- sembly, V ol. 40, PSB.1–14–14
2014
-
[43]
Urry, C. M. & Padovani, P. 1995, PASP, 107, 803
1995
-
[44]
C., Soffitta, P., Baldini, L., et al
Weisskopf, M. C., Soffitta, P., Baldini, L., et al. 2022, Journal of Astronomical
2022
-
[45]
& Böttcher, M
Zhang, H. & Böttcher, M. 2013, ApJ, 774, 18
2013
-
[46]
2024, ApJ, 967, 93
Zhang, H., Böttcher, M., & Liodakis, I. 2024, ApJ, 967, 93
2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.