Long-duration Gamma-ray Bursts Research Articles

Unveiling the Progenitors of a Population of Likely Peculiar Gamma-Ray Bursts

Traditionally, gamma-ray bursts (GRBs) are classified as long and short GRBs, with T 90 = 2 s being the threshold duration. Generally, long-duration GRBs (LGRBs; T 90 > 2 s) are associated with the collapse of massive stars, and short-duration (SGRBs; T 90 < 2 s) are associated with compact binary mergers involving at least one neutron star. However, the existence of a population of so-called “peculiar GRBs”—i.e., LGRBs originating from mergers or long Type I GRBs, and SGRBs originating from collapsars or short Type II GRBs—has challenged the traditional paradigm of GRB classification. Finding more peculiar GRBs may help to give more insight into this issue. In this work, we analyze the properties of machine-learning-identified long Type I GRB and short Type II GRB candidates, long GRBs-I and short GRBs-II (the so-called “peculiar GRBs”). We find that long GRBs-I almost always exhibit properties similar to Type I GRBs, which suggests that mergers may indeed produce GRBs with T 90 > 2 s. Furthermore, according to the probability given by the redshift distribution, short GRBs-II almost exhibit properties similar to Type II GRBs. This suggests that the populations of short Type II GRBs are not scarce and that they are hidden in a large number of samples without redshifts, which is unfavorable for the interpretation that the jet progression leads to a missed main emission.

Repeated Partial Disruptions in a White Dwarf–Neutron Star or White Dwarf–Black Hole Merger Modulate the Prompt Emission of Long-duration Merger-type GRBs

The progenitors of gamma-ray bursts (GRBs) have long been an unresolved issue. GRB 230307A stands out as an exceptionally bright event, belonging to the long-duration GRBs but also exhibiting a late-emission component reminiscent of a kilonova. Together with the similar events GRBs 060614 and 211211A, they make up a new subgroup of GRBs with intriguing progenitors. If such long-duration merger-type GRBs originated from the coalescence of a white dwarf (WD) with a neutron star (NS) or a black hole (BH), as proposed in the recent literature, then the larger tidal disruption radius of the WD, together with a nonnegligible residual orbital eccentricity, would make repeated partial tidal disruptions inevitable. This may modulate the mass accretion and jet launching process at the NS or BH, resulting in a quasiperiodic modulation (QPM) in the light curve of the GRB, with a period equal to the orbital period. The detection of potential QPMs during the early episode of prompt emission of these three GRBs supports this scenario, and the relatively slow QPM (> 1 s) suggests that the lighter object cannot be an NS. We propose that the progenitor system of GRBs 230307A, 060614, and 211211A consist of a WD of mass 1.3 M ⊙, 0.9 M ⊙, and 1.4 M ⊙, respectively, and an NS (or BH). After several cycles of modulations, the WD is completely destroyed, and the accretion of the remaining debris dominates the extended emission episode.

Closing the Net on Transient Sources of Ultrahigh-energy Cosmic Rays

The arrival directions of ultrahigh-energy cosmic rays (UHECRs) observed above 4 × 1019 eV provide evidence of localized excesses that are key to identifying their sources. We leverage the 3D matter distribution from optical and infrared surveys as a density model of UHECR sources, which are considered to be transient. Agreement of the sky model with UHECR data imposes constraints on both the emission rate per unit matter and the time spread induced by encountered turbulent magnetic fields. Based on radio measurements of cosmic magnetism, we identify the Local Sheet as the magnetized structure responsible for the kiloyear duration of UHECR bursts for an observer on Earth and find that the turbulence amplitude must be within 0.5–20 nG for a coherence length of 10 kpc. At the same time, the burst-rate density must be above 50 Gpc−3 yr−1 for Local Sheet galaxies to reproduce the UHECR excesses and below 5000 Gpc−3 yr−1 (30,000 Gpc−3 yr−1) for the Milky Way (Local Group galaxies) not to outshine other galaxies. For the transient emissions of protons and nuclei to match the energy spectra of UHECRs, the kinetic energy of the outflows responsible for UHECR acceleration must be below 4 × 1054 erg and above 5 × 1050 erg (2 × 1049 erg) if we consider the Milky Way (or not). The only stellar-sized transients that satisfy both Hillas’ and our criteria are long-duration gamma-ray bursts.

