Pith. sign in

REVIEW 2 major objections 5 minor 80 references

This paper constructs a technically natural axion dark energy model in which thermal freeze-out of WIMP dark matter imprints a nonuniform relic distribution, generating an unsuppressed finite-density potential that holds the axion at its in

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 →

T0 review · deepseek-v4-flash

2026-08-03 00:36 UTC pith:VBEMSKHV

load-bearing objection A clean, honest construction of apparent phantom crossing from axion–WIMP freeze-out memory, with one load-bearing assumption—inter-sector chemical equilibrium through freeze-out—deferred to future portal models; worth refereeing. the 2 major comments →

arxiv 2607.28721 v1 pith:VBEMSKHV submitted 2026-07-30 hep-ph astro-ph.CO

Natural Phantom Crossing from Axion-WIMP Interactions

classification hep-ph astro-ph.CO PACS 95.35.+d98.80.-k
keywords axion dark energyWIMP dark matterZ_N symmetryphantom crossingfreeze-outfinite-density potentialdark sector interactionDESI
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.

The paper attempts to show that an ultralight axion dark energy field can interact with WIMP dark matter without radiative destabilization, and that this interaction naturally yields an apparent late-time phantom crossing. It introduces N fermion WIMP copies under a cyclic Z_N symmetry, so that the leading quantum correction to the axion potential is exponentially suppressed. The crucial step is that thermal freeze-out imprints a nonuniform abundance distribution across the N sectors, breaking the symmetry spontaneously and producing an unsuppressed finite-density potential that holds the axion at its initial value. As the WIMP density redshifts, the axion rolls to its vacuum, increasing the effective WIMP mass and transferring energy from dark energy to dark matter; an observer assuming separate conservation infers w < -1. If true, this would provide a technically natural dark-energy model consistent with DESI hints without ghosts or NEC violation.

Core claim

On the paper's own terms, the central discovery is that the cosmological relic state can spontaneously break the exact Z_N symmetry of the Lagrangian, storing the initial axion value in the WIMP sector. Freeze-out with slow inter-sector conversion produces Boltzmann-suppressed relic fractions that retain a memory of the initial axion angle, generating an unsuppressed finite-density potential V_fd ≈ σ_fd n_χ [1 - cos(Θ - Θ_i)]. This potential traps the axion at Θ_i at early times; once n_χ drops below Λ_D/σ_fd, the vacuum potential releases it, and the relic-weighted WIMP mass increases, transferring energy from DE to DM. The result is an effective equation of state that crosses below -1, com

What carries the argument

The load-bearing mechanism is the cyclic Z_N symmetry acting on N Dirac fermion WIMPs, χ_k → χ_{k+1}, Φ → e^{2πi/N} Φ, combined with a mass matrix M + y_χ e^{2πik/N} Φ. The root-of-unity sums project the Coleman-Weinberg potential onto the Nth harmonic, suppressing radiative corrections by ε^N. During freeze-out, the axion is frozen at Θ_i; the species-dependent masses m_k(Θ_i) create unequal equilibrium abundances, and when inter-sector conversion decouples at x_d, the fractions η_k ≈ e^{-q_d cos(...)}/... become nonuniform, producing the finite-density potential V_fd with amplitude σ_fd = εM I_1(q_d)/I_0(q_d). This potential depends only on Θ - Θ_i and stores the initial field value.

Load-bearing premise

The whole late-time energy transfer relies on inter-sector conversion processes remaining efficient until at least WIMP annihilation freeze-out (x_d ≳ x_f = 25); if conversions decouple earlier, the relic fractions approach uniform, the finite-density potential becomes exponentially suppressed, and the phantom crossing disappears.

What would settle it

Solve the full coupled Boltzmann system for the N sectors with a concrete universal mediator and compute T_d for the benchmark parameters; if T_d > T_f (conversions decouple before annihilation freeze-out), the nonuniformity parameter q_d is reduced, σ_fd drops below the benchmark value, and the phantom crossing vanishes. Observationally, a future growth-rate measurement at z ≈ 0.3-0.5 that shows no percent-level suppression of fσ8 would also pressure the model, since the benchmark predicts such a suppression.

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

If this is right

  • The relic distribution retains a memory of the initial axion angle, and the late-time axion roll transfers energy from DE to DM, producing an effective w crossing below -1.
  • The model is technically natural: the Z_N symmetry keeps radiative corrections to the axion potential exponentially small, so no fine-tuning of the DE potential is needed.
  • The inferred phantom behavior arises within a canonical scalar and respects the null-energy condition; an observer assuming separately conserved components misinterprets the energy exchange.
  • Benchmark cosmological evolution gives a DESI-like phantom crossing at z_c ≈ 0.40 and a percent-level suppression of structure growth, testable with growth-rate data.
  • WIMP multiplicity changes DM signals: direct detection is roughly unchanged, while indirect detection can be enhanced by up to N_f relative to the single-species case.

