Pith. sign in

REVIEW 4 major objections 5 minor 2 cited by

Revision of upper bound on volume-filling intergalactic magnetic fields with LOFAR

T0 review · 4 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The LOFAR Faraday-rotation data that revealed magnetic fields in cosmic filaments also cap the volume-filling intergalactic magnetic field at 70 pG on 1 Mpc scales for scale-invariant spectra, an order of magnitude below the CMB…

desk verdict A genuinely new numerical bound, 70 pG at 1 Mpc for scale-invariant IGMF, but the upper-limit mapping from RRM needs a sign-coherence argument before I'd trust it fully. read the letter →

arxiv 2412.14825 v1 pith:2NWHVTUI submitted 2024-12-19 astro-ph.CO astro-ph.HE

A.Neronov , F.Vazza , S.Mtchedlidze , E.Carretti This is my paper · ORCID
classification astro-ph.COastro-ph.HE
keywords intergalacticmagneticfieldsFaradayrotationmeasureprimordialinflationarymagnetogenesisLOFARlarge-scalestructurecosmicvoids
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that the same LOFAR rotation-measure data that revealed magnetic fields in cosmic filaments also set the tightest known upper limit on the much weaker magnetic field filling the voids of the large-scale structure. The authors treat the measured residual rotation measure (RRM) as an upper bound on the intergalactic contribution, then compare it with magnetohydrodynamic simulations of power-law magnetic spectra to find the maximum allowed field normalization. For a scale-invariant spectrum, the kind inflationary magnetogenesis would produce, the field smoothed over 1 Mpc cannot exceed $B_{\mathrm{Mpc}} = 70\,\mathrm{pG}$. That is an order of magnitude tighter than previous Faraday-rotation limits and an order of magnitude below the CMB anisotropy bound, so radio Faraday data can now probe this class of primordial fields more sharply than CMB data.

What carries the argument

The load-bearing object is the Residual Rotation Measure (RRM), the Faraday rotation left after subtracting the modeled Milky Way contribution and the rotation from massive halos along the line of sight; it is proportional to the line-of-sight integral of the free electron density times the parallel magnetic field component. The argument treats the observed $\mathrm{RRM}(z)$ trend from LoTSS as an upper bound on the IGMF contribution, since the intergalactic field cannot induce more rotation than is actually seen. The comparison side is a set of ENZO cosmological magnetohydrodynamic simulations of $\mathrm{RRM}(z)$ for power-law magnetic spectra $P(k)=A\,(k/k_0)^\alpha$, with the normalization $A$ scaled up until the predicted rotation saturates the observed bound at 90% confidence. Converting the power spectrum to real space gives the field strength $B(\lambda,\alpha)$ averaged over a smoothing scale $\lambda$; the scale-invariant case $\alpha=-3$ becomes a horizontal line in the $B$--$\lambda$ plane, capped at 70 pG.

What would settle it

A decisive check is to simulate the full observation with the intergalactic field set to zero but with all foreground, source-intrinsic, and baryon-fluctuation effects included; if the resulting $\mathrm{RRM}(z)$ distribution matches the LoTSS data as well as the magnetized models do, the $70\,\mathrm{pG}$ cap would not actually constrain the void field. A cheaper observational test is to recompute the bound using only compact, polarization-simple sources and see whether the inferred $B_{\mathrm{Mpc}}$ changes.

Watch

Extended reading notes

Core claim

The central claim is an improved upper bound on the volume-filling intergalactic magnetic field. Using the LoTSS RRM redshift trend from Carretti et al. (2024) after subtracting the Milky Way foreground and massive halo contributions, the paper takes the RRM as an upper limit on the IGMF contribution to Faraday rotation and fits ENZO magnetohydrodynamical simulations for spectra $P(k) = A\,(k/k_0)^\alpha$. For each slope $\alpha$, the maximum allowed normalization $A$ is the value at which the model becomes inconsistent with the data at 90% confidence, converted into a real-space field strength $B(\lambda,\alpha)$. The headline result is that a scale-invariant spectrum ($\alpha = -3$) is capped at $B_{\mathrm{Mpc}} = 70\,\mathrm{pG}$ at the reference smoothing scale $\lambda = 1\,\mathrm{Mpc}$, more than an order of magnitude below the CMB anisotropy limit for such fields; for steeper, causally produced spectra the bound is weaker.

