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

Kilohertz quasiperiodic oscillations in short gamma-ray bursts

Short gamma-ray bursts (GRBs) are associated with binary neutron star mergers, which are multimessenger astronomical events that have been observed both in gravitational waves and in the multiband electromagnetic spectrum¹. Depending on the masses of the stars in the binary and on details of their largely unknown equation of state, a dynamically evolving and short-lived neutron star may be formed after the merger, existing for approximately 10-300 ms before collapsing to a black hole^2,3. Numerical relativity simulations across different groups consistently show broad power spectral features in the 1-5-kHz range in the post-merger gravitational-wave signal^4-14, which is inaccessible by current gravitational-wave detectors but could be seen by future third-generation ground-based detectors in the next decade^15-17. This implies the possibility of quasiperiodic modulation of the emitted gamma rays in a subset of events in which a neutron star is formed shortly before the final collapse to a black hole^18-21. Here we present two such signals identified in the short bursts GRB 910711 and GRB 931101B from archival Burst and Transient Source Experiment (BATSE) data, which are compatible with the predictions from numerical relativity.

Magnetohydrodynamic Simulations of Self-Consistent Rotating Neutron Stars with Mixed Poloidal and Toroidal Magnetic Fields

We perform the first magnetohydrodynamic simulations in full general relativity of self-consistent rotating neutron stars (NSs) with ultrastrong mixed poloidal and toroidal magnetic fields. The initial uniformly rotating NS models are computed assuming perfect conductivity, stationarity, and axisymmetry. Although the specific geometry of the mixed field configuration can delay or accelerate the development of various instabilities known from analytic perturbative studies, all our models finally succumb to them. Differential rotation is developed spontaneously in the cores of our magnetars which, after sufficient time, is converted back to uniform rotation. The rapidly rotating magnetars show a significant amount of ejecta, which can be responsible for transient kilonova signatures. However, no highly collimated, helical magnetic fields or incipient jets, which are necessary for γ-ray bursts, arise at the poles of these magnetars by the time our simulations are terminated.

PIC methods in astrophysics: simulations of relativistic jets and kinetic physics in astrophysical systems

The Particle-In-Cell (PIC) method has been developed by Oscar Buneman, Charles Birdsall, Roger W. Hockney, and John Dawson in the 1950s and, with the advances of computing power, has been further developed for several fields such as astrophysical, magnetospheric as well as solar plasmas and recently also for atmospheric and laser-plasma physics. Currently more than 15 semi-public PIC codes are available which we discuss in this review. Its applications have grown extensively with increasing computing power available on high performance computing facilities around the world. These systems allow the study of various topics of astrophysical plasmas, such as magnetic reconnection, pulsars and black hole magnetosphere, non-relativistic and relativistic shocks, relativistic jets, and laser-plasma physics. We review a plethora of astrophysical phenomena such as relativistic jets, instabilities, magnetic reconnection, pulsars, as well as PIC simulations of laser-plasma physics (until 2021) emphasizing the physics involved in the simulations. Finally, we give an outlook of the future simulations of jets associated to neutron stars, black holes and their merging and discuss the future of PIC simulations in the light of petascale and exascale computing.

Multimessenger Binary Mergers Containing Neutron Stars: Gravitational Waves, Jets, and γ-Ray Bursts