Where Pith is reading between the lines

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

  • If the freeze-out imprint idea generalizes, any ultralight scalar with a discrete symmetry and multi-component thermal relics could use the relic distribution as a natural 'initial condition memory' for late-time scalar dynamics; the phantom crossing epoch would then encode the decoupling temperature of inter-sector conversions.
  • The assumption that inter-sector conversions stay in equilibrium through freeze-out (x_d ≳ x_f) is not realized in a concrete mediator model here; a full Boltzmann solution with a specific portal is the natural next step and could either confirm or suppress the effect.
  • The same mechanism predicts a correlation between the present-day WIMP mass modulation (σ_fd/M) and the timing of the phantom crossing; a measurement of w(z) together with a direct-detection rate characterized by a mass-varying WIMP could test the model.
  • Since the finite-density potential holds the axion at Θ_i, the model also predicts a non-trivial dependence of the DM abundance distribution on the initial axion angle, which in a landscape scenario would make the phantom crossing epoch vary across Hubble patches.

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

2 major / 5 minor

Summary. The paper proposes a technically natural axion–WIMP model in which a cyclic Z_N symmetry suppresses the one-loop Coleman–Weinberg potential to O(ε^N), while the cosmological freeze-out of N WIMP species with axion-dependent masses creates a nonuniform relic distribution that remembers the initial axion value. This asymmetric relic distribution generates a finite-density potential σ_fd n_χ[1−cos(Θ−Θ_i)] that traps the axion at early times; as n_χ dilutes, the axion rolls toward the vacuum minimum, the relic-weighted WIMP mass increases, and an observer assuming separately conserved components infers an effective DE equation of state crossing below −1. The paper derives the freeze-out analysis, finite-temperature and finite-density potentials, presents an illustrative benchmark with a DESI-like phantom crossing at z_c≈0.40 and percent-level suppression of fσ8, and discusses implications for direct and indirect WIMP searches.

Significance. If the mechanism works as presented, this is a notable proof of principle: apparent phantom dark energy from a canonical scalar with no ghosts or NEC violation, with a technically natural protection of the ultralight axion potential. The algebraic derivations in the appendices are careful and internally consistent, and the model produces concrete phenomenological consequences (nonuniform relic fractions, modified effective WIMP multiplicity, enhanced indirect-detection signals). The principal limitation is that the central dynamical premise — inter-sector chemical equilibrium through annihilation freeze-out, x_d ≳ x_f — is explicitly announced but not realized in a concrete portal model. Since the phantom-crossing epoch and its very existence depend on σ_fd, which is computed under that assumption, the headline result is conditional. The paper is transparent about this gap, which is a credit; nevertheless, the gap is load-bearing and should be addressed before the central claim can be regarded as established.

major comments (2)
  1. [Sec. III, Eq. (23), Eq. (38)] The central dynamical premise is the ordering x_d ≳ x_f, stated as the 'principal dynamical assumption' after Eq. (25). Equation (23) for the frozen fractions and Eq. (38) for σ_fd both assume internal chemical equilibrium until annihilation freeze-out. If conversions decouple earlier (x_d < x_f), the relic fractions are the solution of the coupled Boltzmann system (16), not the equilibrium fractions at T_d; they could be substantially more uniform, suppressing σ_fd and possibly erasing the finite-density potential that traps the field and generates the phantom crossing in Sec. V. The manuscript explicitly postpones realization of x_d ≳ x_f to future work. This is a genuine gap in the central claim: without a concrete portal model, or at least a parametric estimate of ⟨σv⟩_{k→j} versus ⟨σv⟩_{k→SM} that guarantees x_d ≳ x_f, the cosmological result is conditional. Please fill this gap or
  2. [Sec. V, Eq. (48), Eq. (31)] The benchmark (48) lists V_0, Λ_D, ε, f, and φ_i but not N. Yet N controls the suppression of the CW potential relative to Λ_D (Eq. (31)) and the width of the relic distribution (Eq. (54)). For ε=0.204, the leading CW harmonic is suppressed by ε^N/N^2; requiring Λ_N ≪ Λ_D with M∼100 GeV forces N≈70–80. With x_f=25, q_f=ε(x_f−3/2)≈4.8, and N_f≈N/√(π q_f)≈15–20, so the per-species annihilation cross section must be enhanced by roughly this factor (Eq. (53)). The benchmark should state N, verify the CW hierarchy, and confirm that the resulting N_f is compatible with the assumed universal thermalization portal. Without this, the benchmark is not a fully specified point in the model parameter space.
