Cornering diphoton resonance models at the LHC
Abstract
We explore the ability of the high luminosity LHC to test models which can explain the 750 GeV diphoton excess. We focus on a wide class of models where a 750 GeV singlet scalar couples to Standard Model gauge bosons and quarks, as well as dark matter. Including both gluon and photon fusion production mechanisms, we show that LHC searches in channels correlated with the diphoton signal will be able to probe wide classes of diphoton models with of data. Furthermore, models in which the scalar is a portal to the dark sector can be cornered with as little as .
B1348 LouvainlaNeuve, BelgiumInstitute of High Energy Physics, Austrian Academy of Sciences,
Nikolsdorfergasse 18, 1050 Vienna, AustriaTheoretische Natuurkunde and IIHE/ELEM, Vrije Universiteit Brussel, and International Solvay Institutes,
Pleinlaan 2, B1050 Brussels, BelgiumNational Centre for Nuclear Research,
Hoża 69, 00681 Warsaw, PolandInstitute for Particle Physics Phenomenology, Department of Physics, Durham University,
DH13LE, UK\preprint
IPPP/16/46, DCPT/16/92,CP31623, HEPHYPUB 967/16
1 Introduction
Both ATLAS and CMS collaborations recently announced an excess in the diphoton spectrum around the invariant mass of . While the excess is not statistically significant to claim a discovery (ATLAS finds a local significance of ATLASCONF2015081 (); ATLASCONF2016018 () and CMS one of CMSPASEXO15004 (); CMSPASEXO16018 ()), it is certainly interesting to entertain the idea that the data points to the existence of a new particle.
If such particle is a singlet under the SM gauge group, it is inevitable that the diphoton excess will be correlated with signals in other channels involving gauge bosons (e.g. , , or ). It has been shown that an excess should appear at least in one of the beforementioned channels, regardless of the underlying model parameters Low:2015qho (); Kamenik:2016tuv (). In the optimistic scenario where the 750 GeV diphoton excess remains as more data comes in, measurements of other final states which are correlated to the diphoton excess will hence become instrumental in both confirming the signal, as well as determining the properties of the new particle. In particular, not observing correlated signals in final states with Standard Model (SM) gauge bosons will have direct implications on many scenarios attempting to explain the excess.
The width of the diphoton excess offers additional crucial information about the nature of the possible new particle. The lineshape of the excess measured by ATLAS indicates a rather broad resonance with a width GeV, which is difficult to account for if it decays only into Standard Model (SM) particles. Large unobserved decay modes can point to interactions between the new resonance and dark matter, leading to collider signatures in channels with large missing energy, as well as signals in direct dark matter detection experiments via scattering off nuclei, and the measurements of galactic ray fluxes Backovic:2015fnp (); Barducci:2015gtd (); Mambrini:2015wyu (); D'Eramo:2016mgv ().
In this paper we explore the reach of LHC run2 searches for the diphoton resonance models. Leading order approximations for the production of channel resonances allow for the use of simple scaling rules to study the constraints of existing and future LHC results on the model parameter space. Specifying the features for the diphoton excess, such as the production mechanism and cross section, defines a hypersurface in the multidimensional parameter space which can explain the excess. Our approach consists in constraining these surfaces further, by imposing collider bounds on correlated final states.
Using concrete examples, we demonstrate the most sensitive channels and relevant bounds, as well as the required integrated luminosity to rule out particular models explaining the diphoton excess. For concreteness, we assume throughout the paper that the resonance is a scalar singlet under the SM gauge group. Hence, its interactions with SM particles are captured at leading order by a set of dimension5 operators suppressed by a new physics scale Franceschini:2015kwy (). We further assume that the new resonance does not mix with the SM Higgs boson, as existing and projected limits from Higgs coupling measurements set strong indirect constraints Falkowski:2015swt ().
We discuss three concrete benchmark scenarios, which serve to encompass a large class of 750 GeV diphoton resonance models. First, we study the “vanilla” scenario, in which a scalar singlet couples only to SM gauge bosons via dimension5 effective interactions. Second, we consider a scenario in which decays of a 750 GeV scalar into an invisible sector ( dark matter) accommodate the potentially large resonance width. Finally, we analyze a scenario in which the scalar is allowed to couple to SM quarks in addition to SM gauge bosons. For the purpose of studying future LHC limits on the three scenarios, we project existing 8 and 13 TeV limits on production of gauge boson, monojet, and final states at various luminosities. We outline the strategy we adopt and the simplified approach we employ to project limits for the LHC in Sec. 2. In Sec. 3 we present our main results, where we confront concrete diphoton scenarios with the existing LHC bounds and our estimated projections for the 13 TeV run. Finally, we briefly summarize our results and conclude in Sec. 4. In Appendix A we provide more technical details about limit projection and in Appendix B we review the analytical forms used here for the calculation of the decay widths.
2 General strategy and LHC limits
We begin with a brief discussion of the possible production modes for the 750 GeV diphoton resonance. We limit our discussion to the case of a pure scalar, however most of the qualitative conclusions in our paper will hold in the case of a pseudoscalar resonance as well. In the most general scenario, the onshell production cross section of the scalar resonance can be approximated by
where are proton constituents (including photons), are the dimensionless parton luminosity factors and are partonic cross sections. Limits from 8 TeV LHC disfavor production via light quarks Franceschini:2015kwy (); Gupta:2015zzs () and we will hence limit ourselves to scenarios in which the new scalar particle is produced via either gluon fusion () or photon fusion () initial state.
The photon fusion production mechanism deserves further discussion. Production of a scalar resonance compatible with the diphoton excess via photon fusion are studied in Refs. Csaki:2015vek (); Csaki:2016raa (); HarlandLang:2016qjy (); Abel:2016pyc (); Fichet:2015vvy ().^{*}^{*}*First coupling constraints for such models using 8 TeV data have been obtained in Jaeckel:2012yz (). However, it is important to note that many subtleties arise in considering the photon fusion channel. The cross section enhancement between 8 and 13 TeV centerofmass energy at the LHC is subject to large uncertainties and can vary between a factor 2 and 4 Csaki:2016raa (); HarlandLang:2016qjy (). Hence, pure photon production is possibly already in tension with 8 TeV data if the ratio is closer to 2. In addition, given the inclusive nature of the diphoton excess measurements in the ATLAS and CMS searches, it is also possible that vector boson fusion (VBF) channels with one or two additional reconstructed jets contribute to the overall production cross section. We estimated the VFB contributions with one or two additional jets for the models we consider in this paper. We found that VBF contributes at most of the inclusive diphoton production cross section in the regions of the parameter space compatible with the observed diphoton excess.^{†}^{†}†The full treatment of multijet merging in electroweak processes is beyond the scope of our paper Alwall:2007fs (). We will thus neglect such VBF contributions in the following.
Continuing, within the narrow width approximation the diphoton cross section at leading order is simply
(1)  
where we have factored out the dependence on couplings to gluons and photons ( and ). are the photon and gluon initiated production cross sections respectively. Note that is an implicit function of all of the theory parameters. Assuming a signal cross section , consistent with the observed excess, Eq. (1) can be solved for as a function of the remaining parameters in a given model, hence defining a slice of the parameter space which can accommodate the excess. Note that in the limit of the branching ratio into photons also vanishes, yielding no viable solution for .
Parameter space slices determined by can then be bound by searches in the complementary final state channels. ATLAS and CMS have recently published the first results from the LHC 13 TeV run, with an integrated luminosity of and respectively, which can be used to constrain existing models. Bounds from resonance searches involving gauge bosons final states are of particular relevance for constraining gauge invariant parameterizations of the diphoton models.
Search  8 TeV limit [fb]  13 TeV limit [fb]  13 TeV limit [fb] (expected)  
(observed)  (observed)  
11 Aad:2014fha ()  30 ATLASCONF2016010 ()  43  14  4.4  1.4  
12 Aad:2015kna ()  180 ATLASCONF2015071 ()  82  27  8.5  2.7  
40 Aad:2015agg ()  400 ATLASCONF2016021 ()  300  98  31  9.8  
460 Chatrchyan:2013lca ()  10000 ATLASCONF2016014 ()  3267  1067  337  107  
MET+  7.2 (SR7) Aad:2015zva ()  61 (IM5) Aaboud:2016tnv ()  51  
19 (IM7) Aaboud:2016tnv ()  15  5  1.5  0.5  

