GRB Jets Research Articles

Collapsar Gamma-Ray Bursts Grind Their Black Hole Spins to a Halt

The spin of a newly formed black hole (BH) at the center of a massive star evolves from its natal value due to two competing processes: accretion of gas angular momentum that increases the spin and extraction of BH angular momentum by outflows that decreases the spin. Ultimately, the final, equilibrium spin is set by a balance between both processes. In order for the BH to launch relativistic jets and power a γ-ray burst (GRB), the BH magnetic field needs to be dynamically important. Thus, we consider the case of a magnetically arrested disk (MAD) driving the spin evolution of the BH. By applying the semianalytic MAD BH spin evolution model of Lowell et al. to collapsars, we show that if the BH accretes ∼20% of its initial mass, its dimensionless spin inevitably reaches small values, a ≲ 0.2. For such spins, and for mass accretion rates inferred from collapsar simulations, we show that our semianalytic model reproduces the energetics of typical GRB jets, L jet ∼ 1050 erg s−1. We show that our semianalytic model reproduces the nearly constant power of typical GRB jets. If the MAD onset is delayed, this allows powerful jets at the high end of the GRB luminosity distribution, L jet ∼ 1052 erg s−1, but the final spin remains low, a ≲ 0.3. These results are consistent with the low spins inferred from gravitational wave detections of binary BH mergers. In a companion paper by Gottlieb et al., we use GRB observations to constrain the natal BH spin to be a ≃ 0.2.

GRB 080710: A narrow, structured jet showing a late, achromatic peak in the optical and infrared afterglow?

We present a possible theoretical interpretation of the observed afterglow emission of long gamma-ray burst GRB 080710. While its prompt GRB emission properties are normal, the afterglow light curves in the optical and infrared bands are exceptional in two respects. One is that the observed light curves of different wavelengths have maximum at the same time, and that the achromatic peak time, 2.2×103 s after the burst trigger, is about an order of magnitude later than typical events. The other is that the observed flux before the peak increases more slowly than theoretically expected so far. Assuming that the angular distribution of the outflow energy is top-hat or Gaussian-shaped, we calculate the observed light curves of the synchrotron emission from the relativistic jets and explore the model parameters that explain the observed data. It is found that a narrowly collimated Gaussian-shaped jet with large isotropic-equivalent energy is the most plausible model for reproducing the observed afterglow behavior. Namely, an off-axis afterglow scenario to the achromatic peak is unlikely. The inferred values of the opening angle and the isotropic equivalent energy of the jet are possibly similar to those of GRB 221009A, but the jet of GRB 080710 has a much smaller efficiency of the prompt gamma-ray emission. Our results indicate a greater diversity of the GRB jet properties than previously thought.

Understanding Binary Systems—a Comparison between COSMIC and MESA

We compare the evolution of binary systems evolved in the MESA stellar evolution code to those in the COSMIC population synthesis code. Our aim is to convey the robustness of the equations that model binary evolution in the COSMIC code, particularly for the cases of high mass stars with closely orbiting compact object companions. Our larger goal is to accurately model the rates of these systems, as they are promising candidates for the progenitor systems behind energetic, longer lasting, radio bright GRB jets. These systems also may be key contributors to the rates of binary black hole mergers throughout our Universe.

TeV emission as an independent probe of GRB jet properties

TeV emission as an independent probe of GRB jet properties

Neutrino search from γ-ray bursts during the prompt and X-ray afterglow phases using 10 years of IceCube public data

Neutrino emission from γ-ray bursts (GRBs) has been sought for a long time, and stringent limits on the most accredited GRB emission models have been obtained from IceCube. Multiwavelength GRB observations of the last decades have improved our knowledge of the GRB emission parameters, such as the Lorentz factor and the luminosity, which can vary from one GRB to another by several orders of magnitude. Empirical correlations among such parameters have been identified during the prompt phase, with direct implications on GRB models. In this work, we use the PSLab open-access code, developed for IceCube data analyses, to search for individual neutrino emission from the prompt and afterglow phases of selected GRBs, and for stacking emission from the ensemble of such GRBs. For the afterglow phase, we focus on GRBs with X-ray flares and plateaus in particular. While past stacking searches assumed the same GRB fluence at Earth, we present a stacking scheme based on physically motivated GRB weights. Moreover, we conceive a new methodology for the prompt phase that uses the empirical correlations to infer the GRB luminosity and Lorentz factor, when redshift measurements are not available. We do not observe any significant neutrino excess. Hence, we set constraints on the GRB neutrino fluxes and on relevant GRB parameters, including the magnetic field in the jet. Notably, the baryon loading is found to be less than ten for typical GRB jets.

