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.

3C 273 is an intensively monitored flat-spectrum radio quasar with both a beamed jet and blue bump together with broad emission lines. The coexistence of the comparably prominent jet and accretion disk leads to complicated variability properties. Recent reverberation mapping monitoring for 3C 273 revealed that the optical continuum shows a distinct long-term trend that does not have a corresponding echo in the Hbeta fluxes. We compile multi-wavelength monitoring data from the Swift archive and other ground-based programs and clearly find two components of emissions at optical wavelength. One component stems from the accretion disk itself and the other component can be ascribed to the jet contribution, which also naturally accounts for the non-echoed trend in reverberation mapping data. We develop an approach to decouple the optical emissions from the jet and accretion disk in 3C 273 with the aid of multi-wavelength monitoring data. By assuming the disk emission has a negligible polarization in consideration of the low inclination of the jet, the results show that the jet contributes a fraction of ~10% at the minimum and up to ~40% at the maximum to the total optical emissions. This is the first time to provide a physical interpretation to the "detrending" manipulation conventionally adopted in reverberation mapping analysis. Our work also illustrates the importance of appropriately analyzing variability properties in cases of coexisting jets and accretion disks.

We performed an intensive accretion disk reverberation mapping campaign on the high accretion rate active galactic nucleus Mrk 142 in early 2019. Mrk 142 was monitored with the Neil Gehrels Swift Observatory for 4 months in X-rays and 6 UV/optical filters. Ground-based photometric monitoring was obtained from the Las Cumbres Observatory, Liverpool Telescope and Dan Zowada Memorial Observatory in ugriz filters and the Yunnan Astronomical Observatory in V. Mrk 142 was highly variable throughout, displaying correlated variability across all wavelengths. We measure significant time lags between the different wavelength light curves, finding that through the UV and optical the wavelength-dependent lags, $\tau(\lambda)$, generally follow the relation $\tau(\lambda) \propto \lambda^{4/3}$, as expected for the $T\propto R^{-3/4}$ profile of a steady-state optically-thick, geometrically-thin accretion disk, though can also be fit by $\tau(\lambda) \propto \lambda^{2}$, as expected for a slim disk. The exceptions are the u and U band, where an excess lag is observed, as has been observed in other AGN and attributed to continuum emission arising in the broad-line region. Furthermore, we perform a flux-flux analysis to separate the constant and variable components of the spectral energy distribution, finding that the flux-dependence of the variable component is consistent with the $f_\nu\propto\nu^{1/3}$ spectrum expected for a geometrically-thin accretion disk. Moreover, the X-ray to UV lag is significantly offset from an extrapolation of the UV/optical trend, with the X-rays showing a poorer correlation with the UV than the UV does with the optical. The magnitude of the UV/optical lags is consistent with a highly super-Eddington accretion rate.

We investigate whether the hard X-ray photon index (${\Gamma}$) versus accretion rate correlation for super-Eddington accreting quasars is different from that for sub-Eddington accreting quasars. We construct a sample of 113 bright quasars from the Sloan Digital Sky Survey Data Release 14 quasar catalog, including 38 quasars as the super-Eddington subsample and 75 quasars as the sub-Eddington subsample. We derive black-hole masses using a simple-epoch virial mass formula based on the ${\rm H\beta}$ lines, and we use the standard thin disk model to derive the dimensionless accretion rates ($\dot{\mathscr{M}}$) for our sample. The X-ray data for these quasars are collected from the Chandra and XMM-Newton archives. We fit the hard X-ray spectra using a single power-law model to obtain ${\Gamma}$ values. We find a statistically significant ($R_{\rm S}=0.43$, $p=7.75\times{10}^{-3}$) correlation between ${\Gamma}$ and $\dot{\mathscr{M}}$ for the super-Eddington subsample. The ${\Gamma}$-$\dot{\mathscr{M}}$ correlation for the sub-Eddington subsample is also significant, but weaker ($R_{\rm S}=0.30$, $p=9.98\times{10}^{-3}$). Linear regression analysis shows that ${\rm \Gamma}=(0.34\pm0.11){\rm log}{\dot{\mathscr{M}}}+(1.71\pm0.17)$ and ${\rm \Gamma}=(0.09\pm0.04){\rm log}{\dot{\mathscr{M}}}+(1.93\pm0.04)$ for the super- and sub-Eddington subsamples, respectively. The ${\Gamma}$-$\dot{\mathscr{M}}$ correlations of the two subsamples are different, suggesting different disk-corona connections in these two types of systems. We propose one qualitative explanation of the steeper ${\Gamma}$-$\dot{\mathscr{M}}$ correlation in the super-Eddington regime that involves larger seed photon fluxes received by the compact coronae from the thick disks in super-Eddington accreting quasars.

At present, the strongest upper limit on $\sum m_{\nu}$, the sum of neutrino masses, is from cosmological measurements. However, this bound assumes that the neutrinos are stable on cosmological timescales, and is not valid if the neutrino lifetime is less than the age of the universe. In this paper, we explore the cosmological signals of theories in which the neutrinos decay into invisible dark radiation on timescales of order the age of the universe, and determine the bound on the sum of neutrino masses in this scenario. We focus on the case in which the neutrinos decay after becoming non-relativistic. We derive the Boltzmann equations that govern the cosmological evolution of density perturbations in the case of unstable neutrinos, and solve them numerically to determine the effects on the matter power spectrum and lensing of the cosmic microwave background. We find that the results admit a simple analytic understanding. We then use these results to perform a Monte Carlo analysis based on the current data to determine the limit on the sum of neutrino masses as a function of the neutrino lifetime. We show that in the case of decaying neutrinos, values of $\sum m_{\nu}$ as large as 0.9 eV are still allowed by the data. Our results have important implications for laboratory experiments that have been designed to detect neutrino masses, such as KATRIN and KamLAND-ZEN.

We explore the cosmological signals of theories in which the neutrinos decay into invisible dark radiation after becoming non-relativistic. We show that in this scenario, near-future large scale structure measurements from the Euclid satellite, when combined with cosmic microwave background data from Planck, may allow an independent determination of both the lifetime of the neutrinos and the sum of their masses. These parameters can be independently determined because the Euclid data will cover a range of redshifts, allowing the growth of structure over time to be tracked. If neutrinos are stable on cosmological timescales, these observations can improve the lower limit on the neutrino lifetime by seven orders of magnitude, from $\mathcal{O}(10)$ years to $2\times 10^8$ years ($95\%$ C.L.), without significantly affecting the measurement of neutrino mass. On the other hand, if neutrinos decay after becoming non-relativistic but on timescales less than $\mathcal{O}(100)$ million years, these observations may allow, not just the first measurement of the sum of neutrino masses, but also the determination of the neutrino lifetime from cosmology.

