[a,b]Srimoy Bhattacharya

Emerging Narrow Resonance at 152 GeV

   Mukesh Kumar    Rachid Mazini    Bruce Mellado
Abstract

The discovery of the Higgs boson at the LHC completed the Standard Model (SM), yet the possibility of additional scalars remains open, provided their contributions to electroweak symmetry breaking are sufficiently small. Recent analyses of LHC data have revealed statistically significant anomalies in multi-lepton final states - events characterized by multiple leptons, missing transverse energy, and (bb-)jets. These anomalies provide intriguing hints of physics beyond the SM. In this work, we present the signature of growing excesses for a new scalar resonance with a mass of 152±1152\pm 1 GeV, observed in the γγ\gamma\gamma, ZγZ\gamma, and WWWW channels. The combined global significance reaches a level that points toward the growing signature of this resonance. The findings align with a simplified model in which a heavy scalar boson decays into two lighter scalars, thus providing a consistent framework explaining the observed multi-lepton anomalies.

These results significantly advance the search for new scalar bosons at the electroweak scale. Future investigations, including precision studies with upcoming HL-LHC data, will be crucial for confirming the nature of this resonance and exploring its implications for extending the SM.

1 Introduction

The Standard Model (SM) of particle physics provides a comprehensive framework describing the fundamental particles and their interactions [ParticleDataGroup:2022pth]. The discovery of the Higgs boson in 2012 by ATLAS and CMS [Aad:2012tfa, Chatrchyan:2012ufa] confirmed the mechanism of electroweak symmetry breaking. Measurements of the 125125 GeV boson show good agreement with SM expectations [Langford:2021osp]. Nonetheless, the SM leaves several key phenomena unexplained – dark matter (DM), neutrino masses, and the baryon asymmetry – while facing theoretical tensions such as the hierarchy problem. These issues suggest the existence of an extended scalar sector, which could also provide explanations for the observational evidence for DM and neutrino masses. Interestingly, over the past decade, several statistically significant anomalies have appeared in LHC data involving final states with multiple leptons (electrons and/or muons) in conjunction with missing energy and with or without (bb-)jets  [Fischer:2021sqw, Crivellin:2023zui], which indicate that the anomalies are not isolated but statistically consistent across channels and datasets. These “multi-lepton anomalies” may hint at additional scalar bosons at the electroweak scale.

Within the framework of a simplified model [vonBuddenbrock:2016rmr, Buddenbrock:2019tua], the multi-lepton anomalies can be explained by the gluon-fusion production of a heavy scalar HH (mH270m_{H}\approx 270 GeV) decaying dominantly into two lighter scalars SS (one off-shell), each primarily decaying to WW bosons and bb-jets. A mass of 150±5150\pm 5 GeV for SS was inferred from the invariant mass of electron–muon pairs in non-resonant WWWW events [vonBuddenbrock:2017gvy], motivating targeted searches for an associated production of narrow resonances around 150 GeV.

2 Multi-lepton anomalies and model motivation

Multi-lepton anomalies first appeared in LHC Run 1 data, including excesses in Higgs pTp_{T} distributions [ATLAS:2014yga, ATLAS:2014xzb, CMS:2015qgt, CMS:2015hja], di-Higgs searches [ATLAS:2015sxd, CMS:2014ipa, CMS:2015uzk, CMS:2014jkv], and associated Higgs production with top quarks [ATLAS:2014ayi, ATLAS:2015xdt, ATLAS:2015utn, CMS:2014tll]. These features motivated a simplified model, as discussed earlier, predicting final states with ++EmissT\ell^{+}\ell^{-}+E^{T}_{\text{miss}}, bb-jets, and photons. Subsequent studies [Buddenbrock:2019tua, Hernandez:2019geu] confirmed that these anomalies persisted in Run 2 data. More recently, excesses in triboson (VVVVVV) final states with V=W,ZV=W,Z have been reported, notably in VVZVVZ, tWZtWZ, and WWWWWW channels [ATLAS:2024nab, CMS:2025hlu, CMS:2020hjs, ATLAS:2022xnu, CMS:2025tre],which aligns with a Real Higgs Triplet model with hypercharge Y=0Y=0 [Ashanujjaman:2024lnr]. These correlated features strengthen the case for an extended Higgs sector involving a new scalar around 150 GeV. Figure 1 shows how anomalies accumulated chronologically, converging toward a predicted scalar around 150±5150\pm 5 GeV.

Refer to caption
Figure 1: Timeline summarizing the evolution from early multi-lepton anomalies to the predicted narrow resonance near 150 GeV [Crivellin:2023zui].

A systematic search for a resonance near mS150±5GeVm_{S}\approx 150\pm 5~\text{GeV} was performed using data from ATLAS and CMS corresponding to proton–proton collisions at s=13TeV\sqrt{s}=13~\text{TeV}. Rather than executing a broad, model-independent scan of all available final states, we adopted a targeted strategy guided by the simplified model introduced in Ref. [vonBuddenbrock:2016rmr]. This approach significantly reduces the “look-elsewhere” effect and potential biases associated with cherry-picking statistically interesting regions in a large dataset. The simplified model assumes that a heavy scalar HH with a mass of mH270GeVm_{H}\simeq 270~\text{GeV} is produced dominantly via gluon fusion and decays into two lighter scalars SS, one of which can be off-shell (HSS()H\to SS^{(*)}). Each SS behaves as a SM like scalar, decaying primarily to WW bosons or bb-jets, while also allowing small branching ratios into photons and invisible states. This topology naturally predicts associated production of photons, leptons, and jets in the final state.

