⚛️ phenomenology

Elastic Scattering of Dirac Dark Matter with Charged Leptons through a Scalar Mediator: Analytical Framework, Numerical Stability, and Direct-Detection Bounds

This paper presents a comprehensive analytical framework for the elastic scattering of Dirac dark matter off charged leptons via a scalar mediator, deriving closed-form cross-section expressions that ensure numerical stability across all mass regimes and utilizing experimental data from multiple direct-detection experiments to establish robust constraints on the model's couplings.

Original authors: Nguyen Nhu Xuan, Vu Toan Thang

Published 2026-10-28

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Nguyen Nhu Xuan, Vu Toan Thang

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The universe is filled with a mysterious substance that we cannot see, touch, or smell, yet its gravity holds galaxies together and shapes the cosmos. Astronomers call this invisible stuff dark matter, and it makes up about a quarter of everything that exists. While we know it is there because of how it pulls on stars and light, we have never directly caught a single particle of it. For decades, scientists have searched for dark matter by looking for heavy, slow-moving particles that might bump into atomic nuclei, much like a bowling ball hitting a pin. However, as these searches have become more sensitive, they have found nothing. This silence has led researchers to consider a different possibility: perhaps dark matter is lighter than we thought, and instead of hitting heavy nuclei, it interacts only with the tiny, fast-moving electrons that orbit atoms.

In this new scenario, dark matter would be a ghost that slips past the heavy parts of an atom and whispers only to the electrons. To test this idea, physicists need to know exactly how often such a whisper would happen and how much energy it would transfer. A team of researchers has now provided a precise map for this interaction. They studied a specific model where a dark matter particle, which they describe as a Dirac fermion, bounces off a charged particle like an electron. In their model, this bounce is not a direct collision but is mediated by a third particle, a real scalar boson, which acts as a messenger carrying the force between them. The researchers wanted to calculate the probability of this event happening across all possible speeds and masses, but they faced a significant hurdle: the standard mathematical tools used to make these calculations break down when the messenger particle is very heavy compared to the energy of the collision.

The team developed a new, more robust way to perform these calculations that avoids the mathematical errors that plague older methods. When the messenger particle is heavy, the standard approach involves subtracting two very large numbers to find a tiny result, a process that causes computers to lose precision and produce garbage data. The researchers rearranged the entire calculation into a form where every term adds a positive amount, ensuring the result remains accurate even when the messenger is extremely massive. They verified this new method against ultra-precise computer arithmetic, confirming that their formulas work correctly across a vast range of conditions, from dark matter that is lighter than an electron to particles as heavy as the heaviest known atoms.

With these reliable formulas in hand, the team explored how the interaction changes as the mass of the messenger particle varies. They found that when the messenger is light, the interaction behaves in a way that is very different from the heavy-messenger scenario often assumed in simpler models. Specifically, they showed that using the simplified "contact" approximation, which treats the interaction as if the messenger were infinitely heavy, can lead to massive errors. If the messenger is light, this simplified model can overestimate the chance of a collision by a factor of two or more, potentially leading scientists to rule out dark matter models that are actually still possible. The researchers established clear rules for when the simple model is safe to use and when the full, complex calculation is required.

Finally, the team applied their new, accurate calculations to data from real-world experiments designed to detect dark matter. They looked at results from several sensitive detectors, including SENSEI, DAMIC, SuperCDMS, XENON10, and XENON1T, which are currently hunting for these electron recoils. By translating their theoretical predictions into the language these experiments use, they were able to set new limits on how strongly dark matter can interact with electrons. They found that for very light messenger particles, the constraints on the interaction strength are different than previously thought, allowing for a wider range of possibilities than the old, simplified models suggested. Their work provides a solid foundation for future searches, ensuring that when scientists interpret data from these delicate experiments, they are using the correct mathematical tools to distinguish between a true signal and a mathematical illusion.

Technical Summary: Elastic Scattering of Dirac Dark Matter with Charged Leptons through a Scalar Mediator

Problem Statement
The paper addresses the phenomenology of leptophilic dark matter (DM), specifically a Dirac fermion $\chi$ interacting with Standard Model charged leptons ( $\ell = e, \mu, \tau$ ) via a real scalar mediator $S$ . While simplified models are preferred over Effective Field Theory (EFT) contact operators to maintain validity across varying momentum transfers, existing analytical treatments often rely on numerical packages (e.g., MADGRAPH5) or suffer from severe numerical instabilities. Specifically, the standard method of calculating the total elastic cross section $\sigma(s)$ by evaluating the difference of an indefinite antiderivative, $I(0) - I(t_{\min})$ , suffers from catastrophic subtractive cancellation in the heavy-mediator regime ( $|t_{\min}| \ll m_S^2$ ). This leads to a complete loss of significance in standard 64-bit floating-point arithmetic, rendering precise predictions for direct-detection bounds unreliable in certain kinematic regions.