Deciphering the ∼18 TeV Photons from GRB 221009A

On 2022 October 9, an extremely powerful gamma-ray burst, GRB 221009A, was detected by several instruments. Despite being obstructed by the Milky Way galaxy, its afterglow outburst outshone all other GRBs seen before. LHAASO detected several thousand very high energy photons extending up to 18 TeV. Detection of such energetic photons is unexpected due to the large opacity of the universe. It is possible that in the afterglow epoch, the intrinsic very high energy photon flux from the source might have increased manifolds, which could compensate for the attenuation by pair production with the extragalactic background light. We propose such a scenario and show that very high energy photons can be observed on the Earth from the interaction of very high energy protons with the seed synchrotron photons in the external forward shock region of the GRB jet.

Constraints on the Physics of the Prompt Emission from Distant and Energetic Gamma-Ray Burst GRB 220101A

The emission region of γ-ray bursts (GRBs) is poorly constrained. The uncertainty on the size of the dissipation site spans over 4 orders of magnitude (1012–1017 cm) depending on the unknown energy composition of the GRB jets. The joint multiband analysis from soft X-rays to high energies (up to ∼1 GeV) of one of the most energetic and distant GRBs, GRB 220101A (z = 4.618), allows us to make an accurate distinction between prompt and early afterglow emissions. The enormous amount of energy released by GRB 220101A (E iso ≈ 3 × 1054 erg) and the spectral cutoff at MeV observed in the prompt emission spectrum constrain the parameter space of the GRB dissipation site. We put stringent constraints on the prompt emission site, requiring 700 < Γ0 < 1160 and R γ ∼ 4.5 × 1013 cm. Our findings further highlight the difficulty of finding a simple self-consistent picture in the electron–synchrotron scenario, favoring instead a proton–synchrotron model, which is also consistent with the observed spectral shape. Deeper measurements of the time variability of GRBs, together with accurate high-energy observations (MeV–GeV), would unveil the nature of the prompt emission.

Gravitational waves from the propagation of long gamma-ray burst jets

ABSTRACT Gamma-ray bursts (GRBs) are produced during the propagation of ultra-relativistic jets. It is challenging to study the jet close to the central source, due to the high opacity of the medium. In this paper, we present numerical simulations of relativistic jets propagating through a massive, stripped envelope star associated to long GRBs, breaking out of the star and accelerating into the circumstellar medium. We compute the gravitational wave (GW) signal resulting from the propagation of the jet through the star and the circumstellar medium. We show that key parameters of the jet propagation can be directly determined by the GW signal. The signal presents a first peak corresponding to the jet duration and a second peak which corresponds to the break-out time for an observer located close to the jet axis (which in turn depends on the stellar size), or to much larger times (corresponding to the end of the acceleration phase) for off-axis observers. We also show that the slope of the GW signal before and around the first peak tracks the jet luminosity history and the structure of the progenitor star. The amplitude of the GW signal is h+D ∼ hundreds to several thousands cm. Although this signal, for extragalactic sources, is outside the range of detectability of current GW detectors, it can be detected by future instruments as BBO, DECIGO, and ALIA. Our results illustrate that future detections of GW associated to GRB jets may represent a revolution in our understanding of this phenomenon.

GeV Signatures of Short Gamma-Ray Bursts in Active Galactic Nuclei

