A supermassive binary black hole with triple disks

Hierarchical structure formation inevitably leads to the formation of supermassive binary black holes (BBHs) with a sub-parsec separation in galactic nuclei. However, to date there has been no unambiguous detection of such systems. In an effort to search for potential observational signatures of supermassive BBHs, we performed high-resolution smoothed particle hydrodynamics (SPH) simulations of two black holes in a binary of moderate eccentricity surrounded by a circumbinary disk. Building on our previous work, which has shown that gas can periodically transfer from the circumbinary disk to the black holes when the binary is on an eccentric orbit, the current set of simulations focuses on the formation of the individual accretion disks, their evolution and mutual interaction, and the predicted radiative signature. The variation in mass transfer with orbital phase from the circumbinary disk induces periodic variations in the light curve of the two accretion disks at ultraviolet wavelengths, but not in the optical or near-infrared. Searches for this signal offer a promising method to detect supermassive BBHs.Comment: Accepted for publication in the Astrophysical Journal, 16 pages, 11 figures. High Resolution Version is Available at http://www2.yukawa.kyoto-u.ac.jp/~kimitake/bbhs.htm

HI Observations of the Supermassive Binary Black Hole System in 0402+379

We have recently discovered a supermassive binary black hole system with a projected separation between the two black holes of 7.3 parsecs in the radio galaxy 0402+379. This is the most compact supermassive binary black hole pair yet imaged by more than two orders of magnitude. We present Global VLBI observations at 1.3464 GHz of this radio galaxy, taken to improve the quality of the HI data. Two absorption lines are found toward the southern jet of the source, one redshifted by 370 +/- 10 km/s and the other blueshifted by 700 +/- 10 km/s with respect to the systemic velocity of the source, which, along with the results obtained for the opacity distribution over the source, suggests the presence of two mass clumps rotating around the central region of the source. We propose a model consisting of a geometrically thick disk, of which we only see a couple of clumps, that reproduces the velocities measured from the HI absorption profiles. These clumps rotate in circular Keplerian orbits around an axis that crosses one of the supermassive black holes of the binary system in 0402+379. We find an upper limit for the inclination angle of the twin jets of the source to the line of sight of 66 degrees, which, according to the proposed model, implies a lower limit on the central mass of ~7 x 10^8 Msun and a lower limit for the scale height of the thick disk of ~12 pc .Comment: 20 pages, 7 figures. Accepted on the Astrophysical Journa

Search of sub-parsec massive binary black holes through line diagnosis II

Massive black hole binaries at sub-parsec separations may display in their spectra anomalously small flux ratios between the MgII and CIV broad emission lines, i.e. F_MgII/F_CIV <~ 0.1, due to the erosion of the broad line region around the active, secondary black hole, by the tidal field of the primary. In Paper I by Montuori et al. (2011), we focussed on broad lines emitted by gas bound to the lighter accreting member of a binary when the binary is at the center of a hollow density region (the gap) inside a circum-binary disc. The main aim of this new study is at exploring the potential contribution to the broad line emission by the circum-binary disc and by gaseous streams flowing toward the black hole through the gap. We carry out a post-process analysis of data extracted from a SPH simulation of a circum-binary disc around a black hole binary. Our main result is that the MgII to CIV flux ratio can be reduced to ~ 0.1 within an interval of sub-pc binary separations of the order of a ~ (0.01-0.2)(f_Edd/0.1)^(1/2) pc corresponding to orbital periods of ~ (20-200) (f_Edd/0.1)^(3/4) years for a secondary BH mass in the range M_2 ~ 10^7-10^9 M_sun and a binary mass ratio of 0.3. At even closer separations this ratio returns to increase to values that are indistinguishable from the case of a single AGN (typically F_MgII/F_CIV ~ 0.3-0.4) because of the contribution to the MgII line from gas in the circum-binary disc.Comment: 7 pages, 3 figure, accepted for publication in MNRA

Accretion vs colliding wind models for the gamma-ray binary LS I +61 303: an assessment

