Cosmic Earthquakes: Algorithm Reveals the Rhythm of Gamma-Ray Bursts ✛ ✦

Original: "An internal shock model calibrated with real gamma-ray burst light curves using a genetic algorithm"
arXiv:2607.12731 · 2026-07-14 · CC BY 4.0 · 3 min
A model of internal shock waves, calibrated by a genetic algorithm, showed: the central engine of gamma-ray bursts follows earthquake statistics.

In the depths of the Universe, when a massive star collapses into a black hole or a magnetar is born, a gamma-ray burst flares up — for an instant, it outshines all the stars in the observable Universe. But how does a dying star turn gravitational collapse into a dazzling dance of gamma photons? Physics offers an elegant model: the central engine ejects hundreds of ultra-relativistic shells — jets. They travel at nearly the speed of light (Lorentz factor Γ>100, as predicted by Einstein). When a fast shell catches up with a slow one, a shock wave is born — like the crack of a whip on a cosmic drum. In these collisions, kinetic energy transforms into radiation, which our telescopes then capture.

But until recently, the internal shock model remained more of a sketch than a precise tool. Astrophysicists decided to approach the problem like breeders: using a genetic algorithm, they 'grew' an optimal population of parameters capable of reproducing real light curves. The algorithm, inspired by Darwinian evolution, sifted through millions of combinations of shell numbers, masses, and delays, crossing and mutating the most successful variants. Training proceeded on three datasets at once: data from Swift/BAT, Fermi/GBM, and BATSE. And then something remarkable happened.

The genetic algorithm selected a distribution of the number of shells that precisely matches the Gutenberg–Richter law — the very one that describes the frequency of earthquakes of different magnitudes. Nature seems to use the same scenario for releasing energy: whether it's an accretion disk or a tectonic fault.

The optimized model perfectly reproduced the observed profiles: the average post-peak decay, the distribution of T20% durations, and even the Hubble law for fluxes — everything matched with an accuracy unattainable by previous semi-empirical fits. The simulated light curves — a superposition of pulses from thousands of collisions — replicate the characteristic 'hiccups' and plateaus of real bursts. It became clear: the central engine is a stochastic generator, ejecting shells with an exponentially decaying frequency. Each shell is a portion of the accretion disk torn from the total mass. All this resembles not an explosion, but a grand percussion, where the rhythm is dictated by internal turbulence.

Analysis showed: the typical decay time of the engine's activity τ is about three on a logarithmic scale. This means that after the burst starts, the probability of a new ejection drops as e^(-t/τ). The cosmic orchestra fades exponentially.

The connection with seismology reveals profound parallels. The Gutenberg–Richter law is a signature of self-organized criticality: the system balances on the edge of stability and releases stress in cascades. The same principle governs the triggering of snow avalanches, earthquakes, and, as it turns out, the birth of gamma-ray bursts. Perhaps accretion disks around supernovae and mergers of compact objects obey the same laws as tectonic plates. This transforms gamma-ray bursts into a unique laboratory: numerical simulations of relativistic magnetohydrodynamics gain a reliable statistical template for calibration.

Ahead lie even bolder plans. Scientists intend to incorporate the microphysics of shock waves into the model: synchrotron radiation and inverse Compton scattering, which depend on the viewing angle. This will allow predicting not only the shape but also the spectra of bursts. Accounting for redshift, the model will become a virtual telescope, predicting which bursts future missions — from SVOM to ATHENA — will see. We are learning to listen to the rhythm of dying stars in order to hear the music of the fabric of spacetime itself.

The shell statistics in gamma-ray bursts exactly reproduce the Gutenberg–Richter law, famous in seismology. Accretion disks around black holes apparently obey the same principles of self-organized criticality as the Earth's crust.

\epsilon = 1 - \frac{(m_r + m_s)\gamma_m}{m_r\gamma_r + m_s\gamma_s}
ϵ is the fraction of energy converted into radiation upon collision of two shells; m_r, m_s are the masses of the fast and slow shells; γ_r, γ_s, γ_m are their Lorentz factors and the factor of the merged shell.
f(t) = \frac{1}{\tau} e^{-t/\tau}
τ is the characteristic decay time of the central engine's activity; the probability of ejecting a new shell drops exponentially over time.
Scientists
Georges LemaîtreFritz ZwickyJacob BekensteinKarl SchwarzschildRalph AlpherStephen Hawking
Concepts
central engine gamma-ray burst gamma ray burst merger genetic algorithm stochastic gene expression afterglow swift satellite xgboost outburst generative model gene expression noise Gaussian distribution random number generator CRISPR gibbs sampling log-normal distribution