We present a summary of the bounds used in this paper in Table 1. In the final state the 95% C.L. 8 TeV ATLAS upper bound on the production cross section times branching ratio Aad:2014fha () reads approximately 11 fb, whereas the bound from the equivalent search in Run 2 ATLASCONF2016010 () yields . While the data from Run 2 is not particularly useful to constrain these scenarios yet, it can nonetheless be used to estimate the reach of these searches for future luminosity. The idea is based on the assumption that, while being model dependent, quantities like cross sections, acceptances and efficiencies do not depend on the integrated luminosity. In the limit of a large number of events, one can obtain the expected 95% C.L. cross section bound at any target luminosity by rescaling the limit with the ratio of corresponding luminosities. Considering, for example, the case in Table 1, rescaling the expected bound of ATLASCONF2016010 () yields the projected values shown in the columns of 30, 300, and .^{‡}^{‡}‡We stress that the limits we obtain in this way are conservative. Datadriven methods can reduce systematic uncertainties when large data samples are available and dedicated reconstruction techniques Abdesselam:2010pt (); Altheimer:2012mn (). exploiting the increased centerofmass energy at TeV and different decay mode scan improve on the limits we extrapolate.
We use the above luminosityrescaling ansatz to obtain the majority of the projections considered in this paper. However, while luminosity rescaling provides conservative estimates in most cases, it does not always reproduce the most realistic expectations. As experience with the large number of search results produced during and after the 8 TeV run has shown, a statistical combination of the data obtained in searches sensitive to different final states often leads to a dramatic improvement in the bounds with respect to searches in single channels. For instance, a direct comparison of the expected 8 TeV bounds on the production cross section of a heavy scalar decaying to in the , , , and a combination thereof Aad:2015kna () shows that the combined limit is at least a factor of two stronger than any of the individual bounds. ATLAS has published results for the 13 TeV resonance searches in the ATLASCONF2015068 () and ATLASCONF2015071 () final states, but at this early stage the combination has not been published. It is reasonable to assume that the final combined limit will be also stronger than the one obtained in Refs. ATLASCONF2015068 () or ATLASCONF2015071 (). Hence, we will adopt the 13 TeV limit extrapolated from the combined 8 TeV LHC limit, using the procedure described in detail in Appendix A. We have verified that the procedure accurately reproduces the existing 13 TeV limits in the and channels, leading us to conclude that our combined limit extrapolation is also accurate (see Appendix A for more details).
Limits on the resonant production both at 8 TeV and 13 TeV exist Aad:2015agg (); ATLASCONF2016021 (), and we adopt the observed limits on the 750 GeV resonance from both LHC runs.
The strongest observed ATLAS limits in the final state with at least one jet and large missing transverse momentum (hereafter MET+) comes from inclusive search bins denominated SR7 (in the 8 TeV search Aad:2015zva ()) and IM5 (at 13 TeV Aaboud:2016tnv ()), which are defined by . The strongest expected limit at 13 TeV comes instead from the inclusive bin IM7 with . Hence, for the purpose of extrapolating the limit to higher luminosities we use the expected limit at 13 TeV in the inclusive bin IM7.
Finally, current experimental searches for resonances at 13 TeV ATLASCONF2016014 () have focused only on the boosted regime, with no publicly available result on searches for resonances in the resolved regime. Boosted top analyses are ill suited for efficient reconstruction of the final states with invariant mass of (assuming the standard fat jet cone of radius ), resulting in 13 TeV limits on a 750 GeV resonance which are far weaker than the extrapolated 8 TeV limits in the resolved jet analysis. For 13 TeV final state, we hence adopt an extrapolated limit from the resolved 8 TeV analysis, obtained with the techniques explained in Appendix A.
3 Diphoton resonance models
In order to illustrate the strategy we have discussed in the previous section, we consider a concrete set of models where the new resonance is represented by a singlet scalar coupled to the SM with dimensionfive operators. Moreover, we also investigate the possibility that the new resonance plays the role of a portal to a dark sector. A wide class of diphoton resonance models can comprehensively be described by the interaction Lagrangian
(2) 
where , , and are the , , and field strength tensors, respectively, indicates SM fermions (of mass ), and is an invisible Dirac fermion which can play the role of dark matter. In the following we will independently study different subsets of this general class of models by switching on and off some of the couplings in Eq. (2).
Note that we assumed that the new scalar resonance does not couple to the SM Higgs boson. The coupling to the Higgs is mainly constrained by the allowed size of the mixing angle, which is bounded by LHC Higgs coupling measurements to be Falkowski:2015swt (). This already puts significant constraints on possible correlated signals of the new resonance in the Higgs final states, and we leave to future studies a detailed investigation of the LHC 13 TeV reach for these signatures.
We point out that the couplings of the scalar are chosen proportional to the quark masses, to respect minimal flavor violation. Since is a singlet of the SM gauge groups, the new couplings to SM fermions should be considered as descending from dimensionfive operators such as , which after electroweak symmetry breaking, generate the couplings in Eq. (2). The couplings with SM fermions in Eq. (2) have an extra suppression factor scaling, , for this reason. Without loss of generality, we have introduced a unique suppression scale for the various operators, which are then weighted by different couplings . For definiteness we will take throughout the paper.
As mentioned before, we will consider a combination of production mechanisms. For couplings of similar size gluonfusion is typically the dominant production mechanism. However we will explore also regions of the parameter space where photonfusion processes, which scale like , are dominating. In the case of quarkinitiated production, the dimensionless Yukawa couplings of to the quarks are suppressed by a factor , as they descend from higher dimensional gauge invariant operators. Hence, the light quark contributions are suppressed by the small quark masses, while the heavy quark ones are suppressed by small proton PDF and by the smallness of (since we are considering couplings and ). In particular, the top loop induced gluon fusion contribution to the production cross section is negligible with respect to other production mechanisms in the range of couplings that we study.
In order to estimate the production cross section for the resonance through the available processes we make use of several tools. We have implemented the model of Eq. (2) in FeynRules Alloul:2013bka () and we simulate the production of at the LHC using MadGraph5_aMC@NLO (MG5_aMC) Alwall:2014hca () with the NN23LO1 Ball:2012cx () PDF set for gluon as well as for photon PDFs. For photonfusion we consider both the inelasticinelastic as well as the elasticinelastic proton scattering processes.
Given the production cross section for the resonance, the cross sections in the various final states are determined by the branching ratios. Analytic formulas for the partial decay widths of in the model (2) are listed in Appendix B.
In exploring the parameter space of the model, our strategy relies on solving the condition (see Eq. (1)) for the coupling . After fixing the couplings and to some representative value, we present the results in the plane. For definiteness we choose but our results are qualitatively robust under change of the required cross section. We will also display the contours necessary to fit the excess, and identify the most relevant production mechanism on each region of the parameter space.
3.1 The “vanilla” model:
We start our analysis by considering the simplest version of the model capable of explaining the diphoton excess, we set the couplings to dark matter and SM fermions to 0. The so called “vanilla” model is then parameterized only by three couplings: and . We explore the parameter space in the range and for every value of we solve the equation fb for , imposing the conservative bound .
Figure 1 shows our first result. In the upper left plot of Fig. 1 we display in solid red the contours of consistent with signal cross section . The values of decrease towards larger values of and since the branching ratio into photons increases. In addition, the photon fusion contribution to the total production cross section also increases with larger and values, requiring a lower gluon fusion contribution to reproduce the signal.
The green dashed contours mark the boundary between regions where gluon fusion dominates and regions where photon fusion dominates instead. Photon fusion can be dominant only for large values of . The shape of the green dashed contour is determined by the competition between the BR() and the branching ratios for the other electroweak bosons (, , and ), which can deplete the signal in .
The gray regions indicate regions where no solutions for resulting into exist. We can identify two distinct gray areas which have different physical interpretations. The almost vertical gray stripe close to the central axis (denoted “Nosoln.”) is located around the straight line . In this regime, the coupling to photons is very small, leading to the fact that no real value of can reproduce the signal strength . The argument can be understood analytically as follows. The coupling to photons is almost vanishing in the central grey region, leading to a gluon fusion dominated production mechanism. We can then write
(3) 
where is the gluon luminosity and is the centre of mass energy. One can impose and solve this equation for the total width of obtaining
(4) 
Given that the total width of the resonance is always larger or equal than the width into gluons, , we arrive to the inequality
(5) 
which implies an absolute lower bound for the partial decay width into photons, necessary to accommodate the diphoton excess. Inserting the explicit expression for the partial width (see Appendix B) we obtain the lines which delimit the vertical gray stripe:
(6) 
The above argument does not depend on the other contributions to and is therefore a generic result for the complete model of Eq. (2), independently of the value of and . We will indeed find the same gray stripe around in all of the other scenarios considered in this paper.
The other gray region, denoted with fb in Fig. 1, are instead characterized by excessively large rates in , completely dominated by photonfusion processes. The internal border of the region identifies the line where the production mechanism is photon fusion, and .
Figure 1 does not show the 8 TeV bound on final states. In the region where gluon fusion dominates, this bound is automatically satisfied since gluon luminosity increases by a factor of , and hence a corresponds to , just below the LHC 8 TeV bound. In the photonfusion dominated regions the argument is less straightforward. Using the NN23LO1 PDF in MG5 the enhancement factor from to TeV in photon fusion is approximately and hence the photonfusion dominated regions would not be compatible with LHC 8 TeV constraints. Given the ongoing discussion in the literature about the exact value of the enhancement factor HarlandLang:2016qjy (), it is still possible that photonfusion is eventually a viable option Csaki:2015vek (); Csaki:2016raa (); HarlandLang:2016qjy (); Abel:2016pyc (). Thus, given the large uncertainties in such estimate, conservatively we do not impose any extra bound on such regions from the LHC 8 TeV final state searches. An ATLAS study of the jet multiplicity distribution in the diphoton events seems to show that the data favors production processes with a small number of accompanying jets ATLASCONF2016018 (), hence consistent with dominant photonfusion. For all of the above reasons, we choose to simply denote the region with a dashed green line, and remain agnostic on whether it is viable or not.
We proceed to investigate the bounds which are imposed by the LHC 8 TeV searches of resonances in the , , and final states. The results for the 8 TeV limits are displayed in the second top panel of Fig. 1, where as usual on every point of the plane we have solved for in order to get fb. The signal cross section in electroweak boson final states, once the signal yield in is imposed, is only a function of the ratio , which controls the relative size of the branching ratios.^{§}^{§}§Note, however, that what we are imposing is a signal cross section in at 13 TeV. In the transition from the gluonfusion to the photonfusion regime, the corresponding 8 TeV signal strength changes since the ratio of the gluon and photon luminosity is different. This effect is not visible in the shape of the regions excluded by the 8 TeV searches since effectively they always lie inside the region dominated by gluonfusion. As a consequence, the excluded region for each signature has the shape of a symmetric triangular angular slice in the plane. The strongest constraints come from the and final states. The limit instead provides inferior exclusion power for regions already bounded by the other searches. The white region is compatible with all existing LHC 8 TeV constraints and fits the 13 TeV diphoton excess.
The remaining panels of Fig. 1 show the LHC 13 TeV reach with increasing luminosity up to . It is interesting to observe that the limit at is essentially equivalent to the 8 TeV bound, while the and the are slightly weaker. Increasing the luminosity reduces the allowed parameter regions, resulting in a tiny remaining portion at . The result suggests that, if the diphoton excess is confirmed, a complementary signature in weak boson final states is highly likely to be discovered in the coming years. Notice that the projection of the 8 TeV combined limit we obtained in Sec. 2 plays a crucial role in closing almost entirely the allowed parameter region at the high luminosity LHC.
3.2 The dark matter model:
Current ATLAS results favor the interpretation of the diphoton excess in terms of a resonance with a relatively large width ( ). The large width cannot be explained by decays to gluons and photons alone. Unitarity and the existing dijet bounds exclude the coupling sizes necessary to generate the large width Knapen:2015dap (), suggesting that a wide 750 GeV resonance would have to decay to other states as well. As no new charged particles with mass have been observed at the LHC, it is reasonable to consider that the large resonance width can be explained by decays to new invisible particles. Decays of the 750 GeV resonance to neutral states are conceptually very interesting, as nonSM massive particles with no electric charge are natural candidates for dark matter.
Reference Backovic:2015fnp (); Barducci:2015gtd (); Mambrini:2015wyu (); D'Eramo:2016mgv () already considered scenarios in which a scalar with mass of 750 GeV is allowed to decay to dark matter. A generic feature appears in most models which explain the large width of via decays to dark matter: once the values for the decay width and dark matter relic density are fixed, the parameters of the dark sector (, ) are fully determined. For instance, in cases where dark matter is a Dirac fermion coupling to a pure scalar , a large width and dark matter relic density predict , .^{¶}^{¶}¶The values of the fixed parameter point can change based on the assumptions on the spin and CP properties of dark matter and .
As an illustration of the LHC prospects to probe the class of the dark matter models for the 750 GeV resonance, here we will consider a benchmark point from Ref. Backovic:2015fnp () which is allowed by the current astrophysical and collider constraints:
The first panel of Fig. 2 shows in solid red contours the values necessary to explain the diphoton excess in the dark matter model. The required values of at a fixed are significantly higher compared to the vanilla scenario of the previous section. The reason for a larger stems from the fact that in our dark matter model , requiring larger couplings to compensate for a smaller . Notice also that the photon fusion contribution to the production becomes dominant only for .
The remaining panels of Fig. 2 show the results of the current LHC exclusion of the dark matter model parameter space as well as the future prospects. The main difference compared to the vanilla benchmark model of the previous section is that allowing the 750 GeV scalar decays to dark matter introduces constraints from searches in channels with large missing energy, of which we consider MET+. The , , and results constrain the same regions of the parameter space as in the case of the vanilla model, while the MET+ channel typically provides the strongest limits, except in the corners of large . We find that current 8 TeV and 13 TeV results exclude values in the range of for , up to values of for . Future LHC results at 13 TeV will be able to exclude a majority of the parameter space with as little as of data, while with only the regions of parameter space in which photonfusion dominates will not be ruled out by MET+.
3.3 The “topphilic” model:
As a final concrete example of the diphoton models, we discuss the case in which the new scalar resonance also couples to SM fermions, we set , with nonvanishing in Eq. (2).^{∥}^{∥}∥Note that the coupling of to the SM fermions will generate extra contributions to the effective operator between and the gauge bosons (see for instance Bellazzini:2015nxw () for the case of a pseudoscalar coupled to gauge bosons and top quark). However, since we consider the same suppression scale for all dimension five operators, and all couplings and of order , such loop induced contributions will be typically subleading on the parameter space under study.
Among the various couplings to the SM fermions, the dominant coupling is to the top quark, justifying the title “topphilic”. In particular, the coupling of with the top quark will induce a sizable decay width of into pairs (see Appendix B), which now constitutes the dominant decay mode of the scalar resonance. As a consequence, in order to obtain the desired signal strength in the final state at 13 TeV, the production cross section for should be quite sizeable compared to the vanilla model of Sec. 3.1. The topphilic model is then similar to the dark matter model studied before, where the dominant invisible decay has been now substituted by a dominant decay into topantitop pairs.
We show the results of our analysis in the case of the topphilic model in Fig. 3 for one representative value at the high luminosity LHC. Indeed, note that since the coupling to SM fermions are suppressed by a factor only large values of will induce interesting effects. The only final state which distinguishes the topphilic model from the vanilla scenario of Sec. 3.1 is . The brown shaded regions in Fig. 3 illustrate the regions of the parameter space the future resonance searches will be able to probe. The example we show in Fig. 3 suggests that the topphilic model can be probed with high luminosity LHC only in the regime of . In the large scenario, the addition of the channel to the usual electroweak boson searches essentially covers the entire parameter space that we considered with 3000 of integrated luminosity at 13 TeV LHC. Note that the presence of a large coupling to the top quark, pushes the photonfusion dominated region further to larger values of compared to the vanilla model of Sec. 3.1. The reason is that a large coupling to SM fermions implies a small , resulting in the need of larger production cross section to accommodate the excess, that can essentially be obtained only via gluon fusion in the range of under consideration.
4 Summary
In this paper we have explored the LHC 13 TeV reach for models capable of explaining the diphoton excess at 750 GeV. As illustrative example we have considered a simple model with a scalar resonance coupled to SM gauge bosons, a dark matter candidate, and the SM quarks. We took into account gluonfusion as well as photonfusion as production mechanisms at the LHC. The requirement of generating the correct cross section in final state at 13 TeV imposes relations among the model parameters. We have studied the correlated signatures that can arise in such scenarios, including final states with dibosons, jet plus missing energy, and resonance searches, in order to further constrain the parameter space of the model and establish the exclusion reach of the LHC 13 TeV.
Our findings indicate that correlated LHC searches can exclude most of the relevant parameter space of a broad class of diphoton models during the second run of the LHC. The “vanilla” model (where the scalar resonance is coupled only to SM gauge bosons with dimension five operators) can be almost completely covered by associated signals in dibosons with of integrated luminosity. Concerning models where is a portal to a dark sector, we show that the monojet searches are able to corner the model with as little as . Finally, for models where the scalar resonance couples to SM quarks, the signature in the final state could provide a handle on distinguishing such scenarios from the “vanilla” model. However, in order for the signal in to be accessible, sizeable couplings of quarks to are required, as well as integrated luminosity of at least .
It would be interesting to extend our work to more exotic scenarios that can explain the diphoton excess, including nonresonant production, collimated photons, and models with non trivial coupling with the Higgs boson. The procedure we have adopted in this paper to compare with extrapolated LHC 13 TeV limits could be extended also to such scenarios. If the diphoton excess is confirmed, it becomes of utmost importance to explore the full set of correlated signatures expected to appear in the ongoing run of the LHC.
Note added:
During the final stages of this work, Refs. Sato:2016hls () and No:2016htu () appeared. Both references studied the LHC prospects for exclusion of a simplified diphoton resonance model analogous to the scenario we study in Section 3.1, and obtained results which are in agreement with ours. Compared to Refs. Sato:2016hls () and No:2016htu (), our analysis in Section 3.1 also discusses the production mechanism for the resonance, including gluon and photon fusion.
Acknowledgments
We would like to thank A. Goudelis, K. Kowalska, D. Redigolo, and F. Sala for useful discussions. M.B. is supported by a MOVEIN Louvain Cofund grant. S.K. is supported by the “New Frontiers” program of the Austrian Academy of Sciences. A.M. is supported by the Strategic Research Program High Energy Physics and the Research Council of the Vrije Universiteit Brussel. M.B. and A.M. are also supported in part by the Belgian Federal Science Policy Office through the Interuniversity Attraction Pole P7/37. M.S. is supported in part by the European Commission through the “HiggsTools” Initial Training Network PITNGA2012316704.
Appendix A Limit extrapolation
In order to project the LHC 8 TeV limits to 13 TeV, we employ a simple extrapolation algorithm, similar to Refs. Thamm:2015zwa (); Buttazzo:2015bka (). We begin with the assumption that the 8 TeV and 13 TeV resonance searches are characterized by acceptances and event selection efficiencies which are roughly equal.
The test statistics employed by the experimental collaborations to determine the 95% C.L. upper bounds on the cross section times branching ratio is a constant at different luminosities and center of mass energies:
(7) 
where is the upper bound on the number of signal events and is the expected or observed background, is the center of mass energy, the integrated luminosity, and the invariant mass bin.
Assuming that the background is dominated by a single initial state production mode (which is a decent approximation in most cases) we can write at any given and :
(8) 
where is the parton luminosity ratio and stand for quarks and gluons.
Inserting this in Eq. (7)
(9)  
where in the last line we have used the fact the the is a constant, see Eq. (7).
In the limit of a large number of events the event ditribution becomes well approximated by a Gaussian, so that the equality between the secondtolast and last line of Eq. (9) can be written as
(10) 
Moreover, because of our initial assumption that the efficiencies and acceptances are the same, scales as , so that one solves Eq. (10) to get
(11) 
Equation (11) represents our “master formula” for limit extrapolations. The parton luminosity ratios have been previously calculated in Ref. jsterling (). For completeness, here we give a numerical polynomial fit to the parton luminosity ratios for , and initial states, valid in the range of :
(12)  
where .
We find that when used to extrapolate the expected 8 TeV limits, the extrapolation formula of Eq. (11) gives results which are within from the true expected limits at 13 TeV. In order to validate the procedure, we have compared the results using Eq. (11) to a number of already public ATLAS results from 13 TeV. Table 2 shows the results. The largest error in our limit extrapolation is , in the case of the and searches. This is mostly due to the fact that ATLAS does not provide for those searches the efficiencies for all bins, and a knowledge of the latter is required to extrapolate the cross section bound for the fiducial cross section bound. The average error is about 10%. The uncertainty in the limit extrapolation does not strongly affect our results on the parameter space exclusion. Figure 4 illustrates the result in case of the cross section, extrapolated from the limit to 13 TeV with . The blue, shaded region shows the excluded parameter space, while the dashed regions show where the edge of the exclusion would lie if the maximal cross section was different.
Although Eq. (11) gives reasonably accurate results in many cases, it is important to point out where it fails. If the event reconstruction and selection efficiencies and acceptances differ significantly between 8 TeV and 13 TeV, Eq. (11) can result in errors larger than . The approximation is also not accurate when the 8 TeV expected background is a number of the order of a few units, so that the event distribution is not well approximated by a Gaussian, but rather presents a longer tail.
Another scenario in which the extrapolation of Eq. (11) fails are nonresonance searches ( MET+) or searches for broad resonances. In cases where the signal cross section is not distributed mostly in a narrow range of invariant masses (such as in the case of a narrow resonance), it is inappropriate to use a parton luminosity ratio evaluated at a single . Instead, an integral value over the parton luminosities is more appropriate, as the signal cross section will be distributed over a wider range of .
F.S.  Ref.  [GeV]  I.S.  (13 TeV)  % diff.  

