Gamma-ray Bursts Research Articles

The Dependence of Gamma-Ray Burst Jet Collimation on Black Hole Spin

Gamma-ray bursts (GRBs) are the most luminous events in the Universe and are excellent laboratories to study extreme physical phenomena in the cosmos. Despite a long trajectory of progress in understanding these highly energetic events, there are still many observed features that are yet to be fully explained. Observations of the jet opening angle of long gamma-ray bursts (LGRBs) suggest that LGRB jets are narrower for those GRBs at higher redshift. This phenomenon has been explained in the context of collimation by the stellar envelope, with denser (lower metallicity) stars at higher redshifts able to collimate the jet more effectively. However, until now, the dependence of the jet opening angle on the properties of the central engine has not been explored. We investigate the effect of black hole spin on the jet collimation angle for a magnetically launched jet, using the general relativistic radiation magnetohydrodynamical code ν bhlight. We present 3D results for a range of spin values. The simulations show that higher-spinning black holes tend to create narrower jets. If indeed LGRB progenitors in the early Universe are able to produce black hole central engines with higher spin, this could account for at least some of the observed jet opening angle-redshift correlation.

Self Consistent Modeling of Relativistic Runaway Electron Beams Giving Rise to Terrestrial Gamma‐Rays Flashes

AbstractTerrestrial Gamma Ray Flashes (TGFs) are short bursts of gamma rays occurring during thunderstorms. They are believed to be produced by Relativistic Runaway Electron Avalanches (RREAs). In this paper, we present a new numerical model based on the Particle‐In‐Cell (PIC) method to simulate the interactions between the electromagnetic fields and the electron avalanche self‐consistently. The code uses a cylindrical Yee lattice to numerically solve the electromagnetic fields, a Monte Carlo approach to simulate collisions with air molecules, and a plasma fluid model to calculate the effects of low‐energy electrons and ions. The model is first tested through the reproduction of dispersion relations in a hot plasma. RREAs propagating under a homogeneous background electric field are then simulated. Owing to the self‐consistent nature of this description, we report here new physical properties such as saturation processes in the electron density and in the number of high‐energy electrons, detailed dynamical screening of the electric field in the ion trail of the avalanche, and the associated electric currents. We find that the saturation of RREAs is obtained when the numbers of high‐energy electrons and photons is consistent with those believed to be representative of TGF sources.

On the Detection and Characterization of Quasiperiodic Oscillations in Astronomical Time Series: Gamma-Ray Burst X-Ray Light Curves as a Test Case

The study of temporal properties of variable sources can elucidate their physical processes. In this context, we present a critical study comparing three approaches to periodic or quasiperiodic behavior: Gaussian process, power spectrum, and wavelet analysis, using celerite, Lomb–Scargle periodograms, and weighted wavelet Z-transforms, respectively. We use 15 Swift X-ray Telescope light curves of short gamma-ray bursts (sGRBs) as examples. A comprehensive analysis of two sGRB X-ray light curves is performed. The results reveal the importance of artifacts, largely in the form of false quasiperiodic oscillation signals, possibly introduced by preprocessing (such as detrending) or other aspects of the analysis. The exploration described in this paper can be helpful for future studies of variability in gamma-ray bursts, active galactic nuclei, and other astronomical sources.

“Extended Emission” from Fallback Accretion onto Merger Remnants

Using a set of general-relativistic magnetohydrodynamics simulations that include proper neutrino transfer, we assess for the first time the role played by the fallback accretion onto the remnant from a binary neutron star merger over a timescale of hundreds of seconds. In particular, we find that, independently of the equation of state, the properties of the binary, and the fate of the remnant, the fallback material reaches a total mass of ≳10−3 M ⊙, i.e., about 50% of the unbound matter, and that the fallback accretion rate follows a power law in time with slope ∼t −5/3. Interestingly, the timescale of the fallback and the corresponding accretion luminosity are in good agreement with the so-called “extended emission” observed in short gamma-ray bursts (GRBs). Using a simple electromagnetic emission model based on the self-consistent thermodynamical state of the fallback material heated by r-process nucleosynthesis, we show that this fallback material can shine in gamma and X-rays with luminosities ≳1048 erg s−1 for hundreds of seconds, thus making it a good and natural candidate to explain the extended emission in short GRBs. Additionally, our model for fallback emission reproduces well and rather naturally some of the phenomenological traits of extended emission, such as its softer spectra with respect to the prompt emission and the presence of exponential cutoffs in time. Our results clearly highlight that fallback flows onto merger remnants cannot be neglected, and the corresponding emission represents a very promising and largely unexplored avenue to explain the complex phenomenology of GRBs.

