Issue 
A&A
Volume 571, November 2014



Article Number  A60  
Number of page(s)  10  
Section  Extragalactic astronomy  
DOI  https://doi.org/10.1051/00046361/201424776  
Published online  07 November 2014 
Imprints of the quasar structure in timedelay light curves: Microlensingaided reverberation mapping
ArgelanderInstitut für Astronomie,
Auf dem Hügel 71,
53121
Bonn,
Germany
email:
dsluse@astro.unibonn.de
Received: 7 August 2014
Accepted: 12 September 2014
The advent of large area photometric surveys has raised a great deal of interest in the possibility of using broadband photometric data, instead of spectra, to measure the size of the broad line region of active galactic nuclei. We describe here a new method that uses timedelay lensed quasars where one or several images are affected by microlensing due to stars in the lensing galaxy. Because microlensing decreases (or increases) the flux of the continuum compared to the broad line region, it changes the contrast between these two emission components. We show that this effect can be used to effectively disentangle the intrinsic variability of those two regions, offering the opportunity to perform reverberation mapping based on singleband photometric data. Based on simulated light curves generated using a damped random walk model of quasar variability, we show that measurement of the size of the broad line region can be achieved using this method, provided one spectrum has been obtained independently during the monitoring. This method is complementary to photometric reverberation mapping and could also be extended to multiband data. Because the effect described above produces a variability pattern in difference light curves between pairs of lensed images that is correlated with the timelagged continuum variability, it can potentially produce systematic errors in measurement of time delays between pairs of lensed images. Simple simulations indicate that timedelay measurement techniques that use a sufficiently flexible model for the extrinsic variability are not affected by this effect and produce accurate time delays.
Key words: gravitational lensing: micro / quasars: general / gravitational lensing: strong
© ESO, 2014
1. Introduction
Gravitational lensing of distant quasars offers new opportunities in the study of quasars and supermassive black holes. When a galaxy happens to be on the line of sight towards a more distant quasar, multiple images of the quasar are created, with a typical separation of a few arcsecs. The (macro) magnification of the lensed images associated with this phenomenon eases the study of those distant objects and allows one to reconstruct a high resolution image of the quasar’s host (Peng et al. 2006; Claeskens et al. 2006; Suyu et al. 2013). Since the stars in the lensing galaxy also act as many individual gravitational lenses, the macro images are in fact composed of many unresolved microimages which are separated on the sky by a few microarcsec corresponding to the Einstein radius η_{0} of the microlenses. The microlensing effect can be used to probe the source at high resolution because microlensing selectively magnifies the source emission as a function of the size of the emitting region, provided the latter is smaller than ~10 η_{0} (Refsdal & Stabell 1997). Nowadays, quasar microlensing is employed to measure the size and temperature profile of the accretion disc, or the size and geometry of the broad line emitting region (Kochanek 2004; Sluse et al. 2007, 2011; Eigenbrod et al. 2008; Morgan et al. 2010; Blackburne et al. 2011; O’Dowd et al. 2011; Guerras et al. 2013).
The quasar continuum emitting region is more compact than η_{0}, and is therefore significantly microlensed. The more extended broad line region (BLR) is generally less affected, with typically 10–20% of its flux being microlensed (Sluse et al. 2012). Consequently, microlensing effectively modifies the contrast between the flux of the continuum and the flux from the broad lines. Because a lensed system is composed of several (two to four) lensed images of the quasar, we observe multiple realizations of the same intrinsic light curve with different amount of microlensing of the continuum and broad lines. The proposed technique of microlensingaided reverberation mapping aims at taking advantage of this effect to measure the time lag between the continuum and broad line variations. Our method is not conceptually very different from photometric reverberation mapping (Haas et al. 2011; Chelouche & Daniel 2012; Chelouche & Zucker 2013; Pozo Nuñez et al. 2012; Edri et al. 2012; Zu et al. 2014; Bachev et al. 2014). In that case, multiband photometry is used to disentangle the flux of the continuum and of the BLR, while our technique can already be applied to singleband data. The same data can now be used to derive the size of the continuum emission based on the variability of the microlensing signal (Kochanek 2004), and if multiband data are available, to measure the temperature profile of the accretion disc (Anguita et al. 2008). This opens the possibility of studying the properties of the accretion disc and of the BLR in the same population of AGNs.
In Sect. 2, we present a fiducial example of lensed quasar light curves that demonstrates that imprints from the continuum and broad line variations are present in the difference light curve between pairs of lensed images under the simplifying assumption of nonvariable microlensing. We explain how we proceed to simulate daily sampled lensed quasar light curves, and show how to use them to measure the time lag between the continuum and the broad line variations. We also discuss how our results depend on the properties of the source, and show that this signal should be detected in timedelay light curves. In Sect. 3, we increase the complexity of the simulations, producing mock light curves with irregular sampling and seasonal gaps, as well as microlensing signal drawn from microlensing simulations. We propose a technique to detrend the time variable microlensing signal with a Bspline model in order to recover the time lag. In the next section (Sect. 4) we search for a bias on the measurement of the time delay in mock light curves, which include continuum and BLR flux. Finally, Sect. 5 summarizes our main results and outlines future work necessary to turn this technique into a robust probe of the quasar structure at intermediate and high redshift.
2. Fiducial case and timelag measurement
In this section, we explain how we simulate realistic light curves of lensed quasars, and test simple techniques to measure the time lag between the continuum and the response from the BLR. We start with the ideal situation where a lensed quasar has been observed for nine years on a daily basis. The duration of this mock light curve is chosen to mimic existing light curves of lensed quasars as obtained by the COSMOGRAIL collaboration (Courbin et al. 2011; Tewes et al. 2013b), and light curves provided by future surveys like the Large Synoptic Survey Telescope (LSST). We first assume a noisefree light curve, and a constant amount of microlensing. Section 3 is dedicated to the simulation of more realistic light curves. Finally, for simplicity, we assume that the time delay from the pair of lensed images has been obtained independently and corrected for, or is naturally close to zero. This situation occurs for lensed systems where the four lensed images have a cross like configuration around the lensing galaxy, e.g. the Einstein Cross Q2237+0305 (Huchra et al. 1985), or where two (three) of the four lensed images are “merging”, a situation that occurs when the source is located close to a “fold” (“cusp”) caustic, e.g. WFI 2033−4723 and RXS J1131−1231 (Morgan et al. 2004; Sluse et al. 2003). Preliminary investigation of the impact of the presence of multiple sources of emission on timedelay measurements is discussed in Sect. 4.
2.1. Intrinsic variability
First, we describe how we simulate the intrinsic variability of the lensed quasar. Despite our limited understanding of the detailed processes governing quasar variability, it has been shown by several authors (e.g. Kelly et al. 2009; Zu et al. 2013; Graham et al. 2014) that Gaussian processes, and in particular damped random walk (DRW), provides a satisfying mathematical description of the AGN variability. Deviations from this model on timescales smaller than five days, as well as possible turnover in the properties of the signal for timescales above 54 days have been suggested (Mushotzky et al. 2011; Zu et al. 2013; Graham et al. 2014). However, those deviations are relatively subtle and overall the DRW process provide a good proxy of the AGN variability. In the following, we generate mock AGN light curves using the JAVELIN code (Zu et al. 2011, 2013, 2014).
The continuum variability is described by: (1)where the mean function of the DRW is (constant over time), and its associated covariance function between two epochs t and t′ is κ(t,t′). Following several authors, we use an exponential covariance function of the form κ(t,t′) = σ^{2}exp(− t − t′  /τ_{d}) where σ^{2} and τ_{d} are the variance and characteristic timescale of the process (Zu et al. 2014; Graham et al. 2014). Note that instead of σ, various variability studies (MacLeod et al. 2010; Butler & Bloom 2011) of AGNs use , which is the amplitude of the DRW. It is related to σ through the relation (Kozlowski et al. 2010). On long timescales, the variance of the light curve is and on short timescales .
The variations in the BLR is modelled as the variation of the continuum convolved with a timelagging transfer function Ψ(t): (2)Following (Zu et al. 2011, 2014), we use a tophat transfer function centered on the time lag τ, with width w and amplitude A, so that (3)Chelouche & Daniel (2012); Chelouche & Zucker (2013) suggest that photometric reverberation mapping using timelag measurements based on crosscorrelation methods are sensitive to the choice of Ψ(t). However, the transfer function is observationally poorly constrained, as reverberation mapping studies generally concentrate on the measurement of the timelag between the continuum and broad line variations, but not on recovering Ψ(t). Since this work is a proofofconcept of the technique, we will only use Ψ(t) as defined in (3), which corresponds to a thin shell geometry of the BLR (Peterson 1993; Pancoast et al. 2011). On the other hand, we also choose the properties of the intrinsic variability light curve to match qualitatively the variability of the lensed quasar RXS J1131−1231 (Sluse et al. 2003, 2006; Tewes et al. 2013b). Future works should investigate in more details whether results are biased by those particular choices.
In the following, the total signal is assumed to be the sum of the continuum and of a single emission line. Furthermore, we fix τ = 100 days. Since the bulk of time lags observed for H β in local AGNs is observed in the range 10–100 days (Kaspi et al. 2000; Bentz et al. 2009), a lag of 100 days may be representative of expected (cosmologically dilated) lags for high ionization UV lines in lensed AGNs at intermediate redshifts.
2.2. Simulating the light curves
The light curves of a pair of lensed images (at a given wavelength and to a good approximation in a given band; assuming that differential extinction due to the lensing galaxy is negligible, as commonly observed) can be expressed as: (4)where F_{1,2} is the flux of the pair of lensed images, M is the relative macromagnification ratio between the pair of images^{1}, μ_{1}(t) (resp. μ_{2}(t)) is the amplitude of microlensing of the most compact region (i.e. the continuum) in image 1 (resp. 2), and F_{Mμ}(t) (F_{M}(t)) the part of the flux that is (not) affected by microlensing. For the fiducial case, all the microlensed flux is emitted by the continuum as defined in (1), hence F_{Mμ}(t) = c(t). On the other hand, the flux that is not microlensed originates from the broad line only, such that F_{M}(t) = f_{BLR}l(t), with f_{BLR} being the flux ratio between the line and the continuum. We should emphasize that in known lensed quasars, a small fraction (in general 10–20%) of the flux from the BLR is microlensed (see Sluse et al. 2012). Therefore, the above identification of the microlensed flux F_{Mμ} to the continuum (only) and of the non microlensed flux to the BLR is an approximation whose impact is discussed in Sect. 3.2.2. In Eq. (4), and in the following, we assume that the macromagnification ratio M = M_{1}/M_{2}, which corresponds to the continuum flux ratio that would be measured in absence of microlensing, has been derived independently from e.g. flux ratios measured from spectroscopy in narrow emission lines, or at longer wavelengths such as midinfrared (MIR) or radio wavelengths (however see Sluse et al. 2013 regarding presence of microlensing in the MIR). For simplicity, we set M = 1.
We start by using f_{BLR} = 0.2. Although arbitrary, this choice may be representative of a large population of quasars. Indeed, the median equivalent width of the main UV emission lines (Mg ii, C iii], C iv) is of the order of 40 Å (Croom et al. 2002; Shen et al. 2011), while the width of red optical filters is typically 120 Å. One may also note that some planned surveys will use narrow band filters (Benitez et al. 2014; Martì et al. 2014), which would significantly increase f_{BLR}. In addition, we assume a constant microlensing in image 1, μ_{1}(t) = 0.5 at all epoch, and no microlensing in image 2. Simulations including more realistic microlensing signal are presented in Sect. 3.2. Figure 1 shows the continuum and emission line light curves in the top panel, and the simulated light curves of the pair of images in the bottom panel. The latter panel also shows in magnitude the differential light curve between images 1 and 2 (solid black line), and the same signal if only the continuum was varying (dashed line). Two important messages have to be kept from this figure. First, the presence of a fraction of the flux (time variable or not) that is not microlensed, leads to differential light curve that does not only contain a signal from microlensing, as commonly assumed, but is modulated at a detectable level (modulation of about ±0.05 mag in Fig. 1) by a signal associated with the quasar intrinsic variability. Second, the differential light curve shows a different shape for a variable and a non variable BLR. In the next section, we cross correlate the differential signal with the intrinsic signal to unveil the imprint of the timelagged signal from the BLR.
Fig. 1
Upper panel: fiducial realization of the continuum (black) and emission line (red) light curves of a quasar (Sect. 2). The light curve corresponding to the emission line has been divided by f_{BLR} to show more clearly the time lag τ ~ 100 days between the continuum and the line. A negligible amount of noise has been considered to ease legibility. Bottom panel: corresponding light curves of the two lensed images in magnitude (thick blue and red lines) and the corresponding difference light curve (solid black). The dotted grey line shows the differential light curve that would be observed if the broad line was not varying, the solid green line a model of the extrinsic variations, and the thick dashed grey line is Â(t). 
2.3. Cross correlation function
Figure 1 demonstrates that the difference light curve between two lensed quasar images is not the same when the contribution associated with the BLR is constant or responds to the continuum variations. The cross correlation function (CCF) of the ratio light curve F_{2}/F_{1} with the microlensed signal F_{1} of the lensed image^{2} displayed in Fig. 2, confirms that variations of the BLR are imprinted in the differential signal. The CCF nicely peaks at 100 days, namely the input time lag τ, but it also shows a pronounced peak at 0 days, as the CCF obtained for a nonvariable emission line. In fact, if one calculates the CCF between F_{2}/F_{1} and F_{2}, this secondary peak becomes the main peak. We discuss in Sect. 2.4 how the linear combination of F_{1} and F_{2} can be used to reduce the power at a zero day lag, and more robustly retrieve the peak at the time lag τ.
Fig. 2
Cross correlation function (normalized to peak intensity) of F_{2}/F_{1} with F_{1} in the case of a timelagged BLR (thin blue line) and of a nonvariable BLR (dashed black). The third curve (thick solid black line) shows the CCF between the fraction of the quasar flux affected by (free of) microlensing F_{Mμ} (F_{M}), as obtained by linearly combining F_{1} and F_{2} (Sect. 2.4, Eq. (5)). The dotteddashed vertical line shows the fiducial time lag τ = 100 days between the variation of the continuum and of the BLR in the simulation. 
2.4. Macromicro Decomposition (MmD)
We describe here a method to deblend the continuum and the BLR signal from the intrinsic light curves. This technique is inspired by a similar method devised in Sluse et al. (2007, see also Hutsemékers et al. 2010 and Sluse et al. 2012) but applied to quasar spectra instead of time series. Simple linear combinations of the signal of a pair of images, as expressed in Eq. (4), allows us to isolate F_{M} using the observed fluxes F_{1} and F_{2}: (5)where we have introduced A(t) = M × μ(t) (with μ(t) = μ_{1}(t) /μ_{2}(t) ≠ 1). This quantity is preferred to M and μ(t) because it is more closely related to observations, with A(t) = F_{1}(t) /F_{2}(t) when there is no flux from the emission lines.
In order to perform this decomposition, it is necessary to know M and to have a proxy to A(t). The acquisition of a spatially resolved spectrum of the lensed images at an epoch t_{1} (ideally part of the photometric monitoring period), with at least the same wavelength range as the broadband data, allows one to derive those two quantities. First, M is obtained based on the flux ratio of the narrow emission lines, or of the broad lines if the latter are at least partly unaffected by microlensing (Sluse et al. 2007, 2012; Hutsemékers et al. 2010; Braibant et al. 2014; Nierenberg et al. 2014). Second, the flux ratio measured in the continuum of the spectra gives us A(t_{1}). Third, it is possible to model F_{1}(t) /F_{2}(t) with a smooth model , which encodes the large scale extrinsic variations of F_{1}(t) /F_{2}(t). For the example depicted in Fig. 1, the solid horizontal green line is such a model. We can then define an empirical estimate of A(t) such that: (6)In general Â(t) ≠ A(t), but the difference may not be large as far as the flux from the continuum is much larger than the flux from the BLR, and microlensing variations remain modest over the time of the monitoring. This is further discussed in Sect. 3.2.
In summary, the MmD allows one to empirically deblend the signal that is microlensed, and mostly originating from the continuum emission, from the signal that is not affected by microlensing. The cross correlation of those two signals, is used in the following to measure the time lag τ. Alternative methodologies may be developed, but we focus in this paper on the use of a zeropadded CCF applied to F_{M} and F_{Mμ} as derived with the MmD.
Fig. 3
Average timelag measurement (left) and associated standard deviation (right) for various fractional flux contribution from the line f_{BLR}, and six different amplitude of microlensing of the continuum μ. For each pair (f_{BLR}, μ), 500 different quasar light curves have been generated, with a time lag τ = 100 days between the variation of the continuum and of the BLR. Note that due to the definition of f_{BLR}, the case f_{BLR} = 0.9 corresponds to an average contribution of the line of 0.9/1.9 = 0.47 to the total flux. 
2.5. Modification of microlensing and BLR contributions
Several properties of the signal might hamper the detectability of a lag, such as the relative contribution of flux from the line, f_{BLR}, the amount of microlensing from the continuum, μ, and large photometric errors. In order to test those effects, we have simulated light curves in the same way as our fiducial light curves (i.e. time delay of 0 days, regular sampling of 1 point per day) for nine different values of f_{BLR} uniformly distributed in the range [ 0.1,0.9 ], and for six values of amplitudes μ, chosen such that −2.5 log (μ) uniformly covers the range [–0.75, 0.75] mag (μ = 1 excluded). First, we consider a noiseless situation. For each couple (f_{BLR}, μ), we have generated 500 different light curves, and for each light curve we have measured the time lag using the peak of the CCF between F_{M} and F_{Mμ} as described in Sect. 2.4. Following this procedure we retrieve a median time lag τ ~ 99 days, with a standard deviation σ_{τ} ~ 2.2 days. The median time lag is 1% smaller than 100 days because the distribution of time lags is asymmetric. The peak of the distribution is in fact found at 100 days. We have not identified the reason for this asymmetry. The transfer function does not seem to be the reason, as we obtained exactly the same average time lag when using a deltafunction for Ψ(t) in (2). A possibility could be lowfrequency variations of the quasar, known to produce similar biases in CCF analysis (Chelouche & Daniel 2012). Because timelag measurements will be affected by larger errors than this bias at the percent level, we do not investigate it further as it will produce only a second order effect.
Fig. 4
Impact of the photometric errors (xaxis) on the distribution of retrieved timelag measurements. Upper panel: μ = 0.5. Three different values of f_{BLR} are shown as squares, stars, and diamonds. Lower panel: f_{BLR} = 0.2. Three different values of μ are shown as stars, squares and diamonds. The dashed black line shows the input time lag τ = 100 days and the solid black line the median lag measured on light curves free of noise. The points with the same photometric uncertainty (i.e. squares, stars, and diamonds) have been slightly shifted between each other towards larger uncertainty to ease legibility. In both panels, the solid blue squares correspond to our fiducial case for different photometric accuracies. 
Finally, we generate a second set of simulations, identical to the above ones, but to which we have added uniform Gaussian noise σ_{f}/f = 0.01. Such a small photometric uncertainty is aimed by upcoming surveys and is reached for good timedelay light curves currently obtained by monitoring programmes such as COSMOGRAIL. Figure 3 shows the results of this procedure. At this noise level, the mean time lag agrees perfectly with the time lag retrieved by noiseless simulations. Uncertainties (at 68.2% confidence level) smaller than ten days can be obtained. However, as expected, the distribution of time lags broadens significantly (as revealed by the larger σ_{τ}) when the amplitude of microlensing is small (i.e. −2.5 log (μ) ~ ± 0.25 mag) or for low values of f_{BLR} (i.e., 0.1). The impact of the photometric accuracy on the results is discussed in the next section.
2.6. Photometric accuracy
The signal to noise of the input light curves is expected to limit the ability to measure τ. For various pairs (f_{BLR}, μ), we generate 500 noiseless simulated light curves. Then, we generate realizations of those light curves for different amounts of noise σ_{f}/f ∈ [ 0.01,0.09 ], and measure the time lag as described earlier. As shown above, in such a situation, the uncertainty in the results depends on the amplitude of microlensing μ and on the relative contribution of flux from the emission line, f_{BLR}. We study two cases: i) we fix f_{BLR} = 0.2, and vary μ = 0.5,0.63,0.79; ii) we fix μ = 0.5 and vary the amplitudes of f_{BLR} = 0.2,0.4,0.7. The median value of τ, and 1σ interval (i.e. interval containing 68.2% of the distribution) are displayed in Fig. 4 as a function of the noise. Three conclusions can be drawn from that figure. First, the median of the distribution agrees with the noiseless case, and is biased low by 1.5% compared to the true time lag, because of the asymmetry of the lag distribution. Second, for an amplitude of microlensing μ = 0.5, corresponding to a demagnification of the continuum by 0.75 mag, the uncertainty on the time lag increases by typically a factor 4 when the photometric uncertainty is increased from 0.01 mag to 0.09 mag. Third, the ability to measure a time lag at low signal to noise strongly depends on the amplitude of microlensing. While a time lag can still be measured in most cases for (μ,f_{BLR}) = (0.79,0.2) if the photometric accuracy is better than 0.03 mag, this measurement becomes highly uncertain for lower signal to noise. We should emphasize that these trends should also depend on the amplitude of intrinsic variability of the quasar (and to some extent of the macro magnification M), and therefore are only representative of the variability properties assumed for our fiducial quasar.
3. Simulations of more realistic light curves
In the previous section, we have demonstrated that light curves of multiple imaged quasars can be used, in presence of microlensing, to perform reverberation measurement of the size of the BLR. The fiducial light curves we presented are however highly idealized. Real light curves will be sampled on a less regular baseline and with gaps between seasons, while microlensing will not be constant but vary with time. In addition, a small amount of microlensing of the emission line may be expected. We consider the impact of all these features hereafter. This may however not cover all possible complications encountered in nature, and future work is needed to quantify how the method behaves when multiple emission lines fall in the same broadband filter, or when emission which is not time variable (e.g. host galaxy flux, narrow emission lines) is present.
3.1. Sampling and gaps in the light curve
Simulating sparsely sampled light curves with seasonal gaps can be performed in a simple way. For a given mock light curve, we first create a new light curve with regular gaps every year to mimic an ensemble of observing seasons, and then we reduce the sampling of each season by keeping a given fraction f of the points per season. In practice, we create light curves with gaps of 130 days and keep only f = 34% of the daily sampled points, which corresponds to a mean frequency of observation of 1 point every three days. This kind of sampling is representative of the best light curves currently obtained for lensed quasars (Vuissoz et al. 2008; Courbin et al. 2011; Tewes et al. 2013b).
Fig. 5
Impact of gaps and sampling on the cross correlation. Upper panel: realistic light curve (black circles) with irregular sampling and seasonal observing gaps, generated based on the fiducial “continuous” light curve (light solid grey). Bottom panel: DCF obtained for the fiducial “continuous” light curve (black band), when data are obtained with a mean sampling of 1 point every three days (green band), when data are obtained with seasonal gaps of 130 days (blue band), and with a sparse sampling and seasonal gaps as shown on the upper panel (red points with error bars). 
The main complication when one works with sparsely sampled light curves comes from the use of the CCF, that cannot be applied to irregularly sampled time series. This problem is common to all reverberation mapping measurements and several techniques have been introduced to address it. We use hereafter the most simple techniques, namely the interpolated cross correlation function (ICCF, Gaskell & Peterson 1987), which consists of applying the cross correlation to the data set after interpolation, and the discrete correlation function (DCF) introduced by Edelson & Krolik (1988). An example of sparsely sampled light curve is shown in Fig. 5, with the corresponding cross correlation, calculated with a bin of three days (i.e. average seasonal sampling of the light curves) for the DCF. The DCF is shown for i) the fiducial “continuous” light curve; ii) with a sampling rate of 1 point every three days; iii) with a daily sampling but seasonal gaps; and iv) in case of sparse sampling and seasonal gaps. A fifth case, with the same number of points as (iv) but no gaps (i.e. a sampling rate of about 1 point every five days) has also been tested but is not shown as the DCF is similar to (ii). In all the cases, the amplitude of the correlation function is oscillating around zero, with a main peak around the true time lag. The DCF obtained with only one third of the points is not very different from the fiducial DCF, showing that the irregular sampling has little impact on the time lag measurement for observing rates of a few days. More critical however is the presence of gaps in the light curves that lead to peaks with distorted shapes. Once light curves with both irregular sampling and seasonal gaps are considered, a secondary peak appears superimposed to the main peak, at τ ~ 220 days. Interestingly, this secondary peak is not visible if we calculate the ICCF instead of the DCF. As we show hereafter, this secondary peak has important consequences for measuring the time lag when seasonal gaps are present.
Instead of generating a large sample of different realizations of the continuum variability c(t), we have repeated the procedure outlined above with different gap locations and different sampling of the input fiducial light curve. We have used the maximum of the correlation function to estimate the time lag. This quantity is well defined for the CCF studied until now, but is more sensitive to noise fluctuations when we calculate the DCF. Therefore, we have compared different techniques to measure the time lag^{3}: 1) we search for the maximum of the correlation function, assuming that the latter is positive and smaller than 1000 days; 2) we fit a Gaussian to the main peak of the correlation function; and 3) we measure the centroid of the correlation function. In the last two cases, we search for the peak after an automatic identification of the main peak, assumed to be the signal located between the two minima of the correlation function for τ ∈ [ −200, + 400 ] days. Figure 5 displays the distribution of lags derived using those three methods for the ICCF and the DCF. The distribution of lags of the ICCF is well described by a Gaussian with a width of ~12 days, centered on the input lag. Hence, the degradation of the observing conditions mostly introduces noise, but does not bias the timelag measurement if the ICCF is used. This is not true when we use the DCF. In that case the distribution becomes broader and multimodal, and the lag can either be biased low or high depending of the method used to measure it. This behaviour is caused by the secondary peak visible in the DCF at τ ~ 220 days. As suggested by Fig. 5, this peak seems to be associated with the gaps in the light curve, but it is not obvious that its location can be predicted a priori. Indeed, this peak is only marginally detected in the DCF of the continuous light curve with gaps (blue band in Fig. 5) but for a smaller lag.
Fig. 6
Distribution of lags τ measured on the fiducial light curve for different realizations of the sampling and location of the gaps. The left panel shows the result for the ICCF and the right panel the results for the DCF. The distributions are shown for three different estimators of the lag (cf. Sect. 3.1): 1) the maximum of the correlation function for lags shorter than 1000 days (thick black); 2) the maximum of a Gaussian function fitted on the main peak; 3) The centroid of the main peak. 
3.2. Time variable microlensing signal
Until now, we have assumed that μ is perfectly known and does not vary with time. Although this assumption may be valid in some realistic situations, the amplitude of microlensing generally varies on timescales of several months to several years (Mosquera & Kochanek 2011). Rapid microlensing variations may be difficult to deblend from intrinsic variability when both variations occur on the same timescale, as observed in Q 01584325 (Morgan et al. 2012). Such a situation is rare however and in most of the known lensed systems, the microlensing signal modulates the intrinsic variability with a lower frequency signal (e.g. Vuissoz et al. 2008; Hutsemékers et al. 2010). In practice, extrinsic variability can often be modelled as a low order polynomial or a spline (Kochanek et al. 2006; Tewes et al. 2013a), corresponding to in (6). In the following, we use spline models of and follow the prescription of Sect. 2.4 to derive F_{M} and F_{Mμ} with the MmD (Eq. (5)). Then, we look whether the cross correlation of F_{M} and F_{Mμ} lead to an accurate estimate of τ.
3.2.1. Generating the microlensing signal
We generate mock light curves similar to the fiducial light curve (i.e. f_{BLR} = 0.2, τ_{fid} = 100 days, σ_{f}/f = 0.01 and same intrinsic signal F_{M},F_{Mμ} as Sect. 2), but now with time variable microlensing. The synthetic microlensing light curve is obtained by drawing 500 random trajectories in two different microlensing magnification maps constructed with the inverse rayshooting code developed by Wambsganss (1990, 2001). We used two maps representative of microlensing occurring in a saddlepoint image with macro magnification μ ~ 20 ((κ,γ) = (0.47,0.57)), and a minimum image characterized by μ ~ 12 ((κ,γ) = (0.42,0.50)). This arbitrary choice, which simulates microlensing for the images A & B−C of the lensed quasar RXS J1131−1231, should be representative of microlensing in many lensed AGNs (Vernardos et al. 2014,their Fig. 2). However, because of the large magnification of the images, the fraction of large amplitude microlensing events over a period of ten years may be larger than commonly observed. In addition, we set the fraction of objects in compact form towards the lensed images to f_{∗} ~ 7%, typical of the stellar fraction at a galactocentric distance of a few effective radii, i.e. where lensed images are commonly located (Mediavilla et al. 2009; Pooley et al. 2012; but see JiménezVicente 2014 who derived f_{∗} ~ 0.2). We assume that the continuum arises from a disc with halflight radius R_{1/2} = 0.06 η_{0} (where η_{0} is the Einstein radius of a microlens), and a track length of 0.45 η_{0}, which corresponds to a transverse velocity of 0.05 η_{0}/ year. Those estimates match expectations for known lensed quasars (Mosquera & Kochanek 2011).
3.2.2. Microlensing models and timelag measurement
We use the PyCS package (Tewes et al. 2013a) to construct, from the pair of simulated light curves, an empirical model of the variability (see Tewes et al. 2013a,Sect. 4). The intrinsic variability signal, common to the pair of light curves, and the extrinsic signal, are simultaneously fitted to the data with separate freeknot spline models^{4}. In addition to the parameters of the spline model, the magnitude shift between the curves and the time delay are free to vary. We employ that technique because it is now commonly used for timedelay measurements (Courbin et al. 2011; Tewes et al. 2013b; Eulaers et al. 2013; Kumar et al. 2013). It provides very good fit to our synthetic light curves, for which we also retrieve our input delay of 0 day. Figure 7 displays various splines reproducing the differential light curve F_{1}/F_{2}, and compares this signal to the input microlensing signal. We compare three splines for the extrinsic variations, differing by their number of “free knots”. The spline with the largest number of knots (32) reproduces best the extrinsic variations but includes variations that are not those of the microlensed continuum. With a lower number of knots (5), the modeling of the light curves is poorer but the model better represents the microlensing variations.
Fig. 7
Two examples of differential light curves and associated spline models . The grey points represent the simulated flux ratio with associated error bars. The thin solid lines show the spline models. The thick solid black curve shows the true microlensing variability A(t) and the dashed black curve shows Â(t) (Eq. (6)) when is a spline with 5 knots. 
Figure 8 shows the distribution of time lags derived with the three spline models described above. The results for the fiducial light curve, namely when microlensing does not vary with time, are also shown. On the left panel, we show the lag as measured using in the MmD, while the right panel assumes that Â(t) (Eq. (6)) has been used. If we focus on the fiducial case, we see that the distribution of lags peaks at τ ~ 150 days when is used. This bias is expected because F_{M} contains a fraction of the microlensed flux. Once we use Â(t), we recover the true lag for the fiducial distribution. The latter is the reference towards which lags derived in presence of time variable microlensing have to be compared. The spline model with 5 knots leads to a distribution of measured lags similar to the fiducial distribution, while more flexible spline models remove a fraction of the intrinsic signal and bias, or even preclude, the lag measurement. When using in the MmD, we find for a significant fraction^{5} of the light curves lags with τ< 50 days or τ> 200 days. Those incorrect lags are derived in two cases: when the microlensing signal varies so quickly that its variations are not adequately reproduced by the spline model, and when both magnification and demagnification occur during the monitoring. In the latter case, the MmD fails once Â(t) ~ M (i.e. μ ~ 1). The use of the DCF instead of the CCF, reduces the weight of those regions in the correlation function, but generally still produces a peak at τ ~ 0 day, or occasionally at τ> 200 days. The correct time lag can be recovered for most of those light curves if the lag at τ ~ 0 day is ignored or if the error bars on F_{M} − F_{Mμ} are artificially increased for the time range where A(t) ~ (1 ± 0.3) M. Similarly, increasing the number of knots in the spline allows for the recovery of the time lag for most of the light curves where F_{1}/F_{2} is poorly approximated by a 5knots spline. When we use Â(t), we find a distribution of lags compatible with the fiducial distribution. Although this procedure properly shifts the first peak of the DCF to the correct lag, it often increases the height of the second peak at τ> 500 days. Therefore, our automatic identification of the peak of the DCF sometimes fails to identify the peak at τ ~ 100 days, artificially decreasing the number of retrieved lags. A visual inspection of the DCF alleviates this problem.
Fig. 8
Distribution of lags τ measured on the fiducial light curve for different realizations of time variable microlensing (cf. Sect. 3.2). Left: distribution of lags τ using (Eq. (6)). Right: distribution of lags using Â(t), and assuming that the lag τ< 210 days as suggested by the left panel. Each panel shows the result for three different spline models of the extrinsic variability, and applying the same procedure to the fiducial case where microlensing is constant over time. 
The distribution of recovered lags does not perfectly match the fiducial distribution, suggesting that a bias at a few percent level might exist in the lag measurement. In addition, we find that time lags are more likely to be correctly retrieved when the average amplitude of microlensing over the monitoring period is larger than typically 0.75 mag. Among all events, time lags are more efficiently measured when the continuum is strongly demagnified.
We have currently assumed that the amplitude of microlensing of the BLR is negligible. However, it has been shown that 10–20% of the flux of broad line is typically microlensed (Sluse et al. 2012; Guerras et al. 2013). To test for this effect, we have assumed that the whole BLR responds to the flux from the continuum, and is microlensed according to its halflight radius size . We have chosen , typical of the size of the BLR in known lensed quasars (Mosquera & Kochanek 2011, their Fig. 3). Because the BLR size is in units of Einstein radius (and therefore depends on the source and lens redshifts) we do not scale the time lag according to the BLR size. Considering only reliable lags (i.e. such that τ ∈ [ 50,200 ] days), we find a distribution of lags in statistical agreement with those measured for a nonmicrolensed BLR.
3.2.3. Summary
The above results demonstrate that it is possible to correct empirically for the time variable microlensing signal without removing the imprinted timelagged signal from the BLR. As expected, time lags are more easily retrieved when microlensing evolves almost linearly over the period of monitoring, as often observed in lensed quasars light curves (Vuissoz et al. 2007, 2008; Courbin et al. 2011; Eulaers et al. 2013). When the microlensing variations are large, some fine tuning may be necessary to detrend the light curves with a spline model. On one hand, the spline used to model extrinsic variations should not be too flexible, since it can then remove signal that is not associated with differential microlensing. On the other hand, regions of the light curves where there is no differential microlensing (i.e. μ ~ 1) lead to an artificial peak at τ ~ 0 day in the DCF. Ignoring this peak, or artificially increasing the error bars of the concerned points in F_{M} and F_{Mμ}, often solves the problem and allows one to identify the peak produced by the lagged BLR signal.
4. Impact on timedelay measurements
The presence of an intrinsic variability signal in the difference light curve between pairs of lensed images could be a source of systematic errors for the measurement of the (stronglensing) time delay between the light curves. To test for this effect, we use PyCS to measure time delays of a set of mock light curves. The measurement of the time delay is a complex problem by itself (see, e.g., Eigenbrod et al. 2005; Tewes et al. 2013a), which depends on the length of the light curve, on the shape and amplitude of the variability signal, on the time sampling, and presence of observing gaps. Therefore, we will not quantify biases possibly taking place depending of these properties but limit ourselves to simple and well controlled situations.
Our intrinsic signal is the same as the fiducial signal, with f_{BLR} = 0.2 and τ = 100 days and Gaussian noise with σ_{f}/f = 0.01. The light curve is shortened to 3099 days, and the observing rate is chosen to be one point every three days on average. Two groups of light curves, one with a slowly varying microlensing signal and one with a large amplitude microlensing, as depicted in Fig. 7, are generated. For each group, three different delays Δt = 20,100, and 200 days, shorter, similar, or larger than the time lag τ, are assumed. We hope this ensemble of light curves to be sufficiently representative of real light curves. For each situation, we generate 500 light curves that have different sampling and noise realizations.
To measure the time delay, we follow the prescriptions of Tewes et al. (2013a), and fit a sum of two freeknot spline models, one for the intrinsic and one for the extrinsic light curve (cf. Sect. 3.2). In addition, we also apply the same technique to a set of reference light curves with the same characteristics as above but containing only flux from the continuum (i.e. there is no flux from the BLR). Following that procedure, we find that timedelay measurements are unaffected by the lagged signal from the BLR, provided the spline used to model the extrinsic variability is flexible enough. A comparison of the three spline models shown in Fig. 7 reveals that a small bias of a fraction of a day may take place when the spline used for microlensing has few knots. This suggests that methods employing an insufficiently flexible microlensing model could suffer from biases due to the quasar structure, in addition to the biases due to the poorly fitted microlensing variability. A detailed comparison of timedelay measuring techniques in presence of quasar structure is beyond the scope of the present paper.
5. Summary and conclusions
Nowadays, high quality photometric monitoring data of gravitationally lensed quasars are obtained for more then 30 systems. In the next decades, owing to the advent of large surveys like the LSST, the time domain will be accessible for an increasing number of astrophysical phenomena, and in particular for quasar and gravitational lensing studies. Current analysis of optical lensed quasars light curves implicitly assume that the flux originates only from the accretion disc. Under this assumption, the differential signal between pairs of light curves shifted by the time delay, yields to the time variable microlensing signal produced by the stars in the lensing galaxy. However, the hypothesis that the entire broadband flux originates from a single emitting region does not hold in general. Broadband quasar emission includes flux from the broad emission lines, and to a lesser extent from narrow emission lines, and other sources of emission (e.g. Balmer continuum, flux from the host galaxy). Because the various emitting regions have different sizes, they are affected differently by microlensing. The latter modifies the contrast between the continuum and the other sources of emission, in particular the broad lines (Fig. 1). Since broad lines respond to continuum variations, their light curve can be cross correlated to that of the continuum to derive their size. This is exactly what is performed by the reverberation mapping technique, which is one of our most powerful probe of the structure of the BLR, and a robust proxy to the mass of the central black hole. That technique, originally designed for spectrophotometric data (e.g. Peterson 1993; Kaspi et al. 2000; Bentz et al. 2009; Denney et al. 2010; Pancoast et al. 2011), is currently extended to multiband photometric data (e.g. Chelouche & Daniel 2012; Zu et al. 2014). The signal produced by microlensing on quasar light curves is conceptually very similar to photometric reverberation mapping, but is potentially applicable to singleband data. In addition, since microlensing provides an independent probe of the accretion disc and BLR size (Kochanek 2004; Eigenbrod et al. 2008; Morgan et al. 2010; Blackburne et al. 2011; Sluse et al. 2011; Garsden et al. 2011; Guerras et al. 2013), it offers a potentially more complete picture of the same sample of objects.
We have studied the modulation of the differential microlensing signal between pairs of lensed quasar light curves when broadband emission originates from two regions: the accretion disc and the BLR. Assuming intrinsic quasar variations of about 0.5 mag, and ~20% of the emission originating from the broad line, we have shown that modulation of the microlensing signal as large as 0.05 mag could be detected for differential microlensing larger than 0.5 mag. We have introduced a technique to combine pairs of lensed quasars light curves which enables one to disentangle the flux of the continuum and of the broad line, provided a spectra of the lensed images have been obtained once during the monitoring to derive the macrolensing flux ratio. The measurement of the time lag τ between the continuum and BLR variations obtained this way, is what we named microlensingaided reverberation mapping. This technique has been applied to several sets of mock light curves, under the simplifying hypothesis of constant microlensing, to test the ability of the method to retrieve τ under various observational conditions. We found that unbiased time lags could be retrieved for fractional flux from the BLR as small as ~10% of the continuum flux, and amplitude of microlensing as small as 0.25 mag. The precision on the timelag measurement depends little on the fraction of flux from the BLR (provided it is typically >10%), but requires sufficiently large amplitude of microlensing, and photometric uncertainties typically better than 0.04 mag, to measure a time lag with small uncertainties. The precision depends significantly on the rate at which data points are observed, and more crucially on the absence/presence of seasonal observing gaps. Gaps and sparse sampling might bias the measurement of τ, depending of the cross correlation technique used. We suggest that more advanced methods, such as the damped random walk models designed by Zu et al. (2011) combined with microlensing simulations, may offer a framework to tackle that problem.
We have tested if time variable microlensing is a strong obstacle to timelag measurements. We found the use of a spline to model the extrinsic variability to be efficient in detrending light curves from microlensing, but some fine tuning in the spline model is often necessary to obtain optimal results. More problematic have been the light curves showing both deamplification and amplification of the continuum during the observational period. In that case, the decomposition method we introduced fails to deblend the continuum and the BLR signal over the monitoring period, and leads to spurious peaks in the cross correlation. Ignoring the time range where there is no differential microlensing between the quasar images allows one to generally solve the problem and recover the correct lag.
Finally, we have performed a preliminary investigation of the impact of the above discussed effect on the measurement of the time delay between lensed quasar images. Multiple techniques exist to measure time delays. We have focused on the one introduced by Tewes et al. (2013a), which uses a freeknots spline function to model the intrinsic and extrinsic variability of quasar light curves. This method, applied to mock light curves with time delays larger, equal or shorter than τ, robustly retrieves the input delay provided the spline modelling the extrinsic variability is flexible enough. In such a situation, the spline function does not only model the microlensing, but also the intrinsic signal superimposed to it. We anticipate that methods which do not account for extrinsic variations in a sufficiently flexible way may lead to a biased estimate of the time delay. A timedelay challenge has been recently set up to test the ability of current techniques to measure accurate time delays from the thousands of lensed quasars light curves that should be monitored with LSST (Dobler et al. 2013; Liao et al. 2014). We suggest that future timedelay challenges account for the effect outlined in this paper.
Microlensingaided reverberation mapping is a promising technique to study the quasar structure from light curves of timedelay lensed quasars up to high redshift. The natural magnification of the lensed images, which happens as a consequence of strong lensing, offers a natural boost of signal to noise. The small image separation, the time delay between the lensed images and the microlensing produced by the stars in the lensing galaxy were in the past a strong limitation in the use of lensed systems to study quasars. Current observational techniques and analysis methods allow one to tackle those difficulties, and use those systems as powerful astrophysical laboratories. While the analysis of the microlensing signal can be used to derive the size and temperature profile of the continuum emission (Anguita et al. 2008; Bate et al. 2008; Floyd et al. 2009; Eigenbrod et al. 2008; Poindexter & Kochanek 2010), the intrinsic variability can be studied in the same systems to derive the size of the BLR. Although, we have conceptually demonstrated the feasibility of this technique, more exhaustive set of simulations need to be carried out in the future. For example, it is necessary to estimate the impact on timelag measurements of subtle effects likely to take place in real data, such as contribution of non variable flux (from an intermediate line region, from the narrow lines and/or from the host galaxy) to the broadband signal. It also has to be seen how those results depend on the absolute value of the time lag, on the properties of the intrinsic variability signal, and on multiple observational constraints that might refrain us from accurately modelling the time variable microlensing signal.
For a given light curve, a proper estimate of the noise associated with the time lag, is discussed in Alexander (1997) and Peterson et al. (1998).
Acknowledgments
D.S. acknowledges support from the Deutsche Forschungsgemeinschaft, reference SL172/11, and M.T. acknowledges support by the DFG grant Hi 1495/21. This research made use of Astropy, a communitydeveloped core Python package for Astronomy (The Astropy Collaboration et al. 2013).
References
 Alexander, T. 1997, in Astronomical Time Series, eds. D. Maoz, A. Sternberg, & E. M. Leibowitz, Astrophys. Space Sci. Lib., 218, 163 [Google Scholar]
 Anguita, T., Schmidt, R. W., Turner, E. L., et al. 2008, A&A, 480, 327 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bachev, R., Strigachev, A., Semkov, E., et al. 2014, Bulg. Astron. J., 2013, 20, 26 [Google Scholar]
 Bate, N. F., Floyd, D. J. E., Webster, R. L., & Wyithe, J. S. B. 2008, MNRAS, 391, 1955 [NASA ADS] [CrossRef] [Google Scholar]
 Benitez, N., Dupke, R., Moles, M., et al. 2014, submitted [arXiv:1403.5237] [Google Scholar]
 Bentz, M. C., Peterson, B. M., Netzer, H., Pogge, R. W., & Vestergaard, M. 2009, ApJ, 697, 160 [Google Scholar]
 Blackburne, J. A., Pooley, D., Rappaport, S., & Schechter, P. L. 2011, ApJ, 729, 34 [NASA ADS] [CrossRef] [Google Scholar]
 Braibant, L., Hutsemékers, D., Sluse, D., Anguita, T., & GarcíaVergara, C. J. 2014, A&A, 565, L11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Butler, N. R., & Bloom, J. S. 2011, AJ, 141, 93 [NASA ADS] [CrossRef] [Google Scholar]
 Chelouche, D., & Daniel, E. 2012, ApJ, 747, 62 [NASA ADS] [CrossRef] [Google Scholar]
 Chelouche, D., & Zucker, S. 2013, ApJ, 2013, 124 [NASA ADS] [CrossRef] [Google Scholar]
 Claeskens, J., Sluse, D., Riaud, P., & Surdej, J. 2006, A&A, 451, 865 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Courbin, F., Chantry, V., Revaz, Y., et al. 2011, A&A, 536, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Croom, S. M., Rhook, K., Corbett, E. A., et al. 2002, MNRAS, 337, 275 [NASA ADS] [CrossRef] [Google Scholar]
 Denney, K. D., Peterson, B. M., Pogge, R. W., et al. 2010, ApJ, 721, 715 [NASA ADS] [CrossRef] [Google Scholar]
 Dobler, G., Fassnacht, C., Treu, T., et al. 2013, ApJ, submitted [arXiv:1310.4830] [Google Scholar]
 Edelson, R. A., & Krolik, J. H. 1988, ApJ, 333, 646 [NASA ADS] [CrossRef] [Google Scholar]
 Edri, H., Rafter, S. E., Chelouche, D., Kaspi, S., & Behar, E. 2012, ApJ, 756, 73 [NASA ADS] [CrossRef] [Google Scholar]
 Eigenbrod, A., Courbin, F., Vuissoz, C., et al. 2005, A&A, 436, 25 [Google Scholar]
 Eigenbrod, A., Courbin, F., Meylan, G., et al. 2008, A&A, 490, 933 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eulaers, E., Tewes, M., Magain, P., et al. 2013, A&A, 553, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Floyd, D. J. E., Bate, N. F., & Webster, R. L. 2009, MNRAS, 398, 233 [NASA ADS] [CrossRef] [Google Scholar]
 Garsden, H., Bate, N. F., & Lewis, G. F. 2011, MNRAS, 418, 1012 [NASA ADS] [CrossRef] [Google Scholar]
 Gaskell, C. M., & Peterson, B. M. 1987, ApJS, 65, 1 [Google Scholar]
 Graham, M. J., Djorgovski, S. G., Drake, A. J., et al. 2014, MNRAS, 439, 703 [NASA ADS] [CrossRef] [Google Scholar]
 Guerras, E., Mediavilla, E., JimenezVicente, J., et al. 2013, ApJ, 764, 160 [NASA ADS] [CrossRef] [Google Scholar]
 Haas, M., Chini, R., Ramolla, M., et al. 2011, A&A, 535, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Huchra, J., Gorenstein, M., Kent, S., et al. 1985, AJ, 90, 691 [NASA ADS] [CrossRef] [Google Scholar]
 Hutsemékers, D., Borguet, B., Sluse, D., Riaud, P., & Anguita, T. 2010, A&A, 519, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 JiménezVicente, J., Mediavilla, E., Kochanek, C. S., & Muñoz, J. A. M. 2014, ApJ, submitted [arXiv:1407.3955] [Google Scholar]
 Kaspi, S., Smith, P. S., Netzer, H., et al. 2000, ApJ, 533, 631 [NASA ADS] [CrossRef] [Google Scholar]
 Kelly, B. C., Bechtold, J., & Siemiginowska, A. 2009, ApJ, 698, 895 [Google Scholar]
 Kochanek, C. S. 2004, ApJ, 605, 58 [Google Scholar]
 Kochanek, C. S., Morgan, N. D., Falco, E. E., et al. 2006, ApJ, 640, 47 [NASA ADS] [CrossRef] [Google Scholar]
 Kozlowski, S., Kochanek, C. S., Udalski, A., et al. 2010, ApJ, 708, 927 [NASA ADS] [CrossRef] [Google Scholar]
 Kumar, S. R., Tewes, M., Stalin, C. S., et al. 2013, A&A, 57, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Liao, K., Treu, T., Marshall, P., et al. 2014, ApJ, submitted [arXiv:1409.1254] [Google Scholar]
 MacLeod, C. L., Ivezić, Ž., Kochanek, C. S., et al. 2010, ApJ, 721, 1014 [NASA ADS] [CrossRef] [Google Scholar]
 Martì, P., Miquel, R., Castander, F. J., et al. 2014, MNRAS, 442, 92 [NASA ADS] [CrossRef] [Google Scholar]
 Mediavilla, E., Munoz, J. A., Falco, E., et al. 2009, ApJ, 706, 1451 [NASA ADS] [CrossRef] [Google Scholar]
 Morgan, N. D., Caldwell, J. A. R., Schechter, P. L., et al. 2004, AJ, 127, 2617 [NASA ADS] [CrossRef] [Google Scholar]
 Morgan, C. W., Kochanek, C. S., Morgan, N. D., & Falco, E. E. 2010, ApJ, 712, 1129 [NASA ADS] [CrossRef] [Google Scholar]
 Morgan, C. W., Hainline, L. J., Chen, B., et al. 2012, ApJ, 756, 52 [NASA ADS] [CrossRef] [Google Scholar]
 Mosquera, A. M., & Kochanek, C. S. 2011, ApJ, 738, 96 [NASA ADS] [CrossRef] [Google Scholar]
 Mushotzky, R. F., Edelson, R., Baumgartner, W., & Gandhi, P. 2011, ApJ, 743, L12 [NASA ADS] [CrossRef] [Google Scholar]
 Nierenberg, A. M., Treu, T., Wright, S. A., Fassnacht, C. D., & Auger, M. W. 2014, MNRAS, 4420, 2434 [NASA ADS] [CrossRef] [Google Scholar]
 O’Dowd, M., Bate, N. F., Webster, R. L., Wayth, R., & Labrie, K. 2011, MNRAS, 415, 1985 [NASA ADS] [CrossRef] [Google Scholar]
 Pancoast, A., Brewer, B. J., & Treu, T. 2011, ApJ, 730, 139 [NASA ADS] [CrossRef] [Google Scholar]
 Peng, C. Y., Impey, C. D., Rix, H., et al. 2006, ApJ, 649, 616 [Google Scholar]
 Peterson, B. M. 1993, PASP, 105, 247 [NASA ADS] [CrossRef] [Google Scholar]
 Peterson, B. M., Wanders, I., Horne, K., et al. 1998, PASP, 110, 660 [NASA ADS] [CrossRef] [Google Scholar]
 Poindexter, S., & Kochanek, C. S. 2010, ApJ, 712, 668 [NASA ADS] [CrossRef] [Google Scholar]
 Pooley, D., Rappaport, S., Blackburne, J. A., Schechter, P. L., & Wambsganss, J. 2012, ApJ, 744, 111 [NASA ADS] [CrossRef] [Google Scholar]
 Pozo Nuñez, F., Ramolla, M., Westhues, C., et al. 2012, A&A, 545, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Refsdal, S., & Stabell, R. 1997, A&A, 325, 877 [NASA ADS] [Google Scholar]
 Shen, Y., Hall, P. B., Richards, G. T., et al. 2011, ApJS, 194, 45 [NASA ADS] [CrossRef] [Google Scholar]
 Sluse, D., Surdej, J., Claeskens, J.F., et al. 2003, A&A, 406, L43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sluse, D., Claeskens, J.F., Altieri, B., et al. 2006, A&A, 449, 539 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sluse, D., Claeskens, J.F., Hutsemékers, D., & Surdej, J. 2007, A&A, 468, 885 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sluse, D., Schmidt, R., Courbin, F., et al. 2011, A&A, 528, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sluse, D., Hutsemékers, D., Courbin, F., Meylan, G., & Wambsganss, J. 2012, A&A, 544, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Sluse, D., Kishimoto, M., Anguita, T., Wucknitz, O., & Wambsganss, J. 2013, A&A, 553, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Suyu, S. H., Auger, M. W., Hilbert, S., et al. 2013, ApJ, 766, 70 [NASA ADS] [CrossRef] [Google Scholar]
 Tewes, M., Courbin, F., & Meylan, G. 2013a, A&A, 553, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tewes, M., Courbin, F., Meylan, G., et al. 2013b, A&A, 556, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 The Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33 [Google Scholar]
 Vernardos, G., Fluke, C. J., Bate, N. F., & Croton, D. 2014, ApJS, 211, 16 [NASA ADS] [CrossRef] [Google Scholar]
 Vuissoz, C., Courbin, F., Sluse, D., et al. 2007, A&A, 464, 845 [Google Scholar]
 Vuissoz, C., Courbin, F., Sluse, D., et al. 2008, A&A, 488, 481 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wambsganss, J. 2001, PASA, 18, 207 [NASA ADS] [CrossRef] [Google Scholar]
 Wambsganss, J. 1990, Ph.D. Thesis, Thesis LudwigMaximiliansUniv., Munich (Germany), Fakultät für Physik [Google Scholar]
 Zu, Y., Kochanek, C. S., & Peterson, B. M. 2011, ApJ, 735, 80 [Google Scholar]
 Zu, Y., Kochanek, C. S., Kozłowski, S., & Udalski, A. 2013, ApJ, 765, 106 [NASA ADS] [CrossRef] [Google Scholar]
 Zu, Y., Kochanek, C. S., Kozlowski, S., & Peterson, B. M. 2014, submitted [arXiv:1310.6774] [Google Scholar]
