**Gravitational Lensing in Astronomy**

Joachim Wambsganss

Astrophysikalisches Institut Potsdam

An der Sternwarte 16

14482 Potsdam, Germany

jwambsganss@aip.de

Review article for *Living Reviews in Relativity*

http://www.livingreviews.org

Abstract

Deflection of light by gravity was predicted by General Relativity and observationaly confirmed in 1919. In the following decades various aspects of the gravitational lens effect were explored theoretically, among them the possibility of multiple or ring-like images of background sources, the use of lensing as a gravitational telescope on very faint and distant objects, and the possibility to determine Hubble’s constant with lensing. Only relatively recently gravitational lensing became an observational science after the discovery of the first doubly imaged quasar in 1979. Today lensing is a booming part of astrophysics.

In addition to multiply-imaged quasars, a number of other aspects of lensing have been discovered since, e.g. giant luminous arcs, quasar microlensing, Einstein rings, galactic microlensing events, arclets, or weak gravitational lensing. By now literally hundreds of individual gravitational lens phenomena are known.

Although still in its childhood, lensing has established itself as a very useful astrophysical tool with some remarkable successes. It has contributed significant new results in areas as different as the cosmological distance scale, the large scale matter distribution in the universe, mass and mass distribution of galaxy clusters, physics of quasars, dark matter in galaxy halos, or galaxy structure. Looking at these successes in the recent past we predict an even more luminous future for gravitational lensing.

1

Introduction
============

Within the last 20 years gravitational lensing has changed from being considered a geometric curiosity to a helpful and in some ways unique tool of modern astrophysics. Although the deflection of light at the solar limb was very successfully hailed as the first experiment to confirm a prediction of Einstein’s theory of General Relativity in 1919, it took more than half a century to establish this phenomenon observationally in some other environment. By now almost a dozen different realizations of lensing are known and observed, and surely more will show up.

Gravitational lensing – the attraction of light by matter – displays a number of attractive features as an academic discipline. Its principles are very easy to understand and to explain due to its being a geometrical effect. Its ability to produce optical illusions is fascinating to scientists and laypeople alike. And – most importantly of course – its usefulness for a number of astrophysical problems makes it an attractive tool in many branches of astronomy. All three aspects will be considered below.

In the not quite two decades of its existence as an observational branch of astrophysics, the field of gravitational lensing has been continuously growing. Every few years a new realisation of the phenomenon was discovered: multiple quasars, giant luminous arcs, quasar microlensing, Einstein rings, galactic microlensing, weak lensing, galaxy-galaxy lensing opened up very different regimes for the gravitational telescope. This trend is reflected in a growing number of people working in the fi2ld. In Figure \[fig-lens-papers\] the number of publications in scientific journals that deal with gravitational lensing is plotted over time: It is obvious that lensing is booming.

(15.0,12.0) (2.0,-0.5) [=12.0cm ]{}

Although there had been a slight sense of disappointment in the astronomical community a few years ago about the fact that lensing had not yet solved all the big problems of astrophysics right away (e.g. determination of the Hubble constant; nature of dark matter; physics/size of quasars) this feeling has apparently been reversed by now. With its many applications and quantitative results lensing starts to fulfill its astrophysical promises.

We shall start with a brief look back in time and mention some historic aspects of light deflection and lensing in Chapter \[sec-hist\]. We then attempt to explain the basic features of gravitational lensing quantitatively, deriving some of the relevant equations (Chapter \[sec-basi\]). A whole variety of lensing observations and phenomena which curved space-time provides for us is presented in Chapter \[sec-phen\], e.g. multiple versions of quasars, gigantically distorted images of galaxies, and highly magnified stars. Along that we explain and discuss the astrophysical applications of lensing which show the use of this tool. This chapter will be the most detailed one. Finally in the concluding Chapter \[sec-futu\] we try to extrapolate and speculate about the future development of the field.

This article cannot cover the whole field of lensing and its applications. Therefore we list here a number of books, proceedings and other (partly complementary) review articles on the subject of gravitational lensing.

The textbook by Schneider, Ehlers & Falco [@SEF] contains the most comprehensive presentation of gravitational lensing. A new edition is underway. The book by Bliokh & Minakov [@BM89] on gravitational lensing is still only available in Russian. A new book currently in press by Petters, Levine & Wambsganss [@PLW] treats mainly the mathematical aspects of lensing, in particular its applications to singularity theory.

The contributions to the most important conferences on gravitational lensing in the last few years were all published: Swings [@liege83] edited the Proceedings on the first conference on lensing in Liège in 1983; Moran et al. [@MIT] are the editors of the MIT workshop on lensing in 1988; Mellier et al. [@toulouse] of the Toulouse conference in 1989; Kayser et al. [@hamburg] of the Hamburg meeting in 1991; Surdej et al. [@liege93] of the Liège conference in 1993; Kochanek & Hewitt [@melbourne] of the IAU Symposium 173 in Melbourne in 1995; Jackson [@Jac97] of the Jodrell Bank Meeting “Golden Lenses" in 1997.

There exist a number of excellent reviews on gravitational lensing. Blandford & Kochanek [@BK87] give a nice introduction on the theory of lensing. The optical aspects of lensing are derived elegantly in [@Bla89]. The presentation of Blandford & Narayan [@BN92] emphasizes in particular the cosmological applications of gravitational lensing. The review by Refsdal & Surdej [@RS94] contains a section on optical model lenses that simulate the lensing action of certain astrophysical objects. A recent review article by Narayan & Bartelmann [@NB96] summarizes in a very nice and easy-to-understand way the basics and the latest in the gravitational lens business. In the sections below some more specific review articles will be mentioned.