Load-bearing premise

The bound rests on treating the residual rotation measure as an upper limit on the intergalactic magnetic field's contribution, which holds only if no other sign-coherent effect, such as source-intrinsic Faraday rotation, intergalactic baryon fluctuations, or foreground-model error, dominates the residual after the Milky Way and massive halo contributions are subtracted.

Editorial extensions

If this is right

  • A scale-invariant intergalactic field originating from inflation must be no stronger than $70\,\mathrm{pG}$ at $1\,\mathrm{Mpc}$, an order of magnitude below what CMB anisotropy studies allow.
  • If such a scale-invariant field is responsible for easing the Hubble tension through its effect on recombination, it should be detectable through Faraday rotation in the upcoming LOFAR, ASKAP, and SKA data, because the new bound is close to the recombination-based limit.
  • For causally produced fields with steeper spectra ($\alpha > -3$), the CMB anisotropy and clumping constraints remain stronger than the LOFAR bound, so Faraday detection of such fields would favor a non-primordial origin and challenge galactic-outflow or AGN-jet magnetisation scenarios, which require roughly Mpc coherence lengths.
  • The new bound improves the previous Faraday-rotation upper limits by an order of magnitude, making the Faraday-rotation probe competitive with CMB data for a significant range of primordial field configurations.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The 70 pG cap is only as clean as the residual rotation measure: if source-intrinsic Faraday rotation, intergalactic baryon fluctuations, or foreground-model errors contribute sign-coherently to the RRM, the true IGMF upper limit could be higher than quoted; the paper's conservative 'cannot exceed RRM' choice brackets this uncertainty from one side.
  • A direct extension would be to split the LoTSS source sample by polarization morphology and by redshift and test whether the inferred $B_{\mathrm{Mpc}}$ limit is stable; instability would flag contamination rather than a genuine intergalactic signal.
  • If the bound holds and is combined with the gamma-ray lower bounds near $10^{-17}$ G, the allowed window for a scale-invariant void field is narrowed to roughly four orders of magnitude, a range the next generation of Faraday surveys could close entirely.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper uses LoTSS-derived residual rotation measures (RRM) and ENZO cosmological MHD simulations from companion papers by the same group to set upper limits on the normalization of a power-law intergalactic magnetic field (IGMF) power spectrum. For the scale-invariant case (α=-3) it claims an upper limit BMpc ≤ 70 pG at the reference smoothing scale λ = 1 Mpc, which is more than an order of magnitude tighter than previous Faraday-rotation limits and than the CMB-anisotropy bound. The method treats the observed RRM(z) trend as an upper bound on the IGMF contribution, fits simulated RRM predictions to the binned data, and reports 90% confidence upper limits for different spectral slopes and coherence lengths.

Significance. If the upper bound is robust, the result is significant: it exploits an external data set (LoTSS), improves on earlier Faraday constraints by an order of magnitude, and places the scale-invariant inflationary IGMF normalization at a level where it can be confronted with CMB-anisotropy and Hubble-tension-motivated predictions. The paper also makes a concrete, falsifiable statement about detectability with future SKA/ASKAP data. The main caveat is that the core inequality connecting RRM to the IGMF is not fully justified; with that issue addressed, this would be a compact and useful contribution to the IGMF literature.

