---
title: Sommerfeld Enhancement in Dark Matter Physics
url: https://www.emergentmind.com/topics/sommerfeld-enhancement
type: topic
---

# Sommerfeld Enhancement in Dark Matter Physics

Sommerfeld enhancement is the non-perturbative modification of a short-distance reaction rate by a long-range interaction acting on a non-relativistic two-body state. In dark-matter physics it usually denotes the change in annihilation or coannihilation rates caused by the distortion of the incoming wave function relative to a free plane wave, so that the probability density at short distance is enhanced or suppressed as the relative velocity decreases; related constructions also appear for final-state interactions and in above-gap absorption in two-dimensional Dirac materials [2403.01682][2603.02647][2104.12569].

## 1. Definition, physical content, and scope

The central physical picture is that slowly moving particles spend enough time in a long-range potential for their scattering state to differ substantially from a plane wave. If the potential is attractive, the wave function can be enhanced near the origin and the annihilation probability rises; if the potential is repulsive, the wave function is depleted and the rate is suppressed. Several of the cited works stress that the effect is strongest at low relative velocity, which is why it enters both early-universe freeze-out and late-time indirect detection [2209.14343][2507.12522].

In the most common formulation, the long-distance dynamics modify a hard annihilation process multiplicatively rather than replacing it. A representative expression is
$$
\langle \sigma v\rangle = \sum_l S_l \cdot \langle \sigma v\rangle_{(l)}^{\rm pert},
\qquad
\langle \sigma v\rangle_{(l)}^{\rm pert}\simeq a_l v^{2l},
$$
where the partial-wave dependence of the perturbative rate is preserved and the Sommerfeld factor $S_l$ encodes the wave-function distortion [2512.24662]. For asymmetric dark matter, the same structure is often written as
$$
\langle \sigma v \rangle = a\,\langle S_s\rangle + b\,\langle v^2 S_p\rangle,
$$
making explicit the separate $s$-wave and $p$-wave contributions [2403.01682].

A common misconception is that the effect is always an enhancement and always belongs to the initial state. The literature supplied here shows both statements to be too narrow. Repulsive loop-induced quantum forces can suppress the rate rather than enhance it [2507.12522], and for annihilation into heavier unstable particles the same formalism can act in the final state, where the slow outgoing pair experiences its own long-range distortion near threshold [2603.02647].

## 2. Non-relativistic formulation

The standard treatment is a Schrödinger problem for the two-body relative coordinate. In electroweak dark matter, for example, the non-relativistic dynamics are written as
$$
\left( - \frac{1}{2 \mu}\nabla^2 + V(r) - i \frac{1}{2}\Gamma(r) \right)\psi(r) = E\psi(r),
$$
followed in the usual Sommerfeld treatment by neglect of the imaginary part and solution of the radial equation in partial waves [2606.16442]. For a light force carrier the canonical potential is Yukawa,
$$
V(r) = - \frac{\alpha_X}{r} e^{-m_X r},
$$
or equivalently $V(r)=-\alpha e^{-m_\phi r}/r$ in the scalar-mediator notation [2209.14343][2512.24662].

The partial-wave enhancement factor is extracted from the short-distance behavior of the distorted wave function. One explicit definition used for general orbital angular momentum is
$$
S_l = \left| \frac{(2l+1)!!}{|\mathbf{p}|^l\,l!} \left. \frac{\partial^l R_{El}(r)}{\partial r^l} \right|_{r=0} \right|^2,
$$
which reduces to the familiar wave-function-at-the-origin expression in the $s$-wave channel [2512.24662]. In coupled-channel electroweak problems, the same logic is implemented through matrix Schrödinger equations for nearly degenerate neutral and charged states [2606.16442].

For massless or Coulombic force carriers, analytic expressions are available. In the asymmetric-dark-matter treatment with a massless mediator,
$$
S_s = \frac{2 \pi\alpha/v}{1-e^{-2\pi\alpha/v}},
\qquad
S_p = \left[1 + \left(\frac{\alpha}{v}\right)^2\right] \frac{2 \pi\alpha/v}{1-e^{-2\pi\alpha/v}},
$$
so the enhancement grows strongly at low velocity [2001.08404]. The thermal average is then performed over the non-relativistic velocity distribution, and in the massless force-carrier limit rational-function interpolations reproduce the exact thermal averages with less than $0.3\%$ error for $s$-wave and less than $0.9\%$ error for $p$-wave, allowing the standard analytic freeze-out method to remain accurate to within $1\%$ [1008.2905].