minor comments (5)
  1. [Sec. V, Eq. (40)] The phrase 'phantom crossing is obtained for A(ϕ)<1, that is, when the effective WIMP mass increases at late times' is easy to misread. Since A is normalized to the present mass, A<1 at earlier times is precisely the signature of a mass that increases toward the present. Please rephrase to avoid the apparent contradiction.
  2. [Sec. IV B, Eqs. (C2)-(C6)] The approximation leading to Eq. (38) neglects O(ε^2 x_d) terms. For the benchmark ε=0.204 and x_d=25, ε^2 x_d≈1.0, so these corrections could shift σ_fd at the tens-of-percent level. The authors should either use the unexpanded expression or explicitly quantify the error for the benchmark.
  3. [Introduction, reference [24]] The citation '[24? –27]' contains a stray question mark and should be corrected.
  4. [Sec. VI, Eq. (58)] The statement that the total direct-detection rate is unchanged for equal per-sector scattering cross sections is correct, but the subsequent discussion of enhancement assumes that annihilation and scattering cross sections share the same coupling dependence. This model-dependence should be flagged explicitly in the text.
  5. [Sec. V, Figures 3–6] The numerical solutions are presented without enough detail for exact reproduction (e.g., the precise form of the coupled equations solved, the treatment of radiation and neutrinos, and the implementation of Eq. (40)). For an illustrative benchmark, this is acceptable, but a brief statement of the numerical procedure or a reproducibility note would strengthen the paper.

Circularity Check

0 steps flagged

No significant circularity: the central mechanism is derived from the stated freeze-out equilibrium assumption and is not fitted to the target w(z) nor imported via load-bearing self-citation.

full rationale

The derivation is self-contained. The Z_N suppression of the Coleman-Weinberg potential (Eq. 30) follows from standard root-of-unity projection, with external citations [48,49] to Hook and Brzeminski et al., not to the authors' own work. The relic fractions (Eq. 23) are obtained from the explicit assumption that inter-sector conversions maintain internal chemical equilibrium until T_d; this is an openly stated modeling assumption (after Eq. 25), not a quantity fitted to dark-energy data. The finite-density potential V_fd (Eqs. 36-38) is then computed from those fractions, and its minimum at Θ=Θ_i is a mathematical consequence of the equilibrium distribution, not a separately imposed condition. The phantom-crossing condition (Eqs. 40-42) is a known, cited mapping from DM-mass variation to apparent w_eff (Das et al. [28]); the sign of the energy transfer is fixed by σ_fd>0 and the location of the vacuum, not chosen to match DESI. The benchmark in Eq. (48) is explicitly illustrative and is matched only to the Planck angular acoustic scale, not to the reconstructed w(z), so there is no fitted-input-called-prediction reduction. Reference [37], by one of the present authors, appears only in a broad citation list for interacting dark-sector models and is not load-bearing for any step of the derivation. The principal dynamical assumption x_d≳x_f is flagged by the authors as postponed to future work; this is a robustness/correctness concern, not circularity.

Axiom & Free-Parameter Ledger

8 free parameters · 5 axioms · 3 invented entities

The central claims rest on: (i) the imposed Z_N cyclic symmetry, (ii) the freeze-out assumption x_d ≳ x_f that guarantees a nonuniform relic distribution, (iii) the assumption that the dark confining sector is unpopulated, and (iv) standard statistical mechanics and CW/Boltzmann formulas. Free parameters are the benchmark cosmological/particle numbers (V_0, Λ_D, ε, f, φ_i, N, M, x_d); none are fitted to w(z) data, but they are chosen to produce the DESI-like outcome.

free parameters (8)
  • V_0 = 0.784 H_0^2 M_Pl^2 (benchmark)
    Field-independent vacuum energy offset; set by hand to match the observed late-time DE density; the paper explicitly delegates its origin to the cosmological constant problem (Sec. IV.B, Eq. (39)).
  • Λ_D = 0.784 H_0^2 M_Pl^2 (benchmark)
    Dark topological susceptibility; sets the axion mass m_ϕ ~ H0 and the DE scale; chosen so the effective potential reproduces the observed expansion and θ_* (Sec. V.A).
  • ε = m_χ/M = 0.204 (benchmark)
    Ratio controlling axion modulation of WIMP masses; sets the finite-density potential amplitude and the required N; chosen to give a DESI-like crossing.
  • f = 0.327 M_Pl (benchmark)
    Axion decay constant; chosen so m_ϕ ~ H0 and to match the Planck angular acoustic scale θ_*.
  • φ_i = 0.784π f (benchmark)
    Initial axion displacement from inflation; the phantom crossing requires a nonzero displacement from the vacuum; this value places the field near the maximum of the confining potential, delaying the roll to z_c ~ 0.4.
  • N = not stated for benchmark; N ~ 25–90 for ε ~ 0.01–0.3
    Number of WIMP sectors; must be large enough that the CW potential Λ_N ≪ Λ_D; affects the multiplicity factors N_f, N_d.
  • M = not specified (≈100 GeV used in estimates)
    WIMP mass scale; enters σ_fd, the relic abundance, and DM-search cross sections.
  • x_d (conversion decoupling) = x_d = x_f = 25 (benchmark)
    Temperature of inter-sector conversion decoupling in units of m/T; sets q_d = ε(x_d − 3/2) and hence the nonuniformity of the relic distribution; treated as an assumption, not derived from a portal model (Sec. III).
