No Arabic abstract
We study the evolution of a massive black hole pair in a rotationally supported nuclear disc. The distributions of stars and gas mimic the nuclear region of a gas-rich galaxy merger remnant. Using high-resolution SPH simulations, we follow the black hole dynamics and trace the evolution of the underlying background, until the black holes form a binary. We find that the gravitational perturbation of the pair creates a core in the disc density profile, hence decreasing the gas-dynamical drag. This leads the newly formed binary to stall at a separation of ~5 pc. In the early phases of the sinking, black holes lose memory of their initial orbital eccentricity if they co-rotate with the disc, as rotation of the gaseous background promotes circularization of the black hole orbits. Circularization is efficient until the black holes bind in a binary, though in the latest stages of the simulations a residual eccentricity > 0.1 is still present. Black holes are treated as sink particles, allowing for gas accretion. We find that accretion strongly depends on the dynamical properties of the black holes, and occurs preferentially after circularization.
If massive black holes (BHs) are ubiquitous in galaxies and galaxies experience multiple mergers during their cosmic assembly, then BH binaries should be common albeit temporary features of most galactic bulges. Observationally, the paucity of active BH pairs points toward binary lifetimes far shorter than the Hubble time, indicating rapid inspiral of the BHs down to the domain where gravitational waves lead to their coalescence. Here, we review a series of studies on the dynamics of massive BHs in gas-rich galaxy mergers that underscore the vital role played by a cool, gaseous component in promoting the rapid formation of the BH binary. The BH binary is found to reside at the center of a massive self-gravitating nuclear disc resulting from the collision of the two gaseous discs present in the mother galaxies. Hardening by gravitational torques against gas in this grand disc is found to continue down to sub-parsec scales. The eccentricity decreases with time to zero and when the binary is circular, accretion sets in around the two BHs. When this occurs, each BH is endowed with it own small-size (< 0.01 pc) accretion disc comprising a few percent of the BH mass. Double AGN activity is expected to occur on an estimated timescale of < 1 Myr. The double nuclear point-like sources that may appear have typical separation of < 10 pc, and are likely to be embedded in the still ongoing starburst. We note that a potential threat of binary stalling, in a gaseous environment, may come from radiation and/or mechanical energy injections by the BHs. Only short-lived or sub-Eddington accretion episodes can guarantee the persistence of a dense cool gas structure around the binary necessary for continuing BH inspiral.
We explore the hardening of a massive black hole binary embedded in a circum-binary gas disc when the binary and the gas are coplanar and the gas is counter-rotating. The secondary black hole, revolving in the direction opposite to the gas, experiences a drag from gas-dynamical friction and from direct accretion of part of it. Using two-dimensional (2D) hydrodynamical grid simulations we investigate the effect of changing the accretion prescriptions on the dynamics of the secondary black hole which in turn affect the binary hardening and eccentricity evolution. We find that realistic accretion prescriptions lead to results that differ from those inferred assuming accretion of all the gas within the Roche Lobe of the secondary black hole. Different accretion prescriptions result in different discs surface densities which alter the black holes dynamics back. Full 3D SPH realizations of a number of representative cases, run over a shorter interval of time, validate the general trends observed in the less computationally demanding 2D simulations. Initially circular black hole binaries increase only slightly their eccentricity which then oscillates around small values (<0.1) while they harden. By contrast, initially eccentric binaries become more and more eccentric. A semi-analytical model describing the black holes dynamics under accretion only explores the late evolution stages of the binary in an otherwise unperturbed retrograde disc to illustrate how eccentricity evolves with time in relation to the shape of the underlying surface density distribution.
We want to test if self-similar magneto-hydrodynamic (MHD) accretion-ejection models can explain the observational results for accretion disk winds in BHBs. In our models, the density at the base of the outflow, from the accretion disk, is not a free parameter, but is determined by solving the full set of dynamical MHD equations without neglecting any physical term. Different MHD solutions were generated for different values of (a) the disk aspect ratio ($varepsilon$) and (b) the ejection efficiency ($p$). We generated two kinds of MHD solutions depending on the absence (cold solution) or presence (warm solution) of heating at the disk surface. The cold MHD solutions are found to be inadequate to account for winds due to their low ejection efficiency. The warm solutions can have sufficiently high values of $p (gtrsim 0.1)$ which is required to explain the observed physical quantities in the wind. The heating (required at the disk surface for the warm solutions) could be due to the illumination which would be more efficient in the Soft state. We found that in the Hard state a range of ionisation parameter is thermodynamically unstable, which makes it impossible to have any wind at all, in the Hard state. Our results would suggest that a thermo-magnetic process is required to explain winds in BHBs.
The community may be on the verge of detecting low-frequency gravitational waves from massive black hole binaries (MBHBs), but no examples of binary active galactic nuclei (AGN) have been confirmed. Because MBHBs are intrinsically rare, the most promising detection methods utilize photometric data from all-sky surveys. Recently DOrazio & Di Stefano 2018 (arXiv:1707.02335) suggested gravitational self-lensing as a method of detecting AGN in close separation binaries. In this study we calculate the detectability of lensing signatures in realistic populations of simulated MBHBs. Within our model assumptions, we find that VROs LSST should be able to detect 10s to 100s of self-lensing binaries, with the rate uncertainty depending primarily on the orientation of AGN disks relative to their binary orbits. Roughly a quarter of lensing detectable systems should also show detectable Doppler boosting signatures. If AGN disks tend to be aligned with the orbit, lensing signatures are very nearly achromatic, while in misaligned configurations the bluer optical bands are lensed more than redder ones. Whether substantial obscuring material (e.g.~a dusty torus) will be present in close binaries remains uncertain, but our estimates suggest that a substantial fraction of systems would still be observable in this case.
Observations of accreting systems often show significant variability (10-20 percent of accretion luminosity) on timescales much longer than expected for the disc regions releasing most of the luminosity. We propose an explicit physical model for disc variability, consistent with Lyubarskiis (1997) general scheme for solving this problem. We suggest that local dynamo processes can affect the evolution of an accretion disc by driving angular momentum loss in the form of an outflow (a wind or jet). We model the dynamo as a small-scale stochastic phenomenon, operating on roughly the local dynamical timescale. We argue that large-scale outflow can only occur when the small-scale random processes in neighbouring disc annuli give rise by chance to a coherent large-scale magnetic field. This occurs on much longer timescales, and causes a bright large-amplitude flare as a wide range of disc radii evolve in a coherent fashion. Most of the time, dynamo action instead produces small-amplitude flickering. We reproduce power spectra similar to those observed, including a 1/f power spectrum below a break frequency given by the magnetic alignment timescale at the inner disc edge. However the relation between the black hole mass and the value of the break frequency is less straightforward than often assumed in the literature. The effect of an outer disc edge is to flatten the spectrum below the magnetic alignment frequency there. We also find a correlation between the variability amplitude and luminosity, similar to that found in some AGN.