Black Hole Binaries Research Articles

Inferring the spin distribution of binary black holes using deep learning

The spin characteristics of black holes offer valuable insights into the evolutionary pathways of their progenitor stars. This is crucial for understanding the broader population properties of black holes. Traditional hierarchical Bayesian inference techniques employed to discern these properties often demand substantial time, and consensus regarding the spin distribution of binary black hole (BBH) systems remains elusive. In this study, leveraging observations from GWTC-3, we adopted a machine learning approach to infer the spin distribution of black holes within BBH systems. Specifically, we developed a deep neural network (DNN) and trained it using data generated from a Beta distribution. Our training strategy, involving the segregation of data into 10 bins, not only expedites model training but also enhances the versatility and adaptability of the DNN to accommodate the growing volume of gravitational wave observations. Utilizing Monte Carlo-bootstrap (MC-bootstrap) to generate observation-simulated samples, we derived spin distribution parameters: for the larger BH sample and for the smaller BH sample. Within our constraints, the distributions of component spin magnitudes suggest the likelihood of both black holes in the BBH merger possessing non-zero spin.

Self-Interacting Dark Matter Solves the Final Parsec Problem of Supermassive Black Hole Mergers.

Evidence for a stochastic gravitational wave (GW) background, plausibly originating from the merger of supermassive black holes (SMBHs), is accumulating with observations from pulsar timing arrays. An outstanding question is how inspiraling SMBHs get past the "final parsec" of separation, where they have a tendency to stall before GW emission alone can make the binary coalesce. We argue that dynamical friction from the dark matter (DM) spike surrounding the black holes is sufficient to resolve this puzzle, if the DM has a self-interaction cross section of order cm^{2}/g. The same effect leads to a softening of the GW spectrum at low frequencies as suggested by the current data. For collisionless cold DM, the friction deposits so much energy that the spike is disrupted and cannot bridge the final parsec, while for self-interacting DM, the isothermal core of the halo can act as a reservoir for the energy liberated from the SMBH orbits. A realistic velocity dependence, such as generated by the exchange of a massive mediator like a dark photon, is favored to give a good fit to the GW spectrum while providing a large enough core. A similar velocity dependence has been advocated for solving the small-scale structure problems of cold DM.

Perturbative frequency expansion for nearly monochromatic binary black holes detectable with LISA

Perturbative frequency expansion for nearly monochromatic binary black holes detectable with LISA

The X-ray rise and fall of the symbiotic recurrent nova system T CrB

ABSTRACT We present the analysis of publicly available NuSTAR, Suzaku, and XMM–Newton observations of the symbiotic recurrent nova T CrB covering the 2006.77–2022.66 yr period. The X-ray spectra are analysed by adopting a model that includes a reflection component produced by the presence of a disc that mimics the accretion disc and the immediate surrounding medium. Our best-fitting model requires this disc to have a radius of 1 au, effective thickness of 0.1 au, averaged column density 10$^{25}$ cm$^{-2}$ and orientation of 50$^{\circ }$ with respect to the line of sight. This disc is about a factor of two larger than recent estimations for the accretion disc and its presence contributes significantly via reflection to the total X-ray flux detected from T CrB, which naturally produces the emission of the 6.4 keV Fe line. Our analysis suggests that the temperature of the boundary layer evolved from 14.8 keV in the steady-state phase (before 2016), to 2.8 keV in the 2017.24 epoch, to finally stabilize to about $\sim$8 keV in the subsequent epochs. These variations in the plasma temperature of the boundary layer are attributed to the evolution of the mass accretion rate ($\dot{M}_\mathrm{acc}$), which is estimated to have an averaged value of $\dot{M}_\mathrm{acc}$ = 2.6$\times \,10^{-8}$ M$_\odot$ yr$^{-1}$ for the current active phase. The presence of emission lines in the XMM–Newton Reflection Grating Spectrometer spectrum of 2017.24 prevents from adopting a blackbody emission model to fit the soft X-ray range. Instead, we use plasma emission models that suggest the presence of adiabatically shocked gas produced by gas velocities of 110–200 km s$^{-1}$, very likely tracing jet-like ejections similar to what is found in other symbiotic systems. The analysis of X-ray and optical data together show that T CrB has a similar evolution as black hole binaries, accreting neutron stars and active galactic nuclei in the hardness–intensity diagram.

Global analysis of LISA data with Galactic binaries and massive black hole binaries

Global analysis of LISA data with Galactic binaries and massive black hole binaries

