Long-duration Gamma-ray Bursts Research Articles

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.

Black hole growths in gamma-ray bursts driven by the Blandford–Znajek mechanism

ABSTRACT The Blandford–Znajek (BZ) mechanism in stellar-mass black hole (BH) hyperaccretion systems is generally considered to power gamma-ray bursts (GRBs). Based on observational GRB data, we use the BZ mechanism driven by the BH hyperaccretion disc to investigate the evolution of the BH mass and spin after the jets break out from the progenitors. We find that the BH growths are almost independent of initial BH masses. Meanwhile, the BH growths will be more efficient with smaller initial spin parameters. We conclude that (i) the BZ mechanism is efficient for triggering BH growths for only 1 of 206 typical long-duration GRBs; (ii) the mean BH mass growths of ultra-long GRBs are marginal for all 7 samples collected; (iii) for the short-duration GRBs, the results that BHs show minimal growths is consistent with the mass supply limitation in the scenario of compact object mergers.

Solving the inverse cosmological calibration problem of gamma-ray bursts

ABSTRACT We have received a new physical characteristics fitting based on actual observational data from the Swift mission’s long-duration gamma-ray bursts (LGRBs). We considered such characteristics as the Amati parameters for linear correlation (Eiso–Ep,i) and the k-correction for gravitational lensing and Malmquist bias (GLMB) effect. We used the Pantheon SN Ia catalogue and the standard Lambda cold dark matter model with a fixed Hubble constant of H0 = 70 km s−1 Mpc−1 as the baseline for the Hubble function μ(z). In our paper, we formulated the inverse cosmological calibration problem (ICCP) in the non-parametric statistics framework. The ICCP involves fitting non-observable physical characteristics while assuming a fixed cosmological model. To solve this problem, we developed a new method that is resistant to non-Gaussian processes. This method is based on error propagation through the Monte Carlo method and the Theil–Sen method for linear regression estimate. We have demonstrated the stability and robustness of this assessment method. The parameter estimates are as follows: $a=0.92^{+0.12}_{-0.12}$, $b=50.32^{+0.33}_{-0.32}$ without considering the GLMB effect, and $a=0.63^{+0.13}_{-0.14}$, $b=50.12^{+0.33}_{-0.31}$, and $k=1.98^{+0.25}_{-0.24}$ with the effect included. The proposed method can be applied to any other calibration sample of known standard candles, a calibrated sample of LGRBs, and the Hubble function μ(z). In the future, the ICCP idea can be used as an alternative cosmological test for estimating cosmological parameters, including the GLMB effect, or even for the selection of models, providing new information about the Universe. This can be done by analysing the residual values of observational data within the Bayesian statistics paradigm.

Heavy-element production in a compact object merger observed by JWST.

The mergers of binary compact objects such as neutron stars and black holes are of central interest to several areas of astrophysics, including as the progenitors of gamma-ray bursts (GRBs)1, sources of high-frequency gravitational waves (GWs)2 and likely production sites for heavy-element nucleosynthesis by means of rapid neutron capture (the r-process)3. Here we present observations of the exceptionally bright GRB 230307A. We show that GRB 230307A belongs to the class of long-duration GRBs associated with compact object mergers4-6 and contains a kilonova similar to AT2017gfo, associated with the GW merger GW170817 (refs. 7-12). We obtained James Webb Space Telescope (JWST) mid-infrared imaging and spectroscopy 29 and 61 days after the burst. The spectroscopy shows an emission line at 2.15 microns, which we interpret as tellurium (atomic mass A = 130) and a very red source, emitting most of its light in the mid-infrared owing to the production of lanthanides. These observations demonstrate that nucleosynthesis in GRBs can create r-process elements across a broad atomic mass range and play a central role in heavy-element nucleosynthesis across the Universe.

AGILE Gamma-Ray Detection of the Exceptional GRB 221009A

Gamma-ray emission in the MeV–GeV range from explosive cosmic events is of invaluable relevance to understanding physical processes related to the formation of neutron stars and black holes. Here we report on the detection by the AGILE satellite in the MeV–GeV energy range of the remarkable long-duration gamma-ray burst GRB 221009A. The AGILE onboard detectors have good exposure to GRB 221009A during its initial crucial phases. Hard X-ray/MeV emission in the prompt phase lasted hundreds of seconds, with the brightest radiation being emitted between 200 and 300 s after the initial trigger. Very intense GeV gamma-ray emission is detected by AGILE in the prompt and early afterglow phase up to 10,000 s. Time-resolved spectral analysis shows time-variable MeV-peaked emission simultaneous with intense power-law GeV radiation that persists in the afterglow phase. The coexistence during the prompt phase of very intense MeV emission together with highly nonthermal and hardening GeV radiation is a remarkable feature of GRB 221009A. During the prompt phase, the event shows spectrally different MeV and GeV emissions that are most likely generated by physical mechanisms occurring in different locations. AGILE observations provide crucial flux and spectral gamma-ray information regarding the early phases of GRB 221009A during which emission in the TeV range was reported.

