Issue 
A&A
Volume 580, August 2015



Article Number  A38  
Number of page(s)  10  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201423977  
Published online  23 July 2015 
H_{0} from ten wellmeasured time delay lenses
^{1} Indian Institute of Astrophysics, II Block, Koramangala, 560 034 Bangalore, India
email: rathna@iiap.res.in
^{2} Physical Research Laboratory, Navrangpura, 380 009 Ahmedabad, India
Received: 10 April 2014
Accepted: 27 May 2015
In this work, we present a homogeneous curveshifting analysis using the differencesmoothing technique of the publicly available light curves of 24 gravitationally lensed quasars, for which time delays have been reported in the literature. The uncertainty of each measured time delay was estimated using realistic simulated light curves. The recipe for generating such simulated light curves with known time delays in a plausible range around the measured time delay is introduced here. We identified 14 gravitationally lensed quasars that have light curves of sufficiently good quality to enable the measurement of at least one time delay between the images, adjacent to each other in terms of arrivaltime order, to a precision of better than 20% (including systematic errors). We modeled the mass distribution of ten of those systems that have known lens redshifts, accurate astrometric data, and sufficiently simple mass distribution, using the publicly available PixeLens code to infer a value of H_{0} of 68.1 ± 5.9 km s^{1} Mpc^{1} (1σ uncertainty, 8.7% precision) for a spatially flat universe having Ω_{m} = 0.3 and Ω_{Λ} = 0.7. We note here that the lens modeling approach followed in this work is a relatively simple one and does not account for subtle systematics such as those resulting from lineofsight effects and hence our H_{0} estimate should be considered as indicative.
Key words: gravitational lensing: strong / methods: numerical / cosmological parameters / quasars: general
© ESO, 2015
1. Introduction
The Hubble constant at the present epoch (H_{0}), the current expansion rate of the universe, is an important cosmological parameter. All extragalactic distances, as well as the age and size of the universe depend on H_{0}. It is also an important parameter in constraining the dark energy equation of state and it is used as input in many cosmological simulations (Freedman & Madore 2010; Planck Collaboration XVI 2014). Therefore, precise estimation of H_{0} is of utmost importance in cosmology.
Estimates of H_{0} available in the literature cover a wide range of uncertainties from ~2% to ~10% and the value ranges between 60 and 75 km s^{1} Mpc^{1}. The most reliable measurements of H_{0} known to date include

the Hubble Space Telescope (HST) Key Project(72 ± 8 km s^{1} Mpc^{1}; Freedman et al. 2001),

the HST Program for the Luminosity Calibration of Type Ia Supernovae by Means of Cepheids (62.3 ± 5.2 km s^{1} Mpc^{1}; Sandage et al. 2006),

Wilkinson Microwave Anisotropy Probe (WMAP) (70.0 ± 2.2 km s^{1} Mpc^{1}; Hinshaw et al. 2013),

Supernovae and H_{0} for the Equation of State (SH0ES) Program (73.8 ± 2.4 km s^{1} Mpc^{1}; Riess et al. 2011),

Carnegie Hubble Program (CHP) (74.3 ± 2.6 km s^{1} Mpc^{1}; Freedman et al. 2012),

the Megamaser Cosmology Project (MCP) (68.9 ± 7.1 km s^{1} Mpc^{1}; Reid et al. 2013; Braatz et al. 2013),

Planck measurements of the cosmic microwave background (CMB) anisotropies (67.3 ± 1.2 km s^{1} Mpc^{1}; Planck Collaboration XVI 2014), and

