Engine Of Gamma-ray Bursts Research Articles

Effects of Fallback Accretion on Protomagnetar Outflows in Gamma-Ray Bursts and Superluminous Supernovae

Abstract Rapidly spinning, strongly magnetized protoneutron stars (“millisecond protomagnetars”) are candidate central engines of long-duration gamma-ray bursts (GRBs), superluminous supernovae (SLSNe), and binary neutron star mergers. Magnetar birth may be accompanied by the fallback of stellar debris, lasting for seconds or longer following the explosion. Accretion alters the magnetar evolution by (1) providing an additional source of rotational energy (or a potential sink, if the propeller mechanism operates), (2) enhancing the spin-down luminosity above the dipole rate by compressing the magnetosphere and expanding the polar cap region of open magnetic field lines, and (3) supplying an additional accretion-powered neutrino luminosity that sustains the wind baryon loading, even after the magnetar’s internal neutrino luminosity has subsided. The more complex evolution of the jet power and magnetization of an accreting magnetar more readily accounts for the high 56Ni yields of GRB SNe and the irregular time evolution of some GRB light curves (e.g., bursts with precursors followed by a long quiescent interval before the main emission episode). Additional baryon loading from accretion-powered neutrino irradiation of the polar cap lengthens the time frame over which the jet magnetization is in the requisite range σ ≲ 103 for efficient gamma-ray emission, thereby accommodating GRBs with ultralong durations. Though accretion does not significantly raise the maximum energy budget from the limit of ≲ few × 1052 erg for an isolated magnetar, it greatly expands the range of magnetic field strengths and birth spin periods capable of powering GRB jets, reducing the differences between the magnetar properties normally invoked to explain GRBs versus SLSNe.

The X-Ray Light Curve in GRB 170714A: Evidence for a Quark Star?

Abstract Two plateaus and a following bump in the X-ray light curve of GRB 170714A have been detected by the Swift/X-ray Telescope, which could be very significant for the central engine of gamma-ray bursts (GRBs), implying that the origin of this burst might be different from those of other ultra-long GRBs. We propose that merging two neutron stars into a hyper-massive quark star (QS) and then collapsing into a black hole (BH), with a delay time around 104 s, could be responsible for these X-ray components. The hyper-massive QS is initially in a fluid state, being turbulent and differentially rotating, but would solidify and release its latent heat, injecting it into the GRB fireball (lasting about 103 s during the liquid–solid phase transition). A magnetic field as high as ∼1015 G can be created by dynamo action of the newborn liquid QS, and a magnetar-like central engine (after solidification) supplies significant energy for the second plateau. More energy could be released during a fall-back accretion after the post-merger QS collapses to a BH, and the X-ray bump forms. This post-merger QS model could be tested by future observations, with either advanced gravitational wave detectors (e.g., advanced LIGO and VIRGO) or X-ray/optical telescopes.

External Shock in a Multi-bursting Gamma-Ray Burst: Energy Injection Phase Induced by the Later Launched Ejecta

Abstract Central engines of gamma-ray bursts (GRBs) may be intermittent and launch several episodes of ejecta separated by a long quiescent interval. In this scenario, an external shock is formed due to the propagation of the first launched ejecta into the circum-burst medium and the later launched ejecta may interact with the external shock at a later period. Owing to the internal dissipation, the later launched ejecta may be observed at a later time (t jet). In this paper, we study the relation of t b and t jet, where t b is the collision time of the later launched ejecta with the formed external shock. It is found that the relation of t b and t jet depends on the bulk Lorentz factor (Γjet) of the later launched ejecta and the density (ρ) of the circum-burst medium. If the value of Γjet or ρ is low, the t b would be significantly larger than t jet. However, the t b ∼ t jet can be found if the value of Γjet or ρ is significantly large. Our results can explain the large lag of the optical emission relative to the γ-ray/X-ray emission in GRBs, e.g., GRB 111209A. For GRBs with a precursor, our results suggest that the energy injection into the external shock and thus more than one external-reverse shock may appear in the main prompt emission phase. According to our model, we estimate the Lorentz factor of the second launched ejecta in GRB 160625B.

