Interpreting the 750 GeV diphoton excess in the Minimal Dilaton Model

We try to interpret the 750 GeV diphoton excess in the Minimal Dilaton Model, which extends the SM by adding one linearized dilaton field and vector-like fermions. We first show by analytic formulae in this framework that the production rates of the $\gamma \gamma$, $gg$, $Z\gamma$, $ZZ$, $WW^\ast$, $t\bar{t}$ and $hh$ signals at the $750 {\rm GeV}$ resonance are only sensitive to the dilaton-Higgs mixing angle $\theta_S$ and the parameter $\eta \equiv v N_X/f$, where $f$ is the dilaton decay constant and $N_X$ denotes the number of the fermions. Then we scan the two parameters by considering various theoretical and experimental constraints to find the solutions to the diphoton excess. We conclude that the model can predict the central value of the diphoton rate without conflicting with any constraints. The signatures of our explanation at the LHC Run II and the vacuum stability at high energy scale are also discussed.


I. INTRODUCTION
About four years ago, the hint of a 125 GeV Higgs boson was reported in the diphoton channel by both the ATLAS and CMS collaborations based on about 5 fb −1 data for each collaboration at the 7-TeV LHC [1,2], and this led to the great discovery of the Higgs boson in July 2012 [3,4].
Recently another excess in the diphoton channel was reported by the first 3.2 fb −1 data at the 6]. This time the invariant mass of the signal locates around 750GeV, and its local and global significances are about 3.6σ and 2.3σ respectively for the ATLAS analysis, and 2.6σ and 2σ for the CMS analysis. Interestingly, although there exists a ostensible inconsistence in the width of the resonance 1 , both the analyses favored the diphoton production rate at about 4fb in the narrow width approximation. Such a rate is about 10 4 times larger than the prediction of the Standard Model (SM) with a 750GeV Higgs boson [8]. Obviously, if this excess is confirmed in near future, it points undoubtedly to the existence of new physics.
So far more than one hundred theoretical papers have appeared to interpret the excess in new physics models [9][10][11][12][13][14][15][16][17][18][19], and most of them employed the process gg → S → γγ with S denoting a scalar particle with mass around 750 GeV to fit the data. From these studies, one can infer two essential ingredients of the explanations. One is that there must exist other charged and colored particles to generate by loop effects sufficiently large Sγγ and Sgg interactions. The other is, given the fact that no excess was observed in the channels such as ZZ, W W * and tt at the LHC Run I, the particle S is preferred to be gauge singlet dominated so that the branching ratios of S → ZZ, W W * , tt are not much larger than that of S → γγ. These requirements guide us in seeking for the explanations of the excess.
In this work, we consider interpreting the diphoton excess in the Minimal Dilaton Model (MDM), which extends the SM by one gauge singlet field called dilaton [20][21][22]. Just like the traditional dilaton theories [23], the dilaton in this model arises from a strong interaction theory with approximate scale invariance at a certain high energy scale. The breakdown of the invariance then triggers the electroweak symmetry breaking, and during this process, the dilaton as the pseudo Nambu-Goldstone particle of the broken invariance can be naturally light in comparison with the high energy scale. Furthermore, this model assumes that all SM particles except for the Higgs field do not interact with the dynamics sector, and consequently the dilaton does not couple directly to the fermions and W , Z bosons in the SM. In this sense, the dilaton is equivalent to an electroweak gauge singlet field. The model also consists of massive vector-like fermions acting as the lightest particles in the dynamical sector, to which the dilaton naturally couples in order to recover the scale invariance: M → M e −φ/f . As a result, the interactions between the dilaton and the photons/gluons are induced through loop diagrams of these fermions. These characters enable the MDM as a hopeful theory to explain the diphoton excess through the dilaton production.
Discussing the capability of the MDM in explaining the excess is the aim of this work. This paper is organized as follows. We first introduce briefly the MDM in Section II, and present in Section III some analytical formulae which are used to calculate the diphoton rate. In Section IV, we discuss the constraints on the model, its capability in explaining the excess, and also the related phenomenology at the LHC Run II. For completeness, in section V we turn to discuss the vacuum stability at high energy scale. Finally, we draw our conclusions in Section VI.