major comments (4)
  1. [§2] The mapping RRM(z) as an upper bound on the IGMF contribution is load-bearing but not demonstrated. The statement that 'the RM induced by the IGMF cannot exceed the RRM' holds only if the RRM is dominated by a sign-coherent IGMF contribution along each line of sight, or if all contaminating contributions are rigorously subtracted. The paper itself notes that the RRM scatter has 'unclear origin', and after Galactic and massive-halo subtraction the residual includes source-intrinsic Faraday rotation, baryon density fluctuations in the intergalactic medium, and foreground-model errors. If any of these are sign-coherent in the redshift bins used for the fit, the normalization A saturating the 90% confidence bound is an upper limit on the sum IGMF plus contaminants, not on the IGMF alone, and the BMpc ≤ 70 pG claim in §3 is not established. Please add a sign-coherence or contamination test, for example using the sign distribution of RRM per bin, comparing with contaminant-only null simulations, or deriving the limit from the scatter rather than from the binned mean, or state explicitly the additional assumption and its systematic effect.
  2. [§2] The statistical procedure behind the 'best fit' and 'inconsistent with the data at 90% confidence level' is under-specified. The paper does not define the test statistic, the treatment of the bin-to-bin scatter ('wiggles'), the covariance or independence of the RRM redshift bins, or which simulation parameters are held fixed. Without this information the numerical value BMpc = 70 pG cannot be reproduced or checked. Please report the likelihood or chi-square definition, the number of bins, the treatment of the scatter, and the percentile used for the upper limit, or provide a table of the limits for each α in Fig. 1.
  3. [Fig. 1 and §2] The black curve in Fig. 1 is described as an upper bound 'marginalized over the slope α', but it is the envelope of the individual α lines, not a statistical marginalization. If the envelope is intended, the wording should be changed to 'envelope of upper limits'; if a genuine marginalization is intended, the prior on α and the integration procedure need to be given. This distinction matters because the headline scale-invariant bound is read off one particular α, while the 'marginalized' label may overstate the statistical content of the envelope.
  4. [§2] The simulated RRM predictions used to set the limit are imported from Carretti et al. (2024) and Mtchedlidze et al. (2024) without enough detail to assess their systematics. In particular, the box size, resolution, magnetic field evolution, and the method by which simulated RRM is extracted and scaled to the observed RRM are not summarized. Since the upper limit is set by comparing these simulations to the data, please include a concise summary of these parameters or explicitly state which values from the companion papers are used.
minor comments (5)
  1. [Abstract] Typo: 'to to derive' should read 'to derive'.
  2. [§2] Typo: 'pre-recobmination' should read 'pre-recombination'.
  3. [§2] The sentence describing the RRM scatter has an unclosed parenthesis: '...likely with a physical origin.' needs a closing bracket after 'origin'.
  4. [Fig. 2 caption] The caption contains a stray '[?]' in 'generated during EW or QCD phase transitions [?]'; this should be replaced with proper citations.
  5. [Fig. 2 and §2] The text refers to the 'red dashed line' as a slight improvement, but elsewhere the same bound is called a 'red dotted line' and a 'red solid curve'; the color/line-style labels should be made consistent.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found: the LOFAR upper bound is obtained by comparing forward-modeled simulated rotation measures to external LoTSS-derived RRM data, not by defining or fitting the conclusion into the inputs.

full rationale

