It arises a puzzle in \NGC\, how to secularly maintain the counter-rotating disc from $0.2$ to $7\,$pc unambiguously detected by recent ALMA observations of molecular gas. Upon further analysis of disc dynamics, we find that the Kelvin-Helmholtz (KH) instability (KHI) results in an unavoidable catastrophe of the disc developed at the interface between the reversely rotating parts, and demonstrate that a close binary of supermassive black holes provides tidal torques as the unique external sources to prevent the disc from the KH catastrophe. We are led to the inescapable conclusion that there must be a binary black hole at the center of NGC 1068, to prevent it from the KH catastrophe. The binary is composed of black holes with a separation of $0.1\,$pc from GRAVITY/VLTI observations, a total mass of $1.3\times 10^{7}\:M_{\odot}$ and a mass ratio of $\sim 0.3$ estimated from the angular momentum budge of the global system. The KHI gives rise to forming a gap without cold gas at the velocity interface which overlaps with the observed gap of hot and cold dust regions. Releases of kinematic energies from the KHI of the disc are in agreement with observed emissions in radio and $\gamma$-rays. Such a binary is shrinking with a timescale much longer than the local Hubble time via gravitational waves, however, the KHI leads to an efficient annihilation of the orbital angular momentum and speed up merge of the binary, providing a new paradigm of solving the long term issue of "final parsec problem". Future observations of GRAVITY+/VLTI are expected to be able to spatially resolve the CB-SMBHs suggested in this paper.

Distance measurements for extragalactic objects are a fundamental problem in astronomy and cosmology. In the era of precision cosmology, we urgently need better measurements of cosmological distances to observationally test the increasing $H_{0}$ tension of the Hubble constant measured from different tools. Using spectroastrometry, GRAVITY at The Very Large Telescope Interferometer successfully revealed the structure, kinematics and angular sizes of the broad-line region (BLR) of 3C 273 with an unprecedentedly high spatial resolution. Fortunately, reverberation mapping (RM) of active galactic nuclei (AGNs) reliably provides linear sizes of their BLRs. Here we report a joint analysis of spectroastrometry and RM observations to measure AGN distances. We apply this analysis to 3C 273 observed by both GRAVITY and an RM campaign,and find an angular distance of $551.5_{-78.7}^{+97.3}\, {\rm Mpc}$ and $H_{0}=71.5_{-10.6}^{+11.9}\,{\rm km\,s^{-1}\,Mpc^{-1}}$. Advantages of the analysis are 1) its pure geometrical measurements and 2) it simultaneously yields mass of the central black hole in the BLR. Moreover, we can conveniently repeat measurements of selected AGNs to efficiently reduce the statistical and systematic errors. Future observations of a reasonably sized sample ($\sim 30$ AGNs) will provide distances of the AGNs and hence a new way of measuring $H_{0}$ with a high precision $\left(\lesssim 3\%\right)$ to test the $H_{0}$ tension.

Most large galaxies harbor supermassive black holes (SMBHs) in their centers, and galaxies merge. Consequently, binary SMBHs should be common in galactic nuclei. However, close binaries of SMBH (CB-SMBHs) with sub-parsec separation cannot be imaged directly using current facilities. Some indirect signatures, such as periodic signals in light curves and double peaks in emission-line profile, have been used to find CB-SMBH candidates, but ambiguities still exist and no definitive conclusions can be made. We have recently proposed a new method focusing on kinematic signatures that can be derived from reverberation mapping of CB-SMBHs, one that offers a promising avenue to address this important problem. In this paper, we calculated models for a wide range of parameters, but BLRs of two BHs are close but still not merged. The purpose of this supplementary paper is to provide an atlas of two-dimensional transfer functions of CB-SMBHs with a wide range of orbital and geometrical parameters to aid more efficient identification of CB-SMBH candidates in reverberation mapping data.

In the present paper, we explore opportunities of applying the GRAVITY at the Very Large Telescope Interferometry (VLTI) with unprecedented spatial resolution to identify close binaries of supermassive black holes (CB-SMBHs) in active galactic nuclei (AGNs). Each SMBH is assumed to be separately surrounded by their own broad-line regions (BLRs) composed of clouds with virialized motion. Composition of the binary orbital motion and the virial motion of clouds in each BLR determines the projected velocity fields and hence differential phase curves, which are obviously different from that of a single BLR. We calculate emission line profiles and differential phase curves of CB-SMBHs for the GRAVITY. For the simplest case where angular momentums of two BLRs and orbital motion are parallel, a phase plateau generally appears in the phase curves. For other combinations of the angular momentum, the plateau is replaced by new peaks and valleys variously depending on the situations. Given a combination, phase curves are also sensitive to changes of parameters of CB-SMBHs. All these features are easily distinguished from the well-known $S$-shaped phase curves of a single BLR so that the GRAVITY is expected to reveal signals of CB-SMBH from candidates of AGNs. With joint analysis of observations of reverberation mapping campaigns, we can reliably identify CB-SMBHs, and measure their orbital parameters in the meanwhile. This independent measurement of the orbital parameters also has implications to analysis of Pulsar Timing Array (PTA) observations for properties of low-frequency gravitational waves in future.

