Issue 
A&A
Volume 553, May 2013



Article Number  A121  
Number of page(s)  9  
Section  Catalogs and data  
DOI  https://doi.org/10.1051/00046361/201321140  
Published online  22 May 2013 
COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses^{⋆,}^{⋆⋆}
XII. Time delays of the doubly lensed quasars SDSS J1206+4332 and HS 2209+1914
^{1} Institut d’Astrophysique et de Géophysique, Université de Liège, 17 Allée du 6 Août, Sart Tilman (Bat. B5C), Liège 1, Belgium
email: E.Eulaers@ulg.ac.be
^{2} Laboratoire d’Astrophysique, École Polytechnique Fédérale de Lausanne (EPFL), Observatoire de Sauverny, 1290 Versoix, Switzerland
^{3} Ulugh Beg Astronomical Institute, Uzbek Academy of Sciences, Astronomicheskaya 33, 100052 Tashkent, Uzbekistan
^{4} Instituut voor Sterrenkunde, Katholieke Universiteit Leuven, Celestijnenlaan 200B, 3001 Heverlee, Belgium
^{5} Indian Institute of Astrophysics, II Block, Koramangala, 560 034 Bangalore, India
Received: 21 January 2013
Accepted: 2 March 2013
Aims. Within the framework of the COSMOGRAIL collaboration we present 7 and 8.5yearlong light curves and timedelay estimates for two gravitationally lensed quasars: SDSS J1206+4332 and HS 2209+1914.
Methods. We monitored these doubly lensed quasars in the Rband using four telescopes: the Mercator, Maidanak, Himalayan Chandra, and Euler telescopes, together spanning a period of 7 to 8.5 observing seasons from mid2004 to mid2011. The photometry of the quasar images was obtained through simultaneous deconvolution of these data. The time delays were determined from these resulting light curves using four very different techniques: a dispersion method, a spline fit, a regression difference technique, and a numerical model fit. This minimizes the bias that might be introduced by the use of a single method.
Results. The time delay for SDSS J1206+4332 is Δt_{AB} = 111.3 ± 3 days with A leading B, confirming a previously published result within the error bars. For HS 2209+1914 we present a new time delay of Δt_{BA} = 20.0 ± 5 days with B leading A.
Conclusions. The combination of data from up to four telescopes have led to wellsampled and nearly 9seasonlong light curves, which were necessary to obtain these results, especially for the compact doubly lensed quasar HS 2209+1914.
Key words: gravitational lensing: strong / cosmological parameters / quasars: individual: SDSS J1206+4332 / quasars: individual: HS 2209+1914
Based on observations made with the 1.2m Swiss Euler telescope (La Silla, Chile), the 1.5m AZT22 telescope (Maidanak Observatory, Uzbekistan), the 2.0m HCT telescope (Hanle, India), and the 1.2m Mercator Telescope. Mercator is operated on the island of La Palma by the Flemish Community, at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias.
Numerical values of light curves are only available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/553/A121 and at http://www.cosmograil.org
© ESO, 2013
1. Introduction
Refsdal (1964) was the first to realize that the determination of the time delay between different images of a gravitationally lensed source can be used to derive the Hubble constant H_{0}, the expansion rate of the Universe, if the mass distribution in the lensing galaxy is known. The other way around, if a value for H_{0} is assumed, the time delay between the lensed quasar images allows one to study the mass of the lensing galaxy.
Even if this method presents the advantage of being independent of other calibrations, it has often been subject to controversies on both timedelay values (Press et al. 1992a,b; Pelt et al. 1996; Kundic et al. 1997, see e.g. on QSO 0957+561) and on the Hubble constant H_{0} derived from these time delays, this value seeming lower than the values derived by other techniques (Burud 2001). A decade ago, time delays were measured on relatively short light curves of only a few seasons (e.g. Burud et al. 2002; Patnaik & Narasimha 2001; Corbett et al. 1996), using only one timedelay extracting method at a time. This could lead to biased results and an underestimation of the error bars as pointed out by Eulaers & Magain (2011). Nowadays, the timedelay approach is considered to be competitive and highly complementary to other cosmological probes, and it gives similar results for H_{0} (Oguri 2007; Linder 2011; Suyu 2012).
COSMOGRAIL is a project aimed at measuring time delays in most known gravitationally lensed quasars. The use of 1m to 2m class semidedicated telescopes has allowed us to acquire light curves spanning up to ten years of data. In combination with a homogeneous way of analysing data as described in Tewes et al. (2013a,b) the collaboration is contributing to a significant rise in reliable timedelay measurements of gravitationally lensed quasars. The first results included time delays for SDSS J1650+4251 (Vuissoz et al. 2007), and WFI J20334723 (Vuissoz et al. 2008) and more recently for HE 04351223 (Courbin et al. 2011) and RX J11331231 (Tewes et al. 2013a). In this paper, we present time delays for two doubly lensed quasars: SDSS J1206+4332 () and HS 2209+1914 ().
SDSS J1206+4332 was discovered by Oguri et al. (2005). They found a z = 1.789 quasar whose lensed images are separated by . The redshift of the main lensing galaxy is estimated at z = 0.748. One or two other nearby galaxies might contribute to the lensing potential.
During the Hamburg Quasar Survey Hagen et al. (1999) identified the system HS 2209+1914 as a bright z_{s} = 1.07 quasar having a galaxy neighbour detectable on their Schmidt plates. Nothing else was published on this object until Chantry et al. (2010), who obtained precise astrometry for the two lensed quasar images and the lensing galaxy based on the deconvolution of Hubble Space Telescope images.
In Sect. 2 we summarize the data acquisition and reduction techniques, and briefly explain the methods used to determine the time delays. In Sect. 3, we present the results for the doubly lensed quasar SDSS J1206+4332, including photometry and time delays. In Sect. 4, we take a closer look at HS 2209+1914 for which we give a first timedelay estimate. The last section is a summary of our results.
2. Data acquisition and reduction
2.1. Data acquisition
Both lensed quasars were monitored for more than seven years on different telescopes. The 1.2m Mercator telescope, located at the Roque de los Muchachos Observatory on La Palma, Canary Islands (Spain), is equipped with the MEROPE camera having a field of view of by and a scale of per pixel. COSMOGRAIL observations by the Mercator telescope stopped in December 2008. The remotely operated 2 m Himalayan Chandra Telescope (HCT) started systematic monitoring for the COSMOGRAIL collaboration in 2007 and has continued up to now. It is equipped with the Hanle Faint Object Spectrograph Camera (HFOSC), a 2048 × 4096 pixels CCD camera with a field of view of 10′ by 10′ and a scale of per pixel.
These two telescopes provided the bulk of the data, and were complemented by data from the 1.5m telescope at Maidanak Observatory, Uzbekistan. Most of the data were obtained with the SITE CCD camera, a 800 × 2000 array camera with a scale of per pixel and a field of view of . During two seasons, a second camera (SI) was used with an array of 4096 × 4096 pixels and a field of view of by on the same scale of per pixel.
For the lensed system HS 2209+1914, some data were added from the 1.2m Euler telescope, located at La Silla, Chile. It is equipped with a 2048 × 2048 CCD camera having a field of view of 11′ by 11′ and a scale of per pixel.
At every observing epoch, and as long as there were no technical or meteorological problems, five dithered images were taken in the Rband with an exposure time of six minutes each.
2.2. Data reduction
To analyse COSMOGRAIL data in a consistent and homogeneous way, we use the semiautomated pipeline that was described in Tewes et al. (2013a) to carry out the reduction of the CCD images. We summarize the main points after the bias subtraction and flatfielding of the frames.
Four stars are used to construct a point spread function (PSF) per frame. Photometry of the sources is obtained through simultaneous deconvolution using the MCS algorithm (Magain et al. 1998), where the main feature is the finite final resolution, which avoids artefacts. The deconvolved images have a pixel scale of half the telescope’s pixel scale and have a Gaussian PSF with a two pixel fullwidthathalfmaximum. The flux of the point sources is allowed to vary from one frame to another, unlike the lensing galaxy, which is part of the numerical background that is held fixed throughout the frames. The simultaneous deconvolution of all frames constrains the positions of the quasar images very well, even if not all data have been obtained under optimal seeing conditions.
An estimation of the photometric shot noise σ_{N}, calculated for each quasar image and frame, is given by^{1}(1)with f_{∗} the flux in the quasar image (expressed in number of photons), N_{sky} the sky level of the exposure, R the CCD read noise (both in photons per pixel), and S the PSF sharpness given by (2)i.e. the sum over all pixels of the squared normalized fluxes of the PSF pixels.
The data points in the light curves are the median photometry per night obtained by the simultaneous deconvolution of the images. Their error bars are the maximum of two kinds of error bars: an empirical error bar and an intrinsic one. Indeed, the total error on the night cannot be smaller than the one given by the spread of the measurements in that night, nor can it be smaller than the error bar given by the photon noise combined with the renormalization error, so taking the maximum of both is justified. They are generally very close to one another.
The empirical error reflecting the spread of the measurements during one night is estimated by the standard deviation of the measurements per night normalized by the number of images per night. Since our photometry is a median value per night, the standard deviation is estimated via the median absolute deviation (MAD) estimator (Hoaglin et al. 1983) corrected by a scale factor for normal distributions. The MAD is defined as the median of the absolute values of the residuals from the data’s median: (3)It is related to the standard deviation via σ = K ∗ MAD. In case of a normal distribution, , with φ^{1} being the quantile function of the normal distribution.
On the other hand, the intrinsic or theoretical error, calculated as the quadratic sum of the median shotnoise per night and the median renormalization error per night, gives a theoretical estimate of the error on the combined data point. That neither the empirical nor the theoretical error systematically dominate is a proof of the coherence of both approaches, and it justifies the idea of taking the maximum of both.
Finally, we put together the light curves of the same object coming from different telescopes by optimizing the magnitude shifts between the curves of the same image from different instruments.
Overview of the optical monitoring for SDSS J1206+4332.
2.3. Timedelay analysis
To avoid a bias due to the method, we use four different techniques for our timedelay analysis. Three of them, a spline fit, a dispersion minimization, and a method based on the variability of regression differences, are described and tested extensively in Tewes et al. (2013b). The fourth method, the numerical model fit, was described and applied in Eulaers & Magain (2011). We quickly summarize the main features of each of them:
 1.
