Pith. sign in

REVIEW 3 major objections 3 minor 83 references

A single five-dimensional axion can drive inflation at early times and supply dark matter at late times, with the lightest state stable beyond the age of the Universe.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-01 00:41 UTC pith:VGYOETRB

load-bearing objection The spectral-distortion prediction and the 300-field numerics are the real new content; the headline DM lifetime is asserted rather than derived, and the benchmark has enough internal inconsistencies that it needs major revision before it can be trusted. the 3 major comments →

arxiv 2607.26138 v1 pith:VGYOETRB submitted 2026-07-28 hep-ph

A multi-axion model of inflation and dark matter

Marcos A.G. Garcia , Hansel Gordillo-Ruiz , Raul Henriquez-Ortiz , Saul Ramos-Sanchez This is my paper
classification hep-ph
keywords Kaluza-Klein axionsmulti-axion cosmologydynamical dark matterinflationspectral distortionsmisalignment mechanismreheatingaxion-like particles
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper sets out to show that one five-dimensional axion, compactified on a circle with a brane-localized potential, can be responsible for both cosmic inflation and the observed dark matter. After dimensional reduction the axion becomes a tower of Kaluza-Klein states; the paper argues that the heaviest state acts as the inflaton, while the lighter states behave as a decaying ensemble of dark-matter particles whose lightest member outlives the Universe. For the benchmark parameters the inflationary trajectory is nearly geodesic, so the cosmology is effectively single-field and matches the measured spectral index, tensor-to-scalar ratio, amplitude, and running. The same setup reheats the Universe to about 5 MeV, consistent with nucleosynthesis, and produces percent-level deviations in CMB spectral distortions that future experiments could detect. If correct, inflation and dark matter would share one origin and a concrete observational signature.

Core claim

The central claim is that the Kaluza-Klein tower of a single five-dimensional axion interpolates between inflation and dark matter. In the mass-eigenstate basis the heaviest mode, a_299, carries the inflationary trajectory, with a nearly vanishing turning rate so that the dynamics are effectively single-field and standard slow-roll formulas apply. The predicted observables for the benchmark are n_s = 0.9662, r = 1.4e-8, A_s = 2.105e-9, and alpha_s = -0.0106, consistent with the latest CMB anisotropy data; a small phase shift can also accommodate the larger spectral index favored by newer ACT measurements. After inflation the inflaton decays to photons and reheats the Universe to T_RH about 5

What carries the argument

The central object is the Kaluza-Klein tower generated by compactifying a five-dimensional axion on S1/Z2, with all modes sharing a single decay constant f_a. The mass matrix is diagonalized by a rotation matrix R, yielding mass eigenstates a_i; the hierarchy m_Lambda about 300 M_c selects 300 states, making a_299 the effective inflaton and a_0 the lightest dark-matter candidate. The cascade of decays is controlled by the large ratio Gamma_299/Gamma_0 about 3e26, which allows a reheating temperature of a few MeV and a lightest-state lifetime of 1e25 s at the same time. The geodesic trajectory, quantified by |D_N T| approximately 0, is what licenses the effective single-field treatment that c

Load-bearing premise

The unified story depends on a single five-dimensional axion sector, but the benchmark implements inflation and dark matter with two slightly different potentials by adding a global phase Psi = -3.1 to restore inflationary initial conditions after the misalignment angle theta is fixed to match the dark-matter abundance; if that phase is physical, the two epochs are not described by exactly the same potential.

What would settle it

Recompute the benchmark with one fixed potential, removing the additional global phase Psi, and check whether a single misalignment angle theta can simultaneously reproduce Omega_DM about 0.26 and the approximately 40.6 e-folds required for matching CMB observables; if no common parameter point exists, the unified claim fails. Alternatively, a future spectral-distortion experiment with sensitivity to percent-level deviations that measures a signal matching the standard model exactly would rule out the predicted suppression.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Inflation and dark matter arise from a single higher-dimensional axion sector, so no separate inflaton and dark-matter particle are needed.
  • The dark matter is an ensemble of decaying states; the lightest member is effectively stable on cosmological timescales, so decay searches constrain the model rather than exclude it.
  • The predicted CMB spectral distortion is about 2.5% lower than the standard prediction at frequencies below 100 GHz and above 250 GHz, within the reach of future spectral-distortion experiments.
  • Reheating can proceed to about 5 MeV, just above the BBN bound, and increasing the hierarchy between the cutoff scale and the compactification scale yields higher reheating temperatures while keeping the lightest dark-matter state long-lived.
  • Because all initial misalignments stem from one angle theta, the dark-matter isocurvature fluctuations are correlated with the inflaton fluctuations, a distinguishing feature that the paper leaves for dedicated follow-up.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • A single-potential benchmark without the additional global phase Psi is the real test of unification; the paper's Psi = -3.1 adjustment means the inflation and dark-matter epochs are modeled with slightly different potentials, so the unification claim is not yet fully demonstrated.
  • Stronger gamma-ray limits on decaying dark matter could push the required tau_0 beyond 1e25 s, forcing a larger m_Lambda/M_c hierarchy and shifting the reheating temperature and spectral-distortion predictions.
  • The predicted percent-level spectral-distortion suppression, if confirmed, would distinguish this model from other single-field inflationary scenarios; the paper does not perform that comparison.
  • Computing the linear matter power spectrum from the decaying tower would test whether the multi-component dark-matter structure leaves a cutoff accessible to Lyman-alpha forest surveys.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 3 minor