Expected Gamma-Ray Burst Detection Rates and Redshift Distributions for the BlackCAT CubeSat Mission

We report the results of an extensive set of simulations exploring the sensitivity of the BlackCAT CubeSat to long-duration gamma-ray bursts (GRBs). BlackCAT is a NASA APRA-funded CubeSat mission for the detection and real-time subarcminute localization of high-redshift (z ≳ 3.5) GRBs. Thanks to their luminous and long-lived afterglow emissions, GRBs are uniquely valuable probes of high-redshift star-forming galaxies and the intergalactic medium. In addition, each detected GRB with a known redshift serves to localize a region of high-redshift star formation in three dimensions, enabling deep follow-on searches for host galaxies and associated local and large-scale structures. We explore two distinct models for the GRB redshift distribution and luminosity function, both consistent with Swift observations. We find that, for either model, BlackCAT is expected to detect a mean of 42 bursts per year on orbit, with 6.7% to 10% of these at z > 3.5. BlackCAT bursts will be localized to an r 90 ≲ 55″ precision and reported to the community within seconds. Due to the mission orbit and pointing scheme, bursts will be located in the night sky and well placed for deep multiwavelength follow-up observations. BlackCAT is on schedule to achieve launch readiness in 2025.

The Peculiar Precursor of a Gamma-Ray Burst from a Binary Merger Involving a Magnetar

The milestone discovery of GW170817-GRB 170817A-AT 2017gfo has shown that gravitational waves (GWs) could be produced during the merger of a neutron star–neutron star/black hole and that in electromagnetic (EM) waves, a gamma-ray burst (GRB) and a kilonova (KN) are generated in sequence after the merger. Observationally, however, EM properties before the merger phase are still unclear. Here we report a peculiar precursor in a KN-associated long-duration GRB 211211A, providing evidence of the EM before the merger. This precursor lasts ∼0.2 s, and the waiting time between the precursor and the main burst is ∼1 s, comparable to that between GW170817 and GRB 170817A. The spectrum of the precursor could be well fit with a nonthermal cutoff power-law model instead of a blackbody. In particular, a ∼22 Hz quasiperiodic oscillation candidate (∼3σ) is detected in the precursor. These temporal and spectral properties indicate that this precursor is probably produced by a catastrophic flare accompanied with magnetoelastic or crustal oscillations of a magnetar in a binary compact merger. The strong magnetic field of the magnetar can also account for the prolonged duration of GRB 211211A. However, it poses a challenge to reconcile the rather short lifetime of a magnetar with the rather long spiraling time of a binary neutron star system only by the GW radiation before the merger.

A Comparative Analysis of Two Peculiar Gamma-Ray Bursts: GRB 230307A and GRB 211211A

GRB 211211A is a peculiar long gamma-ray burst (GRB) with very high brightness and short burst properties. Its full light curve consists of three emission episodes, i.e., a precursor, a main burst, and an extended emission. We find that a recently detected long-duration GRB 230307A also includes the three consistent emission episodes. Furthermore, the two bursts have similar redshifts, 0.076 and 0.065, respectively. We perform a detailed temporal and spectral analysis of the two GRBs to compare their temporal and spectral properties. Our analysis shows that the two bursts share great similarities for both the whole emission and the three corresponding emission phases, which are listed as follows: (1) they have near zero spectral lag; (2) they have very short minimum variability timescale (MVT); (3) they lie in the same region of in the MVT–T 90, Amati relation and hardness–T 90 planes; (4) the three phases are quasi-thermal spectra; (5) both the peak energy and the low-energy index track the flux; (6) the time-resolved spectra are much wider than those of the blackbody predicted by the model; (7) there are strong correlations between thermal flux and total flux and the correlation coefficients, and the slopes for the corresponding stages are very consistent; and (8) the photosphere emission properties are very consistent. Other investigations and observations suggest that the two GRBs indeed belong to a short burst with a compact star merger origin. Therefore, we think that GRB 230307A and GRB 211211A are rare and similar GRBs, and the photospheric radiation can interpret their radiation mechanisms.

