Intergalactic Magnetic Field (IGMF)

• IGMF is a weak, pervasive magnetic field hypothesized to exist in low-density cosmic regions like voids and filaments, measurable via gamma-ray cascades and Faraday rotation.
• Gamma-ray observations, including halo detection and cascade time delays, indicate field strengths around 10⁻¹⁵ G with coherence lengths from kiloparsecs to megaparsecs.
• Detection of the IGMF offers crucial insights into cosmic magnetogenesis, affecting cosmic ray propagation, structure formation, and the evolution of large-scale cosmic magnetic fields.

The Intergalactic Magnetic Field (IGMF) is a pervasive, extremely weak magnetic field hypothesized to exist throughout the low-density regions of the universe, especially in cosmic voids and the filaments of large-scale structure (LSS). Its detection and characterization are central to understanding cosmic magnetogenesis, the propagation of ultra-high-energy particles, and the physics of early universe phase transitions. The IGMF is observable primarily through its subtle effects on high-energy astrophysical phenomena, especially electromagnetic cascades from TeV blazars and gamma-ray bursts (GRBs), Faraday rotation of polarized extragalactic radio emission, and, by implication, its cosmological origin and filling factor.

1. Physical Origins and Theoretical Context

The IGMF is expected to originate from one of two broad classes of seed mechanisms: primordial (inflationary or phase transition era) or astrophysical (outflows from galaxies and AGN). Primordial mechanisms predict a volume-filling, weak, and near-uniform field, possibly generated by quantum fluctuations during inflation or by baryogenesis-related processes in cosmic phase transitions. Conversely, astrophysical mechanisms create more patchy magnetization, localized around galaxies, clusters, and filaments, with fields amplified by turbulence and galactic outflows (Tjemsland et al., 2023, Barai et al., 2018).

The topology, strength ($B_\mathrm{IGMF}$), and coherence length ($\lambda_B$) of the IGMF are crucial for distinguishing between these origins. In primordial scenarios, coherence lengths can approach Mpc scales, while astrophysical scenarios predict both shorter coherence scales and lower volumetric filling fractions. Current empirical constraints require a high filling fraction ($f \gtrsim 0.67$) and a minimum uniform field strength ($B_0 > 5.1 \times 10^{-15}$ G for $\Delta t = 10^4$ yr duty cycle in blazars), strongly favoring a primordial origin over astrophysical models unless extremely efficient and widespread dynamo amplification occurred (Tjemsland et al., 2023).

2. Electromagnetic Cascades and Gamma-Ray Halos

The most robust evidence for the IGMF arises from the paper of electromagnetic cascades initiated when very-high-energy (VHE) gamma rays from blazars or GRBs interact with the extragalactic background light (EBL). These interactions produce electron–positron pairs ($\gamma + \gamma_\mathrm{EBL} \rightarrow e^+ + e^-$), which then inverse Compton scatter cosmic microwave background (CMB) photons up to GeV energies. In the absence of a magnetic field, the cascade follows the original gamma-ray direction, leaving the source point-like in gamma rays. In the presence of an IGMF, the charged $e^\pm$ are deflected, broadening the arrival direction of secondary photons and producing a detectable "halo" around the source (Ando et al., 2010).

The angular size of the resulting halo depends sensitively on $B_\mathrm{IGMF}$, $\lambda_B$, the energy of the observed gamma rays ($E_\gamma$), the redshift ($z$), and the optical depth ($\tau$). The analytic scaling for the halo size, $\theta_\mathrm{halo}$, is

$$\theta_\mathrm{halo} = \begin{cases} 5^\circ (1+z)^{-2} \tau^{-1} (E_\gamma / 10\,\mathrm{GeV})^{-1} (B_\mathrm{IGMF} / 10^{-15}\,\mathrm{G}) & \lambda_B \gg D_e \ 0.4^\circ (1+z)^{-1/2} \tau^{-1} (E_\gamma / 10\,\mathrm{GeV})^{-3/4} (B_\mathrm{IGMF} / 10^{-15}\,\mathrm{G}) (\lambda_B / 1\,\mathrm{kpc})^{1/2} & \lambda_B \ll D_e \end{cases}$$

