Black Hole Binaries Research Articles

Connecting low-redshift LISA massive black hole mergers to the nHz stochastic gravitational wave background

Pulsar Timing Array (PTA) experiments worldwide recently reported evidence of a nHz stochastic gravitational wave background (sGWB) compatible with the existence of slowly inspiralling massive black hole (MBH) binaries (MBHBs). The shape of the signal contains valuable information about the evolution of z < 1 MBHs above 108 M⊙, suggesting a faster dynamical evolution of MBHBs towards the gravitational-wave-driven inspiral or a larger MBH growth than usually assumed. In this work, we investigate if the nHz sGWB could also provide constraints on the population of merging lower-mass MBHBs (< 107 M⊙) detectable by LISA. To this end, we use the L-Galaxies semi-analytical model applied to the Millennium suite of simulations. We generate a population of MBHs compatible simultaneously with current electromagnetic and nHz sGWB constraints by including the possibility that, in favourable environments, MBHs can accrete gas beyond the Eddington limit. The predictions of this new model for the sGWB show that the global (integrated up to high-z) LISA detection rate is not significantly affected when compared to a fiducial model whose nHz sGWB signal is ∼2 times smaller. In both cases, the global rate yields ∼12 yr−1 and is dominated by systems of 105 − 6 M⊙. The main differences are limited to low-z (z < 3), high-mass (> 106 M⊙) LISA MBHBs. The model compatible with the latest PTA results predicts up to ∼1.6 times more detections, with a rate of ∼1 yr−1. We find that these LISA MBHB systems have 50% probability of shining with bolometric luminosities > 1043 erg s−1. Hence, in case PTA results are confirmed and given the current MBH modelling, our findings suggest there will be higher chances to perform multimessenger studies with LISA MBHB than previously expected.

Physics-inspired spatiotemporal-graph AI ensemble for the detection of higher order wave mode signals of spinning binary black hole mergers

We present a new class of AI models for the detection of quasi-circular, spinning, non-precessing binary black hole mergers whose waveforms include the higher order gravitational wave modes (ℓ,|m|)={(2,2),(2,1),(3,3),(3,2),(4,4)} , and mode mixing effects in the ℓ=3,|m|=2 harmonics. These AI models combine hybrid dilated convolution neural networks to accurately model both short- and long-range temporal sequential information of gravitational waves; and graph neural networks to capture spatial correlations among gravitational wave observatories to consistently describe and identify the presence of a signal in a three detector network encompassing the Advanced LIGO and Virgo detectors. We first trained these spatiotemporal-graph AI models using synthetic noise, using 1.2 million modeled waveforms to densely sample this signal manifold, within 1.7 h using 256 NVIDIA A100 GPUs in the Polaris supercomputer at the Argonne Leadership Computing Facility. This distributed training approach exhibited optimal classification performance, and strong scaling up to 512 NVIDIA A100 GPUs. With these AI ensembles we processed data from a three detector network, and found that an ensemble of 4 AI models achieves state-of-the-art performance for signal detection, and reports two misclassifications for every decade of searched data. We distributed AI inference over 128 GPUs in the Polaris supercomputer and 128 nodes in the Theta supercomputer, and completed the processing of a decade of gravitational wave data from a three detector network within 3.5 h. Finally, we fine-tuned these AI ensembles to process the entire month of February 2020, which is part of the O3b LIGO/Virgo observation run, and found 6 gravitational waves, concurrently identified in Advanced LIGO and Advanced Virgo data, and zero false positives. This analysis was completed in one hour using one NVIDIA A100 GPU.

Emission signatures from sub-parsec post-Newtonian binaries embedded in circumbinary discs

We studied the dynamical evolution of quasi-circular, equal-mass massive black hole binaries embedded in circumbinary discs from separations of ∼100 Rg down to the merger, following the post merger evolution. The binary orbit evolves owing to the presence of the gaseous disc and the addition of post-Newtonian (PN) corrections up to the 2.5 PN order, therefore including the dissipative gravitational wave back reaction. We investigated two cases of relatively cold and warm circumbinary discs, with aspect ratios of H/R = 0.03, 0.1, respectively, employing 3D hyper-Lagrangian resolution simulations with the GIZMO-MFM code. We extracted spectral energy distributions and light curves in different frequency bands (i.e. X-ray, optical, and UV) from the simulations. We find a clear two orders of magnitude drop in the X-ray flux right before merger if the disc is warm, while we identify a significant increase in the UV flux regardless of the disc temperature. The optical flux shows clear distinctive modulations on the binary orbital period and on the cavity edge period, regardless of the disc temperature. We find that the presence of a cold disc can accelerate the coalescence of the binary by up to 130 s over the last five days of inspiral, implying a phase shift accumulation of about 0.14 radians compared to the binary evolution in vacuum. These differences are triggered by the presence of the gaseous disc and might have implications on the waveforms that can be detected in principle. We discuss the implications that these distinctive signatures might have for existing and upcoming time domain surveys and for multi-messenger astronomy.

