Multiphoton signals of a (96 GeV?) stealth boson

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.


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 SU(3) × SU(2) L × 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]. 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 [7].
In ref. [8] 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 U(1) Y coupling to baryon number up to an arbitrary normalisation constant. To ensure gauge-anomaly cancellation, the SM fermion sector is minimally extended with a set of vector-like leptons containing a doublet (N 2 E 2 ) with SM hypercharge Y = −1/2 plus vector-like singlets N 1 , E 1 with charges 0 and −1, respectively. 1 The scalar sector comprises two complex singlets χ 1,2 that, upon 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.
In the MSBM one can accommodate sizeable branching ratios (BRs) for the H i → γγ decays of the new scalars since the corresponding one-loop amplitudes may be enhanced by the Yukawa couplings with the new charged leptons E 1 , E 2 . (We denote the four scalar mass eigenstates as H i , i = 1 − 4, being H 1 ≡ H the 125 GeV Higgs boson.) This will be shown in sections 2 and 3. After writing down the relevant interactions and the H i → γγ decay widths in section 2, we perform a scan over the parameter space of the model in section 3 to compute the BRs for H i → γγ, as well as for the other decay modes. We remark that although sizeable BRs for H i decays into γγ are naturally accommodated within the MSBM, they are not a straightforward implication of the model. Besides the pp → H i → γγ 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 section 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 → γγ 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 [11] 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σ at M h = 95.3 GeV. The ATLAS Collaboration has also looked for such process with a larger luminosity of 80 fb −1 [12] at 13 TeV, without finding a significant deviation from the background-only expectation at the same mass. However, as noted in ref. [13], 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 section 5. For completeness, we also explore the possibility that the same scalar decays mostly into bb. In such case, an excess observed at the Large Electron-Positron (LEP) collider around the same mass [14] can also be accommodated. This possibility is further explored in appendix A. To conclude, in section 6 we discuss our results in light of current and future experimental searches.

H j → γγ decays in the MSBM
The MSBM [8] extends the SM scalar content with two complex scalar fields χ 1 and χ 2 , 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 U(1) Y hypercharges Y and Y of the various fields are shown in table 1, where Q iL and l iL denote the left-handed (LH) quark and lepton 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 [9]. 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) where Φ = (φ + φ 0 ) is the SM Higgs doublet. The terms in V Z 2 (V Z 2 ) conserve (break) a Z 2 symmetry under which only χ 2 transforms non trivially as , λ 0−3 and λ 5,6 are real, m 2 12 , λ 4 and λ 7−9 can be, in general, complex. We define where ρ i and η i are real fields and with v = 246 GeV. As usual, we will also use A detailed analysis of the scalar potential and scalar mass spectrum can be found in ref. [8]. Upon U(1) Y symmetry breaking, the neutral Z gauge boson acquires the mass where g Z and Y χ are the U(1) Y coupling constant and χ 1,2 hypercharge, respectively. All Z -scalar interactions can be found in ref. [8].
There are two would-be Goldstone bosons in the model, namely G 0 1 = η 0 and G 0 2 = cos β η 1 + sin β η 2 . The orthogonal state A 0 = − sin β η 1 + cos β η 2 is CP-odd, being a mass eigenstate in case the scalar potential parameters are real and ϕ = 0. In general, where O is an orthogonal 4 × 4 real matrix parameterised in terms of 2 × 2 rotations as where O kl corresponds to a rotation in the (k, l) plane by an angle θ kl . We identify the 125 GeV SM Higgs boson H as being H 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 F i and w F i (with i = 1, 2, F = E, N ) are general complex numbers. Taking into account eqs. (3), the corresponding mass terms for the new charged and neutral fermions in the interaction basis F L,R = (F 1 F 2 ) T L,R are with We define the field rotations to the mass-eigenstate basis F L,R = (F 1 F 2 ) T L,R as Here, U F L,R are 2 × 2 unitary complex matrices and the masses of the new leptons m F i are real and positive. To compute U F L,R , one defines the Hermitian matrices , being V F L,R unitary matrices given by: The mixing angles θ F L,R and CP phases δ F L,R satisfy Since, in general, V F L,R 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), (7) and (10), we can write the scalar-fermion interactions in the mass-eigenstate basis as: with (j, k) = 1, 2, a = 1, . . . , 4, and the coefficients again with (j, k) = 1, 2, a = 1, . . . , 4 and the coefficients C F jka , D F jka given by:

H j → γγ decay widths
The Feynman diagrams which contribute at one-loop level to H j → γγ decays are shown in figure 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 thef f H j interactions of eq. (16) in the form: with j = 1, . . . , 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). For the new leptons, Figure 1: One-loop diagrams contributing to H j → γγ. In the fermion loop (a) f stands for any electrically-charged fermion (see table 1).
being the coefficients C E kkj and D E kkj those of (17). As for the W W H j coupling, we have: where g is the SM SU(2) gauge coupling and m W the W boson mass.
Due to electromagnetic gauge invariance, the tensor form of the H j → γγ amplitude can be generically expressed as [15] being p and q the momenta of the outgoing photons, while ε(k) and ε(q) denote the corresponding polarisation vectors. The quantities S X j andS X j 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 figure 1, for which Here, α = e 2 /(4π) is the electromagnetic fine-structure constant, Q f is the fermion electric charge, N f c is the number of colours of f and τ X = M 2 H j /4m 2 X . The function g(τ ) is given by [16] As usual, the decay widths are obtained using where | p | = M H j /2 and a 1/2 factor has been included for identical particles in the final state.

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 λ 7−9 = 0 in eq. (1), which corresponds to having the Z 2 symmetry softly broken by the term m 2 12 = 0. The interested reader is addressed to ref. [8] 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 input values: 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 [8] 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 section 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 section 4 may take place for a wider range of masses.) M H 4 plays little role in our study and we choose it close to M H 3 . 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 The values of the six angles θ ij defining the scalar mixing matrix O in eq. (6), and of the scalar potential parameter λ 2 , are randomly varied in the intervals [0, 2π] and [−1, 1], respectively.
These inputs allow us to determine tan β, m H 2 , λ 0−6 and m 2 12 . As done in ref. [8], we will only keep those points which lead to a global minimum of the potential with v, u 1,2 = 0. Notice also that, since we are considering the case with a softly broken Z 2 symmetry under which χ 2 → −χ 2 , the couplings x 2 and x N,E 1 in eq. (7) vanish. Fermion mixing plays little role in our study and, given that v u 1.6 TeV, we will neglect it in eq. (9) implying θ F L,R = 0 and, thus, U F L,R = 1 1 2×2 . The only relevant couplings among the new fermions and the scalars are then y E,N 1 and y 2 , which we randomly vary in the perturbative interval [0, 4π].
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), V H associated production, ttH, etc.) are all rescaled by the mixing factor O 2 11 ≤ 1, implying for any of those processes. Here, BR(H → γγ) 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 µ γγ values corresponding to the different Higgs production processes, reported by the ATLAS collaboration in ref. [17], 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 H →γγ (withγγ = ZZ, W W, bb, τ τ, gg, no sum over channels), the widths are suppressed by the factor O 2 11 . Then, we have where Γ SM is the SM Higgs total width, Γ new the width into new modes (e.g. lighter scalars) and the sum in the denominator comprises all SM decay channels other than γγ.
which we require in our scan.
In figure 2 we show in ternary plots the results obtained for the decay BRs of the four scalars, corresponding to 6 × 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  (26). All points obey the constraints on the Higgs signal-strength parameters in eqs. (28) and (30). In each plot, the label 'rest' refers to the channels not shown in the other axes.
signal-strength parameters imply for the SM-like Higgs boson BR(H → H i H j ) 0.1 with 0.55 BR(H → bb) 0.6 (the SM value is approximately 0.58 [18]). As expected, the lightest scalar H 2 predominantly decays into γγ and bb (top-right panel), being the probability to decay into other modes (τ τ , cc and gg, labeled as 'rest') around 20%, at most. Notice that BR(H 2 → γγ) 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 γγ. Notice that BR(H 3,4 → γγ) = 1 (apex of the triangle) is possible, as well as BR(H 3,4 → H 2 H 2 ) = 1 (lower-right vertex).
The decay BRs of the Z boson presented in figure 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  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 → H 2 H 4 follows a similar pattern as Z → 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. [8] for model 2.
In summary, figures 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. 2 The scenario corresponding to Z → H 3 H 4 , with H 3,4 → H 2 H 2 followed by H 2 → bb provides the perfect template for the stealth-boson hypothesis, as investigated in refs. [7,8]. 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.
MadGraph5 [19] with the model implemented using Feynrules [20] and the universal Feynrules output [21]. 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 [22] set an upper limit σ(pp → Z → jj) × A ≤ 24.3 fb at the 95% CL, where A is the acceptance for the event selection in the rapidity region |y| ≤ 0.6. We evaluate the acceptance by simulating a Z → jj sample using Pythia [23] and Delphes [24], with jets reconstructed using FastJet [25]. 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 Γ Z ≤ 0.15M Z , which is not always the case. Therefore, the actual limits may be looser.
For M H i M Z (which is the case in our benchmark), and depending on m F , the Z → H i H j branching ratios are for the signals we study in the remainder of this section, with a maximum of 121 fb.
The first signal addressed, boosted diphotons, is produced from H 2−4 → γγ, when the scalars are produced from Z decay. Instead, the second and third signals arise from the (expected) dominant decay H 3,4 → 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 γγ, as it may be the case in our benchmark scenario (see figure 2). Note that in the remainder of this section we discuss the type of 'object' produced by one of the boosted scalars in Z → H i H j ; an actual event will have two of such objects, or one plus a hadronically-decaying scalar, etc.

