Pith. sign in

REVIEW 4 major objections 4 minor 1 cited by

Cosmological Evolution of Gamma Ray Bursts

T0 review · 4 major / 4 minor · reviewed 2026-08-04 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper claims that long gamma-ray bursts track the cosmic star formation rate above redshift 1.5, but exceed it by up to two orders of magnitude at low redshift in the spectroscopic sample, with machine-learning-estimated redshifts dilu

desk verdict Competent EP/C-minus analysis of a larger GRB sample; the low-z excess in the spectroscopic sample is solid, but the ML-driven dilution is not selection-corrected and should be treated as tentative. read the letter →

arxiv 2510.23945 v4 pith:OX7BMY2Y submitted 2025-10-27 astro-ph.HE

Sujay Champati , Vahé Petrosian , Maria G. Dainotti This is my paper · ORCID
classification astro-ph.HE PACS 98.70.Rz
keywords gamma-rayburstslongGRBsstarformationrateluminosityevolutionEfron-PetrosianmethodLynden-BellC-redshiftestimationmachinelearning
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper asks whether long gamma-ray bursts (LGRBs) really form only from collapsing massive stars, tracking the cosmic star formation rate (SFR), or whether an excess of low-redshift LGRBs points to an additional progenitor channel. Using the largest sample to date—695 bursts, merging spectroscopic redshifts with 251 machine-learning-estimated redshifts—and applying non-parametric Efron-Petrosian and Lynden-Bell C- methods to correct for selection bias, the authors find that the LGRB formation rate closely tracks the SFR for z≥1.5. Below z=1.5 the spectroscopically confirmed sample shows a deviation that grows to two orders of magnitude at the lowest redshifts, confirming earlier claims. Adding ML-estimated redshifts, which cluster at 1.5

What carries the argument

The machinery is the Efron-Petrosian (EP) method, a non-parametric test for correlation between luminosity and redshift under one-sided truncation, which the authors use to measure luminosity evolution via a broken power-law g(Z) that renders a transformed luminosity L0 independent of redshift. Once the evolution is removed, the Lynden-Bell C- method reconstructs the local luminosity function and the cumulative formation rate from the rank structure of the truncated sample. The flux limit and averaged K-correction define the truncation boundary that these methods require.

What would settle it

Obtain spectroscopic redshifts for a large, complete sample of satellite-detected long GRBs at z<1 (for instance from a future wide-field mission with immediate IR follow-up) and compare their space density directly to the star formation rate at those redshifts; a measured deviation of two orders of magnitude would confirm the claim, while a rate consistent with the SFR would refute it. Alternatively, show that the ML redshift predictor systematically assigns low-z sources to z>1.5 because of its training-set selection, which would invalidate the combined-catalog comparison.

Watch

Extended reading notes

Core claim

The central claim, stated in Section 5, is that the formation rate of long GRBs closely tracks the cosmic star formation rate for z≥1.5 in both the combined catalog and the spectroscopic-only catalog. In the spectroscopic catalog, the density rate deviates increasingly from the SFR below z<1.5, reaching two orders of magnitude at the lowest redshift—confirming earlier results. When machine-learning-estimated redshifts are added, the concordance with SFR persists until z=1, where the formation rate breaks and increases by a factor of ten; the smaller rise is attributed to the absence of low-redshift bursts in the ML sample. The paper interprets the low-z excess as evidence that a significant

Load-bearing premise

The paper's combined-catalog conclusion that the low-redshift excess is reduced rests on treating machine-learning-estimated redshifts as valid measurements on par with spectroscopic ones; if those estimates are biased toward intermediate redshifts by their training set, the reduction is an artifact.

Editorial extensions