TeV afterglow from GRB 221009A: photohadronic origin?

ABSTRACT Gamma-ray burst (GRB), GRB 221009A, a long-duration GRB, was observed simultaneously by the Water Cherenkov Detector Array (WCDA) and the Kilometer Squared Array (KM2A) of the Large High Altitude Air Shower Observatory (LHAASO) during the prompt emission and the afterglow periods. Characteristic multi-TeV photons up to 13 TeV were observed in the afterglow phase. The observed very high-energy (VHE) gamma-ray spectra by WCDA and KM2A during different time intervals and in different energy ranges can be explained very well in the context of the photohadronic model with the inclusion of extragalactic background light models. In the photohadronic scenario, interaction of high-energy protons with the synchrotron self-Compton (SSC) photons in the forward shock region of the jet is assumed to be the source of these VHE photons. The observed VHE spectra from the afterglow of GRB 221009A are similar to the VHE gamma-ray spectra observed from the temporary extreme high-energy peaked BL Lac (EHBL), 1ES 2344+514 only during the 11th and the 12th of August, 2016. Such spectra are new and have been observed for the first time in a GRB.

GRB Afterglows with Energy Injections in AGN Accretion Disks

Active galactic nucleus (AGN) disks are widely considered potential hosts for various high-energy transients, including gamma-ray bursts (GRBs). The reactivation of GRB central engines can provide additional energy to shocks formed during the interaction of the initially ejected GRB jets with the circumburst material, commonly referred to as energy injections. In this paper, we study GRBs occurring in AGN disks within the context of energy injections. We adopt the standard external forward shock (EFS) model and consider both short- and long-duration GRB scenarios. Light curves for two types of radiation, namely, the radiation from the heated disk material (RHDM) and GRB afterglows, are computed. We find that the energy injection facilitates the EFS to break out from the photosphere of the low-density AGN disk at relativistic velocity. Moreover, the energy injection almost does not affect the RHDM but significantly enhances the peak flux of the GRB afterglows.

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.

GRB 220304A: Another Gamma-Ray Burst Dominated by Thermal Radiation

We report a recently detected long-duration gamma-ray burst (GRB) event by Fermi-GBM, GRB 220304A. The spectral analysis of the burst by the Band function shows that both the time-integrated and time-resolved spectra are very narrow, with low-energy spectral index 〈α〉 = −0.05 ± 0.30 and high-energy spectral index 〈β〉 = −3.53 ± 0.30. It is reminiscent of GRB 090902B, a special GRB with photosphere radiation characteristics. Then, we perform spectral analysis using the Planck function (blackbody, BB) and the multicolor BB model. It is found that the spectra within −1 ∼ 3 s is well fit by the BB model, indicating that the observation within the first 4 s is a pure thermal event. Apart from that most of the spectra can be well modeled as a multicolor BB. We conclude that GRB 220304A is a GRB dominated by thermal radiation. We also find that the spectral widths of the time-resolved spectrum of the burst gradually increase with time. Based on the Amati relation, we infer the redshift to be 0.23, the physical properties of the relativistic outflow, and find that the relationship between the bulk Lorentz factor and the isotropic luminosity Γ–L iso,52 still exists.

Ejecta–Circumstellar Medium Interaction in High-density Environment Contribution to Kilonova Emission: Application to GRB 191019A

