Axion White-Noise Fluctuations

• Axion white-noise fluctuations are scale-independent, stochastic density variations from quantum or cosmological origins in axion dark matter scenarios.
• They exhibit distinctive power spectra with sharp cutoffs and influence minihalo formation, galactic dynamics, and atomic clock precision.
• Observational tests, including galactic heating and microlensing, provide key constraints and potential signatures to differentiate axion dark matter from standard models.

Axion white-noise fluctuations denote the stochastic, scale-independent density inhomogeneities that arise in axion dark-matter scenarios from quantum or cosmological origins. These fluctuations play central roles in the astrophysical impact, detection signatures, and theoretical constraints of both "fuzzy" ultra-light axion dark matter and the post-inflation QCD axion scenario. Across approaches, these fluctuations are described by distinctive power spectra, diffusion coefficients, and observational imprints, underlining their significance in galactic dynamics, minihalo formation, and precision experiments.

1. Physical Origin and Theoretical Context

White-noise fluctuations in axion models emerge from fundamentally different mechanisms depending on the mass regime and cosmological history. In post-inflation scenarios (PQ symmetry breaking after inflation), the axion field’s initial misalignment angle is random and uncorrelated between different causally separated Hubble patches. Upon the onset of mass ($m_a$) and coherent field oscillations at temperature $T_0$, this randomness becomes imprinted as order-unity density fluctuations on the corresponding horizon scale $k_0$, resulting in a scale-independent (white-noise) isocurvature power spectrum up to a sharp cutoff (Dai et al., 2019, Xiao et al., 2021):

$$P_\mathrm{iso}(k, t_0) = \Theta(k_0 - k) \frac{24\pi^2}{5k_0^3}\,,$$

where $k_0 = a(t_0)H(t_0)/c$ and $\Theta$ is the Heaviside step function.

For ultra-light axions ("fuzzy dark matter" or FDM), white-noise-like density fluctuations arise from finite de Broglie wavelength effects in a homogeneous quantum field, governed by the free Schrödinger equation:

$$i\hbar\partial_t\psi = -\frac{\hbar^2}{2m}\nabla^2\psi$$

On scales much larger than the de Broglie length $\lambda_\mathrm{dB} = \hbar/(m\sigma)$ (with $\sigma$ being the velocity dispersion), the density power spectrum becomes flat, transitioning to an exponential cutoff at $k \gtrsim k_\lambda = 2\pi/\lambda_\mathrm{dB}$ (El-Zant et al., 2019):

$$P_\delta(k) = P_0 \exp\left(-\frac{k^2 \lambda_\mathrm{dB}^2}{4\pi^2}\right),\quad P_0 = \left(\frac{2\sqrt{\pi}}{m_\hbar\sigma}\right)^3$$

The scale-independence at low $k$ and sharp suppression at the smallest scales are the defining features of the axion white-noise regime.

2. Mathematical Properties: Power Spectrum and Correlations

The hallmark of axion white-noise fluctuations is the flat density power spectrum up to a cutoff. For post-inflationary axions, this takes the form:

$$P(k) \simeq \mathrm{const}\ \text{for}\ k \ll k_0;\quad P(k) \to 0\ \text{for}\ k \gg k_0$$

With density contrast $\delta(\mathbf r, t) = (|\psi|^2/\rho_0) - 1$, the spatial two-point correlation function $\xi(r)$ is the Fourier transform of $P_\delta(k)$. In the fuzzy dark matter case (El-Zant et al., 2019), this yields:

$$\xi(r) = \exp\left(-\frac{r^2}{\lambda_\mathrm{dB}^2}\right)$$

This exponential decay demonstrates that on scales $r \gtrsim \lambda_\mathrm{dB}$, the density fluctuations are uncorrelated, characteristic of spatially white noise, while on smaller scales they are smoothed out by quantum pressure or coherence.

For minihalo-forming scenarios, the variance on mass scale $M$ is given by integration over the power spectrum:

$$\sigma^2(M,z) = \int \frac{d\ln k}{2\pi^2} k^3 P_m(k,z) |W(kR)|^2$$

where $W(kR)$ is a smoothing window function.

3. Nonlinear Evolution: Formation of Minihalos and Substructure

