Detection of afterglow emission up to 100 GeV through a stacking analysis of gamma-ray bursts

This paper reports the detection of high-energy gamma-ray emission up to 100 GeV from a stacked sample of 330 gamma-ray bursts using Fermi-LAT data, revealing that while individually detected bursts align with standard afterglow models, individually undetected bursts suggest the presence of a previously unobserved energy injection effect in the GeV band.

The Gist

Imagine trying to hear a single whisper in a hurricane. That is essentially what astronomers face when they try to study the faint, high-energy "echoes" (afterglows) of Gamma-Ray Bursts (GRBs). GRBs are the most violent explosions in the universe, like cosmic supernovas that happen in the blink of an eye. While we can see the initial flash easily, the high-energy "afterglow" that follows is often too faint for our telescopes to catch individually.

This paper, written by a team of scientists, describes a clever trick they used to solve this problem: The "Crowd-Sourcing" of Light.

Here is the story of their discovery, broken down into simple concepts.

1. The Problem: Too Faint to See

Think of a GRB as a firework. The initial explosion is bright and easy to see. But the lingering smoke and sparks (the afterglow) fade quickly. For most GRBs, the high-energy gamma rays are so dim that a single telescope (Fermi-LAT) sees nothing but static noise. It's like trying to hear a single person whispering in a crowded stadium; you can't pick out their voice.

2. The Solution: The "Stacking" Technique

Instead of listening to one whisper, the scientists decided to listen to 330 whispers at the same time.

They took data from 330 different GRBs and "stacked" them on top of each other. Imagine taking 330 blurry, low-quality photos of a faint star and layering them perfectly on top of one another. Individually, the photos are useless. But stacked together, the noise cancels out, and the star suddenly becomes bright and clear.

By doing this, they could detect high-energy gamma rays (up to 100 GeV) that were previously invisible. It was like turning up the volume on a choir of faint singers until their combined voice became a roar.

3. The Discovery: Two Different Types of Fireworks

Once they had this clear signal, they split the data into two groups and found something surprising:

• Group A (The Bright Stars): These were the 220 GRBs that the telescope could already see individually. Their "afterglow" behaved exactly as scientists expected. It was like a standard firework: it flared up, faded slowly, and then died out quickly. The physics behind this was the "standard model" of how jets of particles crash into space dust.
• Group B (The Hidden Stars): These were the 110 faint GRBs that the telescope couldn't see on its own. When the scientists stacked these, they saw something weird. The light didn't just fade away; it seemed to get a "second wind."

4. The Mystery: The "Energy Injection"

For the faint group (Group B), the light curve (the graph of brightness over time) showed a strange behavior. After fading for a while, the light stayed bright for much longer than expected before finally dropping off.

The scientists compared this to a car engine.

• Standard GRBs: Like a car running out of gas. It speeds up, coasts for a bit, and then slowly rolls to a stop.
• Faint GRBs: Like a car that seems to run out of gas, but then someone secretly pours a fresh tank of fuel into the engine, making it speed up again before finally stopping.

This suggests that for these weaker explosions, the central engine (the black hole or neutron star left behind) is injecting extra energy into the jet long after the initial explosion. This is a new discovery for high-energy gamma rays, confirming that the "engine" of these cosmic explosions can keep running longer than we thought.

5. Why It Matters

This study is a big deal for three reasons:

1. It proves the technique works: We can now study the "faint end" of the universe by stacking data, opening up a whole new world of research.
2. It reveals the engine's secrets: The discovery of "energy injection" in the faint GRBs tells us that the central engines of these explosions are more complex and active than we realized.
3. It maps the physics: By seeing how the light changes from soft to hard (from low energy to high energy), they confirmed that these explosions involve two types of particle acceleration working together, like a cosmic particle accelerator.

The Bottom Line

The universe is full of faint, high-energy signals that are too quiet for our ears to hear alone. By using a "stacking" technique—listening to hundreds of them at once—this team of scientists turned a whisper into a shout. They discovered that while some cosmic explosions fade away normally, others have a hidden "second wind," proving that the engines driving the universe's most violent events are still firing long after the initial blast.

Technical Summary

Here is a detailed technical summary of the paper "Detection of afterglow emission up to 100 GeV through a stacking analysis of gamma-ray bursts."

1. Problem Statement

Gamma-ray bursts (GRBs) are extreme cosmic explosions where high-energy gamma-ray emission ($>1$ GeV) is crucial for understanding jet evolution, particle acceleration, and magnetic field properties. However, detecting this high-energy emission is challenging because:

• Limited Statistics: Space-based detectors like Fermi-LAT have limited effective areas, resulting in very few high-energy photons for most individual GRBs.
• Observation Constraints: Ground-based experiments suffer from small fields of view or high backgrounds.
• Data Gaps: While a few exceptionally bright GRBs (e.g., GRB 130427A, GRB 221009A) have been well-studied, the vast majority of GRBs lack sufficient photon counts to derive precise spectral and temporal evolution of their afterglows.
• Prompt vs. Afterglow: Distinguishing the prompt emission from the afterglow in the GeV band is difficult for individual events due to low statistics.

The authors aim to overcome these limitations by employing a stacking analysis on a large population of GRBs to detect and characterize high-energy afterglow emission up to 100 GeV.

2. Methodology

The study utilizes data from the Fermi Large Area Telescope (LAT) covering the period from August 2008 to December 2024.