Summary. The manuscript proposes a five-dimensional axion model compactified on S^1/Z_2, whose Kaluza-Klein tower is claimed to provide a unified description of inflation and dark matter. In the benchmark setup, the heaviest eigenstate drives inflation, reheating proceeds via its decay to photons, the lighter eigenstates form a dynamical dark-matter ensemble, and the lightest state is claimed to be stable on timescales τ_0 ≳ 10^25 s. The paper computes inflationary observables, reheating dynamics, spectral distortions, and dark-matter abundances, claiming consistency with Planck and BBN constraints. A phase parameter Ψ is introduced in Sec. 5.3 to reconcile the dark-matter initial condition with the inflationary initial condition.

Significance. If the central claims were correct, the model would be an economical unification of inflation and dark matter in a single higher-dimensional axion sector, with a concrete spectral-distortion signature. The paper also contains a substantial numerical treatment of a 300-field system and a geodesic-trajectory analysis, which is a useful technique. However, the central stability claim is contradicted by the paper's own decay-rate formula, the stationary point used to define the vacuum is not a solution of the potential for the benchmark parameters, and the inflationary and dark-matter calculations rely on two different potentials. These are load-bearing inconsistencies: they undermine the claimed dark-matter lifetime, the mass eigenstate construction, and the single-Lagrangian unification. The numerical machinery is attached to an internally inconsistent benchmark, so the significance of the claimed results is not established.

major comments (3)
  1. [Sec. 3, eqs. (16)–(17) and (45); Table 2] The quoted dark-matter lifetime is inconsistent with the decay-rate formula used elsewhere. For the benchmark, gaγγ = 2.08×10^-27 GeV^-1 and m0 = 7.86×10^5 GeV. Inserting these into the same formula as eq. (45), Γ0 = gaγγ^2 m0^3/(64π), gives Γ0 ≈ 1.0×10^-38 GeV, i.e. τ0 ≈ 6×10^13 s, not Γ0 ≈ 1.52×10^-49 GeV and τ0 ≈ 10^25 s quoted in eqs. (16)–(17). No suppression mechanism is identified anywhere; the mixing factor for the lightest eigenstate is O(1) in this rank-one model. Consequently the abstract's claim of a lightest state 'stable on timescales far exceeding the age of the Universe' fails for the benchmark. Independently, the ratio Γ_299/Γ0 implied by the Table 2 masses and couplings is roughly (m_299/m_0)^3 × (24.5)^2 ≈ 2×10^15, not the ≈10^26 used in eq. (17).
  2. [Sec. 2, eq. (9)] The claimed minimum ⟨ϕ0_min⟩ does not satisfy the stationarity condition of the potential (8). For κ=1, n=3, Θ=7.30×10^-6, the derivative condition is sin(3x) = 3 sin(x+Θ) with x = ϕ0/fa. Eq. (9) with ℓ=1 gives x = π + 1 ≈ 4.14, for which the left-hand side is ≈ −0.14 and the right-hand side is ≈ −2.52; the equation is not close to being satisfied. The actual minimum of the hilltop potential lies near x ≈ π − 0.016 for these parameters. Since eq. (13) and all subsequent mass eigenvalues and couplings are derived by expanding around this claimed minimum, the benchmark's mass spectrum and decay rates are not based on a stationary point of the model potential.
  3. [Sec. 5.3, eq. (77); Secs. 3.2 and 4] The inflation and dark-matter calculations use different potentials. The inflationary observables in Sec. 3.2 and the benchmark in Table 3 are obtained from the potential in eqs. (8) and (34) with phase Θ = 7.30×10^-6. The spectral-distortion analysis in Sec. 4 explicitly refers to 'the inflationary potential V(a_299) defined in eq. (77)', which is introduced only in Sec. 5.3 and contains an additional phase Ψ = −3.1 in both cosine terms. In Sec. 5.3, the dark-matter abundance is first fixed by choosing θ = 6.99×10^-5 using the minima of the Θ-potential, and then the potential is changed by Ψ to preserve the inflationary initial conditions. If Ψ is a physical parameter, the minimum ⟨ϕ_min⟩ entering eq. (58) must be recomputed with Θ+Ψ, so the quoted dark-matter abundance is not tied to a consistent vacuum. If Ψ is a bookkeeping device, no single Lagrangian produces both the inflation obs
minor comments (3)
  1. [Sec. 3, text after eq. (23)] The text states 'The lightest eigenstate has a mass m0 = 7.86×10^6 GeV', while Table 2 lists m0 = 7.86×10^5 GeV. The factor 10 discrepancy changes Γ0 by 10^3 if used in eq. (45), further worsening the lifetime problem.
  2. [Sec. 3, eq. (17)] The arithmetic in eq. (17) is inconsistent: 1.07×10^-23 GeV divided by 1.52×10^-49 GeV is ≈ 7.0×10^25, not ≈ 1.63×10^26.
  3. [Sec. 4] Eq. (77) is cited as the inflationary potential in Sec. 4 before it is defined in Sec. 5.3, and it is not the potential used for the inflationary observables in Sec. 3. This forward reference obscures the two-potential inconsistency.

Circularity Check

3 steps flagged

DM abundance and lightest-state lifetime are fitted inputs, and the DM epoch is computed with a phase-shifted potential, so the claimed inflation+DM unification is partly constructed rather than derived.