Neutron stars (NSs) are extraordinary not only because they are the densest form of matter in the visible Universe but also because they can generate magnetic fields ten orders of magnitude larger than those currently constructed on earth. The combination of extreme gravity with the enormous electromagnetic (EM) fields gives rise to spectacular phenomena like those observed on August 2017 with the merger of a binary neutron star system, an event that generated a gravitational wave (GW) signal, a short γ -ray burst (sGRB), and a kilonova. This event serves as the highlight so far of the era of multimessenger astronomy. In this review, we present the current state of our theoretical understanding of compact binary mergers containing NSs as gleaned from the latest general relativistic magnetohydrodynamic simulations. Such mergers can lead to events like the one on August 2017, GW170817, and its EM counterparts, GRB 170817 and AT 2017gfo. In addition to exploring the GW emission from binary black hole-neutron star and neutron star-neutron star mergers, we also focus on their counterpart EM signals. In particular, we are interested in identifying the conditions under which a relativistic jet can be launched following these mergers. Such a jet is an essential feature of most sGRB models and provides the main conduit of energy from the central object to the outer radiation regions. Jet properties, including their lifetimes and Poynting luminosities, the effects of the initial magnetic field geometries and spins of the coalescing NSs, as well as their governing equation of state, are discussed. Lastly, we present our current understanding of how the Blandford-Znajek mechanism arises from merger remnants as the trigger for launching jets, if, when and how a horizon is necessary for this mechanism, and the possibility that it can turn on in magnetized neutron ergostars, which contain ergoregions, but no horizons.

Gravitational waves from disks around spinning black holes: Simulations in full general relativity

We present fully general-relativistic numerical evolutions of self-gravitating tori around spinning black holes with dimensionless spin a/M = 0.7 parallel or antiparallel to the disk angular momentum. The initial disks are unstable to the hydrodynamic Papaloizou-Pringle instability which causes them to grow persistent orbiting matter clumps. The effect of black hole spin on the growth and saturation of the instability is assessed. We find that the instability behaves similarly to prior simulations with nonspinning black holes, with a shift in frequency due to spin-induced changes in disk orbital period. Copious gravitational waves are generated by these systems, and we analyze their detectability by current and future gravitational wave observatories for a large range of masses. We find that systems of 10 M ⊙-relevant for black hole-neutron star mergers-are detectable by Cosmic Explorer out to ~300 Mpc, while DECIGO (LISA) will be able to detect systems of 1000 M ⊙ (105 M ⊙)-relevant for disks forming in collapsing supermassive stars-out to cosmological redshift of z ~ 5 (z ~ 1). Computing the accretion rate of these systems we find that these systems may also be promising sources of coincident electromagnetic signals.

Black hole-neutron star coalescence: Effects of the neutron star spin on jet launching and dynamical ejecta mass

Black hole-neutron star (BHNS) mergers are thought to be sources of gravitational waves (GWs) with coincident electromagnetic (EM) counterparts. To further probe whether these systems are viable progenitors of short gamma-ray bursts (SGRBs) and kilonovas, and how one may use (the lack of) EM counterparts associated with LIGO/Virgo candidate BHNS GW events to sharpen parameter estimation, we study the impact of neutron star spin in BHNS mergers. Using dynamical spacetime magnetohydrodynamic simulations of BHNSs initially on a quasicircular orbit, we survey configurations that differ in the BH spin (a BH/M BH = 0 and 0.75), the NS spin (a NS/M NS = -0.17, 0, 0.23, and 0.33), and the binary mass ratio (q = M BH:M NS = 3:1 and 5:1). The general trend we find is that increasing the NS prograde spin increases both the rest mass of the accretion disk onto the remnant black hole, and the rest mass of dynamically ejected matter. By a time Δt ~ 3500-5500M ~ 88-138(M NS/1.4 M ⊙) ms after the peak gravitational-wave amplitude, a magnetically driven jet is launched only for q = 3:1 regardless of the initial NS spin. The lifetime of the jets [Δt ~ 0.5-0.8(M NS/1.4 M ⊙) s] and their outgoing Poynting luminosity [L Poyn ~ 1051.5±0.5 erg/s] are consistent with typical SGRBs' luminosities and expectations from the Blandford-Znajek mechanism. By the time we terminate our simulations, we do not observe either an outflow or a large-scale magnetic-field collimation for the other systems we consider. The mass range of dynamically ejected matter is 10-4.5-10-2(M NS/1.4 M ⊙) M ⊙, which can power kilonovas with peak bolometric luminosities L knova ~ 1040-1041.4 erg/s with rise times ≲6.5 h and potentially detectable by the LSST.