Boosted diphotons from H i → γγ
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 η and azimuthal angle φ. 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 ∆η ×∆φ = 0.025×0.0245 for the ATLAS ECAL [26] 3 and 0.0174×0.0174 for the CMS ECAL [27]. Photons are reconstructed from 3×5 clusters in ATLAS, yielding a cluster size of ∆η × ∆φ = 0.075 × 0.124. In CMS 5 × 5 clusters are used, with a size ∆η × ∆φ = 0.087 × 0.087. In both experiments, additional rows up to ∆φ = 0.172 (ATLAS) and ∆φ = 0.609 (CMS) may be added in the case of photon conversion into an e + e − pair, which are curved in opposite directions in φ 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 ∆R ≡ (∆η) 2 + (∆φ) 2 ≤ 0.2 − 0.3 [26,27].
Let us now consider the decay Z → 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 figure 4 we show the normalised distributions of the diphoton separation in ∆φ and ∆η when H i is either H 2 or H 3 with M H 3 = 96 GeV and M H 2 = 30 GeV, together with the sizes of the ATLAS and CMS clusters. For M H 2 = 30 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 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 figure 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 φ 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 [28]. 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 ≤ 0.01, not covering the 30 GeV benchmark mass considered. For example, the limit obtained for M Z = 2 TeV, 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. (34).

Multiphotons from H
The direct production gg → H 3 → H 2 H 2 → 4γ (or, equivalently, H 4 → H 2 H 2 → 4γ) produces a signal with four well-separated photons. An ATLAS search for events with three or more photons [29] 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 → H → aa → 4γ, with m a = 10 GeV. We can conservatively use the same limit for gg → H 3 → H 2 H 2 , bearing in mind that the efficiency will be slightly smaller for M H 3 < M H . Using a cross section σ(gg → h) = 26.5 pb for a 95 GeV SM-like scalar [18], and assuming a 100% branching ratio of H 3 → H 2 H 2 , we can obtain the upper bound With our scan over parameter space we have verified that BR(H 2 → γγ) can be as large as unity fulfilling the above limit.
The cascade decay of a boosted scalar H 3 → H 2 H 2 → 4γ yields collinear multiphotons that are much alike the collinear diphotons examined above. We show in figure 5 the kinematical distributions for the minimum and maximum separation in η and φ 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.  The decay of a boosted H 3 → H 2 H 2 → γγbb produces a conspicuous jet containing hard photons. Such jet is generally narrow, at it can be inferred from the distribution of the ∆R distance between the two b quarks in figure 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

