pith. sign in

arxiv: 2604.03847 · v1 · submitted 2026-04-04 · 🌌 astro-ph.HE

Characterizing Gamma-Radio Delayed Flaring Activity from Blazars

Alina Kochocki , Emma Kun , Sam Hori This is my paper

Pith reviewed 2026-05-13 16:52 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords blazarsgamma-ray flaresradio flarestime delaysjet propagationFermi-LATRATAN-600MOJAVE
0
0 comments X

The pith

Blazars show radio flares delayed by about 180 days after gamma-ray activity in stacked analysis.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper looks for time delays between gamma-ray flares and later radio flares in roughly 100 blazars, motivated by longer delays seen in some neutrino-bright objects. Gaussian process models smooth the Fermi-LAT and RATAN-600 light curves, allowing cross-correlation searches across the full sample and a MOJAVE morphological sub-sample. The stacked result peaks at a radio delay near 180 days, while several individual sources show longer 1-3 year delays and the MOJAVE subset yields a modestly shorter lag. These patterns suggest radio emission often arises farther downstream in the jet than gamma-ray emission. Establishing the typical delay scale would clarify how particle acceleration and jet loading evolve on sub-parsec scales.

Core claim

Using Fermi-LAT gamma-ray data in three bands and RATAN-600 radio observations, the authors apply Gaussian process smoothing to predict smooth light-curve functions and compute cross-correlations. In the broader sample several sources exhibit radio flares delayed by 1-3 years after gamma-ray activity. The stacked correlation across the ensemble reaches its maximum when the radio data are shifted by approximately 180 days. The MOJAVE morphological sub-selection produces a comparable but slightly smaller lag. The authors conclude that delayed radio flares or extended radio emission constitute a notable feature within the general blazar population.

What carries the argument

Stacked cross-correlation of Gaussian-process-smoothed gamma-ray and radio light curves, compared with a MOJAVE morphological sub-sample.

Load-bearing premise

That the measured correlations and lags reflect genuine physical propagation delays in the jet rather than artifacts from Gaussian process smoothing, data sampling gaps, or post-hoc selection of the lag that maximizes the stacked signal.

What would settle it

Re-processing the same light curves with a different smoothing kernel or without selecting the lag that maximizes the stack, and finding no significant correlation peak near 180 days, would indicate the result is likely an artifact.

Figures

Figures reproduced from arXiv: 2604.03847 by Alina Kochocki, Emma Kun, Sam Hori.

Figure 1
Figure 1. Figure 1: Light curves for TXS 0506+056 and PKS 1424+240. In the top plot, we present available RATAN-600 light curves at several frequencies as well as Fermi-LAT data for TXS 0506+056 (Sotnikova et al. 2022; Abdollahi et al. 2023). In the bottom plot, we provide MOJAVE and Fermi-LAT observations for PKS 1424+240 (Lister et al. 2018; Abdollahi et al. 2023). The left axis corresponds to the photon flux of the gamma-r… view at source ↗
Figure 2
Figure 2. Figure 2: Representative light curves and uncertainties predicted through Gaussian process modeling. Here, we provide several light curve examples produced from Gaussian process modeling along with the underlying data. The name of the source is provided as the title in each case. The predicted means are plotted, as well as the 1σ posterior predictive uncertainties as shaded regions. The Fermi photon flux high-energy… view at source ↗
Figure 3
Figure 3. Figure 3: Best-fit time lags and estimated significances for a delayed gamma-ray-radio flare. We find that only a small selection of sources shows some evidence for temporal correlation on this delayed timescale of 0.5–3.5 years. We plot these highest-correlation lags, Tdelay, in the top left corner. Results for all gamma-ray energy bands are provided. We also provide distributions of estimated significances for eac… view at source ↗
Figure 4
Figure 4. Figure 4: Correlation coefficient as a function of time lag. Here we plot the correlation coefficient for the three gamma-ray energy ranges considered in this work. At the top, the correlation result for the set of all sources is shown. At the bottom, the result from the MOJAVE selection is plotted. The shaded bands indicate the 95% confidence regions determined for each light curve selection. The highest correlatio… view at source ↗
read the original abstract