The nearby long-duration GRB 191019A recently detected by Swift lacks an associated supernova and belongs to a host galaxy with little star formation activity, suggesting that the origin of this burst is the result of a merger of two compact objects with dynamical interactions in a high-density medium of an active galactic nucleus. Given the potential motivation of this event, and given that it occurs in such a high-density environment, the ejecta–circumstellar medium (CSM) interaction cannot be ignored as possibly contributing to the kilonova emission. Here, we theoretically calculate the kilonova emission by considering the contribution of the ejecta–CSM interaction in a high-density environment. We find that the contribution to the kilonova emission from the ejecta–CSM interaction will dominate at a later time, and a smaller ejecta mass will have a stronger kilonova emission from the ejecta–CSM interaction. Moreover, we try to apply it to GRB 191019A, but we find that it is difficult to identify the possible kilonova emission from the observations, due to the contribution of the bright host galaxy. On the other hand, less injected mass (less than M ej = 2 × 10−5 M ⊙) will be required if one can detect the kilonova emission associated with a GRB 191019A–like event in the future. The r-process-powered and spin energy contributions from the magnetar are also discussed.

Inferring the Redshift of More than 150 GRBs with a Machine-learning Ensemble Model

Gamma-ray bursts (GRBs), due to their high luminosities, are detected up to a redshift of 10, and thus have the potential to be vital cosmological probes of early processes in the Universe. Fulfilling this potential requires a large sample of GRBs with known redshifts, but due to observational limitations, only 11% have known redshifts (z). There have been numerous attempts to estimate redshifts via correlation studies, most of which have led to inaccurate predictions. To overcome this, we estimated GRB redshift via an ensemble-supervised machine-learning (ML) model that uses X-ray afterglows of long-duration GRBs observed by the Neil Gehrels Swift Observatory. The estimated redshifts are strongly correlated (a Pearson coefficient of 0.93) and have an rms error, namely, the square root of the average squared error 〈Δz 2〉, of 0.46 with the observed redshifts showing the reliability of this method. The addition of GRB afterglow parameters improves the predictions considerably by 63% compared to previous results in peer-reviewed literature. Finally, we use our ML model to infer the redshifts of 154 GRBs, which increase the known redshifts of long GRBs with plateaus by 94%, a significant milestone for enhancing GRB population studies that require large samples with redshift.

Neutron star mergers as the dominant contributor to the production of heavy r-process elements

ABSTRACT The discovery of the radioactively powered kilonova AT2017gfo, associated with the short-duration gamma-ray burst GRB 170817A and the gravitational wave source GW170817, has provided the first direct evidence supporting binary neutron star mergers as crucial astrophysical sites for the synthesis of heavy elements beyond iron through r-process nucleosysthesis in the universe. However, recent identification of kilonovae following long-duration gamma-ray bursts, such as GRB 211211A and GRB 230307A, has sparked discussions about the potential of neutron star–white dwarf mergers to also produce neutron-rich ejecta and contribute to the production of heavy r-process elements. In this work, we estimate the contribution of binary neutron star mergers to the total mass of r-process elements in the Milky Way and investigate the possibility of neutron star–white dwarf mergers as alternative astrophysical sites for r-process nucleosynthesis through an analysis of the total mass of the r-process elements in the Milky Way. Our results reveal that binary neutron star mergers can sufficiently account for the Galactic heavy r-process elements, suggesting that these events are the dominant contributor to the production of heavy r-process elements in the Milky Way. Considering the total mass of r-process elements in the Milky Way and the higher occurrence rate of neutron star–white dwarf mergers, it is unlikely that such mergers can produce a significant amount of neutron-rich ejecta, with the generated mass of r-process elements being lower than $0.005\, {\rm M}_{\odot }$.

Prospects for detecting neutron star–white dwarf mergers with decihertz gravitational-wave observatories