II. THE MINIMAL DILATON MODEL
As introduced in last section, the MDM extends the SM by adding one gauge singlet field S, which represents a linearized dilaton field, and also vector-like fermions X i . The low energy effective Lagrangian is then written as [20,21] where L SM is the SM Lagrangian without Higgs potential, f is the decay constant of the dilaton S, M i is the mass of the fermion X i , and N X is the number of the vector-like fermions. The scalar potential V (S,H) contains terms with explicit breaking of the scale invariance, and its general form is given by where m S , λ S , m H , λ H and λ HS are all free real parameters.
About the Lagrangian in Eq.
(1), one should note following points: • The MDM is actually a low energy theoretical framework describing the breakdown of a UV strong dynamics with approximate scale invariance, and the dilaton in this theory is distinguished from the usual one. Explicitly speaking, in the traditional dilaton models the whole SM sector is usually assumed to be a part of the strong dynamics, and all the fermions and gauge bosons of the SM are composite particles at the weak scale [23]. Under these theoretical assumptions, the couplings of the linearized dilaton S to the SM fields take following form [23] L where T µ µ represents the trace of the energy-momentum tensor of the SM. Through the interactions in Eq. (3), the dilaton couples directly to the fermions and W , Z bosons in the SM with the strengthes proportional to the mass of the involved particle. In this way, the dilaton mimics the properties of the SM Higgs boson. By contrast, in the MDM all SM particles except for the Higgs field are assumed to be the spectators of the strong dynamics, and they are all elementary particles. As a result, the dilaton does not couple directly to these particles.
• In the original version of the MDM, the authors set N X = 1 and chose the quantum numbers of the fermion X i same as those of the right-handed top quark. This setting was motivated by topcolor theory [24], which intended to present a reasonable explanation of the relatively large top quark mass within a minimal framework. However, as we will show below, such a setting is tightly limited by the vacuum stability of the theory at m X i scale in interpreting the diphoton excess. Considering that a strong dynamical theory usually involves rich fermion fields and the assignment on their quantum numbers is somewhat arbitrary, we therefore consider a more general but also simple case, which assumes that all the vector-like fermions are identical, and each of them transforms in the (3, 1, Y = 2Q X ) representation of the SM gauge group SU (3) c SU (2) L U (1) Y . In the following, we vary the number of the fermions N X , their common mass m X , and also their electric charge Q X to discuss the diphoton excess.
If one writes the Higgs field in unitary gauge viaH = 1 √ 2 U (0, H) T , the scalar potential in Eq.(2) can be rewritten asṼ In the following, we consider the most general situation in which both H and S take vacuum expectation values (VEV), H = v and S = f , and they mix to form mass eigenstates h and s: In our scheme for the diphoton excess, h corresponds to the 125 GeV Higgs boson discovered at the LHC, and s is responsible for the 750GeV diphoton excess by the process gg → s → γγ. So in the following, we set m h = 125GeV, m s = 750GeV and v = 246GeV, and for the convenience of our discussion, we choose η ≡ v f N X , sin θ S , Q X , N X and m X as the input parameters of the MDM model. In this case, we have following relations With the assumption that the dilaton is fully responsible for the fermion masses, the Yukawa coupling of X i is given by y X ≡ m X f = ηm X vN X . Obviously y X is inversely proportional to N X for fixed η and m X . As we will show below, the diphoton rate is only sensitive to the parameters η, sin θ S and Q X , and does not depend on y X directly.