History of Gravitational Lensing {#sec-hist}
================================

The first written account of the deflection of light by gravity appeared in the “Berliner Astronomisches Jahrbuch auf das Jahr 1804" in an article entitled: “Ueber die Ablenkung eines Lichtstrals von seiner geradlinigen Bewegung, durch die Attraktion eines Weltkörpers, an welchem er nahe vorbeigeht"[^1] [@Sol1801]. Johann Soldner – a German geodesist, mathematician and astronomer, then working at the Berlin Observatory – explored this effect and inferred that a light ray close to the solar limb would be deflected by an angle $\tilde\alpha = 0.84$ arcsec. It is very interesting to read how carefully and cautiously he investigated this idea and its consequences on practical astronomy.

In the year 1911 – more than a century later – Albert Einstein [@Ein11] directly addressed the influence of gravity on light (“Über den Einflu[ß]{} der Schwerkraft auf die Ausbreitung des Lichtes"[^2]). At this time the General Theory of Relativity was not fully developed yet. This is the reason why Einstein obtained – unaware of the earlier result – the same value for the deflection angle as Soldner had calculated with Newtonian physics. In this paper Einstein found $\tilde\alpha = 2 G M_{\odot} / c^2 R_{\odot}  =  0.83$ arcsec for the deflection angle of a ray grazing the sun (here $M_{\odot}$ and $R_{\odot}$ are the mass and the radius of the sun, $c$ and $G$ are the velocity of light and the gravitational constant, respectively). Einstein emphasized his wish that astronomers investigate this question (“Es wäre dringend zu wünschen, da[ß]{} sich Astronomen der hier aufgerollten Frage annähmen, auch wenn die im vorigen gegebenen Überlegungen ungenügend fundiert oder gar abenteuerlich erscheinen sollten."[^3]). Recently it was discovered that Einstein had derived the lens equation, the possibility of a double image and the magnifications of the images in a notebook in the year 1912 [@Ren97]. In 1913 Einstein even contacted directly the director of the Mt. Wilson Observatory – George Ellery Hale – and asked him whether it would be possible to measure positions of stars near the sun during the day in order to establish the deflection effect of the sun (see [@AIP]). Einstein’s hand-drawn sketch marks the the “wrong" value of 0.84 arcsec. There actually were plans to test Einstein’s wrong prediction of the deflection angle during a solar eclipse in 1914 on the Russian Crimea peninsula. However, when the observers were already in Russia, World War I broke out and they were captured by Russian soldiers [@Bre82]. So – fortunately for Einstein – the measurement of the deflection angle at the solar limb had to be postponed for a few years.

With the completion of the General Theory of Relativity Einstein was the first to derive the correct deflection angle $\tilde \alpha$ of a light ray passing at a distance $r$ from an object of mass $M$ as $$\tilde \alpha  = {4 G M \over c^2}{1 \over  r}.$$ The additional factor of two (compared to the “Newtonian" value) reflects the spatial curvature which is missed if photons are just treated as particles. With the solar values for radius and mass Einstein obtained [@Ein16; @Ein22]: $$\tilde \alpha_{\odot} = { {4 G M_{\odot}} \over  { c^2}} 
            {1 \over  R_{\odot}} =  1.74 {\rm \ arcsec}.$$ It is common wisdom now that the determination of this value to within 20% during the solar eclipse in 1919 by Arthur Eddington and his group was the second observational confirmation of General Relativity [@Dys20] and the basis of Einstein’s huge popularity starting in the 1920s (the first one had been the explanation of Mercury’s perihelion shift). Recently the value predicted by Einstein was confirmed to an accuracy better than 0.02 % [@Leb95].

In the following decades light deflection or gravitational lensing was only very rarely the topic of a research paper: In 1924 Chwolson [@Chw24] mentioned the idea of a “fictitous double star" and the mirror-reversed nature of the secondary image. He also mentioned the symmetric case of star exactly behind star, resulting in a circular image. Einstein also reported in 1936 about the appearance of a “luminous circle" for perfect alignment between source and lens [@Ein36], and of two magnified images for slightly displaced positions[^4]. Today such a lens configuration is called “Einstein-ring", although more correctly it should be called “Chwolson-ring". Influenced by Einstein, Fritz Zwicky [@Zwi37a; @Zwi37b] pointed out in 1937 that galaxies (“extragalactic nebulae") are much more likely to be gravitationally lensed than stars and that one can use the gravitational lens effect as a “natural telescope".

In the 1960s a few partly independent theoretical studies showed the usefulness of lensing for astronomy [@Kli63; @Lie64; @Lin67; @Met63; @Ref64a; @Ref64b]. In particular Sjur Refsdal derived the basic equations of gravitational lens theory and subsequently showed how the gravitational lens effect can be used to determine Hubble’s constant by measuring the time delay between two lensed images. He followed up this work with interesting applications of lensing [@Ref66a; @Ref66b; @Ref70]. The mathematical foundation of how a light bundle is distorted on its passage throught the universe had been derived in the context of gravitational radiation even before [@Sac61].

Originally gravitational lensing was discussed for stars or for galaxies. When in the 1960s quasars were discovered, Barnothy [@Bar65] was the first to connect them with the gravitational lens effect. In the late 60s/early 70s a few groups and individuals explored various aspects of lensing further, e.g.: statistical effects of local inhomogeneities on the propagation of light [@Gun67a; @Gun67b; @PG73]; lensing applied to quasars and clusters of galaxies [@CK75; @Noo71; @San71]; development of a formalism for transparent lenses [@BK75; @Cla72]; or the effect of an inhomogeneous universe on the distance-redshift relations [@DR73].

But only in the year 1979 the whole field received a real boost when the first double quasar was discovered and confirmed to be a real gravitational lens by Walsh, Carswell & Weymann [@WCW]. This discovery and the development of lensing since then will be described in Chapter \[sec-phen\].

There are a few historic accounts of lensing which are more detailed than the one presented here, e.g. in [@Rob83], in [@SEF], or in [@Wam90]. The complete history of gravitational lensing has to be written yet.

Basics of Gravitational Lensing {#sec-basi}
===============================

The path, the size and the cross section of a light bundle propagating through spacetime in principle are affected by all the matter between the light source and the observer. For most practical purposes we can assume that the lensing action is dominated by a single matter inhomogeneity at some location between source and observer. This is usually called the “thin lens approximation": all the action of deflection is thought to take place at a single distance. This approach is valid only if the relative velocities of lens, source and observer are small compared to the velocity of light $v \ll c$ and if the Newtonian potential is small $|\Phi| \ll c^2$. These two assumptions are justified in all astronomical cases of interest. The size of a galaxy, e.g., is of order 50 kpc, even a cluster of galaxies is not much larger than 1 Mpc. This “lens thickness" is small compared to the typical distances of order few Gpc between observer and lens or lens and background quasar/galaxy, respectively. We assume that the underlying spacetime is well described by a perturbed Friedmann-Robertson-Walker metric[^5]:

Lens Equation  {#lens-equation .unnumbered}
--------------

The basic setup for such a simplified gravitational lens scenario involving a point source and a point lens is displayed in Figure \[fig-setup\]. The three ingredients in such a lensing situation are the source S, the lens L, and the observer O. Light rays emitted from the source are deflected by the lens. For a point-like lens, there will always be (at least) two images S$_1$ and S$_2$ of the source. With external shear – due to the tidal field of objects outside but near the light bundles – there can be more images. The observer sees the images in directions corresponding to the tangents to the real incoming light paths.

(14.0,11.0) (0.5,-0.5) [=12.0cm ]{}

In Figure \[fig-setup-eq\] the corresponding angles and angular diameter distances $D_L$, $D_S$, $D_{LS}$ are indicated[^6]. In the thin-lens approximation the hyperbolic paths are approximated by their asymptotes. In the circular-symmetric case the deflection angle is given as

$$\label{eq-angle}
\tilde \alpha (\xi)  =  { { 4 G M(\xi)} \over {c^2} } { 1  \over \xi }.$$

where $M(\xi)$ is the mass inside a radius $\xi$. In this depiction the origin is chosen at the observer. From the diagram it can be seen that the following relation holds: $$\theta D_S = \beta D_S + \tilde \alpha  D_{LS}$$ (for $\theta$, $\beta$, $\tilde \alpha \ll 1$; this condition is fulfilled in practically all astrophysically relevant situations). With the definition of the reduced deflection angle as $\alpha (\theta)  =   ( D_{LS} / D_{S} ) \tilde \alpha (\theta)$, this can be expressed as: $$\label{eq-lenseq}
       \beta = \theta - \alpha (\theta).$$ This relation between the positions of images and source can easily be derived for a non-symmetric mass distribution as well. In that case all angles are vector-valued. The two-dimensional [**lens equation**]{} then reads: $$\label{eq-lenseq-vec}
       \vec\beta = \vec\theta - \vec\alpha (\vec\theta).$$

(14.0,11.0) (2.5,-0.9) [=12.0cm ]{}

Einstein Radius {#einstein-radius .unnumbered}
---------------

For a point lens of mass $M$ the deflection angle is given by equation (\[eq-angle\]). Plugging into equation (\[eq-lenseq\]) and using the relation $\xi = D_L \theta $ (cf. Figure \[fig-setup-eq\]) one obtains: $$\beta(\theta)  = 
       \theta - { D_{LS} \over D_L D_S} { 4 G M \over c^2 \theta}.$$ For the special case in which the source lies exactly behind the lens ($\beta = 0$), due to the symmetry a ring-like image occurs whose angular radius is called [**Einstein radius**]{} $\theta_E$:

$$\theta_E  = \sqrt{               
            { 4 G M \over c^2}  
            {D_{LS} \over D_L D_S}.
                }$$

The Einstein radius defines the angular scale for a lens situation. For a massive galaxy with a mass of $M = 10^{12} M_{\odot}$ at a redshift of $z_L = 0.5$ and a source at redshift $z_S = 2.0$ (we used here $H = 50 $km sec$^{-1}$ Mpc$^{-1}$ as the value of the Hubble constant and an Einstein-deSitter universe) the Einstein radius is $$\label{eq-angle-gal}
       \theta_E  \approx 1.8  \ \sqrt{ M \over 10^{12} M_\odot } \ {\rm arcsec}$$ (note that for cosmological distances in general $D_{LS} \ne D_S - D_L$!). For a galactic microlensing scenario in which stars in the disk of the Milky Way act as lenses for bulge stars close to the center of the Milky Way, the scale defined by the Einstein radius is $$\theta_E  \approx 0.5  \ \sqrt{ M \over M_{\odot} }   \ {\rm milliarcsec}.$$ An application and some illustrations of the point lens case can be found in Section \[sec-gal-micro\] on galactic microlensing.

Critical Surface Mass Density  {#critical-surface-mass-density .unnumbered}
------------------------------

In the more general case of a three-dimensional mass distribution of an extended lens, the density $\rho(\vec r)$ can be projected along the line of sight onto the lens plane to obtain the two-dimensional surface mass density distribution $\Sigma(\vec\xi)$, as $$\Sigma(\vec\xi) = \int_0^{D_S} \rho(\vec r) dz.$$ Here $ \vec r $ is a three-dimensional vector in space, and $ \vec \xi$ is a two-dimensional vector in the lens plane. The two-dimensional deflection angle $\vec{\tilde \alpha}$ is then given as the sum over all mass elements in the lens plane: $$\vec{ \tilde \alpha} (\vec{\xi})=  { {4 G} \over {c^2}}
    \int { { (\vec\xi - \vec\xi') \Sigma(\vec\xi')}  \over 
    { |\vec\xi - \vec\xi'|^2 }  }d^2\xi'.$$ For a finite circle with constant surface mass density $\Sigma$ the deflection angle can be written: $${        \alpha (\xi) } =   { D_{LS} \over D_S} 
            { 4 G \over c^2    } 
            { \Sigma \pi \xi^2 \over     \xi}$$ With $\xi = D_L \theta$ this simplifies to $${        \alpha (\theta) } =   { 4 \pi G \Sigma \over c^2 } 
    { D_L D_{LS} \over D_S} \theta.$$ With the definition of the [**critical surface mass density**]{} $\Sigma_{\rm crit}$ as $$\label{eq-crit}
\Sigma_{\rm crit}  =  { c^2  \over {4 \pi G }} { D_S  \over {D_L D_{LS} }}$$ the deflection angle for a such a mass distribution can be expressed as $${ \tilde \alpha (\theta) } = { \Sigma \over \Sigma_{\rm crit} } \theta.$$ The critical surface mass density is given by the lens mass $M$ “smeared out" over the area of the Einstein ring: $\Sigma_{crit} = M / (R_E^2 \pi)$, where $R_E = \theta_E D_L$. The value of the critical surface mass density is roughly $\Sigma_{\rm crit}  \approx 0.8$ g cm$^{-2}$ for lens and source redshifts of $z_L = 0.5$ and $z_S = 2.0$, respectively. For an arbitrary mass distribution the condition $ \Sigma > \Sigma_{\rm crit}$ at any point is sufficient to produce multiple images.

Image Positions and Magnifications  {#image-positions-and-magnifications .unnumbered}
-----------------------------------

The lens equation (\[eq-lenseq\]) can be re-formulated in the case of a single point lens:

$$\label{eq-lenseq-point}
       \beta = \theta - { \theta_E^2  \over \theta}.$$

Solving this for the image positions $\theta$ one finds that an isolated point source always produces two images of a background source. The positions of the images are given by the two solutions: $$\label{eq-lenseq-images}
       \theta_{1,2} = {1 \over 2 } 
       \left( \beta  \pm \sqrt{\beta^2 + 4 \theta_E^2} \right).$$ The magnification of an image is defined by the ratio between the solid angles of the image and the source, since the surface brightness is conserved. Hence the magnification $\mu$ is given as $$\label{eq-magn}
       \mu  = {\theta \over \beta} {d \theta \over d \beta}.$$ In the symmetric case above the image magnification can be written as (by using the lens equation): $$\mu_{1,2}  
       = 
       \left( 1 - \left[ {\theta_E \over \theta_{1,2}} \right]^4  \right)^{-1}
       = 
           { u^2 + 2 \over 2 u \sqrt{u^2 + 4 } } \pm {1 \over 2}$$ Here we defined $u$ as the “impact parameter", the angular separation between lens and source in units of the Einstein radius: $u = \beta / \theta_E$. The magnification of one image (the one inside the Einstein radius) is negative. This means it has negative parity: it is mirror-inverted. For $\beta \rightarrow 0$ the magnification diverges: in the limit of geometrical optics the Einstein ring of a point source has infinite magnification[^7]! The sum of the absolute values of the two image magnifications is the measurable total magnification $\mu$: $$\label{eq-magn-tot}
       \mu = |\mu_{1}| + |\mu_2| = { u^2 + 2 \over   u \sqrt{u^2 + 4 } }.$$ Note that this value is (always) larger than one[^8]! The difference between the two image magnifications is unity: $$\mu_{1} + \mu_2 = 1.$$

(Non-)Singular Isothermal Sphere  {#non-singular-isothermal-sphere .unnumbered}
---------------------------------

A handy and popular model for galaxy lenses is the singular isothermal sphere with a three-dimensional density distribution of $$\rho(r) = { \sigma_v^2 \over 2 \pi G} {1 \over r^2},$$ where $\sigma_v$ is the one-dimensional velocity dispersion. Projecting the matter on a plane one obtains the circularly-symmetric surface mass distribution $$\Sigma(\xi) = { \sigma_v^2 \over 2 G} {1 \over \xi }.$$ With  $M(\xi) = \int_0^\xi \Sigma(\xi') 2 \pi \xi' d\xi' $   plugged into equation (\[eq-angle\]) one obtains the deflection angle for an isothermal sphere, which is a constant (i.e. independent of the impact parameter $\xi$): $$\tilde \alpha(\xi) = 4 \pi { \sigma_v^2 \over c^2 }.$$ In “practical units" for the velocity dispersion this can be expressed as: $$\tilde \alpha(\xi) 
    = 1.15 \left( \sigma_v \over 200 {\rm km s}^{-1} \right)^2  
    {\rm arcsec} .$$ Two generalizations of this isothermal model are commonly used: Models with finite cores are more realistic for (spiral) galaxies. In this case the deflection angle is modified to (core radius $\xi_c$): $$\tilde \alpha(\xi) 
= 4 \pi { \sigma_v^2 \over c^2 } {\xi \over {(\xi_c^2 + \xi^2)^{1/2}}}.$$ Furthermore, a realistic galaxy lens usually is not perfectly symmetric but is slightly elliptical. Depending on whether one wants an elliptical mass distribution or an elliptical potential, various formalisms have been suggested. Detailed treatments of elliptical lenses can be found in [@bar98; @BK87; @KK93; @KK98; @Kor94; @Schramm].

Lens Mapping  {#lens-mapping .unnumbered}
-------------

In the vicinity of an arbitrary point, the lens mapping as shown in equation (\[eq-lenseq-vec\]) can be described by its Jacobian matrix $\cal A$: $${\cal A} = { \partial \vec \beta \over \partial \vec \theta} = 
         \left( \delta_{ij} 
        - {\partial \alpha_i (\vec \theta) \over \partial \theta_j} \right)
     = 
         \left( \delta_{ij} 
     - {\partial^2 \psi (\vec \theta) 
        \over \partial \theta_i \partial \theta_j} \right).$$ Here we made use of the fact (see [@Bla86; @Sch85]), that the deflection angle can be expressed as the gradient of an effective two-dimensional scalar potential $\psi$:  $\vec \nabla_{\theta} \psi = \vec \alpha$, where $$\psi(\vec \theta) = { D_{LS} \over {D_L D_S} } {2 \over c^2}
    \int \Phi (\vec r) dz$$ and $\Phi (\vec r) $ is the Newtonian potential of the lens.

The determinant of the Jacobian $\cal A$ is the inverse of the magnification: $$\label{eq-magn-det}
    \mu = {1 \over \det {\cal A}}.$$ Let us define $$\psi_{ij} = {\partial^2 \psi  \over \partial \theta_i \partial \theta_j}.$$ The Laplacian of the effective potential $\psi$ is twice the convergence: $$\label{eq-kappa}
 \psi_{11} + \psi_{22}  = 2 \kappa = {\rm tr } \    \psi_{ij} .$$ With the definitions of the components of the external shear $\gamma$: $$\label{eq-gamma1}
  \gamma_1 (\vec \theta) = {1 \over 2} (\psi_{11} - \psi_{22}) = 
  \gamma (\vec \theta) \cos [2 \varphi (\vec \theta)]$$ and $$\label{eq-gamma2}
  \gamma_2 (\vec \theta) =  \psi_{12}   =  \psi_{21} = 
  \gamma (\vec \theta) \sin [2 \varphi (\vec \theta)]$$ (where the angle $\varphi$ reflects the direction of the shear-inducing tidal force relative to the coordinate system) the Jacobian matrix can be written $${\cal A} = 
     \left( 
    \begin{array}{cc} 
            1 - \kappa - \gamma_1 & -\gamma_2 \\
                           -\gamma_2  &  1 - \kappa + \gamma_1 \\
    \end{array} 
    \right)
    =  ( 1 - \kappa )
     \left( 
    \begin{array}{cc} 
            1 & 0 \\
            0 & 1 \\ 
    \end{array} 
    \right)
    - \gamma
     \left( 
    \begin{array}{lr} 
            \cos 2\varphi  &  \sin 2\varphi  \\
            \sin 2\varphi  & -\cos 2\varphi  \\ 
    \end{array} 
    \right).$$ The magnification can now be expressed as a function of the local convergence $\kappa$ and the local shear $\gamma$: $$\mu = ( \det {\cal A})^{-1}  =  { 1 \over (1-\kappa)^2 - \gamma^2}.$$ Locations at which   $\det A = 0$    have formally infinite magnification. They are called [**critical curves**]{} in the lens plane. The corresponding locations in the source plane are the [**caustics**]{}. For spherically symmetric mass distributions, the critical curves are circles. For a point lens, the caustic degenerates into a point. For elliptical lenses or spherically symmetric lenses plus external shear, the caustics can consist of cusps and folds. In Figure \[fig-lens-elliptical\] the caustics and critical curves for an elliptical lens with a finite core are displayed.

(15.0,8.0) (-0.5,-0.5) [= 8.0cm ]{} (7.5,-0.5) [= 8.0cm ]{}

Time delay and “Fermat’s" Theorem {#time-delay-and-fermats-theorem .unnumbered}
---------------------------------

The deflection angle is the gradient of an effective lensing potential $\psi$ (as was first shown by [@Sch85]; see also [@Bla86]). Hence the lens equation can be rewritten as $$(\vec \theta - \vec \beta) - \vec \nabla_\theta \psi  = 0$$ or $$\label{eq-nabla}
\vec \nabla_{\theta} 
    \left( {1 \over 2} (\vec \theta - \vec \beta)^2  - \psi \right) = 0.$$ The term in brackets appears as well in the physical time delay function for gravitationally lensed images: $$\label{eq-time-delay}
\tau(\vec \theta, \vec \beta) =  \tau_{\rm geom} + \tau_{\rm grav} =
    { 1 + z_L \over c } { D_L D_S \over D_{LS}}
  \left( {1 \over 2} (\vec \theta - \vec \beta)^2  - \psi(\theta)  \right).$$ This time delay surface is a function of the image geometry ($\vec \theta$, $\vec \beta$), the gravitational potential $\psi$, and the distances $D_L$, $D_S$, and $D_{LS}$. The first part – the geometrical time delay $\tau_{\rm geom}$ – reflects the extra path length compared to the direct line between observer and source. The second part – the gravitational time delay $\tau_{\rm grav}$ – is the retardation due to gravitational potential of the lensing mass (known and confirmed as Shapiro delay in the solar system). From equations (\[eq-nabla\]), (\[eq-time-delay\]) it follows that the gravitationally lensed images appear at locations that correspond to extrema in the light travel time, which reflects Fermat’s principle in gravitational-lensing optics.

The (angular-diameter) distances that appear in equation (\[eq-time-delay\]) depend on the value of the Hubble constant [@Wei72]; therefore it is possible to determine the latter by measuring the time delay between different images and using a good model for the effective gravitational potential $\psi$ of the lens (see [@Kun97; @Ref64b; @WS97] and Section \[sec-mult\]).

This chapter followed heavily the elegant presentation of the basics of lensing in Narayan & Bartelmann [@NB96]. Many more details can be found there. More complete derivations of the lensing properties are also provided in all the introductory texts mentioned in Chapter 1, in particular in [@SEF]. More on the formulation of gravitational lens theory in terms of time-delay and Fermat’s principle can be found in Blandford & Narayan [@Bla86] and Schneider [@Sch85]. Discussions of the concept of “distance" in relation to cosmology/curved space can be found in chapter 3.5 of [@SEF] or chapter 14.4 of [@Wei72].

Lensing Phenomena {#sec-phen}
=================

In this chapter we describe different groups of gravitational lens observations. The subdivision is pragmatic rather than entirely logical. It is done partly by lensed object, or by lensing object, or by lensing strength. The ordering roughly reflects the chronological appearance of different sub-disciplines to lensing. The following sections are on

-   Multiply-imaged quasars

-   Quasar microlensing

-   Einstein rings

-   Giant luminous arcs and arclets

-   Weak lensing / Statistical lensing

-   Cosmological aspects of (strong) lensing

-   Galactic microlensing

Comprehensive reviews could be written on each separate subject listed here. Hence the treatment here can be only very cursory.

Multiply-imaged quasars {#sec-mult}
-----------------------

In the year 1979 gravitational lensing became an observational science when the double quasar Q0957+561 was discovered - the first example of a lensed object [@WCW]. The discovery itself happened rather by accident; the discoverer Dennis Walsh describes in a nice account how this branch of astrophysics came into being [@Wal88].

It was not entirely clear at the beginning, though, whether the two quasar images really were an illusion provided by curved space-time – or rather physical twins. But intensive observations soon confirmed the almost identical spectra. The intervening “lensing" galaxy was found, and the “supporting" cluster was identified as well. Later very similar lightcurves of the two images (modulo offsets in time and magnitude) confirmed this system beyond any doubt as a bona fide gravitational lens.

By now about two dozen multiply-imaged quasar systems have been found, plus another ten good candidates (updated tables of multiply-imaged quasars and gravitational lens candidates are provided, e.g., by the CASTLE group [@CASTLE]). This is not really an exceedingly big number, considering a 20 year effort to find lensed quasars. The reasons for this “modest" success rate is

1.  Quasars are rare and not easy to find (by now roughly $10^4$ are known).

2.  The fraction of quasars that is lensed is small (less than one percent).

3.  It is not trivial at all to [*identify*]{} the lensed (i.e. multiply-imaged) quasars among the known ones.

Gravitationally lensed quasars come in a variety of classes: double, triple and quadruple systems; symmetric and asymmetric image configurations are known.

A recurring problem connected with double quasars is the question whether they are two images of a single source or rather a physical association of two objects (with three or more images it is more and more likely that it is lensed system). Few systems are as well established as the double quasar Q0957+561; but many are considered “safe" lenses as well. Criteria for “fair", “good", or “excellent" lensed quasar candidates comprise the following:

-   There are two or more point-like images of very similar optical color.

-   Redshifts (or distances) of both quasar images are identical or [*very*]{} similar.

-   Spectra of the various images are identical or [*very*]{} similar to each other.

-   There is a lens (most likely a galaxy) found between the images, with a measured redshift much smaller than the quasar redshift (for a textbook example see Figure \[fig-lens-galaxy\]).

-   If the quasar is intrinsically variable, the fluxes measured from the two (or more) images follow a very similar light curve, except for certain lags – the time delays – and an overall offset in brightness (cf. Figure \[fig-0957-light\]).

For most of the known multiple quasar systems only some of the above criteria are fully confirmed. And there are also good reasons not to require perfect agreement with this list: e.g., the lensing galaxy could be superposed to one quasar image and make the quasar appear extended; color/spectra could be affected by dust absorption in the lensing galaxy and appear not identical; the lens could be too faint to be detectable (or even a real dark lens?); the quasar could be variable on time scales shorter than the time delay; microlensing can affect the lightcurves of the images differently. Hence it is not easy to say how many gravitationally lensed quasar systems exist. The answer depends on the amount of certainty one requires. In a recent compilation Keeton & Kochanek [@KK96] put together 29 quasars as lenses or lens candidates in three probability “classes".

(14.0,8.0) (1.5, 0.1) [=12.0cm ]{}

Gravitationally lensed quasar systems are studied individually in great detail to get a better understanding of both lens and source (so that, e.g., a measurement of the time delay can be used to determine the Hubble constant). As an ensemble, the lens systems are also analysed statistically in order to get information about the population of lenses (and quasars) in the universe, their distribution in distance (i.e. cosmic time) and mass, and hence about the cosmological model (more about that in Section \[sec-cosmological\]). Here we will have a close look on one particularly well investigated system.

### The first Lens: Double Quasar Q0957+561 {#the-first-lens-double-quasar-q0957561 .unnumbered}

The quasar Q0957+561 was originally found in a radio survey, subsequently an optical counterpart was identified as well. After the confirmation of its lens nature [@WCW; @Wal88] this quasar attracted quite some attention. Q0957+561 has been looked at in all available wavebands, from X-rays to radio frequencies. More than 100 scientific papers have appeared on Q0957+561 (cf. [@PSV]), many more than on any other gravitational lens system. Here we will summarize what is known about this system from optical and radio observations.

In the optical light Q0957+561 appears as two point images of roughly 17 mag (R band) separated by 6.1 arcseconds (see Figure \[Q0957+561\_hst\]). The spectra of the two quasars reveal both redshifts to be $z_Q = 1.41$. Between the two images, not quite on the connecting line, the lensing galaxy (with redshift $z_G = 0.36$) appears as a fuzzy patch close to the component B. This galaxy is part of a cluster of galaxies at about the same redshift. This is the reason for the relatively large separation for a galaxy-type lens (typical galaxies with masses of $10^{11-12} M_{\odot}$ produce splitting angles of only about one arcsecond, see equation (\[eq-angle-gal\])). In this lens system the mass in the galaxy cluster helps to increase the deflection angles to this large separation.

(14.0,10.0) (3.0, 0.1) [=10.0cm ]{}

A recent image of Q0957+561 taken with the MERLIN radio telescope is shown in Figure \[Q0957+561\_merlin\]. The positions of the two point-like objects in this radio observation coincide with the optical sources. There is no radio emission detected at the position of the galaxy center, i.e. the lensing galaxy is radio-quiet. But this also tells us that a possible third image of the quasar must be very faint, below the detection limit of all the radio observations[^9]. In Figure \[Q0957+561\_merlin\] a “jet" can be seen emerging from image A (at the top). It is not unusual for radio quasars to have such a “jet" feature. This is most likely matter that is ejected from the central engine of the quasar with very high speed along the polar axis of the central black hole. The reason that this jet is seen only around one image is that it lies outside the caustic region in the source plane, which marks the part that is multiply imaged: only the compact core of the quasar lies inside the caustic and is doubly imaged.

(14.0,10.0) (2.5, 0.1) [=10.0cm ]{}

As stated above, a virtual “proof" of a gravitational lens system is a measurement of the “time delay" $\Delta t$, the relative shift of the light curves of the two or more images, $I_A(t)$ and $I_B(t)$, so that $I_B(t) = const \times I_A(t+\Delta t)$: any intrinsic fluctuation of the quasar shows up in both images, in general with an overall offset in apparent magnitude and an offset in time.

Q0957+561 is the first lens system in which the time delay was firmly established. After a decade long attempt and various groups claiming either of two favorable values [@Pel96; @PRH; @ST95; @Van89], Kundić et al. [@Kun97] confirmed the shorter of the two (cf. Figure \[fig-0957-light\]; see also Oscoz et al. [@Osc97] and Schild & Thomson [@ST97]): $$\Delta t_{Q0957+561} = (417\pm3) {\rm days}$$

(14.0,11.0) (2.5,-0.1) [=10.0cm ]{}

With a model of the lens system, the time delay can be used to determine the Hubble constant[^10]. In Q0957+561 the lensing action is done by an individual galaxy plus an associated galaxy cluster (to which the galaxy belongs). This provides some additional problems, a degeneracy in the determination of the Hubble constant [@GFS]: the appearance of the double quasar system including the time delay could be identical for different partitions of the matter between galaxy and cluster, but the derived value of the Hubble constant could be quite different. However, this degeneracy can be “broken", once the focussing contribution of the galaxy cluster can be determined independently. And the latter has been attempted recently [@Fis97]. The resulting value for the Hubble constant [@Kun97] obtained by employing a detailed lens model [@GN] and the measured velocity dispersion of the lensing galaxy [@Fal97] is $$H_0 = (67 \pm  13) \  {\rm km} \ {\rm sec}^{-1} \ {\rm Mpc}^{-1},$$ where the uncertainty comprises the 95% confidence level.

More time delays are becoming available for other lens systems (e.g., [@Big97; @Sch97a]). Blandford & Kundić [@BK97] provide a nice review in which they explore the potential to get a good determination of the extragalactic distance scale by combining measured time delays with good models; see also [@Sch97b] and [@WS97] for very recent summaries of the current situation on time delays and determination of the Hubble constant from lensing.

Quasar microlensing 
--------------------

Light bundles from “lensed" quasars are split by intervening galaxies. With typical separations of order one arcsecond between center of galaxy and quasar image this means that the quasar light bundle passes through the galaxy and/or the galaxy halo. Galaxies consist at least partly of stars, and galaxy haloes consist possibly of compact objects as well.

Each of these stars (or other compact objects, like black holes, brown dwarfs, or planets) acts as a “compact lens" or “microlens" and produces at least one new image of the source. In fact, the “macro-image" consists of many “micro-images" (Figure \[fig-micro-image\]). But because the image splitting is proportional to the lens mass – see equation (\[eq-angle\]) – these microimages are only of order a microarcsecond apart and can not be resolved. Various aspects of microlensing have been addressed after the first double quasar had been discovered [@CR79; @CR84; @Got81; @KRS; @Pac86a; @SW87; @Wam90].

(14.0,10.0) (2.5,-0.1) [=10.0cm ]{}

The surface mass density in front of a multiply imaged quasar is of order the “critical surface mass density", see equation (\[eq-crit\]). Hence microlensing should be occuring basically all the time. This can be visualized in the following way. If one assigns each microlens a little disk with radius equal to the Einstein ring, then the fraction of sky which is covered by these disks corresponds to the surface mass density in units of the critical density; this fraction is sometimes also called the “optical depth".

The microlenses produce a complicated two-dimensional magnification distribution in the source plane. It consists of many caustics, locations that correspond to formally infinitely high magnification. An example for such a magnification pattern is shown in Figure \[fig-magpat\]. It is determined with the parameters of image A of the quadruple quasar Q2237+0305 (surface mass density $\kappa = 0.36$; external shear $\gamma = 0.44$). The grey scale indicates the magnification: dark grey is relatively low magnification, and white is very high magnification.

Due to the relative motion between observer, lens and source the quasar changes its position relative to this arrangement of caustics, i.e. the apparent brightness of the quasar changes with time. A one-dimensional cut through such a magnification pattern, convolved with a source profile of the quasar, results in a microlensed lightcurve. Examples for microlensed lightcurves taken along the white tracks in Figure \[fig-magpat\] can be seen in Figure \[fig-lightcurve\] for two different quasar sizes.

In particular when the quasar track crosses a caustic (the sharp lines in Figure \[fig-magpat\] for which the magnification formally is infinite, because the determinant of the Jacobian disappears, cf. Equation (\[eq-magn-det\])), a pair of highly magnified microimages appears newly or merges and disappears (see [@Bla86]). Such a microlensing event can easily be detected as a strong peak in the lightcurve of the quasar image.

(14.0,12.0) (1.5, 0.1) [=12.0cm ]{}

(14.0,10.0) (2.5,-0.5) [=10.0cm ]{}

In most simulations it is assumed that the relative positions of the microlenses is fixed and the lightcurves are produced only by the bulk motion between quasar, galaxy and observer. A visualization of a situation with changing microlens positions for three different values of the surface mass density can be found in three video sequences accompanying [@WK95]. This change of caustics shapes due to the motion of individual stars produces additional fluctuations in the lightcurve [@KW93; @WK95].

Microlens-induced fluctuations in the observed brightness of quasars contain information both about the light-emitting source (size of continuum region or broad line region of the quasar, brightness profile of quasar) and about the lensing objects (masses, density, transverse velocity). Hence from a comparison between observed and simulated quasar microlensing (or lack of it) one can draw conclusions about the density and mass scale of the microlenses. It is not trivial, though, to extract this information quantitatively. The reason is that in this regime of optical depth of order one, the magnification is not due to a single isolated microlens, but it rather is a collective effect of many stars. This means individual mass determinations are not just impossible from the detection of a single caustic-crossing microlensing event, but it does not even make sense to try do so, since these events are not produced by individual lenses[^11]. Mass determinations can only be done in a statistical sense, by comparing good observations (frequently sampled, high photometric accuracy) with simulations. Interpreting microlensed lightcurves of multiple quasars allows to determine the size of the continuum emitting region of the quasar and to learn even more about the central engine [@Gou97; @Jar92; @RB91; @WPS90].

So far the “best" example of a microlensed quasar is the quadruple quasar Q2237+0305 [@H85; @Irw89; @Lew97; @Ost96; @WPS90; @Web91; @Wit93]. In Figure \[fig-2237-lewis\] two images of this system are shown which were taken in 1991 and 1994, respectively. Whereas on the earlier observation image B (top) is clearly the brightest, three years later image A (bottom) is at least comparable in brightness. Since the time delay in this system is only a day or shorter (because of the symmetric image arrangement), any brightness change on larger time scales must be due to microlensing. In Figure \[fig-2237-light\] lightcurves are shown for the four images of Q2237+0305 over a period of almost a decade (from [@Lew95]). The changes of the relative brightnesses of these images induced by microlensing are obvious.

(14.0,5.0) (2.0,-0.1) [=11.0cm ]{}

(14.0,8.0) (2.0,-0.1) [=12.0cm ]{}

Einstein rings  {#sec-einstein_rings}
---------------

If a point source lies exactly behind a point lens, a ring-like image occurs. Theorists had recognized early on [@Chw24; @Ein36] that such a symmetric lensing arrangement would result in a ring-image, a so-called “Einstein-ring". Can we observe Einstein rings? There are two necessary requirements for their occurence: the mass distribution of the lens needs to be axially symmetric, as seen from the observer, and the source must lie exactly on top of the resulting degenerate point-like caustic. Such a geometric arrangement is highly unlikely for point-like sources. But astrophysical sources in the real universe have a finite extent, and it is enough if a part of the source covers the point caustic (or the complete astroid caustic in a case of a not quite axial-symmetric mass distribution) in order to produce such an annular image.

In 1988 the first example of an “Einstein ring" was discovered [@Hew88]. With high resolution radio observations, the extended radio source MG1131+0456 turned out to be a ring with a diameter of about 1.75 arcsec. The source was identified as a radio lobe at a redshift of $z_S = 1.13$, whereas the lens is a galaxy at $z_L = 0.85$. Recently a remarkable observation of the Einstein ring 1938+666 was presented [@Kin97]. The infrared HST image shows an almost perfectly circular ring with two bright parts plus the bright central galaxy. The contours agree very well with the MERLIN radio map (see Figure \[fig-1938\]).

(14.0,10.0) (1.5, 0.1) [=12.0cm ]{}

By now about a half dozen cases have been found that qualify as Einstein rings. Their diameters vary between 0.33 and about 2 arcseconds. All of them are found in the radio regime, some have optical or infrared counterparts as well. Some of the Einstein rings are not really complete rings, but they are “broken" rings with one or two interruptions along the circle. The sources of most Einstein rings have both an extended and a compact component. The latter is always seen as a double image, separated by roughly the diameter of the Einstein ring. In some cases monitoring of the radio flux showed that the compact source is variable. This gives the opportunity to measure the time delay and the Hubble constant $H_0$ in these systems.

The Einstein ring systems provide some advantages over the multiply-imaged quasar systems for the goal to determine the lens structure and/or the Hubble constant. First of all the extended image structure provides many constraints on the lens. A lens model can be much better determined than in cases of just two or three or four point-like quasar images. Einstein rings thus help us to understand the mass distribution of galaxies at moderate redshifts. For the Einstein ring MG 1654+561 it was found [@Koc95] that the radially averaged surface mass density of the lens was fitted well with a distribution like $\Sigma(r) \propto r^{\alpha}$, where $\alpha$ lies between $-1.1 \le \alpha \le -0.9$ (an isothermal sphere would have exactly $\alpha = -1$!); there was also evidence found for dark matter in this lensing galaxy.

Second, since the diameters of the observed rings (or the separations of the accompanying double images) are of order one or two arcseconds, the expected time delay must be much shorter than the one in the double quasar Q0957+561 (in fact, it can be arbitrarily short, if the source happens to be very close to the point caustic). This means one does not have to wait so long to establish a time delay (but the source has to be variable intrinsically on even shorter time scales ...).

The third advantage is that since the emitting region of the radio flux is presumably much larger than that of the optical continuum flux, the radio lightcurves of the different images are not affected by microlensing. Hence the radio lightcurves between the images should agree with each other very well.

Another interesting application is the (non-)detection of a central image in the Einstein rings. For singular lenses, there should be no central image (the reason is the discontinuity of the deflection angle). However, many galaxy models predict a finite core in the mass distribution of a galaxy. The non-detection of the central images puts strong constraints on the size of the core radii.

Giant luminous arcs and arclets  {#sec-arcs}
--------------------------------

Zwicky had pointed out the potential use in the 1930s, but nobody had really followed up the idea, not even after the discovery of the lensed quasars: Galaxies can be gravitationally lensed as well. Since galaxies are extended objects, the apparent consequences for them would be far more dramatic than for quasars: galaxies should be heavily deformed once they are strongly lensed.

It came as quite a surprise when in 1986 Lynds & Petrosian [@LP86] and Soucail et al. [@Sou87] independently discovered this new gravitational lensing phenomenon: magnified, distorted and strongly elongated images of background galaxies which happen to lie behind foreground clusters of galaxies.

Rich clusters of galaxies at redshifts beyond $z \approx 0.2$ with masses of order $10^{14} M_{\odot}$ are very effective lenses if they are centrally concentrated. Their Einstein radii are of the order of 20 arcseconds. Since most clusters are not really spherical mass distributions and since the alignment between lens and source is usually not perfect, no complete Einstein rings have been found around clusters. But there are many examples known with spectacularly long arcs which are curved around the cluster center, with lengths up to about 20 arcseconds.

The giant arcs can be exploited in two ways, as is typical for many lens phenomena. Firstly they provide us with strongly magnified galaxies at (very) high redshifts. These galaxies would be too faint to be detected or analysed in their unlensed state. Hence with the lensing boost we can study these galaxies in their early evolutionary stages, possibly as infant or proto-galaxies, relatively shortly after the big bang. The other practical application of the arcs is to take them as tools to study the potential and mass distribution of the lensing galaxy cluster. In the simplest model of a spherically symmetric mass distribution for the cluster, giant arcs form very close to the critical curve, which marks the Einstein ring. So with the redshifts of the cluster and the arc it is easy to determine a rough estimate of the lensing mass by just determining the radius of curvature and interpreting it as the Einstein radius of the lens system.

More detailed modelling of the lensing clusters which allows for the asymmetry of the mass distribution according to the visible galaxies plus an unknown dark matter component provides more accurate determinations for the total cluster mass and its exact distribution. More than once this detailed modelling predicted additional (counter-) images of giant arcs, which later were found and confirmed spectroscopically [@Kne93; @Ebb97].

Gravitational lensing is the third method for the determination of masses of galaxy clusters, complementary to the mass determinations by X-ray analysis and the old art of using the virial theorem and the velocity distribution of the galaxies (the latter two methods use assumptions of hydrostatic or virial equilibrium, respectively). Although there are still some discrepancies between the three methods, it appears that in relaxed galaxy clusters the agreement between these different mass determinations is very good [@All97].

Some general results from the analysis of giant arcs in galaxy clusters are: Clusters of galaxies are dominated by dark matter. The typical “mass-to-light ratios" for clusters obtained from strong (and weak, see below) lensing analyses are   $M/L \ge 100 M_\odot/L_\odot$. The distribution of the dark matter follows roughly the distribution of the light in the galaxies, in particular in the central part of the cluster. The fact that we see such arcs shows that the central surface mass density in clusters must be high. The radii of curvature of many giant arcs is comparable to their distance to the cluster centers; this shows that core radii of clusters – the radii at which the mass profile of the cluster flattens towards the center – must be of order this distance or smaller. For stronger constraints detailed modelling of the mass distribution is required.

In Figures \[fig-2218\] and \[fig-0024\] two of the most spectacular cluster lenses producing arcs can be seen: Clusters Abell 2218 and CL0024+1654. Close inspection of the HST image of Abell 2218 reveals that the giant arcs are resolved (Figure \[fig-2218\]), structure can be seen in the individual components [@Kne96] and used for detailed mass models of the lensing cluster. In addition to the giant arcs, more than 100 smaller “arclets" can be identified in Abell 2218. They are farther away from the lens center and hence are not magnified and curved as much as the few giant arcs. These arclets are all slightly distorted images of background galaxies. With the cluster mass model it is possible to predict the redshift distribution of these galaxies. This has been successfully done in this system with the identification of an arc as a star-forming region, opening up a whole new branch for the application of cluster lenses [@Ebb96]. In another impressive exposure with the Hubble Space Telescope, the galaxy cluster CL0024+1654 (redshift $z = 0.39$) was deeply imaged in two filters [@CTT96]. The combined picture (Figure \[fig-0024\]) shows very nicely the reddish images of cluster galaxies, the brightest of them concentrated around the center, and the bluish arcs. There are four blue images which all have a shape reminiscent of the Greek letter $\Theta$. All the images are resolved and show similar structure (e.g., the bright fishhook-like feature at one end of the arcs), but two of them are mirror inverted, i.e. have different parity! They lie roughly on a circle around the center of the cluster and are tangentially elongated. There is also another faint blue image relatively close to the cluster center, which is extended radially. Modelling reveals that this is a five-image configuration produced by the massive galaxy cluster. All the five arcs are images of the same galaxy, which is far behind the cluster at a much higher redshift and most likely undergoes a burst of star formation. This is a spectacular example of the use of a galaxy cluster as a “Zwicky" telescope.

(14.0,8.0) (1.5, 0.1) [=12.0cm ]{}

(14.0,13.0) (1.5, 0.1) [=12.0cm ]{}

In CL0024+1654 the lensing effect produces a magnification of roughly a factor of ten. Combined with the angular resolution of the HST of 0.1 arcsec, this can be used to yield a resolution that effectively corresponds to 0.01 arcsec (in the tangential direction), unprecedented in direct optical imaging. Colley et al. [@CTT96] map the five images “backward" to the source plane with their model for the cluster lens and hence reconstruct the un-lensed source. They get basically identical source morphology for all arcs, which confirms that the arcs are all images of one source.

Recently, yet another superlative about cluster lenses was found: A new giant luminous arc was discovered in the field of the galaxy cluster CL1358+62 with the HST [@Fra97]. This arc-image turned out to be a galaxy at a redshift of $z = 4.92$. Up to a few months ago this was the most distant object in the universe with a spectroscopically measured redshift! In contrast to most other arcs, this one is very red. The reason is that due to this very high redshift, the Lyman-$\alpha$ emission of the galaxy, which is emitted in the ultra-violet part of the electromagnetic spectrum at a wavelength of 1216 Å is shifted by a factor of $z+1 \approx 6$ to the red part at a wavelength of about $7200$Å!

A review of cluster lensing and giant arcs/arclets can be found in Fort & Mellier [@FM94]. The review by Wu [@Wu96] provides, e.g., an updated table of giant arcs.

Weak/statistical lensing 
-------------------------

In contrast to the phenomena that were mentioned so far, “weak lensing" deals with effects of light deflection that cannot be measured individually, but rather in a statistical way only. As was discussed above “strong lensing" – usually defined as the regime that involves multiple images, high magnifications, and caustics in the source plane – is a rare phenomenon. Weak lensing on the other hand is much more common. In principle, weak lensing acts along each line of sight in the universe, since each photon’s path is affected by matter inhomogeneities along or near its path. It is just a matter of how accurate we can measure (cf. [@Pre95]).

Any non-uniform matter distribution between our observing point and distant light sources affects the measurable properties of the sources in two different ways: the angular size of extended objects is changed and the apparent brightness of a source is affected, as was first formulated in 1967 by Gunn [@Gun67a; @Gun67b].

A weak lensing effect can be a small deformation of the shape of a cosmic object, or a small modification of its brightness, or a small change of its position. In general the latter cannot be observed, since we have no way of knowing the unaffected position[^12].

The first two effects – slight shape deformation or small change in brightness – in general cannot be determined for an individual image. Only when averaging over a whole ensemble of images it is possible to measure the shape distortion, since the weak lensing (due to mass distributions of large angular size) acts as the coherent deformation of the shapes of extended background sources.

The effect on the apparent brightness of sources shows that weak lensing can be both a blessing and a curse for astronomers: the statistical incoherent lens-induced change of the apparent brightness of (widely separated) “standard candles" – like type Ia supernovae – affects the accuracy of the determination of cosmological parameters [@Fri97; @Kan95; @Wam97].

The idea to use the weak distortion and tangential alignment of faint background galaxies to map the mass distribution of foreground galaxies and clusters has been floating around for a long time. The first attempts go back to the years 1978/79, when Tyson and his group tried to measure the positions and orientations of the then newly discovered faint blue galaxies, which were suspected to be at large distances. Due to the not quite adequate techniques at the time (photographic plates), these efforts ended unsuccessfully [@VJMT84; @VTJ83]. Even with the advent of the new technology of CCD cameras, it was not immediately possible to detect weak lensing, since the pixel size originally was relatively large (of order an arcsecond). Only with smaller CCD pixels, improved seeing conditions at the telecope sites and improved image quality of the telescope optics the weak lensing effect could ultimately be measured.

Weak lensing is one of the two sub-disciplines within the field of gravitational lensing with the highest rate of growth in the last couple of years (along with galactic microlensing). There are a number of reasons for that:

1.  the availability of astronomical sites with very good seeing conditions,

2.  the availability of large high resolution cameras with fields of view of half a degree at the moment (aiming for more),

3.  the availability of methods to analyse these coherent small distortions

4.  the awareness of both observers and time allocation committees about the potential of these weak lensing analyses for extragalactic research and cosmology.

Now we will briefly summarize the technique of how to use the weak lensing distortion in order to get the mass distribution of the underlying matter.

### Cluster mass reconstruction {#cluster-mass-reconstruction .unnumbered}

The first real detection of a coherent weak lensing signal of distorted background galaxies was measured in 1990 around the galaxy clusters Abell 1689 and CL1409+52 [@TVW90]. It was shown that the orientation of background galaxies – the angle of the semi-major axes of the elliptical isophotes relative to the center of the cluster – was more likely to be tangentially oriented relative to the cluster than radially. For an unaffected population of background galaxies one would expect no preferential direction. This analysis is based on the assumption that the major axes of the background galaxies are intrinsically randomly oriented.

With the elegant and powerful method developed by Kaiser and Squires [@KS93] the weak lensing signal can be used to quantitatively reconstruct the surface mass distribution of the cluster. This method relies on the fact that the convergence $\kappa (\theta)$ and the two components of the shear $\gamma_1(\theta)$, $\gamma_2(\theta)$ are linear combinations of the second derivative of the effective lensing potential $\Psi(\theta)$ (cf. equations \[eq-kappa\] - \[eq-gamma2\]). After Fourier transforming the expressions for the convergence and the shear one obtains linear relations between the transformed components $\tilde \kappa$, $\tilde \gamma_1$, $\tilde \gamma_2$. Solving for $\tilde \kappa$ and inverse Fourier transforming gives an estimate for the convergence $\kappa$ (details can be found in [@KS93], [@Kai95b], [@NB96], or [@SK96]).

The original Kaiser-Squires method was improved/modified/extended/generalized by various authors subsequently. In particular the constraining fact that observational data are available only in a relatively small, finite area was implemented. Maximum likelihood techniques, non-linear reconstructions as well as methods using the amplification effect rather than the distortion effect complement each other. Various variants of the mass reconstruction technique have been successfully applied to more than a dozen rich clusters by now. Descriptions of various techniques and applications for the cluster mass reconstruction can be found in, e.g., [@Abd97; @Bar95; @Bar96a; @Bar96; @Bon95; @Bro95; @Ham97; @Kai95a; @Nat96; @SS95; @Sei96; @SS96; @Sma95; @Wil96]. In Figure \[fig-hoekstra\] a recent example for the reconstructed mass distribution of galaxy cluster CL1358+62 is shown [@Hoe97].

(14.0,12.3) (1.5, -0.1) [=12.0cm ]{}

We could present here only one weak lensing issue in some detail: the reconstruction of the mass distribution of galaxy clusters from weakly distorted background images. Many more interesting weak lensing applications are under theoretical and observational investigation, though. To name just a few:

-   constraints on the distribution of the faint galaxies from weak lensing (e.g., [@FM97; @LK97]);

-   galaxy-galaxy lensing (e.g., [@Bra96]);

-   lensing by galaxy halos in clusters (e.g., [@Nat97]);

-   weak lensing effect by large scale structure and/or detection of dark matter concentrations; this considers both shear effects as well as magnification effects (e.g. [@barkana; @Bar95a; @Sch96; @Wil97]);

-   determination of the power spectrum of the matter distribution (e.g., [@Ber97]; or

-   the weak lensing effects on the cosmic microwave background (e.g. [@Mar97], [@Met97b], [@Sel96]).

An upcoming comprehensive review on weak lensing by Schneider & Bartelmann [@SB98] treats both theory and applications of weak lensing in great depths.

Cosmological aspects of (strong) lensing  {#sec-cosmological}
-----------------------------------------

Gravitational lenses can be used in two different ways to study the cosmological parameters of the universe. The first is to explore a particular lens system in great detail, determine all possible observational parameters (image positions/brightnesses/shapes; matter/light distribution of lens; time variability etc.) and model both lens and source in as much detail as possible. This way one can in principle determine the amount of dark matter in the lens and – maybe even more importantly – the value of the Hubble constant. A reliable determination of the Hubble constant establishes the extragalactic distance scale, something astronomers have been trying to do for more than 70 years.

The second approach is of statistical nature: find out [*how many*]{} (what fraction of) quasars are multiply imaged by gravitationally lensing, determine their separation and redshift distributions [@TOG] and deduce the value of (or limits to) $\Omega_{\rm compact}$ – matter in clumps of, say, $10^6 \le M/M_{\odot} \le 10^{14}$ – and to $\Omega_{\Lambda}$ – the value of the cosmological constant.

The first approach has already been treated in Section \[sec-mult\]. Here we will concentrate on the statistical approach. In order to determine which fraction of a certain group of objects is affected by strong lensing (i.e. multiply imaged), one first needs a well-defined underlying sample. What needs to be done is the following:

1.  Do a systematic study of a sample of high-redshift objects: quasar surveys.

2.  Identify the gravitational lens systems among them.

3.  Determine the relative frequency of lensed objects, the distribution of splitting angles $\Delta \theta$ as a function of lens and source redshifts $z_L$/$z_S$.

4.  Determine matter content of universe $\Omega_{\rm compact}$, typical mass scale $M_{lens}$, cosmological constant $\Omega_\Lambda$, by comparison with theoretical models/simulations.

Since quasars are rare objects and lensing is a relatively rare phenomenon, steps 1 and 2 are quite difficult and time-consuming. Nevertheless, a number of systematic quasar surveys with the goal to find (many) lens systems with well defined selection criteria have been done in the past and others are underway right now (e.g., [@Bro97; @Mao93; @Mao97; @Web88; @Yee93]).

The largest survey so far, the CLASS survey, has looked at about 7000 radio sources at the moment (the goal is 10000). In total CLASS found 12 new lens systems so far. Interestingly, all the lenses have small separations ($\Delta \theta < 3$arcsec), and all lensing galaxies are detected [@Bro97; @Jac97a]. That leaves little space for a population of dark objects with masses of galaxies or beyond. A detailed discussion of lens surveys and a comparison between optical and radio surveys can be found in [@Koc96a].

The idea for the determination of the cosmological constant $\Omega_{\Lambda} = \Lambda/(3 H_0^2)$ from lens statistics is based on the fact that the relative lens probability for multiple imaging increases rapidly with increasing $\Omega_{\Lambda}$ (cf. Figure 9 of [@CPT92]). This was first pointed out 1990 [@FFK90; @Tur90]. The reason is the fact that the angular diameter distances $D_S$, $D_L$, $D_{LS}$ depend strongly on the cosmological model. And the properties that determine the probability for multiple lensing (i.e. the “fractional volume" that is affected by a certain lens) depend on these distances [@CPT92]. This can be seen, e.g., when one looks at the critical surface mass density required for multiple imaging (cf. equation \[eq-crit\]) which depends on the angular diameter distances.

The consequences of lensing studies on the cosmological constant can be summarized as follows. The analyses of the frequency of lensing are based on lens systems found in different optical and radio surveys. The main problem is still the small number of lenses. Depending on the exact selection criteria, only a few lens systems can be included in the analyses. Nevertheless, one can use the existing samples to put limits on the cosmological constant. Two different studies found 95%-confidence limits of $\Omega_{\Lambda} < 0.66$ [@Koc96b] and $\Omega_{\Lambda} < 0.7$ [@MR93; @Rix96]. This is based on the assumption of a flat universe ($\Omega_{matter} + \Omega_{\Lambda} = 1$). Investigations on the matter content of the universe from (both “macro-" and “micro-") lensing generally conclude that the fractional matter in compact form cannot exceed a few percent of the critical density (e.g. [@Can82; @Dal94; @Nem91; @Sch93]).

Galactic microlensing {#sec-gal-micro}
---------------------

It has been known for more than two decades that halos of galaxies must contain some unknown kind of dark matter. Many different particles/objects had been suggested as constituents of this halo dark matter. The candidates can be divided into the two broad categories “elementary particles" and “astronomical bodies". A conservative candidate for this dark matter are brown dwarfs, objects with masses less than 0.08 $M_{\odot}$ so that the central temperature is not high enough to start helium fusion. These objects are certain to exist, we just do not know how many there are.

In 1986 Paczyński [@Pac86b] suggested a method to test observationally whether the Milky Way halo is made of such brown dwarfs (or other astronomical objects in roughly this mass range). Subsequently this type of dark matter candidate was labelled “Macho" for MAssive Compact Halo Object [@Gri91]. If one could continuously observe the brightness of stars of our neighbouring galaxy Large Magellanic Cloud (LMC) one should see typical fluctuations in some of these stars due to the fact that every now and then one of these compact halo objects passes in front of the star and magnifies its brightness. The only problem with this experiment is the low probability for such an event: only about one out of three million LMC stars would be significantly magnified at any given time.

The underlying scenario is very simple: Due to the relative motion of observer, lensing Macho and source star the projected impact parameter between lens and source changes with time and produces a time dependent magnification. If the impact parameter is smaller than an Einstein radius then the magnification is $\mu_{min} >  1.34$ (cf. equation \[eq-magn-tot\]).

For an extended source such a sequence is illustrated in Figure \[fig-macho-1\]. The separation of the two images is of order two Einstein radii when they are of comparable magnification, which corresponds to only about a milliarcsecond. Hence the two images cannot be resolved individually, we can only observe the brightness of the combined image pair. This is illustrated in Figures \[fig-macho-2a\] and \[fig-macho-2b\] which show the relative tracks and the respective light curves for five values of the minimum impact parameter $u_{min}$.

(14.0,5.0) (1.5,-0.1) [=12.0cm ]{}

(14.0,9.0) (0.5,-2.2) [=12.0cm ]{}

(14.0,11.0) (1.5,-0.5) [=12.0cm ]{}

Quantitatively, the total magnification $\mu = \mu_1 + \mu_2$ of the two images (cf. equation \[eq-magn-tot\]) entirely depends on the impact parameter $u(t) = r(t)/R_E$ between the lensed star and the lensing object, measured in the lens plane (here $R_E$ is the Einstein radius of the lens, i.e. the radius at which a circular image appears for perfect alignment between source, lens and observer, cf. Figure \[fig-macho-1\], rightmost panel): $$\label{eq-light}
\mu(u) = { {u^2 + 2 } \over u \sqrt{ (u^2 + 4) } }.$$ The time scale of such a “microlensing event" is defined as the time it takes the source to cross the Einstein radius: $$\label{eq-macho-time}
t_0    = {R_E \over v_{\perp} } \ \approx \ \ 0.214 \ \  {\rm yr} \ \ \ 
\sqrt{M   \over M_\odot} \ 
\sqrt{D_L   \over 10 {\rm kpc}} \ 
\sqrt{1 - {D_L   \over D_S}            }  \ 
\left({v_\perp   \over 200 {\rm km/sec}}\right)^{-1}.$$ Here $v_{\perp}$ is the (relative) transverse velocity of the lens. We parameterized the time scale by “typical" numbers for the distances of lensed and lensing star and the relative transverse velocity. Note also that here we used the simple relation $D_{LS} = D_S - D_L$ (which is not valid for cosmological distances).

Note that from equation (\[eq-macho-time\]) it is obvious that it is [**not**]{} possible to determine the mass of the lens from one individual microlensing event. The duration of an event is determined by three unknown parameters: the mass of the lens, the transverse velocity and the distances of lens and source. It is impossible to disentangle these for individual events. Only with a model for the spatial and velocity distribution of the lensing objects and comparison with “simulated microlensing events" it is possible to obtain information about the masses of the lensing objects and their density.

What seemed to be an impossible task at the time – namely determine the brightness of millions of stars on an almost nightly basis – was turned into three big observational campaigns within few years (MACHO, EROS, OGLE experiments). These groups looked at millions of stars in the LMC and towards the bulge of the Milky Way, and their first results appeared essentially simultaneously in the fall of 1993 [@macho93; @eros93; @ogle93]. In the meantime more groups have joined this effort, some of them with special emphases: e.g. on covering ongoing microlensing events (PLANET, DUO), or on extending the microlensing search to unresolved stars (“pixel lensing") in the Andromeda galaxy [@Cro96; @Gou96] (AGAPE) or to cover the Magellanic Clouds completely around the year (MOA). Here is a list of groups currently active in the search for microlensing signatures of compact objects in the halo of the Milky Way or elsewhere:

-   MACHO (MAssive Compact Halo Object): [@macho97a; @macho97b]

-   EROS (Experience de Recherche d’Objets Sombres): [@eros96; @eros97]

-   OGLE (Optical Gravitational Lens Experiment): [@ogle94a; @ogle97]

-   AGAPE (Andromeda Galaxy and Amplified Pixels Experiment): [@agape97]

-   MOA (MACHO Observations in Astrophysics): [@moa96]

-   PLANET (Probing Lensing Anomalies NETwork): [@planet96]

-   DUO (Disk Unseen Objects) [@Ala97]

-   GMAN (Global Microlensing Alert Network) [@gman97]

The observations towards the Large Magellanic Cloud show that there are fewer microlensing events than one would expect if the halo of the Milky Way was made entirely of these compact objects. The latest published results from the microlensing experiments that monitor stars in the LMC indicate that the optical depths toward the LMC is about $\tau \approx 3 \times 10^{-7}$. The observations are consistent with 50% of the Milky Way halo made of compact objects with most likely masses of $0.5^{+0.3}_{-0.2}  M_{\odot}$ [@macho97a]. But the number of observed events is still small (in this analysis eight events were used) and hence the uncertainties are large; in fact, it cannot even be excluded that none of the observed events is due to an unknown halo population [@Gat97].

The same type of experiment (searching for microlensing events) is being performed in the direction of the galactic bulge as well, the central part of the Milky Way. By now more than 200 microlensing events have been detected in this direction (for an example see Figure \[fig-bulge-1\]). Among them are a number of “binary lens"-events (which have a very typical signature of at least two caustic crossings, cf. Figure \[fig-bulge-2\]). This is about three times as many microlensing events as were expected/predicted. Several groups try to explain this “over-abundance" of events to have a new look at the stellar content and the dynamics of the bar/bulge of the Galaxy. The latest published results can be found in [@macho97b].

With these microlensing experiments gravitational lensing has established itself as a new tool to study the structure of the Milky Way. This type of microlensing also holds some promise for the future. It can be used, e.g. to study the frequency of binary stars. One of the most interesting possibilities is to detect planets around other stars by extending the sensitivity of the binary lenses to smaller and smaller companion masses [@MP91; @Wam97pla].

(14.0,8.0) (0.0,-0.3) [=15.0cm ]{}

(14.0,8.5) (0.5,-0.1) [=15.0cm ]{}

For a recent comprehensive presentation of galactic microlensing and beyond see [@Pac96]. Various aspects of microlensing in the local group are reviewed in detail. Another review article on the basics and the results of galactic microlensing can be found in [@RM96].

Future Gravitational Lensing {#sec-futu}
============================

Gravitational lensing is an exceptional field in astronomy in the sense that its occurence and many of its features – e.g. multiple images, time delays, Einstein rings, quasar microlensing, galactic microlensing, weak lensing – were [*predicted*]{} (long) before they were actually observed. Although “prediction" or predictability is considered one of the important criteria of modern science, many (astro-)physical phenomena are too complicated for a minute prediction (just think of the weather forecast). The reason why this worked here is that gravitational lensing is a simple geometrical concept which easily allows qualitative estimates and quantitative calculations. Extrapolating from these thoughts, it should be possible to look forward in time once again and predict future applications of gravitational lensing.

However, at any given time it requires very good intuition, some courage and maybe even a bit of ingenuity to predict qualitatively new phenomena. It does not need much of either to envision that the known lensing phenomena will become better, sharper, more. My predictions for the next decade in this sense are humble and modest:

No doubt there will soon be more determinations of accurate time delays in multiply-imaged quasar systems. If the models will get more precise as well, the value of the Hubble constant $H_0$ determined from a number of lens systems will be accurate to a few percent or better and will probably turn out to be as reliable as $H_0$ values obtained with any other method [@WS97].

The frequencies, image separations, redshift distributions of multiply-imaged quasars and their lenses will become a major tool in differentiating between different cosmological models. The Sloan Digital Sky Survey, e.g., will discover a few hundred new lensed quasars with very well defined selection criteria, ideally suited for that purpose. Another angle on the cosmological model and the values of $\Omega$ and $\Lambda$ offer the statistics of arcs. The number of high redshift galaxies seen as arcs depends crucially on the number of rich galaxy clusters at intermediate redshifts. And since different cosmological models predict very different formation redshifts for clusters, this promising road should be followed as well [@Bar98].

The new facilities which become available now or in the near future in the infrared/sub-mm/mm domain – like SCUBA, SIRTF, FIRST, IRAM – will open a completely new window in these wavelength ranges, with supposedly most spectacular results in the arcs and cluster lensing regime.

Quasar microlensing will provide information on the structure of the quasars and the intervening clumped matter. With the new X-ray telescope AXAF with its high spatial resolution it will become possible to obtain X-ray lightcurves which due to the presumably smaller emission region will produce dramatic microlensing events in multiply-imaged quasars. Maybe we can “map" the hot spots of quasars this way.

The largest number of lensing events in the near future will doubtlessly come from the “local" microlensing experiments monitoring galactic bulge stars. The art of taking spectra of highly magnified stars during microlensing events (as pioneered by [@Len96]) will open up the fascinating possibility to investigate the metallicity of bulge stars in detail or even resolve the stellar surfaces and study their center-to-limb variations. In addition of being an excellent tool to study the structure of the Milky Way, galactic microlensing will also provide unbiased statistics on the fraction of binary stars (within certain relative distances). Extending the sensitivity to higher mass ratios between the binary components will naturally lead to the detection of planets around stars (at distances of many kiloparsecs!). Microlensing has the advantage compared to all other Earth-bound planet search techniques that it is able to detect [*Earth-mass planets*]{}! It is also imaginable that before too long such microlensing events could be detected directly by monitoring astrometrically the position of the star very accurately [@miralda].

In due course we should also know quantitatively how much dark compact objects contribute to the mass of the halo of the Milky Way, and what their mass range is. The “pixel lensing" will probe other lines of sight through the Galactic halo by exploring the Andromeda galaxy and other nearby galaxies. This will provide information on the three-dimensional mass distribution of the halo.

Weak lensing will be used to map not just the outskirts of massive galaxy clusters, but also to trace the large scale structure by its effect on the background population of galaxies. If we find good ways to discriminate between source galaxies at various redshifts, this way we can ultimately produce a three-dimensional map of the [*matter*]{} in the universe (rather than a light map)! This will be an utmost useful thing for the understanding of structure formation and evolution; as an aside we will determine the matter content of the universe $\Omega$.

Some other possible applications of lensing will be: The black hole in the Galactic center affects all sources that are near or behind the center. Mapping this effect will be a complementary determination of the black hole mass and will help to study the dynamics near the black hole. The redshift of the most distant object will be pushed beyond $z = 6$, and it is quite likely that it will be magnified by lensing. The next generation of experiments to map the cosmic microwave background will be sensititive enough to detect the gravitational lens signature of the matter “in front".

What about the not-so-predictable or not-so-easily-predictable future of lensing? Ultimately every object in the sky is affected by (ever so slight) lensing effects: this is the not-yet-reached regime of ultra weak lensing. I would like to conclude citing two remarks that Bill Press presented in his lensing outlook at the IAU Symposium 173 in Melbourne (1995). He mentions that “gravitational lens effects ... are present along virtually every line of sight" [@Pre95]. In a not quite so serious extrapolation Press points out that more and more astronomers will (have to) deal with lensing in the next decade, so that lensing will become an “ubiquitous observational technique" and hence – for better or for worse: “gravitational lensing may well disappear as a unique sub-specialty in astronomy".

Acknowledgements {#acknowledgements .unnumbered}
================

It is a pleasure to thank Wes Colley, Frederic Courbin, Emilio Falco, Henk Hoekstra, Neal Jackson, Tomislav Kundić, Geraint Lewis, and Andrzej Udalski for permission to use their figures. I would also like to thank Matthias Bartelmann, Emilio Falco, Jean-Paul Kneib, Bohdan Paczyński, Sjur Refsdal, Robert Schmidt, Liliya Williams and David Woods for their careful reading of (parts of) the manuscript at various stages and their useful comments. Of particular help were the comments of Jürgen Ehlers and an unknown referee which improved the paper considerably.

\#1 [\#1]{} \#1 [ ]{} å\#1,\#2.[[*Astron. Astroph.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Acta Astron.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Astron. Ap. Suppl.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Astron. Ap. Rev.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Astron. J.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Astron. and Geophys.J.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Astron. Nachr.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Ann. Rev. Astr. Ap.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Astrophys. J.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Astrophys. J. Lett.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Ap. J. Suppl.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Astrophys. and Space Sci.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Bull. Am. Astron. Soc.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Comments Astrophys.* ]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Fundamentals Cosm. Phys.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*J. Astrophys. Astr.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*J. Geophys. Res.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*J. Math. Phys.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Liv. Rev. Rel.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Mon. Not. R. Astron. Soc.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Nature*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Observatory*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Publ. Astr. Soc. Pac.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Phys. Rep.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Phil. Mag.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Phys. Rev.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Phys. Rev. B*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Phys. Rev. D*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Phys. Rev. Lett.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Proc. Roy. Soc. London*]{}, [**A\#1**]{}, \#2.]{} \#1,\#2.[[*Q. Jl R. astr. Soc.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Rev. Mod. Phys.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Rep. Prog. Phys.* ]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Sov. Astr.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Sov. Phys. Doklady.*]{}, [**\#1**]{}, \#2.]{} \#1,\#2.[[*Science*]{}, [**\#1**]{}, \#2.]{}

[99]{}

[^1]: “On the Deflection of a Light Ray from its Straight Motion due to the Attraction of a World Body which it Passes Closely"

[^2]: “On the Influence of Gravity on the Propagation of Light"

[^3]: “It would be very desirable that astronomers address the question unrolled here, even if the considerations should seem to be insufficiently founded or entirely speculative."

[^4]: As stated above, only very recently it was shown that Einstein had derived these equations as early as 1912, but did not bother to publish them [@Ren97].

[^5]: A detailed description of optics in curved spacetimes and a derivation of the lens equation from Einstein’s field equations can be found in Chapters 3 and 4 of [@SEF].

[^6]: In cosmology the various methods to define distance diverge (see, e.g. Chapter 14.4 of [@Wei72]) or 3.5 of [@SEF]). The relevant distances for gravitational lensing are the angular diameter distances, see [@NB96].

[^7]: Due to the fact that physical objects have a finite size, and also because at some limit wave optics has to be applied, in reality the magnification stays finite.

[^8]: This does not violate energy conservation, since this is the magnification relative to an “empty" universe and not relative to a “smoothed out" universe. This issue is treated in detail in, e.g., [@Sch84] or in Chapter 4.5 of [@SEF].

[^9]: There exists a theorem that gravitational lenses should produce an odd number of images ( e.g., [@McK85]).

[^10]: This can be seen very simply: Imagine a lens situation like the one displayed in Figure \[fig-setup\]. If now all length scales are reduced by a factor of two and at the same time all masses are reduced by a factor of two, then for an observer the angular configuration in the sky would appear exactly identical. But the total length of the light path is reduced by a factor of two. Now, since the time delay between the two paths is the same fraction of the total lengths in either scenario, a measurement of this [**fractional**]{} length allows to determine the total length, and hence the Hubble constant, the constant of proportionality between distance and redshift.

[^11]: Similarly, one cannot determine the temperature of a black body by measuring the energy of a single photon emitted by the black body, but one needs to measure a large number of them and compare with some underlying theory.

[^12]: A well-known exception is the light deflection at the solar limb, where the difference between the lensed and the unlensed positions of stars was used to confirm General Relativity, see Chapter \[sec-hist\].