Strong lensing time delays (75.2 km s^{1} Mpc^{1}; Suyu et al. 2013).
It is worth noting here that the small uncertainties in H_{0} measurements resulting from WMAP and Planck crucially depend on the assumption of a spatially flat universe.
Although the values of H_{0} obtained from different methods are consistent with each other within 2σ given the current level of precision, all of the above methods of determination of H_{0} suffer from systematic uncertainties. Therefore, as the measurements increase in precision, multiple approaches based on different physical principles need to be pursued so as to be able to identify unknown systematic errors present in any given approach.
The phenomenon of strong gravitational lensing offers an elegant method to measure H_{0}. For gravitationally lensed sources that show variations in flux with time, such as quasars, it is possible to measure the time delay between the various images of the background source. The time delay, which is a result of the travel times for photons being different along the light paths corresponding to the lensed images, has two origins: (i) the geometric difference between the light paths and (ii) gravitational delay due to the dilation of time as photons pass in the vicinity of the lensing mass. Time delays, therefore depend on the cosmology, through the distances between the objects involved, and on the radial mass profile of the lensing galaxies. This was shown theoretically five decades ago by Refsdal (1964) long before the discovery of the first gravitational lens Q0957+561 by Walsh et al. (1979).
Estimation of H_{0} through gravitational lens time delays, although it has its own degeneracies, is based on the wellunderstood physics of General Relativity, and compared to distance ladder methods, is free from various calibration issues. In addition to measuring H_{0}, measurement of time delays between the light curves of a lensed quasar can be used to study the microlensing variations present in the light curves, and to study the structure of the quasar (Hainline et al. 2013; Mosquera et al. 2013). However, these time delay measurements of H_{0} are extremely challenging because of the need of an intensive monitoring program that offers high cadence and goodquality photometric data over a long period of time. This type of program would then be able to cope with the presence of uncorrelated variations present in the lensed quasar light curves, which can interestingly arise due to microlensing by stars in the lensing galaxy (Chang & Refsdal 1979) or for mundane reasons, such as the presence of additive flux shifts in the photometry (Tewes et al. 2013a). Moreover, the estimation of H_{0} from such highquality data is hampered by the uncertainty on lens models. Recently, using time delay measurements from highquality optical and radio light curves, deep and highresolution imaging observations of the lensing galaxies and lensed AGN host galaxy, and the measurement of stellar velocity dispersion of the lens galaxy to perform detailed modeling, Suyu et al. (2013) report a H_{0} of 75.2 km s^{1} Mpc^{1} through the study of two gravitational lenses namely RX J1131−1231 and CLASS B1608+656.
Another approach is to perform simple modeling of a relatively large sample of gravitational lenses with moderateprecision time delay measurements. In this way, it should be possible to obtain a precise determination of the global value of H_{0}, even if the H_{0} measurements from individual lenses have large uncertainties. In addition, when inferring H_{0} from a relatively large sample of lenses, lineofsight effects that bias the H_{0} measurements from individual lenses (see Suyu et al. 2013, Sect. 2) should tend to average out, although a residual systematic error must still remain (Hilbert et al. 2007; Fassnacht et al. 2011). A pixelized method of lens modeling is available in the literature and is also implemented in the publicly available code PixeLens (Saha & Williams 2004). Using this code, Saha et al. (2006) have found km s^{1} Mpc^{1} for a sample of ten time delay lenses. Performing a similar analysis on an extended sample of 18 lenses Paraficz & Hjorth (2010) obtained km s^{1} Mpc^{1}. Here, we present an estimate of H_{0} using the pixellated modeling approach on a sample of carefully selected lensed quasars. So far, time delays have been reported for 24 gravitationally lensed quasars among the hundreds of such strongly lensed quasars known. However, the quality of the light curves and the techniques used to infer these time delays vary between systems. In this work, we apply the differencesmoothing technique, introduced in Rathna Kumar et al. (2013), to the publicly available light curves of the 24 systems in a homogeneous manner, first to crosscheck the previously measured time delays and then to select a subsample of suitable lens systems to determine H_{0}.
The paper is organized as follows. Section 2 describes the technique used for time delay determination and introduces a recipe for creating realistic simulated light curves with known time delays; the simulated light curves are used in this work to estimate the uncertainty of each measured delay. In Sect. 3, the application of the curveshifting procedure to the 24 systems is described. In Sect. 4, we infer H_{0} from the lensmodeling of those systems that have at least one reliably measured time delay, known lens redshift, accurate astrometric data, and sufficiently simple mass distribution. We conclude in Sect. 5.
2. Time delay determination
In this section, we briefly describe the previously reported differencesmoothing technique, which contains one modification to the original version (see Rathna Kumar et al. 2013 for details). We then introduce a recipe for simulating realistic light curves having known time delays in a plausible range around the measured delay in order to estimate its uncertainty. We also present an approach for tuning the free parameters of the differencesmoothing technique for a given dataset.
2.1. Differencesmoothing technique
A_{i} and B_{i} are the observed magnitudes constituting light curves A and B sampled at epochs t_{i} (i = 1,2,3,...,N). Light curve A is selected as the reference. We shift light curve B in time with respect to light curve A by an amount τ. This shifted version B′ of B is given by
We note here that we do not apply any flux shift to light curve B as in Rathna Kumar et al. (2013), since we have found that doing so considerably increases the computational time without significantly changing the results.
For any given estimate of the time delay τ, we form a difference light curve having points d_{i} at epochs t_{i}, (3)where the weights w_{ij} are given by (4)The parameter δ is the decorrelation length and σ_{Bj} denotes the photometric error of the magnitude B_{j}. We calculate the uncertainty of each d_{i} as (5)where w_{ij} are given by Eq. (4).
We now smooth the difference curve d_{i} using a Gaussian kernel to obtain a model f_{i} for the differential extrinsic variability(6)where the weights ν_{ij} are given by (7)The smoothing time scale s is another free parameter of this method. The uncertainty of each f_{i} is computed as (8)We optimize the time delay estimate τ to minimize the residuals between the difference curve d_{i} and the much smoother f_{i}. To quantify the mismatch between d_{i} and f_{i}, we define a normalized χ^{2}, (9)and minimize this using a global optimization.
In the above description, since light curves A and B are not interchangeable, we systematically perform all computations for both permutations of A and B, and minimize the average of the two resulting values of .
2.2. Simulation of light curves
In Rathna Kumar et al. (2013), in order to estimate the uncertainty of the time delay measured using the differencesmoothing technique, we made use of realistic simulated light curves, which were created following the procedure introduced in Tewes et al. (2013a). In this work, we introduce an independent recipe for creating simulated light curves.
We infer the underlying variation A(t) of the light curve A at the epoch t_{i} based on the magnitudes A_{j} for all the epochs as (10)where the value of m is set to equal the mean sampling of the light curves calculated after excluding the large gaps following a 3σ rejection criterion. For those points having the nearest neighboring points on both sides separated by a value less than or equal to m, we compute the values of (A_{i}−A(t_{i})) /σ_{Ai}, the standard deviation of which is multiplied to the error bars σ_{Ai} to obtain the rescaled error bars . We note here that the rescaling is applied for all the epochs and not just the epochs of points used in computing the rescaling factor. Similarly for the B light curves the rescaled error bars are obtained. This rescaling inferred from the local scatter properties of the light curves is done because the magnitudes of the original error bars may suffer from systematic underestimation or overestimation.
We merge light curves A and B by shifting the B light curve by the time delay found (Δt) and subtracting the differential extrinsic variability f_{i} corresponding to the delay from the A light curve. This merged light curve M_{i}, whose errors we denote σ_{Mi}, consists of the magnitudes A_{i}−f_{i} at times t_{i} and having errors and the magnitudes B_{i} at times t_{i} + Δt and having errors . We now model the quasar brightness variation M(t) as (11)We then model the quasar brightness variation using only the A points in M_{i} as (12)and only the B points in M_{i} as (13)The residual extrinsic variations present in the A and B light curves can now be calculated as (14)and (15)We can now simulate light curves and having a time delay of Δt + dt between them by sampling M(t) at appropriate epochs and adding terms for extrinsic variations and noise, (16)and (17)where N^{∗}(0,1) is a random variate drawn from a normal distribution having mean 0 and variance 1. These simulated light curves are then assigned the times t_{i} and the error bars σ_{Ai} and σ_{Bi} for the A and B light curves, respectively. Including the terms f_{Ai} and f_{Bi} in the calculation of and , respectively, ensures that our simulated light curves contain extrinsic variability on all time scales, just as in the real light curves.
Here again in the above description, since light curves A and B are not interchangeable, we systematically perform all computations for both permutations of A and B, and average the corresponding values of and , before adding the noise terms.
2.3. Choice of free parameters
The value chosen for the decorrelation length δ needs to be equivalent to the temporal sampling of the light curves. In this work, we set δ equal to m, the mean sampling of the light curves calculated after excluding the large gaps following a 3σ rejection criterion.
The value chosen for the smoothing time scale s needs to be significantly larger than δ. In this work, its value is optimized such that the larger of the maximum absolute values of and , which quantify the residual extrinsic variations in units of photometric noise for the A and B light curves respectively, is equal to 2. This choice ensures that the value of s is small enough to adequately model the extrinsic variations, so that the extreme values of residual extrinsic variations are not significantly larger than the noise in the data.
Again as in the above description, because light curves A and B are not interchangeable, we systematically perform all the computations for both permutations of A and B, and average the corresponding maximum absolute values.
Fig. 1 Light curves from the Strong Lens Time Delay Challenge file “tdc1_rung3_quad_pair9A.txt”. Light curve A is shown in red and light curve B in blue. 