The dispersionlike technique considers the time delay to be the shift between the light curves that minimizes the dispersion between these data. Microlensing is modelled by polynomials up to the third order, but no model light curve for the overall quasar variability is constructed here.
 2.
The regression difference technique searches for the time delay that minimizes the variability of the differences between the timeshifted continuous regressions of the light curves. Microlensing variability is not explicitly modelled here.
 3.
The free knot spline technique models the intrinsic quasar variability, as well as the microlensing variations, by fitting splines simultaneously to the light curves. The curves are shifted in time to optimize this fit.
 4.
The numerical model fit (NMF) constructs a numerical model for the quasar variability, together with a linear microlensing trend for a given time delay. The optimal time delay is the one that minimizes the difference between the data and this numerical model.
We use a Monte Carlo approach for the timedelay uncertainty calculations. For the first three methods, these are based on simulated light curves with known time delays based on the spline fit as explained in Tewes et al. (2013b). For the last method, we use the numerical model light curve to which we add appropriate Gaussian noise.
3. SDSS J1206+4332
3.1. Deconvolution results and light curves
For SDSS J1206+4332 we use data coming from three telescopes. At Maidanak this lensed quasar was observed from January 2005 until July 2008. During the 2007 season, the second camera SI was used. Mercator data span the period March 2005 till December 2008, and HCT has continued monitoring this lens since May 2007 until now. The combination of these data results in seven complete observing seasons from 2005 until 2011.
The object is unvisible from August until the beginning of November, a period of more or less 100 days. The mean temporal sampling during the visibility window is about one point every fourth night for Mercator and Maidanak data, and one point per week for HCT. In the overlapping 2008 season this gives more than one point every second night. Table 1 gives an overview of the optical monitoring for this quasar. One epoch corresponds to one data point in the light curve, which is the median value for that night as explained in Sect. 2. S is the temporal sampling of the observations expressed as an approximative average in number of days excluding the seasonal gaps.
Astrometry for SDSS J1206+4332 relative to the brighter quasar image A.
The four stars, designated 1 to 4 in Fig. 1, are used to model the PSF needed for the simultaneous deconvolution. If these stars did not satisfy our quality criteria (due to e.g. cosmic ray hits), a second set of four stars was used: star 1 in combination with star 5 to star 7. Five stars, 1 to 5, were also used as flux calibrators.
Fig. 1 Part of the field around SDSS J1206+4332. It is a combination of 29 frames taken at the 1.2m Mercator telescope with a seeing and an ellipticity ≤ 0.12. The stars used to model the PSF are labelled 1 to 4. In case one of these stars did not have sufficient quality, three other stars labelled 5 to 7 were used. North is up and east is left. 
Figure 2 shows that the simultaneous deconvolution of 1109 Rband Mercator images clearly reveals the lens galaxy between the lensed quasar images, as well as two other galaxies close to the system, which have also been identified by Oguri et al. (2005). Table 2 presents the relative astrometry of the system as determined from this deconvolution. It shows a slightly different position for the B image and the galaxy G1 in comparison with the positions in Oguri et al. (2005).
Fig. 2 Result of the simultaneous deconvolution of 1109 Rband Mercator images of J1206+4332. The pixel size is (i.e. half the pixel size of the Merope CCD camera at the Mercator telescope), and the resolution is 2 pixels fullwidthhalfmaximum. In addition to the main lensing galaxy G1 between the lensed images A and B we clearly identify two more galaxies: G2 to the north of the system and G3 to the east. North is up and east is left. 
We obtain an image separation of between the quasar images A and B, which is slightly higher than Oguri et al. (2005), but a lower value for the flux ratio F_{A}/F_{B} = 1.31, for which Oguri et al. (2005) found 1.48. Our value is the result of the deconvolution and has not been corrected for the time delay yet. We come back to this in the next section.
The light curves that have been obtained as described in Sect. 2 are shown in Fig. 3. The 2007 and 2008 seasons include data from all three telescopes and show that they match and complement each other well. We distinguish clear features of quasar variability over periods as short as two months, for example in the 2006 season of the B curve. Visually, the light curves are compatible with a time delay around ∽100 days. Finally, we also observe an evolution of the flux ratio between the quasar images, especially in the 2010 and 2011 seasons.
3.2. Time delay
Based on their lens modelling, Oguri et al. (2005) predicted a time delay for SDSS J1206+4332 of Δt_{AB} = 92.6 or 104.4 days depending on the influence of the secondary galaxy G2. With only one observing season, Paraficz et al. (2009) found a time delay of days.
Our timedelay analysis is based on the light curves presented in Fig. 3, and it uses the four different methods briefly explained in Sect. 2. Since three of these techniques (dispersion, regression difference and spline fit) rely on optimizations of a first guess of the time delay, we test the robustness of the result by running 200 times these techniques on the same data set, while uniformly and randomly varying the initial delay in a range of ± 10 days around the first guess. Figure 4 shows that in the case of SDSS J1206+4332 we have monomodal distributions for the spline and regression techniques, and a bimodal distribution for the dispersion technique, which is inherently more sensitive to the variation of this entry delay. The results of this test confirm that our method parameters are well adapted to the data, so that we can use the mean value of these distributions as our final best time delay.
Summary of time delays for SDSS J1206+4332 obtained by the four different methods.
To obtain realistic error bars, we applied these three methods to 1000 synthetic light curves with known time delays. The error bars for these three techniques consist of a systematic and a random part, which then are summed quadratically to obtain the total error. Figure 5 illustrates both contributions: we plot the mean measurement error as a function of the true delay. The error bars are the standard deviation of these measurement errors. To be conservative, the maximum measurement error is our systematic error, and the maximum standard deviation our random error, which are then summed quadratically.
The numerical model fit (NMF) is completely independent of the three other methods in how it determines the time delay and its associated error bars. We add normally distributed random errors with the appropriate standard deviation to the numerical model light curve and redetermine the time delay. This procedure is repeated 1000 times. The mean value of the timedelay distribution that we obtain in that way is considered to be the final time delay, and its dispersion represents the 1σ error bar. Because the errors added to the model, hence the final error bars on the time delay, depend on the photometric error bars of our data, we checked that our residuals (data minus model) were compatible with a normalized Gaussian distribution.
One of the parameters to be optimized in the NMF method is the flux ratio corrected by the time delay, for which we find F_{A}/F_{B} = 1.45 for the complete data set, which is very close to the Oguri et al. (2005) value of 1.48.
We present an overview of the time delays, along with 1σ error bars, obtained by the different methods in Table 3, first for the Mercator and HCT data, two sets that really complement each other, and then for the whole data set coming from the three telescopes. Adding Maidanak data to the Mercator and HCT light curve adds more scatter, but seems to constrain the model light curves better, which is reflected in larger error bars for the dispersion method that does not involve a model, but smaller ones for the regression difference technique, the spline fit, and the numerical model fit. Image A leads image B by Δt_{AB} days. Our results coming from four basically very different timedelay estimation techniques are very consistent, and confirm the existing estimation from Paraficz et al. (2009), while improving its error bars for the spline fit and the NMF techniques used on the complete data set.
That the error bars for the dispersion technique, the regression difference, and the spline fit are systematically larger than for the NMF is inherent to the way they are calculated. For these first three methods, the error bars are the quadratic sum of two maximum contributions, as described before, whereas the NMF does not take the systematic error into account (generally under one day), and considers the standard deviation of the time delay to be its random error without taking a maximum over models with slightly different time delays. On top of that, the random error contribution of the first three methods is the result of red noise, whereas the NMF only takes white noise into account.
Since the time delays resulting from the four methods all converge towards the same value, we only have to decide on the size of the error bars, which reflect the different sensitivities to various aspects that are present in the data and the method. Given the low spread in timedelay values between the methods, and given the possibility that the error bars from the dispersion methods, the regression difference technique, and the spline fit have been overestimated, while those for the NMF might be underestimated, we think that a final time delay for SDSS J1206+4332 of Δt_{AB} = 111.3 ± 3 is a reasonable and comprehensive value, which can be applied in lens modelling.
Fig. 3 Light curves of SDSS J1206+4332 combining data from Mercator, Maidanak, and HCT, and spanning 7 observing seasons. Photometry is obtained from simultaneous deconvolution, and the 1σ error bars are calculated as described in Sect. 2. 
4. HS 2209+1914
4.1. Deconvolution results and light curves
The doubly lensed quasar HS 2209+1914 was observed by four telescopes involved in COSMOGRAIL: the Euler telescope started observations in May 2004 and was joined by the Mercator telescope, which continued the monitoring until December 2008. In May 2007, observations of this lens also started at the Maidanak Observatory during two seasons, and they were joined by the HCT the same year. This last telescope has continued observations until now. Even without the end of the 2012 season, we dispose of a total of 8.5 seasons of observations from 2004 to mid2012.
The lensed quasar can be observed from May to December suffering from a nonvisibility window of four months. The temporal sampling varies with the seasons and telescopes: combining data from Mercator, HCT, and Maidanak yields the best sampling of more than one point every second night in 2007 and 2008. This drops to one point every fortnight on average in 2010, 2011, and 2012 when we only dispose of HCT data. Table 4 gives an overview of the optical monitoring for this lens. One epoch corresponds to one data point in the light curve, which is the median value for that night, as explained in Sect. 2.
Fig. 4 Distribution of time delays for SDSS J1206+4332 
Fig. 5 Error analysis of the simulated light curves for SDSS J1206+4332 using data from Mercator, HCT, and Maidanak. We plot binned mean measurement errors and their standard deviations against the true delays of the synthetic curves. The systematic error is largest for the dispersion technique, and random errors dominate for the regression and spline methods. 
Figure 6 shows the four stars, labelled 1 to 4, used to model the PSF for the simultaneous deconvolution. Three of these stars, 1 to 3, were also used as flux calibrators.
In order not to degrade our photometry by some bad exposures, we eliminated all images with a seeing >3″ or an ellipticity of the point sources >0.3, as well as images with too high a sky level or deviating normalization coefficient. Simultaneous deconvolution of all remaining 2242 images reveals a faint ringlike structure (see Fig. 7) that is very similar to the deconvolved Hubble Space Telescope images shown in Chantry et al. (2010). A secondary galaxy to the west of the system is also visible.
However, the main lensing galaxy cannot be resolved, given the compactness of the system. Both the image separation of only and the flux ratio after deconvolution of (F_{B}/F_{A}) = 0.82 agree with Chantry et al. (2010), who find and 0.79 ± 0.027, respectively, even if our flux ratio has not been corrected for the time delay yet.
Applying the data reduction described in Sect. 2, we show the complete light curve in Fig. 8. Euler and Mercator data complement each other well in the 2004 and 2005 seasons. The 2007 and 2008 seasons show a very good match between the data from the Mercator, Maidanak, and HCT telescopes. Even if we can observe brightness variations in the light curve up to 0.2 mag over the years, clear features within a season are unfortunately absent, which complicates the timedelay measurement. However, the slight rise in magnitude at the end of 2009 in the B curve pins down the order, with the fluctuations in the B image preceding those observed in A.
Fig. 6 Part of the field of view around HS 2209+1914. It is a combination of 49 frames taken at the Mercator telescope with a seeing . The stars used to model the PSF are labelled 1 to 4. All but the last one are also used to calculate the normalization coefficient. North is up and east is left. 
Overview of the optical monitoring for HS 2209+1914.
Fig. 7 The numerical background common to the simultaneous deconvolution of 2242 Rband images of HS 2209+1914 reveals a faint ringlike structure that is very similar to the deconvolved Hubble Space Telescope images shown in Chantry et al. (2010). A secondary galaxy to the west of the system is also visible. However, the main lensing galaxy cannot be resolved given the compactness of the system. North is up and east is left. 
4.2. Time delay
No theoretical estimations or previous measurements of the time delay are available for this system. Our analysis is based on the light curves shown in Fig. 8, and it makes use of the four techniques of which we recalled the basic principles in Sect. 2. As we explained in more detail in the previous section on SDSS J1206+4332, we first test the sensivity of the results to the initial guess of a timedelay value necessary for the dispersion, regression, and spline fit methods. In Fig. 9 we see that the regression difference technique and, to a lesser extent, the dispersion method show a second solution. We take this into account because our final best time delay is the weighted mean of these different values.
The error bars for the dispersion, regression difference, and spline fit techniques are calculated in the same way as in the previous section. The analysis of the random and systematic errors, based on 1000 synthetic curves with known time delays, and presented in Fig. 10, show very small systematic errors in comparison to the large random errors.
Applying the numerical model fit to our light curve constitutes an independent test for the timedelay results. To have realistic error bars on our time delay, we have to check, as previously explained, whether the residuals from our fit are compatible with a normal Gaussian distribution with appropriate mean and sigma. This is not the case for HS 2209+1914, so we multiply the photometric error bars of the data by the appropriate positive factor to obtain a correct distribution. As mentioned in Tewes et al. (2013b), a rescaling of these error bars does not directly influence the uncertainties on the timedelay measurement for the dispersion, regression difference, and spline fit techniques, but it does of course for the NMF method as the noise added to the model light curve directly depends on the photometric error bars of the data points. Larger photometric errors lead to a larger uncertainty on the time delay obtained by the NMF.
According to the telescopes involved, the NMF method finds a flux ratio corrected by the time delay ranging between F_{B}/F_{A} = 0.82 (Maidanak) and F_{B}/F_{A} = 0.85 (Mercator) and of F_{B}/F_{A} = 0.84 when applied to the complete data set. These values are only slightly higher than the ratio obtained after deconvolution mentioned previously.
We present an overview of our timedelay results from the four different methods on the complete data set in Table 5. Image B leads image A by Δt_{BA} days. We notice that in spite of very large error bars, especially for the dispersion and regression difference techniques, the different timedelay estimators do converge to a common time delay of Δt_{BA} ≈ 20 days with image B leading A.
Fig. 8 Light curves of HS 2209+1914 combining data from Mercator, HCT, Euler, and Maidanak and spanning nearly 9 observing seasons. Photometry is obtained from simultaneous deconvolution, and the 1σ error bars are calculated as described in Sect. 2. 
Summary of time delays for HS 2209+1914 obtained by the four different methods on the complete light curve with data from Euler, Mercator, Maidanak, and HCT.
The large error bars for two of our methods can be partially explained if we take a closer look at the contributions of the data from different telescopes to the time delay. We applied the NMF to five subsets of data: 1) Mercator only; 2) Maidanak only; 3) HCT only; 4) Mercator + HCT; 5) Euler + Mercator + HCT. Comparing the results of a single run of the method on these subsets and a Monte Carlo approach of 1000 runs on perturbed model curves, both presented in Table 6, where image B leads image A by Δt_{BA} days, reveals interesting information. Mercator and HCT data do not converge to the same time delay: Mercator data point to one around ≈23 days, whereas HCT data allow a possible delay around ≈33 days. With Maidanak data, which cover the two overlapping seasons between Mercator and HCT exactly, we find both delays. With bare Maidanak data, the lower Mercator value is preferred, but once we perturb the model curve, we find a bimodal distribution with peaks around ≈22 and ≈34 days. The more data we add to the light curve, the more the lower delay dominates.
Fig. 9 Distribution of time delays for HS 2209+1914 obtained by running three of the timedelay estimation techniques on the complete data set 200 times, while uniformly and randomly varying the initial delay in a range of ± 10 days around our first guess for the time delay. 
Fig. 10 Error analysis of the simulated light curves for HS 2209+1914 using all data from Euler, Mercator, Maidanak and HCT. We plot binned mean measurement errors and their standard deviations against the true delays of the synthetic curves. The random errors largely dominate the systematic ones, and are the smallest for the spline fit method. 
The remarks we made at the end of Sect. 3.2 on the comparison of the results are also valid for the time delay in HS 2209+1914. However, the presence of a possible second solution in the data and the question of whether the error bars should reflect this presence complicate the picture. Again, basically different methods converge towards the same value of Δt_{BA} ≈ 20 days, so a final result of Δt_{BA} = 20 ± 5 is a reasonable and intermediate solution between possibly too optimistic and too pessimistic error bars, while not completely excluding higher values since these uncertainties are 1σ error bars. However, given the lack of obvious fast intrinsic quasar variability features in these light curves, and the possible presence of microlensing, this timedelay measurement should be handled with caution.
Summary of time delays for HS 2209+1914 obtained by the NMF on five subsets of data and on the complete set.
5. Conclusions
In this paper we present sevenseasonlong light curves of two doubly lensed quasars based on COSMOGRAIL observations combining data from Mercator, Maidanak, and HCT telescopes for SDSS J1206+4332, as well as with data from Euler for HS 2209+1914. All data were reduced in a homogeneous way described in Tewes et al. (2013a). We applied four timedelay estimation methods, based on different principles, on these light curves.
For SDSS J1206+4332 simultaneous deconvolution reveals the three galaxies that possibly contribute to the lensing potential of the system. We consider the results of the spline fit method as our final best time delay of 111 days and pin down the error bars to 3.5%, thus improving the uncertainty over a previously published result (Paraficz et al. 2009).
Hardly anything has been published on the lensed quasar HS 2209+1914 since its discovery. Our work confirms the image separation, the flux ratio, and a faint ringlike structure as mentioned in Chantry et al. (2010). We are the first to measure a time delay of Δt_{BA} ≈ 20 days, which was only possible by combining data from four telescopes into a 8.5yearlong light curve.
Even if quadruply lensed quasar systems allow a more straightforward estimation of the Hubble constant H_{0}, Suyu (2012) has shown that doubly lensed systems can augment the number of useful time delays when high resolution images of the source galaxy in extended arcs can break the radial profile slope degeneracy. HS 2209+1914 is promising in this respect even if the compactness of the system and the short time delay make it difficult to reach sufficiently high precision.
Acknowledgments
We are deeply grateful to the numerous observers who contributed to the data acquisition at the FlemishBelgian 1.2m Mercator telescope (La Palma, Spain), the Swiss 1.2m Euler telescope (La Silla, Chile), the Uzbek 1.5m AZT22 telescope (Maidanak, Uzbekistan), and the Indian 2.0m HCT telescope (Hanle, India). COSMOGRAIL is financially supported by the Swiss National Science Foundation (SNSF). This work was partially supported by ESA and the Belgian Federal Science Policy (BELSPO) in the framework of the PRODEX Experiment Arrangement C90312.
References
 Burud, I. 2001, Ph.D. Thesis, Institute of Astrophysics and Geophysics, Liège [Google Scholar]
 Burud, I., Courbin, F., Magain, P., et al. 2002, A&A, 383, 71 [Google Scholar]
 Chantry, V., Sluse, D., & Magain, P. 2010, A&A, 522, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Corbett, E. A., Browne, I. W. A., Wilkinson, P. N., & Patnaik, A. 1996, in Astrophysical Applications of Gravitational Lensing, eds. C. S. Kochanek, & J. N. Hewitt, IAU Symp., 173, 37 [Google Scholar]
 Courbin, F., Chantry, V., Revaz, Y., et al. 2011, A&A, 536, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eulaers, E., & Magain, P. 2011, A&A, 536, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hagen, H.J., Engels, D., & Reimers, D. 1999, A&AS, 134, 483 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hoaglin, D. C., Mosteller, F., & Tukey, J. W. 1983, Understanding robust and exploratory data analysis (New York: Wiley) [Google Scholar]
 Linder, E. V. 2011, Phys. Rev. D, 84, 123529 [NASA ADS] [CrossRef] [Google Scholar]
 Magain, P., Courbin, F., & Sohy, S. 1998, ApJ, 494, 472 [NASA ADS] [CrossRef] [Google Scholar]
 Oguri, M. 2007, ApJ, 660, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Oguri, M., Inada, N., Hennawi, J. F., et al. 2005, ApJ, 622, 106 [NASA ADS] [CrossRef] [Google Scholar]
 Paraficz, D., Hjorth, J., & Elíasdóttir, Á. 2009, A&A, 499, 395 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Patnaik, A. R., & Narasimha, D. 2001, MNRAS, 326, 1403 [NASA ADS] [CrossRef] [Google Scholar]
 Refsdal, S. 1964, MNRAS, 128, 307 [NASA ADS] [CrossRef] [Google Scholar]
 Suyu, S. H. 2012, MNRAS, 426, 868 [NASA ADS] [CrossRef] [Google Scholar]
 Tewes, M., Courbin, F., Meylan, G., et al. 2013a, A&A, submitted [arXiv:1208.6009] [Google Scholar]
 Tewes, M., Courbin, F., & Meylan, G. 2013b, A&A, in press (Paper XI) [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]