where $D_e$ is the cascade electron energy-loss length (Ando et al., 2010). Observed halo sizes of a few tenths of a degree around AGN in stacked Fermi-LAT images yield $B_\mathrm{IGMF} \sim 10^{-15}$ G for reasonable parameter choices.

3. Faraday Rotation and Polarimetry

The IGMF also contributes to the rotation measure (RM) of linearly polarized radio sources via the Faraday effect. The expected RM from filaments is modeled as a random walk process due to the relatively short coherence length compared to total path length: $$\mathrm{RM}(z_s) \approx 812 \int_0^{z_s} \frac{n_e(z) B_\parallel(z)}{(1+z)^2} \frac{dl}{dz} dz \quad \mathrm{rad}\,\mathrm{m}^{-2}$$ with $n_e(z)$ the thermal electron density and $B_\parallel(z)$ the line-of-sight magnetic field. In hydrodynamical simulations with turbulence-dynamo amplification, filaments achieve $\langle B \rangle \sim 10$ nG and a coherence length of several hundred kpc, yielding rms RM values $\sim 1$ rad m$^{-2}$ (Akahori et al., 2010). The predicted RM distribution is log-normal, with the power spectrum peaking at a scale of $\sim 1\,h^{-1}\,\mathrm{Mpc}$. Detecting such low RM signals in the IGMF context requires high-sensitivity wide-band radio polarimetric surveys (e.g., SKA, LOFAR), advanced RM-synthesis methods, and careful foreground subtraction (Akahori et al., 2014, O'Sullivan et al., 2018).

Faraday tomography and statistical techniques applied to dense RM grids enable disentangling the IGMF contribution from Galactic, intrinsic, and intervening galaxy RM, especially using high-pass spatial filters and source selection criteria (high redshift, no depolarization, no optical absorption) (Akahori et al., 2014). The observational threshold for detecting IGMF-induced RM is generally $\gtrsim 10$ rad m$^{-2}$ in favorable cases; for single filaments, the expected signal remains at or just below the threshold of current surveys.

4. Constraints from Gamma-Ray Cascades and Time Delay/Echo Effects

Observations of secondary gamma-ray emission produced by cascade processes—both point-like suppression/extension and delayed "pair echo" signatures—impose stringent lower bounds on $B_\mathrm{IGMF}$. The absence of extended, bright GeV emission in Fermi-LAT data around extreme TeV blazars places lower limits: recent analyses with Fermi-LAT and H.E.S.S. data yield $B_\mathrm{IGMF} > 7.1 \times 10^{-16}$ G for $\lambda_B = 1$ Mpc, improving to $\sim 10^{-14}$ G for longer blazar duty cycles ($10^4$–$10^7$ yr) (S. et al., 2023). Conservative assumptions about blazar activity intervals produce somewhat lower bounds, $B_\mathrm{IGMF} \gtrsim 10^{-18}$–$10^{-15}$ G, depending on cascade modeling details and EBL uncertainties (Dermer et al., 2010, Finke et al., 2015, Venters et al., 2012). The result is sensitive to the source's TeV duty cycle, EBL model choice, and the adopted spectral cutoff or break energy.

Detection of delayed GeV "echo" emission from GRBs sets independent IGMF limits. Observations of GRB 221009A with Fermi/LAT and LHAASO established $B_\mathrm{IGMF} \gtrsim 10^{-19}$ G ($\lambda_B > 1$ Mpc), with the limit scaling as $$B_\mathrm{IGMF} > 10^{-19}(l/1\,\mathrm{Mpc})^{-1/2}$$ G for $l < 1$ Mpc (Vovk et al., 2023, Wang et al., 2020, Veres et al., 2017). Similar logic applies to time delays in the cascade from TeV blazars, with the delay relation

$$\Delta t \approx \frac{\lambda_{\gamma\gamma}}{2c} \theta_\mathrm{dfl}^2$$

where $\theta_\mathrm{dfl}$ is the total electron deflection angle determined by $B_\mathrm{IGMF}$ and $\lambda_B$ (Dermer et al., 2010).

5. Detection of Extended Gamma-Ray Halos: Direct IGMF Evidence

Robust direct detection of extended secondary gamma-ray emission, consistent with the cascade halo hypothesis, was reported around the nearby blazar Mkn 501. Analysis of 14 years of Fermi-LAT data finds an extended component with high statistical significance ($\sim 5\sigma$), fitted by an IGMF with $B_\mathrm{rms}=1.5^{+1.6}_{-0.6} \times 10^{-15}$ G and a coherence length $\ell_c=10 \pm 3$ kpc, assuming a source activity timescale exceeding 45,000 years (Webar et al., 15 Sep 2025). The analysis employs the ELMAG Monte Carlo code for cascade simulations, jointly fitting the spectral and angular profiles of the emission to extract the best-fit IGMF parameters. This measurement aligns with the parameter range inferred via the suppression of "cascade bump" GeV emission in more distant blazars and constitutes direct empirical evidence for extremely weak ($\sim10^{-15}$ G) magnetic fields permeating cosmic voids.

6. Astrophysical Implications and Magnetogenesis Constraints

The IGMF’s measured strength and coherence length exclude many astrophysical seeding models unless an unrealistically high magnetized filling factor can be achieved ($f \gtrsim 0.67$), which is difficult via galactic winds or outflows but natural for primordial mechanisms (Tjemsland et al., 2023). The lower bound on $B_\mathrm{IGMF}$ rules out non-helical, weak seed fields arising only from astrophysical origin or passive phase transition turbulence, unless very strong and widespread amplification occurred. The measured IGMF persists to the present epoch and has survived the formation of large-scale structure, providing a unique probe of early universe physics and possible inflationary or baryogenesis processes.

Moreover, the IGMF modifies charged particle propagation (notably ultra-high energy cosmic rays), affects the formation and thermodynamics of cosmic filaments, and influences the anisotropy of the extragalactic gamma-ray background (Venters et al., 2012). Its detection constrains the overall power spectrum and morphology of cosmic magnetism and provides critical input for simulations of structure formation and galaxy cluster evolution.

7. Future Prospects and Methodological Developments

Sensitivity to faint, extended gamma-ray halos or delayed cascade emission will improve with next-generation telescopes such as the Cherenkov Telescope Array (CTA) and with further Fermi-LAT observations (Barai et al., 2018, Keita et al., 19 Sep 2025). GRBs, due to their transient nature, offer sharp temporal windows for echo searches, enabling constraints with minimal contamination from persistent or variable foregrounds. In the radio domain, future SKA RM surveys will push the detection threshold for IGMF-induced Faraday rotation well below 1 rad m$^{-2}$, probing the structure and topology of the field in LSS filaments and voids (Akahori et al., 2014, O'Sullivan et al., 2018). Statistical and tomographic techniques, including high-pass spatial filtering, RM-synthesis, and advanced decomposition methods (e.g., RMCLEAN, QU-fitting), will be required to unambiguously isolate the IGMF signatures from dominant galactic and source-associated foregrounds (Akahori et al., 2014).

A full understanding of the IGMF will require the synthesis of gamma-ray, radio, and possibly CMB probes, along with refined cosmological MHD simulations incorporating different primordial and astrophysical magnetogenesis scenarios. Plasma instability-induced energy loss effects, though generally subdominant in most blazar sightlines, should be considered especially for exceptionally luminous or compact sources (Yan et al., 2018, Webar et al., 15 Sep 2025).

---

In summary, the Intergalactic Magnetic Field is now measured to have a nominal strength of $10^{-15}$ G (femtogauss) with coherence lengths from kpc to Mpc, strongly favoring a primordial, volume-filling origin. Empirical constraints are derived from an array of observational proxies—gamma-ray halos, Faraday rotation measures, cascade suppression, and temporal delays—cementing the IGMF as a key observable linking high-energy astrophysics with the microphysical processes of the early universe, cosmic structure formation, and the evolution of large-scale cosmic magnetism.