Binary Mergers Research Articles

Imprints of massive black-hole binaries on neighbouring decihertz gravitational-wave sources

The most massive black holes in our Universe form binaries at the centre of merging galaxies. The recent evidence for a gravitational-wave (GW) background from pulsar timing may constitute the first observation that these supermassive black-hole binaries (SMBHBs) merge. Yet, the most massive SMBHBs are out of reach of interferometric GW detectors and are exceedingly difficult to resolve individually with pulsar timing. These limitations call for unexplored strategies to detect individual SMBHBs in the uncharted frequency band ≲10−5 Hz to establish their abundance and decipher the coevolution with their host galaxies. Here we show that SMBHBs imprint detectable long-term modulations on GWs from stellar-mass binaries residing in the same galaxy at a distance d ≲ 1 kpc. We determine that proposed decihertz GW interferometers sensitive to numerous stellar-mass binaries could uncover modulations from ~O(10−1–104) SMBHBs with masses ~O(107–108) M⊙out to redshift z ≈ 3.5. This offers a unique opportunity to map the population of SMBHBs through cosmic time, which might remain inaccessible otherwise.

Prospects for Time-Domain and Multi-Messenger Science with AXIS

The Advanced X-ray Imaging Satellite (AXIS) promises revolutionary science in the X-ray and multi-messenger time domain. AXIS will leverage excellent spatial resolution (<1.5 arcsec), sensitivity (80× that of Swift), and a large collecting area (5–10× that of Chandra) across a 24-arcmin diameter field of view at soft X-ray energies (0.3–10.0 keV) to discover and characterize a wide range of X-ray transients from supernova-shock breakouts to tidal disruption events to highly variable supermassive black holes. The observatory’s ability to localize and monitor faint X-ray sources opens up new opportunities to hunt for counterparts to distant binary neutron star mergers, fast radio bursts, and exotic phenomena like fast X-ray transients. AXIS will offer a response time of <2 h to community alerts, enabling studies of gravitational wave sources, high-energy neutrino emitters, X-ray binaries, magnetars, and other targets of opportunity. This white paper highlights some of the discovery science that will be driven by AXIS in this burgeoning field of time domain and multi-messenger astrophysics. This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at the AXIS website.

Revisiting the Tertiary-induced Binary Black Hole Mergers: The Role of Superthermal Wide Tertiary Eccentricity Distributions

Recent studies show that the eccentricity distribution of wide binaries (semimajor axis ≳103 au) observed by Gaia tends to favor large eccentricities more strongly than the canonical thermal distribution (P(e) ∝ e)—such distributions are termed “superthermal.” Motivated by this observation, we revisit the formation channel of black hole (BH) binary mergers in triple stellar systems and study the impact of superthermal eccentricity distributions in the outer binaries. We explore the persistence of the highly eccentric outer orbits after each component in a stellar triple has undergone mass loss due to supernova explosions. We find that the outer eccentricity distribution can remain significantly superthermal for modestly hierarchical BH triples satisfying a in/a out ≳ 0.005 (where a in and a out are the semimajor axes of the inner and outer orbits), and are otherwise shaped by mass-loss induced kicks and dynamical instability. We then study the impact of these different outer eccentricity distributions of the remaining BH triples on mergers via the tertiary-induced channel. Of interest, we find that mergers can sometimes be produced even when the initial stellar orbits are near alignment (not subject to the von-Zeipel–Lidov–Kozai effect; ZLK effect) as long as the system is sufficiently hierarchical. On the other hand, although the impact of the octupole-order ZLK effect is much greater when the outer binary is more eccentric, we find that the merger fraction only changes modestly for extreme outer eccentricity distributions because the largest eccentricities tend to lead to dynamical instability.

Did Binary Neutron Star Merger GW170817 Leave Behind a Long-lived Neutron Star?