All Tables
Summary of time delays for SDSS J1206+4332 obtained by the four different methods.
Summary of time delays for HS 2209+1914 obtained by the four different methods on the complete light curve with data from Euler, Mercator, Maidanak, and HCT.
Summary of time delays for HS 2209+1914 obtained by the NMF on five subsets of data and on the complete set.
All Figures
Fig. 1 Part of the field around SDSS J1206+4332. It is a combination of 29 frames taken at the 1.2m Mercator telescope with a seeing and an ellipticity ≤ 0.12. The stars used to model the PSF are labelled 1 to 4. In case one of these stars did not have sufficient quality, three other stars labelled 5 to 7 were used. North is up and east is left. 

In the text 
Fig. 2 Result of the simultaneous deconvolution of 1109 Rband Mercator images of J1206+4332. The pixel size is (i.e. half the pixel size of the Merope CCD camera at the Mercator telescope), and the resolution is 2 pixels fullwidthhalfmaximum. In addition to the main lensing galaxy G1 between the lensed images A and B we clearly identify two more galaxies: G2 to the north of the system and G3 to the east. North is up and east is left. 

In the text 
Fig. 3 Light curves of SDSS J1206+4332 combining data from Mercator, Maidanak, and HCT, and spanning 7 observing seasons. Photometry is obtained from simultaneous deconvolution, and the 1σ error bars are calculated as described in Sect. 2. 

