User:TallJimbo/WeakLensing/Full

While the presence of any mass bends the path of light passing near it, this effect rarely produces the giant arcs and multiple images associated with strong gravitational lensing. Most lines of sight in the universe are thoroughly in the weak lensing regime, in which the deflection is impossible to detect in a single background source. However, even in these cases, the presence of the foreground mass can be detected in a systematic alignment of background sources around the lensing mass. Weak gravitational lensing is thus an intrinsically statistical measurement, but it provides a way to measure the masses of astronomical objects without requiring assumptions about their composition or dynamical state.

Model

Gravitational lensing acts as a coordinate transformation that distorts the images of background objects (usually galaxies) near a foreground mass. The transformation can be split into two terms, the convergence and shear. The convergence term magnifies the background objects by increasing their size, and the shear term stretches them tangentially around the foreground mass.

To measure this tangential alignment, it is necessary to measure the ellipticities of the background galaxies and construct a statistical estimate of their systematic alignment. The fundamental problem is that galaxies are not intrinsically circular, so their measured ellipticity is a combination of their intrinsic ellipticity and the gravitational lensing shear. Typically, the intrinsic ellipticity is much greater than the shear (by a factor of 3-300, depending on the foreground mass). The orientation of intrinsic ellipticities of galaxies should be almost entirely random, so any systematic alignment between multiple galaxies can generally be assumed to be caused by lensing. However, the measurements of many background galaxies must be combined to average down this "shape noise."

Distortions of the type produced by lensing, acting on circles and a distribution of ellipses similar to that of real galaxies. Note that the distortion shown here is greatly exaggerated relative to real astronomical systems.

Another major challenge for weak lensing is correction for the point spread function due to instrumental and atmospheric effects, which causes the observed images to be smeared relative to the "true sky". This smearing tends to make small objects more round, destroying some of the information about their true ellipticity. As a further complication, the PSF typically adds a small level of ellipticity to objects in the image, which is not at all random, and can in fact mimic a true lensing signal. Even for the most modern telescopes, this effect is usually at least the same order of magnitude as the gravitational lensing shear, and is often much larger. Correcting for the PSF requires building a model for how it varies across the field. Stars in our own galaxy provide a direct measurement of the PSF, and these can be used to construct such a model, usually by interpolating between the points where stars appear on the image. This model can then be used to reconstruct the "true" ellipticities from the smeared ones. Ground-based and space-based data typically undergo distinct reduction procedures due to the differences in instruments and observing conditions.

Angular diameter distances to the lenses and background sources are important for converting the lensing observables to physically meaningful quantities. These distances are often estimated using photometric redshifts when spectroscopic redshifts are unavailable. Redshift information is also important to separate the background source population from other galaxies in the foreground or associated with the mass responsible for the lensing. With no redshift information, the foreground and background populations can be split by an apparent magnitude or a color cut, but this is much less accurate and can only benefit analyses of lensing by clusters or galaxies.

Weak lensing by clusters of galaxies

Galaxy clusters are among the largest gravitationally bound structures in the Universe, surpassed only by superclusters; approximately 80% of cluster content is in the form of dark matter; their gravitational fields deflect light-rays traveling near the cluster. As seen from Earth, this effect can cause dramatic distortions of a background source object detectable by eye such as multiple images, arcs, and rings (cluster strong lensing). More generally, the effect causes small, but statistically coherent, distortions of background sources on the order of 10% (cluster weak lensing). Abell 1689, CL0024+17, and the Bullet Cluster are among the most prominent examples of lensing clusters.

History

The effects of cluster strong lensing were first detected by Roger Lynds of the National Optical Astronomy Observatories and Vahe Petrosian of Stanford University who discovered giant luminous arcs in a survey of galaxy clusters in the late 1970s. Lynds and Petrosian published their findings in 1986 without knowing the origin of the arcs. In 1987, Genevieve Soucail of the Toulouse Observatory and her collaborators presented data of a blue ring-like structure in Abell 370 and proposed a gravitational lensing interpretation. The first cluster weak lensing analysis was conducted in 1990 by J. Anthony Tyson of Bell Laboratories and collaborators. Tyson et al. detected a coherent alignment of the ellipticities of the faint blue galaxies behind both Abell 1689 and CL 1409+52. Lensing has been used as a tool to investigate a tiny fraction of the thousands of known galaxy clusters.