8 TeV Aad:2015wra ()  300  220 fb, 20.3  1250 fb, 3.2  952 fb  27 %  
400  92 fb, 20.3  500 fb, 3.2  423 fb  17 %  
13 TeV ATLASCONF2016015 ()  750  16 fb, 20.3  73 fb, 3.2  87 fb  +17 %  
1000  10 fb, 20.3  50 fb, 3.2  61 fb  +20 %  
8 TeV Aad:2014fha ()  400  0.5 fb, 20.3  2.3 fb, 3.2  1.8 fb  24 %  
750  0.2 fb, 20.3  1.2 fb, 3.2  0.9 fb  29 %  
13 TeV ATLASCONF2016010 ()  1600  0.1 fb, 20.3  0.6 fb, 3.2  0.6 fb  0 %  
8 TeV Aad:2014cka ()  500  3.2 fb, 20.4  11 fb, 3.2  12 fb  +9 %  
750  1.2 fb, 20.4  4.8 fb, 3.2  4.9 fb  +2 %  
13 TeV ATLASCONF2015070 ()  1500  0.4 fb, 20.4  1.6 fb, 3.2  1.9 fb  +17 %  
8 TeV Aad:2014xka ()  750  48 fb, 20.3  200 fb, 3.2  197 fb  2 %  
1000  19 fb, 20.3  105 fb, 3.2  85 fb  21 %  
13 TeV ATLASCONF2015071 ()  2000  6.0 fb, 20.3  38 fb, 3.2  41 fb  +8 %  
8 TeV Aad:2015uka ()  600  22 fb, 19.5  110 fb, 3.2  110 fb  0 %  
800  9 fb, 19.5  60 fb, 3.2  49 fb  20 %  
13 TeV ATLASCONF2016017 ()  1400  3.9 fb, 19.5  22 fb, 3.2  28 fb  +24 %  
8 TeV Aad:2015mna ()  500  4.1 fb, 20.3  22 fb, 3.2  20 fb  10 %  
750  2.0 fb, 20.3  8.2 fb, 3.2  10.9 fb  +28 %  
13 TeV ATLASCONF2015081 ()  1500  0.5 fb, 20.3  3.9 fb, 3.2  3.8 fb  3 % 
Appendix B Analytical form of the decay widths
In this appendix we report the analytic formulas for the partial decay widths of the resonance. The following expressions were used for the analysis discussed in the main body of the paper.
(13)  
(14)  
(15)  
(16)  
(17)  
(18)  
(19) 
References
 (1) Search for resonances decaying to photon pairs in 3.2 fb of collisions at = 13 TeV with the ATLAS detector, Tech. Rep. ATLASCONF2015081, CERN, Geneva, Dec, 2015.
 (2) Search for resonances in diphoton events with the ATLAS detector at = 13 TeV, Tech. Rep. ATLASCONF2016018, CERN, Geneva, Mar, 2016.
 (3) CMS Collaboration Collaboration, Search for new physics in high mass diphoton events in protonproton collisions at TeV, Tech. Rep. CMSPASEXO15004, CERN, Geneva, 2015.
 (4) CMS Collaboration Collaboration, Search for new physics in high mass diphoton events in of protonproton collisions at and combined interpretation of searches at and , Tech. Rep. CMSPASEXO16018, CERN, Geneva, 2016.
 (5) I. Low and J. Lykken, Implications of Gauge Invariance on a Heavy Diphoton Resonance, arXiv:1512.09089.
 (6) J. F. Kamenik, B. R. Safdi, Y. Soreq, and J. Zupan, Comments on the diphoton excess: critical reappraisal of effective field theory interpretations, arXiv:1603.06566.
 (7) M. Backovic, A. Mariotti, and D. Redigolo, Diphoton excess illuminates Dark Matter, JHEP 03 (2016) 157, [arXiv:1512.04917].
 (8) D. Barducci, A. Goudelis, S. Kulkarni, and D. Sengupta, One jet to rule them all: monojet constraints and invisible decays of a 750 GeV diphoton resonance, arXiv:1512.06842.
 (9) Y. Mambrini, G. Arcadi, and A. Djouadi, The LHC diphoton resonance and dark matter, Phys. Lett. B755 (2016) 426–432, [arXiv:1512.04913].
 (10) F. D’Eramo, J. de Vries, and P. Panci, A 750 GeV Portal: LHC Phenomenology and Dark Matter Candidates, JHEP 05 (2016) 089, [arXiv:1601.01571].
 (11) R. Franceschini, G. F. Giudice, J. F. Kamenik, M. McCullough, A. Pomarol, R. Rattazzi, M. Redi, F. Riva, A. Strumia, and R. Torre, What is the resonance at 750 GeV?, JHEP 03 (2016) 144, [arXiv:1512.04933].
 (12) A. Falkowski, O. Slone, and T. Volansky, Phenomenology of a 750 GeV Singlet, JHEP 02 (2016) 152, [arXiv:1512.05777].
 (13) R. S. Gupta, S. JÃ¤ger, Y. Kats, G. Perez, and E. Stamou, Interpreting a 750 GeV Diphoton Resonance, arXiv:1512.05332.
 (14) C. Csaki, J. Hubisz, and J. Terning, Minimal model of a diphoton resonance: Production without gluon couplings, Phys. Rev. D93 (2016), no. 3 035002, [arXiv:1512.05776].
 (15) C. Csaki, J. Hubisz, S. Lombardo, and J. Terning, Gluon vs. Photon Production of a 750 GeV Diphoton Resonance, arXiv:1601.00638.
 (16) L. A. HarlandLang, V. A. Khoze, and M. G. Ryskin, The production of a diphoton resonance via photonphoton fusion, JHEP 03 (2016) 182, [arXiv:1601.07187].
 (17) S. Abel and V. V. Khoze, Photoproduction of a 750 GeV diphoton resonance mediated by KaluzaKlein leptons in the loop, JHEP 05 (2016) 063, [arXiv:1601.07167].
 (18) S. Fichet, G. von Gersdorff, and C. Royon, Scattering light by light at 750 GeV at the LHC, Phys. Rev. D93 (2016), no. 7 075031, [arXiv:1512.05751].
 (19) J. Jaeckel, M. Jankowiak, and M. Spannowsky, LHC probes the hidden sector, Phys. Dark Univ. 2 (2013) 111–117, [arXiv:1212.3620].
 (20) J. Alwall et al., Comparative study of various algorithms for the merging of parton showers and matrix elements in hadronic collisions, Eur. Phys. J. C53 (2008) 473–500, [arXiv:0706.2569].
 (21) ATLAS Collaboration, G. Aad et al., Search for new resonances in and final states in collisions at TeV with the ATLAS detector, Phys. Lett. B738 (2014) 428–447, [arXiv:1407.8150].
 (22) Search for heavy resonances decaying to a boson and a photon in collisions at TeV with the ATLAS detector, Tech. Rep. ATLASCONF2016010, CERN, Geneva, Mar, 2016.
 (23) ATLAS Collaboration, G. Aad et al., Search for an additional, heavy Higgs boson in the decay channel at in collision data with the ATLAS detector, Eur. Phys. J. C76 (2016), no. 1 45, [arXiv:1507.05930].
 (24) Search for diboson resonances in the llqq final state in pp collisions at = 13 TeV with the ATLAS detector, Tech. Rep. ATLASCONF2015071, CERN, Geneva, Dec, 2015.
 (25) ATLAS Collaboration, G. Aad et al., Search for a highmass Higgs boson decaying to a boson pair in collisions at TeV with the ATLAS detector, JHEP 01 (2016) 032, [arXiv:1509.00389].
 (26) Search for a highmass Higgs boson decaying to a pair of W bosons in pp collisions at sqrt(s)=13 TeV with the ATLAS detector, Tech. Rep. ATLASCONF2016021, CERN, Geneva, Apr, 2016.
 (27) CMS Collaboration, S. Chatrchyan et al., Searches for new physics using the invariant mass distribution in pp collisions at =8ââTeV, Phys. Rev. Lett. 111 (2013), no. 21 211804, [arXiv:1309.2030]. [Erratum: Phys. Rev. Lett.112,no.11,119903(2014)].
 (28) Search for heavy particles decaying to pairs of highlyboosted top quarks using leptonplusjets events in proton–proton collisions at TeV with the ATLAS detector, Tech. Rep. ATLASCONF2016014, CERN, Geneva, Mar, 2016.
 (29) ATLAS Collaboration, G. Aad et al., Search for new phenomena in final states with an energetic jet and large missing transverse momentum in pp collisions at 8 TeV with the ATLAS detector, Eur. Phys. J. C75 (2015), no. 7 299, [arXiv:1502.01518]. [Erratum: Eur. Phys. J.C75,no.9,408(2015)].
 (30) ATLAS Collaboration, M. Aaboud et al., Search for new phenomena in final states with an energetic jet and large missing transverse momentum in collisions at TeV using the ATLAS detector, arXiv:1604.07773.
 (31) A. Abdesselam et al., Boosted objects: A Probe of beyond the Standard Model physics, Eur. Phys. J. C71 (2011) 1661, [arXiv:1012.5412].
 (32) A. Altheimer et al., Jet Substructure at the Tevatron and LHC: New results, new tools, new benchmarks, J. Phys. G39 (2012) 063001, [arXiv:1201.0008].
 (33) Search for diboson resonances in the final state in collisions at 13 TeV with the ATLAS detector, Tech. Rep. ATLASCONF2015068, CERN, Geneva, Dec, 2015.
 (34) A. Alloul, N. D. Christensen, C. Degrande, C. Duhr, and B. Fuks, FeynRules 2.0  A complete toolbox for treelevel phenomenology, Comput. Phys. Commun. 185 (2014) 2250–2300, [arXiv:1310.1921].
 (35) J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer, H. S. Shao, T. Stelzer, P. Torrielli, and M. Zaro, The automated computation of treelevel and nexttoleading order differential cross sections, and their matching to parton shower simulations, JHEP 07 (2014) 079, [arXiv:1405.0301].
 (36) R. D. Ball et al., Parton distributions with LHC data, Nucl. Phys. B867 (2013) 244–289, [arXiv:1207.1303].
 (37) S. Knapen, T. Melia, M. Papucci, and K. Zurek, Rays of light from the LHC, Phys. Rev. D93 (2016), no. 7 075020, [arXiv:1512.04928].
 (38) B. Bellazzini, R. Franceschini, F. Sala, and J. Serra, Goldstones in Diphotons, JHEP 04 (2016) 072, [arXiv:1512.05330].
 (39) R. Sato and K. Tobioka, LHC Future Prospects of the 750 GeV Resonance, arXiv:1605.05366.
 (40) J. M. No, Is it or just ? GeV diphoton probes of the electroweak nature of new states, arXiv:1605.05900.
 (41) A. Thamm, R. Torre, and A. Wulzer, Future tests of Higgs compositeness: direct vs indirect, JHEP 07 (2015) 100, [arXiv:1502.01701].
 (42) D. Buttazzo, F. Sala, and A. Tesi, Singletlike Higgs bosons at present and future colliders, JHEP 11 (2015) 158, [arXiv:1505.05488].
 (43) J. W. Sterling. personal communication.
 (44) ATLAS Collaboration, G. Aad et al., Search for a CPodd Higgs boson decaying to Zh in pp collisions at TeV with the ATLAS detector, Phys. Lett. B744 (2015) 163–183, [arXiv:1502.04478].
 (45) Search for a CPodd Higgs boson decaying to Zh in pp collisions at âs = 13 TeV with the ATLAS detector, Tech. Rep. ATLASCONF2016015, CERN, Geneva, Mar, 2016.
 (46) ATLAS Collaboration, G. Aad et al., Search for highmass dilepton resonances in pp collisions at ??TeV with the ATLAS detector, Phys. Rev. D90 (2014), no. 5 052005, [arXiv:1405.4123].
 (47) Search for new phenomena in the dilepton final state using protonproton collisions at â s = 13 TeV with the ATLAS detector, Tech. Rep. ATLASCONF2015070, CERN, Geneva, Dec, 2015.
 (48) ATLAS Collaboration, G. Aad et al., Search for resonant diboson production in the final state in collisions at TeV with the ATLAS detector, Eur. Phys. J. C75 (2015) 69, [arXiv:1409.6190].
 (49) ATLAS Collaboration, G. Aad et al., Search for Higgs boson pair production in the final state from pp collisions at TeVwith the ATLAS detector, Eur. Phys. J. C75 (2015), no. 9 412, [arXiv:1506.00285].
 (50) Search for pair production of Higgs bosons in the final state using protonproton collisions at TeV with the ATLAS detector, Tech. Rep. ATLASCONF2016017, CERN, Geneva, Mar, 2016.
 (51) ATLAS Collaboration, G. Aad et al., Search for highmass diphoton resonances in collisions at TeV with the ATLAS detector, Phys. Rev. D92 (2015), no. 3 032004, [arXiv:1504.05511].