{"id":"b36d4a57-0375-4be4-b709-85e5097af686","arxiv_id":"2506.18968","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Faraday complexity in 191 background radio sources is dominated by sub-2.5-arcsecond magneto-ionic structures in the Milky Way, likely linked to anisotropic turbulent magnetic fields.","lead":"Using radio observations of 191 distant galaxies behind the Milky Way's disk, the authors find that the twisting of radio waves varies sharply across each galaxy's image, more than expected from the galaxies themselves. This points to tiny, unseen magnetic structures in our Galaxy's gas, opening a new way to map small-scale Galactic magnetic fields.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No direct RM SF at sub-arcminute separations; the two-order discrepancy that rules out conventional turbulence is an extrapolated power-law result.","rationale":"I agree with the reader that the RM SF extrapolation is the load-bearing assumption. The size-independence of FD spread (Figure 5) is suggestive of small scales, but it does not by itself quantify the amplitude anomaly; the 'two orders of magnitude' statement is what rules out the conventional turbulent cascade. That statement rests on a power law measured between source pairs at separations much larger than 2.5 arcsec and extrapolated without a demonstrated break. The reader's verdict of CONDITIONAL is appropriate; no change needed. The unmatched comparison samples are a secondary issue; the RM SF extrapolation is the more direct support for the scale claim.","tokens_in":51128,"tokens_out":21626,"duration_ms":247264,"concrete_test":"Recompute Figure 6 using only RM SF bins with separation below 0.5 deg (or the smallest decade with enough pairs), refit the power-law slope, and extrapolate to 2.5 arcsec. If the refitted slope is flatter than +1.07 by more than its uncertainty, or if the extrapolated RM SF at 2.5 arcsec moves to within one order of magnitude of 2 x (FD Spread)^2, the two-order discrepancy is not robust. Additionally, fit a broken power law with a free break below 0.01 deg and compare Bayesian evidence against the no-break model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative inference (Section 3.1.5, Figures 6-7) is that 2 x (FD Spread)^2 lies up to two orders of magnitude above the RM structure function at angular scales of 2.5-300 arcsec, ruling out the conventional supernova-driven isotropic turbulent field. This comparison depends entirely on an extrapolation. The RM SF is built from integrated RMs of 191 sources observed with a ~50-arcsec VLA D-array beam; source pairs closer than roughly one beam cannot be separated, so there are no RM SF measurements at sub-arcminute separations. The plotted small-scale slope (+1.07 +/- 0.10) is fit over separations from about 0.01 deg (36 arcsec) up to the break at 3 deg, then extrapolated about three orders of magnitude downward to 2.5 arcsec. If the RM SF flattens or breaks below the smallest measured pair separation, the extrapolated amplitude at 2.5 arcsec could be much higher than the power-law value, and the claimed two-order discrepancy would shrink or disappear. Because that discrepancy is the main quantitative basis for excluding the conventional isotropic turbulent magnetic field (Section 4.2.3), the conclusion that the structures are <2.5 arcsec is not uniquely established by the amplitude comparison. The reader's condition is appropriate: either direct RM SF measurements at sub-arcminute separations, or a physical argument for why no break occurs, is required.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper re-uses VLA L-band spectro-polarimetric data from Ma et al. (2020) for 191 polarized extragalactic sources at Galactic latitude |b| <= 5 deg and longitude 20-52 deg, applies Stokes QU-fitting with eight astrophysical models, and defines a new 'FD spread' parameter (Eq. 12) to quantify intra-source Faraday-depth variance. The authors report an enhancement of FD spread near the Galactic mid-plane with a fitted exponential scale height of about 5 deg, hints of longitudinal modulation around the Scutum arm tangent, no dependence of FD spread on source angular size over 2.5 arcsec to 300 arcsec, and no correlation with spectral index or H-alpha intensity. By comparing 2 x (FD spread)^2 with an RM structure function constructed from the same sources (Figs. 6-7), they conclude that the supernova-driven isotropic turbulent magnetic field cannot explain the amplitude, and that the Faraday complexity is dominated by <2.5 arcsec-scale magneto-ionic structures in the Milky Way, possibly the anisotropic turbulent magnetic field or stellar-wind-driven turbulence.","tokens_in":51346,"tokens_out":8702,"duration_ms":87234,"significance":"If the central conclusion is correct, the paper opens a genuinely new observational window: broadband polarimetry of background EGSs could map sub-arcsecond to sub-parsec scale magneto-ionic structure in the Milky Way using unresolved sources. The paper is methodologically transparent and unusually thorough in its robustness checks: it publishes machine-readable QU-fitting results and spectra, tests foreground subtraction (App. C), Delta-BIC filtering (App. D), and exclusion of sinc-component models (App. E), and it provides three concrete, falsifiable predictions in Section 4.2.4. These strengths make the study a useful contribution even if the <2.5 arcsec scale interpretation is later refined. However, the strongest quantitative inference depends on an extrapolation of the RM structure function to angular scales where no direct measurements exist, and the extragalactic comparison is based on heterogeneous samples, so the paper's headline claim needs additional support or careful qualification.","major_comments":[{"comment":"The quantitative case against the conventional isotropic turbulent field is carried by comparing 2 x (FD spread)^2 with the RM structure function extrapolated down to 2.5 arcsec. The RM SF is computed from source pairs whose smallest angular separation is limited by the roughly 50 arcsec D-array beam, and the text does not state the minimum pair separation or the number of pairs in the innermost bins. The small-scale fitted slope (+1.07 +/- 0.10) is therefore extrapolated about a decade or more in angular scale down to 2.5 arcsec. If the RM SF flattens or breaks below the smallest measured separations, as expected near an inner scale of magneto-ionic turbulence, the predicted amplitude of 2 x sigma_FD^2 at 2.5 arcsec could be much higher and the claimed two-order discrepancy could largely disappear. Because this discrepancy is the main basis for excluding the supernova-driven isotropic turbulent field in Section 4.2.3, the conclusion is not uniquely established. Please report the minimum separation and fit range explicitly, test the sensitivity to removing the smallest bins, and either provide direct RM SF measurements at sub-arcminute separations or a physical argument bounding the inner scale and break behavior.","section":"Section 3.1.5, Figs. 6-7, Eq. (17)"},{"comment":"The claim in the abstract that the FD spread amplitude is higher than expected from extragalactic structures is only weakly supported by the statistical tests as presented. The KS tests comparing the full sample with Anderson et al. (2015) and O'Sullivan et al. (2017) give p = 0.13 and p = 0.33, respectively, meaning the full-sample distributions are not significantly different; only the |b| < 3 deg subset gives p = 0.058 and p = 8.6 x 10^-3 against those two samples. The comparison samples also differ in frequency band, angular resolution, sky region, and in the use of the second moment M2 versus the newly defined FD spread, so the quantitative distribution comparison is not apples-to-apples. Please either repeat the comparison on matched subsamples with comparable frequency coverage and resolution, or explicitly qualify the amplitude claim as applying only to the low-latitude subset.","section":"Section 4.1.2, Fig. 10"},{"comment":"The inference that the responsible magneto-ionic structures are smaller than 2.5 arcsec uses the total-intensity angular size of each EGS as the scale sampled by the polarization measurement. The VLASS and RACS images trace all radio emission, whereas the polarized emission may be concentrated in compact cores or hotspots embedded in larger lobes. If the polarized emission is compact, the FD spread would be nearly independent of total source size even if the foreground magneto-ionic structures have scales of tens of arcseconds or more. The flat relation in Figure 5 therefore does not by itself place the structure scale below 2.5 arcsec. Please quantify the size of the polarized emitting regions, for example with high-resolution polarimetric data for a subsample, or by splitting the sample according to morphology and fractional polarization, before asserting the sub-2.5-arcsec scale in the conclusions.","section":"Section 2.4, Section 3.1.5, Fig. 5"},{"comment":"FD spread is a model-dependent composite: depending on the best-fit Stokes QU model, it reduces to sigma_phi, Delta_phi, or a mixture of component separation and dispersion terms. The paper provides signal-to-noise cuts and robustness checks against Delta BIC and sinc-component exclusion, but no end-to-end injection-recovery simulations showing that the FD spread estimator is unbiased across the full range of input models, signal-to-noise, and frequency coverage. Such simulations would also test whether measurement noise can produce spurious multi-component fits and hence a positive FD-spread floor that varies with sky position. Given that FD spread is the central new observable, adding at least a targeted noise-injection study would materially strengthen the quantitative latitude and longitude results.","section":"Eq. (12), Section 2.3"}],"minor_comments":[{"comment":"The axis labels in the provided text extract appear garbled (e.g. '2 x (FD )p(ea )2'); please ensure the typeset version has correct mathematical notation.","section":"Figs. 6 and 7"},{"comment":"Please state explicitly the smallest angular separation used in the RM SF, the number of pairs per bin, and the range over which the broken power-law slopes are fitted, so that the reader can assess the extrapolation.","section":"Section 3.1.5"},{"comment":"The phrases 'ruled out' and 'impossible' are stronger than the evidence warrants given the extrapolation discussed in the first major comment; consider phrasing such conclusions as inconsistent with an unbroken power-law extrapolation of the RM SF.","section":"Section 4.2.3"},{"comment":"The factor 1/N applies only to the first term in the definition of the FD spread; making this explicit with clearly separated summands would avoid ambiguity.","section":"Eq. (12)"}],"recommendation":"major_revision","confidential_remarks":"This is a well-executed and potentially influential paper that is clearly within the scope of MNRAS. My main reservation is that the abstract's strongest claim (that the Faraday complexity is dominated by <2.5 arcsec-scale Milky Way structures) is not uniquely established by the RM SF extrapolation or by total-intensity source sizes. I would be comfortable with publication after the authors either provide direct high-resolution RM SF or polarization-size constraints, or substantially soften and carefully qualify the headline claim. The heterogeneous extragalactic comparison and the absence of injection-recovery tests for the new FD spread parameter should also be acknowledged in the text."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper is worth your time. It introduces the FD spread parameter and uses it to claim that Faraday complexity seen toward 191 EGSs in the northern Galactic plane comes from magneto-ionic structures smaller than 2.5 arcsec in the Milky Way, not from the sources or instrumental effects. The new parameter is transparently defined and the analysis is refreshingly careful: foreground subtraction, BIC filtering, exclusion of sinc models, and appendix checks all support the main statistical claims. The latitude enhancement (stronger within |b|<3°) is solid, and the constancy of FD spread with source size is an interesting new observational fact.\n\nThe soft spots are real but less severe than a first pass suggests. The RM structure function comparison (Figure 6) extrapolates a power law fit from separations above ~36 arcsec down to 2.5 arcsec. That is about an order of magnitude, not the 'three orders' claimed in the stress-test note, and the mismatch between 2×(FD spread)² and the extrapolated RM SF is so large that even a flattened RM SF at small scales would not erase it entirely. Still, the conclusion that conventional supernova-driven turbulence is excluded depends on that extrapolation, and the paper would be stronger with a physical argument for no break or with higher-resolution RM SF data. The extragalactic comparison is also heterogeneous — different frequencies, resolutions, and sample selections — and the KS tests are not significant for the full sample against two of the three comparison works. The paper acknowledges this and does not oversell the extragalactic exclusion, but the abstract's central claim leans on it.\n\nWho is this for? Anyone working on Faraday rotation, Galactic magnetic fields, or polarized source populations. The paper is honest, the new parameter is useful, and the robustness checks set a good standard. It deserves a serious referee and likely publication after the authors temper the small-scale conclusion or add supporting evidence.\n\nRecommendation: engage, and send the authors back to the RM SF extrapolation and the comparison-sample matching.","headline":"A careful, honest paper that likely opens a new observational window into small-scale Galactic magnetic structure, but the central amplitude comparison rests on an extrapolated RM structure function and a heterogeneous extragalactic benchmark.","tokens_in":51974,"tokens_out":2888,"would_cite":true,"duration_ms":32006,"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":"The Faraday complexity seen toward 191 extragalactic radio sources near the Galactic plane comes primarily from Milky Way magneto-ionic structures smaller than 2.5 arcseconds, not from the sources or the telescope.","keywords":["Faraday complexity","Faraday depth","Galactic magnetic fields","interstellar medium","radio polarimetry","rotation measure","magneto-ionic structures","turbulence"],"falsifier":"A direct measurement of the rotation-measure structure function below an arcminute, for example from high-resolution Faraday depth maps of a compact source field at about one arcsecond resolution, would show whether the RM variance actually reaches the level implied by the extrapolation; if the power law breaks and rises before 2.5 arcseconds, the two-order-of-magnitude gap disappears. Alternatively, a survey toward the Galactic anticentre testing the paper's prediction that no enhanced FD spread should appear there would distinguish the anisotropic-field explanation from a stellar-feedback origin.","tokens_in":50900,"feed_emoji":"🧲","tokens_out":7167,"duration_ms":68455,"temperature":0.7,"pith_summary":"This paper tries to establish where the 'Faraday complexity' of background radio sources near the Milky Way's plane comes from: the sources themselves, the instrument, or gas and magnetic fields inside our Galaxy. Using 191 polarised extragalactic sources observed at 1–2 GHz, the authors define a new quantity, the FD spread, that captures the spatial scatter of Faraday depth within each source. They find the spread is largest near the Galactic mid-plane, roughly constant for source sizes from 2.5 to 300 arcseconds, and larger than extrapolations of the rotation-measure structure function predict. Together the evidence points to magneto-ionic structures smaller than 2.5 arcseconds in the Milky Way—most plausibly an anisotropic turbulent magnetic field shaped by spiral-arm shocks and shear, or small-scale turbulence from stellar winds—as the primary cause.","feed_headline":"Milky Way's sub-2.5-arcsec structures dominate Faraday complexity","feed_subtitle":"A new measure of Faraday-depth jitter shows tiny Galactic magnetic structures, not the background radio sources, drive the signal.","key_machinery":"The load-bearing object is the FD spread, a single number per source constructed from the best-fit Stokes QU-fitting model: $$\\text{FD spread}=\\sqrt{\\frac{1}{N}\\sum_i(\\phi_i-\\bar{\\phi})^2+\\sum_i\\sigma_{\\phi,i}^2+\\sum_i\\$\\Delta$\\$phi_i^{2}$}.$$ It converts the multi-component Faraday depth information of each source into a uniform measurement of spatial FD fluctuation, deliberately not weighting by polarised intensity as the RM-spectrum second moment does. The argument then leans on comparing twice the square of the FD spread against the RM structure function, exploiting the identity that a Gaussian FD distribution with variance $\\sigma^2$ produces an RM structure-function amplitude of $2\\sigma^2$. The paper's conclusion follows from the mismatch between the two quantities at angular scales of 2.5–300 arcseconds, combined with the absence of any dependence of FD spread on source size.","core_discovery":"On the paper's own terms, the central discovery is that the Faraday complexity seen in the sample is dominated by sub-2.5-arcsecond magneto-ionic structures in the Milky Way foreground, not by the extragalactic sources or by instrumental artifacts. The new FD spread statistic has a mean of about 42 rad m$^{-2}$ (65 rad m$^{-2}$ among Faraday-complex sources), enhances toward $|b|<3^\\circ$, shows a hint of excess around the Scutum arm tangent near $\\ell=28$–$32^\\circ$, and is independent of source angular size between 2.5 and 300 arcseconds. The decisive comparison is that twice the square of the FD spread lies up to two orders of magnitude above the rotation-measure structure function extrapolated to the same angular scales, so the conventional supernova-driven isotropic turbulent field cannot account for the signal.","pith_inferences":["If the sub-2.5-arcsecond structures are real and widespread, the effective small-scale Faraday sky is much rougher than turbulence-cascade models predict, implying an additional energy-injection or anisotropy mechanism at physical scales below about 0.1 parsec at typical Galactic-plane distances.","The FD spread statistic could be applied to polarised sources behind other spiral galaxies or gas-rich systems, where the foreground structure function is measured at parsec scales, to test whether such ultra-compact magneto-ionic structure is generic rather than unique to the Milky Way.","A direct statistical test would compare the FD spread distribution of sources behind the Galactic plane with a matched sample at high latitude at the same angular resolution, isolating the Galactic contribution without relying on the structure-function extrapolation."],"forward_implications":["Galactic magnetism studies that use compact extragalactic sources as Faraday probes must treat each source's FD as potentially contaminated by sub-2.5-arcsecond Milky Way structures; selecting spatially extended sources and discarding Faraday-simple compact sources may be the safer strategy in the surveyed region.","High-resolution spectro-polarimetric follow-up at about one arcsecond should resolve the responsible structures into patches of different Faraday depth, and if the anisotropic turbulent field is the cause, the coherence length should be shorter along Galactic longitude than latitude.","The Scutum-arm-tangent excess, where only 15 percent of sources are Faraday simple compared with 41 percent elsewhere, predicts that Faraday complexity can serve as a tracer of spiral-arm compressions.","Omitting foreground diffuse-emission subtraction inflates the measured FD spread by about 50 percent, so future RM surveys in this region need that subtraction before quantifying Faraday complexity.","The two proposed mechanisms make opposite predictions toward the Galactic anticentre: no enhanced FD spread if anisotropic fields aligned with ring tangents dominate, versus possibly still elevated spread from stellar-wind bubbles, so observations there can separate the two."],"supporting_citations":[{"why":"Supplies the VLA L-band spectro-polarimetric image cubes and the source sample whose Stokes IQU spectra are re-extracted here.","marker":"Ma et al. 2020"},{"why":"Provides the high-latitude Faraday-complexity distribution used as an upper-limit comparison for the extragalactic contribution.","marker":"Anderson et al. 2015"},{"why":"High-latitude Stokes QU-fitting sample used to estimate the expected FD spread from source-intrinsic Faraday complexity.","marker":"O'Sullivan et al. 2017"},{"why":"Galactic-centre Faraday-complexity sample showing where the complexity is enhanced by known turbulent structures.","marker":"Livingston et al. 2021"},{"why":"Southern Galactic plane study that found enhanced Faraday complexity near spiral-arm tangents, supporting the longitude modulation seen here.","marker":"Ranchod et al. 2024"},{"why":"Supplies coherence lengths and RM structure-function behaviour for supernova-driven isotropic turbulence against which the FD spread is compared.","marker":"Haverkorn et al. 2008"},{"why":"Establishes the 2-sigma-squared relation between Faraday-depth variance and RM structure-function amplitude used in the comparison figures.","marker":"Seta et al. 2023"},{"why":"Derives the sinc modulation produced by a plane-of-sky Faraday-depth gradient, the basis for treating Delta-phi as a foreground probe.","marker":"Sokoloff et al. 1998"},{"why":"Provides the RM structure function in the inertial range that sets the expected slope and amplitude at arcminute-to-degree scales.","marker":"Stil et al. 2011"},{"why":"RM-Synthesis method that underlies the Faraday spectra and sets the FD resolution and maximum-scale limits for these data.","marker":"Brentjens & de Bruyn 2005"}],"fun_headline_variants":["Sub-2.5-arcsec Milky Way structures dominate Faraday complexity","Tiny Galactic magnetic knots, not quasars, drive Faraday complexity","New FD-spread probe ties Faraday complexity to Milky Way small scales","Sub-arcsecond Galactic magnetic clumps control Faraday complexity","Milky Way's second-scale magnetic wiggles rule Faraday sky"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument that the FD spread amplitude exceeds the expected turbulence signal depends on the rotation-measure structure function power law measured between source pairs at larger separations continuing down to 2.5 arcseconds without a break or flattening.","fun_headline_variants_meta":{"raw":{"variants":["Sub-2.5-arcsec Milky Way structures dominate Faraday complexity","Tiny Galactic magnetic knots, not quasars, drive Faraday complexity","New FD-spread probe ties Faraday complexity to Milky Way small scales","Sub-arcsecond Galactic magnetic clumps control Faraday complexity","Milky Way's second-scale magnetic wiggles rule Faraday sky"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000617,"raw_usage":{"total_tokens":2913,"prompt_tokens":1046,"completion_tokens":1867,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":662,"completion_tokens_details":{"reasoning_tokens":1777}},"tokens_in":662,"tokens_out":1867,"duration_ms":15177,"temperature":1.0,"reasoning_tokens":1777,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:40:50.310356+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct measurement of the rotation-measure structure function below an arcminute, for example from high-resolution Faraday depth maps of a compact source field at about one arcsecond resolution, would show whether the RM variance actually reaches the level implied by the extrapolation; if the power law breaks and rises before 2.5 arcseconds, the two-order-of-magnitude gap disappears. Alternatively, a survey toward the Galactic anticentre testing the paper's prediction that no enhanced FD spread should appear there would distinguish the anisotropic-field explanation from a stellar-feedback origin.","supporting_citations":[],"review_version":2}