Recognition: unknown
Detectability of Polarized Gamma-ray Emission from Blazar Flares with COSI
Pith reviewed 2026-05-09 18:11 UTC · model grok-4.3
The pith
COSI can detect MeV polarization from up to six blazar flares in two years.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
With baseline background rates and flare statistics in the MeV band assumed comparable to GeV observations, COSI is projected to measure polarization in approximately six blazar flares with MDP99 below 50 percent over its two-year mission. Only the strongest flares reach MDP99 below 20 percent. Detection prospects improve by analyzing shorter segments around flare peaks, and the most promising targets are predominantly flat-spectrum radio quasars.
What carries the argument
Calculation of the 99 percent minimum detectable polarization (MDP99) using the instrument's response functions for different assumed spectra and backgrounds, applied to each identified flare to determine detectability.
If this is right
- Flat-spectrum radio quasars make up most of the detectable events.
- Focusing on peak intervals within flares boosts the expected detections.
- COSI will provide the first direct MeV-band polarization measurements from blazars.
- These data will constrain models of jet geometry and particle acceleration.
Where Pith is reading between the lines
- Polarization measurements could distinguish between competing models for the high-energy emission in blazar jets.
- If the MeV flare population differs from GeV, the number of detections may vary from the estimate.
- Combining with other wavelengths could give a fuller picture of flare evolution.
Load-bearing premise
That the rate and characteristics of flares observed in blazars at GeV energies apply similarly in the lower MeV energy band.
What would settle it
Observing far fewer than six polarized flares in the first two years of COSI data, or finding no polarization in the highest-ranked candidates, would indicate that the assumptions about MeV flare statistics or polarization levels do not hold.
Figures
read the original abstract
We investigate the detectability of polarized gamma-ray emission from blazar flares with the Compton Spectrometer and Imager (COSI). Using 17 years of Fermi Large Area Telescope observations, we analyze light curves for 1413 blazars and identify a maximum of 787 sources with flaring episodes through Bayesian block analysis. For each flare, we estimate the minimum detectable polarization MDP99 in the COSI energy band (0.2-5 MeV) using instrument response functions under a range of spectral assumptions and background conditions. Under baseline background levels (1 counts/s), and assuming that blazar flare statistics in the MeV band are comparable to those observed at GeV energies, we find that COSI can realistically detect polarization in up to ~6 flares with MDP99<50% over its two-year prime mission depending on different spectral and flare identification assumptions, with only a few most powerful ones reaching MDP99<20%. These expectations are shown to improve when shorter intervals around bright peaks within long flares are considered. We provide a ranked list of the most promising targets, finding that flat-spectrum radio quasars dominate the population of polarization-detectable events. Through its continuous all-sky monitoring in the largely unexplored MeV band, COSI will open a new observational window on blazar variability and deliver the first direct measurements of MeV polarization, offering unique insights into jet geometry and high-energy emission processes.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that by analyzing 17 years of Fermi-LAT observations of 1413 blazars and identifying flaring episodes with Bayesian blocks (yielding up to 787 flaring sources), and then estimating the minimum detectable polarization (MDP99) in the 0.2-5 MeV band for COSI using instrument response functions under various spectral and background assumptions, COSI can detect polarized emission in up to approximately 6 blazar flares over its two-year prime mission with MDP99 < 50%, and a few with <20%, assuming that MeV flare statistics are similar to those at GeV energies. A ranked list of the most promising targets, dominated by flat-spectrum radio quasars, is provided.
Significance. If the central assumption regarding flare statistics holds, this work would be significant for providing the first quantitative predictions for detecting polarized gamma-ray emission from blazar flares in the MeV band with COSI. It leverages extensive public Fermi data and detailed instrument modeling to offer actionable insights for the mission, including a target list that could guide observations and advance understanding of jet geometry and high-energy processes in blazars.
major comments (2)
- [Abstract] The headline number of up to ~6 flares with MDP99<50% (and a few <20%) is derived under the assumption that blazar flare statistics in the MeV band are comparable to those observed at GeV energies. This assumption is load-bearing for the central claim, as the paper does not provide validation, error analysis, or sensitivity tests for the GeV-to-MeV extrapolation of flare rates, amplitudes, and duty cycles. Given known energy-dependent variations in blazar SEDs and flare hardness, this could substantially alter the predicted yield.
- [Methods (flare identification and MDP99 estimation)] The procedure folds GeV-derived peak flux, duration, and spectrum through COSI response matrices to obtain MDP99, but lacks quantitative assessment of uncertainties in the spectral extrapolation or how different flare identification assumptions quantitatively affect the final count of six flares.
minor comments (2)
- [Abstract] The phrase 'depending on different spectral and flare identification assumptions' is vague; specifying the range of assumptions and their impact on the ~6 number would improve clarity.
- Consider adding a reference to prior studies on MeV blazar emission or polarization expectations to contextualize the novelty.
Simulated Author's Rebuttal
We thank the referee for the constructive feedback and for acknowledging the potential significance of our work. We address the major comments point by point below. We agree that additional quantitative sensitivity analyses strengthen the paper and have revised the manuscript to include them.
read point-by-point responses
-
Referee: [Abstract] The headline number of up to ~6 flares with MDP99<50% (and a few <20%) is derived under the assumption that blazar flare statistics in the MeV band are comparable to those observed at GeV energies. This assumption is load-bearing for the central claim, as the paper does not provide validation, error analysis, or sensitivity tests for the GeV-to-MeV extrapolation of flare rates, amplitudes, and duty cycles. Given known energy-dependent variations in blazar SEDs and flare hardness, this could substantially alter the predicted yield.
Authors: We agree the GeV-to-MeV extrapolation assumption is central and that more explicit sensitivity testing is valuable. Direct validation is not feasible without MeV flare catalogs, but we have tested variations in flare amplitudes, durations, and duty cycles informed by multi-wavelength studies. The revised manuscript adds a new subsection and Figure 8 showing the detectable flare count ranges from 2 to 12 for scaling factors 0.5x–2x the GeV rate, and from 3 to 9 under plausible amplitude/duration changes. The abstract has been updated to emphasize this range. revision: yes
-
Referee: [Methods (flare identification and MDP99 estimation)] The procedure folds GeV-derived peak flux, duration, and spectrum through COSI response matrices to obtain MDP99, but lacks quantitative assessment of uncertainties in the spectral extrapolation or how different flare identification assumptions quantitatively affect the final count of six flares.
Authors: We acknowledge the original presentation focused on baseline cases. The revised manuscript now includes explicit quantification: we tested three Bayesian block priors and two significance thresholds, producing 4–9 flares with MDP99<50%. For spectral extrapolation we varied photon indices (1.8–2.5) and added cutoff energies (10 MeV, 100 MeV), with results summarized in new Table 3. These show the count changes by at most 2–3 sources, confirming the headline result is stable within the explored range while transparently displaying the uncertainties. revision: yes
- Full observational validation of the assumption that MeV flare rates, amplitudes, and duty cycles match GeV statistics, which cannot be performed without comprehensive existing MeV flare monitoring data.
Circularity Check
No significant circularity; forecast uses external data and explicit assumption
full rationale
The paper derives its ~6-flare prediction by running Bayesian blocks on external 17-year Fermi-LAT light curves of 1413 blazars to extract flares, then folding GeV peak fluxes, durations and spectra through COSI response matrices to obtain MDP99 values. The single scaling step (MeV flare statistics assumed identical to GeV) is stated explicitly as an assumption rather than fitted or derived inside the work. No equation reduces to a self-definition, no fitted parameter is relabeled as a prediction, and no load-bearing premise rests on a self-citation chain. The derivation chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- Baseline background rate =
1 counts/s
axioms (1)
- domain assumption Blazar flare statistics in the MeV band are comparable to those observed at GeV energies
Reference graph
Works this paper leans on
-
[1]
Aartsen, M., Ackermann, M., Adams, J., et al. 2018, Science, 361, eaat1378, doi: 10.1126/science.aat1378
-
[2]
Abdollahi, S., Acero, F., Ackermann, M., et al. 2020, ApJS, 247, 33, doi: 10.3847/1538-4365/ab6bcb 18 T able 4.List of BB-selected flares identified below 50% MDP 99% for COSI for our baseline background case (η= 0.3, Rbkg = 1 ct/s) and a softer X-ray spectrum within Fig. 3, with their associated MDP 99%, average and peak flare flux, and flare duration wi...
-
[3]
Abdollahi, S., Ajello, M., Baldini, L., et al. 2023, ApJS, 265, 31, doi: 10.3847/1538-4365/acbb6a
-
[4]
Agudo, I., Liodakis, I., Otero-Santos, J., et al. 2025, ApJL, 985, L15, doi: 10.3847/2041-8213/adc572
-
[5]
W., & Darling, D
Anderson, T. W., & Darling, D. A. 1954, Journal of the American Statistical Association, 49, 765
1954
-
[6]
2016, MNRAS, 463, 3365, doi: 10.1093/mnras/stw2217
Angelakis, E., Hovatta, T., Blinov, D., et al. 2016, MNRAS, 463, 3365, doi: 10.1093/mnras/stw2217
-
[7]
Atwood, W. B., Abdo, A. A., Ackermann, M., et al. 2009, ApJ, 697, 1071, doi: 10.1088/0004-637X/697/2/1071
-
[8]
H., Lott, B., & collaboration, T
Ballet, J., Bruel, P., Burnett, T. H., Lott, B., & The Fermi-LAT collaboration. 2023, arXiv e-prints, arXiv:2307.12546, doi: 10.48550/arXiv.2307.12546
-
[9]
Ballet, J., Burnett, T. H., Digel, S. W., & Lott, B. 2020, arXiv e-prints, arXiv:2005.11208, doi: 10.48550/arXiv.2005.11208
-
[10]
Blandford, R. D., & K¨ onigl, A. 1979, ApJ, 232, 34, doi: 10.1086/157262
-
[11]
Blinov, D., Pavlidou, V., Papadakis, I., et al. 2018, MNRAS, 474, 1296, doi: 10.1093/mnras/stx2786 19 Figure 11.Example of a Long-flare (gray shaded region) identified in the Fermi-LAT light curve of source 4FGL J2253.9+1609. Manually-selected time windows with bright peaks are highlighted in orange and blue shades. Dotted red lines represent the average ...
-
[12]
2013, ApJ, 769, 70, doi: 10.1088/0004-637X/769/1/70
Chang, Z., Jiang, Y., & Lin, H.-N. 2013, ApJ, 769, 70, doi: 10.1088/0004-637X/769/1/70
-
[13]
Connolly, S. D. 2015, arXiv e-prints, arXiv:1503.06676, doi: 10.48550/arXiv.1503.06676
-
[14]
Dermer, C. D., & Schlickeiser, R. 1993, ApJ, 416, 458, doi: 10.1086/173251
-
[15]
Emmanoulopoulos, D., McHardy, I. M., & Papadakis, I. E. 2013, MNRAS, 433, 907, doi: 10.1093/mnras/stt764
-
[16]
Finke, J. D., & Becker, P. A. 2014, ApJ, 791, 21, doi: 10.1088/0004-637X/791/1/21
-
[17]
A., Ramsay, G., Andronov, I., et al
Ghisellini, G. 1998, MNRAS, 299, 433, doi: 10.1046/j.1365-8711.1998.01828.x
-
[18]
2025, arXiv e-prints, arXiv:2510.25304, doi: 10.48550/arXiv.2510.25304
Gallego, S., Oberlack, U., Lommler, J., et al. 2025, arXiv e-prints, arXiv:2510.25304, doi: 10.48550/arXiv.2510.25304
-
[19]
2017, MNRAS, 469, 255, doi: 10.1093/mnras/stx806
Ghisellini, G., Righi, C., Costamante, L., & Tavecchio, F. 2017, MNRAS, 469, 255, doi: 10.1093/mnras/stx806
-
[20]
Giommi, P., Perri, M., Capalbi, M., et al. 2021, MNRAS, 507, 5690, doi: 10.1093/mnras/stab2425
-
[21]
Hovatta, T., Pavlidou, V., King, O. G., et al. 2014, MNRAS, 439, 690, doi: 10.1093/mnras/stt2494
- [22]
-
[23]
Lister, M. L., Aller, M. F., Aller, H. D., et al. 2015, ApJL, 810, L9, doi: 10.1088/2041-8205/810/1/L9
-
[24]
B., & Whitney, D
Mann, H. B., & Whitney, D. R. 1947, The Annals of Mathematical Statistics, 18, 50
1947
-
[25]
1993, A&A, 269, 67, doi: 10.48550/arXiv.astro-ph/9302006
Mannheim, K. 1993, A&A, 269, 67, doi: 10.48550/arXiv.astro-ph/9302006
-
[26]
Maraschi, L., Ghisellini, G., & Celotti, A. 1992, ApJL, 397, L5, doi: 10.1086/186531
-
[27]
Marcotulli, L., Ajello, M., Urry, C. M., et al. 2022, ApJ, 940, 77, doi: 10.3847/1538-4357/ac937f
-
[28]
Martinez, I. 2023, in Proceedings of 38th International Cosmic Ray Conference — PoS(ICRC2023), ICRC2023 (Sissa Medialab), 858, doi: 10.22323/1.444.0858
-
[29]
2023, arXiv e-prints, arXiv:2308.11436, doi: 10.48550/arXiv.2308.11436
Martinez-Castellanos, I., Gallego, S., Huang, C.-Y., et al. 2023, arXiv e-prints, arXiv:2308.11436, doi: 10.48550/arXiv.2308.11436
-
[30]
Meyer, M., Scargle, J. D., & Blandford, R. D. 2019, ApJ, 877, 39, doi: 10.3847/1538-4357/ab1651
-
[31]
Middei, R., Liodakis, I., Perri, M., et al. 2023, ApJL, 942, L10, doi: 10.3847/2041-8213/aca281
-
[32]
2013, MNRAS, 430, 1324, doi: 10.1093/mnras/sts711
Nalewajko, K. 2013, MNRAS, 430, 1324, doi: 10.1093/mnras/sts711
-
[33]
Oh, K., Koss, M., Markwardt, C. B., et al. 2018, ApJS, 235, 4, doi: 10.3847/1538-4365/aaa7fd
-
[34]
A., Becerra Gonz´ alez, J., et al
Otero-Santos, J., Acosta-Pulido, J. A., Becerra Gonz´ alez, J., et al. 2023, MNRAS, 523, 4504, doi: 10.1093/mnras/stad1722
-
[35]
Otero-Santos, J., Raiteri, C. M., Acosta-Pulido, J. A., et al. 2024, A&A, 686, A228, doi: 10.1051/0004-6361/202449647
-
[36]
Peirson, A. L., Negro, M., Liodakis, I., et al. 2023, ApJL, 948, L25, doi: 10.3847/2041-8213/acd242
-
[37]
Raiteri, C. M., Villata, M., Jorstad, S. G., et al. 2023, MNRAS, 522, 102, doi: 10.1093/mnras/stad942
-
[38]
Rajput, B., Stalin, C. S., & Rakshit, S. 2020, A&A, 634, A80, doi: 10.1051/0004-6361/201936769
-
[39]
L., Siemiginowska A., Sobolewska M
Grindlay, J. 2019, ApJ, 885, 12, doi: 10.3847/1538-4357/ab426a 20
-
[40]
Sacahui, J. R., Penacchioni, A. V., Marinelli, A., et al. 2021, RMxAA, 57, 251, doi: 10.22201/ia.01851101p.2021.57.02.01
-
[41]
Scargle, J. D., Norris, J. P., Jackson, B., & Chiang, J. 2013, arXiv e-prints, arXiv:1304.2818, doi: 10.48550/arXiv.1304.2818
-
[42]
2025, arXiv e-prints, arXiv:2511.09420
Sciaccaluga, A., Tavecchio, F., Sbarrato, T., & Bonnoli, G. 2025, arXiv e-prints, arXiv:2511.09420. https://arxiv.org/abs/2511.09420
-
[43]
Smith, P. S. 1996, in Astronomical Society of the Pacific Conference Series, Vol. 110, Blazar Continuum Variability, ed. H. R. Miller, J. R. Webb, & J. C. Noble, 135
1996
-
[44]
Pascual-Granado, J. 2020, ApJS, 250, 1, doi: 10.3847/1538-4365/aba2c7 The Compton Spectrometer and Imager (COSI) Collaboration. 2025, 3.1 Zenodo, doi: 10.5281/zenodo.15126188
-
[45]
Tomsick, J., Boggs, S., Zoglauer, A., et al. 2024, in 38th International Cosmic Ray Conference, 745, doi: 10.48550/arXiv.2308.12362
-
[46]
2021, ApJ, 916, 28, doi: 10.3847/1538-4357/ac0341
Tsuji, N., Yoneda, H., Inoue, Y., et al. 2021, ApJ, 916, 28, doi: 10.3847/1538-4357/ac0341
-
[47]
Urry, C. M., & Padovani, P. 1995, PASP, 107, 803, doi: 10.1086/133630
work page internal anchor Pith review doi:10.1086/133630 1995
-
[48]
M., Burd, P., Dorner, D., et al
Wagner, S. M., Burd, P., Dorner, D., et al. 2022, in 37th International Cosmic Ray Conference, 868, doi: 10.22323/1.395.0868
-
[49]
Weisskopf, M. C., Elsner, R. F., & O’Dell, S. L. 2010, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7732, Space Telescopes and Instrumentation 2010: Ultraviolet to Gamma Ray, ed. M. Arnaud, S. S. Murray, & T. Takahashi, 77320E, doi: 10.1117/12.857357
-
[50]
Journal of Astronomical Telescopes, Instruments, and Systems , keywords =
Weisskopf, M. C., Soffitta, P., Baldini, L., et al. 2022, Journal of Astronomical Telescopes, Instruments, and Systems, 8, 026002, doi: 10.1117/1.JATIS.8.2.026002
-
[51]
2023, ApJ, 954, 194, doi: 10.3847/1538-4357/acea74
Yoshida, K., Petropoulou, M., Murase, K., & Oikonomou, F. 2023, ApJ, 954, 194, doi: 10.3847/1538-4357/acea74
-
[52]
2013, ApJ, 774, 18, doi: 10.1088/0004-637X/774/1/18
Zhang, H., & B¨ ottcher, M. 2013, ApJ, 774, 18, doi: 10.1088/0004-637X/774/1/18
-
[53]
2024, ApJ, 967, 93, doi: 10.3847/1538-4357/ad4112
Zhang, H., B¨ ottcher, M., & Liodakis, I. 2024, ApJ, 967, 93, doi: 10.3847/1538-4357/ad4112
-
[54]
2016, ApJ, 829, 69, doi: 10.3847/0004-637X/829/2/69
Zhang, H., Diltz, C., & B¨ ottcher, M. 2016, ApJ, 829, 69, doi: 10.3847/0004-637X/829/2/69
-
[55]
2019, ApJ, 876, 109, doi: 10.3847/1538-4357/ab158d
Zhang, H., Fang, K., Li, H., et al. 2019, ApJ, 876, 109, doi: 10.3847/1538-4357/ab158d
-
[56]
Zoglauer, A., Andritschke, R., Boggs, S. E., et al. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7011, Space Telescopes and Instrumentation 2008: Ultraviolet to Gamma Ray, ed. M. J. L. Turner & K. A. Flanagan, 70113F, doi: 10.1117/12.789537
-
[57]
2021, arXiv e-prints, arXiv:2102.13158, doi: 10.48550/arXiv.2102.13158
Zoglauer, A., Siegert, T., Lowell, A., et al. 2021, arXiv e-prints, arXiv:2102.13158, doi: 10.48550/arXiv.2102.13158
-
[58]
Zoglauer, A. C. 2006, PhD thesis, -
2006
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.