If this is right

  • If the low-redshift excess is real, long GRBs cannot be exclusively collapsars at low z; a merger channel would increase the predicted rate of gravitational-wave sources and could explain the kilonova associations seen in GRB 211211A and GRB 230307A.
  • The formation rate tracking the SFR at z≥1.5 supports the collapsar model for the high-redshift LGRB population, where lower metallicity favors envelope retention and jet production.
  • The dilution of the excess when ML redshifts are added implies that redshift-estimation methods with selection functions concentrated at intermediate z can mask genuine low-z features; future population studies should treat ML and spectroscopic samples carefully.
  • The broken power-law luminosity evolution index (k≈2.8 for the full catalog, 3.7 for the spectroscopic) is consistent with earlier estimates, validating the non-parametric approach on a larger sample.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the low-z excess is treated as astrophysical, a direct cross-check is to compare the implied local rate of compact mergers from LGRBs with the LIGO/Virgo binary neutron star merger rate; an inconsistency would force a rethink of the progenitor interpretation.
  • The Anderson-Darling test reported (p=0.001) shows the ML and spectroscopic redshift distributions are not drawn from the same parent population. A joint-likelihood analysis that models both selection functions simultaneously would be a natural next step, rather than pooling the samples as done here.
  • The same EP + C- pipeline could be applied to short GRBs with increasing redshift samples to test whether their formation rate follows a delayed SFR, providing an independent handle on merger delay-time distributions.
  • If the ML redshift distribution's concentration at 1.5<z<3 reflects training-set bias rather than a true dearth of low-z bursts, then the combined-catalog reduction is an artifact; this could be tested by constructing an ML estimator trained on a redshift-complete sample and checking where the predicted redshifts fall.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper analyzes the cosmological evolution of long GRBs using a sample of Swift LGRBs with spectroscopic redshifts (444 after cuts) augmented by 251 ML-estimated redshifts from Dainotti et al. (2025), for a working sample of 499 bursts with fluxes and spectral fits. The authors apply the nonparametric Efron–Petrosian method to remove luminosity evolution and the Lynden-Bell C^- method to derive luminosity functions and formation rates for both the full (spectroscopic+ML) and the spectroscopic-only catalogs. They find that the spectroscopic sample shows a low-redshift excess of the formation rate relative to the Madau–Dickinson SFR, growing to about two orders of magnitude at the lowest redshifts, while the full catalog shows a smaller rise because the ML sample contributes few low-z bursts. The main claimed result is that the LGRB formation rate closely tracks the SFR for z>=1.5 in both catalogs, with the low-z excess confirming earlier work.

Significance. If the ML-augmented comparison were robust, this would be a valuable step forward in using larger, less incomplete GRB samples for population studies. The paper provides a detailed description of the data pipeline, flux limits, K-corrections, and spectral fits, and it explicitly documents the distributional difference between ML and spectroscopic redshifts (Anderson–Darling p=0.001), which is commendable transparency. The spectroscopic-only analysis reproduces and strengthens previously reported low-z LGRB excess with a larger sample. However, the novel full-catalog comparison is not yet convincing because the ML sample is not corrected for its z-dependent inclusion probability, and because the analysis lacks error propagation and uncertainty estimates on key outputs.

major comments (4)
  1. [Section 2, Fig. 1, Section 5]
  2. [Section 4.2, 4.3, Tables 1-2, Fig. 7]
  3. [Section 3, Eqs. (6), (10), (11)]
  4. [Section 2, Fig. 2]
minor comments (4)
  1. [Section 2]
  2. [Section 4.1]
  3. [Fig. 6 caption]
  4. [General]

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the EP/C^- derivation is self-contained, and the low-z excess is an empirical comparison against the external Madau-Dickinson SFR, not an identity.

full rationale

The paper's central derivation chain is not circular. The luminosity-evolution slope k is fitted within the Efron-Petrosian procedure to make L0 = L/g(Z) uncorrelated with Z (Eqs. 6-7, Fig. 5), and this fitted k is a nuisance parameter, not the claimed result. The subsequent Lynden-Bell C^- reconstruction of the luminosity function and formation rate (Eqs. 8-14) is an independent nonparametric inversion of the de-evolved data. The headline claim is a comparison of the resulting formation-rate shape with the externally compiled Madau & Dickinson (2014) SFR; the GRB rate is not normalized to, or defined in terms of, the SFR curve. Thus the low-z excess is an empirical contrast, not an identity. The authors' use of prior same-group work (Petrosian et al. 2015 for Zcr ~ 3.5; Petrosian & Dainotti 2024 for earlier excess claims; Narendra et al. 2025 / Dainotti et al. 2025 for ML redshifts) is contextual or supplies input data, but does not by itself force the conclusions: the spectroscopic-only excess is obtained from the same pipeline without ML redshifts, and the reduction of the excess in the combined catalog is explicitly attributed to the absence of low-z ML bursts (Sec. 5), which is a selection caveat rather than a derivation of the result from the input. The Anderson-Darling p = 0.001 between the ML and spectroscopic redshift distributions is disclosed, and the difference is attributed to training-set selection, which is a correctness risk, not circularity. One minor self-citation issue is that Zcr is imported from prior same-author work rather than independently re-fit, but because k is fit and the conclusions are robust to the reported k ranges, this is not load-bearing.

Assumptions & free parameters 7 free parameters · 6 assumptions · 0 invented entities

The central claim rests on standard nonparametric truncation statistics plus two fragile data assumptions: the truncation boundary is unbiased, and the ML-redshift sample can be pooled with the spectroscopic sample. No new physical entities are introduced. The model parameters in Tables 1 and 2 are fits to the derived distributions, not independent measurements.

