Source author record

Yan-Rong Li

Yan-Rong Li appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

35works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2024arXiv

AGN STORM 2: V. Anomalous Behavior of the CIV Light Curve in Mrk 817

An intensive reverberation mapping campaign on the Seyfert 1 galaxy Mrk817 using the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST) revealed significant variations in the response of the broad UV emission lines to fluctuations in the continuum emission. The response of the prominent UV emission lines changes over a $\sim$60-day duration, resulting in distinctly different time lags in the various segments of the light curve over the 14 months observing campaign. One-dimensional echo-mapping models fit these variations if a slowly varying background is included for each emission line. These variations are more evident in the CIV light curve, which is the line least affected by intrinsic absorption in Mrk817 and least blended with neighboring emission lines. We identify five temporal windows with distinct emission line response, and measure their corresponding time delays, which range from 2 to 13 days. These temporal windows are plausibly linked to changes in the UV and X-ray obscuration occurring during these same intervals. The shortest time lags occur during periods with diminishing obscuration, whereas the longest lags occur during periods with rising obscuration. We propose that the obscuring outflow shields the ultraviolet broad lines from the ionizing continuum. The resulting change in the spectral energy distribution of the ionizing continuum, as seen by clouds at a range of distances from the nucleus, is responsible for the changes in the line response.

preprint2022arXiv

Accretion Disk Size Measurements of Active Galactic Nuclei Monitored by the Zwicky Transient Facility

We compile a sample of 92 active galactic nuclei (AGNs) at z<0.75 with $gri$ photometric light curves from the archival data of the Zwicky Transient Facility and measure the accretion disk sizes via continuum reverberation mapping. We employ Monte Carlo simulation tests to assess the influences of data sampling and broad emission lines and select out the sample with adequately high sampling cadences (3 days apart in average) and minimum contaminations of broad emission lines. The inter-band time delays of individual AGNs are calculated using the interpolated cross-correlation function and then these delays are fitted with a generalized accretion disk model, in which inter-band time delays are a power function of wavelength, black hole mass, and luminosity. A Markov-chain Monte Carlo method is adopted to determine the best parameter values. Overall the inter-band time delays can be fitted with the $τ\ \propto λ^{4/3}$ relation as predicted from a steady-state, optically thick, geometrically thin accretion disk, however, the yielded disk size is systematically larger than expected, although the ratio of the measured to theoretical disk sizes depend on using the emissivity -- or responsivity -- weighted disk radius. These results are broadly consistent with previous studies, all together raising a puzzle about the "standard" accretion disk model.

preprint2022arXiv

Monitoring AGNs with H$β$ Asymmetry. III. Long-term Reverberation Mapping Results of 15 Palomar-Green Quasars

In this third paper of the series reporting on the reverberation mapping (RM) campaign of active galactic nuclei with asymmetric H$β$ 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$β$ 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.

preprint2022arXiv

Revisiting the Continuum Reverberation Lags in the AGN PKS 0558-504

We present a revised analysis of the photometric reverberation mapping campaign of the narrow-line Seyfert 1 galaxy PKS 0558-504 carried out with the Swift Observatory during 2008--2010. Previously, Gliozzi et al.\ found using the Discrete Correlation Function (DCF) method that the short-wavelength continuum variations lagged behind variations at longer wavelengths, the opposite of the trend expected for thermal reprocessing of X-rays by the accretion disc, and they interpreted their results as evidence against the reprocessing model. We carried out new DCF measurements that demonstrate that the inverted lag-wavelength relationship found by Gliozzi et al.\ resulted from their having interchanged the order of the driving and responding light curves when measuring the lags. To determine the inter-band lags and uncertainties more accurately, we carried out new measurements with four independent methods. These give consistent results showing time delays increasing as a function of wavelength, as expected for the disc reprocessing scenario. The slope of the re-analysed delay spectrum appears to be roughly compatible with the predicted $τ\propto λ^{4/3}$ relationship for reprocessing by an optically thick and geometrically thin accretion disc, although the data points exhibit a large scatter about the fitted power-law trend.

preprint2022arXiv

Spectroastrometry and Reverberation Mapping: the Mass and Geometric Distance of the Supermassive Black Hole in the Quasar 3C 273

