Knocking on New Physics' door with a Scalar Resonance

We speculate about the origin of the recent excess at ~750 GeV in diphoton resonance searches observed by the ATLAS and CMS experiments using the first 13 TeV data. Its interpretation as a new scalar resonance produced in gluon fusion and decaying to photons is consistent with all relevant exclusion bounds from the 8 TeV LHC run. We provide a simple phenomenological framework to parametrize the properties of the new resonance and show in a model-independent way that, if the scalar is produced in gluon fusion, additional new colored and charged particles are required. Finally, we discuss some interpretations in various concrete setups, such as a singlet (pseudo-) scalar, composite Higgs, and the MSSM.


Introduction
Very recently, the ATLAS and CMS collaborations presented first results from 13 TeV proton-proton collisions at LHC Run-II [1, 2]. Intriguingly, both experiments found a resonance-like excess in the diphoton invariant mass spectrum around 750 GeV.
The CMS collaboration reported a 95% CL upper limit of 13.7 fb on the cross section times branching ratio of a narrow spin-2 resonance decaying into two photons, compared with an expected exclusion of 6.3 fb (Fig. 6 of Ref. [1]), which corresponds to an excess with a local significance of 2.6 σ [1]. When interpreting the excess in terms of the signal strength for a narrow scalar (or pseudo-scalar) resonance, based on the expected and observed exclusion limits, the CMS measurement in the Gaussian approximation reads µ CMS 13TeV = σ(pp → S) 13 TeV × B(S → γγ) = (5.6 ± 2.4) fb . (1) Moreover, the ATLAS collaboration reported the observed exclusion limit in the fiducial region for a narrow-width scalar resonance µ ATLAS 13TeV,fid < 11.5 fb, compared with an expected exclusion of 2.6 fb (Fig. 3 of Ref. [2]), showing an excess of 3.6 σ significance [2]. Using Monte Carlo simulation we estimate the acceptance of the fiducial region for scalar production via gluon fusion to be ∼ 60%. In this case, the Gaussian approximation cannot be used to estimate of the signal strength, as is clear also from the very large value of the observed limit, compared to the expected one. We therefore parameterize the likelihood with a Poissonian function, requiring the correct observed exclusion limit and local significance to be reproduced, resulting in µ ATLAS 13 TeV = σ(pp → S) 13 TeV × B(S → γγ) = 10 +4 −3 fb . (2) The CMS search for a diphoton scalar resonance [3] performed during the Run-I phase at proton-proton collision energy of 8 TeV, sets a 95% CL observed upper limit of σ(pp → S) 8 TeV × B(S → γγ) < 1.32 fb, with an expected limit of 0.69 fb (Fig. 10 of Ref. [3]), implying that a ∼ 2σ excess was observed by CMS already at Run-I. The analogous ATLAS search [4] reported an observed upper limit on the RS graviton production cross section times branching ratio of < 2.8 fb at 95% CL, with an expected one of 2.2 fb (Fig. 4 of Ref. [4]). We estimate that the limit improves only by a factor ∼ 1.3 for the scalar resonance case. Based on the expected and observed exclusion limits, the diphoton signals from 8 TeV searches in the Gaussian approximation, for a narrow-width scalar resonance, are µ CMS 8TeV = σ(pp → S) 8 TeV × B(S → γγ) = (0.63 ± 0.35) fb , µ ATLAS 8TeV = σ(pp → S) 8 TeV × B(S → γγ) = (0.46 ± 0.85) fb . ( Before presenting any further discussions, we note that there is a simple (and yet general) way to test the compatibility of 8 and 13 TeV measurements, by fitting to a single parameter: where R pp depends on the production mechanism of the scalar and, to a good approximation, is given by the ratio of the parton luminosities of the relevant initial states at the two collision energies. For instance, this ratio is 4.7, 2.5 and 2.7 for 750 GeV scalar produced via gluon fusion (gg), vector fusion ( qq + qq) and associated production ( qq), respectively [5]. Among these, the most likely possibility to reconcile the measurements at two energies is if the scalar was produced dominantly via gluon fusion. Finally, when interpreting the resonance as a narrow-width (pseudo-) scalar particle produced via gluon fusion, our combination of 8 and 13 TeV measurements leads to µ 13TeV = (4.6 ± 1.2) fb , still showing a significant excess over the SM background. Shown in Fig. 1 are the individual likelihoods at 8 (13) TeV in dashed (dotted-dashed) both for CMS (red) and ATLAS (blue), while the combination is shown in solid-black. Another important experimental input, besides the overall signal yield, is the decay width of the resonance (Γ S ). The typical diphoton invariant mass resolution of the detector at 750 GeV is approximately ∼ 10 GeV [4]. On the one hand, if the width is much smaller than the resolution, as for the SM Higgs, it remains experimentally unmeasurable. On the other hand, for a sizable width a direct measurement from the lineshape is possible. Based on the results presented in [1, 2], the observed width of the excess is easily compatible with Γ S 40 GeV. Indeed, the best-fit in the ATLAS analysis [2] is obtained for Γ S /m S ∼ 6% with 3.9 σ local significance (compared with 3.6 σ in the narrow-width approximation). It is clear that the present data is not yet conclusive and a precise measurement of the width would require further analysis. Until then, we study Δχ 2 Figure 1: Reconstructed likelihoods as a function of µ 13TeV of the diphoton resonance searches at 8 TeV (dashed) and 13 TeV (dotted-dashed) by CMS (red) and ATLAS (blue). The combination (solid-black) is obtained when interpreting the resonance as a narrow-width (pseudo)scalar particle produced via gluon fusion. the phenomenological implications for both cases: negligible and sizable width (compared to the resolution).
In the rest of the paper we entertain the possibility that the excess is due to a new scalar (or pseudo-scalar) resonance. In Sec. 2 we introduce a simplified effective framework to parametrize its couplings to SM particles, while in Sec. 3 we explore all the relevant experimental constraints from other searches performed at Run-I. In Sec. 4 we discuss two phenomenologically different scenarios based on the size of the total decay width. Finally, in Sec. 5 we speculate on possible interpretations of the excess in some concrete new physics models, after which we conclude.

Simplified characterization framework
Let us add to the Standard Model (SM) a neutral scalar resonance S with mass m S . The relevant effective operators for gg → S → γγ are where G a µν and F µν are the SU(3) C and U(1) QED field strength tensors, respectively, while α s and α are the strong and electromagnetic coupling constants. Analogously, if the resonance is a pseudo-scalar, the effective Lagrangian describing its interactions with SM particles remains the same as in Eq. (6) after the substitution c G,γγ →c G,γγ and X µν X µν → 1 2 µνρσ X µν X ρσ , where (X = G, F ). We assume that the dominant production mechanism of the scalar S at hadron collider is gluon fusion (gg → S), induced by c G . As anticipated in the literature, the leading order cross section is expected to receive large higher order QCD correction factor [6]. We benefit from the precise computations of the SM Higgs boson production at the LHC. In fact, the QCD correction factor in the infinite top mass limit (to match our calculation) is a good approximation of the full top mass dependence even for a relatively heavy Higgs boson [6]. In the numerical analysis below, we fix m S = 750 GeV and estimate the total cross section using the available NNLO QCD predictions for the SM Higgs [6,7] and rescaling these for the heavy m top limit (since we assume heavy mediators in the loops generating Eq. (6)). Finally, the production cross sections for a 750 GeV scalar are For the pseudo-scalar case, we note that the N 3 LO QCD corrections, computed recently in Ref. [8] for 13 TeV, provide a K-factor very similar to the one of the scalar [6], well within the scale uncertainty of the NNLO computation. The decay width of S into two photons or two gluons following from Eq. (6) is given by where we include the NLO QCD correction factor K F = 1 + 67α S 4π from Ref. [6], and the running coupling constant at the appropriate scale, α s ≡ α s (m S ). Note that the leading order partial decay widths for the pseudo-scalar case are the same, with the substitution c X →c X . Since the c γγ and c G effective couplings are expected to be generated at loop level (c γγ,G ∼ O(1)), any tree-level coupling of S to lighter states is typically expected to substantially increase the total decay width to ∼ O(GeV), thus strongly reducing the branching fraction in two photons, in which case a larger value of both c G and c γγ is needed to fit the excess. In the Sec. 4.2 we quantify this statement.
In order to assess the allowed parameter space for a generic model with a 750 GeV neutral scalar, we parametrize its tree-level couplings to SM particles as follows: The partial decay widths which follow from these couplings are where β X = 4m 2 X /m 2 S . Also,m f is the MS mass that should be evaluated at the scale m S and ∆ QCD are the higher order QCD corrections [6]. In the W W and ZZ decay formulas given in Eq. (12) we neglected the contributions from loop-induced couplings, since these are expected to be subleading whenever the tree-level couplings are present. The tree-level couplings to top and W induce at 1-loop the effective Sγγ and Sgg interactions. Their contribution to the decay widths and production cross section can be obtained with the substitution Let us analyze the question if the excess can be accommodated by top coupling only, without invoking further contributions to c G and c γγ . Using the expressions derived in this section, we find that the required coupling is c t 50. This, however, would imply an unphysical partial decay width Γ(S → tt) ∼ 8 TeV. We also note that the situation cannot be ameliorated by considering non-zero values of c V , nor by introducing a large coupling c b to the bottom quark as well. We therefore exclude this possibility and conclude that new colored and charged particles inside the loop are necessary in order to explain the excess. As we will see later on, these particles are expected to be light in order to accommodate the excess, potentially within the reach of LHC.
In addition to the decay channels listed in Eq. (12), we also consider a possible invisible decay width Γ inv . This could be due to decays into dark matter particles, or into particles which escape detection, or for which no experimental bound is present.

Experimental constraints
The framework introduced above can be employed to analyze other potentially relevant experimental constraints in a model-independent way. In Table 1 we summarize the LHC Run-I limits on σ × B for a given decay channel for a narrow 750 GeV neutral scalar (pseudo-scalar) resonance (Γ S /m S < few %). The extended discussion is given below.
The CMS search for a dijet resonance [9] at 8 TeV and 18.8 fb −1 , optimized for the mass window between 500-800 GeV, imposes a 95% CL upper limit on the production of a RS graviton decaying to gg: where A is the acceptance. We conservatively assume A = 1. The ATLAS [11,13] and CMS [10] searches for a scalar resonance decaying to ZZ and W W with the full data set, combining all the relevant Z and W decay channels, impose a 95% CL upper limits of On the other hand, the ATLAS collaboration [12] has performed a search for resonance decaying to γ and Z(→ + − ). From Monte Carlo simulations we estimate that ∼ 70% of all events fall into the fiducial region, and we interpret the search as a 95% CL upper limit on the inclusive σ × B, The decay into a Higgs boson pair is also subject to constraints. Several searches have been performed during the LHC Run-I by ATLAS [15] and CMS [14]; for a resonance mass of 750 GeV, the most stringent constraints come from the 4b channel, where each of the 125 GeV Higgs bosons decays into a bb pair. One has This constrains the B(S → hh) at a level similar to the branching ratio into vector bosons, i.e. a few 10 −2 for a production cross cross sectionsection of the order of a few picobarn. Decays into quarks are experimentally less constrained. In particular, ATLAS [17] and CMS [16] searches for resonance decaying into a pair of top quarks result in the following constraints: Other fermionic channels, such as bb, τ + τ − , or light quarks and leptons, are less relevant if one assumes Yukawa-like couplings to the new resonance as in Eq. (11).
An interesting possibility is that S predominantly decays to invisible particles that might constitute (part of) the observed dark matter in the universe. This scenario, on the other hand, might lead to sizable mono-jet (missing energy plus jet) signatures at the LHC. The present ATLAS search [18] sets a bound on σ(pp → H) × B(H → inv) for a heavy Higgs-like particle, but only up to a mass of 300 GeV. We perform a Monte Carlo simulation of the signal for 750 GeV applying the same cuts as in the analysis. We find that the acceptance times efficiency improves by a factor of ∼ 3 with respect to 300 GeV case. We use this to estimate the bound on a 750 GeV scalar Negative results from 8 TeV resonance searches, together with the positive signal in γγ channel from Eq. (5), imply model-independent upper limits on B(S → XX)/B(S → γγ), for a given channel XX. These are illustrated in Fig. 2 for XX = Zγ, ZZ, W W , hh, tt, gg and invisible.

Phenomenological scenarios
With this experimental results in mind, we face two different classes of scenarios, depending on whether or not the resonance has extra tree-level decay channels, that would lead to a sizable decay width. Since this is an open issue in the present analyses [1, 2], we study both cases separately.
• Only loop-induced decays: The dominant decay channels of S are into the SM gauge bosons (gg, γγ, Zγ, ZZ,pseudo-scalar and W + W − ) through loops of heavy states of the new sector, described effectively by in SU(2) L × U(1) Y invariant way. In terms of the mass eigenstates, the relevant part of the above Lagrangian reads where Analogous operators can be written in the pseudo-scalar case using the substitutions in Eq. (7), In this scenario, the relevant phenomenology depends entirely on c G , c W and c B (or c G ,c W andc B ).
• Sizable extra decay channels: S is allowed to couple at tree level to other lighter particles (either SM or BSM ones) and, therefore, Γ S is expected to be dominated by extra channels, rendering it largely independent on c G and c γγ . In this case, the phenomenological parameters relevant to the observed excess are c G , c γγ , and Γ S .

Phenomenology of the "loop-only" scenario
The loop-induced decay channels generated by the effective couplings in Eq. (21) are dominant only if any tree-level couplings to lighter states are strongly suppressed. Starting with Eq. (21), we compute the partial decay widths for S → V V , where V V = γγ, ZZ, W W , Zγ and gg, as a function of c G , c B , and c W . The total width is then simply given by the sum of the partial widths, for m S = 750 GeV. The branching ratio into two photons, B(S → γγ), is also a function of the same three coefficients only. In Fig. 3, we show 68 and 95% CL preferred regions from the combined diphoton signal in Eq. the effective Sgg coupling, |c G | 12. On the other hand, this scenario also predicts a correlated signal in ZZ, Zγ, and W + W − channels: where R W B = c W /c B . The best present limit is due to the Zγ channel (see Fig. 2) with a lower limit R W B −1.7 and a loose upper limit. The correlation induced via the single parameter R W B is a striking prediction for future searches at LHC Run-II. In Sec. 5 we discuss the interpretation of the excess in some concrete models which fall into this class.

Phenomenology of the "extra-width" scenario
In this section we entertain the possibility that the decay width of the new scalar is dominated by other decay channels than the loop generated ones described previously, in which case the total decay width Γ S becomes an independent free parameter. This implies that the branching fraction in two photons, B(S → γγ), is a function of c γγ and Γ S only. Within this class of scenarios, two possibilities could be realized: • either the width is bigger than the experimental resolution (Γ S 10 GeV) and could then be measured directly,  Figure 4: Preferred region at 68% (green) and 95% (yellow) CL by the diphoton excess from Eq. (5) and constraints from dijet production in Eq. (14) in (c G , c γγ ) plane for Γ S = 1 GeV (upper plot) and Γ S = 20 GeV (lower plot). Shown in black-dotted is the gluon fusion production cross section at 8 TeV. The blue line is the prediction for a composite vector-like fermion bidoublet (2, 2) 2/3 . See text for details.
• or it is still smaller than the resolution, Γ S 10 GeV, but larger than the "looponly" contribution, Γ loop S < Γ S .
As a benchmark point, we assume the total decay width to be either Γ S = 20 GeV or Γ S = 1 GeV, and repeat the survey of the relevant phenomenological constraints. In this case, the diphoton signal strength at 13 TeV is where one can notice the simple scaling with Γ S . We show in Fig. 4 in green (yellow) the 68 % (95 %) CL region preferred by the combined fit to diphoton resonance searches from Eq. (5). The main difference with respect to the previous case, is that the c G,γγ couplings have to be large in order to fit the diphoton excess. This should also be confronted with the dijet bound from Eq. (14), which corresponds to a constraint |c G | 13 if Γ S = 1 GeV as shown in Fig. 4 (top) and |c G | 28 if Γ S = 20 GeV, as shown in Fig. 4 (bottom). According to Eq. (8), the large values of c G enhance the production cross section, therefore it is important and non-trivial to point out to which final states S is allowed to decay and compare with the present experimental constraints discussed in Sec. 3. The total production cross section σ(pp → S) isolines at 8 TeV are shown with black-dashed vertical lines in Fig. 4. The 95% CL constraints on tree-level decay channels from 8 TeV resonance searches are Comparing these bounds (Table 1)

Explicit models and insights
In this section we provide interpretations of the generic scenarios described above, in the context of some more concrete models.

Scalar singlet with loop decays only
In this class of models the scalar resonance is a SU(2) L × U(1) Y singlet, which decays predominantly in SM gauge bosons via the couplings in Eq. (21). This feature can be naturally achieved by assuming that S is the lightest state of a new heavy sector, coupled to the SM only via the SM gauge interactions, and that it does not participate at tree level to the electroweak symmetry breaking, i.e. there is no mixing with the SM Higgs boson. In this case the operators in Eq. (21) are generated via loops of (colored and charged) heavy states of the new sector. Given that the main phenomenological constraints concern only the gg → S → γγ process, the relevant properties of the new particles inside the loop are the spin and the quantum numbers under SU(3) C and U(1) QED . As a simple benchmark scenario, we assume S is coupled to a set of vector-like heavy fermions Ψ i , triplets (or singlets) of SU(3) C with electric charge Q i and mass M i , via marginal (Yukawa type) operators, At one loop they generate the necessary couplings of S with gluon and photon pairs. In the limit of heavy fermions (τ Ψ i = m 2 S /4M 2 i 1), we can match the c G , c γγ and c Zγ coefficients to the model parameters, where N c i = 3 (1) for the color triplets (singlets) and T 3 i is the SU(2) L isospin quantum number. For completeness, c Zγ parameterizes the effective scalar coupling to Zγ. Using the relations in Eq. (23), the matching in Eq. (30) can easily be translated to c B and c W .
Based on the analysis in Sec. 4.1, the "loop-only" scenario points to relatively small effective couplings c G , c B and c W , which can in principle be due to a single (or few) particle(s) in the loop. For example, consider a single vector-like quark representation which is a singlet under SU(2) L and has electric charge Q f = Y f : In Fig. 3 (top), we show in solid-blue the predictions for electric charges 5/3, 2/3 and 1/3. The first two, in particular, can nicely explain the observed excess for reasonable values of g * ∼ O(1) and M f ∼ O(TeV). If this is correct, one can expect vector-like quark signatures to show up at the LHC.

Singlet mixed with the Higgs
The previous scenario can be generalized allowing also for other decay channels of the singlet into SM particles, which can arise at a renormalizable level through mixing with the Higgs. Indeed, the cubic term S|H| 2 will in general be present in the scalar potential, if no particular assumption is made in order to suppress it, and it will give rise to a mass mixing between the two CP even states after electroweak symmetry breaking. A situation of this kind can be found in the context of many most natural extensions of the SM. The effective Lagrangian for the singlet is then given by Eq. (11), with where θ is the singlet-Higgs mixing angle. For sufficiently high masses, the main decay widths of S are into W , Z, and Higgs bosons, in an approximate ratio dictated by the equivalence theorem. The exact value of Γ(S → hh), which is determined by c hm , depends on the details of the scalar potential; see e.g. [19]. If no other relevant decay modes are present, the branching ratios into W + W − and ZZ are close to 0.5 and 0.25, respectively. The total width in this case is Γ S = sin 2 θ Γ SM , where Γ SM 250 GeV is the width of a SM-like Higgs of 750 GeV. One can then see from Fig. 2 that, due to the bound on S → ZZ, a branching ratio into γγ of at least about 2% is needed in order to reproduce the observed signal if B(S → ZZ) ∼ 25%. Independently of the total decay width, this can be recast as a bound on the mixing angle sin θ 2 × 10 −3 c γγ . The strongest experimental constraint on the mixing angle from the LHC Run-I Higgs coupling analysis, sin 2 θ < 0.2 [20], is several orders of magnitude weaker for physically motivated values of c γγ .

Pseudo-scalar singlet
Similar couplings to photons and gluons can be generated also for a pseudo-scalar singlet S, with a coupling to the heavy fermions of the type L Y uk ⊃ −g * i SΨ i iγ 5 Ψ i . In this case, the matching toc G andc γγ is similar to Eq. (30) with g * i → 3/2g * i . A nice feature of the pseudo-scalar scenario is that the CP symmetry automatically forbids a mixing of this particle with the SM Higgs, thus also forbidding tree-level decays to SM gauge bosons. Axion-like particles fall in this class of models, an early analysis of collider bounds in this context can be found e.g. in Ref. [21].
In addition to this, if the (pseudo-scalar) singlet is one of the pseudo-Nambu-Goldstone bosons (pNGB) of a non-minimal composite Higgs scenario in a symmetry-breaking pattern G/H, one expects that a Wess-Zumino-Witten term in the low-energy effective theory is generated; see e.g. Refs. [22,23]. This effectively provides couplings of the pseudo-scalar singlet to SM gauge bosons exactly as in Eq. (24), with the matching where f is the scale of the spontaneous symmetry breaking in the strong sector and n G,W,B are O(1) coefficients which depend on the symmetries and fermion content of the underlying UV theory. Using Eq. (23) for the pseudo-scalar one obtainsc γγ = (n B + n W )m S /f . This contribution, together with those from loops of heavy fermions discussed above, could easily match the observed excess for O(1) values of the n B,W,G parameters (see Fig. 3). In this context, measuring the n B,W,G parameters could offer an insight into the UV structure of the strong sector [23], in the same way as measuring the π 0 → γγ decay width offered insights on the structure of QCD. A simple model which provides a singlet pNGB, as well as a solution to the electroweak naturalness problem, can be found in the context of non-minimal composite Higgs models, for example those based on the spontaneous symmetry-breaking pattern SO(6)/SO(5) [23][24][25][26]. In particular, it has been shown [25,26] that a ∼ 750 GeV singlet pNGB can be accommodated in such models. Even though in this case the UV anomaly is such that n G = 0 and n B = −n W , 1 which gives vanishingc γγ andc G , the necessary non-zero contributions to these coefficients in order to explain the excess can be obtained from loops of SM fermions or vector-like top partners, generically predicted in these setups.

Composite resonance
The "extra-width" scenario points to large effective couplings which, based on perturbativity arguments, require several fermion representations to coherently contribute to gg → S. Interestingly enough, we find this to be plausible in the context of composite Higgs models if S is a composite scalar singlet resonance of the strong sector. In this example we focus on the minimal composite Higgs model SO(5)/SO(4) [27], where such a resonance could play a role in the unitarization of W W scattering at high energy [28]. Its mass is naively expected to be near the strong coupling scale Λ ∼ (few) TeV, unless it is protected by a symmetry. Nevertheless, the resonance could still be light due to some accidental cancellations or peculiarities of non-perturbative dynamics. For example, in QCD the σ meson (or f 0 (500)) is much lighter than the typical mass scale of the other resonances, even though it is not a pNGB like the pions.
In order to generate sizable values of c G consider, for example, a color triplet vector-like fermion resonance that transforms as a bidoublet (2, 2) 2/3 under the SU(2) L × SU(2) R × U(1) X global symmetry of the strong sector. The four mass eigenstates have electric charges (Q = T 3 L + T 3 R + X): 5/3, 2/3, 2/3 and −1/3. The composite scalar S couples to the bi-doublets via strong sector coupling g * as described in Sec. 5.1. Using Eq. (30), we find In Fig. 4, we show in solid-blue the correlation in the (c G , c γγ ) plane. The excess can easily be accounted for reasonable values of the parameters, e.g. for Γ S = 20 GeV, M f ∼ 1 TeV, and g * ∼ 3. Such a composite singlet resonance is expected to couple strongly with the pNGBs of the model, in particular, the Higgs and the longitudinal polarizations of the W and Z bosons. In terms of the four pNGBs π a , the coupling reads L S,CH ⊃ a S (∂ µ π a ) 2 S/f . In terms of the physical states this interaction corresponds to [28] where f is the scale of the spontaneous breaking of the global symmetry in the composite sector, f ∼ 1 TeV, and a S is an O(1) parameter. In particular, notice that the derivative coupling with the Higgs is given as c h∂ = c V = a S m S f . With this matching, and using Eq. (12), the decay width into the pNGBs (h, W , and Z) is The strongest constraint from LHC Run-I resonance searches comes from the ZZ decay channel, Eq. (15). It imposes a bound of For c G = 5 (10), as suggested in Fig. 4 top (bottom) for the bidoublet representation, this corresponds to B(S → pNGB) 16 (4)%. This implies that, for reasonable values of the parameters (i.e. assuming no extremely large contributions to c γγ ), a sizable total width can only be obtained via decays to other channels, such as tt or invisible. Therefore, one should expect to observe the signal in these channels soon if this scenario is indeed realized in Nature. In other words, the upper limit on Γ(S → ZZ)/Γ(S → γγ) from Fig. 2, corresponds to a S m S f < 0.6 × 10 −2 c γγ . The fact that the coupling to pNGBs has to be suppressed with respect to the naive expectation a S ∼ 1, puts in some tension this scenario as a natural interpretation of the excess.

A second doublet: the MSSM and beyond
Extra Higgs bosons below the TeV scale are naturally predicted in supersymmetric models. A prime example is the second doublet of the Minimal Supersymmetric SM (MSSM), on which we now focus.
The mass matrix of the CP-even Higgs system of the MSSM contains three free parameters -two masses and one mixing angle, or equivalently m A , tan β, and the well-known top-stop radiative correction ∆ t . Identifying the 750 GeV resonance with the CP-even component of the heavier doublet, the masses, mixing, and couplings of all the Higgs states are determined as functions of tan β alone, which remains the only free parameter of the model.
The mixing angle between the two doublets, in the basis where one of the states takes all the vacuum expectation value, reads (see e.g. [29]) where m H = 750 GeV and m h = 125 GeV are the masses of the two physical states, and the radiative correction ∆ t is determined as a function of the masses and tan β. The mixing is largest, δ 0.3, for tan β = 1, which is close to the edge of future sensitivity of the high-luminosity LHC to modified Higgs couplings [30]. Neglecting loop effects due to new (supersymmetric) particles coupled to the Higgs bosons, the production cross section and branching ratios of H are also determined. We have already shown in Sec. 2 that with modified couplings to SM particles alone it is not possible to reproduce the diphoton excess. In this simple case, one finds the highest values of the γγ signal strength, at the level of only ∼ 10 −2 fb, for very low values of tan β, where also the production cross section is the largest.
Contributions from additional (supersymmetric) particles are required in order to further enhance the γγ rate. For low tan β the width of the 750 GeV state reaches 10-20 GeV, and Fig. 4 (bottom) shows that very large effective couplings to gluons and photons are needed in this case. Furthermore, since the branching fraction into tt is close to 1, also direct constraints from tt resonance searches are relevant, and require c G 7 and therefore c γγ 40. On the other hand, the width of H reaches its minimum of around a GeV for tan β 6-8. In this case, as can be seen in Fig. 4 (top), smaller values of c G and c γγ are needed. However, loops of top squarks can increase the sizable only marginally [31].
B(H → γγ) can get a more significant contribution from light charginos [31], but still at a level which is not enough to reproduce the observed signal strength.
It therefore looks difficult to accommodate the observed excess in the MSSM. A generic two Higgs doublet model of type II, on the other hand, with the addition of new light charged and colored states coupled to the Higgs bosons, can easily accommodate a diphoton signal compatible with the excess. In this case one simply has to satisfy all the constraints as described in Sec. 4.2, and the c G and c γγ coefficients can be estimated as in Eq. (30).

Conclusion
After analyzing the very first data at 13 TeV collision energy, both ATLAS and CMS collaborations have recently reported a tantalizing excess in the diphoton invariant mass spectrum around ∼ 750 GeV. When interpreting the excess as a scalar (or pseudo-scalar) resonance, produced dominantly via gluon-gluon fusion, no tension with 8 TeV analyses is found. On the contrary, as shown in Fig. 1, a slight excess observed at 8 TeV in CMS nicely dovetail with the recent excess, contributing to the combined signal strength We introduced, in Sec. 2, an effective parameterization of scalar and pseudo-scalar interactions with the SM fields, and computed the relevant production cross sections and decay widths in terms of the effective couplings. Assuming the production cross section to be dominated by gluon fusion, we showed that the top quark and W contributions in the loop are not sufficient to explain the excess, requiring new colored and charged particles to exist. In Sec. 3, we did a survey of all potentially relevant resonance searches for a neutral scalar at the LHC Run-I, summarizing the limits on σ × B in Table 1. These, on the other hand, imply upper limits on the size of potential decay modes as shown in Fig. 2. Working in the effective framework, we identified, in general terms, two phenomenologically distinct scenarios based on the assumptions as regards the total decay width, to which we refer as "loop-only" and "extra-width" scenario.
In the "loop-only" scenario, the main assumption is that the new resonance mainly couples to the SM gauge bosons at loop level, and thus the total decay width is in the MeV range. The excess, in this case, can easily be explained for O(1) effective couplings c G and c γγ , while remaining in agreement with all other data. In Sec. 5, we give two specific examples that match to this scenario, namely: a model with a single vector-like quark generating the effective ggS and γγS couplings, and a setup in which the pseudoscalar singlet, along with the Higgs, arises as a pseudo-Nambu-Goldstone boson of some spontaneous symmetry breaking, in which case the required couplings are generated by a combination of Wess-Zumino-Witten terms and loops of composite fermions.
The "extra-width" scenario assumes that there are additional tree-level decay channels that dominate over the loop induced ones. We investigate the phenomenological implications specifying Γ S = 1 or 20 GeV. In the latter case, employing the limits from Table 1, we argue that the resonance can not dominantly decay to the SM gauge bosons, and identify the tt or monojet signatures at the LHC as a possible way out. This scenario can be realized in composite Higgs models (see Sec. 5.4), where large c G and c γγ couplings are obtained via interactions with composite resonances in the large representations (or large multiplicity) of the global symmetry of the strong sector.
In the intermediate regime (Γ S = 1 GeV), somewhat smaller photon and gluon couplings are required to fit the excess. However, also in this case strong constraints on the other couplings apply, making tt and invisible again the channels in which a sizable branching ratio is allowed. An example of a scenario of this kind could arise in a type-II two Higgs doublet model, like the MSSM. Here, knowing the mass of the heavier doublet leaves only one free parameter that determines the phenomenology at tree level. The total width of the heavier doublet is in the few GeV range for moderate values of tan β, but larger widths can also be attained. Even in the optimal case of smallest total width, we estimate that loop corrections from supersymmetric particles are too small to explain the observed excess, calling for an interpretation beyond the MSSM.
Even though it is still too early to draw definite conclusions about the existence of a new resonance, our analysis shows that both experiments consistently point to a sizable excess at an invariant mass of around 750 GeV. Moreover, such a resonance could fit well in many reasonable scenarios beyond the SM. In all cases other light particles are required and interesting signatures are predicted to show up during the rest of Run-II at the LHC.