## 3. Interaction structures and mediator sectors

The literature no longer treats Sommerfeld enhancement as a phenomenon tied only to a single tree-level Yukawa potential. The supplied works span single mediators, mediator towers, loop-induced quantum forces, neutrino-exchange forces, finite-density backgrounds, and finite-size composites. The interaction structure determines not only the size of the effect but also whether the rate is enhanced or suppressed, whether resonances occur, and which observables are most sensitive.

| Setting | Interaction structure | Reported feature |
|---|---|---|
| Light vector or scalar mediator | Yukawa potential $-\alpha e^{-m_\phi r}/r$ | Standard low-velocity enhancement [2403.01682] |
| Narrow resonance plus long-range force | Product $S_{\rm SF}(v)\,R(v)$ at leading order | Complete factorization [2209.14343] |
| Tower of mediators | $V(r)=-\sum_n \frac{\alpha_n}{2r}e^{-m_n r}$ | Off-resonant enhancement increased by about $15$–$20\%$ or $\sim 20\%$ in the studied setup [1203.6341] |
| Loop-induced quantum force | Vacuum potentials scaling as $1/r^3$ or $1/r^5$ | Enhancement or suppression depending on operator structure [2507.12522] |
| Two-neutrino exchange | Naive $-1/r^5$ potential with UV-completion dependence | Strong sensitivity to short-distance completion [2203.05455] |
| Finite-density dark background | Density-dependent Yukawa-like effective force | Stronger enhancement in dense environments [2512.24188] |

For mediator towers, the potential becomes a sum of Yukawa contributions,
$$
V(r)=-\sum_n \frac{\alpha_n}{2r}e^{-m_n r},
$$
and in the regime $\epsilon_\phi\ll\epsilon_v$ the combined effect is an effective larger coupling. In the confining-CFT / warped-dual construction analyzed in the cited work, off-resonant enhancement is larger than the standard single-mediator result by about $15$–$20\%$ for $\epsilon_\phi\sim 10^{-3}$ when the UV cutoff is near the Planck scale [1203.6341].

Loop-induced quantum forces are more singular than Yukawa exchange. In the examples studied, two-scalar exchange yields an attractive potential behaving as $-1/r^3$ at short distance, attractive two-fermion exchange behaves as $-1/r^5$, and vector-vector two-fermion exchange is repulsive and behaves as $+1/r^5$ [2507.12522]. The neutrino-force study reaches a parallel conclusion: the long-range two-neutrino-exchange potential scales as $-1/r^5$, but the physical effect on the Sommerfeld factor depends crucially on the UV completion because the short-distance behavior becomes $1/r$ for a $t$-channel completion or $1/r^2$ for an $s$-channel completion [2203.05455].

## 4. Resonances, saturation, and unitarization

The enhancement becomes especially large when the long-range potential supports a shallow bound state or a near-threshold resonance. In the Yukawa/Hulthén analyses used for halo phenomenology, the low-velocity behavior is commonly organized into perturbative, Sommerfeld, and saturation regimes. For the $p$-wave example driven by a light CP-even scalar mediator, the intermediate regime gives $S_{l=1}\propto v^{-3}$ and hence $\langle \sigma v\rangle\propto v^{-1}$, while below
$$
v_{\rm sat}\sim |\alpha_{\rm eff}-\alpha_{\rm crit}| \equiv \Delta\alpha,
$$
the enhancement stops growing and the cross section returns to the underlying $p$-wave suppression; the cited paper gives a $p$-wave resonance near
$$
\alpha_{\rm crit}\approx 9.08\,\frac{m_\varphi}{m_\chi}.
$$
[2512.24662]

