---
title: 'Axion-Like Early Dark Energy: Theory, Methods, Results'
url: https://www.emergentmind.com/topics/axion-like-early-dark-energy-ede
type: topic
---

# Axion-Like Early Dark Energy: Theory, Methods, Results

Axion-like early dark energy (EDE) is a class of pre-recombination dark-energy models in which a pseudo-Nambu–Goldstone scalar field is initially frozen by Hubble friction, contributes a transient vacuum-like energy component near matter–radiation equality, and then becomes dynamical at a critical redshift \(z_c\), after which its energy density redshifts away faster than matter. In the standard implementation, the field is governed by an axion-like potential \(V(\phi)=m^2f^2[1-\cos(\phi/f)]^n\), usually with \(n=3\), so that after oscillations the effective equation of state is \(w\simeq (n-1)/(n+1)=1/2\) [2202.13379]. Its cosmological role is to raise the pre-recombination expansion rate, reduce the sound horizon \(r_s\), and thereby permit a larger CMB-inferred \(H_0\) at fixed acoustic angle. Subsequent work has treated axion-like EDE both as a candidate resolution of the Hubble tension and as a target of increasingly stringent tests from Planck, ACT, SPT, BAO, supernovae, CMB lensing, and full-shape large-scale-structure likelihoods, while also extending the framework to dark-sector triggers, cosmic birefringence, and multi-field realizations [2311.00524].

## 1. Field-theoretic structure and background dynamics

A canonical formulation writes the total action as
\[
S=\int d^4x\sqrt{-g}\,\Bigl[\tfrac12R-\tfrac12(\partial\phi)^2-V(\phi)\Bigr]+S_{\rm m}[g_{\mu\nu},\Psi_{\rm m}],
\]
with the axion-like potential
\[
V(\phi)=m^2f^2\bigl[1-\cos(\phi/f)\bigr]^n,
\]
and \(n=3\) in the main observational studies [2312.04007]. The homogeneous field obeys the Klein–Gordon equation
\[
\ddot\phi+3H\dot\phi+V_{,\phi}=0,
\]
while the Friedmann equation includes \(\rho_\phi=\tfrac12\dot\phi^2+V(\phi)\) and \(p_\phi=\tfrac12\dot\phi^2-V(\phi)\) alongside the standard radiation, baryon, dark-matter, and late-time dark-energy sectors [2404.16805].

The standard phenomenological picture is that \(\phi\) is frozen at early times when \(H\gg m\), so \(\rho_\phi\) is approximately constant. The field unfreezes when \(H(z_c)\simeq m\), or more precisely \(H(z_c)\simeq m_{\rm eff}(\theta_i)\) when the initial misalignment angle is important, and then rolls and oscillates. For \(n=3\), the post-oscillation energy density redshifts faster than matter; this rapid dilution is one of the reasons the \(n=3\) potential is repeatedly adopted in data analyses [2202.13379].

Operationally, most cosmological fits trade the microscopic parameters \(\{m,f\}\) for the phenomenological triplet \(\{f_{\rm EDE},z_c,\theta_i\}\), where \(\theta_i\equiv\phi_i/f\) is the initial displacement and
\[
f_{\rm EDE}\equiv \frac{\rho_\phi(z_c)}{\rho_{\rm tot}(z_c)}.
\]
This parameterization is designed to isolate the quantities that directly control the transient energy injection, its timing, and its perturbation dynamics [2312.04007].

## 2. Cosmological mechanism and perturbative signatures

The principal effect of axion-like EDE is a temporary enhancement of \(H(z)\) before recombination. Because
\[
r_s(z_*)=\int_{z_*}^{\infty}\frac{c_s(z)}{H(z)}\,dz,
\]
a larger pre-recombination expansion rate reduces the sound horizon. At fixed measured acoustic angle \(\theta_s=r_s/D_A\), this permits a larger inferred \(H_0\) than in \(\Lambda\)CDM [2109.04451]. This is the core mechanism by which EDE addresses the Hubble tension.

The same injection also perturbs the diffusion scale, damping tail, and early-ISW contribution. Accordingly, axion-like EDE is not a purely background-level modification. In synchronous gauge, the field perturbation satisfies
\[
\delta\phi_k''+2\mathcal H\,\delta\phi_k'+\bigl[k^2+a^2V_{,\phi\phi}\bigr]\delta\phi_k
= -\frac12 h'\phi',
\]
and these perturbations feed back into the CMB TT, TE, and EE spectra [2310.16800]. Several analyses therefore emphasize that accurate constraints require solving the full field dynamics rather than using a purely effective fluid description [2109.04451].

