REVIEW 2 major objections 2 minor 58 references
Detectability of Polarized Gamma-ray Emission from Blazar Flares with COSI
T0 review · 2 major / 2 minor · reviewed 2026-05-25 · grok-4.3
Pith's one-line read COSI is predicted to detect MeV polarization from up to six blazar flares over its two-year mission.
desk verdict The paper forecasts that COSI might catch polarized MeV emission from up to ~6 blazar flares in two years, but only if GeV flare rates and durations transfer directly to the MeV band. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
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.
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.
Extended reading notes
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.
Load-bearing premise
Blazar flare statistics and properties in the MeV band match those observed at GeV energies by Fermi LAT.
Editorial extensions
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.
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.
Signed reviews
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.
Assumptions & free parameters
free parameters (2)
- baseline background level
- spectral assumptions
assumptions (1)
- domain assumption blazar flare statistics in the MeV band are comparable to those observed at GeV energies
Cite this review
Pith. "Pith review of Detectability of Polarized Gamma-ray Emission from Blazar Flares with COSI." pith.science (2026). https://pith.science/paper/PLEN5G2U
@misc{pith2026260501139,
author = {Pith},
title = {Pith review of: Detectability of Polarized Gamma-ray Emission from Blazar Flares with COSI},
year = {2026},
howpublished = {\url{https://pith.science/paper/PLEN5G2U}},
note = {Machine review of arXiv:2605.01139}
}
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.
Figures
Figures from the paper (8 more)
Reference graph
Works this paper leans on
-
[1]
Science , archivePrefix = "arXiv", eprint =
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]
, archivePrefix = "arXiv", eprint =
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 Pith review arXiv doi:10.48550/arxiv.2307.12546 2023
Show all 58 references
-
[9]
H., Digel, S
Ballet, J., Burnett, T. H., Digel, S. W., & Lott, B. 2020, arXiv e-prints, arXiv:2005.11208, doi: 10.48550/arXiv.2005.11208
2020 doi
- [10]
-
[11]
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
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 ar...
2018 doi
-
[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
2013 doi
- [13]
-
[14]
D., & Schlickeiser, R
Dermer, C. D., & Schlickeiser, R. 1993, ApJ, 416, 458, doi: 10.1086/173251
1993 doi
-
[15]
M., & Papadakis, I
Emmanoulopoulos, D., McHardy, I. M., & Papadakis, I. E. 2013, MNRAS, 433, 907, doi: 10.1093/mnras/stt764
2013 doi
-
[16]
D., & Becker, P
Finke, J. D., & Becker, P. A. 2014, ApJ, 791, 21, doi: 10.1088/0004-637X/791/1/21
2014 doi
-
[17]
1998, MNRAS, 299, 433, doi: 10.1046/j.1365-8711.1998.01828.x
Ghisellini, G. 1998, MNRAS, 299, 433, doi: 10.1046/j.1365-8711.1998.01828.x
1998 doi
-
[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
2026 doi
-
[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
2017 doi
-
[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
2021 doi
-
[21]
G., et al
Hovatta, T., Pavlidou, V., King, O. G., et al. 2014, MNRAS, 439, 690, doi: 10.1093/mnras/stt2494
2014 doi
-
[22]
1981, ApJ, 243, 700, doi: 10.1086/158638
Konigl, A. 1981, ApJ, 243, 700, doi: 10.1086/158638
1981 doi
-
[23]
L., Aller, M
Lister, M. L., Aller, M. F., Aller, H. D., et al. 2015, ApJL, 810, L9, doi: 10.1088/2041-8205/810/1/L9
2015 doi
-
[24]
B., & Whitney, D
Mann, H. B., & Whitney, D. R. 1947, The Annals of Mathematical Statistics, 18, 50
1947
- [25]
-
[26]
1992, ApJL, 397, L5, doi: 10.1086/186531
Maraschi, L., Ghisellini, G., & Celotti, A. 1992, ApJL, 397, L5, doi: 10.1086/186531
1992 doi
-
[27]
M., et al
Marcotulli, L., Ajello, M., Urry, C. M., et al. 2022, ApJ, 940, 77, doi: 10.3847/1538-4357/ac937f
2022 doi
-
[28]
2023, in Proceedings of 38th International Cosmic Ray Conference — PoS(ICRC2023), ICRC2023 (Sissa Medialab), 858, doi: 10.22323/1.444.0858
Martinez, I. 2023, in Proceedings of 38th International Cosmic Ray Conference — PoS(ICRC2023), ICRC2023 (Sissa Medialab), 858, doi: 10.22323/1.444.0858
2023 doi
-
[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
2023 doi
-
[30]
D., & Blandford, R
Meyer, M., Scargle, J. D., & Blandford, R. D. 2019, ApJ, 877, 39, doi: 10.3847/1538-4357/ab1651
2019 doi
-
[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
2023 doi
-
[32]
2013, MNRAS, 430, 1324, doi: 10.1093/mnras/sts711
Nalewajko, K. 2013, MNRAS, 430, 1324, doi: 10.1093/mnras/sts711
2013 doi
-
[33]
B., et al
Oh, K., Koss, M., Markwardt, C. B., et al. 2018, ApJS, 235, 4, doi: 10.3847/1538-4365/aaa7fd
2018 doi
-
[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
2023 doi
-
[35]
M., Acosta-Pulido, J
Otero-Santos, J., Raiteri, C. M., Acosta-Pulido, J. A., et al. 2024, A&A, 686, A228, doi: 10.1051/0004-6361/202449647
2024 doi
-
[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
2023 doi
-
[37]
M., Villata, M., Jorstad, S
Raiteri, C. M., Villata, M., Jorstad, S. G., et al. 2023, MNRAS, 522, 102, doi: 10.1093/mnras/stad942
2023 doi
-
[38]
S., & Rakshit, S
Rajput, B., Stalin, C. S., & Rakshit, S. 2020, A&A, 634, A80, doi: 10.1051/0004-6361/201936769
2020 doi
-
[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
2019 doi
-
[40]
R., Penacchioni, A
Sacahui, J. R., Penacchioni, A. V., Marinelli, A., et al. 2021, RMxAA, 57, 251, doi: 10.22201/ia.01851101p.2021.57.02.01
2021 doi
- [41]
-
[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
2025
-
[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]
2020, ApJS, 250, 1, doi: 10.3847/1538-4365/aba2c7 The Compton Spectrometer and Imager (COSI) Collaboration
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
2020 doi
-
[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
2024 doi
-
[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
2021 doi
- [47]
-
[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
2022 doi
-
[49]
C., Elsner, R
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.1...
2010 doi
-
[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
2022 doi
-
[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
2023 doi
-
[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
2013 doi
-
[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
2024 doi
-
[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
2016 doi
-
[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
2019 doi
-
[56]
E., et al
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/...
2008 doi
-
[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
2021 doi
-
[58]
Zoglauer, A. C. 2006, PhD thesis, -
2006
Reviewed May 25, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.