Yan-Rong Li 李彦荣Astrophysics - Scientific Programming
About
This is Yan-Rong Li (李彦荣)'s homepage.
I was born in Gansu Province, China and got my Ph.D in
astrophysics at the Institute of High Energy Physics (IHEP), Chinese Academy of Sciences (CAS).
I am currently a researcher staff in astrophysics at IHEP, CAS.
My research focuses on super-massive black holes in active galacit nuclei.
I developed a number of software packages: dynamical modeling of broad-line regions (BRAINS),
reverberation mapping analysis (MICA2), time series reconstruction (RECON),
light-curve inter-calibration (PyCALI), diffusive nested sampling (CDNest),
a Pixon-based reverberation mapping analysis method (PIXON), and a package
for astronomical analysis tools.
Measuring Geometric Distances of Quasars with SARM
The geometric distances of active galactic nuclei (AGNs) are challenging to measure because of their exceptionally compact structure yet vast cosmic distances. A combination of spectroastrometry and reverberation mapping (SARM) of broad-line regions (BLRs) constitutes a novel means to probe the geometric distance of AGNs, which has recently become practically feasible owing to successful interferometric observations with VLTI/GRAVITY. Here, we perform SARM analysis of four nearby quasars, Mrk 509, PDS 456, 3C 273, and NGC 3783. Results for the former two are reported for the first time and the latter two are revisited using our improved BLR dynamical modeling that includes the radial-dependent responsivity of BLRs. This allows us to self-consistently account for the emissivity weighting of the BLR in spectroastrometry and responsivity weighting in reverberation mapping. We obtain angular-diameter distances of the four quasars, from which we derive a Hubble constant of H0=69-10+12 km/s/Mpc. Although this consititutes a large uncertainty for a measurement of H0 , it is anticipated that the precision will improve to a competitive level once a greater number of AGNs are accessible following the upgrade of GRAVITY in the near future. From SARM analysis, the black hole masses of the four quasars are also measured with the statistical uncertainty ranging from 0.06 to 0.23 dex, consistent with the correlations between black hole masses and properties of the host bulges.
Spectroastrometry and Reverberation Mapping of Active Galactic Nuclei
Spectroastrometry measures source astrometry as a function of wavelength/velocity. Reverberations of spectroastrometric signals naturally arise in broad-line regions (BLRs) of active galactic nuclei (AGNs) as a result of the continuum variations that drive responses of the broad emission lines with time delays. Such signals provide a new diagnostic for mapping BLR kinematics and geometry, complementary to the traditional intensity reverberation mapping (RM) technique. We present a generic mathematical formalism for spectroastrometric RM and show that under realistic parameters of a phenomenological BLR model, the spectroastrometric reverberation signals vary on a level of several to tens of microarcseconds, depending on the BLR size, continuum variability, and angular-size distance. We also derive the analytical expressions of spectroastrometric RM for an inclined ring-like BLR. We develop a Bayesian framework with a sophisticated Monte Carlo sampling technique to analyze spectroastrometric data and infer the BLR properties, including the central black hole mass and angular-size distance. We demonstrate the potential of spectroastrometric RM in spatially resolving BLR kinematics and geometry through a suite of simulation tests. The application to realistic observation data of 3C 273 obtains tentative, but enlightening results, reinforcing the practical feasibility of conducting spectroastrometric RM experiments on bright AGNs with the operating Very Large Telescope Interferometer as well as possibly with the planned next-generation 30 m class telescopes.
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.
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 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 and Press and Zu et al., 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 continuum variations as a damped random walk process and takes into account long-term secular variations which are irrelevant to RM signals. The application to RM data shows the fidelity of our approach.
Cosmological Evolution of Super-massive Black Holes
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, 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), while the trend reverses at lower redshifts, such that the highest efficiencies are attained by the lowest mass black holes. Assuming that the standard accretion disk model applies, we suggest that the accretion history of SMBHs and their accompanying spins evolves 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.
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 1 billion yr ago, exhibits long-term variability with a period of ~14 yr in the optical continuum and broad Hβ emission line. Remarkably, the double-peaked profile of Hβ shows systematic velocity changes with a similar period. These pieces of observations plausibly indicate that an 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 semimajor 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.
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.
Dynamical Modeling of Broad-line Regions in Active Galactic Nuclei
We employ the recently developed dynamical modeling approach for broad-line regions (BLRs) based on the method of Pancoast et al. to analyze the RM data set 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 from BLR clouds, nonlinear 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 the mass of the central BH. We apply three BLR models with different prescriptions of BLR cloud 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 $\mathrm{log}({M}_{\bullet }/{M}_{\odot })={6.23}_{-0.45}^{+0.26}$, resulting in a virial factor of $\mathrm{log}f=-{0.36}_{-0.54}^{+0.33}$ for the full width at half maximum of the Hβ line measured from the mean spectrum. The virial factors for the other measures of the Hβ line width are also presented.