Pith. Sign in

REVIEW

New cosmological mass limit on thermal relic axions

High Energy Physics - Phenomenology · Submitted 7 Apr 2005 · 47 references

Pith's verdict: A solid, honest recomputation of axion freeze-out that yields a real 1.05 eV bound and closes the hadronic axion window; the only soft spot is the model-dependence of the one-dimensional collapse, which the authors themselves flag.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on November 11, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 10 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-ph/0504059 (version 1)
Pith IDTFQTRW7G, Pith's permanent record of this paper and its review
CitationAarhus), A.Mirizzi (Bari U. & INFN, Bari & Munich, G.Raffelt (Max Planck Inst., Max Planck Inst.), Munich), S.Hannestad (Aarhus U.. "New cosmological mass limit on thermal relic axions." arXiv:hep-ph/0504059 (2005). https://pith.science/paper/TFQTRW7G
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Cosmological data push the thermal-relic axion below about one electronvolt, closing a window that supernova arguments left open.

The reading

The paper asks whether the same large-scale-structure argument that yields neutrino mass limits can be turned on axions that were once in thermal equilibrium in the early universe. Axions that interact with pions strongly enough stay in equilibrium until late times; after they decouple, they free-stream and smooth out small-scale density fluctuations exactly as light neutrinos do, so a hot-dark-matter component of them can be bounded by galaxy-survey and microwave-background data. The authors compute the axion freeze-out temperature from the axion-pion interaction rather than assuming it, then feed the resulting abundance and velocity distribution into a likelihood analysis of WMAP, SDSS, supernova and Lyman-$α$ data. For hadronic axions obeying the standard mass-decay-constant relation they obtain $m_a < 1.05$ eV at 95% CL, i.e. $f_a > 5.7×10^{6}$ GeV, which would close the hadronic axion window near $10^6$ GeV that survived the supernova 1987A arguments. The same machinery gives $∑ m_ν < 0.65$ eV for neutrinos, a useful calibration of how much constraining power the data set carries.

What carries the argument

The load-bearing object is the axion-pion absorption rate $⟨Γ_a⟩_T = A (C_{aπ}/(f_a f_π))^2 T^5 h(m_π/T)$ with $A = 0.215$ and $h(0)=1$, obtained by integrating the three pion channels $a+π^0 → π^+ + π^-$ and its charged partners over thermal phase space using the squared matrix element $∑|M|^2 = (9/4)(C_{aπ}/(f_a f_π))^2 (s^2+t^2+u^2-3m_π^4)$. Setting this absorption rate equal to the Hubble rate $H = [4π^3 g_*(T)/45]^{1/2} T^2/m_{Pl}$ defines the decoupling temperature, and the value of $g_*$ there fixes the relic number density through entropy conservation, $n_a = [g_{*S}({today})/g_{*S}({decoupling})] (n_γ/2)$. The resulting density and velocity distribution are then inserted into a modified CMBFAST Boltzmann code and marginalized over seven cosmological parameters with a simulated-annealing sampler, giving the exclusion contours in the $m_a$–$g_*$ and $m_a$–$f_a$ planes.

Extended reading notes

Core claim

The central claim is a new upper bound on the mass of a thermally produced axion: $m_a < 1.05$ eV at 95% CL, equivalently $f_a > 5.7×10^{6}$ GeV, for hadronic axions that obey the standard relation $m_a = 6.0$ eV $(f_a/10^{6}$ GeV$)^{-1}$ and the hadronic axion-pion coupling $C_{aπ} = (1-z)/[3(1+z)]$ with $z = 0.56$. The bound comes from structure formation, not from stellar energy loss or from a laboratory signal: an axion that was once in thermal equilibrium and decoupled while the universe was still hot behaves as hot dark matter, suppressing power on small scales. Because the axion-pion interaction is the dominant thermalization channel in the post-QCD epoch, the authors can compute the decoupling temperature as a function of $f_a$ (ranging from about 13 MeV to 145 MeV over the plotted range) and hence the present-day axion number density, which falls from roughly 74 cm$^{-3}$ to 38 cm$^{-3}$ as $f_a$ grows. That abundance, together with the free-streaming length set by the decoupling epoch, is what the cosmological data constrain. The authors state explicitly that the result closes the hadronic axion window of roughly $3×10^{5}$ GeV $≲ f_a ≲ 3×10^{6}$ GeV at generic coupling values, and they present the two-dimensional exclusion in the $m_a$–$f_a$ plane so that other axion models can be recovered by rescaling the vertical axis.

Load-bearing premise

The 1.05 eV number assumes both the standard axion mass-decay-constant relation and the specific hadronic axion-pion coupling with $z = 0.56$; if the axion-pion coupling is suppressed, the same data allow a heavier axion.

Editorial extensions

If this is right

  • Hadronic axions with masses of a few eV, previously tolerated in the supernova 1987A window around $f_a ∼ 10^{6}$ GeV, are excluded at 95% CL unless the axion-pion coupling is suppressed relative to the hadronic value.
  • The bound is independent of the axion-photon and axion-electron couplings, so it constrains models that stellar cooling arguments cannot reach.
  • The upper reach of the CAST helioscope search, roughly 1 eV, coincides with this cosmological limit, so a CAST detection in that range would be in tension with the structure-formation bound.
  • The same analysis pipeline reproduces $∑ m_ν < 0.65$ eV for neutrinos, giving a direct comparison of how strongly the data constrain a bosonic versus a fermionic hot relic of the same mass.
  • For non-hadronic models the exclusion contour shifts only through the rescaling of $C_{aπ}$, so the bound weakens in a way that can be read off directly from the two-dimensional plot.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Because the bound scales roughly as $C_{aπ}^2$ through the decoupling temperature, a model with $C_{aπ}$ smaller by a factor of a few would move the 1.05 eV limit up by a comparable factor, so the number is best read as a constraint on the combination $m_a$ and axion-pion coupling rather than on mass alone.
  • The axion density quoted here, 40 to 75 cm$^{-3}$, is below the neutrino density of about 112 cm$^{-3}$ per species, so an axion of the same mass contributes less hot dark matter than a neutrino; the tighter axion mass limit therefore reflects the bosonic distribution and the earlier decoupling as much as the abundance.
  • A future measurement of the axion-nucleon coupling, for instance from a next-generation supernova neutrino signal or a direct detection experiment, would fix the remaining freedom in the model and turn the contour plot into a single testable point.
  • If small-scale structure data improve below the $k = 0.15\ h$ Mpc$^{-1}$ cutoff used here, the Lyman-$α$ sensitivity noted for $g_* ≳ 60$ could extend the reach to lighter axions or to more weakly coupled ones.
Original abstract

Observations of the cosmological large-scale structure provide well-established neutrino mass limits. We extend this argument to thermal relic axions. We calculate the axion thermal freeze-out temperature and thus their cosmological abundance on the basis of their interaction with pions. For hadronic axions we find a new mass limit $m_a<1.05$ eV (95% CL), corresponding to a limit on the axion decay constant of $f_a>5.7\times 10^6$ GeV. For other models this constraint is significantly weakened only if the axion-pion coupling is strongly suppressed. For comparison we note that the same approach leads to $\sum m_\nu<0.65$ eV (95% CL) for neutrinos.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 4,877 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.11 · states limitation 0.78 · human participants 0.01 · controlled comparison 0.13 · randomized assignment 0.01 · states quantitative result 0.99 · compares with prior results 0.91 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

Observations of the cosmological large-scale structure provide well-established neutrino mass limits. We extend this argument to thermal relic axions. We calculate the axion thermal freeze-out temperature and thus their cosmological abundance on the basis of their interaction with pions. For hadronic axions we find a new mass limit $m_a<1.05$ eV (95% CL), corresponding to a limit on the axion decay constant of $f_a>5.7\times 10^6$ GeV. For other models this constraint is significantly weakened only if the axion-pion coupling is strongly suppressed. For comparison we note that the same approach leads to $\sum m_\nu<0.65$ eV (95% CL) for neutrinos.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Observations of the cosmological large-scale structure provide well-established neutrino mass limits. background 1.00 0.05 0.02 asserted 1.00
s1 We extend this argument to thermal relic axions. contribution 0.95 0.77 0.05 asserted 0.92
s2 We calculate the axion thermal freeze-out temperature and thus their cosmological abundance on the basis of their interaction with pions. method 0.92 0.96 0.04 asserted 1.00
s3 For hadronic axions we find a new mass limit $m_a<1.05$ eV (95% CL), corresponding to a limit on the axion decay constant of $f_a>5.7\times 10^6$ GeV. result 1.00 0.96 0.90 asserted 1.00
s4 For other models this constraint is significantly weakened only if the axion-pion coupling is strongly suppressed. result 0.39 0.76 0.07 asserted 0.75
s5 For comparison we note that the same approach leads to $\sum m_\nu<0.65$ eV (95% CL) for neutrinos. result 0.93 0.87 0.05 asserted 0.99

What kind of paper the classifier reports this is: theoretical result 1.00 · evidence from proof or derivation 0.58.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Dimensionless axion absorption rate h(µ) as defined by Eq. (9)
Figure 2. Effective number of thermal degrees of freedom in the early universe during the post-QCD epoch, assuming the particle content listed in
Figure 3. Average axion absorption rate for fa = 107 GeV and cosmic expansion rate as a function of the cosmic temperature. Rates are in units of T 2/mPl
Figure 4. Axion freeze-out, assuming the a-π coupling of Eq. (5). Top: Decoupling temperature. Middle: Corresponding g∗. Bottom: Present-day axion density
Figure 5. Likelihood contours for the allowed axion parameters. Everything to the right of the dark shaded region is excluded at the 68% CL, and everything to the right of the light shaded region is excluded at the 95% CL. Contours of equal axion mass with the indicated ma values in eV are shown as dashed lines. Upper panel: All data included. Lower panel: Lyman-α data excluded. In
Figure 6. Likelihood contours including all data, equivalent to the upper panel of
Figure 7. Likelihood contours including all data, equivalent to the upper panel of
Figure 8. One-dimensional likelihood analysis for the hadronic axion model obeying Eq. (1). Shown is ∆χ 2 = χ 2 − χ 2 0 as a function of ma where χ 2 0 refers to the best-fit model. In this 1D analysis the 95% excluded region is at ∆χ 2 > 4, and the 68% excluded region at ∆χ 2 > 1, as indicated by the horizontal dashed lines

Figures reproduced from arXiv:hep-ph/0504059 by the authors.

References
How later papers used it

60 later papers cite it

Evidence: substantial. 51 of the citing papers share no author with it, and 11 of those carry a confident label from Pith's citation classifier.

Pith read 16 of the 108 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 3 papers · 1 by the same authors
  • The crucial quantity that is needed to extract the HDM bound is the axion decoupling temperature, $T_D$, obtained via the freeze-out condition (following the same criterium as in [ref] ) [equation] Here, $H(T) = √(4π^3 g_⋆(T) / 45) T^2 / m_pl$ is the Hubble rate (assuming a radiation dominated Universe) in terms of the Planck mass $m_pl = 1.22 ×10^19$ GeV and the effective number of relativistic degrees of freedom, $g_⋆(T)$, while $Γ_a$ is the axion thermalization rate entering the Boltzmann equation [equation] where $n_a^eq = (ζ_3 / π^2) T^3$ and $f_i = 1/(e^E_i/T - 1)$.

    Breakdown of chiral perturbation theory for the axion hot dark matter bound (2021) · Introduction · 0.94

  • For the perturbative calculation up to NLO, we expand the squared amplitude as $∑|ℳ|^2≃∑|ℳ^(2)|^2+∑2ℳ^(2)Re(ℳ^(4))$, and then the perturbative axion rate will be casted into the following form [ref] [equation] where the dimensionless functions $h_lo$ and $h_nlo$ denote the contributions from $∑|ℳ^(2)|^2$ and $∑2Re(ℳ^(2)ℳ^(4))$, respectively.

    Axion-pion scattering at finite temperature in chiral perturbation theory and its influence in axion thermalization (2023) · Brief discussions on the axion thermalization rate and its cosmology implication · 0.96

By the same authors

Uses its data or values 1 of 5 papers with this label
Compares with it 3 of 5 papers with this label · 2 by the same authors

By the same authors

  • Said that its latest limit on the axion mass is almost identical to one that some of its authors derived several years ago in this paper.

    Axion hot dark matter bounds (2008) (same authors) · 0.97 · paraphrased: the licence does not allow quoting

  • Said that its axion mass limit of 1.4 eV in the absence of neutrino masses is significantly weaker than the 1.05 eV limit found earlier in this paper.

    Cosmological constraints on neutrino plus axion hot dark matter (2007) (same authors) · 1.00 · paraphrased: the licence does not allow quoting

How the literature uses it

55 arXiv papers cite this in 100 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 62 sentences without a confident label

Builds on it 2 papers
Uses its method or tool 3 papers
  • Therefore, we limit ourselves to the traditional treatment of axion thermalization presented in [ref] and comment later on this important point.

    Thermal axions with multi-eV masses are possible in low-reheating scenarios · Axion decoupling in standard cosmology · 0.94

  • The crucial quantity that is needed to extract the HDM bound is the axion decoupling temperature, $T_D$, obtained via the freeze-out condition (following the same criterium as in [ref] ) [equation] Here, $H(T) = √(4π^3 g_⋆(T) / 45) T^2 / m_pl$ is the Hubble rate (assuming a radiation dominated Universe) in terms of the Planck mass $m_pl = 1.22 ×10^19$ GeV and the effective number of relativistic degrees of freedom, $g_⋆(T)$, while $Γ_a$ is the axion thermalization rate entering the Boltzmann equation [equation] where $n_a^eq = (ζ_3 / π^2) T^3$ and $f_i = 1/(e^E_i/T - 1)$.

    Breakdown of chiral perturbation theory for the axion hot dark matter bound · Introduction · 0.94

  • The explicit freeze out condition advocated in Ref. [ref] to determine the axion decoupling temperature is [equation] where $H(T)=T^2√(4π^3g_*(T)/45)/m_Pl$ denotes the Hubble expansion parameter, with $g_*(T)$ the effective number of relativistic thermal degrees of freedom (d.o.f) and $m_Pl$ the Plank mass. $Γ_a(T)$ in Eq. [ref] stands for the averaged axion thermalization rate in the pion thermal bath, and it can be calculated via [equation] where the sum runs over all the possible $aπ$ reaction channels, the axion number density in equilibrium is $n_a^eq=ζ(3) T^3/π^2$ with $ζ$ the Riemann Zeta function, and the phase space integral is [equation] The recipe given in Ref [ref] will be used to calculate the above phase space integral.

    Axion-pion scattering at finite temperature in chiral perturbation theory and its influence in axion thermalization · Brief discussions on the axion thermalization rate and its cosmology implication · 0.97

Uses its data or values 5 papers
  • The thermally averaged rates for these processes were first computed by [ref] and updated by [ref] , which found [equation] where $x_π ≡m_π/T$, $A=0.215$ and $h(x)$ is a numerical function normalized to $h(0)=1$ that decays exponentially below $T=m_π$ (see Appendix [ref]).

    The DFSZ axion in the CMB · Interactions with pions · 0.83

  • For instance, for the hadronic axion upper mass bound quoted in Ref. [ref] , i.e. $m_a∼1.05$ eV, the axion decouples at $T_D∼64$ MeV, at which $g_⋆S(T_D)≃15.24$ and $Ω_ah^2 ≃0.0053$.

    An improved cosmological bound on the thermal axion mass · Cosmological Constraints · 0.87

  • The processes to consider are $π^0π^±→a π^±$ and $π^+π^-→a π^0$, with an approximate total scattering rate for both given by [ref] :

    New cosmological bounds on axions in the XENON1T window · Background · 0.84

Compares with it 4 papers
Cites it as support 1 paper
Cites it as background 18 papers
  • The LO calculation was performed in Refs. [ref] , while Ref. [ref] considered the impact of NLO corrections in order to assess the convergence of the chiral expansion.

    Axion-pion thermalization rate in unitarized NLO chiral perturbation theory · Introduction · 0.89

  • From the standpoint of the axion solution to the strong CP problem, an unavoidable process arises from the model-independent coupling to gluons, $α_s/8π a/f_a G G̃$.[Other thermalization channels arise from model-dependent axion couplings to photons [ref] , SM quarks [ref] and leptons [ref] .] For $T_D ≳1$ GeV thermal axion production proceeds via its scatterings with gluons in the quark-gluon plasma [ref] , while for $T_D ≲1$ GeV processes involving pions and nucleons must be considered [ref] .

    Breakdown of chiral perturbation theory for the axion hot dark matter bound · Introduction · 0.89

  • Quantifying the production rate across the QCD phase transition (QCDPT) in the strongly coupled regime is far from being straightforward, and axion production was analyzed either exclusively above the confinement scale [ref] or below [ref] .

    Thermal axion production at low temperatures: a smooth treatment of the QCD phase transition · Introduction · 0.82

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.