The mystery of fast radio bursts (FRBs) has puzzled astrophysicists since the first event was discovered in 2007. These flashes last milliseconds but release energy comparable to the Sun's monthly output. Most FRBs come from distant galaxies, indicating a cosmological origin: their radio signal traveled billions of years to reach us, experiencing delays due to relativistic time dilation and dispersion in the intergalactic medium. This makes them useful probes for studying the expansion of the universe and the distribution of baryonic matter. Historically, neutron stars, predicted by Fritz Zwicky as remnants of supernovae, became prime candidates after Jocelyn Bell Burnell discovered pulsars and Subrahmanyan Chandrasekhar established their stability limit. Today, models split into two camps: one posits starquakes in the crust, the other magnetospheric flares. Neither provides a complete explanation.
The researchers assembled an extensive collection of event sequences: for five repeating FRBs (including the record-breaker FRB 20240114A with over 10,000 bursts in 8 months), magnetar X-ray flares, pulsar glitches, solar flares, and earthquakes. For each sequence, waiting times and energy changes were computed. Then, using the Pincus index (PI), a measure of entropy-based complexity, the degree of stochasticity was assessed, and through the Lyapunov exponent (LE) — chaoticity. Since radio signals travel at the speed of light, precise timing is critical. By combining PI and LE, they constructed a two-dimensional 'stochasticity–chaos' diagram (PLD), where each phenomenon is marked by a point. A key advantage of the method is that the diagnosis does not depend on absolute energetics but captures relative patterns.
On the PLD diagram, all studied classes distributed themselves in a patterned way. Earthquakes occupied a region of high chaos and low stochasticity — their sequences 'remember' past shocks. Solar and magnetar flares fell in the middle, consistent with their magnetic reconnection nature. Pulsar glitches showed high stochasticity and weak chaos. The most interesting part: repeating FRBs clustered into a compact cloud, markedly isolated from all others. A statistical test based on energy distance with permutations confirmed the significance of the difference between FRBs and magnetar flares (p=0.031) and glitches (p=0.020). At the same time, the internal consistency of the FRB class is high: splitting the sample into parts yields p=0.984, indicating homogeneity. The most active repeater FRB 20240114A deserves special attention. Its monitoring with the FAST telescope revealed colossal rate variations — from tens to almost 730 bursts per hour. Yet, its position on the PLD remained remarkably stable: over eight months of observations, PI hovered around 0.90, LE around 0.10, with no significant trends. The Dickey–Fuller test rejected non-stationarity with p-values around 10^{-3}. This means that the statistical 'handwriting' of the source does not depend on its instantaneous activity.
The isolation of FRBs in the 'stochasticity–chaos' phase space suggests that their physical mechanism is more complex than previously thought. High PI (close to unity) indicates an almost complete absence of deterministic memory: each burst is practically independent of the previous ones. Low LE points to weak chaos. Together, this paints a picture of a noisy, high-entropy source — perhaps a neutron star with a highly turbulent magnetosphere or multiple independent emission sites. This aligns with the lack of stable rotational periodicity, which has been searched for unsuccessfully in repeating FRBs. In the extreme gravitational fields of such objects, effects of time dilation according to Einstein become noticeable, which may influence the observed pattern. Traditional scenarios — singular magnetospheric flares or discrete starquake-glitches — do not reproduce the observed statistics.
The proposed PLD method opens a new path for classifying transients. In the future, it can be applied to other types of cosmic bursts, such as those from accreting black holes in quasars or events detected by gravitational wave detectors (LIGO and Virgo). As data accumulate from new repeaters and more sensitive instruments, it will be possible to refine the statistical patterns and possibly identify subclasses within the FRB population.
The work connects seemingly distant fields: neutron star astrophysics, plasma physics, seismology, and nonlinear dynamics. The results are important for understanding extreme states of matter in ultra-strong magnetic fields and may find resonance in laboratory experiments with high-power lasers.
The next step is to detect more bursts from a single source across different spectral ranges to test whether wavelength affects the position on the PLD. It is also important to include non-repeaters in the analysis to understand whether they are driven by different physics or by the same mechanism with infrequent repetition.
The nature of fast radio bursts remains one of the key unsolved problems of modern astrophysics, alongside the mysteries of dark matter and dark energy. This study narrows the range of possible models, but the fundamental question — exactly how the colossal energy is released in a millisecond radio pulse — remains open. Furthermore, the connection to pulsar glitches and sudden magnetic field reconfigurations sheds light on the internal structure of neutron stars, the equation of state of nuclear matter.
🎯 In a single second, the most active repeater FRB 20240114A belts out up to 12 bursts! And its total number of recorded events rivals the number of stars in some globular clusters. Meanwhile, each burst carries the energy the Sun radiates over several days.
🎬 In Carl Sagan's novel 'Contact,' radio signals from the constellation Vega become a bridge to interstellar dialogue. Fast radio bursts, discovered decades later, initially also sparked hypotheses of artificial origin, though today scientists lean toward natural, but no less exciting, explanations.