Open with DEXTER 
Fig. 2 Light curves A and B from Fig. 1 have been merged, with light curve A as reference, after shifting light curve B by the measured time delay of Δt = −20.5 days and subtracting the differential extrinsic variability from A. M_{A}(t) sampled at the epochs t_{i} and M_{B}(t) sampled at the epochs t_{i} + Δt are connected by red and blue lines, respectively. M(t) sampled at the epochs t_{i} and t_{i} + Δt are connected by black lines. The optimum free parameters for this pair of light curves were found to be δ = 3.1 days and s = 139.0 days. The magnitudes at those epochs corresponding to maximum absolute values of and have been circled. The negative value of time delay implies that light curve A leads light curve B. The magnitudes are shown without error bars for convenience of display. 

Open with DEXTER 
2.4. Estimation of uncertainty
We create 200 simulated light curves having a true delay of Δt between them. The differencesmoothing technique is applied on each of them to obtain 200 delay values. The standard deviation of the 200 delay values gives us the random error, and the systematic error is obtained by the difference between the mean of the 200 delay values and the true delay. The total error Δτ_{0} is obtained by adding the random error and the systematic error in quadrature.
However, as noted by Tewes et al. (2013a), it is important to simulate light curves that have not only the time delay Δt found, but also other time delays in a plausible range around Δt, so as to obtain a reliable estimate of the uncertainty (see also Sect. 3.2 in Rathna Kumar et al. 2013). To this end, we also simulate 200 light curves for each true delay that differs from Δt by ± Δτ_{0}, ± (Δτ_{0} + Δτ_{1}),..., ± (Δτ_{0} + Δτ_{1} + ... + Δτ_{n−1}), in each step updating the total error Δτ_{n} by adding the maximum obtained value of the random error and the maximum obtained absolute value of the systematic error in quadrature. The value of n is chosen to be the smallest integer for which (18)This ensures that we have simulated light curves over a range of delay values that is at least as wide as or wider than the 95.4% confidence interval implied by the stated final error Δτ_{n}.
Fig. 3 Error analysis of the time delay measurement based on delay estimations on simulated light curves that mimic the light curves displayed in Fig. 1. The horizontal axis corresponds to the value of the true time delay used in these simulated light curves. The gray colored rods and 1σ error bars show the systematic biases and random errors, respectively. Our measured time delay of Δt = −20.5 ± 1.0 days is discrepant with the true time delay of 22.75 days listed in the TDC1 truth files at the level of 2.25σ. The difference in sign of time delay is simply a matter of convention. 

Open with DEXTER 
2.5. Testing the robustness of the procedure
In order to test the robustness of our procedure for estimating the time delay and its uncertainty, we made use of synthetic light curves from the TDC1 stage of the Strong Lens Time Delay Challenge^{1} (Liao et al. 2015), which are arranged in five rungs having different sampling properties (see Liao et al. 2015, Table 1). We applied our procedure on a sample of 250 light curves, 50 from each rung, selected such that we were able to reliably measure time delays from them. Comparing our results with the truth files, we found that all the measured delays agreed with the true delays to within twice the estimated uncertainties, except in one case. For the exceptional case, the discrepancy between the measured delay and the true delay was found to be 2.25σ. This is still a reasonably good level of agreement, thus demonstrating the robustness of our procedure. We note here that this property of robustness also depends on the careful choice of free parameters as presented here. For instance, setting δ equal to the mean sampling of the light curves computed without excluding the large gaps was found to lead to biased time delay measurements, which was especially noticeable for light curves having shorter seasons and larger cadence. We show some plots for the pair of TDC1 light curves corresponding to the exceptional case mentioned above in Figs. 1–3.
3. Time delays of 24 gravitationally lensed quasars
Time delays have been reported for 24 gravitationally lensed quasars. However, the quality of the data and the curveshifting procedure followed differs from system to system. In this section, we present a homogeneous analysis of their publicly available light curves following the procedure described in the previous section, with the aim of identifying those systems that have reliable time delay measurements. In the case of systems with more than two images, we measured the time delays between all pairs of light curves. The results are summarized in Table 1. All quoted uncertainties are 1σ error bars, unless stated otherwise. Additional information on some systems listed in Table 1 and discussion on the possible reasons for our inability to reliably measure some of the time delays follow. For all the other systems, our time delay measurements agree with the previously reported values to within 2σ.

