Mapping Dark Matter: Gravitational Lensing

Learning Goal: Analyzing the physics of strong and weak gravitational lensing to map the distribution of dark matter in galaxy clusters and individual galaxies.

• Prerequisites: Multivariable calculus, introductory classical mechanics (Newtonian gravity), and basic electromagnetism (wave propagation).
• Estimated Total Study Time: 20 Hours

---

Module 1: Foundations of General Relativity and Spacetime

Overview

To understand why light bends across the universe, we must abandon the Newtonian concept of gravity as an instantaneous, straight-line force. This module introduces Albert Einstein's General Theory of Relativity, focusing on how mass and energy warp the four-dimensional fabric of spacetime, forcing light paths to follow curved trajectories called null geodesics.

Recommended Videos

Why this video

Delivered by theoretical physicist Leonard Susskind, this classic Stanford lecture uses the equivalence principle to show why gravity bends light. By comparing a stationary frame in a gravitational field to an accelerating elevator, Susskind walks through the kinematic derivation of light bending, showing that the path of a horizontal light ray must curve downward in an accelerating frame.

---

Why this video

This video provides an intuitive visual transition from Newton's gravitational model to Einstein's geometric model. It explains how massive bodies deform spacetime coordinates, visually demonstrating why objects (including photons) traveling in straight paths actually follow curved paths (geodesics) through warped space.

---

Why this video

In this conversational yet deep explanation, theoretical physicist Sean Carroll explains the core concept of General Relativity: matter tells spacetime how to curve, and spacetime tells matter how to move. This video provides a conceptual bridge between the mathematics of the metric tensor and observable astrophysical phenomena.

Knowledge Checkpoint

• Understand the Equivalence Principle and how it demands the bending of light in a gravitational field.
• Define a "geodesic" and explain why light follows curved paths despite having zero rest mass.
• Explain how Einstein's field equations relate spacetime curvature to mass-energy density.

---

Module 2: Cosmology and the Evidence for Dark Matter

Overview

This module establishes the cosmological context of dark matter. Students will explore observational anomalies—primarily galactic rotation curves and galaxy cluster velocities—that cannot be explained by visible matter, showing why dark matter is an essential pillar of our standard model of cosmology ($\Lambda$CDM).

Recommended Videos

Why this video

This comprehensive, academically rigorous video explains why dark matter is treated as an observed physical phenomenon rather than a loose hypothesis. It details the history of galaxy rotation curves, the contributions of Vera Rubin and Fritz Zwicky, and how alternative theories like MOND (Modified Newtonian Dynamics) struggle to explain the broader range of observations that dark matter easily accounts for.

---

Why this video

This lecture focuses on the mathematical derivation of galactic rotation curves. It shows how Newtonian dynamics predicts that outer stellar orbital speeds should fall off as $v \propto r^{-1/2}$ based on visible mass distributions, and contrasts this prediction with the observed flat rotation curves ($v \approx \text{constant}$), proving the existence of a massive, invisible dark matter halo.

---

Why this video

This video profiles Fritz Zwicky’s 1933 discovery of mass discrepancies in the Coma Cluster and Vera Rubin’s subsequent work on spiral galaxies. It serves as an excellent introduction to why baryonic (normal) matter is insufficient to bind large-scale structures, setting the stage for gravitational lensing as the ultimate mapping tool.

Knowledge Checkpoint

• Derive the expected Newtonian rotation velocity $v(r)$ for stars outside a galaxy's optical radius and explain the flat curves observed in real data.
• Explain how Fritz Zwicky used the Virial Theorem to conclude that the Coma Cluster contains far more mass than is visible.
• Summarize the primary components of the $\Lambda$CDM model and the estimated percentages of dark matter, dark energy, and baryonic matter.

---

Module 3: Physics of Strong Gravitational Lensing

Overview