As one of the most interesting Seyfert 1 galaxies, PG 2130+099 has been the target of several reverberation mapping (RM) campaigns over the years. However, its measured broad H$\beta$ line responses have been inconsistent, with time lags of $\sim$200 days, $\sim$25 days, and $\sim$10 days being reported for different epochs while its optical luminosity changed no more than 40%. To investigate this issue, we conducted a new RM-campaign with homogenous and high cadence (about $\sim$3 days) for two years during 2017--2019 to measure the kinematics and structure of the ionized gas. We successfully detected time lags of broad H$\beta$, He II, He I, and Fe II lines with respect to the varying 5100\AA continuum, revealing a stratified structure that is likely virialized with Keplerian kinematics in the first year of observations, but an inflow kinematics of the broad-line region from the second year. With a central black hole mass of $0.97_{-0.18}^{+0.15}\times 10^7~M_{\odot}$, PG 2130+099 has an accretion rate of $10^{2.1\pm0.5}L_{\rm Edd}c^{-2}$, where $L_{\rm Edd}$ is the Eddington luminosity and $c$ is speed of light, implying that it is a super-Eddington accretor and likely possesses a slim, rather than thin, accretion disk. The fast changes of the ionization structures of the three broad lines remain puzzling.

We developed a spectroscopic monitoring project to investigate the kinematics of the broad-line region (BLR) in active galactic nuclei (AGN) with ultra-fast outflows (UFOs). Mrk~79 is a radio-quiet AGN with UFOs and warm absorbers, had been monitored by three reverberation mapping (RM) campaigns, but its BLR kinematics is not understood yet. In this paper, we report the results from a new RM-campaign of Mrk~79, which was undertaken by Lijiang 2.4-m telescope. Mrk~79 is seeming to come out the faint state, the mean flux approximates a magnitude fainter than historical record. We successfully measured the lags of the broad emission lines including H$\beta~\lambda4861$, H$\gamma~\lambda4340$, He II $\lambda4686$ and He I $\lambda5876$ with respect to the varying AGN continuum. Based on the broad H$\beta~\lambda4861$ line, we measured black hole (BH) mass of $M_{\bullet}=5.13^{+1.57}_{-1.55}\times10^{7}M_{\odot}$, estimated accretion rates of ${\dot{M}_{\bullet}}=(0.05\pm0.02)~L_{\rm Edd}~c^{-2}$, indicating that Mrk~79 is a sub-Eddington accretor. We found that Mrk~79 deviates from the canonical Radius$-$Luminosity relationship. The marginal blueshift of the broad He II $\lambda4686$ line detected from rms spectrum indicates outflow of high-ionization gas. The velocity-resolved lag profiles of the broad H$\gamma~\lambda4340$, H$\beta~\lambda4861$, and He I $\lambda5876$ lines show similar signatures that the largest lag occurs in the red wing of the lines then the lag decreases to both sides. These signatures should suggest that the BLR of Keplerian motion probably exists the outflow gas motion. All findings including UFOs, warm absorbers, and the kinematics of high- and low-ionization BLR, may provide an indirect evidence that the BLR of Mrk~79 probably originates from disk wind.

We study the long-term variability in the optical monitoring database of Ark~120, a nearby radio-quiet active galactic nucleus (AGN) at a distance of 143 Mpc (z=0.03271). We compiled the historical archival photometric and spectroscopic data since 1974 and conducted a new two-year monitoring campaign in 2015-2017, resulting in a total temporal baseline over four decades. The long-term variations in the optical continuum exhibit a wave-like pattern and the Hbeta integrated flux series varies with a similar behavior. The broad Hbeta profiles have asymmetric double peaks, which change strongly with time and tend to merge into a single peak during some epochs. The period in the optical continuum determined from various period-search methods is about 20 yr and the estimated false alarm probability with null hypothesis simulations is about 1*10^-3. The overall variations of the broad Hbeta profiles also follow the same period. However, the present database only covers two cycles of the suggested period, which strongly encourages continued monitoring to track more cycles and confirm the periodicity. Nevertheless, in light of the possible periodicity and the complicated Hbeta profile, Ark~120 is one candidate of the nearest radio-quiet AGNs with possible periodic variability, and it is thereby a potential candidate host for a sub-parsec supermassive black hole binary.

The radius-luminosity (R-L) relationship of active galactic nuclei (AGNs) established by the reverberation mapping (RM) observations has been widely used as a single-epoch black hole mass estimator in the research of large AGN samples. However, the recent RM campaigns discovered that the AGNs with high accretion rates show shorter time lags by factors of a few comparing with the predictions from the R-L relationship. The explanation of the shortened time lags has not been finalized yet. We collect 8 different single-epoch spectral properties to investigate how the shortening of the time lags correlate with those properties and to understand what is the origin of the shortened lags. We find that the flux ratio between Fe II and H$\beta$ emission lines shows the most prominent correlation, thus confirm that accretion rate is the main driver for the shortened lags. In addition, we establish a new scaling relation including the relative strength of Fe II emission. This new scaling relation can provide less biased estimates of the black hole mass and accretion rate from the single-epoch spectra of AGNs.

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.

We compiled a sample of 73 active galactic nuclei (AGNs) with reverberation mapping (RM) observations from RM campaigns including our ongoing campaign of monitoring super-Eddington accreting massive black holes (SEAMBHs). This sample covers a large range of black hole (BH) mass $(M_{\bullet}=10^{6-9}~M_\odot)$, dimensionless accretion rates $(\dot{\mathscr{M}}=10^{-2.7}-10^{2.7})$ and 5100~\AA~luminosity $(L_{5100}=10^{42-46}~\rm erg~s^{-1})$, allowing us to systematically study the AGN variability and their relations with BH mass, accretion rates, and optical luminosity. We employed the damped random walk (DRW) model to delineate the optical variability of continuum at 5100~\AA~and obtained damped variability timescale ($\tau_{\rm d}$) and amplitude ($\sigma_{\rm d}$) using a Markov Chain Monte Carlo (MCMC) method. We also estimated the traditional variability amplitudes ($F_{\rm var}$), which provide a model-independent measure and therefore are used to test the DRW results. We found that AGN variability characteristics are generally correlated with $(M_{\bullet},\dot{\mathscr{M}},L_{5100})$. These correlations are smooth from sub-Eddington to super-Eddington accretion AGNs, probably implying that the AGN variability may be caused by the same physical mechanism.

