Secondary Black Hole Research Articles

Being Careful with the Field Formation Interpretation of GW190412

Abstract The LIGO/Virgo Scientific Collaboration recently announced the detection of a compact object binary merger, GW190412, as the first asymmetric binary black hole (BBH) merger with mass ratio q ≈ 0.25. Other than the mass ratio, this BBH has been shown to have a positive effective spin of around χ eff ≈ 0.28. Assuming a field formation channel, associating this effective spin to either the primary or the secondary black hole (BH) each have implications: if the spin of the BBH comes from the primary BH, it has consequences for the efficiency of angular momentum transport in the formation of the BH. If, on the other hand, the spin is due to the secondary BH through tidal spin-up processes, it should be noted that (i) such processes have very short delay times, and (ii) subsequently, their local merger rate is determined by local star formation rate at the assumed formation metallicity of the BBH. We show that the predicted merger rate density from this channel is , and is in tension with the rather high local merger rate of such systems, which we estimate from this single event to be (90% confidence interval, and assuming 50 days of observing time). Large natal kicks (v ≳ 500 km s−1) would be required to get such an BBHs with an in-plane spin component to account for the marginal detection of precession in GW190412. However, this would only exacerbate the tension as the estimated local merger rate would be further decreased. Similarly, the formation of such systems through the dynamical assembly is exceedingly rare; consequently, it is difficult to account for this system with the currently accepted paradigms of BBH formation.

Spitzer Observations of the Predicted Eddington Flare from Blazar OJ 287

Abstract Binary black hole (BH) central engine description for the unique blazar OJ 287 predicted that the next secondary BH impact-induced bremsstrahlung flare should peak on 2019 July 31. This prediction was based on detailed general relativistic modeling of the secondary BH trajectory around the primary BH and its accretion disk. The expected flare was termed the Eddington flare to commemorate the centennial celebrations of now-famous solar eclipse observations to test general relativity by Sir Arthur Eddington. We analyze the multi-epoch Spitzer observations of the expected flare between 2019 July 31 and 2019 September 6, as well as baseline observations during 2019 February–March. Observed Spitzer flux density variations during the predicted outburst time display a strong similarity with the observed optical pericenter flare from OJ 287 during 2007 September. The predicted flare appears comparable to the 2007 flare after subtracting the expected higher base-level Spitzer flux densities at 3.55 and 4.49 μm compared to the optical R-band. Comparing the 2019 and 2007 outburst lightcurves and the previously calculated predictions, we find that the Eddington flare arrived within 4 hr of the predicted time. Our Spitzer observations are well consistent with the presence of a nano-Hertz gravitational-wave emitting spinning massive binary BH that inspirals along a general relativistic eccentric orbit in OJ 287. These multi-epoch Spitzer observations provide a parametric constraint on the celebrated BH no-hair theorem.

On Post-starburst Galaxies Dominating Tidal Disruption Events

Abstract A starburst induced by a galaxy merger may create a relatively thin central stellar disk at radius ≤100 pc. We calculate the rate of tidal disruption events (TDEs) by the inspiraling secondary supermassive black hole (SMBH) through the disk. With a small enough stellar velocity dispersion ( ) in the disk, it is shown that 105–106 TDEs of solar-type main-sequence stars per post-starburst galaxy (PSB) can be produced to explain their dominance in producing observed TDEs. Although the time it takes to bring the secondary SMBH to the disk apparently varies in the range of ∼0.1–1 Gyr since the starburst, depending on its landing location and subsequently due to dynamical friction with stars exterior to the central stellar disk in question, the vast majority of TDEs by the secondary SMBH in any individual PSB occur within a time span shorter than ∼30 Myr. Five unique testable predictions of this model are suggested.

Coalescence of Kerr Black Holes—Binary Systems from GW150914 to GW170814