Evidence of a Strong 19.5 Hz Flux Oscillation in Swift BAT and Fermi GBM Gamma-Ray Data from GRB 211211A

The gamma-ray burst (GRB) GRB 211211A is believed to have occurred due to the merger of two neutron stars or a neutron star and a black hole, despite its duration of more than a minute. Subsequent analysis has revealed numerous interesting properties including the possible presence of a ∼22 Hz quasiperiodic oscillation (QPO) during precursor emission. Here we perform timing analysis of Fermi and Swift gamma-ray data on GRB 211211A and, although we do not find a strong QPO during the precursor, we do find an extremely significant 19.5 Hz flux oscillation, which has higher fractional amplitude at higher energies, in a ∼0.2 s segment beginning ∼1.6 s after the start of the burst. After presenting our analysis we discuss possible mechanisms for the oscillation.

Distribution of the number of peaks within a long gamma-ray burst

Context. The variety and complexity of long duration gamma-ray burst (LGRB) light curves (LCs) encode a wealth of information about the way LGRB engines release their energy following the collapse of the progenitor massive star. Thus far, attempts to characterise GRB LCs have focused on a number of properties, such as the minimum variability timescale and power density spectra (both ensemble average and individual), or considering different definitions of variability. In parallel, a characterisation as a stochastic process has been pursued by studying the distributions of waiting times, peak flux, and fluence of individual peaks that can be identified within GRB time profiles. However, an important question remains as to whether the diversity of GRB profiles can be described in terms of a common stochastic process. Aims. Here, we address this issue by extracting and modelling, for the first time, the distribution of the number of peaks within a GRB profile. Methods. We analysed four different GRB catalogues: CGRO/BATSE, Swift/BAT, BeppoSAX/GRBM, and Insight-HXMT. The statistically significant peaks were identified by means of well tested and calibrated algorithm MEPSA and further selected by applying a set of thresholds on the signal-to-noise ratio. We then extracted the corresponding distributions of number of peaks per GRB. Results. Among the different models considered (power-law, simple or stretched exponential), we find that only a mixture of two exponentials was able to model all the observed distributions. This suggests the existence of two distinct behaviours: (i) an average number of 2.1 ± 0.1 peaks per GRB (“peak-poor”), accounting for about 80% of the observed population of GRBs; and (ii) an average number of 8.3 ± 1.0 peaks per GRB (“peak-rich”), accounting for the remaining 20% of the observed population. Conclusions. We associate the class of peak-rich GRBs with the presence of sub-second variability, which appears to be surprisingly absent among peak-poor GRBs. The two classes could result from two distinct regimes in which the inner engines of GRBs release their energy or otherwise dissipate that energy as gamma rays.

A proposal for a fast infrared bursts detector

