Engine Of Gamma-ray Bursts Research Articles

Probing Blackbody Components in Gamma-Ray Bursts from Black Hole Neutrino-dominated Accretion Flows

A stellar-mass black hole (BH) surrounded by a neutrino-dominated accretion flow (NDAF) is generally considered to be the central engine of gamma-ray bursts (GRBs). Neutrinos escaping from the disk will annihilate outside the disk to produce the fireball that could power GRBs with blackbody (BB) components. The initial GRB jet power and fireball launch radius are related to the annihilation luminosity and annihilation height of the NDAFs, respectively. In this paper, we collect seven GRBs with known redshifts and identified BB components to test whether the NDAF model works. We find that, in most cases, the values of the accretion rates and the central BH properties are all in the reasonable range, suggesting that these BB components indeed originate from the neutrino annihilation process.

Neutrino Oscillation Effects on the Luminosity of Neutrino-dominated Accretion Flows around Black Holes

A stellar-mass black hole surrounded by the neutrino-dominated accretion flow (NDAF) has been proposed as the central engine of gamma-ray bursts. In this work, we investigate the neutrino/antineutrino luminosity and annihilation luminosity from the NDAF and pair annihilation model, taking into account the neutrino oscillation above the accretion disk. The disk's hydrodynamical properties are modeled using an empirical solution previously derived with boundary conditions, including the effect of electron degeneracy and neutrino trapping. Our key parameters are the mass accretion rate and the black hole spin given in the ranges of Ṁ=0.1 –10 M ⊙ s−1 and 0 ≤ a < 1, respectively. Without neutrino oscillation, the obtained neutrino/antineutrino luminosity is ∼1051–1053 erg s−1, while the neutrino annihilation luminosity is found to be ∼1046–1051 erg s−1. In the presence of neutrino oscillation in the vacuum limit, the electron neutrino annihilation luminosity decreases by ≲22% through the flavor transformation, while the muon and tau neutrino luminosity can increase by up to ∼45% and 60%, respectively. As a result, the total annihilation luminosity can be reduced by up to ∼19% due to the oscillation process above the disk. Finally, we also investigate the case whereby the charge-conjugation-parity-symmetry-violating (CP-violating) phase is changed from δCP = 0° to 245°. However, our results reveal that the CP-violating phase has minimal impact on the neutrino annihilation luminosity.

GRB 211211A: The Case for an Engine-powered over r-process-powered Blue Kilonova

The recent gamma-ray burst (GRB) GRB 211211A provides the earliest (∼5 hr) data of a kilonova (KN) event, displaying bright (∼1042 erg s−1) and blue early emission. Previously, this KN was explained using simplistic multicomponent fitting methods. Here, in order to understand the physical origin of the KN emission in GRB 211211A, we employ an analytic multizone model for r-process-powered KNe. We find that r-process-powered KN models alone cannot explain the fast temporal evolution and the spectral energy distribution (SED) of the observed emission. Specifically, (i) r-process models require high ejecta mass to match early luminosity, which overpredicts late-time emission, while (ii) red KN models that reproduce late emission underpredict early luminosity. We propose an alternative scenario involving early contributions from the GRB central engine via a late low-power jet, consistent with plateau emission in short GRBs and GeV emission detected by Fermi-LAT at ∼104 s after GRB 211211A. Such late central engine activity, with an energy budget of ∼a few percent of that of the prompt jet, combined with a single red KN ejecta component, can naturally explain the light curve and SED of the observed emission, with the late-jet–ejecta interaction reproducing the early blue emission and r-process heating reproducing the late red emission. This supports claims that late low-power engine activity after prompt emission may be common. We encourage early follow-up observations of future nearby GRBs and compact binary merger events to reveal more about the central engine of GRBs and r-process events.

Distributions of Energy, Luminosity, Duration, and Waiting Times of Gamma-Ray Burst Pulses with Known Redshift Detected by Fermi/GBM