The derivation chain is: (i) RRM(z) is taken from Carretti et al. (2024), which is based on LoTSS data (O'Sullivan et al. 2023) after subtracting the Milky Way foreground and massive-halo contributions; (ii) RRM predictions are generated from ENZO cosmological MHD simulations of IGMF models with power spectra P(k) = A(k/k0)^alpha (Carretti et al. 2024; Mtchedlidze et al. 2024); (iii) for each simulated model the normalization A is scaled until the predicted RRM becomes inconsistent with the observed RRM at 90% confidence; (iv) the resulting maximal A is converted to the real-space amplitude BMpc at 1 Mpc. None of these steps defines the target bound from itself: the simulated RRM is a forward prediction of the model, and the LoTSS-derived RRM is external to the simulation. The step 'we consider the RRM(z) trend ... as an upper bound on the contribution from IGMF, noticing that the RM induced by the IGMF cannot exceed the RRM' is a physical interpretation of the residual data, not a circular reduction; it is a conservative upper bound only if non-IGMF residuals are sign-coherent, and the paper itself flags the 'significant scatter of unclear origin,' so this is a robustness caveat rather than a circularity. The paper also notes that the 'constrained RRM' comes from Carretti et al. (2024) and the simulations from Mtchedlidze et al. (2024), both with overlapping authorship; however, these are used as data products and model templates with external LoTSS anchoring, not as an appeal to authority for the conclusion. The final claim BMpc <= 70 pG is an upper limit computed from the data and models, not a quantity already contained in the inputs by construction. No fitted parameter is relabeled as a prediction.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The central bound rests on the assumed mapping from RRM to IGMF (which relies on simulations of the baryon distribution) and on the input spectral slopes and damping scale. No new particles or forces are introduced. The only fitted quantity is the magnetic field normalization, which is the target parameter of an upper-limit measurement.

free parameters (2)
  • Magnetic field normalization BMpc at 1 Mpc = ≤70 pG for alpha = -3; slope-dependent for other alpha
    Normalization of the magnetic field power spectrum P(k) = A(k/k0)^alpha. The upper limit is found by fitting simulated RRM to the LoTSS constrained RRM data and finding the value where the model becomes inconsistent at 90% CL.
  • Spectral slope alpha = -3 to 2
    The spectral index of the magnetic field power spectrum is scanned over a range of values, chosen by hand to cover scale-invariant inflationary fields (alpha=-3) to causal phase-transition fields (alpha=2).
assumptions (4)
  • domain assumption The RRM(z) trend from Carretti et al. (2024) is an upper bound on the IGMF contribution to Faraday rotation.
    Section 2 states this directly. If unmodeled contributions contaminate RRM, the bound could be biased; the paper acknowledges scatter of unclear origin.
  • domain assumption ENZO cosmological MHD simulations accurately predict RRM as a function of magnetic field strength and free electron density.
    Section 2 uses the simulated RRM models from Carretti et al. (2024) and Mtchedlidze et al. (2024) to map P(k) normalization to observable RM. The mapping depends on the simulated baryon distribution in voids.
  • domain assumption Small-scale damping scale lD lies in the 10-100 kpc range.
    Section 2 states this is 'usually assumed' and does not affect the relevant scales; it limits extrapolation of the (B,lambda) relation.
  • standard math The real-space field strength is related to the power spectrum by B = BMpc (lambda/1 Mpc)^(-(alpha+3)/2).
    Section 2 uses this to convert the power spectrum normalization to field strength at a given smoothing scale. This follows from the definitions of P(k) and the spectral index convention.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Revision of upper bound on volume-filling intergalactic magnetic fields with LOFAR." pith.science (2026). https://pith.science/paper/2NWHVTUI

@misc{pith2026241214825,
  author       = {Pith},
  title        = {Pith review of: Revision of upper bound on volume-filling intergalactic magnetic fields with LOFAR},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2NWHVTUI}},
  note         = {Machine review of arXiv:2412.14825}
}
read the original abstract

Magnetic fields present in the Large Scale Structure (LSS) of the Universe change polarization of radio waves arriving from distant extragalactic sources through the effect of Faraday rotation. This effect has been recently used to detect magnetic field in the LSS filaments based on the Rotation Measure data of the LOFAR Two-Meter Sky Survey (LoTSS). We notice that the same data also constrain the strength of the volume-filling magnetic field in the voids of the LSS. We use the LoTSS data to to derive an improved upper bound on the volume-filling field. The new upper bound provides an order of magnitude improvement on the previous Faraday rotation bounds. The new Faraday Rotation bound on the scale-invariant field that may originate from the epoch of inflation is also an order of magnitude lower than the bound on such field derived from the anisotropy analysis of the Cosmic Microwave Background.

Figures

Figures reproduced from arXiv: 2412.14825 by the authors.

