Properties of the Binary Black Hole Merger GW150914

Testing the No-Hair Theorem with GW150914

We analyze gravitational-wave data from the first LIGO detection of a binary black-hole merger (GW150914) in search of the ringdown of the remnant black hole. Using observations beginning at the peak of the signal, we find evidence of the fundamental quasinormal mode and at least one overtone, both associated with the dominant angular mode (ℓ=m=2), with 3.6σ confidence. A ringdown model including overtones allows us to measure the final mass and spin magnitude of the remnant exclusively from postinspiral data, obtaining an estimate in agreement with the values inferred from the full signal. The mass and spin values we measure from the ringdown agree with those obtained using solely the fundamental mode at a later time, but have smaller uncertainties. Agreement between the postinspiral measurements of mass and spin and those using the full waveform supports the hypothesis that the GW150914 merger produced a Kerr black hole, as predicted by general relativity, and provides a test of the no-hair theorem at the ∼10% level. An independent measurement of the frequency of the first overtone yields agreement with the no-hair hypothesis at the ∼20% level. As the detector sensitivity improves and the detected population of black-hole mergers grows, we can expect that using overtones will provide even stronger tests.

Lyman-α Forest Constraints on Primordial Black Holes as Dark Matter

The renewed interest in the possibility that primordial black holes (PBHs) may constitute a significant part of dark matter has provided motivation for revisiting old observational constraints, as well as developing new ones. We present new limits on the PBH abundance, from a comprehensive analysis of high-resolution high-redshift Lyman-α forest data. Poisson fluctuations in the PBH number density induce a small-scale power enhancement which departs from the standard cold dark matter prediction. Using a grid of hydrodynamic simulations exploring different values of astrophysical parameters, we obtain a marginalized upper limit on the PBH mass of f_{PBH}M_{PBH}∼60M_{⊙} at 2σ, when a Gaussian prior on the reionization redshift is imposed, preventing its posterior distribution from peaking on very high values, which are disfavored by the most recent estimates obtained both through cosmic microwave background and intergalactic medium observations. Such a bound weakens to f_{PBH}M_{PBH}∼170M_{⊙} when a conservative flat prior is instead assumed. Both limits significantly improve on previous constraints from the same physical observable. We also extend our predictions to nonmonochromatic PBH mass distributions, ruling out large regions of the parameter space for some of the most viable PBH extended mass functions.

Reduced-Order Modeling with Artificial Neurons for Gravitational-Wave Inference

Gravitational-wave data analysis is rapidly absorbing techniques from deep learning, with a focus on convolutional networks and related methods that treat noisy time series as images. We pursue an alternative approach, in which waveforms are first represented as weighted sums over reduced bases (reduced-order modeling); we then train artificial neural networks to map gravitational-wave source parameters into basis coefficients. Statistical inference proceeds directly in coefficient space, where it is theoretically straightforward and computationally efficient. The neural networks also provide analytic waveform derivatives, which are useful for gradient-based sampling schemes. We demonstrate fast and accurate coefficient interpolation for the case of a four-dimensional binary-inspiral waveform family and discuss promising applications of our framework in parameter estimation.

Approaching the Black Hole by Numerical Simulations

Black holes represent extreme conditions of physical laws. Predicted about a century ago, they are now accepted as astrophysical reality by most of the scientific community. Only recently has more direct evidence of their existence been found—the detection of gravitational waves from black hole mergers and of the shadow of a supermassive black hole in the center of a galaxy. Astrophysical black holes are typically embedded in an active environment which is affected by the strong gravity. When the environmental material emits radiation, this radiation may carry imprints of the black hole that is hosting the radiation source. In order to understand the physical processes that take place in the close neighborhood of astrophysical black holes, numerical methods and simulations play an essential role. This is simply because the dynamical evolution and the radiative interaction are far too complex in order to allow for an analytic solution of the physical equations. A huge progress has been made over the last decade(s) in the numerical code development, as well as in the computer power that is needed to run these codes. This review tries to summarize the basic questions and methods that are involved in the undertaking of investigating the astrophysics of black holes by numerical means. It is intended for a non-expert audience interested in an overview over this broad field. The review comes along without equations and thus without a detailed expert discussion of the underlying physical processes or numerical specifics. Instead, it intends to illustrate the richness of the field and to motivate further reading. The review puts some emphasis on magneto-hydrodynamic simulations but also touches radiation transfer and merger simulations, in particular pointing out differences in these approaches.

