🔭 astrophysics

Probing Cosmic Magnetism with Rotation Measure-Squared-Galaxy Cross-Correlations

This paper proposes a new tomographic method to detect and reconstruct the evolution of large-scale cosmic magnetic fields by cross-correlating the squared Faraday rotation measures of background sources with foreground galaxy densities, a technique validated by Illustris-TNG simulations and forecasted to be detectable with current surveys and future SKA data.

Original authors: Zekai Zhang, Adam Lidz

Published 2026-06-08

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Zekai Zhang, Adam Lidz

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe is a giant, invisible ocean. We know there are currents and waves in this ocean (magnetic fields), but they are so faint and spread out that we can't see them directly. We only know they exist because they twist the light from distant lighthouses (galaxies and quasars) as that light travels through the ocean. This twisting is called the Faraday Rotation.

However, there's a problem. When we look at a lighthouse, the light has traveled through many different parts of the ocean: the water right next to the lighthouse, the water right next to our own telescope (in the Milky Way), and the vast, empty spaces in between. It's like trying to hear a whisper in a crowded room; you hear the whisper, but you also hear the chatter of the crowd and the noise of the room itself. It's very hard to tell which part of the sound came from where.

The New Idea: Squaring the Signal
The authors of this paper propose a clever new trick to isolate the "whisper" from the "crowd." Instead of just measuring how much the light twists (the Rotation Measure, or RM), they suggest squaring that number (RM²) and then looking for patterns.

Think of it like this:

The Old Way: If you try to measure the average twist, positive twists cancel out negative twists, and the noise from the crowd (our galaxy) and the lighthouse itself gets in the way. It's like trying to find a specific person in a crowd by asking everyone to shout "Hello!"—you just get a wall of noise.
The New Way (RM²): By squaring the twist, everything becomes positive. It's like turning all the voices in the crowd into a single, steady hum. Now, if you look for where this hum is loudest, you can find the specific spots where the "crowd" (the magnetic fields) is thickest.

The "Stacking" Trick
To find these patterns, the scientists used a method called stacking. Imagine you have a map of the universe with millions of galaxies.

You take a tiny piece of the magnetic "hum" map for every single galaxy.
You line them all up so the galaxy is in the exact center of each piece.
You add them all together.

If the magnetic fields are just random noise, the pieces will cancel each other out when you stack them. But, if the magnetic fields are actually stronger near galaxies (like a halo of invisible energy surrounding them), the "hum" will get louder and louder as you stack more and more galaxies. The paper shows that this "loudness" does indeed build up, proving that magnetic fields are concentrated around galaxies.

What They Found in the Simulations
The authors tested this idea using a super-computer simulation called Illustris-TNG, which is like a video game that creates a fake universe with realistic physics, including magnetic fields.

Time Travel: They looked at this fake universe at different times in its history (from when it was young to now). They found that the magnetic "hum" gets much, much louder as the universe gets older.
Why? It's because galaxies act like cosmic generators. Over billions of years, the movement of gas and stars inside galaxies (dynamo processes) and powerful winds blowing out from them have amplified the magnetic fields, making them stronger and more widespread.