III. USEFUL FORMULAE IN GETTING THE DIPHOTON EXCESS
In the MDM, the particle s may decay into gg, γγ, Zγ, ZZ, W W * , ff and hh. In this section, we list the formulae for the widths of these decays, which are needed to get the diphoton rate. As we will show below, these formulae are helpful to understand our results.
• The widths of φ → γγ, gg, Zγ with φ = h, s: where the I φ g , I φ γ and I φ Zγ are given by In above expressions, About these formulae, one should note that the terms proportional to cos θ S in the expressions of I s i are contributed by the dilaton component of s, while those proportional to sin θ S come from the H-component of s. One should also note that in the case of sin θ S ∼ 0, which is required by the null excess in the channels such as ZZ and hh at the 750GeV invariant mass (see below) and also by the 125GeV Higgs data, I s γ , I s g and I s Zγ are all dominated by the contribution from the vector-like fermions, and consequently they are correlated. Explicitly speaking, we have I s γ : I s g : correlation may sever as a test of the model at future LHC experiments.
• The widths of the decays s → V V * with V = W, Z.
If one parameterizes the effective sV V * interaction as then the decay width of s → V V * is given by [13] Γ where In the MDM, we have Note that in the expressions of A s Z and B s Z , we have included the one-loop corrections. This is because in case of sin θ S ∼ 0, the corrections are not always smaller than the tree level contributions. Also note that in getting A s Z and B s Z , to a good approximation we have neglected the Z boson mass appeared in the loop functions, and that is why we can express the corrections in term of the simple function A 1 2 (τ X ).
• The width of the tree-level decay s → ff with f denoting any of the fermions in the SM: Note that for this kind of decays, the widths are proportional to sin 2 θ S .
• The width of the tree level decay s → hh: where In getting the final expression of C shh , we have used the relation m 2 s m 2 h and sin θ S ∼ 0 to neglect some unimportant terms. Just like the decays s → W W * and s → tt, Γ s→hh is proportional to sin 2 θ S .
With these formulae, the total width of the scalar s and the s-induced diphoton rate can be written as where the Γ new in Eq. (19) represents the contribution from the exotic decays of s, which may exist if the MDM is embedded in a more complex theoretical framework, Γ SM H→gg denotes the decay width of the SM Higgs H into gg with m H = 750GeV, and σ SM √ s=13TeV (H) = 735fb is the NNLO production rate of the H at the 13 TeV LHC [29]. Obviously, if Γ tot is determined mainly by Γ gg , the rate can be approximated by while if Γ tot takes a fixed value, we have where the normalized cross section σ norm is equal to 0.019 fb (0.018 fb) for m X = 1TeV (1.5TeV).
From the discussion in this section, one can get following important conclusions: • The widths of s → gg, γγ, Zγ or the production rates of the gg, γγ and Zγ signals at the LHC are correlated by Γ s→gg : Γ s→γγ : Γ s→Zγ 1 : • The widths listed from Eq. (7) to Eq. (18) depend on the number of the vector-like fermions N X only through the parameter η ≡ vN X f . As a result, explaining the diphoton excess puts non-trivial requirements on the combination vN X f , instead of on the individual parameter N X or y X = ηm X vN X .
• Since the recent LHC searches for right-handed heavy quarks have required m X 900 GeV [26][27][28] and thus τ X ≡ m 2 s /(4m 2 X ) < 0.2, the loop functions appeared in the widths change slightly with the further increase of m X . This implies that the widths and also the cross section have a very weak dependence on the value of m X . As a result, the results obtained in this work are only sensitive to the parameters η, sin θ S and Q X .
At this stage, one can infer that the parameter N X may also be understood as the total number of the vector-like fermions with the electric charge Q X in the strong dynamics because the contributions of the fermions to the diphoton rate are roughly identical. Since the particle content of a strong dynamics is usually rich, N X is naturally larger than 1.
We remind that the second and third conclusions depend on the assumption that the dilaton is fully responsible for the masses of the vector-like fermions, and within our knowledge, they were not paid attention to in previous literatures. In this section, we discuss the diphoton excess in the MDM. In order to get the favored parameter space for the excess, we fix Q X = 2 3 , 5 3 and m X = 1TeV, 1.5TeV at each time, and scan following parameter space During the scan, we consider following theoretical and experimental constraints: • The vacuum stability at the scale of m s = 750GeV for the scalar potential, which corresponds to the requirement 4λ H λ S − λ 2 HS > 0 [20].
• Constraints from the perturbativity at the scale of m s = 750GeV, which requires λ S , λ H , λ HS 4π, and y X 4π/ √ N c [18].
• Constraints from the electroweak precision data. We calculate the Peskin-Takeuchi S and T parameters [30] with the formulae presented in [20], and construct χ 2 ST by following experi- In our calculation, we require that the samples satisfy χ 2 ST ≤ 6.18.
• Experimental constraints from the 125 GeV Higgs data, which include the updated exclusive signal rates for γγ, ZZ * , W W * , bb and ττ channels [32,33]. We perform the fits like our previous paper [22,34], and require the samples to coincide with the combined data at 2σ level.
• Experimental constraints from the null results in the search for the 750 GeV resonance through other channels such as s → ZZ, hh at Run I, just like what we did in [13]. The upper bounds on these channels at 95% C.L. are listed in Table I.  Table I, which are projected on the σ 13T eV γγ − Γ tot planes for Q X = 2/3 (left panel) and Q X = 5/3 (right panel) respectively. The regions filled by the colors from gray to deep blue represent the parameter spaces that can fit the diphoton data within 3σ, 2σ and 1σ level respectively, and by contrast the regions covered by straw color are excluded by the constraints. The boundaries for the hh, ZZ and W W * channels are also plotted, which correspond to blue lines, red lines and brown lines respectively, and the other constraints listed in Table I are too weak to be drawn on the panels. In each panel, the green line represents the best-fit samples. In getting this panel, we have set Γ new = 0 and m X = 1TeV, and we checked that m X = 1.5TeV predicts roughly same results, which reflects that our results are insensitive to m X .
For each sample surviving the constraints, we perform a fit to the 750GeV diphoton data collected at the 8 TeV and the 13 TeV LHC. In doing this, we use the method introduced in [9], where the data were given by and the χ 2 γγ function was given by [9,13] for the 13 TeV ATLAS data, for the other three sets of data, with µ i denoting the theoretical prediction of the diphoton rate. In the following, we only consider the samples surviving the first four constraints. In Fig.1, we project these samples on the σ 13T eV γγ − Γ tot planes for Q X = 2/3 (left panel) and Q X = 5/3 (right panel) respectively. The details of this figure are explained in its caption. From this figure, one can get following facts: • The central value of the diphoton rate is 3.9fb at the 13TeV LHC from the fit, and the 1σ, 2σ and 3σ ranges of the rate are (2.5 ∼ 5.3) fb, (1.5 ∼ 6.3) fb, (0.2 ∼ 7.9) fb respectively.
Note that this conclusion is independent of the value of Q X .
• For both Q X = 2 3 and Q X = 5 3 cases, the diphoton excess can be well explained. The difference of the two options comes from the fact that for Q X = 2 3 case, Γ tot 0.15GeV if one wants to explain the excess at 2σ level, while for Q X = 5 3 case, Γ tot 1.6GeV. The reason for such a difference is that in the Q X = 5 3 case, sin θ S can take a larger value (see discussion below).
• Among the channels listed in Table I, the hh channel puts the tightest constraints on the parameter space regardless the value of Q X .
Next we illustrate the favored parameter regions for the excess. For this purpose, we project the samples used in Fig.1 on the η − tan θ S planes, which are shown in Fig.2. This figure indicates following facts: • In order to explain the diphoton excess at 2σ level, 0.65 ≤ η ≤ 1.55 and | tan θ S | ≤ 0.012 are preferred for Q X = 2 3 case, and by contrast 0.15 ≤ η ≤ 0.8 and | tan θ S | ≤ 0.06 are preferred for Q X = 5 3 case. Note that in the Q X = 5 3 case, a smaller η as well as a wider range of tan θ S are favored to explain the excess in comparison with the Q X = 2 3 case. The reason is that a larger Q X can increase greatly the width and also the branching ratio of s → γγ, which in return needs a smaller s production rate to explain the excess.
• The channels listed in Table I exclude the parameter space characterized by a large η and/or a large | tan θ S |. For these cases, the production rates of the channels are usually enhanced, which can be inferred from the expressions of the widths.
• In case of tan θ S 0, the Zγ channel may impose upper bounds on η, which is shown in the right panel of Fig.2.
• The favored parameter space is not symmetric if the sign for tan θ S is reversed, and this asymmetry turns out to be more obvious for larger Q X and | tan θ S |. The source of such a asymmetry comes from the expressions of Γ s→gg , Γ s→γγ , Γ s→Zγ and Γ s→ZZ , which are presented from Eq. (7) to Eq. (16).
In Table II, we show the detailed information for one of the best points in the left and right panels of Fig.2 respectively. In the following, we label the two points by P 1 and P 2 respectively.
From this table, one can learn that to explain the diphoton excess in the MDM, the branching ratio of s → γγ is usually at 1% level, which is significantly larger than that of the Higgs boson in the SM. One can also learn that for the best points, s → gg may be either dominant or subdominant decay channel of the s.
Finally, we study the correlations between the diphoton rate at the 13TeV LHC with the rates of the ZZ, W W * , hh and tt signals respectively. The results are presented in Fig.3  • Since for a moderately small sin θ S , the sZZ, sW W , shh and stt couplings are roughly proportional to sin θ S tan θ S , the constant contours of the signal rates exhibit similar behaviors on the η − tan θ S plane. Obviously, if the diphoton excess persists at future LHC experiments and meanwhile none of the other signals is observed, a small tan θ S is preferred.
• More important, if more than one type of the signals are measured at the future LHC experiments, one can decide the parameters of the MDM. For example, given that σ γγ and σ jj are precisely known, one can get the value of Q X , and if σ γγ and σ ZZ are also measured, FIG. 4. Similar to Fig.3, but for the Q X = 5 3 case. Here the scale µ is in unit of GeV, and the points P 1 and P 2 correspond to the two benchmark points in Table II. We checked that for the point P2 with N X = 5, 6, the vacuum keeps stable before λ H reaches its Landau poles, which are roughly at 5.6 × 10 11 GeV and 3.8 × 10 10 GeV respectively. We also checked that for the P 2 with N X = 4, the Landau pole of λ H is roughly at 2.2 × 10 13 GeV.
Point µ(N X = 1) µ(N X = 2) µ(N X = 3) µ(N X = 4) µ(N X = 5) µ(N X = 6) In Fig.4, we show the correlations of the different signals for the Q X = 5 3 case. The features of this figure are quite similar to those of Fig.3 except that: i) now the diphoton rate becomes more sensitive to η and sin θ S , So to extract the values of the two parameters in this case, a more precise measurement of the diphoton signal is needed. ii) the asymmetry between ± tan θ S on the rates at 13TeV LHC becomes more obvious.

