⚛️ phenomenology

Inelastic Dark Matter from a Higher-Dimensional Operator: The LUX-ZEPLIN High-Recoil Event and Thermal Relic

This paper proposes a fermionic inelastic dark matter model mediated by a dimension-6 operator that naturally explains a high-energy nuclear-recoil event observed by LUX-ZEPLIN and the correct thermal relic density for dark matter masses between 0.2 and 1 without requiring fine-tuning.

Original authors: Debajyoti Choudhury, Jaydeb Das, Divya Sachdeva

Published 2026-10-02

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Debajyoti Choudhury, Jaydeb Das, Divya Sachdeva

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 known as dark matter. We cannot see it, and it does not emit light, yet we know it is there because its gravity holds galaxies together and shapes the cosmos. For decades, scientists have searched for this invisible material by building massive detectors deep underground, hoping to catch a rare moment when a dark matter particle bumps into an atom in the detector. Most of these searches have come up empty, leading researchers to refine their ideas about what dark matter might be and how it behaves. The prevailing theory suggests that dark matter particles are heavy and move relatively slowly, colliding with atoms in a way that produces a small, predictable kick. However, the universe does not always follow the most common script, and sometimes a single, unusual event can force scientists to rethink the entire story.

Recently, the LUX-ZEPLIN experiment, a state-of-the-art detector filled with liquid xenon, recorded a single, striking event that did not fit the standard pattern. The detector observed a nucleus recoiling with an energy of approximately 248 keV, a value that is significantly higher than what is typically expected from the slow-moving dark matter particles that most theories predict. In the world of particle physics, such a high-energy hit usually requires a particle moving at an extraordinary speed, far faster than the average speed of dark matter in our galaxy. Furthermore, standard theories struggle to explain why this single high-energy event would appear without a flood of lower-energy events accompanying it. This anomaly has prompted a fresh look at the nature of dark matter, suggesting that the particle responsible might not be a simple, static object, but rather something that changes its state when it interacts with ordinary matter.

A team of researchers has proposed a new explanation for this puzzling signal, suggesting that the dark matter particle undergoes an "inelastic" collision. In this scenario, the dark matter particle is not a single, unchanging entity but exists in two slightly different forms, like a person wearing two different hats that are nearly identical in weight but distinct in nature. The lighter version, which makes up the dark matter in our galaxy, travels through space until it hits a xenon nucleus. Upon impact, it absorbs a tiny bit of energy from the collision to transform into its heavier, excited version. This process is similar to a ball bouncing off a wall and losing some of its speed to change its shape, but in this case, the energy cost of the transformation is so specific that it filters out the slower, more common dark matter particles. Only the fastest particles in the galactic halo have enough speed to pay the energy price for this transformation, which explains why the detector saw a high-energy event while missing the lower-energy ones that usually dominate the data.

To test this idea, the researchers built a simple mathematical model that adds just one new type of particle and one specific rule for how it interacts with the known particles of the Standard Model. They found that this minimal setup naturally creates the two different states of dark matter and allows for the energy-absorbing collision required to explain the LUX-ZEPLIN event. Crucially, the same model also accounts for the total amount of dark matter we observe in the universe today. In the early universe, these particles would have been created in equal numbers of their light and heavy forms. As the universe cooled, the heavier particles would have decayed into the lighter, stable ones, leaving behind the correct amount of dark matter to match cosmological observations. The researchers calculated that this scenario works perfectly for dark matter particles with a mass between 0.2 and 1 TeV, a range that is heavy but within the realm of possibility for modern physics.

The study demonstrates that this specific type of interaction can simultaneously explain the strange high-energy event and the overall abundance of dark matter without requiring any arbitrary adjustments to the numbers. The model predicts that the dark matter particles must be split by a very precise amount of energy, roughly between 240 and 270 keV, to align the speed requirements of the collision with the observed data. If the split were any smaller, the detector would have seen too many low-energy events; if it were larger, the particles would need to move faster than the galaxy allows, making the event impossible. The researchers emphasize that their work is a proof of concept, showing that a very simple extension of our current understanding of physics can resolve a complex mystery. While the model is purely theoretical at this stage and relies on a specific type of interaction that has not yet been proven, it offers a robust and economical path forward for interpreting the data. As more data is collected from underground detectors, scientists will be able to test whether this specific mechanism is indeed the key to unlocking the nature of the dark matter that surrounds us.

Technical Summary: Inelastic Dark Matter from a Higher-Dimensional Operator

Problem Statement
The LUX-ZEPLIN (LZ) experiment recently reported a single nuclear-recoil event with a reconstructed energy of approximately $E_R \simeq 248$ keV. Interpreting this event within the framework of standard elastic scattering of non-relativistic Weakly Interacting Massive Particles (WIMPs) is challenging. Conventional contact interactions predict a recoil spectrum that decreases with increasing energy, making an isolated high-energy event difficult to accommodate without predicting a larger excess of events at lower energies. Furthermore, producing such a large recoil typically requires particles from the high-velocity tail of the Galactic dark matter (DM) distribution, which is suppressed by the nuclear form factor at high momentum transfers. This motivates the exploration of scenarios where scattering kinematics or momentum dependence differ from standard elastic WIMP scattering, specifically through inelastic dark matter (iDM) models.