Methodology
The authors construct a rigorous analytical framework for the process $\chi + \ell \to \chi + \ell$ in the $t$ -channel:

Theoretical Setup: The model is defined by a low-energy Yukawa interaction Lagrangian in the mass eigenbasis, derived from a gauge-invariant UV completion involving a dimension-five operator or scalar-Higgs mixing. The interaction terms are $-g_\chi S \bar{\chi}\chi - g_\ell S \bar{\ell}\ell$ .
Analytical Derivation: Closed-form tree-level expressions are derived for the invariant amplitude, the spin-averaged squared matrix element $|\mathcal{M}|^2$ , and the differential cross sections $d\sigma/dt$ and $d\sigma/d\Omega$ .
Numerical Reformulation: To resolve the numerical instability, the authors rearrange the phase-space integral. Instead of subtracting antiderivatives, they transform the integration variable to positive-definite recoil contributions. The total cross section is expressed as a sum of positive terms supplemented by convergent power-series expansions for auxiliary recoil integrals $\Phi_1(x)$ and $\Phi_2(x)$ (where $x = |t_{\min}|/m_S^2$ ).
Validation: The reformulated expressions are benchmarked against 80-digit arbitrary-precision arithmetic across a broad kinematic grid ( $m_\chi \in [10^{-4}, 10^2]$ GeV, $v_{DM} \in [10^{-4}, 10^{-1}]c$ , $m_S \in [10^{-8}, 10^2]$ GeV).
Experimental Mapping: The microscopic differential cross section is mapped to the reference electron-scattering cross section $\bar{\sigma}_e$ used in experimental literature. This mapping accounts for the mediator form factor $F_{med}(q)$ and distinguishes between fixed-reference ( $q_0 = \alpha m_e$ ) and detector-response ( $q_{eff}$ ) projections.

Key Contributions and Results

Closed-Form Analytical Expressions: The paper provides exact tree-level formulas for the scattering amplitude and cross sections without relying on EFT contact approximations, valid from the light-mediator to the heavy-mediator regimes.
Numerical Stability: The proposed positive-term reformulation eliminates subtractive cancellation. The authors demonstrate that the new method maintains a relative error of $\le 1.5 \times 10^{-15}$ across all tested regimes, whereas the naive antiderivative subtraction fails completely (underflowing to zero or producing errors $>10^2$ ) when $m_S \gtrsim 10$ GeV.
EFT Validity Criteria: The study derives the exact propagator suppression ratio $R_{EFT} = \sigma_{exact}/\sigma_{EFT}$ $R_{E F T} = σ_{e x a c t} / σ_{E F T}$ . It establishes quantitative criteria for the validity of the contact approximation:
- At the transition scale $m_S = q_{\max}$ , the EFT approximation overestimates the cross section by 100%.
- To achieve $\le 1\%$ accuracy, the mediator mass must satisfy $m_S \ge 10 q_{\max} \approx 20 \mu_{\chi\ell} v_{DM}$ .
Direct-Detection Bounds: Using experimental data from SENSEI, DAMIC, SuperCDMS, XENON10, and XENON1T, the authors map the theoretical cross sections to constraints on the effective coupling $\sqrt{g_\chi g_\ell}$ $g_{χ} g_{ℓ}$ .
- They highlight that for light mediators ( $m_S \ll \alpha m_e$ ), standard EFT extrapolations significantly underestimate the allowed coupling strength.
- The dynamical propagator leads to a saturation of the coupling bound at a constant plateau for ultra-light mediators, a feature missed by contact-interaction assumptions.
- The paper distinguishes between "fixed- $q_0$ " projections and "detector-response $q_{eff}$ -recast" methods, showing that the latter relaxes the bounds on ultra-light mediators by a factor of $\kappa_{det} \approx 1.6\text{--}2.8$ .

Significance and Claims
The paper claims to provide a robust, numerically stable analytical tool for leptophilic dark matter searches that bridges the gap between light and heavy mediator regimes without resorting to EFT approximations where they break down. By resolving the numerical loss of significance in standard arithmetic, the work ensures that quantitative validity criteria for the contact limit are established and that experimental exclusion limits are interpreted correctly, particularly for scenarios involving light scalar mediators where the momentum transfer is comparable to or exceeds the mediator mass. The authors position this framework as a necessary benchmark for numerical implementations and a precise tool for interpreting current and future direct-detection data in the sub-GeV to GeV mass range.

More like this