ABSTRACT Based on different neutron star–white dwarf (NS–WD) population models, we investigate the prospects of gravitational-wave (GW) detections for NS–WD mergers, with the help of early warnings from two space-borne decihertz GW observatories, DO-Optimal and DECIGO. We not only give quick assessments of the GW detection rates for NS–WD mergers with the two decihertz GW detectors, but also report systematic analyses on the characteristics of GW-detectable merger events using the method of Fisher matrix. With a sufficient 1-d early-warning time, the yearly GW detection number for DO-Optimal is in the range of (1.5–1.9) × 103, while it is (3.3–4.6) × 104 for DECIGO. More importantly, our results show that most NS–WD mergers can be localized with an uncertainty of $\mathcal {O}(10^{-2})\, \mathrm{deg}^2$. Given the NS–WD merger as a possible origin for a peculiar long-duration gamma-ray burst, GRB 211211A, followed with kilonova-like emissions, we further suggest that the GW early-warning detection would allow future electromagnetic telescopes to get prepared to follow up transients after some special NS–WD mergers. Based on our analyses, we emphasize that such a feasible ‘wait-for’ pattern can help to firmly identify the origin of GRB 211211A-like events in the future and bring excellent opportunities for the multimessenger astronomy.

Collapsars as Sites of r-process Nucleosynthesis: Systematic Photometric Near-infrared Follow-up of Type Ic-BL Supernovae

One of the open questions following the discovery of GW170817 is whether neutron star (NS) mergers are the only astrophysical sites capable of producing r-process elements. Simulations have shown that 0.01–0.1 M ⊙ of r-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both NS mergers and collapsing massive stars associated with long-duration gamma-ray bursts (collapsars). The hallmark signature of r-process nucleosynthesis in the binary NS merger GW170817 was its long-lasting near-infrared (NIR) emission, thus motivating a systematic photometric study of the light curves of broad-lined stripped-envelope (Ic-BL) supernovae (SNe) associated with collapsars. We present the first systematic study of 25 SNe Ic-BL—including 18 observed with the Zwicky Transient Facility and 7 from the literature—in the optical/NIR bands to determine what quantity of r-process material, if any, is synthesized in these explosions. Using semi-analytic models designed to account for r-process production in SNe Ic-BL, we perform light curve fitting to derive constraints on the r-process mass for these SNe. We also perform independent light curve fits to models without the r-process. We find that the r-process-free models are a better fit to the light curves of the objects in our sample. Thus, we find no compelling evidence of r-process enrichment in any of our objects. Further high-cadence infrared photometric studies and nebular spectroscopic analysis would be sensitive to smaller quantities of r-process ejecta mass or indicate whether all collapsars are completely devoid of r-process nucleosynthesis.

The Progenitor and Central Engine of a Peculiar GRB 230307A

Recently, a lack of supernova-associated with long-duration gamma-ray burst (GRB 230307A) at such a low redshift z = 0.065, but associated with a possible kilonova emission, has attracted great attention. Its heavy element nucleosynthesis and the characteristic of soft X-ray emission suggest that the central engine of GRB 230307A is a magnetar that is originated from a binary compact star merger. The calculated lower value of ε ∼ 0.05 suggests that GRB 230307A seems to have an ambiguous progenitor. The lower value of f eff = 1.23 implies that GRB 230307A is not likely to be from the effect of “tip of iceberg.” We adopt the magnetar central engine model to fit the observed soft X-ray emission with varying efficiency and find that the parameter constraints of the magnetar falls into a reasonable range, i.e., B < 9.4 × 1015 G and P < 2.5 ms for Γsat = 103, and B < 3.6 × 1015 G and P < 1.05 ms for Γsat = 104. Whether the progenitor of GBR 230307A is from the mergers of neutron star–white dwarf (NS–WD) or neutron star–neutron star (NS–NS) remains unknown. The difference of GW radiation between NS–NS merger and NS–WD merger may be a probe to distinguish the progenitor of GRB 230307A-like events in the future.

A Photohadronic Interpretation of H.E.S.S. Afterglow Observations of GRB 221009A