Discovered more than 50 years ago, gamma-ray burst (GRB) prompt emission remains the most puzzling aspect of GRB physics. Its complex and irregular nature should reveal how newborn GRB engines release their energy. In this respect, the possibility that GRB engines could operate as self-organized critical (SOC) systems has been put forward. Here, we present the energy, luminosity, waiting time, and duration distributions of individual pulses of GRBs with known redshift detected by the Fermi Gamma-ray Burst Monitor. This is the first study of this kind in which selection effects are accounted for. The compatibility of our results with the framework of SOC theory is discussed. We found evidence for an intrinsic break in the power-law models that describe the energy and the luminosity distributions.

Black Hole Activity Imprints on the Internal Plateau and the Subsequent Sharp Decay

A stellar-mass black hole (BH) or a millisecond magnetar is believed to be born as the central engine of gamma-ray bursts (GRBs). The presence of plateaus in the X-ray extended emission or afterglow of GRBs is widely accepted as an indicator of a magnetar central engine, particularly those with a sharp decay (faster than t −3), the so-called internal plateau. However, an alternative model, by taking the evolution of the magnetic flux at the BH horizon into account, suggests that an internal plateau can also arise from a jet powered by the Blandford–Znajek (BZ) mechanism (hereafter, a BZ jet). In this study, we propose that a precessional BZ jet would manifest a quasiperiodic oscillation (QPO) signature on the internal plateau and the subsequent sharp decay. Such lightcurves cannot be readily explained by the activity of a short-lived, supermassive magnetar, thus favoring a Kerr BH as the central engine. The X-ray afterglow of GRB 050904, comprising nine flares, is characterized by a QPO-modulated plateau and sharp decay, which can be well reproduced by a precessional BZ jet model. Therefore, one potential clue for distinguishing between these two engines lies in whether the QPO signature is present throughout the entire plateau and the subsequent sharp decay, as the magnetar scenario suggests a collapse at the end of the plateau.

The Fundamental Plane of GRBs

Gamma-ray bursts (GRBs) exhibit powerful radiation and relativistic jets similar to blazars. However, the central engine of GRBs remains unknown. In this paper, we use the fundamental plane to analyze a sample of GRBs with measured mass. We extend, over ∼12 orders of magnitude, the correlation analysis and fundamental plane with a sample of X-ray binaries, active galactic nuclei (AGNs) including blazars, and brightest cluster galaxies. The fundamental plane of our de-beamed sample, with a measured mass ( = (0.60 ± 0.03)+(0.78 ± 0.02)+7.23 ± 0.95), closely aligns with the findings of previous work on AGNs. This finding suggests that GRBs adhere to the fundamental plane of AGNs and supports the theory proposed in previous work that the central engine of GRBs may be black holes. This observation provides a plausible explanation for the striking similarities between GRBs and AGNs.

Magnetar Flare-driven Bumpy Declining Light Curves in Hydrogen-poor Superluminous Supernovae

Recent observations indicate that hydrogen-poor superluminous supernovae (SLSNe-I) often display bumpy declining light curves. However, the cause of these undulations remains unclear. In this paper, we have improved the magnetar model, which includes flare activities. We present a systematic analysis of a well-observed SLSN-I sample with bumpy light curves in their late phase. These SLSNe-I were identified from multiple transient surveys, such as the Pan-STARRS1 Medium Deep Survey and the Zwicky Transient Facility. Our study provides a set of magnetar-powered model light curve fits for five SLSNe-I, which accurately reproduce the observed light curves using reasonable physical parameters. By extracting the essential characteristics of both the explosions and central engines, these fits provide valuable insights into investigating their potential association with gamma-ray burst engines. We found that the SLSN flares tend to be the dim and long extension of GRB flares in the peak luminosity versus peak time plane. Conducting large-scale, high-cadence surveys in the near future could enhance our comprehension of both SLSN undulating light curve properties and their potential relationship with GRBs by modeling their light curve characteristics.

Scale invariance in x-ray flares of gamma-ray bursts

