🔭 astrophysics

Simulating the LOcal Web (SLOW): VII. Intergalactic magnetic field models for multi-messenger applications

This paper presents and validates new intergalactic magnetic field models derived from the constrained SLOW cosmological simulation, including a novel algorithm for optimizing line-of-sight extraction toward galaxies, to enhance the interpretation of multi-messenger observations involving ultra-high-energy cosmic rays and gamma-rays.

Original authors: Johannes Stoiber, Klaus Dolag, Francesca Capel, Benjamin Seidel, Michael Kachelrieß, Ludwig M. Böss, Jenny G. Sorce

Published 2026-07-09

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Johannes Stoiber, Klaus Dolag, Francesca Capel, Benjamin Seidel, Michael Kachelrieß, Ludwig M. Böss, Jenny G. Sorce

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 as a vast, invisible ocean. In this ocean, there are currents and winds that we can't see, but they affect everything that swims through them. In the world of astronomy, these "currents" are magnetic fields that stretch across the empty spaces between galaxies. Scientists call this the Intergalactic Magnetic Field (IGMF).

This paper is about creating a better, more accurate map of these invisible currents to help us understand how high-energy particles (like cosmic rays and gamma rays) travel through the universe.

Here is a breakdown of what the researchers did, using simple analogies:

1. The Problem: The "Blind" Map

For a long time, scientists trying to track these high-energy particles had to guess what the magnetic fields looked like.

The Old Way: Some simulations were like drawing a map of a city without ever looking at the actual city. They used random guesses for how the "buildings" (galaxies) were arranged.
The Size Problem: Other maps were too small. If you want to trace a path from New York to London, but your map only covers New York, you have to just guess what the ocean looks like in between. This leads to errors.
The "Ghost" Problem: Sometimes, the specific galaxy a particle comes from is too small or faint to be perfectly matched in a computer simulation. It's like trying to find a specific house in a neighborhood, but the map only shows the major streets.

2. The Solution: The "SLOW" Simulation

The researchers used a super-advanced computer simulation called SLOW (Simulating the LOcal Web).

The "Constrained" Advantage: Unlike the old "blind" maps, SLOW is a constrained simulation. Think of it like a 3D puzzle where the pieces are forced to fit the actual shape of our local neighborhood in the universe. It matches the real locations of major galaxy clusters (like the Coma cluster) that we can see through telescopes.
The Result: This gives a much more realistic "terrain" for the magnetic fields to flow through.

3. The New Trick: Finding the "Ideal Position"

This is the paper's biggest innovation. Sometimes, a source of high-energy particles (like a blazar, a super-bright galaxy) is too small to be perfectly pinned down in the simulation.

The Analogy: Imagine you are trying to find a specific tree in a forest, but your map only clearly shows the big oak trees. You know the tree you want is near a big oak, but the map isn't precise enough to show the small tree.
The Algorithm: The team invented a new math trick called "Fuzzy Triangulation." They look at the three biggest, most clearly mapped galaxy clusters near the target. They then ask the computer: "Where in the simulation does the distance between these three big clusters match the real-world distances?"
The "Ideal Position": Once they find that spot, they treat it as the "ideal" location for the small target galaxy. It's like using the positions of three big lighthouses to triangulate the exact location of a small boat, even if the boat itself is hard to see.

4. Testing the Map: The "Gamma-Ray Cascade"

How do they know their new map is good? They test it using a phenomenon called a gamma-ray cascade.

The Analogy: Imagine throwing a bright flashlight beam (a high-energy gamma ray) into the foggy ocean. If the water is calm, the beam goes straight. But if there are strong magnetic currents (the fog), the beam gets scattered and twisted.
The Test: When these high-energy rays hit the "fog" of the universe, they create a shower of particles that eventually turn into lower-energy light (GeV gamma rays). If the magnetic field is strong, it bends the path so much that we see less light than expected.
The Finding: The researchers ran their new magnetic field maps through a computer program (ELMAG) to simulate this scattering. They found that their simulated magnetic field (the one that comes directly from the physics of the SLOW simulation) matched the real-world observations best. It predicted the right amount of "scattering" to explain what telescopes actually see.

5. Why This Matters

Better Accuracy: By using the "Ideal Position" trick, they reduced the uncertainty in their predictions. It's like switching from a blurry, hand-drawn sketch to a high-definition GPS route.
Multi-Messenger Science: This helps scientists who study "multi-messenger" astronomy (using light, particles, and gravitational waves together). With a better map of the magnetic currents, they can trace high-energy particles back to their true origins more accurately.

Summary

The authors built a high-definition, realistic map of the magnetic "wind" in our local universe. They created a clever method to find the exact location of small, hard-to-see galaxies within this map. When they tested this map against real observations of how light gets scattered in space, it worked better than previous models. This gives scientists a sharper tool to understand the invisible forces shaping our universe.

Technical Summary: Simulating the LOcal Web (SLOW) VII. Intergalactic magnetic field models for multi-messenger applications

Problem Statement
The propagation of ultra-high-energy cosmic rays (UHECRs) and ultra-high-energy gamma-rays is heavily influenced by the intergalactic magnetic field (IGMF). Accurately characterizing the IGMF across cosmic large-scale structure is essential for interpreting multi-messenger observations and constraining magnetogenesis scenarios. However, existing magnetic field models for multi-messenger studies suffer from two primary limitations:

Unconstrained Initial Conditions: Many simulations do not reproduce the observed local large-scale structure, leading to inaccurate line-of-sight profiles for specific astrophysical sources.
Scale and Fidelity Issues: Even constrained simulations are often too small to encompass distances to relevant sources (e.g., blazars) without resorting to artificial techniques like mirroring or concatenating volumes, which degrade the fidelity of the magnetic field profile. Furthermore, many studies rely on simplified representations (e.g., Gaussian random fields) rather than physical models extracted from simulations. Additionally, identifying simulated counterparts for individual galaxy-mass halos (potential UHECR/TeV sources) is difficult because initial conditions are linearly reconstructed and only constrain scales above ~3–4 Mpc.

Methodology
The authors present magnetic field models derived from the SLOW (Simulating the LOcal Web) suite of constrained cosmological simulations, specifically utilizing two components:

SLOW-AGN15363: A high-resolution simulation including Active Galactic Nuclei (AGN) and supermassive black holes (SMBH) for identifying galaxy and cluster positions, but lacking magnetic fields.
SLOW-CR30723: A pure non-radiative MHD simulation with higher resolution, containing magnetic fields and cosmic ray evolution, but lacking AGN feedback.
These simulations share initial conditions based on the Cosmicflows-2 catalog (peculiar velocities), allowing for cross-matching of halos.

The study analyzes six magnetic field models:

$B_{sim}$ : The raw magnetic field from the MHD simulation (seeded by a uniform primordial magnetic field of $10^{-14}$ G).
Rescaled Models: Five models rescaled based on physical assumptions to address resolution limits in filaments:
- $B_{\beta}$ : Constant plasma- $\beta$ (thermal to magnetic pressure ratio).
- $B_{F}$ : Equilibrium with turbulent pressure.
- $B_{ff}$ : Scaled via ideal MHD flux-freezing (electron density).
- $B_{dyn, \downarrow}$ and $B_{dyn, \uparrow}$ : Scaled based on turbulent dynamo amplification scenarios (lower and upper limits).

Novel Algorithm: "Ideal Position" Determination
To address the inability of initial conditions to constrain individual galaxy-mass halos, the authors introduce a novel algorithm to determine an "ideal position" ( $\mathbf{x}_{ideal}$ ) within the simulation for a given observed source:

Triangulation: The algorithm identifies the three nearest observed galaxy clusters to the source.
Fuzzy Triangulation: It calculates a "shifted centroid" by adding the observed relative position vectors to the simulated cluster positions.
Optimization: A parameter $R$ (a sum of Gaussian-weighted differences between observed and simulated distances) is minimized on a grid around the centroid. This "fuzzy" approach allows for a global minimum that accounts for large-scale drifts between the simulated and observed structures, rather than enforcing strict equality.
Uncertainty Estimation: The authors estimate uncertainty by calculating magnetic field profiles along lines-of-sight to random positions offset from the ideal position (using 2D sky-plane shifts, 3D spatial shifts, and selecting realistic subhalos hosting SMBHs).

Key Results

Magnetic Field Properties: The models span a wide range of filling factors. While all models converge in the cores of galaxy clusters (matching observations of the Coma cluster), they diverge significantly in filaments and voids. The raw simulation ( $B_{sim}$ ) predicts significantly lower field strengths in filaments compared to rescaled models, often dropping below 0.1 nG.
Gamma-Ray Cascade Analysis: Using the Monte Carlo code ELMAG3.03, the authors simulated the electromagnetic $\gamma$ $γ$ -ray cascade for a source at the "ideal position" of Mrk 501.
- Models with high filling factors in voids ( $B_{dyn, \uparrow}$ , $B_F$ , $B_{ff}$ , $B_{\beta}$ ) deflect the cascade particles entirely, suppressing GeV emission.
- The raw simulation ( $B_{sim}$ ) and the lower-limit dynamo model ( $B_{dyn, \downarrow}$ ) allow for significant secondary emission, consistent with the lower limits derived from $\gamma$ -ray cascade observations (e.g., $B > 1.3 \times 10^{-17}$ G).
- The results suggest that $B_{sim}$ best reproduces the IGMF derived from electromagnetic $\gamma$ -ray cascade constraints, implying that the filling factor required to deflect particles is not as high as some rescaled models suggest.
Uncertainty Reduction: The "ideal position" method yields magnetic field profiles with lower variance (measured by mean quartile distance in cascade spectra) compared to random lines-of-sight. Specifically, using lines-of-sight ending in realistic galaxy halos ("subhalos" method) provided the most accurate and least uncertain cascade spectra.

Significance and Claims
The paper claims to provide the astroparticle physics community with accurate, physically motivated IGMF models derived from a constrained simulation of the local Universe.

Methodological Advance: The authors claim to be the first to propose an algorithm for identifying "ideal positions" for galaxies below the constraining scale of initial conditions, thereby enabling accurate line-of-sight magnetic field extraction for specific sources.
Robustness: The study demonstrates that using these constrained "ideal position" models reduces the uncertainty in $\gamma$ -ray cascade simulations compared to using random or unconstrained magnetic field profiles.
Future Utility: The large-scale structure drift observed between simulated and real clusters, which the algorithm exploits, is suggested as a potential route for refining the simulation's constrained initial conditions in future work.
Limitations: The authors modestly note that current models are based on a uniform primordial magnetic field seed. They acknowledge that future work must explore astrophysical seeding scenarios and stochastic primordial seeds (with/without helicity) to fully capture the topology of the IGMF, which is crucial for observables like $\gamma$ -ray halos. The models are currently limited to the local Universe (box size 500 $h^{-1}$ Mpc) and do not yet cover all magnetogenesis scenarios.