Comparison of Gravitational Waves from Central Engines of Gamma-Ray Bursts: Neutrino-dominated Accretion Flows, Blandford–Znajek Mechanisms, and Millisecond Magnetars

Abstract Neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH) is one of the plausible candidates for the central engines of gamma-ray bursts (GRBs). Because the time-variant and anisotropic emission of neutrinos from NDAFs leads to GRB variability, NDAFs can be regarded as the sources of strong gravitational waves (GWs). We calculate the dependences of the GW strains on both the BH spin and the accretion rate. We demonstrate that, for typical GRBs with either a single pulse or multiple pulses, the GWs from NDAFs might be detected at a distance of ∼100 kpc/∼1 Mpc by the advanced LIGO/Einstein Telescope with a typical frequency of ∼ 10 – 100 Hz . Besides NDAFs, the other two competitive candidates for the GRB central engine are the Blandford–Znajek (BZ) mechanism and millisecond magnetars. We explore the GW signals from these two as well, and compare the corresponding results with NDAFs’. We find that for a certain GRB, the possible detected distance from NDAFs is about two orders of magnitude higher than that from the BZ mechanism, but at least two orders of magnitude lower than that from magnetars. The typical GW frequency for the BZ mechanism is the same as that of NDAFs, ∼ 10 – 100 Hz , while the typical frequency for magnetars is ∼2000 Hz. Therefore, the GWs released by the central engines of adjacent GRBs might be used to determine whether there is an NDAF, a BZ jet, or a magnetar in the GRB center.

Hyperaccreting Black Hole as Gamma-Ray Burst Central Engine. II. Temporal Evolution of the Central Engine Parameters during the Prompt and Afterglow Phases

Abstract A hyperaccreting stellar-mass black hole (BH) has been proposed as the candidate central engine of gamma-ray bursts (GRBs). The rich observations of GRBs by Fermi and Swift make it possible to constrain the central engine model by comparing the model predictions against data. This paper is dedicated to studying the temporal evolution of the central engine parameters for both the prompt emission and afterglow phases. We consider two jet-launching mechanisms, i.e., ν ν ¯ annihilations and the Blandford–Znajek (BZ) process, and obtain analytical solutions to these two models. We then investigate the BH central engine parameters, such as the jet power, the dimensionless entropy η, and the central engine parameter μ 0 = η (1 + σ 0) (where σ 0 is the initial magnetization of the engine) at the base of the jet. The BH may be spun up by accretion or spun down by the BZ process, leaving imprints in the GRB light curves. Usually, a BZ jet is more powerful and is likely responsible for the late-time central engine activities. However, an initially non-spinning BH central engine may first launch a thermal “fireball” via neutrino annihilations, and then launch a Poynting-flux-dominated jet via the BZ process. Multiple flares, giant bumps, and plateaus in GRB afterglows can be produced as the result of late-time accretion onto the BH.

Constraints on millisecond magnetars as the engines of prompt emission in gamma-ray bursts

