{"id":"191528b1-664b-4e5b-83ba-4d9ac2dc8faf","arxiv_id":"2505.18619","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":5,"one_line_summary":"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.","lead":"This paper compares radio observations of cosmic filaments with new computer simulations to test where the Universe's magnetic fields came from. It concludes that the data favor magnetic fields born in the early Universe, with an average filament field of a few tens of nanoGauss, rather than fields created mostly by galaxies.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline 43±7 nG and the 'rapidly decreasing astrophysical component' rest on the ad hoc γ=5 shape in Eq. (9); the LOFAR data alone do not constrain γ, and γ=1 gives B_f,0≈14 nG.","rationale":"The paper makes a qualitative claim (primordial magnetic fields dominate; astrophysical-only disfavored) and a quantitative claim (best combined model gives 43±7 nG with a rapidly decreasing astrophysical component). The qualitative claim has independent support: the astro B4 simulation falls short of the high-redshift RRM data, and the synchrotron stacking provides a separate observable that the astrophysical-only model fails to match. The quantitative claim, however, is set by an assumed functional form for the astrophysical RM component. The only shapes tested are γ=1 and γ=5 in Eq. (9), with γ=5 preferred because it matches the integrated 21±4% astrophysical fraction; that constraint does not determine the redshift evolution of the component. The fitted B_f,0 varies from ~14 nG to ~43–60 nG across the two shapes, so the headline value is not pinned down by the data. This is the single most load-bearing concern: if the true astrophysical component evolves differently (e.g., rises toward z~1–2 as expected from AGN activity), both the 'rapidly decreasing' claim and the 43±7 nG value would be wrong. A secondary circularity is that the nB=-1 seed amplitude was calibrated on the same LOFAR RRM data used for model selection (Section 3, Section 4.1); this affects the 'nB=-1 favored' claim but not the primordial-vs-astrophysical dichotomy. The paper itself notes that independent knowledge of the astrophysical component shape is 'much-needed' (Section 6), which is an implicit admission of the limitation. The reader's conditional accept is appropriate; no change in verdict is needed. The proposed free-γ refit would settle whether the data actually constrain γ or merely tolerate it.","tokens_in":28199,"tokens_out":10649,"duration_ms":95465,"concrete_test":"Re-fit Eq. (9) to the same LoTSS RRM rms data (Section 4.1) with γ as a free parameter, or with a physically motivated astrophysical shape (e.g., ∝ star-formation-rate density or AGN luminosity density), and compute the profile likelihood of B_f,0. If the 68% interval for B_f,0 spans both ~14 nG (γ=1) and ~43 nG (γ=5), or the likelihood is flat in γ, then the headline 43±7 nG is not supported and the model-independent 10–60 nG range should be the primary result. If the data strongly prefer γ≈5 and exclude γ≈1, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative headline (B_f,0 = 43±7 nG with an astrophysical RRM component rapidly decreasing with z) is produced by fitting Eq. (9) with the astrophysical component fixed to A_rrm/(1+z)^γ and γ=5. Only two shapes are tested, γ=1 and γ=5; γ=1, which the authors' previous work preferred, yields B_f,0 = 11–14±4 nG, while γ=5 yields 43–60±13 nG. The paper chooses γ=5 partly because it reproduces the observed 21±4% astrophysical fraction, but that fraction is an integrated constraint and does not fix the redshift dependence of the astrophysical term. Because the astrophysical and filament terms are both power-laws in (1+z) added in quadrature, the RRM(z) data alone are degenerate in γ. If the true astrophysical component declines more slowly (or rises) with z, the inferred filament field shifts by roughly a factor of four, and the 'best-matching combined model' as well as the 'rapidly decreasing' claim are not robust. The qualitative conclusion that a solely astrophysical simulation fails at high z is less affected, since the astro B4 model falls short for both tested shapes, but the paper's distinctive quantitative result and its model selection are conditional on an unconstrained functional form. The abstract's broad 10–60 nG range is honest, but it is not the claim highlighted in the conclusions, and the paper itself calls independent knowledge of this shape 'much-needed' (Section 6).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compares three radio observables of cosmic filaments—LOFAR RRM rms versus redshift, stacked synchrotron emission, and POSSUM MAD RRM radial profiles around galaxy groups—with predictions from cosmological MHD simulations of astrophysical and primordial magnetic seeding scenarios. Using an analytical fit with a power-law filament field and an ad hoc astrophysical RM component, the authors find that a combined primordial (nB = -1) plus astrophysical model best matches the first two probes, yielding B_f,0 = 43 ± 7 nG at z = 0 for a rapidly decreasing astrophysical component (γ = 5), while a solely astrophysical model fails at high redshift. The third probe is inconclusive, and the paper recommends independent measurement of the astrophysical RM shape.","tokens_in":28542,"tokens_out":8086,"duration_ms":64000,"significance":"If the main inference holds, radio observations of the cosmic web would provide competitive constraints on primordial magnetic fields, potentially more sensitive than CMB studies for inflationary-like spectra, and would support a dominant primordial component in filament magnetization. The paper's strengths include the use of recent high-quality LOFAR data, an updated astrophysical simulation (B4) calibrated to multiple galaxy observables, careful treatment of line-of-sight selection and high-density flagging, and an honest discussion of the third probe's inconclusiveness. The data availability statement provides access to the simulated LOS data. However, the quantitative headline and the model-selection step are conditional on an unconstrained shape of the astrophysical RM component and on a seed normalization calibrated to the same LOFAR data, so the specific 43 ± 7 nG claim is currently not as robust as the abstract suggests.","major_comments":[{"comment":"The headline result B_f,0 = 43 ± 7 nG (abstract; §6) is obtained only for γ = 5, one of two arbitrarily chosen shapes for the astrophysical component. The data do not constrain γ: with γ = 1, the same fits give B_f,0 = 11–14 ± 4 nG (Table 2). The choice of γ = 5 is justified partly by matching the observed 21 ± 4% astrophysical fraction, but that integrated fraction does not fix the redshift dependence of the astrophysical term; since both terms in Eq. (9) are power laws in (1+z) added in quadrature, the RRM(z) data alone are degenerate in γ. Thus the abstract's 'best-matching combined model' and the 'rapidly decreasing astrophysical component' are not robust. The paper's own §6 calls independent knowledge of the astrophysical RM shape 'much-needed,' which correctly identifies the load-bearing uncertainty.","section":"§4.1, Eq. (9), Tables 2–3; §6"},{"comment":"The nB = -1 seed amplitude, <B>_1Mpc = 0.37 nG, is set from comparison with LOFAR RRMs (Section 3, citing Ref. [18]), not from CMB constraints. The same LOFAR RRM-z data are then used in Section 4.1 to fit Eq. (9) and in Figures 4–5 to identify this model as the unique match. This is a circular validation: the model is normalized to the very data set used to test it, so the statement in Section 5 that this model is 'the only scenario that is consistent with both the z ≤ 3 cosmic web and CMB limits' is true partly by construction. A non-circular test would require leaving the seed amplitude free (or using a CMB-only prior) and showing that the other models cannot be brought into agreement without violating CMB limits; the paper only performs this approximately for the stacked-emission test in Section 4.2.","section":"§3 (seed normalization) and §4.1/Figs. 4–5"}],"minor_comments":[{"comment":"The reference list contains duplicates that should be merged: [32] and [68] are the same Vernstrom et al. 2017 paper, [36] and [66] are the same Locatelli et al. 2021 paper, and [31] and [41] are the same Vernstrom et al. 2023 paper.","section":"Reference list"},{"comment":"The citation for the splashback radius is missing; the text contains a placeholder '[? ]' that should be replaced with the proper reference (likely Diemer et al. 2017, which appears in the bibliography).","section":"§4.3"},{"comment":"The criterion for flagging high-density points along the LOS is not specified; please state the quantitative threshold or give an explicit pointer to the exact procedure in Ref. [18] so that the analysis is reproducible.","section":"§4.1"},{"comment":"The text states an upper limit of 1.5 nG from Ref. [50] for the RM pair-difference experiment, while Table 1 lists ≤ 9 nG for the same entry; please reconcile this apparent inconsistency (e.g., proper versus comoving values or different filtering choices).","section":"§2.2 and Table 1"},{"comment":"The phrase 'Independently of the scenario and the shape of the astrophysical component RM' overstates the support from the analysis, because only two shapes (γ = 1 and γ = 5) are tested; the range 10–60 nG should be qualified as covering the two adopted shapes.","section":"Abstract and §6"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is essentially an extension of the authors' previous work [18] with an updated simulation and an additional probe, and the central quantitative claim is heavily tied to the authors' own earlier calibration of the nB = -1 seed amplitude. The editor may wish to consider whether the incremental advance over Ref. [18] is sufficient for the journal's scope, and whether the circularity in the model-selection step should be addressed more prominently in the revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take on Carretti & Vazza. The genuinely new thing is the joint comparison: RRM-z, synchrotron stacking, and group RM profiles, all run through the same suite of MHD simulations including the new B4 astrophysical model. That three-probe consistency check is real progress, and the paper is honest about the third probe being inconclusive. The qualitative conclusion—an astrophysical-only scenario cannot explain the high-redshift RRM—looks robust across the tested shapes, and the 10–60 nG range quoted for the field strength is a fair summary of the systematic spread. The nB=−1 primordial model is also individually interesting because its best-fit amplitude sits about five times below the CMB bound, which is the kind of claim radio data can actually make.\n\nThe soft spot is the quantitative headline. The 43±7 nG and the 'rapidly decreasing astrophysical component' come from fitting Eq. (9) with the astrophysical term fixed to A_rrm/(1+z)^5. Only two shapes are tested. With γ=1 you get B_f,0≈11–14 nG; with γ=5 you get 43–60 nG. The paper chooses γ=5 partly because it reproduces the observed 21±4% astrophysical fraction, but that integrated fraction does not pin down the redshift dependence. The two terms are added in quadrature and are degenerate in γ, so the RRM(z) data alone cannot select the shape. That means the 'best-matching combined model' is conditional on an assumption the data do not constrain. The paper actually says this—Section 6 calls for independent knowledge of the astrophysical component shape—so the limitation is stated, not hidden. But the abstract and conclusions still lead with the γ=5 number, and a busy reader could easily take 43±7 nG as the measured field.\n\nThere's also a milder circularity: the nB=−1 seed amplitude is calibrated on LOFAR RRMs in earlier work, and the same RRM-z relation is then used to select that model. That is a real concern, though the synchrotron stacking does provide some independent leverage.\n\nNone of this kills the paper. The qualitative message is probably right, and the joint framework is worth having. But the headline number needs to be framed as model-dependent, and the γ-shape issue deserves a serious treatment before this becomes the canonical value. I'd send it to peer review. It deserves referee time, and the discussion of the degeneracy would be constructive. I'd probably recommend major revision, not rejection, with the ask to reframe the conclusions and ideally add a systematic scan over γ.","headline":"Useful update with a genuinely new joint test, but the 43±7 nG headline is soft because it rests on an assumed γ=5 astrophysical shape that the LOFAR data do not constrain.","tokens_in":29125,"tokens_out":5208,"would_cite":true,"duration_ms":36490,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Low-frequency radio observations of cosmic filaments favour a primordial origin for cosmic magnetism, and place the average filament field at 10–60 nG at z=0.","keywords":["magnetic fields","intergalactic medium","large-scale structure of the Universe","methods: statistical","cosmic filaments","primordial magnetic fields","Faraday rotation","synchrotron cosmic web"],"falsifier":"Measure the redshift dependence of the astrophysical RM contribution directly at gigahertz frequencies, where this component dominates over the cosmic web. If an independent survey finds the astrophysical rms to be roughly proportional to $(1+z)^{-1}$ instead of $(1+z)^{-5}$, the best-fit filament field at $z=0$ shifts from $43\\pm7$ nG to roughly $11$–$14$ nG, and the RRM-vs-redshift argument for a dominant primordial seed loses its quantitative support.","tokens_in":27886,"feed_emoji":"🌌","tokens_out":17231,"duration_ms":126731,"temperature":0.7,"pith_summary":"The paper updates the comparison between radio observations of the cosmic web and recent cosmological MHD simulations to test which class of magnetogenesis scenarios is favoured: primordial fields seeded in the early Universe, astrophysical fields injected later by galaxy feedback, or a combination. It reports that two probes — the low-frequency rotation-measure scatter of filaments as a function of redshift and the stacked synchrotron emission of filaments — favour a dominant primordial magnetic field and disfavour a solely astrophysical origin, while a third probe, the rotation-measure radial profile around galaxy groups, does not yet give an unambiguous answer. If this is right, the magnetic fields filling the cosmic web preserve information from the early Universe, and low-frequency radio surveys can discriminate among primordial models more sensitively than CMB observations for inflationary-like spectra. The paper also gives an average filament field at $z=0$ of 10–60 nG independently of the model, and $43\\pm7$ nG for the combined model that best matches the simulations.","feed_headline":"Radio data favour a primordial origin for cosmic magnetism","feed_subtitle":"Rotation measures and stacked radio emission put filament fields at 10–60 nG, 43 ± 7 nG in the best-fit model.","key_machinery":"The machinery is the residual rotation measure (RRM) that remains after subtracting the Milky Way contribution, treated as a redshift-dependent statistic: the observed rms is fit to $\\langle RRM^2\\rangle^{1/2}=A_{rrm}/(1+z)^\\gamma+\\langle RRM_f^2\\rangle^{1/2}$, where the astrophysical term is a power law and the filament term is computed by integrating $0.812\\, n_e B_\\parallel/(1+z)^2$ along simulated lines of sight. The filament field itself is parametrized as $B_f=B_{f,0}(1+z)^\\alpha$, with electron densities drawn from the cosmological MHD simulations. The second, independent probe is the stacked synchrotron surface brightness of massive simulated filaments at about 118 MHz, computed with the shock-acceleration formula (Equation (1)), which responds approximately as $B^2$ and therefore breaks degeneracies that the RM alone leaves open. The third probe, the median-absolute-deviation radial profile of RRM around group-mass halos, is simulated and compared but is not yet conclusive.","core_discovery":"The paper's central claim is that the low-frequency radio data already select a magnetogenesis scenario: a dominant primordial field, specifically a stochastic seed with spectral slope $n_B=-1$ and $\\langle B\\rangle_{\\rm 1Mpc}=0.37$ nG combined with the most realistic astrophysical feedback model, while a purely astrophysical origin is disfavoured. In this combined model the average proper magnetic field in cosmic filaments is $B_{f,0}=43\\pm7$ nG at $z=0$ with a redshift slope $\\alpha=0.8\\pm0.5$, and the astrophysical term in the rotation-measure budget declines quickly with redshift. The paper also claims a scenario-independent range of $10$–$60$ nG for the filament field at $z=0$, and argues that the same radio observables can discriminate among primordial models more sensitively than current CMB analyses for inflationary-like spectra. The third probe, the rotation-measure radial profile around group-mass halos, is still inconclusive because the simulations fall short of the observed amplitude by roughly a factor of two.","pith_inferences":["Beyond the paper: if the steep $\\gamma=5$ astrophysical shape is right, astrophysical magnetization of the intergalactic medium is essentially a low-redshift phenomenon, which would imply that the magnetization of cosmic voids probed by gamma-ray pair-echo limits is also almost entirely primordial; the paper does not pursue this connection.","Beyond the paper: the paper's exclusion of the causal $n_B=2$ model is weakened by finite numerical resolution at about $41.5$ kpc cells; a higher-resolution rerun could determine whether artificially damped small-scale field energy accounts for the required factor of $\\sim100$–$200$ in magnetic energy, which would be a cleaner test than rescaling seed amplitudes.","Beyond the paper: the persistent factor-of-two shortfall near galaxy groups suggests the simulated feedback bubbles are too sparse or too weak, so observing the same group RM profiles with a full treatment of depolarization at low frequency could independently calibrate the astrophysical contribution and might shift the primordial/astrophysical balance at low redshift.","Beyond the paper: a direct cross-check would be to apply the same joint test to fast-radio-burst rotation and dispersion measures, whose simultaneous measurement maps $n_e B_\\parallel$ along individual lines of sight through filaments rather than statistically."],"forward_implications":["If the preferred model is correct, the bulk of the magnetic field energy in cosmic filaments at $z \\lesssim 3$ is a relic of an early-universe seed, and galaxy feedback contributes at most about a quarter of the observed low-frequency rotation-measure scatter.","Low-frequency radio surveys would then already constrain inflationary-like primordial field spectra roughly five times more tightly, in amplitude terms, than present CMB analyses.","The causal ($n_B=2$) phase-transition scenario would be excluded by the joint RM-plus-synchrotron test unless its seed amplitude exceeds CMB limits by about an order of magnitude, an exclusion that sharper CMB constraints can confirm.","The best-matching combined model predicts a rapidly decreasing astrophysical RM component, so high-redshift observations ($z>1$) should show almost purely primordial filament signal; this is a direct prediction of the model.","Pinpointing the redshift shape of the astrophysical RM component with gigahertz-frequency surveys is the single most informative next observation, because the current factor-of-four spread between the two shapes dominates the uncertainty in the filament field."],"supporting_citations":[{"why":"Provides the low-frequency rotation-measure rms versus redshift data and the measured 21 ± 4% astrophysical fraction that the fit in Equation (9) targets.","marker":"[18]"},{"why":"Supplies the updated astrophysical seeding model B4 and the calibrated cosmological MHD simulations used to build synthetic rotation measures and synchrotron stacks.","marker":"[13]"},{"why":"Gives the stacked synchrotron detection from hundreds of thousands of cluster-pair filaments that the simulated stacking must match.","marker":"[16]"},{"why":"Defines the stochastic primordial seed-field models with power-law spectra and CMB-compatible normalizations, including the preferred nB = -1 case.","marker":"[7]"},{"why":"Provides the radial rotation-measure profile around group-mass halos that is replicated for the third, inconclusive probe.","marker":"[63]"},{"why":"Establishes the earlier low-frequency radio constraints on filament magnetic field strength and evolution that motivate the power-law B(z) model.","marker":"[17]"},{"why":"Sets the CMB upper limits on primordial magnetic fields that fix the seed amplitudes of most simulated models and define the sensitivity comparison.","marker":"[2]"},{"why":"Supports the claim that radio RM probes of the cosmic web now constrain inflationary-like primordial fields more tightly than current CMB analysis.","marker":"[21]"}],"fun_headline_variants":["Radio probes point to primordial cosmic magnetism","Filament radio data favor primordial magnetic seeds","Rotation measures support primordial cosmic magnetic fields","Radio observations weigh in on cosmic magnetogenesis","Primordial fields win in cosmic filament radio data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is an unmeasured assumption about the redshift evolution of the astrophysical part of the rotation-measure signal: the paper prefers the steep form $A_{rrm}/(1+z)^5$ because it matches the observed $21\\%$ astrophysical fraction, but if the true astrophysical component evolves more mildly the inferred filament field drops from roughly $43$ nG to about $11$–$14$ nG and the preferred combined model is no longer forced by the data.","fun_headline_variants_meta":{"raw":{"variants":["Radio probes point to primordial cosmic magnetism","Filament radio data favor primordial magnetic seeds","Rotation measures support primordial cosmic magnetic fields","Radio observations weigh in on cosmic magnetogenesis","Primordial fields win in cosmic filament radio data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000424,"raw_usage":{"total_tokens":2184,"prompt_tokens":964,"completion_tokens":1220,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":580,"completion_tokens_details":{"reasoning_tokens":1155}},"tokens_in":580,"tokens_out":1220,"duration_ms":7720,"temperature":1.0,"reasoning_tokens":1155,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:28:47.681468+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the redshift dependence of the astrophysical RM contribution directly at gigahertz frequencies, where this component dominates over the cosmic web. If an independent survey finds the astrophysical rms to be roughly proportional to $(1+z)^{-1}$ instead of $(1+z)^{-5}$, the best-fit filament field at $z=0$ shifts from $43\\pm7$ nG to roughly $11$–$14$ nG, and the RRM-vs-redshift argument for a dominant primordial seed loses its quantitative support.","supporting_citations":[],"review_version":1}