The joint detection of gravitational waves and the gamma-ray counterpart of a binary neutron star merger event, GW170817, unambiguously validates the connection between short gamma-ray bursts and compact binary object (CBO) mergers. We focus on a special scenario where short gamma-ray bursts produced by CBO mergers are embedded in disks of active galactic nuclei (AGNs), and we investigate the γ-ray emission produced in the internal dissipation region via synchrotron, synchrotron self-Compton, and external inverse Compton (EIC) processes. In this scenario, isotropic thermal photons from the AGN disks contribute to the EIC component. We show that a low-density cavity can be formed in the migration traps, leading to the embedded mergers producing successful GRB jets. We find that the EIC component would dominate the GeV emission for typical CBO mergers with an isotropic-equivalent luminosity of L j,iso = 1048.5 erg s−1 that are located close to the central supermassive black hole. Considering a long-lasting jet of duration T dur ∼ 102–103 s, we find that the future Cherenkov Telescope Array (CTA) will be able to detect its 25–100 GeV emission out to a redshift z = 1.0. In the optimistic case, it is possible to detect the on-axis extended emission simultaneously with GWs within one decade using MAGIC, H.E.S.S., VERITAS, CTA, and LHAASO-WCDA. Early diagnosis of prompt emissions with Fermi-GBM and HAWC can provide valuable directional information for the follow-up observations.

Very High-energy Afterglow Emission of GRB 190829A: Evidence for Its Hadronic Origin?

The detection of multi-TeV gamma rays from the afterglow phase of GRB 190829A by the High Energy Stereoscopic System telescope is an addition to the already existing list of two GRBs observed in very high-energy (VHE) gamma rays in recent years. Jets of blazars and GRBs have many similarities and the photohadronic model is very successful in explaining the VHE gamma-ray spectra from the high-energy blazars. Recently, the photohadronic model has been successfully applied to study the sub-TeV gamma rays from the afterglow phases of GRB 180720B and GRB 190114C. We employed this model again to explain the VHE spectra observed for the two consecutive nights from GRB 190829A. We show that the spectra of GRB 190829A can be due to the interactions of high-energy protons with the synchrotron self-Compton photons in the forward shock region of the GRB jet, similar to the low emission state of the VHE flaring events of high-energy blazars. We speculate that, if in the future, it is possible to observe the VHE gamma-ray spectra from nearby GRBs in their afterglow phases, then some of them could only be explained by employing two different spectral indices. If confirmed, such VHE spectra could be interpreted as a result of the interactions of the high-energy protons with the photons, both from the synchrotron background and the synchrotron self-Compton background in the forward shock region.

Fast Particle Acceleration in Three-dimensional Relativistic Reconnection

Abstract Magnetic reconnection is invoked as one of the primary mechanisms to produce energetic particles. We employ large-scale 3D particle-in-cell simulations of reconnection in magnetically dominated (σ = 10) pair plasmas to study the energization physics of high-energy particles. We identify an acceleration mechanism that only operates in 3D. For weak guide fields, 3D plasmoids/flux ropes extend along the z-direction of the electric current for a length comparable to their cross-sectional radius. Unlike in 2D simulations, where particles are buried in plasmoids, in 3D we find that a fraction of particles with γ ≳ 3σ can escape from plasmoids by moving along z, and so they can experience the large-scale fields in the upstream region. These “free” particles preferentially move in z along Speiser-like orbits sampling both sides of the layer and are accelerated linearly in time—their Lorentz factor scales as γ ∝ t, in contrast to γ ∝ t in 2D. The energy gain rate approaches ∼eE rec c, where E rec ≃ 0.1B 0 is the reconnection electric field and B 0 the upstream magnetic field. The spectrum of free particles is hard, dN free / d γ ∝ γ − 1.5 , contains ∼20% of the dissipated magnetic energy independently of domain size, and extends up to a cutoff energy scaling linearly with box size. Our results demonstrate that relativistic reconnection in GRB and AGN jets may be a promising mechanism for generating ultra-high-energy cosmic rays.

Acceleration of ultrahigh-energy cosmic rays in the early afterglows of gamma-ray bursts: Concurrence of jet dynamics and wave-particle interactions

