Pseudo-Dirac Dark Matter Candidate

• Pseudo-Dirac dark matter candidates are models featuring a Dirac fermion that acquires small Majorana masses, resulting in two nearly degenerate mass eigenstates.
• They interpolate between Dirac and Majorana phenomenology, influencing s-wave coannihilation and thermal freeze-out mechanisms that determine relic abundance.
• Experimental signatures include inelastic direct detection, distinctive collider displaced vertices, and unique astrophysical signals that test the model's parameter space.

A pseudo-Dirac dark matter candidate consists of a new Dirac fermion whose two Weyl components acquire small but nonzero Majorana mass terms, resulting in a pair of nearly degenerate Majorana mass eigenstates separated by a small splitting. Such models arise in a broad class of dark sector constructions, including minimal extensions with new singlet and doublet fermions, dark-photon portals, Dirac-gaugino supersymmetry, and radiative neutrino mass models. Pseudo-Dirac dark matter generically interpolates between the phenomenology of Dirac (s-wave annihilation, unsuppressed vector-like couplings) and Majorana (suppressed vector interactions, p-wave annihilation) particles, depending on the size of the splitting and mediation mechanism. This structure underpins a variety of robust signatures and model-building strategies for satisfying relic abundance, direct- and indirect-detection, collider, and astrophysical constraints.

1. Theoretical Structure and Mass Spectrum

The defining feature of pseudo-Dirac dark matter is a new Dirac spinor—typically denoted χ or Ψ—composed of two Weyl fermions (ξ, η), endowed with a Dirac mass term $m_D$ and small Majorana masses $m_{χ_L}$, $m_{χ_R}$ (or $\mu_{L,R}$, depending on basis): $$\mathcal{L} \supset -m_D \bar{\chi} \chi - \frac{1}{2} m_{χ_L} \bar{\chi}_L^c \chi_L - \frac{1}{2} m_{χ_R} \bar{\chi}_R^c \chi_R +\text{h.c.}$$ In the limit $m_{χ_{L,R}} \ll m_D$, diagonalization yields two Majorana mass eigenstates $χ_1$ and $χ_2$ with

$$m_{1,2} = m_D \mp m_M,\qquad m_M \equiv \frac{m_{χ_L}+m_{χ_R}}{2}$$

so that the splitting $\delta m = m_2 - m_1 \ll m_D$ defines the "pseudo-Dirac" regime (Simone et al., 2010, Konar et al., 2020).

In scenarios with SM gauge charges or mixing, pseudo-Dirac dark matter may arise from singlet-doublet mixing (Konar et al., 2020), Dirac gauginos with R-symmetry breaking (Goodsell et al., 2015), or Higgsino–singlino mixing in high-scale SUSY (Mummidi et al., 2018), always leading to a similar quasi-degenerate pair structure.

2. Relic Abundance and Thermal Mechanisms

The relic abundance of pseudo-Dirac dark matter is fixed by thermal freeze-out in most constructions, typically via s-wave coannihilation of the two quasi-degenerate Majorana states. The annihilation cross section retains essentially Dirac-like efficiency, since processes such as $m_{χ_L}$0 SM are unsuppressed, while self-annihilations $m_{χ_L}$1, $m_{χ_L}$2 are p-wave or suppressed by $m_{χ_L}$3 (Simone et al., 2010, Jordan et al., 2018). The effective cross section is (Simone et al., 2010): $m_{χ_L}$4 with

$m_{χ_L}$5

and $m_{χ_L}$6 at freeze-out.

In models with additional mixing, such as singlet-doublet or Higgsino–singlino frameworks, the Boltzmann equations must account for the thermal populations and interactions of all near-degenerate states (charged and neutral), but coannihilation generally remains dominant as long as the mass splittings remain below $m_{χ_L}$7 GeV (Konar et al., 2020, Goodsell et al., 2015, Mummidi et al., 2018).

Freeze-in production (rather than freeze-out) is possible if the coupling is extremely suppressed; in dipole-interaction models with keV splittings and MeV–GeV masses, the observed relic abundance can be achieved via freeze-in enabled by transition dipole moments (Chatterjee et al., 2022).

3. Direct and Indirect Detection Signatures

Spin-Independent and Inelastic Scattering

Direct detection via spin-independent (SI) elastic scattering is generically suppressed for pseudo-Dirac fermions with pure Majorana interactions, as diagonal Z-boson vector couplings cancel at leading order (Simone et al., 2010, Konar et al., 2020, Mummidi et al., 2018). Off-diagonal couplings induce inelastic up-scattering $m_{χ_L}$8, but this is kinematically inaccessible if $m_{χ_L}$9, where $m_{χ_R}$0–$m_{χ_R}$1 keV is the typical momentum exchange in terrestrial detectors. Constraints require $m_{χ_R}$2 keV for tree-level Z-exchange to be forbidden kinematically at XENON1T (Mummidi et al., 2018).

