Model Of GRBs Research Articles

Stripping Model for Short GRBs: The Impact of Nuclear Data

We investigate the impact of forthcoming nuclear data on the predictions of the neutron star (NS) stripping model for short gamma-ray bursts. The main area to which we pay attention is the NS crust. We show that the uncertain properties of the NS equation of state can significantly influence the stripping time tstr, the main dynamical parameter of the model. Based on the known time delay (tstr≈1.7 s) between the peak of the gravitational wave signal GW170817 and the detection of gamma photons from GRB170817A, we obtain new restrictions on the nuclear matter parameters, in particular, the symmetry energy slope parameter: L<114.5MeV. In addition, we study the process of nucleosynthesis in the outer and inner crusts of a low-mass NS. We show that the nucleosynthesis is strongly influenced by both the forthcoming nuclear data and the equation of state of the NS matter.

On the nature of the ultrarelativistic prompt emission phase of GRB 190114C and GRB 180720B

A profound difference has occurred in the analysis of GRBs initially analyzed in the domain of gamma-ray astronomy with the Compton Gamma Ray Observatory and the BATSE Detector with the extension to X-ray astronomy and optical astronomy introduced by Beppo-SAX and the KEK and VLT optical observatory, followed by the AGILE, Fermi mission, Niels Gehrels SWIFT mission as well as HEHSS and MAGIC observations. An authentic ``Copernican revolution" has occurred in the transition from the traditional model of GRBs originating in a single star progenitor (collapsar) to GRBs with a binary progenitor composed of a CO core and a companion NS. We evidence from 24 SN observations related to GRBs that all of them, once analyzed with the general relativistic corrections, can be identified with I Bc Sn with a common value of the luminosity and common time of occurrence of the peak of the optical emission. This can be understood in terms of a precursor composed of a CO core and a binary NS companion. By contrast, the GRBs differ profoundly in their energetics which can be expressed in terms of 3 different classes of BDHNe: BDHN I, BDHN II, BDHN III, also originating from precursors composed of a CO core and a binary companion. It is pointed out how the analysis of these systems has profoundly modified the concept of the BH and the associated fundamental physics necessary for their description since the publication of ``Introducing the Black Hole" by Ruffini and Wheeler in 1971. At first, it is illustrated how the most radical change has occurred in the introduction of the effective BH charge, with an electric field only function of the BH dimensionless spin and a background magnetic field $B0$. This definitely indicates the abandonment of the Kerr--Newman solution in favour of the Kerr metric, for energy extraction processes in BH physics. We then enter in the detailed description of the ``seven Seals", the seven Episodes characterizing GRB 190114c and the associated SN 2019jrj. Based on the article by Aimuratov et al., in preparation.

Reverse Shock Emission from Short GRBs

We investigate the expected radio emission from the reverse shock of short GRBs, using the fitted afterglow parameters. In light of recent results suggesting that in some cases the radio afterglow is due to emission from the reverse shock, we examine the extent to which this component is detectable for short GRBs. In some GRBs, the standard synchrotron shock model predicts detectable radio emission from the reverse shock when none was seen. Many physical parameters play a role in these estimates, and our results highlight the need to explore the fundamental processes involved in GRB particle acceleration and emission more deeply. However, with a more rapid follow-up, we can test our standard model of GRBs, which predicts an early, radio bright reverse shock in many cases.

On the lack of a radio afterglow from some gamma-ray bursts – insight into their progenitors?

We investigate the intrinsic properties of a sample of bright (E_iso > 10^52 erg) gamma-ray bursts, comparing those with and without radio afterglows. We find that the sample of bursts with no radio afterglows has a significantly shorter mean intrinsic duration of the prompt gamma-ray radiation, and the distribution of this duration is significantly different from those bursts with a radio afterglow. Although the sample with no radio afterglow has on average lower isotropic energy, the lack of radio afterglow does not appear to be a result of simply energetics of the burst, but a reflection of a separate physical phenomenon likely related to the circumburst density profile. We also find a weak correlation between the isotropic gamma-ray energy and intrinsic duration in the sample with no radio afterglow, but not in the sample which have observed radio afterglows. We give possible explanations for why there may exist a sample of GRBs with no radio afterglow depending on whether the radio emission comes from the forward or reverse shock, and why these bursts appear to have intrinsically shorter prompt emission durations. We discuss how our results may have implications for progenitor models of GRBs.