The quasar 3C 273 has been observed with infrared spectroastrometry (SA) on broad Pa$α$ line and optical reverberation mapping (RM) on broad H$β$ line. SA delivers information about the angular size and structure of the Pa$α$ broad-line region (BLR), while RM delivers information about the physical size and structure of the H$β$ 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.

preprint2020arXiv

A parallax distance to 3C 273 through spectroastrometry and reverberation mapping

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.

preprint2020arXiv

Broad-line Region of the Quasar PG 2130+099 from a Two-Year Reverberation Mapping Campaign with High Cadence

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$β$ 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$β$, He II, He I, and Fe II lines with respect to the varying 5100Å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.

preprint2020arXiv

Constraints on individual supermassive binary black holes using observations of PSR J1909-3744

We perform a search for gravitational waves (GWs) from several supermassive binary black hole (SMBBH) candidates (NGC 5548, Mrk 231, OJ 287, PG 1302-102, NGC 4151, Ark 120 and 3C 66B) in long-term timing observations of the pulsar PSR J1909$-$3744 obtained using the Parkes radio telescope. No statistically significant signals were found. We constrain the chirp masses of those SMBBH candidates and find the chirp mass of NGC 5548 and 3C 66B to be less than $2.4 \times 10^9\,\rm M_{\odot}$ and $2.5 \times 10^9\,\rm M_{\odot}$ (with 95% confidence), respectively. Our upper limits remain a factor of 3 to 370 above the likely chirp masses for these candidates as estimated from other approaches. The observations processed here provide upper limits on the GW strain amplitude that improve upon the results from the first Parkes Pulsar Timing Array data release by a factor of 2 to 7. We investigate how information about the orbital parameters can help improve the search sensitivity for individual SMBBH systems. Finally, we show that these limits are insensitive to uncertainties in the Solar System ephemeris model.

preprint2020arXiv

