Spin-Dependent WIMP Capture

• Spin-dependent WIMP capture is the process where WIMPs scatter off nucleons via axial-vector interactions in astrophysical bodies, becoming gravitationally bound.
• This mechanism uses detailed nuclear structure inputs, combining chiral effective field theory and shell-model calculations to derive cross sections and capture rates.
• Experimental bounds from direct detection and neutrino telescopes tightly connect theory with astrophysical observations, refining dark matter interaction constraints.

Spin-dependent WIMP (Weakly Interacting Massive Particle) capture refers to the process by which WIMPs, through axial-vector (spin-dependent) interactions, scatter off nucleons inside astrophysical bodies—primarily the Sun and Earth—and lose sufficient kinetic energy to become gravitationally bound. This mechanism directly links the limits from direct detection experiments with indirect searches, such as neutrino telescopes, and is sensitive to the details of both nuclear structure and the particle physics underlying WIMP interactions.

1. Fundamental Formalism for Spin-Dependent WIMP Capture

The capture rate $C$ of WIMPs via spin-dependent (SD) elastic scattering in a body of radius $R$ is given, in the optically thin (single-scattering) limit, by integrating the local dark matter flux times differential cross section over the phase space, spatial distribution, and composition of the capturing body. The master equation, valid for arbitrary nuclear content and halo velocity distribution, is: $$C = \frac{\rho_\chi}{m_\chi} \int_0^R 4\pi r^2 dr \sum_A n_A(r) \int_0^\infty du\, \frac{f(u)}{u} \, w(u,r) \int_{E_{\min}}^{E_{\max}} dE_R \frac{d\sigma_A^{SD}}{dE_R}(w)$$ where:

• $\rho_\chi$ is the local DM density, $m_\chi$ the WIMP mass,
• $f(u)$ is the WIMP speed distribution at infinity,
• $w(u,r) = \sqrt{u^2 + v_{\rm esc}^2(r)}$ is the WIMP speed at radius $r$,
• $n_A(r)$ is the number density of nuclear species $A$ at $r$,
• $E_{\min} = \frac{1}{2} m_\chi u^2$ (the minimum recoil energy for capture),
• $E_{\max} = 2 \mu_{A}^2 w^2/m_A$ (the kinematic maximum for a WIMP–nucleus interaction; $\mu_A$ is the reduced mass),
• $d\sigma_A^{SD}/dE_R$ is the differential spin-dependent cross section.

The cross section encodes the particle physics and nuclear structure inputs essential for meaningful predictions and constraints.

2. Spin-Dependent WIMP–Nucleus Cross Section and Nuclear Structure Functions

The spin-dependent WIMP–nucleus cross section depends critically on nuclear spin structure and WIMP couplings. The zero–momentum-transfer differential cross section is typically expressed as: $$\frac{d\sigma_A^{SD}}{dq^2} = \frac{8 G_F^2}{(2J+1) v^2} S_A(q)$$ where:

• $G_F$ is the Fermi constant,
• $J$ is the spin of the target nucleus,
• $v$ is the WIMP–nucleus relative speed,
• $S_A(q)$ is the nuclear spin structure function.

The structure function decomposes as: $$S_A(q) = a_0^2 S_{00}(q) + a_0 a_1 S_{01}(q) + a_1^2 S_{11}(q)$$ with $a_0$ (isoscalar) and $a_1$ (isovector) WIMP–nucleon axial-vector couplings. For practical application, calculations such as those in "Large-scale nuclear structure calculations for spin-dependent WIMP scattering with chiral effective field theory currents" (Klos et al., 2013) provide $S_{ij}(q)$ for relevant isotopes, incorporating both state-of-the-art shell-model wavefunctions and leading long-range two-body chiral EFT currents. Two-body currents induce a universal suppression of the isovector coupling by $\sim$15–30%, reducing the total SD cross section by $(1-\delta)^2 \simeq 0.5$–$0.6$ (Divari et al., 2013, Menendez et al., 2012).

At zero momentum transfer ($q\to0$), $S_A(0)$ reduces to a form involving the proton and neutron spin expectation values $\langle S_p \rangle$ and $\langle S_n \rangle$, leading to cross sections: $$\sigma_A^{SD}(0) = \frac{4 \mu_A^2}{\pi} \frac{J+1}{J} [a_p \langle S_p \rangle + a_n \langle S_n \rangle]^2$$ with $a_p$ (proton) and $a_n$ (neutron) couplings.