RELATIVISTIC MHD SIMULATIONS OF COLLISION-INDUCED MAGNETIC DISSIPATION IN POYNTING-FLUX-DOMINATED JETS/OUTFLOWS

We perform 3D relativistic ideal MHD simulations to study the collisions between high-$\sigma$ (Poynting-flux-dominated) blobs which contain both poloidal and toroidal magnetic field components. This is meant to mimic the interactions inside a highly variable Poynting-flux-dominated jet. We discover a significant electromagnetic field (EMF) energy dissipation with an Alfv\'enic rate with the efficiency around 35\%. Detailed analyses show that this dissipation is mostly facilitated by the collision-induced magnetic reconnection. Additional resolution and parameter studies show a robust result that the relative EMF energy dissipation efficiency is nearly independent of the numerical resolution or most physical parameters in the relevant parameter range. The reconnection outflows in our simulation can potentially form the multi-orientation relativistic mini-jets as needed for several analytical models. We also find a linear relationship between the $\sigma$ values before and after the major EMF energy dissipation process. Our results give support to the proposed astrophysical models that invoke significant magnetic energy dissipation in Poynting-flux-dominated jets, such as the internal collision-induced magnetic reconnection and turbulence (ICMART) model for GRBs, and reconnection triggered mini-jets model for AGNs.

THE MILLISECOND MAGNETAR CENTRAL ENGINE IN SHORT GRBs

One favored progenitor model for short duration gamma-ray bursts (SGRBs) is the coalescence of two neutron stars (NS-NS). One possible outcome of such a merger would be a rapidly spinning, strongly magnetized neutron star (known as a millisecond magnetar). These magnetars may be "supra-massive", implying they would collapse to black holes after losing centrifugal support due to magnetic dipole spindown. By systematically analyzing the BAT-XRT light curves of all short GRBs detected by {\em swift}, we test how well the data are consistent with this central engine model of short GRBs. We find that the so-called "extended emission" observed with BAT in some short GRBs are fundamentally the same component as the "internal X-ray plateau" as observed in many short GRBs, which is defined as a plateau in the lightcurve followed by a very rapid drop. Based on how likely a short GRB hosts a magnetar, we characterize the entire {\em Swift} short GRB sample into three categories: the "internal plateau" sample, the "external plateau" sample, and the "no plateau" sample. The derived magnetar surface magnetic field $B_{\rm p}$ and the initial spin period $P_0$ fall into the reasonable range. No GRBs in the internal plateau sample have the total energy exceeding the maximum energy budget of a millisecond magnetar. Assuming that the rapid fall time at the end of the internal plateau is the collapse time of a supra-massive magnetar to a black hole, and applying the measured mass distribution of NS-NS systems in our Galaxy, we constrain the neutron star equation of state (EOS). The data suggest that the NS equation of state is close to the GM1 model, which has a maximum non-rotating NS mass $M_{\rm TOV} \sim 2.37 M_\odot$.

Search for gravitational waves associated with γ-ray bursts detected by the interplanetary network.