Flaring activity from the jets of active galactic nuclei has been studied for several decades, closely related to the loading and evolution of the jet. In this work, we focus on the sub-hundred parsec jet region, well traced by non-thermal radio and gamma-ray emission. Only in recent years have light curves capturing the decade-long behavior of such sources become available for a large ensemble of objects. While previous studies have focused on a direct correlation or few-month lag between gamma-ray and radio activity, recent neutrino-bright blazars observed by the IceCube Neutrino Observatory present multi-year delays between initial gamma-ray activity and subsequent radio flares. In this work, we search for similar-timescale correlations between Fermi-LAT gamma-ray data and RATAN-600 radio data from ~100 blazars. We consider two gamma-ray bands, 100 MeV-1 GeV and 1 GeV-500 GeV, as well as the integral band, to compare correlations between potential opaque and unabsorbed regions of the jet. Gaussian process modeling is used for smooth light curve function prediction. We also analyze morphological AGN core data from the MOJAVE survey, forming a sub-selection to better illustrate potential dependence on location. In the broader selection, several sources exhibit delayed flares on the order of 1-3 years. In the stacked analysis, we find the highest correlation for a radio delay on the order of 180 days. The stacked correlation resulting from the MOJAVE sub-selection corresponds to a slightly smaller time lag. Delayed radio flares or extended radio emission appear to be notable features within the general blazar population.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 1 minor

Summary. The manuscript analyzes Fermi-LAT gamma-ray and RATAN-600 radio light curves for ~100 blazars, using Gaussian process interpolation to search for delayed radio flares following gamma-ray activity. It reports individual sources with 1-3 year delays and, in a stacked correlation analysis, identifies the strongest signal at a radio lag of ~180 days; a MOJAVE morphological sub-sample yields a slightly shorter lag. The work compares results across gamma-ray energy bands and links the findings to jet physics and neutrino-bright blazars.

Significance. If the reported lag survives proper statistical controls, the result would be of moderate significance for AGN jet studies, extending prior work on gamma-radio correlations to longer timescales and a larger sample while incorporating multi-band gamma-ray data and MOJAVE core information. The absence of quantitative error bars, control samples, or robustness tests in the current version limits the strength of the claim.

major comments (3)
  1. [Abstract / stacked analysis] Abstract and stacked-analysis description: the highest-correlation lag of ~180 days is identified after scanning delays on GP-interpolated light curves, yet no Monte Carlo null distribution is presented that preserves the observed sampling, variability, and GP kernel while destroying physical correlations. Without this, the peak cannot be distinguished from a look-elsewhere artifact induced by lag scanning and smoothing.
  2. [Methods / Results] Methods and results sections: the manuscript provides no quantitative assessment of how the measured lags or stacked correlation change under variations in GP kernel choice, hyperparameter priors, or data-gap handling, nor does it report control samples of randomized or shuffled light curves. These omissions directly affect the central claim that the 180-day feature reflects physical propagation rather than methodological bias.
  3. [Results] Results section: individual sources are stated to show 1-3 year delays, but no per-source significance, error bars on the lags, or comparison to the distribution expected from uncorrelated series with the same cadence are supplied. This weakens the supporting evidence for the stacked result.
minor comments (1)
  1. [Abstract] The abstract would benefit from a brief statement of the number of sources, the exact lag range scanned, and the correlation metric used in the stack.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We agree that additional statistical controls are needed to support the central claims and have revised the paper to incorporate Monte Carlo null distributions, robustness tests, and per-source significance assessments.