X-ray flares are generally believed to be produced by the reactivation of the central engine, and may have the same energy dissipation mechanism as the prompt emission of gamma-ray bursts (GRBs). X-ray flares can therefore provide important clues to understanding the nature of the central engines of GRBs. In this work, we study for the first time the physical connection between differential size and return distributions of X-ray flares of GRBs with known redshifts. We find that the differential distributions of duration, energy, and waiting time can be well fitted by a power-law function. In particular, the distributions for the differences of durations, energies, and waiting times at different times (i.e., the return distributions) well follow a $q$-Gaussian form. The $q$ values in the $q$-Gaussian distributions remain nearly steady for different temporal interval scales, implying a scale-invariant structure of GRB X-ray flares. Moreover, we verify that the $q$ parameters are related to the power-law indices $\alpha$ of the differential size distributions, characterized as $q=(\alpha+2)/\alpha$. These statistical features can be well explained within the physical framework of a self-organizing criticality system.

Magnetically arrested accretion disks launching structured jets in application to GRB and AGN engines

Context. We explore the formation, energetics, and geometry of relativistic jets along with the variability of their central engine. We study both fast and slowly rotating black holes and address our simulations to active galaxy centers as well as gamma ray burst engines. Aims. The structured jets are postulated to account for emission properties of high energy sources across the mass scale, launched from stellar mass black holes in gamma ray bursts (GRBs) and from supermassive black holes in active galactic nuclei (AGNs). Their active cores contain magnetized accretion disks and the rotation of the Kerr black hole provides a mechanism for launching relativistic jets. This process works most effectively if the mode of accretion turns out to be magnetically arrested. In this mode, the modulation of jets launched from the engine is related to internal instabilities in the accretion flow that operate on smallest time and spatial scales. As these scales are related to the light-crossing time and the black hole gravitational radius, the universal model of jet-disk connection is expected to scale with the mass of the black hole. Methods. We investigated the jet-disk connection by means of 3D general relativistic magneto-hydrodynamical simulations of the magnetically arrested disk accretion in Kerr geometry. We also quantified the variability of the disk by means of a Fourier analysis. Results. We found that the system evolution is governed by the physical parameters of the engine, such as the black hole spin and disk size, as well as disk magnetization, and we applied our scenarios to typical types of sources in AGN and GRB classes. We found that the magnetically arrested disk (MAD) scenario is applicable to AGN engines and supports persistent jet emissions. It can also be applied to GRBs, as it gives the variability pattern roughly consistent with observations. However, in some cases, strong magnetic fields may lead to jet quenching, and this effect is found to be important mainly for GRB jets. We speculate that it may be related to the strength of magnetically driven winds from the GRB engines.

“Super-kilonovae” from Massive Collapsars as Signatures of Black Hole Birth in the Pair-instability Mass Gap

The core collapse of rapidly rotating massive ∼ 10M ⊙ stars (“collapsars”), and the resulting formation of hyperaccreting black holes, comprise a leading model for the central engines of long-duration gamma-ray bursts (GRBs) and promising sources of r-process nucleosynthesis. Here, we explore the signatures of collapsars from progenitors with helium cores ≳ 130M ⊙ above the pair-instability mass gap. While the rapid collapse to a black hole likely precludes prompt explosions in these systems, we demonstrate that disk outflows can generate a large quantity (up to ≳ 50M ⊙) of ejecta, comprised of ≳ 5–10M ⊙ in r-process elements and ∼ 0.1–1M ⊙ of 56Ni, expanding at velocities ∼0.1 c. Radioactive heating of the disk wind ejecta powers an optical/IR transient, with a characteristic luminosity ∼ 1042 erg s−1 and a spectral peak in the near-IR (due to the high optical/UV opacities of lanthanide elements), similar to kilonovae from neutron star mergers, but with longer durations ≳1 month. These “super-kilonovae” (superKNe) herald the birth of massive black holes ≳ 60M ⊙, which—as a result of disk wind mass loss—can populate the pair-instability mass gap “from above,” and could potentially create the binary components of GW190521. SuperKNe could be discovered via wide-field surveys, such as those planned with the Roman Space Telescope, or via late-time IR follow-up observations of extremely energetic GRBs. Multiband gravitational waves of ∼ 0.1–50 Hz from nonaxisymmetric instabilities in self-gravitating massive collapsar disks are potentially detectable by proposed observatories out to hundreds of Mpc; in contrast to the “chirp” from binary mergers, the collapsar gravitational-wave signal decreases in frequency as the disk radius grows (“sad trombone”).

