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 formulas in this framework that the production rates of the γγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma \gamma $$\end{document}, gg, Zγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z\gamma $$\end{document}, ZZ, WW∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$WW^*$$\end{document}, tt¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t\bar{t}$$\end{document}, and hh signals at the 750GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$750\, \mathrm{GeV}$$\end{document} resonance are only sensitive to the dilaton–Higgs mixing angle θS\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta _S$$\end{document} and the parameter η≡vNX/f\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta \equiv v N_X/f$$\end{document}, where f is the dilaton decay constant and NX\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_X$$\end{document} 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.


Introduction
About four years ago, the hint of a 125 GeV Higgs boson was reported in the diphoton channel by both the ATLAS and the 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 13-TeV LHC [5,6]. This time the invariant mass of the signal is located around 750 GeV, 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 e-mail: zhujy@itp.ac.cn an ostensible inconsistence in the width of the resonance, 1 both analyses favored the diphoton production rate at about 4 fb 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 100 theoretical papers have appeared to interpret the excess in new physics models , 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 Z Z, 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 → Z Z, 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 the dilaton [165][166][167]. Just like the traditional dilaton theories [168][169][170][171][172][173][174], 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 → Me −φ/ f . As a result, the interactions between the dilaton and the photons/gluons are induced through loop diagrams of these fermions. These features render the MDM 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 Sect. 2, and we present in Sect. 3 some analytical formulas which are used to calculate the diphoton rate. In Sect. 4, 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 Sect. 5 we turn to a discussion of the vacuum stability at high energy scale. Finally, we draw our conclusions in Sect. 6.

The minimal dilaton model
As introduced in the 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 [165,166] 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 λ H S are all free real parameters. About the Lagrangian in Eq.
(1), one should note the 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 [168][169][170][171][172][173][174]. Under these theoretical assumptions, the couplings of the linearized dilaton S to the SM fields take the following form [168][169][170][171][172][173][174]: 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 strengths 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 top-color theory [175], 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 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 = (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 750 GeV diphoton excess by the process gg → s → γ γ . So, in the following, we set m h = 125 GeV, m s = 750 GeV, and v = 246 GeV, 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 the following relations: With the assumption that the dilaton is fully responsible for the fermion masses, the Yukawa coupling of X i is given by 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 it does not depend on y X directly.

Useful formulas in getting the diphoton excess
In the MDM, the particle s may decay into gg, γ γ , Z γ , Z Z, W W * , ff , and hh. In this section, we list the formulas for the widths of these decays, which are needed to get the diphoton rate. As we will show below, these formulas 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 the above expressions, About these formulas, 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 Z Z and hh at the 750 GeV invariant mass (see below) and also by the 125 GeV 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, we have This correlation may serve as a test of the model at future LHC experiments.
If one parameterizes the effective sV V * interaction as then the decay width of s → V V * is given by [158] s→V 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 the 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 terms 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 formulas, the total width of the scalar s and the s-induced diphoton rate can be written as where the new in Eq. (19) 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 the 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 Eqs. (7) to (18) depend on the number of the vector-like fermions N X only through the parameter η ≡ v N X f . As a result, explaining the diphoton excess puts non-trivial requirements on the combination • Since the recent LHC searches for right-handed heavy quarks have required m X 900 GeV [178-180] 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 vectorlike 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 the strong dynamics is usually rich, N X is naturally larger than 1.
We recall 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 as far as we know were not paid attention to in previous literature.

Numerical results and discussions
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 = 1 TeV, 1.5 TeV at each time, and scan the following parameter space: During the scan, we consider the following theoretical and experimental constraints: • The vacuum stability at the scale of m s = 750 GeV for the scalar potential, which corresponds to the requirement 4λ H λ S − λ 2 H S > 0 [165]. • Constraints from the perturbativity at the scale of m s = 750 GeV, which requires λ S , λ H , λ H S 4π , and y X 4π/ √ N c [163]. • Constraints from the electroweak precision data. We calculate the Peskin-Takeuchi S and T parameters [182] with the formulas presented in [165], and construct χ 2 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 γ γ , Z Z * , W W * , bb, and ττ channels [184,185]. We perform the fits like our previous paper [167,186], 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 → Z Z, hh at Run I, just like what we did in [158]. The upper bounds on these channels at 95 % C.L. are listed in Table 1.
For each sample surviving the constraints, we perform a fit to the 750 GeV 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,158] ] 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 σ 13TeV γ γ 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 the following facts: • The central value of the diphoton rate is 3.9fb at the 13-TeV 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.15 GeV if one wants to explain the excess at 2σ level, Table 1 Upper limits on various 750 GeV resonant signals at 8-TeV LHC set by either ATLAS or CMS collaboration [158] Channel j j [187,188] hh [189][190][191][192] W W * [193,194] Z Z [193,195] Z γ [196] tt [197,198] 95 % CL limits 1800 fb 35 fb 37 fb 12 fb 3.6 fb 450 fb Fig. 1 The fit results of the MDM to the 750GeV diphoton data together with the LHC Run I constraints listed in Table 1, 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, Z Z, 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 1 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 Fig. 2 Same samples as those in Fig. 1, but projected on the η-tan θ S planes. Although we take m X = 1 TeV in getting this figure, we check that setting m X = 1.5 TeV produces indistinguishable difference on the figure due to the comments below Eq. (23) while for Q X = 5 3 case, tot 1.6 GeV. 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 1, 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 the 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 1 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. Table 2 Detailed information for one of the best points in the left and right panels of Fig. 2 (labeled by P1 and P2 hereafter), respectively. We checked that all these points predict χ 2 γ γ = 2.32, which corresponds to a p-value of 0.68  . 3 Correlations of the diphoton rate at the 13 TeV with those of Z Z, W W * , hh, and tt signals respectively for the Q X = 2 3 case, which are shown on the η-tan θ S planes. Colors in this figure have same meanings as those in Fig. 2, and from the left to right and upper to lower panels, the constant contours (red lines) of the production rates for Z Z, W W * , hh, and tt signals are shown, respectively. The numbers on the red lines represent the corresponding production rates at the 13-TeV LHC. Note that the correlations of the diphoton rate with those of the gg and Z γ signals are presented in Eq. (23) • 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→Z Z , which are presented from Eqs. (7) to (16).
In Table 2, we show the detailed information for one of the best points in the left and right panels of Fig. 2 respec-tively. 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 13-TeV LHC with the rates of the Z Z, W W * , hh, and tt signals, respectively. The results are presented in Fig. 3 for the Q X = 2 3 case with the implication of the figure • Current LHC data have put upper limits on the rates of the different signals at the 13-TeV LHC, which are σ Z Z 48 fb, σ W W 96 fb, σ hh 190 fb, and σ tt 19 fb. • Since for a moderately small sin θ S , the s Z Z, 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 σ j j are precisely known, one can get the value of Q X , and if σ γ γ and σ Z Z are also measured, one can pin down the favored regions of η and sin θ S .
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 , thus, 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 13-TeV LHC becomes more obvious.

