# A Comparison of Cosmological Models Using Time Delay Lenses

## Abstract

The use of time-delay gravitational lenses to examine the cosmological expansion introduces a new standard ruler with which to test theoretical models. The sample suitable for this kind of work now includes 12 lens systems, which have thus far been used solely for optimizing the parameters of CDM. In this paper, we broaden the base of support for this new, important cosmic probe by using these observations to carry out a one-on-one comparison between competing models. The currently available sample indicates a likelihood of that the Universe is the correct cosmology versus for the standard model. This possibly interesting result reinforces the need to greatly expand the sample of time-delay lenses, e.g., with the successful implementation of the Dark Energy Survey, the VST ATLAS survey, and the Large Synoptic Survey Telescope. In anticipation of a greatly expanded catalog of time-delay lenses identified with these surveys, we have produced synthetic samples to estimate how large they would have to be in order to rule out either model at a confidence level. We find that if the real cosmology is CDM, a sample of time-delay lenses would be sufficient to rule out at this level of accuracy, while time-delay lenses would be required to rule out CDM if the real Universe is instead . This difference in required sample size reflects the greater number of free parameters available to fit the data with CDM.

^{7}

## 1 Introduction

The idea of using gravitational lenses with time delays between the various images of a background quasar as a cosmological probe was first suggested by Refsdal (1964). Null geodesics originating with distant variable sources have different optical paths and pass through dissimilar gravitational potentials. Their deflection angles and time delays can thus be used to infer lens-system dependent angular-diameter distances, which can then be compared to theoretical predictions from general relativity to test the background cosmological expansion and offer the possibility of testing competing models (see, e.g., Petters et al. 2001; Schneider et al. 2006).

As of today, time delays have been observed from 21 lensed quasars, a relatively small
subset of the several hundred known strong-lens systems. But this is only the beginning.
In the near future, observational programmes, such as the COSmological MOnitoring of
GRAvItational Lenses (COSMOGRAIL; Eigenbrod et al. 2005) and perhaps also the
International Liquid Mirror Telescope (ILMT) project (Jean et al. 2001), which is now
in the final phases of construction in the Kumaun region of the Himalayas (Sagar et al.
2013), should increase this sample considerably. New strong gravitational lens systems
will also be discovered by cosmic structure surveys, including the Dark Energy
Survey^{8}

Several attempts have already been made to demonstrate the usefulness of these data for constraining the cosmological parameters in the standard model, CDM (see, e.g., Paraficz & Hjorth 2009, 2010; Balmès & Corasaniti 2013; and Suyu et al. 2013). In this paper, we broaden the base of support for this cosmic probe by demonstrating its usefulness in testing competing cosmological models. In particular, we show that the currently available sample of time-delay lensing systems favors the Universe with a likelihood of of being correct, versus for CDM. Though this result is still only marginal, it nonetheless calls for a significant increase in the sample of suitable lensing systems in order to carry out more sophisticated and higher precision measurements.

In § 2, we describe the key theoretical steps used in the application of time-delay lenses for cosmological testing, and we apply this procedure to the currently known sample of 12 systems in § 3. We discuss the results of our one-on-one comparison between CDM and in § 4. As we shall see, model selection tools favor the latter, but the distinction, given the relatively small number of lenses, is still not strong enough to completely rule out either model. We therefore estimate the sample size required from future surveys to reach likelihoods of and , for a confidence limit, and we present our conclusions in § 6.

## 2 Strong Lensing

The measurement of time delays in strong gravitational lenses is not straightforward, due in part to the uncertainty associated with the lens mass distribution and the possible presence of other perturbing masses along the line-of-sight. To this point, two principal methods have been employed to model the lens itself, which may be characterized as either simple-parametric (see, e.g., Oguri et al. 2002; Keeton et al. 2003) or grid-based parametric (see, e.g., Warren & Dye 2003; Bradac et al. 2008; Coles 2008; Suyu et al. 2013) approaches. The former uses simply-parametrized forms for the mass distribution of the deflector, while the latter uses as parameters a grid of pixels, to describe either the potential or the mass distribution of the deflector, and/or the source surface brightness distribution. Some have also used a hybrid approach, in which pixellated corrections were made to a simply parametrized mass model (Suyu et al. 2010; Vegetti et al. 2010).