RisingTides: An Analytic Modeling Code of Tidal Effects in Binary Neutron Star Mergers

Gravitational waves produced by binary neutron star mergers offer a unique window into matter behavior under extreme conditions. In this context, we analytically model the effect of matter on gravitational waves from binary neutron star mergers. We start with a binary black hole system, leveraging the post-Newtonian formalism for the inspiral and the Backwards-one-Body model for the merger. We combine the two methods to generate a baseline waveform and we validate our results against numerical relativity simulations. Next, we integrate tidal effects in phase and amplitude to account for matter and spacetime interaction using the NRTidal model and test its accuracy against numerical relativity predictions for two equations of state, finding a mismatch around the merger. Subsequently, we lift the restriction on the coefficients to be independent of the tidal deformability and recalibrate them using the numerical relativity predictions. We obtain better fits for phase and amplitude around the merger and are able to extend the phase modeling beyond the merger. We implement our method in new open-source, user-friendly Python code, steered by a Jupyter Notebook, named RisingTides. Our research offers new perspectives on analytically modeling the effect of tides on the gravitational waves from binary neutron star mergers.

Inferring Small Neutron Star Spins with Neutron Star–Black Hole Mergers

The precise measurement of neutron star (NS) spins can provide important insight into the formation and evolution of compact binaries containing NSs. While traditional methods of NS spin measurement rely on pulsar observations, gravitational-wave detections offer a complementary avenue. However, determining component spins with gravitational waves is hindered by the small dimensionless spins of the NSs and the degeneracy in the mass and spin parameters. This degeneracy can be addressed by the inclusion of higher-order modes in the waveform, which are important for systems with unequal masses. This study shows the suitability of NS–black hole mergers, which are naturally mass-asymmetric, for precise NS spin measurements. We explore the effects of the black hole masses and spins, higher-mode content, inclination angles, and detector sensitivity on the measurement of NS spin. We find that networks with next-generation observatories like the Cosmic Explorer and the Einstein Telescope can distinguish NS dimensionless spin of 0.04 (0.1) from zero at 1σ confidence for events within ∼350 (∼1000) Mpc. Networks with A+ and A♯ detectors achieve similar distinction within ∼30 (∼70) Mpc and ∼50 (∼110) Mpc, respectively.

Forming off-center massive black hole binaries in dwarf galaxies through Jacobi capture

It is well established that black holes reside in the central regions of virtually all types of known galaxies. Recent observational and numerical studies however challenge this picture, suggesting that intermediate-mass black holes in dwarf galaxies may be found on orbits far from the center. In particular, constant-density cores minimize orbital energy losses due to dynamical friction, and allow black holes to settle on stable off-center orbits. Using controlled simulations, we study the dynamics of off-center black holes in dwarf galaxies with such cores. We propose a new scenario to describe off-center mergers of massive black holes, starting with a Jacobi capture. We focus on initially circular and co-planar black hole orbits and explore a large parameter space of black hole masses and orbital parameters. We find that Jacobi captures are a complex and chaotic phenomenon that occurs in about 13% of cases in this simplified setup, and we quantify how the likelihood of capture depends on the simulation parameters. We note that this percentage is likely an upper limit of the general case. Nevertheless, we show that Jacobi captures in cored dwarf galaxies can facilitate the formation of off-center black hole binaries, and that this process is sufficiently common to have a substantial effect on black hole growth. While our setup only allows for temporary captures, we expect dissipative forces from baryons and post-Newtonian corrections to maintain the captures over time and to lead to the formation of stable binary systems. This motivates future studies of the effectiveness of such dissipative forces, within stripped nuclei or globular clusters, in forming stable bound systems.

Viscous torque in turbulent magnetized active galactic nucleus accretion disks and its effects on the gravitational waves of extreme mass ratio inspirals

The merger of supermassive black holes produces millihertz gravitational waves (GWs), which are potentially detectable by the future Laser Interferometer Space Antenna (LISA). Such binary systems are usually embedded in an accretion disk environment at the center of an active galactic nucleus (AGN). Recent studies suggest the plasma environment imposes measurable imprints on the GW signal if the mass ratio of the binary is around q ∼ 10−4 − 10−3. The effect of the gaseous environment on the GW signal is strongly dependent on the disk’s parameters; therefore, it is believed that future low-frequency GW detections will provide us with precious information about the physics of AGN accretion disks. We investigated this effect by measuring the viscous torque via modeling of the evolution of magnetized tori around the primary massive black hole. Using the general relativistic magnetohydrodynamic HARM-COOL code, we performed 2D and 3D simulations of weakly magnetized, thin accretion disks, with a possible truncation and transition to advection-dominated accretion flow. We studied the angular momentum transport and turbulence generated by the magnetorotational instability. We quantified the disk’s effective alpha viscosity and its evolution over time. We applied our numerical results to quantify the relativistic viscous torque on a hypothetical low-mass secondary black hole via a 1D analytical approach, and we estimated the GW phase shift due to the gas environment.