When a massive object (like a galaxy or cluster) aligns closely with a background light source, it acts as a cosmic lens. This module covers the core geometric optics of strong gravitational lensing, the derivation of the lens equation, the Einstein radius, and the formation of highly distorted images, including Einstein rings, arcs, and multiple images.

Recommended Videos

Why this video

This video explains the thin-lens approximation and the fundamental lens equation, which relates the observed image position ($\theta$) to the true source position ($\beta$) through the deflection angle ($\alpha$). It offers a great, concise mathematical visualization of the geometric configuration of the observer, lens, and source planes.

---

Why this video

This video offers clean 3D animations and astronomical observations of Einstein rings and gravitational arcs. It illustrates how perfect alignment yields rings, while small offsets produce multiple discrete images or stretched arcs.

---

Why this video

This video presents a real-world case study of strong lensing: the discovery of Earendel, a star magnified thousands of times by an intervening galaxy cluster. It illustrates "critical curves" and "caustics"—boundaries in the source and lens planes where magnification theoretically diverges to infinity.

---

Supplemental Mathematics: Derivation of the Lens Equation

Curriculum Gap Note: Existing video content focuses on qualitative explanations of strong lensing. To master this topic, work through the following mathematical derivation:

Source (S) o [True Position: beta] / \ / \ / \ Deflection Angle (alpha) / \ / (L) Lens (Mass M) / |

/ | [Observed Position: theta] / | o------------+--------------------o Observer (O) <- - - D_L - ->[Lens Plane] < - - - - - - D_S - - - - - - - - >[Source Plane]

1. Light Deflection Angle ($\alpha$): For a point mass $M$ with an impact parameter $\xi$ (where $\xi = D_L \theta$), General Relativity predicts a deflection angle of: $\alpha(\theta) = \frac{4GM}{c^2 D_L \theta}$
2. The Lens Equation (Thin Lens Approximation): Using the geometry of the diagram above, the relationship between the true source angular position $\beta$ and the observed angular position $\theta$ is: $\beta = \theta - \alpha'(\theta)$ where the scaled deflection angle is: $\alpha'(\theta) = \frac{D_{LS}}{D_S} \alpha(\theta)$ Therefore: $\beta = \theta - \frac{D_{LS}}{D_L D_S} \frac{4GM}{c^2 \theta}$
3. The Einstein Radius ($\theta_E$): When the source is perfectly aligned behind the lens ($\beta = 0$), the lensed image forms a symmetric ring. Solving for $\theta$ under this condition yields the Einstein Radius: $\theta_E = \sqrt{\frac{4GM}{c^2} \frac{D_{LS}}{D_L D_S}}$

Knowledge Checkpoint

• Derive the lens equation from the geometric relationship between the observer, lens, and source planes.
• Calculate the Einstein radius (in arcseconds) for a galaxy of mass $10^{12} M_{\odot}$ acting as a lens at $D_L = 1 \text{ Gpc}$ with a source at $D_S = 2 \text{ Gpc}$.
• Explain how multiple images are produced by solving the lens equation for non-zero source offsets ($\beta \neq 0$).

---

Module 4: Physics of Weak Gravitational Lensing

Overview

Unlike strong lensing, which requires near-perfect alignment and produces dramatic arcs, weak gravitational lensing is a statistical phenomenon. It slightly distorts the shapes of almost all background galaxies. This module introduces the mathematical framework of weak lensing, detailing how the gravitational potential field creates shear (stretching) and convergence (magnification) across the sky.

Recommended Videos

Why this video

This academic lecture explains the physics of the weak lensing regime. It demonstrates how an individual galaxy's light is slightly distorted when passing through a mass distribution, showing how astronomers use statistical averages of multiple background galaxies to reconstruct the shape profile of a lens.

---

Why this video

This lecture provides the mathematical foundation of cosmic shear. Alonso presents the key observational equation: the observed ellipticity ($e_{\text{obs}}$) is the sum of the true unlensed ellipticity ($e_{\text{true}}$) and the gravitational shear ($\gamma$). Since galaxy orientations are randomly distributed, averaging $e_{\text{obs}}$ over many galaxies allows astronomers to isolate the shear term.