We report the discovery of extreme X-ray variability in a type 1 quasar: SDSS J$075101.42+291419.1$. It has a black hole mass of $1.6\times 10^7~\rm M_\odot$ measured from reverberation mapping (RM), and the black hole is accreting with a super-Eddington accretion rate. Its XMM-Newton observation in 2015 May reveals a flux drop by a factor of $\sim 22$ with respect to the Swift observation in 2013 May when it showed a typical level of X-ray emission relative to its UV/optical emission. The lack of correlated UV variability results in a steep X-ray-to-optical power-law slope ($\alpha_{\rm OX}$) of -1.97 in the low X-ray flux state, corresponding to an X-ray weakness factor of 36.2 at rest-frame 2 keV relative to its UV/optical luminosity. The mild UV/optical continuum and emission-line variability also suggest that the accretion rate did not change significantly. A single power-law model modified by Galactic absorption describes well the $0.3-10$ keV spectra of the X-ray observations in general. The spectral fitting reveals steep spectral shapes with $\Gamma\approx3$. We search for active galactic nuclei (AGNs) with such extreme X-ray variability in the literature and find that most of them are narrow-line Seyfert 1 galaxies and quasars with high accretion rates. The fraction of extremely X-ray variable objects among super-Eddington accreting AGNs is estimated to be $\approx 15-24\%$. We discuss two possible scenarios, disk reflection and partial covering absorption, to explain the extreme X-ray variability of SDSS J$075101.42+291419.1$. We propose a possible origin for the partial covering absorber, which is the thick inner accretion disk and its associated outflow in AGNs with high accretion rates.

We report results of the first reverberation mapping campaign of I Zwicky 1 during $2014$-$2016$, which showed unambiguous reverberations of the broad H$\beta$ line emission to the varying optical continuum. From analysis using several methods, we obtain a reverberation lag of $\tau_{\rm H\beta}=37.2^{+4.5}_{-4.9}\,$ days. Taking a virial factor of $f_{_{\rm BLR}}=1$, we find a black hole mass of $M_{\bullet}=9.30_{-1.38}^{+1.26}\times 10^6 M_{\odot}$ from the mean spectra. The accretion rate is estimated to be $203.9_{-65.8}^{+61.0}\,L_{\rm Edd}c^{-2}$, suggesting a super-Eddington accretor, where $L_{\rm Edd}$ is the Eddington luminosity and $c$ is the speed of light. By decomposing {\it Hubble Space Telescope} images, we find that the stellar mass of the bulge of its host galaxy is $\log (M_{\rm bulge}/M_{\odot}) = \rm 10.92\pm 0.07$. This leads to a black hole to bulge mass ratio of $\sim 10^{-4}$, which is significantly smaller than that of classical bulges and elliptical galaxies. After subtracting the host contamination from the observed luminosity, we find that I Zw 1 follows the empirical $R_{\rm BLR}\propto L_{5100}^{1/2}$ relation.

Despite many decades of study, the kinematics of the broad-line region of 3C~273 are still poorly understood. We report a new, high signal-to-noise, reverberation mapping campaign carried out from November 2008 to March 2018 that allows the determination of time lags between emission lines and the variable continuum with high precision. The time lag of variations in H$\beta$ relative to those of the 5100 Angstrom continuum is $146.8_{-12.1}^{+8.3}$ days in the rest frame, which agrees very well with the Paschen-$\alpha$ region measured by the GRAVITY at The Very Large Telescope Interferometer. The time lag of the H$\gamma$ emission line is found to be nearly the same as for H$\beta$. The lag of the Fe II emission is $322.0_{-57.9}^{+55.5}$ days, longer by a factor of $\sim$2 than that of the Balmer lines. The velocity-resolved lag measurements of the H$\beta$ line show a complex structure which can be possibly explained by a rotation-dominated disk with some inflowing radial velocity in the H$\beta$-emitting region. Taking the virial factor of $f_{\rm BLR} = 1.3$, we derive a BH mass of $M_{\bullet} = 4.1_{-0.4}^{+0.3} \times 10^8 M_{\odot}$ and an accretion rate of $9.3\,L_{\rm Edd}\,c^{-2}$ from the H$\beta$ line. The decomposition of its $HST$ images yields a host stellar mass of $M_* = 10^{11.3 \pm 0.7} M_\odot$, and a ratio of $M_{\bullet}/M_*\approx 2.0\times 10^{-3}$ in agreement with the Magorrian relation. In the near future, it is expected to compare the geometrically-thick BLR discovered by the GRAVITY in 3C 273 with its spatially-resolved torus in order to understand the potential connection between the BLR and the torus.

Most results of the reverberation monitoring of active galaxies showed a universal scaling of the time delay of the Hbeta emission region with the monochromatic flux at 5100 A, with very small dipersion. Such a scaling favored the dust-based formation mechanism of the Broad Line Region (BLR). Recent reverberation measurements showed that actually a significant fraction of objects exhibits horter lags than the previously found scaling. Here we demonstrate that these shorter lags can be explained by the old concept of scaling of the BLR size with the ionization parameter. Assuming a universal value of this parameter and universal value of the cloud density reproduces the distribution of observational points in the time delay vs. monochromatic flux plane, provided that a range of black hole spins is allowed. However, a confirmation of the new measurements for low/moderate Eddington ratio sources is strongly needed before the dust-based origin of the BLR can be excluded.

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.

We have started a long-term reverberation mapping project using the Wyoming Infrared Observatory 2.3 meter telescope titled "Monitoring AGNs with H\beta\ Asymmetry" (MAHA). The motivations of the project are to explore the geometry and kinematics of the gas responsible for complex H\beta\ emission-line profiles, ideally leading to an understanding of the structures and origins of the broad-line region (BLR). Furthermore, such a project provides the opportunity to search for evidence of close binary supermassive black holes. We describe MAHA and report initial results from our first campaign, from December 2016 to May 2017, highlighting velocity-resolved time lags for four AGNs with asymmetric H\beta\ lines. We find that 3C 120, Ark 120, and Mrk 6 display complex features different from the simple signatures expected for pure outflow, inflow, or a Keplerian disk. While three of the objects have been previously reverberation mapped, including velocity-resolved time lags in the cases of 3C 120 and Mrk 6, we report a time lag and corresponding black hole mass measurement for SBS 1518+593 for the first time. Furthermore, SBS 1518+593, the least asymmetric of the four, does show velocity-resolved time lags characteristic of a Keplerian disk or virialized motion more generally. Also, the velocity-resolved time lags of 3C 120 have significantly changed since previously observed, indicating an evolution of its BLR structure. Future analyses of the data for these objects and others in MAHA will explore the full diversity of H\beta\ lines and the physics of AGN BLRs.