Figure 1
Figure 1. shows the upper bounds on the field strength as a function of distance scale for different assumptions about the slope of the power spectrum α and the overall upper bound marginalized over α. The magnetic field energy density is re￾lated to the power spectrum as EM = R 4πk 2P(k)dk. The mini￾mal possible value of α = −3 corresponds to the scale-invariant magnetic field (which only inflationary processes can produce) … view at source ↗
Figure 2
Figure 2. LOFAR RM bounds on IGMF strength and correlation length, compared with other known constraints. Green area shows possible lo￾cations of endpoints of the evolution of PMF Hosking & Schekochihin (2023); Brandenburg et al. (2024); Banerjee & Jedamzik (2004) gener￾ated during EW or QCD phase transitions [?].. The red semi-transparent arrow shows the level of other recent limits from the Faraday rota￾tion, obtained by co… view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Radio Observations as a Probe of Cosmic Web Magnetism

    astro-ph.CO 2025-05 conditional novelty 6.0 of 10

    Radio observations of cosmic filaments favor a dominant primordial magnetic field over an astrophysical-only origin, with best-fit filament field 43 ± 7 nG at z=0.

  2. Imprints of primordial magnetic fields in gravitational collapse during early structure formation

    astro-ph.CO 2026-02 conditional novelty 5.0 of 10

    Gravitational collapse can trigger a small-scale dynamo that amplifies magnetic fields below the Jeans scale, potentially erasing primordial spectral features unless the turbulent inertial range is resolved.

Reference graph

Works this paper leans on

38 extracted references · 24 canonical work pages · cited by 2 Pith papers

  1. [1]

    Acciari, V . A. et al. 2023, Astron. Astrophys., 670, A145

    work page 2023

  2. [2]

    Ade, P. A. R. et al. 2016, Astron. Astrophys., 594, A19

    2016

  3. [3]

    Aharonian, F. et al. 2023, Astrophys. J. Lett., 950, L16

    2023

  4. [4]

    2022, Mon

    Aramburo-Garcia, A., Bondarenko, K., Boyarsky, A., et al. 2022, Mon. Not. Roy. Astron. Soc., 515, 5673

    work page 2022

  5. [5]

    & Jedamzik, K

    Banerjee, R. & Jedamzik, K. 2004, Phys. Rev. D, 70, 123003

    2004

  6. [6]

    2006, Mon

    Bertone, S., V ogt, C., & Ensslin, T. 2006, Mon. Not. Roy. Astron. Soc., 370, 319

    work page 2006

  7. [7]

    Blasi, P., Burles, S., & Olinto, A. V . 1999, Astrophys. J. Lett., 514, L79

    1999

  8. [8]

    & Neronov, A

    Blunier, J. & Neronov, A. 2024, Astron. Astrophys., 691, A34

    work page 2024

