Reverberation Mapping of High-z, High-luminosity Quasars

We present Reverberation Mapping results after monitoring a sample of 17 high-z, high-luminosity quasars for more than 10 years using photometric and spectroscopic capabilities. Continuum and line emission flux variability is observed in all quasars. Using cross-correlation analysis we successfully determine lags between the variations in the continuum and broad emission lines for several sources. Here we present a highlight of our results and the determined radius--luminosity relations for Ly_alpha and CIV.


INTRODUCTION
Reverberation Mapping (RM) has been an extremely successful technique used to study the innermost regions of Active Galactic Nuclei (AGN). The determination of lags between variations in the continuum emission coming from the accretion disk near the central Black Hole (BH), and the response from the emission lines produced in the Broad Line Region (BLR) have shown that the BLR is an extended, virialized and ionizion-stratified structure. Furthermore, the determination of the radius-luminosity relation between the BLR radius at which Hβ is produced and the luminosity of the central source has open, through cross-calibration to other BLR lines, the possibility to measure BH masses in hundred of thousands of sources. The cross-calibration is, however, subject to many uncertainties due to the extrapolations necessary to apply the radius-luminosity relation to sources of very different luminosities to those actually probed with RM experiments, and to the use of emission lines produced by regions of the BLR that can be far from that producing Hβ. This is the motivation to conduct RM campaigns in high-z, high-luminosity quasars for those emission lines available in the observed optical domain. selected from the SDSS (Schneider et al., 2005) and Calán-Tololo samples (Maza et al., 1996). We started with a purely R-band imaging (corresponding to rest frame wavelengths ∼ 1700 − 1800Å, depending on the exact redshift of the source) and two years later triggered the first spectroscopic observations of those quasars with the largest photometric variations. Over the years the campaign was narrowed down to a final sample of 17 quasars which have good quality R-band and emission line light curves.
Line fluxes were measured using spectral windows tailored to each line and each quasar. We avoided regions where the lines were contaminated by variable absorbing features, but did not attempt to correct for the contribution of other (weaker) emission lines, either narrow or broad. For further details see Lira et al. (2018, submitted).
Over the campaign we found that most quasars showed a substantial degree of variability in the continuum and line emission line fluxes, with typical normalized variability amplitudes f var (Rodríguez-Pascual et al., 1997) of ∼ 10%. Most quasars also showed the expected behavior, where the emission lines followed the trends seen in the continuum as traced by the R-band light curves. Many sources, unfortunately, did not shown enough structure in their light curves (which in the rest frame only map the continuum and line flux variations during ∼ 3 years), to allow for statically significant lag determinations.
Two sources showed unexpected Lyα light curves, where the line fluxes depart from the behavior shown by the continuum and the remaining emission lines. One example (J224743) is shown in Figure 1, where we also include a source that presents the expected line response to the continuum variations (CT650).

TIME SERIES ANALYSIS
Cross-correlation analysis was conducted using the ICCF and ZCCF methods (Gaskell and Sparke, 1986;Alexander, 1997). Errors were determined using Monte Carlo simulations were the light curve fluxes were randomized and bootstrapped to contain about 70% of the original data points (Peterson et al., 2004). We determined statistically significant lags with respect to that of the R-band continuum for 3 Lyα, 5 CIV, 1 SiIV, 1 CIII], and 1 MgII emission line light curves.
In Figure 2 we reproduce the light curves for two of our sources, CT286 and J221516. The emission line light curves have been shifted according to the calculated lags, while all curves have been taken to a mean of zero and a standard deviation of one. Figure 2 illustrates that emission line light curves can closely match the observed UV continuum, like in the case of J221516, or can show rather different trends, as seen in CT286, where the line light curves correspond to a heavily smoothed version of the continuum light curve.
Two reasons can be invoked to explain such differences: 1) the response of the BLR to continuum variations differs from object to object, and it can be nonlinear and show variations with time; 2) the observed UV continuum might not be a good representation of the ionizing continuum, which is actually responsible of driving the emission line changes.
The recent monitoring of NGC5548 by the STORM consortium displays several of these traits during the 170 days of monitoring (Goad et al., 2016). Emission line light curves follow closely that of the continuum during the first 1/3 of the campaign, to then disengage from it showing a decorrelated behavior for ∼ 60 − 70 days, to finally going back to the original state. Besides, while some emission lines show a smoother light curve than that of the observed continuum during the last segment of the monitoring, SiIV stands out for showing larger amplitude in its peaks and troughs than that of the continuum.