NGC 5548 has been well spectroscopically monitored for reverberation mapping of the central kinematics by 19 campaigns. Using the maximum entropy method in this Letter, we build up a high-quality velocity-delay map of the H$\beta$ emission line in the light curves of the continuum and the line variations observed between 2015-2016. The map shows the response strength and lags of the velocity fields of the H$\beta$ emitting regions. The velocity-delay structure of the map is generally symmetric, with strong red and blue wings at time lag $\tau \leq 15$ days, a narrower velocity distribution at $\tau \geq 15$ days, and a deficit of response in the core. This is suggestive of a disk geometry of the broad-line region (BLR). The relatively weaker H$\beta$ response at the longer lags in the red side indicates anisotropic emission from the outer part of the BLR. We also recover the velocity-delay maps of NGC 5548 from the historical data of 13 years to investigate the long-term variability of its BLR. In general, the BLR of NGC 5548 was switching between the inflow and virialized phases in the past years. The resultant maps of seven years reveal inflow signatures and show decreasing lags, indicating that the changes in the BLR size are related to the infalling BLR gas. The other four maps show potential disk signatures which are similar to our map.

As one of the series of papers reporting on a large reverberation mapping campaign, we apply the maximum entropy method (MEM) to 9 narrow-line Seyfert 1 galaxies with super-Eddington accretion rates observed during 2012-2013 for the velocity-delay maps of their H\beta\ and H\gamma\ emission lines. The maps of 6 objects are reliably reconstructed using MEM. The maps of H\beta\ and H\gamma\ emission lines of Mrk 335 indicate that the gas of its broad-line region (BLR) is infalling. For Mrk 142, its H\beta\ and H\gamma\ lines show signatures of outflow. The H\beta\ and H\gamma\ maps of Mrk 1044 demonstrate complex kinematics -- a virialized motion accompanied by an outflow signature, and the H\beta\ map of IRAS F12397+3333 is consistent with a disk or a spherical shell. The H\beta\ maps of Mrk 486 and MCG +06-26-012 suggest the presence of an inflow and outflow, respectively. These super-Eddington accretors show diverse geometry and kinematics. Brief discussions of their BLRs are provided for each individual object.

Narrow-line Seyfert 1 galaxies are defined on the basis of their line widths, and they are generally considered to be high Eddington ratio sources. But in the context of the Quasar Main Sequence, high Eddington rate sources are those which have weak [O III] lines and strong Fe II lines. There is an overlap between the two populations, but they are not identical. Thus these two selection criteria give a different view on which objects are actually high Eddington ratio sources. We discuss this issue in the context of the broad band spectral energy density, emission line shape modeling, Fe II pseudo-continuum strength, and the level of X-ray variability. We also discuss the issue of the viewing angle and the insight one can gain from spectropolarimetric observations. We conclude that there is no single driver behind the Quasar Main Sequence, and the Eddington ratio of the source cannot be determined from the source location on the optical plane alone. On the other hand, an expected range of AGN parameters combined with a simple model of the Fe II production represent well the observed coverage pattern of the plane, with not much effect needed from the dispersion due to the viewing angle.

As one of the series of papers reporting on a large reverberation mapping campaign of super-Eddington accreting massive black holes (SEAMBHs) in active galactic nuclei (AGNs), we present the results of 10 SEAMBHs monitored spectroscopically during 2015-2017. Six of them are observed for the first time, and have generally higher 5100 \AA\ luminosities than the SEAMBHs monitored in our campaign from 2012 to 2015; the remaining four are repeat observations to check if their previous lags change. Similar to the previous SEAMBHs, the H$\beta$ time lags of the newly observed objects are shorter than the values predicted by the canonical $R_{\mathrm{H\beta}}$-$L_{5100}$ relation of sub-Eddington AGNs, by factors of $\sim2-6$, depending on the accretion rate. The four previously observed objects have lags consistent with previous measurements. We provide linear regressions for the $R_{\mathrm{H\beta}}$-$L_{5100}$ relation, solely for the SEAMBH sample and for low-accretion AGNs. We find that the relative strength of Fe II and the profile of the H$\beta$ emission line can be used as proxies of accretion rate, showing that the shortening of H$\beta$ lags depends on accretion rates. The recent SDSS-RM discovery of shortened H$\beta$ lags in AGNs with low accretion rates provides compelling evidence for retrograde accretion onto the black hole. These evidences show that the canonical $R_{\mathrm{H\beta}}$-$L_{5100}$ relation holds only in AGNs with moderate accretion rates. At low accretion rates, it should be revised to include the effects of black hole spin, whereas the accretion rate itself becomes a key factor in the regime of high accretion rates.

Optical and near-infrared photometry, optical spectroscopy, and soft X-ray and UV monitoring of the changing-look active galactic nucleus NGC 2617 show that it continues to have the appearance of a type-1 Seyfert galaxy. An optical light curve for 2010$-$2017 indicates that the change of type probably occurred between 2010 October and 2012 February and was not related to the brightening in 2013. In 2016 and 2017 NGC 2617 brightened again to a level of activity close to that in 2013 April. However, in 2017 from the end of the March to end of July 2017 it was in very low level and starting to change back to a Seyfert 1.8. We find variations in all passbands and in both the intensities and profiles of the broad Balmer lines. A new displaced emission peak has appeared in H$\beta$. X-ray variations are well correlated with UV$-$optical variability and possibly lead by $\sim$ 2$-$3 d. The $K$ band lags the $J$ band by about 21.5 $\pm$ 2.5 d and lags the combined $B + J$ bands by $\sim$ 25 d. $J$ lags $B$ by $\sim$ 3 d. This could be because $J$-band variability arises predominantly from the outer part of the accretion disc, while $K$-band variability is dominated by thermal re-emission by dust. We propose that spectral-type changes are a result of increasing central luminosity causing sublimation of the innermost dust in the hollow bi-conical outflow. We briefly discuss various other possible reasons that might explain the dramatic changes in NGC 2617.

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.

The hidden broad-line regions (BLRs) in Seyfert 2 galaxies, which display broad emission lines (BELs) in their polarized spectra, are a key piece of evidence in support of the unified model for active galactic nuclei (AGNs). However, the detailed kinematics and geometry of hidden BLRs are still not fully understood. The virial factor obtained from reverberation mapping of type 1 AGNs may be a useful diagnostic of the nature of hidden BLRs in type 2 objects. In order to understand the hidden BLRs, we compile six type 2 objects from the literature with polarized BELs and dynamical measurements of black hole masses. All of them contain pseudobulges. We estimate their virial factors, and find the average value is 0.60 and the standard deviation is 0.69, which agree well with the value of type 1 AGNs with pseudobulges. This study demonstrates that (1) the geometry and kinematics of BLR are similar in type 1 and type 2 AGNs of the same bulge type (pseudobulges), and (2) the small values of virial factors in Seyfert 2 galaxies suggest that, similar to type 1 AGNs, BLRs tend to be very thick disks in type 2 objects.