Modeling the Gamma-Ray Burst Jet Properties with 3D General Relativistic Simulations of Magnetically Arrested Accretion Flows

We investigate the dependence of the gamma-ray burst (GRB) jet structure and its evolution on the properties of the accreting torus in the central engine. Our models numerically evolve the accretion disk around a Kerr black hole using three-dimensional general relativistic magnetohydrodynamic simulations. We use two different analytical hydrodynamical models of the accretion disk, based on the Fishbone–Moncrief and Chakrabarti solutions, as our initial states for the structure of the collapsar disk and the remnant after a binary neutron star (BNS) merger, respectively. We impose poloidal magnetic fields of two different geometries upon the initial stable solutions. We study the formation and evolution of the magnetically arrested disk state and its effect on the properties of the emitted jet. The jets produced in our models are structured and have a relatively hollow core and reach higher Lorentz factors at an angle ≳9° from the axis. The jet in our short GRB model has an opening angle of up to ∼25° while our long GRB engine produces a narrower jet, of up to ∼11°. We also study the time variability of the jets and provide an estimate of the minimum variability timescale in our models. The application of our models to the GRB jets in the BNS postmerger system and to the ultrarelativistic jets launched from collapsing stars are briefly discussed.

Relativistic global solutions of neutrino-dominated accretion flows with magnetic coupling

ABSTRACT A Kerr black hole (BH) surrounded by a neutrino-dominated accretion flow (NDAF) is one of plausible candidates of the central engine in gamma-ray bursts. The accretion material might inherit and restructure strong magnetic fields from the compact object mergers or massive collapsars. The magnetic coupling (MC) process between a rapid rotating BH and an accretion disc is one of possible magnetic configurations that transfers the energy and angular momentum from the BH to the disc. In this paper, we investigate one-dimensional global solutions of NDAFs with MC (MCNDAFs), taking into account general relativistic effects, detailed neutrino physics, different MC geometries, and reasonable nucleosynthesis processes. Six cases with different accretion rates and power-law indices of magnetic fields are presented and compared with NDAFs without MC. Our results indict that the MC process can prominently impact the structure, thermal properties, and microphysics of MCNDAFs, increase luminosities of neutrinos and their annihilations, result in the changing of radial distributions of nucleons, and push the region of heavy nuclei synthesis to a larger radius than counterparts in NDAFs.

X-ray flares raising upon magnetar plateau as an implication of a surrounding disk of newborn magnetized neutron star

The X-ray flares have usually been ascribed to long-lasting activities of the central engine of gamma-ray bursts (GRBs), e.g., fallback accretion. The GRB X-ray plateaus, however, favor a millisecond magnetar central engine. The fallback accretion can be significantly suppressed due to the propeller effect of a magnetar. Therefore, if the propeller regime cannot resist the mass flow onto the surface of the magnetar efficiently, the X-ray flares raising upon the magnetar plateau would be expected. In this work, such peculiar cases are connected to the accretion process of the magnetars, and an implication for magnetar-disc structure is given. We investigate the repeated accretion process with multi-flare GRB 050730, and give a discussion for the accretion-induced variation of the magnetic field in GRB 111209A. Two or more flares exhibit in the GRB 050730, 060607A and 140304A; by adopting magnetar mass M = 1.4 M ⊙ and radius R = 12 km, the average mass flow rates of the corresponding surrounding disk are 3.53 × 10−4 M ⊙ s−1, 4.23 × 10−4 M ⊙ s−1, and 4.33 × 10−4 M ⊙ s−1, and the corresponding average sizes of the magnetosphere are 5.01 × 106 cm, 6.45 × 106 cm, and 1.09 × 107 cm, respectively. A statistic analysis that contains eight GRBs within 12 flares shows that the total mass loading in single flare is ∼ 2 × 10−5 M ⊙. In the lost mass of a disk, there are about 0.1% used to feed a collimated jet.