Historically, lensing analyses were conducted on galaxy clusters detected via their baryon content (e.g. from optical or X-ray surveys). The sample of galaxy clusters studied with lensing was thus subject to various selection effects; for example, the most luminous clusters were investigated. In 2006, David Wittman of the University of California at Davis and collaborators published the first sample of galaxy clusters detected via their lensing signals, completely independent of their baryon content. However, clusters discovered through lensing are subject to mass selection effects because the more massive clusters produce lensing signals with higher signal-to-noise.

Observational products

The projected mass density can be recovered from the measurement of the ellipticities of the lensed background galaxies through techniques that can be classifed into two types: direct reconstruction and inversion. However, a mass distribution reconstructed without knowledge of the magnification suffers from a limitation known as the mass sheet degeneracy, where the cluster surface mass density κ can be determined only up to a transformation $\kappa \rightarrow \kappa ^{\prime }=\lambda \kappa +(1-\lambda )$ where λ is an arbitrary constant. This degeneracy can be broken if an independent measurement of the magnification is available because the magnification is not invariant under the aforementioned degeneracy transformation.

Given a centroid for the cluster, which can determined by using a reconstructed mass distribution or optical or X-ray data, a parametric model such as a singular isothermal sphere (SIS) profile or a Navarro-Frenk-White (NFW) profile can be fit to the shear profile as a function of clustrocentric radius. Knowledge of the lensing cluster redshift and the redshift distribution of the background galaxies is also necessary for estimation of the mass and size from a model fit; these redshifts can be measured precisely using spectroscopy or estimated using photometry. Individual mass estimates from weak lensing can only be derived for the most massive clusters, and the accuracy of these mass estimates are limited by projections along the line of sight.

Scientific implications

Photo of the Bullet Cluster from the Hubble Space Telescope with total mass contours (dominated by dark matter) from a lensing analysis overlaid.

Cluster mass estimates determined by lensing are valuable because the method requires no assumption about the dynamical state or star formation history of the cluster in question. Lensing mass maps can also potentially reveal "dark clusters," clusters containing overdense concentrations of dark matter, but relatively insignificant amounts of baryonic matter. Comparison of the dark matter distribution mapped using lensing with the distribution of the baryons using optical and X-ray data reveals the interplay of the dark matter with the stellar and gas components. A notable example of such a joint analysis is the so-called Bullet Cluster. The Bullet Cluster data provide constraints on models relating light, gas, and dark matter distributions such as Modified Newtonian dynamics (MOND) and Λ-Cold Dark Matter (Λ-CDM).

In principle, since the number density of clusters as a function of mass and redshift is sensitive to the underlying cosmology, cluster counts derived from large weak lensing surveys should be able to investigate cosmological parameters. In practice, however, projections along the line of sight cause many false positives. Weak lensing can also be used to calibrate the mass-observable relation via a stacked weak lensing signal around an ensemble of clusters, although this relation is expected to have an intrinsic scatter. In order for lensing clusters to be a precision probe of cosmology in the future, the projection effects and the scatter in the lensing mass-observable relation need to be thoroughly characterized and modeled.

Galaxy-galaxy lensing

Galaxy-galaxy lensing is a specific type of weak (and occasionally strong) gravitational lensing, in which the foreground object responsible for distorting the shapes of background galaxies is itself an individual field galaxy (as opposed to a galaxy cluster or the large-scale structure of the cosmos). Of the three typical mass regimes in weak lensing, galaxy-galaxy lensing produces a “mid-range” signal (shear correlations of ~1%) that is weaker than the signal due to cluster lensing, but stronger than the signal due to cosmic shear.

History

J.A. Tyson and collaborators first postulated the concept of galaxy-galaxy lensing in 1984, though the observational results of their study were inconclusive. It was not until 1996 that evidence of such distortion was tentatively discovered, with the first statistically significant results not published until the year 2000. Since those initial discoveries, the construction of larger, high resolution telescopes and the advent of dedicated wide field galaxy surveys have greatly increased the observed number density of both background source and foreground lens galaxies, allowing for a much more robust statistical sample of galaxies, making the lensing signal much easier to detect. Today, measuring the shear signal due to galaxy-galaxy lensing is a widely used technique in observational astronomy and cosmology, often used in parallel with other measurements in determining physical characteristics of foreground galaxies.

Stacking

