Gamma-ray bursts as collimated jets from neutron star/black hole mergers

The Rate of Short-Duration Gamma-Ray Bursts in the Local Universe

Following the faint gamma-ray burst, GRB 170817A, coincident with a gravitational wave-detected binary neutron star merger at d ∼ 40 Mpc, we consider the constraints on a local population of faint short duration GRBs (defined here broadly as T 90 < 4 s). We review proposed low-redshift short-GRBs and consider statistical limits on a d ≲ 200 Mpc population using Swift/Burst Alert Telescope (BAT), Fermi/Gamma-ray Burst Monitor (GBM), and Compton Gamma-Ray Observatory (CGRO) Burst and Transient Source Experiment (BATSE) GRBs. Swift/BAT short-GRBs give an upper limit for the all-sky rate of < 4 y − 1 at d < 200 Mpc, corresponding to < 5% of SGRBs. Cross-correlation of selected CGRO/BATSE and Fermi/GBM GRBs with d < 100 Mpc galaxy positions returns a weaker constraint of ≲ 12 y − 1 . A separate search for correlations due to SGR giant flares in nearby ( d < 11 Mpc) galaxies finds an upper limit of < 3 y − 1 . Our analysis suggests that GRB 170817A-like events are likely to be rare in existing SGRB catalogues. The best candidate for an analogue remains GRB 050906, where the Swift/BAT location was consistent with the galaxy IC 0327 at d ≈ 132 Mpc. If binary neutron star merger rates are at the high end of current estimates, then our results imply that at most a few percent will be accompanied by detectable gamma-ray flashes in the forthcoming LIGO/Virgo science runs.

GW170817, general relativistic magnetohydrodynamic simulations, and the neutron star maximum mass

Recent numerical simulations in general relativistic magnetohydrodynamics (GRMHD) provide useful constraints for the interpretation of the GW170817 discovery. Combining the observed data with these simulations leads to a bound on the maximum mass of a cold, spherical neutron star (the TOV limit): M max sph ≲ 2.74 / β , where β is the ratio of the maximum mass of a uniformly rotating neutron star (the supramassive limit) over the maximum mass of a nonrotating star. Causality arguments allow β to be as high as 1.27, while most realistic candidate equations of state predict β to be closer to 1.2, yielding M max sph in the range 2.16-2.28M⊙. A minimal set of assumptions based on these simulations distinguishes this analysis from previous ones, but leads a to similar estimate. There are caveats, however, and they are enumerated and discussed. The caveats can be removed by further simulations and analysis to firm up the basic argument.

General relativistic magnetohydrodynamics simulations of prompt-collapse neutron star mergers: The absence of jets

Inspiraling and merging binary neutron stars are not only important source of gravitational waves, but also promising candidates for coincident electromagnetic counterparts. These systems are thought to be progenitors of short gamma-ray bursts (sGRBs). We have shown previously that binary neutron star mergers that undergo delayed collapse to a black hole surrounded by a weighty magnetized accretion disk can drive magnetically powered jets. We now perform magnetohydrodynamic simulations in full general relativity of binary neutron stars mergers that undergo prompt collapse to explore the possibility of jet formation from black hole-light accretion disk remnants. We find that after t - tBH ~26(MNS/1.8 M⊙) ms (MNS is the ADM mass) following prompt black hole formation, there is no evidence of mass outflow or magnetic field collimation. The rapid formation of the black hole following merger prevents magnetic energy from approaching force-free values above the magnetic poles, which is required for the launching of a jet by the usual Blandford-Znajek mechanism. Detection of gravitational waves in coincidence with sGRBs may provide constraints on the nuclear equation of state (EOS): the fate of an NSNS merger-delayed or prompt collapse, and hence the appearance or nonappearance of an sGRB-depends on a critical value of the total mass of the binary, and this value is sensitive to the EOS.

Gamma-ray bursts and their use as cosmic probes

Since the launch of the highly successful and ongoing Swift mission, the field of gamma-ray bursts (GRBs) has undergone a revolution. The arcsecond GRB localizations available within just a few minutes of the GRB alert has signified the continual sampling of the GRB evolution through the prompt to afterglow phases revealing unexpected flaring and plateau phases, the first detection of a kilonova coincident with a short GRB, and the identification of samples of low-luminosity, ultra-long and highly dust-extinguished GRBs. The increased numbers of GRB afterglows, GRB-supernova detections, redshifts and host galaxy associations has greatly improved our understanding of what produces and powers these immense, cosmological explosions. Nevertheless, more high-quality data often also reveal greater complexity. In this review, I summarize some of the milestones made in GRB research during the Swift era, and how previous widely accepted theoretical models have had to adapt to accommodate the new wealth of observational data.