In this work, which is a continuation of Castell\'o-Mor et al. (2016), we present new X-ray and infrared (IR) data for a sample of active galactic nuclei (AGN) covering a wide range in Eddington ratio over a small luminosity range. In particular, we rigorously explore the dependence of the optical-to-X-ray spectral index $\alpha_{OX}$ and the IR-to-optical spectral index on the dimensionless accretion rate, $\dot{\mathcal{M}}=\dot{m}/\eta$ where $\dot{m}=L_{AGN}/L_{Edd}$ and $\eta$ is the mass-to-radiation conversion efficiency, in low and high accretion rate sources. We find that the SED of the faster accreting sources are surprisingly similar to those from the comparison sample of sources with lower accretion rate. In particular: i) the optical-to-UV AGN SED of slow and fast accreting AGN can be fitted with thin AD models. ii) The value of $\alpha_{OX}$ is very similar in slow and fast accreting systems up to a dimensionless accretion rate $\dot{\mathcal{M}}_{c}\sim$10. We only find a correlation between $\alpha_{OX}$ and $\dot{\mathcal{M}}$ for sources with $\dot{\mathcal{M}} > \dot{\mathcal{M}}_{c}$. In such cases, the faster accreting sources appear to have systematically larger $\alpha_{OX}$ values. iii) We also find that the torus in the faster accreting systems seems to be less efficient in reprocessing the primary AGN radiation having lower IR-to-optical spectral slopes. These findings, failing to recover the predicted differences between the SEDs of slim and thin ADs within the observed spectral window, suggest that additional physical processes or very special geometry act to reduce the extreme UV radiation in fast accreting AGN. This may be related to photon trapping, strong winds, and perhaps other yet unknown physical processes.

The origin of the Broad Line Region (BLR) in active galaxies remains unknown. It seems to be related to the underlying accretion disk but an efficient mechanism is required to rise the material from the disk surface without giving too strong signatures of the outflow in the case of the low ionization lines. We discuss in detail two proposed mechanisms: (i) radiation pressure acting on dust in the disk atmosphere creating a failed wind (ii) the gravitational instability of the underlying disk. We compare the predicted location of the inner radius of the BLR in those two scenarios with the observed position obtained from the reverberation studies of several active galaxies. The failed dusty outflow model well represents the observational data while the predictions of the self-gravitational instability are not consistent with observations. The issue remains why actually we do not see any imprints of the underlying disk instability in the BLR properties.

We propose a method for the flux calibration of reverberation mapping spectra based on accurate measurement of [O III] $\lambda 5007$ emission by spectral fitting. The method can achieve better accuracy than the traditional method of van Groningen & Wanders (1992), allowing reverberation mapping measurements for object with variability amplitudes as low as $\sim$ 5%. As a demonstration, we reanalyze the data of the Seyfert 1 galaxy MCG--6-30-15 taken from the 2008 campaign of the Lick AGN Monitoring Project, which previously failed to obtain a time lag for this weakly variable object owing to a relatively large flux calibration uncertainty. We detect a statistically significant rest-frame time lag of $6.38_{-2.69}^{+3.07}$ days between the H$\beta$ and $V$-band light curves. Combining this lag with FWHM(H$\beta$) = $1933\pm81$ $\rm km~s^{-1}$ and a virial coefficient of $f$ = 0.7, we derive a virial black hole mass of $3.26_{-1.40}^{+1.59}\times10^6$ $M_{\odot}$, which agrees well with previous estimates by other methods.

NGC 5548 is the best-observed reverberation-mapped active galactic nucleus with long-term, intensive monitoring. Here we report results from a new observational campaign between January and July, 2015. We measure the centroid time lag of the broad H$\beta$ emission line with respect to the 5100 \AA continuum and obtain $\tau_{\rm cent} = 7.20^{+1.33}_{-0.35}$ days in the rest frame. This yields a black hole mass of $M_{\bullet}=8.71^{+3.21}_{-2.61} $x$ 10^{7}M_{\odot}$ using a broad H$\beta$ line dispersion of $3124\pm302$ km s$^{-1}$ and a virial factor of $f_{_{\rm BLR}}=6.3\pm1.5$ for the broad-line region (BLR), consistent with the mass measurements from previous H$\beta$ campaigns. The high-quality data allow us to construct a velocity-binned delay map for the broad H$\beta$ line, which shows a symmetric response pattern around the line center, a plausible kinematic signature of virialized motion of the BLR. Combining all the available measurements of H$\beta$ time lags and the associated mean 5100 {\AA} luminosities over 18 campaigns between 1989 and 2015, we find that the H$\beta$ BLR size varies with the mean optical luminosity, but, interestingly, with a possible delay of $2.35_{-1.25}^{+3.47}$ yrs. This delay coincides with the typical BLR dynamical timescale of NGC 5548, indicating that the BLR undergoes dynamical changes, possibly driven by radiation pressure.

This paper reports results of the third-year campaign of monitoring super-Eddington accreting massive black holes (SEAMBHs) in active galactic nuclei (AGNs) between 2014-2015. Ten new targets were selected from quasar sample of Sloan Digital Sky Survey (SDSS), which are generally more luminous than the SEAMBH candidates in last two years. H$\beta$ lags ($\tau_{_{\rm H\beta}}$) in five of the 10 quasars have been successfully measured in this monitoring season. We find that the lags are generally shorter, by large factors, than those of objects with same optical luminosity, in light of the well-known $R_{_{\rm H\beta}}-L_{5100}$ relation. The five quasars have dimensionless accretion rates of $\dot{\mathscr{M}}=10-10^3$. Combining measurements of the previous SEAMBHs, we find that the reduction of H$\beta$ lags tightly depends on accretion rates, $\tau_{_{\rm H\beta}}/\tau_{_{R-L}}\propto\dot{\mathscr{M}}^{-0.42}$, where $\tau_{_{R-L}}$ is the H$\beta$ lag from the normal $R_{_{\rm H\beta}}-L_{5100}$ relation. Fitting 63 mapped AGNs, we present a new scaling relation for the broad-line region: $R_{_{\rm H\beta}}=\alpha_1\ell_{44}^{\beta_1}\,\min\left[1,\left(\dot{\mathscr{M}}/\dot{\mathscr{M}}_c\right)^{-\gamma_1}\right]$, where $\ell_{44}=L_{5100}/10^{44}\,\rm erg~s^{-1}$ is 5100 \AA\ continuum luminosity, and coefficients of $\alpha_1=(29.6_{-2.8}^{+2.7})$ lt-d, $\beta_1=0.56_{-0.03}^{+0.03}$, $\gamma_1=0.52_{-0.16}^{+0.33}$ and $\dot{\mathscr{M}}_c=11.19_{-6.22}^{+2.29}$. This relation is applicable to AGNs over a wide range of accretion rates, from $10^{-3}$ to $10^3$. Implications of this new relation are briefly discussed.

