{"id":"8343ef45-1da0-43b9-becc-7d1737b48d52","arxiv_id":"2501.01489","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Dark inverse Compton scattering of cosmic-ray electrons and positrons on ultralight dark photon dark matter produces a sub-MHz radio background that, compared with IMP-6, RAE-2, and Parker Solar Probe data, constrains the kinetic mixing below about 2e-6 for masses below 2e-17 eV.","lead":"A new calculation shows that ultralight dark photons scattering off cosmic-ray electrons would produce a faint radio glow below 1 MHz. Comparing this predicted glow to spacecraft radio surveys, the authors set some of the strongest limits yet on dark photon dark matter in this mass range.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Absorption-model uncertainty at 0.1–0.5 MHz is exponential, and the quoted O(1) variation in the epsilon limits may underestimate the shift; the 'leading constraints' claim needs an uncertainty band.","rationale":"I read the paper as a conservative upper-limit analysis: the dark Compton cross-section in Appendix A reproduces the standard inverse-Compton spectrum in the Thomson limit, the cosmic-ray fluxes are anchored by AMS-02/Voyager, and the inequality limit procedure avoids background fitting. The central derivation appears internally consistent. The weakest link is the absorption optical depth in Eq. (8), which controls the normalization of the sub-MHz signal. At 0.1 MHz the optical depth reaches unity within ~14 pc for nominal YMW16 parameters, so the predicted flux is set by the least well-mapped part of the free-electron distribution, and the exponential sensitivity amplifies the model uncertainty far more than the linear 50% statement suggests. A concrete propagation of the YMW16 parameter uncertainties and the local temperature range would settle whether the constraints remain leading. This is an addressable systematic rather than a fatal flaw, so I agree with the CONDITIONAL verdict and recommend verifying with an uncertainty band.","tokens_in":13064,"tokens_out":50881,"duration_ms":462534,"concrete_test":"Recompute the exclusion curves in Fig. 4 with (i) n_e globally rescaled by 0.5 and 2.0 in Eq. (8), and (ii) T = 5e3 K and 2e4 K, producing a shaded band for each dataset. If the band at every mA' <~ 2e-17 eV lies below the SuperMAG and other astrophysical constraints plotted in Fig. 4, the central claim survives; if the band crosses those constraints, the 'leading constraints' claim needs qualification or the conditions under which it holds must be stated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is the free-free absorption factor P_nu = exp(-tau_nu) in Eq. (8), which sets the normalization of the predicted sub-MHz flux. At the constraining frequencies 0.1–0.5 MHz, the optical depth per parsec is tau/s = 0.65 (nu/MHz)^-2.1 (T/10^4 K)^-1.35 (n_e/0.03 cm^-3)^2. For nu = 0.1 MHz this is about 0.074 pc^-1 at the nominal local density, so the emission region is the first ~14 pc. The YMW16 model has a 50% uncertainty in n_e, and the local temperature may be as low as ~7000 K in the Local Interstellar Cloud; together these change tau at the emission horizon by a factor of several. Because the signal is the line-of-sight integral of emissivity times exp(-tau), a factor of 2 in n_e (factor of 4 in tau) changes the signal by exp(-Delta tau) ~ 50 for tau ~ 1, shifting the epsilon limit by roughly a factor of 7, not the quoted O(1). The cancellation with the observed synchrotron background is only partial, since at a fixed frequency the dark-Compton signal is produced by few-GeV electrons while the synchrotron at 0.1 MHz is from ~200 MeV electrons; these populations have different spatial distributions and energy-loss lengths. Without an uncertainty band on the limits in Fig. 4, the 'leading constraints' claim for mA' <~ 2e-17 eV is not fully established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a new observable for ultralight dark photon dark matter (DPDM): dark inverse Compton scattering of cosmic-ray electrons and positrons off the DPDM background produces a diffuse radio background, most prominently below 1 MHz. The authors derive the relevant cross section from first principles in Appendix A, compute the sky flux using the GALPROP cosmic-ray model SSZ4R20T150C2 and the YMW16 free-electron density map, and compare the predicted sky-averaged brightness with IMP-6, RAE-2, and Parker Solar Probe measurements. They derive upper limits on the kinetic mixing parameter, claiming leading constraints for masses mA' ≲ 2×10^-17 eV, with ε ≲ (2-10)×10^-6 depending on the dataset. The paper also discusses the dominant systematic from free-free absorption and mentions future sub-MHz radio missions.","tokens_in":13328,"tokens_out":8986,"duration_ms":96396,"significance":"If the limits survive a careful treatment of systematic uncertainties, this is a genuinely new and competitive probe of an experimentally difficult mass range. The central derivation is clean and self-contained, the input models are independent published products, and the limit procedure is a conservative inequality rather than a fit. The paper is also explicit about the linear scaling with the dark matter fraction and about the sub-dominant role of cosmic-ray protons. The main weakness is that the normalization of the predicted flux at the constraining frequencies is exponentially sensitive to the local free-electron density and temperature, and the final limits are presented without an uncertainty band; this is the key issue that needs to be addressed before the leading-constraint claim is fully established.","major_comments":[{"comment":"The normalization of the predicted flux in the constraining frequency range is controlled by the absorption factor P_nu = exp(-tau_nu) in Eq. (8). At 0.1-0.5 MHz the absorption horizon is only tens to hundreds of parsecs, so the limits are set almost entirely by the local electron density and temperature assumed in the YMW16 model. The text acknowledges a 50% uncertainty in the electron density and quotes an O(1) variation in the limits, but because tau_nu is proportional to n_e^2 T^-1.35 and the signal is an integral of emissivity times exp(-tau_nu), a 50% density change or a plausible change in the local temperature can alter the effective integration depth by substantially more than a factor of two. The exclusion curves in Fig. 4 are drawn without any band, so the claim that these are leading constraints for mA' <~ 2×10^-17 eV is not yet fully supported. Please propagate the plausible range of local n_e and T into the limits and show them as a band in Fig. 4, or otherwise quantify how the plotted curves map onto the assumed absorption parameters.","section":"III, Eq. (8), Fig. 4"},{"comment":"The treatment of the observational data is not fully specified. For IMP-6 the paper states that the 'probable maximum spectrum observed' is used, but it does not define the digitization procedure, the frequency binning, or how the upper envelope is constructed. For RAE-2 and PSP it is not stated whether the published measurement uncertainties are propagated into the limits. Since the constraints are driven by a small number of frequency bins near 0.1-0.5 MHz, the exact choice of sky-averaging and upper-envelope treatment can affect the quoted epsilon values. Please describe the data-reduction procedure in enough detail to reproduce it, and show the measurement uncertainties in Fig. 3 if they are available.","section":"III, Figs. 3 and 4"}],"minor_comments":[{"comment":"The notation for the integrated cosmic-ray flux is inconsistent: Eq. (C1) defines Phi, while the text immediately after uses the barred symbol \\bar{Phi}; please make the notation uniform.","section":"Appendix C, Eqs. (C1)-(C2)"},{"comment":"The legend entries 'SuperMAG (1s)' and 'SuperMAG (1m)' are not defined in the caption or the text; please state what these two curves represent.","section":"Fig. 4 caption and legend"},{"comment":"After Eq. (A8) the text refers to the spin-averaged matrix element square but writes \\overline{|M|^2} = (1/6) \\sum |M|^2; please specify explicitly which initial and final spin states are being summed and averaged.","section":"Appendix A, Eq. (A8)"},{"comment":"The text says the electron and positron fluxes are consistent with diffuse synchrotron emission [42]; since this consistency check involves assumptions about the magnetic field, it would be helpful to state the level of agreement or the reference for that comparison.","section":"III, paragraph on cosmic-ray models"},{"comment":"The phrase 'decameter wavelength regime' is imprecise for the sub-MHz band discussed in the paper; consider using 'hectometer' or a quantitative frequency range to avoid confusion with the decameter band at ~3-30 MHz.","section":"IV, outlook"}],"recommendation":"major_revision","confidential_remarks":"This is a well-executed and interesting paper: the cross-section derivation is first-principles, the inputs are independent, and the limit procedure is conservative. The main risk is not circularity or internal inconsistency but external systematic uncertainty in the free-free absorption model, which sets the normalization of the predicted sub-MHz flux. The stress-test concern about the exponential sensitivity of the limits to local n_e and T lands on a real gap: the manuscript acknowledges the uncertainty but does not propagate it into a band on the final exclusion curves. A revision that adds such a band and documents the data treatment for IMP-6, RAE-2, and PSP would make the leading-constraint claim solid. The paper is within scope for the journal and does not require rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a genuinely new probe. The idea is that cosmic-ray electrons and positrons scatter off the background of ultralight dark photons and up-scatter them to radio frequencies. The cross-section in Appendix A is derived from first principles, the flux integral is standard, and the comparison with IMP-6, RAE-2, and PSP data is handled conservatively—they demand the predicted brightness not exceed the observed spectrum, rather than fitting backgrounds. For m_A' below about 1e-17 eV this appears to give the strongest current bound on the kinetic mixing, epsilon <~ few x 1e-6.\n\nThe paper does its homework: the GALPROP cosmic-ray model is consistent with AMS-02 and Fermi, the YMW16 electron density map is used for absorption, and they explicitly check that cosmic-ray model uncertainty is at the percent level within the first ~100 pc that matters at these frequencies. They also verify that stopping power on DPDM is negligible and that proton contributions are negligible. The limiting procedure is a conservative inequality, not a fit, so there is no circularity issue.\n\nSoft spots. The main one is the free-free absorption model. At 0.1-0.5 MHz the optical depth is so large that only the first few tens of parsecs contribute. The YMW16 local density has a stated 50% uncertainty, and the warm ISM temperature is an assumption. The stress-test worry about exponential amplification sounds scary, but once you integrate over the line of sight the flux scales roughly as 1/n_e^2, so a 50% density error moves the epsilon limit by a factor of about 1.5-2, consistent with the paper's own O(1) estimate. Still, they should show that as a band on the exclusion curves; they already do the calculation, and the absence of a band undercuts the 'leading constraints' claim. A second, more minor issue is reproducibility: they don't release the digitized 1970s spectra or the code, so an independent check requires some archaeology.\n\nWho is this for? Dark photon DM aficionados and anyone interested in low-frequency radio astronomy as a DM probe. It's a solid, honest first calculation of a new channel. It deserves a serious referee: the cross-section derivation should be checked, the absorption treatment scrutinized, and the data handling verified, but this is exactly the kind of paper that should go through review rather than be desk-rejected.","headline":"New and honestly derived sub-MHz radio probe of ultralight dark photons; the constraints likely survive, but the missing absorption uncertainty band should be added before this is accepted as 'leading.'","tokens_in":13903,"tokens_out":3453,"would_cite":true,"duration_ms":34122,"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":"Ultralight dark photons would inverse-Compton scatter off cosmic-ray electrons into a diffuse sub-MHz radio glow, and existing sky-brightness data already constrain the kinetic mixing parameter to about 2e-6 to 1e-5 for masses below about…","keywords":["ultralight dark photon dark matter","dark inverse Compton scattering","sub-MHz radio background","kinetic mixing parameter","cosmic-ray electrons and positrons","Galactic radio spectrum","radio absorption","interstellar medium"],"falsifier":"Measure the free-electron density within the first ~100 pc along several high-latitude sightlines using pulsar dispersion measures toward the nearest pulsars, recompute $\\tau_\\nu$ at 0.1 to 0.5 MHz from Eq. (8), and re-evaluate the predicted sky brightness; if the flux at 0.2 MHz changes by more than a factor of two relative to the YMW16-based estimate, the exclusion curves shift by the same factor. Alternatively, a balloon or lunar sub-MHz telescope that maps the 0.1 to 1 MHz sky and finds no frequency-dependent component with the morphology shown in Figure 2 would directly test whether this scattering channel produces the claimed background.","tokens_in":12818,"feed_emoji":"📡","tokens_out":14424,"duration_ms":121899,"temperature":0.7,"pith_summary":"Ultralight dark photons—a well-motivated dark-matter candidate—are usually hard to detect because laboratory shields and haloscope setups lose sensitivity at tiny masses. The paper proposes a new route: cosmic-ray electrons and positrons scattering off the dark-photon field will up-scatter the dark photons into ordinary radio photons, producing a diffuse Galactic glow that peaks below 1 MHz. The glow's brightness grows linearly with the dark-matter density and quadratically with the kinetic-mixing parameter $\\epsilon$, and its typical frequency is set by the dark-photon mass times the square of the cosmic-ray energy. Comparing the predicted sky brightness with archived sub-MHz measurements from IMP-6, RAE-2, and the Parker Solar Probe, the paper derives the strongest existing limits on $\\epsilon$ for masses below about $2\\times10^{-17}$ eV, reaching $\\epsilon \\sim 2\\times10^{-6}$ with IMP-6 and about $1\\times10^{-5}$ with the other two datasets. The main caveat, stated in Section III, is that the sub-MHz absorption depends on the Local Bubble's electron density through the YMW16 model—quoted at up to 50% uncertainty with an order-one effect on the limits—and the final exclusion curves are drawn without an uncertainty band.","feed_headline":"Dark photons would glow below 1 MHz, and old maps already bound them","feed_subtitle":"Three space missions' sky maps constrain the mixing parameter to a few parts per million for the lightest dark photons.","key_machinery":"The load-bearing object is the dark inverse Compton cross-section, Eq. (4), $\\frac{d\\sigma_{\\rm DC}}{dE_\\gamma} \\simeq \\frac{\\epsilon^2 e^4}{12\\pi m_{A'} p_e^2}\\left[1 - 2x + 2x^2\\right]$ with $x = E_\\gamma/(2E_\\gamma^{\\max})$, together with the kinematic relation $E_\\gamma^{\\max} = 2\\gamma_e^2 m_{A'}$ that maps a cosmic-ray energy to an emitted photon frequency. This relation gives the paper its characteristic prediction, $\\nu \\simeq 0.93\\,{\\rm MHz}\\, (E_e/10\\,{\\rm GeV})^2 (m_{A'}/10^{-17}\\,{\\rm eV})$, so each dark-photon mass is probed by a specific cosmic-ray energy band. The flux integral in Eq. (6) sums electron, positron, and secondary-electron spectra over the line of sight and multiplies by the absorption factor $P_\\nu = \\exp(-\\tau_\\nu)$, with $\\tau_\\nu \\simeq 0.65 (\\nu/{\\rm MHz})^{-2.1} (T/10^4\\,{\\rm K})^{-1.35} \\int n_e^2\\, ds$. Because $\\tau_\\nu$ scales so steeply with frequency, the 0.1 to 0.5 MHz measurements that drive the constraints see only the first ~100 pc of the Galaxy, making the local electron density the controlling input and reducing cosmic-ray and dark-matter-profile uncertainties to the percent level.","core_discovery":"The paper establishes that the dark inverse Compton process $e^- + A' \\to e^- + \\gamma$ (the same for positrons) converts a small but observable fraction of cosmic-ray electron and positron energy into a diffuse, almost isotropic background of sub-MHz radio photons. For a dark-photon mass around $10^{-17}$ eV and a 10 GeV cosmic-ray electron, the average emitted photon lands near 0.93 MHz, so the signal sits below the frequencies normally used in radio astronomy. When the production rate is integrated through the Milky Way with absorption modeled by the YMW16 electron-density map, the resulting sky brightness is large enough that the measured 0.1 to 1 MHz spectra from IMP-6, RAE-2, and the Parker Solar Probe already exclude kinetic mixing down to $\\epsilon \\sim (2\\text{--}10)\\times10^{-6}$ for dark-photon masses below about $2\\times10^{-17}$ eV. The limits are conservative in that no astrophysical radio background is subtracted, and they scale linearly with the dark-photon abundance if dark photons are only a sub-component of the dark matter.","pith_inferences":["Because the signal scales linearly with the dark-photon fraction, the same sky-brightness data can be recast as upper limits on the fractional abundance of dark photons at a fixed mixing; for $m_{A'}$ near $10^{-18}$ eV and $\\epsilon$ at the current exclusion, the allowed fraction would be well below unity, a quantitative bound the paper leaves implicit.","Beyond what the paper computes, the 0.1 MHz glow's near-isotropy and the 1 MHz glow's sensitivity to local electron-density structure suggest a future sub-MHz survey could fit $\\epsilon$ and the local free-electron density jointly, breaking the absorption degeneracy and strengthening the limits.","The same flux pipeline, applied to axion-like particles and dilatons with photon and electron couplings, would predict background shapes that differ in frequency and anisotropy; matching those shapes could eventually distinguish which ultralight boson is responsible for any detected sub-MHz excess."],"forward_implications":["For dark-photon masses below about $2\\times10^{-17}$ eV, sub-MHz radio sky-brightness data now provide the strongest constraint on the kinetic mixing parameter, reaching $\\epsilon \\sim 2\\times10^{-6}$ with the IMP-6 dataset and about $1\\times10^{-5}$ with RAE-2 and Parker Solar Probe data.","If dark photons are only a fraction of the dark matter, the excluded $\\epsilon$ values scale linearly with that fraction, so the same observations continue to produce useful limits in sub-component scenarios.","The predicted glow is almost isotropic at 0.1 MHz and becomes anisotropic near 1 MHz, where the electron-density structure of the Milky Way is resolved, so frequency-dependent anisotropy is a handle for future searches.","The same cosmic-ray up-scattering mechanism should generate a comparable sub-MHz background for other ultralight bosons, such as axion-like particles and dilatons with photon and electron couplings, giving future sub-MHz space missions a generic light-dark-matter signature to look for."],"supporting_citations":[{"why":"introduced inverse Compton-like scattering of ultrahigh-energy cosmic-ray protons against dark photons, the approach this paper extends to cosmic-ray electrons and positrons.","marker":"[24]"},{"why":"provides the Explorer 43 (IMP-6) Galactic radio spectrum between 130 and 2600 kHz that sets the tightest of the paper's constraints.","marker":"[25]"},{"why":"provides the Radio Astronomy Explorer 2 spectrum below 10 MHz used as a second sub-MHz sky-brightness dataset.","marker":"[26]"},{"why":"provides the Parker Solar Probe sky-brightness measurements from 0.5 to 7 MHz used as the third dataset.","marker":"[27]"},{"why":"supplies the free-free absorption optical-depth formula $P_\\nu = \\exp(-\\tau_\\nu)$ used to attenuate the predicted sub-MHz flux.","marker":"[33]"},{"why":"supplies the YMW16 Galactic free-electron density model that sets the optical depth and, through the Local Bubble, the overall normalization of the signal.","marker":"[34]"},{"why":"supplies the cosmic-ray electron, positron, and secondary-electron spectra (model SSZ4R20T150C2) used as the scattering target flux.","marker":"[35]"},{"why":"provides the numerical cosmic-ray propagation machinery used to generate the adopted particle spectra.","marker":"[36]"}],"fun_headline_variants":["Sub-MHz radio glow from dark photons pinned by old space data","Dark photon scattering off cosmic rays yields sub-MHz signal","Old maps constrain ultralight dark photons via sub-MHz glow","Ultralight dark photons revealed by sub-MHz cosmic-ray glow","Sub-MHz background from dark photons bound by three probes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the YMW16 local free-electron density and an assumed interstellar temperature of $10^4\\,$K set the sub-MHz absorption; if the gas within the first ~100 pc is denser than modeled, the predicted glow is absorbed more, the predicted flux drops, and the quoted $\\epsilon$ bounds weaken, while a less dense medium would strengthen them.","fun_headline_variants_meta":{"raw":{"variants":["Sub-MHz radio glow from dark photons pinned by old space data","Dark photon scattering off cosmic rays yields sub-MHz signal","Old maps constrain ultralight dark photons via sub-MHz glow","Ultralight dark photons revealed by sub-MHz cosmic-ray glow","Sub-MHz background from dark photons bound by three probes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000574,"raw_usage":{"total_tokens":2694,"prompt_tokens":910,"completion_tokens":1784,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":1697}},"tokens_in":526,"tokens_out":1784,"duration_ms":12282,"temperature":1.0,"reasoning_tokens":1697,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:28:23.876914+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the free-electron density within the first ~100 pc along several high-latitude sightlines using pulsar dispersion measures toward the nearest pulsars, recompute $\\tau_\\nu$ at 0.1 to 0.5 MHz from Eq. (8), and re-evaluate the predicted sky brightness; if the flux at 0.2 MHz changes by more than a factor of two relative to the YMW16-based estimate, the exclusion curves shift by the same factor. Alternatively, a balloon or lunar sub-MHz telescope that maps the 0.1 to 1 MHz sky and finds no frequency-dependent component with the morphology shown in Figure 2 would directly test whether this scattering channel produces the claimed background.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the Explorer 43 (IMP-6) Galactic radio spectrum between 130 and 2600 kHz that sets the tightest of the paper's constraints."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the Radio Astronomy Explorer 2 spectrum below 10 MHz used as a second sub-MHz sky-brightness dataset."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the Parker Solar Probe sky-brightness measurements from 0.5 to 7 MHz used as the third dataset."},{"cited_title":"An Ultra-long Wavelength Sky Model with Absorption Effect","cited_arxiv_id":"2104.03170","evidence_quote":"supplies the free-free absorption optical-depth formula $P_\\nu = \\exp(-\\tau_\\nu)$ used to attenuate the predicted sub-MHz flux."}],"review_version":1}