A 95 GeV light Higgs in the top-pair-associated diphoton channel at the LHC in the Minimal Dilaton Model

Whether there is another Higgs scalar remains a fundamental and unresolved question, especially following the discovery of the 125 GeV Higgs boson [1, 2]. Beyond the Standard Model (SM), numerous new physics theories naturally suggest the existence of additional Higgs scalars, alongside a $ 125 \; {\rm{GeV}} $ SM-like Higgs boson. In recent years, more and more data have demonstrated that the discovered $ 125 \; {\rm{GeV}} $ Higgs boson closely resembles the SM, which narrows down the scope of new physics theories, although several theoretical models remain viable. Direct detection of additional scalar particles, particularly in the low-mass region, is crucial for testing these new physics models. Experimentalists have conducted extensive direct searches at colliders such as LEP, Tevatron, and LHC, yielding some intriguing hints.

The first indication of a light Higgs boson emerged in 2003 when the LEP experiment reported a 2.3σ local excess in the $ e^+ e^- \to Z(H\to b \bar{b}) $ channel around 98 GeV [3]. Subsequent LHC experiments provided further clues in 2015, with the CMS collaboration observing a 2.0σ local significance of a diphoton excess at around 97 GeV with 8 TeV 19.7 $ {{\rm{fb}}}^{-1} $ data [4]. This was followed in 2018 by a 2.8σ local (1.3σ global) excess in the 95.3 GeV diphoton invariant mass, using a combination of 13 TeV 35.9 $ {{\rm{fb}}}^{-1} $ data [5]. In 2022, observations by CMS of a 2.6σ local (2.3σ global) $ \tau^+\tau^- $ excess in the 95-100 GeV range with 13 TeV 138 $ {{\rm{fb}}}^{-1} $ data further supported these hints [6]. Most recently, in 2023, CMS and ATLAS reported diphoton excesses around 95 GeV with 2.9σ local (1.9σ global) [7] and 1.7σ local excesses [8], respectively, using 13 TeV 132.2 $ {{\rm{fb}}}^{-1} $ and 140 $ { {\rm{fb}}}^{-1} $ data.

The recurring excesses around 95 GeV suggest the potential existence of an undiscovered light scalar particle at this mass. Numerous recent studies have attempted to account for this 95 GeV excess through various new physics models [9–53].

We focus on the Minimal Dilaton Model (MDM) [54–58], a straightforward extension of the SM. The MDM introduces a singlet scalar field S and a fermion top partner, which is inspired by concepts from strong interaction theory and topcolor theory. This model aims to clarify the mechanisms underlying electroweak symmetry breaking and the role of scale invariance at high energy levels. The singlet does not directly couple with fermions and gauge bosons in SM, thus maintaining its status as a pseudo Nambu-Goldstone particle of broken scale invariance. In scenarios where the top partner has considerable mass, the model essentially converges to the SM plus one scalar. However, a distinct characteristic of the MDM emerges in such contexts: even when the mixing angle between the singlet and the SM-like Higgs boson is zero, the light scalar remains capable of being produced through gluon-gluon fusion processes and subsequently decaying into diphoton channels via loops involving the top partner.

Our previous analyses indicate that the MDM is capable of producing a 95 GeV scalar particle [59]. In this work, we explore the detectability of this scalar at the LHC. To minimize the QCD background, we focus on the production of the $ 95 \; {\rm{GeV}} $ scalar associated with a top pair, decaying to diphoton, with the top quark decaying into $ W^+b $, and the W boson cascade decaying into leptons accompanied by missing transverse energy $ \not E_T $.

The remainder of the paper is structured as follows: In Sec. II, we introduce the Minimal Dilaton Model and delineate the parameter space under theoretical and experimental constraints. In Sec. III, we use Madgraph for detailed collider simulations and then discuss the results. Finally, in Sec. IV, we present our conclusions.

The Minimal Dilaton Model introduces a dilaton-like singlet field, denoted as S, and a vector-like fermion field, represented by T. The Lagrangian of the model is defined as follows [54, 55]:

where $ {\cal{L}}_{{\rm{SM}}}^{\mathrm{without\; }V(H)} $ represents the SM Lagrangian excluding the scalar potential. M is the scale of strong dynamics, and f is the vacuum expectation value (VEV) of the singlet field S. The term $ q_{3 L} $ refers to the quark $ SU(2)_L $ doublet of the third generation. The modified scalar potential incorporating the scalar S, given by $ \tilde{V}(S, H) $, is formulated as [54, 55]:

where $ m_{H} $, $ \lambda_{H} $, $ m_{S} $, $ \lambda_{S} $ and κ are free real parameters. For convenience, we introduce a dimensionless parameter [54, 55]:

Here, v and f represent the VEVs of the Higgs field H and the singlet field S, respectively, with $ v=246 \; {\rm{GeV}} $. The variable $ N_{T} $ denotes the number of the field T, where $ N_{T}=1 $ in the minimal model. To achieve a larger value of η, one may increase $ N_{T} $, thereby avoiding the necessity for an excessively small f.

In the Minimal Dilaton Model, the fields S and H mix, resulting in the formation of two CP-even mass eigenstates: the SM-liked Higgs boson h and the light Higgs s (also called dilaton). The mixing angle, denoted as $ \theta_{S} $, is defined by the following equations:

where $ H^{0} $ represents the neutral component of the Higgs field, and v and f are the VEVs of the Higgs and the singlet field, respectively.

Additionally, the mixing angle $ \theta_{L} $ for the top quark involving the mass eigenstates t and $ t' $ can be defined as follows:

In the conditions where $ m_{t'} \gg m_{t} $ and $ \tan \theta_{L} \ll m_{t'}/m_{t} $, the normalized couplings of the Higgs boson h and the light higgs s can be derived as [54, 55]:

where V represents either the $ W^{\pm} $ or Z boson, and f signifies the fermions, excluding the top quark. Here, $ C_{hVV}/SM $ and $ C_{hff}/{SM} $ represent the normalized couplings of the Higgs boson h to vector bosons and fermions, respectively, relative to their SM predictions. Similarly, $ C_{sVV}/{SM} $ and $ C_{sff}/{SM} $ denote the normalized couplings of the light Higgs s to vector bosons and fermions. The normalized couplings between the light Higgs s and the top quark states t and top partner $ t' $ are specified as follows [54, 55]:

Here, the terms $ C_{st\bar{t}}/SM $ and $ C_{st'\bar{t'}}/SM $ represent the normalized couplings of the light Higgs s to the top quark pair and the top partner pair, respectively. When the mass of the light Higgs, $ m_s $, is set to 95 GeV, the normalized loop-induced couplings with gluons and photons can be approximated as follows [59]:

Here, $ C_{sgg}/SM $ and $ C_{s\gamma\gamma}/SM $ denote the normalized couplings of the light Higgs s to gluons and photons, respectively.