The origin of ultra-high energy cosmic rays $($UHECRs$)$ remains a mystery. It has been suggested that UHECRs can be produced by the stochastic acceleration in relativistic jets of gamma-ray bursts $($GRBs$)$ at the early afterglow phase. Here, we develop a time-dependent model for proton energization by cascading compressible waves in GRB jets considering the concurrent effect of the jet's dynamics and the mutual interactions between turbulent waves and particles. Considering fast mode of magnetosonic wave as the dominant particle scatterer and assuming interstellar medium $($ISM$)$ for the circumburst environment, our numerical results suggest that protons can be accelerated up to $\textrm{10}^{\textrm{19}}\,$eV during the early afterglow. An estimation shows ultra-high energy nuclei can easily survive photodisintegration in the external shocks in most cases, thus allowing the acceleration of $\textrm{10}^{\textrm{20}}\,$eV cosmic-ray nuclei in the proposed frame. The spectral slope can be as hard as $\textrm{d}\textit{N}/\textrm{d}\textit{E} \propto \textit{E}^\textrm{0}$, which is consistent with the requirement for the interpretation of intermediate-mass composition of UHECR as measured by the Pierre Auger Observatory.

GRB 101225A as Orphan Dipole Radiation of a Newborn Magnetar with Precession Rotation in an Off-axis Gamma-ray Burst

Abstract The unusual multiwavelength lightcurves of GRB 101225A are revisited by assuming that they are from an off-axis GRB powered by a newborn magnetar. We show that GRB 101225A’s optical afterglow lightcurve is fitted with the forward shock model by parameterizing its jet structure as a Gaussian function with a half-opening angle of the jet core as 1.°67. The derived initial Lorentz factor (Γ0) is 120, and the viewing angle to the jet axis is θ v = 3.°7. Tentative QPO signatures of P = 488 s and P = 250 ∼ 300 s are found with a confidence level of 90% by analyzing its X-ray flares observed in the time interval of [4900, 7500] s. Its global gamma-ray/X-ray lightcurve and the QPO signatures are represented with the magnetar dipole radiation (DR) model by considering the magnetar precession motion, assuming that the magnetar spindown is dominated by GW emission. The bulk Lorentz factor of the DR ejecta is limited to 8, being much lower than Γ0. Comparing GRB 101225A with the extremely off-axis GRB 170817A, we suspect that the nature of the two-component jet in GRB 170817A is a combination of a co-axial GRB jet and a DR ejecta. GRB 101225A would be among the brightest ones of the CDF-S XT2-like X-ray transient population driven by newborn magnetars. A discussion of the detectability of its gravitational wave emission is also presented.

How to probe GRB jets

Abstract Núria Jordana-Mitjans explains how she uses linear polarization to unveil strong magnetic fields in gamma-ray burst outflows

High-energy Neutrinos from Choked Gamma-Ray Bursts in Active Galactic Nucleus Accretion Disks

Abstract Both long-duration gamma-ray bursts (LGRBs) from the core collapse of massive stars and short-duration GRBs (SGRBs) from mergers of a binary neutron star or a neutron star–black hole are expected to occur in the accretion disk of active galactic nuclei (AGNs). We show that GRB jets embedded in the migration traps of AGN disks are promised to be choked by the dense disk material. Efficient shock acceleration of cosmic rays at the reverse shock is expected, and high-energy neutrinos would be produced. We find that these sources can effectively produce detectable TeV–PeV neutrinos through pγ interactions. From a choked LGRB jet with isotropic equivalent energy of 1053 erg at 100 Mpc, one expects ∼2(7) neutrino events detectable by IceCube (IceCube-Gen2). The contribution from choked LGRBs to the observed diffuse neutrino background depends on the unknown local event rate density of these GRBs in AGN disks. For example, if the local event rate density of choked LGRBs in an AGN disk is ∼5% that of low-luminosity GRBs (∼10 Gpc−3 yr−1), the neutrinos from these events would contribute to ∼10% of the observed diffuse neutrino background. Choked SGRBs in AGN disks are potential sources for future joint electromagnetic, neutrino, and gravitational wave multimessenger observations.

The Fate of the Merger Remnant in GW170817 and Its Imprint on the Jet Structure

Abstract The first neutron star binary merger detected in gravitational waves, GW170817, and the subsequent detection of its emission across the electromagnetic spectrum showed that these systems are viable progenitors of short γ-ray bursts (sGRB). The afterglow signal of GW170817 has been found to be consistent with a structured GRB jet seen off-axis, requiring significant amounts of relativistic material at large angles. This trait can be attributed to the interaction of the relativistic jet with the external wind medium. Here we perform numerical simulations of relativistic jets interacting with realistic wind environments in order to explore how the properties of the wind and central engine affect the structure of successful jets. We find that the angular energy distribution of the jet depends primarily on the ratio between the lifetime of the jet and the time it takes the merger remnant to collapse. We make use of these simulations to constrain the time it took for the merger remnant in GW170817 to collapse into a black hole based on the angular structure of the jet as inferred from afterglow observations. We conclude that the lifetime of the merger remnant in GW170817 was ≈1–1.7 s, which, after collapse, triggered the formation of the jet.

