{"id":"635bf73e-1d73-4758-b3bf-c8a0e339e598","arxiv_id":"2411.13701","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"Resonant conversion of CMB photons into axions in a primordial magnetic field creates a characteristic spectral distortion, and the updated calculations tighten axion-photon coupling constraints and reach previously unexplored parameter space.","lead":"This paper updates limits on axions, hypothetical particles that could steal tiny amounts of light from the cosmic microwave background and distort its frequency spectrum. The new calculations show that, if large-scale magnetic fields exist near current upper limits, spectral distortion measurements can probe axion parameter space that other experiments cannot reach.","discovery_kind":"extension","skeptic_critique":null,"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits constraints on resonant CMB photon-to-axion conversions from CMB spectral distortions. It improves the modeling of the photon plasma mass by adding neutral hydrogen and helium contributions to the free-electron term, derives an analytic spectral template for the axion distortion, and uses the CosmoTherm code to evolve the distortion through the thermalization epochs. The resulting limits are compared with COBE/FIRAS data and with a PIXIE-like forecast, and are mapped into the gaγγ–ma plane for several assumed primordial magnetic field amplitudes. The paper also explores the large-distortion regime, where a significant fraction of CMB photons are converted.","tokens_in":24826,"tokens_out":18149,"duration_ms":197622,"significance":"If the central derivation holds, the paper provides a useful update to spectral-distortion constraints on axion-photon couplings, with quantitative improvements from the full spectral template rather than energy-density-only limits, and a physically better motivated plasma-mass treatment. The analytic framework in Eqs. (17)–(21) is transparent and cross-checked against CosmoTherm, the broad agreement with Mirizzi et al. is a valuable consistency check, and the PIXIE forecasts give a concrete, falsifiable target. The authors are also explicit that the constraints are conditional on the existence of a large-scale primordial magnetic field near the current Planck upper bound, and they show how the reach degrades with B0. The main caveats are associated with the normalization convention for the magnetic-field fudge factor, the self-consistency of the low-mass constraints with the recombination history, and the completeness of the large-distortion treatment.","major_comments":[{"comment":"The angle-averaging argument in footnote 12 gives ⟨B_T²⟩ = (1/3) B0², so the rms transverse field entering Eq. (3b) linearly is B_T = B0/√3, not f = 1/3 as stated in the text. Since κ is linear in B_T and the conversion signal scales as κ², using f = 1/3 in Fig. 13 normalizes the gaγγ constraints differently from the stated convention by a factor √3 (the limits are weaker by √3 for fixed B0). Please correct the definition of f or recompute Fig. 13, and state the convention unambiguously.","section":"Sec. V B and Eq. (3b)"},{"comment":"The low-mass constraints (ma ≲ few × 10⁻¹⁰ eV) are computed with the standard CosmoRec recombination history and assume that the conversion does not alter that history. At the FIRAS 95% limit in this mass range, Eq. (21) gives |Δργ/ργ| ≈ 6 × 10⁻⁵, corresponding to |ϵρ| ≈ 10⁻³ at the conversion redshift; this is not obviously negligible for the ionization balance and would feed back into Xe and hence into mγ and the resonance locations. Please quantify this feedback or restrict the mass range over which the constraints are quoted.","section":"Sec. II B and Sec. V (Figs. 1 and 11)"},{"comment":"The treatment of the large-distortion regime is presented as complete (\"we treat for the first time the large-distortion regime\"), but the CosmoTherm integrations are stopped at a scattering y-parameter of about 0.3 because a photon shock renders the numerical treatment insufficient (Sec. IV A), and the Landau-Zener expression is used up to γcon ≃ 1 despite caveats from Refs. [48–50]. The conclusion that no sign flip occurs and the interpretation of Fig. 12 therefore rest on an incomplete numerical evolution. Please either complete the evolution with a method that handles the shock or soften the claims to describe an exploratory treatment of the onset of the large-distortion regime.","section":"Sec. IV A and Sec. V A (Figs. 4–6, 12)"}],"minor_comments":[{"comment":"The sentence \"This assumption is not crucial when considering considering axion masses\" contains a duplicated word; please edit.","section":"Sec. II A, after Eq. (3)"},{"comment":"The index set \"i = e, HI, HeI, HeII\" is not consistent with the displayed term Xe having no κi; please clarify that κe ≡ 1 and define Xp explicitly.","section":"Sec. II B, Eq. (15)"},{"comment":"The ordinate label \"f ga\" should be \"f gaγγ\" and the f convention should be stated in the caption, especially in light of the normalization issue raised above.","section":"Fig. 13"},{"comment":"The multiple-conversion implementation replaces δ(z − zcon,i) with a Gaussian of width Δz/z ≈ 10⁻²; a sentence justifying that this width is small compared with the thermalization timescale would improve reproducibility.","section":"Sec. IV C"},{"comment":"The paper would benefit from a data/code availability statement, since the FIRAS likelihood is described only by reference to Ref. [45].","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper builds substantially on the authors' own dark-photon paper (Ref. [45]), including the FIRAS likelihood. A referee with access to that paper should verify that the reanalysis is not circular and that the reported ~30% tightening from the full spectral shape is robust. The f-normalization issue in Fig. 13 should be fixed before publication; it does not change the qualitative conclusions but affects the quantitative limits. The low-mass back-reaction concern is, in my view, the most scientifically important point to address, as it bears directly on the claimed reach to ma ≲ 10⁻¹⁰ eV."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one. It's a careful, honest reanalysis of resonant axion-photon conversion limits from CMB spectral distortions, and it holds up better than most papers in this space. The genuinely new bits are the corrected atomic-hydrogen refractive index coefficient (the old literature value comes from molecular hydrogen, H2, not HI; they show the coefficient drops from ~7.7e-3 to ~5.0e-3 eV^-2), the inclusion of helium in the photon plasma mass, and the explicit axion distortion template A(x). They validate the analytic small-distortion framework with CosmoTherm and find the full-shape analysis tightens the COBE/FIRAS and PIXIE constraints by up to ~30% relative to the energy-density-only estimate. The broad agreement with Mirizzi, Redondo, and Sigl (2009) is a good consistency check. The central constraints are built on clean logic: Landau-Zener conversion probability, no fitted normalization, external COBE/FIRAS and Planck data. I did not find a load-bearing error.\n\nThe soft spots are real but not fatal. The 'previously unexplored regions' claim in the abstract and Fig. 13 is conditional on a ~1 nG comoving primordial magnetic field at ~1 Mpc coherence. That's a Planck upper bound, not a detection, and the authors say the limits degrade roughly linearly with B0, becoming suboptimal below ~1e-2 nG. They flag it, but the visibility of that caveat matters. The large-distortion regime is only partially treated: the numerical run develops a photon shock at low frequencies and they stop at scattering y~0.3, and they note Landau-Zener corrections may matter at gamma_con ~1. That's fine because the Neff bound keeps them out of that regime, but the 'first treatment' is more of a first look. No code or data release accompanies the numerical constraint curves, and the FIRAS likelihood itself is deferred to the companion dark-photon paper. Minor: the paper leans on that companion, Ref [45], for machinery, which is legitimate but makes this paper less self-contained.\n\nBottom line: this is a solid refinement of an existing constraint program, with one genuinely useful correction (the HI vs H2 coefficient) and a credible numerical upgrade. Send it to a serious referee. I'd bring it to reading group and would cite it for the plasma mass treatment.","headline":"A careful, honest reanalysis of axion-photon conversion limits from CMB spectral distortions; the corrected atomic-hydrogen plasma mass is the real gem and the constraints hold up, so it deserves serious review.","tokens_in":25368,"tokens_out":2099,"would_cite":true,"duration_ms":63896,"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":"Resonant CMB photon-to-axion conversion in a ~1 nG primordial magnetic field would imprint a characteristic spectral distortion, and COBE/FIRAS residuals already exclude axion-photon couplings that other experiments have not reached.","keywords":["axions","axion-like particles","CMB spectral distortions","resonant photon-axion conversion","primordial magnetic fields","photon plasma mass","COBE/FIRAS","PIXIE"],"falsifier":"A measurement that places an upper limit $B_0^{\\rm rms} \\lesssim 10^{-2}$ nG on a $\\sim$1 Mpc-scale cosmological magnetic field would falsify the core claim, since the distortion scales roughly as $(g_{a\\gamma\\gamma} B_0)^2$ and would fall below the COBE/FIRAS threshold. A second decisive test would be a direct axion detection at a coupling and mass inside the region the paper's Fig. 13 rules out under the assumed 1 nG field.","tokens_in":24684,"feed_emoji":"🌌","tokens_out":10958,"duration_ms":108167,"temperature":0.7,"pith_summary":"This paper argues that the near-perfect blackbody spectrum measured by COBE/FIRAS can act as a laboratory for axions and axion-like particles: when the plasma mass of a CMB photon matches the axion mass during cosmic expansion, the photon can resonantly convert into an axion, leaving a frequency-dependent deficit in the CMB spectrum. The authors improve the photon plasma mass by including neutral hydrogen and helium, derive the spectral template of the resulting distortion, and evolve it through the mu- and y-eras with CosmoTherm. For a large-scale primordial magnetic field near the current CMB upper bound, about 1 nG, they find that COBE/FIRAS residuals already exclude axion-photon couplings in mass regions that other experiments have not reached, and a PIXIE-like experiment would extend the reach by about two orders of magnitude. They also treat the large-distortion regime for the first time, finding that stimulated Compton scattering creates a photon excess at low frequencies that prevents the naively expected sign flip of the distortion. A reader should care because CMB spectral distortions are one of the few probes sensitive to axion masses around $10^{-13}$ to $10^{-4}$ eV, a window that other searches cover poorly if at all.","feed_headline":"COBE/FIRAS data reaches axion masses other probes miss","feed_subtitle":"A ~1 nG primordial magnetic field would make CMB spectral distortions an axion probe; COBE/FIRAS already sets limits.","key_machinery":"The load-bearing object is the resonant conversion condition $m_a \\simeq m_\\gamma(z,\\omega)$, where the effective photon mass is computed from the refractive index of the pre- and post-recombination plasma. The paper uses $m_\\gamma^2 = \\tilde{\\omega}_p^2 [X_e - (\\omega/\\mathrm{eV})^2 \\sum_i \\kappa_i X_i]$ with $\\kappa_{\\mathrm{HI}} \\simeq 5.0\\times 10^{-3}\\,\\mathrm{eV}^{-2}$, $\\kappa_{\\mathrm{HeI}} \\simeq 1.8\\times 10^{-3}\\,\\mathrm{eV}^{-2}$, and $\\kappa_{\\mathrm{HeII}} \\simeq 3.1\\times 10^{-4}\\,\\mathrm{eV}^{-2}$, so that high-frequency photons can have a negative effective mass-squared and stop converting. Conversion probabilities are handled with the Landau-Zener form $P = 1 - \\exp(-\\gamma_{\\mathrm{con}} x)$, where $\\gamma_{\\mathrm{con}} \\propto (g_{a\\gamma\\gamma}B_0)^2$ divided by the derivative of $\\ln m_\\gamma^2$ at the resonance; this expression determines which photon frequencies convert at which redshifts. The paper also defines the axion distortion template $A(x) = (G_3/(3G_2))G(x) - x n_{bb}(x)$, the photon-number-conserving spectral shape left by the conversion before Comptonization drives it toward a $\\mu$-type distortion.","core_discovery":"The central claim is that resonant CMB photon-to-axion conversion, computed with an updated plasma mass and a full spectral-distortion treatment, produces competitive new constraints on the axion-photon coupling $g_{a\\gamma\\gamma}$. For axion masses roughly $10^{-13}$ to $10^{-4}$ eV and a comoving magnetic field of order 1 nG, the paper finds that COBE/FIRAS data explore regions of parameter space not yet accessed by other experimental probes, on both the high- and low-mass ends, with a PIXIE-type experiment improving the reach by about two orders of magnitude. This improvement is attributed to three changes: the refractive index of neutral atomic hydrogen is used instead of a molecular-hydrogen value, helium contributions are included, and the full shape of the distortion rather than a total-energy proxy is compared to data. For multiple conversions at low masses, the analytic conversion probability can change sign near $10^{-12}$ to $10^{-11}$ eV, but the CosmoTherm solutions keep the final distortion negative because photon condensation at low frequencies supplies the extra photons. The resulting constraints and forecasts appear in the paper's Figs. 11 and 13.","pith_inferences":["Editorial inference: the corrected atomic-hydrogen and helium refractive indices should also change predictions for axion dark matter converting into photons; the helium-induced conversion spikes in the Wien tail would be a source of signal rather than a negligible correction in the $\\Omega_a \\simeq \\Omega_{\\rm dm}$ scenario.","Editorial inference: the new template $A(x)$ is a generic signature of high-frequency photon removal, so future spectral-distortion fits could include it as an independent component alongside $\\mu$ and $y$; because its null lies near the $\\mu$-distortion null, multi-band data will be needed to separate them.","Editorial inference: the limits in Fig. 13 are conditioned on the assumed $B_0$; a null PIXIE measurement would tighten $g_{a\\gamma\\gamma}$ only for a fixed field, so the analytic framework could be used to present joint exclusions on $(m_a, g_{a\\gamma\\gamma}, B_0)$ rather than separate curves."],"forward_implications":["COBE/FIRAS residuals already exclude some $g_{a\\gamma\\gamma}$ values near $m_a \\sim 10^{-13}$ to $10^{-4}$ eV that other laboratory and astrophysical searches do not reach, if a ~1 nG large-scale magnetic field exists.","A PIXIE-like experiment would improve the $g_{a\\gamma\\gamma}$ reach by about two orders of magnitude, allowing a present-day field strength near $10^{-4}$ nG to produce a detectable distortion.","The analytic distortion template and energy estimate reproduce the CosmoTherm constraints well, and the full numerical shape tightens the PIXIE forecast by roughly 30%.","Multiple conversions for $m_a \\lesssim 10^{-10}$ eV make the distortion sign and shape sensitive to the details of recombination and reionization, so spectral-distortion data becomes a probe of those epochs as well.","In the large-distortion regime, strong axion conversions do not produce the sign-flipped positive distortion one might expect; stimulated Compton scattering and low-frequency photon destruction keep the net distortion negative, weakening constraints relative to naive estimates."],"supporting_citations":[{"why":"Supplies the COBE/FIRAS residual limits that set the distortion constraints.","marker":"[3]"},{"why":"The CosmoTherm code used to evolve the initial conversion spectrum through the Comptonization and thermalization eras.","marker":"[9]"},{"why":"Earlier analytic resonant-conversion constraint that this work updates and cross-checks.","marker":"[40]"},{"why":"Prior axion spectral-distortion treatment with multiple conversions; the paper compares its Boltzmann treatment with this probability formula.","marker":"[42]"},{"why":"Dark-photon conversion machinery, likelihood code, and Green's function formalism that the axion calculation builds on.","marker":"[45]"},{"why":"Planck-based upper bound on present-day large-scale magnetic field strength used as the benchmark $B_0^{\\rm rms} \\sim 1$ nG.","marker":"[51]"},{"why":"Updated Planck-related primordial-magnetic-field constraint used together with [51] to set the field-strength benchmark.","marker":"[52]"},{"why":"CosmoRec recombination history supplies the hydrogen and helium ionization fractions used in the plasma mass.","marker":"[64]"},{"why":"Catalog of existing axion limits used for the comparison in Fig. 13.","marker":"[35]"}],"fun_headline_variants":["COBE/FIRAS probes new axion masses via CMB distortions","CMB spectral distortions extend axion reach at ~1 nG field","Axion-photon conversion in CMB sets stringent limits","CMB distortion analysis opens fresh axion parameter space"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a large-scale primordial magnetic field with present-day amplitude near 1 nG, coherence length near 1 Mpc, and redshift scaling $B(z)=B_0(1+z)^2$ exists; if the field is absent or weaker than about $10^{-2}$ nG today, the claimed unexplored parameter-space regions from COBE/FIRAS disappear.","fun_headline_variants_meta":{"raw":{"variants":["COBE/FIRAS probes new axion masses via CMB distortions","CMB spectral distortions extend axion reach at ~1 nG field","Axion-photon conversion in CMB sets stringent limits","CMB distortion analysis opens fresh axion parameter space"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000236,"raw_usage":{"total_tokens":1545,"prompt_tokens":1029,"completion_tokens":516,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":645,"completion_tokens_details":{"reasoning_tokens":444}},"tokens_in":645,"tokens_out":516,"duration_ms":5672,"temperature":1.0,"reasoning_tokens":444,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:00:03.789522+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement that places an upper limit $B_0^{\\rm rms} \\lesssim 10^{-2}$ nG on a $\\sim$1 Mpc-scale cosmological magnetic field would falsify the core claim, since the distortion scales roughly as $(g_{a\\gamma\\gamma} B_0)^2$ and would fall below the COBE/FIRAS threshold. A second decisive test would be a direct axion detection at a coupling and mass inside the region the paper's Fig. 13 rules out under the assumed 1 nG field.","supporting_citations":[{"cited_title":"Hecht, Optics (Pearson, 2017)","cited_arxiv_id":null,"evidence_quote":"CosmoRec recombination history supplies the hydrogen and helium ionization fractions used in the plasma mass."}],"review_version":1}