Determining the Hubble constant with AGN-assisted black hole mergers

ABSTRACT Gravitational waves from neutron star mergers have long been considered a promising way to measure the Hubble constant, H0, which describes the local expansion rate of the Universe. While black hole mergers are more abundantly observed, their expected lack of electromagnetic emission and poor gravitational-wave localization make them less well suited for measuring H0. Black hole mergers within the discs of Active Galactic Nuclei (AGN) could be an exception. Accretion from the AGN disc may produce an electromagnetic signal, pointing observers to the host galaxy. Alternatively, the low number density of AGNs could help identify the host galaxy of $1{\!-\!}5~{{\ \rm per\ cent}}$ of mergers. Here we show that black hole mergers in AGN discs may be a sensitive way to determine H0 with gravitational waves. If 1 per cent (10 per cent) of LIGO’s observations occur in AGN discs with identified host galaxies, we could measure H0 with 12 per cent (4 per cent) uncertainty in five years, possibly comparable to the sensitivity of neutron star mergers and set to considerably improve current gravitational wave measurements.

The Metallicity Dependence and Evolutionary Times of Merging Binary Black Holes: Combined Constraints from Individual Gravitational-wave Detections and the Stochastic Background

The advent of gravitational-wave astronomy is now allowing for the study of compact binary merger demographics throughout the Universe. This information can be leveraged as tools for understanding massive stars, their environments, and their evolution. One active question is the nature of compact binary formation: the environmental and chemical conditions required for black hole birth and the time delays experienced by binaries before they merge. Gravitational-wave events detected today, however, primarily occur at low or moderate redshifts due to current interferometer sensitivity, therefore limiting our ability to probe the high-redshift behavior of these quantities. In this work, we circumvent this limitation by using an additional source of information: observational limits on the gravitational-wave background from unresolved binaries in the distant Universe. Using current gravitational-wave data from the first three observing runs of LIGO–Virgo–KAGRA, we combine catalogs of directly detected binaries and limits on the stochastic background to constrain the time-delay distribution and metallicity dependence of binary black hole evolution. Looking to the future, we also explore how these constraints will be improved at the Advanced LIGO A+ sensitivity. We conclude that, although binary black hole formation cannot be strongly constrained with today’s data, the future detection (or a nondetection) of the gravitational-wave background with Advanced LIGO A+ will carry strong implications for the evolution of binary black holes.

Two Repeated Transient Quasiperiodic Oscillations in the γ-Ray Emission from the Blazar 3C 279

In this work, we report for the first time two repeated transient quasiperiodic oscillations (QPOs) in the γ-ray light curve of the TeV blazar 3C 279. We search for the periodicity in the light curve and estimate its confidence level using the weighted wavelet Z-transform, the Lomb–Scargle periodogram, and the REDFIT techniques. The main results are as follows: (1) a QPO of ∼33 days (>2.5σ) is found during the flare of 117 days (MJD 55008–55125) from 2009 June to November. Interestingly, the same QPO (∼39 days) reappeared in the flaring duration from MJD 59430 to 59585, with the confidence level of >4σ. (2) Another transient QPO of ∼91 days with a significance of >3.8σ is found during a period with 455 days (MJD 58430–58985) from 2019 February to 2020 May. Under the premise of considering the QPOs reported in the literature, the QPO of ∼40 days is repeated three times and the QPO of ∼91 days is repeated twice. We discuss several physical models explaining the QPOs of this blazar. Our study may suggest that the two QPOs originate from the twin jets of the binary black holes at the center of 3C 279. The repeated occurrence of QPOs of a similar scale strongly supports the geometric scenario of a blob spiraling within the jet. Furthermore, the hypothesis of a sheath in the jet may also be a potential explanation for the repetitive γ-ray flare patterns observed in the light curve.

Can Supercooled Phase Transitions Explain the Gravitational Wave Background Observed by Pulsar Timing Arrays?