This is the eighth in a series of papers reporting on a large reverberation mapping campaign to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). We employ the recently developed dynamical modeling approach for broad-line regions (BLRs) based on the method of Pancoast et al. to analyze the reverberation mapping dataset of Mrk 142 observed in the first monitoring season. In this approach, continuum variations are reconstructed using a damped random walk process, and BLR structure is delineated using a flexible disk-like geometry, in which BLR clouds move around the central BH with Keplerian orbits or inflow/outflow motion. The approach also includes the possibilities of anisotropic emission of BLR clouds, non-linear response of the line emission to the continuum, and different long-term trends in the continuum and emission-line variations. We implement the approach in a Bayesian framework that is apt for parallel computation and use a Markov Chain Monte Carlo technique to recover the parameters and uncertainties for the modeling, including mass of the central BH. We apply three BLR models with different prescriptions of BLR clouds distributions and find that the best model for fitting the data of Mrk 142 is a two-zone BLR model, consistent with the theoretical BLR model surrounding slim accretion disks. The best model yields a BH mass of $\log (M_\bullet/M_\odot)=6.23_{-0.45}^{+0.26}$, resulting in a virial factor of $\log f=-0.36_{-0.54}^{+0.33}$ for the full width at half maximum of the H$\beta$ line measured from the mean spectrum. The virial factors for the other measures of the H$\beta$ line width are also presented.

Close binaries of supermassive black holes (CB-SMBHs) with separations of $\lesssim 0.1$pc as the final stage of galaxy mergers are sources of low frequency gravitational waves (GW), however, they are still elusive observationally because they are not spatially resolved. Fortunately, reverberation as echoes of broad emission lines to ionizing continuum conveys invaluable information of the dynamics of broad-line regions (BLRs) governed by supermassive black holes in the central regions of active galactic nuclei (AGNs). In this paper, we demonstrate how to composite the hybrid 2-dimensional transfer functions of binary BLRs around the CB-SMBHs in AGNs, providing an opportunity of identifying them from reverberation mapping (RM) data. It is found that there are variation-coupling effects in the transfer functions, arising from the coupling of CB-SMBH light curves in the Fourier space. We provide semi-analytical formulations of the transfer functions for kinematic maps of the gas. For cases with the simplest variation-coupling effects, we make calculations for several BLR models and reveal significant distinctions from those of single active black holes. In principle, the difference is caused by the orbital motion of the CB-SMBH systems. In order to search for CB-SMBHs in time-domain space, selection of target candidates should focus on local AGNs with H$\beta$ double-peaked profiles and weaker near-infrared emission. High-fidelity RM-campaigns of monitoring the targets in future will provide opportunities to reveal these kinematic signatures of the CB-SMBHs and hence for measurements of their orbital parameters.

Type 1 active galactic nuclei display broad emission lines, regarded as arising from photoionized gas moving in the gravitational potential of a supermassive black hole. The origin of this broad-line region gas is unresolved so far, however. Another component is the dusty torus beyond the broad-line region, likely an assembly of discrete clumps that can hide the region from some viewing angles and make them observationally appear as Type 2 objects. Here we report that these clumps moving within the dust sublimation radius, like the molecular cloud G2 discovered in the Galactic center, will be tidally disrupted by the hole, resulting in some gas becoming bound at smaller radii while other gas is ejected and returns to the torus. The clumps fulfill necessary conditions to be photoionized. Specific dynamical components of tidally disrupted clumps include spiral-in gas as inflow, circularized gas, and ejecta as outflow. We calculate various profiles of emission lines from these clouds, and find they generally agree with H$\beta$ profiles of Palomar-Green quasars. We find that asymmetry, shape and shift of the profiles strongly depend on [O III], luminosity, which we interpret as a proxy of dusty torus angles. Tidally disrupted clumps from the torus may represent the source of the broad-line region gas.

Black hole (BH) mass of Type I active galactic nuclei (AGN) can be measured or estimated through either reverberation mapping (RM) or empirical $R-L$ relation, however, both of them suffer from uncertainties of the virial factor ($f_{\rm BLR}$), thus limiting the measurement accuracy. In this letter, we make an effort to investigate $f_{\rm BLR}$ through polarised spectra of the broad-line regions (BLR) arisen from electrons in the equatorial plane. Given the BLR composed of discrete clouds with Keplerian velocity around the central BH, we simulate a large number of spectra of total and polarised flux with wide ranges of parameters of the BLR model and equatorial scatters. We find that the $f_{\rm BLR}$-distribution of polarised spectra is much narrower than that of total ones. This provides a way of n accurately estimating BH mass from single spectropolarimetric observations of type I AGN whose equatorial scatters are identified.