We examine millisecond magnetars as central engines of Gamma Ray Bursts' (GRB) prompt emission. Using the proto-magnetar wind model of Metzger et al. 2011, we estimate the temporal evolution of the magnetization and power injection at the base of the GRB jet and apply these to different prompt emission models to make predictions for the GRB energetics, spectra and lightcurves. We investigate both shock and magnetic reconnection models for the particle acceleration, as well as the effects of energy dissipation across optically thick and thin regions of the jet. The magnetization at the base of the jet, $\sigma_0$, is the main parameter driving the GRB evolution in the magnetar model and the emission is typically released for $100\lesssim \sigma_0 \lesssim 3000$. Given the rapid increase in $\sigma_0$ as the proto-magnetar cools and its neutrino-driven mass loss subsides, the GRB duration is typically limited to $\lesssim 100$ s. This low baryon loading at late times challenges magnetar models for ultra-long GRBs, though black hole models likely run into similar difficulties without substantial entrainment from the jet walls. The maximum radiated gamma-ray energy is $\lesssim 5 \times 10^{51}$erg, significantly less than the magnetar's total initial rotational energy and in strong tension with the high end of the observed GRB energy distribution. However, the gradual magnetic dissipation model (Beniamini & Giannios 2017) applied to a magnetar central engine, naturally explains several key observables of typical GRBs, including energetics, durations, stable peak energies, spectral slopes and a hard to soft evolution during the burst.

A common central engine for long gamma-ray bursts and Type Ib/c supernovae

Long-duration, spectrally-soft Gamma-Ray Bursts (GRBs) are associated with Type Ic Core Collapse (CC) Supernovae (SNe), and thus arise from the death of massive stars. In the collapsar model, the jet launched by the central engine must bore its way out of the progenitor star before it can produce a GRB. Most of these jets do not break out, and are instead "choked" inside the star, as the central-engine activity time, $t_{\rm e}$, is not long enough. Modelling the long-soft GRB duration distribution assuming a power-law distribution for their central-engine activity times, $\propto t_{\rm e}^{-\alpha}$ for $t_{\rm e}>t_{\rm b}$, we find a steep distribution ($\alpha\sim4$) and a typical GRB jet breakout time of $t_{\rm b}\sim 60\text{ s}$ in the star's frame. The latter suggests the presence of a low-density, extended envelope surrounding the progenitor star, similar to that previously inferred for low-luminosity GRBs. Extrapolating the range of validity of this power law below what is directly observable, to $t_{\rm e}<t_{\rm b}$, by only a factor of $\sim$4-5 produces enough events to account for all Type Ib/c SNe. Such extrapolation is necessary to avoid fine-tuning the distribution of central engine activity times with the breakout time, which are presumably unrelated. We speculate that central engines launching relativistic jets may operate in all Type Ib/c SNe. In this case, the existence of a common central engine would imply that (i) the jet may significantly contribute to the energy of the SN; (ii) various observational signatures, like the asphericity of the explosion, could be directly related to jet's interaction with the star.

Self-gravity in Magnetized Neutrino-dominated Accretion Disks

Abstract In the present work we study self-gravity effects on the vertical structure of a magnetized neutrino-dominated accretion disk as a central engine for gamma-ray bursts (GRBs). Some of the disk physical timescales that are supposed to play a pivotal role in the late-time evolutions of the disk, such as viscous, cooling, and diffusion timescales, have been studied. We are interested in investigating the possibility of the occurrence of X-ray flares, observed in late-time GRBs’ extended emission through the “magnetic barrier” and “fragmentation” processes in our model. The results lead us to interpret self-gravity as an amplifier for Blandford–Payne luminosity (BP power) and the generated magnetic field, but a suppressor for neutrino luminosity and magnetic barrier processes via highlighting the fragmentation mechanism in the outer disk, especially for the higher mass accretion rates.

The status of general relativistic simulations of compact binary mergers as engines of short gamma-ray bursts

Black hole - neutron star (BHNS) and neutron star - neutron star (NSNS) binaries are perhaps the most popular progenitors for short-hard gamma ray bursts. After about two decades of numerical relativity simulations of binary compact objects we are beginning to understand the necessary ingredients for jets to emerge from these systems following merger. We report on the latest development of this field summarizing the results from state-of-the-art numerical relativity (magnetohydrodynamic) simulations of compact binary mergers as progenitors of short-hard gamma-ray bursts.

The Signature of the Blandford-Znajek Mechanism in GRB Light Curves