White-noise axion fluctuations drive early, nonlinear collapse in the cosmic density field, resulting in the formation of ultra-compact minihalos. The characteristic mass is set by the horizon crossing scale at the onset of axion oscillations:

$$M_0 = \frac{4\pi}{3} \left(\frac{\pi}{k_0}\right)^3 \bar\rho_{a0}$$

yielding $M_0 \sim 10^{-10} M_\odot$ for typical QCD axion parameters (Dai et al., 2019, Xiao et al., 2021). The minimum halo mass is numerically found to be $M_\mathrm{min} \sim 0.01 M_0$.

N-body simulations confirm that these minihalos rapidly virialize at $z \sim 50-200$, adopting universal Navarro-Frenk-White (NFW) density profiles with extreme concentrations ($c \sim 10^3$–$10^4$) (Xiao et al., 2021):

$\rho(r) = \frac{\rho_s}{(r/r_s)(1 + r/r_s)^2}$

with mass-dependent scale radii and concentrations.

The abundance of such minihalos follows a Press-Schechter or Sheth-Tormen formalism, calibrated to simulations:

$$\nu f_\mathrm{ST}(\nu) = A[1 + (q\nu)^{-p}] \sqrt{\frac{q\nu}{2\pi}} e^{-q\nu/2}$$

with specific parameter fits reproducing the simulation results for the white-noise spectrum.

4. Stochastic Effects in Stellar and Atomic Dynamics

On galactic scales, white-noise axion fluctuations induce stochastic forces on embedded test particles (e.g., stars in disks) (El-Zant et al., 2019). The power spectrum of the potential fluctuations, via Poisson's equation, is:

$P_\Phi(k) = (4\pi G \rho_0)^2 k^{-4} P_\delta(k)$

leading to a force power spectrum $P_F(k) = k^2 P_\Phi(k)$. The velocity variance acquired by a test star traversing these fluctuations grows diffusively:

$$\langle (\Delta v)^2 \rangle = 8\pi G^2 \rho_0 m_\mathrm{eff} \ln\Lambda \int \frac{f_\mathrm{eff}(\mathbf v) d^3 v}{|\mathbf{v}_p - \mathbf{v}|} T$$

where $m_\mathrm{eff}$ encodes the effective mass arising from interference:

$$m_\mathrm{eff} = \left(\frac{4\pi}{m_\hbar}\right)^3 \frac{\int f^2(\mathbf v) d^3 v}{\rho_0}$$

and $\ln\Lambda$ is a Coulomb logarithm regulated by the de Broglie wavelength. The resulting relaxation time,

$$t_r \simeq \frac{v^3}{8\pi G^2 \rho_0 m_\mathrm{eff} \ln\Lambda}$$

mirrors Chandrasekhar's two-body relaxation, but crucially replaces the particle mass $m$ with $m_\mathrm{eff}$ controlled by quantum statistics.

White-noise axion field fluctuations similarly impact atomic frequencies. Treating the axion as a classical field with coherence time $\tau_\mathrm{coh}$, the observable shifts in clock or energy transitions are characterized by their power spectral density (PSD), autocorrelation, and amplitude distributions (Flambaum et al., 2023). The fractional frequency shift autocorrelation

$$R_x(\tau) = \langle x(t+\tau) x(t) \rangle \simeq \sigma_x^2 e^{-|\tau|/\tau_\mathrm{coh}}$$

yields a flat (white) one-sided PSD over frequencies $\omega \ll 1/\tau_\mathrm{coh}$.

5. Observational Implications and Constraints

Astrophysical and precision measurement implications of axion white-noise fluctuations are multifaceted.

Astrophysical Constraints

• Galactic disk heating: Assuming that observed vertical or radial velocity dispersions in stellar disks (e.g., Milky Way) arise entirely from FDM-induced fluctuations sets a lower bound:

$m_a \gtrsim 2 \times 10^{-22}\,\mathrm{eV}$

based on the empirical requirement that FDM-induced heating not overshoot observations (El-Zant et al., 2019).