free parameters (7)
  • luminosity evolution index k = Full: 2.8 (+2.3/-3.3); non-ML: 3.7 (+2.9/-4.2)
    Chosen so Kendall's tau equals zero in Eq. 6; defines the de-evolved luminosity L0 = L/g(Z), so all downstream LF and formation-rate results depend on it.
  • break redshift Zcr in g(Z) = 3.5
    Adopted from Petrosian et al. (2015), not refit here; sets where luminosity evolution flattens.
  • flux limit f_lim = 4e-8 erg s^-1 cm^-2 keV^-1
    Conservative truncation limit used to construct associated sets; changing it alters the bias correction and the derived rates.
  • K-correction power-law parameters = amplitude 0.93, slope 0.69
    Fit to the windowed running average of K values and used in Lmin(Z); affects the truncation boundary for every burst.
  • moving window size for K average = 10 bursts
    Smoothing parameter chosen by hand for constructing Kbar(z).
  • luminosity function fit parameters = Full: phi0=308, L0=3.56e50, delta1=0.27, delta2=1.35; non-ML: phi0=172, L0=1.56e50, delta1=0.31, delta2=1.22
    Fitted to the cumulative luminosity function from Eq. 10; reported without uncertainties.
  • formation rate double-break parameters = Full: alpha=16.76, Zc1=1.20, beta1=4.77, Zc2=2.98, beta2=0.09, N0=1.94; non-ML: alpha=11.70, Zc1=1.39, beta1=3.57, Zc2=2
    Fitted to the cumulative rate from Eq. 11 and used for the comparison to SFR in Figure 7.
assumptions (6)
  • domain assumption After correcting luminosity by g(Z), L0 and Z are statistically independent, so the bivariate distribution factorizes as Psi(L,Z) = psi(L) rho(Z).
    Invoked in Section 3 to justify the Lynden-Bell C-minus method; residual correlation would bias the luminosity function and formation rate.
  • domain assumption The chosen f_lim yields a complete, unbiased truncation boundary that is more conservative than the survey limit.
    Section 2 assumes this removes selection bias, but completeness is not independently verified.
  • domain assumption ML-estimated redshifts from Dainotti et al. (2025) / Narendra et al. (2025) are accurate enough to pool with spectroscopic redshifts in population analysis.
    Section 2; the paper's own Anderson-Darling test shows the ML and spectroscopic distributions differ (p=0.001), making this assumption fragile and load-bearing for the full-sample result.
  • domain assumption Standard flat Lambda-CDM cosmology with Omega_m = 0.3 and H0 = 70 km/s/Mpc for luminosity distances.
    Used in Eq. 1 to convert flux to luminosity; a different cosmology would shift luminosities and truncation limits.
  • domain assumption Power-law and cutoff-power-law spectral models, selected by the Sakamoto criterion, correctly describe the bursts for K-corrections.
    Section 2; K(z) depends on these spectral fits, and 74 Swift bursts lack reliable model fits.
  • standard math The Efron-Petrosian and Lynden-Bell C-minus estimators are unbiased for one-sided truncation in this setting.
    Section 3; relies on established statistical results from Efron & Petrosian (1992) and Lynden-Bell (1971).

how reviews work

0 comments
Cite this review

Pith. "Pith review of Cosmological Evolution of Gamma Ray Bursts." pith.science (2026). https://pith.science/paper/OX7BMY2Y

@misc{pith2026251023945,
  author       = {Pith},
  title        = {Pith review of: Cosmological Evolution of Gamma Ray Bursts},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OX7BMY2Y}},
  note         = {Machine review of arXiv:2510.23945}
}
abstract

Gamma-ray bursts (GRBs) are classified as long (LGRBs) and short (SGRBs), with collapsars and compact-object mergers (NS-NS or NS-Black Holes) as progenitors, respectively. LGRBs are expected to follow the cosmic star formation rate (SFR), while SGRBs follow a delayed version of the SFR. However, this division has come under question, most prominently by observational evidence of an excess of LGRBs at low redshifts by several investigations, summarized in \cite{Petrosian_2024}. Two recent observations of low-redshift LGRBs show associations with kilonovae. Both of these indicate compact mergers as a potential source of LGRBs as well. Most results showing this separation are based on analyses of small (less than 200) samples of LGRBs with measured redshifts. The aim of this paper is to use a larger sample of LGRBs. The number of LGRBs with measured redshifts has increased by more than a factor of 2 over the last decade. To this data set we add a sample of LGRBs whose redshifts are estimated using a machine learning (ML) method (\cite{Narendra_2025}). To account for the observational selection bias due to redshift measurements, we use the non-parametric, non-binning Efron-Petrosian method to establish the degree of correlation between luminosity and redshift, \textit{the luminosity evolution}, and then use the Lynden-Bell $C^-$ method to obtain the luminosity function. We find a low redshift excess for the larger sample with measured redshifts. Adding the sources with ML-estimated redshifts, which shows overabundance of the mid-range redshifts, the excess is reduced.