The Second-Order Correction to the Energy and Momentum in Plane Symmetric Gravitational Waves Like Spacetimes

This research provides second-order approximate Noether symmetries of geodetic Lagrangian of time-conformal plane symmetric spacetime. A time-conformal factor is of the form e ϵ f ( t ) which perturbs the plane symmetric static spacetime, where ϵ is small a positive parameter that produces perturbation in the spacetime. By considering the perturbation up to second-order in ϵ in plane symmetric spacetime, we find the second order approximate Noether symmetries for the corresponding Lagrangian. Using Noether theorem, the corresponding second order approximate conservation laws are investigated for plane symmetric gravitational waves like spacetimes. This technique tells about the energy content of the gravitational waves.

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.

Tracking Black Hole Kicks from Gravitational-Wave Observations

Coalescing binary black holes emit anisotropic gravitational radiation. This causes a net emission of linear momentum that produces a gradual acceleration of the source. As a result, the final remnant black hole acquires a characteristic velocity known as recoil velocity or gravitational kick. The symmetries of gravitational wave emission are reflected in the interactions of the gravitational wave modes emitted by the binary. In this Letter, we make use of the rich information encoded in the higher-order modes of the gravitational wave emission to infer the component of the kick along the line of sight (or radial kick). We do this by performing parameter inference on simulated signals given by numerical relativity waveforms for nonspinning binaries using numerical relativity templates of aligned-spin (nonprecessing) binary black holes. We find that for suitable sources, namely those with mass ratio q≥2 and total mass M∼100  M_{⊙}, and for modest radial kicks of 120  km/s, the 90% credible intervals of our posterior probability distributions can exclude a zero kick at a signal-to-noise ratio of 15, using a single Advanced LIGO detector working at its early sensitivity. The measurement of a nonzero radial kick component would provide the first observational signature of net transport of linear momentum by gravitational waves away from their source.

GW170817: Measurements of Neutron Star Radii and Equation of State

On 17 August 2017, the LIGO and Virgo observatories made the first direct detection of gravitational waves from the coalescence of a neutron star binary system. The detection of this gravitational-wave signal, GW170817, offers a novel opportunity to directly probe the properties of matter at the extreme conditions found in the interior of these stars. The initial, minimal-assumption analysis of the LIGO and Virgo data placed constraints on the tidal effects of the coalescing bodies, which were then translated to constraints on neutron star radii. Here, we expand upon previous analyses by working under the hypothesis that both bodies were neutron stars that are described by the same equation of state and have spins within the range observed in Galactic binary neutron stars. Our analysis employs two methods: the use of equation-of-state-insensitive relations between various macroscopic properties of the neutron stars and the use of an efficient parametrization of the defining function p(ρ) of the equation of state itself. From the LIGO and Virgo data alone and the first method, we measure the two neutron star radii as R_{1}=10.8_{-1.7}^{+2.0}  km for the heavier star and R_{2}=10.7_{-1.5}^{+2.1}  km for the lighter star at the 90% credible level. If we additionally require that the equation of state supports neutron stars with masses larger than 1.97  M_{⊙} as required from electromagnetic observations and employ the equation-of-state parametrization, we further constrain R_{1}=11.9_{-1.4}^{+1.4}  km and R_{2}=11.9_{-1.4}^{+1.4}  km at the 90% credible level. Finally, we obtain constraints on p(ρ) at supranuclear densities, with pressure at twice nuclear saturation density measured at 3.5_{-1.7}^{+2.7}×10^{34}  dyn cm^{-2} at the 90% level.

Erratum: GW170104: Observation of a 50-Solar-Mass Binary Black Hole Coalescence at Redshift 0.2 [Phys. Rev. Lett. 118, 221101 (2017)]

This corrects the article DOI: 10.1103/PhysRevLett.118.221101.

Erratum: Tests of General Relativity with GW150914 [Phys. Rev. Lett. 116, 221101 (2016)]