We investigate the energy of the gravitational wave from a binary black hole merger by the coalescence of two Kerr black holes with an orbital angular momentum. The coalescence is constructed to be consistent with particle absorption in the limit in which the primary black hole is sufficiently large compared with the secondary black hole. In this limit, we analytically obtain an effective gravitational spin–orbit interaction dependent on the alignments of the angular momenta. Then, binary systems with various parameters including equal masses are numerically analyzed. According to the numerical analysis, the energy of the gravitational wave still depends on the effective interactions, as expected from the analytical form. In particular, we ensure that the final black hole obtains a large portion of its spin angular momentum from the orbital angular momentum of the initial binary black hole. To estimate the angular momentum released by the gravitational wave in the actual binary black hole, we apply our results to observations at the Laser Interferometer Gravitational-Wave Observatory: GW150914, GW151226, GW170104, GW170608 and GW170814.

Accretion Disk Parameters Determined from the Great 2015 Flare of OJ 287

Abstract In the binary black hole model of OJ 287, the secondary black hole orbits a much more massive primary, and impacts on the primary accretion disk at predictable times. We update the parameters of the disk, the viscosity, α, and the mass accretion rate, . We find α = 0.26 ± 0.1 and in Eddington units. The former value is consistent with Coroniti, and the latter with Marscher & Jorstad. Predictions are made for the 2019 July 30 superflare in OJ 287. We expect that it will take place simultaneously at the Spitzer infrared channels, as well as in the optical, and that therefore the timing of the flare in optical can be accurately determined from Spitzer observations. We also discuss in detail the light curve of the 2015 flare, and find that the radiating volume has regions where bremsstrahlung dominates, as well as regions that radiate primarily in synchrotron radiation. The former region produces the unpolarized first flare, while the latter region gives rise to a highly polarized second flare.

Tidal disruption event rates in galaxy merger remnants

ABSTRACT The rate of tidal disruption events (TDEs) depends sensitively on the stellar properties of the central galactic regions. Simulations show that galaxy mergers cause gas inflows, triggering nuclear starbursts, increasing the central stellar density. Motivated by these numerical results, and by the observed overrepresentation of post-starburst galaxies among TDE hosts, we study the evolution of the TDE rate in high-resolution hydrodynamical simulations of a galaxy merger, in which we capture the evolution of the stellar density around the massive black holes (BHs). We apply analytical estimates of the loss-cone theory, using the stellar density profiles from simulations, to estimate the time evolution of the TDE rate. At the second pericentre, a nuclear starburst enhances the stellar density around the BH in the least massive galaxy, leading to an enhancement of the TDE rate around the secondary BH, although the magnitude and the duration of the increase depend on the stochasticity of star formation on very small scales. The central stellar density around the primary BH remains instead fairly constant, and so is its TDE rate. After the formation of the binary, the stellar density decreases, and so does the TDE rate.

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.

Evidence for a TDE origin of the radio transient Cygnus A-2

In 2015, a radio transient named Cygnus A-2 was discovered in Cygnus A with the Very Large Array. Because of its radio brightness ($\nu F_{\nu} \approx 6 \times 10^{39}$ erg s$^{-1}$), this transient likely represents a secondary black hole in orbit around the AGN. Using {\it Chandra} ACIS observations from 2015 to 2017, we have looked for an X-ray counterpart to Cygnus A-2. The separation of 0.42 arcsec means that Cygnus A-2 can not be spatially resolved, but by comparing the data with simulated \texttt{marx} data, we put an upper limit to the 2-10 keV X-ray luminosity of Cygnus A-2 of $1 \times 10^{43}$ erg s$^{-1}$. Using the Fundamental Plane for accreting black holes, we find that our upper limit to the X-ray flux of Cygnus A-2 in 2015-2017 disfavours the interpretation of Cygnus A-2 as a steadily accreting black hole. We suggest instead that Cygnus A-2 is the radio afterglow of a tidal disruption event (TDE), and that a peak in the 2-10 keV luminosity of the nuclear region in 2013, when it was observed by {\it Swift} and {\it NuSTAR}, is X-ray emission from the TDE. A TDE could naturally explain the X-ray light curve of the nuclear region, as well as the appearance of a short-lived, fast, and ionized outflow previously detected in the 2013 {\it NuSTAR} spectrum. Both the radio and X-ray luminosities fall in between typical luminosities for 'thermal' and 'jetted' TDE types, suggesting that Cygnus A-2 would be unlike previously seen TDE's.

