WIMP-Nucleon SI Cross-Section Sensitivity

• The topic defines SI sensitivity as the lowest elastic scattering cross section measurable by dark matter detectors, crucial for evaluating experimental reach.
• Methodologies include differential cross-section modeling using nuclear form factors and statistical analyses incorporating detector responses and background estimates.
• Experimental techniques, ranging from ZEPLIN-III to XENONnT, constrain particle physics models such as supersymmetry and minimal dark matter through precise sensitivity metrics.

Weakly Interacting Massive Particle (WIMP)–nucleon spin–independent (SI) cross–section sensitivity quantifies the lowest SI elastic cross section that direct dark matter detection experiments can exclude or discover as a function of WIMP mass. This sensitivity is a fundamental metric for evaluating experimental reach and directly constrains particle physics models such as supersymmetry and minimal dark matter, where WIMPs are generic dark matter candidates. The extraction and interpretation of SI sensitivity depend crucially on the interplay between detector technology, background modeling, nuclear response, and theoretical scattering formalism.

1. Theoretical Framework for SI WIMP–Nucleon Scattering

The SI WIMP–nucleus differential cross section at zero momentum transfer is parameterized by the per–nucleon cross section, $\sigma_{\text{SI}}^0 \equiv \sigma_N$. For isospin–conserving couplings ($f_p = f_n$), the WIMP–nucleus differential rate in nuclear recoil energy $E_R$ is

$$\frac{d\sigma}{dE_R} = \frac{m_A}{2\mu_A^2\,v^2}\,\sigma_{\text{SI}}^0\,A^2\,F^2(q)$$

where $m_A$ is the nuclear mass, $\mu_A = m_\chi m_A/(m_\chi+m_A)$ is the WIMP–nucleus reduced mass, $v$ is the lab–frame WIMP speed, $A$ is the atomic mass number, $q = \sqrt{2m_A E_R}$ is the momentum transfer, and $F(q)$ is the normalized scalar nuclear form factor ($f_p = f_n$0). The structure factor $f_p = f_n$1, in terms of multipole transitions, can be written as

$f_p = f_n$2

with the leading $f_p = f_n$3 multipole dominating the SI response. The per–nucleus cross section at $f_p = f_n$4 becomes

$f_p = f_n$5

where $f_p = f_n$6 is the WIMP–nucleon reduced mass. Hence, experimental limits or prospective sensitivities on $f_p = f_n$7 are converted to $f_p = f_n$8 by dividing by $f_p = f_n$9 and reduced mass scaling (Vietze et al., 2014).

2. Nuclear Structure and Form Factor Modeling

The extraction of accurate SI cross–section sensitivities depends on robust nuclear structure calculations. Large–scale shell model fits for isotopes (e.g., xenon) model the structure factor as

$E_R$0

with $E_R$1 the oscillator length and $E_R$2 isotope–specific coefficients. These fits, validated for recoil energies up to $E_R$3 keV, agree with the phenomenological Helm form factor

$E_R$4

within a few percent for $E_R$5 (Vietze et al., 2014). Nuclear–structure uncertainty in extracted $E_R$6 from Xe detectors is thus at the level of $E_R$7 in the most relevant recoil window.

3. Experimental Techniques, Analysis, and Sensitivity Metrics

The sensitivity to the SI WIMP–nucleon cross section in direct detection is established by folding the theoretical recoil spectrum, including form factor and detector response, with measured data and background models. Experiments such as ZEPLIN-III, CDMS-II/SuperCDMS, XENONnT, CDEX-50, and TEXONO use different target nuclei, exposures, thresholds, and analysis protocols.

Experimental upper limits or sensitivities are set using likelihood or counting–based statistical methods. The expected number of signal events is calculated as

$E_R$8

where $E_R$9 is the fiducial mass, $$\frac{d\sigma}{dE_R} = \frac{m_A}{2\mu_A^2\,v^2}\,\sigma_{\text{SI}}^0\,A^2\,F^2(q)$$0 exposure, $$\frac{d\sigma}{dE_R} = \frac{m_A}{2\mu_A^2\,v^2}\,\sigma_{\text{SI}}^0\,A^2\,F^2(q)$$1 detection efficiency, and $$\frac{d\sigma}{dE_R} = \frac{m_A}{2\mu_A^2\,v^2}\,\sigma_{\text{SI}}^0\,A^2\,F^2(q)$$2 the predicted rate. Limits on $$\frac{d\sigma}{dE_R} = \frac{m_A}{2\mu_A^2\,v^2}\,\sigma_{\text{SI}}^0\,A^2\,F^2(q)$$3 are then derived using profile likelihood ratios, binned Poisson methods, Feldman–Cousins, or optimum interval methods, incorporating background estimates and systematic uncertainties from detector response, energy calibration, nuclear structure, and astrophysical inputs (Akimov et al., 2011, Bruch, 2010, collaboration et al., 2020, Geng et al., 2023, Collaboration, 2013).

Representative Sensitivities