This corrects the article DOI: 10.1103/PhysRevLett.116.221101.

Fundamentals of numerical relativity for gravitational wave sources

Einstein's theory of general relativity affords an enormously successful description of gravity. The theory encodes the gravitational interaction in the metric, a tensor field on spacetime that satisfies partial differential equations known as the Einstein equations. This review introduces some of the fundamental concepts of numerical relativity-solving the Einstein equations on the computer-in simple terms. As a primary example, we consider the solution of the general relativistic two-body problem, which features prominently in the new field of gravitational wave astronomy.

Measuring the Hubble Constant with Neutron Star Black Hole Mergers

The detection of GW170817 and the identification of its host galaxy have allowed for the first standard-siren measurement of the Hubble constant, with an uncertainty of ∼14%. As more detections of binary neutron stars with redshift measurement are made, the uncertainty will shrink. The dominating factors will be the number of joint detections and the uncertainty on the luminosity distance of each event. Neutron star black hole mergers are also promising sources for advanced LIGO and Virgo. If the black hole spin induces precession of the orbital plane, the degeneracy between luminosity distance and the orbital inclination is broken, leading to a much better distance measurement. In addition, neutron star black hole sources are observable to larger distances, owing to their higher mass. Neutron star black holes could also emit electromagnetic radiation: depending on the black hole spin and on the mass ratio, the neutron star can be tidally disrupted, resulting in electromagnetic emission. We quantify the distance uncertainty for a wide range of black hole mass, spin, and orientations and find that the 1σ statistical uncertainty can be up to a factor of ∼10 better than for a nonspinning binary neutron star merger with the same signal-to-noise ratio. The better distance measurement, the larger gravitational-wave detectable volume, and the potentially bright electromagnetic emission imply that spinning black hole neutron star binaries can be the optimal standard-siren sources as long as their astrophysical rate is larger than O(10)  Gpc^{-3} yr^{-1}, a value allowed by current astrophysical constraints.

Prospects for observing and localizing gravitational-wave transients with Advanced LIGO, Advanced Virgo and KAGRA

We present possible observing scenarios for the Advanced LIGO, Advanced Virgo and KAGRA gravitational-wave detectors over the next decade, with the intention of providing information to the astronomy community to facilitate planning for multi-messenger astronomy with gravitational waves. We estimate the sensitivity of the network to transient gravitational-wave signals, and study the capability of the network to determine the sky location of the source. We report our findings for gravitational-wave transients, with particular focus on gravitational-wave signals from the inspiral of binary neutron star systems, which are the most promising targets for multi-messenger astronomy. The ability to localize the sources of the detected signals depends on the geographical distribution of the detectors and their relative sensitivity, and [Formula: see text] credible regions can be as large as thousands of square degrees when only two sensitive detectors are operational. Determining the sky position of a significant fraction of detected signals to areas of 5-[Formula: see text] requires at least three detectors of sensitivity within a factor of [Formula: see text] of each other and with a broad frequency bandwidth. When all detectors, including KAGRA and the third LIGO detector in India, reach design sensitivity, a significant fraction of gravitational-wave signals will be localized to a few square degrees by gravitational-wave observations alone.

First Higher-Multipole Model of Gravitational Waves from Spinning and Coalescing Black-Hole Binaries

Gravitational-wave observations of binary black holes currently rely on theoretical models that predict the dominant multipoles (ℓ=2,|m|=2) of the radiation during inspiral, merger, and ringdown. We introduce a simple method to include the subdominant multipoles to binary black hole gravitational waveforms, given a frequency-domain model for the dominant multipoles. The amplitude and phase of the original model are appropriately stretched and rescaled using post-Newtonian results (for the inspiral), perturbation theory (for the ringdown), and a smooth transition between the two. No additional tuning to numerical-relativity simulations is required. We apply a variant of this method to the nonprecessing PhenomD model. The result, PhenomHM, constitutes the first higher-multipole model of spinning and coalescing black-hole binaries, and currently includes the (ℓ,|m|)=(2,2),(3,3),(4,4),(2,1),(3,2),(4,3) radiative moments. Comparisons with numerical-relativity waveforms demonstrate that PhenomHM is more accurate than dominant-multipole-only models for all binary configurations, and typically improves the measurement of binary properties.

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.