The grid-based methods are regularized, often by imposing physical priors, otherwise they would just fit the noise. The simply parametrized methods can even be computationally more intensive, depending on the choice of the parameters. If an appropriate sub-sample of homogeneous systems can be identified—meaning a set of lenses whose properties provide evidence that the same lens model description may be used with the same level of statistical significance—the simply parametrized method can serve as an ideal first attempt at gauging whether the image-inversion effort is warranted with follow-up analysis. This is the method we will be using in this paper.

For a given image at angular position , with the source position at angle , the time delay, , is the combined effect of the difference in path length between the straight and deflected rays, and the gravitational time dilation for the ray passing through the effective gravitational potential of the lens, :

(1) |

(see, e.g., Blandford & Narayan 1986, and references cited therein), where and are the lens and source redshifts, respectively, and is the angular-diameter distance between redshifts and . If the lens geometry and the lens potential are known, the time delay measures the ratio

(2) |

also known as the time-delay distance, which depends on the cosmological model.

It has been known for over a decade that lens spiral and elliptical galaxies have a mass distribution that is well approximated by power-law density profiles (Witt et al. 2000; Rusin et al. 2003), for which the lens potential may be written

(3) |

in terms of the deflection scale and index . The single isothermal sphere (SIS) is the special case corresponding to , for which , where is the velocity dispersion of the lensing galaxy. Though some systems have shallow profiles with , measurements of galaxy density distributions suggest that is generally close to the isothermal value. Thus, in addition to the SIS model being convenient for its simplicity, it is actually also a surprisingly useful and accurate model for lens galaxies (Guimaraes & Sodré 2009; Koopmans et al. 2009). And for such systems with only two images at and , the time delay is given by the expression

(4) |

Treu et al. (2006) found that the ratio is very close to unity, where includes systematic errors in the rms deviation of the velocity dispersion and a softened isothermal sphere potential (see additional details below). Note that if the velocity dispersion of the lensing galaxy could also be observed, and assuming a ratio , two of the angular-diameter distances appearing in equation (2) could be replaced with the measured value of and the Einstein radius . An alternative approach would be to impose some prior (see, e.g., Oguri 2007), or to compute from a dynamical model (see, e.g., Treu & Koopmans 2002). As shown by Paraficz & Hjorth (2009), the quantity is more sensitive to the cosmological parameters than or separately, so this additional datum would improve the reliability with which this approach could distinguish between competing models. As of today, however, there are simply too few time-delay lenses with the corresponding measure of , so all of the analysis we carry out in this paper will be based solely on the use of equation (4). Even looking to the future, velocity dispersions are particularly difficult to measure for these systems precisely because they are crowded by quasars that make the time delay measurement possible.

In CDM, the angular-diameter distance depends on several parameters, including and the mass fractions , , and , defined in terms of the current matter (), radiation (), and dark energy () densities, and the critical density . Assuming zero spatial curvature, so that , the angular-diameter distance between redshifts and () is given by the expression

(5) |

where is the dark-energy equation of state. Thus, the essential free parameters in flat CDM are , and (since radiation is insignificant at gravitational lensing redshifts). In the Universe (Melia 2007; Melia & Shevchuk 2012), the angular-diameter distance depends only on . In this cosmology,

(6) |

## 3 Sample of Time-Delay (Two-image) Lensing Systems

In their careful Bayesian approach to constraining within the framework of CDM, Balmès & Corasaniti (2013) collected a sample of time-delay lenses for which Bayesian selection techniques can identify the lens mass model with the highest probability of describing the lens system. Rather than attempting to model individual lenses in all their complexity, the goal was to identify the model whose parameters significantly influence the time-delay, allowing them to average over individual mass model parameter uncertainties on a homogeneous mass sample. The first selection criterion in such an approach is therefore a restriction to two-image lenses, listed in Table 1, which seem to be more likely than other lens systems to be consistent with a simple power-law (or even SIS) profile. Paraficz & Hjorth (2010) followed the alternative method of using inversion techniques with each individual intensity image to map the mass distribution in each individual lens system, and produced a very useful comparison of their mass profiles, shown in Figure 1 of that paper. Indeed, the two-image lenses are significantly more symmetric than the rest. But though the object constituting the lens has been identified unambiguously in all the entries listed in Table 1, it is not yet clear whether these are part of a group or whether perturbators contribute along the line-of-sight. Thus, at this stage, an important caveat to our conclusions is that external shear may yet be contributing to some selection bias.