We present the results of a search for gravitational waves associated with 223 γ-ray bursts (GRBs) detected by the InterPlanetary Network (IPN) in 2005-2010 during LIGO's fifth and sixth science runs and Virgo's first, second, and third science runs. The IPN satellites provide accurate times of the bursts and sky localizations that vary significantly from degree scale to hundreds of square degrees. We search for both a well-modeled binary coalescence signal, the favored progenitor model for short GRBs, and for generic, unmodeled gravitational wave bursts. Both searches use the event time and sky localization to improve the gravitational wave search sensitivity as compared to corresponding all-time, all-sky searches. We find no evidence of a gravitational wave signal associated with any of the IPN GRBs in the sample, nor do we find evidence for a population of weak gravitational wave signals associated with the GRBs. For all IPN-detected GRBs, for which a sufficient duration of quality gravitational wave data are available, we place lower bounds on the distance to the source in accordance with an optimistic assumption of gravitational wave emission energy of 10(-2)M⊙c(2) at 150 Hz, and find a median of 13 Mpc. For the 27 short-hard GRBs we place 90% confidence exclusion distances to two source models: a binary neutron star coalescence, with a median distance of 12 Mpc, or the coalescence of a neutron star and black hole, with a median distance of 22 Mpc. Finally, we combine this search with previously published results to provide a population statement for GRB searches in first-generation LIGO and Virgo gravitational wave detectors and a resulting examination of prospects for the advanced gravitational wave detectors.

CONSTRAINTS ON THE SYNCHROTRON EMISSION MECHANISM IN GAMMA-RAY BURSTS

We reexamine the general synchrotron model for GRBs' prompt emission and determine the regime in the parameter phase space in which it is viable. We characterize a typical GRB pulse in terms of its peak energy, peak flux and duration and use the latest Fermi observations to constrain the high energy part of the spectrum. We solve for the intrinsic parameters at the emission region and find the possible parameter phase space for synchrotron emission. Our approach is general and it does not depend on a specific energy dissipation mechanism. Reasonable synchrotron solutions are found with energy ratios of $10^{-4}<\epsilon_B/\epsilon_e<10$, bulk Lorentz factor values of $300<\Gamma<3000$, typical electrons' Lorentz factor values of $3\times 10^3<\gamma_e<10^5$ and emission radii of the order $10^{15}cm<R<10^{17}$cm. Most remarkable among those are the rather large values of the emission radius and the electron's Lorentz factor. We find that soft (with peak energy less than 100KeV) but luminous (isotropic luminosity of $~1.5 \times 10^{53}$) pulses are inefficient. This may explain the lack of strong soft bursts. In cases when most of the energy is carried out by the kinetic energy of the flow, such as in the internal shocks, the synchrotron solution requires that only a small fraction of the electrons are accelerated to relativistic velocities by the shocks. We show that future observations of very high energy photons from GRBs by CTA, could possibly determine all parameters of the synchrotron model or rule it out altogether.

THE ELECTROMAGNETIC MODEL OF SHORT GRBs, THE NATURE OF PROMPT TAILS, SUPERNOVA-LESS LONG GRBs, AND HIGHLY EFFICIENT EPISODIC ACCRETION

Many short GRBs show prompt tails lasting up to hundreds of seconds that can be energetically dominant over the initial sub-second spike. In this paper we develop an electromagnetic model of short GRBs that explains the two stages of the energy release, the prompt spike and the prompt tail. The key ingredient of the model is the recent discovery that an isolated black hole can keep its open magnetic flux for times much longer than the collapse time and, thus, can spin-down electromagnetically, driving the relativistic wind. First, the merger is preceded by an electromagnetic precursor wind. If a fraction of the wind power is converted into pulsar-like coherent radio emission, this may produce an observable radio burst of few milliseconds. At the active stage of the merger, two neutron stars produces a black hole surrounded by an accretion torus in which the amplified magnetic field extracts the rotational energy of the black hole and drives an axially-collimated electromagnetic wind. For observers nearly aligned with the orbital normal this is seen as a classical short GRB. After the accretion of the torus, the isolated black hole keeps the open magnetic flux and drives the equatorially (not axially) collimated outflow, which is seen by an observer at intermediate polar angles as a prompt tail. The tail carries more energy than the prompt spike, but its emission is de-boosted for observers along the orbital normal. Observers in the equatorial plane miss the prompt spike and interpret the prompt tail as an energetic supernova-less long GRB. We also demonstrate that episodic accretion onto the BH of magnetized clouds that carry non-zero magnetic flux can be highly efficient in extracting the spin energy of the BH, producing the outflows with the power exceeding the average accretion power.

