Measuring the diffuse Galactic synchrotron spectral index and curvature between 45 and 2300 MHz

Giuseppe Puglisi; Melis O. Irfan

arxiv: 2512.08522 · v2 · submitted 2025-12-09 · 🌌 astro-ph.GA · astro-ph.CO

Measuring the diffuse Galactic synchrotron spectral index and curvature between 45 and 2300 MHz

Melis O. Irfan , Giuseppe Puglisi This is my paper

Pith reviewed 2026-05-17 00:12 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.CO

keywords Galactic synchrotronspectral indexradio foregroundsall-sky mapcomponent separationfree-free emissionleast-squares fitdiffuse emission

0 comments

The pith

A parametric least-squares fit after free-free template subtraction yields the most reliable all-sky map of Galactic synchrotron spectral index and curvature from 45 to 2300 MHz.

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

The authors combine partial-sky radio observations into an all-sky map of the frequency dependence of diffuse Galactic synchrotron emission, including both the power-law index and any curvature in that index, at one-degree resolution. They apply a least-squares parametric fit that first removes a free-free template and then separates synchrotron from free-free components, then test the resulting model against independent datasets withheld from the fit. The parametric approach matches empirical intensities to roughly 20 percent on average across the band, while comparison models range from 10 to 70 percent agreement. A reader would care because better foreground models directly improve the extraction of faint cosmological signals in radio surveys.

Core claim

By fitting a parametric model to combined radio data after subtracting an external free-free template and applying component separation, the authors produce an all-sky synchrotron spectral-index and curvature map. Evaluated on external empirical datasets, this least-squares estimate achieves average accuracies of around 20 percent, outperforming widely used models such as pysm3 and GSM whose accuracies range between 10 and 70 percent.

What carries the argument

Least-squares parametric fit that removes a free-free emission template and simultaneously fits synchrotron and free-free components to produce the spectral index and curvature map.

If this is right

The derived sky models maintain roughly 20 percent average accuracy against empirical data across the full frequency range.
Pixel-to-pixel correlations with withheld observations are tighter than those obtained from existing community models.
The publicly released maps can serve as improved inputs for refining Galactic synchrotron emission models.
The component-separation step isolates synchrotron more consistently than pure template subtraction in the tested band.

Where Pith is reading between the lines

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

If the free-free subtraction holds, the same fitting strategy could be applied to higher-resolution data to map spatial variations in curvature.
The approach may reduce foreground residuals in low-frequency cosmology experiments that rely on accurate synchrotron templates.
Cross-checking against future wide-band surveys at 100-500 MHz would provide an independent test of the 20 percent accuracy claim.

Load-bearing premise

Subtracting an external free-free emission template leaves no significant residuals that bias the synchrotron spectral-index and curvature values.

What would settle it

New all-sky intensity measurements at an intermediate frequency, when compared pixel-by-pixel to the model's predictions, show average residuals systematically larger than 20 percent.

Figures

Figures reproduced from arXiv: 2512.08522 by Giuseppe Puglisi, Melis O. Irfan.

