{"id":"0ef43af1-e8a9-49ac-952a-ccfd99fb72dc","arxiv_id":"2505.21276","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":13,"one_line_summary":"Axion-kaon scattering, previously omitted from axion thermalization calculations, exceeds the axion-pion channel near T=130 MeV and tightens the Planck Delta N_eff bound on f_a by about 30%.","lead":"This paper calculates axion production in the early universe and finds that axion-kaon scattering, previously ignored, dominates over axion-pion scattering at temperatures above about 130 MeV. This makes the axion hot dark matter bound from Planck data about 30% stronger, meaning axions must be lighter or more weakly coupled than previously thought.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ~30% tightening rests on Eq. (4)'s unvalidated dressing of the LO aK→πK vertex; no error propagation from the hadronic fit is provided.","rationale":"The paper's central number is f_a >= 3.18e7 GeV, a ~30% tightening from the pion-only bound. The calculation has two layers: measured meson-meson scattering data fix the unitarized T-matrix, and Eq. (4) appends the axion to that T-matrix. The first layer is well supported—Fig. 1 shows good fits and the pion cross-section agrees with the IAM result of Ref. [33]. The second layer is where the claim is least secure. Eq. (4) is not derived from QCD; it is an on-shell resummation ansatz, and the K*(892) enhancement that drives the effect is precisely the part of the result that depends on the unconstrained residue of the production amplitude at the resonance pole. The paper gives no estimate of how the 30% number changes if the off-shell prescription or the axion vertex at next order is varied. Table I's parameter errors are not propagated, so even within the model there is no error band. I agree with the reader that the direction is plausible and the pion-channel crosscheck is good, but the precise bound should not be used as firm until Eq. (4) is checked with an independent unitarization for aK→πK. The missing code link in [59] also prevents an immediate independent numerical check. Hence I recommend no change to the reader's CONDITIONAL verdict.","tokens_in":12197,"tokens_out":9491,"duration_ms":111247,"concrete_test":"Obtain the aK→πK partial waves from the same LO χPT vertices using an independent unitarization, e.g. the IAM prescription of Ref. [33] (M_IAM = M^(2) [1 - G T^(2)]^{-1} in the single-channel case, with the same G and LECs fitted to the data in Fig. 1), and recompute Γ_{aK→πK}(T) and the resulting f_a bound. If at T=130 MeV the rate changes by more than ~20% or the f_a bound moves by more than ~10%, Eq. (4)'s model dependence materially affects the headline 30% tightening; if it changes by less, the conditional is resolved in the paper's favor.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (4) is the single bridge from measured ππ/πK phase shifts to the claimed aK→πK enhancement. It assumes the axion production amplitude is the LO χPT vertex M^(2) resummed by the same on-shell two-meson T-matrix used for scattering, i.e. M_uni = T_uni (T^(2))^{-1} M^(2). This is an ansatz, not a QCD derivation; the K*(892) and K0*(700) peaks in Fig. 2 are inherited as poles of T_uni, and the residues at those poles are not constrained by the scattering data used to fit G and the LECs. A different off-shell prescription for G, a different production-vertex matching, or NLO axion couplings can move the resonance contribution. The paper's aπ→ππ cross-check against IAM (green dotted, Fig. 2 left) validates the method for pions but has no analogue for aK→πK. Moreover, Table I reports fit errors but no errors are propagated to Γ(T) or to f_a=3.18e7 GeV; the +0.04/-0.03 uncertainties in Fig. 3 reflect only the Planck ΔN_eff band, not the hadronic model. Footnote [59] states the Mathematica code 'can be downloaded from here' but no link is given, so the amplitudes cannot currently be independently re-evaluated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reassesses the thermalization rate of QCD axions in the hadronic phase below T_c. Using SU(3) chiral perturbation theory amplitudes for aP1→P2P3 processes, the authors unitarize both the two-meson scattering amplitudes and the axion production amplitudes via the same chiral-unitarization kernel, Eq. (4). They find that the aK→πK channel, which contains the K*(892) and K0*(700) resonances, has much larger cross sections than aπ→ππ for √s around 0.8–1 GeV (Fig. 2). Consequently, the axion thermalization rate Γ_a(T) is substantially enhanced for T≳100 MeV, and the aK→πK contribution exceeds the pion-only contribution around T≈130 MeV. Confronting the Planck 2018 ΔN_eff constraint, they obtain a tightened lower bound f_a ≥ 3.18×10^7 GeV (versus 2.45×10^7 GeV with pions only), a tightening of about 30%.","tokens_in":12610,"tokens_out":6815,"duration_ms":76000,"significance":"If the central result holds, the paper identifies a previously overlooked channel that materially changes the hot-dark-matter bound on axions. The approach is a standard chiral-unitarization framework, and the paper includes two valuable robustness checks: a crosscheck of the aπ→ππ amplitude against the inverse-amplitude-method result of Ref. [33] (green dotted curve in Fig. 2, left panel), and a cutoff analysis showing that for T<150 MeV more than 96% of the thermal rate comes from √s≤1.2 GeV, where the unitarized amplitudes are pinned to scattering data. These are real strengths. The main significance depends on Eq. (4), which is an ansatz for the unitarized production amplitude, and on the numerical reliability of the kaon-channel enhancement; both need further scrutiny before the 30% tightening can be taken as a firm quantitative prediction.","major_comments":[{"comment":"The unitarized axion production amplitude is constructed as M_uni = T^(2) [T^(2) − T^(4)LECs − T^(2) G T^(2)]^{-1} M^(2), i.e., the LO axion vertex dressed by the same T-matrix used for meson-meson scattering. This is an ansatz, not a derivation from QCD. The K*(892) and K0*(700) peaks in Fig. 2 are poles of the two-meson T-matrix, and the axion-production residues at these poles are not constrained by the phase-shift and inelasticity data that determine G and the LECs. The pion-channel crosscheck against the IAM result (green dotted, Fig. 2 left) validates the procedure for aπ→ππ, but there is no analogous crosscheck for aK→πK. Since the central claim of a 30% tightening rests on this unvalidated dressing of the kaon vertex, please quantify the sensitivity of Γ_a(T) and of f_a to (i) alternative off-shell prescriptions for the loop function G, and (ii) possible next-to-leading-order axion-meson couplings.","section":"Eq. (4) and the paragraph beginning 'Regarding the aP1→P2P3 reaction'"},{"comment":"Table I reports fit uncertainties for the subtraction constants and LECs, but these uncertainties are not propagated to Γ_a(T), ΔN_eff, or the final bound f_a ≥ 3.18^{+0.04}_{−0.03}×10^7 GeV. The error bars shown in Fig. 3 reflect only the Planck ΔN_eff band, not the hadronic-model uncertainty. Because the 30% tightening relative to 2.45×10^7 GeV is the central quantitative result, the absence of any hadronic error estimate makes the significance of the shift unclear. Please propagate the parameter covariance matrix through the thermal-rate calculation, or at least provide a conservative estimate of how the 30% shift and the aK/aπ crossover temperature depend on the hadronic-model uncertainties.","section":"Table I and the paragraph 'Significant enhancement...' / Fig. 3"}],"minor_comments":[{"comment":"The statement 'The mathematica code for the amplitudes can be downloaded from here' contains no URL, DOI, or ancillary-file identifier. Without a working link, the amplitudes cannot be independently re-evaluated; please provide a full link and, ideally, a version/checksum.","section":"Footnote [59]"},{"comment":"The partial-wave projection uses a factor 1/2(√2)^{N_i+N_f} with the rule N_i(N_f)=1 for identical particles, but the values of N_i and N_f for each of the four independent production types (especially aπ0→π0π0 and the charged kaon modes) are not explicitly listed. Please state them or give a table.","section":"Supplement, Eq. (A.7)"},{"comment":"The y-axis label is 'Γ_ch(T)/Γ_a(T)' with axis values from 0 to 1, while the text and caption describe 'relative contributions in percentage.' Please make the units consistent (fraction vs. percent) to avoid confusion about the 40% and 'exceeds' statements.","section":"Fig. 2, right panel"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and addresses an important question in axion cosmology. The central claim, however, hinges on an unvalidated production-vertex unitarization in the kaon channel and on the absence of hadronic error propagation. These issues are fixable in a revision, so I do not recommend rejection, but the 30% tightening should not be presented as firm until the sensitivity analysis is provided. The missing code link is also a reproducibility concern that should be fixed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"New result: the aK -> pi K channel, previously omitted, grows to dominate axion thermalization above about 130 MeV and tightens the Planck 2018 f_a bound by ~30% (3.18e7 vs 2.45e7 GeV). That is a real, quantitative new result, not a repackaging of pion-only calculations. The chiral unitarization machinery is standard, the fit to phase shifts in Fig. 1 looks reasonable, and the IAM crosscheck for the pion channel is a sensible validation. The cutoff stability check (96% of the rate below 1.2 GeV for T<150 MeV) is also welcome.\n\nThe soft spot is Eq. (4). The axion production amplitude is taken to be the LO vertex resummed by the same T-matrix that fits pi-pi and pi-K scattering. That's an ansatz, not a QCD derivation. The K*(892) and K0*(700) peaks in the cross sections come from the poles of that T-matrix, and nothing in the scattering data fixes the residues at those poles for the production process. A different off-shell prescription for the loop function G or NLO axion couplings could shift the resonance contribution. The pion case has the IAM cross-check; the kaon case has no analogue, so you're relying on the ansatz more heavily there.\n\nThe second issue is error propagation. Table I gives fit uncertainties for the LECs and subtraction constants, but they aren't propagated to Gamma(T) or to the final f_a bound. The +0.04/-0.03 in Fig. 3 is only the Planck band, so the hadronic model uncertainty on the 30% number is unquantified. That doesn't invalidate the direction, but it means the precise bound shouldn't be quoted as firm.\n\nMinor: footnote [59] says the Mathematica code can be downloaded 'from here' but no link is given. And 'systematical' overstates the channel set; you don't include a eta -> pi pi or heavier mesons, though the kaon is clearly the most plausible next channel.\n\nOverall: the direction is likely right, and the paper deserves a serious referee. I'd send it to peer review and ask for (i) a propagation of the hadronic fit errors to the final bound and (ii) a discussion of the model dependence of Eq. (4), ideally with a sensitivity estimate. The paper is useful for anyone relying on the axion hot dark matter bound.","headline":"The kaon channel plausibly dominates axion thermalization near 130 MeV and tightens the f_a bound by ~30%, but the central number rests on an unvalidated production-vertex ansatz and no error propagation from the hadronic fit.","tokens_in":13179,"tokens_out":2487,"would_cite":false,"duration_ms":24859,"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":"The paper claims that above $T\\simeq100$ MeV the $aK\\leftrightarrow\\pi K$ channel, not $a\\pi\\leftrightarrow\\pi\\pi$, dominates axion thermalization, tightening the hot-dark-matter bound on $f_a$ by about 30%.","keywords":["axion","axion thermalization rate","hot dark matter","extra relativistic degrees of freedom","chiral unitarization","meson-meson scattering","K*(892) resonance","QCD axion"],"falsifier":"A lattice-QCD computation of the $aK\\to\\pi K$ amplitude at center-of-mass energies near the $K^*(892)$ resonance, or a next-to-leading-order chiral calculation of the axion-kaon vertex, would settle the claim: if the resonance-enhanced kaon cross section comes out close to the pion-only one, the reported 30% tightening of $f_a$ would not survive.","tokens_in":11964,"feed_emoji":"🌌","tokens_out":11507,"duration_ms":115683,"temperature":0.7,"pith_summary":"The paper argues that the standard way of computing how axions stay in thermal contact with the early-universe bath—keeping only the $a\\pi\\leftrightarrow\\pi\\pi$ reaction—misses a much stronger process. By unitarizing chiral perturbation theory amplitudes, the authors find that $aK\\leftrightarrow\\pi K$ scattering, pumped up by the $K^*(892)$ and $K_0^*(700)$ resonances, has far larger cross sections, and above $T\\simeq100$ MeV it substantially raises the total axion thermalization rate. Around $T\\simeq130$ MeV the kaon channel overtakes the pion channel. If this calculation is right, the hot-dark-matter bound on the axion gets roughly 30% tighter: the axion decay constant must satisfy $f_a\\geq 3.18\\times 10^7$ GeV rather than $2.45\\times 10^7$ GeV under the measured bound on extra relativistic degrees of freedom.","feed_headline":"Kaon scattering tightens the axion hot-dark-matter bound","feed_subtitle":"Above 100 MeV, aK→πK overtakes pion scattering and raises the f_a lower bound by about 30 percent.","key_machinery":"The central object is the unitarized axion-meson partial-wave amplitude, Eq. (4): the leading-order chiral perturbation theory ($\\chi$PT) axion amplitude $\\vec M^{(2)}_{IJ}$ is dressed by the same inverse-amplitude combination $T^{(2)}_{IJ}\\,[T^{(2)}_{IJ}-T^{(4)\\,{\\rm LECs}}_{IJ}-T^{(2)}_{IJ}\\,G\\,T^{(2)}_{IJ}]^{-1}$ used for the meson-meson $T$-matrix, with $G$ the two-point loop function and subtraction constants and $O(p^4)$ low-energy constants fit to phase-shift and inelasticity data. This construction enforces two-body unitarity, generates the $f_0(500)$, $f_0(980)$, $\\rho(770)$, $K^*(892)$, and $K_0^*(700)$ resonances in the relevant partial waves, and couples $\\pi K$ with $\\eta K$ in the $I=1/2$ channels. The same dressed amplitudes feed the 12-dimensional phase-space integral, reduced to five numerical variables, that defines the axion thermalization rate.","core_discovery":"On the paper's own terms, the claim is that the dominant axion thermalization channel below the QCD crossover is not the pion-only channel but $aK\\leftrightarrow\\pi K$, once two-meson rescattering is resummed. The axion amplitude is built from the leading-order chiral axion vertex dressed by the same unitarized partial-wave $T$-matrix that describes $\\pi\\pi$ and $\\pi K$ scattering data; the resulting amplitudes carry the light resonances, especially $K^*(892)$ and $K_0^*(700)$ in the kaon channels. The integrated rates then show the kaon contribution exceeds 40% of the total for $T\\gtrsim110$ MeV and surpasses the pion contribution around $T\\gtrsim130$ MeV. Using the measured $\\Delta N_{\\rm eff}$, the lower bound on $f_a$ rises from $2.45\\times 10^7$ GeV (pion only) to $3.18\\times 10^7$ GeV (pion plus kaon), and the corresponding bound on the axion mass tightens by roughly the same factor.","pith_inferences":["Reading beyond the paper, the same resonance-dressing logic could be applied to other strange and heavier channels, such as axion-$\\eta$ or axion-nucleon couplings; if those channels show comparable enhancements, the $f_a$ bound would tighten further than the 30% reported here.","A direct lattice-QCD calculation of the $aK\\to\\pi K$ amplitude across the $K^*(892)$ region would provide a model-independent check of the unitarized dressing prescription, since the predicted enhancement is large enough to be visible.","The temperature-dependent shape of the thermalization rate changes above 100 MeV, which could alter not only the central $\\Delta N_{\\rm eff}\\to f_a$ mapping but also the decoupling history used in future cosmic microwave background forecasts.","One implicit consequence is that earlier axion and axion-like-particle constraints that assumed pion-only thermalization may need re-evaluation whenever the gluonic axion coupling is the operative interaction, not only for the QCD axion."],"forward_implications":["Below the QCD crossover, axion decoupling is governed by a total rate that is substantially larger than the pion-only rate for $T\\gtrsim100$ MeV, so pion-only calculations understate axion hot-dark-matter production in that window.","The lower bound $f_a\\geq 3.18\\times 10^7$ GeV (pion plus kaon channels) tightens the previously quoted $2.45\\times 10^7$ GeV by about 30%, and the corresponding upper bound on the axion mass moves downward by a comparable amount.","The long-standing assumption that $a\\pi\\leftrightarrow\\pi\\pi$ is the only dominant thermalization channel below $T_c$ is invalidated in the temperature range where the kaon channel contributes more than 40% of the total rate.","For $T$ below 150 MeV, more than 96% of the computed rate comes from $\\sqrt{s}\\leq 1.2$ GeV, so the central result does not rely on the behavior of the amplitudes above the fitted energy region."],"supporting_citations":[{"why":"Supplies the unitarization resummation recipe used to generate the resonances in the two-meson amplitudes.","marker":"[38]"},{"why":"Provides the companion derivation of the inverse-amplitude method that dresses the axion-meson vertex.","marker":"[39]"},{"why":"Supplies the $O(p^4)$ chiral Lagrangian whose terms define the leading-order axion-meson couplings and the low-energy constants.","marker":"[36]"},{"why":"Defines the two-point loop function $G$ with subtraction constants used in the unitarized amplitudes.","marker":"[40]"},{"why":"Gives the earlier inverse-amplitude calculation of $a\\pi\\to\\pi\\pi$ used as a crosscheck of the present pion-channel result.","marker":"[33]"},{"why":"Sets out the isospin-breaking suppression of the pion channel that motivates adding the kaon channel.","marker":"[34]"},{"why":"Supplies the phase-space reduction of the thermalization rate integral used in the numerical evaluation.","marker":"[35]"},{"why":"Provides the cosmic microwave background measurement of $\\Delta N_{\\rm eff}$ from which the $f_a$ bound is derived.","marker":"[22]"},{"why":"Supplies the cosmological inputs (Hubble rate and entropy) that convert the thermalization rate into $\\Delta N_{\\rm eff}$.","marker":"[57]"},{"why":"Supplies the lattice-averaged quark mass ratios that set the values of $r$ and $z$ in the axion couplings.","marker":"[60]"}],"fun_headline_variants":["Kaon scattering raises axion mass bound by 30 percent","Axion thermalization rate enhanced by kaon interactions","Kaon channel dominates axion thermalization above 130 MeV","Pion-only axion rate misses kaon enhancement at 100 MeV"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on assuming that the axion's interaction with pions and kaons is the simplest low-energy form prescribed by chiral symmetry, amplified by the same mathematical resummation that fits ordinary meson scattering data; that resummation is not derived from quantum chromodynamics and has never been tested by any axion-scattering measurement.","fun_headline_variants_meta":{"raw":{"variants":["Kaon scattering raises axion mass bound by 30 percent","Axion thermalization rate enhanced by kaon interactions","Kaon channel dominates axion thermalization above 130 MeV","Pion-only axion rate misses kaon enhancement at 100 MeV"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00072,"raw_usage":{"total_tokens":3275,"prompt_tokens":1033,"completion_tokens":2242,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":649,"completion_tokens_details":{"reasoning_tokens":2183}},"tokens_in":649,"tokens_out":2242,"duration_ms":21837,"temperature":1.0,"reasoning_tokens":2183,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:31:10.934938+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A lattice-QCD computation of the $aK\\to\\pi K$ amplitude at center-of-mass energies near the $K^*(892)$ resonance, or a next-to-leading-order chiral calculation of the axion-kaon vertex, would settle the claim: if the resonance-enhanced kaon cross section comes out close to the pion-only one, the reported 30% tightening of $f_a$ would not survive.","supporting_citations":[{"cited_title":"Non-perturbative Approach to effective chiral Lagrangians and Meson Interactions","cited_arxiv_id":"hep-ph/9803242","evidence_quote":"Supplies the unitarization resummation recipe used to generate the resonances in the two-meson amplitudes."}],"review_version":1}