Search for: All records

The final state of a binary black hole merger is predicted with high precision by numerical relativity, but could there be a simple thermodynamic principle within general relativity that governs the selection of the remnant? Using post-Newtonian relations between the mass 𝑀 (including the binding energy) and angular momentum 𝐽 of quasicircular, nonspinning binaries, we uncover a puzzling result: When the binary’s instantaneous 𝑀 and 𝐽 are mapped to those of a hypothetical Kerr black hole, the corresponding entropy exhibits a maximum during the evolution. This maximum occurs at values of 𝑀 and 𝐽, strikingly close to those of the final remnant predicted by numerical relativity. Consistent behavior is observed when using the relation between 𝑀 and 𝐽 obtained from numerical relativity evolution. Although this procedure is somewhat ad hoc, the agreement between the masses and spins of the final state obtained from numerical relativity and the results of this maximum entropy procedure is remarkable, with agreement to within a few percent when using either post-Newtonian or numerical relativity results for 𝑀 and 𝐽. These findings allow us to propose an entropy maximization conjecture for binary black hole mergers, hinting that thermodynamic principles may govern the selection of the final black hole state. 

Black holes (BHs) with masses between $\sim 3–5M_{\odot }$, produced by a binary neutron star (BNS) merger, can further pair up with a neutron star or BH and merge again within a Hubble time. However, the astrophysical environments in which this can happen and the rate of such mergers are open questions in astrophysics. Gravitational waves may play an important role in answering these questions. In this context, we discuss the possibility that the primary of the recent LIGO-Virgo-KAGRA binary GW230529_181500 (GW230529, in short) is the product of a previous BNS merger. Invoking numerical relativity (NR)-based fitting formulas that map the binary constituents’ masses and tidal deformabilities to the mass, spin, and kick velocity of the remnant BH, we investigate the potential parents of GW230529’s primary. Our calculations using NR fits based on BNS simulations reveal that the remnant of a high-mass BNS merger similar to GW190425 is consistent with the primary of GW230529. This argument is further strengthened by the gravitational wave-based merger rate estimation of GW190425-like and GW230529-like populations. We show that around 18% (median) of the GW190425-like remnants could become the primary component in GW230529-like mergers. The dimensionless tidal deformability parameter of the heavier neutron star in the parent binary is constrained to ${67}_{-61}^{+163}$at 90% credibility. Using estimates of the gravitational-wave kick imparted to the remnant, we also discuss the astrophysical environments in which these types of mergers can take place and the implications for their future observations. 

The production of black holes with masses between ∼50⁢𝑀⊙−130⁢𝑀⊙ is believed to be prohibited by stellar processes due to (pulsational) pair-instability supernovae. Hierarchical mergers of black holes in dense star clusters are proposed as a mechanism to explain the observations of binary black holes with component masses in this range by LIGO/Virgo. We study the efficiency with which hierarchical mergers can produce higher and higher masses using a simple model of the forward evolution of binary black hole populations in gravitationally bound systems like stellar clusters. The model relies on pairing probability and initial mass functions for the black hole population, along with numerical relativity fitting formulas for the mass, spin, and kick speed of the merger remnant. We carry out an extensive comparison of the predictions of our model with clusterBHBdynamics (cBHBD) model, a fast method for the evolution of star clusters and black holes therein. For this comparison, we consider three different pairing functions of black holes and consider simulations from high- and low-metallicity cluster environments from cBHBD. We find good agreements between our model and the cBHBD results when the pairing probability of binaries depends on both total mass and mass ratio. We also assess the efficiency of hierarchical mergers as a function of merger generation and derive the mass distribution of black holes using our model. We find that the multimodal features in the observed binary black hole mass spectrum—revealed by the nonparametric population models—can be interpreted by invoking the hierarchical merger scenario in dense, metal-rich, stellar environments. Further, the two subdominant peaks in the GWTC-3 component mass spectrum are consistent with second and third-generation mergers in metal-rich, dense environments. With more binary black hole detections, our model could be used to infer the black hole initial mass function and pairing probability exponents. 

Abstract This study investigates the origins of GW230529, delving into its formation from massive stars within isolated binary systems. Utilizing population-synthesis models, we present compelling evidence that the neutron star component forms second. However, the event’s low signal-to-noise ratio introduces complexities in identifying the underlying physical mechanisms driving its formation. Augmenting our analysis with insights from numerical relativity, we estimate the final black hole mass and spin to be approximately 5.3M_⊙and 0.53, respectively. Furthermore, we employ the obtained posterior samples to calculate the ejecta mass and kilonova light curves resulting fromr-process nucleosynthesis. We find the ejecta mass to be within 0–0.06M_⊙, contingent on the neutron star equation of state. The peak brightness of the kilonova light curves indicates that targeted follow-up observations with a Rubin-like observatory may have detected this emission. 

