Massive Black Hole Binaries Research Articles

LISA parameter estimation and source localization with higher harmonics of the ringdown

LISA can detect higher harmonics of the ringdown gravitational-wave signal from massive black-hole binary mergers with large signal-to-noise ratio. The most massive black-hole binaries are more likely to have electromagnetic counterparts, and the inspiral will contribute little to their signal-to-noise ratio. Here we address the following question: can we extract the binary parameters and localize the source using LISA observations of the ringdown only? Modulations of the amplitude and phase due to LISA's motion around the Sun can be used to disentangle the source location and orientation when we detect the long-lived inspiral signal, but they can not be used for ringdown-dominated signals, which are very short-lived. We show that (i) we can still measure the mass ratio and inclination of high-mass binaries by carefully combining multiple ringdown harmonics, and (ii) we can constrain the sky location and luminosity distance by relying on the relative amplitudes and phases of various harmonics, as measured in different LISA channels.

The Most Massive Binary Black Hole Detections and the Identification of Population Outliers

Abstract Advanced LIGO and Virgo detected 10 binary black holes (BBHs) in their first two observing runs (O1 and O2). Analysis of these events found evidence for a dearth of BBHs with component masses greater than ∼45 M ⊙, as would be expected from pair-instability supernovae. Meanwhile, a standalone analysis of the merger GW170729 found its primary mass to be , which appears to be in contradiction with the existence of a limit at ∼45 M ⊙. In this work, we argue that the masses of individual events can only be evaluated with reference to the full population. When GW170729 is analyzed jointly with the remaining detections, its inferred primary mass tightens considerably, to . In the presence of noise, apparent outliers in the detected distribution are inevitable. We discuss methods of distinguishing between statistical fluctuations and true population outliers using posterior predictive tests. Applying these tests to O1 and O2, we find that the 10 detections are consistent with even the simplest power-law plus maximum-mass model considered by the LIGO–Virgo Collaboration, supporting the claim that GW170729 is not a population outlier. We also provide non-parametric constraints on the rate of high-mass mergers and conservatively bound the rate of mergers with m 1 > 45 M ⊙ at of the total merger rate. After 100 detections like those of O1 and O2 from a population with a maximum primary mass of 45 M ⊙, it would be common for the most massive system to have an observed maximum-likelihood mass m 1 ≳ 70 M ⊙.

Dynamic signatures of black hole binaries with superradiant clouds

Superradiant clouds may develop around a rotating black hole, if there is a bosonic field with Compton wavelength comparable to the size of the black hole. In this paper, we investigate the effects of the cloud on the orbits of nearby compact objects. In particular, we consider the dynamical friction and the backreaction due to level mixing. Under these interactions, the probability of a black hole dynamically capturing other compact objects, such as stellar mass black holes and neutron stars, is generally enhanced with the presence of the cloud. For extreme mass ratio inspirals and binary stellar mass binary black holes, the cloud-induced orbital modulation may be detected by observing the gravitational waveform using space borne gravitational wave detectors, such as LISA. Interestingly within certain range of boson Compton wavelength, the enhanced capture rate of stellar mass black holes could accelerate hierarchical mergers, with higher-generation merger product being more massive than the mass threshold predicted by supernova pair instability. These observational signatures provide promising ways of searching light bosons with gravitational waves.

On the eccentricity evolution of massive black hole binaries in stellar backgrounds

ABSTRACT We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three-body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of q = m2/m1 = 10−4, where m1 and m2 are the masses of the primary and secondary hole, respectively. Contrary to a simple extrapolation from higher mass ratios, we find that for q ≲ 10−3 the eccentricity growth rate becomes negative, i.e. the binary circularizes as it shrinks. This behaviour is due to the subset of interacting stars captured in metastable counter-rotating orbits; those stars tend to extract angular momentum from the binary, promoting eccentricity growth for q > 10−3, but tend to inject angular momentum into the binary driving it towards circularization for q < 10−3. The physical origin of this behaviour requires a detailed study of the orbits of this subset of stars and is currently under investigation. Our findings might have important consequences for intermediate massive black holes (IMBHs) inspiralling on to MBHs (e.g. a putative $10^3\,\rm M_{\odot }$ black hole inspiralling on to Sgr A*).

No tension between pulsar timing array upper limits on the nano-Hertz gravitational wave background and assembly models of massive black hole binaries