We consider the observational implications of the binary neutron star merger GW170817, which may have left behind a rapidly rotating massive neutron star that launches a relativistic, equatorial outflow as well as a jet. We show that if the equatorial outflow (ring) is highly beamed in the equatorial plane, its luminosity can be “hidden” from view until late times, even if carrying a significant fraction of the spin-down energy of the merger remnant. This hidden ring reveals itself as a rebrightening in the light curve once it slows down enough for the Earth to be within the ring’s relativistic beaming solid angle. We compute semianalytic light curves using this model and find they are in agreement with the observations thus far, and we provide predictions for the ensuing afterglow.

Impact of Higher Harmonics of Gravitational Radiation on the Population Inference of Binary Black Holes

Templates modeling just the dominant mode of gravitational radiation are generally sufficient for the unbiased parameter inference of near-equal-mass compact binary mergers. However, neglecting the subdominant modes can bias the inference if the binary is significantly asymmetric, very massive, or has misaligned spins. In this work, we explore if neglecting these subdominant modes in the parameter estimation of nonspinning binary black hole mergers can bias the inference of their population-level properties such as mass and merger redshift distributions. Assuming the design sensitivity of the advanced LIGO-Virgo detector network, we find that neglecting subdominant modes will not cause a significant bias in the population inference, although including them will provide more precise estimates. This is primarily because asymmetric binaries are expected to be rarer in our detected sample, due to their intrinsic rareness and the observational selection effects. The increased precision in the measurement of the maximum black hole mass can help in better constraining the upper mass gap in the mass spectrum.

GW230529_181500: a potential primordial binary black hole merger in the mass gap

During the fourth observing run of the LIGO-Virgo-KAGRA detector network, the LIGO Livingston observatory detected a coalescing compact binary, GW230529_181500, with component masses of 2.5–4.5 M ⊙ and 1.2–2.0 M ⊙ at the 90% credible level. The gravitational-wave data alone is insufficient to determine whether the components are neutron stars or black holes. In this paper, we propose that GW230529_181500 originated from the merger of two primordial black holes (PBHs). We estimate a merger rate of 5.0+47.0 -4.9 Gpc-3 yr-1 for compact binary coalescences with properties similar to GW230529_181500. Assuming the source is a PBH-PBH merger, GW230529_181500-like events lead to approximately 1.7+36.2 -1.5 × 10-3 of the dark matter in the form of PBHs. The required abundance of PBHs to explain this event is consistent with existing upper limits derived from microlensing, cosmic microwave background observations and the null detection of gravitational-wave background by LIGO-Virgo-KAGRA.

Standard siren cosmology in the era of the 2.5-generation ground-based gravitational wave detectors: bright and dark sirens of LIGO Voyager and NEMO

The 2.5-generation (2.5G) ground-based gravitational wave (GW) detectors LIGO Voyager and NEMO are expected to be operational in the late 2020s and early 2030s. In this work, we explore the potential of GW standard sirens observed by the 2.5G GW detectors in measuring cosmological parameters, especially for the Hubble constant. Using GWs to measure cosmological parameters is inherently challenging, especially for 2.5G detectors, given their limited capability, which results in weaker constraints on cosmological parameters from the detected standard sirens. However, the measurement of the Hubble constant using standard siren observations from Voyager and NEMO is still promising. For example, using bright sirens from Voyager and NEMO can measure the Hubble constant with a precision of about 2% and 6% respectively, and using the Voyager-NEMO network can improve the precision to about 1.6%. Moreover, bright sirens can be used to break the degeneracy of cosmological parameters generated by CMB data, and to a certain extent, 2.5G detectors can also play a role in this aspect. Observations of dark sirens by 2.5G detectors can achieve relatively good results in measuring the Hubble constant, with a precision of within 2%, and if combining observations of bright and dark sirens, the precision of the Hubble constant measurement can reach about 1.4%. Finally, we also discussed the impact of the uncertainty in the binary neutron star merger rate on the estimation of cosmological parameters. We conclude that the magnificent prospect for solving the Hubble tension is worth expecting in the era of the 2.5G ground-based GW detectors.