Investigating the Candidate Displaced Active Galactic Nucleus in NGC 3115

The nearby galaxy NGC 3115 contains a known radio-emitting, low-luminosity active galactic nucleus (AGN), and was recently claimed to host a candidate AGN displaced 14.3 pc from the galaxy's optical photocenter. Our goal is to understand whether this represents a single offset AGN, an AGN in orbit around a central black hole, or something else. We present a new, sensitive (RMS = 4.4 $\mu$Jy beam$^{-1}$) 10 GHz image, which finds evidence for only one AGN. We place a stringent limit on the radio luminosity of any secondary supermassive black hole of $L_{10~\rm{GHz}}<5.8\times10^{33}$ ergs/s. An analysis of the relative positioning of the radio core, X-ray nucleus, and stellar bulge in this galaxy indicate that the radio source is centrally located, and not offset from the galactic bulge. This provides an argument against a single offset AGN in NGC 3115, however does not provide conclusive evidence against the purported offset AGN as an in-spiralling secondary black hole.

Multi‐messenger observations of supermassive black holes binary mergers

When a circumprimary gaseous disk is pushed into the primary supermassive black hole (SMBH) by the tidal force of the decaying secondary black hole, it can produce an increase in the luminosity. In the case of a disk that lies in the same plane of the SMBHs binary, that luminosity exceeds the Eddington limit. In this paper, we extend these previous investigations to the case of systems misaligned with different mass values for the binary components. Our aim is to probe which mass range can produce stronger electromagnetic counterparts in the last stages of coalescence. We concentrate in geometrically thin disk models and binary SMBHs with mass ratio q = 10−3. We compare our numerical results with previously obtained analytical results in the literature. Disks misaligned surrounding a primary SMBH can lead to an electromagnetic signature of the gravitational wave signal emitted from final stages of the SMBH's binary orbital decay.

What is the final state of a black hole merger?

In this short paper, we discuss the possibility of testing the nature of astrophysical black holes using the recently observed black hole mergers. We investigate the possibility that a secondary black hole is created in the merger of two astrophysical black holes and discuss potential astrophysical signatures. We point out that black hole mergers are a possible astrophysical mechanism for the creation of quantum black holes with masses close to the Planck mass.

The complex jet- and bar-perturbed kinematics in NGC 3393 as revealed with ALMA and GEMINI–GMOS/IFU

NGC 3393, a nearby Seyfert 2 galaxy with nuclear radio jets, large-scale and nuclear bars, and a posited secondary super massive black hole, provides an interesting laboratory to test the physics of inflows and outflows. Here we present and analyse the molecular gas (ALMA observations of CO J:2-1 emission over a field of view (FOV) of 45\arcsec \times 45\arcsec, at 0\farcs56 (143 pc) spatial and 5 km/s spectral resolution), ionised gas and stars (GEMINI-GMOS/IFU; over a FOV of 4\arcsec \times 5\arcsec, at 0\farcs62 (159 pc) spatial and 23 km/s spectral resolution) in NGC 3393. The ionised gas emission, detected over the complete GEMINI-GMOS FOV, has three identifiable kinematic components. A narrow (\sigma < 115 km/s) component present in the complete FOV, which is consistent with rotation in the galaxy disk. A broad (\sigma > 115 km/s) redshifted component, detected near the NE and SW radio lobes; which we interpret as a radio jet driven outflow. And a broad (\sigma > 115 km/s) blueshifted component that shows high velocities in a region perpendicular to the radio jet axis; we interpret this as an equatorial outflow. The CO J:2-1 emission is detected in spiral arms on 5\arcsec - 20\arcsec scales, and in two disturbed circumnuclear regions. The molecular kinematics in the spiral arms can be explained by rotation. The highly disturbed kinematics of the inner region can be explained by perturbations induced by the nuclear bar and interactions with the large scale bar. We find no evidence for, but cannot strongly rule out, the presence of the posited secondary black hole.