Pulsar timing arrays provide a means to observe the nano-Hertz gravitational wave background from the population of merging massive black hole binaries. Observations are placing increasingly stringent upper limits on the gravitational wave background. Upper limits and future detections will enable the study of the properties of the merging population. Recent upper limits have cast doubt on current predictions of the gravitational wave background. Here we perform a Bayesian analysis comparing upper limits to astrophysical prediction. So far models are consistent with observation. These proceedings summarise previous work in Ref. [1].

Formation of binary black holes in star clusters as gravitational wave sources

Direct detections of gravitational wave suggest that ~ 30 solar mass binary black holes (BBHs) commonly exist in the Universe. One possible formation scenario of such BBHs is dynamical three-body encounters in a dense core of globular clusters, which consist of millions of stars. Compared to globular clusters, open clusters are less dense and less massive but more populous. Because of their shallow potential, they have not been considered to be a formation site of merging BBHs. However, we found a new channel for the formation of BBHs, which is dominant in open clusters. We performed direct N-body simulations of open clusters with a mass of 2500 solar mass for metallicity of Z/Z ⊙ = 0.1. The core-collapse time of open clusters is much shorter than that of globular clusters. Therefore, massive main-sequence stars can form binaries before they evolve to BHs. These main-sequence binaries experience mass transfer evolution, and some of them evolve to BBHs merging within the Hubble time. From our simulations, we estimated the merger rates of BBHs originated from open clusters.

NGC 6240: A triple nucleus system in the advanced or final state of merging

Aims. NGC 6240 is a well-studied nearby galaxy system in the process of merging. Based on optical, X-ray, and radio observations, it is thought to harbor two active nuclei. We carried out a detailed optical 3D spectroscopic study to investigate the inner region of this system in connection with existing MERLIN and VLBA data. Methods. We observed NGC 6240 with very high spatial resolution using the MUSE instrument in the Narrow-Field Mode with the four-laser GALACSI adaptive optics system on the ESO VLT under seeing conditions of 0″.49. Our 3D spectra cover the wavelength range from 4725 to 9350 Å at a spatial resolution of ∼75 mas. Results. We report the discovery of three nuclei in the final state of merging within a region of only 1 kpc in the NGC 6240 system. Thanks to MUSE we are able to show that the formerly unresolved southern component actually consists of two distinct nuclei separated by only 198 pc. In combination with Gaia data we reach an absolute positional accuracy of only 30 mas that is essential to compare optical spectra with MERLIN and VLBA radio positions. Conclusions. The verification and detailed study of a system with three nuclei, two of which are active and each with a mass in excess of 9 × 107 M⊙, is of great importance for the understanding of hierarchical galaxy formation via merging processes since multiple mergers lead to a faster evolution of massive galaxies in comparison to binary mergers. So far it has been suggested that the formation of galactic nuclei with multiple supermassive black holes (SMBHs) is expected to be rare in the local universe. Triple massive black hole systems might be of fundamental importance for the coalescence of massive black hole binaries in less than a Hubble time leading to the loudest sources of gravitational waves in the millihertz regime.

OUP accepted manuscript

Gravitational wave (GW) radiations from numerous cosmic stellar-compact-binaries form a stochastic GW background (GWB), which is expected to be detected by ground and space GW detectors in future. Theoretical predictions of this GWB were mostly obtained by assuming either circular binaries and/or a specific channel for binary formation, which may have some uncertainties. In this paper, we estimate the GWB and its spectrum by using simple models for the formation of both stellar mass binary black holes (sBBHs) and binary neutron stars (BNSs). We consider that the dynamically originated sBBHs have relatively larger masses and higher eccentricities compared with those from field binary stars and its possible effect on the GWB spectrum. We find that the GWB spectrum may have a turnover in the low-frequency (Laser Interferometer Space Antenna; LISA) band and may be better described by a broken double power-law than a single power-law with the unique index $2/3$, and the low-frequency slope depends on the significance of the dynamically originated sBBHs with high eccentricities. We further generate mock samples of sBBHs and BNSs, and simulate the mock GWB strain in the time domain. We find that GWB can be detected with signal-to-noise ratio (SNR) $\gtrsim274/255/21$ by LISA/Taiji/TianQin over $5$-years' observation and $\gtrsim3$ by LIGO over $2$-years' observation. Furthermore, we estimate that the number of sBBHs that may be detected by LISA/Taiji/TianQin individually with SNR $\gtrsim8$ is $\sim5$-$221$/$7$-$365$/$3$-$223$ over $5$-years' observation.

