🔭 astrophysics

Probing Evolution of Long Gamma-Ray Burst Properties through Their Cosmic Formation History

This study analyzes Swift and ground-based data to demonstrate that the cosmic Long Gamma-Ray Burst rate density cannot be fully explained by the star formation rate density through single-property evolution models, indicating that a combination of multiple evolving LGRB properties is required to reconcile the two rate densities across cosmic history.

Original authors: Nikita S. Khatiya, Maria Giovanna Dainotti, Aditya Narendra, Dhruv S. Bal, Aleksander Ł. Lenart, Dieter H. Hartmann

Published 2026-05-13

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Nikita S. Khatiya, Maria Giovanna Dainotti, Aditya Narendra, Dhruv S. Bal, Aleksander Ł. Lenart, Dieter H. Hartmann

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Cosmic Detective Story: Why Long Gamma-Ray Bursts Don't Match the Star-Formation Map

Imagine the universe as a giant, bustling construction site. For decades, astronomers have been trying to figure out exactly how many new "buildings" (stars) are being built at any given time. They call this the Star Formation Rate Density (SFRD).

Then, there are these spectacular, blinding flashes of light called Long Gamma-Ray Bursts (LGRBs). Think of these as the universe's most dramatic fireworks displays, which we believe happen when massive stars collapse and die.

For a long time, scientists thought these fireworks were a perfect, 1-to-1 map of the construction site. The logic was simple: "More stars being born = more massive stars dying = more fireworks."

But this new paper by Khatiya and colleagues says: "Not so fast."

Here is what they found, explained simply:

1. The "Fireworks" Don't Match the "Blueprint"

The researchers gathered data on 242 of these cosmic fireworks (specifically those with a unique "X-ray plateau" feature, which makes them a very consistent group to study). They tried to plot where and when these explosions happened against the map of where stars are being born.

The Result: The maps didn't line up.

In the distant past (High Redshift): There were more fireworks than the star-formation map predicted.
In the recent past (Low Redshift): There were also way more fireworks than expected, creating a strange "bump" in the data that shouldn't be there.

It's like trying to predict how many car accidents happen on a highway by counting how many cars are driving. If you see a massive pile-up in a quiet neighborhood (low redshift) or a sudden spike in a desert (high redshift) that doesn't match the traffic flow, you know something else is going on.

2. The "No Evolution" Theory is Dead

The team tested a simple idea: What if the fireworks have always been exactly the same, and we just need to count them better? They called this the "no evolution" case.

The Verdict: This idea was ruled out immediately. The data simply doesn't fit. The fireworks aren't static; they are changing as the universe ages.

3. Trying to Fix the Map with "One-Size-Fits-All" Adjustments

The scientists tried to fix the mismatch by applying simple rules, like:

The "Beam" Theory: Maybe the fireworks are just pointing at us more often now than before? (Changing the "beaming angle").
The "Power Law" Theory: Maybe the brightness of the fireworks changes in a simple, predictable mathematical way over time?

The Verdict: Neither of these simple fixes worked.

Adjusting the "beam" helped explain the distant past but made the recent past look even worse.
A simple mathematical curve could fix the recent past but broke the distant past.

The Analogy: Imagine trying to fit a square peg into a round hole. You try sanding off one corner (fixing the high-z data), but then the other side doesn't fit (the low-z data breaks). You try sanding the other side, and now the first side is loose. You can't fix the shape with just one simple cut.

4. The Real Conclusion: It's a Complex Recipe

The paper concludes that to make the "Fireworks Map" match the "Star Formation Map," you can't just tweak one knob. You need to change multiple properties at the same time.

The universe isn't just turning a dial up or down. The properties of the stars that cause these explosions (like their metal content, how they spin, or how they are grouped in binary systems) are all evolving together in a complex dance.

Why the mismatch?
- At the beginning of the universe: Stars were made of "pure" ingredients (low metallicity), which might make them more likely to explode as fireworks.
- In the recent universe: There might be a hidden population of "low-luminosity" fireworks (fainter ones) that we are accidentally counting as the big ones, or perhaps some of these explosions are caused by a different mechanism entirely (like merging black holes) rather than a dying star.

The Bottom Line

The authors are essentially saying: "We have a great map of where stars are born, and we have a great list of where the big explosions happen. But they don't match because the explosions are more complicated than we thought. We can't explain the difference with a single simple rule; we need to understand a whole new set of evolving rules that govern how these stars die."

They didn't find a new application for this data (like medical uses or new technology); they simply proved that our current understanding of the "recipe" for these cosmic explosions is incomplete and needs a much more sophisticated cookbook.

Technical Summary: Probing Evolution of Long GRB Properties through Their Cosmic Formation History

Problem Statement
The astrophysics of Long Gamma-Ray Burst (LGRB) progenitors and the potential cosmological evolution of their properties remain subjects of significant debate. While previous studies suggest that the LGRB rate density (LGRB-RD) tracks the cosmic star formation rate density (SFRD) at high redshifts ( $z$ ), attributing this to metallicity evolution, the relationship at low redshifts is contentious. Some studies indicate an "uptick" or excess in LGRB-RD at low $z$ compared to SFRD, potentially caused by a distinct class of low-luminosity GRBs (LLGRBs), sample contamination, or different progenitor subclasses. Furthermore, existing literature offers conflicting views on whether LGRBs are unbiased tracers of star formation or biased tracers influenced by metallicity and the initial mass function (IMF). This study aims to resolve these discrepancies by investigating the LGRB-RD and its mapping to the SFRD using a homogeneous sample of LGRBs with specific X-ray plateau features.

