Global Study of the Simplest Scalar Phantom Dark Matter Model

We present a global study of the simplest scalar phantom dark matter model. The best fit parameters of the model are determined by simultaneously imposing (i) relic density constraint from WMAP, (ii) 225 live days data from direct experiment XENON100, (iii) upper limit of gamma-ray flux from Fermi-LAT indirect detection based on dwarf spheroidal satellite galaxies, and (iv) the Higgs boson candidate with a mass about 125 GeV and its invisible branching ratio no larger than 40% if the decay of the Higgs boson into a pair of dark matter is kinematically allowed. The allowed parameter space is then used to predict annihilation cross sections for gamma-ray lines, event rates for three processes mono-b jet, single charged lepton and two charged leptons plus missing energies at the Large Hadron Collider, as well as to evaluate the muon anomalous magnetic dipole moment for the model.

💡 Research Summary

The paper presents a comprehensive global analysis of the simplest scalar‑phantom dark‑matter model, in which a real singlet scalar field S is stabilized by a Z₂ symmetry and interacts with the Standard Model solely through the Higgs‑portal term λ S² H†H. The authors treat the dark‑matter mass m_S and the portal coupling λ as the only free parameters and confront the model with four independent experimental/observational constraints: (i) the relic‑density measurement from WMAP (Ω_c h² = 0.1123 ± 0.0035), (ii) the 225‑day XENON100 direct‑detection limit on the spin‑independent scattering cross‑section σ_SI, (iii) the Fermi‑LAT upper limits on gamma‑ray fluxes from dwarf spheroidal galaxies, which constrain the annihilation cross‑sections into γγ and γZ final states, and (iv) the LHC observation of a Higgs‑like boson with a mass around 125 GeV together with the requirement that its invisible branching ratio does not exceed 40 % when the decay h → SS is kinematically allowed.

Using a numerical solution of the Boltzmann equation, the authors compute the thermally averaged annihilation cross‑section ⟨σv⟩ required to reproduce the observed relic density. They evaluate σ_SI analytically (σ_SI ∝ λ² f_N² μ_N² / m_h⁴ m_S²) and implement the latest nuclear form‑factor values. For indirect detection they employ PPPC4DMID and DarkSUSY to obtain the gamma‑ray spectra and compare them with the Fermi‑LAT dwarf limits. The invisible Higgs width is calculated from Γ(h → SS) = (λ² v² / 8π m_h) √(1 − 4 m_S²/m_h²) and fed into HiggsSignals/HiggsBounds to enforce the BR_inv ≤ 0.4 constraint.

Scanning over m_S from 1 GeV to 1 TeV and λ from 10⁻⁴ to 1, the authors identify two distinct viable regions. The first is a low‑mass “resonance” band (m_S ≲ 55 GeV) where 2 m_S ≈ m_h enhances annihilation through an s‑channel Higgs pole; here the relic density can be achieved with very small λ (10⁻³–10⁻²). However, XENON100 imposes a stringent σ_SI bound, forcing λ to be ≲ 10⁻², and the invisible‑Higgs limit further restricts λ to ≤ 0.01. The second region lies at higher masses (m_S ≳ 200 GeV). In this regime larger couplings (λ ≈ 0.1–0.3) are allowed because direct‑detection limits are weaker, but the annihilation cross‑section must remain below the Fermi‑LAT dwarf limits, which requires a careful balance between λ and m_S. The authors find that masses around 300–600 GeV with λ ≈ 0.1–0.2 provide the widest allowed parameter space.

Having delineated the allowed parameter space, the paper proceeds to predict collider signatures at the 13 TeV LHC. Three processes are simulated with MadGraph5_aMC@NLO, Pythia 8, and Delphes: (1) mono‑b‑jet + missing transverse energy (MET) from pp → h* → SS + b‑jet, (2) single charged lepton + MET from pp → W* → ℓ ν + SS, and (3) dilepton + MET from pp → Z* → ℓ⁺ℓ⁻ + SS. Existing ATLAS/CMS search strategies are recast to estimate signal efficiencies. In the low‑mass resonance region, the mono‑b channel offers the best sensitivity, though current data only set weak limits (signal efficiencies of order 1–2 %). In the high‑mass region, the dilepton + MET channel is most promising because of its low Standard‑Model background; with an integrated luminosity of 300 fb⁻¹ the authors anticipate a 3σ excess for m_S ≈ 400 GeV and λ ≈ 0.15. The single‑lepton channel provides complementary coverage but is less powerful due to larger W+jets backgrounds.

The authors also evaluate the model’s contribution to the muon anomalous magnetic moment (g‑2)μ. The dominant effect arises from two‑loop Barr‑Zee diagrams involving the Higgs portal. For the allowed λ–m_S values, the calculated shift Δa_μ lies in the range (0.5–1.5) × 10⁻¹¹, far below the current discrepancy between experiment and the Standard Model (≈ 2.7 × 10⁻⁹). Consequently, the scalar‑phantom model does not address the (g‑2)μ anomaly.

In summary, the paper demonstrates that the minimal scalar‑phantom dark‑matter model survives all present constraints only within a narrowly defined parameter space: a low‑mass Higgs‑resonance strip with very small portal coupling, and a higher‑mass window around a few hundred GeV with moderate coupling. Future direct‑detection experiments such as XENONnT, LZ, and DARWIN, which aim for σ_SI sensitivities down to 10⁻⁴⁸ cm², will likely probe and possibly exclude the low‑mass region entirely. At the LHC, high‑luminosity runs (up to 3 ab⁻¹) will improve the reach of dilepton + MET searches, offering a realistic chance to test the high‑mass window. The model’s inability to explain the muon (g‑2) anomaly suggests that any simultaneous solution to dark matter and (g‑2) would require additional fields or interactions beyond the simple Higgs portal. The work thus provides a clear roadmap for experimental tests of one of the most economical dark‑matter scenarios.