Tidal disruption of stars in a supermassive black hole binary system: the influence of orbital properties on fallback and accretion rates

The disruption of a star by a supermassive black hole generates a sudden bright flare. Previous studies have focused on the disruption by single black holes, for which the fallback rate decays as~$\propto t^{-5/3}$. In this paper, we generalise the study to the case of a supermassive black hole binary (SMBHB), using both analytical estimates and hydrodynamical simulations, looking for specific observable signatures. The range of binary separation for which it is possible to distinguish between the disruption created by a single or a binary black hole concerns typically separations of order a few milliparsecs for a primary of mass $\sim 10^6M_{\odot}$. When the fallback rate is affected by the secondary, it undergoes two types interruptions, depending on the initial inclination $\theta$ of the orbit of the star relative to the plane of the SMBHB. For $\theta \lesssim 70^\circ$, periodic sharp interruptions occur and the time of first interruption depends on the distance of the secondary black hole with the debris. If $\theta \gtrsim 70^\circ$, a first smooth interruption occurs, but not always followed by a further recovery of the fallback rate. This implies that most of the TDEs around a SMBHB will undergo periodic sharp interruptions of their lightcurve.

Outflow-driven Transients from the Birth of Binary Black Holes. I. Tidally Locked Secondary Supernovae

Abstract We propose a new type of electromagnetic transient associated with the birth of binary black holes (BBHs), which may lead to merger events accompanied by gravitational waves in . We consider the newborn BBHs formed through the evolution of isolated massive stellar binaries. For a close massive binary, consisting of a primary black hole (BH) and a secondary Wolf–Rayet (WR) star that are orbiting around each other, the spin period of the WR star can be tidally synchronized to its orbital period. Then the angular momentum of the outer material of the WR star is large enough to form an accretion disk around a newborn, secondary BH, following its core-collapse. This disk may produce an energetic outflow with a kinetic energy of and an outflow velocity of , resulting in an optical transient with an absolute magnitude from approximately −14 to approximatley −17 with a duration of around a day. This type of transient also produces detectable radio signals years after the birth of BBHs, via synchrotron emission from nonthermal electrons accelerated at external shocks. The predicted optical transients have a shorter duration than ordinary core-collapse supernovae. Dedicated optical transient surveys can detect them and distinguish them from ordinary SNe using the different features of its light curve and late-time spectrum. In this paper (Paper I), we investigate disk-driven outflows from the secondary BH, whereas possible signals from the primary BH will be discussed in Paper II.

Probabilistic eccentricity bifurcation for stars around shrinking massive black hole binaries

Based on the secular theory, we discuss the orbital evolution of stars in a nuclear star cluster to which a secondary massive black hole is infalling with vanishing eccentricity. We find that the eccentricities of the stars could show sharp transitions, depending strongly on their initial conditions. By examining the phase-space structure of an associated Hamiltonian, we show that these characteristic behaviours are partly due to a probabilistic bifurcation at a separatrix crossing, resulting from the retrograde apsidal precession by the cluster potential. We also show that separatrix crossings are closely related to realization of a large eccentricity and could be important for astrophysical phenomena such as tidal disruption events or gravitational wave emissions.

VLBA 24 and 43 GHz observations of massive binary black hole candidate PKS 1155 + 251

