DRAFT VERSION MAY 15, 218 Preprint typeset using L A TEX style emulateapj v. 5/2/11 MEASUREMENT OF THE BROAD LINE REGION SIZE IN A LUMINOUS MACHO QUASAR DORON CHELOUCHE 1, ELIRAN DANIEL 2, & SHAI KASPI 2,3 Draft version May 15, 218 arxiv:124.1862v1 [astro-ph.co] 9 Apr 212 ABSTRACT We measure the broad emission line region (BLR) size of a luminous, L 1 47 ergs 1, high-z quasar using broadband photometric reverberation mapping. To this end, we analyze 7.5 years of photometric data for MACHO 13.685.324 (z 1.72) in the B and R MACHO bands and find a time delay of 18±4 days in the rest frame of the object. Given the spectral-variability properties of high-z quasars, we associate this lag with the rest-uv iron emission blends. Our findings are consistent with a simple extrapolation of the BLR size-luminosity relation in local active galactic nuclei to the more luminous, high-z quasar population. Longterm spectroscopic monitoring of MACHO 13.685.324 may be able to directly measure the line-to-continuum time-delay and test our findings. Subject headings: galaxies: active methods: data analysis quasars: emission lines quasars: individual (MACHO 13.685.324) techniques: photometric 1. INTRODUCTION Our understanding of galaxy and black hole (BH) formation and co-evolution has progressed significantly in recent years (Wandel 1999; Ferrarese & Merritt 2; Häring & Rix 24; Peng et al. 26; Bennert et al. 21; Cisternas et al. 211) with new observations being able to probe the first epoch of quasar activity, and place interesting constraints on the mechanisms responsible for BH growth and galaxy formation (Hughes & Blandford 23; Trakhtenbrot et al. 211). Gaining further physical insight requires that the BH masses be determined with good accuracy in a large sample of objects. Presently, the best means for weighing BHs in quasars is via the reverberation mapping (RM) technique (Peterson 1993), which measures the size of the broad line region (BLR). This, combined with a measure of the velocity dispersion of the BLR, can be used to estimate the BH mass (Peterson et al. 24). To date, BH masses in 45 low-z objects have been measured in this way leading to various scaling laws with other quasar properties. These relations are often extrapolated to high-z objects to allow for the indirect estimate of their BH mass (Netzer 23). Nevertheless, it is not clear that such extrapolations are meaningful, and a more direct measure of the BH mass in such objects is highly desirable. Here we report new measurements for the BLR size in a z 1.72 luminous quasar, using photometric RM (Chelouche & Daniel 212), and by analyzing 7.5 years of data from the MACHO survey (Geha et al. 23). Section 2 presents the analysis and results, with a follow-up discussion in 3. 2. ANALYSIS & RESULTS To date, some 2 quasars have been confirmed behind the Magellanic clouds (Dobrzycki et al. 22; Geha et al. 23; Kozłowski et al. 212). In the course of analyzing their photometric data, we report our findings for MACHO 13.685.324, 1 Department of Physics, Faculty of Natural Sciences, University of Haifa, Haifa 3195, Israel; doron@sci.haifa.ac.il 2 School of Physics & Astronomy and the Wise Observatory, Tel- Aviv University, Tel-Aviv 69978, Israel; shai@wise.tau.ac.il, elirandviv@gmail.com 3 Department of Physics, Technion, Haifa 32, Israel an m R = 18.66 quasar at z 1.72 having a monochromatic luminosity 4 λ L λ (135Å) 3 1 46 ergs 1. We focus on this object as it has the best-sampled B and R light curves, as well as the highest fractional variability (Kaspi et al. 27, and references therein), F var.16, of a high-z quasar in the Geha et al. (23) sample. This object is consistent with being radio-quiet down to a flux limit of 6 mjy (5 mjy) at 21 cm (13 cm) (Marx et al. 1997; Mauch et al. 23). A more complete analysis of the MACHO quasar sample is deferred to a forthcoming publication. We choose to work with the original MACHO bandpasses rather than transform into standard B and R bands. This reduces the spectral overlap as the photometric transformations require linear combinations of the MACHO magnitudes. Light curves were filtered against poor data: all points with errors larger than ζ times the mean error or those which deviate from the light curve mean by more than ζ standard deviations were discarded (Fig. 1). 2.1. Spectral decomposition To effectively use the photometric RM technique of Chelouche & Daniel (212), it is helpful to identify the band with the larger contribution of emission lines to its variance. If the emission line contribution to the flux may be used as a proxy to its variable component (e.g., as would be the case if the relative flux variations in all emission lines were similar; see 3), then by spectral decomposition and knowledge of the instrumental throughput, one can identify the line-rich and line-poor bands. Unfortunately, the quality of available spectra of MACHO 13.685.324 is low, and the data are not 4 To estimate the monochromatic luminosity of MACHO 13.685.324 at 135Å, we take the reported R-magnitudes from Geha et al. (23), employ a standard K-correction, and account for the extinction behind the LMC using the results of Zaritsky et al. (24, see: http://ngala.as.arizona.edu/dennis/lmcext.html). To this end, we average over 5 stars of all types in the direction of the quasar, and identify the upper turn-over in the extinction distribution, at A V 1.3 mag, with the effective extinction of the background source. We conservatively estimate the uncertainty on the luminosity to be a factor 3, to account for the potentially patchy structure of LMC disk. Our estimate is consistent with m B -m R for MA- CHO 13.685.324 being larger by.67 mag than typical, as concluded by integrating the Vanden Berk et al. (21) composite quasar spectrum over the MACHO bands, and assuming a Galactic extinction law. Using concordance (.3,.7,.3) cosmology, we find λl λ (135Å) 3 1 46 ergs 1.
2 Chelouche et al. fr (mag) fb (mag) 18 18.5 19 19.5 18.5 19 19.5 2 914 pts 998 pts 91 pts 12 pts 5 1 15 2 25 3 Days Figure 1. R (upper panel) and B (lower panel) light curves for MA- CHO 13.685.324 cover 7.5 years in the observed frame ( 2.8 years in the quasar frame). The number of points in each band is of order 1 3, and depends on the filtering scheme used: darker shades correspond to effectively unfiltered data (ζ = 1) while brighter colors use a ζ = 2 filter. The fractional variability measure, F var.14 (.16) when ζ = 2 (ζ = 1) filter is used. The mean photometric error is 5.5% (7.6)% for ξ = 2 (ξ = 1). flux-calibrated. Nevertheless, their broad consistency with the composite quasar spectrum of Vanden Berk et al. (21, see Fig. 2) motivates us to use the latter for our purpose. Prominent emission lines at z = 1.72 include the C III] λ 199, Mg II λ 2799, as well as the iron complexes (Fig. 2). To account for the iron-blend contribution to the spectrum, we use the iron template of Vestergaard & Wilkes (21). An underlying powerlaw continuum model was assumed, F λ λ 1.63, which is somewhat different from the one used by Vanden Berk et al. (21), and better accounts for the flux level at 6Å (observed), where little iron blend and Balmer continuum emission is present. Emission lines and blends were individually convolved with a single Gaussian kernel of some amplitude and width, and a qualitative agreement was sought with the Vanden Berk et al. (21) composite spectrum (Fig. 2). The solution is not strictly unique, nor very physical, and its sole purpose is to qualitatively assess the relative contribution of the various emission components to the broadband flux. Weighing the decomposed spectrum by the instrumental throughput (Alcock et al. 1999), we find that the R-band has the greater relative contribution of emission lines to its flux. Specifically, the iron emission blends contribute a net of 1% to the flux in the bands, with the relative contribution of the Mg II λ 2799 and C III] λ 199 to their respective bands being 4 times smaller. 2.2. Uncovering the time-lag Having identified the line-rich band, we now proceed to measure the line-to-continuum time-delay. Following Chelouche & Daniel (212), we compute two statistical estimators for the line-to-continuum cross-correlation function: ξ CA (δt) = f R f B f B f B, and ξ AA (δt)= f R f R f B f B ( denotes convolution and f B, f R are the light curves in the B and R bands, respectively; see Fig. 3). With these definitions, a peak at δt > indicates an emission component in the R-band, which lags behind the B-band. To calculate these statistical estimators, we use two independent schemes: the interpolated cross-correlation function (Peterson et al. 24, ICCF) and the z-transformed discrete correlation function (Alexander 1997, ZDCF). The former has the advantage of being somewhat more sensitive, while the latter is less af- Rest wavelength (Å) 16 18 2 22 24 26 28 1.9.8.7.6.5.4.3.2.1 HeII OIII] NIII] AlIII SiIII] CIII] MACHO B-band Fλ λ 1.63 instrumental MACHO R-band 45 5 55 6 65 7 75 8 Observed wavelength (Å) Figure 2. A quasar composite spectrum at z = 1.72 (Vanden Berk et al. 21) overlaid with the throughput curves of the MACHO B- and R-filters [the somewhat erratic behavior is real and results from an interference pattern (C. Stubbs; private communication)]. A qualitative spectral decomposition into the prominent emission components is shown including the iron blend (gray curve), the C III] λ199, and Mg II λ2799 lines (magenta curves), and a powerlaw continuum (dashed line). Designated minor emission lines were also included in the spectral decomposition (green curve). Integrating over the emission line flux, and taking into account the wavelengthdependent throughput, the largest contribution to line emission is from the iron complexes in the R band. The contribution of the C III and the Mg II lines is roughly 4 times smaller than iron s and their relative contribution to their respective bands is similar. Also shown are the two (non fluxcalibrated, slightly shifted in wavelength range for clarity) optical spectra for MACHO 13.685.324, which were multiplied by a power-law function of wavelength so that their trend roughly follows that of the composite spectrum (dashed blue-shaded curves). The spectral shape is qualitatively consistent with the quasar template (note the instrumental feature at around 73Å) and suggests a significant contribution of the iron emission blend to the flux in the R-band. Note the prominent atmospheric absorption feature which happens to coincide with the Mg II emission line at this particular redshift. fected by sampling. Results for ξ CA and ξ AA using the ICCF and ZDCF schemes, and for two levels of light curve filtering, are presented in figure 3. All analyses give consistent results: ξ CA and ξ AA peak at δt 4 6 days. Unsurprisingly, the peak is more significant when modest filtering (ζ = 2) is applied, resulting in ξ CA, ξ AA being significant at the 2.8σ, 5σ (σ is the standard deviation) level, respectively. To verify that the observed signal is indeed associated with a lagging component in the R-band, we reversed our choice of line-rich and line-poor bands, and repeated the analysis: no significant peak was detected at δt > (instead, a highly significant trough was detected at δt 4 6 days; not shown). In addition, we carried out sets of Monte Carlo simulations (Chelouche & Daniel 212) wherein the B-band light curve was convolved with a Gaussian kernel having prescribed time-delay and width (the latter was taken to be half of the former but the results are not very sensitive to the particular kernel chosen). The resulting light curve was then scaled down and added to the B-band signal so that the line contribution to the flux is 1%. The combined light curve was then sampled with the cadence of the R-band to create a synthetic light curve, fr s. The ξ -estimators were calculated for the light curve pair [ f B, fr s ], and the results for several input time-lags are shown in figure 3. Clearly, the observed signal is qualitatively reproduced by our simulations for lags of the order of 4 days. This a) confirms that the observed signal is indeed associated with a lagging emission component in the R-band and is not spurious, and b) that the time-lag constraints deduced below are meaningful and are not significantly biased. We estimate the time-lag and its uncertainty in the fol- MgII
Correlation coefficient The BLR size of a Luminous MACHO Quasar 3.8.6.4.2.2 2 4 6 8 1 f R f B f B f B fr s(3days) f B fr s(4days) f B fr s(5days) f B ζ = 2 ζ = 1 3 days 5 days 4 days 2 4 6 8 2 4 6 8 Figure 3. Photometric RM analysis. Top panel: an excess power is seen at around 5 days in the cross-correlation signal compared to the auto-correlation signal (thick lines). The simulated signal qualitatively reproduces the excess power in the cross-correlation term, at the right timescales, for three input time lags (see text and legend). Bottom panels: the calculated ξ-functions for two levels of filtering: ζ = 2 (left) and ζ = 1 (right). Solid curves show the results of the ICCF algorithm with shaded regions indicating the uncertainty. The ZDCF calculations are shown as points with error bars in dotted lines. Blue (red) colors correspond to ξ CA (ξ AA ). All methods give consistent results and a modest level of light curve filtering yields a more significant result. We constrain the time scale at which the ξ-functions peak to be 5 ± 1 days in the observed frame. Repeating the analysis using model-dependent synthetic R-band light curves (blue-shaded curves; see text) shows an agreement with a time lag of 4 days; considerably shorter or longer lags are excluded by the data. lowing way: the peak in ξ CA and ξ AA is identified with the time-lag for each Monte Carlo realization (a total of 1 realizations are used for each numerical scheme). A timelag distribution is then obtained whose mean is identified with the time lag and its uncertainty with the standard deviation. We find that the time-lag is 48 ± 1/53 ± 7 days (ZDCF), or 49±8/49±5days (ICCF) for the ξ CA /ξ AA statistical estimator. We conservatively estimate the lag to be 5±1 days, which corresponds to a rest-frame lag, τ 18±4 days. Geha Kozlowski C IV HeII AlIII SiIII] C III] 3. DISCUSSION & CONCLUSIONS With only photometric data at hand, interpreting our findings is not straightforward: while we have a statistically robust time lag measurement, it is difficult to associate it with a specific emission line 5. That being said, it is very likely that the observed time-lag is associated with the iron emission blends for the following reasons: 1) a positive peak in ξ CA, ξ AA indicates that emission lines contributing to the R- band are responsible for the signal, 2) the iron contribution to the R-band typically exceeds that of the Mg II 2799Å line by a factor 4(Fig. 2), and 3) the implied flux variation of the iron blend is consistent with that seen in low-z objects (Maoz et al. 1993; Vestergaard & Peterson 25), and indirectly deduced for high-z quasars (Meusinger et al. 211) 6. 5 Another possibility is that the signal originates in an unidentified continuum emission component. While this cannot be excluded based on available data, we find it rather unlikely given that continuum emission, by definition, extends over a broad wavelength range, and the fact that there is no clear reason why the lag should be similar to that which characterizes emission lines. 6 In a recent work, Meusinger et al. (211) analyzed Stripe 82 SDSS quasars and found that the combined flux from all emission lines and blends varies by 1%, which is consistent with the relative variations measured FWHM=42-92 kms 1 FWHM=59-92kms 1 Fe emission 4 42 44 5 52 54 Observed wavelength (Å) Figure 4. The arbitrarily normalized (non flux-calibrated) spectra of MA- CHO 13.685.324 (Geha et al. 23, and S. Kozłowski, private communication; see legend) show two prominent emission lines. A cubic spline was fit far from the prominent line locations to normalize the spectrum (not shown), and single Gaussian models were overlaid to constrain the FWHM range, as indicated by the shaded areas. Only the red wing of the C IV λ1549 line was used for fitting purposes due to potential absorption features just blue-ward of its peak (left panel). Similarly, only the red wing of the C III] λ199 line was used for FWHM estimation to avoid blending with other emission lines just blueward of its center (right panel). We note that FWHM measurements for this line should be treated with caution since iron blend emission may be substantial (see also Fig. 2), and is dififcult to estimate given the poor S/N of the non flux-calibrated data. Little is known with confidence about the physical properties of the UV iron-emitting region, and a statistically significant time-lag, reflecting on its size, exists for only one additional low-luminosity source, NGC5548 (Maoz et al. 1993). for individual lines by Kaspi et al. (27). This implies that similar flux variations of the iron blend are characteristics of quasars [see also Wilhite et al. (25) for some additional constraints on the time variability of emission lines in high-z quasars].
4 Chelouche et al. BLR size (light days) 1 3 1 2 1 1 1 1 1 NGC 4395 MACHO 13.685.324 S5 836+71 CIV λ1549 rest-uv iron blend 3C 39.3 MgII λ28 CIII] λ199 NGC 5548 NGC 3783 NGC 4151 NGC 7469 Kaspi et al. (27): FITEXY Peterson et al. (25) Kaspi et al. (27): BCES Mass (M ) 1 1 1 9 1 8 1 7 1 6 1 3 LEdd.1LEdd.1LEdd LEdd MACHO 13.685.324 1 2 1 4 1 41 1 42 1 43 1 44 1 45 1 46 1 47 λl λ(135å) (erg s 1 ) 1 5 1 4 1 41 1 42 1 43 1 44 1 45 1 46 1 47 λl λ(135å) (erg s 1 ) Figure 5. The BLR size-luminosity relation (left; adopted from Kaspi et al. (27); see their paper for a discussion of the various fitted relations) and the massluminosity relation [right; adopted from Peterson et al. (24, and augmented by up-tp-date data from Peterson et al. (25); Denney et al. (26); Kaspi et al. (27); Grier et al. (28); Bentz et al. (29); Denney et al. (21); Grier et al. (212)) using the composite spectrum of Vanden Berk et al. (21) to convert monochromatic luminosities to λl λ (135Å); uncertainties are not shown]. Left: plotted are data from Peterson et al. (24, 25); Kaspi et al. (25, 27) as well as our results for MACHO 13.685.324, which are in agreement with both relations. Also shown, for comparison, are reliable BLR size measurements for the iron blend (Maoz et al. 1993) and the C III λ199 and Mg II λ2799 emission lines (Metzroth et al. 26); see legend. Evidently, the BLR size measured here lies on the extrapolation of the R L relation, as obtained for low luminosity, low-z objects. Right: Interestingly, MACHO 13.685.324 is consistent with emitting at 1% of its Eddington rate, in agreement with the efficiency of other bright quasars [note that the mass reported here for S5 836+71 is 1.8 times larger than the value reported in Kaspi et al. (27) and follows from the normalization of Onken et al. (24)]. Curves of constant L/L Edd. are also plotted (dashed lines) assuming the bolometric correction used by Kaspi et al. (27). Taken together, the results imply that the size-luminosity scaling for this region is consistent with that of other broad lines (Kaspi et al. 25, 27, see our Fig. 5). Interestingly, the inferred sizes of the iron- and C IV λ 1549-emitting regions are comparable in NGC 5548 (Maoz et al. 1993). If true also for MACHO 13.685.324 then we find consistency with a simple extrapolation of size-luminosity relation for the C IV-emitting region, as obtained for lower luminosity sources, to the more luminous quasar population (Fig. 5). A meaningful comparison with the size-luminosity relations for other rest-uv lines is currently limited by small number statistics. Once the size of the line-emitting region is known, the black hole mass may be determined via M BH 2.7 1 5 (FWHM/1 3 kms 1 ) 2 (τ/days)m (Peterson et al. 24, and using the normalization of Onken et al. (24)), where FWHM is the full width at half maximum of the relevant emission line. This expression assumes that the gas is virialized and uses the FWHM as a proxy for its velocity dispersion in the black hole s potential well. Nevertheless, the FWHM of the iron blend in MACHO 13.685.324 cannot be determined from the available spectra, and it is not clear that the iron-emitting gas is virialized in quasars (Hu et al. 28). With these uncertainties in mind, we note that currently available data point to a roughly comparable size for the iron, C IV λ 1549, and C III λ 199 emitting regions (Maoz et al. 1993; Metzroth et al. 26, see our Fig. 5) 7. If true also for MACHO 13.685.324 then we may use the FWHM of the carbon lines and our iron line region size measurement to estimate M BH 8. Taking FWHM=65±25km s 1 (Fig. 4), we estimate the black hole mass to be 2 1 9 M with a formal uncertainty of.3 dex (see Fig. 5). We note, however, that this mass estimate does not account for potentially important, but poorly understood, systematic effects concerning the physics and geometry of the iron-blend emitting gas. It is worth noting that, being>1 quieter in the radio band than S5 836+71 (Kaspi et al. 27), MACHO 13.685.324 is not a blazar, and may well be a radio-quiet object. As such, it may have the largest black hole mass ever measured using the RM technique in a representative member of the quasar population. Our findings also indicate that MACHO 13.685.324 is probably shining at 1% of its Eddington luminosity (Fig. 5), and in agreement with similar quasars in its class (Shemmer et al. 26). A more reliable interpretation of our results, as well as a more accurate BH mass determination, require a higher signal-to-noise spectrum. A sparse (one visit every two months) spectroscopic-monitoring campaign of MA- CHO 13.685.324, over a period of 7 years, will determine the variability amplitude of the iron blend as well as potentially (independently) measure the time-lag in this object. Furthermore, by analyzing large enough samples of MACHO quasars, with a similar redshift and luminosity range to MA- CHO 13.685.324, it may be possible to statistically corroborate the results presented here (Chelouche & Daniel 212). With these future tests and improvements in mind, our present findings already demonstrate the feasibility of long-term photometric surveys in determining the BLR size (and BH mass) in luminous high-z quasars. 7 The situation concerning the size of the C III λ199 region is far from being clear: two low-luminosity objects having only loose constraints on the time-lag indicate that it may be twice the size of the C IV λ1549 region (Peterson 1994), while a tentative time-lag measurement in a high-z quasar indicates that the line traces the continuum level with no apparent delay (Kaspi et al. 27). 8 There are uncertainties associated with the use of carbon lines for measuring M BH : as noted by Baskin & Laor (25), the C IV λ1549 FWHM may be affected by non-virial motions, especially for bright objects emitting close to their Eddington limit. Also, line de-blending may be required to properly estimate the FWHM of the C III λ199 line. Lastly, a generic problem We are grateful to Marla Geha for providing us with the MACHO quasar light curves and spectra in electronic form, and for commenting on an earlier version of this letter. We thank S. Kozłowski for providing us with a recently acquired spectrum of MACHO 13.685.324, and M. Vestergaard for supplying us with iron emission templates in electronic form. with single-epoch FWHM estimates concerns the fact that the FWHM of the varying component of the line may differ from that measured from the mean spectrum (Kaspi et al. 2).
The BLR size of a Luminous MACHO Quasar 5 Fruitful discussions with H. Netzer and O. Shemmer are greatly appreciated, as well as helpful comments by the referee. This research has been supported in part by a FP7/IRG PIRG-GA-29-256434 grant as well as by grant 927/11 from the Israeli Science Foundation. REFERENCES Alcock, C., et al. 1999, PASP, 111, 1539 Alexander, T. 1997, Astronomical Time Series, 218, 163 Baskin, A., & Laor, A. 25, MNRAS, 356, 129 Bennert, V. N., Treu, T., Woo, J.-H., Malkan, M. A., Le Bris, A., Auger, M. W., Gallagher, S., & Blandford, R. D. 21, ApJ, 78, 157 Bentz, M. C., Peterson, B. M., Netzer, H., Pogge, R. W., & Vestergaard, M. 29, ApJ, 697, 16 Vanden Berk, D. E., et al. 21, AJ, 122, 549 Cisternas, M., et al. 211, ApJ, 726, 57 Chelouche, D., & Daniel E. 212, ApJ, 747, 62 Denney, K. D., Bentz, M. C., Peterson, B. M., et al. 26, ApJ, 653, 152 Denney, K. D., Peterson, B. M., Pogge, R. W., et al. 21, ApJ, 721, 715 Dobrzycki, A., Groot, P. J., Macri, L. M., & Stanek, K. Z. 22, ApJ, 569, L15 Ferrarese, L., & Merritt, D. 2, ApJ, 539, L9 Geha, M., et al. 23, AJ, 125, 1 Grier, C. J., Peterson, B. M., Bentz, M. C., et al. 28, ApJ, 688, 837 Grier, C. J., Peterson, B. M., Pogge, R. W., et al. 212, ApJ, 744, L4 Häring, N., & Rix, H.-W. 24, ApJ, 64, L89 Hu, C., Wang, J.-M., Ho, L. C., et al. 28, ApJ, 687, 78 Hughes, S. A., & Blandford, R. D. 23, ApJ, 585, L11 Kaspi, S., Brandt, W. N., Maoz, D., Netzer, H., Schneider, D. P., & Shemmer, O. 27, ApJ, 659, 997 Kaspi, S., Maoz, D., Netzer, H., Peterson, B. M., Vestergaard, M., & Jannuzi, B. T. 25, ApJ, 629, 61 Kaspi, S., Smith, P. S., Netzer, H., Maoz, D., Jannuzi, B. T., & Giveon, U. 2, ApJ, 533, 631 Kozłowski, S., Kochanek, C. S., Jacyszyn, A. M., et al. 212, ApJ, 746, 27 Maoz, D., et al. 1993, ApJ, 44, 576 Mauch, T., Murphy, T., Buttery, H. J., Curran, J., Hunstead, R. W., Piestrzynski, B., Robertson, J. G., & Sadler, E. M. 23, MNRAS, 342, 1117 Marx, M., Dickey, J. M., & Mebold, U. 1997, A&AS, 126, 325 Metzroth, K. G., Onken, C. A., & Peterson, B. M. 26, ApJ, 647, 91 Meusinger, H., Hinze, A., & de Hoon, A. 211, A&A, 525, A37 Netzer, H. 23, ApJ, 583, L5 Onken, C. A., Ferrarese, L., Merritt, D., et al. 24, ApJ, 615, 645 Page, K. L., Reeves, J. N., O Brien, P. T., & Turner, M. J. L. 25, MNRAS, 364, 195 Peng, C. Y., Impey, C. D., Ho, L. C., Barton, E. J., & Rix, H.-W. 26, ApJ, 64, 114 Peterson, B. M. 1993, PASP, 15, 247 Peterson, B. M. 1994, Reverberation Mapping of the Broad-Line Region in Active Galactic Nuclei, 69, 1 Peterson, B. M., et al. 24, ApJ, 613, 682 Peterson, B. M., et al. 25, ApJ, 632, 799 Shemmer, O., Brandt, W. N., Netzer, H., Maiolino, R., & Kaspi, S. 26, ApJ, 646, L29 Trakhtenbrot, B., Netzer, H., Lira, P., & Shemmer, O. 211, ApJ, 73, 7 Vestergaard, M., & Peterson, B. M. 25, ApJ, 625, 688 Vestergaard, M., & Wilkes, B. J. 21, ApJS, 134, 1 Wandel, A. 1999, ApJ, 519, L39 Wilhite, B. C., Vanden Berk, D. E., Kron, R. G., et al. 25, ApJ, 633, 638 Zaritsky, D., Harris, J., Thompson, I. B., & Grebel, E. K. 24, AJ, 128, 166