Reverberation Mapping Data Research Articles

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.

ARE OUTFLOWS BIASING SINGLE-EPOCH C IV BLACK HOLE MASS ESTIMATES?

We use a combination of reverberation mapping data and single-epoch spectra of the CIV emission line in a sample of both low and high-redshift active galactic nuclei (AGNs) to investigate sources of the discrepancies between CIV- and Hbeta-based single-epoch black hole mass estimates. We find that for all reverberation mapped sources, there is a component of the line profile that does not reverberate, and the velocity characteristics of this component vary from object-to-object. The differing strength and properties of this non-variable component are responsible for much of the scatter in CIV-based black hole masses compared to Hbeta masses. The CIV mass bias introduced by this non-variable component is correlated with the shape of the CIV line, allowing us to make an empirical correction to the black hole mass estimates. Using this correction and accounting for other sources of scatter such as poor data quality and data inhomogeneity reduces the scatter between the CIV and Hbeta masses in our sample by a factor of ~2, to only ~0.2 dex. We discuss the possibility that this non-variable CIV component originates in an orientation-dependent outflow from either the proposed broad line region (BLR) disk-wind or the intermediate line region (ILR), a high-velocity inner extension of the narrow line region (NLR).

THE LICK AGN MONITORING PROJECT 2011: DYNAMICAL MODELING OF THE BROAD-LINE REGION IN Mrk 50

We present dynamical modeling of the broad line region (BLR) in the Seyfert 1 galaxy Mrk 50 using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. We model the reverberation mapping data directly, constraining the geometry and kinematics of the BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or virial coefficient. We find that the geometry of the BLR in Mrk 50 is a nearly face-on thick disk, with a mean radius of 9.6(+1.2,-0.9) light days, a width of the BLR of 6.9(+1.2,-1.1) light days, and a disk opening angle of 25\pm10 degrees above the plane. We also constrain the inclination angle to be 9(+7,-5) degrees, close to face-on. Finally, the black hole mass of Mrk 50 is inferred to be log10(M(BH)/Msun) = 7.57(+0.44,-0.27). By comparison to the virial black hole mass estimate from traditional reverberation mapping analysis, we find the normalizing constant (virial coefficient) to be log10(f) = 0.78(+0.44,-0.27), consistent with the commonly adopted mean value of 0.74 based on aligning the M(BH)-{\sigma}* relation for AGN and quiescent galaxies. While our dynamical model includes the possibility of a net inflow or outflow in the BLR, we cannot distinguish between these two scenarios.

GEOMETRIC AND DYNAMICAL MODELS OF REVERBERATION MAPPING DATA

We present a general method to analyze reverberation mapping data that provides both estimates for the black hole mass and for the geometry and dynamics of the broad line region (BLR) in active galactic nuclei (AGN). Our method directly infers the spatial and velocity distribution of the BLR from the data, allowing us to easily derive a velocity-resolved transfer function and allowing for a self-consistent estimate of the black hole mass without a virial coefficient. We obtain estimates and reasonable uncertainties of the BLR model parameters by implementing a Markov Chain Monte Carlo algorithm using the formalism of Bayesian probability theory. We use Gaussian Processes to interpolate the the continuum light curve data and create mock light curves that can be fitted to the data. We test our method by creating simulated reverberation mapping data-sets with known true parameter values and by trying to recover these parameter values using our models. We are able to recover the parameters with realistic uncertainties that depend upon the variability of the AGN and the quality of the reverberation mapping campaign. With a geometry model we can recover the mean radius of the BLR to within ~0.1dex random uncertainty for simulated data with an integrated line flux uncertainty of 1.5%, while with a dynamical model we can recover the black hole mass and the mean radius to within ~0.05dex random uncertainty, for simulated data with a line profile average signal to noise ratio of 4 per spectral pixel. These uncertainties do not include modeling errors, which are likely to be present in the analysis of real data, and should therefore be considered as lower limits to the accuracy of the method.

Diagnostics for the structure of AGNs’ broad line regions with reverberation mapping data: confirmation of the two-component broad line region model

We re-examine the ten Reverberation Mapping (RM) sources with public data based on the two-component model of the Broad Line Region (BLR). In fitting their broad H-beta lines, six of them only need one Gaussian component, one of them has a double-peak profile, one has an irregular profile, and only two of them need two components, i.e., a Very Broad Gaussian Component (VBGC) and an Inter-Mediate Gaussian Component (IMGC). The Gaussian components are assumed to come from two distinct regions in the two-component model; they are Very Broad Line Region (VBLR) and Inter-Mediate Line region (IMLR). The two sources with a two-component profile are Mrk 509 and NGC 4051. The time lags of the two components of both sources satisfy $t_{IMLR}/t_{VBLR}=V^2_{VBLR}/V^2_{IMLR}$, where $t_{IMLR}$ and $t_{VBLR}$ are the lags of the two components while $V_{IMLR}$ and $V_{VBLR}$ represent the mean gas velocities of the two regions, supporting the two-component model of the BLR of Active Galactic Nuclei (AGN). The fact that most of these ten sources only have the VBGC confirms the assumption that RM mainly measures the radius of the VBLR; consequently, the radius obtained from the R-L relationship mainly represent the radius of VBLR. Moreover, NGC 4051, with a lag of about 5 days in the one component model, is an outlier on the R-L relationship as shown in Kaspi et al. (2005); however this problem disappears in our two-component model with lags of about 2 and 6 days for the VBGC and IMGC, respectively.

The Effect of Radiation Pressure on Virial Black Hole Mass Estimates and the Case of Narrow‐Line Seyfert 1 Galaxies

We consider the effect of radiation pressure from ionizing photons on black hole (BH) mass estimates based on the application of the virial theorem to broad emission lines in AGN spectra. BH masses based only on the virial product Δ V2R and neglecting the effect of radiation pressure can be severely underestimated, especially in objects close to the Eddington limit. We provide an empirical calibration of the correction for radiation pressure, and we show that it is consistent with a simple physical model in which BLR clouds are optically thick to ionizing radiation and have average column densities of NH ∼ 1023 cm −2. This value is remarkably similar to what is required in standard BLR photoionization models to explain observed spectra. With the inclusion of radiation pressure, the discrepancy between virial BH masses based on single-epoch spectra and on reverberation mapping data drops from 0.4 to 0.2 dex rms. The use of single-epoch observations as surrogates of reverberation mapping campaigns can thus provide more accurate BH masses than previously thought. Finally, we show that narrow-line Seyfert 1 (NLS1) galaxies have apparently low BH masses because they are radiating close to their Eddington limit. After the radiation pressure correction, NLS1 galaxies have BH masses similar to other broad-line AGNs and follow the same MBH − σe/Lsph relations as other active and normal galaxies. Radiation forces arising from ionizing photon momentum deposition constitute an important physical effect which must be taken into account when computing virial BH masses.

Black Hole Masses in Three Seyfert Galaxies

We analyze published reverberation mapping data for three Seyfert galaxies (NGC 3227, NGC 3516, and NGC 4593) to refine the mass estimate for the supermassive black hole in the center of each object. Treatment of the data in a manner more consistent with other large compilations of such masses allows us to more securely compare our results to wider samples of data, e.g., in the investigation of the M(sub bh-sigma*) relationship for active and quiescent galaxies.