specific steps
  1. fitted input called prediction [Sec. 5.3, after eq. (76)]
    "Our choice of initial conditions is fixed such that the total relic abundance of the ensemble matches this observed value at late times. For this we require that θ = 6.99×10^(-5) in eq. (58)."

    The observed DM abundance is imposed by hand through the misalignment angle θ, and the subsequent statement that the ensemble reproduces the observed relic abundance is a restatement of this input. Ω_DM is normalized to the Planck value, so the abundance consistency check cannot fail and provides no independent confirmation of the model.

  2. other [Sec. 5.3, eq. (77)]
    "This election of course can change the initial conditions that we chose for inflation, for these we change the original potential by a global phase in [eq. (77)]."

    The claimed unification relies on one 5D Lagrangian governing both epochs. Here, after fixing θ from the DM abundance, the paper changes the global phase Ψ of the brane potential for the DM computation, so the same potential is not used for inflation and DM. The successful coexistence of inflation and DM is not derived from a single model; it is manufactured by switching the potential between epochs. The central unification claim is thus equivalent to allowing different potentials for the two epochs.

  3. fitted input called prediction [Sec. 3, eqs. (16)-(17) and Table 2]
    "we adopt the conservative requirement τ0 = 1/Γ0 ≳ 10^25 s ... we set the cutoff scale in eq. (8) to mΛ ≃ 300Mc ... corresponding to a lifetime τ0 = 10^25 s."

    The lightest-state lifetime quoted in the abstract is the benchmark requirement itself: eq. (16) imposes τ0 ≳ 10^25 s, and then the cutoff hierarchy mΛ ≃ 300 Mc is chosen so that eq. (17) outputs τ0 = 10^25 s. The 'result' is an input constraint used to select the parameters, not a parameter-free prediction of the model.

full rationale

Two headline outputs are set by hand: the DM abundance is fixed by choosing θ, and the cutoff scale is chosen so that the lightest lifetime meets the adopted constraint. The abundance and stability claims are therefore parameter-selection checks rather than free predictions. More seriously, the DM computation uses eq. (77), where the global phase is shifted relative to the inflation computation, so the same 5D sector is not carried through both epochs; the unification claim is partly manufactured. I do not count the self-citations ([48], [67], [74], [75]) as load-bearing, and there is no imported uniqueness theorem or ansatz hidden behind a citation. The numerical multifield evolution, the geodesic check, and the spectral-distortion estimate provide real independent content, which prevents the score from being maximal. But because the central unification result depends on fitted θ, a tuned cutoff, and an epoch-dependent phase, the partial-circularity score is 6.

Axiom & Free-Parameter Ledger

10 free parameters · 7 axioms · 0 invented entities

The paper introduces no genuinely new particle species beyond the KK tower of a 5D axion, which is standard. The central claims rest on a large set of tuned parameters (fa, Λ, Θ, κ, n, Mc, mode count, θ, Ψ, and the axion-photon coupling) and on several domain assumptions about the brane potential, misalignment, and instantaneous thermalization. The free-parameter count is high, and two of the central outputs (ΩDM and ns/As) are explicitly fit to data.

free parameters (10)
  • fa = 6.92×10^18 GeV
    Effective four-dimensional decay constant; trans-Planckian; tuned together with Λ and Θ to match As and ns.
  • Λ = 3.04×10^14 GeV
    Scale of the brane-localized potential; tuned to set inflation energy scale and mass matrix.
  • Θ = 7.30×10^-6 (Planck) or 8.21×10^-6 (ACT)
    Phase in the hilltop potential; fitted to the scalar spectral index ns.
  • κ = 1
    Coefficient of the second cosine in the potential; chosen benchmark value.
  • n = 3
    Integer harmonic of the second cosine; chosen benchmark value.
  • Mc = 1.57×10^6 GeV
    Compactification scale; chosen so that mΛ ≈ 300 Mc and to satisfy experimental bounds.
  • Number of modes = M+1 = 300
    Truncated KK tower; cutoff mΛ/Mc = 300 chosen by hand; changing it changes lifetimes and reheating.
  • θ = 6.99×10^-5
    Misalignment angle; explicitly fixed so the total DM abundance matches the observed ΩDM.
  • Ψ = -3.1
    Global phase added to the inflaton potential in eq. (77) to restore the inflationary initial conditions after θ is set for DM.
  • gaγγ / cγ = 2.08×10^-27 GeV^-1 (implies cγ ~ 10^-5)
    Axion-photon coupling; the implied EM anomaly coefficient is tiny and unexplained; effectively tuned to suppress decay rates.
axioms (7)
  • standard math KK decomposition of a massless 5D scalar on S1/Z2 with cosine modes
    Eq. (3) uses the standard expansion; no independent verification is provided but it is a conventional mathematical result.
  • domain assumption SM fields are confined to the brane at y=0 and the axion couples only to photons
    Eq. (4) restricts interactions to the brane and to electromagnetism; this drives the reheating and DM decay phenomenology.
  • domain assumption A non-perturbatively generated brane-localized hilltop potential with scale Λ and phase Θ, absent at T ≫ Λ
    Eqs. (6)-(8) define the potential; Sec. 5.1 assumes it is absent at high temperatures for the misalignment mechanism.
  • domain assumption Misalignment mechanism generates the initial field displacements
    Eqs. (57)-(60) parametrize the initial VEV through θ; this sets all DM abundances and the inflaton initial condition.
  • domain assumption Instantaneous thermalization and simplified piecewise Hubble evolution during reheating and later eras
    Eqs. (22), (64), and (71) assume radiation and matter epochs transition abruptly and that decay products thermalize rapidly.
  • domain assumption Slow-roll single-field formulas are valid because the trajectory is geodesic
    Sec. 3.2 uses the standard slow-roll expressions (29)-(32) after numerically demonstrating |D_N T| ≪ 1 for the benchmark.
  • domain assumption The lightest state must satisfy τ0 ≳ 10^25 s (conservative requirement)
    Eq. (16) is imposed rather than derived; it is used to select parameters and decay rates.

