Extreme-Value Distribution Analysis of the Second CHIME/FRB Catalog: Assessing the Rarity of the One-off FRB 20250316A

We present a statistical analysis of the extreme brightness of the fast radio burst FRB 20250316A, a luminous, apparently non-repeating event detected by CHIME/FRB. Employing a model-agnostic framework based on the Generalized Extreme Value (GEV) distribution applied to the second CHIME/FRB catalog, we quantify its rarity within the current population. Bayesian fitting of GEV models to block maxima of peak flux and fluence data reveals FRB 20250316A to be a profound statistical outlier. For the peak flux, the analysis yields a return period of $\sim$ 600 years ($1σ$ credible level), with the underlying distribution being of the heavy-tailed, unbounded Fréchet type ($ξ> 0$). The fluence analysis indicates greater complexity: while the full sample suggests a Fréchet-type distribution with a $\sim50$-year return period in $1σ$ credible level, the removal of three other notable outliers points toward a light-tailed Weibull-type distribution ($ξ< 0$) with a finite upper bound far exceeded by FRB 20250316A. This dichotomy highlights the challenge in characterizing the tail of the FRB luminosity function with limited data. Although less extreme in recurrence time than the Brightest Of All Time'' gamma-ray burst GRB 221009A, FRB 20250316A constitutes a similarly exceptional event (a potential FRB BOAT’’) within the short observational history of wide-field radio surveys. Our results underscore the existence of rare, highly luminous events that may probe the upper limits or distinct sub-populations of the FRB luminosity distribution.

💡 Research Summary

The paper presents a rigorous statistical investigation of the exceptionally bright, apparently non‑repeating fast radio burst FRB 20250316A, discovered by CHIME/FRB in March 2025. Using the second CHIME/FRB catalog (spanning roughly five years of observations), the authors apply a model‑agnostic extreme‑value framework based on the Generalized Extreme Value (GEV) distribution to quantify how rare such an event is within the current FRB population.

First, the catalog of non‑repeating FRBs is divided into contiguous 30‑day blocks. Within each block the maximum peak flux and the maximum fluence are extracted, yielding two independent series of block‑maxima. The choice of a 30‑day block balances the need for enough events per block (to approach the asymptotic GEV limit) against the desire for a sufficient number of block maxima for reliable inference.

The GEV distribution is characterized by three parameters: location (µ), scale (σ), and shape (ξ). The shape parameter determines the tail class: ξ > 0 corresponds to a Fréchet (heavy‑tailed, unbounded) distribution, ξ = 0 to a Gumbel (exponential) distribution, and ξ < 0 to a Weibull (light‑tailed, finite upper bound). A Bayesian inference scheme is employed, with uniform priors on µ and σ and a weakly informative normal prior (mean 0, σ = 0.5) on ξ. Posterior sampling is performed using the affine‑invariant MCMC ensemble sampler from the Python package emcee, running 200 000 steps with 100 walkers and confirming convergence (Gelman‑Rubin R̂ ≈ 1.01).

For the peak‑flux series, the posterior mean of ξ is 0.32 ± 0.14, clearly indicating a Fréchet‑type tail. The fitted GEV model predicts a 1‑σ (68 % credible) return period of ≈ 593 years for a burst as bright as FRB 20250316A, and a 2‑σ (≈ 95 % credible) return period of about 73 years. This places the event far beyond the typical extremes observed in the five‑year CHIME dataset, confirming its status as a statistical outlier. QQ‑plots show that the bulk of the block‑maxima lie within the 95 % confidence band, supporting the adequacy of the model for the peak‑flux data.

The fluence analysis is more nuanced. When the full sample—including three other high‑fluence outliers (FRB 20200723B, 20220222B, 20210922C)—is used, the posterior shape parameter is ξ = 0.37 ± 0.18, again suggesting a heavy‑tailed Fréchet distribution. The corresponding return periods are ≈ 14 years (1‑σ) and ≈ 51 years (2‑σ). However, the QQ‑plot reveals systematic deviations beyond the 95 % band, indicating that the model does not fully capture the data’s tail.

To test robustness, the authors repeat the fluence fit after removing the three additional outliers, treating them as distinct rare events. In this reduced sample the posterior mean of ξ becomes –0.24 ± 0.12, pointing to a Weibull‑type distribution with a finite upper bound estimated at µ – σ/ξ ≈ 174.5 Jy ms. FRB 20250316A’s fluence of 1.7 kJy ms exceeds this bound by roughly an order of magnitude, making it an extreme outlier even under a light‑tailed model. Notably, the posterior for ξ still spans both positive and negative values, reflecting the limited statistical power of the current catalog to definitively determine the tail class.

The authors conclude that FRB 20250316A is a profound outlier in both peak flux and fluence. Its peak‑flux behavior aligns with a heavy‑tailed Fréchet distribution, implying a return period of several centuries, while the fluence tail is ambiguous: it could be heavy‑tailed (return period of a few decades) or light‑tailed with a hard upper limit far surpassed by this burst. This dichotomy mirrors findings for the “Brightest Of All Time” gamma‑ray burst GRB 221009A, suggesting that the brightest transients may arise either from the extreme tail of a single population or from distinct sub‑populations with different physical origins.

Methodologically, the work demonstrates the utility of GEV extreme‑value analysis for radio transients, while also highlighting challenges: the choice of block length, the small number of extreme events, and the sensitivity of the shape parameter to a handful of outliers. The authors advocate for larger, longer‑duration FRB catalogs from CHIME, the upcoming Square Kilometre Array, and other wide‑field facilities to refine tail‑parameter estimates, resolve whether a single luminosity function suffices, and ultimately improve our understanding of the most energetic FRB phenomena.