Abstract We propose a Bayesian inference framework to predict the merger history of LIGO-Virgo binary black holes (BHs), whose binary components may have undergone hierarchical mergers in the past. The framework relies on numerical relativity predictions for the mass, spin, and kick velocity of the remnant BHs. This proposed framework computes the masses, spins, and kicks imparted to the remnant of the parent binaries, given the initial masses and spin magnitudes of the binary constituents. We validate our approach by performing an “injection study” based on a constructed sequence of hierarchically formed binaries. Noise is added to the final binary in the sequence, and the parameters of the “parent” and “grandparent” binaries in the merger chain are then reconstructed. This method is then applied to three GWTC-3 events: GW190521, GW200220_061928, and GW190426_190642. These events were selected because at least one of the binary companions lies in the putative pair-instability supernova mass gap, in which stellar processes alone cannot produce BHs. Hierarchical mergers offer a natural explanation for the formation of BHs in the pair-instability mass gap. We use the backward evolution framework to predict the parameters of the parents of the primary companion of these three binaries. For instance, the parent binary of GW190521 has masses ${72}_{-22}^{+32}{M}_{\odot }$and ${31}_{-23}^{+24}{M}_{\odot }$within the 90% credible interval. Astrophysical environments with escape speeds ≥100 km s⁻¹are preferred sites to host these events. Our approach can be readily applied to future high-mass gravitational wave events to predict their formation history under the hierarchical merger assumption. 

We present a comprehensive assessment of multiparameter tests of general relativity (GR) in the inspiral regime of compact binary coalescences using principal component analysis (PCA). Our analysis is based on an extensive set of simulated gravitational-wave (GW) signals, including both general relativistic and non-GR sources, injected into zero-noise data colored by the noise power spectral densities of the LIGO and Virgo GW detectors at their designed sensitivities. We evaluate the performance of PCA-based methods in the context of two established frameworks: and . For GR-consistent signals, we find that PCA enables stringent constraints on potential deviations from GR, even in the presence of multiple free parameters. Applying the method to simulated signals that explicitly violate GR, we demonstrate that PCA is effective at identifying such deviations. We further test the method using numerical relativity waveforms of eccentric binary black hole systems and show that missing physical effects—such as orbital eccentricity—can lead to apparent violations of GR if not properly included in the waveform models used for analysis. Finally, we apply our PCA-based test to selected real gravitational-wave events from GWTC-3, including GW190814 and GW190412. We present joint constraints from selected binary black hole events in GWTC-3, finding that the 90% credible bound on the most informative PCA parameter is ${0.03}_{-0.08}^{+0.08}$in the framework and $-0.01_{-0.04}^{+0.05}$in the framework, both of which are consistent with GR. These results highlight the sensitivity and robustness of the PCA-based approach and demonstrate its readiness for application to future observational data from the fourth observing runs of LIGO, Virgo, and KAGRA. 

Principal component analysis (PCA) is an efficient tool to optimize multiparameter tests of general relativity (GR), wherein one looks for simultaneous deviations in multiple post-Newtonian phasing coefficients. This is accomplished by introducing non-GR deformation parameters in the phase evolution of the gravitational-wave templates used in the analysis. A PCA is performed to construct the “best-measured” linear combinations of the deformation parameters. This helps to set stringent limits on deviations from GR and to more readily detect possible beyond-GR physics. In this paper, we study the effectiveness of this method with the proposed next-generation gravitational-wave detectors, Cosmic Explorer (CE) and Einstein Telescope (ET). For compact binaries at a luminosity distance of 500 Mpc and the detector-frame total mass in the range 20–200M⊙, CE can measure the most dominant linear combination with a 1-σ uncertainty ∼0.1% and the next two subdominant linear combinations with a 1-σ uncertainty of ≤ 10%. For a specific range of masses, constraints from ET are better by a factor of a few than CE. This improvement is because of the improved low frequency sensitivity of ET compared to CE (between 1–7 Hz). In addition, we explain the sensitivity of the PCA parameters to the different post-Newtonian deformation parameters and discuss their variation with total mass. We also discuss a criterion for quantifying the number of most dominant linear combinations that capture the information in the signal up to a threshold.