---

Why this video

This long-form, highly advanced academic lecture goes deep into the observational challenges of weak lensing. Bernstein explains how systematic errors, telescope optics (Point Spread Function), and detector noise distort galaxy shapes, and discusses the statistical correction methods required to extract clean lensing signals.

---

Supplemental Mathematics: Shear vs. Convergence

Curriculum Gap Note: To bridge the gap between qualitative overviews and Bernstein's highly advanced lecture, study this mathematical breakdown of the Jacobian matrix mapping background sources to observed images.

The distortion of an image is described by the Jacobian matrix $A(\theta)$, which maps coordinates in the source plane to coordinates in the image plane:

$A(\theta) = \frac{\partial \beta}{\partial \theta} = \begin{pmatrix} 1-\kappa-\gamma_1 & -\gamma_2 \\ -\gamma_2 & 1-\kappa+\gamma_1 \end{pmatrix}$

This matrix is decomposed into two main physical components:

1. Convergence ($\kappa$):

   • Definition: A scalar field representing the isotropic magnification of the source (changes size and flux, but preserves shape).
   • Physics: Proportional to the projected local surface mass density $\Sigma(\theta)$ divided by the critical surface mass density $\Sigma_{\text{crit}}$: $\kappa(\theta) = \frac{\Sigma(\theta)}{\Sigma_{\text{crit}}}$

2. Shear ($\gamma = \gamma_1 + i\gamma_2$):

   • Definition: A complex trace-free tensor describing anisotropic stretching (stretches circular sources into ellipses).

   • Physics: Directly related to the tidal gravitational forces of the lens. $\gamma_1$ represents stretching along the x/y axes, while $\gamma_2$ represents stretching at $45^\circ$ angles.

     Unlensed Convergence (kappa) Shear (gamma) [O] ( O ) ( O ) Circular Source Isotropic expansion Anisotropic stretching

Knowledge Checkpoint

• Write down the deformation matrix $A(\theta)$ and explain the geometric impact of convergence ($\kappa$) versus shear ($\gamma_1, \gamma_2$).
• Explain why the average unlensed ellipticity $\langle e \rangle$ of a large population of galaxies is assumed to be zero.
• Explain how the Point Spread Function (PSF) of a telescope can mimic weak lensing shear, and describe one method used to correct for it.

---

Module 5: Dark Matter Mass Reconstruction and Mapping

Overview

This capstone module synthesizes strong and weak lensing data into a unified framework for mapping dark matter. Students will learn the step-by-step pipeline used by cosmologists to convert observed galaxy ellipticities into 2D and 3D density maps of galaxy clusters, with a focus on computational modeling using software like lenstronomy and Lenstool.

Recommended Videos

Why this video

This computational walkthrough shows how researchers reconstruct mass distributions in practice. The tutorial demonstrates how to use the open-source software Lenstool to model galaxy-galaxy lensing using Singular Isothermal Sphere (SIS) and Ellipsoid (SIE) mass profiles, showing how observational data is fitted to physical parameters.

---

Why this video

Using recent James Webb Space Telescope data, this video explains how weak lensing measurements of minor 1% distortions are combined over millions of background galaxies to reconstruct the "cosmic web" of dark matter filaments.

---

Why this video

Though brief, this clip displays the overlay mapping of the Bullet Cluster (1E 0657-568). It shows the critical offset between baryonic matter (hot gas seen in pink by Chandra X-ray) and total mass (mapped in blue by weak lensing), which serves as direct proof of dark matter.

---

Computational Pipeline: Mass Mapping Walkthrough

Curriculum Gap Note: To address the lack of technical walkthroughs for lensing reconstruction, review this Python/lenstronomy pipeline overview.