When a narrow annihilation resonance is present in addition to the long-range force, the leading-order structure simplifies rather than becoming more entangled. The effective-field-theory analysis of resonant dark-matter annihilation shows that Sommerfeld enhancement and Breit-Wigner enhancement factorize completely at leading order in the non-relativistic and narrow-width expansions,
$$
\sigma v_{\rm rel} =  \frac{C^2 Y}{4m_\chi^2} \,S_{\rm SF}(v)\, R(v),
$$
with all non-factorizable long-distance effects canceling at leading power [2209.14343].

Another recurring issue is the apparent violation of partial-wave unitarity near zero-energy resonances. In the Jost-function formulation, the conventional Sommerfeld factor is
$$
S_\ell(p)=\frac{1}{|\mathscr{J}_\ell(p)|^2},
$$
and if the long-range potential produces a zero of $\mathscr{J}_\ell$ close to threshold, the naive enhancement can become too singular. The renormalization-group treatment identifies the pathology as a secular failure of the perturbative expansion in the short-range annihilation sector and replaces the naive factor by a unitarity-consistent form with an improved Jost function, so that the pole acquires an imaginary part reflecting the decay width of the annihilating bound state [2508.09511]. A subsequent comparison of unitarization prescriptions shows that cutoff-based, RG-improved, and Green’s-function methods coincide to leading order when the unitarity-preserving corrections are large, and that the regulated cross sections are independent of the UV regulator to a good approximation [2605.05309].

The same resonance logic extends to unstable final states. For annihilation into heavier unstable particles, the width is incorporated through a complex energy argument in the final-state Schrödinger equation, and narrow bound states of the final-state pair generate resonant enhancements, while a large width washes them out. The formulation simultaneously includes off-shell final states below the nominal threshold [2603.02647].

## 5. Freeze-out, relic density, and coannihilation

In cosmology, Sommerfeld enhancement enters the Boltzmann equation through the thermally averaged annihilation cross section. For asymmetric dark matter with comoving abundances $Y_{\chi(\bar\chi)}=n_{\chi(\bar\chi)}/s$ and $x=m/T$, the standard equation is
$$
\frac{dY_{\chi(\bar\chi)}}{dx} = -\frac{\lambda\langle\sigma v\rangle}{x^2} \left(Y_\chi Y_{\bar\chi}-Y_{\chi,\rm eq}Y_{\bar\chi,\rm eq}\right),
\qquad
\lambda = 1.32\,mM_{\rm Pl}\sqrt{g_*},
$$
with the conserved asymmetry written as $Y_\chi-Y_{\bar\chi}=\eta$ [2403.01682]. In this setting the enhanced annihilation rate reduces the relic abundance by enlarging the integrated annihilation term, and the depletion is especially visible in the antiparticle abundance $Y_{\bar\chi}$ rather than in the majority species [2001.08404].

For asymmetric dark matter with a light but massive mediator, the Planck-normalized analysis yields lower bounds on the coupling and upper bounds on the mass. For the benchmark $m_\varphi=0.25$ GeV, $g_\chi=2$, and $g_*=90$, the smallest couplings in the near-symmetric limit are $\alpha_{\min}\simeq 0.0043$ for $s$-wave and $\alpha_{\min}\simeq 0.0188$ for $p$-wave, and for asymmetries $\eta = 3\times 10^{-12}$ and $\eta = 8\times 10^{-13}$ the reported maximal masses are about $143$ GeV and $535$ GeV in the $s$-wave case and about $143$ GeV and $532$ GeV in the $p$-wave case [2403.01682]. The kinetic-decoupling analysis further shows that plain $s$-wave annihilation is essentially unaffected, whereas $p$-wave and Sommerfeld-enhanced $s$- and $p$-wave annihilations remain sensitive after chemical freeze-out because the interaction rate continues to depend on velocity [2001.08404].

Sommerfeld enhancement also competes with modifications of the expansion history. In non-standard cosmology with an extra energy density $\rho_D$, faster expansion leads to earlier decoupling and therefore raises the relic abundance, while Sommerfeld enhancement lowers it. The two effects oppose each other during freeze-out, and the cited study concludes that the allowed region in coupling and perturbative annihilation cross section is enlarged, even though the maximum allowed mass is essentially unchanged relative to the standard case [2501.18209].