Jiang and Piao found that the compensating parameter shifts induced by pre-recombination injection empirically satisfy
\[
\delta n_s \approx 0.3\,\frac{\delta H_0}{H_0},
\]
so that increasing \(H_0\) toward \(\sim 73\) km/s/Mpc pushes \(n_s\) toward unity [2202.13379]. In their combined CMB runs, \(n_s\) becomes consistent with \(1\) at approximately the \(1\sigma\) level. This connection ties axion-like EDE not only to the Hubble tension but also to the question of whether a Harrison–Zeldovich primordial spectrum can be made compatible with modern data.

The detailed observational leverage comes disproportionately from the small-scale CMB. SPT-3G and ACT analyses stress the role of the damping scale \(\theta_D\), the shape of the high-\(\ell\) TT spectrum, and the interplay between TT and polarization. One SPT-3G study reports that ACT DR4 prefers \(100\,\theta_D=0.16327\pm0.00051\), whereas Planck gives \(0.16161\pm0.00019\), while SPT-3G 2018 agrees with Planck with \(100\,\theta_D=0.16202\pm0.00051\) [2309.03265]. This suggests that the viability of the model is controlled at least as much by damping-tail and polarization residuals as by the overall shrinkage of \(r_s\).

## 3. Inference pipelines and data combinations

Axion-like EDE has been constrained with several closely related but non-identical inference pipelines. The dominant pattern is the use of modified Boltzmann solvers derived from CLASS, including CLASS_EDE and AxiCLASS, occasionally with CAMB-based implementations, coupled to samplers such as Cobaya, MontePython, or CosmoMC [2202.13379]. Reported convergence thresholds vary across analyses, including Gelman–Rubin criteria \(R-1<0.1\), \(R-1<0.05\), \(R-1<0.03\), and \(R-1<0.01\), reflecting differing posterior geometries and precision targets [2309.03265].

The basic EDE priors are typically flat in
\[
f_{\rm EDE},\qquad \log_{10}z_c,\qquad \theta_i,
\]
with representative choices such as \(f_{\rm EDE}\in[0,0.3]\) or \([0,0.5]\), \(\log_{10}z_c\in[3,4]\) or \([3,4.5]\), and \(\theta_i\in[0,\pi]\) [2408.09521]. Some analyses also vary neutrino mass, dark-matter properties, or late-time dark-energy parameters, but the core parameterization remains the same.

The data combinations span a wide range of CMB and low-redshift inputs: Planck 2018 or Planck NPIPE primary spectra, Planck lensing, ACT DR4 or DR6 TT/TE/EE or lensing, SPT-3G 2018 temperature and polarization, BAO from 6dFGS, SDSS MGS, BOSS DR12, eBOSS DR16, or DESI Y1, Pantheon or Pantheon+ supernovae, full-shape BOSS likelihoods within the EFT of LSS, and optional local-distance-ladder priors from SH0ES or H0DN [2311.00524].

Methodological details materially affect the resulting posteriors. Jiang and Piao truncated Planck high-\(\ell\) TT to \(\ell_{\rm TT}\le 1000\) to avoid the known over-smoothing anomaly at \(\ell>1000\), and found best-fit points with BOBYQA and posterior visualizations with GetDist [2202.13379]. By contrast, the NPIPE analysis used a new high-\(\ell\) CamSpec-style TT, TE, EE likelihood cleaned with a 545 GHz dust template and explicitly compared Bayesian posteriors with profile likelihoods to assess prior-volume effects [2311.00524].

A notable methodological development is the use of neural-network emulators. One DESI-era analysis employed the TensorFlow package “cosmopower,” trained on \(\sim 2\times 10^5\) CLASS_EDE runs, to reproduce observables at sub-percent accuracy while running \(\sim 100\times\) faster than full Boltzmann solutions, enabling convergence in \(\sim 10\) hours rather than weeks [2404.16805].

## 4. Observational status and dataset dependence

The empirical status of axion-like EDE is not uniform across datasets. Some ACT-centered or Pantheon\(+\)-inclusive analyses find a nonzero EDE fraction and \(H_0\) in the low \(70\)s, whereas SPT-, NPIPE-, DESI-, and EFT-BOSS-centered analyses pull the model much closer to the \(\Lambda\)CDM limit.

