Volume 568, August 2014
|Number of page(s)||4|
|Published online||14 August 2014|
Can one determine cosmological parameters from multi-plane strong lens systems?
Argelander-Institut für Astronomie, Universität Bonn, Auf dem Hügel 71, 53121 Bonn, Germany
Received: 23 June 2014
Accepted: 18 July 2014
Strong gravitational lensing of sources with different redshifts has been used to determine cosmological distance ratios, which in turn depend on the expansion history. Hence, such systems are viewed as potential tools for constraining cosmological parameters. Here we show that in lens systems with two distinct source redshifts, of which the nearest one contributes to the light deflection toward the more distant one, there exists an invariance transformation that leaves all strong-lensing observables unchanged (except for the product of time delay and Hubble constant), generalizing the well-known mass-sheet transformation in single-plane lens systems. The transformation preserves the relative location of mass and light. All time delays (from sources on both planes) scale with the same factor – time-delay ratios are therefore invariant under the mass-sheet transformation. Changing cosmological parameters, and thus distance ratios, is essentially equivalent to such a mass-sheet transformation. As an example, we discuss the double-source plane system SDSSJ0946+1006, which has recently been studied by Collett and Auger, and show that variations of cosmological parameters within reasonable ranges lead to only a weak mass-sheet transformation in both lens planes. Hence, the ability to extract cosmological information from such systems depends heavily on the ability to break the mass-sheet degeneracy.
Key words: gravitational lensing: strong / cosmological parameters
© ESO, 2014
Strong gravitational lensing by galaxies is a powerful tool for cosmological studies, in particular regarding precise mass estimates of the inner region of galaxies, the angular structure of the mass distribution (e.g., ellipticity and orientation), and mass substructure (see, e.g., Kochanek 2006; Treu 2010; Bartelmann 2010, and references therein). In his pioneering paper, Refsdal (1964) pointed out the possibility of determining cosmological parameters from lensing, specifically the Hubble constant, from measurements of time delays in lens systems. Whereas time delays have been determined in some 20 lens systems by now, the accuracy of the corresponding values of H0 is difficult to judge, because of the difficulty of reliably constraining the mass distribution of the lens: on the one hand, the number of observational constraints can be insufficient to constrain the mass distribution to sufficient accuracy. On the other hand, Falco et al. (1985) have found that there exists a transformation of the mass distribution of the lens that leaves all observables invariant, except for the product of time delay and H0. Thus, although very detailed studies of some lens systems have been conducted (see, e.g., Suyu et al. 2013), this mass-sheet transformation (MST) poses a fundamental limitation of the accuracy of the derived value for H0 from gravitational lensing alone. Although the degeneracy caused by the MST can be broken with additional observations, such as stellar dynamics in the lens or some additional information about the properties of the source (such as its luminosity), the current accuracy on these quantities leads to uncertainties of H0 larger than those from other methods, and the method may be biased. In practice, the degeneracy due to the MST is broken by assumptions about the mass distribution, e.g., that it follows a power law. A recent discussion of the impact of the MST on H0 determination can be found in Schneider & Sluse (2013).
A different approach to cosmology by strong-lens systems involves the relative lensing strength for two sources at two different redshifts. It is known that the classical MST is an exact transformation of the lens mass distribution only for sources at a single distance, and that it can be broken in principle when sources at several redshifts are employed (Bradač et al. 2004). However, it must be stressed that this degeneracy-breaking assumes that there is no second deflector along the line of sight to the sources. In reality, the lower-redshift sources have masses and consequently act as additional lenses for higher-redshift sources.
The recent discovery of a strong galaxy-scale lens with two extended multiple-image arc structures from two sources at vastly different redshifts (SDSSJ0946+1006; Gavazzi et al. 2008) opened up the possibility of studying a multi-plane lens system in great detail. Collett & Auger (2014; hereafter CA14) performed a detailed analysis of this system by constructing a model that involves both the main lens in the foreground of the two sources, and a smaller-mass deflector associated with the lower-redshift source. The relative lens strength β of the main lens on the two sources depends on distance ratios, which in turn depend on the redshifts of lenses and sources involved as well as on the distance-redshift relation in the Universe. Since the latter is sensitive to the density parameters and the equation of state parameter w of dark energy, CA14 were able to obtain constraints on w from this lens system.
In this Letter, we investigate whether this method can yield reliable results. Specifically, we show in Sect.2 that an analog of the MST also exists for the case of two lens and two source planes, as is the case for SDSSJ0946+1006. We then demonstrate in Sect.3 that a change of cosmological parameters that leads to a change of the expected value of β is equivalent to such a generalized MST. We then study the amplitude of this equivalent MST for a range of values of w, concluding that the transformation amplitude across a plausible range of w is indeed very small. Since the MST leads to a shape of the transformed density profile that is different from the original, we conclude that this method heavily relies on assumptions made for the shape of the mass profiles of the lenses.
In this section we first summarize the lensing equations for the multi-plane case (Sect.2.1), specify our notation, and recall the MST in the single-plane case (Sect.2.2) before we derive the new MST for two lens and two source planes.
We consider lenses and sources distributed along nearly the same line of sight at N different distances from us, characterized by their redshifts zi, or angular-diameter distances Di from us, 1 ≤ i ≤ N (we largely follow the notation of Schneider et al. 1992, where more details of the derivation are given). Perpendicular to the line of sight, we consider planes in which the sources at these distances are located, and onto which the mass distribution of the deflecting masses are projected (lens/source planes). We denote by Dij the angular-diameter distance of the jth plane as seen from the i-th plane, with 1 ≤ i<j ≤ N. The projected mass distribution in the i-th plane is characterized by its surface mass density Σi(ξ) and gives rise to a deflection angle , where (1)and ξi = Diθi is a transverse separation vector in the ith plane, with the corresponding unlensed angular position θi. The resulting propagation equations of a light ray follows solely from geometry and the definition of angular-diameter distances, (2)where we set θ ≡ θ1, scaled the deflection angles to the final source plane at i = N, i.e., (3)as coefficients of relative distance ratios. Accordingly, we define the dimensionless surface mass densities κi(θi) = 4πGDiDiN Σi(Diθi)/(c2DN), which satisfy ∇·αi = 2κi.
We briefly recall the MST for the single-lens plane, for which N = 2 (see Falco et al. 1985; Schneider & Sluse 2013, and references therein). Hence, we have a single lens plane at redshift z1 with deflection α(θ) ≡ α1(θ1) and dimensionless surface mass density κ(θ) ≡ κ1(θ1), and a single source plane at redshift z2.
If a mass distribution κ(θ) can explain all the observed lensing features, such as image positions, flux ratios, relative image shapes, and time-delay ratios, of a source with unlensed brightness profile Is(θ2), then the whole family of mass distributions and corresponding scaled deflection angles (4)explains the lensed features equally well for a source of brighness profile . Hence, the transformed brightness distribution in the source plane is rescaled by a factor λ. This rescaling of the source plane implies that the magnification of images is changed, μλ = μ/λ2, but that can only be observed if the luminosity or physical size of the source is known (standard candle or standard rod). In most cases, this information is unavailable, so that the magnification cannot be obtained from observations, whereas magnification (and thus flux) ratios are invariant under the MST. The time delay between pairs of images is changed under the MST, so that a time delay measurement can be employed to break the degeneracy caused by the MST, provided the Hubble constant (and the density parameters of our Universe) are assumed to be known.
We now study whether a similar invariance transformation exists for sources at two different redshifts z2 and z3, with lenses at redshift z1 and z2. In particular, we consider that the closer source at z2 is associated with mass that deflects light rays from the more distant source. By specializing the equations of Sect.2.1 to the case N = 3, we obtain the pair of relations (5)where we defined β ≡ β12 (see Eq. (3)). A mass-sheet transformation is defined as a change of the deflection angles (and correspondingly the lensing mass distributions) such that the lens equations remain invariant, with only a uniform isotropic scaling in the source planes. Considering the first of Eq. (5), a MST for the first source plane is obtained by applying the single-lens plane MST to the effective deflection βα1, i.e., by setting (here and in the following, a prime denotes transformed quantities) (6)after which the transformed position at z2 becomes (7)i.e., the required uniform isotropic scaling. The question now is whether we can find a transformation of the deflection angle α2 such that the second of Eq. (5) also remains invariant up to a scaling of θ3. If denotes the transformed deflection angle, then the transformed lens equation reads (8)Requiring that , corresponding to uniform scaling with the factor ν3, we obtain (9)To satisfy this equation, the term on the l.h.s. of Eq. (9) must contain the term ν3α2(θ2). Therefore, we set (10)where in the second step we used and the first of Eq. (5). Here, K2 is a constant, to be constrained later. Inserting Eq. (10) into Eq. (9), we see that the terms involving α2 vanish. Comparing the remaining terms proportional to α1 and θ, we find the pair of constraints λ − K2λβ = ν3 and β + λ − 1 − K2λβ = βν3, which have the unique solution (11)Hence we obtain a solution of Eq. (9) and accordingly a transformation of the mass distributions in both lens planes, which leads to at most a uniform isotropic scaling of the source planes.
The transformation of the first lens plane is a normal MST, in that the deflection angle (and thus the surface mass density) is scaled by an overall factor λ, and a uniform density K1 = (1 − λ) /β is added, leading to a scaling of the first source plane by a factor λ. This scaling equally applies to the light distribution and the mass distribution, as seen by Eq. (10). Thus, after the MST, both the light and the mass in the plane i = 2 are transformed in exactly the same way. In particular, this means that if the original model has a mass component centered on a light component, the same remains true after the MST, but both are located at a position that differs by a factor of λ.
Condition (10) furthermore implies that the transformed deflection in the second lens plane is a scaled version of the original deflection, plus a contribution from a uniform mass sheet K2. The implication of the scaling of the first term in Eq. (10) on the corresponding mass distribution can be seen as follows: if (12)is the deflection caused by the mass distribution κ, then (13)Hence, the transformed mass distribution is a scaled version of the original one, multiplied by a factor λ-1, plus a uniform mass sheet. It must be stressed here that the phrase “adding a uniform mass sheet” corresponds to a global interpretation of the MST; however, only a relatively small inner region of the lens is probed by strong lensing, and hence the MST needs to apply only locally. Its main effect is the change of the local slope of the mass profile near the Einstein radius of a lens, making it flatter (steeper) for λ< 1 (λ> 1). The global interpretation of the mass sheet as an “external convergence”, relating it to the large-scale environment of the lens which can be probed by the observed density field of galaxies (e.g., Wong et al. 2011; Collett et al. 2013; Greene et al. 2013), constitutes an extrapolation over a vast range of scales.
Surprisingly, the second source plane remains unscaled under this MST, because of ν3 = 1. Hence, the MST implies no change in the mapping of the second source plane, including no change in the magnification matrix. Whereas simple algebra has straightforwardly led to this result, the geometrical reason for this appears unclear.
If the MST is such that the scaling in the first lens plane corresponds to an additional focusing, which means λ< 1, so that the uniform mass sheet has positive convergence, then the mass sheet in the second lens plane has negative convergence, since the sign of K2 is opposite to that of 1 − λ. In particular, λ = 1 implies K2 = 0, so there is no MST that leaves the first lens plane invariant and only affects the second one. It must be stressed that the MST leaves the global lens mapping invariant and thus applies to sources of (in principle) arbitrary extent. This is quite different from other modifications of the lens mass distribution (e.g., Coe et al. 2008; Liesenborgs et al. 2008; Liesenborgs & De Rijcke 2012, and references therein), which apply to a discrete set of isolated small images of sources.
We now consider the impact of the MST on the time delay. For sources at z2, the MST is a normal single-plane MST, and the time delays are changed by a factor λ. For sources at z3, we use the expression for the light travel-time from θ3 via θ2 and θ to the observer, as given in Schneider et al. (1992), (14)where τi,i + 1 = (θi − θi + 1)2/ 2 − βi,i + 1ψi(θi) is the Fermat potential corresponding to neighboring planes, with β12 ≡ β and β23 = 1, and ψi(θi) is the deflection potential with ∇ψi = αi. We now consider how T behaves under an MST. The scalings (6) and (10) of αi imply that (15)Together with and , we obtain (16)Hence, the light travel-time function transforms as
The second term in Eq. (17) only depends on the source position θ3; therefore, this term does not contribute to the time delay, which is obtained as the difference of T between images. The third term of Eq. (17) vanishes, since the expression in the bracket is zero because of relations between the D’s – see Schneider et al. (1992). We thus find that all time delays in a two-source plane lens scale with λ under an MST, for sources on both source planes. Hence, time-delay ratios do not break the degeneracy of the MST; on the other hand, any time delay from a source in either source plane breaks it, provided H0 is assumed to be known.
In the light of the new MST, we now consider the possibility of constraining cosmological parameters from two-source plane lenses. In the approach of CA14, the sensitivity to cosmology is due to the distance ratio parameter β, which depends on the density parameters and the dark energy e.o.s. parameter w. Different cosmological models yield different distance-redshift relations, and thus different β. Thus let β correspond to a fiducial cosmological model, and β′ to a model with different parameters. We examine below whether we can find deflection angles such that the lens mappings are unchanged between the fiducial and modified models, up to a uniform scaling (by a factor νi) in the lens/source planes. Thus we require for the first lens equation (18)With the ansatz , we obtain by comparing terms in Eq. (18) proportional to α1 and θ the two relations β′λ = ν2β and 1 − β′K1 = ν2, with solutions ν2 = β′λ/β and K1 = 1 /β′ − λ/β. The requirement for the second lens equation reads (19)Inserting the ansatz
into Eq. (19), the comparison of terms propotional to α1 and θ yields the equations λ − K2ν2β = ν3 and 1 − K1 − K2ν2 = ν3, which have the solutions (20)Hence we see that a change of cosmology changes the lens equations in a similar way as an MST. Even with a change in β, we retain the freedom of a one-parameter family of mass models that leave the lens mapping invariant, up to a uniform scaling. In contrast to the MST discussed in Sect.2.3, here a non-trivial scaling of the second lens plane is implied, where ν3 solely depends on change of β.
MST parameters for modified cosmological parameters.
Three special choices of λ are worth to be discussed separately: (1) No mass sheet in the first lens plane (NMS1): for λ = β/β′, K1 = 0 and ν2 = 1, so that there is no mass sheet in the first lens and no scaling in the first source/second lens plane is implied. The mass sheet in the second lens plane then becomes K2 = (β′ − β)/ [β′(1 − β)]. For four different variations around a fiducial cosmological model, all with zero spatial curvature, we give the values of β′, K2 and ν3 in Table1, using the redshifts for the system SDSSJ0946+1006 (setting the uncertain redshift of the second source to be z3 = 2.4). (2) No mass sheet in the second lens plane (NMS2): setting λ = ν3, we derive K2 = 0, ν2 = (1 − β′)/(1 − β), and K1 = (β′ − β)/ [β′(1 − β)], so that λ + K1 = 1. Hence, in this case, there is no mass sheet in the second lens plane, and the one in the first lens plane is the same as K2 was for NMS1. (3) Equal-mass scaling (EMS): choosing sheets of equal density in both lens planes, we derive , , and . For these choices of λ, we also list some quantities for the different cosmological parameters in Table1.
From the table, we see that variations of the cosmological parameters within generous plausible ranges yield only small deviations of β′ from β. Correspondingly, the required MSTs are small; for the three special choices discussed above they imply mass sheets with | Ki | < 0.06. In particular, for EMS, the mass-sheet densities are smaller than 0.03. Such low values of Ki are perfectly acceptable, even if fairly accurate measurements of the velocity dispersion of the lenses were available (see Schneider & Sluse 2013, for a more detailed discussion). In particular, if the original mass distribution were a power law, the transformed ones would deviate only very little from a power law over the range where multiple images are formed.
We have shown that a lens system with two lens and two source planes admits a mass-sheet transformation that leaves the lens mapping invariant up to a uniform scaling in the source and lens plane(s). Furthermore, we demonstrated in the previous section that a change of cosmological parameters is essentially equivalent to an MST. A mass sheet acts like a magnifying glass, changing apparent distances, much in the same way as a different expansion history changes the distance-redshift relation. For the particular case of SDSSJ0946+1006 (Gavazzi et al. 2008), we have shown that the changes of the mass distribution in the second lens plane are very small when the cosmological parameters are changed within currently acceptable ranges.
As explicitly stated in CA14, the cosmological constraints they obtained are based on the assumption of a power-law mass distribution in either lens planes. This assumption about the functional form of the density profile formally breaks the degeneracy implied by the MST, but there are no good reasons to assume that the profiles of galaxies are exactly described by a power law (see Schneider & Sluse 2013). Indeed, the case NMS1 discussed above implies no added mass sheet in the first lens plane, thus preserving the power-law property, whereas the case EMS corresponds to only small Ki, yielding only marginal modifications of the power-law mass profile. The MST also for multiple source and lens planes provides a degeneracy in the determination of lens-mass distributions and cosmological distance ratios in lens systems, and needs to be accounted for in future studies.
The author thanks Dominique Sluse, Tom Collett, Sherry Suyu, and an anonymous referee for helpful discussions and comments on the manuscript. This work was supported in part by the Deutsche Forschungsgemeinschaft under the TR33 “The Dark Universe”.
- Bartelmann, M. 2010, Class. Quant. Gravity, 27, 3001 [NASA ADS] [CrossRef] [Google Scholar]
- Bradač, M., Lombardi, M., & Schneider, P. 2004, A&A, 424, 13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Coe, D., Fuselier, E., Benítez, N., et al. 2008, ApJ, 681, 814 [NASA ADS] [CrossRef] [Google Scholar]
- Collett, T. E., & Auger, M. W. 2014, MNRAS, 443, 969 [NASA ADS] [CrossRef] [Google Scholar]
- Collett, T. E., Marshall, P. J., Auger, M. W., et al. 2013, MNRAS, 432, 679 [NASA ADS] [CrossRef] [Google Scholar]
- Falco, E. E., Gorenstein, M. V., & Shapiro, I. I. 1985, ApJ, 289, L1 [NASA ADS] [CrossRef] [Google Scholar]
- Gavazzi, R., Treu, T., Koopmans, L. V. E., et al. 2008, ApJ, 677, 1046 [NASA ADS] [CrossRef] [Google Scholar]
- Greene, Z. S., Suyu, S. H., Treu, T., et al. 2013, ApJ, 768, 39 [NASA ADS] [CrossRef] [Google Scholar]
- Kochanek, C. S. 2006, in Saas-Fee Advanced Course 33: Gravitational Lensing: Strong, Weak and Micro, eds. G. Meylan, P. Jetzer, P. North, et al., 91 [Google Scholar]
- Liesenborgs, J., & De Rijcke, S. 2012, MNRAS, 425, 1772 [NASA ADS] [CrossRef] [Google Scholar]
- Liesenborgs, J., de Rijcke, S., Dejonghe, H., & Bekaert, P. 2008, MNRAS, 386, 307 [NASA ADS] [CrossRef] [Google Scholar]
- Refsdal, S. 1964, MNRAS, 128, 307 [NASA ADS] [CrossRef] [Google Scholar]
- Schneider, P., & Sluse, D. 2013, A&A, 559, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schneider, P., Ehlers, J., & Falco, E. E. 1992, Gravitational Lenses (New York, Heidelberg: Springer Verlag) [Google Scholar]
- Suyu, S. H., Auger, M. W., Hilbert, S., et al. 2013, ApJ, 766, 70 [NASA ADS] [CrossRef] [Google Scholar]
- Treu, T. 2010, ARA&A, 48, 87 [NASA ADS] [CrossRef] [Google Scholar]
- Wong, K. C., Keeton, C. R., Williams, K. A., Momcheva, I. G., & Zabludoff, A. I. 2011, ApJ, 726, 84 [NASA ADS] [CrossRef] [Google Scholar]
Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.