Why This Matters (and Why It's Better)
The paper highlights two main reasons why this new "squared" method is superior to previous attempts:

No Noise Bias: Previous methods tried to measure the "absolute size" of the twist. The problem is that even if there is no real magnetic field, the instrument's own "static" (noise) will make the number look bigger than it is. It's like a broken scale that always adds 5 pounds to your weight. The new "squared" method is immune to this; the static doesn't trick it.
Clearer Picture: Because they can measure this "hum" at different distances (redshifts), they can create a 3D movie of how magnetic fields have grown over the history of the universe.

The Future
The authors predict that with upcoming giant radio telescopes (like the SKA), which will measure the magnetic twist of millions of galaxies, this method will work beautifully. They estimate that we will be able to detect this signal with very high confidence, allowing us to finally map out the "magnetic cosmic web" and understand how the universe became magnetized.

In short: They found a way to turn the confusing noise of cosmic magnetic fields into a clear, measurable signal by squaring the data and stacking it, revealing that magnetic fields have been growing stronger around galaxies for billions of years.

Technical Summary: Probing Cosmic Magnetism with Rotation Measure-Squared-Galaxy Cross-Correlations

Problem Statement
Cosmic magnetic fields permeate the universe, from planets to the intergalactic medium (IGM), yet their origin, evolution, and spatial structure remain poorly understood. A primary observational tool, Faraday rotation, measures the Rotation Measure (RM) of background sources. However, RM provides only an integral constraint along the line of sight, mixing contributions from the source, the Milky Way, and intervening large-scale structure (LSS). Furthermore, direct cross-correlations between RM and galaxy distributions are problematic because the RM field can be positive or negative, leading to cancellations that wash out the signal. Previous attempts to circumvent this by correlating the absolute value of RM ($|RM|$) with galaxies introduce non-differentiability at zero and, crucially, suffer from noise bias, where the expectation value of $|RM + n|$ (where $n$ is noise) differs from $|RM|$.

Methodology
The authors propose a new statistic: the cross-correlation between the squared Rotation Measure field ( $RM^2$ ) and the projected galaxy overdensity ( $g$ ), denoted as $\langle RM^2 \times g \rangle$ . This approach is analogous to the "projected fields" estimator used for the kinetic Sunyaev-Zel'dovich (kSZ) effect ( $\langle kSZ^2 \times g \rangle$ ).

Theoretical Framework: The authors derive that the angular cross-power spectrum of $RM^2$ and $g$ relates to a bispectrum involving two copies of the electron-density-weighted magnetic field ( $\tilde{B}_{\perp,\parallel}$ ) and one copy of the galaxy overdensity. In real space, the correlation function $w_{RM^2,g}(r_\perp)$ is shown to be a line-of-sight projection of the correlation between the squared magnetic field and galaxies.
Analytic Approximation: An approximate analytic model is developed, expressing the correlation as the product of the electron-galaxy cross-power spectrum and an integral over the power spectrum of the electron-density-weighted magnetic field. This model suggests the signal's shape is governed by electron clustering, while its amplitude traces the projected magnetic field strength.
Simulations: The authors utilize the Illustris-TNG cosmological magnetohydrodynamic (MHD) simulations (specifically the TNG300-3 run) to construct mock RM maps. They generate $RM^2$ maps by integrating the electron density and magnetic field along lines of sight, then apply a stacking method around dark matter halos to estimate the cross-correlation function $w_{RM^2,g}(r_\perp; z)$ .
Noise Analysis: The paper explicitly compares the $\langle RM^2 \times g \rangle$ estimator against the $\langle |RM| \times g \rangle$ estimator. It demonstrates mathematically and via simulation that $\langle RM^2 \times g \rangle$ is immune to noise bias (assuming noise is uncorrelated with galaxies), whereas $\langle |RM| \times g \rangle$ suffers from systematic suppression due to the non-linear nature of the absolute value operation on noisy data.

Key Results

Signal Detection in Simulations: Using Illustris-TNG, the authors detect a strong excess $RM^2$ signal around halos. The signal increases by approximately three orders of magnitude from redshift $z \sim 3$ to $z \sim 0$ .
Redshift Evolution: The evolution is driven by the amplification of magnetic fields in overdense regions via turbulent dynamos and the transport of magnetized plasma by AGN and galactic outflows. The signal is dominated by highly magnetized regions within dark matter halos rather than the diffuse IGM.
Analytic Model Validation: The approximate analytic model (Eq. 3.3) successfully reproduces the scale dependence and redshift evolution of the simulated signals, confirming that the amplitude is primarily set by the electron-density-weighted magnetic field strength.
Smoothing Scale Dependence: The signal is highly sensitive to the smoothing scale of the RM measurements. Because the magnetic field power spectrum contains strong small-scale fluctuations, coarser smoothing significantly suppresses the signal amplitude.
Resolution Dependence: Higher-resolution simulations (TNG300-1) yield lower signal amplitudes compared to lower-resolution runs (TNG300-3). This is attributed to better resolution of magnetic field reversals, which cause cancellations in the line-of-sight projection, a feature missed in lower-resolution runs.
Forecasts: Forecasts for current and future surveys (e.g., SKA, POSSUM, NVSS) suggest that $\langle RM^2 \times g \rangle$ $⟨ R M^{2} \times g ⟩$ is detectable with high significance.
- An SKA-like survey with $\sim 10^7$ RM sources and $\sim 10^6$ galaxies per redshift bin could detect the signal across a wide range of scales ($0.5 - 50$ Mpc/h) and redshifts.
- Even current datasets (e.g., NVSS combined with WISE) may allow for initial detections, particularly at low redshifts, though with lower signal-to-noise ratios.
- The $\langle RM^2 \times g \rangle$ statistic yields a slightly higher signal-to-noise ratio than $\langle |RM| \times g \rangle$ and avoids the latter's noise bias.

Significance and Claims
The paper claims to provide a robust, tomographic method for probing the evolution of large-scale cosmic magnetic fields. By utilizing $\langle RM^2 \times g \rangle$ , the method:

Isolates Redshift Evolution: It cleanly extracts contributions from magnetic fields near the redshifts of foreground tracer galaxies, filtering out intrinsic source and Milky Way contributions which act only as noise variance.
Avoids Systematic Biases: Unlike $\langle |RM| \times g \rangle$ , the squared estimator is unbiased by measurement noise, making it more reliable for quantitative analysis.
Connects to Fundamental Physics: The statistic links directly to the underlying bispectra of magnetic fields and electron densities, offering a clearer path to interpreting the physical properties of the magnetized cosmic web.
Enables Future Science: The authors argue that with upcoming data from the Square Kilometer Array (SKA) and large galaxy surveys (e.g., DESI), this technique will enable a census of cosmic magnetic energy density across cosmic time, addressing open questions regarding magnetic field seeding and amplification.

The paper remains modest regarding the interpretation of the signal's origin, noting that while the signal traces the evolution of magnetic energy density, it is currently dominated by magnetized plasma within halos. Future work is required to refine galaxy mocks, match observational smoothing scales precisely, and account for sample variance and non-Gaussian covariance terms.

More like this