On the properties of the massive binary black hole merger GW170729

We present a detailed investigation into the properties of GW170729, the gravitational wave with the most massive and distant source confirmed to date. We employ an extensive set of waveform models, including new improved models that incorporate the effect of higher-order waveform modes which are particularly important for massive systems. We find no indication of spin-precession, but the inclusion of higher-order modes in the models results in an improved estimate for the mass ratio of $(0.3-0.8)$ at the 90\% credible level. Our updated measurement excludes equal masses at that level. We also find that models with higher-order modes lead to the data being more consistent with a smaller effective spin, with the probability that the effective spin is greater than zero being reduced from $99\%$ to $94\%$. The 90\% credible interval for the effective spin parameter is now $(-0.01-0.50)$. Additionally, the recovered signal-to-noise ratio increases by $\sim0.3$ units compared to analyses without higher-order modes. We study the effect of common spin priors on the derived spin and mass measurements, and observe small shifts in the spins, while the masses remain unaffected. We argue that our conclusions are robust against systematic errors in the waveform models. We also compare the above waveform-based analysis which employs compact-binary waveform models to a more flexible wavelet- and chirplet-based analysis. We find consistency between the two, with overlaps of $\sim 0.9$, typical of what is expected from simulations of signals similar to GW170729, confirming that the data are well-described by the existing waveform models. Finally, we study the possibility that the primary component of GW170729 was the remnant of a past merger of two black holes and find this scenario to be indistinguishable from the standard formation scenario.

Probing Massive Black Hole Binary Populations with LISA

Abstract ESA and NASA are moving forward with plans to launch LISA around 2034. With data from the Illustris cosmological simulation, we provide analysis of LISA detection rates accompanied by characterization of the merging massive black hole population. Massive black holes of total mass ∼105 − 1010M⊙ are the focus of this study. We evolve Illustris massive black hole mergers, which form at separations on the order of the simulation resolution (∼kpc scales), through coalescence with two different treatments for the binary massive black hole evolutionary process. The coalescence times of the population, as well as physical properties of the black holes, form a statistical basis for each evolutionary treatment. From these bases, we Monte Carlo synthesize many realizations of the merging massive black hole population to build mock LISA detection catalogs. We analyze how our massive black hole binary evolutionary models affect detection rates and the associated parameter distributions measured by LISA. With our models, we find massive black hole binary detection rates with LISA of ∼0.5 − 1 yr−1 for massive black holes with masses greater than 105M⊙. This should be treated as a lower limit primarily because our massive black hole sample does not include masses below 105M⊙, which may significantly add to the observed rate. We suggest reasons why we predict lower detection rates compared to much of the literature.

Search for intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network

Gravitational wave astronomy has been firmly established with the detection of gravitational waves from the merger of ten stellar mass binary black holes and a neutron star binary. This paper reports on the all-sky search for gravitational waves from intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network. The search uses three independent algorithms: two based on matched filtering of the data with waveform templates of gravitational wave signals from compact binaries, and a third, model-independent algorithm that employs no signal model for the incoming signal. No intermediate mass black hole binary event was detected in this search. Consequently, we place upper limits on the merger rate density for a family of intermediate mass black hole binaries. In particular, we choose sources with total masses $M=m_1+m_2\in[120,800]$M$_\odot$ and mass ratios $q = m_2/m_1 \in[0.1,1.0]$. For the first time, this calculation is done using numerical relativity waveforms (which include higher modes) as models of the real emitted signal. We place a most stringent upper limit of $0.20$~Gpc$^{-3}$yr$^{-1}$ (in co-moving units at the 90% confidence level) for equal-mass binaries with individual masses $m_{1,2}=100$M$_\odot$ and dimensionless spins $\chi_{1,2}= 0.8$ aligned with the orbital angular momentum of the binary. This improves by a factor of $\sim 5$ that reported after Advanced LIGO's first observing run.

The Impact of Pair-instability Mass Loss on the Binary Black Hole Mass Distribution

