Unveiling the Merger Structure of Black Hole Binaries in Generic Planar Orbits.

The precise modeling of binary black hole coalescences in generic planar orbits is a crucial step to disentangle dynamical and isolated binary formation channels through gravitational-wave observations. The merger regime of such coalescences exhibits a significantly higher complexity compared to the quasicircular case, and cannot be readily described through standard parametrizations in terms of eccentricity and anomaly. In the spirit of the effective one body formalism, we build on the study of the test-mass limit, and introduce a new modeling strategy to describe the general-relativistic dynamics of two-body systems in generic orbits. This is achieved through gauge-invariant combinations of the binary energy and angular momentum, such as a dynamical "impact parameter" at merger. These variables reveal simple "quasi-universal" structures of the pivotal merger parameters, allowing us to build an accurate analytical representation of generic (bounded and dynamically bounded) orbital configurations. We demonstrate the validity of these analytical relations using 311 numerical simulations of bounded noncircular binaries with progenitors from the RIT and SXS catalogs, together with a custom dataset of dynamical captures generated using the Einstein Toolkit, and test-mass data in bound orbits. Our modeling strategy lays the foundations of accurate and complete waveform models for systems in arbitrary orbits, bolstering observational explorations of dynamical formation scenarios and the discovery of new classes of gravitational wave sources.

Neutrino Trapping and Out-of-Equilibrium Effects in Binary Neutron-Star Merger Remnants.

We study out-of-thermodynamic-equilibrium effects in neutron-star mergers with 3D general-relativistic neutrino-radiation large-eddy simulations. During mergers, the cores of the neutron stars remain cold (T∼ a few MeV) and out of thermodynamic equilibrium with trapped neutrinos originating from the hot collisional interface between the stars. However, within ∼2 to 3 ms matter and neutrinos reach equilibrium everywhere in the remnant massive neutron star. Our results show that dissipative effects, such as bulk viscosity, if present, are only active for a short window of time after the merger.

Constraints on the Maximum Densities of Neutron Stars from Postmerger Gravitational Waves with Third-Generation Observations.

Using data from 289 numerical relativity simulations of binary neutron star mergers, we identify, for the first time, a robust quasiuniversal relation connecting the postmerger peak gravitational-wave frequency and the value of the density at the center of the maximum mass nonrotating neutron star. This relation offers a new possibility for precision equation-of-state constraints with next-generation ground-based gravitational-wave interferometers. Mock Einstein Telescope observations of fiducial events indicate that Bayesian inferences can constrain the maximum density to ∼15% (90% credibility level) for a single signal at the minimum sensitivity threshold for a detection. If the postmerger signal is included in a full-spectrum (inspiral-merger-postmerger) analysis of such a signal, the pressure-density function can be tightly constrained up to the maximum density, and the maximum neutron star mass can be measured with an accuracy better than 12% (90% credibility level).

Probing the Incompressibility of Nuclear Matter at Ultrahigh Density through the Prompt Collapse of Asymmetric Neutron Star Binaries.

Using 250 neutron star merger simulations with microphysics, we explore for the first time the role of nuclear incompressibility in the prompt collapse threshold for binaries with different mass ratios. We demonstrate that observations of prompt collapse thresholds, either from binaries with two different mass ratios or with one mass ratio but combined with the knowledge of the maximum neutron star mass or compactness, will constrain the incompressibility at the maximum neutron star density K_{max} to within tens of percent. This otherwise inaccessible measure of K_{max} can potentially reveal the presence of hyperons or quarks inside neutron stars.

Black-Hole Remnants from Black-Hole-Neutron-Star Mergers.

Observations of gravitational waves and their electromagnetic counterparts may soon uncover the existence of coalescing compact binary systems formed by a stellar-mass black hole and a neutron star. These mergers result in a remnant black hole, possibly surrounded by an accretion disk. The mass and spin of the remnant black hole depend on the properties of the coalescing binary. We construct a map from the binary components to the remnant black hole using a sample of numerical-relativity simulations of different mass ratios q, (anti)aligned dimensionless spins of the black hole a_{BH}, and several neutron star equations of state. Given the binary total mass, the mass and spin of the remnant black hole can therefore be determined from the three parameters (q,a_{BH},Λ), where Λ is the tidal deformability of the neutron star. Our models also incorporate the binary black hole and test-mass limit cases and we discuss a simple extension for generic black-hole spins. We combine the remnant characterization with recent population synthesis simulations for various metallicities of the progenitor stars that generated the binary system. We predict that black-hole-neutron-star mergers produce a population of remnant black holes with masses distributed around 7 M_{⊙} and 9 M_{⊙}. For isotropic spin distributions, nonmassive accretion disks are favored: no bright electromagnetic counterparts are expected in such mergers.

Gravitational-Wave Luminosity of Binary Neutron Stars Mergers.

We study the gravitational-wave peak luminosity and radiated energy of quasicircular neutron star mergers using a large sample of numerical relativity simulations with different binary parameters and input physics. The peak luminosity for all the binaries can be described in terms of the mass ratio and of the leading-order post-Newtonian tidal parameter solely. The mergers resulting in a prompt collapse to black hole have the largest peak luminosities. However, the largest amount of energy per unit mass is radiated by mergers that produce a hypermassive neutron star or a massive neutron star remnant. We quantify the gravitational-wave luminosity of binary neutron star merger events, and set upper limits on the radiated energy and the remnant angular momentum from these events. We find that there is an empirical universal relation connecting the total gravitational radiation and the angular momentum of the remnant. Our results constrain the final spin of the remnant black hole and also indicate that stable neutron star remnant forms with super-Keplerian angular momentum.

Modeling the Complete Gravitational Wave Spectrum of Neutron Star Mergers.

In the context of neutron star mergers, we study the gravitational wave spectrum of the merger remnant using numerical relativity simulations. Postmerger spectra are characterized by a main peak frequency f2 related to the particular structure and dynamics of the remnant hot hypermassive neutron star. We show that f(2) is correlated with the tidal coupling constant κ(2)^T that characterizes the binary tidal interactions during the late-inspiral merger. The relation f(2)(κ(2)^T) depends very weakly on the binary total mass, mass ratio, equation of state, and thermal effects. This observation opens up the possibility of developing a model of the gravitational spectrum of every merger unifying the late-inspiral and postmerger descriptions.

Modeling the Dynamics of Tidally Interacting Binary Neutron Stars up to the Merger.

The data analysis of the gravitational wave signals emitted by coalescing neutron star binaries requires the availability of an accurate analytical representation of the dynamics and waveforms of these systems. We propose an effective-one-body model that describes the general relativistic dynamics of neutron star binaries from the early inspiral up to the merger. Our effective-one-body model incorporates an enhanced attractive tidal potential motivated by recent analytical advances in the post-Newtonian and gravitational self-force description of relativistic tidal interactions. No fitting parameters are introduced for the description of tidal interaction in the late, strong-field dynamics. We compare the model energetics and the gravitational wave phasing with new high-resolution multiorbit numerical relativity simulations of equal-mass configurations with different equations of state. We find agreement within the uncertainty of the numerical data for all configurations. Our model is the first semianalytical model that captures the tidal amplification effects close to merger. It thereby provides the most accurate analytical representation of binary neutron star dynamics and waveforms currently available.