read point-by-point responses

  1. Referee: [Abstract / stacked analysis] Abstract and stacked-analysis description: the highest-correlation lag of ~180 days is identified after scanning delays on GP-interpolated light curves, yet no Monte Carlo null distribution is presented that preserves the observed sampling, variability, and GP kernel while destroying physical correlations. Without this, the peak cannot be distinguished from a look-elsewhere artifact induced by lag scanning and smoothing.

    Authors: We agree that a Monte Carlo null distribution is required to evaluate the significance of the 180-day peak. In the revised manuscript we have added a dedicated subsection describing 1000 surrogate datasets generated by shuffling the radio fluxes while preserving the exact sampling times, variability amplitudes, and GP kernel hyperparameters. The resulting null distribution places the observed stacked correlation at the 180-day lag above the 99th percentile, and we have updated both the abstract and results text to report this significance level. revision: yes

  2. Referee: [Methods / Results] Methods and results sections: the manuscript provides no quantitative assessment of how the measured lags or stacked correlation change under variations in GP kernel choice, hyperparameter priors, or data-gap handling, nor does it report control samples of randomized or shuffled light curves. These omissions directly affect the central claim that the 180-day feature reflects physical propagation rather than methodological bias.

    Authors: We have expanded the Methods section with quantitative robustness checks. We re-ran the full analysis using both squared-exponential and Matérn-3/2 kernels, varied the hyperparameter priors over a factor of three in length-scale and variance, and tested two gap-handling approaches (linear interpolation versus masking). The 180-day peak remains within ±25 days across all variants. We also include shuffled-light-curve control samples in the Results; these controls show no comparable peak, supporting that the feature is not an artifact of the GP procedure. revision: yes

  3. Referee: [Results] Results section: individual sources are stated to show 1-3 year delays, but no per-source significance, error bars on the lags, or comparison to the distribution expected from uncorrelated series with the same cadence are supplied. This weakens the supporting evidence for the stacked result.

    Authors: We now report per-source lag uncertainties derived from the GP posterior covariance and compare each lag to a null distribution obtained from 500 randomized radio series that retain the original cadence and GP kernel. For the sources previously noted as having 1–3 year delays we list the measured lag, its 1σ uncertainty, and the empirical p-value from the null distribution; most exceed 2σ. These values and the corresponding null histograms are added to the revised Results section and supplementary figures. revision: yes

Circularity Check

0 steps flagged

No significant circularity; lags measured directly from cross-correlation of observed data

full rationale

The paper reports an observational analysis of gamma-ray and radio light curves from ~100 blazars using Gaussian process interpolation followed by cross-correlation to identify time lags. No equations, derivations, or first-principles results are presented that reduce to fitted inputs by construction. The central claim (highest stacked correlation at ~180 days) is a direct numerical output from the data processing pipeline rather than a self-referential prediction or renamed fit. Any self-citations concern prior observational work and are not load-bearing for the reported lags. The analysis is self-contained against external benchmarks (Fermi-LAT, RATAN-600, MOJAVE data) with no ansatz smuggling or uniqueness theorems invoked.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The claim rests on the domain assumption that Gaussian process interpolation faithfully recovers flare timing without introducing spurious lags, plus the implicit assumption that radio and gamma-ray variability trace the same jet disturbances.

free parameters (1)
  • radio delay lag
    The 180-day value is the lag that maximizes the stacked correlation coefficient and is therefore fitted to the data.
axioms (1)
  • domain assumption Gaussian processes provide unbiased smooth reconstructions of blazar light curves
    Invoked to justify using GP modeling for flare timing extraction.

pith-pipeline@v0.9.0 · 5598 in / 1299 out tokens · 49056 ms · 2026-05-13T16:52:00.483280+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