V. VACUUM STABILITY AT HIGH ENERGY SCALE
About one week before we finish this work, several papers appeared to discuss the vacuum stability in a theoretical framework which is quite similar to the MDM [16][17][18]. The main argument of these papers was that, in order to explain the diphoton excess, the Yukawa coupling y X must be so large that the vacuum becomes unstable at a certain high energy scale 2 . In our opinion, the MDM may be free of this problem due to following two reasons. One is that the MDM is actually a low energy effective theory describing the breakdown of a strong dynamics with approximate scale invariance. This means that the physics beyond the MDM must appear at a certain high energy scale. The other is that, as we emphasized in Section III, the diphoton excess actually imposes non-trivial requirements on the parameter η ≡ vN X f , instead of on the Yukawa coupling y X ≡ ηm X vN X directly. For a given value of η, one may increase N X to suppress the Yukawa coupling y X , and thus alleviate the problem. In order to verify our speculation, we assume that there are no particles in the strong interaction sector other than the vector-like fermions, and consider the two benchmark points presented in Table II. We repeat the analysis in [18], i.e. we use the same RGEs as those in [18] to run all parameters in the MDM, and also consider the threshold correction to λ S at the scale m X . In Table III, we present the scale where the vacuum becomes unstable for different choices of N X . This table indicates that moderately large N X and Q X are helpful to stabilize the vacuum state.
Finally, we remind that, although large Q X and/or N X are welcomed to explain the excess, they can not be arbitrarily large in the extension of the SM by one gauge singlet scalar and the vector-like fermions. The reason is that the β function of the gauge coupling g 1 is given by β g 1 = ( 41 10 + N X Q 2 X 12 5 )g 3 1 [18], and consequently g 1 increases rapidly with the RGE energy scale for large N X and Q X . In this case, the β function of λ H is dominated by the term proportional to g 4 1 , and consequently, λ H may reach its Landau pole at an energy scale not far above the weak scale.