Black holes, disks, and jets following binary mergers and stellar collapse: The narrow range of electromagnetic luminosities and accretion rates

We have performed magnetohydrodynamic simulations in general relativity of binary neutron star and binary black hole-neutron star mergers, as well as the magnetorotational collapse of supermassive stars. In many cases the outcome is a spinnng black hole (BH) immersed in a magnetized disk, with a jet emanating from the poles of the BH. While their formation scenarios differ and their BH masses, as well as their disk masses, densities, and magnetic field strengths, vary by orders of magnitude, these features conspire to generate jet Poynting luminosities that all lie in the same, narrow range of ~1052±1 erg s-1. A similar result applies to their BH accretion rates upon jet launch, which is ~0.1-10 M⊙ s-1. We provide a simple model that explains these unanticipated findings. Interestingly, these luminosities reside in the same narrow range characterizing the observed luminosity distributions of over 400 short and long GRBs with distances inferred from spectroscopic redshifts or host galaxies. This result, together with the GRB lifetimes predicted by the model, supports the belief that a compact binary merger is the progenitor of an SGRB, while a massive, stellar magnetorotational collapse is the progenitor of an LGRB.

BINARY NEUTRON STAR MERGERS: A JET ENGINE FOR SHORT GAMMA-RAY BURSTS

We perform magnetohydrodynamic simulations in full general relativity (GRMHD) of quasi-circular, equal-mass, binary neutron stars that undergo merger. The initial stars are irrotational, n = 1 polytropes and are magnetized. We explore two types of magnetic-field geometries: one where each star is endowed with a dipole magnetic field extending from the interior into the exterior, as in a pulsar, and the other where the dipole field is initially confined to the interior. In both cases the adopted magnetic fields are initially dynamically unimportant. The merger outcome is a hypermassive neutron star that undergoes delayed collapse to a black hole (spin parameter a/MBH ~ 0.74) immersed in a magnetized accretion disk. About 4000M ~ 60(MNS/1.625M⊙) ms following merger, the region above the black hole poles becomes strongly magnetized, and a collimated, mildly relativistic outflow-an incipient jet-is launched. The lifetime of the accretion disk, which likely equals the lifetime of the jet, is Δ t ~ 0.1 (MNS/1.625M⊙) s. In contrast to black hole-neutron star mergers, we find that incipient jets are launched even when the initial magnetic field is confined to the interior of the stars.

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.

The Evolution of Compact Binary Star Systems

We review the formation and evolution of compact binary stars consisting of white dwarfs (WDs), neutron stars (NSs), and black holes (BHs). Binary NSs and BHs are thought to be the primary astrophysical sources of gravitational waves (GWs) within the frequency band of ground-based detectors, while compact binaries of WDs are important sources of GWs at lower frequencies to be covered by space interferometers (LISA). Major uncertainties in the current understanding of properties of NSs and BHs most relevant to the GW studies are discussed, including the treatment of the natal kicks which compact stellar remnants acquire during the core collapse of massive stars and the common envelope phase of binary evolution. We discuss the coalescence rates of binary NSs and BHs and prospects for their detections, the formation and evolution of binary WDs and their observational manifestations. Special attention is given to AM CVn-stars - compact binaries in which the Roche lobe is filled by another WD or a low-mass partially degenerate helium-star, as these stars are thought to be the best LISA verification binary GW sources.

X-ray Flares from Postmerger Millisecond Pulsars

Recent observations support the suggestion that short-duration gamma-ray bursts are produced by compact star mergers. The x-ray flares discovered in two short gamma-ray bursts last much longer than the previously proposed postmerger energy-release time scales. Here, we show that they can be produced by differentially rotating, millisecond pulsars after the mergers of binary neutron stars. The differential rotation leads to windup of interior poloidal magnetic fields and the resulting toroidal fields are strong enough to float up and break through the stellar surface. Magnetic reconnection-driven explosive events then occur, leading to multiple x-ray flares minutes after the original gamma-ray burst.

An origin in the local Universe for some short γ-ray bursts

Gamma-ray bursts (GRBs) divide into two classes: 'long', which typically have initial durations of T90 > 2 s, and 'short', with durations of T90 < 2 s (where T90 is the time to detect 90% of the observed fluence). Long bursts, which on average have softer gamma-ray spectra, are known to be associated with stellar core-collapse events-in some cases simultaneously producing powerful type Ic supernovae. In contrast, the origin of short bursts has remained mysterious until recently. A subsecond intense 'spike' of gamma-rays during a giant flare from the Galactic soft gamma-ray repeater, SGR 1806-20, reopened an old debate over whether some short GRBs could be similar events seen in galaxies out to approximately 70 Mpc (refs 6-10; redshift z approximately 0.016). Shortly after that, localizations of a few short GRBs (with optical afterglows detected in two cases) have shown an apparent association with a variety of host galaxies at moderate redshifts. Here we report a correlation between the locations of previously observed short bursts and the positions of galaxies in the local Universe, indicating that between 10 and 25 per cent of short GRBs originate at low redshifts (z < 0.025).

