Sub-MeV Electrophilic Dark Matter

• Sub-MeV electrophilic dark matter is defined as a class of light dark-sector particles below 1 MeV that preferentially interact with electrons, influencing both cosmological and experimental outcomes.
• Direct detection strategies employ low-threshold experiments, including semiconductor superlattices, Dirac materials, and doped semiconductors, to probe electron recoils with high sensitivity.
• A combination of cosmological, astrophysical, and laboratory constraints narrows the viable parameter space and motivates alternative models with light mediators and absorption mechanisms.

Sub-MeV electrophilic dark matter refers to hypothesized dark-sector particles with mass below 1 MeV that couple preferentially—or exclusively—to electrons rather than nuclei. This interaction channel is of particular theoretical and experimental interest given both the suppressed nuclear-recoil rates for such light dark matter (DM) and the unique cosmological implications of their possible electromagnetic couplings. The sub-MeV regime is characterized by strong cosmological, astrophysical, and laboratory constraints, and is the subject of active research in both direct detection methodologies and model building.

1. Effective Field Theory and Basic Model Structures

Sub-MeV electrophilic DM is typically described within an effective field theory (EFT) framework. The canonical scenario involves a single DM particle χ, stabilized by a Z₂ symmetry and singlet under the Standard Model (SM) gauge group. DM–electron interactions are assumed to be mediated by a heavy new particle (mass $m_\phi \gg \mathrm{MeV}$), which, when integrated out, yields a leading four-fermion operator: $$\mathcal{L}_\mathrm{int} = \frac{g_\chi g_e}{M^2} \, \bar\chi\chi \, \bar e e + \ldots$$ where $g_\chi$ and $g_e$ are the DM– and electron–mediator couplings, $M \simeq m_\phi$, and the ellipsis denotes other possible operator structures (scalar, pseudoscalar, vector, etc.) (Lehmann et al., 2020). The EFT description is valid for momentum transfers and plasma temperatures $T$, $q \ll m_\phi$.