Q0142−100 (UM673): we were unable to make a reliable timedelay measurement using the light curves presented in Koptelovaet al. (2012). This is not surprisinggiven that the light curves are characterized by large seasonal gapsand there are no clear variability features that could be matchedbetween the A and B light curves.

JVAS B0218+357: from 8 GHz and 15 GHz VLA observations reported by Cohen et al. (2000), we measured time delays of 10.4 ± 1.0 days and 11.4 ± 1.5 days, respectively. Taking the weighted average of the two results, we find the time delay to be 10.7 ± 0.8 days. We note here that Biggs et al. (1999) monitored this system using VLA during the same period as Cohen et al. (2000) at the same two frequencies and report a time delay of 10.5 ± 0.4 days (95% CI).

HE 0435−1223: we made use of the light curves presented in Courbin et al. (2011) spanning seven seasons using data from Euler, Mercator, Maidanak, and SMARTS and the light curves presented in Blackburne et al. (2014) spanning eight seasons using data from SMARTS. The SMARTS data used by Courbin et al. (2011) is the same as the first two seasons of data presented in Blackburne et al. (2014). Hence we excluded the SMARTS data points from the light curves of Courbin et al. (2011) to make it independent of the light curves of Blackburne et al. (2014). Owing to the differences in the approaches followed by these two teams of authors to derive photometry and also the photometric uncertainties, we avoided merging the two datasets. Our time delay measurements listed in Table 1 are the weighted averages of the time delays measured from the two independent sets of light curves. The reported time delay values in Table 1 are from Courbin et al. (2011). The bestfit time delay values reported without uncertainties by Blackburne et al. (2014) are consistent with the values of Courbin et al. (2011) to within 1σ. In Table 2, we present our measurements of the time delays of HE 0435−1223 from the two independent sets of light curves and the resulting weighted averages. For each pair of quasar images, we see that the time delays measured from the two datasets agree to within 2σ.

SBS 0909+532: for our analysis, we used only the rband data points obtained using the Liverpool Robotic Telescope between 2005 January and 2007 January presented in Goicoechea et al. (2008) and Hainline et al. (2013), based on homogeneity and sampling considerations.

RX J0911.4+0551: we used the light curves presented in Hjorth et al. (2002), which were made publicly available by Paraficz et al. (2006).

FBQ 0951+2635: we used the light curves presented in Jakobsson et al. (2005), which were made publicly available by Paraficz et al. (2006).

Q0957+561: from the rband and gband light curves presented in Shalyapin et al. (2012), we measured time delays of 420.6 ± 1.8 days and 419.2 ± 2.2 days, respectively. Taking the weighted average of the two results, we find the time delay to be 420.0 ± 1.4 days. The reported delay listed is the weighted average of the two delays found by Shalyapin et al. (2012).

RX J1131−1231: Tewes et al. (2013b) measured time delays between all pairs of light curves using three different numerical techniques. The time delay value listed in the table for each pair of light curves is for the technique that resulted in the smallest uncertainty.

H1413+117: the light curves presented in Goicoechea & Shalyapin (2010) span less than one season and display poor variability. Hence our time delay measurements for the pairs AB, AC, and AD although in good agreement with the reported values, are of low precision and we could not reliably measure time delays for the pairs BC and BD.

CLASS B1600+434: from both the optical light curves presented in Burud et al. (2000; and made publicly available by Paraficz et al. 2006) and the radio light curves presented in Koopmans et al. (2000), we were unable to make a reliable time delay measurement. Although the optical light curves show good variability, they suffer from poor sampling and thus exclude the possibility of convincingly matching the variability features between light curves A and B. The radio light curves spanning one season is well sampled; however, light curve A displays short time scale fluctuations that are not seen in light curve B, thus making it difficult to measure the time delay unambiguously.

