Perturbative and numerical study of nonlinear relativistic effects in weak lensing

This paper demonstrates that standard weak lensing formalisms are incomplete beyond linear order by deriving and numerically quantifying second-order relativistic corrections—specifically Sachs-basis rotation and frame-dragging effects—which break the degeneracy between rotation and shear B-modes and dominate the B-mode power spectrum on large scales, despite their overall small impact on observed galaxy ellipticity.

The Gist

The Big Picture: Looking at the Universe Through a Distorted Lens

Imagine you are standing in a vast, dark forest (the Universe) trying to look at a distant campfire (a galaxy). Between you and the fire, there are trees, bushes, and uneven ground (dark matter and gravity).

In Weak Lensing, astronomers use the distortion of the campfire's shape to map out the invisible trees and bushes in between. If the fire looks stretched into an oval, it means gravity pulled on the light.

For decades, scientists have used a "standard rulebook" to interpret these distortions. This rulebook assumes that gravity acts like a simple, flat sheet of rubber: it stretches the image in one direction and squeezes it in another, but it doesn't twist it or rotate it.

This paper says: "That rulebook is incomplete."

The authors, Matteo Magi, Francesca Lepori, and Julian Adamek, show that when you look closely enough (using the full power of Einstein's General Relativity), the "rubber sheet" isn't just stretching; it's also twisting, rotating, and dragging the light in ways the old rulebook missed.

---

Analogy 1: The Moving Walkway (The "Standard" vs. "Real" View)

The Old Way (Standard Formalism):
Imagine you are walking on a straight, flat moving walkway at an airport. You look at a sign on the wall. If the walkway vibrates slightly, the sign might look a bit blurry or stretched. The old rulebook says, "Okay, the sign is stretched. Let's measure how much it's stretched and calculate the vibration." It assumes the sign stays perfectly upright relative to the floor.

The New Way (Jacobi Map & Parallel Transport):
Now, imagine the walkway isn't just vibrating; it's also twisting slightly as it moves, like a corkscrew.

• If you hold a pen pointing "North" at the start of the walkway, and you walk to the end, the pen might now be pointing "Northeast" because the walkway twisted underneath you.
• The old rulebook ignores this twist. It assumes the pen still points North.
• The new rulebook (the Jacobi Map) accounts for the twist. It says, "Hey, the pen rotated because the path itself twisted."

Why does this matter?
In the old rulebook, "stretching" (Shear) and "twisting" (Rotation) are mathematically linked in a specific way. If you see a twist, you expect a specific amount of stretch. But because the path actually twists (due to Frame Dragging and Parallel Transport), that link breaks. The twist and the stretch become independent things.

---

Analogy 2: The Spinning Ice Skater (Frame Dragging)

One of the most exciting discoveries in this paper is about Frame Dragging.

Imagine a massive, spinning ice skater (a cluster of galaxies) spinning in a pool of water (space-time).

• Newtonian Physics (Old view): The water stays still; only the skater moves.
• Einstein's Relativity (New view): The spinning skater actually drags the water around with them. The water near the skater starts swirling.

In the Universe, massive objects spinning or moving fast drag the fabric of space-time with them. This is called Frame Dragging.

• The Paper's Finding: This "swirling water" creates a specific type of distortion in the light from distant galaxies. It creates a "twist" (B-mode) that the old rulebook completely missed.
• The Result: On very large scales (looking at huge chunks of the sky), this frame-dragging effect is actually the dominant source of these twists, even more than the standard stretching effects!

---

Analogy 3: The Messy Kitchen Counter (Scalar vs. Vector)

The authors break down the distortions into two types of ingredients:

1. Scalar Perturbations (The Flour): These are the standard "clumps" of matter (like piles of flour). They create the main stretching (Shear E-modes) that we usually measure.
2. Vector Perturbations (The Whisk): These are the "swirls" or currents (like the whisk moving through the flour). These are usually tiny and ignored.

The Surprise:
The paper shows that even if you start with just "Flour" (scalar matter), the act of mixing it (non-linear gravity) creates "Whisk" movements (vector modes) at the second order.

• The Twist: These "Whisk" movements (Frame Dragging) are the main culprit behind the "B-mode" distortions (the odd-parity twists) on large scales.
• The Catch: While these effects are huge in theory, they are still very small in practice. The paper calculates that for a galaxy at a moderate distance, these new effects change the measurement by about 5%.

---

The "Reduced Shear" Problem: The Onion Layers

There is one final complication.
When we look at a galaxy, we don't measure the "Shear" (the distortion) directly. We measure the Ellipticity (how oval the galaxy looks).

• Linear Theory (Simple Onion): The oval shape is directly proportional to the distortion. Easy peasy.
• Non-Linear Theory (Complex Onion): At high precision, the oval shape is a messy mix of the distortion and the magnification (how much bigger the galaxy looks).

The authors show that this "messy mix" (called Reduced Shear) creates a huge amount of "fake" twisting noise.

• The Problem: The real "Frame Dragging" signal (the cool new physics) is about 1% of the total signal.
• The Noise: The "Reduced Shear" noise is about 100 times larger than the Frame Dragging signal.

The Conclusion:
It's like trying to hear a whisper (Frame Dragging) in a room where someone is screaming (Reduced Shear).

• The authors have successfully calculated exactly what the whisper sounds like.
• They have built a super-accurate map (using N-body simulations) to show where the whisper is.
• But: To actually hear it, future telescopes (like Euclid or the Rubin Observatory) will need to be incredibly good at "turning down the volume" on the screaming noise.

Summary for the General Public

1. The Old Map was Wrong: The standard way we calculate how gravity bends light missed some subtle, real-world effects like twisting and dragging.
2. Space-Time Drags: Massive objects don't just sit there; they drag space-time around them, creating a "swirl" in the light.
3. The Twist is Real: This swirling creates a specific pattern (B-modes) that is different from the stretching pattern.
4. It's Hard to See: While these effects are mathematically significant (changing results by ~5% on large scales), they are tiny in the real world.
5. The Future: We now have the correct formulas and computer simulations to find these effects. The challenge for the next generation of telescopes is to filter out the "noise" of galaxy shapes to finally hear the "whisper" of Einstein's frame-dragging.

In short: The authors have updated the rulebook for reading the Universe's distortions. They found that space-time is more dynamic and "twisty" than we thought, but catching that twist will require the most precise measurements humanity has ever attempted.

Technical Summary

1. Problem Statement

Standard weak lensing formalism relies on the deflection angle to map observed galaxy images to their intrinsic shapes. This approach assumes that the lensing map depends solely on the deflection angle, leading to specific degeneracies in the power spectra of lensing observables:

• The convergence ($\kappa$) and shear E-modes ($\gamma_E$) share the same power spectrum.
• The image rotation ($\omega$) and shear B-modes ($\gamma_B$) share the same power spectrum.

However, the deflection angle is not a gauge-invariant observable. The authors argue that this standard description is incomplete beyond linear perturbation theory, even when only scalar perturbations exist at first order. At second order, scalar perturbations generate vector modes (frame-dragging) and tensor modes, and the parallel transport of the Sachs basis (the local screen basis used to define image orientation) introduces corrections that break the standard degeneracies. The paper aims to quantify these missing relativistic corrections and validate them against numerical simulations.

2. Methodology

A. Theoretical Framework: Jacobi Map Formalism

Instead of the deflection angle, the authors utilize the Jacobi map formalism, which is fully relativistic and coordinate-independent.

• Geodesic Deviation: They solve the geodesic deviation equation for a bundle of null geodesics connecting the source to the observer.
• Sachs Basis: They explicitly track the parallel transport of an orthonormal screen basis (Sachs basis) along the light path. This ensures that image distortions are defined relative to a local rest frame rather than global coordinates.
• Conformal Transformation: To simplify calculations, they work in a conformally rescaled metric where the background is flat FLRW. They derive the evolution of the conformal Jacobi map ($\hat{D}_{IJ}$) and relate it to physical observables (convergence, shear, rotation) via the amplification matrix.
• Perturbation Theory: The analysis is carried out up to second order in perturbation theory within the Poisson gauge. They include:
  • Scalar perturbations (Weyl potential $\Psi$).
  • Vector perturbations (frame-dragging potential $B_\alpha$), sourced by the curl of momentum density (product of first-order density and velocity).
  • Tensor perturbations are neglected as they are subdominant in $\Lambda$CDM.

B. Analytical Derivations

• Second-Order Expressions: They derive explicit full-sky expressions for the image rotation ($\omega$) and shear components ($\gamma_1, \gamma_2$) in terms of metric perturbations.
• E/B-Mode Decomposition: They decompose the shear into E-modes (parity-even) and B-modes (parity-odd).
• Power Spectra: Using the flat-sky and Limber approximations, they compute the angular power spectra ($C_\ell$) for:
  • Rotation ($C_{\omega\omega}$).
  • Shear B-modes ($C_{\gamma_B\gamma_B}$).
  • Ellipticity B-modes (accounting for the non-linear relation between shear and observed ellipticity, known as "reduced shear").
• Key Insight: They derive the analytical difference between the rotation spectrum and the shear B-mode spectrum, showing that the standard relation $\omega \approx \gamma_B$ is broken by terms arising from the parallel transport of the Sachs basis and frame-dragging effects.

C. Numerical Validation

• Simulation: The authors use gevolution, a relativistic N-body code, to generate a high-resolution simulation of a cosmological volume ($1440 \, h^{-1}\text{Mpc}$ box, $1920^3$ particles).
• Ray Tracing: They perform non-perturbative ray tracing backwards from the observer to the source, solving the geodesic and geodesic deviation equations directly.
• Map Generation: They construct full-sky maps (HEALPix $N_{side}=16384$) of the convergence, shear, rotation, and ellipticity.
• Isolation of Effects: By setting specific potentials ($\Psi$ or $B_\alpha$) to zero in different simulation runs, they isolate scalar-induced effects from vector-induced (frame-dragging) effects.

3. Key Contributions

1. Breaking the Degeneracy: The paper proves analytically and numerically that the standard relation $C_{\omega\omega} = C_{\gamma_B\gamma_B}$ (valid in the standard formalism) is violated in the fully relativistic Jacobi map formalism.
2. Frame-Dragging Impact: This is the first analytical derivation and numerical study of the impact of frame-dragging (vector modes) on weak lensing power spectra.
3. Reduced Shear Correction: They quantify the impact of the "reduced shear" (non-linear mixing of convergence and shear) on the observed ellipticity B-modes, showing it dominates the signal compared to relativistic corrections.
4. Numerical Benchmark: They provide the first numerical validation of frame-dragging contributions to lensing convergence and shear power spectra using relativistic N-body simulations.

4. Key Results

• Scalar-Induced Corrections:

  • On large angular scales ($\ell \sim 5$) for sources at $z_s = 0.5$, the difference between the rotation spectrum and the shear B-mode spectrum is approximately 5%.
  • The standard formalism underestimates the convergence power spectrum on large scales compared to the relativistic Jacobi map result.
  • The analytical predictions agree well with numerical simulations for $\ell < 100$.

• Frame-Dragging (Vector) Contributions:

  • Frame dragging does not contribute to image rotation at leading order.
  • However, it dominates the shear B-mode signal on very large scales ($\ell \lesssim 10$), becoming comparable to or larger than the scalar-induced B-modes at low redshift.
  • The contribution of frame dragging to convergence and shear E-modes is negligible (seven orders of magnitude smaller than scalar contributions).

• Ellipticity vs. Shear:

  • While relativistic corrections significantly affect the theoretical shear B-modes, their impact on the observed galaxy ellipticity is suppressed to the ~1% level.
  • The dominant source of ellipticity B-modes is the "reduced shear" effect (convolution of linear convergence and shear), which acts as a contaminant. Detecting the subtle relativistic frame-dragging signal would require precise modeling and subtraction of this dominant term.

5. Significance and Implications

• Precision Cosmology: As next-generation surveys (Euclid, LSST) aim for sub-percent precision, understanding and modeling these nonlinear relativistic effects is crucial to avoid biases in cosmological parameter inference.
• Systematics vs. Physics: The breakdown of the standard degeneracy between rotation and B-modes provides a new diagnostic. While B-modes are often used to identify systematics, this work shows that a portion of the B-mode signal is a genuine General Relativistic effect (frame-dragging).
• Testing General Relativity: The detection of frame-dragging on cosmological scales would provide a unique test of General Relativity in the weak-field, large-scale regime.
• Future Directions: The authors suggest that while current surveys may struggle to detect these effects due to the dominance of the reduced shear and astrophysical contaminants (like intrinsic alignments), future techniques analogous to CMB delensing could be developed to isolate these signals.

In summary, the paper establishes that the standard weak lensing formalism is insufficient for high-precision cosmology at second order. It provides the necessary theoretical framework and numerical tools to account for parallel transport effects and frame-dragging, revealing that while these effects are small (~1-5%), they are physically significant and distinct from standard scalar lensing.