Show all 38 references
  1. [9]

    2022, Astron

    Bondarenko, K., Boyarsky, A., Korochkin, A., et al. 2022, Astron. Astrophys., 660, A80

    2022

  2. [10]

    2017, Phys

    Brandenburg, A., Kahniashvili, T., Mandal, S., et al. 2017, Phys. Rev. D, 96, 123528

    2017

  3. [11]

    2024, Astron

    Brandenburg, A., Neronov, A., & Vazza, F. 2024, Astron. Astrophys., 687, A186

    2024

  4. [12]

    P., Vacca, V ., et al

    Carretti, E., O’Sullivan, S. P., Vacca, V ., et al. 2023, MNRAS, 518, 2273

    2023

  5. [13]

    P., et al

    Carretti, E., Vazza, F., O’Sullivan, S. P., et al. 2024 [arXiv:2411.13499]

    2024 arXiv

  6. [14]

    2009, Mon

    Donnert, J., Dolag, K., Lesch, H., & Müller, E. 2009, Mon. Not. Roy. Astron. Soc., 392, 1008

    2009

  7. [15]

    & Neronov, A

    Durrer, R. & Neronov, A. 2013, Astron. Astrophys. Rev., 21, 62

    2013

  8. [16]

    Hogan, C. J. 1983, Phys. Rev. Lett., 51, 1488

    1983

  9. [17]

    Hosking, D. N. & Schekochihin, A. A. 2023, Nature Commun., 14, 7523

    2023

  10. [18]

    2023, arXiv e-prints, arXiv:2312.11448

    Jedamzik, K., Abel, T., & Ali-Haimoud, Y . 2023, arXiv e-prints, arXiv:2312.11448

    2023 arXiv

  11. [19]

    & Pogosian, L

    Jedamzik, K. & Pogosian, L. 2020, Phys. Rev. Lett., 125, 181302

    2020

  12. [20]

    & Saveliev, A

    Jedamzik, K. & Saveliev, A. 2019, Phys. Rev. Lett., 123, 021301

    2019

  13. [21]

    Marinacci, F. et al. 2018, Mon. Not. Roy. Astron. Soc., 480, 5113

    2018

  14. [22]

    2024, The Astrophysi- cal Journal, 977, 128

    Mtchedlidze, S., Domínguez-Fernández, P., Du, X., et al. 2024, The Astrophysi- cal Journal, 977, 128

    2024

  15. [23]

    2023, Phys

    Neronov, A., Semikoz, D., & Kalashev, O. 2023, Phys. Rev. D, 108, 103008

    2023

  16. [24]

    & V ovk, I

    Neronov, A. & V ovk, I. 2010, Science, 328, 73 O’Sullivan, S. P., Shimwell, T. W., Hardcastle, M. J., et al. 2023, MNRAS, 519, 5723

    2010

  17. [25]

    Paoletti, D., Chluba, J., Finelli, F., & Rubiño-Martín, J. A. 2019, MNRAS, 484, 185

    2019

  18. [26]

    & Finelli, F

    Paoletti, D. & Finelli, F. 2019, JCAP, 2019, 028

    2019

  19. [27]

    S., Tinyakov, P

    Pshirkov, M. S., Tinyakov, P. G., & Urban, F. R. 2016, Phys. Rev. Lett., 116, 191302

    2016

  20. [28]

    M., Loeb, A., & Spergel, D

    Quashnock, J. M., Loeb, A., & Spergel, D. N. 1989, AP.J.Lett, 344, L49

    1989

  21. [29]

    1992, AP.J.Lett, 391, L1

    Ratra, B. 1992, AP.J.Lett, 391, L1

    1992

  22. [30]

    2016, Reports on Progress in Physics, 79, 076901

    Subramanian, K. 2016, Reports on Progress in Physics, 79, 076901

    2016

  23. [31]

    Tajima, T., Cable, S., Shibata, K., & Kulsrud, R. M. 1992, ApJ, 390, 309

    1992

  24. [32]

    G., Kisslinger, L., Brandenburg, A., & Kahniashvili, T

    Tevzadze, A. G., Kisslinger, L., Brandenburg, A., & Kahniashvili, T. 2012, ApJ, 759, 54

    2012

  25. [33]

    2018, MNRAS, 481, 3401

    Trivedi, P., Reppin, J., Chluba, J., & Banerjee, R. 2018, MNRAS, 481, 3401

    2018

  26. [34]

    Turner, M. S. & Widrow, L. M. 1988, Phys. Rev. D, 37, 2743

    1988

  27. [35]

    2021, Rept

    Vachaspati, T. 2021, Rept. Prog. Phys., 84, 074901

    2021

  28. [36]

    2017, Classical and Quantum Gravity V ovk, I., Korochkin, A., Neronov, A., & Semikoz, D

    Vazza, F., Brueggen, M., Gheller, C., et al. 2017, Classical and Quantum Gravity V ovk, I., Korochkin, A., Neronov, A., & Semikoz, D. 2024, Astron. Astrophys., 683, A25

    2017

  29. [37]

    C., Li, S., & Norman, M

    Xu, H., Li, H., Collins, D. C., Li, S., & Norman, M. L. 2011, ApJ, 739, 77

    2011

  30. [38]

    2017, Phys

    Zucca, A., Li, Y ., & Pogosian, L. 2017, Phys. Rev. D, 95, 063506 Article number, page 4 of 4

    2017

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.