RADIUS-LUMINOSITY RELATIONS
As well as providing extremely valuable information about the innerworks of the BLR (see next section), reverberation mapping has provided a huge scientific legacy with the determination of the so called radius-luminosity relations. This tight correlations between the distance at which one particular emission line is produced and the continuum luminosity of the central engine allows for the determination of BH masses by applying these calibrations to a virialized BLR (i.e., M BH ∝ Rv 2 , where R comes from the radius-luminosity relation and v is measured from the width of the broad emission lines).
So far, RM of the Hβ line for nearby (z < 0.3) AGN has produced a solid radius-luminosity relation for this line (Wandel et al., 1999;Kaspi et al., 2000Kaspi et al., , 2005Bentz et al., 2006Bentz et al., , 2009Bentz et al., , 2013. Cross-calibration of the correlation to other lines has been a significant enterprise which has allowed to determine BH masses of high-z quasars whose Balmer lines are redshifted into the infrared realm. In particular, MgII has proven to be a safe line to be used as BH mass estimator (McLure and Dunlop, 2004;Trakhtenbrot and Netzer, 2012), while it has been extensively shown that CIV yields unreliable results (Baskin and Laor, 2005;Netzer et al., 2007;Shen et al., 2008;Shen and Liu, 2012;Mejía-Restrepo et al., 2016).
Our monitoring effort has provided a sizeable number of Lyα and CIV lags at the high-luminosity end of the quasar distribution. This, together with other measurements found in the literature for lower luminosity AGN allows us now to determine radius-luminosity relations for these lines. These are presented in Figure  3, while the analytical expressions are as follows: where R line is the measured lag and λL λ (1345Å) corresponds to the product λ × L λ as measured at 1345 A from each spectra. As can be seen from Equation 1, the zero point and slope of the Lyα radius-luminosity correlation are poorly constrained. The reason is clear after an inspection of Figure 3 which reveals that Seyfert type sources show a large dispersion around the best fitted correlation, in contrast with the situation for CIV. In fact, the new CIV radius-luminosity relation is very close to that reported by (Kaspi et al., 2007), who found a zero point of 0.24 ± 0.06 and a slope of 0.55 ± 0.04. Notice that the Lyα and CIV relations are consistent within the errors.

BLR STRATIFICATION
One crucial result that emerged early during the reverberation campaigns of nearby AGN was that the BLR was compact, dense and stratified (Peterson, 1994), in contrast with previous photoionization results that attempted to explain all emission lines as produced by one set of physical parameters (or one single cloud). RM results made clear that different regions, with different properties, and located at different distances from the central engine, were producing the observed set of emission lines. To build a consistent picture of the BLR, therefore, it is important to determine where different lines are produced.
Our results allows to put constraints on the distance at which Lyα and CIV are produced by obtaining R Lyα /R CIV for all source for which both lags have been measured. This is presented in Figure 4, where  . Lyα (left) and CIV (right) radius-luminosity relations. Blue triangles correspond to results from our high-z, high-luminosity sample, while red circles are taken from the literature. These are NGC3783, NGC5548, NGC7469, F9 and 3C390.3 for the Lyα relation, and NGC4395, NGC5548, NGC3783, NGC7469, and 3C390.3 for the CIV relation De Rosa et al., 2015;Lira et al., submitted). Taken from Lira et al. (submitted), c AAS. three sources at high-luminosities come from our lag determinations. Figure 4 clearly supports that R Lyα /R CIV ∼ 1 and that this ratio is independent of luminosity.

SUMMARY AND CONCLUSIONS
We have presented selected results from the RM campaign of 17 high-z, high-luminosity quasars, which lasted more than 10 years. For several sources lags between the continuum and BLR line emission were determined, allowing us to extend radius-luminosity relationships up to λL λ (1345Å) ∼ 10 47 ergs/s. Continuum and line light curves for all sources can be found in Lira et al. (submitted).