Detection, sky localization, and early warning for binary neutron star mergers by detectors located in China of different configurations in a third generation detector network

Detection, sky localization, and early warning for binary neutron star mergers by detectors located in China of different configurations in a third generation detector network

Constraining the Dense Matter Equation of State with New NICER Mass–Radius Measurements and New Chiral Effective Field Theory Inputs

Pulse profile modeling of X-ray data from the Neutron Star Interior Composition Explorer is now enabling precision inference of neutron star mass and radius. Combined with nuclear physics constraints from chiral effective field theory (χEFT), and masses and tidal deformabilities inferred from gravitational-wave detections of binary neutron star mergers, this has led to a steady improvement in our understanding of the dense matter equation of state (EOS). Here, we consider the impact of several new results: the radius measurement for the 1.42 M ⊙ pulsar PSR J0437−4715 presented by Choudhury et al., updates to the masses and radii of PSR J0740+6620 and PSR J0030+0451, and new χEFT results for neutron star matter up to 1.5 times nuclear saturation density. Using two different high-density EOS extensions—a piecewise-polytropic (PP) model and a model based on the speed of sound in a neutron star (CS)—we find the radius of a 1.4 M ⊙ (2.0 M ⊙) neutron star to be constrained to the 95% credible ranges 12.28−0.76+0.50 km ( 12.33−1.34+0.70 km) for the PP model and 12.01−0.75+0.56 km ( 11.55−1.09+0.94 km) for the CS model. The maximum neutron star mass is predicted to be 2.15−0.16+0.14 M ⊙ and 2.08−0.16+0.28 M ⊙ for the PP and CS models, respectively. We explore the sensitivity of our results to different orders and different densities up to which χEFT is used, and show how the astrophysical observations provide constraints for the pressure at intermediate densities. Moreover, we investigate the difference R 2.0 − R 1.4 of the radius of 2 M ⊙ and 1.4 M ⊙ neutron stars within our EOS inference.

GRB 211211A: The Case for an Engine-powered over r-process-powered Blue Kilonova

The recent gamma-ray burst (GRB) GRB 211211A provides the earliest (∼5 hr) data of a kilonova (KN) event, displaying bright (∼1042 erg s−1) and blue early emission. Previously, this KN was explained using simplistic multicomponent fitting methods. Here, in order to understand the physical origin of the KN emission in GRB 211211A, we employ an analytic multizone model for r-process-powered KNe. We find that r-process-powered KN models alone cannot explain the fast temporal evolution and the spectral energy distribution (SED) of the observed emission. Specifically, (i) r-process models require high ejecta mass to match early luminosity, which overpredicts late-time emission, while (ii) red KN models that reproduce late emission underpredict early luminosity. We propose an alternative scenario involving early contributions from the GRB central engine via a late low-power jet, consistent with plateau emission in short GRBs and GeV emission detected by Fermi-LAT at ∼104 s after GRB 211211A. Such late central engine activity, with an energy budget of ∼a few percent of that of the prompt jet, combined with a single red KN ejecta component, can naturally explain the light curve and SED of the observed emission, with the late-jet–ejecta interaction reproducing the early blue emission and r-process heating reproducing the late red emission. This supports claims that late low-power engine activity after prompt emission may be common. We encourage early follow-up observations of future nearby GRBs and compact binary merger events to reveal more about the central engine of GRBs and r-process events.

A NICER View of the Nearest and Brightest Millisecond Pulsar: PSR J0437–4715