GRB 160625B: Evidence for a Gaussian-shaped Jet

Abstract We present multiwavelength modeling of the afterglow from the long γ-ray burst (GRB) 160625B using Markov Chain Monte Carlo techniques of the afterglowpy Python package. GRB 160625B is an extremely bright burst with a rich set of observations spanning from radio to γ-ray frequencies. These observations range from ∼0.1 days to &gt;1000 days, thus making this event extremely well suited to such modeling. In this work we compare top-hat and Gaussian jet structure types in order to find best-fit values for the GRB jet collimation angle, viewing angle, and other physical parameters. We find that a Gaussian-shaped jet is preferred (2.7σ–5.3σ) over the traditional top-hat model. Our estimate for the opening angle of the burst ranges from 1.°26 to 3.°90, depending on jet-shape model. We also discuss the implications that assumptions on jet shape, viewing angle, and particularly the participation a fraction of electrons have on the final estimation of GRB intrinsic energy release and the resulting energy budget of the relativistic outflow. Most notably, allowing the participation fraction to vary results in an estimated total relativistic energy of ∼1053 erg. This is two orders of magnitude higher than when the total fraction is assumed to be unity; thus, this parameter has strong relevance for placing constraints on long GRB central engines, details of the circumburst media, and host environment.

Jet propagation in expanding medium for gamma-ray bursts

ABSTRACT The binary neutron star (BNS) merger event GW170817 clearly shows that a BNS merger launches a short gamma-ray burst (sGRB) jet. Unlike collapsars, where the ambient medium is static, in BNS mergers the jet propagates through the merger ejecta that is expanding outward at substantial velocities (∼0.2c). Here, we present semi-analytical and analytical models to solve the propagation of GRB jets through their surrounding media. These models improve our previous model by including the jet collimation by the cocoon self-consistently. We also perform a series of 2D numerical simulations of jet propagation in BNS mergers and in collapsars to test our models. Our models are consistent with numerical simulations in every aspect (the jet head radius, the cocoon’s lateral width, the jet opening angle including collimation, the cocoon pressure, and the jet–cocoon morphology). The energy composition of the cocoon is found to be different depending on whether the ambient medium is expanding or not; in the case of BNS merger jets, the cocoon energy is dominated by kinetic energy, while it is dominated by internal energy in collapsars. Our model will be useful for estimating electromagnetic counterparts to gravitational waves.

Extracting high-level information from gamma-ray burst supernova spectra

ABSTRACT Radiation transport codes are often used in astrophysics to construct spectral models. In this work, we demonstrate how producing these models for a time series of data can provide unique information about supernovae (SNe). Unlike previous work, we specifically concentrate on the method for obtaining the best synthetic spectral fits, and the errors associated with the preferred model parameters. We demonstrate how varying the ejecta mass, bolometric luminosity (Lbol) and photospheric velocity (vph), affects the outcome of the synthetic spectra. As an example we analyse the photospheric phase spectra of the GRB-SN 2016jca. It is found that for most epochs (where the afterglow subtraction is small) the error on Lbol and vph was ∼5 per cent. The uncertainty on ejecta mass and Ekin was found to be ∼20 per cent, although this can be expected to dramatically decrease if models of nebular phase data can be simultaneously produced. We also demonstrate how varying the elemental abundance in the ejecta can produce better synthetic spectral fits. In the case of SN 2016jca it is found that a decreasing 56Ni abundance as a function of decreasing velocity produces the best-fitting models. This could be the case if the 56Ni was synthesized at the side of the GRB jet, or dredged up from the centre of the explosion. The work presented here can be used as a guideline for future studies on SNe which use the same or similar radiation transfer code.

Erratum: On the existence of a luminosity threshold of GRB jets in massive stars

Erratum: On the existence of a luminosity threshold of GRB jets in massive stars