From Seeds to Supermassive Black Holes: Capture, Growth, Migration, and Pairing in Dense Protobulge Environments

The origins and mergers of supermassive black holes (SMBHs) remain a mystery. We describe a scenario from a novel multiphysics simulation featuring rapid (≲1 Myr) hyper-Eddington gas capture by a ∼1000 M ⊙ “seed” black hole (BH) up to supermassive (≳106 M ⊙) masses in a massive, dense molecular cloud complex typical of high-redshift starbursts. Due to the high cloud density, stellar feedback is inefficient, and most of the gas turns into stars in star clusters that rapidly merge hierarchically, creating deep potential wells. Relatively low-mass BH seeds at random positions can be “captured” by merging subclusters and migrate to the center in ∼1 freefall time (vastly faster than dynamical friction). This also efficiently produces a paired BH binary with ∼0.1 pc separation. The centrally concentrated stellar density profile (akin to a “protobulge”) allows the cluster as a whole to capture and retain gas and build up a large (parsec-scale) circumbinary accretion disk with gas coherently funneled to the central BH (even when the BH radius of influence is small). The disk is “hypermagnetized” and “flux-frozen”: dominated by a toroidal magnetic field with plasma β ∼ 10−3, with the fields amplified by flux-freezing. This drives hyper-Eddington inflow rates ≳1 M ⊙ yr−1, which also drive the two BHs to nearly equal masses. The late-stage system appears remarkably similar to recently observed high-redshift “little red dots.” This scenario can provide an explanation for rapid SMBH formation, growth, and mergers in high-redshift galaxies.

Constraining Possible γ-Ray Burst Emission from GW230529 Using Swift-BAT and Fermi-GBM

GW230529 is the first compact binary coalescence detected by the LIGO–Virgo–KAGRA collaboration with at least one component mass confidently in the lower mass gap, corresponding to the range 3–5 M ⊙. If interpreted as a neutron star–black hole merger, this event has the most symmetric mass ratio detected so far and therefore has a relatively high probability of producing electromagnetic (EM) emission. However, no EM counterpart has been reported. At the merger time t 0, Swift-BAT and Fermi-GBM together covered 100% of the sky. Performing a targeted search in a time window [t 0 − 20 s, t 0 + 20 s], we report no detection by the Swift-BAT and Fermi-GBM instruments. Combining the position-dependent γ-ray flux upper limits and the gravitational-wave posterior distribution of luminosity distance, sky localization, and inclination angle of the binary, we derive constraints on the characteristic luminosity and structure of the jet possibly launched during the merger. Assuming a top-hat jet structure, we exclude at 90% credibility the presence of a jet that has at the same time an on-axis isotropic luminosity ≳1048 erg s−1 in the bolometric band 1 keV–10 MeV and a jet opening angle ≳15°. Similar constraints are derived by testing other assumptions about the jet structure profile. Excluding GRB 170817A, the luminosity upper limits derived here are below the luminosity of any GRB observed so far.

Discrete orbit effect lengthens merger times for inspiraling binary black holes

The inspiral merger time for two black holes captured into a nonrelativistic bound orbit by gravitational radiation emission has been often calculated by a formula of Peters that assumes the adiabatic approximation that the changes per orbit are small. However, initially this is not true for the semimajor axis and period of most of the initially highly eccentric orbits, which change significantly during closest approach and much less elsewhere along the orbit. This effect can make the merger time much longer (using other formulas from Peters that do not assume the adiabatic approximation) than that calculated by the adiabatic formula of Peters.

The comparison of an optical and X-ray counterpart of subparsec supermassive binary black holes