Supermassive black holes with high accretion rates in active galactic nuclei. XI. Accretion disk reverberation mapping of Mrk 142

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, $τ(λ)$, generally follow the relation $τ(λ) \propto λ^{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 $τ(λ) \propto λ^{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_ν\proptoν^{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.

preprint2020arXiv

Untangling Optical Emissions of the Jet and Accretion Disk in the Flat-Spectrum Radio Quasar 3C 273 with Reverberation Mapping Data

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.

preprint2019arXiv

Active Galactic Nuclei with Ultra-fast Outflows Monitoring Project: The Broad-line Region of Mrk 79 as a Disk Wind

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$β~\lambda4861$, H$γ~\lambda4340$, He II $\lambda4686$ and He I $\lambda5876$ with respect to the varying AGN continuum. Based on the broad H$β~\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$γ~\lambda4340$, H$β~\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.

preprint2016arXiv

A Non-parametric Approach to Constrain the Transfer Function in Reverberation Mapping

Broad emission lines of active galactic nuclei stem from a spatially extended region (broad-line region, BLR) that is composed of discrete clouds and photoionized by the central ionizing continuum. The temporal behaviors of these emission lines are blurred echoes of the continuum variations (i.e., reverberation mapping, RM) and directly reflect the structures and kinematic information of BLRs through the so-called transfer function (also known as the velocity-delay map). Based on the previous works of Rybicki & Press (1992) and Zu et al. (2011), we develop an extended, non-parametric approach to determine the transfer function for RM data, in which the transfer function is expressed as a sum of a family of relatively displaced Gaussian response functions. Therefore, arbitrary shapes of transfer functions associated with complicated BLR geometry can be seamlessly included, enabling us to relax the presumption of a specified transfer function frequently adopted in previous studies and to let it be determined by observation data. We formulate our approach in a previously well-established framework that incorporates the statistical modeling of the continuum variations as a damped random walk process and takes into account the long-term secular variations which are irrelevant to RM signals. The Application to RM data shows the fidelity of our approach.

preprint2016arXiv

A note on periodicity of long-term variations of optical continuum in active galactic nuclei

Graham et al. found a sample of active galactic nuclei (AGNs) and quasars from the Catalina Real-time Transient Survey (CRTS) that have long-term periodic variations in optical continuum, the nature of the periodicity remains uncertain. We investigate the periodic variability characteristics of the sample by testing the relations of the observed variability periods with AGN optical luminosity, black hole mass and accretion rates, and find no significant correlations. We also test the observed periods in several different aspects related to accretion disks surrounding single black holes, such as the Keplerian rotational periods of 5100~Å photon-emission regions and self-gravity dominated regions and the precessing period of warped disks. These tests shed new lights on understanding AGN variability in general. Under the assumption that the periodic behavior is associated with SMBHB systems in particular, we compare the separations ($\mathscr{D}_{\bullet}$) against characteristic radii of broad-line regions ($R_{\rm BLR}$) of the binaries and find $\mathscr{D}_{\bullet}$$\approx$ 0.05 $R_{\rm BLR}$. This interestingly implies that these binaries have only circumbinary BLRs.

preprint2016arXiv

Improving the Flux Calibration in Reverberation Mapping by Spectral Fitting: Application to the Seyfert Galaxy MCG--6-30-15

We propose a method for the flux calibration of reverberation mapping spectra based on accurate measurement of [O III] $λ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$β$ and $V$-band light curves. Combining this lag with FWHM(H$β$) = $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.

preprint2016arXiv

Reverberation Mapping of the Broad-line Region in NGC 5548: Evidence for Radiation Pressure?

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$β$ emission line with respect to the 5100 Åcontinuum and obtain $τ_{\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$β$ 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$β$ campaigns. The high-quality data allow us to construct a velocity-binned delay map for the broad H$β$ 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$β$ time lags and the associated mean 5100 Å luminosities over 18 campaigns between 1989 and 2015, we find that the H$β$ 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.

preprint2016arXiv

Spectroscopic Indication for a Centi-parsec Supermassive Black Hole Binary in the Galactic Center of NGC 5548

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.

preprint2016arXiv

Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. V. A New Size-Luminosity Scaling Relation for the Broad-Line Region

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$β$ lags ($τ_{_{\rm Hβ}}$) 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β}}-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$β$ lags tightly depends on accretion rates, $τ_{_{\rm Hβ}}/τ_{_{R-L}}\propto\dot{\mathscr{M}}^{-0.42}$, where $τ_{_{R-L}}$ is the H$β$ lag from the normal $R_{_{\rm Hβ}}-L_{5100}$ relation. Fitting 63 mapped AGNs, we present a new scaling relation for the broad-line region: $R_{_{\rm Hβ}}=α_1\ell_{44}^{β_1}\,\min\left[1,\left(\dot{\mathscr{M}}/\dot{\mathscr{M}}_c\right)^{-γ_1}\right]$, where $\ell_{44}=L_{5100}/10^{44}\,\rm erg~s^{-1}$ is 5100 Å continuum luminosity, and coefficients of $α_1=(29.6_{-2.8}^{+2.7})$ lt-d, $β_1=0.56_{-0.03}^{+0.03}$, $γ_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.

preprint2016arXiv

Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. VI. Velocity-resolved Reverberation Mapping of H$β$ Line

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$β$ 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$β$ 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.

preprint2016arXiv

The Fundamental Plane of the Broad-line Region in Active Galactic Nuclei

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β}}={\rm FWHM}/σ_{_{\rm Hβ}}$, the ratio of full width at half maximum to the dispersion of broad H$β$, 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$β$ line flux emission. Assembling all AGNs with reverberation mapping measurements of broad H$β$, 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})=α+β{\cal D}_{_{\rm Hβ}}+γ{\cal R}_{\rm Fe},$ where $α=(2.47,0.31)$, $β=-(1.59,0.82)$ and $γ=(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.

preprint2015arXiv

Alignments Of Black Holes With Their Warped Accretion Disks And Episodic Lifetimes Of Active Galactic Nuclei

Warped accretion disks have attracted intensive attention because of their critical role on shaping the spin of supermassive massive black holes (SMBHs) through the Bardeen-Petterson effect, a general relativistic effect that leads to final alignments or anti-alignments between black holes and warped accretion disks. We study such alignment processes by explicitly taking into account the finite sizes of accretion disks and the episodic lifetimes of AGNs that delineate the duration of gas fueling onto accretion disks. We employ an approximate global model to simulate the evolution of accretion disks, allowing to determine the gravitomagnetic torque that drives the alignments in a quite simple way. We then track down the evolutionary paths for mass and spin of black holes both in a single activity episode and over a series of episodes. Given with randomly and isotropically oriented gas fueling over episodes, we calculate the spin evolution with different episodic lifetimes and find that it is quite sensitive to the lifetimes. We therefore propose that spin distribution of SMBHs can place constraints on the episodic lifetimes of AGNs and vice versa. Applications of our results on the observed spin distributions of SMBHs and the observed episodic lifetimes of AGNs are discussed, although both the measurements at present are yet ambiguous to draw a firm conclusion. Our prescription can be easily incorporated into semi-analytic models for black hole growth and spin evolution.

preprint2015arXiv

Coexistence and interconvertibility of ferromagnetic and antiferromagnetic phases in the single crystal of Mn$_3$ single-molecule magnet

We report the coexistence of ferromagnetic and antiferromagnetic phases in the single crystal of Mn$_3$ single-molecule magnet. The coexistent state appears within a certain period of time in the progress of either oxidation or reduction during a reversible oxidation-reduction process, when the sample is exposed in the air (oxygen) or the methyl vapor. We noticed an apparent change in the molecular structure, which is also reversible in terms of that the methyl group is dropped or added to the molecules during the oxidation or reduction. The absence or presence of the methyl group in the molecules exert an essential impact upon the intermolecular exchange interaction, and the ferromagnetic phase comes from the heterogenous intermolecular bonds between pairs of molecules of which one molecule has a methyl group whereas the other has lost the methyl group. The reversible change in molecular structure suggests the magnetic structure of Mn$_3$ might be designed and modulated at molecular scale, which implies Mn$_3$ has a great application potential.

preprint2015arXiv

Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. III. Detection of Fe II Reverberation in Nine Narrow-Line Seyfert 1 Galaxies

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.

preprint2015arXiv

Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. IV. H$β$ Time Lags and Implications for Super-Eddington Accretion

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$β$ 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$β$ 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β}}$) and optical luminosity at 5100Å, $R_{_{\rm Hβ}}-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β}}$ 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.

