Axion Dark Matter Distribution

• Axion dark matter distribution is the spatial and velocity arrangement of axion particles shaped by production mechanisms, self-interactions, and gravitational assembly.
• It encompasses models from smooth halo distributions to fine-grained streams, caustic rings, dark disks, and compact clumps, guiding direct-detection strategies.
• Recent analyses integrate analytic methods, simulations, and experimental data to refine likelihood frameworks for identifying axion signals in haloscope experiments.

Axion dark matter distribution refers to the spatial and momentum-space arrangement of axions or @@@@1@@@@ as the dominant component of the cosmological dark matter (DM). Unlike conventional cold dark matter (CDM), axion DM exhibits phase-space structure determined by its production mechanism, self-interactions, coherence properties, and gravitational assembly history. Recent decades have seen the emergence of diverse theoretical models and observational strategies that elucidate the expected axion dark matter distribution, ranging from smooth halos and fine-grained streams to compact clumps, solitons, and complex multi-axion “axiverse” scenarios. This article provides a comprehensive technical overview, synthesizing analytic results, simulations, and phenomenological implications from leading research.

1. Smooth Halo Models and Standard Likelihood Formalism

The bulk of galactic axion dark matter is often modeled as a smooth, collisionless field following a quasi-equilibrium distribution, parameterized by the Standard Halo Model (SHM). In direct-detection experiments—e.g., axion haloscopes—this motivates the use of a likelihood-based framework for analyzing time-series data $\Phi(t)$ proportional to the local axion field $a(t)$. After Fourier transformation and forming the power spectral density (PSD), the signal expectation at each frequency bin $\omega_k$ is given by

$$\lambda_k(\theta) = A\, \left(\pi\, f(v)/[m_a v]\right)_{v=\sqrt{2\omega_k/m_a-2}} + \lambda_B(\omega_k)$$

where $$A\propto g_{a\gamma\gamma}^2 B^2 V^2 \rho_\mathrm{DM}$$ encodes the axion-photon coupling, detector geometry, and local DM density, $f(v)$ is the local axion speed distribution, and $\lambda_B$ is the background noise PSD (Foster et al., 2017).

The canonical SHM is the isotropic, boosted Maxwell-Boltzmann:

$$f_\mathrm{SHM}(v|v_0, v_\mathrm{obs}) = \frac{v}{\sqrt{\pi}\, v_0 v_\mathrm{obs}}\,\exp\left[-\frac{v^2+v_\mathrm{obs}^2}{v_0^2}\right]\,\sinh\left(\frac{2v v_\mathrm{obs}}{v_0^2}\right)$$

with dispersion $\sigma_v = v_0/\sqrt{2}$, lab speed $v_\mathrm{obs}$, and sometimes truncated at escape speed $v_\mathrm{esc}$. Substituting $f_\mathrm{SHM}$ into the signal model yields the standard “Maxwellian-folded” line in frequency/energy, with width $\sim m_a v_0^2$ (Foster et al., 2017).

Annual modulation and gravitational focusing by the Sun introduce time-dependent perturbations to $f(v)$; solar gravity induces a correction $\delta f_\mathrm{GF}(v, t)$ computed by Liouville’s theorem, generating time-dependent effective $f_\mathrm{eff}(v, t)$ that can be directly fit from data (Foster et al., 2017).

2. Fine-Grained Phase-Space Structure and Streams

On solar-system scales, the local axion dark matter phase-space is generically a superposition of an extremely large number of cold streams, resulting from the folding of the initial CDM sheet under gravity. The distribution is

$$\rho(x, v) = \sum_{i=1}^N \rho_i\, \delta^3(v - v_i)$$