Evidence of X-Ray Plateaus Driven by the Magnetar Spindown Winds in Gamma-Ray Burst Afterglows

Abstract The central engine of gamma-ray bursts (GRBs) remains an open and cutting-edge topic in the era of multimessenger astrophysics. X-ray plateaus appear in some GRB afterglows, which are widely considered to originate from the spindown of magnetars. According to the stable magnetar scenario of GRBs, an X-ray plateau and a decay phase ∼t −2 should appear in X-ray afterglows. Meanwhile, the “normal” X-ray afterglow is produced by the external shock from a GRB fireball. We analyze the Neil Gehrels Swift GRB data, then find three gold samples that have an X-ray plateau and a decay phase ∼t −2 superimposed on the jet-driven normal component. Based on these features of the lightcurves, we argue that the magnetars should be the central engines of these three GRBs. Future joint multimessenger observations might further test this possibility, which can then be beneficial to constrain GRB physics.

Ultra-delayed Neutrino-driven Explosion of Rotating Massive-star Collapse

Abstract Long-term neutrino-radiation hydrodynamics simulations in full general relativity are performed for the collapse of rotating massive stars that are evolved from He-stars with initial masses of 20 and 32 M ⊙. It is shown that if the collapsing stellar core has sufficient angular momentum, the rotationally supported proto-neutron star (PNS) survives for seconds accompanying the formation of a massive torus of mass larger than 1 M ⊙. Subsequent mass accretion onto the central region produces a massive and compact central object, and eventually enhances the neutrino luminosity beyond 1053 erg s−1, resulting in a very delayed neutrino-driven explosion, in particular toward the polar direction. The kinetic energy of the explosion can be appreciably higher than 1052 erg for a massive progenitor star and compatible with that of energetic supernovae like broad-line type-Ic supernovae. By the subsequent accretion, the massive PNS collapses eventually into a rapidly spinning black hole, which could be a central engine for gamma-ray bursts if a massive torus surrounds it.

Constraining properties of GRB central engines with X-ray flares

ABSTRACT X-ray flares in gamma-ray bursts (GRBs) are believed to be generated by the late activities of central engine, and thus provide an useful tool to diagnose the properties of central objects. In this paper, we work on a GRB X-ray flare sample whose bulk Lorentz factors are constrained by two different methods and the jet opening angles are determined by the jet breaks in afterglow light curves. Considering a hyperaccreting stellar-mass black hole (BH) as the central engine of GRBs and the Blandford &amp; Znajek process (BZ) as the jet production mechanism, we constrain the parameters of central engine by using the X-ray flare data. We find that the BZ mechanism is so powerful making it possible to interpret both GRB prompt emissions and bright X-ray flares. The wind parameter (p) and accreted mass (Md) fall into reasonable ranges. Our result is also applied to GRB 170817A. The late X-ray flare in GRB 170817A, if it is true, might not be a BH origin.

Two X-Ray Plateaus of Gamma-Ray Bursts: Energy Injection from Nascent Magnetars with an Evolving Magnetic Inclination Angle

Abstract As the central engine of gamma-ray bursts (GRBs), a nascent magnetar can release its rotational energy by magnetic dipole radiation. This energy is injected into the outflow of the GRB jet. With the rapid evolution of the magnetic inclination angle, the perpendicular component of the magnetar’s surface magnetic field decreases rapidly. The energy loss rate of the magnetar is characterized by two different timescales. This peculiar two-stage energy injection model will affect the shape of the X-ray afterglow light curve. It is shown that a new shallow decay phase will emerge following the initial shallow decay phase. As a result, there could be two plateaus in the X-ray light curve. The effects of various physical parameters on the light curves are studied in detail. It is found that the observed afterglows of GRB 161217A and GRB 190114A, both having two plateaus in the light curve, can be well fitted in our framework.

Toward observing neutron star collapse with gravitational wave detectors