A short γ-ray burst apparently associated with an elliptical galaxy at redshift z = 0.225

Gamma-ray bursts (GRBs) come in two classes: long (> 2 s), soft-spectrum bursts and short, hard events. Most progress has been made on understanding the long GRBs, which are typically observed at high redshift (z approximately 1) and found in subluminous star-forming host galaxies. They are likely to be produced in core-collapse explosions of massive stars. In contrast, no short GRB had been accurately (< 10'') and rapidly (minutes) located. Here we report the detection of the X-ray afterglow from--and the localization of--the short burst GRB 050509B. Its position on the sky is near a luminous, non-star-forming elliptical galaxy at a redshift of 0.225, which is the location one would expect if the origin of this GRB is through the merger of neutron-star or black-hole binaries. The X-ray afterglow was weak and faded below the detection limit within a few hours; no optical afterglow was detected to stringent limits, explaining the past difficulty in localizing short GRBs.

Numerical Hydrodynamics in Special Relativity

This review is concerned with a discussion of numerical methods for the solution of the equations of special relativistic hydrodynamics (SRHD). Particular emphasis is put on a comprehensive review of the application of high-resolution shock-capturing methods in SRHD. Results of a set of demanding test bench simulations obtained with different numerical SRHD methods are compared. Three applications (astrophysical jets, gamma-ray bursts and heavy ion collisions) of relativistic flows are discussed. An evaluation of various SRHD methods is presented, and future developments in SRHD are analyzed involving extension to general relativistic hydrodynamics and relativistic magneto-hydrodynamics. The review further provides FORTRAN programs to compute the exact solution of a 1D relativistic Riemann problem with zero and nonzero tangential velocities, and to simulate 1D relativistic flows in Cartesian Eulerian coordinates using the exact SRHD Riemann solver and PPM reconstruction.

ELECTRONIC SUPPLEMENTARY MATERIAL

Supplementary material is available for this article at 10.12942/lrr-2003-7 and is accessible for authorized users.

The signature of supernova ejecta in the X-ray afterglow of the gamma-ray burst 011211

Now that gamma-ray bursts (GRBs) have been determined to lie at cosmological distances, their isotropic burst energies are estimated to be as high as 1054 erg (ref. 2), making them the most energetic phenomena in the Universe. The nature of the progenitors responsible for the bursts remains, however, elusive. The favoured models range from the merger of two neutron stars in a binary system to the collapse of a massive star. Spectroscopic studies of the afterglow emission could reveal details of the environment of the burst, by indicating the elements present, the speed of the outflow and an estimate of the temperature. Here we report an X-ray spectrum of the afterglow of GRB011211, which shows emission lines of magnesium, silicon, sulphur, argon, calcium and possibly nickel, arising in metal-enriched material with an outflow velocity of the order of one-tenth the speed of light. These observations strongly favour models where a supernova explosion from a massive stellar progenitor precedes the burst event and is responsible for the outflowing matter.

Relativistic Jets from Collapsars

Using a collapsar progenitor model of MacFadyen & Woosley, we have simulated the propagation of an axisymmetric jet through a collapsing rotating massive star with the GENESIS multidimensional relativistic hydrodynamic code. The jet forms as a consequence of an assumed (constant or variable) energy deposition in the range of 1050-1051 ergs s-1 within a 30 degrees cone around the rotation axis. The jet flow is strongly beamed (approximately less than a few degrees), spatially inhomogeneous, and time dependent. The jet reaches the surface of the stellar progenitor (R*=2.98x1010 cm) intact. At breakout, the maximum Lorentz factor of the jet flow is 33. After breakout, the jet accelerates into the circumstellar medium, whose density is assumed to decrease exponentially and then become constant, rhoext=10-5 g cm-3. Outside the star, the flow begins to expand laterally also (v approximately c), but the beam remains very well collimated. At a distance of 2.54 R*, where the simulation ends, the Lorentz factor has increased to 44.

Numerical Hydrodynamics in Special Relativity

This review is concerned with a discussion of numerical methods for the solution of the equations of special relativistic hydrodynamics (SRHD). Particular emphasis is put on a comprehensive review of the application of high-resolution shock-capturing methods in SRHD. Results obtained with different numerical SRHD methods are compared, and two astrophysical applications of SRHD flows are discussed. An evaluation of the various numerical methods is given and future developments are analyzed.

ELECTRONIC SUPPLEMENTARY MATERIAL

Supplementary material is available for this article at 10.12942/lrr-1999-3.