pith. sign in

arxiv: 2606.25630 · v1 · pith:GIQ26ZQUnew · submitted 2026-06-24 · 🌌 astro-ph.IM

Application of Bayesian Statistical Tools to SKA Telescopes Polarization Surveys to Study Magnetization of the Large-scale Structure of the Universe

Valentina Vacca , Sebastian Hutschenreuter , Andrea Cabriolu , Torsten A. Ensslin , Philipp Frank , Jakob Roth , Martin Reinecke , Jens Jasche
show 13 more authors
Florent Leclercq Ettore Carretti Luigina Feretti Chiara Ferrari Federica Govoni Cathy Horellou Melanie Johnston-Hollitt Francesca Loi Matteo Murgia Shane O'Sullivan Rosita Paladino Tessa Vernstrom Gianni Fenu
This is my paper

Pith reviewed 2026-06-25 20:23 UTC · model grok-4.3

classification 🌌 astro-ph.IM
keywords Faraday rotationBayesian statisticspolarization surveyscosmic webgalaxy clustersfilamentsmagnetizationlarge-scale structure
0
0 comments X

The pith

Bayesian analysis of Faraday rotation data shows 50,000 mid-frequency measurements with high-precision redshifts are needed to constrain magnetization in galaxy clusters.

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

The paper develops a Bayesian statistical method to extract information on magnetic fields from Faraday rotation observations of background radio sources across different parts of the cosmic web. It calculates the sample sizes required from planned polarization surveys to place limits on magnetization levels in dense versus sparse environments. Roughly 50,000 mid-frequency measurements paired with accurate redshifts would allow constraints on fields inside galaxy clusters, while low-frequency data and spectroscopic redshifts for at least 17,000 sources could yield first limits on weaker fields along filaments. These numbers matter because they indicate what is realistically achievable for mapping how magnetic fields are distributed on the largest scales.

Core claim

A Bayesian framework applied to large samples of Faraday rotation measures can separate contributions from galaxy clusters, filaments, and voids, with the result that 50,000 mid-frequency measurements and high-precision redshifts suffice to constrain magnetization in dense environments while 17,000 spectroscopic redshifts at low frequencies can provide initial constraints on filaments.

What carries the argument

Bayesian statistical approach to separating Faraday rotation contributions from different large-scale structure environments

If this is right

  • Magnetization inside galaxy clusters becomes measurable once 50,000 mid-frequency Faraday rotation values are combined with precise redshifts.
  • Initial limits on magnetization along filaments become feasible with low-frequency observations and spectroscopic redshifts for at least 17,000 sources.
  • Magnetic fields in sheets and voids stay largely unconstrained even with the largest planned samples.
  • Redshift accuracy directly determines whether environmental signals can be isolated without mixing.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Future optical or spectroscopic surveys will need to be coordinated with radio observations to meet the redshift requirements.
  • The remaining difficulty for filaments suggests that stacking or cross-correlation techniques could be tested as ways to boost sensitivity.
  • If the separation works, the same data set could be reanalyzed for other statistical properties of the magnetic field distribution.

Load-bearing premise

The spatial distribution of background radio sources, the coherence scales of magnetic fields, and the accuracy of redshift information will allow the Bayesian framework to cleanly separate contributions from clusters, filaments, and voids without significant contamination or model mismatch.

What would settle it

If real survey data produce magnetization estimates for clusters or filaments that disagree with independent probes such as X-ray or synchrotron observations at a level exceeding the expected statistical uncertainties, the required sample sizes would need revision.

Figures

Figures reproduced from arXiv: 2606.25630 by Andrea Cabriolu, Cathy Horellou, Chiara Ferrari, Ettore Carretti, Federica Govoni, Florent Leclercq, Francesca Loi, Gianni Fenu, Jakob Roth, Jens Jasche, Luigina Feretti, Martin Reinecke, Matteo Murgia, Melanie Johnston-Hollitt, Philipp Frank, Rosita Paladino, Sebastian Hutschenreuter, Shane O'Sullivan, Tessa Vernstrom, Torsten A. Ensslin, Valentina Vacca.