Coannihilation and electroweak multiplet dark matter provide additional model-specific realizations. In neutralino–sfermion coannihilation, the Sommerfeld correction is of the order of several per cent in stau regions, while in stop regions it can reach a factor of $5$ because of strong QCD attraction [1102.4295]. In the renormalizable spin-1 electroweak model studied in the supplied literature, inclusion of long-range electroweak Sommerfeld effects shifts the thermal relic window to
$$
3.6~\mathrm{TeV} \lesssim m_V \lesssim 9.2~\mathrm{TeV},
$$
within a perturbative regime, substantially above the familiar $\sim 3$ TeV scale associated with scalar or fermionic electroweak triplets [2606.16442].

## 6. Extensions, phenomenology, and nonstandard realizations

Indirect-detection phenomenology often relies on the separation between freeze-out and present-day velocities. In the halo-excess model with a light CP-even scalar mediator, the three-regime structure of perturbative behavior, Sommerfeld scaling, and saturation permits a $p$-wave annihilation model to remain thermal at freeze-out, grow to halo-scale rates, and then saturate in dwarf spheroidals. The benchmark
$$
(m_\chi,m_\varphi,f_\varphi)=(575,\ 11.6,\ 873)\ {\rm GeV}
$$
gives
$$
\langle \sigma v\rangle_{\rm Halo}\simeq 1.09\times 10^{-24}\ {\rm cm}^3/{\rm s},
$$
with dwarf constraints avoided because the enhancement shuts off in the colder systems [2512.24662].

A different mechanism appears when the long-range force itself is medium induced. In the two-component model with a dominant fermionic dark matter particle $\chi$ and an ultralight pseudoscalar background $\phi$, a finite-density loop effect produces an effective Yukawa-like potential whose coupling is proportional to the local background energy density,
$$
\alpha_{\rm bkg} \equiv \frac{4\pi \alpha_\phi \alpha_\eta \rho_\phi}{m_\chi^2 m_\phi^2}.
$$
The corresponding inner-Galaxy Sommerfeld factor can reach $\mathcal{O}(10^2)$, opening parameter space capable of fitting the Fermi-LAT Galactic Center GeV excess [2512.24188]. The broader quantum-force analysis reaches a related conclusion: coherent background corrections can generate temperature-induced resonances for massless mediators, and in a non-thermal bosonic background the effective annihilation rate can be either enhanced or suppressed because the background both creates a Coulomb-like force and shifts the effective dark-matter mass [2507.12522].

Finite-size dark matter changes the standard point-particle picture qualitatively. For puffy dark matter the potential is softened inside the overlap region, the relevant size parameter is $y=R_\chi m_\phi$, resonance peaks shift, and increasing the radius suppresses the enhancement. In the model-independent analysis summarized in the supplied material, finite size turns a sharp zero-energy resonance into a finite-width band, while in the nugget-type example the band collapses back to a sharp resonance line once the internal structure is fixed [2606.29710]. The self-interacting puffy-dark-matter study reports the same qualitative tendency and finds that for a large ratio between $R_{\chi}$ and $1/m_{\phi}$ the Sommerfeld factor approaches $1$ [2308.02170].

The concept is also not restricted to dark matter. In two-dimensional Dirac materials the absorbance above the gap is written as
$$
\eta(\omega)=\eta_0(\omega)\,S(\omega),
$$
where $S(\omega)$ is the Sommerfeld factor produced by electron–hole attraction. The cited work finds strong enhancement near the interband edge, followed by rapid reduction at higher frequencies due to single-particle energy renormalization; in graphene, by contrast, the interaction correction is small and the Sommerfeld factor stays close to $1$ [2104.12569].

Taken together, these studies present Sommerfeld enhancement not as a single formula attached to a Yukawa potential, but as a general low-velocity framework for incorporating long-distance wave-function distortion into short-distance reaction rates. The supplied literature shows that its magnitude, sign, resonance pattern, and phenomenological meaning are controlled by mediator content, background occupation, channel structure, cosmological history, and, in composite scenarios, the spatial extent of the interacting objects [2209.14343][2508.09511][2606.29710].

Source: https://www.emergentmind.com/topics/sommerfeld-enhancement