The gravitational wave GW170817 from a binary neutron star merger and the simultaneous electromagnetic detection of the GRB170717A by Fermi Gamma-Ray Space Telescope, opened a new era in the multi-messenger astronomy. Furthermore, the GRBs (Gamma-Ray Bursts) and the mysterious FRBs (Fast Radio Bursts) have sparked interest in the development of new detectors and telescopes dedicated to the time-domain astronomy across the electromagnetic spectrum. Time-domain astronomy aims to acquire fast astronomical bursts in temporal range between a few seconds down to 1 ns. Fast InfraRed Bursts (FIRB's) have been relatively understudied, often due to the lack of appropriate tools for observation and analysis. In this scientific scenario, the present contribution proposes a new detection system for ground-based reflecting telescopes working in the mid-infrared (mid-IR) range to search for astronomical FIRB's. Experience developed in the diagnostics for lepton circular accelerators can be used to design temporal devices for astronomy. Longitudinal diagnostic instruments acquire bunch-by-bunch particle shifts in the direction of flight, that is equivalent to temporal. Transverse device integrates the beam signal in the horizontal and vertical coordinates, as standard telescopes. The proposed instrument aims to work in temporal mode. Feasibility study tests have been carried out at SINBAD, the infrared beam line of DAFNE, the e+e- collider of INFN. SINBAD releases pulsed infrared synchrotron light with 2.7 ns separation. The front-end detector system has been evaluated to detect temporal fast infrared signals with 2–12 μm wavelengths and 1 ns rise times. The present contribute aims to be a step toward a feasibility study report.

Constraints on the z ∼ 5 Star-forming Galaxy Luminosity Function From Hubble Space Telescope Imaging of an Unbiased and Complete Sample of Long Gamma-Ray Burst Host Galaxies

We present rest-frame UV Hubble Space Telescope imaging of the largest and most complete sample of 23 long-duration gamma-ray burst (GRB) host galaxies between redshifts 4 and 6. Of these 23, we present new WFC3/F110W imaging for 19 of the hosts, which we combine with archival WFC3/F110W and WFC3/F140W imaging for the remaining four. We use the photometry of the host galaxies from this sample to characterize both the rest-frame UV luminosity function (LF) and the size–luminosity relation of the sample. We find that when assuming the standard Schechter-function parameterization for the UV LF, the GRB host sample is best fit with α=−1.30−0.25+0.30 and M*=−20.33−0.54+0.44 mag, which are consistent with results based on z ∼ 5 Lyman-break galaxies. We find that ∼68% of our size–luminosity measurements fall within or below the same relation for Lyman-break galaxies at z ∼ 4. This study observationally confirms expectations that at z ∼ 5 Lyman-break and GRB host galaxies should trace the same population and demonstrates the utility of GRBs as probes of hidden star formation in the high-redshift Universe. Under the assumption that GRBs unbiasedly trace star formation at this redshift, our nondetection fraction of 7/23 is consistent at the 95% confidence level with 13%–53% of star formation at redshift z ∼ 5 occurring in galaxies fainter than our detection limit of M 1600Å ≈ −18.3 mag.

A high-redshift calibration of the [O I]-to-H I conversion factor in star-forming galaxies

The assembly and build-up of neutral atomic hydrogen (H I) in galaxies is one of the most fundamental processes in galaxy formation and evolution. Studying this process directly in the early universe is hindered by the weakness of the hyperfine 21-cm H I line transition, impeding direct detections and measurements of the H I gas masses (MHI). Here we present a new method to infer MHI of high-redshift galaxies using neutral, atomic oxygen as a proxy. Specifically, we derive metallicity-dependent conversion factors relating the far-infrared [O I]-63 μm and [O I]-145 μm emission line luminosities and MHI in star-forming galaxies at z ≈ 2 − 6 using gamma-ray bursts (GRBs) as probes. We calibrate the [O I]-to-H I conversion factor relying on a sample of local galaxies with direct measurements of MHI and [O I]-63 μm and [O I]-145 μm line luminosities in addition to the SIGAME hydrodynamical simulation framework at similar epochs (z ≈ 0). We find that the [O I]63 μm-to-H I and [O I]145 μm-to-H I conversion factors, here denoted β[OI]−63 μm and β[OI]−145 μm, respectively, universally appear to be anti-correlated with the gas-phase metallicity. The GRB measurements further predict a mean ratio of L[OI]−63 μm/L[OI]−145 μm = 1.55 ± 0.12 and reveal generally less excited [C II] over [O I] compared to the local galaxy sample. The z ≈ 0 galaxy sample also shows systematically higher β[OI]−63 μm and β[OI]−145 μm conversion factors than the GRB sample, indicating either suppressed [O I] emission in local galaxies likely due to their lower hydrogen densities or more extended, diffuse H I gas reservoirs traced by the H I 21-cm. Finally, we apply these empirical calibrations to the few detections of [O I]-63 μm and [O I]-145 μm line transitions at z ≈ 2 from the literature and further discuss the applicability of these conversion factors to probe the H I gas content in the dense, star-forming interstellar medium (ISM) of galaxies well into the epoch of reionization.

GRB 221009A: Spectral Signatures Based on ALPs Candidates

GRB 221009A has posed a significant challenge to our current understanding of the mechanisms that produce TeV photons in gamma-ray bursts (GRBs). On one hand, the Klein–Nishina (KN) effect of the inverse Compton scattering leads to less efficient energy losses of high-energy electrons. On the other hand, at a redshift of 0.151, the TeV spectrum of GRB 221009A undergoes significant absorption by the extragalactic background light (EBL). Therefore, the observation of a 13 TeV photon in this event implies the presence of enormous photon fluxes at the source, which the synchrotron self-Compton mechanism in external shocks cannot easily generate. As an alternative, some authors have suggested the possibility of converting the TeV photons into axion-like particles (ALPs) at the host galaxy, in order to avoid the effects of EBL absorption, and then reconverting them into photons within the Milky Way. While this solution relaxes the requirement of very high photon fluxes, the KN effect still poses a challenge. Previously, we have shown that the injections of ALPs could explain the observation of 13 TeV photons. Here, we include the energy dependence of the probability of survival and the amount of energy carried to determine the ALP candidates, which could potentially explain the TeV photons observed by the Large High Altitude Air Shower Observatory and their hard spectrum. We found that the allowed candidates are generally clustered around masses of 10−7 eV. We also considered different EBL models, for the one predicting larger attenuation tends to reject ALP candidates with the lowest coupling factor. For some hypotheses of the EBL model, these candidates are found below a region of the parameter space in which, if detected, ALPs could account for all of the cold dark matter in the Universe.

Model-independent test of the running Hubble constant from the Type Ia supernovae and the Hubble parameter data

ABSTRACT In this study, we model-independently investigate the behaviour of running Hubble constant, characterized by the fit function $H_{\rm 0}(z)=\tilde{H_{\rm {0}}}/(1+z)^{\alpha }$, where α represents the evolutionary parameter and ${\tilde{H_{\rm {0}}}}$ corresponds to the current value of Hubble constant. Our analysis utilizes the expansion rate E(z) data points measured from the Pantheon + Multi-Cycle Treasury compilation of Type Ia supernova data, the measurements of H0 obtained by Riess et al., and the Hubble parameter H(z) data obtained from the differential ages of passive galaxies [known as cosmic chronometer (CC) method] and from the baryon acoustic oscillation (BAO) in the radial direction of galaxy clustering. To resolve the redshift mismatch problem between the E(z) and H(z) data sets, we adopt the Hubble parameter data obtained via CC or BAO along with the measurements of H0 obtained by Riess et al. to reconstruct the H(z) function using the Gaussian process. Our constraint yields α values of 0.125 ± 0.063 or 0.095 ± 0.052 when combining six pairs of the E(z) data and the reconstructed H(z) points via CC or BAO. These findings reveal that the Hubble constant may evolve with redshift, exhibiting a slowly decreasing trend, with α coefficients consistent with zero only at 2.0σ or 1.8σ. Therefore, the running Hubble constant might offer a promising resolution to the Hubble tension, and its reliability should be further tested through high-precision measurement at higher redshifts, such as the upcoming gamma-ray bursts and quasars.

Gamma-ray burst studies with wide-field X-ray detectors on a SmallSat platform

Gamma-ray burst studies with wide-field X-ray detectors on a SmallSat platform

Scaling relations for gamma-ray burst afterglow light curves and centroid motion independent of jet structure and dynamics

ABSTRACT Models for gamma-ray burst afterglow dynamics and synchrotron spectra are known to exhibit various scale invariances, owing to the scale-free nature of fluid dynamics and the power-law shape of synchrotron spectra. Since GRB 170817A, off-axis jet models including a lateral energy structure in the initial outflow geometry have gained in prominence. Here, we demonstrate how the scale invariance for arbitrary jet structure and dynamical stage can be expressed locally as a function of jet temporal light-curve slope. We provide afterglow flux expressions and demonstrate their use to quickly assess the physical implications of observations. We apply the scaling expressions to the Swift X-ray Telescope sample, which shows a spread in observed fluxes, binned by light-curve slope at time of observation, that increases with increasing light-curve slope. According to the scaling relations, this pattern is inconsistent with a large spread in environment densities if these were the dominant factor determining the variability of light curves. We further show how the late deep Newtonian afterglow stage remains scale-invariant but adds distinct spectral scaling regimes. Finally, we show that for given jet structure a universal curve can be constructed of the centroid offset, image size, and ellipticity (that can be measured using very large baseline interferometry) versus observer angle, in a manner independent of explosion energy and circumburst density. Our results apply to any synchrotron transient characterized by a release of energy in an external medium, including supernova remnants, kilonova afterglows, and soft gamma-repeater flares.

A magnetar giant flare in the nearby starburst galaxy M82.

Magnetar giant flares are rare explosive events releasing up to 1047 erg in gamma rays in less than 1 second from young neutron stars with magnetic fields up to 1015-16 G (refs. 1,2). Only three such flares have been seen from magnetars in our Galaxy3,4 and in the Large Magellanic Cloud5 in roughly 50 years. This small sample can be enlarged by the discovery of extragalactic events, as for a fraction of a second giant flares reach luminosities above 1046 erg s-1, which makes them visible up to a few tens of megaparsecs. However, at these distances they are difficult to distinguish from short gamma-ray bursts (GRBs); much more distant and energetic (1050-53 erg) events, originating in compact binary mergers6. A few short GRBs have been proposed7-11, with different amounts of confidence, as candidate giant magnetar flares in nearby galaxies. Here we report observations of GRB 231115A, positionally coincident with the starburst galaxy M82 (ref. 12). Its spectral properties, along with the length of the burst, the limits on its X-ray and optical counterparts obtained within a few hours, and the lack of a gravitational wave signal, unambiguously qualify this burst as a giant flare from a magnetar in M82.

Naked Forward Shock Seen in the TeV Afterglow Data of GRB 221009A

We explore the implications of the light curve of the early TeV gamma-ray afterglow of GRB 221009A reported by the LHAASO Collaboration. We show that the reported offset of the reference time, T *, allows the determination of the activation time of the main jet component, which occurs approximately 200 s after the Gamma-ray Burst Monitor (GBM) trigger time and closely precedes the moment at which GBM was saturated. We find that while the LHAASO data do not exclude the homogeneous circumburst medium scenario, the progenitor wind scenario looks preferable, finding excellent agreement with the expected size of the stellar bubble. We conclude that the initial growth of the light curve is dominated by processes internal to the jet or by gamma-gamma attenuation on the photons emitted during the prompt phase. Namely, either the activation of the acceleration process or the decrease of internal gamma-gamma absorption can naturally explain the initial rapid flux increase. The subsequent slow flux growth phase observed up to T * + 18 s can be related to the process of particle accumulation in the production region. The duration of this phase allows an almost parameter-independent determination of the jet’s initial Lorentz factor, Γ0 ≈ 600, and magnetic field strength, B′∼0.3G . These values appear to match well those previously revealed through spectral modeling of the GRB emission.

Ultrafast Compact Binary Mergers

The duration of orbital decay induced by gravitational waves (GWs) is often the bottleneck of the evolutionary phases going from star formation to a merger. We show here that kicks imparted to the newly born compact object during the second collapse generically result in a GW merger time distribution behaving like dN/dlogt∝t2/7 at short durations, leading to ultrafast mergers. Namely, a nonnegligible fraction of neutron star binaries, formed in this way, will merge on a timescale as short as 10 Myr, and a small fraction will merge even on a timescale less than 10 kyr. The results can be applied to different types of compact binaries. We discuss here the implications for binary neutron star mergers. These include unique short gamma-ray bursts (GRBs), eccentric and misaligned mergers, r-process enrichment in the very early Universe and in highly star-forming regions, and possible radio precursors. Interestingly, we conclude that among the few hundred short GRBs detected so far, a few must have formed via this ultrafast channel.

Application of Deep-learning Methods for Distinguishing Gamma-Ray Bursts from Fermi/GBM Time-tagged Event Data

To investigate gamma-ray bursts (GRBs) in depth, it is crucial to develop an effective method for identifying GRBs accurately. Current criteria, e.g., onboard blind search, ground blind search, and target search, are limited by manually set thresholds and perhaps miss GRBs, especially for subthreshold events. We proposed a novel approach that utilizes convolutional neural networks (CNNs) to distinguish GRBs and non-GRBs directly. We structured three CNN models, plain-CNN, ResNet, and ResNet-CBAM, and endeavored to exercise fusing strategy models. Count maps of NaI detectors on board Fermi/Gamma-ray Burst Monitor were employed, as the input samples of data sets and models were implemented to evaluate their performance on different timescale data. The ResNet-CBAM model trained on the 64 ms data set achieves high accuracy overall, which includes residual and attention mechanism modules. The visualization methods of Grad-CAM and t-SNE explicitly displayed that the optimal model focuses on the key features of GRBs precisely. The model was applied to analyze 1 yr data, accurately identifying approximately 98% of GRBs listed in the Fermi burst catalog, eight out of nine subthreshold GRBs, and five GRBs triggered by other satellites, which demonstrated that the deep-learning methods could effectively distinguish GRBs from observational data. Besides, thousands of unknown candidates were retrieved and compared with the bursts of SGR J1935+2154, for instance, which exemplified the potential scientific value of these candidates indeed. Detailed studies on integrating our model into real-time analysis pipelines thus may improve their accuracy of inspection and provide valuable guidance for rapid follow-up observations of multiband telescopes.

Gamma-Ray Burst Pulses and Lateral Jet Motion

We propose that gamma-ray burst (GRB) pulses are produced when highly relativistic jets sweep across an observer’s line of sight. We hypothesize that axisymmetric jet profiles, coupled with special relativistic effects, produce the time-reversed properties of GRB pulses. Curvature resulting from rapid jet expansion is responsible for much of the observed pulse asymmetry and hard-to-soft evolution. The relative obliqueness with which the jet crosses the line of sight explains the known GRB pulse morphological types. We explore two scenarios: one in which a rigid/semirigid jet moves laterally and another in which a ballistic jet sprays material from a laterally moving nozzle. The ballistic jet model is favored based upon its consistency with standard emission mechanisms.

Future Perspectives for Gamma-ray Burst Detection from Space

Since their first discovery in the late 1960s, gamma-ray bursts have attracted an exponentially growing interest from the international community due to their central role in the most highly debated open questions of the modern research of astronomy, astrophysics, cosmology, and fundamental physics. These range from the intimate nuclear composition of high-density material within the core of ultra-dense neuron stars, to stellar evolution via the collapse of massive stars, the production and propagation of gravitational waves, as well as the exploration of the early universe by unveiling the first stars and galaxies (assessing also their evolution and cosmic re-ionization). GRBs in the past ∼50 years have stimulated the development of cutting-edge technological instruments for observations of high-energy celestial sources from space, leading to the launch and successful operations of many different scientific missions (several of them still in data-taking mode currently). In this review, we provide a brief description of the GRB-dedicated missions from space being designed and developed for the future. The list of these projects, not meant to be exhaustive, shall serve as a reference to interested readers to understand what is likely to come next to lead the further development of GRB research and the associated phenomenology.

A comparative study of outflow structures of two classes of gamma-ray bursts

ABSTRACT The outflow structures of gamma-ray bursts (GRBs) can provide insights into the origins and radiation mechanisms of these cosmic explosions. We systematically study the GRB outflow structures by modelling their afterglow light curves and check if the structures of long gamma-ray bursts (LGRBs) and short gamma-ray bursts (SGRBs) are different. The sample consists of Swift-XRT afterglows with sufficient coverage and known redshift, which includes 195 well-fit LGRBs and 13 well-fit SGRBs. The model we use is a two-parameter ‘boosted fireball’ model, which consists of a family of outflows, with shapes varying smoothly from a quasi-spherical outflow to a highly collimated jet. We use the jetfit package to fit afterglow light curves and obtain the jet parameters. We find that there are no statistical differences in the distributions of jet parameters between LGRBs and SGRBs by performing K–S test and 74 per cent of the ratios of the observer angle to jet opening angle are in the range of 0.2 to 1. Our analysis indicates that the majority of GRB afterglows are viewed off-axis and there has no statistical difference between LGRBs and SGRBs. We also find that both the LGRBs and SGRBs exhibit two similar correlations: the jet opening angle is positively correlated with the observer angle, with the correlation coefficient 0.61 for LGRBs and 0.63 for SGRBs; the circumburst density is inversely correlated with the explosion energy with the correlation coefficient −0.89 for LGRBs and −0.69 for SGRBs. Our results suggest that the outflow structures are similar for the LGRBs and SGRBs.