Neutron star merger remnants

Binary neutron star mergers observations are a unique way to constrain fundamental physics and astrophysics at the extreme. The interpretation of gravitational-wave events and their electromagnetic counterparts crucially relies on general-relativistic models of the merger remnants. Quantitative models can be obtained only by means of numerical relativity simulations in 3 + 1 dimensions including detailed input physics for the nuclear matter, electromagnetic and weak interactions. This review summarizes the current understanding of merger remnants focusing on some of the aspects that are relevant for multimessenger observations.

Binary Neutron Star Merger Simulations with a Calibrated Turbulence Model

Magnetohydrodynamic (MHD) turbulence in neutron star (NS) merger remnants can impact their evolution and multi-messenger signatures, complicating the interpretation of present and future observations. Due to the high Reynolds numbers and the large computational costs of numerical relativity simulations, resolving all the relevant scales of the turbulence will be impossible for the foreseeable future. Here, we adopt a method to include subgrid-scale turbulence in moderate resolution simulations by extending the large-eddy simulation (LES) method to general relativity (GR). We calibrate our subgrid turbulence model with results from very-high-resolution GRMHD simulations, and we use it to perform NS merger simulations and study the impact of turbulence. We find that turbulence has a quantitative, but not qualitative, impact on the evolution of NS merger remnants, on their gravitational wave signatures, and on the outflows generated in binary NS mergers. Our approach provides a viable path to quantify uncertainties due to turbulence in NS mergers.

Magnetohydrodynamic simulations of binary neutron star mergers in general relativity: Effects of magnetic field orientation on jet launching

Binary neutron star mergers can be sources of gravitational waves coincident with electromagnetic counterpart emission across the spectrum. To solidify their role as multimessenger sources, we present fully 3D, general relativistic, magnetohydrodynamic simulations of highly spinning binary neutrons stars initially on quasicircular orbits that merge and undergo delayed collapse to a black hole. The binaries consist of two identical stars modeled as Γ = 2 polytropes with spin χ NS = 0.36 aligned along the direction of the total orbital angular momentum L. Each star is initially threaded by a dynamical unimportant interior dipole magnetic field. The field is extended into the exterior where a nearly force-free magnetosphere resembles that of a pulsar. The magnetic dipole moment μ is either aligned or perpendicular to L and has the same initial magnitude for each orientation. For comparison, we also impose symmetry across the orbital plane in one case where μ in both stars is aligned along L. We find that the lifetime of the transient hypermassive neutron star remnant, the jet launching time, and the ejecta (which can give rise to a detectable kilonova) are very sensitive to the magnetic field orientation. By contrast, the physical properties of the black hole + disk remnant, such as the mass and spin of the black hole, the accretion rate, and the electromagnetic (Poynting) luminosity, are roughly independent of the initial magnetic field orientation. In addition, we find imposing symmetry across the orbital plane does not play a significant role in the final outcome of the mergers. Our results suggest that, as in the black hole-neutron star merger scenario, an incipient jet emerges only when the seed magnetic field has a sufficiently large-scale poloidal component aligned to the initial orbital angular momentum. The lifetime [Δt ≳ 140(M NS/1.625 M ⊙) ms] and Poynting luminosities [L EM ≃ 1052 erg/s] of the jet, when it forms, are consistent with typical short gamma-ray bursts, as well as with the Blandford-Znajek mechanism for launching jets.

Dynamically Stable Ergostars Exist: General Relativistic Models and Simulations