Acceleration of Magnetized Collapsar Jets After Breakout

In the collapsar model of long GRBs the jet is formed at the center of the progenitor star, propagates in its interior, and produces the observed gamma rays much after its breakout from the star. The loss of pressure support during breakout induces a strong rarefaction wave that propagates inside the jet and causes its bulk acceleration. This mechanism has been already studied using axisymmetric magnetohydrodynamic (MHD) simulations assuming a prescribed shape for the surface between the jet and its environment, as well as using simple rarefaction waves in planar geometry. Trying to improve over these works, we solve the steady-state, axisymmetric, relativistic MHD equations using the method of characteristics. In this way the jet boundary is found self-consistently and the rarefaction wave is studied in the axisymmetric geometry. In this poster we present our first results and a comparison with previous works.

Two-component jets from 3-dimensional magnetohydrodynamic jet simulations of disk winds at sub-parsec scales

We explore the effect of large scale magnetic field on the formation of two-component jets in magnetohydrodynamic disk winds simulations. Our simulations show a one-component and two-component jets develop depending on the magnetic field distribution along the surface of the accretion disk. Magnetic field configurations with the least steep gradient along the disk lead to a well defined two-component jet with the self-similar (Blandford-Payne) configuration separating the two regimes. Our results have direct implications to jets models of AGN and GRBs if indeed two-component jets emanate directly from the accretion disk. Our findings imply that a three-component jets may exist in AGN jets if one takes into account a Blandford-Znajek component in the innermost, relativistic, regions.

Magnetic jet model for GRBs and the delayed arrival of &amp;gt;100 MeV photons

ABSTRACT Photons of energy larger than 100 MeV from long GRBs arrive a few seconds after &amp;lt;10 MeV photons do. We show that this delay is a natural consequence of a magnetic-dominated relativistic jet. The much slower acceleration of a magnetic jet with radius (compared with a hot baryonic outflow) results in high-energy γ-ray photons to be converted to electron–positron pairs out to a larger radius, whereas lower energy γ-rays of energy less than ∼10 MeV can escape when the jet crosses the Thomson photosphere. The resulting delay for the arrival of high-energy photons is found to be similar to the value observed by the Fermi satellite for a number of GRBs. A prediction of this model is that the delay should increase with photon energy (E) as ∼E0.17 for E &amp;gt; 100 MeV. The delay depends almost linearly on burst redshift, and on the distance from the central compact object where the jet is launched (R0). Therefore, the delay in arrival of &amp;gt;102 MeV photons can be used to estimate burst redshift if the magnetic jet model for γ-ray generation is correct and R0 is roughly the same for long GRBs.

Testing Two-Component Jet Models of GRBs with Orphan Afterglows

Abstract In the Swift era, two-component jet models were introduced to explain the complex temporal profiles and the diversity of early afterglows. In this paper, we concentrate on the two-component jet model: the first component is the conventional afterglow, and second is emission due to late internal dissipation, such as late-prompt emission. We suggest herein that the two-component jet model can be probed by the existence of two optical peaks for orphan GRB afterglows. Each peak is caused by its respective jet as its relativistic beaming cone widens to encompass the off-axis line of sight. Typically, the first peak appears at 10$^{4}$–10$^{5}\ $s, and the second at 10$^{5}\ $–10$^{6}\ $s. Furthermore, we expect to observe a single, bright X-ray peak at the same time as the first optical peak. Because orphan afterglows do not have prompt emission, it is necessary to monitor the entire sky every 10$^{4}$s in the X-ray regime. We can test the model with orphan afterglows through the X-ray all-sky survey collaboration, and by using ground-based optical telescopes.

Investigating γ- and X-ray properties of GRBs using multivariate statistics

