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

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.

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.

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.

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.

Maximally rotating supermassive stars at the onset of collapse: the perturbative effects of gas pressure, magnetic fields, dark matter, and dark energy

The discovery of quasars at increasingly large cosmological redshifts may favour 'direct collapse' as the most promising evolutionary route to the formation of supermassive black holes. In this scenario, supermassive black holes form when their progenitors - supermassive stars - become unstable to gravitational collapse. For uniformly rotating stars supported by pure radiation pressure and spinning at the mass-shedding limit, the critical configuration at the onset of collapse is characterized by universal values of the dimensionless spin and radius parameters J/M2 and R/M, independent of mass M. We consider perturbative effects of gas pressure, magnetic fields, dark matter, and dark energy on these parameters, and thereby determine the domain of validity of this universality. We obtain leading-order corrections for the critical parameters and establish their scaling with the relevant physical parameters. We compare two different approaches to approximate the effects of gas pressure, which plays the most important role, find identical results for the above dimensionless parameters, and also find good agreement with recent numerical results.

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.

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.