No Arabic abstract
In this paper, we explore the mechanisms that regulate the formation and evolution of stellar black hole binaries (BHBs) around supermassive black holes (SMBHs). We show that dynamical interactions can efficiently drive in-situ BHB formation if the SMBH is surrounded by a massive nuclear cluster (NC), while orbitally segregated star clusters can replenish the BHB reservoir in SMBH-dominated nuclei. We discuss how the combined action of stellar hardening and mass segregation sculpts the BHB orbital properties. We use direct N-body simulations including post-Newtonian corrections up to 2.5 order to study the BHB-SMBH interplay, showing that the Kozai-Lidov mechanism plays a crucial role in shortening binaries lifetime. We find that the merging probability weakly depends on the SMBH mass in the $10^6-10^9{rm ~M}_odot$ mass range, leading to a merger rate $Gamma simeq 3-8$ yr$^{-1}$ Gpc$^{-3}$ at redshift zero. Nearly $40%$ of the mergers have masses in the BH mass gap, $50-140{rm ~M}_odot$, thus indicating that galactic nuclei are ideal places to form BHs in this mass range. We argue that gravitational wave (GW) sources with components mass $m_1>40{rm ~M}_odot$ and $m_2<30{rm ~M}_odot$ would represent a strong indicator of a galactic nuclei origin. The majority of these mergers could be multiband GW sources in the local Universe: nearly $40%$ might be seen by LISA as eccentric sources and, a few years later, as circular sources by LIGO and the Einstein Telescope, making decihertz observatories like DECIGO unique instruments to bridge the observations during the binary inspiral.
Galaxy mergers produce supermassive black hole binaries, which emit gravitational waves prior to their coalescence. We perform three-dimensional hydrodynamic simulations to study the tidal disruption of stars by such a binary in the final centuries of its life. We find that the gas stream of the stellar debris moves chaotically in the binary potential and forms accretion disks around both black holes. The accretion light curve is modulated over the binary orbital period owing to relativistic beaming. This periodic signal allows to detect the decay of the binary orbit due to gravitational wave emission by observing two tidal disruption events that are separated by more than a decade.
We study the collapse of rapidly rotating supermassive stars that may have formed in the early Universe. By self-consistently simulating the dynamics from the onset of collapse using three-dimensional general-relativistic hydrodynamics with fully dynamical spacetime evolution, we show that seed perturbations in the progenitor can lead to the formation of a system of two high-spin supermassive black holes, which inspiral and merge under the emission of powerful gravitational radiation that could be observed at redshifts z>10 with the DECIGO or Big Bang Observer gravitational-wave observatories, assuming supermassive stars in the mass range 10^4-10^6 Msol. The remnant is rapidly spinning with dimensionless spin a^*=0.9. The surrounding accretion disk contains ~10% of the initial mass.
We investigate the evolution of supermassive binary black holes (BBHs) in galaxies with realistic property distributions and the gravitational-wave (GW) radiation from the cosmic population of these BBHs. We incorporate a comprehensive treatment of the dynamical interactions of the BBHs with their environments by including the effects of galaxy triaxial shapes and inner stellar distributions, and generate a large number of BBH evolution tracks. By combining these BBH evolution tracks, galaxy mass functions, galaxy merger rates, and supermassive black hole-host galaxy relations into our model, we obtain the statistical distributions of surviving BBHs, BBH coalescence rates, the strength of their GW radiation, and the stochastic GW background (GWB) contributed by the cosmic BBH population. About ~1%-3% (or ~10%) of supermassive BHs at nearby galactic centers are expected to be binaries with mass ratio >1/3 (or >1/100). The characteristic strain amplitude of the GWB at frequency 1/yr is estimated to be ~$2.0^{+1.4}_{-0.8}times 10^{-16}$, and the upper bound of its results obtained with the different BH-host galaxy relations can be up to $5.4times 10^{-16}$, which await testing by future experiments (e.g., the Square Kilometer Array, FAST, Next-Generation Very Large Array). The turnover frequency of the GWB spectrum is at ~0.25nHz. The uncertainties on the above estimates and prospects for detecting individual sources are also discussed. The application of the cosmic BBH population to the Laser Interferometer Space Antenna (LISA) band provides a lower limit to the detection rate of BBHs by LISA, ~0.9/yr.
We introduce a technique for gravitational-wave analysis, where Gaussian process regression is used to emulate the strain spectrum of a stochastic background using population-synthesis simulations. This leads to direct Bayesian inference on astrophysical parameters. For PTAs specifically, we interpolate over the parameter space of supermassive black-hole binary environments, including 3-body stellar scattering, and evolving orbital eccentricity. We illustrate our approach on mock data, and assess the prospects for inference with data similar to the NANOGrav 9-yr data release.
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project currently observes 43 pulsars using the Green Bank and Arecibo radio telescopes. In this work we use a subset of 17 pulsars timed for a span of roughly five years (2005--2010). We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Within the timing data, we perform a search for continuous gravitational waves from individual supermassive black hole binaries in circular orbits using robust frequentist and Bayesian techniques. We find that there is no evidence for the presence of a detectable continuous gravitational wave; however, we can use these data to place the most constraining upper limits to date on the strength of such gravitational waves. Using the full 17 pulsar dataset we place a 95% upper limit on the sky-averaged strain amplitude of $h_0lesssim 3.8times 10^{-14}$ at a frequency of 10 nHz. Furthermore, we place 95% emph{all sky} lower limits on the luminosity distance to such gravitational wave sources finding that the $d_L gtrsim 425$ Mpc for sources at a frequency of 10 nHz and chirp mass $10^{10}{rm M}_{odot}$. We find that for gravitational wave sources near our best timed pulsars in the sky, the sensitivity of the pulsar timing array is increased by a factor of $sim$4 over the sky-averaged sensitivity. Finally we place limits on the coalescence rate of the most massive supermassive black hole binaries.