A population of binary black hole mergers has now been observed in gravitational waves by Advanced LIGO and Virgo. The masses of these black holes appear to show evidence for a pile-up between $30$--$45$ $M_\odot$ and a cut-off above $\sim 45$ $M_\odot$. One possible explanation for such a pile-up and subsequent cut-off are pulsational pair-instability supernovae (PPISNe) and pair-instability supernovae (PISNe) in massive stars. We investigate the plausibility of this explanation in the context of isolated massive binaries. We study a population of massive binaries using the rapid population synthesis software COMPAS, incorporating models for PPISNe and PISNe. Our models predict a maximum black hole mass of $40$ $M_\odot$. We expect $\sim 10$\% of all binary black hole mergers at redshift z = 0 will include at least one component that went through a PPISN (with mass $30$--$40$ $M_\odot$), constituting $\sim 20$--$50$\% of binary black hole mergers observed during the first two observing runs of Advanced LIGO and Virgo. Empirical models based on fitting the gravitational-wave mass measurements to a combination of a power law and a Gaussian find a fraction too large to be associated with PPISNe in our models. The rates of PPISNe and PISNe track the low metallicity star formation rate, increasing out to redshift $z = 2$. These predictions may be tested both with future gravitational-wave observations and with observations of superluminous supernovae.

Binary Black Hole Population Properties Inferred from the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

Abstract We present results on the mass, spin, and redshift distributions with phenomenological population models using the 10 binary black hole (BBH) mergers detected in the first and second observing runs completed by Advanced LIGO and Advanced Virgo. We constrain properties of the BBH mass spectrum using models with a range of parameterizations of the BBH mass and spin distributions. We find that the mass distribution of the more massive BH in such binaries is well approximated by models with no more than 1% of BHs more massive than 45 M ⊙ and a power-law index of α = 1.3 − 1.7 + 1.4 (90% credibility). We also show that BBHs are unlikely to be composed of BHs with large spins aligned to the orbital angular momentum. Modeling the evolution of the BBH merger rate with redshift, we show that it is flat or increasing with redshift with 93% probability. Marginalizing over uncertainties in the BBH population, we find robust estimates of the BBH merger rate density of R = 53.2 − 28.2 + 55.8 Gpc−3 yr−1 (90% credibility). As the BBH catalog grows in future observing runs, we expect that uncertainties in the population model parameters will shrink, potentially providing insights into the formation of BHs via supernovae, binary interactions of massive stars, stellar cluster dynamics, and the formation history of BHs across cosmic time.

Science with the TianQin observatory: Preliminary results on testing the no-hair theorem with ringdown signals

We study the capability of the space-based gravitational wave observatory TianQin to test the no-hair theorem of General Relativity, using the ringdown signal from the coalescence of massive black hole binaries. We parameterize the ringdown signal by the four strongest quasinormal modes and estimate the signal to noise ratio for various source parameters. We consider constraints both from single detections and from all the events combined throughout the lifetime of the observatory, for different astrophysical models. We find that at the end of the mission, TianQin will have constrained deviations of the frequency and decay time of the dominant 22 mode from the general relativistic predictions to within 0.2% and 1.5% respectively, the frequencies of the subleading modes can be also constrained within 0.3%. We also find that TianQin and LISA are highly complementary, by virtue of their different frequency windows. Indeed, LISA can best perform ringdown tests for black hole masses in excess of $\sim 3\times 10^6 M_\odot$, while TianQin is best suited for lower masses.

Compact Object Binary Mergers Driven By Cluster Tides: A New Channel for LIGO/Virgo Gravitational-wave Events

Abstract The detections of gravitational waves (GWs) produced in mergers of binary black holes (BHs) and neutron stars (NSs) by LIGO/Virgo have stimulated interest in the origin of the progenitor binaries. Dense stellar systems—globular and nuclear star clusters—are natural sites of compact object binary formation and evolution toward merger. Here we explore a new channel for the production of binary mergers in clusters, in which the tidal field of the cluster secularly drives the binary to high eccentricity (even in the absence of a central massive BH) until GW emission becomes important. We employ the recently developed secular theory of cluster tide-driven binary evolution to compute present day merger rates for BH–BH, NS–BH, and NS–NS binaries, varying cluster potential and central concentration of the binary population (but ignoring cluster evolution and stellar flybys for now). Unlike other mechanisms, this new dynamical channel can produce a significant number of mergers out to cluster-centric distances of several parsecs. For NS–NS binaries we find merger rates in the range of 0.01–0.07 Gpc−3 yr−1 from globular clusters and 0.1–0.2 Gpc−3 yr−1 from cusped nuclear clusters. For NS–BH and BH–BH binaries we find small merger rates from globular clusters, but a rate of 0.1–0.2 Gpc−3 yr−1 from cusped nuclear clusters, contributing to the observed LIGO/Virgo rate at the level of several percent. Therefore, cluster tide-driven mergers constitute a new channel that can be further explored with current and future GW detectors.