In this paper, we study and compare the optical and X-ray counterparts of subparsec supermassive black hole binaries (SMBHBs). With that aim, we simulated the profiles of optical spectral lines emitted from the broad line region (BLR) as well as X-ray spectral lines emitted from the relativistic accretion disks around both black holes and compared them with each other. The obtained results showed that SMBHBs could cause a specific, but different variability of the lines from the optical part and Fe Kα line, leaving potentially detectable imprints in their profiles. Since these imprints depend on the orbital phase of the system, they could be used for reconstructing the Keplerian orbits of the components in the observed SMBHBs. Moreover, such signatures in the optical and X-ray line profiles of the observed SMBHBs could be used as a tool for the detection of these objects as well as for studying their properties.

Consistent Eccentricities for Gravitational-wave Astronomy: Resolving Discrepancies between Astrophysical Simulations and Waveform Models

Detecting imprints of orbital eccentricity in gravitational-wave (GW) signals promises to shed light on the formation mechanisms of binary black holes. To constrain the formation mechanisms, distributions of eccentricity derived from numerical simulations of astrophysical formation channels are compared to the estimates of eccentricity inferred from GW signals. We report that the definition of eccentricity typically used in astrophysical simulations is inconsistent with the one used while modeling GW signals, with the differences mainly arising due to the choice of reference frequency used in both cases. We also posit a prescription for calculating eccentricity from astrophysical simulations, by evolving ordinary differential equations obtained from post-Newtonian theory and using the dominant (ℓ = m = 2) mode’s frequency as the reference frequency; this ensures consistency in the definitions. On comparing the existing eccentricities of the binaries present in the Cluster Monte Carlo catalog of globular cluster simulations with the eccentricities calculated using the prescription presented here, we find a significant discrepancy at e ≳ 0.2; this discrepancy becomes worse with increasing eccentricity. We note the implications this discrepancy has for existing studies and recommend that care be taken when comparing data-driven constraints on eccentricity to expectations from astrophysical formation channels.

Long term variability of Cygnus X-1

The Neutron Star Interior Composition Explorer (NICER) monitoring campaign of Cyg X-1 allows us to study its spectral-timing behavior at energies < 1 keV across all states. The hard state power spectrum can be decomposed into two main broad Lorentzians with a transition at around 1 Hz. The lower-frequency Lorentzian is the dominant component at low energies. The higher-frequency Lorentzian begins to contribute significantly to the variability above 1.5 keV and dominates at high energies. We show that the low- and high-frequency Lorentzians likely represent individual physical processes. The lower-frequency Lorentzian can be associated with a (possibly Comptonized) disk component, while the higher-frequency Lorentzian is clearly associated with the Comptonizing plasma. At the transition of these components, we discover a low-energy timing phenomenon characterized by an abrupt lag change of hard (≳2 keV) with respect to soft (≲1.5 keV) photons, accompanied by a drop in coherence, and a reduction in amplitude of the second broad Lorentzian. The frequency of the phenomenon increases with the frequencies of the Lorentzians as the source softens and cannot be seen when the power spectrum is single-humped. A comparison to transient low-mass X-ray binaries shows that this feature does not only appear in Cyg X-1, but that it is a general property of accreting black hole binaries. In Cyg X-1, we find that the variability at low and high energies is overall highly coherent in the hard and intermediate states. The high coherence shows that there is a process at work which links the variability, suggesting a physical connection between the accretion disk and Comptonizing plasma. This process fundamentally changes in the soft state, where strong red noise at high energies is incoherent to the variability at low energies.

Atacama Large Aperture Submillimeter Telescope (AtLAST) science: Probing the transient and time-variable sky

The study of transient and variable events, including novae, active galactic nuclei, and black hole binaries, has historically been a fruitful path for elucidating the evolutionary mechanisms of our universe. The study of such events in the millimeter and submillimeter is, however, still in its infancy. Submillimeter observations probe a variety of materials, such as optically thick dust, which are hard to study in other wavelengths. Submillimeter observations are sensitive to a number of emission mechanisms, from the aforementioned cold dust, to hot free-free emission, and synchrotron emission from energetic particles. Study of these phenomena has been hampered by a lack of prompt, high sensitivity submillimeter follow-up, as well as by a lack of high-sky-coverage submillimeter surveys. In this paper, we describe how the proposed Atacama Large Aperture Submillimeter Telescope (AtLAST) could fill in these gaps in our understanding of the transient universe. We discuss a number of science cases that would benefit from AtLAST observations, and detail how AtLAST is uniquely suited to contributing to them. In particular, AtLAST’s large field of view will enable serendipitous detections of transient events, while its anticipated ability to get on source quickly and observe simultaneously in multiple bands make it also ideally suited for transient follow-up. We make theoretical predictions for the instrumental and observatory properties required to significantly contribute to these science cases, and compare them to the projected AtLAST capabilities. Finally, we consider the unique ways in which transient science cases constrain the observational strategies of AtLAST, and make prescriptions for how AtLAST should observe in order to maximize its transient science output without impinging on other science cases.

