Abstract
Cascade decays of new scalars into final states with multiple photons and possibly quarks may lead to distinctive experimental signatures at high-energy colliders. Such signals are even more striking if the scalars are highly boosted, as when produced from the decay of a much heavier resonance. We study this type of events within the framework of the minimal stealth boson model, an anomaly-free \(\text {U}(1)_{Y'}\) extension of the Standard Model with two complex scalar singlets. It is shown that, while those signals may have cross sections that might render them observable with LHC Run 2 data, they have little experimental coverage. We also establish a connection with a CMS excess observed in searches for new scalars decaying into diphoton final states near 96 GeV. In particular, we conclude that the predicted multiphoton signatures are compatible with such excess.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.Avoid common mistakes on your manuscript.
1 Introduction
The scalar sector of the Standard Model (SM) contains one scalar doublet which spontaneously breaks the gauge symmetry via the Brout–Englert–Higgs mechanism [1,2,3], predicting the existence of the so-called Higgs boson. A particle compatible with the SM predictions for the Higgs boson and with a mass of approximately 125 GeV was discovered by the ATLAS and CMS Collaborations [4, 5] at the Large Hadron Collider (LHC), thereby culminating several decades of searches. Although in the SM the scalar content is the minimal one required to break the symmetry, several additional scalar particles may be present in SM extensions. For instance, models with an augmented gauge symmetry require extra scalar fields to break that symmetry down to the SM group \(\text {SU}(3) \times \text {SU}(2)_L \times \text {U}(1)_Y\). Thus, scenarios with heavier scalars decaying into lighter ones (or into pairs of SM weak bosons \(V=W,Z\)) may be naturally envisaged [6] as, for example, in the context of left-right models [7, 8]. When the former are produced from decays of a very heavy resonance, and the latter decay hadronically, the experimental signature is a multi-pronged fat jet. In this case, the heavier decaying scalar has been dubbed as ‘stealth boson’ for its elusive character [9].
In Ref. [10] we have proposed the simplest model that accounts for the said cascade decays – the minimal stealth boson model (MSBM) – in which the SM gauge symmetry is enlarged with an extra \(\text {U}(1)_{Y'}\) coupling to baryon number up to an arbitrary normalisation constant. The scalar sector comprises two complex SM singlets \(\chi _{1,2}\) that, upon \(\text {U}(1)_{Y'}\) symmetry breaking, provide masses to the \(Z'\) boson and to the new fermions. We point out that, in order for the \(Z'\) to decay into scalar pairs, two complex singlets are required, since one of the degrees of freedom is ‘eaten’ by the \(Z'\) boson, and also because the required coupling involves one CP-even and one CP-odd scalar weak eigenstate. To ensure gauge-anomaly cancellation, the SM fermion sector is minimally extended with extra leptons.Footnote 1 In the MSBM one can accommodate sizeable branching ratios (BRs) for the \(H_i \rightarrow \gamma \gamma \) decays of the new scalars since the corresponding one-loop amplitudes may be enhanced by the Yukawa couplings with the new charged leptons. (We denote the four scalar mass eigenstates as \(H_i\), \(i=1-4\), being \(H_1 \equiv H\) the 125 GeV Higgs boson.) This will be shown in Sects. 2 and 3. After writing down the relevant interactions and the \(H_i \rightarrow \gamma \gamma \) decay widths in Sect. 2, we perform a scan over the parameter space of the model in Sect. 3 to compute the BRs for \(H_i \rightarrow \gamma \gamma \), as well as for the other decay modes. We remark that although sizeable BRs for \(H_i\) decays into \(\gamma \gamma \) are naturally accommodated within the MSBM, they are not a straightforward implication of the model.
Besides the \(pp \rightarrow H_i \rightarrow \gamma \gamma \) signals from direct production of the new scalars, the presence of a heavy \(Z'\) resonance opens up the possibility of several conspicuous signals from cascade decays, like those involving collimated photons and/or jets containing hard photons. Their features will be discussed in Sect. 4 where we will also show that, remarkably, those signals have little experimental coverage, in the sense that the efficiency for such signals in current searches is marginal. Direct \(H_i \rightarrow \gamma \gamma \) decays are also interesting by themselves, and may reach detectable levels in some regions of the parameter space. In this regard, it is worth noting that the CMS Collaboration has found [14] an excess in the searches for new scalars h decaying into photon pairs with Run 1 (8 TeV) and Run 2 (13 TeV) data, which reaches a statistical significance of \(2.8\sigma \) at \(M_h = 95.3~ \text {GeV}\). The ATLAS Collaboration has also looked for such process with a larger luminosity of \(80 \, {\hbox {fb}}^{-1}\) [15] at 13 TeV, without finding a significant deviation from the background-only expectation at the same mass. However, as noted in Ref. [16], the limits from the ATLAS Collaboration are weaker and are not able to exclude the CMS excess as due to a new particle. Whether this excess can be due to a stealth boson (dominantly decaying into lighter scalars) will be examined in Sect. 5. For completeness, we also explore the possibility that the same scalar decays mostly into \(b {{\bar{b}}}\). In such case, an excess observed at the Large Electron-Positron (LEP) collider around the same mass [17] can also be accommodated. This possibility is further explored in Appendix A. To conclude, in Sect. 6 we discuss our results in light of current and future experimental searches.
2 \(H_j \rightarrow \gamma \gamma \) decays in the MSBM
The MSBM [10] extends the SM scalar content with two complex scalar fields \(\chi _1\) and \(\chi _2\), which are SM singlets and are equally charged under a gauged U(1)\(_{Y'}\) symmetry. The anomaly-cancellation conditions require extra matter fields which, in this work, we consider to be a set of two vector-like lepton singlets \(N'_{1}\), \(E'_{1}\) and a doublet \((N'_{2},E'_{2})\) (we use primes on fermion weak eigenstates in order to distinguish them from unprimed mass eigenstates). The SM and \({\hbox {U(1)}}_{Y'}\) hypercharges Y and \(Y^\prime \) of the various fields are shown in Table 1, where \(Q_{iL}\) and \(l_{iL}\) denote the left-handed (LH) quark and lepton \(\text {SU}(2)_L\) doublets of the SM, respectively. The right-handed (RH) quarks and lepton singlets are \(u_{iR}\), \(d_{iR}\) and \(e_{iR}\). A model with the same fermion content and hypercharge assignments as those in Table 1, but with only one scalar singlet, has been previously considered [11, 12]. However, the decays of the \(Z'\) boson into two scalars require two or more singlets, opening also the possibility of pure scalar cascade decays, and implying a richer phenomenology.
The most general scalar potential invariant under SU(2)\(_L\times \)U(1)\(_Y \times \)U(1)\(_{Y'}\) is \(V = V_{Z_2} + V_{\not Z_2}\), with
where \(\Phi =(\phi ^+\;\phi ^0)\) is the SM Higgs doublet. The terms in \(V_{Z_2}\) (\(V_{\not Z_2}\)) conserve (break) a \(Z_2\) symmetry under which only \(\chi _2\) transforms non trivially as \(\chi _2\rightarrow - \chi _2\). While \(m_0^2\), \(m_{11}^2\), \(m_{22}^2\), \(\lambda _{0-3}\) and \(\lambda _{5,6}\) are real, \(m_{12}^2\), \(\lambda _{4}\) and \(\lambda _{7-9}\) can be, in general, complex. We define
where \(\rho _i\) and \(\eta _i\) are real fields and
with \(v = 246\,\mathrm{GeV}\). As usual, we will also use
A detailed analysis of the scalar potential and scalar mass spectrum can be found in Ref. [10]. Upon U(1)\(_{Y'}\) symmetry breaking, the neutral \(Z'\) gauge boson acquires the mass
where \(g_{Z'}\) and \(Y'_{\chi }\) are the U(1)\(_{Y'}\) coupling constant and \(\chi _{1,2}\) hypercharge, respectively. All \(Z'\)-scalar interactions can be found in Ref. [10].
There are two would-be Goldstone bosons in the model, namely \(G_1^0 = \eta _0\) and \(G_2^0 = \cos \beta \, \eta _1 + \sin \beta \, \eta _2\). The orthogonal state \(A^0 = - \sin \beta \, \eta _1 + \cos \beta \, \eta _2\) is CP-odd, being a mass eigenstate in case the scalar potential parameters are real and \(\varphi = 0\). In general, \(H'_i = (\rho _0 \; \rho _1 \; \rho _2 \; A^0)\) are related to the mass eigenstates \(H_i = (H_1 \; H_2 \; H_3 \; H_4)\) by
where O is an orthogonal \(4\times 4\) real matrix parameterised in terms of \(2\times 2\) rotations as \(O = {\widehat{O}}_{34} {\widehat{O}}_{24} {\widehat{O}}_{14} {\widehat{O}}_{23} {\widehat{O}}_{13} {\widehat{O}}_{12}\), where \({\widehat{O}}_{kl}\) corresponds to a rotation in the (k, l) plane by an angle \(\theta _{kl}\). We identify the 125 GeV SM Higgs boson H as being \(H_1\). The scalar interactions with the \(Z'\) boson field originate from the term
which in the mass-eigenstate basis reads
with \(i < j\) and mixing factors
Notice that \(R_{ij}\) are anti-symmetric and therefore \(R_{ii} = 0\), reflecting the fact that \(Z' \rightarrow H_i H_i\) is forbidden. Also, it can be shown that \(\sum _{i<j} R_{ij}^2 = 1\) due to the orthogonality of the mixing matrix O.
2.1 Scalar-fermion interactions in the mass basis
Given their relevance for the computation of scalar diphoton decay amplitudes, we now obtain the scalar-fermion interactions in the mass basis. The masses of the new fermions and their interactions with scalars are determined by the gauge-invariant Lagrangian
where the Yukawa couplings \(x_i\), \(y_{ i}^F\) and \(w_{i}^F\) (with \(i=1,2\), \(F=E,N\)) are general complex numbers. Taking into account Eq. (3), the corresponding mass terms for the new charged and neutral fermions in the interaction basis \(F'_{L,R}=(F'_1\;F'_2)_{L,R}^T\) are
with
We define the field rotations to the mass-eigenstate basis \(F_{L,R}=(F_1\;F_2)_{L,R}^T\) as
Here, \(U_{L,R}^F\) are \(2\times 2\) unitary complex matrices and the masses of the new leptons \(m_{F_i}\) are real and positive. To compute \(U_{L,R}^F\), one defines the Hermitian matrices \(H_L^F=M_F M_F^{\dag }\) and \(H_R^F=M_F^{\dag }M_F\) diagonalised by \(V_{L,R}^{F\dagger }\,H_{L,R}^F\, V_{L,R}^F=\mathrm{diag}(m_{F_1}^2,m_{F_2}^2)\), being \(V_{L,R}^F\) unitary matrices given by:
The mixing angles \(\theta _{L,R}^F\) and CP phases \(\delta _{L,R}^F\) satisfy
Since, in general, \(V_{L,R}^F\) diagonalise \(M_F\) up to diagonal phases, an additional phase transformation must be performed in either the LH or RH fields to express the interactions in the physical mass basis. Namely, we follow the convention
Considering Eqs. (6), (10) and (13), we can write the scalar-fermion interactions in the mass-eigenstate basis as:
with \((j,k)=1,2\), \(a=1,\dots ,4\), and the coefficients
with \(y_2^F \equiv y_2\), \(x_2^F \equiv x_2\). In terms of \(F = F_L+F_R\), Eq. (17) reads
again with \((j,k)=1,2\), \(a=1,\dots ,4\) and the coefficients \(C_{jka}^F\), \(D_{jka}^F\) given by:
2.2 \(H_j \rightarrow \gamma \gamma \) decay widths
The Feynman diagrams which contribute at one-loop level to \(H_j \rightarrow \gamma \gamma \) decays are shown in Fig. 1. In diagram (a) f is any electrically-charged fermion in the model, namely \(f=e_i,u_i,d_i,E_i\). Due to scalar mixing, all charged fermions enter the \(H_j\) diphoton decay loop, although in practice some of the contributions are suppressed due to experimental constraints on the SM Higgs couplings. Diagrams (b) and (c) stand for the W-boson contributions. It is convenient to write the \({\bar{f}} f H_j\) interactions of Eq. (19) in the form:
with \(j=1,\dots ,4\). The scalar and pseudoscalar couplings \(a_f\) and \(b_f\), respectively, are given by
for SM fermions, where O is the scalar mixing matrix defined in Eq. (6). Notice that \(O_{1j}\) is the admixture between the \(\text {SU}(2)_L\) doublet \(\Phi \) and the j-th scalar eigenstate; the SM fermions do not couple to the \(\text {SU}(2)_L\) singlets. For the new leptons,
being the coefficients \(C_{kkj}^E\) and \(D_{kkj}^E\) those of (20). As for the \(WWH_j\) coupling, we have:
where g is the SM SU(2) gauge coupling and \(m_W\) the W boson mass.
One-loop diagrams contributing to \(H_j\rightarrow \gamma \gamma \). In the fermion loop (a) f stands for any electrically-charged fermion (see Table 1)
Due to electromagnetic gauge invariance, the tensor form of the \(H_j\rightarrow \gamma \gamma \) amplitude can be generically expressed as [18]
being p and q the momenta of the outgoing photons, while \(\varepsilon (p)\) and \(\varepsilon (q)\) denote the corresponding polarisation vectors. The quantities \(S_j^X\) and \({\tilde{S}}_j^X\) with \(X=f,W\), are scalar form factors computed by considering the relevant contributions. In the present case, these stem from the one-loop diagrams shown in Fig. 1, for which
Here, \(\alpha =e^2/(4\pi )\) is the electromagnetic fine-structure constant, \(Q_{f}\) is the fermion electric charge, \(N_c^{f}\) is the number of colours of f and \(\tau _X=M_{H_j}^2 / 4m_X^2\). The function \(g(\tau )\) is given by [19]
As usual, the decay widths are obtained using
where \(|\vec {p}\,|=M_{H_j}/2\) and a 1/2 factor has been included for identical particles in the final state.
3 Parameter-space analysis
In order to investigate the possible multiphoton signals from \(H_j\) decays in the MSBM, we proceed with a parameter-space analysis. For simplicity, we restrict ourselves to the case with \(\lambda _{7-9}=0\) in Eq. (1), which corresponds to having the \(Z_2\) symmetry softly broken by the term \(m_{12}^2\ne 0\). The interested reader is addressed to Ref. [10] for details on the full analysis of the scalar potential, and on the reconstruction of its parameters in terms of the scalar masses \(M_{H_i}\) and mixings \(O_{ij}\). Here, we briefly summarise the general procedure. We will consider the following values for scalar masses and for the \(Z'\) mass and coupling:
The TeV-scale mass for \(Z'\) boson is chosen so as to have the scalars from its decay highly boosted. For better comparison with previous work [10] we choose the value of 2.2 TeV. The chosen coupling \(g_{Z'} Y_q'\) fulfills the limits from dijet production (see next section), while \(M_{H_3}\) is set to 96 GeV in order to establish a connection with Sect. 5, where we investigate if a CMS excess can be due to one of the scalars present in the spectrum. (Of course, the signals studied in Sect. 4 may take place for a wider range of masses.) We choose \(M_{H_4}\) of the same order. The mass of the lightest scalar \(H_2\) is an output parameter which in our scan we require to be within the interval [20, 40] GeV in order to allow \(H_3 \rightarrow H_2 H_2\) and \(H_4 \rightarrow H_2 H_2\). The values of the six angles \(\theta _{ij}\) defining the scalar mixing matrix O in Eq. (6), and of the scalar potential parameter \(\lambda _2\), are randomly varied in the intervals \([0,2\pi ]\) and \([-1,1]\), respectively. The latter range is enough to obtain the predictions later presented.
These inputs allow us to determine \(\tan \beta \), \(m_{H_2}\), \(\lambda _{0-6}\) and \(m_{12}^2\). As done in Ref. [10], we will only keep those points which lead to a global minimum of the potential with \(v,u_{1,2}\ne 0\) and \(\varphi =0\). Notice also that, since we are considering the case with a softly broken \(Z_2\) symmetry under which \(\chi _2 \rightarrow -\chi _2\), the couplings \(x_2\) and \(x_1^{N,E}\) in Eq. (10) vanish. We have verified that this simplification does not affect the predictions for the \(H_i\) decays, as there is still enough freedom for the coefficients in Eqs. (18). Fermion mixing plays little role in our study and, given that \(v\ll u \simeq 1.6\,\mathrm{TeV}\), we will neglect it in Eq. (12) implying \(\theta _{L,R}^{F}=0\) and, thus, \(U_{L,R}^F={\mathbb {1}}_{2\times 2}\). The only relevant couplings among the new fermions and the scalars are then \(y_1^{E,N}\) and \(y_2\), which we randomly vary in the perturbative interval \([0,4\pi ]\). The resulting heavy leptons masses range from 300 GeV to 10 TeV. The \(\lambda \) couplings are small because the four scalars are relatively light, \(m_{H_i} \ll u\), which guarantees perturbativity.
Experimental constraints on the signal-strength parameters for each production and decay mode combination of the SM Higgs H must be taken into account. Since there are no new coloured particles in our model, the SM Higgs production cross sections for the various processes (gg fusion, vector boson fusion (VBF), VH associated production, \(t {\bar{t}} H\), etc.) are all rescaled by the mixing factor \(O_{11}^2\le 1\), implying
for any of those processes. Here, \(\mathrm{BR}(H\rightarrow \gamma \gamma )\) is computed considering the results obtained in the previous section for the diphoton Higgs decay in the MSBM. The subindex ‘SM’ obviously refers to the quantities within the SM. In our scan we use the naive weighted average (without including possible correlations) of the \(\mu _{\gamma \gamma }\) values corresponding to the different Higgs production processes, reported by the ATLAS collaboration in Ref. [20],
We use as constraint the agreement of the model prediction with this value within one standard deviation.
For the remaining SM decay modes of the Higgs boson, which we generically denote with the shorthand (with
, no sum over channels), the widths are suppressed by the factor \(O_{11}^2\). Then, we have
![](http://media.springernature.com/lw377/springer-static/image/art%3A10.1140%2Fepjc%2Fs10052-020-7952-4/MediaObjects/10052_2020_7952_Equ32_HTML.png)
where \(\Gamma _\text {SM}\) is the SM Higgs total width, \(\Gamma _\text {new}\) the width into new modes (e.g. lighter scalars) and the sum in the denominator comprises all SM decay channels other than \(\gamma \gamma \). Note that, as defined above, is the same for all the individual channels and production processes. Computing the naive weighted average of the
signal strengths for non-diphoton decays given in Ref. [20], one obtains
. Restricting ourselves to the
interval that has physical meaning within this model, the constraint translates into the \(1\sigma \) lower limit
![](http://media.springernature.com/lw230/springer-static/image/art%3A10.1140%2Fepjc%2Fs10052-020-7952-4/MediaObjects/10052_2020_7952_Equ33_HTML.png)
which we require in our scan. In addition, we include in our scan the limit
from direct searches for new scalars at LEP [17]. The limits on \(H_4\) from LEP are weaker than the indirect constraints from Higgs measurements. On the other hand, there are limits on its decay into \(\gamma \gamma \). Using \(\text {BR}(H_4 \rightarrow \gamma \gamma )_\text {SM} = 1.77 \times 10^{-3}\) [21], the limits in Ref. [14] translate into the constraint
Direct limits from the production of \(H_{3,4}\) and cascade decay \(H_{3,4} \rightarrow H_2 H_2 \rightarrow 4b\) were reviewed in Ref. [10] (see also Sect. 6). In Sect. 4 we mention, when relevant, the limits from processes involving the decays \(H_2 \rightarrow \gamma \gamma \).
Indirect limits from the S and T parameters [22, 23] do not lead to further constraints in the parameter space. The contribution of the new scalars is quite small because it is suppressed by the small mixings \(O_{1j}^2\), and the new scalars have masses that are not far away from \(M_H\). Using the expressions in Ref. [24], we find that the extra contribution to T is at the level of \(2 \times 10^{-3}\), and the extra contribution to S is of the order of \(8 \times 10^{-3}\), in good agreement with the latest determinations (assuming \(U=0\)) of \(T = 0.06 \pm 0.06\), \(S = 0.02 \pm 0.07\) [25]. The contributions of the lepton singlets to S and T vanish, and so does the contribution of the vector-like lepton doublet in the limit considered of no mixing, since in that case the two mass eigenstates are degenerate (see Ref. [25] for a review).
The purely scalar contributions to electric dipole moments vanish because the pseudoscalar component of the coupling to SM fermions vanishes, see Eq. (22). The only contributions arise from two-loop Barr-Zee [26] diagrams with closed loops of heavy leptons and exchange of a scalar \(H_j\) and a Z boson or photon. The amplitudes are proportional to
where we have simplified Eq. (18) for the case of vanishing heavy lepton mixing. Therefore, the individual amplitudes are already suppressed by small mixing factors, and in most of the parameter space fulfiling the rest of contstraints \(\tan \beta \le 1\). Moreover, the contributions of the two heavy leptons have opposite sign and there is a (partial) cancellation between them, which is exact when \(m_{E_1} = m_{E_2}\). Furthermore, there is a partial cancellation between the contributions of the several scalars, since \(\sum _{j=1}^4 O_{1j} O_{4j} = 0\) by unitarity. As a result, the contributions to electric dipole moments are below experimental bounds, as we have explicitly verified for the case of the electron.
In Fig. 2 we show in ternary plots the results obtained for the decay BRs of the four scalars, corresponding to \(6 \times 10^4\) allowed points in parameter space. As one can see from the top-left panel of this figure, the aforementioned constraints on the Higgs signal-strength parameters imply for the SM-like Higgs boson \(\mathrm{BR}(H\rightarrow H_i H_j) \lesssim 0.1\) with \(0.55 \lesssim \mathrm{BR}(H\rightarrow b {\bar{b}}) \lesssim 0.6\) (the SM value is approximately 0.58 [27]). As expected, the lightest scalar \(H_2\) predominantly decays into \(\gamma \gamma \) and \(b{\bar{b}}\) (top-right panel), being the probability to decay into other modes (\(\tau \tau \), \(c{{\bar{c}}}\) and gg, labeled as ‘rest’) around 20%, at most. Notice that \(\text {BR}(H_2\rightarrow \gamma \gamma )\) can even reach unity. For the remaining two scalars \(H_{3,4}\) (bottom panels) the trend is less clear, though they are expected to predominantly decay into \(H_2 H_2\) and \(\gamma \gamma \). Notice that \(\text {BR}(H_{3,4} \rightarrow \gamma \gamma ) = 1\) (apex of the triangle) is possible, as well as \(\text {BR}(H_{3,4} \rightarrow H_2 H_2) = 1\) (lower-right vertex).
The decay BRs of the \(Z'\) boson presented in Fig. 3 show little dependence on scalar mixing. In the left panel the decays into pairs of quarks, heavy leptons and scalars are compared. The relative BRs basically depend on the masses of the new leptons. When these are heavy, the corresponding \(Z'\) decay modes are kinematically forbidden (base of the triangle), while if \(m_F \ll M_{Z'}\) they dominate the \(Z'\) decays with \(\text {Br}(Z' \rightarrow F {{\bar{F}}}) = 36/53 \simeq 0.68\). In the right panel we show \(\text {BR}(Z' \rightarrow H_3 H_4)\) and \(\text {BR}(Z' \rightarrow H_2 H_3)\), which are the most interesting ones for the discussion in the next section. (The decays involving the SM Higgs H are very suppressed by mixing, and \(Z' \rightarrow H_2 H_4\) follows a similar pattern as \(Z' \rightarrow H_2 H_3\).) As it can be seen from the same plot, these decays can have a BR up to 50%. All these results are in agreement with those in Ref. [10] for model 2.
In summary, Figs. 2 and 3 reflect the fact that there is a variety of interesting signals coming from the cascade decays of the heavy resonance \(Z'\) into scalar pairs, allowed by present constraints. The scenario corresponding to \(Z' \rightarrow H_3 H_4\), with \(H_{3,4} \rightarrow H_2 H_2\) followed by \(H_2 \rightarrow b \bar{b}\) provides the perfect template for the stealth-boson hypothesis, as investigated in Refs. [9, 10]. On the other hand, final states with multi-pronged jets and/or multiphotons have been overlooked in the literature. In the following section we will discuss some interesting features of these signals in the context of the MSBM.
4 Multiphoton signals
When scalars \(H_i\) are produced either directly or from the decay of a heavy resonance, their decays may lead to multiphoton signals. We will review their main features in this section, with results obtained at the parton level. For the signal generation we use MadGraph5 [28] with the model implemented using Feynrules [29] and the universal Feynrules output [30]. We do not aim at a sensitivity estimate for these signals, as their detection relies on tools (e.g. the identification of collimated photons and photons within jets) that are not available at the level of fast simulation.
The possible cross section for these signals is determined in first place by the \(Z'\) production cross section which is
for the \(M_{Z'} = 2.2\) TeV benchmark point considered. For this mass, dijet constraints [31] set an upper limit \(\sigma (pp \rightarrow Z' \rightarrow jj) \times A \le 24.3\) fb at the 95% CL, where A is the acceptance for the event selection in the rapidity region \(|y| \le 0.6\). We evaluate the acceptance by simulating a \(Z' \rightarrow jj\) sample using Pythia [32] and Delphes [33], with jets reconstructed using FastJet [34]. The resulting limits on the coupling depend on the masses of the heavy leptons \(F=E,N\), and range between
Note however that the ATLAS limits are obtained under the assumption of a \(Z'\) width \(\Gamma _{Z'} \le 0.15 M_{Z'}\), which is not always the case. Therefore, the actual limits may be looser. For \(M_{H_i} \ll M_{Z'}\) (which is the case in our benchmark), and depending on \(m_F\), the \(Z' \rightarrow H_i H_j\) branching ratios are
Since the SM scalar doublet \(\Phi \) does not couple to the \(Z'\) boson, the decay modes involving the Higgs boson H turn out to be very suppressed by mixing. The decays into the remaining scalar pairs \(H_2 H_3\), \(H_2 H_4\) and \(H_3 H_4\) may be comparable, or the decay into one pair may dominate over the others, as seen in Fig. 3. If we assume for definiteness \(m_{F} > M_{Z'}/2\), we then have
for the signals we study in the remainder of this section, with a maximum of 121 fb.
Parton-level distributions for \(Z' \rightarrow H_{2,3} X, H_{2,3} \rightarrow \gamma \gamma \), being X some additional particle. Left: azimuthal angle difference between the two photons. Right: rapidity separation between the two photons. For comparison, the size of the ECAL clusters used for photon reconstruction in the CMS and ATLAS detectors is shown with vertical solid lines. The vertical dashed lines indicate the maximum extension in \(\Delta \phi \) in case of photon conversion
The first signal addressed, boosted diphotons, is produced from \(H_{2-4} \rightarrow \gamma \gamma \), when the scalars are produced from \(Z'\) decay. Instead, the second and third signals arise from the (expected) dominant decay \(H_{3,4} \rightarrow H_2 H_2\), when one or the two lighter scalars decay into photons. These signals are possible if the lightest scalar has a sizeable branching ratio into \(\gamma \gamma \), as it may be the case in our benchmark scenario (see Fig. 2). Note that in the remainder of this section we discuss the type of ‘object’ produced by one of the boosted scalars in \(Z' \rightarrow H_i H_j\); an actual event will have two of such objects, or one plus a hadronically-decaying scalar, etc.
4.1 Boosted diphotons from \(H_i \rightarrow \gamma \gamma \)
The scalars \(H_i\) can be directly produced, and their decay into two photons gives signals such at the CMS excess investigated in the next section. On the other hand, when \(H_i\) originate from the decay of a heavy resonance (in our case the \(Z'\) boson), they produce two photons that in the laboratory frame are relatively close in pseudo-rapidity \(\eta \) and azimuthal angle \(\phi \). These diphotons usually fail to pass the isolation criteria for prompt photons imposed in experimental analyses. To clarify this point, we briefly review how photons are identified.
Photons are reconstructed from energy depositions in the electromagnetic calorimeters (ECAL) of the ATLAS and CMS detectors. (The discrimination between photons and electrons is based on tracks in the inner detectors.) In the central region, the calorimeter cells have a size \(\Delta \eta \times \Delta \phi = 0.025 \times 0.0245\) for the ATLAS ECAL [35]Footnote 2 and \(0.0174 \times 0.0174\) for the CMS ECAL [36]. Photons are reconstructed from \(3 \times 5\) clusters in ATLAS, yielding a cluster size of \(\Delta \eta \times \Delta \phi = 0.075 \times 0.124\). In CMS \(5 \times 5\) clusters are used, with a size \(\Delta \eta \times \Delta \phi = 0.087 \times 0.087\). In both experiments, additional rows up to \(\Delta \phi = 0.172\) (ATLAS) and \(\Delta \phi = 0.609\) (CMS) may be added in the case of photon conversion into an \(e^+ e^-\) pair, which are curved in opposite directions in \(\phi \) due to the magnetic field. For the photons so identified, there are several ‘reference’ sets of photon isolation conditions, from tighter to looser, which require, among other, the absence of additional radiation above some given thresholds in a cone of \(\Delta R \equiv \sqrt{(\Delta \eta )^2 + (\Delta \phi )^2} \le 0.2-0.3\) [35, 36].
Parton-level distributions for \(Z' \rightarrow H_3 X, H_3 \rightarrow H_2 H_2 \rightarrow 4\gamma \), being X some additional particle. Left: minimum and maximum azimuthal angle difference among the photons. Right: the same, for the rapidity difference. For comparison, the size of the ECAL clusters used for photon reconstruction in the CMS and ATLAS detectors is shown. The dashed lines indicate the maximum extension in \(\Delta \phi \) in case of photon conversion
Let us now consider the decay \(Z' \rightarrow H_i X\) with \(M_{Z'} = 2.2\) TeV, being X an additional particle (e.g. another scalar of our model with an unspecified decay mode). In Fig. 4 we show the normalised distributions of the diphoton separation in \(\Delta \phi \) and \(\Delta \eta \) when \(H_i\) is either \(H_2\) or \(H_3\) with \(M_{H_3} = 96 \, \text {GeV}\) and \(M_{H_2} = 30 \, \text {GeV}\), together with the sizes of the ATLAS and CMS clusters. For \(M_{H_2} = 30 \, \text {GeV}\), the two photons lie within the same cluster most of the time. Thus, one expects that an event selection for prompt photons may be sensitive to this decay. When \(M_{H_3} = 96 \, \text {GeV}\), the two photons are more separate and seldom hit the ECAL within the same cluster. The probability for that to happen can be estimated from the distributions shown in Fig. 4, being 0.046 for CMS and 0.057 for ATLAS. But the efficiency for the detection of a diphoton as an isolated photon, taking into account the possible extension in \(\phi \) of the clusters in case of photon conversion, cannot be simply estimated at the level of fast simulation.
An ATLAS search further attempts to discriminate between a single photon and more than one photon by using the fine granularity of the first layer of the ECAL [37]. Although this search looks for signatures with a pair of highly-collinear photons, the event selection is only sensitive to \(M_{H_i} / M_{Z'} \le 0.01\), not covering the 30 GeV benchmark mass considered. For example, the limit obtained for \(M_{Z'} = 2\) TeV, \(M_{H_i,H_j} \lesssim 10\) GeV is \(\sigma (pp \rightarrow Z') \times \text {BR}(Z' \rightarrow H_i H_j \rightarrow 4\gamma ) \lesssim 0.15\) fb. Therefore, it is safe to say that a dedicated search into pairs of collimated diphotons (or a diphoton plus a jet) using the information from adjacent clusters and applying ad-hoc isolation criteria will be sensitive to the cross sections expected from Eq. (40).
4.2 Multiphotons from \(H_3 \rightarrow H_2 H_2 \rightarrow 4\gamma \)
The direct production \(gg \rightarrow H_3 \rightarrow H_2 H_2 \rightarrow 4 \gamma \) (or, equivalently, \(H_4 \rightarrow H_2 H_2 \rightarrow 4 \gamma \)) produces a signal with four well-separated photons. An ATLAS search for events with three or more photons [38] in Run 1 at 8 TeV did not observe any significant excess over the SM background expectation, setting an upper limit of 171 fb on the cross section for \(gg \rightarrow H \rightarrow aa \rightarrow 4 \gamma \), with \(m_a = 10\) GeV. We can conservatively use the same limit for \(gg \rightarrow H_3 \rightarrow H_2 H_2\), bearing in mind that the efficiency will be slightly smaller for \(M_{H_3} < M_H\). Using a cross section \(\sigma (gg \rightarrow h) = 26.5\) pb for a 95 GeV SM-like scalar [27], and assuming a 100% branching ratio of \(H_3 \rightarrow H_2 H_2\), we can obtain the upper bound
With our scan over parameter space we have verified that \(\text {BR}(H_2 \rightarrow \gamma \gamma )\) can be as large as unity fulfilling the above limit.
The cascade decay of a boosted scalar \(H_3 \rightarrow H_2 H_2 \rightarrow 4 \gamma \) yields collinear multiphotons that are much alike the collinear diphotons examined above. We show in Fig. 5 the kinematical distributions for the minimum and maximum separation in \(\eta \) and \(\phi \) among the four photons, taking \(M_{Z'} = 2.2\) TeV, \(M_{H_3} = 96\) GeV, \(M_{H_2} = 30\) GeV. It follows from the previous discussion that some of the photons may be caught within the same cluster, but there are often additional photons that produce energy depositions outside it, thereby spoiling the standard isolation requirements. Such signal also has little or no coverage at all by current searches. We point out that, as in the previous case of boosted diphotons, modified isolation criteria may be applied in dedicated searches for collimated multiphotons, that is, with the radiation contributions to the isolation cone from other photons subtracted before applying the isolation requirements.
4.3 Photons within jets: \(H_i \rightarrow H_j H_j \rightarrow \gamma \gamma b {{\bar{b}}}\)
In analogy with the previous case, \(gg \rightarrow H_3 \rightarrow H_2 H_2 \rightarrow \gamma \gamma b {{\bar{b}}}\) produces a signal with two photons and two b quarks that are generally isolated. Searches for heavy scalars h decaying into Higgs boson pairs, \(gg \rightarrow h \rightarrow H H \rightarrow \gamma \gamma b {\bar{b}}\), are carried out at the LHC [39, 40] but they obviously do not cover the range of interest \(M_{H_3} \le M_H\). On the other hand, searches for exotic decays of the Higgs boson into a pair of lighter particles \(H \rightarrow aa\) are investigated in a variety of final states, but not with a photon pair. Besides, those analyses are restricted to a narrow mass window around \(M_{H}\).
The decay of a boosted \(H_3 \rightarrow H_2 H_2 \rightarrow \gamma \gamma b {\bar{b}}\) produces a conspicuous jet containing hard photons. Such jet is generally narrow, at it can be inferred from the distribution of the \(\Delta R\) distance between the two b quarks in Fig. 6 (left). The two photons are at a medium distance from the b quarks, as shown in the right panel. A quite distinctive feature of these jets is the large energy fraction carried by the photons. We can define it as
where the subindices 1 and 2 label the different photons and b quarks, and the energies are considered in the laboratory frame. The kinematical distribution of this quantity, depicted in Fig. 7, shows that a large fraction of the jet energy may be carried by the photons. As antecipated above, this is a quite peculiar signature. Searches for such objects seem difficult but feasible, indeed the CMS Collaboration has searched for events with jets that contain one photon and two gluons [41].
5 The 96 GeV case
The diphoton signals from direct production of new scalars bear a special interest in view of some CMS measurements which we now summarise. The CMS Collaboration has found [14] two nearby excesses when searching for new scalars h decaying into photon pairs in Run 1 and Run 2 data,
-
Run 1: \(2\sigma \) local significance at a diphoton invariant mass \(m_{\gamma \gamma } = 97.6\) GeV, with 19.7 \({\hbox {fb}}^{-1}\) at 8 TeV;
-
Run 2: \(3\sigma \) local significance at a slightly lower mass \(m_{\gamma \gamma } = 95.3\) GeV, with 35.9 \({\hbox {fb}}^{-1}\) at 13 TeV.
They also provide limits on two groups of processes mediating the production of h: (a) \(gg \rightarrow h\) plus \(t {\bar{t}}h\) production; (b) VBF plus Vh production. The excess is present in both, though more prominent in the former group. The reported combination of Run 1 and Run 2 data has a maximum \(2.8\sigma \) local significance at \(m_{\gamma \gamma } = 95.3\) GeV. This significance is slightly smaller than that of Run 2 data alone because the two excesses are not located precisely at the same mass, and the mass shift is larger than the detector resolution. However, one cannot discard the possibility that the slight mass difference arises from some statistical or systematic effect: one can also notice a shift in the same direction between the peaks at \(M_Z\), coming from \(Z \rightarrow e^+ e^-\), with the electrons misidentified as photons. In Run 1, the peak is located at 91.0 GeV, while in Run 2 it is at 90.0 GeV. In this regard, it is also worth remembering that the final combined Run 1 measurements of the Higgs boson mass by the ATLAS Collaboration in the cleanest decay modes \(H \rightarrow \gamma \gamma \) and \(H \rightarrow ZZ^*\) yielded masses \(M_H = 125.98 \pm 0.50\) GeV and \(M_H = 124.51 \pm 0.51\) GeV, respectively, with uncertainties dominated by statistics [42]. These two measurements were compatible only at the \(2\sigma \) level even if the collected samples were large enough to have \(5.2\sigma \) and \(8.1\sigma \) significances for the Higgs signal, respectively, in the \(\gamma \gamma \) and \(ZZ^*\) channel.
The CMS excess corresponds to a signal strength \(\mu _\text {CMS} = 0.6 \pm 0.2\) [49], defined as
with the subindex ‘SM’ indicating that the quantity is referred to a SM-like scalar h, of course evaluated for its corresponding mass \(M_h\). This excess can be accommodated with an eigenstate that is mainly a singlet, by either enhancing the production via gluon fusion or the decay into \(\gamma \gamma \). As the CMS excess is also seen in the VBF and Vh category, the latter possibility seems quite reasonable. Since in our model the new particles do not contribute to the \(gg \rightarrow H_i\) amplitudes, then the cross-section ratio in the r.h.s. of Eq. (43) reduces to the mixing factor \(O_{13}^2\). (As mentioned, we consider \(H_3\) as being the new scalar candidate to explain the excess.) Taking \(\text {BR}(H_3 \rightarrow \gamma \gamma )_\text {SM} = 1.43 \times 10^{-3}\) [27], we find
From Fig. 2 it is clear that the agreement with (44) is possible even for small \(O_{13}\) mixing, and in the presence of non-SM decay modes. Also, we see that it is quite possible that one scalar gives an observable signal \(pp \rightarrow H_3 \rightarrow \gamma \gamma \) while the others do not, if the mixings \(O_{12}\) and \(O_{14}\) are smaller.
The CMS excess has sparked some interest in the framework of SM extensions with a scalar singlet [43,44,45,46,47,48,49,50,51,52,53], also in connection with a \(2.3\sigma \) local excess observed at a mass of 98 GeV in LEP searches for the SM Higgs boson (\(e^+ e^- \rightarrow Zh\), \(h \rightarrow b {\bar{b}} \)) [17]. The fact that the two excesses are located at nearly the same mass, suggests that they could originate from \(\gamma \gamma \) and \(b {\bar{b}}\) decays of the same particle. The signal strength of the LEP excess, defined as
has been computed in Ref. [54], yielding \(\mu _\text {LEP} = 0.117 \pm 0.057\). The subindex ‘SM’ refers to quantities corresponding to a SM-like scalar with mass \(M_h\). In our case, the cross section ratio in the r.h.s. of Eq. (45) is again \(O_{13}^2\). Taking \(\text {BR}(h \rightarrow b {\bar{b}})_\text {SM} = 0.8\) [27], Eq. (45) implies
As shown in Fig. 8, both the CMS and LEP excesses can be accommodated in the MSBM. The value of \(\mu _\text {CMS}\) can be either larger or smaller than the best-fit value, in agreement with the results in Fig. 2. On the other hand, \(\mu _\text {LEP}\) cannot be much larger than the best-fit value, as it is limited by the Higgs constraints on \(O_{13}^2\) — see Eq. (46). Note that in our benchmark scenario, the decay \(H_3 \rightarrow H_2 H_2\) is allowed; for heavier \(H_2\) this channel would be closed and the lower left side of the plot with both \(\mu _\text {CMS}\) and \(\mu _\text {LEP}\) small would be less populated.
The scenarios we are interested in are those in which \(H_3 \rightarrow H_2 H_2\) (as well as \(H_4 \rightarrow H_2 H_2\)) have sizeable branching ratios, thus producing the multiphoton signals described in the previous section. Therefore, we drop the interpretation of the LEP excess, which requires large branching ratio \(H_3 \rightarrow b {\bar{b}}\) (see the appendix). When we require \(\mu _\text {CMS} = 0.6 \pm 0.2\), the phenomenology remains mostly the same even if the parameter space is reduced. This is shown in Fig. 9. On the top-left panel we observe that the pattern of \(Z'\) decays into scalars remains practically unaltered. The same happens for the lightest scalar \(H_2\) (top right). For \(H_3\) (bottom left) there is a stronger preference for \(H_3 \rightarrow H_2 H_2\) decays, and for \(H_4\) (bottom right) the behaviour is similar. The gray points in these plots correspond to the points in Figs. 3 and 2, while the black points are those which are also compatible with the CMS excess, with \(\mu _\text {CMS} \in [0.4,0.8]\). Therefore, one can conclude that the multiphoton signals studied in Sect. 4 are compatible with the CMS excess too.
6 Summary and discussion
New scalars are thoroughly searched for at the LHC. Still, there is a number of theoretically and phenomenologically well-motivated signatures that are yet to be explored, even for scalars lighter than the Higgs boson. The measurement of the couplings of the SM Higgs [20] suggests that new scalars, if they exist, have small couplings to the weak bosons.Footnote 3 Therefore, provided these new scalars have small Yukawa coupling, their most visible signals may arise when produced from the decay of a new heavy resonance. (In case of very small mixing, these could even be the signals with larger cross section.) On the other hand, the cascade decay into lighter scalars may be the dominant decay mode, if kinematically allowed. In this context, we have investigated a number of collider signals involving diphotons from heavy resonance decays into scalars.
We have worked within the framework of the MSBM [10], an anomaly-free leptophobic \(\text {U}(1)_{Y'}\) extension of the SM, in which anomaly cancellation is ensured by a set of vector-like leptons. Additionally to the SM Higgs doublet, two complex singlets charged under \(\text {U}(1)_{Y'}\) are introduced with the purpose of providing mass to the \(Z'\) boson and the new leptons, and of allowing for \(Z'\) decays into scalars. Overall, the model contains four physical scalar particles \(H_i\), being \(H_1 \equiv H\) the 125 GeV Higgs boson. To provide an interpretation for a CMS diphoton excess, we choose the mass of \(H_3\) to be 96 GeV. Moreover, a lighter scalar (\(H_2\)) with mass in the \(20-40\) GeV range has also been considered to allow for the decay \(H_3 \rightarrow H_2 H_2\). Since the mass of \(H_4\) does not play a significant role in the signals involving \(H_3\) and \(H_2\), we take \(M_{H_4}=105\) GeV (relatively close to \(M_{H_3}\)). Its signals are analogous to the ones with \(H_3\).
We have performed a scan over the MSBM parameter space in order to explore the allowed range of the mixing and the decay modes of the new scalars, compatible with existing constraints from several SM Higgs production and decay processes [20]. We have found that the allowed \(\text {BR}(H_{2-4} \rightarrow \gamma \gamma )\) span over the whole range from zero to unity, while keeping \(\text {BR}(H \rightarrow \gamma \gamma )\) for the SM Higgs boson in agreement with current measurements. This is possible due to the fact that the scalar singlets couple to the new leptons (they actually provide their mass), while the SM Higgs doublet does not. Thus, the \(H_{2-4} \rightarrow \gamma \gamma \) loop amplitudes can be enhanced, while those for \(H \rightarrow \gamma \gamma \) and, of course, \(gg \rightarrow H\) are not modified.
Direct production of \(H_{2-4}\) is suppressed by the small mixing with the SM Higgs doublet. Still, we have found that a CMS excess at 96 GeV may be easily accommodated by one of the \(H_i\) (which we have labelled as \(H_3\)) since, as mentioned, the decay width to \(\gamma \gamma \) can be enhanced. If the dominant decay mode is \(H_3 \rightarrow b {\bar{b}}\), then the same scalar can also explain a LEP excess found around the same mass. However, the most intriguing consequences arise when the main decay channel of this particle is \(H_3 \rightarrow H_2 H_2\), which is otherwise the most natural situation if there is a lighter scalar \(H_2\) and available phase space. From Fig. 9, we observe that in most of the parameter space where \(H_3 \rightarrow \gamma \gamma \) fits the CMS excess:
-
(i)
\(H_3\) predominantly decays into \(H_2 H_2\), with \(\text {BR}(H_3 \rightarrow \gamma \gamma )\) at the percent level. (Note however that there are points where \(H_3 \rightarrow \gamma \gamma \) dominates.)
-
(ii)
\(H_2 \rightarrow b {\bar{b}}\) is expected to dominate.
-
(iii)
The main decay modes of \(H_4\) are \(H_2 H_2\) and \(\gamma \gamma \).
We can then have \(pp \rightarrow Z' \rightarrow H_3 H_4\), with \(H_3 \rightarrow \gamma \gamma \) and \(H_4 \rightarrow H_2 H_2 \rightarrow 4b\). (If, instead, we consider the decay \(H_3 \rightarrow H_2 H_2 \rightarrow 4b\), the signal contains two four-pronged jets.) This signal yields a pair of collinear photons from \(H_3\), plus a multi-pronged jet from \(H_4\) which is quite elusive and requires specific tools to be efficiently separated from the QCD background (see for example Refs. [55,56,57]). So does the boosted diphoton which, as seen in Sect. 4.1, will usually fail the isolation criteria for prompt photons and requires a specific selection. The decay \(pp \rightarrow Z' \rightarrow H_3 H_2\) produces a collinear photon pair plus a two-pronged jet of very small mass, which is hard to separate from the QCD background.
It is worth recalling that the CMS Collaboration has reported another excess when searching for heavy resonances decaying into \(Z\gamma \), with 35.9 \({\hbox {fb}}^{-1}\) of data at 13 TeV [58]. The excess reaches a local significance of \(3.6\sigma \) at a mass of 2 TeV in the hadronic channel, without any excess in the leptonic channel. In a search with similar sensitivity [59], the ATLAS Collaboration did not find any excess at that mass. The aforementioned process \(pp \rightarrow Z' \rightarrow H_{3} H_4\), with \(H_{3} \rightarrow \gamma \gamma \) and \(H_4 \rightarrow H_2 H_2 \rightarrow 4b\), is an obvious candidate to explain an excess in the hadronic channel without a counterpart in the leptonic channel. (Note that the pruned jet mass distribution of the \(H_4\) jet with \(M_{H_4} = 105\) GeV falls within the window selected in the analysis, as the pruning significantly decreases the mass of multi-pronged jets [60].) However, this interpretation faces two difficulties. First and foremost, the fact that the CMS Collaboration does not see any excess in an event category with b-tagged jets. A multi-pronged jet with four b quarks should in principle yield an excess in this category, although the b tagging efficiency for such a jet can only be verified with a full simulation of the CMS detector. Second, the fact that the excess is not seen in the ATLAS analysis does not have an obvious explanation: as we have mentioned in Sect. 4.1, the efficiency for a boosted diphoton to be identified as a photon depends on details of the photon reconstruction that are not available without the full detector simulation. Besides, with a fast simulation we have estimated that the efficiency for the fat jet from \(H_4 \rightarrow 4b\) to pass the jet substructure selection in the ATLAS and CMS analyses is similar.
Dropping the connection with the CMS 96 GeV excess, the cascade decay of the heavier scalars with diphoton decay of the lighter one may give two interesting signals:
-
(i)
\(gg \rightarrow H_3 \rightarrow H_2 H_2 \rightarrow 4 \gamma \). This signal is covered by a search for events with at least three photons by the ATLAS Collaboration using 8 TeV data [38]. To the best of our knowledge, there are no searches at 13 TeV in this final state.
-
(ii)
\(gg \rightarrow H_3 \rightarrow H_2 H_2 \rightarrow \gamma \gamma b {\bar{b}}\). There are searches for heavy scalars decaying into the Higgs boson in this final state [39, 40] but for lighter \(H_3\) this is not covered.
Note that the same signals can appear with the scalar \(H_4\) replacing \(H_3\), but for brevity we only refer to the latter. The fully hadronic signals \(H_3 \rightarrow H_2 H_2 \rightarrow 4b\) have quite large backgrounds and small sensitivity. In WH/ZH associated production, a search for \(H \rightarrow aa \rightarrow b {\bar{b}} b {\bar{b}} \) (where a is a lighter particle, e.g. an axion) by the ATLAS Collaboration [61] has a sensitivity one order of magnitude smaller than would be required, and does not cover \(M_{H_3} = 96\) GeV. At the Tevatron, the CDF Collaboration performed a search for pair production of new particles Y, each decaying into two jets, \(p {\bar{p}} \rightarrow YY \rightarrow jjjj\) [62]. The mass range explored \(M_{Y} \ge 50\) GeV does not cover \(M_{H_2} = 30\) GeV and, in addition, the sensitivity is two orders of magnitude below the possible size of the signals.
The same cascade decays, when \(H_3\) is produced in the decay of the \(Z'\) boson, yield boosted and collinear final state particles that form complex objects:
-
(i)
\(Z' \rightarrow H_3 X \rightarrow 4 \gamma + X\), with X an additional particle. The four photons are produced closely (see Fig. 5) and require a special object selection.
-
(ii)
\(Z' \rightarrow H_3 X \rightarrow \gamma \gamma b {\bar{b}}+ X\). The two b quarks are produced close in \(\Delta R\) and merge into a single jet; the same jet contains two very energetic photons. The identification of such an object requires a specific tagging.
This type of final-state objects (collinear multiphotons and jets with photons) are not currently searched for. Notice that the ‘photon jets’ searched for by the ATLAS Collaboration [37] are much more collinear, and that search does not cover the multiphotons produced at intermediate distances, such as those shown in Fig. 5.
To conclude, cascade decays including multiphotons in the final state yield intriguing signals that are not covered by current experimental searches. The CMS 96 GeV excess, whether corresponds to a new particle or not, provides an extra motivation to explore these final states. In direct connection with this excess, a dedicated search for heavy resonances decaying into collimated diphotons plus a fat (multi-pronged) jet would be of great interest. Especially, because it might clarify the situation regarding the 2 TeV excess in the search for heavy resonances decaying into \(Z\gamma \).
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: This manuscript has no associated data.]
Notes
This model is denoted by ‘model 2’ in ref. [10]; ‘model 1’ is similar, differing only from the fact that vector-like quarks are introduced instead of leptons. The heavy lepton sector has also been considered in refs. [11,12,13] in the context of dark-matter phenomenology. Indeed, if \(N_1\) is the lightest among the new leptons, then it is a natural dark-matter candidate.
The ATLAS ECAL has a first layer with a higher granularity, this size refers to the second layer where most of the photon energy is deposited.
This is for example the case of two Higgs doublet models, where in the alignment limit the new neutral scalars \(H^0\), \(A^0\) do not couple to WW and ZZ, and of models with additional singlets, where the coupling of the mass eigenstate is acquired via a small doublet component.
References
F. Englert, R. Brout, Broken symmetry and the mass of gauge vector mesons. Phys. Rev. Lett. 13, 321 (1964)
P.W. Higgs, Broken symmetries, massless particles and gauge fields. Phys. Lett. 12, 132 (1964)
P.W. Higgs, Broken symmetries and the masses of gauge bosons. Phys. Rev. Lett. 13, 508 (1964)
G. Aad et al., [ATLAS Collaboration], Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC. Phys. Lett. B 716, 1 (2012). arXiv:1207.7214 [hep-ex]
S. Chatrchyan et al., [CMS Collaboration], Observation of a New Boson at a Mass of 125 GeV with the CMS Experiment at the LHC. Phys. Lett. B 716, 30 (2012). arXiv:1207.7235 [hep-ex]
J.A. Aguilar-Saavedra, F.R. Joaquim, Multiboson production in \(W^{\prime }\) decays. JHEP 1601, 183 (2016). arXiv:1512.00396 [hep-ph]
R.N. Mohapatra, J.C. Pati, A natural left–right symmetry. Phys. Rev. D 11, 2558 (1975)
G. Senjanovic, R.N. Mohapatra, Exact left–right symmetry and spontaneous violation of parity. Phys. Rev. D 12, 1502 (1975)
J.A. Aguilar-Saavedra, Stealth multiboson signals. Eur. Phys. J. C 77(10), 703 (2017). [arXiv:1705.07885 [hep-ph]]
J.A. Aguilar-Saavedra, F.R. Joaquim, The minimal stealth boson: models and benchmarks. JHEP 1910, 237 (2019). arXiv:1905.12651 [hep-ph]
M. Duerr, P. Fileviez Perez, Theory for Baryon number and dark matter at the LHC. Phys. Rev. D 91(9), 095001 (2015). arXiv:1409.8165 [hep-ph]
S. Caron, J.A. Casas, J. Quilis, R. Ruiz de Austri, Anomaly-free dark matter with harmless direct detection constraints. JHEP 1812, 126 (2018). arXiv:1807.07921 [hep-ph]
J.A. Aguilar-Saavedra, J.A. Casas, J. Quilis, R. Ruiz de Austri, Multilepton dark matter signals. JHEP 04, 069 (2020)
A.M. Sirunyan et al., [CMS Collaboration], Search for a standard model-like Higgs boson in the mass range between 70 and 110 GeV in the diphoton final state in proton-proton collisions at \(\sqrt{s}=\) 8 and 13 TeV. Phys. Lett. B 793, 320 (2019). arXiv:1811.08459 [hep-ex]
ATLAS Collaboration, Search for resonances in the 65 to 110 GeV diphoton invariant mass range using \(80~{\text{fb}}^{-1}\) of \(pp\) collisions collected at \(\sqrt{s}=13 ~\text{ TeV }\) with the ATLAS detector, ATLAS-CONF-2018-025
S. Heinemeyer, T. Stefaniak, A Higgs Boson at 96 GeV?!. PoS CHARGED 2018, 016 (2019). arXiv:1812.05864 [hep-ph]
R. Barate et al., [ALEPH and DELPHI and L3 and OPAL Collaborations and LEP Working Group for Higgs boson searches], Search for the standard model Higgs boson at LEP. Phys. Lett. B 565, 61 (2003). arXiv:hep-ex/0306033
V. Ilisie, A. Pich, Low-mass fermiophobic charged Higgs phenomenology in two-Higgs-doublet models. JHEP 1409, 089 (2014). arXiv:1405.6639 [hep-ph]
J.F. Gunion, H.E. Haber, G.L. Kane, S. Dawson, The Higgs Hunter’s Guide. Front. Phys. 80, 1 (2000)
G. Aad et al., [ATLAS Collaboration], Combined measurements of Higgs boson production and decay using up to \(80 ~ {\text{ fb }}^{-1}\) of proton-proton collision data at \(\sqrt{s}= 13~\text{ TeV }\) collected with the ATLAS experiment. Phys. Rev. D 101(1), 012002 (2020). arXiv:1909.02845 [hep-ex]
S. Heinemeyer et al., [LHC Higgs Cross Section Working Group], Handbook of LHC Higgs Cross Sections: 3. Higgs Properties. arXiv:1307.1347 [hep-ph]
M.E. Peskin, T. Takeuchi, A new constraint on a strongly interacting Higgs sector. Phys. Rev. Lett. 65, 964 (1990)
M.E. Peskin, T. Takeuchi, Estimation of oblique electroweak corrections. Phys. Rev. D 46, 381 (1992)
J.A. Aguilar-Saavedra, J. Alonso-González, L. Merlo, J.M. No, Exotic vectorlike quark phenomenology in the minimal linear model. Phys. Rev. D 101(3), 035015 (2020). arXiv:1911.10202 [hep-ph]
M. Tanabashi et al., [Particle Data Group], Review of Particle Physics. Phys. Rev. D 98(3), 030001 (2018)
S.M. Barr, A. Zee, Electric Dipole Moment of the Electron and of the Neutron. Phys. Rev. Lett. 65, 21 (1990). Erratum: [Phys. Rev. Lett. 65, 2920 (1990)]
D. de Florian et al. [LHC Higgs Cross Section Working Group], Handbook of LHC Higgs Cross Sections: 4. Deciphering the nature of the higgs sector. arXiv:1610.07922 [hep-ph]. See https://twiki.cern.ch/twiki/bin/view/LHCPhysics/CrossSections
J. Alwall et al., The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations. JHEP 1407, 079 (2014). arXiv:1405.0301 [hep-ph]
A. Alloul, N.D. Christensen, C. Degrande, C. Duhr, B. Fuks, FeynRules 2.0—a complete toolbox for tree-level phenomenology. Comput. Phys. Commun. 185, 2250 (2014). arXiv:1310.1921 [hep-ph]
C. Degrande, C. Duhr, B. Fuks, D. Grellscheid, O. Mattelaer, T. Reiter, UFO - the universal FeynRules output. Comput. Phys. Commun. 183, 1201 (2012). arXiv:1108.2040 [hep-ph]
G. Aad et al. [ATLAS Collaboration], Search for new resonances in mass distributions of jet pairs using \(139~{\text{ fb }}^{-1}\) of \(pp\) collisions at \(\sqrt{s}=13~ \text{ TeV }\) with the ATLAS detector. arXiv:1910.08447 [hep-ex]
T. Sjostrand, S. Mrenna, P.Z. Skands, A brief introduction to PYTHIA 8.1. Comput. Phys. Commun. 178, 852 (2008). arXiv:0710.3820 [hep-ph]
J. de Favereau et al., [DELPHES 3 Collaboration], DELPHES 3, A modular framework for fast simulation of a generic collider experiment. JHEP 1402, 057 (2014). arXiv:1307.6346 [hep-ex]
M. Cacciari, G.P. Salam, G. Soyez, FastJet user manual. Eur. Phys. J. C 72, 1896 (2012). arXiv:1111.6097 [hep-ph]
M. Aaboud et al. [ATLAS Collaboration], Measurement of the photon identification efficiencies with the ATLAS detector using LHC Run-1 data. Eur. Phys. J. C 76(12), 666 (2016). arXiv:1606.01813 [hep-ex]
V. Khachatryan et al. [CMS Collaboration], Performance of photon reconstruction and identification with the CMS detector in proton–proton collisions at \(\sqrt{s} = 8 ~\text{ TeV }\). JINST 10(08), P08010 (2015). arXiv:1502.02702 [physics.ins-det]
M. Aaboud et al. [ATLAS Collaboration], A search for pairs of highly collimated photon-jets in \(pp\) collisions at \(\sqrt{s} = 13 ~\text{ TeV }\) with the ATLAS detector, Phys. Rev. D 99(1), 012008 (2019). arXiv:1808.10515 [hep-ex]
G. Aad et al. [ATLAS Collaboration], Search for new phenomena in events with at least three photons collected in \(pp\) collisions at \(\sqrt{s} = 8 ~\text{ TeV }\) with the ATLAS detector. Eur. Phys. J. C 76(4), 210 (2016). arXiv:1509.05051 [hep-ex]
A.M. Sirunyan et al., [CMS Collaboration], Search for Higgs boson pair production in the \(\gamma \gamma \text{ b } \overline{\text{ b }}\) final state in pp collisions at \(\sqrt{s}=\) 13 TeV. Phys. Lett. B 788, 7 (2019). arXiv:1806.00408 [hep-ex]
M. Aaboud et al., [ATLAS Collaboration], Search for Higgs boson pair production in the \(\gamma \gamma b\bar{b}\) final state with 13 TeV \(pp\) collision data collected by the ATLAS experiment. JHEP 1811, 040 (2018). arXiv:1807.04873 [hep-ex]
A. M. Sirunyan et al. [CMS Collaboration], Search for physics beyond the standard model in events with overlapping photons and jets. Phys. Rev. Lett. 123(24), 241801 (2019)
G. Aad et al. [ATLAS Collaboration], Measurement of the Higgs boson mass from the \(H\rightarrow \gamma \gamma \) and \(H \rightarrow ZZ^{*} \rightarrow 4\ell \) channels with the ATLAS detector using \(25~ {\text{ fb }}^{-1}\) of \(pp\) collision data. Phys. Rev. D 90(5), 052004 (2014). arXiv:1406.3827 [hep-ex]
P.J. Fox, N. Weiner, Light signals from a lighter higgs. JHEP 1808, 025 (2018). arXiv:1710.07649 [hep-ph]
T. Biekötter, S. Heinemeyer, C. Muñoz, Precise prediction for the Higgs-boson masses in the \(\mu \nu \)SSM. Eur. Phys. J. C 78(6), 504 (2018). arXiv:1712.07475 [hep-ph]
D. Liu, J. Liu, C.E.M. Wagner, X.P. Wang, A light Higgs at the LHC and the B-Anomalies. JHEP 1806, 150 (2018). arXiv:1805.01476 [hep-ph]
F. Domingo, S. Heinemeyer, S. Paßehr, G. Weiglein, Decays of the neutral Higgs bosons into SM fermions and gauge bosons in the \({\cal{C}}{\cal{P}}\)-violating NMSSM. Eur. Phys. J. C 78(11), 942 (2018). arXiv:1807.06322 [hep-ph]
W.G. Hollik, S. Liebler, G. Moortgat-Pick, S. Paßehr, G. Weiglein, Phenomenology of the inflation-inspired NMSSM at the electroweak scale. Eur. Phys. J. C 79(1), 75 (2019). arXiv:1809.07371 [hep-ph]
L. Liu, H. Qiao, K. Wang, J. Zhu, A light scalar in the minimal dilaton model in light of LHC constraints. Chin. Phys. C 43(2), 023104 (2019). arXiv:1812.00107 [hep-ph]
T. Biekötter, M. Chakraborti, S. Heinemeyer, A 96 GeV Higgs Boson in the N2HDM. Eur. Phys. J. C 80(1)2, (2020)
J.M. Cline, T. Toma, Pseudo–Goldstone dark matter confronts cosmic ray and collider anomalies. Phys. Rev. D 100(3), 035023 (2019). arXiv:1906.02175 [hep-ph]
K. Choi, S.H. Im, K.S. Jeong, C.B. Park, A 96 GeV Higgs boson in the general NMSSM. Eur. Phys. J. C 79(11), 956 (2019)
T. Biekötter, S. Heinemeyer, C. Muñoz, Precise prediction for the Higgs-Boson masses in the \(\mu \nu \)SSM with three right-handed neutrino superfields. Eur. Phys. J. C 79(8), 667 (2019). arXiv:1906.06173 [hep-ph]
J. Cao, X. Jia, Y. Yue, H. Zhou, P. Zhu, The 96 GeV diphoton excess in the seesaw extensions of the natural NMSSM. Phys. Rev. D 101(5), 055008 (2020)
J. Cao, X. Guo, Y. He, P. Wu, Y. Zhang, Diphoton signal of the light Higgs boson in natural NMSSM. Phys. Rev. D 95(11), 116001 (2017). arXiv:1612.08522 [hep-ph]
J.A. Aguilar-Saavedra, J.H. Collins, R.K. Mishra, A generic anti-QCD jet tagger. JHEP 1711, 163 (2017). arXiv:1709.01087 [hep-ph]
T. Heimel, G. Kasieczka, T. Plehn, J.M. Thompson, QCD or what? Sci. Post Phys. 6(3), 030 (2019). arXiv:1808.08979 [hep-ph]
M. Farina, Y. Nakai, D. Shih, Searching for new physics with deep autoencoders. Phys. Rev. D 101(7), 075021 (2020)
A.M. Sirunyan et al., [CMS Collaboration], Search for Z\(\gamma \) resonances using leptonic and hadronic final states in proton–proton collisions at \(\sqrt{s}= 13~ \text{ TeV }\). JHEP 1809, 148 (2018). arXiv:1712.03143 [hep-ex]
M. Aaboud et al. [ATLAS Collaboration], Search for heavy resonances decaying to a photon and a hadronically decaying \(Z/W/H\) boson in \(pp\) collisions at \(\sqrt{s}=13 ~ \text{ TeV }\) with the ATLAS detector. Phys. Rev. D 98(3), 032015 (2018). arXiv:1805.01908 [hep-ex]
J.A. Aguilar-Saavedra, Running bumps from stealth bosons. Eur. Phys. J. C 78(3), 206 (2018). arXiv:1801.08129 [hep-ph]
M. Aaboud et al., [ATLAS Collaboration], Search for the Higgs boson produced in association with a vector boson and decaying into two spin-zero particles in the \(H \rightarrow aa \rightarrow 4b\) channel in \(pp\) collisions at \(\sqrt{s} = 13 ~ \text{ TeV }\) with the ATLAS detector. JHEP 1810, 031 (2018). arXiv:1806.07355 [hep-ex]
T. Aaltonen et al. [CDF Collaboration], Search for pair production of strongly interacting particles decaying to pairs of jets in \(p {\bar{p}}\) collisions at \(\sqrt{s} = 1.96 ~ \text{ TeV }\). Phys. Rev. Lett. 111(3), 031802 (2013). arXiv:1303.2699 [hep-ex]
Acknowledgements
J.A.A.S. thanks J. Alcaraz, A. Casas, M.C. Fouz, C. Muñoz and J. Terrón for useful discussions. F.R.J. is grateful for the warm hospitality and financial support from IFT/UAM-CSIC (Madrid) where most of this work was done. He also acknowledges Fundação para a Ciência e a Tecnologia (FCT, Portugal) for support through the projects CFTP-FCT Unit 777 (UID/FIS/00777/2013 and UID/FIS/00777/2019), CERN/FIS-PAR/0004/2017 and PTDC/FIS-PAR/29436/2017, which are partly funded through POCTI (FEDER), COMPETE, QREN and EU. The work of J.A.A.S. is supported by Spanish Agencia Estatal de Investigación through the grant ‘IFT Centro de Excelencia Severo Ochoa SEV-2016-0597’ and by MINECO project FPA 2013-47836-C3-2-P (including ERDF).
Author information
Authors and Affiliations
Corresponding author
Additional information
J. A. Aguilar-Saavedra: On leave of absence from Universidad de Granada, E-18071 Granada, Spain.
Scalar decays and interpretation of the CMS and LEP excesses
Scalar decays and interpretation of the CMS and LEP excesses
As seen in Fig. 8, the CMS and LEP excesses can be simultaneously fitted within our model. We present in Fig. 10 the allowed BRs of the three scalars in such case. Besides the fact that \(\text {BR}(H_3 \rightarrow b {\bar{b}})\) has to be large, as it can easily inferred from Eq. (46), the lightest scalar does not decay into \(\gamma \gamma \) and, on the other hand, \(\text {BR}(H_4 \rightarrow \gamma \gamma )\) is close to unity. Therefore, as anticipated, the cascade decays \(H_{3,4} \rightarrow H_2 H_2\) that produce the signals studied in Sect. 4 have small BRs.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3
About this article
Cite this article
Aguilar-Saavedra, J.A., Joaquim, F.R. Multiphoton signals of a (96 GeV?) stealth boson. Eur. Phys. J. C 80, 403 (2020). https://doi.org/10.1140/epjc/s10052-020-7952-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-020-7952-4