with $N_\mathrm{streams}\sim10^{14}$ and $N_\mathrm{dense}\sim10^6$ for $\rho_i/\rho_\mathrm{DM}\gtrsim10^{-8}$ at the solar radius (O'Hare et al., 18 Sep 2025). The local velocity distribution in the laboratory frame is then

$$f(\mathbf v) = (1-\epsilon)f_\mathrm{bg}(\mathbf v+\mathbf v_\mathrm{lab}) + \sum_{i=1}^N \frac{\rho_i}{\rho_\mathrm{DM}}\, \delta^3(\mathbf v - (\mathbf v_i - \mathbf v_\mathrm{lab}))$$

where $f_\mathrm{bg}$ is the SHM or alternative background, and $\epsilon$ is the net substructure fraction.

Fine-grained and minicluster streams manifest in ultra-high-spectral-resolution haloscope data as narrow frequency spikes, with width set by the stream velocity dispersion $\sigma_\mathrm{str}$ (for CDM, $\sigma_\mathrm{str}\ll10^{-3}$ km/s) (O'Hare et al., 18 Sep 2025). The lineshape and statistical weight of these streams critically influence discovery significance, with single streams of $\rho_i/\rho_\mathrm{DM}\sim10^{-3}$ greatly enhancing sensitivity for sufficiently long integration times (O'Hare et al., 18 Sep 2025).

3. Substructure: Caustic Rings, Dark Disks, and Clumps

Axion dark matter may display additional substructure, including caustic rings, thick disks, and solitonic or self-gravitating clumps:

• Caustic rings: In tidal-torque-induced Bose–Einstein condensates (BECs), axions thermalize into the lowest energy state consistent with acquired angular momentum, yielding a velocity field corresponding to rigid rotation on the turnaround sphere. The resulting inner caustics have ring (D$_{-4}$ elliptic umbilic) geometry, with predicted radii

$$a_n \simeq 40\,\mathrm{kpc}\,\frac{1}{n}\,\frac{v_\mathrm{rot}}{220\,\mathrm{km/s}}\,\frac{j_\mathrm{max}}{0.18}$$

where $n$ labels the inflow cycle and $j_\mathrm{max}\simeq0.18$ (Sikivie, 2010, Banik et al., 2013). Observational evidence appears as narrow features in rotation curves, IRAS ring-signatures, and stellar overdensities.

• Dark disks: A fraction $x_\mathrm{dd}$ of the local density in a co-rotating disk, with lag velocity $v_\mathrm{lag}\sim50$ km/s and dispersion $v_\mathrm{dd}$, can dominate the axion signal in certain models (Foster et al., 2017).
• Clumps and miniclusters: For QCD axions born after inflation, misalignment and network nonuniformity induce minicluster formation, a fraction of which survive as "axion stars" or dense clumps with mass up to $M_\mathrm{max}\sim10^{-15}M_\odot$ and radius $R_{99}\sim10^{2}$ m (Barranco et al., 2012, Schiappacasse et al., 2017, Vicens et al., 2018). Spherically-symmetric solutions exhibit a "soliton" ground state and a corona with $\rho(r)\propto r^{-5/3}$.

A summary table of substructure types and key parameters:

Structure	Typical Mass Scale	Fraction of DM
Caustic rings	Galactic-scale, diffuse	$\sim$100%
Dense clumps	$10^{17}$–$10^{19}$ kg	$\lesssim10^{-4}$
Minicluster streams	$10^{-13}$–$10^{-12}M_\odot$	Model-dependent
Dark disk	--	$\lesssim20\%$

Additional astrophysical substructures (e.g., Sagittarius streams) as low as a few percent of the local density are detectable at $\sim 2\sigma$ significance in haloscope experiments (Foster et al., 2017).

4. Nonlinear Structures: Axion Quark Nuggets and "Drops"

Nonlinear QCD-era physics can produce baryon- or anti-baryon-rich compact objects—axion quark nuggets (AQNs)—via domain wall collapse. The resultant size-spectrum is a power law

$$\frac{dN}{dB} \propto B^{-\alpha},\quad \alpha = 1 - \frac{\beta\delta + 1}{3(\beta+1)}$$

with $B$ the baryon number per nugget, parameters $\beta\simeq3.925$, $\delta$ phenomenological, and $B \gtrsim 10^{24}$ (Ge et al., 2019). Survival analysis shows this spectrum remains largely intact to present day, providing up to $\Omega_\mathrm{dark}\sim\Omega_\mathrm{visible}$ without fine-tuning (Ge et al., 2019). Constraints from direct-detection, astrophysical, and solar bounds localize the AQN mass window and abundance.

Similarly, "rotating drops" are non-relativistic axion field configurations stabilized by self-gravity, gradient pressure, and self-interaction, with mass up to $M_\mathrm{max}^{(\ell)} \sim 10^{-12}M_\odot$ for high angular momentum (Davidson et al., 2016). The spatial and local velocity distribution of such drops mirrors the underlying macrohalo potential, and they are subject to microlensing, femtolensing, and stellar stream-heating constraints.

5. Multi-Axion and Stochastic Landscape Distributions

In string theory-inspired “axiverse” scenarios, the axion DM distribution is a statistical ensemble of many fields. Each axion’s mass and decay constant are modeled as random variables, often with log-flat or Wishart-distributed eigenvalues. The cumulative dark matter density then depends on a handful of “hyperparameters” (e.g., mean mass scale, variance, spectral width), which can be constrained by Bayesian inference with CMB and large-scale-structure data (Stott et al., 2017).

Specifically,

$$\Omega_a h^2 \sim \sum_i f_{a,i}^2 m_{a,i}^{1/2} \theta_i^2$$

where $f_{a,i}$ and $m_{a,i}$ are decay constants and masses, $\theta_i$ are initial misalignments, and the sum runs over all $n_\mathrm{ax}\sim10$–$100$ fields. The probability of producing the observed $\Omega_\mathrm{DM}$ for given model priors quantifies the likelihood of specific compactification scenarios (Stott et al., 2017).

In stochastic scenarios during inflation, the initial misalignment angle distribution is governed by the Bunch–Davies equilibrium of the Fokker–Planck equation. For a quadratic potential,

$$P_\mathrm{stat}(\theta) \propto \exp\left[-\frac{4\pi^2 m^2 f^2}{3H_I^4} \theta^2\right]$$

leading to today's axion density being a function of the inflationary Hubble scale $H_I$, axion mass $m$, and decay constant $f$ (Nakagawa et al., 2020). In more general “flat-bottomed” potentials, suppressed curvature delays oscillations and allows for a broader parameter regime.

6. Observational and Experimental Probes

The phase-space distribution of axion DM—spanning smooth backgrounds, fine-grained streams, and nonlinear clumps—enters directly into the expected signal in laboratory searches. Haloscope experiments must account for annual modulation and substructure via the likelihood formalism based on $f(v)$ (Foster et al., 2017, O'Hare et al., 18 Sep 2025). High-frequency resolution enables the detection of discrete phase-space spikes, and the presence of fine-grained substructure can significantly boost sensitivity for sufficiently long integration (O'Hare et al., 18 Sep 2025). Detection strategies are sensitive to the lineshape (PSD template) of the axion signal, requiring joint analysis over $g_{a\gamma\gamma}$, $m_a$, and $f(v)$ parameters (Foster et al., 2017).

Outside the laboratory, constraints and possible discovery channels include:

• Microlensing and femtolensing: Compact axion clumps and nuggets are constrained by non-detections in optical and gamma-ray burst surveys (Barranco et al., 2012, Davidson et al., 2016).
• Stellar kinematics and rotation curves: Caustic ring structures are potentially observable as sharp features in galactic rotation curves and local stellar kinematics (Sikivie, 2010, Banik et al., 2013).
• Astrophysical transient signals: Axion–photon conversion in neutron-star magnetospheres and the passage of drops through strong magnetic fields are possible, albeit challenging, indirect detection channels (Barranco et al., 2012, Davidson et al., 2016).

Systematic uncertainties originate from the assumed background and substructure density ($\rho_\mathrm{DM}$, $x$, $x_\mathrm{dd}$), escape speed, Galoactic potential, and primordial axion velocity dispersion. Accurate modeling of axion DM’s spatial and velocity distribution remains a critical input for next-generation axion discovery and exclusion analyses.

7. Summary Table: Principal Axion DM Distributions

Distribution Type	Physics Origin	Spatial Profile	Velocity Structure	Observational Signature
SHM (Maxwellian)	Virialized halo, coll. regime	NFW/Isothermal on $\gtrsim$kpc scales	$\sim$Gaussian ($\sigma_v\sim150$–220 km/s)	Broad lineshape
Fine-grained Streams	Phase-space sheet folding	Overlapping streams, $\sim10^{14}$	Narrow spikes, $\Delta v\ll 1$ km/s	Power-spectrum spikes, modulation
Caustic Rings	BEC, tidal torque, self-similar	Ring caustics at $a_n \propto 1/n$	Rigid-rotation on turnaround sphere	Rotation curve bumps, IRAS features
Clumps/Solitons	Minicluster/soliton collapse	Compact, $R\sim10$–$10^3$ m	Bound-state structure	Microlensing, transient events
Axion Quark Nuggets	QCD domain-wall collapse	Gramscale, $R\sim10^{-5}$ cm	Clustered as cold DM	Direct detection, solar/lithium
Multi-axion ensemble	Axiverse/statistical ensemble	Sum of individual $f_{a,i}, m_{a,i}$	Mixture-of-distributions	Relic density constraints

The distribution of axion dark matter is thus a multi-scale, compositionally complex, and dynamically evolving structure. Its predictive modeling relies on the synthesis of cosmological initial conditions, nonlinear gravitational evolution, particle physics inputs, and high-precision observational and experimental data (Foster et al., 2017, O'Hare et al., 18 Sep 2025, Barranco et al., 2012, Sikivie, 2010, Banik et al., 2013, Vicens et al., 2018, Ge et al., 2019, Davidson et al., 2016, Stott et al., 2017, Nakagawa et al., 2020, Schiappacasse et al., 2017).