Revisiting Singlet Fermion Dark Matter with a Scalar Portal: Connecting Higgs Phenomenology and Strong Electroweak Phase Transition

We investigate a minimal extension of the Standard Model with a real singlet scalar and a singlet Dirac fermion acting as dark matter. Unlike a conventional singlet scalar setup, we assume that the singlet scalar does not acquire a vacuum expectation value at zero temperature. This decouples the scalar mixing angle from the Higgs-portal quartic coupling responsible for the strong first-order electroweak phase transition, allowing it to coexist with current collider and direct-detection constraints. The Higgs-singlet mixing is generated independently through a trilinear portal interaction. We check theoretical consistency conditions, various LHC limits on heavy scalar resonances, dark matter relic abundance, and direct detection bounds to delineate the viable parameter space. We perform a detailed analysis of the electroweak phase transition and show that a strong first-order transition is realized for a selected set of benchmark points. We further compute the resulting stochastic gravitational wave spectra and find that several scenarios yield signals potentially observable at future space-based interferometers. Our results establish a unified and testable framework that connects collider phenomenology, first-order electroweak phase transition, and the resulting production of gravitational waves, along with the dark matter phenomenology, all within a simple renormalizable extension of the Standard Model.

💡 Research Summary

In this work the authors present a minimal renormalizable extension of the Standard Model (SM) that simultaneously addresses three major shortcomings of the SM: the nature of dark matter (DM), the origin of the baryon asymmetry of the Universe, and the possibility of observable stochastic gravitational waves (GWs) from the early Universe. The model adds a real gauge‑singlet scalar field s and a Dirac fermion χ, both odd under an imposed Z₂ symmetry that guarantees the stability of χ, which therefore plays the role of a weakly‑interacting massive particle (WIMP) dark matter candidate. Crucially, the scalar is assumed not to acquire a vacuum expectation value (VEV) at zero temperature. Consequently the Higgs–singlet mixing angle θ is generated solely by a dimension‑three trilinear portal term μₕₛ |H|² s, while the quartic portal coupling λₕₛ remains free to be large. This decouples the mixing angle from the coupling that controls the finite‑temperature effective potential, allowing a strong first‑order electroweak phase transition (SFOEWPT) without violating current Higgs‑signal‑strength measurements or direct‑detection limits.

The scalar potential is written as
V(H,s)=μ²ₕ|H|²+λₕ|H|⁴+½μ²ₛ s²+¼λₛ s⁴+¼λₕₛ|H|² s²+½μₕₛ|H|² s+⅓μ³ₛ s³+a₁ s,
with the linear term a₁ fixed by the condition that the minimum at zero temperature lies at (v,0). After electroweak symmetry breaking the CP‑even mass matrix mixes h and s, giving two mass eigenstates h₁ (identified with the 125 GeV Higgs) and h₂ (a heavier scalar). The independent input parameters are {v, m_{h₁}=125 GeV, m_{h₂}, sin θ, g_χ, m_χ, λₕₛ, μ₃, λₛ}.

The authors impose a series of theoretical and experimental constraints:

• Vacuum stability – positivity of the quartic couplings (λₕ, λₛ ≥ 0, λₕₛ ≥ −2√(λₕλₛ)).
• Perturbative unitarity – bounds on |λₕ|, |λₛ|, |λₕₛ| from 2→2 scalar scattering.
• Electroweak precision observables – the S and T parameters constrain |sin θ| ≲ 0.2 for m_{h₂} ≳ 800 GeV.
• LHC heavy‑scalar searches – resonant production σ(pp→h₂)≈sin²θ σ_SM(m_{h₂}) is compared with ATLAS/CMS limits in the h₂→hh, ZZ→4ℓ, and WW channels. The most stringent bound comes from ZZ→4ℓ, excluding sin θ ≳ 0.13–0.20 for 250–700 GeV, while di‑Higgs searches exclude sin θ ≳ 0.36 in the 460–600 GeV window. Overall, sin θ ≲ 0.1 for m_{h₂} ≲ 600 GeV.
• Higgs signal‑strength – the universal scaling μ_sig = cos²θ yields μ_ATLAS = 1.04 ± 0.06, implying sin θ ≲ 0.2, compatible with the collider limits.

Dark‑matter phenomenology is governed by the Yukawa interaction L_{χs}=g_χ \barχχ s. The relic density is computed with MicrOMEGAs, taking into account s‑channel χχ→SM via h₁/h₂ exchange and t‑channel χχ→ss annihilation when kinematically allowed. The authors scan over (m_χ, g_χ) and find viable regions that reproduce the observed Ω h²≈0.12 while respecting the LUX‑2022 (and newer XENONnT) spin‑independent limits. Because the direct‑detection cross section depends on the mixing angle (through h₁/h₂ exchange) and not on λₕₛ, the model can accommodate relatively large λₕₛ values needed for a strong phase transition without being excluded. Future collider prospects (mono‑X, invisible Higgs decays) are discussed, showing that the parameter space can be probed at the HL‑LHC and future lepton colliders.

The finite‑temperature effective potential V_eff(h,s,T) is built at one‑loop order, including the Daisy resummation of bosonic Matsubara zero modes. The authors evaluate the critical temperature T_c and the order‑parameter v_c (the Higgs VEV at T_c). A strong first‑order transition requires v_c/T_c > 1. By varying λₕₛ and μₕₛ while keeping sin θ small, they identify benchmark points (BP1–BP4) with m_{h₂}=300–800 GeV, sin θ≈0.05–0.1, λₕₛ≈0.6–1.2 that satisfy v_c/T_c > 1. The phase‑transition dynamics (bubble nucleation temperature T_n, inverse duration β/H, latent heat α) are computed using CosmoTransitions.

Gravitational‑wave spectra are then derived from the three standard sources: bubble collisions, sound‑waves in the plasma, and magnetohydrodynamic turbulence. For the selected benchmarks the peak frequencies lie in the millihertz to decihertz range, with peak amplitudes Ω_{GW}h²≈10⁻¹²–10⁻⁹. The authors overlay the sensitivity curves of planned space‑based interferometers (LISA, DECIGO, BBO, Taiji, TianQin) and show that BP2 (m_{h₂}=500 GeV, λₕₛ≈0.9) and BP4 (m_{h₂}=800 GeV, λₕₛ≈1.2) produce signals that could be observed by LISA (optimistic configuration) or DECIGO.

In summary, the paper demonstrates that a simple SM extension with a non‑VEV singlet scalar and a singlet Dirac fermion can simultaneously:

1. Provide a viable WIMP dark‑matter candidate compatible with relic‑density and direct‑detection data.
2. Generate a strong first‑order electroweak phase transition, a prerequisite for electroweak baryogenesis, by exploiting a large Higgs‑portal quartic coupling while keeping the Higgs‑singlet mixing angle small.
3. Predict stochastic gravitational‑wave signals within the reach of upcoming space‑based detectors, thereby linking collider physics, cosmology, and GW astronomy in a testable framework.

The analysis is thorough, covering theoretical consistency, current LHC bounds, dark‑matter phenomenology, finite‑temperature phase‑transition dynamics, and GW prospects. The decoupling of the mixing angle from λₕₛ is the key innovation that resolves the tension present in many singlet‑scalar models between a strong phase transition and direct‑detection limits. The work therefore offers a compelling, minimal, and experimentally accessible scenario for new physics beyond the SM.