In 1977, Blandford and Znajek showed how the spin energy of a rotating black hole may be extracted electromagnetically through a magnetic field that threads the black hole horizon. A characteristic feature of this mechanism is that, under certain fairly general conditions, the energy loss rate decays exponentially. We looked precisely for such behavior in the X-ray light curves of Long and Ultra Long duration Gamma-Ray Bursts (GRBs) observed with the XRT instrument on board the Swift satellite, and found that almost 30% of XRT light curves show an exponential decay before they reach the afterglow plateau. A similar behavior (Fast Rise Exponential Decay-FRED) was observed in γ -rays with the BATSE instrument aboard the CGRO satellite. We consider both of these findings as the signature of the Blandford-Znajek mechanism in action in the central engine of GRBs.

What Can We Learn about GRB from the Variability Timescale Related Correlations?

Abstract Recently, two empirical correlations related to the minimum variability timescale (MTS) of the light curves are discovered in gamma-ray bursts (GRBs). One is the anti-correlation between MTS and Lorentz factor Γ, and the other is the anti-correlation between the MTS and gamma-ray luminosity L γ . Both of the two correlations might be used to explore the activity of the central engine of GRBs. In this paper, we try to understand these empirical correlations by combining two popular black hole central engine models (namely, the Blandford &amp; Znajek mechanism (BZ) and the neutrino-dominated accretion flow (NDAF)). By taking the MTS as the timescale of viscous instability of the NDAF, we find that these correlations favor the scenario in which the jet is driven by the BZ mechanism.

General relativistic simulations of compact binary mergers as engines for short gamma-ray bursts

Black hole—neutron star (BHNS) and neutron star—neutron star (NSNS) binaries are among the favored candidates for the progenitors of the black hole—disk systems that may be the engines powering short-hard gamma ray bursts. After almost two decades of simulations of binary NSNSs and BHNSs in full general relativity we are now beginning to understand the ingredients that may be necessary for these systems to launch incipient jets. Here, we review our current understanding, and summarize the surprises and lessons learned from state-of-the-art (magnetohydrodynamic) simulations in full general relativity of BHNS and NSNS mergers as jet engines for short-hard gamma-ray bursts. We also propose a new approach to probing the nuclear equation of state by virtue of multimessenger observations.

Evidence for Magnetar Formation in Broad-lined Type Ic Supernovae1998bw and 2002ap

Abstract Broad-lined type Ic supernovae (SNe Ic-BL) are peculiar stellar explosions that are distinct from ordinary SNe. Some SNe Ic-BL are associated with long-duration ( ) gamma-ray bursts (GRBs). Black holes and magnetars are two types of compact objects that are hypothesized to be central engines of GRBs. In spite of decades of investigations, no direct evidence for the formation of black holes or magnetars has yet been found for GRBs. Here we report the finding that the early peak ( ) and late-time ( ) slow decay displayed in the light curves of SNe 1998bw (associated with GRB 980425) and 2002ap (not GRB-associated) can be attributed to magnetar spin-down with an initial rotation period , while the intermediate-time ( ) exponential decline is caused by the radioactive decay of 56Ni. The connection between the early peak and late-time slow decline in the light curves is unexpected in alternative models. We thus suggest that GRB 980425 and SN 2002ap were powered by magnetars.

NUMERICAL AND ANALYTICAL SOLUTIONS OF NEUTRINO-DOMINATED ACCRETION FLOWS WITH A NON-ZERO TORQUE BOUNDARY CONDITION AND ITS APPLICATIONS IN GAMMA-RAY BURSTS