pith-pipeline@v1.3.0-alltime-deepseek · 23367 in / 30089 out tokens · 277741 ms · 2026-08-01T00:41:57.110773+00:00 · methodology

0 comments
read the original abstract

Models with extra dimensions include a tower of massive states that unavoidably contribute to cosmological dynamics. In particular, after compactification, a higher-dimensional axion gives rise to a set of states whose dynamics at early times exhibit the properties of inflation and, at late times, those of dynamical dark matter. Interestingly, natural choices of parameters imply that the inflationary trajectory in field space is approximately geodesic, meaning that inflation is effectively driven by a single inflaton corresponding to the heaviest state. In this scenario, all inflationary observables, including the emerging spectral distortions, are consistent with current observational constraints, and exact numerical solutions of the equations of motion show that reheating can consistently proceed down to temperatures of a few MeV. Furthermore, all states lighter than the inflaton only come to dominate the energy density of the Universe during the matter-dominated era, thus behaving as dark matter at late times. These states undergo a cascade of decays into radiation, while maintaining a sufficient dark-matter abundance to account for the observed relic density. The lightest dark-matter candidate remains stable on timescales far exceeding the age of the Universe, thereby satisfying even the most conservative observational constraints.

Figures

Figures reproduced from arXiv: 2607.26138 by Hansel Gordillo-Ruiz, Marcos A.G. Garcia, Raul Henriquez-Ortiz, Saul Ramos-Sanchez.

