skip to main content

Title: An estimate of the stochastic gravitational wave background from the MassiveBlackII simulation

A population of supermassive black hole (SMBH) binaries is expected to generate a stochastic gravitational wave background (SGWB) in the pulsar timing array (PTA) frequency range of 10−9 to $10^{-7}\, {\rm Hz}$. Detection of this signal is a current observational goal and so predictions of its characteristics are of significant interest. In this work, we use SMBH binary mergers from the MassiveBlackII simulation to estimate the characteristic strain of the stochastic background. We examine both a gravitational wave (GW) driven model of binary evolution and a model which also includes the effects of stellar scattering and a circumbinary gas disc. Results are consistent with PTA upper limits and similar to estimates in the literature. The characteristic strain at a reference frequency of $1\, {\rm yr}^{-1}$ is found to be $A_{\rm {yr}^{-1}}= 6.9 \times 10^{-16}$ and $A_{\rm {yr}^{-1}}= 6.4 \times 10^{-16}$ in the GW-driven and stellar scattering/gas disc cases, respectively. Using the latter approach, our models show that the SGWB is mildly suppressed compared to the purely GW-driven model as frequency decreases inside the PTA frequency band.

; ; ; ; ;
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 5241-5250
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT Massive black hole (MBH) binary inspiral time-scales are uncertain, and their spins are even more poorly constrained. Spin misalignment introduces asymmetry in the gravitational radiation, which imparts a recoil kick to the merged MBH. Understanding how MBH binary spins evolve is crucial for determining their recoil velocities, their gravitational wave (GW) waveforms detectable with Laser Interferometer Space Antenna, and their retention rate in galaxies. Here, we introduce a sub-resolution model for gas- and gravitational wave (GW)-driven MBH binary spin evolution using accreting MBHs from the Illustris cosmological hydrodynamic simulations. We also model binary inspiral via dynamical friction, stellar scattering, viscous gas drag, and GW emission. Our model assumes that the circumbinary disc always removes angular momentum from the binary. It also assumes differential accretion, which causes greater alignment of the secondary MBH spin in unequal-mass mergers. We find that 47 per cent of the MBHs in our population merge by z = 0. Of these, 19 per cent have misaligned primaries and 10 per cent have misaligned secondaries at the time of merger in our fiducial model with initial eccentricity of 0.6 and accretion rates from Illustris. The MBH misalignment fraction depends strongly on the accretion disc parameters, however. Reducing accretion rates by a factor ofmore »100, in a thicker disc, yields 79 and 42 per cent misalignment for primaries and secondaries, respectively. Even in the more conservative fiducial model, more than 12 per cent of binaries experience recoils of >500 km s−1, which could displace them at least temporarily from galactic nuclei. We additionally find that a significant number of systems experience strong precession.« less

    Rotating neutron stars (NSs) are promising sources of gravitational waves (GWs) in the frequency band of ground-based detectors. They are expected to emit quasi-monochromatic, long-duration GW signals, called continuous waves (CWs), due to their deviations from spherical symmetry. The degree of such deformations, and hence the information about the internal structure of an NS, is encoded in a dimension-less parameter ε called ellipticity. Searches for CW signals from isolated Galactic NSs have shown to be sensitive to ellipticities as low as $\varepsilon \sim \mathcal {O}(10^{-9})$. These searches are optimal for detecting and characterizing GWs from individual NSs, but they are not designed to measure the properties of NSs as population, such as the average ellipticity εav. These ensemble properties can be determined by the measurement of the stochastic gravitational-wave background (SGWB) arising from the superposition of GW signals from individually undetectable NSs. In this work, we perform a cross-correlation search for such a SGWB using the data from the first three observation runs of Advanced LIGO and Virgo. Finding no evidence for an SGWB signal, we set upper limits on the dimension-less energy density parameter Ωgw(f). Using these results, we also constrain the average ellipticity of Galactic NSs andmore »five NS ‘hotspots’, as a function of the number of NSs emitting GWs within the frequency band of the search Nband. We find $\varepsilon _{\mathrm{av}} \lesssim 1.8 \times 10^{-8}$, with Nband = 1.6 × 107, for Galactic NSs, and $\varepsilon _{\mathrm{av}} \lesssim [3.5-11.8]\times 10^{-7}$, with Nband = 1.6 × 1010, for NS hotspots.

    « less
  3. ABSTRACT We present X-ray and radio observations of what may be the closest Type Iax supernova (SN) to date, SN 2014dt (d = 12.3–19.3 Mpc), and provide tight constraints on the radio and X-ray emission. We infer a specific radio luminosity $L_R\lt (1.0\!-\!2.4)\times 10^{25}\, \rm {erg\, s^{-1}\, Hz^{-1}}$ at a frequency of 7.5 GHz and a X-ray luminosity $L_X\lt 1.4\times 10^{38}\, \rm {erg\, s^{-1}}$ (0.3–10 keV) at ∼38–48 d post-explosion. We interpret these limits in the context of Inverse Compton (IC) emission and synchrotron emission from a population of electrons accelerated at the forward shock of the explosion in a power-law distribution $N_e(\gamma _e)\propto \gamma _e^{-p}$ with p = 3. Our analysis constrains the progenitor system mass-loss rate to be $\dot{M}\lt 5.0 \times 10^{-6} \rm {M_{\odot }\, yr^{-1}}$ at distances $r\lesssim 10^{16}\, \rm {cm}$ for an assumed wind velocity $v_w=100\, \rm {km\, s^{-1}}$, and a fraction of post-shock energy into magnetic fields and relativistic electrons of ϵB = 0.01 and ϵe = 0.1, respectively. This result rules out some of the parameter space of symbiotic giant star companions, and it is consistent with the low mass-loss rates expected from He-star companions. Our calculations also show that the improved sensitivity of the next-generation Very Largemore »Array (ngVLA) is needed to probe the very low-density media characteristic of He stars that are the leading model for binary stellar companions of white dwarfs giving origin to Type Iax SNe.« less
  4. Abstract

    The nanohertz gravitational wave background (GWB) is believed to be dominated by GW emission from supermassive black hole binaries (SMBHBs). Observations of several dual-active galactic nuclei (AGN) strongly suggest a link between AGN and SMBHBs, given that these dual-AGN systems will eventually form bound binary pairs. Here we develop an exploratory SMBHB population model based on empirically constrained quasar populations, allowing us to decompose the GWB amplitude into an underlying distribution of SMBH masses, SMBHB number density, and volume enclosing the GWB. Our approach also allows us to self-consistently predict the number of local SMBHB systems from the GWB amplitude. Interestingly, we find the local number density of SMBHBs implied by the common-process signal in the NANOGrav 12.5-yr data set to be roughly five times larger than previously predicted by other models. We also find that at most ∼25% of SMBHBs can be associated with quasars. Furthermore, our quasar-based approach predicts ≳95% of the GWB signal comes fromz≲ 2.5, and that SMBHBs contributing to the GWB have masses ≳108M. We also explore how different empirical galaxy–black hole scaling relations affect the local number density of GW sources, and find that relations predicting more massive black holes decrease the localmore »number density of SMBHBs. Overall, our results point to the important role that a measurement of the GWB will play in directly constraining the cosmic population of SMBHBs, as well as their connections to quasars and galaxy mergers.

    « less
  5. ABSTRACT Active galactic nuclei (AGN) are powered by the accretion of discs of gas on to supermassive black holes (SMBHs). Stars and stellar remnants orbiting the SMBH in the nuclear star cluster (NSC) will interact with the AGN disc. Orbiters plunging through the disc experience a drag force and, through repeated passage, can have their orbits captured by the disc. A population of embedded objects in AGN discs may be a significant source of binary black hole mergers, supernovae, tidal disruption events, and embedded gamma-ray bursts. For two representative AGN disc models, we use geometric drag and Bondi–Hoyle–Littleton drag to determine the time to capture for stars and stellar remnants. We assume a range of initial inclination angles and semimajor axes for circular Keplerian prograde orbiters. Capture time strongly depends on the density and aspect ratio of the chosen disc model, the relative velocity of the stellar object with respect to the disc, and the AGN lifetime. We expect that for an AGN disc density $\rho \gtrsim 10^{-11}{\rm g\, cm^{-3}}$ and disc lifetime ≥1 Myr, there is a significant population of embedded stellar objects, which can fuel mergers detectable in gravitational waves with LIGO-Virgo and LISA.