• Sample Selection:
  • Sample A (LAT Detected): 220 GRBs individually detected by Fermi-LAT (triggered by LAT or GBM), excluding the two brightest events (130427A and 221009A) to avoid bias.
  • Sample B (LAT Undetected): 110 GRBs triggered by Swift/BAT but not individually detected by Fermi-LAT. A filtering condition required at least one event within a $1^\circ$ radius of the burst location, representing the faintest end of the GRB population accessible to Fermi-LAT.
• Data Processing:
  • Stacking: Photons in the $0.1-100$GeV band were collected from a$10^\circ \times 10^\circ$region centered on each GRB for a duration of 50,000 seconds post-trigger ($T_0$).
  • Background Subtraction: The background model included Galactic diffuse emission, isotropic diffuse background, and sources from the 4FGL catalog.
  • Time Binning: Light curves and Spectral Energy Distributions (SEDs) were analyzed in specific time intervals (e.g., $T_0+[0, 10]$s to $T_0+[10000, 50000]$s).
  • Prompt Emission Removal: To isolate the afterglow, photons within the $T_{90}$ duration (90% of prompt gamma-ray photons) were subtracted for specific analyses.
• Statistical Analysis:
  • Likelihood analysis was performed to calculate the Test Statistic (TS).
  • Significance was determined via pseudo-experiments (shifting time windows by 3 months) to fit a $\chi^2_\nu$ distribution.
• Modeling:
  • Phenomenological: SEDs were fitted with an Exponential Cutoff Power-Law (ECPL) model.
  • Physical: A forward shock afterglow model (synchrotron + Synchrotron Self-Compton/SSC) was used. The model accounted for a population of GRBs with a log-normal distribution of initial Lorentz factors ($\Gamma_0$).
  • Parameters: Sample A was modeled with high $\Gamma_0$ ($130-800$) and high kinetic energy ($E_0 \sim 4 \times 10^{53}$erg); Sample B was modeled with lower$\Gamma_0$($60-250$) and lower$E_0$($4 \times 10^{52}$ erg).

3. Key Results

A. Detection Significance

• Sample A (220 LAT Detected): High-significance detection of excess emission with a significance of $99.7\sigma$.
• Sample B (110 LAT Undetected): Significant detection of excess emission with a significance of **$17.9\sigma$** (in the$1-10$ GeV band).
• Angular Distribution: The spatial distribution of excess photons matches the Fermi-LAT Point Spread Function (PSF), confirming the emission originates from the GRB locations.

B. Temporal Evolution (Light Curves)

• Sample A: The light curves exhibit a three-segment evolution:
  1. Rise: Slope $\sim 0.7$ (coasting phase).
  2. Shallow Decay: Slope $\sim -0.2$ to $-0.3$ (deceleration phase).
  3. Steep Decay: Slope $\sim -1.5$ starting around $t \sim 100$ s (jet break).
  • Note: The shallow decay slope is flatter than typical individual GRBs due to the superposition of GRBs with different $\Gamma_0$.
• Sample B: The light curve shows a shallow decay ($\sim -0.2$ to $-0.3$) followed by a steep decay ($\sim -3.7$) at a much later time ($t \sim 10^4$ s). The rising phase is not detected due to low photon counts.

C. Spectral Evolution

• SEDs: The spectra evolve from soft to hard over time.
• ECPL Fitting: The spectral index ($\alpha$) remains roughly constant at $\sim 1.8$ when prompt emission ($T_{90}$) is removed. The cutoff energy ($E_{cut}$) increases with time.
• Radiation Mechanism: Both synchrotron and SSC components contribute. Synchrotron dominates early, while SSC becomes increasingly important at later times and higher energies, explaining the rising $E_{cut}$.

D. The "Weak" Sample Anomaly

• The standard forward shock model (with jet breaks) successfully reproduces the light curves of Sample A.
• However, for Sample B (LAT undetected), the standard model fails to reproduce the late-time emission. The observed break time ($\sim 10^4$ s) aligns with the duration of X-ray plateaus in canonical GRBs.
• Cross-Check: Analysis of Swift/XRT light curves for Sample B shows a higher fraction of "Canonical" and "One break" light curves (indicative of energy injection) compared to Sample A.
• Conclusion: The weak GeV sample likely exhibits late-time energy injection from the central engine, a feature not prominent in the brighter GeV sample.

4. Key Contributions

1. First Detection up to 100 GeV via Stacking: Demonstrated that stacking analysis can detect high-energy afterglow emission up to 100 GeV for populations of GRBs that are individually undetectable.
2. Characterization of Faint GRBs: Provided the first precise light curves and spectra for a sample of 110 GRBs that were individually undetected by Fermi-LAT, revealing distinct temporal behaviors compared to bright GRBs.
3. Evidence for Energy Injection: Identified a discrepancy in the weak GRB sample that suggests the presence of energy injection effects in the GeV band, linking GeV afterglows to X-ray plateaus.
4. Validation of SSC Mechanism: Confirmed the transition from synchrotron to SSC dominance in the GeV band and the evolution of the cutoff energy over time.

5. Significance

• Jet Dynamics: The detection of jet breaks in the stacked light curves provides robust constraints on jet opening angles and initial Lorentz factors.
• Particle Acceleration: The results verify that GRB jets accelerate particles to energies capable of producing GeV photons via SSC processes, constraining the energy allocation between electrons and magnetic fields.
• Central Engine Physics: The identification of energy injection signatures in the faint sample offers new insights into the longevity and activity of the GRB central engine (e.g., magnetar spin-down or black hole accretion).
• Methodological Advance: Establishes stacking as a powerful tool for studying high-energy astrophysics when individual source statistics are insufficient, paving the way for future studies with next-generation telescopes (e.g., AMEGO-X, CTA).