System | Refs. | ||||||||
---|---|---|---|---|---|---|---|---|---|

(arcsec) | (arcsec) | (days) | (Gpc) | (Gpc) | (Gpc) | ||||

B0218+357 | 0.685 | 0.944 | 5.268 | 5.361 | 1–3 | ||||

B1600+434 | 0.414 | 1.589 | 1.403 | 1.435 | 4,5 | ||||

FBQ0951+2635 | 0.26 | 1.246 | 0.917 | 0.956 | 6 | ||||

HE1104-1805 | 0.729 | 2.319 | 2.202 | 2.170 | 2,7,8 | ||||

HE2149-2745 | 0.603 | 2.033 | 1.909 | 1.908 | 6,9 | ||||

PKS1830-211 | 0.89 | 2.507 | 2.620 | 2.546 | 10,11 | ||||

Q0142-100 | 0.49 | 2.719 | 1.428 | 1.431 | 6,12 | ||||

Q0957+561 | 0.36 | 1.413 | 1.256 | 1.294 | 6, 13,14 | ||||

SBS 0909+532 | 0.83 | 1.377 | 4.085 | 4.072 | 6,15 | ||||

SBS 1520+530 | 0.717 | 1.855 | 2.435 | 2.419 | 6,16 | ||||

SDSS J1206+4332 | 0.748 | 1.789 | 2.632 | 2.612 | 17 | ||||

SDSS J1650+4251 | 0.577 | 1.547 | 2.079 | 2.098 | 6,18 |

^{9}

Note, however, that the two-image criterion may not be sufficient to guarantee a simple power-law lens model. Balmès and Corasaniti (2013) concluded from this sample that nine have Bayes factors favoring such a mass profile, though six of these are somewhat inconclusive. Thus, for a more balanced assessment in our analysis, we will consider two sample cuts, one with the full set of 12 two-image lenses listed in Table 1, and the second with just these nine: B1600+434, SBS 1520+530, SDSS J1650+4251, B0218+357, FBQ 0951+2635, HE 2149-2745, PKS 1830-211, Q0142-100, and SBS 0909+532.

For each model, we find the optimized fit by maximizing the joint likelihood function

(7) |

where ‘th’ stands for either CDM or , is the theoretical time-delay distance calculated from , and the model specific parameters , is the measured value, and is the dispersion of . The measured time-delay distance is

(8) |

so the propagated error in is

(9) |

The dispersion in the measured redshifts and (which enter through the angular distances in ) will be ignored here because a careful analysis of SDSS quasar spectra shows that (Hewett & Wild 2010), much smaller than the other errors appearing in equation (9).

However, we must include another source of error, in addition to , that we will call ; this takes into account at least several effects that apparently give rise to the observed scatter of individual lenses about the pure SIS profile. These include: systematic errors in the rms deviation of the velocity dispersion and a softened isothermal sphere potential, which tends to decrease the typical image separations. Moreover, Koopmans et al. (2009) showed that the mean mass density profile is slightly steeper than SIS and has significant scatter, not to mention that the line of sight contribution is generally non-zero on average (Suyu et al. 2010). According to Cao et al. (2012), may be as big as , depending on how many such factors actually contribute to this scatter. In our approach, we will adopt the additional free parameter to relate the dispersion to the measured effective lensing distance , according to

(10) |

We will add and in quadrature, and optimize the parameter for each fit individually though, as we shall see, the value of appears to be quite independent of the model itself. Thus, the total uncertainty in is given by the expression .

## 4 Results and Discussion

We have used the data shown in Table 1 to directly compare CDM with the Universe. The parameters in each model were individually optimized by maximizing the likelihood estimation, as described above. We will use two well established priors associated with the concordance CDM model: (i) dark energy is a cosmological constant, so ; and the spatial curvature constant is zero, i.e., the Universe is flat, so that . But to allow for added flexibility in the optimization of the model fit, we will employ three free parameters: , , and . We have restricted the fraction to the range . With the full sample of 12 time-delay lenses, CDM fits the data with a maximum likelihood for () km , () and (). The best fit with the Universe has only two free parameters: () km and (). The entries in column 7 of Table 1 are calculated from the observed angles and time delays. By comparison, columns 8 and 9 show the entries for and , respectively, corresponding to these best-fit parameters using all 12 lens systems.