ABSTRACT A stellar-mass black hole (BH) surrounded by a neutrino-dominated accretion flow (NDAF) has been discussed in a number of works as the central engine of gamma-ray bursts (GRBs). It is widely believed that NDAF cannot liberate enough energy for bright GRBs. However, these works have been based on the assumption of a “no torque” boundary condition, which is invalid when the disk is magnetized. In this paper, we present both numerical and analytical solutions for NDAFs with non-zero boundary stresses and reexamine their properties. We find that an NDAF with such a boundary torque can be powerful enough to account for those bright short GRBs, energetic long GRBs, and ultra-long GRBs. The disk becomes viscously unstable, which makes it possible to interpret the variability of GRB prompt emission and the steep decay phase in the early X-ray afterglow. Finally, we study the gravitational waves radiated from a processing BH-NDAF. We find that the effects of the boundary torque on the strength of the gravitational waves can be ignored.

Revisiting the Correlations of Peak Luminosity with Spectral Lag and Peak Energy of the Observed Gamma-ray Bursts

An analysis of light curves and spectra of observed gamma-ray bursts in gamma-ray ranges is frequently demanded because the prompt emission contains immediate details regarding the central engine of gamma-ray bursts (GRBs). We have revisited the relationship between the collimation-corrected peak luminosity and the spectral lag, investigating the lag-luminosity relationships in great detail by focusing on spectral lags resulting from all possible combinations of channels. Firstly, we compiled the opening angle data and demonstrated that the distribution of opening angles of 205 long GRBs is represented by a double Gaussian function having maxima at ~ 0.1 and ~ 0.3 radians. We confirmed that the peak luminosity and the spectral lag are anti-correlated, both in the observer frame and in the source frame. We found that, in agreement with our previous conclusion, the correlation coefficient improves significantly in the source frame. It should be noted that spectral lags involving channel 2 (25-50 keV) yield high correlation coefficients, where Swift/Burst Alert Telescope (BAT) has four energy channels (channel 1: 15-25 keV, channel 2: 25-50 keV, channel 3: 50-100 keV, channel 4: 100-200 keV). We also found that peak luminosity is positively correlated with peak energy.

General relativistic magnetohydrodynamic simulations of binary neutron star mergers with the APR4 equation of state

We present new results of fully general relativistic magnetohydrodynamic simulations of binary neutron star (BNS) mergers performed with the Whisky code. All the models use a piecewise polytropic approximation of the APR4 equation of state for cold matter, together with a ‘hybrid’ part to incorporate thermal effects during the evolution. We consider both equal and unequal-mass models, with total masses such that either a supramassive NS or a black hole is formed after merger. Each model is evolved with and without a magnetic field initially confined to the stellar interior. We present the different gravitational wave (GW) signals as well as a detailed description of the matter dynamics (magnetic field evolution, ejected mass, post-merger remnant/disk properties). Our simulations provide new insights into BNS mergers, the associated GW emission and the possible connection with the engine of short gamma-ray bursts (both in the ‘standard’ and in the ‘time-reversal’ scenarios) and other electromagnetic counterparts.

A METHOD TO CONSTRAIN MASS AND SPIN OF GRB BLACK HOLES WITHIN THE NDAF MODEL

ABSTRACT Black holes (BHs) hide themselves behind various astronomical phenomena and their properties, i.e., mass and spin, are usually difficult to constrain. One leading candidate for the central engine model of gamma-ray bursts (GRBs) invokes a stellar mass BH and a neutrino-dominated accretion flow (NDAF), with the relativistic jet launched due to neutrino-anti-neutrino annihilations. Such a model gives rise to a matter-dominated fireball, and is suitable to interpret GRBs with a dominant thermal component with a photospheric origin. We propose a method to constrain BH mass and spin within the framework of this model and apply the method to the thermally dominant GRB 101219B, whose initial jet launching radius, r 0, is constrained from the data. Using our numerical model of NDAF jets, we estimate the following constraints on the central BH: mass M BH ∼ 5–9 M ⊙, spin parameter a * ≳ 0.6, and disk mass 3 M ⊙ ≲ M disk ≲ 4 M ⊙. Our results also suggest that the NDAF model is a competitive candidate for the central engine of GRBs with a strong thermal component.

Testing black hole neutrino-dominated accretion discs for long-duration gamma-ray bursts

