Ultralight Scalar Dark Matter Overview

• Ultralight Scalar Dark Matter is a class of models employing sub-eV bosonic fields that exhibit macroscopic quantum phenomena, including solitonic cores and interference patterns in galactic halos.
• These models are governed by the Schrödinger–Poisson equations and encompass variants like Fuzzy Dark Matter, Self-Interacting FDM, and axion-like scenarios with distinct astrophysical signatures.
• Observational and experimental probes, from galactic rotation curves to precision interferometry, tightly constrain ULDM parameters while offering new detection avenues for its unique quantum effects.

Ultralight Scalar Dark Matter (ULDM) denotes a class of dark matter models in which the dark-sector component consists of bosonic fields with masses far below the eV scale, typically $$10^{-24}~\mathrm{eV} \lesssim m \lesssim 1~\mathrm{eV}$$. These models are motivated by both theoretical considerations, such as the ubiquity of light scalar fields in string-theoretic compactifications and axion-like scenarios, and the phenomenological tension between standard cold dark matter (CDM) and observed small-scale galactic structure. ULDM exhibits macroscopic quantum behavior on scales that can reach kiloparsecs, manifesting in wave-like phenomena such as solitonic cores, interference patterns, and a granular structure in the cores of galactic dark matter halos.

1. Model Structure and Fundamental Dynamics

At the field-theoretic level, the prototypical ULDM scenario features a real or complex scalar field $\phi$ minimally coupled to gravity, with a canonical or weakly self-interacting potential. The nonrelativistic regime, relevant for halo dynamics, is governed by the Schrödinger–Poisson equations: $$i\hbar\,\partial_t\psi = -\frac{\hbar^2}{2m}\nabla^2\psi + m\,\Phi\,\psi,\qquad \nabla^2\Phi = 4\pi G\,m\,|\psi|^2,$$ where $\psi$ is the field amplitude, $m$ the ULDM mass, and $\Phi$ the Newtonian gravitational potential (Ferreira, 2020). The interplay of "quantum pressure," gravity, and possible self-interactions regulates the formation and structure of halos.

Key classes within the scalar ULDM paradigm include:

• Fuzzy Dark Matter (FDM): free scalar field, $m\sim10^{-22}$ eV, forming Bose–Einstein condensates and solitonic cores.
• Self-Interacting FDM (SIFDM): quartic self-coupling, modifying core size and stability.
• Axion-like ULDM: pseudo-Nambu–Goldstone bosons, possibly with nontrivial cosmological production and a periodically modulated potential (Delaunay et al., 16 Jul 2025).

The characteristic de Broglie wavelength,

$\lambda_{\mathrm{dB}} = h/(mv),$

can reach $\mathcal{O}(\mathrm{kpc})$ for $m\sim10^{-22}$ eV and $v\sim100\,\mathrm{km\,s^{-1}}$, leading to coherent quantum phenomena on galactic scales.

2. Cosmological Production, Redshift, and Abundance

Scalar ULDM can be produced through several mechanisms:

• Vacuum misalignment (standard for axions): field starts displaced from the minimum of its potential in the early universe, oscillates as $H<m$ and redshifts as matter.
• Non-adiabatic cosmological production: for a minimally coupled, free ultra-light scalar in its Bunch-Davies vacuum during inflation followed by instantaneous reheating, the resulting momentum distribution is IR-enhanced, $\mathcal{N}_k\propto 1/k^3$. This process yields a cold, nonthermal spectrum with an equation of state $w\sim10^{-14}$. The observed dark matter abundance can be saturated for $m\sim 1.5\times10^{-5}$ eV at the inflationary energy scale limit (Herring et al., 2019).
• Coupled scenarios: Interactions with other sectors (e.g., right-handed neutrinos) can induce an asymmetric effective potential, yielding nonstandard redshift evolution and imposing constraints from cosmology and lab probes due to time-varying masses (Plestid et al., 2024).

The free-streaming length is negligible except for extremely small masses; e.g., $\lambda_{\textrm{fs}}\sim 70$ pc for $m=1.5\times10^{-5}$ eV. Thus, scalar ULDM with $m\gtrsim 10^{-24}$ eV is "cold" for structure formation.