As a natural consequence of cosmological hierarchical structure formation, sub-parsec supermassive black hole binaries (SMBHBs) should be common in galaxies but thus far have eluded spectroscopic identification. Based on four decades of optical spectroscopic monitoring, we report that the nucleus of NGC 5548, a nearby Seyfert galaxy long suspected to have experienced a major merger about one billion years ago, exhibits long-term variability with a period of 14 years in the optical continuum and broad Hbeta emission line. Remarkably, the double-peaked profile of Hbeta shows systematic velocity changes with a similar period. These pieces of observations plausibly indicate that a SMBHB resides in the center of NGC 5548. The complex, secular variations in the line profiles can be explained by orbital motion of a binary with equal mass and a semi-major axis of 22 light-days (corresponding to 18 milli-parsec). At a distance of 75 Mpc, NGC 5548 is one of the nearest sub-parsec SMBHB candidates that offers an ideal laboratory for gravitational wave detection.

In the sixth of the series of papers reporting on a large reverberation mapping (RM) campaign of active galactic nuclei (AGNs) with high accretion rates, we present velocity-resolved time lags of H$\beta$ emission lines for nine objects observed in the campaign during 2012$-$2013. In order to correct the line-broadening caused by seeing and instruments before the analysis of velocity-resolved RM, we adopt Richardson-Lucy deconvolution to reconstruct their H$\beta$ profiles. The validity and effectiveness of the deconvolution are checked out by Monte Carlo simulation. Five among the nine objects show clear dependence of time delay on velocity. Mrk 335 and Mrk 486 show signatures of gas inflow whereas the clouds in the broad-line regions (BLRs) of Mrk 142 and MCG +06-26-012 tend to be radial outflowing. Mrk 1044 is consistent with the case of virialized motions. The lags of the rest four are not velocity-resolvable. The velocity-resolved RM of super-Eddington accreting massive black holes (SEAMBHs) shows that they have diversity of the kinematics in their BLRs. Comparing with the AGNs with sub-Eddington accretion rates, we do not find significant differences in the BLR kinematics of SEAMBHs.

Broad emission lines in active galactic nuclei (AGNs) mainly arise from gas photoionized by continuum radiation from an accretion disk around a central black hole. The shape of the broad-line profile, described by ${\cal D}_{_{\rm H\beta}}={\rm FWHM}/\sigma_{_{\rm H\beta}}$, the ratio of full width at half maximum to the dispersion of broad H$\beta$, reflects the dynamics of the broad-line region (BLR) and correlates with the dimensionless accretion rate ($\dot{\mathscr{M}}$) or Eddington ratio ($L_{\rm bol}/L_{\rm Edd}$). At the same time, $\dot{\mathscr{M}}$ and $L_{\rm bol}/L_{\rm Edd}$ correlate with ${\cal R}_{\rm Fe}$, the ratio of optical Fe II to H$\beta$ line flux emission. Assembling all AGNs with reverberation mapping measurements of broad H$\beta$, both from the literature and from new observations reported here, we find a strong bivariate correlation of the form $\log(\dot{\mathscr{M}},L_{\rm bol}/L_{\rm Edd})=\alpha+\beta{\cal D}_{_{\rm H\beta}}+\gamma{\cal R}_{\rm Fe},$ where $\alpha=(2.47,0.31)$, $\beta=-(1.59,0.82)$ and $\gamma=(1.34,0.80)$. We refer to this as the fundamental plane of the BLR. We apply the plane to a sample of $z < 0.8$ quasars to demonstrate the prevalence of super-Eddington accreting AGNs are quite common at low redshifts.

We have completed two years of photometric and spectroscopic monitoring of a large number of active galactic nuclei (AGNs) with very high accretion rates. In this paper, we report on the result of the second phase of the campaign, during 2013--2014, and the measurements of five new H$\beta$ time lags out of eight monitored AGNs. All five objects were identified as super-Eddington accreting massive black holes (SEAMBHs). The highest measured accretion rates for the objects in this campaign are $\dot{\mathscr{M}}\gtrsim 200$, where $\dot{\mathscr{M}}= \dot{M}_{\bullet}/L_{\rm Edd}c^{-2}$, $\dot{M}_{\bullet}$ is the mass accretion rates, $L_{\rm Edd}$ is the Eddington luminosity and $c$ is the speed of light. We find that the H$\beta$ time lags in SEAMBHs are significantly shorter than those measured in sub-Eddington AGNs, and the deviations increase with increasing accretion rates. Thus, the relationship between broad-line region size ($R_{_{\rm H\beta}}$) and optical luminosity at 5100\AA, $R_{_{\rm H\beta}}-L_{5100}$, requires accretion rate as an additional parameter. We propose that much of the effect may be due to the strong anisotropy of the emitted slim-disk radiation. Scaling $R_{_{\rm H\beta}}$ by the gravitational radius of the black hole, we define a new radius-mass parameter ($Y$) and show that it saturates at a critical accretion rate of $\dot{\mathscr{M}}_c=6\sim 30$, indicating a transition from thin to slim accretion disk and a saturated luminosity of the slim disks. The parameter $Y$ is a very useful probe for understanding the various types of accretion onto massive black holes. We briefly comment on implications to the general population of super-Eddington AGNs in the universe and applications to cosmology.

This is the third in a series of papers reporting on a large reverberation-mapping campaign aimed to study the properties of active galactic nuclei (AGNs) with high accretion rates. We present new results on the variability of the optical Fe II emission lines in 10 AGNs observed by the Yunnan Observatory 2.4m telescope during 2012--2013. We detect statistically significant time lags, relative to the AGN continuum, in nine of the sources. This accurate measurement is achieved by using a sophisticated spectral fitting scheme that allows for apparent flux variations of the host galaxy, and several narrow lines, due to the changing observing conditions. Six of the newly detected lags are indistinguishable from the Hbeta lags measured in the same sources. Two are significantly longer and one is slightly shorter. Combining with Fe II lags reported in previous studies, we find a Fe II radius--luminosity relationship similar to the one for Hbeta, although our sample by itself shows no clear correlation. The results support the idea that Fe II emission lines originate in photoionized gas which, for the majority of the newly reported objects, is indistinguishable from the Hbeta-emitting gas. We also present a tentative correlation between the lag and intensity of Fe II and Hbeta and comment on its possible origin.