Cosmological inference using gravitational waves and normalizing flows

We present a machine learning approach using normalizing flows for inferring cosmological parameters from gravitational wave events. Our methodology is general to any type of compact binary coalescence event and cosmological model and relies on the generation of training data representing distributions of gravitational wave event parameters. These parameters are conditional on the underlying cosmology and incorporate prior information from galaxy catalogues. We provide an example analysis inferring the Hubble constant using binary black holes detected during the O1, O2, and O3 observational runs conducted by the advanced LIGO/VIRGO gravitational wave detectors. We obtain a Bayesian posterior on the Hubble constant from which we derive an estimate and 1σ confidence bounds of H0=74.51−13.63+14.80 km s−1 Mpc−1. We are able to compute this result in O(1) s using our trained normalizing flow model. Published by the American Physical Society 2024

Reconstruction of Binary Black Hole Harmonics in LIGO Using Deep Learning

Gravitational-wave signals from coalescing compact binaries in the LIGO and Virgo interferometers are primarily detected by the template-based matched filtering method. While this method is optimal for stationary and Gaussian data scenarios, its sensitivity is often affected by nonstationary noise transients in the detectors. Moreover, most of the current searches do not account for the effects of precession of black hole spins and higher-order waveform harmonics, focusing solely on the leading-order quadrupolar modes. This limitation impacts our search for interesting astrophysical sources, such as intermediate-mass black hole binaries and hierarchical mergers. Here we show, for the first time, that deep learning can be used for accurate waveform reconstruction of precessing binary black hole signals with higher-order modes. This approach can also be adapted into a rapid trigger generation algorithm to enhance online searches. Our model, tested on simulated injections in real LIGO noise from the third observing run (2019–2020) achieved a high degree of overlap with injected signals. This accuracy was consistent across a wide range of black hole masses and spin configurations chosen for this study. When applied to real gravitational-wave events, our model's reconstructions achieved between 85% and 98% overlap with those obtained by Coherent WaveBurst (unmodeled) and LALInference (modeled) analyses. These results suggest that deep learning is a potent tool for analyzing signals from a diverse catalog of compact binaries.

Gravitational-wave matched filtering on a quantum computer

State-of-the-art quantum computers have very limited applicability for accurate calculations. Here, we report the first experimental demonstration of qubit-based matched filtering for a detection of the gravitational-wave signal from a binary black hole merger. With our implementation on noisy superconducting qubits, we obtained a similar signal-to-noise ratio for the binary black hole merger as achievable with classical computation, providing evidence for the utility of qubits for practically relevant tasks. The algorithm we invented for this application is a Monte Carlo algorithm which uses quantum and classical computation together. It provides a quasi-quadratic speed-up for time-domain convolution, similar to achievable with fast Fourier transform.

Discovering Love numbers through resonance excitation during extreme mass ratio inspirals

General Relativity predicts that black holes (BHs) do not possess an internal structure and consequently cannot be excited. This leads to a specific prediction about the waveform of gravitational waves (GWs) which they emit during a binary BH inspiral and to the vanishing of their Love numbers. However, if astrophysical BHs do possess an internal structure, their Love numbers would no longer vanish, and they could be excited during an inspiral by the transfer of orbital energy. This would affect the orbital period and lead to an observable imprint on the emitted GWs waveform. The effect is enhanced if one of the binary companions is resonantly excited. We discuss the conditions for resonant excitation of a hypothetical internal structure of BHs and calculate the phase change of the GWs waveform that is induced due to such resonant excitation during intermediate- and extreme-mass-ratio inspirals. We then relate the phase change to the electric quadrupolar Love number of the larger companion, which is resonantly excited by its smaller companion. We discuss the statistical error on measuring the Love number by LISA and show that, because of this phase change, the statistical error is small even for values of the Love number as small as 10−4 for moderate values of the spin parameter. Our results indicate that, for extreme-mass-ratio inspirals with moderate spin parameter, the Love number could be detected by LISA with an accuracy which is higher by up to two orders of magnitude than what can be achieved via tidal deformation effects. Thus, our results indicate that resonant excitation of the central BH during an extreme- or intermediate-mass-ratio inspirals is the most promising effect for putting bounds on, or detecting, non-vanishing tidal Love numbers of BHs.