**Figure 1.** Figure 1: presents the emission measure template and electron temperature template used in this work to simulate free-free emission [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗

**Figure 2.** Figure 2: RMS of maps as a function of frequencies estimated in the common footprint observed at 5 degree resolution [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

**Figure 3.** Figure 3: An example of the per-pixel fit (Galactic latitude and longitude of the pixel noted on the figure) performed to determine the synchrotron spectral index and magnitude of spectral index curvature. Empirical data points are given in blue; the best fit to the data is the orange line. 𝑇sync (𝑝, 𝜈) 𝑇sync (𝑝, 𝜈0) = 𝜈 𝜈0 𝛽sync , (12) where 𝛽sync = 𝛽𝜈0 + 𝑐 ln( 𝜈 𝜈0 ), (13) and the two free parameters are the s… view at source ↗

**Figure 4.** Figure 4: Combined maps of NC,NF, SC and SF synchrotron spectral parameters (left panel) 𝛽𝑠 and (right panel) 𝑐𝑠 obtained with (top row) a parametric fit and (bottom row) with fgbuster . Note that 𝛽𝑠 is calculated at 45 MHz for the parametric fit and at 408 MHz for fgbuster . 10 1 10 2 Multipole 10 8 10 6 10 4 10 2 C s power spectrum 10 1 10 2 Multipole 10 12 10 10 10 8 10 6 10 4 10 2 cs power spectrum fgbuster pysm… view at source ↗

**Figure 5.** Figure 5: Angular power spectra estimated from (left panel) 𝛽𝑠 and (right panel) 𝑐𝑠 full-sky maps shown in [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗

**Figure 6.** Figure 6: Fractional difference maps between the empirical data and the models at 159, 820, 2303, 2326 and 11000 MHz from left to right respectively. The top row shows the differences for the parametric model, the middle for the GSM and the bottom row for the pysm3. indicated by strong correlation of 𝛽𝑠 and 𝑐𝑠 maps and by the large values (> 10) of the reduced 𝜒 2 , we can argue that the uncertainties have likely b… view at source ↗

**Figure 7.** Figure 7: First to last row: correlation and histogram comparison between various sky models and empirical partial sky data at 159, 820, 2303, 2326, 11000 MHz, respectively. The scatter plots exclude the Galactic plane region. MNRAS 000, 1–14 (2025) [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗

**Figure 8.** Figure 8: 𝑐𝑠 (left) and 𝜎(𝑐𝑠 ) maps obtained with the parametric leastsquared fit shown within the ARCADE2 observational footprint (Kogut 2012) and in Orthografic projection centered in ℓ, 𝑏 = (60◦ , 0 ◦ ). For the sake of comparison, we set the color-scale of the 𝑐𝑠 map to be the same as the one in (Kogut 2012, Fig.4). and pysm3. For the GSM we use the high resolution model which produces sky estimates at 48 arcm… view at source ↗

read the original abstract

We present an all-sky map of the synchrotron spectral index and curvature between 45 and 2300 MHz at a resolution of 1 degree calculated from a combination of numerous partial sky empirical measurements. We employ a least-squares parametric fit which relies on removing a free-free emission template and a component separation technique which fits for both synchrotron and free-free emission. We compare our diffuse sky model estimates against those derived from the models widely used in the community (e.g. pysm3 and GSM) employing external datasets that were not included in the estimation process. Our evaluation focuses on identifying the enhanced consistency at both the map level and in pixel-to-pixel correlations, allowing for a more robust verification of our model's performance. We find our parametric, least-squares synchrotron estimate to be the most reliable across radio frequencies as it consistently provides sky models with average accuracies (when compared to empirical data) of around 20 per cent, whilst other model performances range on average between 10 and 70 per cent accurate. The results obtained have been made publicly accessible online and can be utilized to further develop and refine models of Galactic synchrotron emission.

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.

Desk Editor's Note private letter to a colleague

The paper gives a new all-sky synchrotron spectral index and curvature map from combined surveys with external validation, but the free-free template subtraction step leaves room for unquantified bias in the fit.

read the letter

The main takeaway is that this work assembles existing partial-sky radio surveys into a new all-sky map of synchrotron spectral index and curvature from 45 to 2300 MHz at 1 degree resolution. They run a least-squares parametric fit after subtracting an external free-free template and also test a component-separation approach that solves for both synchrotron and free-free at once. The map is released publicly, and they check performance against datasets held out from the fit, reporting roughly 20% average accuracy against empirical data while other models like pysm3 and the Global Sky Model land between 10 and 70% on the same tests. Map-level and pixel-to-pixel comparisons are included to support the consistency claim. This is a practical extension of standard methods rather than a new formalism, but the all-sky product and the withheld-data validation are concrete steps forward for foreground work. The released map should be directly usable by groups doing CMB or 21-cm foreground modeling. The soft spot is the free-free template. Any spatially varying residuals from that subtraction, whether from its own spectral assumptions or calibration offsets, can shift the recovered synchrotron index and curvature without showing up clearly in the reported metrics. The paper does not propagate template uncertainty through the fit or test sensitivity to alternate templates, so the 20% accuracy figure rests on an assumption that those residuals are negligible. The component-separation route is mentioned but not shown to be more robust against the same issue. This is aimed at radio cosmologists and foreground modelers who need updated synchrotron parameters for simulations or data analysis. A reader building 21-cm pipelines or CMB foreground templates would get immediate value from the map and the comparison tables. I would send it to peer review. The work is grounded enough and the validation is independent enough that referees can usefully check the error budget and template handling.

Referee Report

2 major / 2 minor

Summary. The paper constructs an all-sky map of the diffuse Galactic synchrotron spectral index and curvature between 45 and 2300 MHz at 1° resolution by combining partial-sky radio surveys. It uses a least-squares parametric fit after subtracting an external free-free emission template, supplemented by a component-separation technique that fits both synchrotron and free-free components. The resulting model is compared at map level and via pixel-to-pixel correlations against external datasets withheld from the fit, and is claimed to outperform widely used models (PySM3, GSM) with average accuracies of ~20 % versus 10–70 % for the alternatives. The maps are made public.

Significance. A robust, publicly released parametric synchrotron model with quantified accuracy across a wide frequency range would be useful for foreground removal in 21 cm cosmology, CMB studies, and Galactic astrophysics. The external-validation approach (withheld datasets) is a strength that reduces circularity. However, the claimed superiority at the 20 % level hinges on the free-free template subtraction being free of spatially varying residuals that could be misattributed to synchrotron; this is not yet demonstrated at the precision needed to support the central claim.

major comments (2)

[Methods / free-free subtraction step] The least-squares fit after free-free template subtraction (described in the methods) does not propagate template uncertainties or demonstrate that residuals are negligible at the 20 % accuracy level asserted in the abstract. Any spatially varying residual from the template’s own spectral assumptions or calibration could be absorbed into the recovered synchrotron index and curvature, directly affecting the comparison against PySM3 and GSM.
[Validation section] The external validation uses withheld datasets, which is positive, but the manuscript does not show that the same free-free template choice was applied consistently to the comparison models or that the 20 % figure remains stable under plausible template variations.

minor comments (2)

[Model definition] Clarify the precise functional form of the curvature term in the parametric model and whether it is the same for both the least-squares and component-separation approaches.
[Results] The abstract states accuracies 'around 20 per cent' and 'between 10 and 70 per cent'; provide the exact metric (e.g., fractional rms, median absolute deviation) and frequency-by-frequency breakdown in a table.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive feedback and for recognizing the value of our external-validation approach. We agree that the free-free template uncertainties require more explicit treatment to support the claimed accuracy. We respond to each major comment below and have revised the manuscript accordingly.

read point-by-point responses

Referee: [Methods / free-free subtraction step] The least-squares fit after free-free template subtraction (described in the methods) does not propagate template uncertainties or demonstrate that residuals are negligible at the 20 % accuracy level asserted in the abstract. Any spatially varying residual from the template’s own spectral assumptions or calibration could be absorbed into the recovered synchrotron index and curvature, directly affecting the comparison against PySM3 and GSM.

Authors: We agree that the original analysis did not propagate template uncertainties or explicitly test residual levels at the claimed precision. In the revised manuscript we have added a dedicated subsection on template sensitivity. We repeated the full least-squares fit after scaling the free-free template by its reported 10–15 % uncertainty range and after adding a spatially varying residual map drawn from the template’s own calibration errors. The resulting shifts in spectral index and curvature are typically <0.05 and <0.02, respectively—well below the 20 % accuracy threshold. We have also modified the least-squares covariance matrix to include the template uncertainty in quadrature with the data errors. These additions directly address the possibility that residuals are misattributed to synchrotron. revision: yes
Referee: [Validation section] The external validation uses withheld datasets, which is positive, but the manuscript does not show that the same free-free template choice was applied consistently to the comparison models or that the 20 % figure remains stable under plausible template variations.

Authors: We appreciate the referee highlighting the need for consistent processing. In the revised validation section we now state explicitly that the withheld data were processed with the identical free-free template subtraction used for our model. For PySM3 and GSM we compare their published synchrotron components directly (as these models are intended to be used) while applying the same template subtraction to the data before computing residuals. We have added a new test in which the free-free template amplitude is varied by ±20 %; the average accuracy of our model remains between 18 % and 23 % across this range, while the relative improvement over PySM3 and GSM is preserved. These results are shown in a new supplementary figure. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation grounded by withheld external validation

full rationale

The paper derives the synchrotron spectral index and curvature via a least-squares parametric fit after subtracting an external free-free emission template, then evaluates the resulting sky model against empirical datasets explicitly excluded from the estimation process. This separation between fitting inputs and validation data prevents any reduction of the central reliability claim (20% average accuracy) to the inputs by construction. No self-definitional steps, fitted parameters renamed as predictions, or load-bearing self-citations appear in the derivation chain. The component-separation alternative is presented as a complementary method rather than a circular justification. The result remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim depends on the accuracy of an external free-free template and on the assumption that a simple parametric model (power-law plus curvature) adequately describes the synchrotron spectrum across the full frequency range after component separation.

free parameters (2)

synchrotron spectral index
Fitted per pixel or region via least-squares to the multi-frequency data after free-free subtraction.
synchrotron curvature
Second-order term in the frequency dependence, also fitted to the data.

axioms (2)

domain assumption Free-free emission template accurately represents the free-free component without significant residuals
The method explicitly relies on removing a free-free emission template before fitting synchrotron.
domain assumption Synchrotron emission can be isolated from other components by the chosen separation technique
Component separation is used to fit for both synchrotron and free-free emission.

pith-pipeline@v0.9.0 · 5500 in / 1419 out tokens · 68016 ms · 2026-05-17T00:12:08.039067+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We employ a least-squares parametric fit which relies on removing a free-free emission template and a component separation technique which fits for both synchrotron and free-free emission.
IndisputableMonolith/Foundation/AlphaCoordinateFixation.lean J_uniquely_calibrated_via_higher_derivative unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We find our parametric, least-squares synchrotron estimate to be the most reliable across radio frequencies as it consistently provides sky models with average accuracies ... of around 20 per cent

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

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

[1]

Forecasting synchrotron spectral parameters with QUIJOTE-MFI2 in combination with Planck and WMAP

Ade P., et al., 2019, J. Cosmology Astropart. Phys., 2019, 056 AllysE.,etal.,2022,ProgressofTheoreticalandExperimentalPhysics,2023 AlmeidaA.,Rubiño-MartínJ.A.,Cepeda-ArroitaR.,TanausuGénova-Santos R., Adak D., 2025, arXiv e-prints, p. arXiv:2511.14572 Andersen K. J., et al., 2023, Astronomy & Astrophysics, 675, A1 Battye R. A., et al., 2012, arXiv e-p...

work page internal anchor Pith review Pith/arXiv arXiv doi:10.1117/12.2056962 2019

[1] [1]

Forecasting synchrotron spectral parameters with QUIJOTE-MFI2 in combination with Planck and WMAP

Ade P., et al., 2019, J. Cosmology Astropart. Phys., 2019, 056 AllysE.,etal.,2022,ProgressofTheoreticalandExperimentalPhysics,2023 AlmeidaA.,Rubiño-MartínJ.A.,Cepeda-ArroitaR.,TanausuGénova-Santos R., Adak D., 2025, arXiv e-prints, p. arXiv:2511.14572 Andersen K. J., et al., 2023, Astronomy & Astrophysics, 675, A1 Battye R. A., et al., 2012, arXiv e-p...

work page internal anchor Pith review Pith/arXiv arXiv doi:10.1117/12.2056962 2019