Supermassive black holes in active galactic nuclei (AGNs) undergo a wide range of accretion rates, which lead to diversity of appearance. We consider the effects of anisotropic radiation from accretion disks on the broad-line region (BLR), from the Shakura-Sunyaev regime to slim disks with super-Eddington accretion rates. The geometrically thick funnel of the inner region of slim disks produces strong self-shadowing effects that lead to very strong anisotropy of the radiation field. We demonstrate that the degree of anisotropy of the radiation fields grows with increasing accretion rate. As a result of this anisotropy, BLR clouds receive different spectral energy distributions depending on their location relative to the disk, resulting in diverse observational appearance of the BLR. We show that the self-shadowing of the inner parts of the disk naturally produces two dynamically distinct regions of the BLR, depending on accretion rate. These two regions manifest themselves as kinematically distinct components of the broad H$\beta$ line profile with different line widths and fluxes, which jointly account for the Lorentzian profile generally observed in narrow-line Seyfert 1 galaxies. In the time domain, these two components are expected reverberate with different time lags with respect to the varying ionizing continuum, depending on the accretion rate and the viewing angle of the observer. The diverse appearance of the BLR due to the anisotropic ionizing energy source can be tested by reverberation mapping of H$\beta$ and other broad emission lines (e.g., \feii), providing a new tool to diagnose the structure and dynamics of the BLR. Other observational consequences of our model are also explored.

A tight relation between the size of the broad-line region (BLR) and optical luminosity has been established in about 50 active galactic nuclei studied through reverberation mapping of the broad Hbeta emission line. The R_blr-L relation arises from simple photoionization considerations. Using a general relativistic model of an optically thick, geometrically thin accretion disk, we show that the ionizing luminosity jointly depends on black hole mass, accretion rate, and spin. The non-monotonic relation between the ionizing and optical luminosity gives rise to a complicated relation between the BLR size and the optical luminosity. We show that the reverberation lag of Hbeta to the varying continuum depends very sensitively to black hole spin. For retrograde spins, the disk is so cold that there is a deficit of ionizing photons in the BLR, resulting in shrinkage of the hydrogen ionization front with increasing optical luminosity, and hence shortened Hbeta lags. This effect is specially striking for luminous quasars undergoing retrograde accretion, manifesting in strong deviations from the canonical R_blr-L relation. This could lead to a method to estimate black hole spins of quasars and to study their cosmic evolution. At the same time, the small scatter of the observed R_blr-L relation for the current sample of reverberation-mapped active galaxies implies that the majority of these sources have rapidly spinning black holes.

This is the second in a series of papers reporting on a large reverberation mapping (RM) campaign to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). The goal is to identify super-Eddington accreting massive black holes (SEAMBHs) and to use their unique properties to construct a new method for measuring cosmological distances. Based on theoretical models, the saturated bolometric luminosity of such sources is proportional to the BH mass which can be used to obtain their distance. Here we report on five new RM measurements and show that in four of the cases we can measure the BH mass and three of these sources are SEAMBHs. Together with the three sources from our earlier work, we now have six new sources of this type. We use a novel method based on a minimal radiation efficiency to identify nine additional SEAMBHs from earlier RM-based mass measurements. We use a Bayesian analysis to determine the parameters of the new distance expression, and the method uncertainties, from the observed properties of the objects in the sample. The ratio of the newly measured distances to the standard cosmological ones has a mean scatter of 0.14 dex, indicating that SEAMBHs can be use as cosmological distance probes. With their high luminosity, long period of activity and large numbers at high redshifts, SEAMBHs have a potential to extend the cosmic distance ladder beyond the range now explored by type Ia supernovae.

Flux calibration of spectra in reverberation mapping (RM) is most often performed by assuming the flux constancy of some specified narrow emission lines, which stem from an extended region that is sometimes partially spatially resolved, in contrast to the point-like broad-line region and the central continuum source. The inhomogeneous aperture geometries used among different observation sets in a joint monitoring campaign introduce systematic deviations to the fluxes of broad lines and central continuum, and intercalibration over these data sets is required. As an improvement to the previous empirical correction performed by comparing the (nearly) contemporaneous observation points, we describe a feasible Bayesian method that obviates the need for (nearly) contemporaneous observations, naturally incorporates physical models of flux variations, and fully takes into account the measurement errors. In particular, it fits all the data sets simultaneously regardless of samplings and makes use of all of the information in the data sets. A Markov Chain Monte Carlo implementation is employed to recover the parameters and uncertainties for intercalibration. Application to the RM data sets of NGC 5548 with joint monitoring shows the high fidelity of our method.

This is the first paper in a series devoted to systematic study of the size and structure of the broad-line region (BLR) in active galactic nuclei (AGNs) using reverberation mapping (RM) data. We employ a recently developed Bayesian approach that statistically describes the variabibility as a damped random walk process and delineates the BLR structure using a flexible disk geometry that can account for a variety of shapes, including disks, rings, shells, and spheres. We allow for the possibility that the line emission may respond non-linearly to the continuum, and we detrend the light curves when there is clear evidence for secular variation. We use a Markov Chain Monte Carlo implementation based on Bayesian statistics to recover the parameters and uncertainties for the BLR model. The corresponding transfer function is obtained self-consistently. We tentatively constrain the virial factor used to estimate black hole masses; more accurate determinations will have to await velocity-resolved RM data. Application of our method to RM data with Hbeta monitoring for about 40 objects shows that the assumed BLR geometry can reproduce quite well the observed emission-line fluxes from the continuum light curves. We find that the Hbeta BLR sizes obtained from our method are on average ~20% larger than those derived from the traditional cross-correlation method. Nevertheless, we still find a tight BLR size-luminosity relation with a slope of alpha=0.55\pm0.03 and an intrinsic scatter of ~0.18 dex. In particular, we demonstrate that our approach yields appropriate BLR sizes for some objects (such as Mrk 142 and PG 2130+099) where traditional methods previously encountered difficulties.

The metallicity of active galactic nuclei (AGNs), which can be measured by emission line ratios in their broad and narrow line regions (BLRs and NLRs), provides invaluable information about the physical connection between the different components of AGNs. From the archival databases of the International Ultraviolet Explorer, the Hubble Space Telescope and the Sloan Digital Sky Survey, we have assembled the largest sample available of AGNs which have adequate spectra in both the optical and ultraviolet bands to measure the narrow line ratio [N II]/H{\alpha} and also, in the same objects, the broad-line N V/C IV ratio. These permit the measurement of the metallicities in the NLRs and BLRs in the same objects. We find that neither the BLR nor the NLR metallicity correlate with black hole masses or Eddington ratios, but there is a strong correlation between NLR and BLR metallicities. This metallicity correlation implies that outflows from BLRs carry metal-rich gas to NLRs at characteristic radial distances of ~ 1.0 kiloparsec. This chemical connection provides evidence for a kinetic feedback of the outflows to their hosts. Metals transported into the NLR enhance the cooling of the ISM in this region, leading to local star formation after the AGNs turn to narrow line LINERs. This post-AGN star formation is predicted to be observable as an excess continuum emission from the host galaxies in the near infrared and ultraviolet, which needs to be further explored.