Several pulsar timing array collaborations recently reported evidence of a stochastic gravitational wave background (SGWB) at nHz frequencies. While the SGWB could originate from the merger of supermassive black holes, it could be a signature of new physics near the 100MeV scale. Supercooled first-order phase transitions (FOPTs) that end at the 100MeV scale are intriguing explanations, because they could connect the nHz signal to new physics at the electroweak scale or beyond. Here, however, we provide a clear demonstration that it is not simple to create a nHz signal from a supercooled phase transition, due to two crucial issues that could rule out many proposed supercooled explanations and should be checked. As an example, we use a model based on nonlinearly realized electroweak symmetry that has been cited as evidence for a supercooled explanation. First, we show that a FOPT cannot complete for the required transition temperature of around 100MeV. Such supercooling implies a period of vacuum domination that hinders bubble percolation and transition completion. Second, we show that even if completion is not required or if this constraint is evaded, the Universe typically reheats to the scale of any physics driving the FOPT. The hierarchy between the transition and reheating temperature makes it challenging to compute the spectrum of the SGWB.

Tracking Supermassive Black Hole Mergers from kpc to sub-pc Scales with AXIS

We present an analysis showcasing how the Advanced X-ray Imaging Satellite (AXIS), a proposed NASA Probe-class mission, will significantly increase our understanding of supermassive black holes undergoing mergers—from kpc to sub-pc scales. In particular, the AXIS point spread function, field of view, and effective area are expected to result in (1) the detection of hundreds to thousands of new dual AGNs across the redshift range 0<z<5 and (2) blind searches for binary AGNs that are exhibiting merger signatures in their light curves and spectra. AXIS will detect some of the highest-redshift dual AGNs to date, over a large range of physical separations. The large sample of AGN pairs detected by AXIS (over a magnitude more than currently known) will result in the first X-ray study that quantifies the frequency of dual AGNs as a function of redshift up to z=4.

Pre-supernova evolution and final fate of stellar mergers and accretors of binary mass transfer

The majority of massive stars are expected to exchange mass or merge with a companion during their lives. This immediately implies that most supernovae (SNe) are from such post-mass-exchange objects. Here, we explore how mass accretion and merging affect the pre-SN structures of stars and their final fates. To this end, we modelled these complex processes by rapid mass accretion onto stars of different evolutionary stages and followed their evolution up to iron core collapse. We used the stellar evolution code MESA and inferred the outcome of core-collapse using a neutrino-driven SN model. Our models cover initial masses from 11 to 70 M⊙ and the accreted mass ranges from 10−200% of the initial mass. All models are non-rotating and for solar metallicity. The rapid accretion model offers a systematic way to approach the landscape of mass accretion and stellar mergers. It is naturally limited in scope and serves as a clean zeroth order baseline for these processes. We find that mass accretion, in particular onto post-main-sequence (post-MS) stars, can lead to a long-lived blue supergiant (BSG) phase during which stars burn helium in their cores. In comparison to genuine single stars, post-MS accretors have small core-to-total mass ratios, regardless of whether they end their lives as BSGs or cool supergiants (CSGs), and they can have genuinely different pre-SN core structures. As in single and binary-stripped stars, we find black-hole (BH) formation for the same characteristic CO core masses MCO of ≈7 M⊙ and ≳13 M⊙. In models with the largest mass accretion, the BH formation landscape as a function of MCO is shifted by about 0.5 M⊙ to lower masses, that is, such accretors are more difficult to explode. We find a tight relation between our neutron-star (NS) masses and the central entropy of the pre-SN models in all accretors and single stars, suggesting a universal relation that is independent of the evolutionary history of stars. Post-MS accretors explode both as BSGs and CSGs, and we show how to understand their pre-SN locations in the Hertzsprung-Russell (HR) diagram. Accretors exploding as CSGs can have much higher envelope masses than single stars. Some BSGs that avoid the luminous-blue-variable (LBV) regime in the HR diagram are predicted to collapse into BHs of up to 50 M⊙, while others explode in SNe and eject up to 40 M⊙, greatly exceeding ejecta masses from single stars. Both the BH and SN ejecta masses increase to about 80 M⊙ in our models when allowing for multiple mergers, for example, in initial triple-star systems, and they can be even higher at lower metallicities. Such high BH masses may fall into the pair-instability-SN mass gap and could help explain binary BH mergers involving very massive BHs as observed in GW190521. We further find that some of the BSG models explode as LBVs, which may lead to interacting SNe and possibly even superluminous SNe.

Constraining the LVK AGN channel with black hole spins