preprint2014arXiv

A Bayesian Approach to Estimate the Size and Structure of the Broad-line Region in Active Galactic Nuclei Using Reverberation Mapping Data

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.

preprint2014arXiv

A Bayesian Method for the Intercalibration of Spectra in Reverberation Mapping

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.

preprint2014arXiv

A New Approach to Constrain Black Hole Spins in Active Galaxies Using Optical Reverberation Mapping

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.

preprint2014arXiv

Equal-Interval Splitting of Quantum Tunneling in Single-Molecule Magnets with Identical Exchange Coupling

The equal-interval splitting of quantum tunneling observed in simple-Ising-model systems of Ni$_{4}$ (3D) and Mn$_3$ (2D) single-molecule magnets (SMMs) is reported. The splitting is due to the identical exchange coupling in the SMMs, and is simply determined by the difference between the two numbers of the spin-down $n_{\downarrow}$ and spin-up $n_{\uparrow}$ molecules neighboring to the tunneling molecule. The splitting may be presented as $(n_{\downarrow}-n_{\uparrow})JS/{gμ_{0}μ_{B}}$, and the number of the splittings follows $n+1$ where $n=n_{\downarrow}+n_{\uparrow}$ is the coordination number. Besides, since the quantum tunneling is heavily dependent on local spin environment, the manipulation of quantum tunneling may become feasible for this kind of system, which may shed new light on novel applications of SMMs.

preprint2014arXiv

Spin-Pair Tunneling in Mn$_3$ Single-Molecule Magnet

We report spin-pair tunneling observed in Mn$_3$ single-molecule magnet, which is a crystal with 2D network of identical exchange coupling. We observed a series of extra quantum tunnelings by the ac susceptibility measurements, and demonstrated these are concerted tunnelings of two spins taking place from the same initial state to the same final state simultaneously. The resonant field of spin-pair tunneling can be expressed as $H_z=lD/gμ_{0}μ_{B}+(n_{\downarrow}-n_{\uparrow})JS/2{gμ_{0}μ_{B}}$, and the splitting interval ($|J|S/{gμ_{0}μ_{B}}$) is half of that of the single-spin tunneling ($2|J|S/{gμ_{0}μ_{B}}$), which is analogous to the relationship between the magnetic flux quantum in superconductor ($h/2e$) and common metal ($h/e$).