Methodology
The authors utilize a sample of 242 LGRBs detected between March 2005 and December 2023, characterized by X-ray plateaus in their afterglow lightcurves observed by the Neil Gehrels Swift Observatory. The selection of this specific subsample is motivated by the homogeneity of the data and reduced scatter in correlations compared to prompt emission studies.

The analysis proceeds through the following steps:

Luminosity and Flux Limits: Observed luminosities at the end of the plateau ( $L_{a,obs}$ ) and minimum observable luminosities ( $L_{min,obs}$ ) are calculated. A flux limit cut ( $\log F_{lim} = -11.9$ ) is applied, removing 17 bursts to mitigate Malmquist bias.
Non-Parametric Correction: To address intrinsic luminosity evolution and selection biases, the study employs the non-parametric Efron & Petrosian (EP) method. This technique tests the independence of luminosity ( $L$ ) and redshift ( $1+z$ ) using the Kendall $\tau$ statistic.
Evolution Factor Determination: The analysis determines a luminosity evolution factor, $(1+z)^k$ , where $k = 2.1^{+0.5}_{-0.7}$ . This factor is used to correct the cumulative rate density ( $\sigma$ ) and luminosity function ( $\Phi$ ), yielding corrected distributions ( $\sigma_{corr}$ and $\Phi_{corr}$ ).
Rate Density Calculation: The differential LGRB-RD ( $\rho(z)$ ) is derived from the corrected cumulative distribution. The authors apply correction factors ( $F_1$ and $F_2$ ) to account for excluded bursts (due to flux limits or zero associated sets) and instrumental limitations (Swift-BAT field of view, uptime, and operation years) to estimate the true LGRB-RD.
Theoretical Comparison: The corrected LGRB-RD ( $\rho_{corr}$ ) is compared against a theoretical LGRB-RD ( $\rho_{theor}$ ) derived from the SFRD (Madau & Dickinson 2014) using a "Drake equation-like" approach. This theoretical model incorporates factors for the fraction of massive stars, core-collapse supernovae, Type Ic supernovae, beaming, and metallicity.

Key Contributions

Updated Luminosity Function: The study updates the LGRB luminosity function (LF) for LGRBs with X-ray plateaus using a more recent and larger sample than previous works (e.g., Dainotti et al. 2013).
First RD Analysis for X-ray Plateau Sample: This is the first study to derive the LGRB-RD specifically for the X-ray plateau sample and compare it directly with SFRD models.
Evolutionary Constraints: The paper provides a quantitative constraint on the luminosity evolution factor ( $k \approx 2.1$ ) required to make the luminosity and redshift distributions independent, correcting for Malmquist bias.
Comparative Analysis: The work compares the derived $\rho_{corr}$ with both theoretical models and previous observational results (including an optical plateau sample from Dainotti et al. 2024) to isolate selection effects and evolutionary trends.

Results

Luminosity Evolution: The analysis confirms significant luminosity evolution, with the best-fit evolution factor being $(1+z)^{2.1}$ . Without this correction, the intrinsic nature of the LGRB-RD is obscured.
Rate Density Behavior: The corrected LGRB-RD ( $\rho_{corr}$ ) exhibits a relatively flat distribution across most redshifts with a subtle peak around $(1+z) \approx 2.5$ .
Mismatch with SFRD:
- Low- $z$ Excess: The study confirms an excess of LGRBs at low redshifts ( $z < 1$ ) compared to the SFRD, consistent with findings by Yu et al. (2015), Petrosian et al. (2015), and others. This excess diminishes at $z \sim 0.1$ , likely due to small number statistics.
- High- $z$ Discrepancy: At high redshifts, the $\rho_{corr}$ shows a mismatch with most SFRD models, though it partially aligns with the Hopkins & Beacom (2006) model. The SFRD peak occurs at a slightly higher redshift than the LGRB-RD peak.
Testing Evolution Hypotheses: The authors test three scenarios to map $\rho_{theor}$ $ρ_{t h eor}$ to $\rho_{corr}$ $ρ_{cor r}$ :
1. No Evolution: Ruled out as the distributions are incompatible.
2. Beaming Angle Evolution: A factor of $(1+z)^{1.6}$ flattens the high- $z$ slope of the theoretical curve but fails to explain the low- $z$ behavior.
3. Simple Power-Law Evolution: A single power-law factor $(1+z)^\delta$ cannot simultaneously match both the low- $z$ and high- $z$ ends of the observed rate density. Negative $\delta$ values fit the low- $z$ end, while positive values fit the high- $z$ end, but neither works for the entire range.

Significance and Claims
The paper concludes that the "no evolution" hypothesis is insufficient to explain the relationship between LGRB formation and star formation. Furthermore, simple evolutionary models involving only the beaming angle or a single power-law evolution of properties are inadequate to map the LGRB-RD to the SFRD across the entire redshift range.

The authors claim that matching the LGRB-RD and SFRD in their entirety likely requires the inclusion of multiple evolving properties of LGRBs simultaneously. The observed discrepancies suggest that selection biases (such as LLGRB contamination or merger-driven GRBs misclassified as LGRBs) and complex evolutionary factors (such as metallicity dependence, evolving IMFs, or progenitor system evolution) must be accounted for. The study emphasizes that while the derived $\rho_{corr}(z)$ is a powerful tool for probing progenitor systems and stellar mass black hole populations, a more elaborate redshift evolution function is necessary to fully understand the cosmic formation history of LGRBs. The authors modestly suggest that future data, potentially aided by machine learning for redshift inference, will be required to further address these complex issues.