33 extracted references · 33 canonical work pages

  1. [1]

    Neutrino emission from the direction of the blazar TXS 0506+056 prior to the IceCube-170922A alert

    Aartsen, M., Ackermann, M., Adams, J., et al. 2018, Science, 361, 147, 10.1126/science.aat2890

    work page doi:10.1126/science.aat2890 2018

  2. [2]

    and others

    ---. 2020a, Physical Review Letters, 124, 10.1103/physrevlett.124.051103

    work page doi:10.1103/physrevlett.124.051103

  3. [3]

    2023, ApJS, 265, 31, doi: 10.3847/1538-4365/acbb6a

    Abdollahi, S., Ajello, M., Baldini, L., et al. 2023, The Astrophysical Journal Supplement Series, 265, 31, 10.3847/1538-4365/acbb6a

    work page doi:10.3847/1538-4365/acbb6a 2023

  4. [4]

    T., Bosch-Ramon, V., & Romero, G

    Araudo, A. T., Bosch-Ramon, V., & Romero, G. E. 2010, Astronomy & Astrophysics, 522, A97, 10.1051/0004-6361/201014660

    work page doi:10.1051/0004-6361/201014660 2010

  5. [5]

    B., Abdo, A

    Atwood , W. B., Abdo , A. A., Ackermann , M., et al. 2009, , 697, 1071, 10.1088/0004-637X/697/2/1071

    work page doi:10.1088/0004-637x/697/2/1071 2009

  6. [6]

    D., & K¨ onigl, A

    Blandford , R. D., & K \"o nigl , A. 1979, , 232, 34, 10.1086/157262

    work page doi:10.1086/157262 1979

  7. [7]

    D., & Znajek, R

    Blandford, R. D., & Znajek, R. L. 1977, Mon. Not. Roy. Astron. Soc., 179, 433, 10.1093/mnras/179.3.433

    work page doi:10.1093/mnras/179.3.433 1977

  8. [8]

    Reviews of Modern Physics , year = 1970, month = jan, volume =

    Blumenthal, G. R., & Gould, R. J. 1970, Rev. Mod. Phys., 42, 237, 10.1103/RevModPhys.42.237

    work page doi:10.1103/revmodphys.42.237 1970

  9. [9]

    2018, Galaxies, 7, 3, 10.3390/galaxies7010003

    Boula, S., Petropoulou, M., & Mastichiadis, A. 2018, Galaxies, 7, 3, 10.3390/galaxies7010003

    work page doi:10.3390/galaxies7010003 2018

  10. [10]

    2020, Galaxies, 8, 72, 10.3390/galaxies8040072

    Cerruti, M. 2020, Galaxies, 8, 72, 10.3390/galaxies8040072

    work page doi:10.3390/galaxies8040072 2020

  11. [11]

    2018, Research Notes of the American Astronomical Society, 2, 31, doi: 10.3847/2515-5172/aaaf6c

    Foreman-Mackey , D. 2018, Research Notes of the American Astronomical Society, 2, 31, 10.3847/2515-5172/aaaf6c

    work page doi:10.3847/2515-5172/aaaf6c 2018

  12. [12]

    2014, Monthly Notices of the Royal Astronomical Society, 441, 1899–1909, 10.1093/mnras/stu540

    Fuhrmann, L., Larsson, S., Chiang, J., et al. 2014, Monthly Notices of the Royal Astronomical Society, 441, 1899–1909, 10.1093/mnras/stu540

    work page doi:10.1093/mnras/stu540 2014

  13. [13]

    2019, New A Rev., 87, 101541, doi: 10.1016/j.newar.2020.101541

    Hovatta , T., & Lindfors , E. 2019, , 87, 101541, 10.1016/j.newar.2020.101541

    work page doi:10.1016/j.newar.2020.101541 2019

  14. [14]

    G., Pushkarev, A

    Kramarenko, I. G., Pushkarev, A. B., Kovalev, Y. Y., et al. 2021, Monthly Notices of the Royal Astronomical Society, 510, 469–480, 10.1093/mnras/stab3358

    work page doi:10.1093/mnras/stab3358 2021

  15. [15]

    L., & Gergely, L

    Kun, E., Biermann, P. L., & Gergely, L. \'A . 2019, Monthly Notices of the Royal Astronomical Society: Letters, 483, L42–L46, 10.1093/mnrasl/sly216

    work page doi:10.1093/mnrasl/sly216 2019

  16. [16]

    \'E ., et al

    Kun , E., Frey , S., Gab \'a nyi , K. \'E ., et al. 2015, , 454, 1290, 10.1093/mnras/stv2049

    work page doi:10.1093/mnras/stv2049 2015

  17. [17]

    \'E ., Karouzos , M., Britzen , S., & Gergely , L

    Kun , E., Gab \'a nyi , K. \'E ., Karouzos , M., Britzen , S., & Gergely , L. \'A . 2014, , 445, 1370, 10.1093/mnras/stu1813

    work page doi:10.1093/mnras/stu1813 2014

  18. [18]

    2023, Symmetry, 15, 270, 10.3390/sym15020270

    Kun , E., & Medveczky , A. 2023, Symmetry, 15, 270, 10.3390/sym15020270

    work page doi:10.3390/sym15020270 2023

  19. [19]

    L., Aller, M

    Lister , M. L., Aller , M. F., Aller , H. D., et al. 2018, , 234, 12, 10.3847/1538-4365/aa9c44

    work page doi:10.3847/1538-4365/aa9c44 2018

  20. [20]

    T., & Bai , J

    Liu , H. T., & Bai , J. M. 2006, , 653, 1089, 10.1086/509097

    work page doi:10.1086/509097 2006

  21. [21]

    Murase, K., & Stecker, F. W. 2023, High-Energy Neutrinos from Active Galactic Nuclei (WORLD SCIENTIFIC), 10.1142/9789811282645_0010

    work page doi:10.1142/9789811282645_0010 2023

  22. [22]

    2022, Monthly Notices of the Royal Astronomical Society, 511, 4697–4701, 10.1093/mnras/stac376

    Padovani, P., Boccardi, B., Falomo, R., & Giommi, P. 2022, Monthly Notices of the Royal Astronomical Society, 511, 4697–4701, 10.1093/mnras/stac376

    work page doi:10.1093/mnras/stac376 2022

  23. [23]

    , author Oikonomou, F

    Padovani, P., Oikonomou, F., Petropoulou, M., Giommi, P., & Resconi, E. 2019, Monthly Notices of the Royal Astronomical Society: Letters, 484, L104–L108, 10.1093/mnrasl/slz011

    work page doi:10.1093/mnrasl/slz011 2019

  24. [24]

    2018, The Astrophysical Journal Letters, 854, L32, 10.3847/2041-8213/aaad5e

    Paiano, S., Falomo, R., Treves, A., & Scarpa, R. 2018, The Astrophysical Journal Letters, 854, L32, 10.3847/2041-8213/aaad5e

    work page doi:10.3847/2041-8213/aaad5e 2018

  25. [25]

    2014, Monthly Notices of the Royal Astronomical Society, 447, 36–48, 10.1093/mnras/stu2364

    Petropoulou, M., & Mastichiadis, A. 2014, Monthly Notices of the Royal Astronomical Society, 447, 36–48, 10.1093/mnras/stu2364

    work page doi:10.1093/mnras/stu2364 2014

  26. [26]

    and Shibahashi , H

    Potter, W. J., & Cotter, G. 2012, Monthly Notices of the Royal Astronomical Society, 423, 756–765, 10.1111/j.1365-2966.2012.20918.x

    work page doi:10.1111/j.1365-2966.2012.20918.x 2012

  27. [27]

    2013, Monthly Notices of the Royal Astronomical Society, 431, 1840–1852, 10.1093/mnras/stt300

    ---. 2013, Monthly Notices of the Royal Astronomical Society, 431, 1840–1852, 10.1093/mnras/stt300

    work page doi:10.1093/mnras/stt300 2013

  28. [28]

    2015, Monthly Notices of the Royal Astronomical Society, 453, 4071–4089, 10.1093/mnras/stv1657

    ---. 2015, Monthly Notices of the Royal Astronomical Society, 453, 4071–4089, 10.1093/mnras/stv1657

    work page doi:10.1093/mnras/stv1657 2015

  29. [29]

    I., Medina-Carrillo, B., et al

    Sahu, S., Páez-Sánchez, D. I., Medina-Carrillo, B., et al. 2024, Monthly Notices of the Royal Astronomical Society, 533, 2156–2161, 10.1093/mnras/stae1847

    work page doi:10.1093/mnras/stae1847 2024

  30. [30]

    Shepherd , M. C. 1997, in Astronomical Society of the Pacific Conference Series, Vol. 125, Astronomical Data Analysis Software and Systems VI, ed. G. Hunt & H. Payne , 77

    work page 1997

  31. [31]

    V., Mufakharov, T

    Sotnikova, Y. V., Mufakharov, T. V., Mingaliev, M. G., et al. 2022, Astrophysical Bulletin, 77, 361–371, 10.1134/s1990341322040149

    work page doi:10.1134/s1990341322040149 2022

  32. [32]

    T., Buson , S., & Kocevski , D

    Tanaka , Y. T., Buson , S., & Kocevski , D. 2017, The Astronomer's Telegram, 10791, 1

    work page 2017

  33. [33]

    2022, Physical Review D, 105, 10.1103/physrevd.105.023005

    Wang, Z.-R., Liu, R.-Y., Petropoulou, M., et al. 2022, Physical Review D, 105, 10.1103/physrevd.105.023005

    work page doi:10.1103/physrevd.105.023005 2022