For inter-state splittings $m_{χ_R}$3 in the keV range, down-scattering of the subdominant excited component $m_{χ_R}$4 can proceed, depositing a monoenergetic electron-recoil $m_{χ_R}$5 in semiconductor or xenon detectors. This is a target for ongoing and planned low-threshold experiments and can fit potential excesses such as the Xenon1T 2–3 keV range (Chao et al., 2020, González et al., 2021, Chatterjee et al., 2022, Brahma et al., 2023). In models with only transition dipole couplings (MDM or EDM), the electron channel dominates and the predicted rates are sharply peaked at low $m_{χ_R}$6 (Chatterjee et al., 2022).

Indirect Detection

Pseudo-Dirac dark matter preserves efficient s-wave annihilation via $m_{χ_R}$7 SM particles, so indirect detection bounds from gamma-ray, $m_{χ_R}$8, or neutrino fluxes apply, subject to the late-time abundance of the excited state. For small splittings (sub-keV), excited-state depletion is often efficient, and present-day annihilation is minimal; but in the resonant regime or for weak kinetic decoupling, a non-negligible $m_{χ_R}$9 abundance persists, leading to indirect detection signals both from annihilations and characteristic de-excitation/decay signatures (Brahma et al., 2023, González et al., 2021). CMB bounds are satisfied as long as down-scattering processes and late-time annihilations are sufficiently suppressed by the small excited fraction or Boltzmann factors (Brahma et al., 2023, Garcia et al., 2024).

4. Realizations and Model Embeddings

Minimal Effective Theory

A generic effective operator realization features a SM-singlet Dirac fermion with dimension-6 interactions with SM fermion currents, where the mass-splitting arises from Majorana corrections (Simone et al., 2010). The hallmark signatures are efficient freeze-out through off-diagonal (coannihilation) operators, suppressed direct detection, and the collider displaced-vertex signature from $\mu_{L,R}$0 decay to $\mu_{L,R}$1.

Singlet-Doublet/Scotogenic Models

In singlet-doublet extensions of the SM, the $\mu_{L,R}$2 sector mixes via Yukawa couplings, and Majorana mass terms induce splitting, radiative neutrino mass (via a loop mechanism), and link DM to neutrino phenomenology. The diagonalization yields lightest Majorana DM $\mu_{L,R}$3 with sub-EV–consistent neutrino mass and all direct, indirect, and collider bounds are respected for moderate mixing angles $\mu_{L,R}$4 and splittings $\mu_{L,R}$5 GeV (Konar et al., 2020).

Supersymmetry and Dirac Gaugino Sectors

In Dirac gaugino models and high-scale SUSY, neutralino or Higgsino–singlino mixing leads to pseudo-Dirac fermions. Small mass splittings $\mu_{L,R}$6 from R-symmetry breaking (e.g., via the $\mu_{L,R}$7 term) ensure direct detection safety while preserving correct relic density via (co)annihilation and compatibility with collider and vacuum-stability constraints (Goodsell et al., 2015, Mummidi et al., 2018).

Dark Photon and Dipole Portals

Dark sector U(1)′ models assign off-diagonal gauge couplings between the Majorana mass eigenstates, with kinetic mixing ($\mu_{L,R}$8) to the SM photon/hypercharge. In this setup, annihilation and scattering are mediated by a dark photon or through dipole operators (González et al., 2021, Brahma et al., 2023, Mohlabeng et al., 2024). Dipole models realize keV-scale splittings and enable direct detection via both solar-upscattered and halo components (Chatterjee et al., 2022).

5. Collider and Intensity Frontier Probes

Pseudo-Dirac dark matter predicts distinctive collider signatures stemming from the decay $\mu_{L,R}$9, where the decay length is of order centimeters for typical weak-scale masses and GeV splittings: $$\mathcal{L} \supset -m_D \bar{\chi} \chi - \frac{1}{2} m_{χ_L} \bar{\chi}_L^c \chi_L - \frac{1}{2} m_{χ_R} \bar{\chi}_R^c \chi_R +\text{h.c.}$$0 Observing a displaced vertex associated with missing energy and reconstructing the mass edge allows extraction of both the DM mass and splitting, providing a direct collider-cosmology test (Simone et al., 2010).