Figures

Figures reproduced from arXiv: 2510.23945 by the authors.

Figure 1
Figure 1. Comparison of fractional binned redshift distri￾butions of ML (green) and non-ML (red) samples. The union of these two sets gives us a total of 695 bursts. In [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Flowchart of Data Processing Pipeline 1 2 3 4 5 6 7 8 9 10 1 2 3 4 Z = 1 + z K-correction (Window=10) Power-Law Fit to K-corrections Running Average (Window=10) K(z) = 0.93 (1 + z) 0.69 [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Power-law fit overlaid on running-average K￾corrections. However, as emphasized in Petrosian (1992), these nonparametric methods assume statistical independence between L and Z, i.e., Ψ(L, Z) = ψ(L)ρ(Z), and thus cannot capture luminosity evolution. To address this, Efron and Petrosian (EP) developed a method Efron & Petrosian (1992) that tests for correlations between variables under one-sided truncation.4 We adopt… view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: (Left): Redshifts and Luminosities for combined GRB catalog. The black truncation line is obtained from flim and the average K. The power law fit for the bursts is shown in green. (Right): Redshifts and Luminosities for the Non-ML Catalog de-evolved L0 interpretable as…
Figure 5
Figure 5. Figure 5: Kendall’s τ statistic as a function of k for the full GRB catalog (blue) and the non-ML subset (red). Solid curves show the measured τ (k). Dashed lines indicate the central (fit) guides for each dataset, while the lighter dashed bands denote the corresponding ±1σ rang…
Figure 6
Figure 6. Figure 6: (Left): Formation rates derived from the C − method for both the full (blue) and non-ML (red) GRB catalogs. The solid curves indicate the best-fit double-break power law models. Cumulative luminosity functions for the GRB catalogs. Blue points with a solid fit line rep…
Figure 7
Figure 7. Figure 7: Comparison of the LGRB-derived formation rates with the star formation rate density from Madau & Dickin￾son (2014). Green points with black error bars denote the ob￾servational SFR data, while the dashed blue and red curves correspond to the best-fit double-break model…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Gamma-Ray Bursts: Evidence for a Common Origin of X-ray Plateaus with Diverse Temporal Decay Index

    astro-ph.HE 2026-01 conditional novelty 5.0 of 10

    Gamma-ray burst X-ray plateaus with rising, flat, and decaying slopes come from one statistically uniform population, not distinct subclasses.

Reference graph

Works this paper leans on

3 extracted references · 1 linked inside Pith · cited by 1 Pith paper

  1. [1]

    2014, ARA&A, 52, 43, doi: 10.1146/annurev-astro-081913-035926 Bromberg, O., Nakar, E., & Piran, T

    Berger, E. 2014, ARA&A, 52, 43, doi: 10.1146/annurev-astro-081913-035926 Bromberg, O., Nakar, E., & Piran, T. 2011, ApJL, 739, L55, doi: 10.1088/2041-8205/739/2/L55 Cucchiara, A., Levan, A. J., Fox, D. B., et al. 2011, The Astrophysical Journal, 736, 7, doi: 10.1088/0004-637X/736/1/7 Dainotti, M. G., Bhardwaj, S., Cook, C., et al. 2025, GRB Redshift Class...

    arXiv 2014

  2. [163]

    http://www.jstor.org/stable/2241151 Woosley, S. E. 1993, ApJ, 405, 273, doi: 10.1086/172359 Woosley, S. E., & Bloom, J. S. 2006, ARA&A, 44, 507, doi: 10.1146/annurev.astro.43.072103.150558 Yang, Y.-H., Troja, E., O’Connor, B., et al. 2024, Nature, 626, 742, doi: 10.1038/s41586-023-06979-5 9 Yu, Y.-W., Gao, H., Wang, F.-Y., & Zhang, B.-B. 2022, Gamma-Ray B...

    arXiv 1993

  3. [1597]

    1985, The Annals of Statistics, 13,

    http://www.jstor.org/stable/2241492 Woodroofe, M. 1985, The Annals of Statistics, 13,

    arXiv 1985

Pith tools

Reviewed August 4, 2026 · model on record in the stance chip above.