To facilitate a direct comparison between CDM and , we show in Figure 1 the 12 observed values of compared with those predicted by these two theoretical models. The optimized values of and the maximum likelihood are quite similar for these two cases. However, these models formulate their observables (such as the angular diameter distances in Equations 5 and 6) differently, and do not have the same number of free parameters. Therefore a comparison of the likelihoods for either being closer to the ‘true’ model must be based on model selection tools.

Several information criteria commonly used in cosmology (see, e.g., Melia & Maier 2013, and references cited therein) include the Akaike Information Criterion, , where is the maximum likelihood, is the number of free parameters (Liddle 2007), the Kullback Information Criterion, (Cavanaugh et al. 2004), and the Bayes Information Criterion, , where is the number of data points (Schwarz et al. 1978). With characterizing model , the unnormalized confidence that this model is true is the Akaike weight . Model has likelihood

(11) |

of being the correct choice in this one-on-one comparison. Thus, the difference determines the extent to which is favoured over . For Kullback and Bayes, the likelihoods are defined analogously. For the two optimized fits discussed above, the magnitude of the difference , indicates that is to be preferred over . According to Equation (11), the likelihood of (i.e. ) being the correct choice is . For CDM (i.e. ), the corresponding value is . With the alternatives KIC and BIC, the magnitude of the differences and , indicates that is favored over CDM by a likelihood of versus .

We also carried out a one-on-one comparision using the reduced sample of only 9 two-image lens systems. In this case, the best CDM fit has a maximum likelihood for () km , () and (). For , the best fit corresponds to () km and (). Figure 2 is similar to Figure 1, except now for the reduced sample of 9 lenses. In this case, the magnitude of the differences , , and , indicates that is preferred over CDM with a likelihood of versus using AIC, versus using KIC, and versus using BIC.

## 5 Monte Carlo Simulations with a Synthetic Sample

Our results in this paper have shown that time-delay lenses can in fact be used to select one model over another in a one-on-one comparison. But though the likelihood of being closer to the correct cosmology than CDM is or higher, depending on the choice of information criterion, the outcome AIC AIC (and analogously for KIC and BIC) is judged ‘positive’ in the range , and ‘strong’ for . The constraints based on the currently known 12 lens systems should therefore be characterized as ‘positive,’ though not yet strong. In this section, we will estimate the sample size required to significantly strengthen the evidence in favor of or CDM, by conservatively seeking an outcome even beyond , i.e., we will see what is required to produce a likelihood versus , corresponding to .

We will consider two cases: one in which the background cosmology is assumed to be CDM, and a second in which it is , and we will attempt to estimate the number of time-delay lenses required in each case in order to rule out the alternative (incorrect) model at a confidence level. The synthetic time-delay lenses are each characterized by a set of parameters denoted as (, , , ), where , and are generated using the following procedure:

1. Since the 12 observed lens redshifts all fall in the range , and the source redshifts are , with a time delay (days), we assign uniformly between and , uniformly between and , and uniformly between and days.

2. We first infer from Equation 4 corresponding either to the Universe with km (§ 5.1), or to CDM with , and km (§ 5.2). We then assign a random deviation () to the value within the confidence interval, i.e., we put , where is a uniform random variable between and , and . This value of is taken from the current observed sample, which shows a median deviation . The same relative error is assigned to .

3. Since the observed is about of , we will also assign dispersion to the synthetic sample.

This sequence of steps is repeated for each lens system in the sample, which is enlarged until the likelihood criterion discussed above is reached. As with the real 12-lens sample, we optimize the model fits by maximizing the joint likelihood function in equation (7). We employ Markov-chain Monte Carlo techniques. In each Markov chain, we generate samples according to the likelihood function. Then we derive the coefficients and the cosmological parameters from a statistical analysis of the sample.

### 5.1 Assuming as the Background Cosmology

We have found that a sample of at least 1,000 time-delay lenses is required in order to rule out CDM at the confidence level. The optimized parameters corresponding to the best-fit CDM model for these simulated data are displayed in figure 3. To allow for the greatest flexibility in this fit, we relax the assumption of flatness, and allow to be a free parameter, along with . Figure 3 shows the 1-D probability distribution for each parameter (, , , ), and 2-D plots of the and confidence regions for two-parameter combinations. The best-fit values for CDM using the simulated sample with 1,000 lens systems in the Universe are , , km , and .