All Figures
Fig. 1
Upper panel: fiducial realization of the continuum (black) and emission line (red) light curves of a quasar (Sect. 2). The light curve corresponding to the emission line has been divided by f_{BLR} to show more clearly the time lag τ ~ 100 days between the continuum and the line. A negligible amount of noise has been considered to ease legibility. Bottom panel: corresponding light curves of the two lensed images in magnitude (thick blue and red lines) and the corresponding difference light curve (solid black). The dotted grey line shows the differential light curve that would be observed if the broad line was not varying, the solid green line a model of the extrinsic variations, and the thick dashed grey line is Â(t). 

In the text 
Fig. 2
Cross correlation function (normalized to peak intensity) of F_{2}/F_{1} with F_{1} in the case of a timelagged BLR (thin blue line) and of a nonvariable BLR (dashed black). The third curve (thick solid black line) shows the CCF between the fraction of the quasar flux affected by (free of) microlensing F_{Mμ} (F_{M}), as obtained by linearly combining F_{1} and F_{2} (Sect. 2.4, Eq. (5)). The dotteddashed vertical line shows the fiducial time lag τ = 100 days between the variation of the continuum and of the BLR in the simulation. 

In the text 
Fig. 3
Average timelag measurement (left) and associated standard deviation (right) for various fractional flux contribution from the line f_{BLR}, and six different amplitude of microlensing of the continuum μ. For each pair (f_{BLR}, μ), 500 different quasar light curves have been generated, with a time lag τ = 100 days between the variation of the continuum and of the BLR. Note that due to the definition of f_{BLR}, the case f_{BLR} = 0.9 corresponds to an average contribution of the line of 0.9/1.9 = 0.47 to the total flux. 

