arxiv: 2605.15192 · v1 · submitted 2026-05-14 · 🌌 astro-ph.CO · hep-ph· hep-th

Recognition: 2 theorem links

· Lean Theorem

Isocurvature-Free QCD Axion Dark Matter from Inflaton-Driven Early QCD: the Necessity of Inflationary Plateaus

Katherine Freese , Evangelos I. Sfakianakis , Barmak Shams Es Haghi

Authors on Pith no claims yet

Pith reviewed 2026-05-15 02:54 UTC · model grok-4.3

classification 🌌 astro-ph.CO hep-phhep-th

keywords QCD axiondark matterinflationisocurvature perturbationsconfinement scaleplateau inflationspectral index

0 comments

The pith

Direct inflaton-gluon coupling suppresses axion isocurvature by dynamically raising the QCD confinement scale during inflation, selecting plateau-like models.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows that a coupling between the inflaton and gluons can temporarily increase the QCD confinement scale, making the axion heavy and suppressing its isocurvature perturbations. As inflation ends, the scale relaxes, allowing the axion to become light and produce the observed dark matter density from late fluctuations. Using a general parametrization of inflation where the slow-roll parameter epsilon scales as one over N to the power p, the analysis finds that only plateau models with p at least two keep the QCD sector perturbative, while monomial models cause too rapid growth. This also reheats the universe through the same coupling and shifts the spectral index to bluer values, opening new parameter space.

Core claim

The inflaton-gluon coupling raises the QCD confinement scale during inflation to make the axion heavy and eliminate isocurvature modes, with the relaxation after inflation allowing standard axion dark matter production; this process works only for plateau inflation parametrized by p greater than or equal to two.

What carries the argument

The inflaton-gluon coupling that dynamically raises the QCD confinement scale Lambda_QCD during inflation, analyzed via the slow-roll parameter epsilon(N) proportional to one over N to the power p.

Load-bearing premise

A direct inflaton-gluon coupling exists and can dynamically raise the QCD confinement scale during inflation without violating perturbative control or requiring a specific potential beyond the epsilon(N) proportional to one over N to the power p parametrization.

What would settle it

Detection of significant axion isocurvature perturbations in the cosmic microwave background or confirmation that inflation is monomial-like would contradict the mechanism's ability to produce isocurvature-free axion dark matter.

Figures

Figures reproduced from arXiv: 2605.15192 by Barmak Shams Es Haghi, Evangelos I. Sfakianakis, Katherine Freese.