ABSTRACT Merging black holes (BHs) are expected to produce remnants with large dimensionless spin parameters (aspin ∼ 0.7). However, gravitational wave (GW) observations with LIGO–Virgo–Kagra (LVK) suggest that merging BHs are consistent with modestly positive but not high spin (aspin ∼ 0.2), causing tension with models suggesting that high-mass mergers are produced by hierarchical merger channels. Some BHs also show evidence for strong in-plane spin components. Here, we point out that spin-down of BHs due to eccentric prograde post-merger orbits within the gas of an active galactic nucleus (AGN) disc can yield BHs with masses in the upper mass gap, but only modestly positive aspin, and thus observations of BHs with low spin do not rule out hierarchical models. We also point out that the fraction of binary black hole (BBH) mergers with significant in-plane spin components is a strong test of interactions between disc BBHs and nuclear spheroid orbiters. Spin magnitude and spin tilt constraints from LVK observations of BBHs are an excellent test of dynamics of BHs in AGN discs, disc properties, and the nuclear clusters interacting with AGNs.

Effect of Wave Dark Matter on Equal Mass Black Hole Mergers

For dark matter to be detectable with gravitational waves from binary black holes, it must reach higher than average densities in their vicinity. In the case of light (wavelike) dark matter, the density of dark matter between the binary can be significantly enhanced by accretion from the surrounding environment. Here we show that the resulting dephasing effect on the last ten orbits of an equal mass binary is maximized when the Compton wavelength of the scalar particle is comparable to the orbital separation, 2π/μ∼d. The phenomenology of the effect is different from the channels that are usually discussed, where dynamical friction (along the orbital path) and radiation of energy and angular momentum drive the dephasing, and is rather dominated by the radial force (the spacetime curvature in the radial direction) towards the overdensity between the black holes. While our numerical studies limit us to scales of the same order, this effect may persist at larger separations and/or particle masses, playing a significant role in the merger history of binaries. Published by the American Physical Society 2024

Testing space-time non-commutativity with TianQin

The direct detection of gravitational waves offers a powerful tool to explore the nature of gravity and the structure of space-time. This paper focuses on the capabilities of space-based gravitational wave detectors in testing space-time non-commutativity. Our findings indicate that TianQin has the potential to impose constraints on the non-commutative scale at a sub-Planckian level using massive black hole binaries. Additionally, we have developed a pipeline tailored to this specific topic.

Multi-messenger astrophysics of black holes and neutron stars as probed by ground-based gravitational wave detectors: from present to future

The ground-based gravitational wave (GW) detectors LIGO and Virgo have enabled the birth of multi-messenger GW astronomy via the detection of GWs from merging stellar-mass black holes (BHs) and neutron stars (NSs). GW170817, the first binary NS merger detected in GWs and all bands of the electromagnetic spectrum, is an outstanding example of the impact that GW discoveries can have on multi-messenger astronomy. Yet, GW170817 is only one of the many and varied multi-messenger sources that can be unveiled using ground-based GW detectors. In this contribution, we summarize key open questions in the astrophysics of stellar-mass BHs and NSs that can be answered using current and future-generation ground-based GW detectors, and highlight the potential for new multi-messenger discoveries ahead.

Shadows and Photon Rings of Binary Black Holes

Shadows and Photon Rings of Binary Black Holes

Analytically separating the source of the Teukolsky equation

Recent gravitational wave detections from black hole mergers have underscored the critical role black hole perturbation theory and the Teukolsky equation play in understanding the behavior of black holes. The separable nature of the Teukolsky equation has long been leveraged to study the vacuum linear Teukolsky equation; however, as theory and measurements advance, solving the sourced Teukolsky equation is becoming a frontier of research. In particular, second-order calculations, such as in quasinormal mode and self-force problems, have extended sources. This paper presents a novel method for analytically separating the Teukolsky equation’s source, aimed to improve efficiency. Separating the source is a nontrivial problem due to the angular and radial mixing of generic quantities in Kerr spacetime. We provide a proof-of-concept demonstration of our method and show that it is accurate, separating the Teukolsky source produced by the stress-energy tensor of an ideal gas cloud surrounding a Kerr black hole. The detailed application of our method is provided in an accompanying notebook. Our approach opens up a new avenue for accurate black hole perturbation theory calculations with sources in both the time and frequency domain. Published by the American Physical Society 2024

Generating higher order modes from binary black hole mergers with machine learning

Generating higher order modes from binary black hole mergers with machine learning

A JVLA, LOFAR, e-Merlin, VLBA, and EVN study of RBS 797: can binary supermassive black holes explain the outburst history of the central radio galaxy?