3. Effective Theory and Operator Structure in Spin-Dependent Scattering

Spin-dependent capture is controlled by a set of operators within the non-relativistic effective field theory (NREFT) framework, which systematically enumerates all allowed Galilean-invariant interactions for a WIMP with spin $j_\chi \leq 1/2$ (Kang et al., 2022, Widmark, 2017). Key SD operators include:

• $\hat O_4 = \vec S_\chi \cdot \vec S_N$ (axial–axial, "standard" SD),
• $$\hat O_6 = (\vec S_\chi \cdot \vec q/m_N)(\vec S_N \cdot \vec q/m_N)$$ (momentum-suppressed SD),
• $\hat O_7 = \vec S_N \cdot \vec v^\perp$ (spin-velocity),
• $\hat O_9$, $\hat O_{10}$, etc. (spin–momentum and mixed operators).

Each operator gives rise to a unique combination of nuclear response functions ($W_{\Sigma'}, W_{\Sigma''},$ etc.), with the canonical SD interaction mediated by $\Sigma'/\Sigma''$ channels. In the case of a massless mediator ($1/q^2$ propagator), cross sections are IR-enhanced and require regularization (see Section 5).

The NREFT capture rate can be written as a bilinear in Wilson coefficients $c_j^\tau$, integrated over response and halo functions: $$C_{SD} = \sum_{j,k} \sum_{\tau,\tau'} c_j^\tau c_k^{\tau'} \int_0^\infty du\, \eta(u)\, W_{jk}^{\tau\tau'}(u)$$ as implemented in flexible codes such as WimPyC (Kang et al., 24 Oct 2025).

4. Experimental Direct Detection Constraints and Their Folding into Capture

Direct detection experiments (e.g., LUX, PICO, XENONnT) set upper limits on SD cross sections, usually quoted for pure proton ($a_p=1,a_n=0$) and neutron ($a_p=0,a_n=1$) couplings. For instance, LUX Run 3 reports 90% CL bounds at $m_\chi = 33$ GeV/c$^2$: $\sigma_n = 9.4\times10^{-41}$ cm$^2$ and $\sigma_p = 2.9\times10^{-39}$ cm$^2$ (Collaboration et al., 2016). These constraints, together with nuclear spin values, map onto capture cross sections for target nuclei in astrophysical bodies.

To connect direct search results to capture, one translates per-nucleon cross section bounds into per-nucleus quantities using: $$\sigma_{A}^{SD}(q) = \frac{4 \mu_A^2}{\pi} \frac{J+1}{J} [a_p \langle S_p \rangle + a_n \langle S_n \rangle]^2 S_A(q)$$ and incorporates these into the full capture integrals.

The capture rate for hydrogen (dominant in the Sun) thus scales generically as: $$C_\odot \propto \frac{\rho_0}{m_\chi} \, \sigma_p \langle \Phi(u) \rangle$$ where $\langle \Phi(u) \rangle$ encodes the solar potential and velocity integrals (Collaboration et al., 2016, Baum et al., 2016).

5. Astrophysical Capture Phenomenology and Halo-Independent Approaches

Capture by the Sun is extremely sensitive to spin-dependent WIMP–proton couplings due to the Sun's high $^1$H abundance, and by neutron couplings in the case of neutron-rich elements in the Earth or in direct detection targets (Baum et al., 2016). Capture rates for the Sun with $\sigma_p^{SD} = 10^{-40}$ cm$^2$ and $m_\chi = 100$ GeV yield $C_\odot^{SD} \sim 10^{24}$ s$^{-1}$. For the Earth, including all relevant isotopes enhances the rate by a factor of three over previous tabulations (Baum et al., 2016).

Uncertainties in the local halo WIMP velocity distribution translate into uncertainties in capture rates and derived cross section limits. The Standard Halo Model (SHM)—a truncated Maxwell-Boltzmann velocity distribution—can yield bounds that are up to orders of magnitude stronger than those obtained using strictly halo-independent (single-stream) methods. In the latter, the most conservative bound on coupling constants $c_i^p$ at fixed $m_\chi$ is derived by maximizing over all possible $f(u)$ subject to normalization (Kang et al., 2022, Choi et al., 2024); for SD WIMP-proton couplings, the relaxation factor compared to SHM is generally $\lesssim 10$ outside the 10–200 GeV region, where it can reach $\sim10^2$.

Massless mediator scenarios ($1/q^2$ enhancement in the amplitude) lead to UV-finite but IR-divergent capture rates unless the maximum WIMP aphelion is limited (the "Jupiter cut"). The single-stream halo-independent method circumvents this sensitivity as it is independent of low-$u$ tails (Choi et al., 2024).

6. Spin-Dependent Capture, Thermalization, and Indirect Detection Probes

After capture, WIMPs thermalize with the stellar nucleus distribution, settling to a density profile described by the solar core temperature and gravitational potential. Thermalization times range from $10^4$ to $10^8$ years depending on operator structure and $m_\chi$; these are nearly always much less than the Solar age of $4.5 \times 10^9$ yr for cross sections at current experimental limits (Widmark, 2017). Thus, instantaneous thermalization is an excellent approximation except for severely suppressed (e.g., momentum-suppressed) cases.

Captured WIMPs may annihilate to Standard Model products (notably neutrinos for indirect searches). The relationship between capture rate $C$ and annihilation rate $\Gamma_A$ depends on the equilibrium condition: $$\Gamma_A = \frac{C}{2} \tanh^2 \left( t/\tau \right)$$ with $\tau = (CC_A)^{-1/2}$, $C_A$ parameterizing the annihilation cross section and effective volume. For $t \gg \tau$ (equilibrium), $\Gamma_A \approx C/2$, and neutrino flux limits translate directly into constraints on $C$ and thereby $\sigma_p$, $\sigma_n$, or general NREFT couplings (Baum et al., 2016, Liang et al., 2013).

7. Nuclear Structure, Operator Uncertainties, and Renormalization

The nuclear physics underlying spin-dependent capture is characterized by:

• Large-scale shell-model wavefunctions,
• Chiral EFT-derived one-body and two-body axial currents (Menendez et al., 2012, Klos et al., 2013),
• Operator truncations and uncertainty bands arising from nuclear density $\rho$, low-energy constants $c_3$, $c_4$, and shell-model choices.

The dominant uncertainty on computed $S_A(q)$ and thus on capture rates is $\sim$15–35%, rooted both in nuclear-structure and current-operator truncation. Two-body currents suppress the effective isovector SD coupling by a nearly universal factor $(1-\delta) \sim 0.7$–$0.8$, yielding a $\sim$0.5–0.6 overall suppression in cross section and capture rate, independent of the nuclear mass number $A$ (Divari et al., 2013, Klos et al., 2013).

Summary Table: Spin-Dependent Capture—Core Quantities

Quantity	Physical Content	Reference Equation(s)
$d\sigma_A^{SD}/dq^2$	SD differential cross section	$8G_F^2/(2J+1)v^2 S_A(q)$ (Collaboration et al., 2016, Klos et al., 2013)
$S_A(q)$	Nuclear spin structure function	$a_0^2 S_{00} + a_0 a_1 S_{01} + a_1^2 S_{11}$ (Menendez et al., 2012, Klos et al., 2013)
$\sigma_A^{SD}(0)$	Total SD cross section at $q=0$	$$4\mu_A^2/\pi(J+1)/J [a_p\langle S_p \rangle + a_n\langle S_n \rangle]^2$$
$C_{SD}$	SD capture rate	See master equations above; e.g., Eq. (7) in (Collaboration et al., 2016)
Renorm. factor $R$	Two-body-current suppression	$R = (1-\delta) \sim 0.75$ (Divari et al., 2013, Menendez et al., 2012)

Spin-dependent WIMP capture theory tightly integrates nuclear physics, particle physics operator structure, and astrophysical modeling. The most stringent constraints on $\sigma_p$ arise from solar capture (and thus solar neutrino fluxes), especially for contact and massless mediator interactions, while neutron cross section constraints are currently dominated by direct detection with neutron-odd targets. Two-body current effects and halo uncertainties are essential for robust interpretation. These capture rates underpin all indirect probes of WIMP annihilation in astrophysical bodies and constitute the bridge between experimental non-observations and theoretical models of DM scattering.