We propose to study the interrelation between the γ- (Fluence, 1sec Peakflux, duration) and X-ray (early X flux, 24 h X flux, X decay index, X spectral index, X HI column density) properties using the canonical correlation method. Computing the canonical correlations and variables we show that there is a significant interrelation between the γ- and X-ray data. Using the canonical variables from the analysis, we computed their correlations (canonical loadings) with the original ones. The canonical loadings revealed that the γ-ray fluence and the early X-ray flux give the strongest contribution to the correlation in contrast to the X-ray decay index and spectral index. An interesting new result appears to be the strong contribution of the HI column density to the correlation. Accepting the collapsar model of long GRBs this effect may be interpreted as an indication for the ejection of an HI envelope by the progenitor in the course of producing the GRB.

Gamma-ray bursts as a result of an interaction of a supernova shock with the star-companion in a binary system

Interaction of a fast shock wave generated during a supernova explosion with a magnetized star-companion of the supernova precursor produces a current sheet. We consider an evolution of this current sheet and show that a singularity (shock) is formed in finite time within the ideal magnetohydrodynamics framework. Charged particles (electrons) are accelerated in the vicinity of the singularity, and their distribution function has a plateau up to the energies of the order of 104mc2. These fast particles radiate in the γ-range in the strong magnetic field of the current sheet (B ≃ 106 G). Radiation is concentrated within a narrow angle around the current sheet, Δθ ≃ 3 × 10−4, and its spectrum has the maximum at several hundreds of keV. Presented calculations confirm the model of cosmological GRBs proposed earlier by Istomin & Komberg.

OPTICAL VARIABILITY OF THE QUASAR 3C 279: THE SIGNATURE OF A DECELERATING JET?

A recent optical monitoring campaign on the prominent quasar 3C279 revealed a period of a remarkably clean exponential decay of BVRI fluxes with time, with a time constant of 12.8 d, over about 14 days. This is clearly too long to be associated with radiative cooling. Here we propose that this may be the signature of deceleration of the synchrotron emitting jet component. We develop a model analogous to the relativistic blast wave model for GRBs, including radiative energy losses and radiation drag, to simulate the deceleration of a relativistically moving plasmoid in the moderately dense AGN environment. Synchrotron, SSC and external Compton emission are evaluated self-consistently. We show that the observed optical light curve decay can be successfully reproduced with this model. The decelerating plasmoid model predicts a delayed X-ray flare, about 2 - 3 weeks after the onset of the quasi-exponential light curve decay in the optical. A robust prediction of this model, which can be tested with Fermi and simultaneous optical monitoring, is that the peak in the gamma-ray light curve at ~ 100 MeV is expected to be delayed by a few days with respect to the onset of the optical decay, while the VHE gamma-rays are expected to track the optical light curve closely with a delay of at most a few hours.

A new general relativistic magnetohydrodynamics code for dynamical spacetimes

We present a new numerical code which solves the general relativistic magneto-hydrodynamics (GRMHD) equations coupled to the Einstein equations for the evolution of a dynamical spacetime within the conformally-flat approximation. This code has been developed with the main objective of studying astrophysical scenarios in which both, high magnetic fields and strong gravitational fields appear, such as the magneto-rotational collapse of stellar cores, the collapsar model of GRBs, and the evolution of neutron stars. The code is based on an existing and thoroughly tested purely hydrodynamics code and on its extension to accommodate weakly magnetized fluids (passive magnetic field approximation). The numerical code we present here is based on high-resolution shock-capturing schemes to solve the GRMHD equations together with the flux constraint transport method to ensure the solenoidal condition of the magnetic field. Since the astrophysical applications envisaged do not deviate much from spherical symmetry, the conformal flatness condition approximation is used for the formulation of the Einstein equations. In addition, the code can handle several equations of state, from simple analytical expressions to microphysical tabulated ones. In this paper we present stringent tests of our new GRMHD numerical code, which show its ability to handle all aspects appearing in the astrophysical scenarios for which the code is intended, namely relativistic shocks, highly magnetized fluids, and equilibrium configurations of magnetized neutron stars. As an application, magneto-rotational core collapse simulations of a realistic progenitor are presented, comparing the results with our previous finding in the passive magnetic field approximation.