In the text 
Fig. 4 Distribution of time delays for SDSS J1206+4332 

In the text 
Fig. 5 Error analysis of the simulated light curves for SDSS J1206+4332 using data from Mercator, HCT, and Maidanak. We plot binned mean measurement errors and their standard deviations against the true delays of the synthetic curves. The systematic error is largest for the dispersion technique, and random errors dominate for the regression and spline methods. 

In the text 
Fig. 6 Part of the field of view around HS 2209+1914. It is a combination of 49 frames taken at the Mercator telescope with a seeing . The stars used to model the PSF are labelled 1 to 4. All but the last one are also used to calculate the normalization coefficient. North is up and east is left. 

In the text 
Fig. 7 The numerical background common to the simultaneous deconvolution of 2242 Rband images of HS 2209+1914 reveals a faint ringlike structure that is very similar to the deconvolved Hubble Space Telescope images shown in Chantry et al. (2010). A secondary galaxy to the west of the system is also visible. However, the main lensing galaxy cannot be resolved given the compactness of the system. North is up and east is left. 

In the text 
Fig. 8 Light curves of HS 2209+1914 combining data from Mercator, HCT, Euler, and Maidanak and spanning nearly 9 observing seasons. Photometry is obtained from simultaneous deconvolution, and the 1σ error bars are calculated as described in Sect. 2. 

In the text 
Fig. 9 Distribution of time delays for HS 2209+1914 obtained by running three of the timedelay estimation techniques on the complete data set 200 times, while uniformly and randomly varying the initial delay in a range of ± 10 days around our first guess for the time delay. 

In the text 
Fig. 10 Error analysis of the simulated light curves for HS 2209+1914 using all data from Euler, Mercator, Maidanak and HCT. We plot binned mean measurement errors and their standard deviations against the true delays of the synthetic curves. The random errors largely dominate the systematic ones, and are the smallest for the spline fit method. 

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.