Figure 1
Figure 1. Figure 1: Synthetic Galactic Faraday sky. Concerning the Galactic term, we produced synthetic rotation measure values by multiplying the dispersion measure image by Hutschenreuter et al. (2024) by a Gaussian random magnetic field with zero-mean and standard deviation equal to one. This magnetic field has been generated by using the following power spectrum |𝐵ℓ | 2 as a function of wave number ℓ, |𝐵ℓ | 2 = 𝑃0 1.0 + … view at source ↗
Figure 2
Figure 2. Figure 2: Top panel: histogram of the uncertainties in the Faraday rotation measurements from the LoTSS DR2 RM catalogue (O’Sullivan et al., 2023) in red colors and from the MeerKAT data (Loi et al., 2025) in blue colors. Bottom panel: histogram of the uncertainties in the Faraday rotation measurements from the MeerKAT data (Loi et al., 2025) at 0-9-1.4 GHz in blue colors with no ionospheric corrections and in black… view at source ↗
Figure 3
Figure 3. Figure 3: Results obtained assuming an extragalactic Faraday rotation standard deviation of ≈10 rad/m2 and using a rotation measure catalogue based on observations at mid-frequencies. Top panels: Reconstructed mean and uncertainty of the Galactic Faraday sky. Bottom panels 2-dimensional marginalization (in colours) and 1-dimensional marginalization (grey histograms) of the posterior distribution of the extragalactic… view at source ↗
Figure 4
Figure 4. Figure 4: As in [PITH_FULL_IMAGE:figures/full_fig_p017_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Left panel: distribution of observed photometric (blue) and spectroscopic (red) redshifts. The purple colour shows the overlapping region. Right panel: distribution of the luminosity of radio sources. Data have been taken from the LoTSS DR2 RM catalogue by O’Sullivan et al. (2023), see text for more details. 18 [PITH_FULL_IMAGE:figures/full_fig_p018_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Uncertainty in Faraday rotation standard deviation as a function of redshift (left) and luminosity (right) of radio sources for spectroscopic (red) and photometric (blue) redshifts. On top panels, the uncertainty in Faraday rotation has been estimated considering an overall extragalactic Faraday rotation standard deviation of about 1.5 rad/m2 while on bottom panels of about 10 rad/m2 . Redhsifts and lumino… view at source ↗
read the original abstract

Understanding cosmological magnetic fields requires a detailed knowledge of magnetism in the different environments of the large-scale structure of the Universe. Magnetic fields are well known to inhabit galaxy clusters, and recently their presence has been detected between galaxy clusters, along filaments extending up to 10-15 Mpc. Beyond that, there is limited information on the existence of magnetic fields in sheets and voids of the cosmic web. We propose a Bayesian statistical approach to study magnetic fields on large scales through observations of the Faraday rotation effect in large samples of polarized point-like background radio sources. We present the expectations to detect magnetization in environments of the large-scale structure with the SKA-Mid polarization survey planned by the SKAO Magnetism Science Working Group and with SKA-Low with AA4 telescopes, and discuss the required level of accuracy on the redshifts of the host galaxies for such a study. We find that about 50,000 mid-frequency Faraday rotation measurements complemented by high-precision redshifts are needed to constrain magnetization of dense environments as galaxy clusters. Investigation of magnetization in weakly-magnetized low-density enviroments, as filaments, will remain challenging, but low frequencies radio observations and spectroscopic redhifts for at least 17,000 will allow us to put first constraints.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The manuscript proposes a Bayesian hierarchical statistical framework for analyzing Faraday rotation measures (RMs) from polarized background sources in planned SKA-Mid and SKA-Low surveys. It derives numerical forecasts showing that approximately 50,000 mid-frequency RM measurements paired with high-precision redshifts are required to constrain magnetization in dense environments such as galaxy clusters, while low-frequency observations combined with at least 17,000 spectroscopic redshifts could yield initial constraints on weakly magnetized filaments.

Significance. If the underlying simulation pipeline, priors on coherence lengths, and source-density assumptions hold, the work supplies concrete observational targets that can directly inform SKAO Magnetism Science Working Group survey design and data-analysis strategies. The explicit forward-modeling of contributions from clusters, filaments, and voids, together with the hierarchical Bayesian separation, represents a strength that allows falsifiable predictions for required sample sizes.

minor comments (2)
  1. Abstract: the text contains typographical errors ('enviroments', 'redhifts') that should be corrected prior to publication.
  2. The presentation of the numerical thresholds (50,000 and 17,000) would be strengthened by an explicit statement, perhaps in the methods or results section, of how variations in the assumed spatial distribution of sources affect the final counts.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their supportive review, positive assessment of the significance of our Bayesian hierarchical framework for SKA RM surveys, and recommendation for minor revision. No specific major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity; forward-modeling forecasts are self-contained

full rationale

The manuscript performs forward-modeling of expected Faraday RM contributions from clusters, filaments and voids under an explicit Bayesian hierarchical model. The quoted thresholds (50 000 mid-frequency RMs with high-precision z; 17 000 spectroscopic redshifts at low frequency) are numerical outputs of that pipeline given stated priors on coherence lengths and source densities. No equation reduces a claimed prediction to a fitted parameter from the same dataset, no load-bearing self-citation chain is invoked, and the derivation does not rename known results or smuggle ansatzes. The work is therefore internally consistent and externally falsifiable against independent RM surveys and simulations.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The estimates rest on standard domain assumptions about Faraday rotation physics and source populations rather than new free parameters or invented entities.

axioms (1)
  • domain assumption Magnetic field properties and source distributions in clusters, filaments, and voids follow models that permit statistical separation via Bayesian inference.
    The numerical thresholds depend on these prior models of magnetization and radio source statistics.

pith-pipeline@v0.9.1-grok · 5852 in / 1201 out tokens · 22876 ms · 2026-06-25T20:23:55.725802+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

29 extracted references · 28 canonical work pages · 1 internal anchor

  1. [1]

    doi: 10.1093/mnras/133.1.67. E. Carretti et al.MNRAS, 512(1):945–959, May

    work page doi:10.1093/mnras/133.1.67

  2. [2]

    doi: 10.1093/mnras/stac384. E. Carretti et al.MNRAS, 518(2):2273–2286, Jan

    work page doi:10.1093/mnras/stac384

  3. [3]

    doi: 10.1093/mnras/stac2966. E. Carretti et al.A&A, 693:A208, Jan

    work page doi:10.1093/mnras/stac2966

  4. [4]

    doi: 10.1051/0004-6361/202451333. T. E. Clarke.Journal of Korean Astronomical Society, 37(5):337–342, Dec

    work page doi:10.1051/0004-6361/202451333

  5. [5]

    doi: 10.5303/ JKAS.2004.37.5.337. T. A. Enßlin, M. Frommert, and F. S. Kitaura.Physical Review D, 80(10):105005, Nov

    2004

  6. [6]

    doi: 10.1103/PhysRevD.80.105005. P. Frank, R. Leike, and T. A. Enßlin.Entropy, 23(7):853, July

    work page doi:10.1103/physrevd.80.105005

  7. [7]

    GeometricVariationalInference.Entropy, 23(7):853, July 2021

    doi: 10.3390/e23070853. F.GovoniandL.Feretti.InternationalJournalofModernPhysicsD,13(8):1549–1594,Jan.2004. doi: 10.1142/S0218271804005080. F. Govoni et al.Science, 364(6444):981–984, June

    work page doi:10.3390/e23070853 2004

  8. [8]

    doi: 10.1126/science.aat7500. A. M. Hammond, T. Robishaw, and B. M. Gaensler.arXiv e-prints, art. arXiv:1209.1438, Sept

    work page doi:10.1126/science.aat7500

  9. [9]

    doi: 10.48550/arXiv.1209.1438. G. Heald et al.Galaxies, 8(3):53, July

    work page doi:10.48550/arxiv.1209.1438

  10. [10]

    doi: 10.3390/galaxies8030053. B. Hugo and R. Perley. art. SSA-0004E-001, Jan

    work page doi:10.3390/galaxies8030053

  11. [11]

    doi: 10.1051/0004-6361/ 201935479. S. Hutschenreuter et al.A&A, 657:A43, Jan

    work page doi:10.1051/0004-6361/

  12. [12]

    doi: 10.1051/0004-6361/202140486. S. Hutschenreuter et al.A&A, 690:A314, Oct

    work page doi:10.1051/0004-6361/202140486

  13. [13]

    doi: 10.1051/0004-6361/202346740. P. P. Kronberg et al.ApJ, 676(1):70–79, Mar

    work page doi:10.1051/0004-6361/202346740

  14. [14]

    J.Linetal.RASTechniquesandInstruments,3(1):737–747,Jan.2024

    doi: 10.1086/527281. J.Linetal.RASTechniquesandInstruments,3(1):737–747,Jan.2024. doi: 10.1093/rasti/rzae051. F. Loi et al.A&A, 694:A125, Feb

    work page doi:10.1086/527281 2024

  15. [15]

    doi: 10.1051/0004-6361/202451711. V. Mainieri et al.arXiv e-prints, art. arXiv:2403.05398, Mar

    work page doi:10.1051/0004-6361/202451711

  16. [16]

    doi: 10.48550/arXiv.2403. 05398. N. Mauri et al. InJournal of Physics Conference Series, volume 1342 ofJournal of Physics Conference Series, page 012122. IOP, Jan

    work page doi:10.48550/arxiv.2403

  17. [17]

    doi: 10.1088/1742-6596/1342/1/012122. J. A. Newman and D. Gruen.ARA&A, 60:363–414, Aug

    work page doi:10.1088/1742-6596/1342/1/012122

  18. [18]

    doi: 10.1051/0004-6361/201423995. S. P. O’Sullivan et al.MNRAS, 519(4):5723–5742, Mar

    work page doi:10.1051/0004-6361/201423995

  19. [19]

    S.P.O’Sullivanetal

    doi: 10.1093/mnras/stac3820. S.P.O’Sullivanetal. InAdvancingAstrophysicswiththeSKA–II(AASKAII).2026. arXivsearch: Report number AASKAII/OSullivan01. S. Piras et al.A&A, 687:A267, July

    work page doi:10.1093/mnras/stac3820 2026

  20. [20]

    doi: 10.1051/0004-6361/202349085. L. Rudnick and F. N. Owen.ApJ, 785(1):45, Apr

    work page doi:10.1051/0004-6361/202349085

  21. [21]

    doi: 10.1088/0004-637X/785/1/45. D. H. F. M. Schnitzeler.MNRAS, 409(1):L99–L103, Nov

    work page doi:10.1088/0004-637x/785/1/45

  22. [22]

    doi: 10.1111/j.1745-3933.2010. 00957.x. J. M. Stil.ApJ, 987(2):173, July

    work page doi:10.1111/j.1745-3933.2010 2010

  23. [23]

    doi: 10.3847/1538-4357/add69e. J. M. Stil, B. W. Keller, S. J. George, and A. R. Taylor.ApJ, 787(2):99, June

    work page doi:10.3847/1538-4357/add69e

  24. [24]

    doi: 10.3847/1538-3881/ad4150. A. R. Taylor, J. M. Stil, and C. Sunstrum.ApJ, 702(2):1230–1236, Sept

    work page doi:10.3847/1538-3881/ad4150

  25. [25]

    doi: 10.3847/1538-4357/ ad22dd. V. Vacca et al. InAdvancing Astrophysics with the Square Kilometre Array (AASKA14), page 114, Apr

    work page doi:10.3847/1538-4357/

  26. [26]

    doi: 10.22323/1.215.0114. V. Vacca et al.A&A, 591:A13, June

    work page doi:10.22323/1.215.0114

  27. [27]

    DEFROST: Detecting Excess in Faraday Rotation thrOugh Sophisticated analysis Techniques

    doi: 10.1051/0004-6361/201527291. V.Vaccaetal.arXive-prints,art.arXiv:2605.13605,May2026. doi: 10.48550/arXiv.2605.13605. 14 Application of Bayesian statistical tools to SKA telescopes polarization surveys Vacca et al. C. L. Van Eck et al.ApJSS, 267(2):28, Aug

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/201527291

  28. [28]

    doi: 10.3847/1538-4365/acda24. T. Vernstrom et al.Science Advances, 9(7):eade7233, Feb

    work page doi:10.3847/1538-4365/acda24

  29. [29]

    Science Advances , keywords =

    doi: 10.1126/sciadv.ade7233. 15 Application of Bayesian statistical tools to SKA telescopes polarization surveys Vacca et al. Figure3:ResultsobtainedassuminganextragalacticFaradayrotationstandarddeviationof≈10rad/m 2 and using a rotation measure catalogue based on observations at mid-frequencies.Top panels:Reconstructed meananduncertaintyoftheGalacticFara...

    work page doi:10.1126/sciadv.ade7233 2023