Much like in cluster-scale weak lensing, detection of a galaxy-galaxy shear signal requires one to measure the shapes of background source galaxies, and then look for statistical shape correlations (specifically, source galaxy shapes should be aligned tangentially, relative to the lens center.) In principle, this signal could be measured around any individual foreground lens. In practice, however, due to the relatively low mass of field lenses and the inherent randomness in intrinsic shape of background sources (the “shape noise”), the signal is impossible to measure on a galaxy by galaxy basis. However, by combining the signals of many individual lens measurements together (a technique known as “stacking”), the signal to noise ratio will improve, allowing one to determine a statistically significant signal, averaged over the entire lens set.

Scientific applications

Galaxy-galaxy lensing (like all other types of gravitational lensing) is used to measure several quantities pertaining to mass:

• Mass density profiles
  • Using techniques similar to those in cluster-scale lensing, galaxy-galaxy lensing can provide information about the shape of mass density profiles, though these profiles correspond to galaxy-sized objects, instead of larger clusters or groups. Given a high enough number density of background sources, a typical galaxy-galaxy mass density profile can cover a wide range of distances (from ~1 to ~100 effective radii). Since the effects of lensing are insensitive to the matter type, a galaxy-galaxy mass density profile can be used to probe a wide range of matter environments: from the central cores of galaxies where baryons dominate the total mass fraction, to the outer halos where dark matter is more prevalent.

• Mass-to-light ratios
  • Comparing the measured mass to the luminosity (averaged over the entire galaxy stack) in a specific filter, galaxy-galaxy lensing can also provide insight into the mass to light ratios of field galaxies. Specifically the quantity measured through lensing is the total (or virial) mass to light ratio – again due to the insensitivity of lensing to matter type. Assuming that luminous matter can trace dark matter, this quantity is of particular importance, since measuring the ratio of luminous (baryonic) matter to total matter can provide information regarding the overall ratio of baryonic to dark matter in the universe.

• Galaxy mass evolution
  • Since the speed of light is finite, an observer on the Earth will see distant galaxies not as they look today, but rather as they appeared at some earlier time. By restricting the lens sample of a galaxy-galaxy lensing study to lie at only one particular redshift, it is possible to understand the mass properties of the field galaxies that existed during this earlier time. Comparing the results of several such redshift-restricted lensing studies (with each study encompassing a different redshift), one can begin to observe changes in the mass features of galaxies over a period of several epochs, leading towards a better understanding of the evolution of mass on the smallest cosmological scales.

• Other mass trends
  • Lens redshift is not the only quantity of interest that can be varied when studying mass differences between galaxy populations, and often there are several parameters used when segregating objects into galaxy-galaxy lens stacks. Two widely used criteria are galaxy color and morphology, which act as tracers of (among other things) stellar population, galaxy age, and local mass environment. By separating lens galaxies based on these properties, and then further segregating samples based on redshift, it is possible to use galaxy-galaxy lensing to see how several different types of galaxies evolve through time.

Cosmic shear

The gravitational lensing by large scale structure also produces an observable pattern of alignments in background galaxies, but this distortion is only ~0.1%-1% - much more subtle than cluster or galaxy-galaxy lensing. The thin lens approximation usually used in cluster and galaxy lensing does not always work in this regime, because structures can be elongated along the line of sight. Instead, the distortion can be derived by assuming that the deflection angle is always small (see gravitational lensing formalism). As in the thin lens case, the effect can be written as a mapping from the unlensed angular position ${\vec {\beta }}$ to the lensed position ${\vec {\theta }}$. The Jacobian of the transform can be written as an integral over the gravitational potential $\Phi$ along the line of sight

${\frac {\partial \beta _{i}}{\partial \theta _{j}}}=\delta _{ij}+\int _{0}^{r_{\infty }}drg(r){\frac {\partial ^{2}\Phi ({\vec {x}}(r))}{\partial x^{i}\partial x^{j}}}$

where $r$ is the comoving distance, $x^{i}$ are the transverse distances, and

$g(r)=2r\int _{r}^{r_{\infty }}\left(1-{\frac {r^{\prime }}{r}}\right)W(r^{\prime })$

is the lensing kernel, which defines the efficiency of lensing for a distribution of sources $W(r)$.

As in the thin-lens approximation, the Jacobian can be decomposed into shear and convergence terms.

Shear correlation functions

Because large scale cosmological structures do not have a well-defined location, detecting cosmological gravitational lensing typically involves the computation of shear correlation functions, which measure the mean product of the shear at two points as a function of the distance between those points. Because there are two components of shear, three different correlation functions can be defined:

$\xi _{++}(\Delta \theta )=\langle \gamma _{+}({\vec {\theta }})\gamma _{+}({\vec {\theta }}+{\vec {\Delta \theta }})\rangle$

$\xi _{\times \times }(\Delta \theta )=\langle \gamma _{\times }({\vec {\theta }})\gamma _{\times }({\vec {\theta }}+{\vec {\Delta \theta }})\rangle$

$\xi _{\times +}(\Delta \theta )=\xi _{+\times }(\Delta \theta )=\langle \gamma _{+}({\vec {\theta }})\gamma _{\times }({\vec {\theta }}+{\vec {\Delta \theta }})\rangle$

where $\gamma _{+~}$ is the component along or perpendicular to ${\vec {\Delta \theta }}$, and $\gamma _{\times }$ is the component at 45°. These correlation functions are typically computed by averaging over many pairs of galaxies. The last correlation function, $\xi _{\times +}$, is not affected at all by lensing, so measuring a value for this function that is inconsistent with zero is often interpreted as a sign of systematic error.

The functions $\xi _{++~}$ and $\xi _{\times \times }$ can be related to projections (integrals with certain weight functions) of the dark matter density correlation function, which can be predicted from theory for a cosmological model through it's Fourier transform, the matter power spectrum.

Because they both depend on a single scalar density field, $\xi _{++~}$ and $\xi _{\times \times }$ are not independent, and they can be decomposed further into E-mode and B-mode correlation functions. In analogy with electric and magnetic fields, the E-mode field is curl-free and the B-mode field is divergence-free. Because gravitational lensing can only produce an E-mode field, the B-mode provides yet another test for systematic errors.

The E-mode correlation function is also known as the aperture mass variance

$\langle M_{ap}^{2}\rangle (\theta )=\int _{0}^{2\theta }{\frac {\phi d\phi }{\theta ^{2}}}\left[\xi _{++}(\phi )+\xi _{\times \times }(\phi )\right]T_{+}\left({\frac {\phi }{\theta }}\right)=\int _{0}^{2\theta }{\frac {\phi d\phi }{\theta ^{2}}}\left[\xi _{++}(\phi )-\xi _{\times \times }(\phi )\right]T_{-}\left({\frac {\phi }{\theta }}\right)$

$T_{+}(x)=576\int _{0}^{\infty }{\frac {dt}{t^{3}}}J_{0}(xt)[J_{4}(t)]^{2}$

$T_{-}(x)=576\int _{0}^{\infty }{\frac {dt}{t^{3}}}J_{4}(xt)[J_{4}(t)]^{2}$

An exact decomposition thus requires knowledge of the shear correlation functions at zero separation, but an approximate decomposition is fairly insensitive to these values because the filters $T_{+}$ and $T_{-}$ are small near $\theta =0$.

Weak lensing and cosmology

The ability of weak lensing to constrain the matter power spectrum makes it a potentially powerful probe of cosmological parameters, especially when combined with other observations such as the cosmic microwave background, supernovae, and galaxy surveys. Detecting he extremely faint cosmic shear signal requires averaging over many background galaxies, so surveys must be both deep and wide, and because these background galaxies are small, the image quality must be very good. Measuring the shear correlations at small scales also requires a high density of background objects (again requiring deep, high quality data), while measurements at large scales push for wider surveys.

Because of these challenges, while weak lensing of large scale structure was discussed as early as 1967, it was not detected until more than 30 years later, when large CCD cameras enabled surveys of the necessary size and quality. In 2000, four independent groups published the first detections of cosmic shear, and subsequent observations have started to put constraints on cosmological parameters (particularly the dark matter density $\Omega _{m}~$ and power spectrum amplitude $\sigma _{8}~$) that are competitive with other cosmological probes.

For current and future surveys, one goal is to use the redshifts of the background galaxies (often approximated using photometric redshifts) to divide the survey into multiple redshift bins. The low-redshift bins will only be lensed by structures very near to us, while the high-redshift bins will be lensed by structures over a wide range of redshift. This technique, dubbed "cosmic tomography", makes it possible to map out the 3D distribution of mass. Because the third dimension involves not only distance but cosmic time, tomographic weak lensing is sensitive not only to the matter power spectrum today, but also to its evolution over the history of the universe, and the expansion history of the universe during that time. This is a much more valuable cosmological probe, and many proposed experiments to measure the properties of Dark Energy and Dark Matter have focused on weak lensing, such as the Dark Energy Survey, Pan-STARRS, and LSST.