Science with the TianQin observatory: Preliminary results on massive black hole binaries

We investigate the prospects of detecting gravitational waves from coalescing massive black hole binaries in the Universe with the TianQin observatory, a space-based gravitational wave interferometer proposed to be launched in the 2030s. To frame the scientific scope of the mission, in this paper, we carry out a preliminary estimation of the signal-to-noise ratio, detection rate and parameter estimation precision of massive black hole binaries detectable by TianQin. In order to make our results as robust as possible, we consider several models of the growth history of massive black holes, exploring the effect of some key astrophysical prescriptions as well the impact of the employed computational methods. In the most optimistic model, TianQin can detect as many as approximately 60 mergers per year. If TianQin detects a merger at redshift of 15, it will be capable of estimating its luminosity distance to within an accuracy of 10%; for a nearby event at redshift approximately 2, TianQin can issue early warnings 24 hours before coalescence, with a timing accuracy of around three hours and a sky localization ability of approximately 80 deg$^2$, thus enabling multi-messenger observations.

Blip glitches in Advanced LIGO data

Blip glitches are short noise transients present in data from ground-based gravitational-wave observatories. These glitches resemble the gravitational-wave signature of massive binary black hole mergers. Hence, the sensitivity of transient gravitational-wave searches to such high-mass systems and other potential short duration sources is degraded by the presence of blip glitches. The origin and rate of occurrence of this type of glitch have been largely unknown. In this paper we explore the population of blip glitches in Advanced LIGO during its first and second observing runs. On average, we find that Advanced LIGO data contains approximately two blip glitches per hour of data. We identify four subsets of blip glitches correlated with detector auxiliary or environmental sensor channels, however the physical causes of the majority of blips remain unclear.

Tidal disruption events from massive black hole binaries: predictions for ongoing and future surveys

ABSTRACTWe compute the expected cosmic rates of tidal disruption events (TDEs) induced by individual massive black holes (MBHs) and by MBH binaries (MBHBs) – with a specific focus on the latter class – to explore the potential of TDEs to probe the cosmic population of sub-pc MBHBs. Rates are computed by combining MBH and MBHB population models derived from large cosmological simulations with estimates of the induced TDE rates for each class of objects. We construct empirical TDE spectra that fit a large number of observations in the optical, UV, and X-ray and consider their observability by current and future survey instruments. Consistent with results in the literature, and depending on the detailed assumption of the model, we find that LSST and Gaia in optical and eROSITA in X-ray will observe a total of 3000–6000, 80–180, and 600–900 TDEs per year, respectively. Depending on the survey, 1 to several per cent of these are prompted by MBHBs. In particular, both LSST and eROSITA are expected to see 150–450 MBHB-induced TDEs in their respective mission lifetimes, including 5–100 repeated flares. The latter provide an observational sample of binary candidates with relatively low contamination and have the potential of unveiling the sub-pc population of MBHBs in the mass range $10^5\lt M\lt 10^7\, \mathrm{M}_\odot$, thus informing future low-frequency gravitational wave observatories.

The Unique Blazar OJ 287 and Its Massive Binary Black Hole Central Engine

The bright blazar OJ 287 is the best-known candidate for hosting a nanohertz gravitational wave (GW) emitting supermassive binary black hole (SMBBH) in the present observable universe. The binary black hole (BBH) central engine model, proposed by Lehto and Valtonen in 1996, was influenced by the two distinct periodicities inferred from the optical light curve of OJ 287. The current improved model employs an accurate general relativistic description to track the trajectory of the secondary black hole (BH) which is crucial to predict the inherent impact flares of OJ 287. The successful observations of three predicted impact flares open up the possibility of using this BBH system to test general relativity in a hitherto unexplored strong field regime. Additionally, we briefly describe an ongoing effort to interpret observations of OJ 287 in a Bayesian framework.

Black Hole dynamics in Young Star Clusters

AbstractYoung star clusters are a promising environment for forming binary black holes. Such binaries may form dynamically or via binary star evolution or through the interplay of these two channels. To study these formation pathways, we have performed high precision direct N-body simulations of low-mass (M < 1000 M⊙) young star clusters. The simulations were carried out with the code Nbody6++GPU coupled with the population synthesis code MOBSE. Our results highlight the importance of dynamics to form massive black hole binaries even in low-mass young star clusters.