Figure 1
Figure 1. Figure 1: Evolution of the fields ϕ i as a function of the number of e-folds N (the superscript is not the power of ϕ). The trajectories were obtained by numerically solving the equations of motion (18) for 300 fields with the parameters shown in [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Evolution of the mass eigenstates a i as a function of the number of e-folds N. The left panel illustrates the evolution of a representative subset of the heavier eigenstates, while the right panel focuses on the lightest eigenstate, which later constitutes the dominant DM component. The evolution is obtained from the numerical solution of the background equations of motion (18) using the benchmark paramet… view at source ↗
Figure 3
Figure 3. Figure 3: Evolution of |DN T| during inflation. The end of inflation occurs at Nend = 72.57. We find |DN T| ≃ 0 throughout the inflationary phase, indicating that the background trajectory follows an (ap￾proximately) geodesic motion. magnitude of the bending is then quantified by |DN T| = p δ IJωIωJ . (28) If |DN T| ≪ 1, the trajectory follows an approximately geodesic path in field space and the dynamics effectivel… view at source ↗
Figure 4
Figure 4. Figure 4: Evolution of the components of the tangent vector [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Evolution of the inflaton energy density [PITH_FULL_IMAGE:figures/full_fig_p017_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Left: total spectral distortion ∆I predicted by our model (red curve), compared with the ΛCDM prediction (dark gray region) and the expected sensitivity of a future PIXIE-like experiment (light gray shaded region). Right: percentage difference between the signal predicted by our model and the ΛCDM prediction. For frequencies outside the range 100 GHz ≲ ν ≲ 250 GHz, the difference with respect to the fiduci… view at source ↗
Figure 7
Figure 7. Figure 7: Evolution of the energy densities ρi associated with the fields a i during the reheating epoch, together with the radiation energy density ρr, indicated by a red dashed line. The figure shows that the energy densities of the DM fields remain subdominant throughout reheating. 5.2 DM during reheating Once the condition 3H(ti,osc) ≃ 2mi is satisfied, the corresponding field a i begins to oscillate around the … view at source ↗
Figure 8
Figure 8. Figure 8: Comparison between the decay source term of the inflaton, Γ [PITH_FULL_IMAGE:figures/full_fig_p023_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Evolution of the individual abundances Ω [PITH_FULL_IMAGE:figures/full_fig_p024_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Evolution of the tower fraction ξ as a function of the scale factor a. The figure shows that ξ ≃ O(0.1) throughout most of the cosmological evolution, indicating a significant departure from the standard single-field DM scenario. and matter-dominated eras at aEQ, as well as the present-day scale factor a0 = 1. At late times, the abundances asymptote to constant values, and their sum reproduces the observe… view at source ↗
Figure 11
Figure 11. Figure 11: Evolution of the effective equation-of-state parameter [PITH_FULL_IMAGE:figures/full_fig_p027_11.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

83 extracted references · 72 linked inside Pith

  1. [1]

    A. H. Guth,The Inflationary Universe: A Possible Solution to the Horizon and Flatness Problems, Phys. Rev. D23(1981), 347–356

    1981

  2. [2]

    A. D. Linde,A New Inflationary Universe Scenario: A Possible Solution of the Horizon, Flatness, Homogeneity, Isotropy and Primordial Monopole Problems, Phys. Lett. B108 (1982), 389–393

    1982

  3. [3]

    A. R. Liddle,An Introduction to cosmological inflation, inICTP Summer School in High- Energy Physics and Cosmology, 1 1999, pp. 260–295

    1999

  4. [4]

    Baumann,Inflation, inTheoretical Advanced Study Institute in Elementary Particle Physics: Physics of the Large and the Small, 2011, pp

    D. Baumann,Inflation, inTheoretical Advanced Study Institute in Elementary Particle Physics: Physics of the Large and the Small, 2011, pp. 523–686

    2011

  5. [5]

    Aghanim et al.,Planck 2018 results

    Planck, N. Aghanim et al.,Planck 2018 results. I. Overview and the cosmological legacy of Planck, Astron. Astrophys.641(2020), A1,arXiv:1807.06205[astro-ph.CO]

    Pith/arXiv arXiv 2018

  6. [6]

    Hu and S

    W. Hu and S. Dodelson,Cosmic Microwave Background Anisotropies, Ann. Rev. Astron. Astrophys.40(2002), 171–216,astro-ph/0110414. 28

    Pith/arXiv arXiv 2002

  7. [7]

    D. H. Lyth and A. Riotto,Particle physics models of inflation and the cosmological density perturbation, Phys. Rept.314(1999), 1–146,hep-ph/9807278

    Pith/arXiv arXiv 1999

  8. [8]

    A. R. Liddle and D. H. Lyth,Cosmological inflation and large scale structure, 2000

    2000

  9. [9]

    Pajer and M

    E. Pajer and M. Peloso,A review of Axion Inflation in the era of Planck, Class. Quant. Grav.30(2013), 214002,arXiv:1305.3557[hep-th]

    Pith/arXiv arXiv 2013

  10. [10]

    J. L. Feng,Dark Matter Candidates from Particle Physics and Methods of Detection, Ann. Rev. Astron. Astrophys.48(2010), 495–545,arXiv:1003.0904[astro-ph.CO]

    Pith/arXiv arXiv 2010

  11. [11]

    Bozorgnia, J

    N. Bozorgnia, J. Bramante, J. M. Cline, D. Curtin, D. McKeen, D. E. Morrissey, A. Ritz, S. Viel, A. C. Vincent, and Y. Zhang,Dark matter candidates and searches, Can. J. Phys. 103(2025), no. 8, 671–703,arXiv:2410.23454[hep-ph]

    Pith/arXiv arXiv 2025

  12. [12]

    R. D. Peccei,The Strong CP problem and axions, Lect. Notes Phys.741(2008), 3–17, hep-ph/0607268

    Pith/arXiv arXiv 2008

  13. [13]

    J. E. Kim and H. P. Nilles,The mu Problem and the Strong CP Problem, Phys. Lett. B138 (1984), 150–154

    1984

  14. [14]

    Randall,Composite axion models and Planck scale physics, Phys

    L. Randall,Composite axion models and Planck scale physics, Phys. Lett. B284(1992), 77–80

    1992

  15. [15]

    L. Covi, J. E. Kim, and L. Roszkowski,Axinos as cold dark matter, Phys. Rev. Lett.82 (1999), 4180–4183,hep-ph/9905212

    Pith/arXiv arXiv 1999

  16. [16]

    K.-Y. Choi, L. Covi, J. E. Kim, and L. Roszkowski,Axino Cold Dark Matter Revisited, JHEP04(2012), 106,arXiv:1108.2282[hep-ph]

    Pith/arXiv arXiv 2012

  17. [17]

    Higaki, K

    T. Higaki, K. S. Jeong, N. Kitajima, and F. Takahashi,Quality of the Peccei-Quinn sym- metry in the Aligned QCD Axion and Cosmological Implications, JHEP06(2016), 150, arXiv:1603.02090[hep-ph]

    Pith/arXiv arXiv 2016

  18. [18]

    Calibbi, F

    L. Calibbi, F. Goertz, D. Redigolo, R. Ziegler, and J. Zupan,Minimal axion model from flavor, Phys. Rev. D95(2017), no. 9, 095009,arXiv:1612.08040[hep-ph]

    Pith/arXiv arXiv 2017

  19. [19]

    Di Luzio, M

    L. Di Luzio, M. Giannotti, E. Nardi, and L. Visinelli,The landscape of QCD axion models, Phys. Rept.870(2020), 1–117,arXiv:2003.01100[hep-ph]

    Pith/arXiv arXiv 2020

  20. [20]

    Svrˇ cek and E

    P. Svrˇ cek and E. Witten,Axions In String Theory, JHEP06(2006), 051,hep-th/0605206

    Pith/arXiv arXiv 2006

  21. [21]

    Arvanitaki, S

    A. Arvanitaki, S. Dimopoulos, S. Dubovsky, N. Kaloper, and J. March-Russell,String Axi- verse, Phys. Rev. D81(2010), 123530,arXiv:0905.4720[hep-th]

    Pith/arXiv arXiv 2010

  22. [22]

    K.-S. Choi, H. P. Nilles, S. Ramos-S´ anchez, and P. K. S. Vaudrevange,Accions, Phys. Lett. B675(2009), 381–386,arXiv:0902.3070[hep-th]. 29

    Pith/arXiv arXiv 2009

  23. [23]

    B. S. Acharya, K. Bobkov, and P. Kumar,An M Theory Solution to the Strong CP Problem and Constraints on the Axiverse, JHEP11(2010), 105,arXiv:1004.5138[hep-th]

    Pith/arXiv arXiv 2010

  24. [24]

    Cicoli, M

    M. Cicoli, M. Goodsell, and A. Ringwald,The type IIB string axiverse and its low-energy phenomenology, JHEP10(2012), 146,arXiv:1206.0819[hep-th]

    Pith/arXiv arXiv 2012

  25. [25]

    Halverson, C

    J. Halverson, C. Long, and P. Nath,Ultralight axion in supersymmetry and strings and cosmology at small scales, Phys. Rev. D96(2017), no. 5, 056025,arXiv:1703.07779[hep- ph]

    Pith/arXiv arXiv 2017

  26. [26]

    Petrossian-Byrne and G

    R. Petrossian-Byrne and G. Villadoro,Open string axiverse, JHEP07(2025), 049, arXiv:2503.16387[hep-ph]

    Pith/arXiv arXiv 2025

  27. [27]

    J. M. Leedom, M. Putti, and A. Westphal,Towards a heterotic axiverse, JHEP05(2026), 271,arXiv:2509.03578[hep-th]

    Pith/arXiv arXiv 2026

  28. [28]

    Gordillo-Ruiz, M

    H. Gordillo-Ruiz, M. Hern´ andez-Segura, I. Portillo-Castillo, S. Ramos-S´ anchez, and I. Zavala,Rolling with modular symmetry: quintessence and de Sitter in heterotic orbifolds, JHEP04(2026), 211,arXiv:2509.22781[hep-th]

    arXiv 2026

  29. [29]

    Freese, J

    K. Freese, J. A. Frieman, and A. V. Olinto,Natural inflation with pseudo - Nambu-Goldstone bosons, Phys. Rev. Lett.65(1990), 3233–3236

    1990

  30. [30]

    F. C. Adams, J. R. Bond, K. Freese, J. A. Frieman, and A. V. Olinto,Natural inflation: Particle physics models, power law spectra for large scale structure, and constraints from COBE, Phys. Rev. D47(1993), 426–455,hep-ph/9207245

    Pith/arXiv arXiv 1993

  31. [31]

    J. E. Kim, H. P. Nilles, and M. Peloso,Completing natural inflation, JCAP01(2005), 005, hep-ph/0409138

    Pith/arXiv arXiv 2005

  32. [32]

    Dimopoulos, S

    S. Dimopoulos, S. Kachru, J. McGreevy, and J. G. Wacker,N-flation, JCAP08(2008), 003, hep-th/0507205

    Pith/arXiv arXiv 2008

  33. [33]

    Czerny and F

    M. Czerny and F. Takahashi,Multi-Natural Inflation, Phys. Lett. B733(2014), 241–246, arXiv:1401.5212[hep-ph]

    Pith/arXiv arXiv 2014

  34. [34]

    Daido, F

    R. Daido, F. Takahashi, and W. Yin,The ALP miracle revisited, JHEP02(2018), 104, arXiv:1710.11107[hep-ph]

    Pith/arXiv arXiv 2018

  35. [35]

    Y. Abe, T. Higaki, F. Kaneko, T. Kobayashi, and H. Otsuka,Moduli inflation from modular flavor symmetries, JHEP06(2023), 187,arXiv:2303.02947[hep-ph]

    Pith/arXiv arXiv 2023

  36. [36]

    H. M. Lee, A. G. Menkara, M.-J. Seong, and J.-H. Song,Peccei-Quinn inflation at the pole and axion kinetic misalignment, JHEP05(2024), 295,arXiv:2310.17710[hep-ph]. 30

    Pith/arXiv arXiv 2024

  37. [37]

    K. V. Berghaus, M. Forslund, and M. V. Guevarra,Warm inflation with a heavy QCD axion, JCAP10(2024), 103,arXiv:2402.13535[hep-ph]

    Pith/arXiv arXiv 2024

  38. [38]

    Huang, R

    R. Huang, R. Zheng, and Q. Yang,Dual axion-like field inflation, Commun. Theor. Phys. 77(2025), no. 12, 125403,arXiv:2508.17224[astro-ph.CO]

    Pith/arXiv arXiv 2025

  39. [39]

    Baumann and L

    D. Baumann and L. McAllister,Advances in Inflation in String Theory, Ann. Rev. Nucl. Part. Sci.59(2009), 67–94,arXiv:0901.0265[hep-th]

    Pith/arXiv arXiv 2009

  40. [40]

    Cicoli and F

    M. Cicoli and F. Quevedo,String moduli inflation: An overview, Class. Quant. Grav.28 (2011), 204001,arXiv:1108.2659[hep-th]

    Pith/arXiv arXiv 2011

  41. [41]

    Marchesano, G

    F. Marchesano, G. Shiu, and A. M. Uranga,F-term Axion Monodromy Inflation, JHEP09 (2014), 184,arXiv:1404.3040[hep-th]

    Pith/arXiv arXiv 2014

  42. [42]

    Blumenhagen, D

    R. Blumenhagen, D. Herschmann, and E. Plauschinn,The Challenge of Realizing F-term Axion Monodromy Inflation in String Theory, JHEP01(2015), 007,arXiv:1409.7075[hep- th]

    Pith/arXiv arXiv 2015

  43. [43]

    E. G. M. Ferreira,Ultra-light dark matter, Astron. Astrophys. Rev.29(2021), no. 1, 7, arXiv:2005.03254[astro-ph.CO]

    Pith/arXiv arXiv 2021

  44. [44]

    Arcadi, M

    G. Arcadi, M. Dutra, P. Ghosh, M. Lindner, Y. Mambrini, M. Pierre, S. Profumo, and F. S. Queiroz,The waning of the WIMP? A review of models, searches, and constraints, Eur. Phys. J. C78(2018), no. 3, 203,arXiv:1703.07364[hep-ph]

    Pith/arXiv arXiv 2018

  45. [45]

    Arcadi, D

    G. Arcadi, D. Cabo-Almeida, M. Dutra, P. Ghosh, M. Lindner, Y. Mambrini, J. P. Neto, M. Pierre, S. Profumo, and F. S. Queiroz,The Waning of the WIMP: Endgame?, Eur. Phys. J. C85(2025), no. 2, 152,arXiv:2403.15860[hep-ph]

    Pith/arXiv arXiv 2025

  46. [46]

    K. R. Dienes and B. Thomas,Dynamical Dark Matter: I. Theoretical Overview, Phys. Rev. D85(2012), 083523,arXiv:1106.4546[hep-ph]

    Pith/arXiv arXiv 2012

  47. [47]

    K. R. Dienes and B. Thomas,Dynamical Dark Matter: II. An Explicit Model, Phys. Rev. D 85(2012), 083524,arXiv:1107.0721[hep-ph]

    Pith/arXiv arXiv 2012

  48. [48]

    Gordillo-Ruiz, C

    H. Gordillo-Ruiz, C. Morales-Mena, and S. Ramos-S´ anchez,The simplest case of dark in- flaxion, (2023),arXiv:2310.06893[hep-ph]

    Pith/arXiv arXiv 2023

  49. [49]

    Particle Data Group, R. L. Workman and Others,Review of Particle Physics, PTEP2022 (2022), 083C01

    2022

  50. [50]

    K. R. Dienes, E. Dudas, and T. Gherghetta,Invisible axions and large radius compactifica- tions, Phys. Rev. D62(2000), 105023,hep-ph/9912455. 31

    Pith/arXiv arXiv 2000

  51. [51]

    D. J. Kapner, T. S. Cook, E. G. Adelberger, J. H. Gundlach, B. R. Heckel, C. D. Hoyle, and H. E. Swanson,Tests of the gravitational inverse-square law below the dark-energy length scale, Phys. Rev. Lett.98(2007), 021101,hep-ph/0611184

    Pith/arXiv arXiv 2007

  52. [52]

    Hannestad and G

    S. Hannestad and G. G. Raffelt,Supernova and neutron star limits on large extra dimensions reexamined, Phys. Rev. D67(2003), 125008,hep-ph/0304029, [Erratum: Phys.Rev.D 69, 029901 (2004)]

    Pith/arXiv arXiv 2003

  53. [53]

    B. D. Fields, K. A. Olive, T.-H. Yeh, and C. Young,Big-Bang Nucleosynthesis after Planck, JCAP03(2020), 010,arXiv:1912.01132[astro-ph.CO], [Erratum: JCAP 11, E02 (2020)]

    Pith/arXiv arXiv 2020

  54. [54]

    Hasegawa, N

    T. Hasegawa, N. Hiroshima, K. Kohri, R. S. L. Hansen, T. Tram, and S. Hannestad, MeV-scale reheating temperature and thermalization of oscillating neutrinos by radiative and hadronic decays of massive particles, JCAP12(2019), 012,arXiv:1908.10189[hep-ph]

    Pith/arXiv arXiv 2019

  55. [55]

    Cohen, K

    T. Cohen, K. Murase, N. L. Rodd, B. R. Safdi, and Y. Soreq,γ-ray Constraints on Decaying Dark Matter and Implications for IceCube, Phys. Rev. Lett.119(2017), no. 2, 021102, arXiv:1612.05638[hep-ph]

    Pith/arXiv arXiv 2017

  56. [56]

    Blanco and D

    C. Blanco and D. Hooper,Constraints on Decaying Dark Matter from the Isotropic Gamma- Ray Background, JCAP03(2019), 019,arXiv:1811.05988[astro-ph.HE]

    Pith/arXiv arXiv 2019

  57. [57]

    Burgess and D

    C. Burgess and D. Roest,Inflation by Alignment, JCAP06(2015), 012,arXiv:1412.1614 [hep-th]

    Pith/arXiv arXiv 2015

  58. [58]

    J. E. Kim, S. Nam, and Y. K. Semetzidis,Fate of global symmetries in the Universe: QCD axion, quintessential axion and trans-Planckian inflaton decay-constant, Int. J. Mod. Phys. A33(2018), no. 03, 1830002,arXiv:1712.08648[hep-ph]

    Pith/arXiv arXiv 2018

  59. [59]

    Aghanim et al.,Planck 2018 results

    Planck, N. Aghanim et al.,Planck 2018 results. VI. Cosmological parameters, Astron. Astro- phys.641(2020), A6,arXiv:1807.06209[astro-ph.CO], [Erratum: Astron.Astrophys. 652, C4 (2021)]

    Pith/arXiv arXiv 2018

  60. [60]

    A. R. Liddle and S. M. Leach,How long before the end of inflation were observable pertur- bations produced?, Phys. Rev. D68(2003), 103503,astro-ph/0305263

    Pith/arXiv arXiv 2003

  61. [61]

    Martin and C

    J. Martin and C. Ringeval,First CMB Constraints on the Inflationary Reheating Tempera- ture, Phys. Rev. D82(2010), 023511,arXiv:1004.5525[astro-ph.CO]

    Pith/arXiv arXiv 2010

  62. [62]

    D. J. Fixsen,The temperature of the cosmic microwave background, The Astrophysical Jour- nal707(2009), no. 2, 916–920,http://dx.doi.org/10.1088/0004-637X/707/2/916

    doi:10.1088/0004-637x/707/2/916 2009

  63. [63]

    Gordon, D

    C. Gordon, D. Wands, B. A. Bassett, and R. Maartens,Adiabatic and entropy perturbations from inflation, Phys. Rev. D63(2000), 023506,astro-ph/0009131. 32

    Pith/arXiv arXiv 2000

  64. [64]

    C. M. Peterson and M. Tegmark,Testing multifield inflation: A geometric approach, Phys. Rev. D87(2013), no. 10, 103507,arXiv:1111.0927[astro-ph.CO]

    Pith/arXiv arXiv 2013

  65. [65]

    Lalak, D

    Z. Lalak, D. Langlois, S. Pokorski, and K. Turzynski,Curvature and isocurvature perturba- tions in two-field inflation, JCAP07(2007), 014,arXiv:0704.0212[hep-th]

    Pith/arXiv arXiv 2007

  66. [66]

    Louis et al.,The Atacama Cosmology Telescope: DR6 power spectra, likelihoods andΛCDM parameters, JCAP11(2025), 062,arXiv:2503.14452 [astro-ph.CO]

    Atacama Cosmology Telescope, T. Louis et al.,The Atacama Cosmology Telescope: DR6 power spectra, likelihoods andΛCDM parameters, JCAP11(2025), 062,arXiv:2503.14452 [astro-ph.CO]

    Pith/arXiv arXiv 2025

  67. [67]

    M. A. G. Garc ´ ıa, K. Kaneta, Y. Mambrini, and K. A. Olive,Inflaton Oscillations and Post- Inflationary Reheating, JCAP04(2021), 012,arXiv:2012.10756[hep-ph]

    Pith/arXiv arXiv 2021

  68. [68]

    Clery, M

    S. Clery, M. A. G. Garc ´ ıa, Y. Mambrini, and K. A. Olive,Bare mass effects on the reheating process after inflation, Phys. Rev. D109(2024), no. 10, 103540,arXiv:2402.16958[hep-ph]

    Pith/arXiv arXiv 2024

  69. [69]

    J. E. Kim and G. Carosi,Axions and the Strong CP Problem, Rev. Mod. Phys.82(2010), 557–602,arXiv:0807.3125[hep-ph], [Erratum: Rev.Mod.Phys. 91, 049902 (2019)]

    Pith/arXiv arXiv 2010

  70. [70]

    W. Hu, D. Scott, and J. Silk,Power spectrum constraints from spectral distortions in the cosmic microwave background, Astrophys. J. Lett.430(1994), L5–L8,astro-ph/9402045

    Pith/arXiv arXiv 1994

  71. [71]

    Chluba et al.,New horizons in cosmology with spectral distortions of the cosmic microwave background, Exper

    J. Chluba et al.,New horizons in cosmology with spectral distortions of the cosmic microwave background, Exper. Astron.51(2021), no. 3, 1515–1554,arXiv:1909.01593[astro-ph.CO]

    Pith/arXiv arXiv 2021

  72. [72]

    Lucca, N

    M. Lucca, N. Sch¨ oneberg, D. C. Hooper, J. Lesgourgues, and J. Chluba,The synergy between CMB spectral distortions and anisotropies, JCAP02(2020), 026,arXiv:1910.04619[astro- ph.CO]

    Pith/arXiv arXiv 2020

  73. [73]

    Chluba and D

    J. Chluba and D. Grin,CMB spectral distortions from small-scale isocurvature fluctuations, Mon. Not. Roy. Astron. Soc.434(2013), 1619–1635,arXiv:1304.4596[astro-ph.CO]

    Pith/arXiv arXiv 2013

  74. [74]

    Henr ´ ıquez-Ortiz, J

    R. Henr ´ ıquez-Ortiz, J. Mastache, and S. Ramos-S´ anchez,Spectral distortions from axion monodromy inflation, JCAP08(2022), no. 08, 054,arXiv:2206.07719[astro-ph.CO]

    Pith/arXiv arXiv 2022

  75. [75]

    A. Baur, M. A. G. Garc ´ ıa, R. Henr ´ ıquez-Ortiz, M. Hern´ andez-Neri, and S. Ramos-S´ anchez, Spectral distortions from promising single and multifield inflationary models, JCAP04 (2024), 090,arXiv:2310.13071[astro-ph.CO]

    Pith/arXiv arXiv 2024

  76. [76]

    M. Viel, G. D. Becker, J. S. Bolton, and M. G. Haehnelt,Warm dark matter as a solution to the small scale crisis: New constraints from high redshift Lyman-αforest data, Phys. Rev. D88(2013), 043502,arXiv:1306.2314[astro-ph.CO]

    Pith/arXiv arXiv 2013

  77. [77]

    V. K. Narayanan, D. N. Spergel, R. Dave, and C.-P. Ma,Constraints on the mass of warm dark matter particles and the shape of the linear power spectrum from the Lyαforest, Astro- phys. J. Lett.543(2000), L103–L106,astro-ph/0005095. 33

    Pith/arXiv arXiv 2000

  78. [78]

    M. Viel, J. Lesgourgues, M. G. Haehnelt, S. Matarrese, and A. Riotto,Constraining warm dark matter candidates including sterile neutrinos and light gravitinos with WMAP and the Lyman-alpha forest, Phys. Rev. D71(2005), 063534,astro-ph/0501562

    Pith/arXiv arXiv 2005

  79. [79]

    J. Baur, N. Palanque-Delabrouille, C. Y` eche, C. Magneville, and M. Viel,Lyman-alpha Forests cool Warm Dark Matter, JCAP08(2016), 012,arXiv:1512.01981[astro-ph.CO]

    Pith/arXiv arXiv 2016

  80. [80]

    Irˇ siˇ c et al.,New Constraints on the free-streaming of warm dark matter from inter- mediate and small scale Lyman-αforest data, Phys

    V. Irˇ siˇ c et al.,New Constraints on the free-streaming of warm dark matter from inter- mediate and small scale Lyman-αforest data, Phys. Rev. D96(2017), no. 2, 023522, arXiv:1702.01764[astro-ph.CO]

    Pith/arXiv arXiv 2017

Showing first 80 references.