Post-Inflationary Axion Dark Matter

• Post-Inflationary Axion Dark Matter is produced when PQ symmetry is restored after inflation, leading to axions via misalignment and the decay of topological defects.
• The topic details how vacuum realignment, cosmic string radiation, and domain-wall collapse quantitatively shape the relic density and axion mass window.
• It underpins critical experimental searches and simulations that probe non-linear dynamics, isocurvature fluctuations, and small-scale structure in dark matter.

Post-inflationary axion dark matter refers to models in which the Peccei–Quinn (PQ) symmetry is restored during or after inflation, resulting in the cosmological production of axion dark matter via topological defect dynamics and misalignment. In these scenarios, the axion—a pseudo-Nambu–Goldstone boson introduced to solve the strong CP problem of QCD—accounts for the observed cold dark matter abundance through a combination of vacuum misalignment, radiative emission from global axion strings, and collapse of the axion string–domain-wall network. The post-inflationary scenario leads to distinctive predictions for the axion mass, relic density, isocurvature perturbations, and small-scale structure, and is a focus of both theoretical modeling and next-generation direct search experiments.

1. Post-Inflationary Peccei–Quinn Symmetry Breaking: Axion Production Mechanisms

When the PQ symmetry is restored after inflation, the axion field acquires random initial values in each causally disconnected Hubble patch. As the Universe cools, the axion mass “switches on” around the QCD scale due to nonperturbative QCD effects. Three production channels contribute to the relic axion dark matter density:

1. Vacuum Realignment (Misalignment Mechanism):
   • Coherent oscillations of the axion field begin when the Hubble parameter drops below the axion mass, resulting in cold axion dark matter.
   • For random initial angles averaged over many domains, the contribution is

   $$\Omega_{a}^{\rm (mis)}\,h^2 \simeq (3.8\pm0.6)\times10^{-3}\left(\frac{f_A}{10^{10}\,\mathrm{GeV}}\right)^{1.165}$$

   where $f_A$ is the axion decay constant and the exponent 1.165 reflects the latest QCD lattice determinations of the topological susceptibility (Ringwald, 2018).

2. Radiation from Global Strings:

   • A global network of axion strings forms at the PQ phase transition. In the scaling regime, these strings radiate axions of energy around the horizon scale.
   • The relic contribution is

   $$\Omega_{a}^{\rm (S)}\,h^2 \approx 7.8^{+6.3}_{-4.5}\times10^{-3}\,N_{\rm DW}^2\,\left(\frac{f_A}{10^{10}\,\mathrm{GeV}}\right)^{1.165}$$

   with $N_{\rm DW}$ the domain-wall number.

3. Collapse of the String–Wall System:

   • The emergence of the QCD potential produces domain walls attached to each string. For the favored case $N_{\rm DW}=1$, the composite network is unstable and collapses, radiating axions.
   • The contribution is

   $$\Omega_{a}^{\rm (C)}\,h^2 \approx 3.9^{+2.3}_{-2.1}\times10^{-3}\left(\frac{f_A}{10^{10}\,\mathrm{GeV}}\right)^{1.165}$$

Summing these, the total axion relic abundance for $N_{\rm DW}=1$ is

$$\Omega_a h^2 \simeq (1.6^{+1.0}_{-0.7})\times10^{-2}\left(\frac{f_A}{10^{10}\,\mathrm{GeV}}\right)^{1.165}$$

and requiring $\Omega_a h^2 \approx 0.12$ yields the post-inflationary axion mass window $m_A \simeq(58\text{--}150)\,\mu\mathrm{eV}$ (Ringwald, 2018).

2. Nonlinear Dynamics: Topological Defects and Simulations

The post-inflationary scenario demands treatment of axion field inhomogeneities and the nonlinear evolution of the axion field network:

• Cosmic Strings and Domain Walls:

  • The defect network’s evolution requires large-scale numerical simulations to capture the scaling behavior, axion radiation spectrum, and ultimate relic density. The number of axions radiated per unit string length, the average string density parameter $\zeta$, and energy spectrum are key simulation outputs.
