--- title: Ultralight Scalar Dark Matter Overview url: https://www.emergentmind.com/topics/ultralight-scalar-dark-matter-ulsdm type: topic --- # Ultralight Scalar Dark Matter Overview Ultralight scalar dark matter (ULSDM) denotes bosonic dark matter candidates with masses far below the eV scale, often described as coherently oscillating classical fields because their occupation numbers are enormous and their de Broglie wavelengths are macroscopic. In the literature, the term covers both “fuzzy” dark-matter scales around \(m_\phi \sim 10^{-22}\,\mathrm{eV}\) and higher ultralight masses extending to the \(\mu\mathrm{eV}\)–\(\mathrm{meV}\) regime in specific portal constructions. Its phenomenology is correspondingly broad: the scalar may be generated by standard misalignment, thermal or environmental mass effects, or non-minimal gravitational couplings; it may couple linearly or quadratically to Standard Model operators, universally through an effective metric, to neutrinos, or through flavor-violating portals; and it can be constrained by Big Bang Nucleosynthesis (BBN), compact objects, pulsar timing, clocks, atom interferometers, neutrino experiments, and other precision probes [2408.12667] [2211.05174] [2306.13799]. ## 1. Coherent-field description and kinematic regime A standard starting point is a real scalar field \(\phi\) with ultralight mass \(m_\phi\), treated as a classical background because of its large occupation number. In halo applications the field is written as \[ \phi(t,\mathbf{x}) \simeq \phi_0 \cos(\omega t - \mathbf{k}\!\cdot\!\mathbf{x} + \theta),\qquad \phi_0=\sqrt{\frac{2\rho_{\rm DM}}{m_\phi^2}}, \] with \(\omega \simeq m_\phi\) and \(k \simeq m_\phi v\) for a nonrelativistic halo. The coherence scales are \[ \lambda_{\rm dB}=\frac{2\pi}{m_\phi v},\qquad \tau_c\sim \frac{2\pi}{m_\phi v^2}, \] and for virialized dark matter \(v\simeq 10^{-3}c\) the quality factor is \(Q\approx v^{-2}\approx 10^6\) [2601.16259] [2110.11564]. Once oscillations begin, the scalar behaves as pressureless matter. In the bare-mass regime, \[ \phi=\Phi(t)\cos(m_\phi t+\varphi),\qquad \Phi\propto a^{-3/2},\qquad \rho_\phi=\frac{m_\phi^2\Phi^2}{2}\propto a^{-3}, \] and in WKB regimes with time-dependent effective mass one instead has \(\phi_{\rm amp}\propto m_{\rm eff}^{-1/2}a^{-3/2}\), but still \(\rho_\phi\propto a^{-3}\) after oscillations start [2006.04820] [2211.09826]. The classical-field description is not the only useful language. For high-energy neutrino propagation through an ultralight background with \(m_\phi \approx 10^{-23}\,\mathrm{eV}\), ULSDM has also been modeled as a cold, nonrelativistic particle bath with \(n_\phi \approx \rho_{\rm DM}/m_\phi \approx 3\times 10^{31}\,\mathrm{cm}^{-3}\), because the relevant observable is scattering with absorption and regeneration rather than coherent forward scattering [1605.09671]. A common misconception is therefore that ULSDM must always be treated in one universal way; in practice, the appropriate description depends on the observable. ## 2. Effective couplings and portal structures A prominent class of models uses universal couplings that preserve the weak equivalence principle by coupling \(\phi\) to Standard Model matter through an effective metric, \[ \bar g_{\mu\nu}=g_{\mu\nu}\bigl(1+2\alpha(\phi)\bigr),\qquad \alpha(\phi)\simeq \pm \phi^2/\Lambda^2. \] In the Einstein frame the leading interaction is \(\mathcal{L}_{\rm int}=-\alpha(\phi)\Theta_{\rm SM}\), so fermion and vector masses rescale as \(m_{f,V}\mapsto m_{f,V}(1+\alpha)\), while WEP preservation also forces \(\Lambda_{\rm QCD}\mapsto \Lambda_{\rm QCD}(1+\alpha)\) and \(M_{\rm pl}^2\mapsto M_{\rm pl}^2(1-2\alpha)\) in the Jordan frame. Physical observables are frame-invariant even though practical implementations differ [2006.04820]. A second broad EFT parameterization employs Damour–Donoghue-type linear and quadratic couplings. For quadratic interactions one writes \[ \mathcal{L}\supset 2\pi \frac{\phi^2}{M_{\rm pl}^2}\left[\frac{d_e^{(2)}}{4e^2}F_{\mu\nu}F^{\mu\nu} -\frac{d_g^{(2)}\beta_3}{2g_3}G^A_{\mu\nu}G^{A\mu\nu} -d_{m_e}^{(2)}m_e\bar ee-\sum_{i=u,d}(d_{m_i}^{(2)}+\gamma_{m_i}d_g^{(2)})m_i\bar\psi_i\psi_i\right], \] with the induced fractional shifts \[ \frac{\delta\alpha}{\alpha}=d_e^{(2)}\frac{\varphi^2}{2},\qquad \frac{\delta\Lambda_{\rm QCD}}{\Lambda_{\rm QCD}}=d_g^{(2)}\frac{\varphi^2}{2},\qquad \frac{\delta m_f}{m_f}=d_{m_f}^{(2)}\frac{\varphi^2}{2}. \] The general phenomenology of linear and quadratic \(d_i^{(1,2)}\) couplings includes oscillations at \(m_\phi\) and \(2m_\phi\), DC offsets for quadratic couplings, and source-dependent screening effects near massive bodies [2211.09826] [2211.05174]. ULSDM portals are not limited to diagonal “dilaton-like” deformations. Neutrino-specific couplings can modulate neutrino masses and mixings through \[ {\cal L}_{\rm eff}=-m_\nu\left(1+y\frac{\phi}{\Lambda}\right)\nu\nu+{\rm h.c.}, \] while right-handed-neutrino portals use \[ \mathcal{L}\supset -(g\phi-m_N)n^cn^c-yLHn^c+{\rm h.c.} \] and make the scalar evolution sensitive to the finite-density neutrino background [1705.06740] [2409.17396]. Flavor-violating constructions provide another qualitatively distinct possibility, \[ \mathcal{L}_{\rm int}=\phi\sum_{i\neq j}\lambda_{ij}\,\bar d_i d_j, \] which induces \(\mathcal{O}(\phi)\) CKM modulation and \(\mathcal{O}(\phi^2)\) down-quark mass shifts [2603.17237]. Gauge-kinetic couplings, \[ I(\phi)F_{\mu\nu}F^{\mu\nu},\qquad I(\phi)=\exp(2\beta \phi/M), \] open yet another sector of ULSDM phenomenology, including late-time magnetogenesis and oscillatory variation of \(\alpha_{\rm eff}=\alpha_0/I\) [2601.00443]. ## 3. Production mechanisms and cosmological evolution The canonical production mechanism is misalignment: a light scalar displaced from its minimum begins oscillating once \(H\) drops below the relevant mass scale and then redshifts like cold matter. In models with quadratic couplings this simple picture is modified because the scalar acquires a finite-temperature or finite-density effective mass, \[ m_{\rm eff}^2(T)=m_\phi^2+\Pi(T), \] and the onset of oscillations satisfies \(H(T_{\rm osc})\simeq m_{\rm eff}(T_{\rm osc})\). For electron couplings, \[ m_{\rm eff}^2(T)=m_\phi^2+\frac{4\pi d_{m_e}^{(2)}}{M_{\rm pl}^2}(\rho_e-3P_e), \] while for photon couplings the leading high-temperature contribution is two-loop and scales as \(\Pi(T)\propto T^4\) [2211.09826]. In universally coupled models the environmental mass is governed by the Standard Model trace, \[ m_{\phi,{\rm eff}}^2=m_\phi^2\pm 2\frac{\Theta_{\rm SM}}{\Lambda^2},\qquad \Theta_{\rm SM}=\rho_{\rm SM}-3p_{\rm SM}. \] During radiation domination relativistic species contribute negligibly, but around BBN the \(e^\pm\) plasma becomes nonrelativistic and generates a large \(\Theta_e\). This produces three characteristic regimes: Hubble friction, a bare-mass oscillatory regime, and an induced-mass regime. For \(\alpha>0\), induced-mass domination can move the onset of oscillations earlier and suppress \(\phi\) by BBN; for \(\alpha<0\), the induced mass becomes tachyonic and can trigger exponential growth unless the initial conditions are extremely fine-tuned [2006.04820]. Several works study departures from ordinary misalignment. “Gravitational misalignment” introduces a non-minimal curvature coupling, \[ I_\phi=\int d^4x\,\frac{\sqrt{-g}}{2}\left[\partial^\mu\phi\partial_\mu\phi-(m_\phi^2+\xi{\cal R})\phi^2-\frac{\lambda}{2}\phi^4\right], \] so that during radiation domination trace-anomaly effects make \({\cal R}\neq 0\), and \(\xi>0\) can generate a tachyonic phase with \[ \langle\phi\rangle=\sqrt{\frac{|\mu_\phi^2|}{\lambda}},\qquad \mu_\phi^2\equiv \xi{\cal R}+m_\phi^2. \] For the benchmark \(\xi=10^3\), \(\epsilon=0.1\), \(g_*=10.75\), and \(m_\phi=10^{-11}\,\mathrm{eV}\), one finds \(T_*\approx 50\,\mathrm{MeV}\), \(\lambda\approx 1.3\times 10^{-68}\), and \(\phi_*\approx 4.6\times 10^{13}\,\mathrm{GeV}\) [2408.12667]. A further departure arises when ULSDM couples to right-handed neutrinos. In that case the relic neutrino density generates an asymmetric finite-density effective potential, producing a “brick-wall” regime in which the scalar energy density redshifts nontrivially and can undergo an \(\mathcal{O}(1)\) drop between recombination and today. This implies that CMB neutrino-mass bounds and laboratory oscillation bounds constrain different regions of parameter space rather than the same naive misalignment amplitude [2409.17396]. ## 4. Big Bang Nucleosynthesis and early-universe limits BBN provides one of the most developed cosmological probes of ULSDM because scalar backgrounds can shift the effective values of Standard Model masses and couplings during neutron freeze-out and the deuterium bottleneck. In universally coupled models, the relevant variations include \[ \frac{\Delta G_F}{G_F}=-2\alpha(\phi),\qquad \frac{\Delta m_{np}}{m_{np}}\simeq \alpha(\phi),\qquad \frac{\Delta \Gamma_n}{\Gamma_n}\simeq \alpha(\phi),\qquad \frac{\Delta B_D}{B_D}\simeq \alpha(\phi). \] The helium-4 mass fraction depends mainly on the neutron-to-proton ratio at freeze-out and on neutron decay before nucleosynthesis, \[ Y_p\simeq 2X_n(t_{\rm nuc}), \] and the observationally allowed fractional deviation is \[ \frac{\Delta Y_p}{Y_p}=-0.008\pm 0.012. \] A full kinetic treatment shows that precision helium-4 measurements set the dominant constraint, deuterium is subdominant, and the lithium effect is insufficient to solve the cosmological lithium problem [2006.04820]. The same conclusion persists in more general quadratic-coupling analyses, but with important refinements. Once the finite-temperature thermal mass and the full weak freeze-out kinetics are included, BBN constrains \(m_\phi\sim 10^{-19}\,\mathrm{eV}\) up to \(\sim 10^{-4}\,\mathrm{eV}\), with qualitatively different scaling regimes depending on whether the field is Hubble-frozen, thermally-induced-mass dominated, or bare-mass dominated during freeze-out. The inclusion of thermal mass and kinetic freeze-out weakens the electron-coupling bound by about two orders of magnitude at \(m_\phi \gtrsim 10^{-14}\,\mathrm{eV}\), strengthens the quark-coupling bound by a factor of about 2 in the same range, and yields new gluon-coupling constraints [2211.09826]. A central technical result is that frame invariance is preserved only if freeze-out is handled beyond the instantaneous approximation. In the Einstein frame ULSDM appears as time-dependent masses and widths, whereas in the Jordan frame it appears as a modified Hubble rate, \[ \bar H(\bar a)=H^{(0)}(\bar a)\left[1+\alpha(\bar a)\left(1+\frac{d\ln\alpha}{d\ln\bar a}\right)\right]. \] The numerical implementation in AlterBBN v2.2 therefore proceeds through a modified “dark density” \[ \rho_{\rm dark}=6M_{\rm pl}^2[H^{(0)}(\bar a)]^2\alpha(\bar a)\left(1+\frac{d\ln\alpha}{d\ln\bar a}\right), \] with conservative cuts imposed once the linearized treatment begins to fail [2006.04820]. The sign of the coupling matters sharply. For \(\alpha=+\phi^2/\Lambda^2\), BBN bounds are strongest at very low \(m_\phi\), weaken for \(m_\phi\gtrsim 10^{-18}\,\mathrm{eV}\) because induced-mass suppression lowers \(\alpha\) at BBN, and develop oscillatory exclusion bands at intermediate masses. For \(\alpha=-\phi^2/\Lambda^2\), a large tachyonic region during BBN requires extreme fine-tuning of the initial field value to avoid dark-matter overproduction and is effectively excluded; outside that region, the BBN effect is negligible, though neutron-star scalarization may remain relevant [2006.04820]. ## 5. Observational searches from clocks to pulsar timing Laboratory and space-based frequency metrology probe ULSDM through periodic shifts of atomic transition frequencies, masses, and effective constants. ISS-based clocks are especially relevant for quadratically coupled ULSDM because atmospheric shielding suppresses the ground-level field for sufficiently strong couplings, whereas at orbital altitudes the orbit-averaged monopole coefficient approaches \(a_0\to 1\) for \(m_\phi\gtrsim 10^{-9}\,\mathrm{eV}\). When \(\lambda_{\rm dB}