We report Bayesian inference of the mass, radius, and hot X-ray emitting region properties—using data from the Neutron Star Interior Composition ExploreR (NICER)—for the brightest rotation-powered millisecond X-ray pulsar, PSR J0437−4715. Our modeling is conditional on informative tight priors on mass, distance, and binary inclination obtained from radio pulsar timing using the Parkes Pulsar Timing Array (PPTA; Reardon et al.), and we use NICER background models to constrain the nonsource background, cross-checking with data from XMM-Newton. We assume two distinct hot emitting regions and various parameterized hot region geometries that are defined in terms of overlapping circles; while simplified, these capture many of the possibilities suggested by detailed modeling of return current heating. For the preferred model identified by our analysis, we infer a mass of M = 1.418 ± 0.037 M ⊙ (largely informed by the PPTA mass prior) and an equatorial radius of R=11.36−0.63+0.95 km, each reported as the posterior credible interval bounded by the 16% and 84% quantiles. This radius favors softer dense matter equations of state and is highly consistent with constraints derived from gravitational wave measurements of neutron star binary mergers. The hot regions are inferred to be nonantipodal and hence inconsistent with a pure centered dipole magnetic field.

Rethinking Thorne–Żytkow Object Formation: Assembly via Common Envelope in Field Binaries

Thorne–Żytkow objects (TŻOs), hypothetical merger products in which a neutron star is embedded in a stellar core, are traditionally considered steady-state configurations. Their assembly, especially through dynamical channels, is not well understood. The predominant focus in the literature has been on the observational signatures related to the evolution and long-term fate of TŻOs, with their initial formation often treated as a given. However, the foundational calculations supporting the existence of TŻOs assume nonrotating spherically symmetric initial conditions that we find to be inconsistent with a binary merger scenario. In this work, we explore the implications of postmerger dynamics in TŻO formation scenarios with field binary progenitors, specifically the role that angular momentum transport during the common envelope phase plays in constraining the possible merger products, using the tools of stellar evolution and three-dimensional hydrodynamics. We also propose an alternative steady-state outcome for these mergers: the thin-envelope TŻO, an equilibrium solution consisting of a low-mass spherical envelope supported by the accretion disk luminosity of a central stellar-mass black hole. These configurations may be of interest to upcoming time-domain surveys as potential X-ray sources that may be preceded by a series of bright transient events.

Higher memory effects in numerical simulations of binary black hole mergers

Higher memory effects in numerical simulations of binary black hole mergers

Observing black hole mergers beyond the pair-instability mass gap with next-generation gravitational wave detectors

Stellar evolution predicts the existence of a mass gap for black hole remnants produced by pair-instability supernova dynamics, whose lower and upper edges are very uncertain. We study the possibility of constraining the location of the upper end of the pair-instability mass gap, which is believed to appear around mmin∼130M⊙, using gravitational wave observations of compact binary mergers with next-generation ground-based detectors. While high metallicity may not allow for the formation of first-generation black holes on the “far side” beyond the gap, metal-poor environments containing population III stars could lead to such heavy black hole mergers. We show that, even in the presence of contamination from other merger channels, next-generation detectors will measure the location of the upper end of the mass gap with a relative precision close to Δmmin/mmin≃4%(Ndet/100)−1/2 at 90% CL, where Ndet is the number of detected mergers with both members beyond the gap. These future observations could reduce current uncertainties in nuclear and astrophysical processes controlling the location of the gap. Published by the American Physical Society 2024

Runaway Eccentricity Growth: A Pathway for Binary Black Hole Mergers in AGN Disks

Binary black holes (BBHs) embedded within the accretion disks that fuel active galactic nuclei (AGN) are promising progenitors for the source of gravitational wave (GW) events detected by LIGO/VIRGO. Several recent studies have shown that when these binaries form, they are likely to be highly eccentric and retrograde. However, many uncertainties remain concerning the orbital evolution of these binaries as they either inspiral toward merger or disassociate. Previous hydrodynamical simulations exploring their orbital evolution have been predominantly two-dimensional or have been restricted to binaries on nearly circular orbits. We present the first high-resolution, three-dimensional local shearing-box simulations of both prograde and retrograde eccentric BBHs embedded in AGN disks. We find that retrograde binaries shrink several times faster than their prograde counterparts and exhibit significant orbital eccentricity growth, the rate of which monotonically increases with binary eccentricity. Our results suggest that retrograde binaries may experience runaway orbital eccentricity growth, which may bring them close enough together at pericenter for GW emission to drive them to coalescence. Although their eccentricity is damped, prograde binaries shrink much faster than their orbital eccentricity decays, suggesting they should remain modestly eccentric as they contract toward merger. Finally, binary precession driven by the AGN disk may dominate over precession induced by the supermassive black hole depending on the binary accretion rate and its location in the AGN disk, which can subdue the evection resonance and von Ziepel–Lidov–Kozai cycles.