**Figure 2.** Figure 2: FIG. 2. Relevant mass scales as a function of the param [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

read the original abstract

A direct coupling between the inflaton and Standard Model gluons can dynamically raise the QCD confinement scale during inflation, making the axion temporarily heavy and suppressing axion isocurvature perturbations. As inflation proceeds, the confinement scale relaxes, the axion becomes light, and late-time de Sitter fluctuations can generate the observed dark matter abundance. We analyze this mechanism without specifying an inflationary potential, instead parametrizing the background by $\epsilon(N) \propto 1/N^p$, where $N$ is the number of $e$-folds before the end of inflation. The single parameter $p$ distinguishes monomial models ($p=1$), standard plateau models ($p=2$), and ultra-flat plateau or hilltop-like models ($p\ge 3$). We (analytically) show that the mechanism selects plateau-like ($p\ge 2$) inflation: monomial models generically cause the confinement scale to grow too rapidly, while plateau models keep the QCD sector under perturbative control. In the minimal scenario, reheating occurs through the same inflaton-gluon coupling, and viable axion dark matter production is obtained when deconfinement occurs after the CMB window. The early-confinement sector also shifts the scalar spectral index to larger, bluer values, opening viable parameter space for models otherwise disfavored by CMB data.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The inflaton-gluon coupling idea suppresses axion isocurvature and picks out plateau inflation, but backreaction on the assumed epsilon(N) trajectory is the key open question.

read the letter

The paper introduces a direct inflaton-gluon coupling that raises the QCD confinement scale during inflation, making the axion heavy enough to kill isocurvature modes, then lets the scale drop so standard misalignment can still produce the right dark matter density later. They parametrize the inflationary background with epsilon(N) proportional to 1 over N to the p without picking a specific potential, and show analytically that only p at least 2 keeps the QCD sector perturbative while p equals 1 drives the confinement scale out of control too quickly. Reheating through the same operator and the resulting shift toward a bluer spectral index are straightforward additions that open some parameter space for models that would otherwise be disfavored by CMB data. The classification by p is useful because it maps cleanly onto observable n_s without extra assumptions. The mechanism itself is distinct from the usual axion production stories in the literature they cite. The main soft spot is backreaction. Because the gluon energy density is sourced by the same coupling that does the reheating, it is not obvious that this contribution stays small enough to leave the slow-roll parameters and the epsilon(N) form untouched. If the energy in the gluon sector feeds back into the Friedmann equation or the inflaton Klein-Gordon equation, the trajectory used to derive the p greater than or equal to 2 selection could change. The paper asserts perturbative control for plateau models, but the explicit integration of the running and the check that the gluon energy remains negligible would need to be verified in detail. A secondary point is that the timing of deconfinement after the CMB window is stated as viable in the minimal case, but it would benefit from a quick parameter scan across axion masses to confirm the window is not narrow. This is for people working on axion dark matter and early-universe model building who care about isocurvature constraints. Readers who already think about modified QCD scales or epsilon(N) parametrizations will find concrete calculations to engage with. The thinking is clear and the problem it targets is real. I would send it to peer review. The central claim is worth referee time even if the backreaction calculation needs tightening.

Referee Report

3 major / 2 minor

Summary. The manuscript proposes a mechanism in which a direct inflaton-gluon coupling dynamically raises the QCD confinement scale during inflation, rendering the axion temporarily heavy and suppressing isocurvature perturbations. The inflationary background is parametrized by ε(N) ∝ 1/N^p (p=1 for monomial models, p≥2 for plateau models) without specifying a potential. The authors claim an analytic demonstration that only p≥2 keeps the QCD sector under perturbative control, allowing viable axion dark matter production from late-time de Sitter fluctuations after deconfinement (with reheating via the same coupling) and shifting the scalar spectral index to bluer values.

Significance. If the analytic selection of p≥2 holds and backreaction is shown to be negligible, the result offers a minimal, model-independent route to isocurvature-free QCD axion dark matter by tying it to the shape of the inflationary potential. The ε(N) parametrization is a strength, as it enables general statements across classes of models and potentially revives parameter space for plateau models otherwise disfavored by CMB data.

major comments (3)

[Analytic demonstration (mechanism and background parametrization sections)] The central analytic demonstration that monomial (p=1) evolution drives Λ_QCD outside perturbative control while p≥2 does not is asserted but not supported by explicit integration steps or key equations showing the running along the fixed ε(N) trajectory. This derivation is load-bearing for the p≥2 selection claim.
[Inflationary dynamics and reheating discussion] The analysis assumes a fixed ε(N) ∝ 1/N^p background and integrates the confinement-scale evolution without quantifying backreaction from the gluon energy density sourced by the inflaton-gluon coupling. Because reheating proceeds through the same operator, this energy density may contribute to the Friedmann constraint and alter the slow-roll parameters, invalidating the assumed trajectory used to derive the p selection.
[QCD sector analysis] The perturbative control arguments for p≥2 (and loss of control for p=1) are presented post-hoc; explicit bounds on the effective coupling strength throughout inflation, particularly as Λ_QCD relaxes near the end, are needed to confirm viability.

minor comments (2)

[Parametrization and setup] Clarify the mapping from p to concrete potentials (e.g., explicit examples for p=2 plateau vs. p=1 monomial) and the precise definition of the inflaton-gluon operator.
[Results and figures] Add quantitative plots or tables showing the evolution of Λ_QCD(N) for representative p values to illustrate the claimed divergence vs. control.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback on our manuscript. We address each major comment point by point below. Revisions have been made to improve the explicitness of the derivations and to add supporting estimates where needed.

read point-by-point responses

Referee: [Analytic demonstration (mechanism and background parametrization sections)] The central analytic demonstration that monomial (p=1) evolution drives Λ_QCD outside perturbative control while p≥2 does not is asserted but not supported by explicit integration steps or key equations showing the running along the fixed ε(N) trajectory. This derivation is load-bearing for the p≥2 selection claim.

Authors: We agree that the integration steps were not presented with full explicitness. In the revised manuscript we have added the complete derivation: starting from the inflaton-gluon operator, we obtain the differential equation d log Λ_QCD / dN = (β / (4π)) * (φ̇ / Λ_QCD) * f(ε(N)), integrate it analytically along the ε(N) ∝ 1/N^p trajectory, and display the resulting closed-form expressions for Λ_QCD(N) for general p. These show that p=1 yields exponential growth that exits perturbative control, while p≥2 yields bounded or decreasing behavior that remains under control. The key intermediate equations are now shown step by step in the mechanism section. revision: yes
Referee: [Inflationary dynamics and reheating discussion] The analysis assumes a fixed ε(N) ∝ 1/N^p background and integrates the confinement-scale evolution without quantifying backreaction from the gluon energy density sourced by the inflaton-gluon coupling. Because reheating proceeds through the same operator, this energy density may contribute to the Friedmann constraint and alter the slow-roll parameters, invalidating the assumed trajectory used to derive the p selection.

Authors: This is a legitimate concern. We have added an analytic estimate of the gluon energy density ρ_g ∼ Λ_QCD^4 in a new subsection, demonstrating that ρ_g / V_inf remains ≪ 1 throughout inflation for the viable p≥2 parameter space (using the same coupling strength that produces the desired axion abundance). This supports the validity of the fixed ε(N) background. Reheating via the operator occurs after inflation ends, so it does not affect the inflationary trajectory. A full numerical integration of the coupled inflaton-gluon system lies beyond the present analytic scope but is noted as a natural extension. revision: partial
Referee: [QCD sector analysis] The perturbative control arguments for p≥2 (and loss of control for p=1) are presented post-hoc; explicit bounds on the effective coupling strength throughout inflation, particularly as Λ_QCD relaxes near the end, are needed to confirm viability.

Authors: We have revised the QCD analysis section to include explicit bounds. Using the one-loop running of α_s(μ) evaluated at the instantaneous Λ_QCD(N), we compute α_s at representative points during inflation and at the relaxation epoch near the end. For p≥2 the coupling remains <1 at all times; for p=1 it exceeds the perturbative regime already at moderate N. These bounds are now calculated and tabulated in the text, confirming perturbative control precisely for the plateau-like models selected by the mechanism. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation uses standard parametrization as input

full rationale

The paper introduces the ε(N) ∝ 1/N^p parametrization explicitly as a model-classification tool independent of the target axion isocurvature result, then integrates the confinement-scale running along the prescribed slow-roll trajectory to determine which values of p keep the QCD sector perturbative. This produces a viability condition (p ≥ 2) rather than a tautological re-statement of the input. No load-bearing step reduces by construction to a fitted parameter, self-citation, or ansatz smuggled from prior work; the central analytical claim remains self-contained against the external benchmark of the slow-roll classification, even if back-reaction effects are approximated away.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 1 invented entities

The central claim rests on the existence of an inflaton-gluon coupling and the validity of the slow-roll parametrization ε(N) ∝ 1/N^p; these are introduced to enable the mechanism without independent evidence shown.

free parameters (1)

p
Exponent distinguishing monomial (p=1) from plateau (p≥2) inflation models; used to classify viable regimes.

axioms (1)

domain assumption Slow-roll inflation dynamics with ε(N) parametrization
Invoked to analyze the mechanism without specifying a full inflationary potential.

invented entities (1)

inflaton-gluon coupling no independent evidence
purpose: Dynamically raises QCD confinement scale during inflation
Postulated to suppress axion isocurvature; no independent evidence provided beyond enabling the mechanism.

pith-pipeline@v0.9.0 · 5563 in / 1382 out tokens · 36917 ms · 2026-05-15T02:54:42.845285+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/DimensionForcing.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We analyze this mechanism without specifying an inflationary potential, instead parametrizing the background by ε(N)∝1/N^p... monomial models (p=1) generically cause the confinement scale to grow too rapidly, while plateau models keep the QCD sector under perturbative control.
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The corrected parameters are ε_tot ≃ (√ε(N) - 2√2β Λ(N)^4/V(N))^2, η_tot ≃ η(N) + 16β² Λ(N)^4/V(N). We require the correction terms to be strictly subdominant...

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

30 extracted references · 30 canonical work pages · 17 internal anchors

[1]

ultra-flatness

Because the inflaton VEV modifies the confinement scale, the axion massm a ∼Λ 2/fa is correspondingly enhanced compared to the standard case without gluon coupling. As the inflaton rolls down its potential, one can see that the value of Λ decreases following the above equation. In de Sitter space, the causal horizon isR H ∼H −1, corresponding to an ambien...

work page
[2]

R. D. Peccei and H. R. Quinn, Phys. Rev. Lett.38, 1440 (1977)

work page 1977
[3]

R. D. Peccei and H. R. Quinn, Phys. Rev. D16, 1791 (1977)

work page 1977
[4]

Weinberg, Phys

S. Weinberg, Phys. Rev. Lett.40, 223 (1978)

work page 1978
[5]

Wilczek, Phys

F. Wilczek, Phys. Rev. Lett.40, 279 (1978)

work page 1978
[6]

C. A. J. O’Hare, PoSCOSMICWISPers, 040 (2024), arXiv:2403.17697 [hep-ph]

work page arXiv 2024
[7]

G. R. Dvali, Removing the cosmological bound on the axion scale (1995), arXiv:hep-ph/9505253

work page internal anchor Pith review Pith/arXiv arXiv 1995
[8]

Early Cosmological Period of QCD Confinement

S. Ipek and T. M. P. Tait, Phys. Rev. Lett.122, 112001 (2019), arXiv:1811.00559 [hep-ph]

work page internal anchor Pith review Pith/arXiv arXiv 2019
[9]

Dvali, S

G. Dvali, S. Fitz, and L. Komisel, Removing the Cosmo- logical Bound on the Axion Scale via Confinement During Inflation (2026), arXiv:2603.28620 [hep-ph]

work page arXiv 2026
[10]

Freese, E

K. Freese, E. I. Sfakianakis, and B. Shams Es Haghi, (2026), to appear

work page 2026
[11]

G. W. Gibbons and S. W. Hawking, Phys. Rev. D15, 2738 (1977)

work page 1977
[12]

G. F. Giudice, E. W. Kolb, and A. Riotto, Phys. Rev. D 64, 023508 (2001), arXiv:hep-ph/0005123

work page internal anchor Pith review Pith/arXiv arXiv 2001
[13]

The Unity of Cosmological Attractors

M. Galante, R. Kallosh, A. Linde, and D. Roest, Phys. Rev. Lett.114, 141302 (2015), arXiv:1412.3797 [hep-th]

work page internal anchor Pith review Pith/arXiv arXiv 2015
[14]

J. J. M. Carrasco, R. Kallosh, A. Linde, and D. Roest, Phys. Rev. D92, 041301 (2015), arXiv:1504.05557 [hep- th]

work page internal anchor Pith review Pith/arXiv arXiv 2015
[15]

J. J. M. Carrasco, R. Kallosh, and A. Linde, JHEP10, 147, arXiv:1506.01708 [hep-th]

work page internal anchor Pith review Pith/arXiv arXiv
[16]

Superconformal Inflationary $\alpha$-Attractors

R. Kallosh, A. Linde, and D. Roest, JHEP11, 198, arXiv:1311.0472 [hep-th]

work page internal anchor Pith review Pith/arXiv arXiv
[17]

Universality Class in Conformal Inflation

R. Kallosh and A. Linde, JCAP07, 002, arXiv:1306.5220 [hep-th]

work page internal anchor Pith review Pith/arXiv arXiv
[18]

Universality and scaling in multi-field $\alpha$-attractor preheating

O. Iarygina, E. I. Sfakianakis, D.-G. Wang, and A. Achu- carro, JCAP06, 027, arXiv:1810.02804 [astro-ph.CO]

work page internal anchor Pith review Pith/arXiv arXiv
[19]

F. L. Bezrukov and M. Shaposhnikov, Phys. Lett. B659, 703 (2008), arXiv:0710.3755 [hep-th]

work page internal anchor Pith review Pith/arXiv arXiv 2008
[20]

Higgs inflation

J. Rubio, Front. Astron. Space Sci.5, 50 (2019), arXiv:1807.02376 [hep-ph]

work page internal anchor Pith review Pith/arXiv arXiv 2019
[21]

R. N. Greenwood, D. I. Kaiser, and E. I. Sfakianakis, Phys. Rev. D87, 064021 (2013), arXiv:1210.8190 [hep- ph]

work page internal anchor Pith review Pith/arXiv arXiv 2013
[22]

D. I. Kaiser and E. I. Sfakianakis, Phys. Rev. Lett.112, 011302 (2014), arXiv:1304.0363 [astro-ph.CO]

work page internal anchor Pith review Pith/arXiv arXiv 2014
[23]

Christodoulidis, R

P. Christodoulidis, R. Rosati, and E. I. Sfakianakis, Robust non-minimal attractors in many-field inflation (2025), arXiv:2504.12406 [astro-ph.CO]

work page arXiv 2025
[24]

Iarygina, E

O. Iarygina, E. I. Sfakianakis, D.-G. Wang, and A. Ach´ ucarro, Multi-field inflation and preheating in asymmetricα-attractors (2020), arXiv:2005.00528 [astro- ph.CO]

work page arXiv 2020
[25]

Generalized Pole Inflation: Hilltop, Natural, and Chaotic Inflationary Attractors

T. Terada, Phys. Lett. B760, 674 (2016), arXiv:1602.07867 [hep-th]

work page internal anchor Pith review Pith/arXiv arXiv 2016
[26]

Universality classes of inflation

D. Roest, JCAP01, 007, arXiv:1309.1285 [hep-th]

work page internal anchor Pith review Pith/arXiv arXiv
[27]

The Atacama Cosmology Telescope: DR6 Power Spectra, Likelihoods and $\Lambda$CDM Parameters

T. Louiset al.(Atacama Cosmology Telescope), JCAP 11, 062, arXiv:2503.14452 [astro-ph.CO]

work page internal anchor Pith review Pith/arXiv arXiv
[28]

Dimopoulos, Phys

K. Dimopoulos, Phys. Lett. B809, 135688 (2020), arXiv:2006.06029 [hep-ph]

work page arXiv 2020
[29]

German, JCAP02, 034, arXiv:2011.12804 [astro- ph.CO]

G. German, JCAP02, 034, arXiv:2011.12804 [astro- ph.CO]

work page arXiv 2011
[30]

Inflation and the Scale Dependent Spectral Index: Prospects and Strategies

P. Adshead, R. Easther, J. Pritchard, and A. Loeb, JCAP 02, 021, arXiv:1007.3748 [astro-ph.CO]. 7 Appendix Axion de Sitter fluctuations:The dynamics can be understood via stochastic inflation. The equation of motion for the coarse-grained axion angleθ=a/f a in de Sitter space is: f2 a ∂ ∂∆N ⟨θ2⟩= H2 4π2 − 2m2 a 3H2 f2 a ⟨θ2⟩.(15) Form a ≫H, the field is d...

work page internal anchor Pith review Pith/arXiv arXiv