We report first results from a large project to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). Such objects may be different from other AGNs in being powered by slim accretion disks and showing saturated accretion luminosities, but both are not yet fully understood. The results are part of a large reverberation mapping (RM) campaign using the 2.4-m Shangri-La telescope at the Yunnan Observatory in China. The goals are to investigate the gas distribution near the BH and the properties of the central accretion disks, to measure BH mass and Eddington ratios, and to test the feasibility of using such objects as a new type of cosmological candles. The paper presents results for three objects, Mrk 335, Mrk 142 and IRAS F12397+3333 with H$\beta$ time lags relative to the 5100\AA\ continuum of $10.6^{+1.7}_{-2.9}$, $6.4^{+0.8}_{-2.2}$ and $11.4^{+2.9}_{-1.9}$ days, respectively. The corresponding BH masses are $(8.3_{-3.2}^{+2.6})\times 10^6M_{\odot}$, $(3.4_{-1.2}^{+0.5})\times 10^6M_{\odot}$ and $(7.5_{-4.1}^{+4.3})\times 10^6M_{\odot}$, and the lower limits on the Eddington ratios 0.6, 2.3, and 4.6 for the minimal radiative efficiency of 0.038. Mrk 142 and IRAS F12397+333 (extinction corrected) clearly deviate from the currently known relation between H$\beta$ lag and continuum luminosity. The three Eddington ratios are beyond the values expected in thin accretion disks and two of them are the largest measured so far among objects with RM-based BH masses. We briefly discuss implications for slim disks, BH growth and cosmology.

We present a spectroscopic catalog of 67082 M dwarfs from the LAMOST pilot survey. For each spectrum of the catalog, spectral subtype, radial velocity, equivalent width of H${\alpha}$, a number of prominent molecular band indices and the metal sensitive parameter $\zeta$ are provided . Spectral subtype have been estimated by a remedied Hammer program (Original Hammer: Covey et al. 2007), in which indices are reselected to obtain more accurate auto-classified spectral subtypes. All spectra in this catalog have been visually inspected to confirm the spectral subtypes. Radial velocities have been well measured by our developed program which uses cross-correlation method and estimates uncertainty of radial velocity as well. We also examine the magnetic activity properties of M dwarfs traced by H${\alpha}$ emission line. The molecular band indices included in this catalog are temperature or metallicity sensitive and can be used for future analysis of the physical properties of M dwarfs. The catalog is available on the website \url{http://sciwiki.lamost.org/MCatalogPilot/}.

Super-Eddington accreting massive black holes (SEAMBHs) reach saturated luminosities above a certain accretion rate due to photon trapping and advection in slim accretion disks. We show that these SEAMBHs could provide a new tool for estimating cosmological distances if they are properly identified by hard X-ray observations, in particular by the slope of their 2-10 keV continuum. To verify this idea we obtained black hole mass estimates and X-ray data for a sample of 60 narrow line Seyfert 1 galaxies that we consider to be the most promising SEAMBH candidates. We demonstrate that the distances derived by the new method for the objects in the sample get closer to the standard luminosity distances as the hard X-ray continuum gets steeper. The results allow us to analyze the requirements for using the method in future samples of active black holes and to demonstrate that the expected uncertainty, given large enough samples, can make them into a useful, new cosmological ruler.

We present the discovery of 6 nebular objects made by amateur astronomers. Four of these discoveries are clearly planetary nebulae (PNe), one is a possible PN, and another is a likely H II region. The bipolar nebula Ou4 presents the largest angular extent ever found : over one degree on the sky! We consider various scenarios that could explain such a nebula. Ou4 could be one of the nearest PNe known, though its possible PN nature will need confirmation.

Recent x-ray observations of Mrk 766 suggest that broad emission line region clouds cross our line of sight and produce variable x-ray absorption. Here we investigate what optical/ultraviolet spectroscopic features would be produced by such "Intervening BLR Clouds" (IBC) crossing our line of sight to the accretion disk, the source of the optical/UV continuum. Although the emission spectrum produced by intervening clouds is identical to the standard BLR model, they may produce absorption features on the optical or UV continuum. Single clouds will have little effect on the optical/UV spectrum because BLR clouds are likely to be much smaller than the accretion disk. This is unlike the X-ray case, where the radiation source is considerably smaller. However, an ensemble of intervening clouds will produce spectroscopic features in the FUV including a strong depression between the Lyman limit and Ly{\alpha}. The amount of the depression will indicate the line-of-sight covering factor of clouds, an unknown quantity that is important for the ionization of the intergalactic medium and the energy budget of AGN. Comparison with observations suggests that the SED of Mrk 766 may be affected by intervening BLR clouds and IBC may exist in most of AGNs.

(abridged) We study the consequence of star formation (SF) in an self-gravity dominated accretion disk in quasars. The warm skins of the SF disk are governed by the radiation from the inner part of the accretion disk to form Compton atmosphere (CAS). The CAS are undergoing four phases to form broad line regions. Phase I is the duration of pure accumulation supplied by the SF disk. During phase II clouds begin to form due to line cooling and sink to the SF disk. Phase III is a period of preventing clouds from sinking to the SF disk through dynamic interaction between clouds and the CAS. Finally, phase IV is an inevitable collapse of the entire CAS through line cooling. This CAS evolution drives the episodic appearance of BLRs. Geometry and dynamics of BLRs can be self-consistently derived from the thermal instability of the CAS during phases II and III by linear analysis. The metallicity gradient of SF disk gives rise to different properties of clouds from outer to inner part of BLRs. We find that clouds have column density N_H < 10^22cm^{-2} in the metal-rich regions whereas they have N_H > 10^22 cm^{-2} in the metal-poor regions. The metal-rich clouds compose the high ionization line (HIL) regions whereas the metal-poor clouds are in low ionization line (LIL) regions. Metal-rich clouds in HIL regions will be blown away by radiation pressure, forming the observed outflows. The LIL regions are episodic due to the mass cycle of clouds with the CAS in response to continuous injection by the SF disk, giving rise to different types of AGNs. Based on SDSS quasar spectra, we identify a spectral sequence in light of emission line equivalent width from Phase I to IV. A key phase in the episodic appearance of the BLRs is bright type II AGNs with no or only weak BLRs. We discuss observational implications and tests of the theoretical predictions of this model.