Methodology and Theoretical Framework
The authors propose a minimal extension of the Standard Model (SM) involving a gauge-singlet Dirac fermion $\chi$ interacting with SM fields via a single dimension-6 effective operator. The Lagrangian is defined as:
$\mathcal{L}_\chi = \bar{\chi}(i\not{\partial}-m_\chi)\chi + \frac{C_V}{\Lambda^2}(\bar{\chi}\gamma^\mu\chi)(H^\dagger i\overleftrightarrow{D}_\mu H)$
where $\Lambda$ is the EFT cutoff scale and $C_V$ is a real Wilson coefficient.

To facilitate inelastic scattering while suppressing elastic scattering, the global $U(1)$ symmetry of the Dirac fermion is softly broken by a mass term $m'(\chi^c\chi + \text{h.c.})$ . This breaks the symmetry to a $Z_2$ , splitting the Dirac fermion into two Majorana states, $\chi_1$ (lighter, DM candidate) and $\chi_2$ (heavier), with a mass splitting $\Delta m = 2m'$ . The interaction current becomes off-diagonal ( $\bar{\chi}_1\gamma^\mu\chi_2$ ), enabling the endothermic transition $\chi_1 + N \to \chi_2 + N$ while forbidding elastic $\chi_1 N \to \chi_1 N$ scattering at tree level.

The analysis proceeds in two main steps:

Relic Density Calculation: The thermal relic abundance is calculated within the freeze-out paradigm. Since $\Delta m \ll T_f$ (freeze-out temperature), the abundances of $\chi_1$ and $\chi_2$ are nearly equal post-freeze-out. The depletion is governed by co-annihilation processes ( $\chi_1\chi_2 \to Zh, W^+W^-, f\bar{f}$ ) rather than direct annihilation. The required coupling $C_V/\Lambda^2$ is determined to match the observed relic density ( $\Omega_\chi h^2 \approx 0.12$ ).
LZ Event Interpretation: The differential nuclear-recoil rate is computed for the inelastic process $\chi_1 N \to \chi_2 N$ . The analysis incorporates the Standard Halo Model (SHM), the Helm form factor for xenon isotopes, and the detector's energy resolution and acceptance. The minimum velocity required to produce a recoil $E_R$ is non-monotonic in $E_R$ for inelastic scattering, peaking at a specific energy $E^*_R$ . The model parameters are tuned such that $E^*_R \approx 248$ keV.

Key Results

Simultaneous Explanation: The model successfully explains both the observed relic density and the LZ high-energy event without fine-tuning. This is achieved for a DM mass range of $m_\chi \sim (0.2, 1)$ TeV.
Mass Splitting Constraints: To reproduce the LZ signal profile (a single event at 248 keV with no excess at lower energies), the mass splitting is constrained to a narrow window: $240 \text{ keV} \lesssim \Delta m \lesssim 270 \text{ keV}$ . This splitting ensures the kinematic threshold selects the high-velocity tail of the DM distribution, suppressing lower-energy recoils.
Coupling Strength: The required Wilson coefficient $C_V/\Lambda^2$ is moderate. At lower DM masses, larger couplings are needed due to phase-space suppression from SM particle masses ( $W, Z, h$ ). At higher masses, the cross-section scales as $(C_V m_\chi/\Lambda^2)^2$ , allowing for smaller couplings.
Consistency Checks:
- Decay Constraints: The heavier state $\chi_2$ decays to $\chi_1$ via off-shell $Z$ bosons or loops ( $\chi_2 \to \chi_1 + \nu\bar{\nu}$ or $n\gamma$ ). The lifetime is calculated to be much shorter than the age of the Universe, ensuring $\chi_2$ decays well before the present epoch. The small mass splitting prevents hadron production or photodisintegration of light elements, rendering cosmological constraints from Big Bang Nucleosynthesis (BBN) and the CMB negligible.
- Direct Detection: The model naturally suppresses elastic scattering, consistent with null results from other direct detection experiments.

Significance and Claims
The paper claims to demonstrate that a very simple, economic extension of the Standard Model—consisting of a single fermionic field and one dimension-6 operator—suffices to explain the puzzling LZ data while simultaneously reproducing the observed relic density. The authors emphasize that the model is phenomenological, remaining agnostic about the specific ultraviolet (UV) completion that would generate the specific operator. The work provides a robust template for future model building, suggesting that the interplay between inelastic kinematics and co-annihilation can resolve the tension between high-energy recoil events and standard WIMP expectations. The authors note that while other phenomenologically acceptable operators exist, switching them on would primarily alter quantitative results rather than the qualitative success of the framework.

Technical Summary: Inelastic Dark Matter from a Higher-Dimensional Operator

More like this