Long Gamma‐Ray Bursts without Visible Supernovae: A Case Study of Redshift Estimators and Alleged Novel Objects

It has been argued that the observational limits on a supernova (SN) associated with GRB060614 convincingly exclude a SN akin to SN1998bw as its originator, and provide evidence for a new class of long-duration GRBs. We discuss this issue in the contexts of indirect 'redshift estimators' and of the fireball and cannonball models of GRBs. The latter explains the unusual properties of GRB060614: at its debated but favoured low redshift (0.125) they are predicted, as opposed to exceptional, if the associated core-collapse SN is of a recently discovered, very faint type. We take the occasion to discuss the 'association' between GRBs and SNe.

Simulations of Magnetically Driven Supernova and Hypernova Explosions in the Context of Rapid Rotation

We present here the first 2D rotating, multi-group, radiation magnetohydrodynamics (RMHD) simulations of supernova core collapse, bounce, and explosion. In the context of rapid rotation, we focus on the dynamical effects of magnetic stresses and the creation and propagation of MHD jets. We find that a quasi-steady state can be quickly established after bounce, during which a well-collimated MHD jet is maintained by continuous pumping of power from the differentially rotating core. If the initial spin period of the progenitor core is $\sles$ 2 seconds, the free energy reservoir in the secularly evolving protoneutron star is adequate to power a supernova explosion, and may be enough for a hypernova. The jets are well collimated by the infalling material and magnetic hoop stresses, and maintain a small opening angle. We see evidence of sausage instabilities in the emerging jet stream. Neutrino heating is sub-dominant in the rapidly rotating models we explore, but can contribute 10$-$25% to the final explosion energy. Our simulations suggest that even in the case of modest or slow rotation, a supernova explosion might be followed by a secondary, weak MHD jet explosion, which, because of its weakness may to date have gone unnoticed in supernova debris. Furthermore, we suggest that the generation of a non-relativistic MHD precursor jet during the early protoneutron star/supernova phase is implicit in both the collapsar and "millisecond magnetar" models of GRBs. The multi-D, multi-group, rapidly rotating RMHD simulations we describe here are a start along the path towards more realistic simulations of the possible role of magnetic fields in some of Nature's most dramatic events.

Optical and X‐Ray Observations of GRB 060526: A Complex Afterglow Consistent with an Achromatic Jet Break

We obtained 98 R-band and 18 B, r', i' images of the optical afterglow of GRB 060526 (z = 3.21) with the MDM 1.3 m, 2.4 m, and the PROMPT telescopes at CTIO over the five nights following the burst trigger. Combining these data with other optical observations reported in GCN and the Swift XRT observations, we compare the optical and X-ray afterglow light curves of GRB 060526. Both the optical and X-ray afterglow light curves show rich features, such as flares and breaks. The densely sampled optical observations provide very good coverage at T > 10(exp 4) s. We observed a break at 2.4 x 10(exp 5) sin the optical afterglow light curve. Compared with the X-ray afterglow light curve, the break is consistent with an achromatic break supporting the beaming models of GRBs. However, the prebreak and postbreak temporal decay slopes are difficult to explain in simple afterglow models. We estimated a jet angle of theta(sub j) approx. 7deg and a prompt emission size of R(sub prompt) approx. 2 x 10(exp 14) cm. In addition, we detected several optical flares with amplitudes of (Delta)m approx. 0.2,0.6, and 0.2 mag. The X-ray afterglows detected by Swift have shown complicated decay patterns. Recently, many well-sampled optical afterglows also show decays with flares and multiple breaks. GRB 060526 provides an additional case of such a complex, well-observed optical afterglow. The accumulated well-sampled afterglows indicate that most of the optical afterglows are complex.