Photons within jets:
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 figure 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 [32].  Figure 7: Parton-level distributions for the photon energy fraction z γ of jets produced in Z → H 3 X, H 3 → H 2 H 2 → γγbb, being X some additional particle.

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 [11] two nearby excesses when searching for new scalars h decaying into photon pairs in Run 1 and Run 2 data, • Run 1: 2σ local significance at a diphoton invariant mass m γγ = 97.6 GeV, with 19.7 fb −1 at 8 TeV; • Run 2: 3σ local significance at a slightly lower mass m γγ = 95.3 GeV, with 35.9 fb −1 at 13 TeV.
They also provides limits on two groups of processes mediating the production of h: (a) gg → h plus tth production; (b) VBF plus V h 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σ local significance at m γγ = 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 → 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 → γγ and H → ZZ * yielded masses M H = 125.98 ± 0.50 GeV and M H = 124.51 ± 0.51 GeV, respectively, with uncertainties dominated by statistics [33]. These two measurements were compatible only at the 2σ level even if the collected samples were large enough to have 5.2σ and 8.1σ significances for the Higgs signal, respectively, in the γγ and ZZ * channel.
The CMS excess corresponds to a signal strength µ CMS = 0.6 ± 0.2 [36], 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 γγ. As the CMS excess is also seen in the VBF and V h category, the latter possibility seems quite reasonable. Since in our model the new particles do not contribute to the gg → H i amplitudes, then the cross-section ratio in the r.h.s. of eq. (37) reduces to the mixing factor O 2 13 . (As mentioned, we consider H 3 as being the new scalar candidate to explain the excess.) Taking BR(H 3 → γγ) SM = 1.43 × 10 −3 [18], we find The CMS excess has sparked some interest in the framework of SM extensions with a scalar singlet [34][35][36][37][38], also in connection with a 2.3σ local excess observed at a mass of 98 GeV in LEP searches for the SM Higgs boson (e + e − → Zh, h → bb) [14]. The fact that the two excesses are located at nearly the same mass, suggests that they could originate from γγ and bb decays of the same particle. The signal strength of the LEP excess, defined as has been computed in ref. [39], yielding µ LEP = 0.117±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. (39) is again O 2 13 . Taking BR(h → bb) SM = 0.8 [18], eq. (39) implies O 2 13 BR(h → bb) = 0.0934 ± 0.045 , As shown in figure 8, both the CMS and LEP excesses can be accommodated in the MSBM. The value of µ CMS can be either larger or smaller than the best-fit value, in Figure 8: Allowed range for µ CMS and µ LEP , as defined in eqs. (37), (39) obtained from the parameter-space scan. The best-fit point corresponds to µ CMS = 0.6 ± 0.2 [36], µ LEP = 0.117 ± 0.057 [39].
agreement with the results in figure 2. On the other hand, µ LEP cannot be much larger than the best-fit value, as it is limited by the Higgs constraints on O 2 13 -see eq. (40). Note that in our benchmark scenario, the decay H 3 → 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 µ CMS and µ LEP small would be less populated.
The scenarios we are interested in are those in which H 3 → H 2 H 2 (as well as H 4 → 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 → bb (see the appendix). When we require µ CMS = 0.6 ± 0.2, the phenomenology remains mostly the same even if the parameter space is reduced. This is shown in figure 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 → 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 figures 3 and 2, while the black points are those which are also compatible with the CMS excess, with µ CMS ∈ [0.4, 0.8]. Therefore, one can conclude that the multiphoton signals studied in section 4 are compatible with the CMS excess too.

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 [17] suggests that new scalars, if they exist, have small couplings to the weak bosons. 4 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 [8], an anomaly-free leptophobic 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 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 ≡ 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 → 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 [17]. We have found that the allowed BR(H 2−4 → γγ) span over the whole range from zero to unity, while keeping BR(H → γγ) 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 → γγ loop amplitudes can be enhanced, while those for H → γγ and, of course, gg → 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 γγ can be enhanced. If the dominant decay mode is H 3 → bb, 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 → H 2 H 2 , which is otherwise the most natural situation if there is a lighter scalar H 2 and available phase space. From figure 9, we observe that in most of the parameter space where H 3 → γγ fits the CMS excess: 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. [40][41][42]). So does the boosted diphoton which, as seen in section 4.1, will usually fail the isolation criteria for prompt photons and requires a specific selection. The decay pp → Z → 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γ, with 35.9 fb −1 of data at 13 TeV [43]. The excess reaches a local significance of 3.6σ at a mass of 2 TeV in the hadronic channel, without any excess in the leptonic channel. In a search with similar sensitivity [44], the ATLAS Collaboration did not find any excess at that mass. The aforementioned process pp → Z → H 3 H 4 , with H 3 → γγ and H 4 → H 2 H 2 → 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 [45].) 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 section 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 → 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 → H 3 → H 2 H 2 → 4γ. This signal is covered by a search for events with at least three photons by the ATLAS Collaboration using 8 TeV data [29]. To the best of our knowledge, there are no searches at 13 TeV in this final state.
(ii) gg → H 3 → H 2 H 2 → γγbb. There are searches for heavy scalars decaying into the Higgs boson in this final state [30,31] 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 → H 2 H 2 → 4b have quite large backgrounds and small sensitivity. In W H/ZH associated production, a search for H → aa → bbbb (where a is a lighter particle, e.g. an axion) by the ATLAS Collaboration [46] 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, pp → Y Y → jjjj [47]. The mass range explored M Y ≥ 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 → H 3 X → 4γ +X, with X an additional particle. The four photons are produced closely (see figure 5) and require a special object selection.
(ii) Z → H 3 X → γγbb+X. The two b quarks are produced close in ∆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 [28] are much more collinear, and that search does not cover the multiphotons produced at intermediate distances, such as those shown in figure 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γ.