The free parameters in this model are η, $ m_s $, $ m_{t'} $, $ \tan\theta_S $, and $ \sin\theta_L $.

In this work, we will investigate the possibility of searching for the light Higgs in $ pp \to t \bar{t}s $ process at the 14 TeV HL-LHC. In this model, the mass of the top partner $ m_{t'} $ is assumed to be very large to satisfy constraints derived from LHC observations, and the mass of the light Higgs $ m_{s} $ is fixed at 95 GeV to be confronted with the CMS result.

From our previous studies [59], which aimed to meet experimental and theoretical constraints such as vacuum stability, avoidance of Landau poles, and results from low-mass Higgs or resonance searches at LEP, Tevatron, and LHC, we have carefully examined the parameter space. Accordingly, we will consider such a parameter space for all constraints and the existence of a 95 GeV light Higgs:

In the MDM, the tree-level Feynman diagrams for the process $ pp \to t\bar{t}s $ at the LHC is illustrated in Fig. 1. This diagram captures the mechanism of light Higgs s radiation off top quarks, a pivotal process for probing the Higgs sector postulated by the MDM. The depicted process bears a close resemblance to the well-known $ t\bar{t}H $ process, with the primary distinction being the substitution of the Higgs-top coupling $ C_{htt} $ with the light Higgs-top coupling $ C_{stt} $. Given this similarity, the cross-section for the $ pp \to t\bar{t}s $ process can be derived analogously to that of the $ t\bar{t}H $ process. Specifically, the cross-section is calculated by modifying the interaction vertex to incorporate the $ C_{stt} $ coupling. So the cross section can be approximately expressed as follows:

Figure 1.  Feynman diagrams for $ pp \to t \bar{t} s $ in the Minimal Dilaton Model (MDM) at the tree level.

The Next-to-Leading Order (NLO) data for $ \sigma(p p \to t \bar{t} H)|^{SM}_{m_H=95 \; {\rm{GeV}}} $ is found to be $ 1268 \; {\rm{fb}} $ [60]. By applying Eq.(10), the resulting cross section as a function of $ \sin\theta_S $ is illustrated in Fig. 2. It is evident that the cross section $ \sigma(p p \to t \bar{t} s) $ primarily depends on $ \sin\theta_S $, and the sign of $ \sin\theta_L $ does not in fact affect the result. This dependency arises because $ |\sin\theta_L| \ll 1 $; in Eq.(10), the term $ \sin^2\theta_L $ can be neglected. Consequently, the coupling ratio $ (C_{st\bar{t}}/SM)^2 $ approximates to $ \cos^4\theta_L\sin^2\theta_S $, showcasing a quadratic relationship with $ \sin\theta_S $.

Figure 2.  (color online) The cross section of the $ p p \to t \bar{t} s $ process versus $ \sin\theta_S $, with $ \eta = 0.2\; (\rm left),\; 0.6\; (middle), \; 1.2\; (right) $, respectively. The blue, orange, and green lines represent $ |\sin\theta_L| = 0.2, 0.1, 0 $, respectively.

The light Higgs is mainly decay into $ b\bar{b} $, $ c\bar{c} $, $ \tau^+\tau^- $, $ W W^* $, $ ZZ^* $, $ gg $, $ \gamma\gamma $. The light Higgs decay branching ratio can be calculated from the SM Higgs decay information:

where the $ xx $ represent the $ b\bar{b} $, $ c\bar{c} $, $ \tau^+\tau^- $, $ W W^* $, $ ZZ^* $, $ gg $ and $ \gamma\gamma $, and the $ \phi = h, s $ in MDM. $ \Gamma_{tot}^{SM} $ and $ \Gamma_{tot}^{\phi} $ are the total decay widths of the 95 GeV SM Higgs H and the light Higgs s, respectively. And the $ \Gamma_{tot}^{\phi} $ can be writen as:

where the 95 GeV SM Higgs decay information is given as [60]:

In Fig. 3, we plot the decay branching ratios of the light Higgs boson s as a function of $ \sin\theta_S $. Our analysis confirms that the decay branching ratios are not significantly influenced by the $ \sin\theta_L $ parameter; therefore, we set $ \sin\theta_L=0 $. It is observed that for samples where $ \sin\theta_S>0 $, the branching ratio $ Br(s\to\gamma\gamma) $ is greater than for those with $ \sin\theta_S<0 $. Additionally, the minimum value of $ Br(s\to\gamma\gamma) $ is influenced by the parameter η. Specifically, when η is set to 0.2, 0.6, and 1.2, the minimum values of $ Br(s\to\gamma\gamma) $ correspond to $ \sin\theta_S $ values of 0.1, 0.2, and 0.4, respectively.

Figure 3.  (color online) Decay branching ratios of the light Higgs boson s as a function of $ \sin\theta_S $ in the MDM, with $ \sin\theta_L = 0 $ and $ \eta = 0.2\; (\rm left),\; 0.6\; (middle), \; 1.2\; (right) $, respectively. The decay channels $ xx $ include $ b\bar{b} $, $ c\bar{c} $, $ \tau^+\tau^- $, $ W W^* $, $ ZZ^* $, $ gg $, and $ \gamma\gamma $.

In Fig. 4, we present the cross section of the $ pp \to t\bar{t}(s \to \gamma\gamma) $ process as a function of $ \sin\theta_S $. The cross section $ \sigma(pp \to t\bar{t}(s \to \gamma\gamma)) $ is calculated as the product of $ \sigma(pp \to t\bar{t}s) $ and $ Br(s\to\gamma\gamma) $, given by:

Figure 4.  (color online) The cross section of the $ p p \to t \bar{t} (s \to \gamma \gamma) $ process versus $ \sin\theta_S $, with $ \eta = 0.2\; (\rm left),\; 0.6\; (middle), \; 1.2\; (right) $, respectively The blue, orange, and green curves represent $ |\sin\theta_L| = 0.2, 0.1, 0 $, respectively.

Our analysis indicates that the cross section $ \sigma(pp \to t\bar{t}(s \to \gamma\gamma)) $ is primarily affected by $ \sin\theta_S $ and η, with the $ \sin\theta_L $ parameter exerting minimal influence. Notably, the cross section becomes significantly small when $ \sin\theta_S \approx 0 $, attributable to $ \sigma(pp \to t\bar{t}s) \approx 0 $ under such conditions. Furthermore, for samples where $ \sin\theta_S<0 $, the cross section $ \sigma(pp \to t\bar{t}(s \to \gamma\gamma)) $ is smaller compared to those with $ \sin\theta_S>0 $, due to the correlation of $ \sin\theta_S<0 $ with lower $ Br(s\to\gamma\gamma) $ values.

To assess observability, we conduct collider Monte Carlo simulations to explore the sensitivity of the $ pp \to t\bar{t}s $ process at the 14 TeV LHC, focusing on the decay channel:

At tree level, the signal produces two top quarks and a light Higgs boson in the final state. The top quarks predominantly decay into b quarks and W bosons. Considering the light Higgs boson's decay into a diphoton, the final state consists of two b quarks, two W bosons, and two photons ($ WWbb\gamma\gamma $), manifesting as a narrow diphoton resonance centered at the light Higgs mass of 95 GeV.

Backgrounds can be categorized based on their behavior in diphoton systems as either resonant or non-resonant. The resonant backgrounds, including $ tth $ and $ tjh $, feature a Higgs boson that can decay into a diphoton in the final state, with the top quark decaying into a b quark and a W boson. In the case of $ tjh $, the light jet can be misidentified as a b jet. The non-resonant backgrounds consist of processes such as $ bb\gamma\gamma $, $ tt\gamma\gamma $, $ tj\gamma\gamma $, $ tt\gamma $, and $ Wjj\gamma\gamma $, with $ tj\gamma\gamma $ and $ Wjj\gamma\gamma $ involving a light jet being misidentified as a b jet.

Signal and background events at the parton level were generated using $\mathrm{MadGraph5_aMC@NLO_v3.4.1}$ [61]. Subsequently, $\mathrm{PYTHIA_v8.2}$ [62] was employed for particle decay, parton shower, and hadronization processes. The simulation of the detector response was conducted with $\mathrm{Delphes_v3.4.1}$ [63], adopting the default values for mistagging efficiencies. Jet clustering was performed with the $\mathrm{anti-kT}$ algorithm [64], and the $\mathrm{NNPDF23LO1}$ [65] parton distribution functions were selected for our simulation. Event analysis was carried out with $\mathrm{MadAnalysis5_v1.9.60}$ [66].

To determine appropriate kinematic cuts, Fig. 5 displays normalized distributions for both the signals and SM backgrounds. The selection criteria for particles and jets are based on quality metrics for leptons, jets, and photons, specified as follows:

Figure 5.  (color online) The normalized distributions of nine kinematic variables, including $ p_T^{\gamma_1} $ (upper left), $ p_T^{\gamma_2} $ (upper middle), $ p_T^{\ell_1} $ (upper right), $ \not E_T $ (middle left), $ \Delta R_{j_1,j_2} $ (middle), $ \Delta R_{\gamma_1, \gamma_2} $ (middle right), $ M [\gamma_1\gamma_2] $ (lower left), $ M_T[\gamma_1\gamma_2] $ (lower middle), and $ M_T[\gamma_1\gamma_2 b_1] $ (lower right), for the signals and SM backgrounds at the 14 TeV LHC.

To enhance the signal-to-background ratio for the $ WWbb\gamma\gamma $ signal amidst the QCD background, it is crucial to include at least one leptonic decay of a W boson into a lepton and a neutrino, manifesting as a lepton with missing transverse energy ($ \not E_T $). To this end, we apply a series of cuts in four stages:

(1) Basic Cut: Selection criteria include at least one lepton ($ N(\ell) \geq 1 $), two photons ($ N(\gamma) \geq 2 $), and one b-jet ($ N(b) \geq 1 $), alongside thresholds for missing transverse energy and spatial separation:

with specified $ p_T $ criteria for photons and leptons:

(2) Diphoton Invariant Mass Cut: Concentrating on the diphoton invariant mass narrowing down around the light Higgs mass ($ 95 \; {\rm{GeV}} $) to improve signal clarity:

(3) Transverse Invariant Mass of Diphoton Cut: Refining the selection further based on the diphoton transverse invariant mass to diminish background:

(4) Transverse Invariant Mass of Diphoton with a b-Jet Cut: Narrowing down events based on the combined transverse invariant mass of the diphoton plus a b-jet:

We apply these cuts to our benchmark points and backgrounds to assess the signal efficiency. The outcomes of these sequential cuts are presented in Table 1. Notably, the predominant background originates from the $ tt\gamma\gamma $ process. After implementing the basic cut, the $ bb\gamma\gamma $ background becomes negligible, attributed to the inclusion of a lepton and missing transverse energy ($ \not E_T $) criteria. Upon applying the cut on the invariant mass around 95 GeV, backgrounds from SM Higgs processes, specifically $ tjh $ and $ tth $, are effectively eliminated. The subsequent cut on the transverse invariant mass of the diphoton system further excludes the background involving missed tagged photons, namely $ tt\gamma $. Finally, the cut on the transverse invariant mass of the diphoton system with a b-jet successfully removes the incorrectly tagged b-jet background from the $ Wjj\gamma\gamma $ process.

Cuts	Signal ($\sigma \times 10^{-3} \;{\rm{fb}} $)		Background ($ \sigma\times 10^{-3} \;{\rm{fb}}$)
	$ tts $		$ bb\gamma\gamma $	$ tj\gamma\gamma $	$ tjh $	$ tt\gamma $	$ tt\gamma\gamma $	$ tth $	$ Wjj\gamma\gamma $
Initial	500		7646000	17510	78	5120000	11180	563	218000
Basic cuts	15.65		0.00	35.02	0.97	1751.04	116.94	22.70	152.60
$ M[\gamma_1 \gamma_2] $	15.00		0.00	1.40	0.00	102.91	6.48	0.02	8.72
$ M_T[\gamma_1 \gamma_2] $	14.40		0.00	1.40	0.00	0.00	5.14	0.02	4.36
$ M_T[\gamma_1 \gamma_2 b] $	12.25		0.00	0.70	0.00	0.00	4.25	0.02	0.00

Table 1.  The cut flow of the cross sections, for a benchmark point with $ \eta=0.2, \sin\theta_S=-0.6, \sin\theta_L = 0.15 $

In Fig. 6, we illustrate the required luminosity to achieve specified statistical significance levels on the η versus $ \sin\theta_S $ plane, assuming $ \sin\theta_L = 0 $. Our analysis confirms that $ \sin\theta_L $ has a minimal impact on the outcomes. To estimate the signal significance, we employ the Poisson formula [67]:

Figure 6.  (color online) Required luminosity to achieve specified statistical significance of $ 2\sigma $ (left), $ 3\sigma $ (middle), and $ 5\sigma $ (right), respectively, on the η versus $ \sin\theta_S $ plane, with $ \sin\theta_L = 0 $. The red curve represents the luminosity threshold of $ L = 300 \, {\rm{fb}}^{-1} $, and the green curve denotes $ L = 3000 \, {\rm{fb}}^{-1} $.

where S and B represent the signal and background cross sections, respectively, and L is the integrated luminosity. The results demonstrate how the statistical significance ($ {\cal{S}} $) values vary with integrated luminosity L. Statistical significance levels of $ 2\sigma $, $ 3\sigma $, and $ 5\sigma $ are illustrated from left to right, respectively. The red line indicates a luminosity threshold of $ L = 300 \, {\rm{fb}}^{-1} $, while the green line marks $ L = 3000 \, {\rm{fb}}^{-1} $. From the left figure, it is evident that regions where $ \eta+4\sin\theta_S<-2 $ or $ \eta+4\sin\theta_S>2 $ can be covered at the $ 2\sigma $ level with an integrated luminosity of $ L = 300 \; {\rm{fb}}^{-1} $. Similarly, regions where $ \eta+4\sin\theta_S<-1.2 $ or $ \sin\theta_S>0.2 $ can be covered at the $ 2\sigma $ level with $ L = 3000 \, {\rm{fb}}^{-1} $. In the middle plane, the region where $ \sin\theta_S>0.4 $ can reach a $ 3\sigma $ significance with $ L = 300 \, {\rm{fb}}^{-1} $. Additionally, areas defined by $ \eta+3.3\sin\theta_S<-1 $ or $ \sin\theta_S>0.2 $ can achieve a $ 3\sigma $ level with $ L = 3000 \, {\rm{fb}}^{-1} $. The right plane demonstrates that regions satisfying $ \eta+3\sin\theta_S<-1.2 $ or $ \sin\theta_S>0.25 $ can achieve a $ 5\sigma $ level at an integrated luminosity of $ L = 3000 \; {\rm{fb}}^{-1} $.

In this study, we explore the detection prospects of a $ 95 \; {\rm{GeV}} $ light Higgs boson in the Minimal Dilaton Model through the top-pair-associated diphoton process $ pp \rightarrow t\bar{t}s (\rightarrow \gamma\gamma) $ at the $ 14 \; {\rm{TeV}} $ LHC. Monte Carlo simulations have been conducted to analyze the signal and background. The final signal state considered is $ WWbb\gamma\gamma $, necessitating at least one leptonic decay of W boson to reduce the QCD background. The production cross section $ \sigma(pp \rightarrow t\bar{t}s) $ demonstrates a quadratic relationship with $ \sin\theta_S $. It's noteworthy that the light Higgs decay branching ratio $ Br(s \rightarrow \gamma\gamma) $ is lower for $ \sin\theta_S < 0 $ compared with that of $ \sin\theta_S > 0 $. Following a meticulous cut-based selection process, we identified that the predominant background is attributed to the SM $ tt\gamma\gamma $ process.

Our results indicate that areas where $ \eta+4\sin\theta_S < -2 $ or $ \eta+4\sin\theta_S > 2 $ can be covered at the $ 2\sigma $ level with an integrated luminosity of $ L = 300 \, {\rm{fb}}^{-1} $. Similarly, regions defined by $ \eta+4\sin\theta_S < -1.2 $ or $ \sin\theta_S > 0.2 $ can be covered at the $ 2\sigma $ level with $ L = 3000 \, {\rm{fb}}^{-1} $. Moreover, regions of $ \eta+3\sin\theta_S < -1.2 $ or $ \sin\theta_S > 0.25 $ can reach a significance level of $ 5\sigma $ with an integrated luminosity of $ L = 3000 \, {\rm{fb}}^{-1} $.

Overall, for areas where $ |\sin\theta_S|>0.2 $, indicative of significant Higgs-dilaton mixing, exclusion or detection at future colliders is feasible with luminosities of $ L = 300 \, {\rm{fb}}^{-1} $ or $ L = 3000 \, {\rm{fb}}^{-1} $. Conversely, the detection of areas where $ |\sin\theta_S|<0.2 $, indicative of minimal Higgs-dilaton mixing, proves challenging to detect through this process, even at an integrated luminosity of $ L = 3000 \, {\rm{fb}}^{-1} $. We aim to further investigate this parameter space in future research, focusing on strategies to identify such scenarios at upcoming collider experiments.