Search

Over the past three decades, multiple reverberation mapping (RM) campaigns conducted for the quasar PG 2130+099 have exhibited inconsistent findings with time delays ranging from $\sim$10 to $\sim$200 days. To achieve a comprehensive understanding of the geometry and dynamics of the broad-line region (BLR) in PG 2130+099, we continued an ongoing high-cadence RM monitoring campaign using the Calar Alto Observatory 2.2m optical telescope for an extra four years from 2019 to 2022. We measured the time lags of several broad emission lines (including He II, He I, H$\beta$, and Fe II) with respect to the 5100 {\AA} continuum, and their time lags continuously vary through the years. Especially, the H$\beta$ time lags exhibited approximately a factor of two increase in the last two years. Additionally, the velocity-resolved time delays of the broad H$\beta$ emission line reveal a back-and-forth change between signs of virial motion and inflow in the BLR. The combination of negligible ($\sim$10%) continuum change and substantial time-lag variation (over two times) results in significant scatter in the intrinsic $R_{\rm H\beta}-L_{\rm 5100}$ relationship for PG 2130+099. Taking into account the consistent changes in the continuum variability time scale and the size of the BLR, we tentatively propose that the changes in the measurement of the BLR size may be affected by 'geometric dilution'., Comment: 21 pages, 13 figures, 7 tables; accepted for publication in ApJ

We conduct a reverberation mapping (RM) campaign to spectroscopically monitor a sample of selected bright active galactic nuclei with large anticipated broad-line region (BLR) sizes adequate for spectroastrometric observations by the GRAVITY instrument on the Very Large Telescope Interferometer. We report the first results for five objects, IC 4329A, Mrk 335, Mrk 509, Mrk 1239, and PDS 456, among which Mrk 1239 and PDS 456 are for the first time spectroscopically monitored. We obtain multi-year monitoring data and perform multi-component spectral decomposition to extract the broad H$\beta$ profiles. We detect significant time lags between the H$\beta$ and continuum variations, generally obeying the previously established BLR size-luminosity relation. Velocity-resolved H$\beta$ time lags illustrate diverse, possibly evolving BLR kinematics. We further measure the H$\beta$ line widths from mean and rms spectra and the resulting virial products show good consistency among different seasons. Adopting a unity virial factor and the full width at half maximum of the broad H$\beta$ line from the mean spectrum as the measure of velocity, the obtained black hole mass averaged over seasons is $\log M_\bullet/M_\odot=8.02_{-0.14}^{+0.09}$, $6.92_{-0.12}^{+0.12}$, $8.01_{-0.25}^{+0.16}$, $7.44_{-0.14}^{+0.13}$, and $8.59_{-0.11}^{+0.07}$ for the five objects, respectively. The black hole mass estimations using other line width measures are also reported (up to the virial factors). For objects with previous RM campaigns, our mass estimates are in agreement with earlier results. In a companion paper, we will employ BLR dynamical modeling to directly infer the black hole mass and thereby determine the virial factors., Comment: 32 pages, 6 tables, 20 figures. To appear in ApJ

We test the potential of Bayesian synthesis of upcoming multi-instrument data to extract orbital parameters and individual light curves of close binary supermassive black holes (CB-SMBH) with subparsec separations. Next generation (ng) interferometers, will make possible the observation of astrometric wobbles in CB-SMBH. Combining them with periodic variable time-domain data from surveys like the Vera C. Rubin Legacy Survey of Space and Time (LSST), allows for a more information on CB-SMBH candidates compared to standalone observational methods. Our method reliably determines binary parameters and component fluxes from binary total flux across long-term, intermediate and short-term binary dynamics and observational configurations, assuming ten annual observations, even in short period "q-accrete" objects. Expected CB-SMBH astrometric wobbles constructed from binary dynamical parameters, might serve in refining observational strategies for CB-SMBH. Combination of inferred mass ratio, light curves of binary components, and observed photocenter wobbles can be a proxy for the activity states of CB-SMBH components., Comment: The article was accepted for publication in the ApJ