In the text 
Fig. 4
Impact of the photometric errors (xaxis) on the distribution of retrieved timelag measurements. Upper panel: μ = 0.5. Three different values of f_{BLR} are shown as squares, stars, and diamonds. Lower panel: f_{BLR} = 0.2. Three different values of μ are shown as stars, squares and diamonds. The dashed black line shows the input time lag τ = 100 days and the solid black line the median lag measured on light curves free of noise. The points with the same photometric uncertainty (i.e. squares, stars, and diamonds) have been slightly shifted between each other towards larger uncertainty to ease legibility. In both panels, the solid blue squares correspond to our fiducial case for different photometric accuracies. 

In the text 
Fig. 5
Impact of gaps and sampling on the cross correlation. Upper panel: realistic light curve (black circles) with irregular sampling and seasonal observing gaps, generated based on the fiducial “continuous” light curve (light solid grey). Bottom panel: DCF obtained for the fiducial “continuous” light curve (black band), when data are obtained with a mean sampling of 1 point every three days (green band), when data are obtained with seasonal gaps of 130 days (blue band), and with a sparse sampling and seasonal gaps as shown on the upper panel (red points with error bars). 

In the text 
Fig. 6
Distribution of lags τ measured on the fiducial light curve for different realizations of the sampling and location of the gaps. The left panel shows the result for the ICCF and the right panel the results for the DCF. The distributions are shown for three different estimators of the lag (cf. Sect. 3.1): 1) the maximum of the correlation function for lags shorter than 1000 days (thick black); 2) the maximum of a Gaussian function fitted on the main peak; 3) The centroid of the main peak. 

In the text 
Fig. 7
Two examples of differential light curves and associated spline models . The grey points represent the simulated flux ratio with associated error bars. The thin solid lines show the spline models. The thick solid black curve shows the true microlensing variability A(t) and the dashed black curve shows Â(t) (Eq. (6)) when is a spline with 5 knots. 

In the text 
Fig. 8
Distribution of lags τ measured on the fiducial light curve for different realizations of time variable microlensing (cf. Sect. 3.2). Left: distribution of lags τ using (Eq. (6)). Right: distribution of lags using Â(t), and assuming that the lag τ< 210 days as suggested by the left panel. Each panel shows the result for three different spline models of the extrinsic variability, and applying the same procedure to the fiducial case where microlensing is constant over time. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.