Detectability of Polarized Gamma-ray Emission from Blazar Flares with COSI
Pith reviewed 2026-05-25 06:00 UTC · model grok-4.3
The pith
COSI is predicted to detect MeV polarization from up to six blazar flares over its two-year mission.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using 17 years of Fermi LAT observations of 1413 blazars, the authors identify a maximum of 787 sources with flaring episodes. They estimate the minimum detectable polarization MDP99 in the COSI energy band for each flare under a range of spectral assumptions and background conditions. Under baseline background levels and assuming MeV flare statistics match GeV observations, COSI can detect polarization in up to approximately 6 flares with MDP99 less than 50% over its two-year prime mission, with only a few reaching below 20%. Flat-spectrum radio quasars dominate the promising targets, and shorter intervals around bright peaks improve prospects.
What carries the argument
The MDP99 calculation using COSI instrument response functions applied to flare light curves and spectra identified from Fermi LAT data via Bayesian blocks.
If this is right
- Flat-spectrum radio quasars make up most of the polarization-detectable flares.
- Shorter time intervals around bright peaks within flares can yield better MDP99 values.
- A small number of the most powerful flares may achieve MDP99 below 20%.
- COSI's continuous monitoring will enable the first direct MeV polarization measurements from blazars.
Where Pith is reading between the lines
- Detections would allow testing whether the same jet regions produce the polarized emission across energy bands.
- If fewer than expected flares are detected, it may indicate that MeV flares have different duty cycles or spectra than assumed.
- This work could guide target selection for future MeV polarimeters or multi-wavelength campaigns.
- Non-detections might point to depolarization effects specific to the MeV regime.
Load-bearing premise
Blazar flare statistics and properties in the MeV band match those observed at GeV energies by Fermi LAT.
What would settle it
If COSI's two-year observations yield a number of polarized flare detections far from the predicted up to six, or if direct MeV flare rates differ markedly from the GeV-based extrapolation.
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 manuscript analyzes 17 years of Fermi LAT data for 1413 blazars, identifying a maximum of 787 flares via Bayesian block analysis. For each flare it computes the minimum detectable polarization (MDP99) in the COSI 0.2-5 MeV band using instrument response functions under varied spectral assumptions and background levels. Under the explicit assumption that MeV-band flare statistics match those observed at GeV energies, and for a baseline background of 1 count/s, the paper concludes that COSI could detect polarization in up to ~6 flares with MDP99 < 50% (a few reaching <20%) over its two-year prime mission, provides a ranked target list dominated by flat-spectrum radio quasars, and notes that shorter intervals around flare peaks improve prospects.
Significance. If the central extrapolation holds, the result supplies the first quantitative forecast for MeV polarization detections with COSI and supplies a practical ranked target list that could guide early observations. The work draws strength from its use of an extensive, publicly available Fermi LAT flare catalog and from direct application of COSI instrument response functions rather than purely theoretical estimates.
major comments (2)
- [Abstract] Abstract and the paragraph containing the ~6-flare claim: the headline number is obtained by transferring the full set of 787 GeV-identified flares (occurrence rate, duration distribution, and peak fluxes) directly to the COSI band. No independent MeV flare catalog, SED-based simulation, or duty-cycle correction is performed to test this transfer; if the MeV flare rate is lower (as expected when the synchrotron or external-Compton peak lies outside 0.2-5 MeV for many sources), the effective sample size and therefore the predicted detection count would decrease substantially.
- [MDP calculation procedure] MDP calculation procedure (the section describing the conversion from Fermi LAT flare properties to COSI MDP99): the final count of ~6 events is shown to vary with the choice of spectral index, background rate, and flare-selection threshold, yet the manuscript presents only a limited set of discrete cases rather than a continuous sensitivity study or error envelope on the extrapolated flare population. This makes it difficult to assess how robust the central prediction remains under plausible variations in the untested MeV extrapolation.
minor comments (2)
- [Flare identification] The description of the Bayesian block parameters (prior, false-positive rate) used to identify the 787 flares should be stated explicitly so that readers can reproduce the flare sample from the public Fermi LAT light curves.
- [Target list figure] Figure captions for the ranked target list should indicate the exact MDP99 threshold and background level adopted for each source so that the ordering can be directly compared with the text.
Simulated Author's Rebuttal
We thank the referee for their insightful comments, which help clarify the strengths and limitations of our analysis. We address the major comments point by point below, maintaining the data-driven nature of the study while improving clarity on assumptions.
read point-by-point responses
-
Referee: [Abstract] Abstract and the paragraph containing the ~6-flare claim: the headline number is obtained by transferring the full set of 787 GeV-identified flares (occurrence rate, duration distribution, and peak fluxes) directly to the COSI band. No independent MeV flare catalog, SED-based simulation, or duty-cycle correction is performed to test this transfer; if the MeV flare rate is lower (as expected when the synchrotron or external-Compton peak lies outside 0.2-5 MeV for many sources), the effective sample size and therefore the predicted detection count would decrease substantially.
Authors: We agree that the central prediction relies on the explicit assumption, stated in the abstract and throughout the manuscript, that MeV-band flare statistics match those at GeV energies. As no MeV flare catalog is available (COSI being a future mission and Fermi-LAT having limited sensitivity in the MeV range for this purpose), we cannot perform an independent verification or SED-based simulation without introducing further model dependencies. This is the strongest data-driven estimate possible with current observations. We will revise the abstract and discussion to more explicitly caution that if the MeV flare rate is lower, the number of detectable events would decrease, and we will reference this as a key uncertainty. revision: partial
-
Referee: [MDP calculation procedure] MDP calculation procedure (the section describing the conversion from Fermi LAT flare properties to COSI MDP99): the final count of ~6 events is shown to vary with the choice of spectral index, background rate, and flare-selection threshold, yet the manuscript presents only a limited set of discrete cases rather than a continuous sensitivity study or error envelope on the extrapolated flare population. This makes it difficult to assess how robust the central prediction remains under plausible variations in the untested MeV extrapolation.
Authors: The manuscript already demonstrates the dependence on spectral assumptions, background levels, and flare identification criteria through multiple discrete scenarios. To better address robustness, we will expand the analysis to include a more continuous exploration of parameter space, such as varying the spectral index over a range and showing the resulting distribution of detectable flares. This will provide a clearer sensitivity envelope without altering the core methodology or conclusions. revision: yes
Circularity Check
No circularity; forward predictions from external Fermi-LAT catalog under explicit extrapolation assumption.
full rationale
The central result (~6 detectable polarized flares) is obtained by applying COSI instrument response functions to flare parameters taken directly from an external 17-year Fermi-LAT analysis of 1413 blazars. The paper states the key transfer assumption explicitly rather than deriving it from its own equations. No self-citation chain, fitted parameter renamed as prediction, or self-definitional step appears in the derivation. The calculation is therefore a standard forward projection whose validity rests on the plausibility of the stated assumption, not on any internal reduction to the paper's own inputs.
Axiom & Free-Parameter Ledger
free parameters (2)
- baseline background level
- spectral assumptions
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]
Multimessenger observations of a flaring blazar coincident with high-energy neutrino IceCube-170922A
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]
2023, ApJS, 265, 31, doi: 10.3847/1538-4365/acbb6a
Abdollahi, S., Ajello, M., Baldini, L., et al. 2023, ApJS, 265, 31, doi: 10.3847/1538-4365/acbb6a
-
[4]
2025, ApJL, 985, L15, doi: 10.3847/2041-8213/adc572
Agudo, I., Liodakis, I., Otero-Santos, J., et al. 2025, ApJL, 985, L15, doi: 10.3847/2041-8213/adc572
-
[5]
Anderson, T. W., & Darling, D. A. 1954, Journal of the American Statistical Association, 49, 765
work page 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]
Fermi Large Area Telescope Fourth Source Catalog Data Release 4 (4FGL-DR4)
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
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2307.12546 2023
-
[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
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1503.06676 2015
-
[14]
Model for the High-Energy Emission from Blazars
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]
Ghisellini, G. 1998, MNRAS, 299, 433, doi: 10.1046/j.1365-8711.1998.01828.x
-
[18]
2026, ApJ, 997, 284, doi: 10.3847/1538-4357/ae32f4
Gallego, S., Oberlack, U., Lommler, J., et al. 2026, ApJ, 997, 284, doi: 10.3847/1538-4357/ae32f4
-
[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]
2021, MNRAS, 507, 5690, doi: 10.1093/mnras/stab2425
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]
1981, ApJ, 243, 700, doi: 10.1086/158638
Konigl, A. 1981, ApJ, 243, 700, doi: 10.1086/158638
-
[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]
Mann, H. B., & Whitney, D. R. 1947, The Annals of Mathematical Statistics, 18, 50
work page 1947
-
[25]
Mannheim, K. 1993, A&A, 269, 67, doi: 10.48550/arXiv.astro-ph/9302006
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.astro-ph/9302006 1993
-
[26]
1992, ApJL, 397, L5, doi: 10.1086/186531
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]
2023, ApJL, 942, L10, doi: 10.3847/2041-8213/aca281
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
work page internal anchor Pith review doi:10.3847/1538-4365/aaa7fd 2018
-
[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]
L., Negro, M., Liodakis, I., et al
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]
2019, ApJ, 885, 12, doi: 10.3847/1538-4357/ab426a 20
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
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1304.2818 2013
-
[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
work page 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]
2024, in 38th International Cosmic Ray Conference, 745, doi: 10.48550/arXiv.2308.12362
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]
Unified Schemes for Radio-Loud Active Galactic Nuclei
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]
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, 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, -
work page 2006
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.