LS I +61 303 is a puzzling Be/X-ray binary with variable gamma-ray emission at up TeV energies. The nature of the compact object and the origin of the high-energy emission are unclear. One family of models invokes particle acceleration in shocks from the collision between the B-star wind and a relativistic pulsar wind, while another centers on a relativistic jet powered by accretion. Recent high-resolution radio observations showing a putative "cometary tail" pointing away from the Be star near periastron have been cited as support for the pulsar-wind model. We wish here to carry out a quantitative assessment of these competing models for this extraordinary source. We apply a 3D SPH code for dynamical simulations of both the pulsar-wind-interaction and accretion-jet models. The former yields a description of the shape of the wind-wind interaction surface. The latter provides an estimation of the accretion rate. The results allow critical evaluation of how the two distinct models confront the data in various wavebands under a range of conditions. When one accounts for the 3D dynamical wind interaction under realistic constraints for the relative strength of the B-star and pulsar winds, the resulting form of the interaction front does not match the putative "cometary tail" claimed from radio observations. On the other hand, dynamical simulations of the accretion-jet model indicate that the orbital phase variation of accretion power includes a secondary broad peak well away from periastron, thus providing a plausible way to explain the observed TeV gamma ray emission toward apastron. We conclude that the colliding-wind model is not clearly established for LS I +61 303, while the accretion-jet model can reproduce many key characteristics of the observed TeV gamma-ray emission.Comment: Accepted for publication in A&A. The resolution of the figures is lower than in the journal paper to minimize file sizes. Seven pages, 5 figure

Dual black holes in merger remnants. I: linking accretion to dynamics

We study the orbital evolution and accretion history of massive black hole (MBH) pairs in rotationally supported circumnuclear discs up to the point where MBHs form binary systems. Our simulations have high resolution in mass and space which, for the first time, makes it feasible to follow the orbital decay of a MBH either counter- or co-rotating with respect to the circumnuclear disc. We show that a moving MBH on an initially counter-rotating orbit experiences an "orbital angular momentum flip" due to the gas-dynamical friction, i.e., it starts to corotate with the disc before a MBH binary forms. We stress that this effect can only be captured in very high resolution simulations. Given the extremely large number of gas particles used, the dynamical range is sufficiently large to resolve the Bondi-Hoyle-Lyttleton radii of individual MBHs. As a consequence, we are able to link the accretion processes to the orbital evolution of the MBH pairs. We predict that the accretion rate is significantly suppressed and extremely variable when the MBH is moving on a retrograde orbit. It is only after the orbital angular momentum flip has taken place that the secondary rapidly "lights up" at which point both MBHs can accrete near the Eddington rate for a few Myr. The separation of the double nucleus is expected to be around ~10 pc at this stage. We show that the accretion rate can be highly variable also when the MBH is co-rotating with the disc (albeit to a lesser extent) provided that its orbit is eccentric. Our results have significant consequences for the expected number of observable double AGNs at separations of <100 pc.Comment: 8 pages, 4 figures, revised version accepted to MNRA

Search of sub-parsec massive binary black holes through line diagnosis

We investigate on the spectral properties of an active black hole, member of a massive (10^7 - 10^9 Msun) sub-parsec black hole binary. We work under the hypothesis that the binary, surrounded by a circum-binary disc, has cleared a gap, and that accretion occurs onto the secondary black hole fed by material closer to the inner edge of the disc. Broad line emission clouds orbit around the active black hole and suffer erosion due to tidal truncation at the Roche Lobe surface, following gap opening and orbital decay. We consider three of the most prominent broad emission lines observed in the spectra of AGNs, i.e. CIV, MgII and H{\beta}, and compute the flux ratios between the lines of MgII and CIV (FMgII/FCIV) and those of MgII and H{\beta} (FMgII/FH{\beta}). We find that close black hole binaries have FMgII/FCIV up to one order of magnitude smaller than single black holes. By contrast FMgII/FH{\beta} may be significantly reduced only at the shortest separations. Peculiarly low values of line flux ratios together with large velocity offsets between the broad and narrow emission lines and/or periodic variability in the continuum (on timescales >= years) would identify genuine sub-pc binary candidates.Comment: 7 pages, 4 figures. Accepted for publication in MNRA

Properties of Accretion Flows Around Coalescing Supermassive Black Holes