| Analysis | Representative constraint | Interpretation |
|---|---|---|
| ACT DR4 + Planck \(\ell<650\) + lensing + BAO [2109.04451] | \(f_{\rm EDE}=0.091^{+0.020}_{-0.036}\), \(H_0=70.9^{+1.0}_{-2.0}\) km/s/Mpc | Preference for EDE at \(>99.7\%\) CL |
| Planck2018 \((\ell_{\rm TT,max}=1000)\) + ACT DR4 + SPT-3G Y1 [2202.13379] | \(f_{\rm EDE}=0.127^{+0.058}_{-0.034}\), \(H_0=72.4^{+2.2}_{-1.7}\), \(n_s=0.9923^{+0.0098}_{-0.0086}\) | \(\Delta\chi^2\simeq -11\), mainly ACT-driven |
| SPT-3G 2018 TT/TE/EE [2309.03265] | \(f_{\rm EDE}<0.172\) at \(95\%\) CL | No preference for axion-like EDE |
| Planck NPIPE + BAO + Pantheon + lensing [2311.00524] | \(f_{\rm EDE}<0.061\), \(H_0=68.11^{+0.47}_{-0.82}\) km/s/Mpc | Residual \(3.7\sigma\) tension with SH0ES |
| Planck + CMB lensing + DESI BAO [2404.16805] | \(f_{\rm EDE}<0.091\) at \(95\%\) CL, \(H_0=69.14^{+0.68}_{-1.10}\) km/s/Mpc | No significant preference for EDE |
| Planck 2018 + EFT-BOSS [2006.11235] | \(f_{\rm EDE}<0.072\) at \(95\%\) CL, \(H_0=68.54^{+0.52}_{-0.95}\) km/s/Mpc | Full-shape LSS sharply tightens constraints |

A further Pantheon\(+\) and EFTofLSS reassessment found that axion-like EDE remained comparatively competitive under the combination Planck+ext-BAO+EFTofBOSS+EFTofeBOSS+Pantheon\(+\)+\(M_b\), reporting \(f_{\rm EDE}(z_c)=0.116^{+0.023}_{-0.021}\), \(\log_{10}(z_c)=3.69^{+0.20}_{-0.16}\), \(H_0=71.67\pm0.77\) km/s/Mpc, \(\Delta\chi^2=-35.1\), \(\Delta{\rm AIC}=-29.1\), and a residual tension \(Q_{\rm DMAP}=2.5\sigma\) [2310.16800]. This stands in marked contrast to the tighter NPIPE and DESI-era exclusions.

Taken together, the literature exhibits a recurring tension between ACT-favoring analyses and those dominated by Planck high-\(\ell\), SPT, DESI BAO, or EFT-based galaxy-clustering information. Several papers identify the decisive modes as high-\(\ell\) TT or TE/EE multipoles, especially around \(\ell\sim 800\)–\(1500\), the damping tail, and the detailed shape of the polarization spectra. The disagreement has been discussed as a possible physical preference, a systematic effect, or a statistical fluctuation, but no single interpretation is established across the cited studies.

A related large-scale-structure conclusion is comparatively stable: when full-shape information is used rather than compressed BAO+\(f\sigma_8\) summaries, the parameter space that allows the \(f_{\rm EDE}\sim 0.1\) injection needed for a full Hubble-tension resolution contracts substantially [2006.11235].

## 5. Extensions, variants, and associated phenomena

A number of extensions have tested whether the core axion-like EDE mechanism can be made more concordant with other observables. Yao and Meng replaced CDM with dark matter having a constant equation of state \(w_{\rm dm}\), so that \(p_{\rm dm}=w_{\rm dm}\rho_{\rm dm}\) and \(\rho_{\rm dm}\propto a^{-3(1+w_{\rm dm})}\). For the combined CMB+BAO+Pantheon+\(H_0\)+\(S_8\) dataset they obtained \(\log_{10}(z_c)=3.67\pm0.22\), \(f_{\rm EDE}=0.066^{+0.025}_{-0.047}\), \(\theta_i=2.48^{+0.53}_{-0.27}\), \(w_{\rm dm}=0.00042^{+0.00069}_{-0.00043}\), \(H_0=71.11^{+0.88}_{-1.00}\) km/s/Mpc, and \(S_8=0.814^{+0.011}_{-0.012}\). They found that the Hubble tension is reduced to \(\sim 1.4\sigma\), but \(w_{\rm dm}\) is driven close to zero, the \(S_8\) tension remains at \(\sim 2.3\sigma\), and AIC mildly favors standard EDE over the \(w_{\rm dm}\neq 0\) extension [2312.04007].

Lin et al. proposed a dark-matter-triggered early dark sector in which \(m_{\rm DM}(\phi)=m_0[1+g\,\phi^2/M_{\rm pl}^2]\), so that the rise of dark matter near equality itself triggers the rolling of the EDE field. In this trigger EDS model the coincidence problem is addressed at the background level without fine tuning of the coupling or the initial conditions, while the best-fit model reaches \(H_0\simeq 70.7\) km/s/Mpc and \(f_{\rm EDE}\simeq 0.11\) with \(1+z_c\simeq 10^{3.85}\). The same work emphasizes that Planck selects a specific range of initial field positions, \(\phi_i/f\approx 1.1\), because of perturbative effects in the high-\(\ell\) acoustic structure [2212.08098].