Hybrid Emission Modeling of GRB 221009A: Shedding Light on TeV Emission Origins in Long GRBs

Observations of long-duration gamma-ray bursts (GRBs) with TeV emission during their afterglow have been on the rise. Recently, GRB 221009A, the most energetic GRB ever observed, was detected by the Large High Altitude Air Shower Observatory experiment in the energy band 0.2–7 TeV. Here, we interpret its afterglow in the context of a hybrid model in which the TeV spectral component is explained by the proton-synchrotron process while the low-energy emission from optical to X-ray is due to synchrotron radiation from electrons. We constrained the model parameters using the observed optical, X-ray, and TeV data. By comparing the parameters of this burst and of GRB 190114C, we deduce that the VHE emission at energies ≥1 TeV in the GRB afterglow requires large explosion kinetic energy, E ≳ 1054 erg and a reasonable circumburst density, n ≳ 10 cm−3. This results in a small injection fraction of particles accelerated to a power law, ∼10−2. A significant fraction of shock energy must be allocated to a near equipartition magnetic field, ϵ B ∼ 10−1, while electrons should only carry a small fraction of this energy, ϵ e ∼ 10−3. Under these conditions required for a proton-synchrotron model, namely ϵ B ≫ ϵ e , the SSC component is substantially subdominant over proton-synchrotron as a source of TeV photons. These results lead us to suggest that proton-synchrotron process is a strong contender for the radiative mechanisms explaining GRB afterglows in the TeV band.

A Luminous Precursor in the Extremely Bright GRB 230307A

GRB 230307A is an extremely bright long-duration GRB with an observed gamma-ray fluence of ≳3 × 10−3 erg cm−2 (10–1000 keV), second only to GRB 221009A. Despite its long duration, it is possibly associated with a kilonova, thus resembling the case of GRB 211211A. In analogy with GRB 211211A, we distinguish three phases in the prompt gamma-ray emission of GRB 230307A: an initial short duration, spectrally soft emission; a main long duration, spectrally hard burst; and a temporally extended and spectrally soft tail. We interpret the initial soft pulse as a bright precursor to the main burst and compare its properties with models of precursors from compact binary mergers. We find that to explain the brightness of GRB 230307A, a magnetar-like (≳1015 G) magnetic field should be retained by the progenitor neutron star. Alternatively, in the postmerger scenario, the luminous precursor could point to the formation of a rapidly rotating massive neutron star.

The peak flux of GRB 221009A measured with GRBAlpha

Context. On 2022 October 9 the brightest gamma-ray burst (GRB) ever observed lit up the high-energy sky. It was detected by a multitude of instruments, attracting the close attention of the GRB community, and saturated many detectors. Aims. GRBAlpha, a nano-satellite with a form factor of a 1U CubeSat, detected this extraordinarily bright long-duration GRB, GRB 221009A, without saturation but affected by pile-up. We present light curves of the prompt emission in 13 energy bands, from 80 keV to 950 keV, and performed a spectral analysis to calculate the peak flux and peak isotropic-equivalent luminosity. Methods. Since the satellite’s attitude information is not available for the time of this GRB, more than 200 incident directions were probed in order to find the median luminosity and its systematic uncertainty. Results. We find that the peak flux in the 80 − 800 keV range (observer frame) was Fphp = 1300−200+1200 ph cm−2 s−1, or Fergp = 5.7−0.7+3.7 × 10−4 erg cm−2 s−1, and the fluence in the same energy range of the first GRB episode, which lasted 300 s and was observable by GRBAlpha, was S = 2.2−0.3+1.4 × 10−2 erg cm−2, or Sbol = 4.9−0.5+0.8 × 10−2 erg cm−2 for the extrapolated range of 0.9 − 8690 keV. We infer the isotropic-equivalent released energy of the first GRB episode to be Eisobol = 2.8−0.5+0.8 × 1054 erg in the 1 − 10 000 keV band (rest frame at z = 0.15). The peak isotropic-equivalent luminosity in the 92 − 920 keV range (rest frame) was Lisop = 3.7−0.5+2.5 × 1052 erg s−1, and the bolometric peak isotropic-equivalent luminosity was Lisop,bol = 8.4−1.5+2.5 × 1052 erg s−1 (4 s scale) in the 1 − 10 000 keV range (rest frame). The peak emitted energy is Ep∗ = Ep(1+z) = 1120 ± 470 keV. Our measurement of Lisop,bol is consistent with the Yonetoku relation. It is possible that, due to the spectral evolution of this GRB and the orientation of GRBAlpha at the peak time, the true values of peak flux, fluence, Liso, and Eiso are even higher.