{"id":"52801745-6ed8-47e9-9a7f-1de4705cbca7","arxiv_id":"2506.14310","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A 154 MHz measurement of the Galactic synchrotron angular power spectrum from six MWA pointing centers, with fitted spectral indices 0.9 to 1.7 that agree with prior work.","lead":"This paper measures the angular power spectrum of diffuse Galactic synchrotron emission at 154 MHz using MWA drift scan observations, fitting a power-law model to six sky positions. The measured spectral indices are consistent with earlier 150 MHz measurements, providing new foreground data for 21-cm cosmology experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Residual C_l may be dominated by source-subtraction systematics rather than DGSE; the paper lacks a point-source-only null simulation showing the shallow l<200 component is real.","rationale":"I read the paper as a careful, incremental measurement whose central claim is that, after subtracting sources above ~430 mJy, the residual angular power spectrum at 65<l<650 contains a genuine diffuse Galactic synchrotron component with fitted slope beta in the range 0.9-1.7. The reader's weakest-assumption identification matches my own: the residual power could be dominated by calibration errors, unmodeled compact sources, or ionospheric effects rather than DGSE. The paper has real strengths: the TGE has been validated in earlier work, source positions and fluxes are checked against GLEAM for one pointing, and the authors are transparent about the limitations of the unsuccessful pointings. However, the step from 'the residual C_l is shallower than l^2 below l~200' to 'this is DGSE' is not secured. The paper's own discussion of the failed pointings (Section 4) concedes that imaging artifacts around bright sources and residual gain errors can affect C_l, and those systematics are not demonstrably absent for the six successful pointings. The missing test is a null simulation containing only point sources and noise; without it, the shallow component could be a signature of imperfect subtraction rather than diffuse emission. This is an addressable concern rather than a fatal flaw, so I do not move the verdict: the paper should remain CONDITIONAL, with the null simulation as the key condition for acceptance.","tokens_in":15560,"tokens_out":5443,"duration_ms":61528,"concrete_test":"Construct an end-to-end null simulation with no diffuse emission: place only the point sources from the fitted PC source catalogues (with flux and position errors matching the GLEAM comparison), add thermal noise at the reported r.m.s., and push these visibilities through the same imaging, CLEAN, 3-sigma source-subtraction, and TGE pipeline used for the data. Compare the residual C_l against the observed UV-Sub C_l over 65<l<650. If the null simulation alone reproduces the observed residual within 1-sigma, the shallow component and fitted beta are not evidence for DGSE. If the observed residual is significantly above the null at l<200, the DGSE interpretation is supported. As a further check, inject a diffuse component with known beta (e.g., 1.5) into the simulation and verify that the full pipeline recovers it without bias.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central inference that the shallow residual component at l<200 (Section 3, after Fig. 1) is genuine DGSE, and that fitting eq. (4) over 65<l<650 recovers its slope beta, has not been separated from source-subtraction and calibration systematics. The paper validates source astrometry/fluxes against GLEAM for one pointing (Appendix 1) and validates TGE normalization with UAPS simulations, but it never runs a null test with a sky that contains only point sources below the subtraction threshold plus noise. This matters because the residual of an imperfectly subtracted bright source is not white: the Appendix itself reports a ~30 arcsec systematic DEC offset and flux deviations up to 25% for some sources, and the paper later attributes flat C_l for the unsuccessful PCs (Fig. 4) to 'imaging artifacts around bright sources, possibly arising from residual gain calibration errors' (Section 4). Such artifacts inject large-scale, shallow-spectrum power that can mimic the A(1000/l)^beta term in eq. (4), particularly since A, beta, and C are strongly degenerate (|r| ~ 0.9, Section 3). Without a demonstration that the observed residual C_l exceeds the point-source-only/systematics expectation, the beta range 0.9-1.7 is not uniquely attributable to DGSE.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes MWA Phase II drift-scan observations at 154.2 MHz to measure the angular power spectrum C_l of diffuse Galactic synchrotron emission (DGSE). After removing point sources above a 3-sigma threshold (~430 mJy), the authors apply the Tapered Gridded Estimator (TGE) to residual visibilities and fit the model C_l = A(1000/l)^beta + C over the multipole range 65 < l < 650. They report successful fits for 6 of 24 pointing centers, with beta in the range 0.9-1.7 and A in the range 155-568 mK^2. The power-law component is interpreted as DGSE and the constant as Poisson fluctuations from residual point sources. The paper validates the source subtraction for one pointing against GLEAM and uses MCMC to estimate parameter uncertainties.","tokens_in":15836,"tokens_out":6832,"duration_ms":67100,"significance":"If the measurement is robust, it provides a new estimate of the DGSE angular power spectrum at 154 MHz in a region relevant for EoR experiments, complementing earlier GMRT and WSRT measurements. The use of drift-scan data and the well-established TGE estimator is appropriate, and the validation of source astrometry and flux against GLEAM for one pointing is a positive feature. However, the scientific impact is currently limited by the fact that only 6 of 24 pointing centers yield a fit, and the interpretation of the residual spectrum as genuine DGSE has not been separated from systematic effects. The paper is candid about the Gaussian-random-field assumption for error bars and the strong parameter degeneracies, but it does not provide the null tests needed to support the central claim.","major_comments":[{"comment":"The conclusion that the shallow component at l<200 is genuine DGSE is not supported by a null test. The source-subtraction validation in Appendix 1 reports a systematic ~30 arcsec DEC offset and flux deviations up to 25% for some sources, and Section 4 attributes the flat C_l of unsuccessful PCs to imaging artifacts and residual gain calibration errors. Since no simulation of a point-source-only sky (with sources below the 430 mJy threshold plus noise) is presented, an imperfectly subtracted bright source could produce a non-white residual that mimics the A(1000/l)^beta term in eq. (4). The authors should demonstrate that the observed residual C_l exceeds the expectation from a point-source-only model with realistic subtraction errors.","section":"Section 3 (after Fig. 1) and Appendix 1"},{"comment":"The strong degeneracies among A, beta, and C (Pearson |r| ~ 0.9) together with the limited l-range (65<l<650, roughly one decade) mean that the individual parameters may not be independently constrained. The posterior distributions shown in Appendix 2 are highly asymmetric and often truncated at the prior boundary, indicating weak constraints. The paper should assess the robustness of the reported beta range, for example by fixing C to the expected Poisson level from the source count model, by using profile likelihoods, or by checking whether the fit is stable when the fitted l-range is varied.","section":"Section 3 (Fig. 3 and Table 1)"},{"comment":"The criterion for being 'able to fit the model' is not defined a priori, and no results are shown for the remaining 18 PCs beyond one representative example (Fig. 4). If the six PCs are selected after inspecting the fit quality, the reported beta range may not be representative of the region and any inference from the scatter of parameters could be biased. The authors should either report the measured C_l for all 24 PCs (including upper limits for the unsuccessful ones) or provide a principled, pre-defined selection rule and justify that the selected subset is unbiased.","section":"Section 3 (selection of the six PCs)"}],"minor_comments":[{"comment":"The caption states that black lines show the total data before point-source subtraction and red lines show C_l after removal, but Section 3 states that red and black lines show the 'No-Sub' and 'UV-Sub' cases respectively; these are opposite and should be corrected.","section":"Figure 1 caption"},{"comment":"The abstract and Figure 2 use the fitting range '65 <= l <= 650', while Section 3 uses '65 < l < 650'; please make these consistent.","section":"Abstract and Section 3"},{"comment":"The abstract states that A varies from 155 to 400 mK^2, but Table 1 and Section 4 report values up to 568.4 mK^2; please correct the abstract.","section":"Abstract"},{"comment":"The phrase 'we consider an angular scale of 0.4 × theta_FWHM' is ambiguous; it likely means a radial cut at 0.4 times the primary beam FWHM, but this should be stated explicitly.","section":"Appendix 1"},{"comment":"The measured C_l values are not provided in a table or as electronic supplementary material, which would facilitate comparison with future measurements and independent re-analysis.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of PASA and uses a standard estimator, but the central claim rests on only 6 out of 24 pointing centers without a null test that separates DGSE from source-subtraction systematics. The authors report the parameter degeneracy but do not address its impact on the claimed beta range. I recommend requiring additional simulations (e.g., point-source-only skies, variation of the subtraction threshold) before the paper can be accepted. The figure-caption error and abstract/table inconsistencies are easily fixed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new content here is straightforward: the authors use the TGE on MWA drift-scan data to measure C_l of diffuse Galactic synchrotron emission after subtracting point sources above 430 mJy, for six pointing centers in the EoR fields. The numerical C_l values are new, and the paper does a reasonable job of validating the source subtraction against GLEAM for one pointing, calibrating TGE with UAPS simulations, and presenting honest MCMC posteriors that show how degenerate A, beta, and C really are. The result is a consistency check with earlier 150 MHz measurements, and the beta values (0.9-1.7) are statistically consistent with those, as the authors say.\n\nWhere I worry is the central inference that the residual C_l at l<200 is genuinely DGSE. The stress-test note is on target: there is no null simulation with only point sources below the subtraction threshold plus noise. Without that, the shallow component at low l could be produced by source-subtraction residuals or gain errors instead of diffuse emission. The paper itself strengthens this concern by attributing the flat C_l in the other 18 PCs to 'imaging artifacts around bright sources' and residual calibration errors, and the appendix reports a 30 arcsec DEC offset and flux deviations up to 25%. Those systematics are exactly the kind that can inject a shallow-spectrum term and mimic the A(1000/l)^beta component, especially given the strong parameter correlations (|r| ~ 0.9). So the claim that beta is a property of the DGSE is plausible but not demonstrated.\n\nThe other soft spots are the selection of only 6 of 24 PCs, all clustered in a narrow RA range, and the Gaussian random field assumption for error bars, which the authors admit is 'quite likely' invalid. The lack of a quantitative comparison with prior measurements is a minor gap; 'consistent at 2 sigma' does not tell the reader whether the new data favor a flatter or steeper spectrum.\n\nThat said, the paper is transparent, the pipeline is standard, and the measurements are new numbers for a region of interest to EoR foreground work. It deserves a serious referee, but the referee should ask for the missing null test and a quantitative comparison with Bernardi et al., Ghosh et al., and Iacobelli et al. before acceptance. I would not let the paper through as is; I would recommend major revision.\n\nFor a reading group, I'd take it along, mostly to discuss how one separates genuine diffuse emission from subtraction systematics in this kind of data. I would not cite it in my own work until the null test is done.","headline":"Useful incremental MWA measurement of the DGSE angular power spectrum at 154 MHz, but the quoted beta values rest on an untested assumption that residual systematics are negligible.","tokens_in":16418,"tokens_out":2398,"would_cite":false,"duration_ms":28886,"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":"Galactic synchrotron spectrum at 154 MHz has beta 0.9–1.7","keywords":["Galactic synchrotron emission","angular power spectrum","Murchison Widefield Array","drift scan","Tapered Gridded Estimator","21-cm foregrounds","point source subtraction","epoch of reionization"],"falsifier":"A decisive test would be to re-observe the same six fields with deeper integration and a lower source-subtraction threshold (e.g., sources above 100 mJy instead of 430 mJy) and check whether the fitted power-law parameters $A$ and $\\beta$ remain within the reported ranges while the constant term $C$ drops proportionally; if the slope or amplitude changes substantially, the claimed diffuse-emission interpretation is contaminated by systematics or unresolved sources.","tokens_in":15379,"feed_emoji":"📡","tokens_out":13615,"duration_ms":113997,"temperature":0.7,"pith_summary":"This paper reports a measurement of the angular power spectrum of diffuse Galactic synchrotron emission (DGSE) at 154.2 MHz using Murchison Widefield Array drift-scan observations along a fixed declination strip. After removing point sources brighter than about 430 mJy from 24 pointing centers, the authors fit the residual power spectrum with a combination of a power law from the diffuse emission and a constant from unresolved point-source Poisson fluctuations. For six of the pointings the fit works, giving power-law slopes $\\beta$ between 0.9 and 1.7 and amplitudes $A$ between about 154 and 568 $\\mathrm{mK}^2$. The slopes are consistent at the $2\\sigma$ level with earlier 150 MHz measurements, supporting the use of such measurements as foreground models for 21-cm epoch-of-reionization experiments and as probes of interstellar turbulence.","feed_headline":"Galactic synchrotron spectrum at 154 MHz has beta 0.9–1.7","feed_subtitle":"These six power-law fits refine 21-cm foreground models and ISM turbulence probes.","key_machinery":"The central tool is the Tapered Gridded Estimator (TGE), a visibility-based estimator that grids the measured visibilities, applies a Gaussian tapering window with $15^\\circ$ FWHM to suppress primary-beam sidelobes, subtracts the noise bias from the self-correlation of visibilities, and normalizes the estimate using simulated unit angular power spectra. The model being fitted is $C_\\ell^M = A (1000/\\ell)^\\beta + C$, where the power-law term represents the diffuse Galactic synchrotron emission and the constant $C$ represents Poisson fluctuations from point sources below the 430 mJy subtraction threshold. Source subtraction is performed by imaging with long baselines ($|u|>50\\lambda$), CLEAN modeling, and subtracting the model visibilities; the analysis is restricted to $65 < \\ell < 650$, where primary-beam convolution effects and noise-dominated high-$\\ell$ bins are avoided.","core_discovery":"The central claim is that the residual angular power spectrum of diffuse Galactic synchrotron emission, measured with the Tapered Gridded Estimator on MWA drift-scan visibilities after subtracting compact sources above $\\sim 430$ mJy, is described in six fields by the model $C_\\ell^M = A (1000/\\ell)^\\beta + C$ over the multipole range $65 < \\ell < 650$. The power-law component is attributed to diffuse Galactic synchrotron emission, the constant to Poisson fluctuations of residual point sources, and the fitted parameters vary with sky position: $A$ from about 154 to 568 $\\mathrm{mK}^2$, $\\beta$ from 0.9 to 1.7, and $C$ from 777 to 4457 $\\mathrm{mK}^2$. The paper argues that the shallow slope at $\\ell \\le 200$, in contrast to the $\\ell^2$ scaling of point-source Poisson noise at $\\ell > 200$, indicates genuine diffuse emission rather than residual sources. The measured $\\beta$ values are consistent at $2\\sigma$ with earlier interferometric measurements at similar frequencies and angular scales.","pith_inferences":["If the measured $\\beta$ range (0.9 to 1.7) is representative of diffuse Galactic synchrotron emission at arc-minute scales, then the foreground is much flatter at $\\ell \\sim 1000$ than the steeper slopes found at degree scales, and locating the break between these regimes would help separate turbulent from Poissonian contributions in 21-cm foreground subtraction.","The constant term $C$ of the fit is effectively a measurement of residual point-source confusion power; comparing $C$ with predictions from source count models at 154 MHz could validate the completeness of the 430 mJy catalog and the accuracy of the source subtraction.","The 18 non-fitted pointings, including the Fornax A region, suggest that bright extended sources and calibration artifacts dominate the residuals, so a re-analysis after improved ionospheric calibration or with the full MWA Phase II configuration could test whether more pointings become fit.","Extending this single-frequency measurement to a wide bandwidth could test the double power-law frequency dependence of the amplitude reported in spectral studies, linking angular and spectral foreground models."],"forward_implications":["The six fitted power spectra provide a foreground model for epoch-of-reionization 21-cm observations in the MWA EoR0/EoR1 region, enabling better foreground avoidance or removal at arc-minute scales.","The pointing-to-pointing variation in $A$ and $\\beta$ constrains fluctuations of the Galactic magnetic field and cosmic-ray electron density in this sky region.","Consistency of $\\beta$ with earlier GMRT and WSRT measurements at 150 MHz strengthens confidence in visibility-based power-spectrum estimation after aggressive source subtraction.","For the 18 pointings that do not fit, the flat residual spectrum indicates that unresolved Poisson sources or imaging artifacts limit shallow drift-scan measurements, pointing to the need for deeper observations."],"supporting_citations":[{"why":"Supplies the MWA drift-scan visibility data (project G0031) analyzed here.","marker":"Patwa et al. (2021)"},{"why":"Provides the Tapered Gridded Estimator (TGE) formalism used to estimate $C_\\ell$ from gridded visibilities.","marker":"Choudhuri et al. (2016)"},{"why":"Validates TGE on simulated MWA observations and establishes that $\\ell \\ge 65$ is unaffected by primary-beam and tapering convolution.","marker":"Chatterjee et al. (2022)"},{"why":"Gives the model that point-source Poisson fluctuations produce $D_\\ell \\propto \\ell^2$ and a constant $C$, used to split the residual spectrum.","marker":"Ali et al. (2008)"},{"why":"Previous WSRT 150 MHz DGSE power spectrum measurement used as a comparison for the fitted $\\beta$ values.","marker":"Bernardi et al. (2009)"},{"why":"Earlier GMRT 150 MHz measurement with $\\beta \\approx 2.34$, used as a consistency check.","marker":"Ghosh et al. (2012)"},{"why":"Earlier arc-minute-scale measurement with $\\beta \\approx 1.8$, consistent with the reported range.","marker":"Iacobelli et al. (2013a)"},{"why":"Earlier TGE-based GMRT measurement at 150 MHz cited as consistent with the fitted $\\beta$ at $2\\sigma$.","marker":"Choudhuri et al. (2017)"},{"why":"TGSS-based all-sky $C_\\ell$ study showing similar amplitude variation with pointing, used to interpret $A$'s sky dependence.","marker":"Choudhuri et al. (2020)"},{"why":"Supplies the Markov Chain Monte Carlo ensemble sampler used for parameter estimation.","marker":"Foreman-Mackey et al. (2013)"}],"fun_headline_variants":["MWA drift scans pin Galactic synchrotron beta 0.9–1.7","Beta 0.9–1.7 for Galactic synchrotron from MWA drift scans","Synchrotron slope beta 0.9–1.7 at 154 MHz from MWA","MWA drift scan power spectrum yields beta 0.9–1.7"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the residual visibilities after subtracting point sources are genuine diffuse Galactic synchrotron emission, not leftover calibration errors, ionospheric distortions, or unmodeled compact sources.","fun_headline_variants_meta":{"raw":{"variants":["MWA drift scans pin Galactic synchrotron beta 0.9–1.7","Beta 0.9–1.7 for Galactic synchrotron from MWA drift scans","Synchrotron slope beta 0.9–1.7 at 154 MHz from MWA","MWA drift scan power spectrum yields beta 0.9–1.7"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001291,"raw_usage":{"total_tokens":5423,"prompt_tokens":1251,"completion_tokens":4172,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":867,"completion_tokens_details":{"reasoning_tokens":4073}},"tokens_in":867,"tokens_out":4172,"duration_ms":27421,"temperature":1.0,"reasoning_tokens":4073,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:17:18.635608+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be to re-observe the same six fields with deeper integration and a lower source-subtraction threshold (e.g., sources above 100 mJy instead of 430 mJy) and check whether the fitted power-law parameters $A$ and $\\beta$ remain within the reported ranges while the constant term $C$ drops proportionally; if the slope or amplitude changes substantially, the claimed diffuse-emission interpretation is contaminated by systematics or unresolved sources.","supporting_citations":[{"cited_title":"K., Sethi, S., & Dwarakanath, K","cited_arxiv_id":null,"evidence_quote":"Supplies the MWA drift-scan visibility data (project G0031) analyzed here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Validates TGE on simulated MWA observations and establishes that $\\ell \\ge 65$ is unaffected by primary-beam and tapering convolution."},{"cited_title":"S., Bharadwaj , S., & Chengalur , J","cited_arxiv_id":null,"evidence_quote":"Gives the model that point-source Poisson fluctuations produce $D_\\ell \\propto \\ell^2$ and a constant $C$, used to split the residual spectrum."},{"cited_title":"G., Brentjens , M","cited_arxiv_id":null,"evidence_quote":"Previous WSRT 150 MHz DGSE power spectrum measurement used as a comparison for the fitted $\\beta$ values."},{"cited_title":"S., & Chengalur , J","cited_arxiv_id":null,"evidence_quote":"Earlier GMRT 150 MHz measurement with $\\beta \\approx 2.34$, used as a consistency check."},{"cited_title":"S., et al","cited_arxiv_id":null,"evidence_quote":"Earlier TGE-based GMRT measurement at 150 MHz cited as consistent with the fitted $\\beta$ at $2\\sigma$."},{"cited_title":"2020, , 494, 1936","cited_arxiv_id":null,"evidence_quote":"TGSS-based all-sky $C_\\ell$ study showing similar amplitude variation with pointing, used to interpret $A$'s sky dependence."}],"review_version":1}