Experiment	Target & Exposure	Threshold	$$\frac{d\sigma}{dE_R} = \frac{m_A}{2\mu_A^2\,v^2}\,\sigma_{\text{SI}}^0\,A^2\,F^2(q)$$4, $$\frac{d\sigma}{dE_R} = \frac{m_A}{2\mu_A^2\,v^2}\,\sigma_{\text{SI}}^0\,A^2\,F^2(q)$$5
ZEPLIN-III	Xe, 1,344 kg⋅days	$$\frac{d\sigma}{dE_R} = \frac{m_A}{2\mu_A^2\,v^2}\,\sigma_{\text{SI}}^0\,A^2\,F^2(q)$$67 keVr	$$\frac{d\sigma}{dE_R} = \frac{m_A}{2\mu_A^2\,v^2}\,\sigma_{\text{SI}}^0\,A^2\,F^2(q)$$7 pb at 51 GeV/$$\frac{d\sigma}{dE_R} = \frac{m_A}{2\mu_A^2\,v^2}\,\sigma_{\text{SI}}^0\,A^2\,F^2(q)$$8 (Akimov et al., 2011)
SuperCDMS 15kg	Ge, 30,000 kg⋅day	few keV	$$\frac{d\sigma}{dE_R} = \frac{m_A}{2\mu_A^2\,v^2}\,\sigma_{\text{SI}}^0\,A^2\,F^2(q)$$9 cm$m_A$0 at 60 GeV/$m_A$1 (Bruch, 2010)
XENONnT	Xe, 20 t⋅y	4–50 keVnr	$m_A$2 cm$m_A$3 at 50 GeV/$m_A$4 (collaboration et al., 2020)
CDEX-50	Ge, 150 kg⋅year	160 eVee	$m_A$5 cm$m_A$6 at 5 GeV/$m_A$7 (Geng et al., 2023)

4. Operator Structure and Generalizations Beyond Standard SI Coupling

In the effective field theory (EFT) context, SI sensitivity encompasses more than the leading $m_A$8 scalar operator. A complete treatment includes isoscalar ($m_A$9), isovector ($\mu_A = m_\chi m_A/(m_\chi+m_A)$0), and two–body couplings ($\mu_A = m_\chi m_A/(m_\chi+m_A)$1, $\mu_A = m_\chi m_A/(m_\chi+m_A)$2). The generalized differential cross section, including all coherently enhanced scalar and vector responses, is (Hoferichter et al., 2016): $\mu_A = m_\chi m_A/(m_\chi+m_A)$3 Two–body effects from pion–exchange diagrams generically shift the total coherent response at the $\mu_A = m_\chi m_A/(m_\chi+m_A)$4 level unless couplings are fine–tuned.

Model–independent EFT analyses further identify subleading operators (e.g., $\mu_A = m_\chi m_A/(m_\chi+m_A)$5, $\mu_A = m_\chi m_A/(m_\chi+m_A)$6, $\mu_A = m_\chi m_A/(m_\chi+m_A)$7) that generate SI–type nuclear responses with only modest suppressions ($\mu_A = m_\chi m_A/(m_\chi+m_A)$8 or $\mu_A = m_\chi m_A/(m_\chi+m_A)$9) relative to the leading $v$0 scalar operator. In particular, derivative couplings can be constrained at the $v$1–fold greater sensitivity than expected from the naive $v$2 suppression alone due to their coupling to distinct nuclear response functions present in the SI channel (Anand et al., 2014).

5. Astrophysical Assumptions and Mass Scaling

Sensitivity projections and extraction of $v$3 require assumptions about the galactic halo. Standard inputs are $v$4 GeV/cm$v$5, Maxwellian WIMP speed distribution with $v$6 km/s, $v$7 km/s. Variation in these parameters can shift exclusion limits by tens of percent, but inter–experiment comparisons typically fix these values for consistency.

The SI cross–section sensitivity exhibits nontrivial dependence on WIMP mass ($v$8). For $v$9, the sensitivity degrades as $A$0; for $A$1, it plateaus. The optimal sensitivity for most heavy–nucleus targets occurs for $A$2 GeV, where the reduced mass is maximized (Bruch, 2010, collaboration et al., 2020).

6. Systematic and Theoretical Uncertainties

The total uncertainty in SI cross–section sensitivity has components from:

• Nuclear structure: For heavy targets (notably Xe), the difference between modern shell–model and Helm form factor is $A$3 below $A$4 keV. Two–body corrections introduce additional $A$5–$A$6 uncertainty in some EFT treatments (Vietze et al., 2014, Hoferichter et al., 2016).
• Astrophysical modeling: Variations in velocity distribution and local density can shift limits by factors of order unity; all leading results assume the Standard Halo Model.
• Experimental response: Uncertainties in energy calibration, detector threshold, exposure, and background modeling are incorporated via nuisance parameters in the limit–setting procedures (Akimov et al., 2011, Geng et al., 2023).
• Hadronic and perturbative QCD inputs: In theory–driven predictions (e.g., pure WIMP multiplets), perturbative and hadronic uncertainties can reach $A$7–$A$8 (Hill et al., 2013).

7. Implications for New Physics and Experimental Programs

The sensitivity of direct detection experiments to $A$9 constrains new physics models such as supersymmetry (neutralino LSPs), minimal dark matter (e.g., wino or higgsino multiplets), and generic WIMP EFTs. Combined with collider searches (jets + MET, monojet at LHC), and astrophysical probes, SI sensitivity enables complementarity that excludes large portions of parameter space. For example, in the MSSM, addition of direct detection (e.g., LUX, XENONnT) exclusion lines markedly expands the coverage beyond that achievable purely by LHC searches in the $q = \sqrt{2m_A E_R}$0 plane (Arbey et al., 2013). Next–generation detectors with ton–scale exposures and sub–keV thresholds (e.g., XENONnT, SuperCDMS, CDEX-50) are projected to reach the neutrino floor, below which sensitivity is limited by irreducible solar and atmospheric neutrino backgrounds (collaboration et al., 2020, Geng et al., 2023).