preprint2014arXiv

Supermassive black holes with high accretion rates in active galactic nuclei: II. the most luminous standard candles in the Universe

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.

preprint2013arXiv

Evolution of Warped Accretion Disks in Active Galactic Nuclei. I. Roles of Feeding at the Outer Boundaries

We investigate the alignment processes of spinning black holes and their surrounding warped accretion disks in a frame of two different types of feeding at the outer boundaries. We consider (1) fixed flows in which gas is continually fed with a preferred angular momentum, and (2) free flows in which there is no gas supply and the disks diffuse freely at their outer edges. As expected, we find that for the cases of fixed flows the black hole disk systems always end up aligning on timescales of several 1e6 yr, irrespective of the initial inclinations. If the initial inclination angles are larger than pi/2, the black hole accretion transits from retrograde to prograde fashion, and the accreted mass onto the black holes during these two phases is comparable. On the other hand, for the cases of free flows, both alignments and anti-alignments can occur, depending on the initial inclinations and the ratios of the angular momentum of the disks to that of the black holes. In such cases, the disks will be consumed within timescales of 1e6 yr by black holes accreting at the Eddington limit. We propose that there is a close connection between the black hole spin and the lifetime for which the feeding persists, which determines the observable episodic lifetimes of active galactic nuclei. We conclude that careful inclusion of the disk feeding at the outer boundaries is crucial for modeling the evolution of the black hole spin.

preprint2012arXiv

Cosmological Evolution of Supermassive Black Holes. II. Evidence for Downsizing of Spin Evolution

The spin is an important but poorly constrained parameter for describing supermassive black holes (SMBHs). Using the continuity equation of SMBH number density, we explicitly obtain the mass-dependent cosmological evolution of the radiative efficiency for accretion, which serves as a proxy for SMBH spin. Our calculations make use of the SMBH mass function of active and inactive galaxies (derived in the first paper of this series), the bolometric luminosity function of active galactic nuclei (AGNs), corrected for the contribution from Compton-thick sources, and the observed Eddington ratio distribution. We find that the radiative efficiency generally increases with increasing black hole mass at high redshifts (z>~1), roughly as η\propto M_bh^0.5, while the trend reverses at lower redshifts, such that the highest efficiencies are attained by the lowest mass black holes. Black holes with M_bh>~10^8.5M_sun maintain radiative efficiencies as high as η~0.3-0.4 at high redshifts, near the maximum for rapidly spinning systems, but their efficiencies drop dramatically (by an order of magnitude) by z~0. The pattern for lower mass holes is somewhat more complicated but qualitatively similar. Assuming that the standard accretion disk model applies, we suggest that the accretion history of SMBHs and their accompanying spins evolve in two distinct regimes: an early phase of prolonged accretion, plausibly driven by major mergers, during which the black hole spins up, then switching to a period of random, episodic accretion, governed by minor mergers and internal secular processes, during which the hole spins down. The transition epoch depends on mass, mirroring other evidence for "cosmic downsizing" in the AGN population; it occurs at z~2 for high-mass black holes, and somewhat later, at z~1, for lower-mass systems.

preprint2011arXiv

A general relativistic external Compton-scattering model for TeV emission from M87

M87 is the first detected non-blazar extragalactic Tera-Electron-Volt (TeV) source with rapid variation and very flat spectrum in the TeV band. To explain the two-peaks in the spectral energy distribution (SED) of the nucleus of M87 which is similar to those of blazars, the most commonly adopted models are the synchrotron self-Compton scattering (SSC) models and the external inverse Compton scattering (EIC) models. Considering that there is no correlated variation in the soft band (from radio to X-ray) matching the TeV variation, and the TeV sources should not suffer from the gamma-gamma absorption due to the flat TeV spectrum, the EIC models are advantageous in modeling the TeV emission from M87. In this paper, we propose a self-consistent EIC model to explain the flat TeV spectrum of M87 within the framework of fully general relativity, where the background soft photons are from the advection-dominated accretion flow (ADAF) around the central black hole, and the high energy electrons are from the mini-jets which are powered by the magnetic reconnection in the main jet (Giannios et al. 2010). In our model, both the TeV flares observed in the years of 2005 and 2008 could be well explained: the gamma-gamma absorption for TeV photons is very low, even inside the region very close to the black hole 20Rg~50Rg; at the same region, the average EIC cooling time (~ 10^2-10^3s) is short, which is consistent with the observed time scale of TeV variation. Furthermore, we also discuss the possibility that the accompanying X-ray flare in 2008 is due to the direct synchrotron radiation of the mini-jets.