At the intensity frontier, models with dark-photon mediation predict signals in fixed-target and beam-dump searches via production and detection of $$\mathcal{L} \supset -m_D \bar{\chi} \chi - \frac{1}{2} m_{χ_L} \bar{\chi}_L^c \chi_L - \frac{1}{2} m_{χ_R} \bar{\chi}_R^c \chi_R +\text{h.c.}$$1, $$\mathcal{L} \supset -m_D \bar{\chi} \chi - \frac{1}{2} m_{χ_L} \bar{\chi}_L^c \chi_L - \frac{1}{2} m_{χ_R} \bar{\chi}_R^c \chi_R +\text{h.c.}$$2 in meson decays followed by scattering or visible $$\mathcal{L} \supset -m_D \bar{\chi} \chi - \frac{1}{2} m_{χ_L} \bar{\chi}_L^c \chi_L - \frac{1}{2} m_{χ_R} \bar{\chi}_R^c \chi_R +\text{h.c.}$$3 decays within the detector. Experiments like NA64, JSNS$$\mathcal{L} \supset -m_D \bar{\chi} \chi - \frac{1}{2} m_{χ_L} \bar{\chi}_L^c \chi_L - \frac{1}{2} m_{χ_R} \bar{\chi}_R^c \chi_R +\text{h.c.}$$4, LSND, and planned runs by Belle II have or will probe large regions of the viable parameter space, particularly for sub-GeV mass windows and splittings $$\mathcal{L} \supset -m_D \bar{\chi} \chi - \frac{1}{2} m_{χ_L} \bar{\chi}_L^c \chi_L - \frac{1}{2} m_{χ_R} \bar{\chi}_R^c \chi_R +\text{h.c.}$$5 GeV (Jordan et al., 2018, Garcia et al., 2024).

6. Cosmological and Astrophysical Implications

Pseudo-Dirac sterile neutrino constructions provide a framework for warm dark matter, evade X-ray constraints via permutation symmetry, and explain the relic abundance via sequential production and decay of the heavy state. The late-time decay injects dark radiation, shifting $$\mathcal{L} \supset -m_D \bar{\chi} \chi - \frac{1}{2} m_{χ_L} \bar{\chi}_L^c \chi_L - \frac{1}{2} m_{χ_R} \bar{\chi}_R^c \chi_R +\text{h.c.}$$6 and potentially ameliorating the Hubble tension at the level of $$\mathcal{L} \supset -m_D \bar{\chi} \chi - \frac{1}{2} m_{χ_L} \bar{\chi}_L^c \chi_L - \frac{1}{2} m_{χ_R} \bar{\chi}_R^c \chi_R +\text{h.c.}$$7 (Chao et al., 2021). In sub-GeV resonant regimes, the decoupling of kinetic and chemical freeze-out leaves an $$\mathcal{L} \supset -m_D \bar{\chi} \chi - \frac{1}{2} m_{χ_L} \bar{\chi}_L^c \chi_L - \frac{1}{2} m_{χ_R} \bar{\chi}_R^c \chi_R +\text{h.c.}$$8 fractional population of excited states, altering the usual relationship between thermal history, CMB limits, and today’s signal rates (Brahma et al., 2023, González et al., 2021).

Self-interactions and exothermic processes in halos are generically suppressed but can open novel signatures for specific parameter choices (Brahma et al., 2023). Solar upscattering is important in models where the ground state composes halo DM, but the excited state is efficiently replenished via interactions in the Sun, enhancing detection prospects (Chatterjee et al., 2022).

7. Parameter Space, Radiative Corrections, and Experimental Prospects

The allowed parameter space of pseudo-Dirac DM candidates spans a wide mass range: MeV–TeV regimes are accessible via different mechanisms and mediators. Direct detection is primarily sensitive through inelastic or exothermic channels, setting bounds on transition dipole couplings or the effective cross section on electrons. Next-generation low-threshold detectors (XENONnT, SuperCDMS, SENSEI, Oscura, DARWIN) and accelerator-based missing-energy searches (NA64, LDMX, Belle II) are poised to fully test the parameter space favored by thermal- or freeze-in production for sub-GeV masses and splittings from eV to keV (Chatterjee et al., 2022, Brahma et al., 2023, Mohlabeng et al., 2024, Garcia et al., 2024).

Radiative corrections can shift the thermal targets for kinetic mixing $$\mathcal{L} \supset -m_D \bar{\chi} \chi - \frac{1}{2} m_{χ_L} \bar{\chi}_L^c \chi_L - \frac{1}{2} m_{χ_R} \bar{\chi}_R^c \chi_R +\text{h.c.}$$9 or cross sections by 10–20%, impacting sensitivity projections and necessitating precise calculations for mapping theoretical predictions to experiment (Mohlabeng et al., 2024). The window for pseudo-Dirac dark matter—combining relic density, suppressed direct detection, testable collider/beam signals, and novel astrophysical signatures—is among the most broadly scrutinized of contemporary WIMP and sub-GeV models.