Calibrating the internal shock model of gamma-ray bursts using a genetic algorithm ✛ ✦
Context
Gamma-ray bursts (GRBs) are the most powerful explosions in the Universe, lasting from milliseconds to thousands of seconds. They are born from the collapse of massive stars into a supernova or the merger of compact objects, leaving behind a black hole or a magnetar. The key puzzle is exactly how the kinetic energy of ultra-relativistic jets is converted into gamma radiation. The internal shock (IS) model provides a mechanism for energy dissipation through collisions of inhomogeneous ejecta from the central engine, but its parameters have never been directly calibrated against observational data.
Methods
The researchers used a genetic algorithm to fit the IS model parameters to three catalogs: BATSE, Swift/BAT, and Fermi/GBM. The optimization was performed against six metrics: average post-peak profile, third moment, autocorrelation function, duration distributions, signal-to-noise ratio, and number of peaks. Each light curve was modeled as a superposition of pulses from shell collisions with realistic distributions of their number (generalized Zipf distribution) and ejection times (exponential), taking into account cosmological redshift and relativistic effects.
Results
The optimized model successfully reproduced key observational properties for all instruments. For example, for Swift/BAT, the average post-peak decay profile, T20% duration distribution, and flux distribution almost perfectly matched the real data. The model predicted that the number of shells per burst follows the Gutenberg–Richter law (similar to earthquake distribution), while ejection times follow an exponential decay with parameter τ ∼ 3 (on a logarithmic scale), pointing to the stochastic nature of the central engine activity. The model performed somewhat worse on the peak count distribution for BATSE, which may be due to instrumental specifics.
Implications
The results show that even a simplified IS model can explain a wide range of empirical properties of gamma-ray bursts. This strengthens the view of the central engine as a system evolving from a non-equilibrium state to equilibrium, where each fragment of the accretion disk has a constant probability of being ejected. The connection with seismic statistics suggests self-organized criticality in accretion processes.
Future development
In the future, it is planned to include the microphysics of relativistic shock waves (synchrotron and inverse Compton emission) and Doppler boosting effects depending on the viewing angle, in order to self-consistently calculate spectra and peak energy. The model can also be applied to catalogs from new missions such as SVOM and Insight-HXMT.
Impact
This work sets a new standard for testing theoretical models of gamma-ray bursts, bringing together big data and machine learning. The optimized IS model will become a tool for predicting GRB populations to be detected by future observatories.
Next steps
The next step will be to incorporate the energy spectrum and temporal evolution of the emission, as well as testing the model on observed samples with measured redshifts.
Key open problems
This research directly addresses a fundamental problem in gamma-ray burst physics — the nature of the central engine and the mechanism for converting energy into radiation. The obtained shell distributions link macroscopic observed variability with microscopic processes in accretion disks or neutron star magnetospheres, echoing unresolved questions in magnetic hydrodynamics and particle acceleration under extreme conditions.
The knowledge machine's view
This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.
What the work does
Researchers directly calibrated the parameters of an internal shock model against gamma-ray burst light curves from three catalogs using a genetic algorithm. The model reproduced key observed properties and also revealed a Gutenberg–Richter distribution of the number of shells and an exponential decay of ejection times.
Where its strength is
The strength of this work lies in directly confronting a simplified physical model with large observational samples, which allowed a quantitative assessment of the stochastic nature of the central engine's activity and linked it to the universal law of self-organized criticality.
Where it could go next
- It might be interesting to see how your internal shock model relates to observations of the record-long gamma-ray burst, where the central engine remained active for several hours—this is an extreme test for the distributions of shell numbers and ejection times.based on: Record Gamma-Ray Burst: 7 Hours of Invisible Light arXiv:2509.22792
- A related study found 34 pulsations in gamma-ray bursts, presumably linked to jet precession or accretion disk oscillations. Your approach to modeling the sequence of shells could help distinguish between these scenarios.based on: Rhythms of Black Holes: 34 Pulsations in Gamma-Ray Bursts arXiv:2510.11136
- There is a study in the archive that estimates jet formation time from the delay between gravitational-wave and gamma-ray signals. Your calibration could be extended by incorporating such multimessenger constraints on shell ejection times.based on: Light and Gravity Lag: The Key to Jet Birth arXiv:2606.29515
Worth trying
- Test whether the model with an exponential distribution of ejection times reproduces long-lasting activity episodes, as in the record-breaking 7-hour burst.
- Compare the predicted distribution of the number of shells with observed periodicity or quasi-periodicity in light curves to constrain the jet modulation mechanism.
- Add metrics based on multimessenger data (e.g., the delay between gravitational waves and gamma rays) to the loss function for joint calibration.
Why this may matter
The connection between gamma-ray burst statistics and self-organized criticality may point to a universal dissipation mechanism in accretion systems, and the proposed calibration method paves the way for systematic comparisons between central engine models and observational data.
Nearby works in our archive
- Record Gamma-Ray Burst: 7 Hours of Invisible Lightsimilarity 0.70
- Gamma-Ray Burst Wings: What Lies Beyond the Edge of the Beamsimilarity 0.70
- Rhythms of Black Holes: 34 Pulsations in Gamma-Ray Burstssimilarity 0.67
- X-ray Beacons in Galactic Centerssimilarity 0.66
- A Cosmic Blast Showed How Its Light Is Bornsimilarity 0.66
- Light and Gravity Lag: The Key to Jet Birthsimilarity 0.65
The archive is crowded around this work: 6 reviews nearby, closest at 0.70. This ground is well walked — and there is plenty to compare against.
Original abstract
as written by the authors arXiv:2607.12731
The origin of gamma-ray burst (GRB) prompt emission remains an open question. The internal shock (IS) model is a leading scenario for converting relativistic ejecta kinetic energy into gamma rays, but its parameters have not yet been fully calibrated against observed GRB light curves (LCs) to reproduce their diversity. We adopt a machine-learning framework to optimise the IS model by comparing simulated and observed LC properties from three GRB catalogues (Swift/BAT, Fermi/GBM, CGRO/BATSE). Assuming a redshift-dependent GRB formation rate, we employ a genetic algorithm to minimise a loss function based on six independent metrics capturing both average behaviours and statistical distributions. The optimised model reproduces several key observational properties, including the average post-peak temporal profile, autocorrelation function, and the distributions of duration, signal-to-noise ratio, number of peaks, peak flux, and fluence. We also derive constraints on the central engine activity: (i) the number of emitted shells is well described by a generalised Zipf distribution, analogous to the Gutenberg-Richter law for earthquakes, and (ii) the rest-frame shell-emission times follow a negative exponential distribution, indicating a stochastic process with a constant ejection probability. This calibrated IS model provides a physically grounded framework for interpreting GRB variability and predicting GRB populations detectable by future missions.
It turns out that the statistics of the number of shells in gamma-ray bursts exactly replicate the law of earthquake magnitude distribution — the famous Gutenberg–Richter law. Perhaps accretion disks behave like tectonic plates, accumulating and releasing stress.
Key numbers
- Number of GRBs in the Swift/BAT sample: 635
- Typical T20% duration: 10–100 s
- Shell distribution index α: ≈ 2.0
- Typical gamma-ray luminosity: 1052–1053 erg/s
- Total gamma-ray energy per event: ~1053 erg