A different extension couples the axion-like EDE field to electromagnetism through a Chern–Simons term,
\[
\mathcal L_{\rm CS}= -\frac{g_{\phi\gamma}}{4}\phi F_{\mu\nu}\tilde F^{\mu\nu},
\]
producing cosmic birefringence. In a joint fit to Planck TT/EE, Planck lensing, a measured EB spectrum, SDSS DR12 BAO, and a local \(H_0\) prior, the \(n=3\) model yielded a marginalized mean \(f_{\rm EDE}=0.180\), \(\log_{10}(z_c)=3.487\), \(g_{\phi\gamma}=0.170\), and derived \(H_0=71.26\) km/s/Mpc, with the best-fit rotation angle \(\alpha_0\approx 0.059^\circ\) and a reduced \(\chi^2\approx 0.94\) for 72 degrees of freedom in the EB fit [2408.09521]. The same analysis notes that current EB measurements are noise-limited and identifies LiteBIRD, PICO, CMB-S4, and AliCPT as decisive future tests.

The most direct modification of the axion-like sector is to add more fields. A 2026 analysis of a two-field model with \(V_i(\phi_i)=m_i^2f_i^2[1-\cos(\phi_i/f_i)]^3\) reported that splitting the EDE injection into two smoother peaks improves the fit to Planck NPIPE high-\(\ell\) data, raises the best-fit \(H_0\) by \(\sim 1.4\sigma\) relative to the one-field case, and reduces the residual tension with H0DN to \(1.5\sigma\), with no substantial improvement beyond two fields [2604.13535]. This suggests that the sharpness of the single-field injection history is itself part of the difficulty.

By contrast, modifications of the late-time dark-energy equation of state appear to have little leverage within axion-like EDE. In analyses that replaced \(w=-1\) by CPL or oscillating \(w(z)\) ansätze, the best fit remained compatible with \(w_0=-1\), \(w_a=0\), and the change in \(H_0\) was \(\Delta H_0\lesssim 0.1\) km/s/Mpc relative to the pure \(w=-1\) case [2201.07079].

## 6. UV realizations, theoretical constraints, and unresolved issues

The microphysical status of axion-like EDE has been investigated most explicitly in string-theoretic constructions. McDonough and Scalisi showed in a supergravity toy model that a potential of the form
\[
V(\phi)=\Lambda+V_0[1-\cos(\phi/f)]^3
\]
can arise from a delicate balance of separate non-perturbative effects. In a KKLT-type compactification with the EDE scalar realized as a two-form axion, however, they found that successful embedding with stabilized moduli requires restrictive assumptions on the Pfaffians and on the exponents of the non-perturbative terms that generate the EDE harmonics [2209.00011].

A more detailed Type IIB construction with full closed-string moduli stabilization considers both \(C_4\) and \(C_2\) axions in KKLT and in the Large Volume Scenario. That analysis argues that the most promising candidates are \(C_2\) axions in LVS with three non-perturbative corrections to the superpotential generated by gaugino condensation on D7-branes with non-zero world-volume fluxes. It reports that the phenomenologically desired values \(f\approx 0.2\,M_P\), \(V_0\approx({\rm few\ meV})^4\), and \(m\sim 10^{-27}\) eV can be achieved without tuning the microscopic prefactors if the relevant axion violates the mild axionic weak gravity conjecture, whereas WGC-respecting candidates require fine tuning [2303.03414].

These UV studies align with a broader theoretical concern already visible at the effective level: the coincidence problem. Standard axion-like EDE assumes that the field becomes dynamical near equality even though it is not, in its minimal form, explicitly linked to the matter or radiation densities. Trigger models seek to address this by coupling the field to dark matter [2212.08098], while string constructions seek a controlled origin for the required ultra-light mass scale and harmonic structure [2209.00011].

Three empirical issues remain persistent across the literature. First, the model is not uniformly favored by current data: some analyses find \(f_{\rm EDE}\sim 0.1\), but others set \(f_{\rm EDE}<0.061\) or \(<0.072\) at \(95\%\) confidence and recover \(H_0\) values close to \(\Lambda\)CDM [2311.00524]. Second, standard axion-like EDE generally raises \(\omega_{\rm cdm}\), \(\sigma_8\), and \(S_8\), so relief of the Hubble tension does not automatically imply improved global concordance [2312.04007]. Third, late-time freedom in \(w(z)\) does not materially change this conclusion within the standard axion-like setup [2201.07079].

The resulting picture is technically well defined but observationally unsettled. Axion-like EDE remains one of the most extensively developed pre-recombination solutions to the Hubble tension, with a clear Lagrangian description, identifiable perturbative signatures, and active UV model-building. At the same time, its status is increasingly controlled by the consistency of high-\(\ell\) CMB spectra, full-shape large-scale-structure data, and the degree to which additional dark-sector structure is admitted.

Source: https://www.emergentmind.com/topics/axion-like-early-dark-energy-ede