We construct the first dynamically stable ergostars (equilibrium neutron stars that contain an ergoregion) for a compressible, causal equation of state. We demonstrate their stability by evolving both strict and perturbed equilibrium configurations in full general relativity for over a hundred dynamical timescales (≳30 rotational periods) and observing their stationary behavior. This stability is in contrast to earlier models which prove radially unstable to collapse. Our solutions are highly differentially rotating hypermassive neutron stars with a corresponding spherical compaction of C=0.3. Such ergostars can provide new insights into the geometry of spacetimes around highly compact, rotating objects and on the equation of state at supranuclear densities. Ergostars may form as remnants of extreme binary neutron star mergers and possibly provide another mechanism for powering short gamma-ray bursts.

Numerical relativity of compact binaries in the 21st century

We review the dramatic progress in the simulations of compact objects and compact-object binaries that has taken place in the first two decades of the twenty-first century. This includes simulations of the inspirals and violent mergers of binaries containing black holes and neutron stars, as well as simulations of black-hole formation through failed supernovae and high-mass neutron star-neutron star mergers. Modeling such events requires numerical integration of the field equations of general relativity in three spatial dimensions, coupled, in the case of neutron-star containing binaries, with increasingly sophisticated treatment of fluids, electromagnetic fields, and neutrino radiation. However, it was not until 2005 that accurate long-term evolutions of binaries containing black holes were even possible (Pretorius 2005 Phys. Rev. Lett. 95 121101, Campanelli et al 2006 Phys. Rev. Lett. 96 111101, Baker et al 2006 Phys. Rev. Lett. 96 111102). Since then, there has been an explosion of new results and insights into the physics of strongly-gravitating system. Particular emphasis has been placed on understanding the gravitational wave and electromagnetic signatures from these extreme events. And with the recent dramatic discoveries of gravitational waves from merging black holes by the Laser Interferometric Gravitational Wave Observatory and Virgo, and the subsequent discovery of both electromagnetic and gravitational wave signals from a merging neutron star-neutron star binary, numerical relativity became an indispensable tool for the new field of multimessenger astronomy.

Simulating the magnetorotational collapse of supermassive stars: Incorporating gas pressure perturbations and different rotation profiles

Collapsing supermassive stars (SMSs) with masses M ≳ 104-6 M ⊙ have long been speculated to be the seeds that can grow and become supermassive black holes (SMBHs). We previously performed general relativistic magnetohydrodynamic (GRMHD) simulations of marginally stable Γ = 4/3 polytropes uniformly rotating at the mass-shedding limit and endowed initially with a dynamically unimportant dipole magnetic field to model the direct collapse of SMSs. These configurations are supported entirely by thermal radiation pressure and reliably model SMSs with M ≳ 106 M ⊙. We found that around 90% of the initial stellar mass forms a spinning black hole (BH) remnant surrounded by a massive, hot, magnetized torus, which eventually launches a magnetically-driven jet. SMSs could be therefore sources of ultra-long gamma-ray bursts (ULGRBs). Here we perform GRMHD simulations of Γ ≳ 4/3, polytropes to account for the perturbative role of gas pressure in SMSs with M ≲ 106 M ⊙. We also consider different initial stellar rotation profiles. The stars are initially seeded with a dynamically weak dipole magnetic field that is either confined to the stellar interior or extended from its interior into the stellar exterior. We calculate the gravitational wave burst signal for the different cases. We find that the mass of the black hole remnant is 90%-99% of the initial stellar mass, depending sharply on Γ - 4/3 as well as on the initial stellar rotation profile. After t ~ 250-550M ≈ 1 - 2 × 103(M/106 M ⊙) s following the appearance of the BH horizon, an incipient jet is launched and it lasts for ~104-105(M/106 M ⊙) s, consistent with the duration of long gamma-ray bursts. Our numerical results suggest that the Blandford-Znajek mechanism powers the incipient jet. They are also in rough agreement with our recently proposed universal model that estimates accretion rates and electromagnetic (Poynting) luminosities that characterize magnetized BH-disk remnant systems that launch a jet. This model helps explain why the outgoing electromagnetic luminosities computed for vastly different BH-disk formation scenarios all reside within a narrow range (~1052±1 erg s-1), roughly independent of M.

