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 U(1)Y′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {U}(1)_{Y'}$$\end{document} 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 heav-J. A. Aguilar-Saavedra: On leave of absence from Universidad de Granada, E-18071 Granada, Spain. a e-mail: jaas@ugr.es (corresponding author) b e-mail: filipe.joaquim@tecnico.ulisboa.pt ier 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 U(1) Y coupling to baryon number up to an arbitrary normalisation constant. The scalar sector comprises two complex SM 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. To ensure gauge-anomaly cancellation, the SM fermion sector is minimally extended with extra leptons. 1 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. (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 Sects. 2 and 3. After writing down the relevant interactions and the H i → γ γ 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 → γ γ , 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 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 → γ γ 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σ at M h = 95.3 GeV. The ATLAS Collaboration has also looked for such process with a larger luminosity of 80 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 bb. 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.

H j → γ γ decays in the MSBM
The MSBM [10] extends the SM scalar content with two complex scalar fields χ 1 and χ 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 U(1) Y hypercharges Y and Y of the various fields are shown in Table 1, where Q i L and l i L denote the left-handed (LH) quark and lepton SU(2) L doublets of the SM, respectively. The righthanded (RH) quarks and lepton singlets are u i R , d i R and e i R . 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 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 χ 2 → −χ 2 . While m 2 0 , m 2 11 , m 2 22 , λ 0−3 and λ 5,6 are real, m 2 12 , λ 4 and λ 7−9 can be, 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. [10]. 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. [10]. 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, H i = (ρ 0 ρ 1 ρ 2 A 0 ) are related to the mass eigenstates which in the mass-eigenstate basis reads with i < j and mixing factors Notice that R i j are anti-symmetric and therefore R ii = 0, reflecting the fact that Z → H i H i is forbidden. Also, it can be shown that i< j R 2 i j = 1 due to the orthogonality of the mixing matrix O.

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 gaugeinvariant 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 Eq. (3), the corresponding mass terms for the new charged and neutral fermions in the interaction basis with We define the field rotations to the mass-eigenstate basis 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), (10) and (13), we can write the scalarfermion 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: 2.2 H j → γ γ decay widths The Feynman diagrams which contribute at one-loop level to H j → γ γ decays are shown in Fig. 1. In diagram (a) f is any electrically-charged fermion in the model, namely − 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). Notice that O 1 j is the admixture between the SU(2) L doublet and the j-th scalar eigenstate; the SM fermions do not couple to the SU(2) L singlets. For the new leptons, being the coefficients C E kk j and D E kk j those of (20). 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 [18] being p and q the momenta of the outgoing photons, while ε( p) 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 Fig. 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 [19] As usual, the decay widths are obtained using Fig. 1 One-loop diagrams contributing to H j → γ γ . In the fermion loop (a) f stands for any electrically-charged fermion (see Table 1) 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. [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 i j . 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 → H 2 H 2 and The values of the six angles θ i j 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. The latter range is enough to obtain the predictions later presented. These inputs allow us to determine tan β, m H 2 , λ 0−6 and m 2 12 . 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 = 0 and ϕ = 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. (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 u 1.6 TeV, we will neglect it in Eq. (12) 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π ]. The resulting heavy leptons masses range from 300 GeV to 10 TeV. The λ couplings are small because the four scalars are relatively light, m H i 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), V H associated production, tt H, 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. [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 H →γ γ (withγ γ = Z Z, 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 γ γ . 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 μγ γ = 1.070 ± 0.096. Restricting ourselves to the μγ γ ≤ 1 interval that has physical meaning within this model, the constraint translates into the 1σ lower limit 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 γ γ . Using BR(H 4 → γ γ ) SM = 1.77 × 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 → H 2 H 2 → 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 → γ γ . 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 2 1 j , 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 × 10 −3 , and the extra contribution to S is of the order of 8 × 10 −3 , in good agreement with the latest determinations (assuming U = 0) of T = 0.06 ± 0.06, S = 0.02 ± 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 β ≤ 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 4 j=1 O 1 j O 4 j = 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 × 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 BR(H → H i H j ) 0.1 with 0.55 BR(H → bb) 0.6 (the SM value is approximately 0.58 [27]). 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 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 M Z they dominate the Z decays with Br(Z → FF) = 36/53 0.68. In the right panel we show BR(Z → H 3 H 4 ) and BR(Z → 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. [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 → H 3 H 4 , with  (29). All points obey the constraints on the Higgs signal-strength parameters in Eqs. (31) and (33). In each plot, the label 'rest' refers to the channels not shown in the other axes 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. [9,10]. On the other hand, final states with multipronged 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.

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 imple-mented 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   31) and (33) tance for the event selection in the rapidity region |y| ≤ 0.6. We evaluate the acceptance by simulating a Z → j j 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 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 Since the SM scalar doublet 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. 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 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 → 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  [36]. 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 [35,36].
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 Fig. 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 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 φ 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 ≤ 0.01, not covering the 30 GeV benchmark mass considered. For example, the limit obtained for 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).

Multiphotons from
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 [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 → 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 [27], 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 Fig. 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.

Photons within jets: H i → H j H j → γ γ bb
In analogy with the previous case, gg → H 3 → H 2 H 2 → γ γ bb 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 → h → H H → γ γ bb, are carried out at the LHC [39,40]  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 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].

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σ 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 provide 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 Fig. 6 Parton-level distributions for Z → H 3 X, H 3 → H 2 H 2 → γ γ bb, being X some additional particle. Left: R separation between the two b quarks. Right: minimum and maximum R separation between the photons and b quarks Fig. 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 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 → Z Z * yielded masses M H = 125.98 ± 0.50 GeV and M H = 124.51 ± 0.51 GeV, respectively, with uncertainties dominated by statistics [42]. 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 Z Z * channel. The CMS excess corresponds to a signal strength μ CMS = 0.6 ± 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 γ γ . 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. (43) 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 [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 → H 3 → γ γ 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σ local excess observed at a mass of 98 GeV in LEP searches for the SM Higgs boson (e + e − → Zh, h → bb) [17]. The fact that the two excesses are located at nearly the same mass, suggests that they could originate from  (45) obtained from the parameter-space scan. The best-fit point corresponds to μ CMS = 0.6 ± 0.2 [49], μ LEP = 0.117 ± 0.057 [54] γ γ and bb decays of the same particle. The signal strength of the LEP excess, defined as As shown in Fig. 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 agreement with the results in Fig. 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. (46). 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 Fig. 9. On the top-left panel we observe that the pattern of Z

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. 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 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 [20]. We have found that the allowed Fig. 9 Branching ratios for the selected Z and H i decays. While all points (gray and black) obey the constraints on the Higgs signal-strength parameters in Eqs. (31) and (33), only the black ones are compatible with the CMS excess 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 Fig. 9, we observe that in most of the parameter space where H 3 → γ γ fits the CMS excess:  [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 → 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 [58]. 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 [59], 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 [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 → 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: . 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 → H 3 → H 2 H 2 → γ γ bb. 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 → H 2 H 2 → 4b have quite large backgrounds and small sensitivity. In W H/Z H associated production, a search for H → aa → bbbb (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, pp → Y Y → j j j j [62]. 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 Fig. 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 [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 (multipronged) 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 γ .

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This manuscript has no associated data.] 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://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

A 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