Dissipative Kinematics in Binary Neutron Star Mergers

We highlight the recent progress in the calculation of transport coefficients pertinent to binary neutron star mergers. Specifically, we analyze the bulk viscosity coefficient driven by both the DURCA and MURCA processes and electron transport coefficients in dense and hot plasma relevant to the merger scenario. The study considers high temperatures (T>6×1010 K) and dense environments (nB∼n0−3n0). Bulk viscosity exhibits resonant behavior, with peak values and Peak positions dependent on particle interaction rates and thermodynamic susceptibilities. Susceptibilities are calculated by modeling the nuclear matter in the density functional approach. The bulk viscosity coefficient peaks at T∼1011 K, with a compression–rarefaction oscillation dissipation time scale of 20–50 ms. Electrical transports incorporate frequency-dependent dynamical screening in quantized electron–ion scattering rates. Consequently, dynamical screening reduces the maxima of electrical and thermal conductivities, shortening corresponding dissipation time scales. These results highlight the crucial role of dissipation coefficients in understanding binary neutron star mergers.

The effect of higher harmonics on gravitational wave dark sirens

ABSTRACT The gravitational wave (GW) signal from the merger of two black holes can serve as a standard sirens for cosmological inference. However, a degeneracy exists between the luminosity distance and the inclination angle between the binary system’s orbital angular momentum and the observer’s line of sight, limiting the precise measurement of the luminosity distance. In this study, we investigate how higher harmonics affect luminosity distance estimation for third-generation (3G) GW detectors in binary black hole mergers. Our findings demonstrate that considering higher harmonics significantly enhances distance inference results compared with using only the (2, 2) mode. This improved accuracy in distance estimates also strengthens constraints on host galaxies, enabling more precise measurements of the Hubble constant. These results highlight the significant influence of higher harmonics on the range estimation accuracy of 3G ground-based GW detectors.

Gravitational waves from binary neutron star mergers with a spectral equation of state

Gravitational waves from binary neutron star mergers with a spectral equation of state

Simulation Study on Constraining Gravitational Wave Propagation Speed by Gravitational Wave and Gamma-ray Burst Joint Observation on Binary Neutron Star Mergers

Theories of modified gravity suggest that the propagation speed of gravitational waves (GW) v g may deviate from the speed of light c. A constraint can be placed on the difference between c and v g with a simple method that uses the arrival time delay between GW and electromagnetic wave simultaneously emitted from a burst event. We simulated the joint observation of GW and short gamma-ray burst signals from binary neutron star merger events in different observation campaigns, involving advanced LIGO (aLIGO) in design sensitivity and Einstein Telescope (ET) joint-detected with Fermi/GBM. As a result, the relative precision of constraint on v g can reach ∼10−17 (aLIGO) and ∼10−18 (ET), which are one and two orders of magnitude better than that from GW170817, respectively. We continue to obtain the bound of graviton mass m g ≤ 7.1(3.2) × 10−20 eV with aLIGO (ET). Applying the Standard-Model Extension test framework, the constraint on v g allows us to study the Lorentz violation in the nondispersive, nonbirefringent limit of the gravitational sector. We obtain the constraints of the dimensionless isotropic coefficients s¯00(4) at mass dimension d = 4, which are −1×10−15<s¯00(4)<9×10−17 for aLIGO and −4×10−16<s¯00(4)<8×10−18 for ET.

Jetlike structures in low-mass binary neutron star merger remnants

Jetlike structures in low-mass binary neutron star merger remnants