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
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
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
- 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.
The abstract, sentence by sentence
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
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| 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 was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:hep-ph/0504059 by the authors.
References
How later papers used it
60 later papers cite it
3 later papers by other authors are listed below: 2 used its method or tool, 1 used its data or values and 1 compared its own results with it. 1 of them appears under more than one label. 3 more papers by the same authors are listed: 1 used its method or tool and 2 compared their own results with it.
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
-
Limited itself to the traditional treatment of axion thermalization presented in this paper.
Thermal axions with multi-eV masses are possible in low-reheating scenarios (2021) (same authors) · 0.94 · paraphrased: the licence does not allow quoting
Uses its data or values 1 of 5 papers with this label
-
Used the hadronic axion upper mass bound of about 1.05 eV quoted in this paper.
An improved cosmological bound on the thermal axion mass (2007) · 0.87 · paraphrased: the licence does not allow quoting
Compares with it 3 of 5 papers with this label · 2 by the same authors
-
We have checked that $h_LO$ reproduces the result of Ref. [ref] within percent accuracy.
Breakdown of chiral perturbation theory for the axion hot dark matter bound (2021) · Introduction · 1.00
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
Papers citing it · 62 sentences without a confident label
Builds on it 2 papers
-
We return to this topic to extend previous studies by some of us [ref] in several ways.
Cosmological constraints on neutrino plus axion hot dark matter · Introduction · 0.99
-
This figure can also be applied to thermally produced axions, like in Refs. [ref] .
Cosmological constraints on a light non-thermal sterile neutrino · Mass/temperature bounds in the thermal case · 0.85
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
-
Numerically our latest limit on $m_a$ [ref] is almost identical to one that some of us have derived several years ago [ref] .
Axion hot dark matter bounds · Introduction · 0.97
-
Our neutrino mass bound $∑m_ν<0.65$ eV (95% C.L.) in the absence of axions is identical to that derived by some of us in Ref. [ref] , whereas the axion mass limit $m_a<1.4$ eV in the absence of neutrino masses found here is significantly weaker than the $1.05$ eV limit found earlier [ref] .
Cosmological constraints on neutrino plus axion hot dark matter · Results · 1.00
-
We have checked that $h_LO$ reproduces the result of Ref. [ref] within percent accuracy.
Breakdown of chiral perturbation theory for the axion hot dark matter bound · Introduction · 1.00
Cites it as support 1 paper
-
Before proceeding let us stress that we can neglect production processes like $ππ→π$ in the primordial hot hadronic gas [ref] because of the $$ limit [ref].
Axions and saxions from the primordial supersymmetric plasma and extra radiation signatures · Thermal axion production · 0.97
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.