The multifaceted central radio galaxy of the cluster RBS 797 shows several episodes of jet activity in multiple directions. We wish to understand the causes behind these dramatic misalignments and measure the timescales of the successive outbursts. We present a multifrequency (144 MHz -- 9 GHz) and multiscale (5 pc -- $50$ kpc) investigation of the central radio galaxy in RBS 797, by means of JVLA, LOFAR (with international stations), e-Merlin, VLBA, and EVN data. We investigate the morphological and spectral properties of the radio lobes, the jets, and the active core. We confirm the co-spatiality of the radio lobes with the four perpendicular X-ray cavities previously discovered. The radiative ages of the east-west lobes ($31.4 Myr) and of the north-south lobes ($32.1 Myr) support a coeval origin of the perpendicular outbursts, which also have similar active phase duration (sim 12 Myr). Based on the analysis of the inner north-south jets (on scales of $ kpc), we (a) confirm the S-shaped jet morphology; (b) show the presence of two hotspots per jet with a similar spectral index; and (c) estimate the age of the twisting north-south jets to be less than $ Myr. Based on these results, we determine that jet precession, with a period sim 9 Myr, half-opening angle sim circ $, and jet advance speed sim 0.01$c$, can explain the properties of the north-south jets. We also find that the synchrotron injection index has steepened from the large, older outbursts ($ to the younger S-shaped jets ($ possibly due to a transition from FR I-like to FR II-like activity. The e-Merlin, VLBA, and EVN data reveal a single, compact core at the heart of RBS 797, surrounded by extended radio emission whose orientation depends on the spatial scale sampled by the data. We explore several engine-based scenarios to explain these results. Piecing together the available evidence, we argue that RBS 797 likely hosts (or hosted) binary active supermassive black holes (SMBHs). The detection of a single component in the VLBA and EVN data is still consistent with this interpretation, since the predicted separation of the binary SMBHs (leq 0.6 pc) is an order of magnitude smaller than the resolution of the available radio data (5 pc).

On the formation of a 33 solar-mass black hole in a low-metallicity binary

A 33M⊙ black hole (BH) was recently discovered in an 11.6-year binary only 590 pc from the Sun. The system, Gaia BH3, contains a 0.8M⊙ low-metallicity giant () that is a member of the ED-2 stellar stream. This paper investigates whether the system could have formed via isolated binary evolution. I construct evolutionary models for metal-poor massive stars with initial masses ranging from 35−55M⊙, which reach maximum radii of 1150−1800R⊙ as red supergiants. I then explore what combinations of initial orbit, mass loss, and natal kick can produce the period and eccentricity of Gaia BH3. Initial orbits wide enough to accommodate the BH progenitor as a red supergiant can match the observed period and eccentricity, but only if the BH formed with a significant natal kick (). These models are disfavored because such a kick would have ejected the binary from the ED-2 progenitor cluster. I conclude that Gaia BH3 likely formed through dynamical interactions, unless the BH progenitor did not expand to red supergiant dimensions. Only about 1 in 10,000 stars in the solar neighborhood have metallicities as low as Gaia BH3. This suggests that BH companions are dramatically over-represented at low-metallicity, though caveats related to small number statistics apply. The fact that the luminous star in Gaia BH3 has been a giant – greatly boosting its detectability – only for $$1% of the time since the system’s formation implies that additional massive BHs remain to be discovered with only moderately fainter companions. Both isolated and dynamically-formed BH binaries with orbits similar to Gaia BH3 are likely to be discovered in Gaia DR4.

GRMHD simulations of accretion flows onto massive binary black hole mergers embedded in a thin slab of gas

GRMHD simulations of accretion flows onto massive binary black hole mergers embedded in a thin slab of gas

Cosmography with next-generation gravitational wave detectors

Advancements in cosmology through next-generation (XG) ground-based gravitational wave (GW) observatories will bring in a paradigm shift. We explore the pivotal role that GW standard sirens will play in inferring cosmological parameters with XG observatories, not only achieving exquisite precision but also opening up unprecedented redshifts. We examine the merits and the systematic biases involved in GW standard sirens utilizing binary black holes, binary neutron stars, and neutron star-black hole mergers. Further, we estimate the precision of bright sirens, golden dark sirens, and spectral sirens for these binary coalescences and compare the abilities of various XG observatories (A ♯ , cosmic explorer, Einstein telescope, and their possible networks). When combining different sirens, we find sub-percent precision over more than 10 billion years of cosmic evolution for the Hubble expansion rate H(z). This work presents a broad view of opportunities to precisely measure the cosmic expansion rate, decipher the elusive dark energy and dark matter, and potentially discover new physics in the uncharted Universe with XG GW detectors.