Scalar Singlet Dark Matter Candidate

• Scalar singlet dark matter is a minimal Standard Model extension where an extra scalar field, protected by a discrete symmetry, acts as a stable DM candidate.
• It exhibits both WIMP and FIMP regimes, enabling precise calculations of relic density and unique detection prospects through Higgs portal interactions.
• Experimental and theoretical constraints, including collider bounds and direct detection limits, rigorously define the viable mass and coupling parameter space.

A scalar singlet dark matter candidate refers to a Standard Model (SM) extension involving a new real or complex scalar field, singlet under the SM gauge group, stabilized by an imposed discrete symmetry (typically $\mathbb{Z}_2$ or variants), with phenomenology determined primarily by the portal coupling(s) to the Higgs sector. This framework produces minimal and predictive dark matter (DM) candidates, offering both weakly interacting massive particle (WIMP) and feebly interacting massive particle (FIMP) regimes and supports efficient calculation of relic density, direct detection rates, and collider signatures. Scalar singlet DM models are among the most constrained—and most extensively studied—single-particle DM frameworks in the literature.

1. Formal Model Definition and Lagrangian Structure

The core scalar singlet DM model extends the SM by a real scalar $S$ (or complex $S$), with a stabilizing symmetry ($\mathbb{Z}_2$, $U(1)$, or similar), yielding the Lagrangian: $$\mathcal{L} = \mathcal{L}_{\rm SM} + \frac12(\partial_\mu S)\partial^\mu S - \frac12 m_S^2 S^2 - \frac{\lambda_S}{4} S^4 - \frac{\lambda_{HS}}{2} S^2 H^\dagger H$$ Here, $H$ is the SM Higgs doublet; $m_S$ the bare singlet mass; $\lambda_{HS}$ the Higgs portal coupling; $\lambda_S$ the singlet self-coupling (Yaguna, 2011, Collaboration et al., 2017). After electroweak symmetry breaking (EWSB), $S$0 ($S$1GeV).

Variants include complex singlet scenarios with $S$2 or $S$3 stabilization (Gonderinger et al., 2012), two-scalar extensions (Bazzocchi et al., 2012), and constructions with extended scalar sectors (e.g., triplet extensions (Campbell et al., 2016), composite Higgs frameworks [(Cai et al., 2020) ]).

2. Relic Density Dynamics: Freeze-out and Freeze-in

Two principal regimes control relic density:

WIMP Regime (Thermal freeze-out):

The number density $S$4 follows the Boltzmann equation,

$S$5

Annihilation proceeds via $S$6-channel Higgs exchange, $S$7 SM SM, with thermally averaged cross section $S$8 computed via standard integrals over phase space (Collaboration et al., 2017, Biswas et al., 2011). The relic abundance is determined by freeze-out at $S$9: $S$0 where $S$1 are the $S$2- and $S$3-wave coefficients extracted from $S$4 expansions.

FIMP Regime (Freeze-in):

For $S$5, $S$6 never thermalizes. Its abundance accrues via out-of-equilibrium 2$S$72 production from the SM plasma: $S$8 with $S$9. The final relic abundance scales as $\mathbb{Z}_2$0—in stark contrast to WIMP models ($\mathbb{Z}_2$1) (Yaguna, 2011). In this regime, direct and indirect detection signals are negligible.

3. Parameter Space, Phenomenology, and Detection Constraints

The model's phenomenology is fixed by $\mathbb{Z}_2$2 and $\mathbb{Z}_2$3 (or generalizations for multi-scalar or multi-portal constructions). Global fits (e.g., GAMBIT (Collaboration et al., 2017)) scan across DM mass $\mathbb{Z}_2$4 (from $\mathbb{Z}_2$51 GeV to multi-TeV) and portal couplings up to $\mathbb{Z}_2$6.

Viable Regions:

• Higgs resonance: $\mathbb{Z}_2$7 (with $\mathbb{Z}_2$8 the physical Higgs mass), tiny $\mathbb{Z}_2$9 ($U(1)$0–$U(1)$1), region is highly fine-tuned but allows $U(1)$2 to saturate all DM.
• High-mass terrace: $U(1)$3 1 TeV with $U(1)$4–3, testable by future ton-scale experiments (Collaboration et al., 2017).
• FIMP window: $U(1)$5, $U(1)$6–$U(1)$7, completely dark in direct and indirect detection (Yaguna, 2011).

Direct Detection:

Spin-independent DM-nucleon cross section (via $U(1)$8-channel Higgs exchange): $U(1)$9 Experimental bounds from LUX, XENON1T, PandaX, etc., exclude much of the $$\mathcal{L} = \mathcal{L}_{\rm SM} + \frac12(\partial_\mu S)\partial^\mu S - \frac12 m_S^2 S^2 - \frac{\lambda_S}{4} S^4 - \frac{\lambda_{HS}}{2} S^2 H^\dagger H$$0–$$\mathcal{L} = \mathcal{L}_{\rm SM} + \frac12(\partial_\mu S)\partial^\mu S - \frac12 m_S^2 S^2 - \frac{\lambda_S}{4} S^4 - \frac{\lambda_{HS}}{2} S^2 H^\dagger H$$1 plane for $$\mathcal{L} = \mathcal{L}_{\rm SM} + \frac12(\partial_\mu S)\partial^\mu S - \frac12 m_S^2 S^2 - \frac{\lambda_S}{4} S^4 - \frac{\lambda_{HS}}{2} S^2 H^\dagger H$$2 GeV at moderate $$\mathcal{L} = \mathcal{L}_{\rm SM} + \frac12(\partial_\mu S)\partial^\mu S - \frac12 m_S^2 S^2 - \frac{\lambda_S}{4} S^4 - \frac{\lambda_{HS}}{2} S^2 H^\dagger H$$3 (Collaboration et al., 2017). For $$\mathcal{L} = \mathcal{L}_{\rm SM} + \frac12(\partial_\mu S)\partial^\mu S - \frac12 m_S^2 S^2 - \frac{\lambda_S}{4} S^4 - \frac{\lambda_{HS}}{2} S^2 H^\dagger H$$4 (FIMP), $$\mathcal{L} = \mathcal{L}_{\rm SM} + \frac12(\partial_\mu S)\partial^\mu S - \frac12 m_S^2 S^2 - \frac{\lambda_S}{4} S^4 - \frac{\lambda_{HS}}{2} S^2 H^\dagger H$$5–$$\mathcal{L} = \mathcal{L}_{\rm SM} + \frac12(\partial_\mu S)\partial^\mu S - \frac12 m_S^2 S^2 - \frac{\lambda_S}{4} S^4 - \frac{\lambda_{HS}}{2} S^2 H^\dagger H$$6 cm$$\mathcal{L} = \mathcal{L}_{\rm SM} + \frac12(\partial_\mu S)\partial^\mu S - \frac12 m_S^2 S^2 - \frac{\lambda_S}{4} S^4 - \frac{\lambda_{HS}}{2} S^2 H^\dagger H$$7, entirely unobservable (Yaguna, 2011).

Collider Constraints:

Invisible Higgs decays provide critical limits: $$\mathcal{L} = \mathcal{L}_{\rm SM} + \frac12(\partial_\mu S)\partial^\mu S - \frac12 m_S^2 S^2 - \frac{\lambda_S}{4} S^4 - \frac{\lambda_{HS}}{2} S^2 H^\dagger H$$8 for $$\mathcal{L} = \mathcal{L}_{\rm SM} + \frac12(\partial_\mu S)\partial^\mu S - \frac12 m_S^2 S^2 - \frac{\lambda_S}{4} S^4 - \frac{\lambda_{HS}}{2} S^2 H^\dagger H$$9 (Collaboration et al., 2017). Direct production of scalar singlet DM is not accessible at current energies except via missing-energy searches and precision measurements of Higgs width.

4. Extensions and Theoretical Variants

Multi-singlet scenarios introduce additional stabilizing symmetries (e.g., $H$0 or $H$1), hence supporting multicomponent DM (Basak et al., 2021, Belanger et al., 2021). These produce new phenomena including semi-annihilations and conversion processes (e.g., $H$2) relaxing direct detection constraints in multi-component frameworks.

Composite models such as the $H$3 pNGB scenario yield singlet DM candidates whose couplings are loop-induced and whose masses are set by vacuum misalignment, offering viable DM in the $H$4 GeV–$H$5 TeV range (Cai et al., 2020). Novel annihilation channels (to heavy exotic scalars) enable relic density saturation for modified couplings and UV completions.

Scenarios addressing additional issues (such as the little hierarchy problem (Bazzocchi et al., 2012), neutrino masses (Bhattacharya et al., 2016), or dark energy (Landim, 2017)) integrate the scalar singlet with extended scalar sectors, yielding altered quartic mixing, mass sum rules, and additional annihilation channels to quadruplet or triplet states, with enhanced parametric freedom.

5. Vacuum Stability, RG Running, and Perturbativity

Vacuum stability imposes nontrivial requirements on the quartic couplings: $H$6 These are enforced up to a high cutoff $H$7 (typically TeV–$H$8 GeV), with perturbativity constraints $H$9 or more conservative bounds (Gonderinger et al., 2012, Landim, 2017). RG running of $m_S$0 and $m_S$1 can induce instability (usually for large negative $m_S$2), reversed by higher-dimension operators ($m_S$3, $m_S$4 terms) (Landim, 2017).

In composite scenarios, the stability is further protected by accidental discrete symmetries inherited from the UV theory (Cai et al., 2020).

6. Indirect Detection and Astrophysical Implications

Indirect detection signals—primarily gamma-ray observations—are sensitive to scalar singlet annihilation to $m_S$5 and rare two-photon final states. FIMP scenarios and regions with suppressed $m_S$6 are not observable, while resonance or heavy territory may be accessible to Fermi-LAT, H.E.S.S., or CTA depending on parameter choices (Basak et al., 2021, Gaitan et al., 2014, Landim, 2017).

Self-interaction cross sections (mainly set by $m_S$7) can be tuned to match astrophysical small-scale structure constraints (e.g., core–cusp, Bullet Cluster), with light ($m_S$8) singlet models providing nonthermal DM and Bose–Einstein condensate scenarios for galactic halos (Matos et al., 2014).

7. Experimental Outlook and Future Probes

A large section of the WIMP parameter space will be tested by XENONnT, LZ, DARWIN and future colliders (HL-LHC, ILC, FCC-ee) (Collaboration et al., 2017, Campbell et al., 2016). FIMP regions are likely to remain inaccessible. Multi-component models and singlet scenarios with suppressed portal couplings will require novel detection strategies, possibly targeting exotic signatures, semi-annihilations, and precision Higgs or electroweak observables.

Composite scenarios and extensions with nontrivial scalar sectors predict direct-detection cross sections near or just below the neutrino floor, as well as rich collider phenomenology including mono-$m_S$9 signatures and invisible decays of non-SM Higgs partners (Cai et al., 2020, Dutta et al., 2022, Bazzocchi et al., 2012).

In sum, scalar singlet dark matter remains an exceptionally active research topic, fully calculable, and testably predictive, with only narrow allowed windows of parameter space persisting under current and projected experimental constraints. The interplay of relic density, direct detection, indirect signals, and theoretical consistency dictates the feasible regimes for this minimal dark sector.