{"id":"2a225625-ec74-44a5-b9cf-060049393c9b","arxiv_id":"2412.03163","paper_version":2,"verdict":"REJECT","confidence":"HIGH","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":10,"one_line_summary":"Fitting LMC low-frequency radio data with a dark matter plus cosmic-ray double power-law model yields dark-matter flux limits of 98.4 to 114.8 Jy at 1.4 GHz, which translate to annihilation cross-section bounds above the thermal relic value.","lead":"This paper fits radio observations of the Large Magellanic Cloud from 19.7 MHz to 1.4 GHz with a two-component spectrum, separating cosmic-ray emission from a possible dark matter annihilation signal, and converts the allowed dark matter component into upper limits on dark matter mass and annihilation rate.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The DM component in the fit is driven by the 19.7/45 MHz historical points, which are inconsistent with GLEAM; a GLEAM-only refit would likely eliminate the claimed constraint.","rationale":"The reader's REJECT is well-founded. Among the listed defects, the most load-bearing for the paper's central claim is the reliance of the DM normalization on the two historical low-frequency data points that are mutually inconsistent with GLEAM at overlapping frequencies. This is not a subtle statistical issue: with αCR fit to 0.40–0.52 and αDM fixed at 0.75, the DM component is the only steep component that can bridge the flat GLEAM spectrum and the high 19.7/45 MHz fluxes. If those old measurements are biased by confusion, ionospheric calibration, or beam-size mismatch, the 98–115 Jy normalization is a systematic artifact, and the exclusion curves in Section III.D do not follow. The 'best-fit instead of statistical upper limit' issue flagged by the reader is real but secondary: recomputing a proper 95% upper limit from the same biased data would still produce an over-strong constraint. The fixed-index choice is theoretically motivated, but the v1 free-index result (0.21–0.66) shows the decomposition is degenerate; the historical-data test is the cleanest way to expose that degeneracy. The proposed GLEAM-only refit is decisive, inexpensive, and directly tests whether the paper's quantitative conclusions survive without the two questionable points. I therefore see no change to the reader's REJECT verdict.","tokens_in":13602,"tokens_out":9423,"duration_ms":86241,"concrete_test":"Refit Eq. (3) with the same fixed α_DM = 0.75 and the same MCMC setup, but using only the GLEAM/For et al. points at 76–277 MHz plus the 408 MHz and 1.4 GHz measurements, excluding the 19.7 and 45 MHz historical points. Compare the resulting posterior for SDM(1.4 GHz) with the 98–115 Jy values used in Fig. 4. If the 95% credible interval includes zero or the best-fit normalization drops to a small fraction of 98 Jy, the claimed DM constraints are driven by the excluded historical data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Table I, the two pre-GLEAM points at 19.7 MHz (5270±1054 Jy) and 45 MHz (2997±450 Jy) are the only data below 76 MHz, and they are much higher than a smooth extrapolation of the GLEAM points. At overlapping or nearby frequencies the disagreement is large: 98.6 MHz gives 2839±600 Jy versus GLEAM 99 MHz at 1451.6±247.0 Jy, and 158 MHz gives 1736±490 Jy versus GLEAM 158 MHz at 1350.4±229.6 Jy. Because the fitted CR spectral index is only 0.40–0.52 while the DM index is fixed to 0.75, the steep DM component in Eq. (3) is the only element in the fit that can account for the high 19.7/45 MHz fluxes. Thus the best-fit SDM = 98.4–114.8 Jy at 1.4 GHz is effectively a fit to the systematic discrepancy of these historical measurements, not a robust limit on DM-induced synchrotron emission. The v1 abstract's free-index result α_DM = 0.21–0.66 reinforces the degeneracy: the data do not independently select a 0.75 DM component. No robustness test against the GLEAM-only dataset or treatment of the historical data's systematic errors is presented. If those two points are biased high, the mχ–⟨σv⟩ exclusion curves in Figs. 4–6 do not follow.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper fits the low-frequency radio spectrum of the Large Magellanic Cloud (19.7 MHz to 1.4 GHz) with a double power-law model consisting of a dark-matter synchrotron component with a fixed spectral index (α_DM = 0.75) and a cosmic-ray synchrotron component with free normalization and spectral index. Using MCMC, the authors obtain best-fit values of S_DM(1.4 GHz) = 114.8 Jy (without thermal emission) and 98.4 Jy (with a thermal component), which they call upper limits on dark-matter-induced synchrotron emission. They then use RX-DMFIT to translate these flux limits into constraints on the dark matter annihilation cross section ⟨σv⟩ as a function of mχ, showing curves for several frequencies and for variations in the diffusion coefficient and magnetic field strength. The central claim is that the LMC radio data exclude cross sections above roughly 10^-23 to 10^-21 cm^3/s for masses 10 to 1000 GeV, with the lowest frequencies most sensitive to low dark matter masses.","tokens_in":13977,"tokens_out":11201,"duration_ms":106920,"significance":"The paper addresses an interesting and timely question: whether low-frequency radio observations of a nearby, massive dwarf galaxy can constrain dark matter annihilation. The use of the updated GLEAM data, the explicit modeling of both cosmic-ray and dark-matter synchrotron components, and the exploration of sensitivity to diffusion and magnetic field parameters are sensible first steps, and the paper makes a falsifiable prediction that future low-frequency surveys could detect a steep dark-matter component. However, the central analysis has two load-bearing problems. First, the quoted 'upper limits' on S_DM are actually MCMC best-fit values, not statistical upper bounds, so the exclusion curves in Figs. 4–6 do not have the claimed meaning. Second, the dark-matter component is driven by pre-GLEAM historical flux measurements that are inconsistent with GLEAM at overlapping frequencies; a GLEAM-only fit could plausibly eliminate the dark-matter component altogether. The fixed α_DM = 0.75 also contradicts the paper's own v1 abstract, which reported a free-index fit of α_DM = 0.21–0.66. Because these issues affect the main result, the paper is not acceptable in its present form.","major_comments":[{"comment":"The quantities quoted as 'upper limits' (S_DM = 114.8 Jy and 98.4 Jy at 1.4 GHz) are the best-fit values from the MCMC, not upper limits. The text reports log10 S_DM = 2.06^{+0.17}_{-0.14} and 1.99^{+0.20}_{-0.17}; these are central estimates with 1σ intervals. A proper upper limit must be a quantile of the marginalized posterior, such as a 95% credibility bound, and it must allow for the possibility that the data are consistent with S_DM = 0. Because the exclusion curves in Figs. 4–6 are computed from the single best-fit normalization, they do not have the statistical meaning claimed in Section III.D. This is a load-bearing error: the paper's main result is a set of upper limits, and those limits are not actually computed.","section":"Section II.B, Figs. 1–2"},{"comment":"The four pre-GLEAM points below or overlapping the GLEAM band (19.7, 45, 85.5, 98.6, and 158 MHz) are inconsistent with GLEAM measurements at the same or nearby frequencies. For example, 98.6 MHz gives 2839±600 Jy versus 1451.6±247.0 Jy at 99 MHz, and 158 MHz gives 1736±490 Jy versus 1350.4±229.6 Jy. Because the fixed dark-matter spectral index (α_DM = 0.75) is steeper than the fitted cosmic-ray index (α_CR = 0.40–0.52), the dark-matter component is the only element of the model that can absorb the excess flux of the historical points. Consequently, the best-fit S_DM is effectively a fit to the systematic offset of the historical data, not a robust constraint on dark-matter-induced synchrotron emission. A fit restricted to the GLEAM data, or a quantitative treatment of the historical data's systematics, is required before any dark-matter limit can be claimed.","section":"Section II.A and Table I"},{"comment":"The fixed value α_DM = 0.75 is adopted from Tasitsiomi et al. [24], but the previous version of this paper on arXiv reported a free-α_DM fit with α_DM = 0.21–0.66. The current manuscript does not mention or discuss this discrepancy. Since a flatter dark-matter spectrum would substantially reduce the need for a dark-matter component at low frequencies, the choice α_DM = 0.75 is not a harmless convention; it is one of the main determinants of the derived constraints. The authors should either justify the fixed value with the data or treat α_DM as a free parameter and show how the limits depend on it.","section":"Section II.B, Eq. (3)"},{"comment":"The likelihood uses only the quoted 1σ statistical errors of the individual flux measurements. The large scatter between overlapping historical and GLEAM points indicates significant unmodeled systematics, so the MCMC error bars are underestimated. The paper does not include a systematic error floor or a covariance between measurements, and it does not test the robustness of the fit to excluding the pre-GLEAM points. The resulting limits are therefore overconfident even if the best-fit-to-upper-limit issue in the first major comment were fixed.","section":"Section II.B, Eq. (5)"}],"minor_comments":[{"comment":"Equation (2) is a polynomial in x, not a pure power law, yet the text states that it leads to a ν^-0.75 power-law dependence. Please state the frequency and energy range over which this approximation is valid and show numerically that it holds across the full 19.7 MHz to 1.4 GHz band and the adopted dark matter mass range.","section":"Section II.B, Eq. (2)"},{"comment":"The 1400 MHz entry cited to For et al. [23] appears to be outside the GLEAM frequency range (76–227 MHz); please verify the provenance of this data point. In addition, the two 1400 MHz measurements (384±30 Jy and 529±30 Jy) differ by nearly 40%, and the paper does not discuss how this systematic difference is handled.","section":"Table I"},{"comment":"The red lines in Figs. 1 and 2 are labeled 'upper limit' in the captions, but they are plotted from the best-fit model. The labels should be corrected to 'best fit' or the plots should show the actual upper-limit curves derived from the posterior.","section":"Section II.B, Figs. 1–2"},{"comment":"The paper does not state the priors used for the MCMC parameters, the chain lengths, or convergence diagnostics. This information is necessary for reproducibility, especially since the authors use interpolation to compute model fluxes.","section":"Section II.B"}],"recommendation":"reject","confidential_remarks":"The previous version of this paper reported α_DM = 0.21–0.66 from a free-index fit, while the current version fixes α_DM = 0.75 without discussing the change; this should be disclosed. The central constraints are not established by the analysis as presented: the 'upper limits' are best-fit values, and the dark-matter component is driven by historical data inconsistent with GLEAM. A new analysis, rather than a minor revision, would be needed, and if the pre-GLEAM points are excluded the claimed dark-matter limits may not survive."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something genuinely useful: it pulls together the GLEAM low-frequency data and the older historical measurements, fits a two-component (DM plus CR) synchrotron model, and then runs the fitted DM flux through RX-DMFIT to get mχ–⟨σv⟩ contours. The modeling is transparent, the authors check D0 and B0 dependence, and they are honest that their limits sit above the thermal relic value and below Fermi-LAT. That part is fine, and the updated data compilation alone has some value.\n\nBut the central quantitative claim does not hold up. Three problems, in order of severity. First, the 'upper limits' on S_DM are simply the MCMC best-fit values (log10 S_DM = 2.06 and 1.99), with uncertainties quoted but then discarded. A best-fit value is not an upper limit; treating it as one makes the constraints artificially strong, not conservative. This is a basic statistical error that propagates directly into the exclusion curves.\n\nSecond, the DM spectral index is fixed to α_DM = 0.75, citing Tasitsiomi et al. But the paper's own v1 abstract reported a free fit giving α_DM = 0.21–0.66. The fact that the free index is much flatter than 0.75 means the data do not independently prefer the steep DM component; the decomposition is degenerate. The fixed index is load-bearing, not a detail.\n\nThird, and most damning, the low-frequency historical points are inconsistent with GLEAM at nearly the same frequencies. The 98.6 MHz point (2839±600 Jy) sits roughly 2σ above GLEAM's 99 MHz point (1451.6±247.0 Jy), and similar tensions appear at 158 MHz. Those high historical points — 19.7 and 45 MHz especially — are what force the steep DM component into the fit. With a CR index of only 0.40–0.52, the ν^-0.75 DM term is the only way to reproduce the 19.7 MHz flux. If those historical measurements are biased high — which the GLEAM comparison suggests — the S_DM values and all subsequent contours evaporate. The paper presents no GLEAM-only robustness test and no treatment of the historical data's systematics.\n\nSo the paper is not publishable as is. The qualitative idea — low-frequency radio can in principle probe low-mass WIMPs — is plausible and already in the literature. What is new here is the specific numerical result, and that result is not reliable. With proper upper limits, a free α_DM fit, and a GLEAM-only cross-check, the analysis could become a solid null result or even a meaningful constraint. As it stands, I would send it back. It deserves referee time only if the authors are willing to redo the statistics and confront the data inconsistency.","headline":"A well-intentioned but statistically flawed attempt to constrain LMC dark matter with low-frequency radio — the claimed limits are driven by suspect historical data points.","tokens_in":747,"tokens_out":807,"would_cite":false,"duration_ms":32266,"reading_group":"maybe","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 the Large Magellanic Cloud place upper limits on dark matter annihilation, excluding cross sections above roughly $10^{-23}$ to $10^{-21}$ cm$^3$ s$^{-1}$ for masses of 10–1000 GeV.","keywords":["dark matter annihilation","Large Magellanic Cloud","synchrotron radiation","low-frequency radio astronomy","indirect dark matter detection","gaugino annihilation","cosmic rays","radio-infrared correlation"],"falsifier":"Rerun the identical fit to the same 19.7 MHz–1.4 GHz data with $\\alpha_{\\rm DM}$ as a free parameter instead of fixing it at 0.75; if the best-fit $\\alpha_{\\rm DM}$ comes out well below 0.75 with a comparable likelihood, the fixed-slope upper limits are not uniquely determined and the cross-section constraints would need revision.","tokens_in":13294,"feed_emoji":"📡","tokens_out":13882,"duration_ms":114554,"temperature":0.7,"pith_summary":"Low-frequency radio emission from the Large Magellanic Cloud carries a faint, mostly nonthermal signal that this paper uses to bound dark matter annihilation. The authors fit 19.7 MHz–1.4 GHz flux measurements with two power-law components: cosmic-ray synchrotron radiation with a free spectral index, and a dark-matter-annihilation component whose spectral index is fixed at $-0.75$ from gaugino-annihilation physics. Treating the fitted dark-matter normalization as an upper limit gives $S_{\\rm DM}(1.4\\,{\\rm GHz}) = 98$--$115$ Jy depending on whether thermal free-free emission is included. Converting that flux limit into particle-physics terms with a diffusion-loss model for $e^\\pm$ yields exclusion curves in the dark-matter mass vs. annihilation-cross-section plane. If the paper is right, the LMC radio data exclude cross sections above roughly $10^{-23}$ to $10^{-21}$ cm$^3$ s$^{-1}$ for dark-matter masses of 10–1000 GeV, with the lowest frequencies probing the lowest masses.","feed_headline":"LMC radio glow sets new dark matter annihilation limits","feed_subtitle":"The limits reach about 10^-23 cm^3/s at low mass, with the lowest frequencies probing deepest.","key_machinery":"The double power-law model $S_{\\rm nth} = S_{\\rm DM}(\\nu/\\nu_\\star)^{-0.75} + S_{\\rm CR}(\\nu/\\nu_\\star)^{-\\alpha_{\\rm CR}}$, with $\\nu_\\star = 1.4$ GHz, is the core object: it separates the radio spectrum into a dark-matter piece with a fixed spectral index and a cosmic-ray piece with a free index. A Markov-chain Monte Carlo fit over $S_{\\rm DM}$, $S_{\\rm CR}$, and $\\alpha_{\\rm CR}$ (plus a thermal free-free component in one variant) returns the normalization that is then read as an upper limit on DM-induced synchrotron emission. The second piece of machinery is the steady-state transport equation for $e^\\pm$ with diffusion and energy losses, solved with a free-escape boundary at 3.5 kpc, which turns the flux limit into a predicted flux for each value of $m_\\chi$ and $\\langle\\sigma v\\rangle$. Comparing predicted to allowed flux over a grid of masses and cross sections draws the exclusion curves.","core_discovery":"The paper's claim is that a two-component spectral decomposition of the LMC's low-frequency radio flux isolates an upper limit on dark-matter-annihilation synchrotron emission. With the dark-matter index fixed at $\\alpha_{\\rm DM}=0.75$, the best-fit dark-matter normalization at 1.4 GHz is 114.8 Jy without a thermal component and 98.4 Jy with one, and the cosmic-ray component comes out flatter ($\\alpha_{\\rm CR}\\approx0.40$--$0.52$) than the canonical 0.8. These normalizations bound the flux a dark-matter signal could contribute at every frequency, because any excess above the fitted CR component is attributed to DM. Using an analytic diffusion-loss Green's function and a synchrotron emissivity calculation, the paper converts these bounds into exclusion curves for $m_\\chi$ vs. $\\langle\\sigma v\\rangle$, finding that lower frequencies give stronger limits on lower-mass dark matter, weaker diffusion gives stronger limits, and stronger magnetic fields give stronger limits.","pith_inferences":["Letting $\\alpha_{\\rm DM}$ float in the same fit would test the fixed $-0.75$ slope; a flatter best-fit index would mean the quoted cross-section limits are an artifact of the assumed spectrum.","A modern re-measurement of the 19.7 and 45 MHz fluxes would check the steep low-frequency excess that drives the DM-limited low-mass constraints, since those historical points are the ones that push the fitted DM component upward.","The same two-component decomposition could be applied to other dwarf irregular galaxies with low-frequency spectra, potentially producing stacked dark-matter limits across a population."],"forward_implications":["For dark-matter masses between 10 and 1000 GeV, the LMC radio limits exclude annihilation cross sections above roughly $10^{-23}$ to $10^{-21}$ cm$^3$ s$^{-1}$.","Including thermal free-free emission tightens the dark-matter normalization by about 17%, from 114.8 Jy to 98.4 Jy at 1.4 GHz.","The fitted cosmic-ray spectral index ($\\alpha_{\\rm CR}\\approx0.4$--$0.5$) is flatter than the canonical 0.8, implying the LMC's radio spectrum is flatter than normal galaxies and that radio-IR correlations calibrated on normal galaxies underpredict the cosmic-ray contribution.","Lower-frequency radio bands are the most sensitive probe of low-mass dark matter, so future very-low-frequency surveys could push the constraints below the current values."],"supporting_citations":[{"why":"Supplies the 76–227 MHz flux measurements that form the bulk of the low-frequency data, and the previous best-fit $\\alpha_{\\rm CR}=0.55$ that motivates leaving $\\alpha_{\\rm CR}$ free.","marker":"[23]"},{"why":"Establishes the $\\nu^{-0.75}$ synchrotron spectrum of gaugino-annihilation $e^\\pm$ that fixes $\\alpha_{\\rm DM}$ in the fit.","marker":"[24]"},{"why":"Derives the steady-state $e^\\pm$ transport equation with diffusion and energy losses, and one of the analytic solutions used to propagate the injected spectrum.","marker":"[37]"},{"why":"Provides the Green's-function solution with free-escape boundary conditions used to compute the equilibrium $e^\\pm$ distribution in the LMC.","marker":"[41]"},{"why":"The prior low-frequency search in the same target whose upper limit this work extends with additional data below 1 GHz.","marker":"[16]"},{"why":"Earlier constraints derived from higher-frequency LMC radio emission that this work extends downward in frequency.","marker":"[18]"},{"why":"The gamma-ray limits to which the radio-derived cross-section bounds are compared and found weaker.","marker":"[44]"}],"fun_headline_variants":["Low-frequency radio tightens dark matter limits","LMC radio data sharpen dark matter annihilation bounds","Lower radio frequencies probe deeper dark matter","Dark matter annihilation constrained by LMC glow","Radio eyes on LMC reveal dark matter ceiling"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the LMC's nonthermal radio spectrum is exactly the sum of a single cosmic-ray power law and a dark-matter component with a fixed $\\nu^{-0.75}$ slope; if the dark-matter spectrum is flatter or the cosmic-ray component is not a single power law, the derived cross-section limits do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Low-frequency radio tightens dark matter limits","LMC radio data sharpen dark matter annihilation bounds","Lower radio frequencies probe deeper dark matter","Dark matter annihilation constrained by LMC glow","Radio eyes on LMC reveal dark matter ceiling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000202,"raw_usage":{"total_tokens":1483,"prompt_tokens":1144,"completion_tokens":339,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":760,"completion_tokens_details":{"reasoning_tokens":272}},"tokens_in":760,"tokens_out":339,"duration_ms":3926,"temperature":1.0,"reasoning_tokens":272,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T22:42:54.068830+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Rerun the identical fit to the same 19.7 MHz–1.4 GHz data with $\\alpha_{\\rm DM}$ as a free parameter instead of fixing it at 0.75; if the best-fit $\\alpha_{\\rm DM}$ comes out well below 0.75 with a comparable likelihood, the fixed-slope upper limits are not uniquely determined and the cross-section constraints would need revision.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the $\\nu^{-0.75}$ synchrotron spectrum of gaugino-annihilation $e^\\pm$ that fixes $\\alpha_{\\rm DM}$ in the fit."},{"cited_title":"& Ullio, P","cited_arxiv_id":null,"evidence_quote":"Derives the steady-state $e^\\pm$ transport equation with diffusion and energy losses, and one of the analytic solutions used to propagate the injected spectrum."},{"cited_title":"2017, JCAP, 9, 027","cited_arxiv_id":null,"evidence_quote":"Provides the Green's-function solution with free-escape boundary conditions used to compute the equilibrium $e^\\pm$ distribution in the LMC."},{"cited_title":"2015, Physical Review Letters, 115, 231301","cited_arxiv_id":null,"evidence_quote":"The gamma-ray limits to which the radio-derived cross-section bounds are compared and found weaker."}],"review_version":1}