Long-duration gamma-ray bursts (LGRBs) are generally considered to originate from the massive collapsars. It is believed that the central engine of gamma-ray bursts (GRBs) is a neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH). The neutrino annihilation above the NDAF is a feasible mechanism to power GRB. In this work, we analyse the distributions of the isotropic gamma-ray radiated energy and jet kinetic energy of 48 LGRBs. According to the NDAF and fireball models, we estimate the mean accreted masses of LGRBs in our sample to investigate whether the NDAFs can power LGRBs with the reasonable BH parameters and conversion efficiency of neutrino annihilation. The results indicate that most of the values of the accreted masses are less than $5~M_\odot$ for the extreme Kerr BHs and high conversion efficiency. It suggests that the NDAFs may be suitable for most of LGRBs except for some extremely high energy sources.

MULTI-WAVELENGTH OBSERVATIONS OF GRB 111228A AND IMPLICATIONS FOR THE FIREBALL AND ITS ENVIRONMENT

ABSTRACT Observations of very early multi-wavelength afterglows are critical to reveal the properties of the radiating fireball and its environment as well as the central engine of gamma-ray bursts (GRBs). We report our optical observations of GRB 111228A from 95 s to about 50 hr after the burst trigger and investigate its properties of the prompt gamma-rays and the ambient medium using our data and the data from the Swift and Fermi missions. Our joint optical and X-ray spectral fits to the afterglow data show that the ambient medium features a low dust-to-gas ratio. Incorporating the energy injection effect, our best fit to the afterglow light curves with the standard afterglow model via the Markov Chain Monte Carlo technique shows that &#x003F5; e = ( 6.9 &PlusMinus; 0.3 ) &times; 10 &minus; 2 ?> , &#x003F5; B = ( 7.73 &PlusMinus; 0.62 ) &times; 10 &minus; 6 , E K = ( 6.32 &PlusMinus; 0.86 ) &times; 10 53 erg ?> , n = 0.100 &PlusMinus; 0.014 ?> cm−3. The low medium density likely implies that the afterglow jet may be in a halo or in a hot ISM. A chromatic shallow decay segment observed in the optical and X-ray bands is well explained with the long-lasting energy injection from the central engine, which would be a magnetar with a period of about 1.92 ms inferred from the data. The E p of its time-integrated prompt gamma-ray spectrum is ∼26 KeV. Using the initial Lorentz factor ( &Gamma; 0 = 476 &minus; 237 &plus; 225 ?> ) derived from our afterglow model fit, it is found that GRB 111228A satisfies the L iso &minus; E p , z &minus; &Gamma; 0 ?> relation and bridges the typical GRBs and low luminosity GRBs in this relation.

POLARIZATION EVOLUTION OF EARLY OPTICAL AFTERGLOWS OF GAMMA-RAY BURSTS

ABSTRACT The central engine and jet composition of gamma-ray bursts (GRBs) remain mysterious. Here we suggest that observations on the polarization evolution of early optical afterglows may shed light on these questions. We first study the dynamics of a reverse shock and a forward shock that are generated during the interaction of a relativistic jet and its ambient medium. The jet is likely magnetized with a globally large-scale magnetic field from the central engine. The existence of the reverse shock requires that the magnetization degree of the jet should not be high (σ ≤ 1), so that the jet is mainly composed of baryons and leptons. We then calculate the light curves and polarization evolution of early optical afterglows and find that when the polarization position angle changes by 90° during the early afterglow, the polarization degree is zero for a toroidal magnetic field but is very likely to be nonzero for an aligned magnetic field. This result would be expected to provide a probe for the central engine of GRBs because an aligned field configuration could originate from a magnetar central engine and a toroidal field configuration could be produced from a black hole via the Blandford–Znajek mechanism. Finally, for such two kinds of magnetic field configurations, we fit the observed data of the early optical afterglow of GRB 120308A equally well.