We report first-time reverberation mapping results for 14 AGNs from the ongoing Monitoring AGNs with H$\beta$ Asymmetry campaign (MAHA). These results utilize optical spectra obtained with the Long Slit Spectrograph on the Wyoming Infrared 2.3m Telescope between 2017 November-2023 May. MAHA combines long-duration monitoring with high cadence. We report results from multiple observing seasons for 9 of the 14 objects. These results include H$\beta$ time lags, supermassive black hole masses, and velocity-resolved time lags. The velocity-resolved lags allow us to investigate the kinematics of the broad-line region., Comment: 35 pages, 19 figures, accepted for publication in ApJ Supplement

Measuring broad emission-line widths in active galactic nuclei (AGN) is not straightforward owing to the complex nature of flux variability in these systems. Line-width measurements become especially challenging when signal-to-noise is low, profiles are narrower, or spectral resolution is low. We conducted an extensive correlation analysis between emission-line measurements from the optical spectra of Markarian 142 (Mrk 142; a narrow-line Seyfert galaxy) taken with the Gemini North Telescope (Gemini) at a spectral resolution of 185.6+\-10.2 km/s and the Lijiang Telescope (LJT) at 695.2+\-3.9 km/s to investigate the disparities in the measured broad-line widths from both telescope data. Mrk~142 posed a challenge due to its narrow broad-line profiles, which were severely affected by instrumental broadening in the lower-resolution LJT spectra. We discovered that allowing the narrow-line flux of permitted lines having broad and narrow components to vary during spectral fitting caused a leak in the narrow-line flux to the broad component, resulting in broader broad-line widths in the LJT spectra. Fixing the narrow-line flux ratios constrained the flux leak and yielded the Hydrogen-beta broad-line widths from LJT spectra $\sim$54\% closer to the Gemini Hydrogen-beta widths than with flexible narrow-line ratios. The availability of spectra at different resolutions presented this unique opportunity to inspect how spectral resolution affected emission-line profiles in our data and adopt a unique method to accurately measure broad-line widths. Reconsidering line-measurement methods while studying diverse AGN populations is critical for the success of future reverberation-mapping studies. Based on the technique used in this work, we offer recommendations for measuring line widths in narrow-line AGN., Comment: 14 pages, 9 figures, accepted for publication in PASP

Pairs of supermassive black holes (SMBHs) at different stages are natural results of galaxy mergers in the hierarchical framework of galaxy formation and evolution. However, identifications of close binaries of SMBHs (CB-SMBHs) with sub-parsec separations in observations are still elusive. Recently, unprecedented spatial resolutions achieved by GRAVITY/GRAVITY+ onboard The Very Large Telescope Interferometer through spectroastrometry (SA) provide new opportunities to resolve CB-SMBHs. Differential phase curves of CB-SMBHs with two independent broad-line regions (BLRs) are found to have distinguished characteristic structures from a single BLR \citep{songsheng2019}. Once the CB-SMBH evolves to the stage where BLRs merge to form a circumbinary BLR, it will hopefully be resolved by the pulsar timing array (PTA) in the near future as sources of nano-hertz gravitational waves. In this work, we use a parameterized model for circumbinary BLRs to calculate line profiles and differential phase curves for SA observations. We show that both profiles and phase curves exhibit asymmetries caused by the Doppler boosting effect of accretion disks around individual black holes, depending on the orbital parameters of the binary and geometries of the BLR. We also generate mock SA data using the model and then recover orbital parameters by fitting the mock data. Degeneracies between parameters contribute greatly to uncertainties of parameters but can be eased through joint analysis of multiple-epoch SA observations and reverberation mappings., Comment: 21 pages, 17 figures, accepted by ApJ