From hadrons to quarks in neutron stars: a review

In recent years our understanding of neutron stars has advanced remarkably, thanks to research converging from many directions. The importance of understanding neutron star behavior and structure has been underlined by the recent direct detection of gravitational radiation from merging neutron stars. The clean identification of several heavy neutron stars, of order two solar masses, challenges our current understanding of how dense matter can be sufficiently stiff to support such a mass against gravitational collapse. Programs underway to determine simultaneously the mass and radius of neutron stars will continue to constrain and inform theories of neutron star interiors. At the same time, an emerging understanding in quantum chromodynamics (QCD) of how nuclear matter can evolve into deconfined quark matter at high baryon densities is leading to advances in understanding the equation of state of the matter under the extreme conditions in neutron star interiors. We review here the equation of state of matter in neutron stars from the solid crust through the liquid nuclear matter interior to the quark regime at higher densities. We focus in detail on the question of how quark matter appears in neutron stars, and how it affects the equation of state. After discussing the crust and liquid nuclear matter in the core we briefly review aspects of microscopic quark physics relevant to neutron stars, and quark models of dense matter based on the Nambu-Jona-Lasinio framework, in which gluonic processes are replaced by effective quark interactions. We turn then to describing equations of state useful for interpretation of both electromagnetic and gravitational observations, reviewing the emerging picture of hadron-quark continuity in which hadronic matter turns relatively smoothly, with at most only a weak first order transition, into quark matter with increasing density. We review construction of unified equations of state that interpolate between the reasonably well understood nuclear matter regime at low densities and the quark matter regime at higher densities. The utility of such interpolations is driven by the present inability to calculate the dense matter equation of state in QCD from first principles. As we review, the parameters of effective quark models-which have direct relevance to the more general structure of the QCD phase diagram of dense and hot matter-are constrained by neutron star mass and radii measurements, in particular favoring large repulsive density-density and attractive diquark pairing interactions. We describe the structure of neutron stars constructed from the unified equations of states with crossover. Lastly we present the current equations of state-called 'QHC18' for quark-hadron crossover-in a parametrized form practical for neutron star modeling.

Disks around merging binary black holes: From GW150914 to supermassive black holes

We perform magnetohydrodynamic simulations in full general relativity of disk accretion onto nonspinning black hole binaries with mass ratio q = 29/36. We survey different disk models which differ in their scale height, total size and magnetic field to quantify the robustness of previous simulations on the initial disk model. Scaling our simulations to LIGO GW150914 we find that such systems could explain possible gravitational wave and electromagnetic counterparts such as the Fermi GBM hard x-ray signal reported 0.4 s after GW150915 ended. Scaling our simulations to supermassive binary black holes, we find that observable flow properties such as accretion rate periodicities, the emergence of jets throughout inspiral, merger and postmerger, disk temperatures, thermal frequencies, and the time delay between merger and the boost in jet outflows that we reported in earlier studies display only modest dependence on the initial disk model we consider here.

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.

What can we learn from GRBs?

We review our recent results on the classification of long and short gamma-ray bursts (GRBs) in different subclasses. We provide observational evidences for the binary nature of GRB progenitors. For long bursts the induced gravitational collapse (IGC) paradigm proposes as progenitor a tight binary system composed of a carbon-oxygen core (COcore) and a neutron star (NS) companion; the supernova (SN) explosion of the COcore triggers a hypercritical accretion process onto the companion NS. For short bursts a NS–NS merger is traditionally adopted as the progenitor. We also indicate additional sub-classes originating from different progenitors: (COcore)–black hole (BH), BH–NS, and white dwarf–NS binaries. We also show how the outcomes of the further evolution of some of these sub-classes may become the progenitor systems of other sub-classes.

Prompt Electromagnetic Transients from Binary Black Hole Mergers

