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.

General Relativistic Simulations of the Quasicircular Inspiral and Merger of Charged Black Holes: GW150914 and Fundamental Physics Implications

We perform general-relativistic simulations of charged black holes targeting GW150914. We show that the inspiral is most efficient for detecting black hole charge through gravitational waves and that GW150914 is compatible with having charge-to-mass ratio as high as 0.3. Our work applies to electric and magnetic charge and to theories with black holes endowed with U(1) (hidden or dark) charges. Using our results, we place an upper bound on the deviation from general relativity in the dynamical strong-filed regime of Moffat's modified gravity.

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.

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.

Rotating stars in relativity

Rotating relativistic stars have been studied extensively in recent years, both theoretically and observationally, because of the information they might yield about the equation of state of matter at extremely high densities and because they are considered to be promising sources of gravitational waves. The latest theoretical understanding of rotating stars in relativity is reviewed in this updated article. The sections on equilibrium properties and on nonaxisymmetric oscillations and instabilities in f-modes and r-modes have been updated. Several new sections have been added on equilibria in modified theories of gravity, approximate universal relationships, the one-arm spiral instability, on analytic solutions for the exterior spacetime, rotating stars in LMXBs, rotating strange stars, and on rotating stars in numerical relativity including both hydrodynamic and magnetohydrodynamic studies of these objects.

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.

Gravitational wave content and stability of uniformly, rotating, triaxial neutron stars in general relativity

Targets for ground-based gravitational wave interferometers include continuous, quasiperiodic sources of gravitational radiation, such as isolated, spinning neutron stars. In this work, we perform evolution simulations of uniformly rotating, triaxially deformed stars, the compressible analogs in general relativity of incompressible, Newtonian Jacobi ellipsoids. We investigate their stability and gravitational wave emission. We employ five models, both normal and supramassive, and track their evolution with different grid setups and resolutions, as well as with two different evolution codes. We find that all models are dynamically stable and produce a strain that is approximately one-tenth the average value of a merging binary system. We track their secular evolution and find that all our stars evolve toward axisymmetry, maintaining their uniform rotation, rotational kinetic energy, and angular momentum profiles while losing their triaxiality.

Black Hole Spectroscopy with Coherent Mode Stacking

The measurement of multiple ringdown modes in gravitational waves from binary black hole mergers will allow for testing the fundamental properties of black holes in general relativity and to constrain modified theories of gravity. To enhance the ability of Advanced LIGO/Virgo to perform such tasks, we propose a coherent mode stacking method to search for a chosen target mode within a collection of multiple merger events. We first rescale each signal so that the target mode in each of them has the same frequency and then sum the waveforms constructively. A crucial element to realize this coherent superposition is to make use of a priori information extracted from the inspiral-merger phase of each event. To illustrate the method, we perform a study with simulated events targeting the ℓ=m=3 ringdown mode of the remnant black holes. We show that this method can significantly boost the signal-to-noise ratio of the collective target mode compared to that of the single loudest event. Using current estimates of merger rates, we show that it is likely that advanced-era detectors can measure this collective ringdown mode with one year of coincident data gathered at design sensitivity.

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.

Binary black-hole mergers in magnetized disks: simulations in full general relativity

We present results from the first fully general relativistic, magnetohydrodynamic (MHD) simulations of an equal-mass black-hole binary (BHBH) in a magnetized, circumbinary accretion disk. We simulate both the pre- and postdecoupling phases of a BHBH-disk system and both "cooling" and "no-cooling" gas flows. Prior to decoupling, the competition between the binary tidal torques and the effective viscous torques due to MHD turbulence depletes the disk interior to the binary orbit. However, it also induces a two-stream accretion flow and mildly relativistic polar outflows from the BHs. Following decoupling, but before gas fills the low-density "hollow" surrounding the remnant, the accretion rate is reduced, while there is a prompt electromagnetic luminosity enhancement following merger due to shock heating and accretion onto the spinning BH remnant. This investigation, though preliminary, previews more detailed general relativistic, MHD simulations we plan to perform in anticipation of future, simultaneous detections of gravitational and electromagnetic radiation from a merging BHBH-disk system.

Binary black-hole mergers in magnetized disks: simulations in full general relativity

We present results from the first fully general relativistic, magnetohydrodynamic (MHD) simulations of an equal-mass black-hole binary (BHBH) in a magnetized, circumbinary accretion disk. We simulate both the pre- and postdecoupling phases of a BHBH-disk system and both "cooling" and "no-cooling" gas flows. Prior to decoupling, the competition between the binary tidal torques and the effective viscous torques due to MHD turbulence depletes the disk interior to the binary orbit. However, it also induces a two-stream accretion flow and mildly relativistic polar outflows from the BHs. Following decoupling, but before gas fills the low-density "hollow" surrounding the remnant, the accretion rate is reduced, while there is a prompt electromagnetic luminosity enhancement following merger due to shock heating and accretion onto the spinning BH remnant. This investigation, though preliminary, previews more detailed general relativistic, MHD simulations we plan to perform in anticipation of future, simultaneous detections of gravitational and electromagnetic radiation from a merging BHBH-disk system.