In figure 4, we show the corresponding 2-D contours in the plane for the Universe. The best-fit values for the simulated sample are km and .

Since the number of data points in the sample is now much greater than one, the most appropriate information criterion to use is the BIC. The logarithmic penalty in this model selection tool strongly suppresses overfitting if is large (the situation we have here, which is deep in the asymptotic regime). With , our analysis of the simulated sample shows that the BIC would favor the Universe over CDM by an overwhelming likelihood of versus only (i.e., the prescribed confidence limit).

### 5.2 Assuming CDM as the Background Cosmology

In this case, we assume that the background cosmology is CDM, and seek the minimum sample size to rule out at the confidence level. We have found that a minimum of 135 time-delay lenses are required to achieve this goal. To allow for the greatest flexibility in the CDM fit, here too we relax the assumption of flatness, and allow to be a free parameter, along with . In figure 5, we show the 1-D probability distribution for each parameter (, , , ), and 2-D plots of the and confidence regions for two-parameter combinations. The best-fit values for CDM using this simulated sample with 135 lens systems are , , km , and . Note that the simulated lenses give a good constraint on , but a weak one on ; only an upper limit of 0.90 can be set at the confidence level.

The corresponding 2-D contours in the plane for the Universe are shown in figure 6. The best-fit values for the simulated sample are km and . These are similar to those in the standard model, but not exactly the same, reaffirming the importance of reducing the data separately for each model being tested. With , our analysis of the simulated sample shows that in this case the BIC would favor CDM over by an overwhelming likelihood of versus only (i.e., the prescribed confidence limit).

## 6 Conclusions

The general agreement between theory and observation displayed in Figures 1 and 2 is promising, particularly since this work was based on the use of a single isothermal sphere profile for the lens mass distribution. It would be helpful to have additional information from which one may extract the lens structure from individual images. Up to now, these approaches have been used to optimize parameters in CDM itself, but not for an actual one-on-one comparison between competing cosmological models. This must be done because the results we have presented here already indicate a strong likelihood of being able to discriminate between models such as CDM and . Such comparisons have already been made using, e.g., cosmic chronometers (Melia & Maier 2013), gamma-ray bursts (Wei et al. 2013), and Type Ia SNe (Wei et al. 2014). The use of time-delay lenses introduces yet another standard ruler, with systematics different from those encountered elsewhere, thus providing an invaluable tool with which to cross-check the outcomes of these other important tests.

But though time-delay lens observations currently tend to favor over CDM, the known sample of such systems is still too small for us to completely rule out either model. We have therefore considered two synthetic samples with characteristics similar to those of the 12 known systems, one based on a CDM background cosmology, the other on . From the analysis of these simulated lenses, we have estimated that a sample of about 150 lenses would be needed to rule out at a confidence level if the real cosmology is in fact CDM, while a sample of at least 1,000 systems would be needed to similarly rule out CDM if the background cosmology were instead . The difference in required sample size results from CDM’s greater flexibility in fitting the data, since it has a larger number of free parameters. Such a level of accuracy may be achievable with the successful implementation of surveys, such as DES, VST ATLAS, and LSST.

### Footnotes

- affiliation: Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, China; jjwei@pmo.ac.cn, xfwu@pmo.ac.cn.
- affiliation: University of Chinese Academy of Sciences, Beijing 100049, China
- affiliation: Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, China; jjwei@pmo.ac.cn, xfwu@pmo.ac.cn.
- affiliation: Chinese Center for Antarctic Astronomy, Nanjing 210008, China.
- affiliation: Joint Center for Particle, Nuclear Physics and Cosmology, Nanjing University-Purple Mountain Observatory, Nanjing 210008, China.
- affiliation: Department of Physics, The Applied Math Program, and Department of Astronomy, The University of Arizona, AZ 85721, USA; fmelia@email.arizona.edu.
- slugcomment:
- One should take note of the fact, however, that the actual image quality in these surveys may be inferior to that expected, somewhat mitigating the possible yield of suitable lens systems for this work. For example, the DES was aiming to get median FWHM in the r, i, and z band images for its wide survey. At the end of the first year, the quality is close to this, though not yet meeting expectations. In addition, g and Y bands are observed in poorer seeing conditions so their quality is even lower (Bernstein 2014). This is an important consideration in any discussion concerning anticipated sample size, given that even SDSS has discovered only of the originally expected lens systems.
- footnotetext: (1) Carilli et al. (1993); (2) Lehár et al. (2000); (3) Wucknitz et al. (2004); (4) Jackson et al. (1995); (5) Dai & Kochanek (2005); (6) Kochanek et al. (2008); (7) Wisotzki et al. (1993); (8) Poindexter et al. (2007); (9) Burud et al. (2002); (10) Lovell et al. (1998); (11) Meylan et al. (2005); (12) Koptelova et al. (2012); (13) Falco et al. (1997); (14) Colley et al. (2003); (15) Dai & Kochanek (2009); (16) Auger et al. (2008); (17) Paraficz et al. (2009); (18) Vuissoz et al. (2007).