Gravitational waves from binary neutron star inspirals have been detected along with the electromagnetic transients coming from the aftermath of the merger in GW170817. However, much is still unknown about the post-merger dynamics that connects these two sets of observables. This includes if, and when, the post-merger remnant star collapses to a black hole, and what are the necessary conditions to power a short gamma-ray burst, and other observed electromagnetic counterparts. Observing the collapse of the post-merger neutron star would shed led on these questions, constraining models for the short gamma-ray burst engine and the hot neutron star equation of state. In this work, we explore the scope of using gravitational wave detectors to measure the timing of the collapse either indirectly, by establishing the shut-off of the post-merger gravitational emission, or---more challengingly---directly, by detecting the collapse signal. For the indirect approach, we consider a kilohertz high-frequency detector design that utilises a previously studied coupled arm cavity and signal recycling cavity resonance. This design would give a signal-to-noise ratio of 0.5\,-\,8.6 (depending on the variation of waveform parameters) for a collapse gravitational wave signal occurring at 10\,ms post-merger of a binary at 50\,Mpc and with total mass $2.7 M_\odot$. For the direct approach, we propose a narrow-band detector design, utilising the sensitivity around the frequency of the arm cavity free spectral range. The proposed detector achieves a signal-to-noise ratio of 0.3\,-\,1.9, independent of the collapse time. This detector is limited by both the fundamental classical and quantum noise with the arm cavity power chosen as 10\,MW.

A Magnetar Engine for Short GRBs and Kilonovae

Abstract We investigate the influence of magnetic fields on the evolution of binary neutron star (BNS) merger remnants via three-dimensional (3D) dynamical-spacetime general-relativistic magnetohydrodynamic (MHD) simulations. We evolve a post-merger remnant with an initial poloidal magnetic field, resolve the magnetoturbulence driven by shear flows, and include a microphysical finite-temperature equation of state. A neutrino leakage scheme that captures the overall energetics and lepton number exchange is also included. We find that turbulence induced by the magnetorotational instability in the hypermassive neutron star (HMNS) amplifies magnetic field to beyond magnetar strength (1015 G). The ultra-strong toroidal field is able to launch a relativistic jet from the HMNS. We also find a magnetized wind that ejects neutron-rich material with a rate of . The total ejecta mass in our simulation is 5 × 10−3 M ⊙. This makes the ejecta from the HMNS an important component in BNS mergers and a promising source of r-process elements that can power a kilonova. The jet from the HMNS reaches a terminal Lorentz factor of ∼5 in our highest-resolution simulation. The formation of this jet is aided by neutrino cooling preventing the accretion disk from protruding into the polar region. As neutrino pair-annihilation and radiative processes in the jet (which were not included in the simulations) will boost the Lorentz factor in the jet further, our simulations demonstrate that magnetars formed in BNS mergers are a viable engine for short gamma-ray bursts.

Deciphering the properties of the central engine in GRB collapsars

ABSTRACT The central engine in long gamma-ray bursts (GRBs) is thought to be a compact object produced by the core collapse of massive stars, but its exact nature (black hole or millisecond magnetar) is still debatable. Although the central engine of GRB collapsars is hidden to direct observation, its properties may be imprinted on the accompanying electromagnetic signals. We aim to decipher the generic properties of central engines that are consistent with prompt observations of long GRBs detected by the Burst Alert Telescope (BAT) on board the Neil Gehrels Swift Observatory. Adopting a generic model for the central engine, in which the engine power and activity time-scale are independent of each other, we perform Monte Carlo simulations of long GRBs produced by jets that successfully breakout from the star. Our simulations consider the dependence of the jet breakout time-scale on the engine luminosity and the effects of the detector’s flux threshold. The two-dimensional (2D) distribution of simulated detectable bursts in the gamma-ray luminosity versus gamma-ray duration plane is consistent with the observed one for a range of parameter values describing the central engine. The intrinsic 2D distribution of simulated collapsar GRBs peaks at lower gamma-ray luminosities and longer durations than the observed one, a prediction that can be tested in the future with more sensitive detectors. Black hole accretors, whose power and activity time are set by the large-scale magnetic flux through the progenitor star and stellar structure, respectively, are compatible with the properties of the central engine inferred by our model.