HE 2149−2745: we used the light curves presented in Burud et al. (2002a), which were made publicly available by Paraficz et al. (2006).
Summary of time delay measurements.
Our measurements of the time delays of HE 0435−1223 from two independent datasets.
4. H_{0} from pixellated modeling of ten gravitational lenses
Of the 24 systems analyzed in the last section, 14 of them had light curves of sufficiently good quality to enable the measurement of at least one time delay between the images, adjacent to each other in terms of arrivaltime order, to a precision of better than 20% (which corresponds to a 5σ detection of time delay). The ten systems which did not satisfy this criterion are Q0142−100 (UM673), FBQ 0951+2635, PG 1115+080, H1413+117, JVAS B1422+231, CLASS B1600+434, SDSS J1650+4251, PKS 1830−211, HE 2149−2745, and HS 2209+1914.
Of the 14 remaining systems, we did not model the mass distribution for four of them for the following reasons. SDSS J1001+5027 and SDSS J1206+4332 do not have accurate astrometric data measured from Hubble Space Telescope (HST) images or groundbased imaging with adaptive optics. Although the astrometry of JVAS B0218+357, which has a small image separation of 0.33″, has been measured from HST images by Sluse et al. (2012), the authors warn about possibly large systematic errors in the published astrometry. SDSS J1029+2623 is a threeimage cluster lens with highly complex mass distribution (see Oguri et al. 2013) and hence not amenable to lensmodeling following the simplistic approach described below.
To perform massmodeling of the remaining ten systems – HE 0435−1223, SBS 0909+532, RX J0911.4+0551, Q0957+561, SDSS J1004+4112, HE 1104−1805, RX J1131− 1231, SBS 1520+530, CLASS B1608+656 and WFI J2033− 4723 – to infer H_{0}, we used the publicly available PixeLens^{2} code (Saha & Williams 2004), which builds an ensemble of pixellated mass maps compatible with the input data for a given system, which is comprised of the redshifts of the quasar and the lensing galaxy, the arrivaltime order of the images, their astrometry relative to the center of the main lensing galaxy, and the known time delays between the images adjacent to each other in terms of arrivaltime order. In case of quadruple lenses in which only some of the time delays are known, it is still possible to guess the arrivaltime order of the images by following certain simple rules (see Saha & Williams 2003).
We model all lenses, except SDSS J1004+4112, such that their mass profiles have inversion symmetry about the lens center, including any companion galaxy to the main lensing galaxy as a point mass. The lensing cluster in SDSS J1004+4112 consists of several galaxies besides the main lensing galaxy (see Inada et al. 2005) and hence was modeled without assuming inversion symmetry about the lens center.
PixeLens builds models such that their projected density profiles are steeper than  θ  ^{− γmin}, where  θ  is the distance from the center of the lens in angular units, the default value of γ_{min} being 0.5. This is based on the observation that the total density distribution in the central regions of elliptical galaxies is close to isothermal (i.e., r^{2}) and also the observation that the total density in the center of our Galaxy scales is r^{1.75} (see Saha & Williams 2004, Sect. 2.2 and references therein). The profiles r^{2} and r^{1.75} correspond to projected density profiles of  θ  ^{1} and  θ  ^{0.75}, respectively, in the special case of spherical symmetry. In this work, we relax the restriction of γ_{min} = 0.5 and set γ_{min} = 0 for those lenses in our sample in which the largest angular separation between the images is greater than 3″. The lenses in our sample that satisfy this criterion are RX J0911.4+0551, Q0957+561, SDSS J1004+4112, HE 1104−1805, and RX J1131−1231. A large image separation implies that there is significant lensing action from the cluster of which the main lensing galaxy is part, in which case the projected density profile can be shallower than  θ  ^{0.5}.
For each system, we build an ensemble of 100 models, corresponding to 100 values of H_{0}. The mean of the 100 values gives the best estimate of H_{0}, the uncertainty of which is the standard deviation of the 100 values. This uncertainty includes only the uncertainty in the mass model. PixeLens assumes that the uncertainty in the input priors to be negligibly small, which is a reasonable assumption for the redshifts, if they are spectroscopically measured, and astrometry, if measured from HST or groundbased adaptive optics imaging. However, the measured time delays have finite uncertainties, which need to be propagated into the uncertainty of the estimated H_{0}. We do this by remodeling each system after perturbing the time delay by its 1σ uncertainty and noticing the deviation of the resulting value of H_{0} from the original value. For highprecision time delays, the deviation in H_{0} was found to be the same whether the delays were perturbed upward or downward. In general, the deviation in H_{0} was found to be slightly larger when the delays were perturbed downward than when they were perturbed upward. Hence in this work, to get a conservative estimate of the contribution of the time delay uncertainty to the uncertainty in H_{0}, we decrease the time delay by its 1σ uncertainty and find the resulting increase in H_{0}. This uncertainty in H_{0} resulting from the time delay uncertainty is added in quadrature to the uncertainty in H_{0} resulting from mass modeling to find the total uncertainty. In the case of quadruple lenses where more than one time delay is known, we perturb each delay individually while leaving the other delays unchanged to infer its uncertainty contribution. The uncertainty contribution from each independent time delay is then added in quadrature to the uncertainty in H_{0} resulting from the uncertainty in the mass model to find the total uncertainty.
Summary of input data to PixeLens and resulting H_{0} estimates.
In order to include the effects of external shear, an approximate direction of the shear axis needs to be specified and PixeLens will search for solutions within 45° of the specified direction. Since there is no simple rule to guess the direction of the external shear for a given system, for each system, we repeated the modeling specifying the approximate direction of the shear axis as 90°, 45°, 0°, and −45° (in this instance, specifying θ and θ+180° are equivalent). We thus obtain four estimates of H_{0} and their uncertainties. In each case, we propagate the uncertainty contributions from the known time delays to the uncertainty in H_{0}, as discussed previously. The final estimate of H_{0} and its uncertainty are found using maximum likelihood analysis, optimizing their values so as to maximize the joint posterior probability of these two parameters for the sample consisting of the four H_{0} values and their uncertainties (see Barnabè et al. 2011, Eq. (7)). In optimizing the value of the uncertainty, we choose the minimum limit to be the smallest of the four uncertainties. We note here that for the system HE 1104−1805, the choices of the approximate direction of the shear axis of 90° and −45° were found to lead to unphysical models involving negative values in the mass pixels. Hence for this system, the maximum likelihood analysis was carried out using only the two H_{0} values resulting for the approximate shear directions of 45° and 0°.
The input priors for each system and the resulting H_{0} estimates are summarized in Table 3. In Fig. 4 we plot the H_{0} estimates from the ten lenses, all of which are seen to agree with each other within their error bars. To combine the ten independent estimates into a best estimate of H_{0}, we again employ maximum likelihood analysis, as described above. However, in this case, in optimizing the value of the uncertainty of the best estimate of H_{0}, the minimum limit is chosen to be the uncertainty of the weighted average of the ten values. We infer a value of H_{0} of 68.1 ± 5.9 km s^{1} Mpc^{1} (1σ uncertainty, 8.7% precision) for a spatially flat universe having Ω_{m} = 0.3 and Ω_{Λ} = 0.7. The reason for employing maximum likelihood analysis in this case, rather than taking a simple weighted average is to detect the presence of any unmodeled uncertainties. However, as can be seen from Fig. 4, the H_{0} estimates from the individual systems all agree with each other within their error bars and hence the H_{0} value inferred above through maximum likelihood analysis is only marginally different from the weighted average. For the source and lens redshifts of the current sample, we find the H_{0} estimate to decrease by 7.1% for the Einsteinde Sitter universe (Ω_{m} = 1.0 and Ω_{Λ} = 0.0) and increase by 2.3% for an open universe having Ω_{m} = 0.3 and Ω_{Λ} = 0.0, thus illustrating the low level of dependence of the inferred value of H_{0} on the precise values of Ω_{m} and Ω_{Λ}. In Table 3, we also list the H_{0} estimates obtained without propagating the time delay uncertainties. We see that the dominant contribution to uncertainty in H_{0} results from the uncertainty in the mass model.
Fig. 4 H_{0} estimates and their 1σ uncertainties for the ten gravitational lenses – (1) HE 0435−1223; (2) SBS 0909+532; (3) RX J0911.4+0551; (4) Q0957+561; (5) SDSS J1004+4112; (6) HE 1104−1805; (7) RX J1131−1231; (8) SBS 1520+530; (9) CLASS B1608+656; and (10) WFI J2033−4723. The best estimate of H_{0} and its 1σ confidence interval, inferred through maximumlikelihood analysis, are represented by the horizontal line and the gray shaded region, respectively. 

