Interpreting the 650 GeV and 95 GeV Higgs anomalies in the next-to-two-Higgs-doublet model

Recent experimental hints from the Large Hadron Collider (LHC) in diphoton and partially in the τ + τ − final states suggest the possible existence of an additional Higgs boson with a mass around 95 GeV. Interestingly, these observations are consistent with earlier results from the Large Electron-Positron collider, which pointed to an excess in b b ¯ final states within a similar mass range. Additionally, CMS has observed an excess in the γ γ b b ¯ final state, indicating a possible resonance near 650 GeV decaying into a pair of Standard Model (SM)-like Higgs bosons or into a SM-like Higgs boson accompanied by a lighter scalar with mass near 95 GeV. In this work, we investigate whether these anomalies can be simultaneously explained within the next-to-two-Higgs-doublet model (N2HDM), an extension of the SM scalar sector featuring two complex Higgs doublets and an additional real singlet. Assuming the existence of a C P -even Higgs state compatible with the 95 GeV excesses (restricted to the γ γ and b b ¯ channels), we analyze the Type-II and Type-Y Yukawa structures, taking the observed 650 GeV resonance to be a C P -even Higgs state. An extensive parameter scan is performed, incorporating the latest constraints from the properties of the observed 125 GeV Higgs boson, direct searches for additional Higgs states, flavor physics data, and electroweak precision observables. Our results show that a heavy C P -even Higgs resonance around 650 GeV, produced predominantly via gluon-gluon fusion and subsequently decaying into a 125 GeV Higgs boson together with another scalar at approximately 95 GeV, can be simultaneously accommodated within both the N2HDM Type-II and Type-Y frameworks in parameter regions that remain consistent with the relevant experimental 2 σ intervals for the reported excesses, once all theoretical and experimental constraints are imposed. This interpretation leads to distinctive and testable predictions for the ongoing LHC Run 3 and the forthcoming high-luminosity LHC phase, in particular through correlated rates in the γ γ b b ¯ , τ + τ − b b ¯ , b b ¯ γ γ , and γ γ τ + τ − final states.