There are increasing interests in binary supermassive black holes (SMBHs), but merging binaries with separations smaller than ~1 light days (~10^2 gravitational radii for 10^8 Msun), which are rapidly evolving under control of gravitational waves, are elusive in observations. In this paper, we discuss fates of mini-disks around component SMBHs for three regimes: 1) low rates (advection-dominated accretion flows: ADAFs); 2) intermediate rates; 3) super-Eddington accretion rates. Mini-disks with intermediate rates are undergoing evaporation through thermal conduction of hot corona forming a hybrid radial structure. When the binary orbital periods are shorter than sound propagation timescales of the evaporated mini-disks, a new instability, denoted as sound instability, arises because the disks will be highly twisted so that they are destroyed. We demonstrate a critical separation of A_{crit}~10^2 Rg from the sound instability of the mini-disks and the cavity is full of hot gas. For those binaries, component SMBHs are accreting with Bondi mode in the ADAF regime, showing periodic variations resulting from Doppler boosting effects in radio from the ADAFs due to orbital motion. In the mean while, the circumbinary disks (CBDs) are still not hot enough (ultraviolet deficit) to generate photons to ionize gas for broad emission lines. For slightly super-Eddington accretion of the CBDs, MgII line appears with decreases of UV deficit, and for intermediate super-Eddington Balmer lines appear, but CIV line never unless CBD accretion rates are extremely high. Moreover, if the CBDs are misaligned with the binary plane, it is then expected to have optical periodical variations with about ten times radio periods., Comment: 10 pages, 3 figures; accepted to MNRAS

In this third paper of the series reporting on the reverberation mapping (RM) campaign of active galactic nuclei with asymmetric H$\beta$ emission-line profiles, we present results for 15 Palomar-Green (PG) quasars using spectra obtained between the end of 2016 to May 2021. This campaign combines long time spans with relatively high cadence. For 8 objects, both the time lags obtained from the entire light curves and the measurements from individual observing seasons are provided. Reverberation mapping of 9 of our targets has been attempted for the first time, while the results for 6 others can be compared with previous campaigns. We measure the H$\beta$ time lags over periods of years and estimate their black hole masses. The long duration of the campaign enables us to investigate their broad line region (BLR) geometry and kinematics for different years by using velocity-resolved lags, which demonstrate signatures of diverse BLR geometry and kinematics. The BLR geometry and kinematics of individual objects are discussed. In this sample, the BLR kinematics of Keplerian/virialized motion and inflow is more common than outflow., Comment: 35 pages, 18 figures, 8 tables, accepted for publication in The Astrophysical Journal Supplement