Binary black hole (BBH) mergers provide a prime source for current and future interferometric GW observatories. Massive BBH mergers may often take place in plasma-rich environments, leading to the exciting possibility of a concurrent electromagnetic (EM) signal observable by traditional astronomical facilities. However, many critical questions about the generation of such counterparts remain unanswered. We explore mechanisms that may drive EM counterparts with magnetohydro-dynamic simulations treating a range of scenarios involving equal-mass black-hole binaries immersed in an initially homogeneous fluid with uniform, orbitally aligned magnetic fields. We find that the time development of Poynting luminosity, which may drive jet-like emissions, is relatively insensitive to aspects of the initial configuration. In particular, over a significant range of initial values, the central magnetic field strength is effectively regulated by the gas flow to yield a Poynting luminosity of 1045 - 1046 ρ -13 M 8 2 ergs-1, with BBH mass scaled to M 8 ≡ M/(108 M ⨀) and ambient density ρ -13 ≡ ρ/(10-13 g cm-3). We also calculate the direct plasma synchrotron emissions processed through geodesic ray-tracing. Despite lensing effects and dynamics, we find the observed synchrotron flux varies little leading up to merger.

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.

Binary neutron star mergers: a review of Einstein's richest laboratory

In a single process, the merger of binary neutron star systems combines extreme gravity, the copious emission of gravitational waves, complex microphysics and electromagnetic processes, which can lead to astrophysical signatures observable at the largest redshifts. We review here the recent progress in understanding what could be considered Einstein's richest laboratory, highlighting in particular the numerous significant advances of the last decade. Although special attention is paid to the status of models, techniques and results for fully general-relativistic dynamical simulations, a review is also offered on the initial data and advanced simulations with approximate treatments of gravity. Finally, we review the considerable amount of work carried out on the post-merger phase, including black-hole formation, torus accretion onto the merged compact object, the connection with gamma-ray burst engines, ejected material, and its nucleosynthesis.

Magnetorotational collapse of supermassive stars: Black hole formation, gravitational waves, and jets

We perform magnetohydrodynamic simulations in full general relativity of uniformly rotating stars that are marginally unstable to collapse. These simulations model the direct collapse of supermassive stars (SMSs) to seed black holes that can grow to become the supermassive black holes at the centers of quasars and active galactic nuclei. They also crudely model the collapse of massive Population III stars to black holes, which could power a fraction of distant, long gamma-ray bursts. The initial stellar models we adopt are Γ = 4/3 polytropes initially with a dynamically unimportant dipole magnetic field. We treat initial magnetic-field configurations either confined to the stellar interior or extending out from the stellar interior into the exterior. We find that the black hole formed following collapse has mass MBH ≃ 0.9M (where M is the mass of the initial star) and dimensionless spin parameter aBH/MBH ≃ 0.7. A massive, hot, magnetized torus surrounds the remnant black hole. At Δt ~ 400-550M ≈ 2000 - 2700(M/106M⊙)s following the gravitational wave peak amplitude, an incipient jet is launched. The disk lifetime is Δt ~ 105(M/106M⊙)s, and the outgoing Poynting luminosity is LEM ~ 1051-52 ergs/s. If > ˜ 1 % - 10 % of this power is converted into gamma rays, Swift and Fermi could potentially detect these events out to large redshifts z ~ 20. Thus, SMSs could be sources of ultra-long gamma-ray bursts (ULGRBs), and massive Population III stars could be the progenitors that power a fraction of the long GRBs observed at redshift z ~ 5-8. Gravitational waves are copiously emitted during the collapse and peak at ~15(106M⊙/M) mHz [~0.15(104 M⊙/M) Hz], i.e., in the LISA (DECIGO/BBO) band; optimally oriented SMSs could be detectable by LISA (DECIGO/BBO) at z < ˜ 3 ( z < ˜ 11 ) .Hence, 104M⊙ SMSs collapsing at z ~ 10 are promising multimessenger sources of coincident gravitational and electromagnetic waves.

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.