What are the properties of accretion flows in the vicinity of coalescing supermassive black holes (SBHs)? The answer to this question has direct implications for the feasibility of coincident detections of electromagnetic (EM) and gravitational wave (GW) signals from coalescences. Such detections are considered to be the next observational grand challenge that will enable testing general relativity in the strong, nonlinear regime and improve our understanding of evolution and growth of these massive compact objects. In this paper we review the properties of the environment of coalescing binaries in the context of the circumbinary disk and hot, radiatively inefficient accretion flow models and use them to mark the extent of the parameter space spanned by this problem. We report the results from an ongoing, general relativistic, hydrodynamical study of the inspiral and merger of black holes, motivated by the latter scenario. We find that correlated EM+GW oscillations can arise during the inspiral phase followed by the gradual rise and subsequent drop-off in the light curve at the time of coalescence. While there are indications that the latter EM signature is a more robust one, a detection of either signal coincidentally with GWs would be a convincing evidence for an impending SBH binary coalescence. The observability of an EM counterpart in the hot accretion flow scenario depends on the details of a model. In the case of the most massive binaries observable by the Laser Interferometer Space Antenna, upper limits on luminosity imply that they may be identified by EM searches out to z~0.1-1. However, given the radiatively inefficient nature of the gas flow, we speculate that a majority of massive binaries may appear as low luminosity AGN in the local universe.Comment: Revised version accepted to Class. Quantum Grav. for proceedings of 8th LISA Symposium. 15 pages, 3 figures, includes changes suggested in referee report

SDSSJ092712.65+294344.0: a candidate massive black hole binary

In this Letter we explore the hypothesis that the quasar SDSSJ092712.65+294344.0 is hosting a massive black hole binary embedded in a circumbinary disc. The lightest, secondary black hole is active, and gas orbiting around it is responsible for the blue-shifted broad emission lines with velocity off-set of 2650 km/s, relative to the galaxy rest frame. As the tidal interaction of the binary with the outer disc is expected to excavate a gap, the blue-shifted narrow emission lines are consistent with being emitted from the low-density inhomogeneous gas of the hollow region. From the observations we infer a binary mass ratio q ~ 0.3, a mass for the primary of M1 ~ 2 billion Msun and a semi-major axis of 0.34 pc, corresponding to an orbital period of 370 years. We use the results of cosmological merger trees to estimate the likely-hood of observing SDSSJ092712.65+294344.0 as recoiling black hole or as a binary. We find that the binary hypothesis is preferred being one hundred times more probable than the ejection hypothesis. If SDSSJ092712.65+294344.0 hosts a binary, it would be the one closest massive black hole binary system ever discovered.Comment: Accepted for publication in MNRAS Letter

Massive binary black holes in galactic nuclei and their path to coalescence

Massive binary black holes form at the centre of galaxies that experience a merger episode. They are expected to coalesce into a larger black hole, following the emission of gravitational waves. Coalescing massive binary black holes are among the loudest sources of gravitational waves in the Universe, and the detection of these events is at the frontier of contemporary astrophysics. Understanding the black hole binary formation path and dynamics in galaxy mergers is therefore mandatory. A key question poses: during a merger, will the black holes descend over time on closer orbits, form a Keplerian binary and coalesce shortly after? Here we review progress on the fate of black holes in both major and minor mergers of galaxies, either gas-free or gas-rich, in smooth and clumpy circum-nuclear discs after a galactic merger, and in circum-binary discs present on the smallest scales inside the relic nucleus.Comment: Accepted for publication in Space Science Reviews. To appear in hard cover in the Space Sciences Series of ISSI "The Physics of Accretion onto Black Holes" (Springer Publisher

Extreme Air Pollution Events in Hokkaido, Japan Traced Back to Early Snowmelt and Large-Scale Wildfires Over East Eurasia: Case Studies

To identify the unusual climate conditions and their connections to air pollutions in a remote area due to wildfires, we examine three anomalous large-scale wildfires in May 2003, April 2008, and July 2014 over East Eurasia, as well as how products of those wildfires reached an urban city, Sapporo, in the northern part of Japan (Hokkaido), significantly affecting the air quality. NASA's MERRA-2 (the Modern-Era Retrospective analysis for Research and Applications, Version 2) aerosol re-analysis data closely reproduced the PM2.5 variations in Sapporo for the case of smoke arrival inJuly 2014. Results show that all three cases featured unusually early snowmelt in East Eurasia, accompanied by warmer and drier surface conditions in the months leading to the fires, inducing long-lasting soil dryness and producing environmental conditions conductive to active wildfires. Due to prevailing anomalous synoptic-scale atmospheric motions, smoke from those fires eventually reached a remote area, Hokkaido, and worsened the air quality in Sapporo. In futurestudies, continuous monitoring of the timing of Eurasian snowmelt and the air quality from the source regions to remote regions, coupled with the analysis of atmospheric and surface conditions, may be essential in more accurately predicting the effects of wildfires on air quality