We conduct a systematic search for quasars with periodic variations from the archival photometric data of the Zwicky Transient Facility by cross-matching with the quasar catalogs of the Sloan Digital Sky Survey and V{\'e}ron-Cetty \& V{\'e}ron. We first select out 184 primitive periodic candidates using the generalized Lomb-Scargle periodogram and auto-correlation function and then estimate their statistical significance of periodicity based on two red-noise models, i.e., damped random walk (DRW) and single power-law (SPL) models. As such, we finally identify 106 (DRW) and 86 (SPL) candidates with the most significant periodic variations out of 143,700 quasars. We further compare DRW and SPL models using Bayes factors, which indicate a relative preference of the SPL model for our primitive sample. We thus adopt the candidates identified with SPL as the final sample and summarize its basic properties. We extend the light curves of the selected candidates by supplying other archival survey data to verify their periodicity. However, only three candidates (with 6-8 cycles of periods) meet the selection criteria. This result clearly implies that, instead of being strictly periodic, the variability must be quasi-periodic or caused by stochastic red-noise. This exerts a challenge to the existing search approaches and calls for developing new effective methods., Comment: Accepted in MNRAS. Comments welcome

Recent studies have proposed that GRAVITY+ instrument is able to trace the circular orbit of the subparsec close-binary supermassive black holes (CB-SMBHs) by measuring the photocentre variation of the hot dust emission. However, the CB-SMBHs orbit may become highly eccentric throughout the evolution of these objects, and the orbital period may be far longer than the observational time baseline. We investigate the problem of detecting the CB-SMBH with hot dust emission and high eccentricity (eCBSMBH, e=0.5) when the observed time baselines of their astrometric data and radial velocities are considerably shorter than the orbital period. The parameter space of the Keplerian model of the eCBSMBH is large for exploratory purposes. We therefore applied the Bayesian method to fit orbital elements of the eCBSMBH to combined radial velocity and astrometric data covering a small fraction of the orbital period. We estimate that a number of potential eCBSMBH systems within reach of GRAVITY + will be similar to the number of the planned circular targets. We show that using observational time baselines that cover ~ 10% of the orbit increases the possibility of determining the period, eccentricity, and total mass of an eCBSMBH. When the observational time baseline becomes too short (~ 5%), the quality of the retrieved eCBSMBH parameters degrades. We also illustrate how interferometry may be used to estimate the photo-centre at the eCBSMBH emission line, which could be relevant for GRAVITY+ successors. Even if the astrometric signal for eCBSMBH systems is reduced by a factor of sqrt{1-e^{2}} compared to circular ones, we find that the hot dust emission of eCBSMBHs can be traced by GRAVITY+ at the elementary level., Comment: accepted for publication in the journal Astronomy & Astrophysics

In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $\sigma_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,\sigma$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $\Omega_m$, and the amplitude or rate of the growth of structure ($\sigma_8,f\sigma_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged], Comment: Contribution to Snowmass 2021. 224 pages, 27 figures. Accepted for publication in JHEAp

The quasar 3C 273 has been observed with infrared spectroastrometry (SA) on broad Pa$\alpha$ line and optical reverberation mapping (RM) on broad H$\beta$ line. SA delivers information about the angular size and structure of the Pa$\alpha$ broad-line region (BLR), while RM delivers information about the physical size and structure of the H$\beta$ BLR. Based on the fact that the two BLRs share the mass of the supermassive black hole (SMBH) and viewing inclination, a combination of SA and velocity-resolved RM (SARM) thereby allows us to simultaneously determine the SMBH mass and geometric distance through dynamically modeling the two BLRs. We construct a suite of dynamical models with different geometric configurations and apply a Bayesian approach to obtain the parameter inferences. Overall the obtained masses and distances are insensitive to specific BLR configurations but more or less depend on parameterizations of the vertical distributions. The most probable model, chosen in light of the Bayes factor, yields an angular-size distance of $\log\,(D_{\rm A}/{\rm Mpc}) = 2.83_{-0.28}^{+0.32}$ and SMBH mass of $\log\,(M_\bullet/M_\odot)=9.06_{-0.27}^{+0.21}$, which agrees with the relationships between SMBH masses and bulge properties. The BLRs have an inclination of $5_{-1}^{+1}$ degrees, consistent with that of the large-scale jet in 3C 273. Our approach reinforces the capability of SARM analysis to measure SMBH mass and distance of AGNs even though SA and RM observations are undertaken with different emission lines and/or in different periods., Comment: 24 pages, 16 figures, 3 tables; ApJ in press; the code BRAINS available at https://github.com/LiyrAstroph/BRAINS

Detecting continuous nanohertz gravitational waves (GWs) generated by individual close binaries of supermassive black holes (CB-SMBHs) is one of the primary objectives of pulsar timing arrays (PTAs). The detection sensitivity is slated to increase significantly as the number of well-timed millisecond pulsars will increase by more than an order of magnitude with the advent of next-generation radio telescopes. Currently, the Bayesian analysis pipeline using parallel tempering Markov chain Monte Carlo has been applied in multiple studies for CB-SMBH searches, but it may be challenged by the high dimensionality of the parameter space for future large-scale PTAs. One solution is to reduce the dimensionality by maximizing or marginalizing over uninformative parameters semi-analytically, but it is not clear whether this approach can be extended to more complex signal models without making overly simplified assumptions. Recently, the method of diffusive nested (DNest) sampling shown the capability of coping with high dimensionality and multimodality effectively in Bayesian analysis. In this paper, we apply DNest to search for continuous GWs in simulated pulsar timing residuals and find that it performs well in terms of accuracy, robustness, and efficiency for a PTA including $\mathcal{O}(10^2)$ pulsars. DNest also allows a simultaneous search of multiple sources elegantly, which demonstrates its scalability and general applicability. Our results show that it is convenient and also high beneficial to include DNest in current toolboxes of PTA analysis., Comment: 20 pages, 10 figures, submitted to ApJ

We report the results of a multi-year spectroscopic and photometric monitoring campaign of two luminous quasars, PG~0923+201 and PG~1001+291, both located at the high-luminosity end of the broad-line region (BLR) size-luminosity relation with optical luminosities above $10^{45}~{\rm erg~s^{-1}}$. PG~0923+201 is for the first time monitored, and PG~1001+291 was previously monitored but our campaign has a much longer temporal baseline. We detect time lags of variations of the broad H$\beta$, H$\gamma$, Fe {\sc ii} lines with respect to those of the 5100~{\AA} continuum. The velocity-resolved delay map of H$\beta$ in PG~0923+201 indicates a complicated structure with a mix of Keplerian disk-like motion and outflow, and the map of H$\beta$ in PG~1001+291 shows a signature of Keplerian disk-like motion. Assuming a virial factor of $f_{\rm BLR}=1$ and FWHM line widths, we measure the black hole mass to be $118_{-16}^{+11}\times 10^7 M_{\odot}$ for PG~0923+201 and $3.33_{-0.54}^{+0.62}\times 10^7 M_{\odot}$ for PG~1001+291. Their respective accretion rates are estimated to be $0.21_{-0.07}^{+0.06} \times L_{\rm Edd}\,c^{-2}$ and $679_{-227}^{+259}\times L_{\rm Edd}\,c^{-2}$, indicating that PG~0923+201 is a sub-Eddington accretor and PG~1001+291 is a super-Eddington accretor. While the H$\beta$ time lag of PG~0923+201 agrees with the size-luminosity relation, the time lag of PG~1001+291 shows a significant deviation, confirming that in high-luminosity AGN the BLR size depends on both luminosity and Eddington ratio. Black hole mass estimates from single AGN spectra will be over-estimated at high luminosities and redshifts if this effect is not taken into account., Comment: 21 pages, 14 figures, accepted

Recently, GRAVITY onboard the Very Large Telescope Interferometer (VLTI) first spatially resolved the structure of the quasar 3C 273 with an unprecedented resolution of $\sim 10\mu$as. A new method of measuring parallax distances has been successfully applied to the quasar through joint analysis of spectroastrometry (SA) and reverberation mapping (RM) observation of its broad line region (BLR). The uncertainty of this SA and RM (SARM) measurement is about $16\%$ from real data, showing its great potential as a powerful tool for precision cosmology. In this paper, we carry out detailed analyses of mock data to study impacts of data qualities of SA observations on distance measurements and establish a quantitative relationship between statistical uncertainties of distances and relative errors of differential phases. We employ a circular disk model of BLR for the SARM analysis. We show that SARM analyses of observations generally generate reliable quasar distances, even for relatively poor SA measurements with error bars of $40\%$ at peaks of phases. Inclinations and opening angles of BLRs are the major parameters governing distance uncertainties. It is found that BLRs with inclinations $\gtrsim 10^{\circ}$ and opening angles $\lesssim 40^{\circ}$ are the most reliable regimes from SARM analysis for distance measurements. Through analysis of a mock sample of AGNs generated by quasar luminosity functions, we find that if the GRAVITY/GRAVITY+ can achieve a phase error of $0.1^{\circ}$ per baseline for targets with magnitudes $K \lesssim 11.5$, the SARM campaign can constrain $H_0$ to an uncertainty of $2\%$ by observing $60$ targets., Comment: 20 pages, 13 figures, accepted by ApJS

We report on the results of a new spectroscopic monitoring campaign of the quasar PG 0026+129 at the Calar Alto Observatory 2.2m telescope from July 2017 to February 2020. Significant variations in the fluxes of the continuum and broad-emission lines, including H$\beta$ and He II, were observed in the first and third years, and clear time lags between them are measured. The broad H$\beta$ line profile consists of two Gaussian components: an intermediate-width H$\beta_{\rm IC}$ with a full width at half-maximum (FWHM) of 1964$\pm$18 $\rm km~s^{-1}$ and another very broad H$\beta_{\rm VBC}$ with a FWHM of 7570$\pm$83 $\rm km~s^{-1}$. H$\beta_{\rm IC}$ has long time lags of $\sim$40--60 days in the rest frame, while H$\beta_{\rm VBC}$ shows nearly zero time delay with respect to the optical continuum at 5100 \AA. The velocity-resolved delays show consistent results: lags of $\sim$30--50 days at the core of the broad H$\beta$ line and roughly zero lags at the wings. H$\beta_{\rm IC}$ has a redshift of $\sim$400 $\rm km~s^{-1}$ which seems to be stable for nearly 30 years by comparing with archived spectra, and may originate from an infall. The root mean square (rms) spectrum of H$\beta_{\rm VBC}$ shows a double-peaked profile with brighter blue peak and extended red wing in the first year, which matches the signature of a thin disk. Both the double-peaked profile and the near-zero lag suggest that H$\beta_{\rm VBC}$ comes from a region associated with the part of the accretion disc that emits the optical continuum. Adopting the FWHM (in the rms spectrum) and the time lag measured for the total H$\beta$ line, and a virial factor of 1.5, we obtain a virial mass of $2.89_{-0.69}^{+0.60} \times10^7 M_{\odot}$ for the central black hole in this quasar., Comment: 21 pages, 17 figures, accepted for publication in The Astrophysical Journal

The Very Large Telescope Interferometer (VLTI), and the Extremely Large Telescope (ELT) will be a robust astrophysics suite offering the opportunity of probing the structure and dynamics of CB-SMBH at high spectral and angular resolution. Here, we explore and illustrate the application of differential interferometry on unresolved the CB-SMBH systems in elliptical orbital configurations and a single SMBH with clouds in elliptical orbital motion. Photocenter displacements between each SMBH and regions in their disc-like broad line regions (BLR) appear as small interferometric differential phase variability. To investigate the application of interferometric phases for the detection of CB-SMBH systems, we simulate a series of differential interferometry signatures, based on our model comprising ensembles of clouds surrounding each of supermassive black hole in a CB-SMBH. Setting model to the parameters of a single SMBH with elliptical cloud motion, we also calculated a series of differential interferometry observables for this case. We found various deviations from the canonical S-shaped of CB-SMBH phases profile for elliptically configured CB-SMBH systems. The amplitude and specific shape of the interferometry observables depend on orbital configurations of the CB-SMBH system. We get distinctive results when considering antialigned angular momenta of cloud orbits regarding total CB-SMBH angular momentum. Some simulated spectral lines from our model closely resemble observations of Pa \alpha line got from near-infrared AGN surveys. We found differences between differential phases zoo of a single SMBH and CB-SMBH systems. The differential phases zoo for a single SMBH comprises deformed S shape. We also showed how their differential phase shape, amplitude, and slope evolve with various sets of cloud orbital parameters and observer position., Comment: accepted by A&A Journal

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., Comment: 24 pages, 14 figures, accepted by ApJS

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., Comment: 52pages, 9 figures

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., Comment: 16 pages, 6 figures, accepted by ApJ

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., Comment: 21 pages, 12 figures, 5 tables; ApJ accepted; The software BRAINS available at https://github.com/LiyrAstroph/BRAINS

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., Comment: 22 pages, 8 figures, Accepted by The Astrophysical Journal

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., Comment: appear in Nature Astronomy on October 9, 59 pages, 9 figures, 3 tables

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., Comment: MNRAS Letters (accepted): 5 pages, 3 figures