Open with DEXTER 
5. Conclusion
We have presented a homogeneous curveshifting analysis of the light curves of 24 gravitationally lensed quasars for which time delays have been reported in the literature so far. Time delays were measured using the differencesmoothing technique and their uncertainties were estimated using realistic simulated light curves; a recipe for creating these light curves with known time delays in a plausible range around the measured delay was introduced in this work. We identified 14 systems to have light curves of sufficiently good quality to enable the measurement of at least one time delay between the images, adjacent to each other in terms of arrivaltime order, to a precision of better than 20% (including systematic errors). Of these 14 systems, we performed pixellated mass modeling using the publicly available PixeLens software for ten of them, which have known lens redshifts, accurate astrometric information, and sufficiently simple mass distributions, to infer the value of H_{0} to be 68.1 ± 5.9 km s^{1} Mpc^{1} (1σ uncertainty, 8.7% precision) for a spatially flat universe having Ω_{m} = 0.3 and Ω_{Λ} = 0.7. We note here that we have followed a relatively simple lens modeling approach to constrain H_{0} and our analysis does not account for biases resulting from lineofsight effects.
Our measurement closely matches a recent estimate of H_{0} = 69.0 ± 6 (stat.) ± 4 (syst.) km s^{1} Mpc^{1} found by Sereno & Paraficz (2014) using a method based on freeform modeling of 18 gravitational lens systems. Our value is also consistent with the recent measurements of H_{0} by Riess et al. (2011), Freedman et al. (2012) and Suyu et al. (2013); however, it has lower precision. Increasing the number of lenses with goodquality light curves, accurate astrometry, and known lens redshift from the current ten used in this study can bring down the uncertainty in H_{0}.
In the future such highprecision time delays will become available from projects such as COSMOGRAIL (Tewes et al. 2012) involving dedicated mediumsized telescopes. In addition, the next generation of cosmic surveys such as the Dark Energy Survey (DES), the Large Synoptic Survey Telescope (LSST; Ivezic et al. 2008), and the Euclid mission will detect a large sample of lenses, and time delays might be available for a large fraction of them, consequently enabling measurement of H_{0} to an accuracy better than 2%. Furthermore, detection of gravitational wave signals from short gammaray bursts associated with neutron star binary mergers in the coming decade could constrain H_{0} to better than 1% (Nissanke et al. 2013).
Acknowledgments
We thank Jim Lovell for providing us with the light curves of PKS 1830−211. We acknowledge useful discussions with Malte Tewes, G. Indu, Leon Koopmans, Prashanth Mohan and Matthias Bartelmann. We thank Frederic Courbin and Georges Meylan for carefully reading the manuscript and offering helpful comments. We thank the organizers of the Strong Lens Time Delay Challenge for enabling a blind test of our algorithm and subsequently providing the truth files, which helped us to refine our curveshifting procedure. We thank the anonymous referee for constructive reports that helped to improve the presentation of this work.
References
 Auger, M. W., Fassnacht, C. D., Wong, K. C., et al. 2008, ApJ, 673, 778 [NASA ADS] [CrossRef] [Google Scholar]
 Bade, N., Siebert, J., Lopez, S., Voges, W., & Reimers, D. 1997, A&A, 317, L13 [NASA ADS] [Google Scholar]
 Barnabè, M., Czoske, O., Koopmans, L. V. E., Treu, T., & Bolton, A. S. 2011, MNRAS, 415, 2215 [NASA ADS] [CrossRef] [Google Scholar]
 Biggs, A. D., Browne, I. W. A., Helbig, P., et al. 1999, MNRAS, 304, 349 [NASA ADS] [CrossRef] [Google Scholar]
 Blackburne, J. A., Kochanek, C. S., Chen, B., Dai, X., & Chartas, G. 2014, ApJ, 789, 125 [NASA ADS] [CrossRef] [Google Scholar]
 Braatz, J., Reid, M., Kuo, C.Y., et al. 2013, in IAU Symp. 289, ed. R. de Grijs, 255 [Google Scholar]
 Burud, I., Hjorth, J., Jaunsen, A. O., et al. 2000, ApJ, 544, 117 [NASA ADS] [CrossRef] [Google Scholar]
 Burud, I., Courbin, F., Magain, P., et al. 2002a, A&A, 383, 71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Burud, I., Hjorth, J., Courbin, F., et al. 2002b, A&A, 391, 481 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chang, K., & Refsdal, S. 1979, Nature, 282, 561 [NASA ADS] [CrossRef] [Google Scholar]
 Chavushyan, V. H., Vlasyuk, V. V., Stepanian, J. A., & Erastova, L. K. 1997, A&A, 318, L67 [NASA ADS] [Google Scholar]
 Cohen, A. S., Hewitt, J. N., Moore, C. B., & Haarsma, D. B. 2000, ApJ, 545, 578 [NASA ADS] [CrossRef] [Google Scholar]
 Courbin, F., Chantry, V., Revaz, Y., et al. 2011, A&A, 536, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eigenbrod, A., Courbin, F., Meylan, G., Vuissoz, C., & Magain, P. 2006, A&A, 451, 759 [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]
 Fadely, R., Keeton, C. R., Nakajima, R., & Bernstein, G. M. 2010, ApJ, 711, 246 [NASA ADS] [CrossRef] [Google Scholar]
 Fassnacht, C. D., Womble, D. S., Neugebauer, G., et al. 1996, ApJ, 460, L103 [NASA ADS] [CrossRef] [Google Scholar]
 Fassnacht, C. D., Pearson, T. J., Readhead, A. C. S., et al. 1999, ApJ, 527, 498 [NASA ADS] [CrossRef] [Google Scholar]
 Fassnacht, C. D., Xanthopoulos, E., Koopmans, L. V. E., & Rusin, D. 2002, ApJ, 581, 823 [NASA ADS] [CrossRef] [Google Scholar]
 Fassnacht, C. D., Koopmans, L. V. E., & Wong, K. C. 2011, MNRAS, 410, 2167 [NASA ADS] [CrossRef] [Google Scholar]
 Fohlmeister, J., Kochanek, C. S., Falco, E. E., et al. 2007, ApJ, 662, 62 [NASA ADS] [CrossRef] [Google Scholar]
 Fohlmeister, J., Kochanek, C. S., Falco, E. E., Morgan, C. W., & Wambsganss, J. 2008, ApJ, 676, 761 [NASA ADS] [CrossRef] [Google Scholar]
 Fohlmeister, J., Kochanek, C. S., Falco, E. E., et al. 2013, ApJ, 764, 186 [NASA ADS] [CrossRef] [Google Scholar]
 Freedman, W. L., & Madore, B. F. 2010, ARA&A, 48, 673 [NASA ADS] [CrossRef] [Google Scholar]
 Freedman, W. L., Madore, B. F., Gibson, B. K., et al. 2001, ApJ, 553, 47 [NASA ADS] [CrossRef] [Google Scholar]
 Freedman, W. L., Madore, B. F., Scowcroft, V., et al. 2012, ApJ, 758, 24 [NASA ADS] [CrossRef] [Google Scholar]
 Goicoechea, L. J., & Shalyapin, V. N. 2010, ApJ, 708, 995 [NASA ADS] [CrossRef] [Google Scholar]
 Goicoechea, L. J., Shalyapin, V. N., Koptelova, E., et al. 2008, New A, 13, 182 [NASA ADS] [CrossRef] [Google Scholar]
 Hainline, L. J., Morgan, C. W., MacLeod, C. L., et al. 2013, ApJ, 774, 69 [NASA ADS] [CrossRef] [Google Scholar]
 Hilbert, S., White, S. D. M., Hartlap, J., & Schneider, P. 2007, MNRAS, 382, 121 [NASA ADS] [CrossRef] [Google Scholar]
 Hinshaw, G., Larson, D., Komatsu, E., et al. 2013, ApJS, 208, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Hjorth, J., Burud, I., Jaunsen, A. O., et al. 2002, ApJ, 572, L11 [NASA ADS] [CrossRef] [Google Scholar]
 Inada, N., Oguri, M., Pindor, B., et al. 2003, Nature, 426, 810 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Inada, N., Oguri, M., Keeton, C. R., et al. 2005, PASJ, 57, L7 [NASA ADS] [Google Scholar]
 Ivezic, Z., Tyson, J. A., Acosta, E., et al. 2008 [arXiv:0805.2366] [Google Scholar]
 Jakobsson, P., Hjorth, J., Burud, I., et al. 2005, A&A, 431, 103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kneib, J.P., Cohen, J. G., & Hjorth, J. 2000, ApJ, 544, L35 [NASA ADS] [CrossRef] [Google Scholar]
 Kochanek, C. S., Falco, E. E., Schild, R., et al. 1997, ApJ, 479, 678 [NASA ADS] [CrossRef] [Google Scholar]
 Koopmans, L. V. E., de Bruyn, A. G., Xanthopoulos, E., & Fassnacht, C. D. 2000, A&A, 356, 391 [NASA ADS] [Google Scholar]
 Koptelova, E., Chen, W. P., Chiueh, T., et al. 2012, A&A, 544, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Liao, K., Treu, T., Marshall, P., et al. 2015, ApJ, 800, 11 [NASA ADS] [CrossRef] [Google Scholar]
 Lidman, C., Courbin, F., Kneib, J.P., et al. 2000, A&A, 364, L62 [NASA ADS] [Google Scholar]
 Lovell, J. E. J., Jauncey, D. L., Reynolds, J. E., et al. 1998, ApJ, 508, L51 [NASA ADS] [CrossRef] [Google Scholar]
 Lubin, L. M., Fassnacht, C. D., Readhead, A. C. S., Blandford, R. D., & Kundić, T. 2000, AJ, 119, 451 [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, N. D., Kochanek, C. S., Pevunova, O., & Schechter, P. L. 2005, AJ, 129, 2531 [NASA ADS] [CrossRef] [Google Scholar]
 Mosquera, A. M., Kochanek, C. S., Chen, B., et al. 2013, ApJ, 769, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Myers, S. T., Fassnacht, C. D., Djorgovski, S. G., et al. 1995, ApJ, 447, L5 [NASA ADS] [CrossRef] [Google Scholar]
 Nissanke, S., Holz, D. E., Dalal, N., et al. 2013 [arXiv:1307.2638] [Google Scholar]
 Oguri, M., Inada, N., Keeton, C. R., et al. 2004, ApJ, 605, 78 [NASA ADS] [CrossRef] [Google Scholar]
 Oguri, M., Schrabback, T., Jullo, E., et al. 2013, MNRAS, 429, 482 [NASA ADS] [CrossRef] [Google Scholar]
 Paraficz, D., & Hjorth, J. 2010, ApJ, 712, 1378 [NASA ADS] [CrossRef] [Google Scholar]
 Paraficz, D., Hjorth, J., Burud, I., Jakobsson, P., & Elíasdóttir, Á. 2006, A&A, 455, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Patnaik, A. R., & Narasimha, D. 2001, MNRAS, 326, 1403 [NASA ADS] [CrossRef] [Google Scholar]
 Planck Collaboration XVI. 2014, A&A, 571, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Poindexter, S., Morgan, N., Kochanek, C. S., & Falco, E. E. 2007, ApJ, 660, 146 [NASA ADS] [CrossRef] [Google Scholar]
 Rathna Kumar, S., Tewes, M., Stalin, C. S., et al. 2013, A&A, 557, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Refsdal, S. 1964, MNRAS, 128, 307 [NASA ADS] [CrossRef] [Google Scholar]
 Reid, M. J., Braatz, J. A., Condon, J. J., et al. 2013, ApJ, 767, 154 [NASA ADS] [CrossRef] [Google Scholar]
 Riess, A. G., Macri, L., Casertano, S., et al. 2011, ApJ, 730, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Saha, P., & Williams, L. L. R. 2003, AJ, 125, 2769 [NASA ADS] [CrossRef] [Google Scholar]
 Saha, P., & Williams, L. L. R. 2004, AJ, 127, 2604 [NASA ADS] [CrossRef] [Google Scholar]
 Saha, P., Coles, J., Macciò, A. V., & Williams, L. L. R. 2006, ApJ, 650, L17 [NASA ADS] [CrossRef] [Google Scholar]
 Sandage, A., Tammann, G. A., Saha, A., et al. 2006, ApJ, 653, 843 [NASA ADS] [CrossRef] [Google Scholar]
 Sereno, M., & Paraficz, D. 2014, MNRAS, 437, 600 [NASA ADS] [CrossRef] [Google Scholar]
 Shalyapin, V. N., Goicoechea, L. J., & GilMerino, R. 2012, A&A, 540, A132 [NASA ADS] [CrossRef] [EDP Sciences] [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., Chantry, V., Magain, P., Courbin, F., & Meylan, G. 2012, A&A, 538, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Smette, A., Robertson, J. G., Shaver, P. A., et al. 1995, A&AS, 113, 199 [NASA ADS] [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., et al. 2012, The Messenger, 150, 49 [NASA ADS] [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]
 Tsvetkova, V. S., Vakulik, V. G., Shulga, V. M., et al. 2010, MNRAS, 406, 2764 [NASA ADS] [CrossRef] [Google Scholar]
 Vuissoz, C., Courbin, F., Sluse, D., et al. 2007, A&A, 464, 845 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vuissoz, C., Courbin, F., Sluse, D., et al. 2008, A&A, 488, 481 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Walsh, D., Carswell, R. F., & Weymann, R. J. 1979, Nature, 279, 381 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Wisotzki, L., Schechter, P. L., Bradt, H. V., Heinmüller, J., & Reimers, D. 2002, A&A, 395, 17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Our measurements of the time delays of HE 0435−1223 from two independent datasets.
All Figures
Fig. 1 Light curves from the Strong Lens Time Delay Challenge file “tdc1_rung3_quad_pair9A.txt”. Light curve A is shown in red and light curve B in blue. 

Open with DEXTER  
In the text 
Fig. 2 Light curves A and B from Fig. 1 have been merged, with light curve A as reference, after shifting light curve B by the measured time delay of Δt = −20.5 days and subtracting the differential extrinsic variability from A. M_{A}(t) sampled at the epochs t_{i} and M_{B}(t) sampled at the epochs t_{i} + Δt are connected by red and blue lines, respectively. M(t) sampled at the epochs t_{i} and t_{i} + Δt are connected by black lines. The optimum free parameters for this pair of light curves were found to be δ = 3.1 days and s = 139.0 days. The magnitudes at those epochs corresponding to maximum absolute values of and have been circled. The negative value of time delay implies that light curve A leads light curve B. The magnitudes are shown without error bars for convenience of display. 

Open with DEXTER  
In the text 
Fig. 3 Error analysis of the time delay measurement based on delay estimations on simulated light curves that mimic the light curves displayed in Fig. 1. The horizontal axis corresponds to the value of the true time delay used in these simulated light curves. The gray colored rods and 1σ error bars show the systematic biases and random errors, respectively. Our measured time delay of Δt = −20.5 ± 1.0 days is discrepant with the true time delay of 22.75 days listed in the TDC1 truth files at the level of 2.25σ. The difference in sign of time delay is simply a matter of convention. 

Open with DEXTER  
In the text 
Fig. 4 H_{0} estimates and their 1σ uncertainties for the ten gravitational lenses – (1) HE 0435−1223; (2) SBS 0909+532; (3) RX J0911.4+0551; (4) Q0957+561; (5) SDSS J1004+4112; (6) HE 1104−1805; (7) RX J1131−1231; (8) SBS 1520+530; (9) CLASS B1608+656; and (10) WFI J2033−4723. The best estimate of H_{0} and its 1σ confidence interval, inferred through maximumlikelihood analysis, are represented by the horizontal line and the gray shaded region, respectively. 

Open with DEXTER  
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.