preprint2011arXiv

Cosmological evolution of supermassive black holes. I. mass function at 0<z<2

We present the mass function of supermassive black holes (SMBHs) over the redshift range z=0-2, using the latest deep luminosity and mass functions of field galaxies to constrain the masses of their spheroids, which we relate to SMBH mass through the empirical correlation between SMBH and spheroid mass (the M_bh-M_sph relation). In ddition to luminosity fading of the stellar content of the spheroids, we carefully consider the variation of the bulge-to-total luminosity ratio of the galaxy populations and the M_bh/M_sph ratio, which, according to numerous recent studies, evolves rapidly with redshift. The SMBH mass functions derived from the galaxy luminosity and mass functions show very good agreement, both in shape and in normalization. The resultant SMBH mass function and integrated mass density for the local epoch (z~0) match well those derived independently by other studies. Consistent with other evidence for cosmic downsizing, the upper end of the mass function remains roughly constant since $z\approx2$, while the space density of lower mass black holes undergoes strong evolution. We carefully assess the impact of various sources of uncertainties on our calculations. A companion paper uses the mass function derived in this work to determine the radiative efficiency of black hole accretion and constraints that can be imposed on the cosmological evolution of black hole spin.

preprint2010arXiv

Episodic activities of supermassive black holes at redshift $z\le 2$: driven by mergers?

It has been suggested for quite a long time that galaxy mergers trigger activities of supermassive black holes (SMBHs) on the grounds of imaging observations of individual galaxies. To quantitatively examine this hypothesis, we calculate quasar luminosity functions (LFs) by manipulating the observed galaxy LFs ($z\lesssim 2$) and theoretical merger rates from semi-analytical formulations. We find that the model reproduces the observed quasar LFs provided that the mass ratio ($q$) of the secondary galaxy to the newly formed one changes with cosmic time. The results show that the fraction of major mergers decreases from $f_{\rm maj}\sim 0.2$ at $z\sim 2$ to $f_{\rm maj}\to 0$ at $z\sim 0$. As a consequence, the newly formed SMBHs from major mergers at $z\sim 2$ may acquire a maximal spin due to the orbital angular momentum of the merging holes. Subsequently, random accretion led by minor mergers rapidly drives the SMBHs to spin down. Such an evolutionary trend of the SMBH spins is consistent with that radiative efficiency of accreting SMBHs strongly declines with cosmic time, reported by Wang et al. (2009). This suggests that minor mergers are important in triggering activities of SMBHs at low redshift while major mergers may dominate at high redshift.

preprint2009arXiv

Constraining spins of supermassive black holes from TeV variability. II. fully general relativistic calculations

The fast variability of energetic TeV photons from the center of M87 has been detected, offering a new clue to estimate spins of supermassive black holes (SMBHs). We extend the study of Wang et al. (2008) by including all of general relativistic effects. We numerically solve the full set of relativistic hydrodynamical equations of the radiatively inefficient accretion flows (RIAFs) and then obtain the radiation fields around the black hole. The optical depth of the radiation fields to TeV photons due to pair productions are calculated in the Kerr metric. We find that the optical depth strongly depends on: (1) accretion rates as $\tautev\propto \dot{M}^{2.5-5.0}$; (2) black hole spins; and (3) location of the TeV source. Jointly considering the optical depth and the spectral energy distribution radiated from the RIAFs, the strong degeneration of the spin with the other free parameters in the RIAF model can be largely relaxed. We apply the present model to M87, wherein the RIAFs are expected to be at work, and find that the minimum specific angular momentum of the hole is $a\sim0.8$. The present methodology is applicable to M87-like sources with future detection of TeV emissions to constrain the spins of SMBHs.