In alternative formulations, the mediator can be a light vector (e.g., a kinetically mixed dark photon $A'$), leading to a form-factor–dependent coupling. In this case, the Lagrangian reads: $$\mathcal{L}_\mathrm{portal} \supset \epsilon F_{\mu\nu}F'^{\mu\nu} + g_\chi A'_\mu \bar\chi\gamma^\mu\chi + e A_\mu \bar e\gamma^\mu e$$ where $\epsilon$ is the kinetic-mixing parameter (Hochberg et al., 2017, Aboubrahim et al., 2021).

Electrophilic refers specifically to the hierarchy $$\mathcal{L}_\mathrm{int} = \frac{g_\chi g_e}{M^2} \, \bar\chi\chi \, \bar e e + \ldots$$0, such that electron scattering dominates over nucleon scattering at achievable thresholds (Zhang et al., 2024).

2. Cosmological and Astrophysical Constraints

Stringent constraints arise from the early universe, predominantly from the following epochs:

A. Big-Bang Nucleosynthesis (BBN). If χ (or its mediator) equilibrates with electrons at $$\mathcal{L}_\mathrm{int} = \frac{g_\chi g_e}{M^2} \, \bar\chi\chi \, \bar e e + \ldots$$1 MeV, its additional energy density alters the Hubble rate, affecting helium and deuterium yields. BBN constraints require the DM–electron coupling satisfy (Lehmann et al., 2020): $$\mathcal{L}_\mathrm{int} = \frac{g_\chi g_e}{M^2} \, \bar\chi\chi \, \bar e e + \ldots$$2 For $$\mathcal{L}_\mathrm{int} = \frac{g_\chi g_e}{M^2} \, \bar\chi\chi \, \bar e e + \ldots$$3 MeV and $$\mathcal{L}_\mathrm{int} = \frac{g_\chi g_e}{M^2} \, \bar\chi\chi \, \bar e e + \ldots$$4 keV: $$\mathcal{L}_\mathrm{int} = \frac{g_\chi g_e}{M^2} \, \bar\chi\chi \, \bar e e + \ldots$$5.

B. Relic Abundance (Freeze-In/Out-of-Equilibrium Production). Overproduction occurs if annihilation $$\mathcal{L}_\mathrm{int} = \frac{g_\chi g_e}{M^2} \, \bar\chi\chi \, \bar e e + \ldots$$6 is too small, restricting the coupling from below. For correct relic density ($$\mathcal{L}_\mathrm{int} = \frac{g_\chi g_e}{M^2} \, \bar\chi\chi \, \bar e e + \ldots$$7), one needs (Lehmann et al., 2020): $$\mathcal{L}_\mathrm{int} = \frac{g_\chi g_e}{M^2} \, \bar\chi\chi \, \bar e e + \ldots$$8 C. Extra Radiation ($$\mathcal{L}_\mathrm{int} = \frac{g_\chi g_e}{M^2} \, \bar\chi\chi \, \bar e e + \ldots$$9). Late-time entropy transfer to/from the electron-photon bath shifts the photon-to-neutrino temperature ratio. Planck and future CMB experiments tightly constrain $g_\chi$0 (2σ), limiting additional light degrees of freedom at $g_\chi$1 MeV (Lehmann et al., 2020, Green et al., 2017). For typical mediator and DM parameter choices, the $g_\chi$2 constraint is as restrictive as BBN.

D. Stellar Cooling. Electron-coupled mediators can lead to excessive energy loss from stars (white dwarfs, red giants). White-dwarf cooling constrains $g_\chi$3 below $g_\chi$4 for $g_\chi$5 keV (Green et al., 2017).

E. Direct Astrophysical Observables. In dense environments (e.g., white dwarfs near high-density DM regions), sub-MeV DM can alter cooling via capture, scattering, and annihilation, providing complementary constraints in $g_\chi$6–$g_\chi$7 space (Zhang et al., 2024).

Combined Parameter Space: Cosmology, BBN, and relic-abundance constraints carve out a narrow viable window: only $g_\chi$8–$g_\chi$9 and $g_e$0 few hundred keV evade all bounds in standard EFT. Lower $g_e$1 values are practically excluded by order-of-magnitude (Lehmann et al., 2020).

3. Direct Detection Approaches and Projected Sensitivities

Sub-MeV electrophilic DM-induced electron recoils are sought in ultralow-threshold experiments, leveraging a variety of condensed-matter targets and detection strategies:

A. Semiconductor Superlattices and Quantum-Cascade Lasers (QCLs). Superlattice superstructures (SSS) engineered for gaps $g_e$2–300 meV enable detection of sub-MeV DM recoils via prompt mid-infrared photon emission, with QCL-based readout reaching practical thresholds $g_e$3 meV and projected $g_e$4 C.L. sensitivities $g_e$5–$g_e$6 cm$g_e$7 for $g_e$8–$g_e$9 MeV—orders of magnitude beyond current semiconductor limits (Bora, 2022).

B. Dirac and Graphene-Based Materials. Three-dimensional Dirac semimetals ($M \simeq m_\phi$0 meV) or voltage-tunable bilayer graphene enable detection thresholds as low as $M \simeq m_\phi$110 meV (Hochberg et al., 2017, Das et al., 2023, Geilhufe et al., 2019). Their unscreened in-medium response allows probing cross sections $M \simeq m_\phi$2–$M \simeq m_\phi$3 cm$M \simeq m_\phi$4 for $M \simeq m_\phi$5 in the 4 keV–1 MeV range, exploiting daily modulation due to anisotropic response and Earth's rotation as a potential discriminant against background.

C. Doped Semiconductors and Skipper-CCD Technology. P-type or n-type semiconductors with dopant-induced shallow energy levels ($M \simeq m_\phi$610–100 meV) enable access to lower DM masses than pure (eV-gap) materials. Projected sensitivities with exposures of $M \simeq m_\phi$7100 g·day and dark counts $M \simeq m_\phi$81/(g·day) can reach the freeze-in target cross sections ($M \simeq m_\phi$9–$T$0 cm$T$1), contingent on improved noise and backgrounds (Du et al., 2022).

D. Plasmon-Enhanced and Quantum Materials Approaches. Plasmon excitations in e.g., silicon CC(D)s, particularly with cosmic-ray (CR)–boosted DM, can yield strong limits on $T$2 down to $T$3 cm$T$4 at $T$5 keV in the light-mediator scenario (Liang et al., 2024).

E. Boosted, Reflected, or Absorption Signals. Additional handles come from CR-boosted DM (Cao et al., 2020, Shang et al., 2024) or solar-reflected DM, which can produce signals well above the kinematic thresholds of halo DM and yield direct-detection constraints for $T$6 as low as several keV (An et al., 2017).

Target Type	Minimum $T$7, ($T$8)	Best Current/Projected Sensitivity ($T$9)	Unique Features
Dirac materials	$q \ll m_\phi$0 keV (few meV)	$q \ll m_\phi$1–$q \ll m_\phi$2 cm$q \ll m_\phi$3 (3–100 events/kg·yr)	No in-medium suppression, daily modulation possible
Superlattice superstructure	$q \ll m_\phi$4 MeV (50 meV)	$q \ll m_\phi$5 cm$q \ll m_\phi$6 (1 kg·yr)	Tunable gap, photon-based readout
Doped semiconductors	30 keV (10–100 meV)	$q \ll m_\phi$7 cm$q \ll m_\phi$8 (100 g·day, DC=0)	Leverages shallow levels, mature technology

No platform currently achieves sensitivity to the couplings allowed by cosmology in the minimal heavy-mediator models; see Section 5.

4. Model Extensions, Absorption Channels, and Non-Minimal Cosmology

While elastic scattering signatures are tightly constrained by cosmology and freeze-in/out-of-equilibrium history, alternative mechanisms can open viable parameter space:

A. Fermionic Absorption by Electrons. Models where light fermionic DM is absorbed (e.g., $q \ll m_\phi$9 + A $A'$0 $A'$1 + A$A'$2 + $A'$3) rather than scattered can evade some cosmological bounds (Dror et al., 2020). For vector-mediated absorption, XENON1T already explores $A'$4 keV $A'$5 MeV with projected reach to lower masses in future liquid-xenon TPCs.

B. Multi-sector or Dark-Sink Cosmologies. Scenarios with multiple hidden sectors or late-time entropy injection can alleviate BBN and $A'$6 constraints, allowing heavier couplings or different thermal histories (Aboubrahim et al., 2021, Bhattiprolu et al., 2024). For instance, the introduction of a "Dark Sink"—a bath of very light fermions interacting with DM—modifies freeze-in and allows present-day cross-sections up to several orders of magnitude above the canonical freeze-in line: $A'$7 for $A'$8 keV (Bhattiprolu et al., 2024).

C. Bosonic Absorption. Absorption of bosonic DM (e.g., dark photons) in targets with low excitation thresholds can provide observable signals independent of velocity distributions and with weaker cosmological model dependence (Hochberg et al., 2017, Geilhufe et al., 2019, Das et al., 2023), achieving sensitivity to kinetic-mixing parameters as low as $A'$9 for sub-eV–MeV dark photon masses.

5. Combined Parameter Space, Experimental Outlook, and Limitations

The convergence of cosmological, astrophysical, and laboratory data imposes a dramatic narrowing of the allowed parameter space for minimal sub-MeV electrophilic DM with heavy mediators. In the ($$\mathcal{L}_\mathrm{portal} \supset \epsilon F_{\mu\nu}F'^{\mu\nu} + g_\chi A'_\mu \bar\chi\gamma^\mu\chi + e A_\mu \bar e\gamma^\mu e$$0, $$\mathcal{L}_\mathrm{portal} \supset \epsilon F_{\mu\nu}F'^{\mu\nu} + g_\chi A'_\mu \bar\chi\gamma^\mu\chi + e A_\mu \bar e\gamma^\mu e$$1) or ($$\mathcal{L}_\mathrm{portal} \supset \epsilon F_{\mu\nu}F'^{\mu\nu} + g_\chi A'_\mu \bar\chi\gamma^\mu\chi + e A_\mu \bar e\gamma^\mu e$$2, $$\mathcal{L}_\mathrm{portal} \supset \epsilon F_{\mu\nu}F'^{\mu\nu} + g_\chi A'_\mu \bar\chi\gamma^\mu\chi + e A_\mu \bar e\gamma^\mu e$$3) plane, combined BBN, $$\mathcal{L}_\mathrm{portal} \supset \epsilon F_{\mu\nu}F'^{\mu\nu} + g_\chi A'_\mu \bar\chi\gamma^\mu\chi + e A_\mu \bar e\gamma^\mu e$$4, and relic-density constraints restrict:

• $$\mathcal{L}_\mathrm{portal} \supset \epsilon F_{\mu\nu}F'^{\mu\nu} + g_\chi A'_\mu \bar\chi\gamma^\mu\chi + e A_\mu \bar e\gamma^\mu e$$5 (effective coupling): $$\mathcal{L}_\mathrm{portal} \supset \epsilon F_{\mu\nu}F'^{\mu\nu} + g_\chi A'_\mu \bar\chi\gamma^\mu\chi + e A_\mu \bar e\gamma^\mu e$$6–$$\mathcal{L}_\mathrm{portal} \supset \epsilon F_{\mu\nu}F'^{\mu\nu} + g_\chi A'_\mu \bar\chi\gamma^\mu\chi + e A_\mu \bar e\gamma^\mu e$$7,
• $$\mathcal{L}_\mathrm{portal} \supset \epsilon F_{\mu\nu}F'^{\mu\nu} + g_\chi A'_\mu \bar\chi\gamma^\mu\chi + e A_\mu \bar e\gamma^\mu e$$8 keV,
• $$\mathcal{L}_\mathrm{portal} \supset \epsilon F_{\mu\nu}F'^{\mu\nu} + g_\chi A'_\mu \bar\chi\gamma^\mu\chi + e A_\mu \bar e\gamma^\mu e$$9 cm$\epsilon$0 (Lehmann et al., 2020, Green et al., 2017).

Proposed next-generation electron-recoil experiments with thresholds of few meV and exposures $\epsilon$1 1 kg·yr target $\epsilon$2 values at least 1–3 orders of magnitude above the cosmological upper bounds. Thus, in a minimal EFT, no unconstrained parameter space remains for observable elastic scattering unless the cosmological background is non-standard (e.g., late-time phase transition, entropy injection, or strong number-changing processes in the dark sector).

However, models with:

• Light mediators (kinetically mixed dark photons) and non-minimal cosmological histories,
• Absorption-based signals (including bosonic or fermionic DM),
• CR- or solar-boosted detection channels,

remain viable and in several cases testable by current or forthcoming low-threshold experiments (Dror et al., 2020, Bhattiprolu et al., 2024, Liang et al., 2024, An et al., 2017).

6. Key Future Directions and Open Issues

Efforts continue on several fronts:

1. Experimental Development. Lowering detection thresholds toward single–electron or single–phonon sensitivity, expanding target materials (e.g., Dirac materials, bilayer graphene, doped semiconductors) (Das et al., 2023, Du et al., 2022, Geilhufe et al., 2019).
2. Theoretical Refinement. Incorporating full atomic and condensed-matter structure factors at low energy transfers, in-medium corrections for light mediators, and precise calculations of backgrounds and daily modulation signatures (Hochberg et al., 2017).
3. Astrophysical Complementarity. Using observations of stellar cooling and white-dwarf pulsation rates in regions of high DM density to extend sensitivity to parameter regions inaccessible to terrestrial detectors (Zhang et al., 2024).
4. Non-Minimal Cosmology. Exploring more complex cosmological histories (multisector, "Dark Sink" depletion, late entropy injection) that enlarge the allowed direct-detection parameter space (Bhattiprolu et al., 2024, Aboubrahim et al., 2021).

In summary, while minimal heavy-mediator sub-MeV electrophilic dark matter is stringently constrained by cosmological considerations, several well-motivated model extensions—especially those involving light mediators, absorption channels, or alternative production histories—remain open and will continue to motivate both theoretical and experimental advances in the field (Lehmann et al., 2020, Green et al., 2017, Bhattiprolu et al., 2024).