Distinguishing spin-aligned and isotropic black hole populations with gravitational waves

The direct detection of gravitational waves from merging binary black holes opens up a window into the environments in which binary black holes form. One signature of such environments is the angular distribution of the black hole spins. Binary systems that formed through dynamical interactions between already-compact objects are expected to have isotropic spin orientations (that is, the spins of the black holes are randomly oriented with respect to the orbit of the binary system), whereas those that formed from pairs of stars born together are more likely to have spins that are preferentially aligned with the orbit. The best-measured combination of spin parameters for each of the four likely binary black hole detections GW150914, LVT151012, GW151226 and GW170104 is the 'effective' spin. Here we report that, if the magnitudes of the black hole spins are allowed to extend to high values, the effective spins for these systems indicate a 0.015 odds ratio against an aligned angular distribution compared to an isotropic one. When considering the effect of ten additional detections, this odds ratio decreases to 2.9 × 10^-7 against alignment. The existing preference for either an isotropic spin distribution or low spin magnitudes for the observed systems will be confirmed (or overturned) confidently in the near future.

Impact of Bayesian Priors on the Characterization of Binary Black Hole Coalescences

In a regime where data are only mildly informative, prior choices can play a significant role in Bayesian statistical inference, potentially affecting the inferred physics. We show this is indeed the case for some of the parameters inferred from current gravitational-wave measurements of binary black hole coalescences. We reanalyze the first detections performed by the twin LIGO interferometers using alternative (and astrophysically motivated) prior assumptions. We find different prior distributions can introduce deviations in the resulting posteriors that impact the physical interpretation of these systems. For instance, (i) limits on the 90% credible interval on the effective black hole spin χ_{eff} are subject to variations of ∼10% if a prior with black hole spins mostly aligned to the binary's angular momentum is considered instead of the standard choice of isotropic spin directions, and (ii) under priors motivated by the initial stellar mass function, we infer tighter constraints on the black hole masses, and in particular, we find no support for any of the inferred masses within the putative mass gap M≲5  M_{⊙}.

Palatini wormholes and energy conditions from the prism of general relativity

Wormholes are hypothetical shortcuts in spacetime that in general relativity unavoidably violate all of the pointwise energy conditions. In this paper, we consider several wormhole spacetimes that, as opposed to the standard designer procedure frequently employed in the literature, arise directly from gravitational actions including additional terms resulting from contractions of the Ricci tensor with the metric, and which are formulated assuming independence between metric and connection (Palatini approach). We reinterpret such wormhole solutions under the prism of General Relativity and study the matter sources that thread them. We discuss the size of violation of the energy conditions in different cases and how this is related to the same spacetimes when viewed from the modified gravity side.

Gravitational Waves from Binary Black Hole Mergers inside Stars

We present results from a controlled numerical experiment investigating the effect of stellar density gas on the coalescence of binary black holes (BBHs) and the resulting gravitational waves (GWs). This investigation is motivated by the proposed stellar core fragmentation scenario for BBH formation and the associated possibility of an electromagnetic counterpart to a BBH GW event. We employ full numerical relativity coupled with general-relativistic hydrodynamics and set up a 30+30  M_{⊙} BBH (motivated by GW150914) inside gas with realistic stellar densities. Our results show that at densities ρ≳10^{6}-10^{7}  g cm^{-3} dynamical friction between the BHs and gas changes the coalescence dynamics and the GW signal in an unmistakable way. We show that for GW150914, LIGO observations appear to rule out BBH coalescence inside stellar gas of ρ≳10^{7}  g cm^{-3}. Typical densities in the collapsing cores of massive stars are in excess of this density. This excludes the fragmentation scenario for the formation of GW150914.

GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral

On August 17, 2017 at 12∶41:04 UTC the Advanced LIGO and Advanced Virgo gravitational-wave detectors made their first observation of a binary neutron star inspiral. The signal, GW170817, was detected with a combined signal-to-noise ratio of 32.4 and a false-alarm-rate estimate of less than one per 8.0×10^{4}  years. We infer the component masses of the binary to be between 0.86 and 2.26  M_{⊙}, in agreement with masses of known neutron stars. Restricting the component spins to the range inferred in binary neutron stars, we find the component masses to be in the range 1.17-1.60  M_{⊙}, with the total mass of the system 2.74_{-0.01}^{+0.04}M_{⊙}. The source was localized within a sky region of 28  deg^{2} (90% probability) and had a luminosity distance of 40_{-14}^{+8}  Mpc, the closest and most precisely localized gravitational-wave signal yet. The association with the γ-ray burst GRB 170817A, detected by Fermi-GBM 1.7 s after the coalescence, corroborates the hypothesis of a neutron star merger and provides the first direct evidence of a link between these mergers and short γ-ray bursts. Subsequent identification of transient counterparts across the electromagnetic spectrum in the same location further supports the interpretation of this event as a neutron star merger. This unprecedented joint gravitational and electromagnetic observation provides insight into astrophysics, dense matter, gravitation, and cosmology.

GW170814: A Three-Detector Observation of Gravitational Waves from a Binary Black Hole Coalescence

On August 14, 2017 at 10∶30:43 UTC, the Advanced Virgo detector and the two Advanced LIGO detectors coherently observed a transient gravitational-wave signal produced by the coalescence of two stellar mass black holes, with a false-alarm rate of ≲1 in 27 000 years. The signal was observed with a three-detector network matched-filter signal-to-noise ratio of 18. The inferred masses of the initial black holes are 30.5_{-3.0}^{+5.7}M_{⊙} and 25.3_{-4.2}^{+2.8}M_{⊙} (at the 90% credible level). The luminosity distance of the source is 540_{-210}^{+130}  Mpc, corresponding to a redshift of z=0.11_{-0.04}^{+0.03}. A network of three detectors improves the sky localization of the source, reducing the area of the 90% credible region from 1160   deg^{2} using only the two LIGO detectors to 60  deg^{2} using all three detectors. For the first time, we can test the nature of gravitational-wave polarizations from the antenna response of the LIGO-Virgo network, thus enabling a new class of phenomenological tests of gravity.

Inferences about Supernova Physics from Gravitational-Wave Measurements: GW151226 Spin Misalignment as an Indicator of Strong Black-Hole Natal Kicks

The inferred parameters of the binary black hole GW151226 are consistent with nonzero spin for the most massive black hole, misaligned from the binary's orbital angular momentum. If the black holes formed through isolated binary evolution from an initially aligned binary star, this misalignment would then arise from a natal kick imparted to the first-born black hole at its birth during stellar collapse. We use simple kinematic arguments to constrain the characteristic magnitude of this kick, and find that a natal kick v_{k}≳50  km/s must be imparted to the black hole at birth to produce misalignments consistent with GW151226. Such large natal kicks exceed those adopted by default in most of the current supernova and binary evolution models.

GW170104: Observation of a 50-Solar-Mass Binary Black Hole Coalescence at Redshift 0.2

We describe the observation of GW170104, a gravitational-wave signal produced by the coalescence of a pair of stellar-mass black holes. The signal was measured on January 4, 2017 at 10∶11:58.6 UTC by the twin advanced detectors of the Laser Interferometer Gravitational-Wave Observatory during their second observing run, with a network signal-to-noise ratio of 13 and a false alarm rate less than 1 in 70 000 years. The inferred component black hole masses are 31.2_{-6.0}^{+8.4}M_{⊙} and 19.4_{-5.9}^{+5.3}M_{⊙} (at the 90% credible level). The black hole spins are best constrained through measurement of the effective inspiral spin parameter, a mass-weighted combination of the spin components perpendicular to the orbital plane, χ_{eff}=-0.12_{-0.30}^{+0.21}. This result implies that spin configurations with both component spins positively aligned with the orbital angular momentum are disfavored. The source luminosity distance is 880_{-390}^{+450}  Mpc corresponding to a redshift of z=0.18_{-0.07}^{+0.08}. We constrain the magnitude of modifications to the gravitational-wave dispersion relation and perform null tests of general relativity. Assuming that gravitons are dispersed in vacuum like massive particles, we bound the graviton mass to m_{g}≤7.7×10^{-23}  eV/c^{2}. In all cases, we find that GW170104 is consistent with general relativity.