3 Methodology

The channels, predicted by the simplified model used to interpret the corresponding experimental results, are organized into three groups according to the time period of analysis and the evolution of the dataset:

  • Group 1 (2021): Initial studies focused on diphoton (γγ\gamma\gamma) and ZγZ\gamma final states with missing energy, bb-jets, or additional leptons. The first combined analysis [Crivellin:2021ubm] yielded a global (local) significance of 3.9σ3.9\sigma (4.3σ4.3\sigma) near 151.5GeV151.5~\text{GeV}.

  • Group 2 (2023): Extended the search to include γγ+(1+1b-jet)\gamma\gamma+\geq(1\ell+1b\text{-jet}), γγ+4jets\gamma\gamma+\geq 4~\text{jets}, and WW+ETmissWW^{*}+E_{T}^{\text{miss}} final states [ATLAS:2023omk, CMS:2022uhn, ATLAS:2022ooq], further strengthening the case for a common underlying resonance.

  • Group 3 (2025): The most recent results [ATLAS:2024lhu] explored the γγ+τ\gamma\gamma+\tau and γγ+2(,τ)\gamma\gamma+2(\ell,\tau) channels.

Each analysis utilized the sidebands of the SM Higgs searches to probe the invariant mass spectra between 140 and 155GeV155~\text{GeV}. The local pp-values were extracted from these spectra using the standard background models provided by the experimental collaborations, ensuring a consistent treatment of statistical uncertainties. For the WW+ETmissWW^{*}+E_{T}^{\text{miss}} final state, the results were derived by recasting the data from Ref. [Coloretti:2023wng] in the context of the simplified model.

For each channel, we determined the best-fit signal yield and corresponding local pp-value as a function of the hypothesized SS-boson mass. The independent pip_{i} values were then combined using Fisher’s method [fisher1925], which defines a global test statistic as:

χ2n2=2i=1nln(pi),\textstyle\chi^{2}_{2n}=-2\sum_{i=1}^{n}\ln(p_{i}), (1)

where nn is the number of independent channels. The resulting statistic follows a χ2\chi^{2} distribution with 2n2n degrees of freedom. In the final combination, five statistically independent channels were included:

  1. 1.

    Sγγ,ZγS\to\gamma\gamma,\,Z\gamma

  2. 2.

    Sγγ+4jetsS\to\gamma\gamma+\geq 4~\text{jets}

  3. 3.

    Sγγ+(1+1b-jet)S\to\gamma\gamma+\geq(1\ell+1b\text{-jet})

  4. 4.

    SWW+ETmissS\to WW^{*}+E_{T}^{\text{miss}}

  5. 5.

    Correlated Sγγ+(1τ,2τ)S\to\gamma\gamma+(1\tau,2\tau) channels.

The correlation between one- and two-τ\tau final states was accounted for by including the ratio of signal efficiencies obtained from simulations of the HSS()H\to SS^{(*)} process. A trial factor of approximately 3.5 was applied to correct for the look-elsewhere effect across the scanned mass range, consistent with the resolution of the diphoton analyses [Gross:2010qma].

4 Results and discussion

Refer to caption
Figure 2: Combined pp-values as a function of the hypothesized scalar mass. The minimum occurs near 152GeV152~\text{GeV}, corresponding to a global significance of 5.3σ5.3\sigma.

Using the methodology explained above, we first revisited the 2021 combination, incorporating additional γγ+b\gamma\gamma+b-jet [ATLAS:2021ifb] and Z(+)γZ(\to\ell^{+}\ell^{-})\gamma [CMS:2022ahq] channels. The inclusion of these datasets slightly reduced the global significance from 3.9σ3.9\sigma to 3.5σ3.5\sigma at 152GeV152~\text{GeV}. When the new WWWW^{*}, 4j4j, and 1+b1\ell+b-jet channels were included, the combined global (local) significance rose to 4.7σ4.7\sigma (4.9σ4.9\sigma[Bhattacharya:2023lmu]. Finally, integrating the 2025 τ\tau-enriched channels led to a further increase, yielding a global (local) significance of 5.3σ5.3\sigma (5.5σ5.5\sigma) at mS=152GeVm_{S}=152~\text{GeV}, as shown in Fig. 2. The consistent excesses across all independent final states strongly indicate the presence of a narrow resonance, compatible with a SM-like scalar. The robustness of this excess was further tested against variations in background modeling and systematic uncertainties. These were found to have only a marginal impact on the significance, typically at the level of 0.1σ0.1\sigma, comparable to the impact observed in the discovery of the 125 GeV Higgs boson [Aad:2012tfa]. This provides additional confidence that the observed structure is unlikely to originate from modeling artifacts. The coherent appearance of signals across multiple independent final states and decay modes provides the strongest indication yet of a new scalar resonance at 152±1152\pm 1 GeV. Confirming these states would mark a breakthrough beyond the SM, opening a new era in particle physics. A future high-precision e+ee^{+}e^{-} collider [FCC:2018evy, CEPCStudyGroup:2018ghi] will be essential to probe their properties and validate their theoretical interpretation.