axioms (5)
  • ad hoc to paper The Z_N cyclic symmetry is an exact symmetry of the perturbative WIMP Lagrangian (Eq. (6)).
    Central ingredient to suppress the CW potential to the N-th harmonic; introduced for this model.
  • domain assumption Inter-sector conversions remain in internal chemical equilibrium until x_d ≳ x_f = 25.
    Principal dynamical assumption of the freeze-out analysis (Sec. III); no concrete portal model realizing it is provided.
  • domain assumption The dark gauge sector is never appreciably populated after inflation and its temperature remains below the confinement scale.
    Stated at the start of Sec. IV.A; needed so dark glueballs contribute negligibly and the confining potential Eq. (27) is the vacuum potential.
  • domain assumption WIMPs are thermalized with the SM bath through a Z_N-symmetric portal and remain in kinetic equilibrium with the SM plasma until T_d.
    Assumed in Sec. III to justify Eqs. (16)–(23); the portal is not specified.
  • domain assumption The axion field is frozen at its inflationary value Θ_i through freeze-out.
    Stated in Sec. III ('as shown in Appendix B') and Sec. V.A; follows from m_ϕ ~ H0 and Hubble friction, but assumes no early-time finite-density kick displaces it.
invented entities (3)
  • N fermion species χ_k with cyclic Z_N symmetry no independent evidence
    purpose: Protect the ultralight axion from radiative corrections and produce a nonuniform relic distribution that stores the initial axion value and generates the finite-density axion potential.
    No independent observable signature; the WIMPs behave as one effective species with rescaled annihilation cross sections; the Z_N multiplicity is not directly detectable.
  • Dark confining gauge group GD = SU(N_c) and its glueballs no independent evidence
    purpose: Generate the axion vacuum potential Λ_D(1 − cos Θ) that sources dark energy.
    The dark sector is assumed unpopulated, so its glueballs are invisible; no external handle.
  • UV completion scalar Σ with U(1)_X no independent evidence
    purpose: Make ε ~ O(1) natural by generating M and m_χ from a common scale (Eqs. (13)–(15)).
    Illustrative UV completion; not needed for the low-energy mechanism and not observable.

pith-pipeline@v1.3.0-alltime-deepseek · 23060 in / 26659 out tokens · 289437 ms · 2026-08-03T00:36:29.911961+00:00 · methodology

0 comments
read the original abstract

We construct a technically natural model in which thermal dark matter (DM) interacts with axion dark energy (DE) and produces an apparent late-time crossing of the phantom divide. A direct axion coupling to weakly-interacting massive particles (WIMPs) would ordinarily radiatively destabilize the ultralight axion potential. We avoid this issue through $N$ fermion species related by a cyclic $\mathbb{Z}_N$ symmetry, which projects the leading Coleman-Weinberg potential onto the exponentially suppressed $N$th harmonic. Although the microscopic theory preserves $\mathbb{Z}_N$, the axion-dependent WIMP masses generate unequal equilibrium abundances that freeze-out imprints on the cosmological relic state, thereby breaking the symmetry spontaneously. The resulting relic distribution retains a memory of the initial axion value and generates an unsuppressed finite-density potential that holds the field fixed at early times. As the WIMP density dilutes, the axion rolls toward the minimum of its confining potential, transferring energy from DE to DM at late times. An observer assuming separately conserved components then infers an effective equation of state that crosses below $-1$, without ghosts or violation of the null-energy condition. We present an illustrative cosmological solution with a DESI-like phantom crossing and percent-level suppression of structure growth, and discuss the implications of the WIMP multiplicity and relic distribution for DM searches.

Figures

Figures reproduced from arXiv: 2607.28721 by Bingrong Yu, C\'edric Delaunay, Seung J. Lee, Yuan Yin.