PKS 1155+251 is a radio-loud quasar source at z=0.203. Observations using very long baseline interferometry (VLBI) at ~2, 5, 8 and 15 GHz show that the structure of the radio source is quite complicated on parsec scales and that the outer hot spots are apparently undergoing a significant contraction. Because these results cannot be fully explained based on the compact symmetric object (CSO) scenario with a radio core located between the northern and southern complexes, we made observations with the Very Long Baseline Array (VLBA) at 24 and 43 GHz to search for compact substructures and alternative interpretations. The results show that the radio core revealed in the previous VLBI observations remains compact with a flat spectrum in our sub-milli-arcsecond--resolution images; the northern lobe emission becomes faint at 24 GHz and is mostly resolving out at 43 GHz; the southern complex is more bright but has been resolved into the brightest southern-end (S1) and jet or tail alike components westwards. Explaining the southern components aligned westward with a standard CSO scenario alone remains a challenge. As for the flatter spectral index of the southern-end component S1 between 24 and 43 GHz in our observations and the significant 15 GHz VLBA flux variability of S1, an alternative scenario is that the southern complex may be powered by a secondary black hole residing at S1. But more sensitive and high-resolution VLBI monitoring is required to discriminate the CSO and the binary black hole scenarios.

Discovery of a Luminous Radio Transient 460 pc from the Central Supermassive Black Hole in Cygnus A

Abstract We report the appearance of a new radio source at a projected offset of 460 pc from the nucleus of Cygnus A. The flux density of the source (which we designate Cygnus A-2) rose from an upper limit of &lt;0.5 mJy in 1989 to 4 mJy in 2016 (ν = 8.5 GHz), but is currently not varying by more than a few percent per year. The radio luminosity of the source is comparable to the most luminous known supernovae, it is compact in Very Long Baseline Array observations down to a scale of 4 pc, and it is coincident with a near-infrared point source seen in pre-existing adaptive optics and HST observations. The most likely interpretation of this source is that it represents a secondary supermassive black hole in a close orbit around the Cygnus A primary, though an exotic supernova model cannot be ruled out. The gravitational influence of a secondary SMBH at this location may have played an important role in triggering the rapid accretion that has powered the Cygnus A radio jet over the past 107 years.

A swirling jet in the quasar 1308+326

Context. Despite numerous and detailed studies of the jets of active galactic nuclei (AGN) on pc-scales, many questions are still debated. The physical nature of the jet components is one of the most prominent unsolved problems, as is the launching mechanism of jets in AGN. The quasar 1308+326 (z = 0.997) allows us to study the overall properties of its jet in detail and to derive a more physical understanding of the nature and origin of jets in general. The long-term data provided by the Monitoring Of Jets in Active galactic nuclei with Very Long Baseline Array (VLBA) experiments (MOJAVE) survey permit us to trace out the structural changes in 1308+326 that we present here. The long-lived jet features in this source can be followed for about two decades.Aims. We investigate the very long baseline interferomety (VLBI) morphology and kinematics of the jet of 1308+326 to understand the physical nature of this jet and jets in general, the role of magnetic fields, and the causal connection between jet features and the launching process.Methods. Fifty VLBA observations performed at 15 GHz from the MOJAVE survey were re-modeled with Gaussian components and re-analyzed (the time covered: 20 Jan. 1995–25 Jan. 2014). The analysis was supplemented by multi-wavelength radio-data (UMRAO, at 4.8, 8.0, and 14.5 GHz) in polarization and total intensity. We fit the apparent motion of the jet features with the help of a model of a precessing nozzle.Results. The jet features seem to be emitted with varying viewing angles and launched into an ejection cone. Tracing the component paths yields evidence for rotational motion. Radio flux-density variability can be explained as a consequence of enhanced Doppler boosting corresponding to the motion of the jet relative to the line of sight. Based on the presented kinematics and other indicators, such as electric-vector polarization position-angle (EVPA) rotation, we conclude that the jet of 1308+326 has a helical structure, meaning that the components are moving along helical trajectories and the trajectories themselves are also experiencing a precessing motion. A model of a precessing nozzle was applied to the data and a subset of the observed jet feature paths can be modeled successfully within this model. The data till 2012 are consistent with a swing period of 16.9 yr. We discuss several scenarios to explain the observed motion phenomena, including a binary black hole model. It seems unlikely that the accretion disk around the primary black hole, which is disturbed by the tidal forces of the secondary black hole, is able to launch a persistent axisymmetric jet.Conclusions. We conclude that we are observing a rotating helix. In particular, the observed EVPA swings can be explained by a shock moving through a straight jet that is pervaded by a helical magnetic field. We compare our results for 1308+326 with other astrophysical scenarios where similar, wound-up filamentary structures are found. They are all related to accretion-driven processes. A helically moving or wound up object is often explained by filamentary features moving along magnetic field lines of magnetic flux tubes. It seems that a “component” comprises plasma tracing the magnetic field, which guides the motion of the radiating radio-band plasma. Further investigations and modeling are in preparation.