Vacuum stability at high energy scale
About one week before we finished this work, several papers appeared to discuss the vacuum stability in a theoretical framework which is quite similar to the MDM [161][162][163]. 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 the following two reasons. One is that 2 The large Yukawa coupling Y X can influence the vacuum stability condition 4λ H λ S − λ 2 H S > 0 by two ways [163]. 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 m X is proportional to −y 4 X , and consequently λ S usually becomes negative after considering the correction. Table 3 The scale where the vacuum becomes unstable for different choice of the vector-like fermion number N X . Here the scale μ is in unit of GeV, and the points P 1 and P 2 correspond to the two benchmark points in Table 2. 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 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 Sect. 3, the diphoton excess actually imposes non-trivial requirements on the parameter η ≡ v N X f , instead of on the Yukawa coupling y X ≡ ηm X v N 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 2. We repeat the analysis in [163], i.e. we use the same RGEs as those in [163] to run all parameters in the MDM, and also consider the threshold correction to λ S at the scale m X . In Table 3, 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 recall that, although large Q X and/or N X are welcomed to explain the excess, they cannot 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 [163], 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.

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 they may be sizable in comparison with the Hgg 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 formulas that the production rates of the γ γ , gg, Z γ , Z Z, W W * , tt, and hh signals at the 750 GeV resonance are only sensitive to the dilaton-Higgs mixing angle θ S and the parameter η ≡ v N 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 125-GeV Higgs data, the LHC searches for exotic quarks, and the upper bounds on the rates of Z Z, 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 the 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 [159,160]. However, after reading these papers, we learned that the paper [159] considered the traditional dilaton model, and the paper [160] 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,161,163]. Compared with these works, our study has the following features (which we consider 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 Sect. 5, 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 formulas 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 η = v N 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 to the best of 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 1), 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 literature.
Before we end this work, we would like to clarify its relation with our previous work [158], 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 models can predict the central value of the excess in their vast parameter space. Since the two works adopted the same χ 2 function for the excess which only depends on the diphoton rate, the χ 2 values for the best points are the 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 [158], 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 the color index (note that there are totally 32 bosonic freedom), so their coupling strengths are proportional to (C sS A * i S A i v)/m 2 S i A 0 (τ S i ) with C sS A * i S A i 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 strengths 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 [176,177], 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 strengths 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 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 [158] the upper limit of the dijet channel in Table  1 has constrained the diphoton rate to be less than about 7.5 fb [158,199], 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 [158], 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.