A Detailed Analysis of GW190521 with Phenomenological Waveform Models

In this paper we present an extensive analysis of the GW190521 gravitational wave event with the current (fourth) generation of phenomenological waveform models for binary black hole coalescences. GW190521 stands out from other events since only a few wave cycles are observable. This leads to a number of challenges, one being that such short signals are prone to not resolving approximate waveform degeneracies, which may result in multimodal posterior distributions. The family of waveform models we use includes a new fast time-domain model (IMRPhenomTPHM), which allows us to extensively test different priors and robustness with respect to variations in the waveform model, including the content of spherical harmonic modes. We clarify some issues raised in a recent paper, Nitz & Capano, associated with possible support for a high-mass-ratio source, but confirm their finding of a multimodal posterior distribution, albeit with important differences in the statistical significance of the peaks. In particular, we find that the support for both masses being outside the pair instability supernova mass gap, and the support for an intermediate-mass-ratio binary are drastically reduced with respect to what Nitz & Capano found. We also provide updated probabilities for associating GW190521 to the potential electromagnetic counterpart from the Zwicky Transient Facility (ZTF) Graham et al.

Post-Newtonian waveforms from spinning scattering amplitudes

We derive the classical gravitational radiation from an aligned spin binary black hole on closed orbits, using a dictionary built from the 5-point QFT scattering amplitude of two massive particles exchanging and emitting a graviton. We show explicitly the agreement of the transverse-traceless components of the radiative linear metric perturbations — and the corresponding gravitational wave energy flux — at future null infinity, derived from the scattering amplitude and those derived utilizing an effective worldline action in conjunction with multipolar post-Minkowskian matching. At the tree-level, this result holds at leading orders in the black holes’ velocities and up to quadratic order in their spins. At sub-leading order in black holes’ velocities, we demonstrate a matching of the radiation field for quasi-circular orbits in the no-spin limit. At the level of the radiation field, and to leading order in the velocities, there exists a one-to-one correspondence between the binary black hole mass and current quadrupole moments, and the scalar and linear-in-spin scattering amplitudes, respectively. Therefore, we show explicitly that waveforms, needed to detect gravitational waves from inspiraling binary black holes, can be derived consistently, to the orders considered, from the classical limit of quantum scattering amplitudes.

Can We Prescribe the Physical Parameters of Multiple Black Holes?

The parabolic-hyperbolic form of the constraints and superposed Kerr-Schild black holes have already been used to provide a radically new initialization of binary black hole configurations. The method generalizes straightforwardly to multiple black hole systems. This paper is to verify that each of the global Arnowitt-Deser-Misner quantities of the constructed multiple black hole initial data can always be prescribed, as desired, in advance of solving the constraints. These global charges are shown to be uniquely determined by the physical parameters of the involved individual Kerr-Schild black holes.

The Eccentric and Accelerating Stellar Binary Black Hole Mergers in Galactic Nuclei: Observing in Ground and Space Gravitational-wave Observatories

We study the stellar binary black holes (BBHs) inspiraling/merging in galactic nuclei based on our numerical method GNC. We find that 3%–40% of all newborn BBHs will finally merge due to various dynamical effects. In a five-year mission, up to 104, 105, and ∼100 of BBHs inspiraling/merging in galactic nuclei can be detected with signal-to-noise ration >8 in Advanced LIGO (aLIGO), Einstein/DECIGO, and TianQin/LISA/TaiJi, respectively. Roughly tens are detectable in both LISA/TaiJi/TianQin and aLIGO. These BBHs have two unique characteristics. (1) Significant eccentricities: 1%–3%, 2%–7%, or 30%–90% of them have e i > 0.1 when they enter into aLIGO, Einstein, or space observatories, respectively. Such high eccentricities provide a possible explanation for that of GW190521. Most highly eccentric BBHs are not detectable in LISA/Tianqin/TaiJi before entering into aLIGO/Einstein, as their strain becomes significant only at f GW ≳ 0.1 Hz. DECIGO becomes an ideal observatory to detect those events, as it can fully cover the rising phase. (2) Up to 2% of BBHs can inspiral/merge at distances ≲103 r SW from the massive black hole, with significant accelerations, such that the Doppler phase drift of ∼10–105 of them can be detected with signal-to-noise ratio >8 in space observatories. The energy density of the gravitational-wave backgrounds (GWBs) contributed by these BBHs deviates from the power-law slope of 2/3 at f GW ≲ 1 mHz. The high eccentricity, significant accelerations, and the different profile of the GWB of these sources make them distinguishable, and thus interesting for future gravitational-wave detections and tests of relativities.

Characterizing the Observation Bias in Gravitational-wave Detections and Finding Structured Population Properties

The observed distributions of the source properties from gravitational-wave (GW) detections are biased due to the selection effects and detection criteria in the detections, analogous to the Malmquist bias. In this work, this observation bias is investigated through its fundamental statistical and physical origins. An efficient semi-analytical formulation for its estimation is derived, which is as accurate as the standard method of numerical simulations, with only a millionth of the computational cost. Then, the estimated bias is used for unmodeled inferences on the binary black hole population. These inferences show additional structures, specifically two peaks in the joint mass distribution around binary masses ∼10 M ⊙ and ∼30 M ⊙. Example ready-to-use scripts and some produced data sets for this method are shared in an online repository.