The Pairing of Accreting Massive Black Holes in Multiphase Circumnuclear Disks: the Interplay Between Radiative Cooling, Star Formation, and Feedback Processes

Abstract We study the orbital decay of a pair of massive black holes (BHs) with masses and 107 , using hydrodynamical simulations of circumnuclear disks (CNDs) with the alternating presence of sub-grid physics, such as radiative cooling, star formation, supernova feedback, BH accretion, and BH feedback. In the absence of such processes, the orbit of the secondary BH decays over timescales of to the center of the CND, where the primary BH resides. When strong dissipation operates in CNDs, fragmentation into massive objects the size of giant molecular clouds with densities in the range 104–107 amu cm−3 occurs, causing stochastic torques and hits that can eject the secondary BH from the midplane. Outside the plane, the low-density medium provides only weak drag, and the BH return is governed by inefficient dynamical friction. In rare cases, clump–BH interactions can lead to a faster decay. Feedback processes lead to outflows, but do not significantly change the overall density of the CND midplane. However, with a spherically distributed BH feedback, a hot bubble is generated behind the secondary, which almost shuts off dynamical friction. We dub this phenomenon “wake evacuation.” It leads to delays in the decay, possibly of . We discuss the non-trivial implications on the discovery space of the eLISA telescope. Our results suggest that the largest uncertainty in predicting BH merger rates lies in the potentially wide variety of galaxy host systems, with different degrees of gas dissipation and heating, yielding decay timescales from to .

Supermassive black hole pairs in clumpy galaxies at high redshift: delayed binary formation and concurrent mass growth

Massive gas-rich galaxy discs at $z \sim 1-3$ host massive star-forming clumps with typical baryonic masses in the range $10^7-10^8$~M$_{\odot}$ which can affect the orbital decay and concurrent growth of supermassive black hole (BH) pairs. Using a set of high-resolution simulations of isolated clumpy galaxies hosting a pair of unequal-mass BHs, we study the interaction between massive clumps and a BH pair at kpc scales, during the early phase of the orbital decay. We find that both the interaction with massive clumps and the heating of the cold gas layer of the disc by BH feedback tend to delay significantly the orbital decay of the secondary, which in many cases is ejected and then hovers for a whole Gyr around a separation of 1--2 kpc. In the envelope, dynamical friction is weak and there is no contribution of disc torques: these lead to the fastest decay once the orbit of the secondary BH has circularised in the disc midplane. In runs with larger eccentricities the delay is stronger, although there are some exceptions. We also show that, even in discs with very sporadic transient clump formation, a strong spiral pattern affects the decay time-scale for BHs on eccentric orbits. We conclude that, contrary to previous belief, a gas-rich background is not necessarily conducive to a fast BH decay and binary formation, which prompts more extensive investigations aimed at calibrating event-rate forecasts for ongoing and future gravitational-wave searches, such as with Pulsar Timing Arrays and the future evolved Laser Interferometer Space Antenna.