To map dark matter, astronomers must invert the observed shear field $\gamma(\theta)$ to find the underlying mass density field (convergence $\kappa(\theta)$). This is done using the Kaiser-Squires Inversion algorithm:

[Background Galaxies] ---> Measure Ellipticities (e_obs) ---> Grid & Average to create Shear Field (gamma_1, gamma_2) | [Dark Matter Map] <--- Project Convergence (kappa) <--- Kaiser-Squires Inversion (Fourier Domain)

Step-by-Step Computational Workflow:

1. Catalog Generation: Collect shapes, positions, and estimated redshifts of background galaxies.
2. PSF Correction: Deconvolve the telescope's point spread function from the galaxy image shapes to extract the pure gravitational shear.
3. Fourier Transformation: In the Fourier domain, shear ($\tilde{\gamma}$) and convergence ($\tilde{\kappa}$) are related by a simple algebraic projection. First, transform the real-space shear fields: $\tilde{\gamma}_1(\vec{k}), \tilde{\gamma}_2(\vec{k}) = \text{FFT}(\gamma_1(\vec{\theta}), \gamma_2(\vec{\theta}))$
4. Kernel Application (Inversion): Compute the convergence in Fourier space: $\tilde{\kappa}(\vec{k}) = \frac{(k_1^2 - k_2^2) \tilde{\gamma}_1(\vec{k}) + 2k_1 k_2 \tilde{\gamma}_2(\vec{k})}{k_1^2 + k_2^2}$
5. Inverse Fourier Transform: Transform $\tilde{\kappa}(\vec{k})$ back to real space to yield the mass density map: $\kappa(\vec{\theta}) = \text{IFFT}(\tilde{\kappa}(\vec{k}))$
6. Visualization: Plot the 2D convergence map $\kappa(\vec{\theta})$ to reveal dark matter overdensities.

Knowledge Checkpoint

• Explain the steps of the Kaiser-Squires inversion algorithm and how it extracts the convergence field from shear measurements.
• Describe the physical significance of the Bullet Cluster separation and why MOND struggles to explain this observation.
• Set up a basic lensing model using Lenstool or Python's lenstronomy package, specifying the lens mass profile (e.g., SIS or NFW) and the source profile.

---

Course Map

---

Key People Index

---

Final Self-Assessment

Complete this comprehensive self-assessment to verify your mastery of the physics of gravitational lensing and dark matter mapping.

• Equivalence Principle: I can explain how Einstein used an accelerating elevator thought experiment to predict that gravity bends light.
• Geodesic Trajectory: I can mathematically define a null geodesic and explain why photons travel along them in curved spacetime.
• Rotation Curves: I can derive the flat velocity curve relationship using a dark matter halo mass profile ($M(r) \propto r$) and contrast it with the Keplerian fall-off ($v \propto r^{-1/2}$).
• Lens Geometry: I can derive the lens equation $\beta = \theta - \alpha'(\theta)$ from a diagram of the observer, lens, and source planes.
• Einstein Radius: I can calculate the angular radius of an Einstein Ring given the mass of a lens and the angular diameter distances $D_L$, $D_S$, and $D_{LS}$.
• Shear vs. Convergence: I can explain the physical and geometric difference between isotropic convergence ($\kappa$) and complex shear ($\gamma = \gamma_1 + i\gamma_2$).
• Jacobian Transformation: I can write out the magnification matrix $A(\theta)$ and compute the magnification factor $\mu = \frac{1}{\det A}$ in terms of $\kappa$ and $\gamma$.
• Kaiser-Squires Inversion: I can explain how the Kaiser-Squires algorithm operates in the Fourier domain to reconstruct convergence fields from measured shear.
• Bullet Cluster Physics: I can explain how the spatial offset between hot gas (X-ray emission) and total mass (weak lensing reconstruction) serves as empirical proof of dark matter.
• Modeling Frameworks: I can describe the setup, input data, and basic parameters required to run a mass reconstruction model in computational tools like lenstronomy or Lenstool.