### References

- Auger, M. W., Fassnacht, C. D., Wong, K. C., Thompson, D., Matthews, K. & Soifer, B. T. 2008, ApJ, 673, 778
- Balmès, I. & Corasaniti, P. S. 2013, MNRAS, 431, 1528
- Banerji, M., Abdalla, F., Lahav, O. & Lin, H. 2008, MNRAS, 386, 1219
- Bernstein, G. 2014, private communication
- Blandford, R. & Narayan, R. 1986, ApJ, 310, 568
- Blandford, R. D., Oguri, M., Marshall, P., Baltz, E. A., Bradac, M., Fassnacht, C. D. & LSST Collaboration 2006, AAS Meeting 209, No. 86.12
- Bradac M., Allen, S. W., Treu, T., Ebeling, H., Massey, R., Morris, R. G., von der Linden, A. & Applegate, D. 2008, ApJ, 687, 959
- Buckley-Geer, E. J., Dark Energy Survey Collaboration 2014, AAS Meeting 223, No. 248.01
- Burud, I. et al. 2002, A&A, 383, 71
- Cao, S., Pan, Y., Biesiada, M., Godlowski, W. & Zhu, Z.-H. 2012, JCAP, issue 3, id. 16
- Carilli, C. L., Rupen, M. P. & Yanny, B. 1993, ApJ, 412, L59
- Cavanaugh, J. E. 2004, Aust. N. Z. J. Stat., 46, 257
- Chang, C., Jarvis, M., Jain, B., Kahn, S. M., Kirkby, D., Connolly, A., Krughoff, S., Peng, E.-H. & Peterson, J. R. 2014, MNRAS, 434, 2121
- Coles, J. 2008, ApJ, 679, 17
- Colley, W. N. et al. 2003, ApJ, 587, 71
- Dai, X. & Kochanek, C. S. 2005, ApJ, 625, 633
- Dai, X. & Kochanek, C. S. 2009, ApJ, 692, 667
- Eigenbrod, A., Courbin, F., Vuissoz, C., Meylan, G., Saha, P. & Dye, S. 2005, A&A, 436, 25
- Falco, E. E., Shapiro, I. I., Moustakas, Leonidas A. & Davis, M. 1997, ApJ, 484, 70
- Guimaraes, A.C.C. & Sodré, L. Jr. 2009, Highlights of Astronomy, 15, 75
- Hewett, P. C. & Wild, V. 2010, MNRAS, 405, 2302
- Jackson, N. et al. 1995, MNRAS, 274, L25
- Jean, C., Claeskens, J.-F. & Surdej, J. 2001, Gravitational Lensing: Recent Progress and Future Goals, ASP Conference Proceedings, Vol. 237. Ed. T. G. Brainerd and C. S. Kochanek (ASP: San Franciso)
- Keeton, C. R., Gaudi, B. S. & Petters, A. O. 2003, ApJ, 598, 138
- Kochanek, C. S., Falco, E. E., Impey, C., Lehár, J., McLeod, B. & Rix, H.-W. 2008, available at http://www.cfa.harvard.edu/castles/
- Koopmans, L.V.E. et al. 2009, ApJ, 703, L51
- Koposov, S. E., Irwin, M., Belokurov, V., Gonzales-Solares, E., Kupcu-Yoldas, A., Lewis, J., Metcalfe, N. and Schanks, T. 2014, MNRAS Letters, submitted (arXiv:1403.3409)
- Koptelova, E. et al. 2012, A&A, 544, A51
- Lehár, J. et al. 2000, ApJ, 536, 584
- Liddle, A. R. 2007, MNRAS, 377, L74
- Lovell, J.E.J., Jauncey, D. L., Reynolds, J. E., Wieringa, M. H., King, E. A., Tzioumis, A. K., McCulloch, P. M. & Edwards, P. G. 1998, ApJ, 508, L51
- Marshall, P. J., Sandford, C. P., Fassnacht, C. D., Meldgin, D. R., Oguri, M., Suyu, S. H., Auger, M. W., LSST Strong Lensing Science Collaboration 2011, AAS Meeting 217, No. 252.22
- Melia, F. 2007, MNRAS, 382, 1917
- Melia, F. 2013, ApJ, 764, 72
- Melia, F. & Maier, R. S. 2013, MNRAS, 432, 2669
- Melia, F., & Shevchuk, A. S. H. 2012, MNRAS, 419, 2579
- Meylan, G., Courbin, F., Lidman, C., Kneib, J.-P. & Tacconi-Garman, L. E. 2005, A&A, 438, L37
- Oguri, M., Taruya, A., Suto, Y. & Turner, E. L. 2002, ApJ, 568, 488
- Oguri, M. 2007, ApJ, 660, 1
- Oguri, M., & Marshall, P. J. 2010, MNRAS, 405, 2579
- Paraficz, D. & Hjorth, J. 2009, A&A Letters, 507, L49
- Paraficz, D. & Hjorth, J. 2010, ApJ, 712, 1378
- Petters, A. O., Levine, H., & Wambsganss, J. 2001, Singularity theory and gravitational lensing / Arlie O. Petters, Harold Levine, Joachim Wambsganss. Boston : Birkhäuser, c2001. (Progress in mathematical physics ; v. 21),
- Poindexter, S., Morgan, N., Kochanek, C. S. & Falco, E. E. 2007, ApJ, 660, 146
- Refsdal, S. 1964, MNRAS, 128, 307
- Rusin, D., Kochanek, C. S. & Keeton, C. R. 2003, ApJ, 595, 29
- Sagar, R., Kumar, B. & Omar, A. 2013, arXiv:1304.0235
- Schneider, P., Kochanek, C. & Wambsgnass, J. 2006, in Gravitational Lensing: Strong, Weak, and Micro: Saas-Fee Advanced Course, eds. Meylan, G., Jetzer, P. & North, P. (Springer: Berlin)
- Schneider, M. D. 2014, PRL, 112, id.061301
- Schwarz, G. 1978, Ann. Statist., 6, 461
- Suyu, S. H., Marshall, P. J., Auger, M. W., et al. 2010, ApJ, 711, 201
- Suyu, S. H. et al. 2013, ApJ, 766, 70
- Treu, T. & Koopmans, L.V.E. 2002, MNRAS Letters, 337, L6
- Treu, T., Koopmans, L. V., Bolton, A. S., Burles, S., & Moustakas, L. A. 2006, ApJ, 640, 662
- Tyson, T., Wittman, D., Hennawi, J. & Spergel, D. 2002, American Physical Society, April Meeting ID: APR02, abstract No. Y6.004
- Vegetti, S., Koopmans, L.V.E., Bolton, A., Treu, T. & Gavazzi, R. 2010, MNRAS, 408, 1969
- Vuissoz, C., Courbin, F., Sluse, D., Meylan, G., Ibrahimov, M., Asfandiyarov, I., Stoops, E., Eigenbrod, A., Le Guillou, L., van Winckel, H. & Magain, P. 2007, A&A, 464, 845
- Warren, S. J., & Dye, S. 2003, ApJ, 590, 673
- Wei, J.-J., Wu, X. & Melia, F. 2013, ApJ, 772, 43
- Wei, J.-J., Wu, X. & Melia, F. 2014, ApJ, submitted
- White, R. L. et al. 2000, ApJS, 126, 133
- Wucknitz, O., Biggs, A. D. & Browne, I.W.A. 2004, MNRAS, 349, 14
- Wisotzki, L., Koehler, T., Kayser, R. & Reimers, D. 1993, A&A, 278, L15
- Witt, H. J., Mao, S. & Keeton, C. R. 2000, ApJ, 544, 98
- York, T., Jackson, N., Browne, I.W.A., Wucknitz, O. & Skelton, J. E. 2005, MNRAS, 357, 124