3. Galactic Halo Structure and Small-Scale Phenomenology

The defining prediction of scalar ULDM is the formation of soliton-like cores at the centers of halos. These exhibit density profiles (Ferreira, 2020, Chan et al., 2021): $$\rho_\mathrm{core}(r) = 1.9 \times 10^9\,M_\odot\,\mathrm{kpc}^{-3}\,\big[(10\,m_{22})^{-2}\, r_c^{-4}\big][1+9.1\times10^{-2}(r/r_c)^2]^{-8},$$ with $m_{22}\equiv m/(10^{-22}\,\mathrm{eV})$. The core size and mass scale with the halo properties as

$$r_c \sim 1.6\, m_{22}^{-1} (M_h/10^9 M_\odot)^{-1/3}\,\mathrm{kpc},\quad M_c \sim 3.1\times10^4 (M_h/M_\odot)^{1/3} m_{22}^{-1} M_\odot,$$

resulting in cored profiles inconsistent with the central cusps of standard CDM.

ULDM's wave nature also leads to:

• Granular or "interference" structure in the halo due to overlapping eigenstates.
• Stochastic heating of stars by density "granules" whose typical scale is set by $\lambda_{\mathrm{dB}}$ (Gosenca et al., 2023).
• Suppression of sub-galactic structure (cut-off in the halo mass function) determined by the linear power-spectrum suppression at the Jeans scale $k_J\sim m^{1/2}$.

4. Laboratory, Astrophysical, and Cosmological Constraints

A range of precision experiments and astrophysical measurements constrain or probe the scalar ULDM parameter space:

• Astrophysical constraints:
  • Stellar kinematics and rotation curves allow model-independent exclusions in the mass range $0.14-3.11\times10^{-22}$ eV for all-DM scenarios (Chan et al., 2021).
  • Limits from central core masses in massive galaxies (e.g., M87 provide $m\lesssim6\times10^{-23}\,$eV for $\lambda=0$; the constraint is sensitive to even extremely weak self-coupling and can reach axion-like couplings $\lambda\sim-10^{-96}$ (Chakrabarti et al., 2022).
  • Pulsar timing around Sgr A* could probe down to $\sim$solar-mass clouds or solitons for $m\sim10^{-18}$ eV (Yu et al., 26 Oct 2025).
• Cosmological probes:
  • The CMB and baryon acoustic oscillations are sensitive to ULDM-induced time-variation of fundamental constants (via quadratic couplings to $\alpha$ or $m_e$) affecting recombination and BBN (Ghosh et al., 18 Nov 2025). For $m_\phi\lesssim 10^{-26}$\,eV, the allowed DM fraction is constrained to $f_\phi\lesssim 0.3$.
• Direct detection and terrestrial experiments:
  • Atom interferometers and atomic clocks search for time-dependent oscillations of fundamental constants or mass ratios induced by scalar ULDM (Deshpande et al., 2024). Differential cavity-length measurements have set limits on the electron-mass coupling $d_{m_e}\sim10^{-6}$ for the SHM scenario, improving by one to two orders of magnitude over previous cavity-based results.
  • Atom gradiometers (both broadband and compact, e.g. AION-10) are sensitive to linearly coupled scalar ULDM through oscillations in atomic transition energies, with optimal reach depending on baseline, atom number, and cycle time (Badurina et al., 2023, Badurina et al., 2021).
  • Pulsar Timing Arrays (PTAs) constrain both the gravitational and direct-coupling signatures in millisecond pulsar arrival times, with EPTA DR2 constraining $d_\gamma\lesssim10^{-8}$, $d_{\hat m}\lesssim2\times10^{-7}$, and $d_g\lesssim3\times10^{-7}$ at $m_\phi\approx10^{-23}$ eV (Wu et al., 2024).
  • LISA and other space-based laser interferometers discriminate ULDM-induced signals (e.g., Doppler modulation) from monochromatic gravitational waves, with reach to $d_g\sim10^{-6}$ at $m_\phi\sim10^{-17}$ eV (Gué et al., 19 Aug 2025).

5. Multifield and Model Extensions

Theoretically motivated UV completions often introduce not a single field, but an ensemble of $N$ ultralight scalars (multifield ULDM). In such multifield models:

• The total density fluctuation amplitude and halo granulation are suppressed by $1/\sqrt{N}$ compared to the single-field scenario.
• Stellar heating rates are reduced, scaling as $1/N$ (for equal-mass, equal-fraction) or as $1/N^2$ if the lightest field dominates (Gosenca et al., 2023).
• This relaxation of granular structure and stochastic heating relaxes lower mass bounds from observation, e.g., if $N=10$, a $10^{-21}$ eV lower bound reduces to $5\times10^{-22}$ eV.
• Particle production and cosmology are altered if ULDM couples to additional sectors, notably to right-handed neutrinos, leading to temperature-dependent effective potentials, nontrivial redshift behavior, and order-unity variations in both DM density and neutrino masses on cosmological timescales (Plestid et al., 2024).

Furthermore, the quadratic twin mechanism provides a solution to radiative instability of the scalar ULDM mass for quadratically coupled fields. By extending the SM with a mirror sector and a $\mathbb{Z}_2$ symmetry, linear radiative corrections to the ULDM mass cancel, leaving only quadratic corrections. This mechanism renders couplings up to $g\sim10^{-10}-10^{-7}$ natural for $m_\phi\sim10^{-20}-10^{-15}$ eV and vastly enlarges the natural parameter space accessible to current and future experiments (Delaunay et al., 16 Jul 2025).

6. Distinctive Experimental and Astrophysical Signatures

ULDM models with linear or quadratic couplings to SM fields predict a suite of experimentally accessible signatures:

• Periodic modulation of atomic transition frequencies, electron or nucleon masses, and the fine-structure constant at the Compton frequency $f=m/(2\pi)$ (Deshpande et al., 2024, Badurina et al., 2023).
• Stochastic and discrete "granular" density structures at the de Broglie scale, leading to observable heating and velocity dispersion in stars and stellar streams (Gosenca et al., 2023).
• Time-averaged and time-resolved modifications of neutrino oscillation probabilities (including unique CP-violating and CP-conserving signatures) via oscillatory shifts in neutrino mass eigenvalues or mixings (Losada et al., 2023, Losada et al., 2021, Delgadillo et al., 20 Dec 2025).
• Suppression of the nanohertz gravitational-wave background from SMBH binaries due to enhanced dynamical friction in soliton cores, testable with PTAs, which places robust constraints on $m\sim(2.3-19)\times 10^{-22}$ eV if ULDM forms a significant DM fraction (Tiruvaskar et al., 17 Dec 2025).
• Prospective detection of soliton cores or gravitational-atom clouds via precise timing of pulsar orbits around supermassive black holes, probing $M_c\sim1-100\,M_\odot$ for $m\sim10^{-20}$–$10^{-17}$ eV (Yu et al., 26 Oct 2025).

The breadth of observable and experimental signatures provides for a multidimensional test of scalar ULDM across both laboratory and astrophysical scales.

7. Current Status, Constraints, and Future Directions

A convergence of astrophysical, cosmological, and laboratory limits has significantly constrained the parameter space for ULDM as the dominant dark matter component. Model-independent analyses of galactic rotation curves exclude $0.14\le m/(10^{-22}\,\mathrm{eV})\le 3.11$ for the total DM (Chan et al., 2021). CMB and Lyman-$\alpha$ forest data push $m\gtrsim 10^{-21}$ eV (Ferreira, 2020, Ghosh et al., 18 Nov 2025). Nonetheless, multifield constructions, environmental screening, and nontrivial self-coupling can substantially alleviate or shift these bounds.

Direct detection efforts are rapidly advancing, with clock-comparison and interferometer-based experiments probing ever deeper into theoretically natural scalar coupling parameter space, especially given the radiative protection afforded by recent model innovations (Delaunay et al., 16 Jul 2025). The next generation of large liquid scintillator detectors, atom interferometers, and SKA-era PTAs, in conjunction with improved astrophysical modeling and broader mass-coupling reach, will enable the mapping of remaining ULDM parameter space with high precision (Delgadillo et al., 20 Dec 2025, Badurina et al., 2021, Wu et al., 2024).

The scalar ULDM framework thus remains at the forefront of combined theoretical and experimental dark matter research, providing both stringent testability and structural innovation within particle cosmology.