VI. CONCLUSION
The MDM extends the SM by adding vector-like fermions and one gauge singlet scalar, which represents a linearized dilaton field. In this theory, the couplings of the dilaton to gg and γγ are induced by the loops of the vector-like fermions, and may be sizable in comparison with the Hgg 2 The large Yukawa coupling YX can influence the vacuum stability condition 4λH λS − λ 2 HS > 0 by two ways [18]. One is that it pulls down the value of λS in its evolution with the energy scale by the renormalization group equation (RGE). The other is that the threshold correction to the λS at the scale mX is proportional to −y 4 X , and consequently λS usually becomes negative after considering the correction. and Hγγ couplings in the SM. On the other hand, due to the singlet nature of the dilaton its decays into the other SM particles are suppressed. These characters make the diphoton signal of the dilaton potentially detectable at the LHC.
In this work, we tried to interpret the diphoton excess recently reported by the ATLAS and CMS collaborations at the 13 TeV LHC in the framework of the MDM. For this purpose, we first showed by analytic formulae that the production rates of the γγ, gg, Zγ, ZZ, W W * , tt and hh signals at the 750GeV resonance are only sensitive to the dilaton-Higgs mixing angle θ S and the parameter η ≡ vN X /f , where N X denotes the number of the vector-like fermions and f is the dilaton decay constant. Then we scanned the two parameters to find the solutions to the excess.
During the scan, we considered various theoretical and experimental constraints, which included the vacuum stability and the perturbativity of the theory at the scale of m s , the electroweak precision data, the 125GeV Higgs data, the LHC searches for exotic quarks, and the upper bounds on the rates of ZZ, W W * , Zγ, tt and hh signals at LHC Run I. We concluded that the model can predict the central value of the diphoton rate without conflicting with any constraints. Moreover, after deciding the parameter space for the excess we discussed the signatures of the theory at the LHC Run II. We showed that the rates of the W W * and hh signals may still reach about 100 fb and 200 fb respectively at the 13 TeV LHC, and thus they provide good prospect for detection in future.
As an indispensable part of this work, we also discussed the vacuum stability of the theory at high energy scales. We showed that, by choosing moderately large N X and Q X , the vacuum in our explanation can retain stable up to 10 11 GeV.
Note added: When we finished this work at the beginning of this January, we noted that two papers had appeared trying to explain the diphoton excess with the dilaton field [14,15]. However, after reading these papers, we learned that the paper [14] considered the traditional dilaton model, and the paper [15] focused on 5D warped models. So their studies are quite different from ours.
We also noted that by then there existed several papers studying the diphoton excess in the model which extends the SM by one gauge singlet scalar field and vector-like fermions [10,16,18].
Compared with these works, our study has following features (improvements): • We considered a generic model which predicts N X vector-like fermions (by contrast, most of the previous studies considered the most economical N X = 1 case). This enables us to explain the diphoton excess without invoking a large Yukawa coupling y X . Such a treatment, as we have discussed in section V, is helpful to retain vacuum stability of the theory at high energy scales.
• More important, by assuming that the dilaton field is fully responsible for the masses of the vector-like fermions, we showed by analytic formulae that the rates for all the signals discussed in this work, such as γγ, gg, Zγ, V V * , ff and hh, are only sensitive to the parameter η = vN X f , the dilaton-Higgs mixing angle θ S and the electric charge of the fermions Q X . This observation can greatly simplify the analysis on the diphoton excess, and within our knowledge, it was not paid due attention in previous studies.
• We considered various constraints on the model, especially those from different observations at the LHC Run I (which were listed in Table I), and we concluded that the hh signal usually puts the tightest constraint on our explanation. This conclusion is rather new. Moreover, we also studied the signatures of our explanation at the LHC Run II, which are helpful to decide the parameters of the model. Such a study was absent in previous literatures.
Before we end this work, we'd like to clarify its relation with our previous work [13], where we utilized the singlet extension of the Manohar-Wise model to explain the diphoton excess. In either of the works, the scalar sector of the considered model contains a doublet and a singlet scalar field, which mix to form a 125 GeV SM-like Higgs h and a 750 GeV new scalar s, and the sγγ and sgg interactions are induced by colored particles through loop effects. In organizing these works, we first introduced the theoretical framework and listed the formula for the partial widths of the scalar s, then we analyzed various constraints on the model and discussed the diphoton signal from the process gg → s → γγ. We concluded that both the models can predict the central value of the excess in their vast parameter space. Since the two works adopted same χ 2 function for the excess which only depends on the diphoton rate, the χ 2 values for the best points are same in the two explanations. In spite of these similarities, we still think that the two works are independent since they are based on different physics. The differences are reflected in following aspects: • The origin of the singlet dominated scalar s. In the work [13], the singlet field is imposed by hand and only for interpreting the excess, while in this work it corresponds to a linearized dilation field, which is well motivated by the broken of a strong dynamic with approximate scale invariance.
• The mechanism to generate sizable sγγ and sgg interactions. In the singlet extension of the Manohar-Wise model, these interactions are induced by color-octet and isospin-doublet scalars S A R , S A I and S A ± with A = 1, · · · 8 denoting color index (Note that there are totally denoting the coupling coefficient for the sS A * i S A i interaction. As a comparison, the couplings in this work are induced by the vector-like fermions, and their strengthes are determined by the factor ηA 1 2 (τ X ). Since the loop function A 0 is usually several times smaller than the function A 1 2 [25], beside the large bosonic freedom, large C sS A * i S A i and meanwhile moderately light S A i are also necessary to get the same sizes of the strengthes as those in this work. By contrast, we only need to tune the value η to get the right couplings for the excess in this work. So the explanation presented in here is rather simple and straightforward.
• The intrinsic features of the explanations. Due to the particle assignments of the models, the two explanations exhibit different features. For example, for the explanation in [13] the upper limit of the dijet channel in Table I has constrained the diphoton rate to be less than about 7.5fb [13,47], while in the present work the constraint from the dijet channel on the rate is rather loose. Another example is that for the explanation in [13], the vacuum stability can never constrain the model parameters, while in this work it acts as a main motivation to consider moderately large N X and Q X to keep the vacuum stability. ACKNOWLEDGEMENT We thank Prof. C. P. Yuan and Fei Wang for helpful discussion, and this work was supported in part by the National Natural Science Foundation of China (NNSFC) under Grant No. 11547103, 11275245, 11547310, 11575053. Dr. Zhu thanks the support of the U.S. National Science Foundation under Grant No. PHY-0855561, while he was working at Michigan State University.