• Relic Abundance Suppression/Enhancement:
  • For a temperature-independent axion mass ($n=0$), the total axion abundance is enhanced by ≈25% compared to misalignment-only estimates; for QCD axions ($n\approx6$), it is suppressed by a factor of ≈6 relative to naïve misalignment (O'Hare et al., 2021). This arises from the full non-linear decay dynamics of the string-wall network.
• Axitons and Minicluster Seeds:
  • Nonlinear phenomena such as “axitons” (high-density, self-bound field lumps) and gravitational collapse of overdensities lead to formation of axion miniclusters—a robust prediction with implications for small-scale dark matter structure (O'Hare et al., 2021).

3. Isocurvature Fluctuations and Small-Scale Structure

Isocurvature perturbations are a defining feature of post-inflationary axion dark matter:

• Order-One Isocurvature:
  • Each Hubble patch’s random axion angle yields density fluctuations of order unity on the horizon scale at the onset of oscillations, mapping to a white-noise isocurvature power spectrum for wavenumbers below the cutoff $k_\text{osc} = a(T_\text{osc}) H(T_\text{osc})$ (Shimabukuro et al., 2020, Feix et al., 2019).
• Minihalo and Minicluster Formation:
  • Enhanced small-scale power translates to an abundance of low-mass axion minihalos and miniclusters. Their characteristic masses, distribution functions, and survival depend on the details of the axion mass’s temperature dependence and on the string-wall network’s decay.
• Observational Probes:
  • The 21 cm forest technique can probe the mass window $10^{-18} \lesssim m_a \lesssim 10^{-12}$ eV for temperature-independent cases, and up to $10^{-6}$ eV for QCD-like scenarios by measuring absorption toward high-redshift radio sources, due to the increased minihalo abundance (Shimabukuro et al., 2020).

4. Theoretical and Model-Dependent Uncertainties

Theoretical predictions and the favored axion mass range in post-inflationary scenarios are sensitive to several uncertainties:

• Misalignment Angle Distribution:
  • The standard estimate assumes a uniform distribution $\langle\theta^2\rangle = \pi^2/3$; deviations or strong anharmonic corrections can impact the yield.
• Domain Wall Number and Explicit PQ Breaking:
  • For $N_{\rm DW}>1$, string-wall networks are stable unless a small explicit PQ-violating potential is present; the dynamics and collapse timing are then sensitive to the magnitude and phase of this term. In concrete model realizations (e.g., DFSZ with $N=10$), small Planck-suppressed operators break the degeneracy, allowing post-inflationary viability for larger wall numbers but requiring minimal tuning (Ringwald et al., 2015).
• String-Network Extrapolation:
  • The extrapolation from numerically accessible string tensions ($\ln(m_\rho t) \sim$ few) to physically relevant values ($\sim 50$) introduces systematic uncertainties.
• Lattice QCD Inputs:
  • The temperature dependence of the QCD topological susceptibility, encoded in the exponent of the misalignment yield, relies on the latest lattice QCD calculations.

These factors collectively determine the mass window for which axions saturate the observed cold dark matter density.

5. Phenomenological Consequences and Experimental Probes

Post-inflationary axion models have sharply defined phenomenological windows and experimental targets:

• Benchmark Mass Window:
  • For QCD axions with $N_{\rm DW}=1$, the consensus prediction is $m_a \simeq 58\text{--}150\,\mu\mathrm{eV}$ from combined misalignment, string, and wall contributions. The lower limit ($m_a>28\,\mu\mathrm{eV}$) ensures the universe does not overclose from misalignment alone (Ringwald, 2018).
• Haloscope Sensitivity:
  • This mass window is the principal target for next-generation haloscope and dielectric-haloscope experiments such as MADMAX, ADMX, and HAYSTAC. MADMAX, in particular, is designed to probe the 40–200 μeV range with a projected sensitivity reaching DFSZ and post-inflationary QCD axion couplings (Hubaut, 8 Sep 2025).
• Complementary Astrophysical and Laboratory Observations:
  • White-dwarf cooling hints and solar axion searches fall within the predicted mass/coupling regions of certain post-inflationary models with accidental PQ symmetries (Ringwald et al., 2015).
  • Constraints from the Planck CMB on isocurvature power exclude post-inflationary ALP dark matter with $m_a < 10^{-20}\text{--}10^{-16}$ eV, with these limits only improvable by future CMB-S4 and 21 cm intensity mapping (Feix et al., 2019).
• Model Variants and Mass Range Extensions:
  • Scenarios with heavy quark domination, late inflation, or explicit PQ-breaking sources can relax the lower bound on the axion mass to $m_a \sim 10^{-8}$ eV by introducing entropy dilution or delayed network reentry (Harigaya et al., 2022, Cheek et al., 7 May 2025, Azatov et al., 21 Oct 2025).
  • In models with composite axions or those where the network enters the horizon only after QCD, the axion relic abundance and allowed decay constant can differ significantly, potentially enabling $f_a$ as low as the astrophysical lower bounds while predicting larger minicluster masses (Harigaya et al., 2022, Azatov et al., 21 Oct 2025).

6. Model-Building, Extensions, and Open Issues

• Model Solutions to Wall Problem:
  • Mechanisms such as the Lazarides–Shafi mechanism or introduction of discrete gauge symmetries can reduce $N_{\rm DW}$ to unity, ensuring a safe cosmology in post-inflationary scenarios (Lazarides et al., 2020, Ringwald et al., 2015).
• Fermionic and Composite Dark Matter:
  • Some GUT-inspired models with post-inflationary PQ breaking allow for an additional stable non-thermal fermionic dark matter component, possibly comparable to the axion’s contribution (Lazarides et al., 2020). Composite axion models offer UV-complete frameworks that naturally realize $N_{\rm DW}=1$ and evade relic overclosure via a secondary inflationary epoch (Azatov et al., 21 Oct 2025).
• Small-Scale Structure as a Probe:
  • The spectrum and distribution of axion miniclusters are sensitive, model-dependent diagnostics. Post-inflationary scenarios generically predict a non-trivial, possibly observable, halo-mass function and streaming fine-grained substructure (Miguel et al., 2024, O'Hare et al., 2021).
• Discriminating Observables:
  • In models where heavy quark decay dilutes the axion abundance, standard observables ($m_a$, $\Delta N_\text{eff}$) may not differentiate between pre- and post-inflationary PQ breaking. Instead, gravitational-wave backgrounds imprinted by early matter domination may provide the needed discriminant for these cosmologies (Cheek et al., 7 May 2025).

7. Summary Table: Key Parameters and Observables

Mechanism/Model Class	Axion Mass Window / $f_a$	Dominant Constraint
Standard post-inflationary QCD axion, $N_{\rm DW}=1$	$m_a\simeq58\text{--}150\,\mu\mathrm{eV}$ <br/>$f_a\simeq(3.8\text{--}9.9)\times10^{10}$ GeV	Relic abundance, cosmology
DFSZ with explicit PQ breaking ($N_{\rm DW}>1$)	$m_a\sim3.5\text{--}4.2\,\mathrm{meV}$ <br/>$f_a\sim1.5\times10^9$ GeV	Need for wall collapse, white-dwarf cooling
Late network reentry (“diluted string-walls”)	$f_a\sim10^8\text{--}10^{12}$ GeV <br/>$m_a\gtrsim10^{-7}$ eV	Network horizon entry, mini-halos
Models with HQ domination/dilution	$m_a\gtrsim10^{-8}$ eV <br/>$f_a\gtrsim10^{14}$ GeV	BBN, GWs, entropy injection (Cheek et al., 7 May 2025)

This structuring of axion dark matter phenomenology in the post-inflationary regime encapsulates the key physics governing the mass window, cosmological constraints, theoretical uncertainties, and experimental probes, as established in comprehensive reviews and targeted analyses (Ringwald, 2018, O'Hare et al., 2021, Feix et al., 2019, Ringwald et al., 2015, Hubaut, 8 Sep 2025, Cheek et al., 7 May 2025, Azatov et al., 21 Oct 2025).