skip to main content
Explore Scholarly Publications and Datasets in the NSF-PAR

This content will become publicly available on April 1, 2023

Title: Electromagnetic Signatures from Supermassive Binary Black Holes Approaching Merger
Abstract We present fully relativistic predictions for the electromagnetic emission produced by accretion disks surrounding spinning and nonspinning supermassive binary black holes on the verge of merging. We use the code Bothros to post-process data from 3D general relativistic magnetohydrodynamic simulations via ray-tracing calculations. These simulations model the dynamics of a circumbinary disk and the mini-disks that form around two equal-mass black holes orbiting each other at an initial separation of 20 gravitational radii, and evolve the system for more than 10 orbits in the inspiral regime. We model the emission as the sum of thermal blackbody radiation emitted by an optically thick accretion disk and a power-law spectrum extending to hard X-rays emitted by a hot optically thin corona. We generate time-dependent spectra, images, and light curves at various frequencies to investigate intrinsic periodic signals in the emission, as well as the effects of the black hole spin. We find that prograde black hole spin makes mini-disks brighter since the smaller innermost stable circular orbit angular momentum demands more dissipation before matter plunges to the horizon. However, compared to mini-disks in larger separation binaries with spinning black holes, our mini-disks are less luminous: unlike those systems, their mass accretion more » rate is lower than in the circumbinary disk, and they radiate with lower efficiency because their inflow times are shorter. Compared to a single black hole system matched in mass and accretion rate, these binaries have spectra noticeably weaker and softer in the UV. Finally, we discuss the implications of our findings for the potential observability of these systems. « less
Authors:
; ; ; ; ; ;
Award ID(s):
2110338 2018420 2031744 1811228 2009330 1912632 1707946 1707826
Publication Date:
NSF-PAR ID:
10331565
Journal Name:
The Astrophysical Journal
Volume:
928
Issue:
2
Page Range or eLocation-ID:
137
ISSN:
0004-637X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ( , The Astrophysical Journal)
    Abstract We perform a full 3D general relativistic magnetohydrodynamical (GRMHD) simulation of an equal-mass, spinning, binary black hole approaching merger, surrounded by a circumbinary disk and with a minidisk around each black hole. For this purpose, we evolve the ideal GRMHD equations on top of an approximated spacetime for the binary that is valid in every position of space, including the black hole horizons, during the inspiral regime. We use relaxed initial data for the circumbinary disk from a previous long-term simulation, where the accretion is dominated by a m = 1 overdensity called the lump. We compare our newmore »spinning simulation with a previous non-spinning run, studying how spin influences the minidisk properties. We analyze the accretion from the inner edge of the lump to the black hole, focusing on the angular momentum budget of the fluid around the minidisks. We find that minidisks in the spinning case have more mass over a cycle than the non-spinning case. However, in both cases we find that most of the mass received by the black holes is delivered by the direct plunging of material from the lump. We also analyze the morphology and variability of the electromagnetic fluxes, and we find they share the same periodicities of the accretion rate. In the spinning case, we find that the outflows are stronger than the non-spinning case. Our results will be useful to understand and produce realistic synthetic light curves and spectra, which can be used in future observations.« less
  2. ; ; ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT Luminous active galactic nuclei and X-ray binaries often contain geometrically thin, radiatively cooled accretion discs. According to theory, these are – in many cases – initially highly misaligned with the black hole equator. In this work, we present the first general relativistic magnetohydrodynamic simulations of very thin (h/r ∼ 0.015–0.05) accretion discs around rapidly spinning (a ∼ 0.9) black holes and tilted by 45°–65°. We show that the inner regions of the discs with h/r ≲ 0.03 align with the black hole equator, though out to smaller radii than predicted by analytic work. The inner aligned and outer misalignedmore »disc regions are separated by a sharp break in tilt angle accompanied by a sharp drop in density. We find that frame dragging by the spinning black hole overpowers the disc viscosity, which is self-consistently produced by magnetized turbulence, tearing the disc apart and forming a rapidly precessing inner sub-disc surrounded by a slowly precessing outer sub-disc. We find that the system produces a pair of relativistic jets for all initial tilt values. At small distances, the black hole launched jets precess rapidly together with the inner sub-disc, whereas at large distances they partially align with the outer sub-disc and precess more slowly. If the tearing radius can be modeled accurately in future work, emission model independent measurements of black hole spin based on precession-driven quasi-periodic oscillations may become possible.« less
  3. ; ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    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,more »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 of 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
  4. ; ; ; ; ; ( , The Astrophysical Journal)
    Abstract For testing different electron temperature ( T e ) prescriptions in general relativistic magnetohydrodynamics (GRMHD) simulations through observations, we propose to utilize linear polarization (LP) and circular polarization (CP) images. We calculate the polarization images based on a semi-magnetically arrested disk GRMHD model for various T e parameters, bearing M87 in mind. We find an LP–CP separation in the images of the low- T e disk cases at 230GHz; namely, the LP flux mainly originates from downstream of the jet, and the CP flux comes from the counter-side jet, while the total intensity is maximum at the jet base.more »This can be understood as follows: although the LP flux is generated through synchrotron emission widely around the black hole, most of the LP flux from the jet base does not reach the observer, since it undergoes Faraday rotation ( ∝ T e − 2 ) when passing through the outer cold disk and is thus depolarized. Hence, only the LP flux from the downstream (not passing the cold dense plasmas) can survive. Meanwhile, the CP flux is generated from the LP flux by Faraday conversion ( ∝ T e ) in the inner hot region. Stronger CP flux is thus observed from the counter-side jet. Moreover, the LP–CP separation is more enhanced at a lower frequency, such as 86 GHz, but is rather weak at 43 GHz, since the media in the latter case is optically thick for synchrotron self-absorption so that all of the fluxes should come from the photosphere. The same is true for cases with higher mass accretion rates and/or larger inclination angles.« less
  5. ; ; ; ; ; ; ; ; ; ; et al ( , Monthly Notices of the Royal Astronomical Society)
    Abstract Binary supermassive black holes (BSBHs) are expected to be a generic byproduct from hierarchical galaxy formation. The final coalescence of BSBHs is thought to be the loudest gravitational wave (GW) siren, yet no confirmed BSBH is known in the GW-dominated regime. While periodic quasars have been proposed as BSBH candidates, the physical origin of the periodicity has been largely uncertain. Here we report discovery of a periodicity (P=1607±7 days) at 99.95% significance (with a global p-value of ∼10−3 accounting for the look elsewhere effect) in the optical light curves of a redshift 1.53 quasar, SDSS J025214.67−002813.7. Combining archival Sloan Digital Skymore »Survey data with new, sensitive imaging from the Dark Energy Survey, the total ∼20-yr time baseline spans ∼4.6 cycles of the observed 4.4-yr (restframe 1.7-yr) periodicity. The light curves are best fit by a bursty model predicted by hydrodynamic simulations of circumbinary accretion disks. The periodicity is likely caused by accretion rate modulation by a milli-parsec BSBH emitting GWs, dynamically coupled to the circumbinary accretion disk. A bursty hydrodynamic variability model is statistically preferred over a smooth, sinusoidal model expected from relativistic Doppler boost, a kinematic effect proposed for PG1302−102. Furthermore, the frequency dependence of the variability amplitudes disfavors Doppler boost, lending independent support to the circumbinary accretion variability hypothesis. Given our detection rate of one BSBH candidate from circumbinary accretion variability out of 625 quasars, it suggests that future large, sensitive synoptic surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time may be able to detect hundreds to thousands of candidate BSBHs from circumbinary accretion with direct implications for Laser Interferometer Space Antenna.« less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email