The High Energy Stereoscopic System (H.E.S.S.) started observing the extremely powerful long-duration gamma-ray burst (GRB) GRB 221009A starting 53 hr after the triggering event. The H.E.S.S. collaboration carried out observations on 2022 October 11, 12, and 17 under poor atmospheric conditions, without detecting significant very-high-energy photons from the source and computed the upper limits of the fluxes for the different nights. We study these flux upper limits by using the photohadronic model and show that the interaction of high-energy protons with synchrotron seed photons in the forward-shock region of the GRB jet exhibits behavior compatible with the upper limits computed by the H.E.S.S. collaboration.

Self-consistent magnetohydrodynamic simulation of jet launching in a neutron star – white dwarf merger

The merger of a white dwarf (WD) and a neutron star (NS) is a relatively common event that produces an observable electromagnetic signal. Furthermore, the compactness of these stellar objects makes them an interesting candidate for gravitational wave (GW) astronomy, potentially being in the frequency range of LISA and other missions. To date, three-dimensional simulations of these mergers have not fully modeled the WD disruption or have used lower resolutions and have not included magnetic fields even though they potentially shape the evolution of the merger remnant. In this work, we simulated the merger of a 1.4 M⊙ NS with a 1 M⊙ carbon oxygen WD in the magnetohydrodynamic moving mesh code AREPO. We find that the disruption of the WD forms an accretion disk around the NS, and the subsequent accretion by the NS powers the launch of strongly magnetized, mildly relativistic jets perpendicular to the orbital plane. Although the exact properties of the jets could be altered by unresolved physics around the NS, the event could result in a transient with a larger luminosity than kilonovae. We discuss possible connections to fast blue optical transients (FBOTs) and long-duration gamma-ray bursts. We find that the frequency of GWs released during the merger is too high to be detectable by the LISA mission, but suitable for deci-hertz observatories such as LGWA, BBO, or DECIGO.

Investigating the Evolution of Amati Parameters with Redshift

Gamma-ray bursts (GRBs) are among the brightest objects in the Universe and, hence, can be observed up to a very high redshift. Properly calibrated empirical correlations between intensity and spectral correlations of GRBs can be used to estimate the cosmological parameters. However, the possibility of the evolution of GRBs with redshift is a long-standing puzzle. In this work, we used 162 long-duration GRBs to determine whether GRBs below and above a certain redshift have different properties. The GRBs are split into two groups, and we fit the Amati relation for each group separately. Our findings demonstrate that estimations of the Amati parameters for the two groups are substantially dissimilar. We perform simulations to investigate whether the selection effects could cause the difference. Our analysis shows that the differences may be either intrinsic or due to systematic errors in the data, and the selection effects are not their true origin. However, in-depth analysis with a new data set comprised of 119 long GRBs shows that intrinsic scatter may partly be responsible for such effects.

The jet composition of GRB 230307A: Poynting-flux-dominated outflow?

ABSTRACT The jet composition of GRB plays an important role in understanding the energy dissipation and radiation mechanisms in GRB physics, but it is poorly constrained from the observational data. Recently, an interesting long-duration GRB 230307A with redshift z = 0.065 has attracted great attention. The lack of detected thermal emission and mini-structure of prompt emission light curve of this burst suggest that the outflow is Poynting-flux-dominated and point towards the ICMART model. In this paper, we invoke two independent methods to investigate the jet composition of GRB 230307A. The high magnetization parameter (σ &amp;gt; 7 or ever large) for R0 = 1010 cm that is used to suppress thermal component, strongly suggests that a significant fraction of the outflow energy is likely in a Poynting flux entrained with the baryonic matter. Moreover, it is found that the radiation efficiency of this burst for typical values ϵe = 0.1 and ϵB = 0.01 can reach as high as $~50{{\ \rm per\ cent}}$ which disfavours the internal shock model, but is consistent with ICMART model. Finally, a possible unified picture to produce GRB 230307A originated from a compact star merger is also discussed.