Figure 1
Figure 1. Figure 1: shows the relic abundance distribution across the N sectors for several values of ϵ and xd. The lightest sectors, located near Θi − ¯θ + 2πk/N ∼ π, are prefer￾entially populated, with this hierarchy becoming more pronounced for larger values of qd due to the stronger Boltzmann suppression of the heavier states. We defer the implications of this nonuniform distribution for DM phenomenology to Section VI. In… view at source ↗
Figure 2
Figure 2. Figure 2: shows the evolution of the effective DE po￾tential for different values of the DM density around the critical value n c χ ≡ ΛD/σfd. For nχ ≫ n c χ, the finite￾density potential dominates, trapping the field to its ini￾tial value Θi . When nχ ≲ n c χ, the confining potential starts to dominate and the minimum of the effective po￾tential moves toward the vacuum at Θ = 0 (modulo 2π). V. APPARENT PHANTOM DARK … view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Effective DE equation of state as a function of redshift [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Evolution of the linear growth observable ∆ [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Background diagnostics for the benchmark solution [PITH_FULL_IMAGE:figures/full_fig_p010_5.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

80 extracted references · 58 linked inside Pith

  1. [4]

    Brout et al.,The Pantheon+ Analysis: Cosmological Constraints,Astrophys

    D. Brout et al.,The Pantheon+ Analysis: Cosmological Constraints,Astrophys. J.938(2022), no. 2 110, [arXiv:2202.04077]

  2. [5]

    Rubin et al.,Union Through UNITY: Cosmology with 2,000 SNe Using a Unified Bayesian Framework, Astrophys

    D. Rubin et al.,Union Through UNITY: Cosmology with 2,000 SNe Using a Unified Bayesian Framework, Astrophys. J.986(2025), no. 2 231, [arXiv:2311.12098]. [6]DESCollaboration, T. M. C. Abbott et al.,The Dark 15 Energy Survey: Cosmology Results with∼1500 New High-redshift Type Ia Supernovae Using the Full 5 yr Data Set,Astrophys. J. Lett.973(2024), no. 1 L14...

  3. [9]

    Efstathiou and S

    G. Efstathiou and S. Gratton,A Detailed Description of the CamSpec Likelihood Pipeline and a Reanalysis of the Planck High Frequency Maps,arXiv:1910.00483

  4. [10]

    Carron, M

    J. Carron, M. Mirmelstein, and A. Lewis,CMB lensing from Planck PR4 maps,JCAP09(2022) 039, [arXiv:2206.07773]

  5. [11]

    Rosenberg, S

    E. Rosenberg, S. Gratton, and G. Efstathiou,CMB power spectra and cosmological parameters from Planck PR4 with CamSpec,Mon. Not. Roy. Astron. Soc.517 (2022), no. 3 4620–4636, [arXiv:2205.10869]. [12]ACTCollaboration, M. S. Madhavacheril et al.,The Atacama Cosmology Telescope: DR6 Gravitational Lensing Map and Cosmological Parameters,Astrophys. J.962(2024)...

  6. [14]

    Jiang, D

    J.-Q. Jiang, D. Pedrotti, S. S. da Costa, and S. Vagnozzi,Nonparametric late-time expansion history reconstruction and implications for the Hubble tension in light of recent DESI and type Ia supernovae data, Phys. Rev. D110(2024), no. 12 123519, [arXiv:2408.02365]. [15]DESICollaboration, K. Lodha et al.,Extended dark energy analysis using DESI DR2 BAO mea...

  7. [16]

    W. J. Wolf, C. Garc ´ ıa-Garc ´ ıa, and P. G. Ferreira, Robustness of dark energy phenomenology across different parameterizations,JCAP05(2025) 034, [arXiv:2502.04929]

  8. [17]

    I. D. Gialamas, G. H¨ utsi, M. Raidal, J. Urrutia, M. Vasar, and H. Veerm¨ ae,Quintessence and phantoms in light of DESI 2025,Phys. Rev. D112(2025), no. 6 063551, [arXiv:2506.21542]

  9. [18]

    M. W. Toomey, G. Montefalcone, E. McDonough, and K. Freese,How theory-informed priors affect DESI evidence for evolving dark energy,Phys. Rev. D113 (2026), no. 12 123532, [arXiv:2509.13318]

  10. [19]

    D. H. Lee, W. Yang, E. Di Valentino, S. Pan, and C. van de Bruck,Shape of dark energy: Constraining its evolution with a general parametrization,Phys. Rev. D 113(2026), no. 6 063554, [arXiv:2507.11432]

  11. [20]

    Ratra and P

    B. Ratra and P. J. E. Peebles,Cosmological Consequences of a Rolling Homogeneous Scalar Field, Phys. Rev. D37(1988) 3406

  12. [21]

    Wetterich,Cosmology and the Fate of Dilatation Symmetry,Nucl

    C. Wetterich,Cosmology and the Fate of Dilatation Symmetry,Nucl. Phys. B302(1988) 668–696, [arXiv:1711.03844]

  13. [22]

    R. R. Caldwell, R. Dave, and P. J. Steinhardt, Cosmological imprint of an energy component with general equation of state,Phys. Rev. Lett.80(1998) 1582–1585, [astro-ph/9708069]

  14. [23]

    R. R. Caldwell,A Phantom menace?,Phys. Lett. B545 (2002) 23–29, [astro-ph/9908168]

  15. [24]

    Qiu, Y.-F

    T. Qiu, Y.-F. Cai, and X.-M. Zhang,Null Energy Condition and Dark Energy Models,Mod. Phys. Lett. A 23(2008) 2787–2798, [arXiv:0710.0115]

  16. [25]

    K. J. Ludwick,The viability of phantom dark energy: A review,Mod. Phys. Lett. A32(2017), no. 28 1730025, [arXiv:1708.06981]

  17. [26]

    Moghtaderi, B

    E. Moghtaderi, B. R. Hull, J. Quintin, and G. Geshnizjani,How much null-energy-condition breaking can the Universe endure?,Phys. Rev. D111 (2025), no. 12 123552, [arXiv:2503.19955]

  18. [27]

    R. R. Caldwell and E. V. Linder,Null impact of the null energy condition in current cosmology,JCAP05(2026) 008, [arXiv:2511.07526]

  19. [28]

    S. Das, P. S. Corasaniti, and J. Khoury, Super-acceleration as signature of dark sector interaction,Phys. Rev. D73(2006) 083509, [astro-ph/0510628]

  20. [29]

    Smith, M

    A. Smith, M. Mylova, P. Brax, C. van de Bruck, C. P. Burgess, and A.-C. Davis,A Minimal Axio-dilaton Dark Sector,arXiv:2410.11099

  21. [30]

    Khoury, M.-X

    J. Khoury, M.-X. Lin, and M. Trodden,Apparent w<-1 and a Lower S8 from Dark Axion and Dark Baryons Interactions,Phys. Rev. Lett.135(2025), no. 18 181001, [arXiv:2503.16415]

  22. [31]

    Bedroya, G

    A. Bedroya, G. Obied, C. Vafa, and D. H. Wu,Evolving Dark Sector and the Dark Dimension Scenario, arXiv:2507.03090

  23. [32]

    Wang, R.-G

    J.-Q. Wang, R.-G. Cai, Z.-K. Guo, and S.-J. Wang, Resolving the Planck-DESI tension by nonminimally coupled quintessence,Phys. Rev. D113(2026), no. 8 083534, [arXiv:2508.01759]

  24. [33]

    S. L. Guedezounme, B. R. Dinda, and R. Maartens, Phantom crossing or dark interaction?,JCAP01 (2026) 062, [arXiv:2507.18274]

  25. [34]

    R. Chen, J. M. Cline, V. Muralidharan, and B. Salewicz, Quintessential dark energy crossing the phantom divide, JCAP03(2026) 044, [arXiv:2508.19101]

  26. [35]

    La Penna, A

    L. La Penna, A. Notari, and M. Redi,Mimicking Phantom Dark Energy with Evolving Dark Matter Mass, arXiv:2601.05235

  27. [36]

    Antusch, S

    S. Antusch, S. F. King, and X. Wang,Coupled Dark Energy and Dark Matter for DESI: An Effective Guide to the Phantom Divide,arXiv:2604.08449

  28. [37]

    Delaunay and A

    C. Delaunay and A. Greljo,Natural Phantom Dark Energy from aZ N –Axion,arXiv:2607.06774

  29. [38]

    Khoury, M.-X

    J. Khoury, M.-X. Lin, and M. Trodden,Cosmological Evidence for Dark Axion-Dark Baryon Interactions from Apparent Phantom Crossing,arXiv:2607.16191

  30. [39]

    Tsujikawa,Quintessence: A Review,Class

    S. Tsujikawa,Quintessence: A Review,Class. Quant. Grav.30(2013) 214003, [arXiv:1304.1961]

  31. [40]

    J. A. Frieman, C. T. Hill, A. Stebbins, and I. Waga, Cosmology with ultralight pseudo Nambu-Goldstone bosons,Phys. Rev. Lett.75(1995) 2077–2080, [astro-ph/9505060]

  32. [41]

    J. E. Kim,Axion and almost massless quark as ingredients of quintessence,JHEP05(1999) 022, [hep-ph/9811509]

  33. [42]

    J. E. Kim,Model dependent axion as quintessence with almost massless hidden sector quarks,JHEP06(2000) 016, [hep-ph/9907528]. 16

  34. [43]

    R. Liu, Y. Zhu, W. Hu, and V. Miranda,Phantom mirage from axion dark energy,Phys. Rev. D113 (2026), no. 8 083506, [arXiv:2510.14957]

  35. [44]

    B. W. Lee and S. Weinberg,Cosmological Lower Bound on Heavy Neutrino Masses,Phys. Rev. Lett.39(1977) 165–168

  36. [45]

    Steigman and M

    G. Steigman and M. S. Turner,Cosmological Constraints on the Properties of Weakly Interacting Massive Particles,Nucl. Phys. B253(1985) 375–386

  37. [46]

    M. W. Goodman and E. Witten,Detectability of Certain Dark Matter Candidates,Phys. Rev. D31 (1985) 3059–3063

  38. [47]

    Cirelli, A

    M. Cirelli, A. Strumia, and J. Zupan,Dark Matter, arXiv:2406.01705

  39. [48]

    Hook,Solving the Hierarchy Problem Discretely, Phys

    A. Hook,Solving the Hierarchy Problem Discretely, Phys. Rev. Lett.120(2018), no. 26 261802, [arXiv:1802.10093]

  40. [49]

    Brzeminski, Z

    D. Brzeminski, Z. Chacko, A. Dev, and A. Hook, Time-varying fine structure constant from naturally ultralight dark matter,Phys. Rev. D104(2021), no. 7 075019, [arXiv:2012.02787]

  41. [50]

    J. E. Kim,Weak Interaction Singlet and Strong CP Invariance,Phys. Rev. Lett.43(1979) 103

  42. [51]

    M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Can Confinement Ensure Natural CP Invariance of Strong Interactions?,Nucl. Phys. B166(1980) 493–506

  43. [52]

    Griest and D

    K. Griest and D. Seckel,Three exceptions in the calculation of relic abundances,Phys. Rev. D43(1991) 3191–3203

  44. [53]

    Gondolo and G

    P. Gondolo and G. Gelmini,Cosmic abundances of stable particles: Improved analysis,Nucl. Phys. B360 (1991) 145–179

  45. [54]

    Witten,Theta dependence in the large N limit of four-dimensional gauge theories,Phys

    E. Witten,Theta dependence in the large N limit of four-dimensional gauge theories,Phys. Rev. Lett.81 (1998) 2862–2865, [hep-th/9807109]

  46. [55]

    Vicari and H

    E. Vicari and H. Panagopoulos,Theta dependence of SU(N) gauge theories in the presence of a topological term,Phys. Rept.470(2009) 93–150, [arXiv:0803.1593]

  47. [56]

    Bonati, M

    C. Bonati, M. D’Elia, P. Rossi, and E. Vicari,Theta dependence of 4D SU(N) gauge theories in the large-N limit,Phys. Rev. D94(2016), no. 8 085017, [arXiv:1607.06360]

  48. [57]

    M. C` e, M. Garc ´ ıa Vera, L. Giusti, and S. Schaefer,The topological susceptibility in the large-N limit of SU(N) Yang–Mills theory,Phys. Lett. B762(2016) 232–236, [arXiv:1607.05939]

  49. [58]

    Bonati, M

    C. Bonati, M. D’Elia, H. Panagopoulos, and E. Vicari, Change of theta dependence in 4D SU(N) gauge theories across the deconfinement transition,Phys. Rev. Lett. 110(2013), no. 25 252003, [arXiv:1301.7640]

  50. [59]

    S. R. Coleman and E. J. Weinberg,Radiative Corrections as the Origin of Spontaneous Symmetry Breaking,Phys. Rev. D7(1973) 1888–1910

  51. [60]

    J. A. Frieman, C. T. Hill, and R. Watkins,Late Time Cosmological Phase Transitions 1: Particle Physics Models and Cosmic Evolution,Phys. Rev. D46(1992) 1226–1238. [61]PlanckCollaboration, N. Aghanim et al.,Planck 2018 results. VI. Cosmological parameters,Astron. Astrophys.641(2020) A6, [arXiv:1807.06209]. [Erratum: Astron.Astrophys. 652, C4 (2021)]

  52. [62]

    Banks and N

    T. Banks and N. Seiberg,Symmetries and Strings in Field Theory and Gravity,Phys. Rev. D83(2011) 084019, [arXiv:1011.5120]

  53. [63]

    Harlow and H

    D. Harlow and H. Ooguri,Symmetries in quantum field theory and quantum gravity,Commun. Math. Phys.383 (2021), no. 3 1669–1804, [arXiv:1810.05338]

  54. [64]

    Choi and S

    K. Choi and S. H. Im,Realizing the relaxion from multiple axions and its UV completion with high scale supersymmetry,JHEP01(2016) 149, [arXiv:1511.00132]

  55. [65]

    D. E. Kaplan and R. Rattazzi,Large field excursions and approximate discrete symmetries from a clockwork axion,Phys. Rev. D93(2016), no. 8 085007, [arXiv:1511.01827]

  56. [66]

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

  57. [67]

    Arkani-Hamed, H.-C

    N. Arkani-Hamed, H.-C. Cheng, P. Creminelli, and L. Randall,Extra natural inflation,Phys. Rev. Lett.90 (2003) 221302, [hep-th/0301218]

  58. [68]

    Choi,A QCD axion from higher dimensional gauge field,Phys

    K. Choi,A QCD axion from higher dimensional gauge field,Phys. Rev. Lett.92(2004) 101602, [hep-ph/0308024]

  59. [69]

    Weinberg,The Cosmological Constant Problem,Rev

    S. Weinberg,The Cosmological Constant Problem,Rev. Mod. Phys.61(1989) 1–23

  60. [70]

    D. J. Fixsen,The Temperature of the Cosmic Microwave Background,Astrophys. J.707(2009) 916–920, [arXiv:0911.1955]

  61. [71]

    Archidiacono, E

    M. Archidiacono, E. Castorina, D. Redigolo, and E. Salvioni,Unveiling dark fifth forces with linear cosmology,JCAP10(2022) 074, [arXiv:2204.08484]

  62. [72]

    Bottaro, E

    S. Bottaro, E. Castorina, M. Costa, D. Redigolo, and E. Salvioni,Unveiling Dark Forces with Measurements of the Large Scale Structure of the Universe,Phys. Rev. Lett.132(2024), no. 20 201002, [arXiv:2309.11496]

  63. [73]

    Bottaro, E

    S. Bottaro, E. Castorina, M. Costa, D. Redigolo, and E. Salvioni,From 100 kpc to 10 Gpc: Dark matter self-interactions before and after DESI observations, Phys. Rev. D112(2025), no. 2 023525, [arXiv:2407.18252]

  64. [74]

    Beutler, C

    F. Beutler, C. Blake, M. Colless, D. H. Jones, L. Staveley-Smith, G. B. Poole, L. Campbell, Q. Parker, W. Saunders, and F. Watson,The 6dF Galaxy Survey: z≈0measurement of the growth rate andσ 8,Mon. Not. Roy. Astron. Soc.423(2012) 3430–3444, [arXiv:1204.4725]

  65. [75]

    Howlett, A

    C. Howlett, A. Ross, L. Samushia, W. Percival, and M. Manera,The clustering of the SDSS main galaxy sample – II. Mock galaxy catalogues and a measurement of the growth of structure from redshift space distortions atz= 0.15,Mon. Not. Roy. Astron. Soc.449(2015), no. 1 848–866, [arXiv:1409.3238]

  66. [76]

    Blake et al.,Galaxy And Mass Assembly (GAMA): improved cosmic growth measurements using multiple tracers of large-scale structure,Mon

    C. Blake et al.,Galaxy And Mass Assembly (GAMA): improved cosmic growth measurements using multiple tracers of large-scale structure,Mon. Not. Roy. Astron. Soc.436(2013) 3089, [arXiv:1309.5556]. [77]WiggleZCollaboration, C. Blake et al.,The WiggleZ Dark Energy Survey: Joint measurements of the expansion and growth history at z<1,Mon. Not. Roy. Astron. Soc...

  67. [78]

    A. Pezzotta et al.,The VIMOS Public Extragalactic Redshift Survey (VIPERS): The growth of structure at 0.5< z <1.2from redshift-space distortions in the clustering of the PDR-2 final sample,Astron. Astrophys.604(2017) A33, [arXiv:1612.05645]

  68. [79]

    Okumura et al.,The Subaru FMOS galaxy redshift 17 survey (FastSound)

    T. Okumura et al.,The Subaru FMOS galaxy redshift 17 survey (FastSound). IV. New constraint on gravity theory from redshift space distortions atz∼1.4,Publ. Astron. Soc. Jap.68(2016), no. 3 38, [arXiv:1511.08083]

  69. [80]

    K. Said, M. Colless, C. Magoulas, J. R. Lucey, and M. J. Hudson,Joint analysis of 6dFGS and SDSS peculiar velocities for the growth rate of cosmic structure and tests of gravity,Mon. Not. Roy. Astron. Soc.497(2020), no. 1 1275–1293, [arXiv:2007.04993]

  70. [81]

    S. S. Boruah, M. J. Hudson, and G. Lavaux,Cosmic flows in the nearby Universe: new peculiar velocities from SNe and cosmological constraints,Mon. Not. Roy. Astron. Soc.498(2020), no. 2 2703–2718, [arXiv:1912.09383]

  71. [82]

    Carrick, S

    J. Carrick, S. J. Turnbull, G. Lavaux, and M. J. Hudson,Cosmological parameters from the comparison of peculiar velocities with predictions from the 2M++ density field,Mon. Not. Roy. Astron. Soc.450(2015), no. 1 317–332, [arXiv:1504.04627]

  72. [83]

    Huterer, D

    D. Huterer, D. Shafer, D. Scolnic, and F. Schmidt, TestingΛCDM at the lowest redshifts with SN Ia and galaxy velocities,JCAP05(2017) 015, [arXiv:1611.09862]

  73. [84]

    R. J. Turner, C. Blake, and R. Ruggeri,A local measurement of the growth rate from peculiar velocities and galaxy clustering correlations in the 6dF Galaxy Survey,Mon. Not. Roy. Astron. Soc.518(2022), no. 2 2436–2452, [arXiv:2207.03707]. [85]EuclidCollaboration, I. Ocampo et al.,Euclid: Forecasts onΛCDM consistency tests with growth rate data,arXiv:2507.22780

  74. [86]

    Cirelli, N

    M. Cirelli, N. Fornengo, and A. Strumia,Minimal dark matter,Nucl. Phys. B753(2006) 178–194, [hep-ph/0512090]

  75. [87]

    Cirelli, A

    M. Cirelli, A. Strumia, and M. Tamburini,Cosmology and astrophysics of minimal dark matter,Nucl. Phys. B 787(2007) 152–175, [arXiv:0706.4071]

  76. [88]

    Cirelli and A

    M. Cirelli and A. Strumia,Minimal Dark Matter: Model and results,New J. Phys.11(2009) 105005, [arXiv:0903.3381]

  77. [89]

    Hook and R

    A. Hook and R. Rattazzi,Softening the UV without new particles,Phys. Rev. D108(2023), no. 11 115019, [arXiv:2306.12489]

  78. [90]

    Dolan and R

    L. Dolan and R. Jackiw,Symmetry behavior at finite temperature,Phys. Rev. D9(1974) 3320–3341

  79. [91]

    Weinberg,Gauge and global symmetries at high temperature,Phys

    S. Weinberg,Gauge and global symmetries at high temperature,Phys. Rev. D9(1974) 3357–3378

  80. [92]

    Quiros,Finite temperature field theory and phase transitions, inICTP Summer School in High-Energy Physics and Cosmology, pp

    M. Quiros,Finite temperature field theory and phase transitions, inICTP Summer School in High-Energy Physics and Cosmology, pp. 187–259, 1, 1999. hep-ph/9901312