{"id":"4b0b195a-b12e-477c-869a-2ba8e6d7e45a","arxiv_id":"2608.07803","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Accurate axion-isocurvature evolution across the full ultralight mass range converts Planck and BICEP bounds into inflation-scale limits and reveals a small discovery window near m_a ~ 10^-25 eV.","lead":"This paper extends the AxiECAMB cosmological code to compute how quantum fluctuations of an ultralight axion field during inflation appear as isocurvature patterns in the cosmic microwave background, across axion masses from dark energy to dark matter scales.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline bound Eq (22) rests on an unresolved factor-of-4 normalization discrepancy with Ref [12] (footnote 5), and the authors do not rerun a Planck MCMC with AxiECAMB to confirm it; this directly affects the mass threshold where isocurvature beats the tensor bound.","rationale":"The paper is careful, honest, and technically substantial: the ETA extension is physically motivated, the accuracy tests in Appendix A2 are extensive, and the analytic derivation of Eq (22) from Eqs (4), (20), and the Planck beta_iso bound is internally consistent. I verified the coefficient 2.41 in Eq (21) and the resulting 0.076 in Eq (22) by direct calculation from F_1/3 and a_osc. The reader's weakest_assumption (frozen-field misalignment) is correctly stated as a condition, and the ETA failure for k >> k_J is openly disclosed; those are limitations, not internal contradictions. The most load-bearing unresolved issue is the factor-of-4 normalization discrepancy with Ref [12], because it directly concerns the numerical value of the headline bound that defines the crossover with tensor bounds. The authors flag it but do not resolve it, and since the old result came from a full MCMC with overlapping authorship, an independent rerun with AxiECAMB is the decisive check. Until that check is done, CONDITIONAL is the right verdict; my read does not move it to ACCEPT or REJECT, so the verdict is UNCHANGED.","tokens_in":25493,"tokens_out":23444,"duration_ms":187240,"concrete_test":"Rerun the Ref [12] Planck likelihood (TT,TE,EE+lowE+lensing) with AxiECAMB v1.1 (movH_switch=20, accuracy_boost=2) for m_a=10^-24 eV, sampling f_dm and r_a, and compare the 95% CL upper limit on r_a f_dm to 0.076 (m_a/10^-27 eV)^-1/2 = 0.0024 at f_dm=1. If the MCMC limit is approximately 0.0024, Eq (22) is confirmed and the Ref [12] result was a code artifact; if it stays near 0.01, the analytic mapping misses a degeneracy or normalization effect and the headline constraint must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Central quantitative claim is Eq (22) from Planck's beta_iso < 0.038, giving r_a f_dm < 0.076 (m_a/10^-27 eV)^-1/2, which beats BICEP for m_a f_dm^2 >~10^-26.4 eV. The normalization is the unresolved factor-of-4 discrepancy with Ref [12] flagged in footnote 5: at m_a=10^-24 eV, f_dm=1, this paper gives r_a < 0.0024 while Ref [12] (overlapping authorship) reported r_a ~< 10^-2 from a full Planck MCMC. The footnote states the old analytic expression was 'improperly normalized by a large factor' and that the old MCMC did not reproduce even the scaling, but the authors 'do not explore this issue further.' Because Eq (22) is derived from the analytic abundance relation rather than a new MCMC, the old MCMC may have captured a degeneracy or transfer-function effect omitted here; this shifts the threshold Eq (23) by up to ~16 in m_a and moves the claimed coexistence window. Separately, the paper's own Appendix A1 limits the 'accurate across the full mass range' claim: for isocurvature k >> k_J, Eq (A5) fails and AxiECAMB reverts to instantaneous values with only order-of-magnitude accuracy, affecting the intermediate-mass signatures though not Eq (22).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript extends the effective time average (ETA) treatment in AxiECAMB v1.1 to inflationary axion isocurvature perturbations, gives a fitting formula for the present axion abundance as a function of the initial field value φ_ini across masses from the dark-energy regime to the dark-matter regime, and uses that mapping to translate the Planck bound on CDM isocurvature, β_iso < 0.038, into the constraint r_a f_dm < 0.076 (m_a/10^-27 eV)^-1/2 (Eq. 22). The paper also derives analytic scalings for Jeans-suppressed intermediate masses and for the quadrupole-dominated dark-energy regime, presents Fisher forecasts, and identifies a coexistence window near m_a ~ 10^-25 eV, f_dm ≳ 0.1 where axion isocurvature and tensor modes could both be detected. The appendices document the code implementation, the criterion (Eq. A5) for when ETA is valid, and accuracy tests against lower-accuracy codes.","tokens_in":25716,"tokens_out":10204,"duration_ms":89225,"significance":"If the central constraint in Eq. (22) is correct, the paper establishes the CMB isocurvature channel as a competitive probe of the inflationary energy scale for ultralight axions, with a bound that is stronger than BICEP/Keck in a well-defined mass range and with falsifiable predictions in three distinct mass regimes. The submission has notable strengths: the code is public, the analytic fitting forms are explicitly calibrated and accompanied by accuracy tests in Appendix A, and the limitations of the ETA approximation are stated honestly, including order-of-magnitude accuracy in the k >> k_J regime. These features make the results reproducible. The significance is conditional on resolving the normalization discrepancy with an existing full Planck MCMC analysis and on qualifying the low-accuracy regime, both discussed below.","major_comments":[{"comment":"The headline bound is not yet settled. The paper reports r_a < 0.0024 f_dm^-1 at m_a = 10^-24 eV while the full Planck MCMC analysis of Ref. [12] with overlapping authorship found r_a ≲ 10^-2 at f_dm = 1, a factor of about four. The footnote attributes this to an improper normalization in Ref. [12] and states that the issue is not explored further. Because Eq. (22) is derived from the analytic abundance mapping and not from a new MCMC, this discrepancy directly affects the claimed threshold where isocurvature beats the tensor bound (Eq. 23) and the coexistence window in Fig. 5; a factor of four in r_a shifts the mass threshold by about a factor of sixteen. I request either a Planck MCMC rerun with AxiECAMB v1.1 or a detailed quantitative reconciliation showing that the old normalization error accounts for the full difference and that no degeneracy or transfer-function effect is missing.","section":"Footnote 5; Eq. (22)"},{"comment":"The abstract's claim that AxiECAMB v1.1 'accurately evolve[s] these perturbations across the full axion mass range' is stronger than the validation supports. Equation (A5) fails for k >> k_J, and the code then reverts to instantaneous field values with only order-of-magnitude accuracy; this is the Jeans-suppressed regime that produces the unique intermediate-mass and lightest dark-energy signatures in Figs. 6 and 7. The authors should either improve the treatment of this regime or explicitly restrict the accuracy claims to the masses and multipoles for which the predictions are better than order of magnitude, and adjust the abstract and Section III accordingly.","section":"Appendix A1; Sec. III B"},{"comment":"The calibrations C_P = 0.32 and C_Q = 0.0097 are fit to AxiECAMB at ℓ = 2 and then used as analytic predictions, while the abundance formula in Eqs. (8)-(10) is calibrated to and validated against the same code. This means the quoted sub-percent accuracy and the plateau/quadrupole estimates are internally consistent but not independent validations of the isocurvature transfer functions. I ask for at least one independent cross-check: e.g., the CDM isocurvature spectrum from CAMB for the m_a ≫ H_eq plateau, or a direct likelihood comparison with Ref. [12], together with a statement of how much the derived r_a bounds change under the calibration uncertainty.","section":"Eqs. (17), (30), Appendix A2"}],"minor_comments":[{"comment":"The double use of τ for conformal time and reionization optical depth is acknowledged in the text, but replacing one of the symbols would reduce confusion.","section":"Sec. II, around Eq. (12)"},{"comment":"The threshold 10^2 in Eq. (A5) is set without a derivation; a brief derivation or a convergence test justifying this value would strengthen the appendix.","section":"Appendix A1, Eq. (A5)"},{"comment":"The caption refers to 'points, blue solid curve' for the AxiECAMB result, but the figure appears to show only a solid curve; please clarify the legend.","section":"Fig. 1 caption"},{"comment":"The statement that the analytic estimates of Ref. [12] shift the relevant mass range by 'orders of magnitude' is hard to reconcile with the factor-of-four normalization discrepancy quoted in footnote 5; please state the resulting mass shift precisely.","section":"Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The key risk is the factor-of-four discrepancy with Ref. [12]. I would not accept the paper until the authors resolve it with a new MCMC or a complete quantitative explanation. The overlapping authorship makes this especially important, and the current footnote leaves the main quantitative claim inadequately defended. The paper is otherwise a solid technical contribution, and the honesty of Appendix A1 is a point in its favor."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is worth engaging. It does something genuinely new: extends the ETA approach to axion isocurvature perturbations, provides analytic abundance fits across a huge mass range, and derives a constraint mapping that shifts the relevant mass window by orders of magnitude. The technical core is plausible. The initial conditions look right, the switch criterion in Eq. (A5) is a sensible guard against the ETA failing for fast Jeans oscillations, and the accuracy tests in Appendix A2 are thorough, including comparisons against exact KG solutions and a documented order-of-magnitude regime. The code is public, and the paper is honest about its own limits. Footnote 5 flags the normalization discrepancy with Ref. [12], and Appendix A1 states where accuracy degrades.\n\nThe soft spots are real but localized. The unresolved factor-of-4 discrepancy with Ref. [12] is the main one. The authors attribute it to an improperly normalized analytic expression in that earlier work, but they do not rerun a Planck MCMC to confirm that their new mapping holds up against actual data. The old MCMC may have captured degeneracies or transfer-function effects that the analytic route misses. Since the headline bound in Eq. (22) comes from the analytic abundance relation plus the Planck beta_iso constraint, this unresolved factor directly shifts the mass threshold where isocurvature beats the tensor bound, by up to a factor of ~16 in m_a. That is not a cosmetic issue; it moves the claimed coexistence window. The paper should either run the MCMC or explain quantitatively why the old chain was wrong. The abstract's \"sub-percent\" accuracy claim also overstates the body, which reports 1-10% accuracy at the lighter end. That is minor, but worth fixing. The k >> k_J failure is a genuine limitation, though it affects the suppressed regimes rather than the main dark-matter bound.\n\nWho is this for? Anyone working on axion cosmology, CMB constraints on inflation, or ultralight dark matter. It deserves a serious referee. My recommendation: send it to peer review, with a requirement that the authors resolve the Ref. [12] normalization discrepancy before publication. The paper does the field a service by providing tools and scaling relations, but the headline number should rest on a verified foundation.","headline":"Careful, honest extension of AxiECAMB to isocurvature modes with a strong headline bound that is conditionally correct pending resolution of a factor-of-4 normalization discrepancy with prior work.","tokens_in":26426,"tokens_out":1790,"would_cite":true,"duration_ms":17566,"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":"Planck's cold dark matter isocurvature bound, translated through the axion abundance–initial-field mapping, implies $r f_{\\rm dm} < 0.076\\,(m_a/10^{-27}\\,{\\rm eV})^{-1/2}$ for dark-matter-like ultralight axions, a constraint stronger than…","keywords":["axion isocurvature","ultralight axions","CMB anisotropy","inflation","tensor-to-scalar ratio","misalignment","dark matter","dark energy"],"falsifier":"Compute the isocurvature CMB power spectrum for $m_a = 10^{-28}\\,{\\rm eV}$, $f_{\\rm dm}=1$ with a switch at $m_a/H_* = 50$ and AccuracyBoost=3; if the low-$\\ell$ plateau or quadrupole moves by more than the quoted percent-level accuracy, the ETA-based numerical claim is not converged. Alternatively, detect isocurvature at $m_a \\sim 10^{-30}\\,{\\rm eV}$ at the level predicted by $r=0.036$; the paper predicts that signal is unobservable, so such a detection refutes its conclusion.","tokens_in":25115,"feed_emoji":"🔭","tokens_out":11197,"duration_ms":93301,"temperature":0.7,"pith_summary":"The paper aims to turn axion isocurvature fluctuations into a quantitative probe of inflation across the entire ultralight axion mass range, from dark-energy-like ($m_a \\lesssim H_0$) to dark-matter-like ($m_a \\gg 10^{-28}\\,{\\rm eV}$) axions. It claims that, once the axion abundance is accurately calibrated to the initial field value, Planck's bound on cold dark matter isocurvature implies $r f_{\\rm dm} < 0.076\\,(m_a/10^{-27}\\,{\\rm eV})^{-1/2}$ in the dark-matter regime, which is stronger than the current BICEP tensor bound for $m_a f_{\\rm dm}^2 \\gtrsim 10^{-26.4}\\,{\\rm eV}$. In the intermediate and dark-energy regimes the isocurvature power is strongly suppressed by Jeans and horizon effects, so the paper concludes that any primary-CMB detection there would signal physics beyond the standard frozen-field misalignment scenario. The payoff is a small region near $m_a \\sim 10^{-25}\\,{\\rm eV}$, $f_{\\rm dm} \\gtrsim 0.1$, where both axion isocurvature and tensor modes could be discovered below current bounds while easing the $S_8$ tension.","feed_headline":"Axion isocurvature beats tensor B-modes as inflation probe","feed_subtitle":"Planck isocurvature data beat tensor limits for most axion masses; the lightest signals would point beyond inflation.","key_machinery":"The load-bearing mechanism is the effective time average (ETA) implemented in the Boltzmann code AxiECAMB, which replaces the exact Klein-Gordon evolution of a light axion with a fluid description after a switch epoch $m_a/H_*$, chosen here as 20 with AccuracyBoost=2 for isocurvature modes. The ETA decomposition and its validity criterion (Eq. A5) determine when the averaging is trustworthy; for isocurvature modes with $k \\gg k_J$ the criterion fails and the code falls back to instantaneous values, a limitation the paper quantifies. The other load-bearing element is the abundance fitting formula (Eqs. 8–10), which connects the initial field value $\\phi_{\\rm ini}$ to the present abundance $\\Omega_a h^2$ across all masses, and through it converts the isocurvature amplitude $A_{\\rm iso}$ into the tensor-to-scalar ratio $r_a$ via $A_{\\rm iso}/A_s = (r_a/2)(M_{\\rm pl}/\\phi_{\\rm ini})^2$.","core_discovery":"The central discovery is that axion isocurvature is not merely a CDM analogue but a mass-dependent CMB observable with three regimes. For $m_a \\gg H_{\\rm eq}$ the axion isocurvature transfer function is indistinguishable from CDM isocurvature, so Planck's $\\beta_{\\rm iso} < 0.038$ bound applies directly and yields the $r f_{\\rm dm}$ inequality of Eq. (22). For $H_0 \\ll m_a \\lesssim H_{\\rm eq}$, Jeans suppression cuts off the spectrum with a $\\sim(\\ell_J/\\ell)^6$ falloff, breaking the degeneracy with CDM and leaving unique scale-dependent signatures, but also making the inflationary signal undetectable even under optimistic cosmic-variance-limited forecasts. For $m_a \\lesssim H_0$, the axion behaves as dark energy and the isocurvature signal peaks at the quadrupole, scaling as $(m_a/H_0)^2$; a detection there would contradict the standard inflationary production mechanism. The paper's quantitative claim is that the crossover where isocurvature constraints beat tensor constraints occurs at $m_a f_{\\rm dm}^2 \\gtrsim 10^{-26.4}\\,{\\rm eV}$, with an open coexistence window near $10^{-25}\\,{\\rm eV}$.","pith_inferences":["Not stated in the paper, but a direct corollary: the same $r f_{\\rm dm}$ calibration means that a deep null search in the coexistence window would push the tensor-to-scalar bound below the current BICEP value using CMB temperature and polarization data already being analyzed.","A testable extension the paper leaves implicit: run an exact (non-time-averaged) field integration over $m_a \\sim 10^{-30}{-}10^{-28}\\,{\\rm eV}$, where its own validity criterion (Eq. A5) fails for $k \\gg k_J$; the result would either confirm the order-of-magnitude accuracy claimed there or sharpen the light-mass bounds.","The abundance fit (Eqs. 8–10) is also a tool for QCD axion relic-density computations, since the same field-to-fluid switching errors that the paper documents for isocurvature modes have historically biased axion abundance estimates by factors of order unity.","If a future CMB experiment detects isocurvature in the light regimes, the paper's logic implies that the amplitude, not just the presence, of the signal could be used to infer how much the axion's effective potential or initial field velocity departs from the frozen quadratic picture; that inversion is not performed here."],"forward_implications":["In the dark-matter-like regime, any measurement of axion mass and abundance plus a null isocurvature search strengthens the inflation probe: Planck data already require $r f_{\\rm dm} < 0.076\\,(m_a/10^{-27}\\,{\\rm eV})^{-1/2}$ (95% CL).","For $m_a f_{\\rm dm}^2 \\gtrsim 10^{-26.4}\\,{\\rm eV}$, the isocurvature bound is the tighter of the two independent inflationary probes, superseding the BICEP tensor limit.","In the intermediate mass range $10^{-32} \\lesssim m_a/{\\rm eV} \\lesssim 10^{-28}$, the Jeans-scale cutoff makes the isocurvature spectrum distinguishable from CDM isocurvature, but even a cosmic-variance-limited CMB experiment cannot detect the inflationary signal; a detection would require non-inflationary generation.","In the dark-energy regime $m_a \\lesssim H_0$, only the CMB quadrupole carries the isocurvature signal, and the BICEP tensor bound forces it far below detectability, so any observed signal would indicate a breakdown of the frozen-field misalignment scenario.","The coexistence window near $m_a \\sim 10^{-25}\\,{\\rm eV}$, $f_{\\rm dm} \\gtrsim 0.1$ is a near-term discovery target: both axion isocurvature and tensor modes could appear just below current limits while simultaneously lowering $S_8$ to the weak-lensing preferred range."],"supporting_citations":[{"why":"Supplies the effective time average method that the paper extends to isocurvature modes.","marker":"[32]"},{"why":"Introduced AxiECAMB and its ETA treatment for adiabatic modes, which this paper extends with new initial conditions and accuracy criteria.","marker":"[34]"},{"why":"Prior calculation of axion abundance and isocurvature constraints whose normalization the paper finds to be off by a large factor.","marker":"[12]"},{"why":"Planck bound on uncorrelated CDM isocurvature, $\\beta_{\\rm iso} < 0.038$, the observational constraint converted into $r_a$ limits.","marker":"[26]"},{"why":"BICEP/Keck upper limit $r_{0.05} < 0.036$ used as the comparison benchmark for the isocurvature bounds.","marker":"[43]"},{"why":"Provided the analytic quadrupole scaling for dark-energy axion isocurvature that the paper recalibrates to $C_Q = 0.0097$.","marker":"[27]"},{"why":"Supplies the isocurvature gravitational-potential and low-$\\ell$ scaling used to calibrate the plateau coefficient $C_P$.","marker":"[41]"}],"fun_headline_variants":["Axion isocurvature outshines tensors for most masses","CMB axion signals: isocurvature beats B-modes","Axion masses split isocurvature into three regimes","Jeans suppression makes axion isocurvature unique","Coexistence window for axion isocurvature and tensors"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole mapping rests on the axion's symmetry being broken before inflation ends with the field frozen in a nearly quadratic potential, so that its present-day density and the size of its quantum fluctuations both follow from one initial field value; change that setup and the derived bounds on the inflationary energy scale no longer follow.","fun_headline_variants_meta":{"raw":{"variants":["Axion isocurvature outshines tensors for most masses","CMB axion signals: isocurvature beats B-modes","Axion masses split isocurvature into three regimes","Jeans suppression makes axion isocurvature unique","Coexistence window for axion isocurvature and tensors"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00066,"raw_usage":{"total_tokens":3168,"prompt_tokens":1244,"completion_tokens":1924,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":860,"completion_tokens_details":{"reasoning_tokens":1834}},"tokens_in":860,"tokens_out":1924,"duration_ms":12451,"temperature":1.0,"reasoning_tokens":1834,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T04:14:33.617165+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the isocurvature CMB power spectrum for $m_a = 10^{-28}\\,{\\rm eV}$, $f_{\\rm dm}=1$ with a switch at $m_a/H_* = 50$ and AccuracyBoost=3; if the low-$\\ell$ plateau or quadrupole moves by more than the quoted percent-level accuracy, the ETA-based numerical claim is not converged. Alternatively, detect isocurvature at $m_a \\sim 10^{-30}\\,{\\rm eV}$ at the level predicted by $r=0.036$; the paper predicts that signal is unobservable, so such a detection refutes its conclusion.","supporting_citations":[{"cited_title":"A Low CMB Quadrupole from Dark Energy Isocurvature Perturbations","cited_arxiv_id":"astro-ph/0406496","evidence_quote":"Supplies the effective time average method that the paper extends to isocurvature modes."},{"cited_title":"In the absence of metric sources, these free field oscillations can become much more rapid than the mass scale and invalidate the effective time average approach","cited_arxiv_id":null,"evidence_quote":"Introduced AxiECAMB and its ETA treatment for adiabatic modes, which this paper extends with new initial conditions and accuracy criteria."}],"review_version":1}