• Ultra-compact cluster survival: In the Eridanus II central cluster, diffusion from axion fluctuations constrains $m_a \gtrsim 8.8 \times 10^{-20}$ eV for cluster expansion over 3 Gyr, although this regime may lie inside the FDM soliton core and requires careful modeling of coherence and locality (El-Zant et al., 2019).
• Minihalo population: The high area-covering fraction ($\gg 1$) of minihalos in clusters renders the projected surface mass density a quasi-Gaussian field with rms fluctuations $\Delta_\kappa \sim 10^{-4} - 10^{-3}$ on au–$10^4$ au scales (Dai et al., 2019, Xiao et al., 2021).
• Microlensing caustics: During microlensing of highly magnified stars in cluster lenses, convergence fluctuations at this amplitude are comparable to $1/\mu$ for $\mu \sim 10^3$–$10^4$, directly producing $O(1)$ irregularities in lightcurves. These features provide a distinctive signature of axion white-noise structure in gravitational lensing events (Dai et al., 2019).

Precision Laboratory Constraints

• Atomic clock comparisons: Stochastic axion-induced frequency fluctuations in atomic clocks are detectable via broadband variance or temporal autocorrelations (Flambaum et al., 2023). Observations with Rb/Cs and H/Si hyperfine clock networks have set limits:

$$f_a \gtrsim 1.8 \times 10^9\,\mathrm{GeV}\cdot(10^{-15}\,\mathrm{eV}/m)\ \text{[Rb/Cs,}\ 2.4\times 10^{-17}~\mathrm{eV} < m < 2.4\times 10^{-12}~\mathrm{eV]}$$

and

$$f_a \gtrsim 5.8 \times 10^9\,\mathrm{GeV}\cdot(10^{-15}\,\mathrm{eV}/m)\ \text{[H/Si,}\ 7.3\times 10^{-16}~\mathrm{eV} < m < 1.9\times 10^{-10}~\mathrm{eV]}$$

filling in otherwise unconstrained windows between astrophysical and laboratory axion limits.

• Frequency-distribution diagnostics: Binning and higher-moment analysis of clock time series facilitate extraction of the axion coherence time $\tau_\mathrm{coh}$ and hence $m$ via features such as skewness~2 and kurtosis~9 in the ideal limit.

6. Simulation and Mass Function Calibration

N-body simulations incorporating the white-noise initial conditions validate the theoretical expectation of early, ultra-compact minihalo formation (Xiao et al., 2021). The relaxed minihalo population at $z\sim 19$–100 closely follows the Sheth–Tormen form, with mass function and concentration scaling laws robust for a wide range of axion parameters:

$$M_0 \sim 10^{-10} - 10^{-12} M_\odot,\quad c(M, z) \sim 10^3 - 10^4$$

Calibration accommodates uncertainties due to axion strings, modeled via a parameter $A_\mathrm{osc}$ controlling the amplitude and cutoff behavior. The mass function delayed after $z \sim 30$ is frozen due to minihalo assimilation into larger structures, but a substantial surviving population is expected, rendering astrophysical observations sensitive to these features.

7. Significance and Limitations

Axion white-noise fluctuations constitute a robust, predictive consequence of post-inflation axion cosmology and are fundamental to the phenomenology of fuzzy dark matter. Their unique power spectrum, correlation structure, and resulting minihalo population distinguish them from adiabatic fluctuations of standard cold dark matter.

Current constraints from galactic heating, cluster dynamics, and atomic-frequency stability collectively bound axion parameter space, sometimes up to several orders of magnitude above the QCD axion line, but all methodologies face intrinsic model-dependent and technical limitations. For example, the white-noise stochastic framework becomes unreliable in regimes dominated by global soliton coherence or when stochasticity is subdominant to systematic core oscillations (El-Zant et al., 2019). Uncertainties in axion string-induced corrections and non-Gaussianity introduce order-unity ambiguities in the cutoff mass and abundance, but do not alter qualitative predictions.

A plausible implication is that upcoming high-sensitivity microlensing, pulsar timing, and atomic clock experiments will be able to probe, detect, or rule out swathes of the axion parameter space by exploiting the stochastic and scale-invariant features of axion white-noise fluctuations, potentially distinguishing axion dark matter from standard cold dark matter via unique signatures in both astrophysical and laboratory domains.