Abstract
Recently, both ATLAS and CMS collaborations at the CERN Large Hadron Collider (LHC) have announced their observations of an excess of diphoton events around the invariant mass of with a local significance of and , respectively. In this paper, we interpret the diphoton excess as the onshell production of a real singlet scalar in the channel. To accommodate the observed production rate, we further introduce a vectorlike fermion , which is carrying both color and electric charges. The viable regions of model parameters are explored for this simple extension of the Standard Model (SM). Moreover, we revisit the problem of electroweak vacuum stability in the same scenario, and find that the requirement for the electroweak vacuum stability up to high energy scales imposes serious constraints on the Yukawa coupling of the vectorlike fermion and the quartic couplings of the SM Higgs boson and the new singlet scalar. Consequently, a successful explanation for the diphoton excess and the absolute stability of electroweak vacuum cannot be achieved simultaneously in this economical setup.
Electroweak Vacuum Stability and Diphoton Excess at 750 GeV
Jue Zhang ^{*}^{*}*Email: zhangjue@ihep.ac.cn Shun Zhou ^{†}^{†}†Email: zhoush@ihep.ac.cn
Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China
Center for High Energy Physics, Peking University, Beijing 100871, China
PACS number(s): 12.38.Bx, 12.60.Fr, 14.80.Bn
1 Motivation
Exactly four years ago, the ATLAS [1] and CMS [2] collaborations announced publicly the experimental indications of a scalar particle in the narrow mass range of both in the diphoton channel and in the fourlepton mode . Only half a year later, with higher integrated luminosities and more accumulated data, both collaborations [3, 4] confirmed the existence of a scalar particle of mass , which is very likely to be the Higgs boson in the SM [5]. The Higgs boson is playing a particularly important role in the SM [6], as it is responsible for the spontaneous gauge symmetry breaking of and for generating the masses of weak gauge bosons, quarks and charged leptons. However, the discovery of the SMlike Higgs boson is not the end of story, since there are several unsolved fundamental problems that are associated with such a scalar particle.
First of all, the radiative corrections to the Higgsboson mass are not protected by any symmetries, so huge corrections from new physics at a highenergy scale, e.g., the scale of Grand Unified Theories (GUT) or the Planck scale , should be cancelled out by an equally large counter term [7]. There already exist possible solutions to this finetuning or naturalness problem, such as the ColemanWeinberg mechanism [8], which attempts to remove all the dimensionful parameters and restore a classically scaleinvariant theory, or supersymmetries [9], which guarantee an elegant cancellation between the radiative corrections from fermions and those from bosons.
Second, the scalar potential of the SM may develop another much deeper minimum at a large field value, destabilizing the electroweak vacuum [10, 11]. This indeed happens when the Higgs quartic coupling becomes negative at a highenergy scale. In particular, for the Higgsboson mass , the renormalizationgroup (RG) evolution of is significantly affected by a large and negative contribution from the topquark Yukawa coupling . In Ref. [12], it has been found that new physics should show up below a cutoff scale of to stabilize the electroweak vacuum, for the Higgsboson mass , the topquark mass and the strong finestructure constant at the scale of the weak boson mass . The problem of vacuum stability has also been investigated by many other authors [13], and the improvements have been made by including highorder corrections [14]. The final verdict is that the absolute stability of the SM vacuum up to the Planck scale is excluded at the confidence level for [14].
In addition, most free parameters in the SM arise from the Yukawa couplings between the Higgs boson and the fermions. The observed strong quark mass hierarchy, flavor mixing angles and CP violation remain unexplained and have been longstanding flavor puzzles [15]. The solutions to those fundamental problems, together with the origin of neutrino masses [16], definitely call for new physics beyond the SM.
Very recently, based on the data collected from the LHC Run II at a centerofmass energy , the ATLAS [17] and CMS [18] collaborations released their observations of a diphoton excess around the invariant mass of with a local significance of and , respectively. If these observations are confirmed by the future data, we can claim the first and direct discovery of new physics beyond the SM. Therefore, it is intriguing to investigate possible explanations for the diphoton excess, as done in a torrent of papers [19]. Different from most of the previous works, we focus on the problem of electroweak vacuum stability, and examine whether the new physics indicated by the diphoton exceess could help stabilize the electroweak vacuum. Although only a minimal extension of the SM with a real singlet scalar and a vectorlike fermion is considered for illustration, the generalization of our analysis to other models proposed in Ref. [19] is straightforward.
The rest of our paper is structured as follows. In Sec. 2, the minimal model is introduced to explain the diphoton excess as a scalar resonance in , where the vectorlike fermion is implemented to mediate the and processes at the oneloop level. The diphoton signal rate and other collider constraints are used to explore the viable parameter space. In Sec. 3, we derive the conditions for the absolute vacuum stability, and present the RG equations (RGEs) for the SM gauge couplings, the Yukawa coupling of the vectorlike fermion and quartic scalar couplings. Then, we have carried out a complete numerical study of model parameters and found that a successful explanation for the diphoton excess at and the vacuum stability up to the Planck scale cannot be achieved simultaneously in this minimal model. Finally, we summarize in Sec. 4.
2 The Diphoton Excess
Motivated by the diphoton excess, we extend the SM with a real singlet scalar and a vectorlike fermion . Under the SM gauge group is singlet while transforms as , where the color, weak isospin and hypercharge assignments are indicated in the conventional form, e.g., or stands for the dimension of the corresponding representation. The general gaugeinvariant Lagrangian can be written as
(1) 
where the covariant derivative is defined as with and being proper representations of the generators of and gauge groups, respectively. The electric charge of is then given by the GellMannNishijima formula .
In addition, the treelevel scalar potential can be parametrized as follows
(2) 
where the vacuum expectation values (vev’s) of the SM Higgs and new singlet scalar fields are and . Note that we have imposed a symmetry on the Lagrangian in Eq. (1), under which all the SM fields are even, while the other fields transform as and and . After the spontaneous symmetry breaking, the vectorlike fermion acquires its mass . Some comments on the model are in order. The spontaneous breaking of the symmetry may lead to the domainwall problem, which however can be solved by adding an explicit symmetrybreaking term into the Lagrangian. For instance, a Dirac mass term suffices for this purpose. In this case, the fermion mass is given by , so we may just take both and as free model parameters. This is also true for a more general theory in which the fermion mass is not subject to spontaneous symmetry breaking but some underlying strong dynamics. In the present work, we take the scalar potential in Eq. (2) as a phenomenological example, which could be realized in various theoretical models with a Higgs portal to new physics [21, 22].
2.1 Signals and Constraints
After specifying the theoretical framework, we are now in a good position to explain the diphoton excess at and examine the existing constraints on the minimal model presented in Eqs. (1) and (2). The signal cross section can be estimated as for ATLAS and for CMS at the centerofmass energy of , which are compatible at the level with the nonobservation of diphoton excess in the LHC Run I data at [23, 24]. Apart from the hint at an invariant mass of , the ATLAS collaboration also suggests a bestfit width of about or equivalently , while the data from CMS favor mostly a narrow width. But a wide resonance is also consistent with the CMS data at the level [25].
Now it becomes evident that the diphoton excess can be interpreted as the channel production of the singlet scalar of mass . Naively, we expect that the production is dominated by the gluongluon fusion and the decay is via . As demonstrated in Ref. [26], just one singlet scalar coupled to the SM particles cannot explain this diphoton excess without violating existing constraints from other channels. It is therefore necessary to add the vectorlike fermion , which meditates both production and decay processes through oneloop diagrams. In order to avoid tight constraints from the diboson searches at the LHC, we take to be an singlet, and its sole contributions to the partial decay widths of and are found as follows [27]
(3)  
(4) 
where the fine structure constants and are both evaluated at . In addition, we have defined , and the function is
(5) 
with given by
(6) 
The partial decay widths to and via the vectorlike fermions are found to correlate with that to , however, they are respectively suppressed by and , where is the Weinberg angle.
Through the interactions with the SM Higgs field, this scalar field acquires additional decay modes. First, according to the scalar potential in Eq. (2), a mixing between and the SM Higgs boson occurs after the spontaneous symmetry breaking. Because of the mixing with , the singlet scalar is coupled to the SM fermions and gauge bosons, although the coupling strengths are suppressed by the mixing angle .^{1}^{1}1As we will show in next section, can be expressed in terms of the model parameters. For the moment, we keep it as general as possible and investigate its viable values. Since a SMlike Higgs boson with a mass of predominantly decays into , and [28], we thus expect contributions from the mixing with the SM Higgs field to the partial decay widths of and . Their contributions can be easily found by multiplying the corresponding partial decay widths of a SMlike Higgs boson of mass with a suppression factor of .
Apart from the above mixing effect, the scalar potential in Eq. (2) also induces the decay of into two SM Higgs bosons through the quartic coupling , when picks up a nonzero vev . A precise study of this decay mode involves changing the basis to the mass eigenstates of both and , i.e., and , and carefully taking into account the singletdoublet mixing effect. See, e.g., Refs. [29, 30], for more details. However, given the fact that the mixing angle tends to be small in order to satisfy the experimental constraints, we can write the trilinear interaction as at leading order, where the coupling is approximately given by . From the later discussion on the diagalization of the scalar mass matrix, we may further write in terms of the mass of the scalar field and the mixing angle when , namely, . Then, according to Ref. [29], we find the partial decay width of as follows
(7) 
where is the mass of the SM Higgs boson.
To summarize, due to the interactions with the vectorlike fermion and the SM Higgs boson , the scalar boson decays in several modes, e.g., and . It should be noted that among these decay modes, and contain interference effects between the contributions from the vectorlike fermion and that from the SM, which have been taken into account in our calculations. Moreover, we also include the nexttoleading order factor of for the decay mode.
Given the partial decay widths of , we are now ready to calculate its diphoton production rate at the LHC. For the resonant production, the signal cross section can be expressed in terms of relevant decay widths corresponding to the production and decay processes. At , the diphoton signal rate can be recast into the form [25]
(8) 
where is the total decay width of , and is the decay width corresponding to the dominant production mechanism of , namely, .
In addition to the condition in Eq. (8) that accounts for the correct diphoton production rate, we should also ensure that the production rates of final states other than two photons are not in contradiction with current experimental observations. To this end, we quote the results from Ref. [25], in which the present constraints from direct searches for have been reinterpreted as ratios of partial decay widths to . Those constraints are shown in Table 1, where the production rates at are scaled from those at by a factor of , and the decay width fits to the central value of the observed diphoton excess.
Final States  Bounds on [25] 

2.2 Viable Parameter Space
Next, using the signal rate in Eq. (8) and the constraints in Table 1, we explore the allowed parameter space of our model. The details of our numerical calculations can be summarized as follows:

For simplicity, we assume that the vectorlike fermion forms a color triplet, belonging to the fundamental representation of the gauge group, i.e., . The generalization to a different color multiplet is straightforward.

We focus on the parameter space of the mixing angle and the Yukawa coupling , which are of crucial importance for our later discussions on the electroweak vacuum stability in the next section. For the values of , we choose a constrained region of , where the upper bound arises from current measurements of Higgs couplings [31]. As for the Yukawa coupling , the range is chosen as , for which the upper bound is required by the perturbativity.

Finally, for a given set of and , we solve Eq. (8) for the fermion mass , and further restrict it into the range of . Here we consider a lower bound of so as to employ the threshold effect presented in to enhance the signal production rate. Vectorlike fermions with such a light mass are actually already excluded by 8 TeV LHC results [32], if assuming that they can decay to the SM top or bottom quark via , or . However, in our scenario the vectorlike fermion has no such decay modes, so the bounds from current experimental searches are evaded. For the upper bound of we choose it for the testability in future collider experiments.
Following the above strategy, we obtain the allowed parameter space for and in Fig. 1, where two possible electric charge assignments of and are assumed for illustration.^{2}^{2}2If the electric charge of new vectorlike fermions is much larger than , the electromagnetic gauge coupling will run quickly into a nonperturbative region at the TeV scale. The requirement for perturbativity has been taken into account in the next section, and also discussed by others [33]. The gray shaded regions indicate the allowed parameter space when only the diphoton production rate condition of Eq. (8) is considered, while further imposing the constraints from other decay modes leads to the colored regions. Notice that vectorlike fermions with exotic electric charges may appear in the composite models [25], and the studies of the scenarios with different electric charges can be performed in a similar way. Two comments on the numerical results are in order:

As one can observed from Fig. 1, to account for the observed diphoton rate, we need a rather large value of and a small mixing angle . This observation is understandable since a large signal rate needs a big Yukawa coupling and a sizable electric charge. Therefore, if a smaller electric charge is taken, one has to increase the Yukawa coupling or multiply the generations of vectorlike fermions in order to explain the diphoton excess.

The color bars in Fig. 1 indicate the total decay width of is in general narrow. In both cases the widths are only on the order of subGeV, which is much narrower than the observed . However, for the time being, this is not a big issue, since a narrow width cannot be excluded. On the theoretical side, this is reasonable for a weaklycoupled theory, whereas a wide resonance is naturally expected in a stronglycoupled theory, such as the composite models.
Finally, it is worthwhile to note that the fermion mass and the Yukawa coupling can be closely related, as we mentioned before. In this case, the model parameters will be more severely constrained.
3 Electroweak Vacuum Stability
If the diphoton excess is further confirmed by future data, one will be convinced that new physics should be introduced to account for the observations, as we have shown. If the minimal scenario of a singlet scalar and a vectorlike fermion works well for this purpose, then an interesting question is whether they help stabilize the electroweak vacuum. The same question can also be asked in any other new physics scenarios. Unfortunately, as we demonstrate below, it is impossible to simultaneously explain the diphoton excess and solve the vacuum stability problem in this minimal scenario. Even worse, the new fermion with a sizable Yukawa coupling will destabilize the electroweak vacuum.
3.1 Conditions for Absolute Stability
Since the SM Higgs boson of mass is relatively light, while the top quark is heavy, the quartic Higgs coupling may run into a region of negative values at a high renormalization scale. As a consequence, the effective potential becomes unbounded from below, or a global minimum that is much deeper than the electroweak vacuum develops at large field values. For , and , it has been found that a negative value of is reached for [12].
In Ref. [20], a simple but intriguing solution to the instability problem of the electroweak vacuum has been proposed. The essential idea is to introduce an extremely heavy scalar of mass , which leads to a threshold effect on the quartic coupling of Higgs boson. From the potential in Eq. (2), one can immediately derive the mass matrix for the scalars
(9) 
from which one can further obtain the scalar masses
(10) 
In Ref. [20], both and are assumed to be much larger than the electroweak scale, so the mixing between singlet and doublet scalar bosons is extremely small (i.e., ) and the scalar masses turn out to be and . Therefore, even for , there exists a constant shift in the quartic Higgs coupling . It is this threshold effect that has been implemented in Ref. [20] to eliminate the instability problem.
Different from the authors of Ref. [20], we fix the scalar masses at those observed from the ATLAS and CMS experiments, namely, and . In this case, the mixing between doublet and singlet scalars can be significant
(11) 
depending on the vev’s and quartic scalar couplings. As shown in the previous section, the mixing angle will be severely constrained from collider experiments.
In order to find out the conditions for the absolute stability, we recast the scalar potential into the following form
(12) 
from which one can see that the condition for requires , and . In the following discussions, we implement the latter three inequalities as the criteria for an absolute electroweak vacuum. Some remarks are helpful. First, the parametrization of ensures that the SM vacuum corresponds to , so the requirement of at large values of scalar fields guarantees that there is no minimum deeper than the SM vacuum. Second, to reliably compute the scalar potential at a large field value, one has to evaluate the quartic couplings at the corresponding energy scale. Therefore, the RGEs for the gauge, Yukawa and quartic couplings should be derived. This is the main task for the next subsection.
3.2 RenormalizationGroup Equations
In this subsection, we derive the oneloop RGEs for the dimensionless couplings in the extension of the SM with a singlet real or complex scalar and a general vectorlike fermion, whose quantum numbers under the SM gauge group can be arbitrary. In Refs. [34, 35, 36, 37], the oneloop RGEs for an additional scalar have been obtained, while no additional vectorlike fermions are considered. A similar scenario to ours has been studied in Ref. [38], where the authors investigated a specific vectorlike fermion model and provided relevant RGEs for the real scalar case. Here we extend the RGE study in Ref. [38] in two aspects. First, both real and complex scalars are examined. Second, the most general assignment of SM gauge quantum numbers to the vectorlike fermions is studied, so that other vectorlike fermion models can be covered as well. In practice, we adopt the general oneloop RGE formulas given in Ref. [39].
Let us start with the real scalar case. First of all, the RGEs of three gauge couplings get modified because of the additional vectorlike fermions that also are carrying SM gauge quantum numbers. Above the mass threshold of these vectorlike fermions, we find^{3}^{3}3A GUT normalization of is used.
(13) 
with the coefficients of given by
where are coefficients in the SM, and denotes the quadratic Dynkin index of the representation of the group . Here we take the convention that for the adjoint representation, while for the fundamental representation, where is the conjugate representation of .
Next, we turn to the Yukawa couplings. Since the SM Higgs field does not couple to the extra vectorlike fermions, as well as no additional Yukawa couplings between SM fermions and the extra scalar field, the RGEs of the SM Yukawa couplings stay the same. For the newly introduced Yukawa coupling , we obtain its RGE as
(14) 
where is the quadratic Casimir invariants for the representation of the group . For example, we have and .
Finally, we give the RGEs of the quartic couplings , and . Neglecting all the SM Yukawa couplings except for that of the top quark , we find
(15) 
As one can see, the RGEs are coupled together, and the running of can now be remarkably affected by . Moreover, the additional Yukawa coupling appears only in the RGEs of and , and it can play an important role in the vacuum stability study, as could be driven to be negative at some high energy scale if were too large.
The complex scalar case can be studied in a similar way, and the corresponding RGEs differ from those in Eqs. (14) and (3.2) only in and three quartic couplings. The Lagrangian takes the same form as in Eq. (1) but with a distinct scalar potential
(16) 
In this case, we find the RGE of as follows
(17) 
while the RGEs of quartic couplings are given by
(18) 
The applications of the RGEs in the case of complex scalar can be made similarly as we shall do in the next subsection for the real case. This indeed makes sense, as recently one also considers the explanation for the diphoton excess with a complex scalar, e.g., see Refs. [40].
3.3 Viable Parameter Space
It is time to examine if the electroweak vacuum can be stabilized. To this end, we explore the parameter space for which the electroweak vacuum can be stable up to a highenergy scale. For illustration, we choose several typical energy scales, such as the instability scale of the SM vacuum, the GUT scale and the Planck scale . The strategy of our numerical analysis is summarized below

Following Ref. [12], we first run the SM couplings from to , by using twoloop SM RGEs and adopting oneloop matching for the Higgs quartic coupling at the topquark mass threshold. In our calculations, the pole mass of top quark is used. As emphasized in Ref. [41], the issue of electroweak vacuum stability is highly sensitive to the topquark mass, whose experimental uncertainty remains large. This uncertainty should be considered in a more precise study, as the running behavior of will also be changed for a different value of in our case.

Then, at , we perform the treelevel matching for the scalar quartic couplings , and . Such a matching is achieved by solving scalar quartic couplings from Eq. (3.1), given the values of and , and the requirements of and .

Finally, for a given value of , we are able to run all parameters in our model to high energy scales to inspect the vacuum stability problem, with oneloop RGEs derived the previous subsection. It is worthwhile to stress that a dedicated study of the vacuum stabiblity should be carried out with twoloop RGEs and oneloop matching conditions, particularly if the diphoton excess is further confirmed. Learning from the SM case [14], we expect that the vacuum stability should be sensitive to highorder radiative corrections. At present, however, the oneloop calculations may be adequate for us to understand how the singlet scalar and vectorlike fermions significantly modify the stability of the electroweak vacuum.
Requiring the vacuum stability conditions below Eq. (12) to be satisfied till some highenergy scale and all couplings to stay within the perturbative regions (namely, , , , and ), we obtain the allowed parameter space for and in Fig. 2. In consideration of the diphoton production rate, we assume the vectorlike fermion to be a color triplet as before. In Fig. 2, we take two different values of , i.e., and , and also assume and as in Fig. 1. Three different highenergy scales, at which a deeper electroweak vacuum would be developed, are discussed. The allowed regions for these three cases are indicated by (dark purple) for the Planck scale, (dark purple + purple) for the GUT scale, and lastly (dark purple + purple + pink) for the SM instability scale. We also show the regions that only satisfy the perturbativity requirements for the above three scales, and those favored by the diphoton excess (gray shaded regions). Several observations from Fig. 2 are in order.

For , we find that the perturbativity and stability requirements till the Planck scale cannot be satisfied at all for the whole parameter space. However, if relaxing such a breakingdown scale to the GUT scale, some allowed parameter space emerges. In contrast, in the case of we do find the parameter space allowed by both perturbativity and stability till the Planck scale.

Comparing with the regions favored by the diphoton excess, one is able to rule out the possibility of with the perturbtivity requirements alone. However, in the case of a large portion of parameter space favored by the diphoton excess is still retained by perturbativity. Thus, requiring perturbativity alone is impossible to exclude all possibilities of explaining the diphoton excess.

Finally, as one can observe from Fig. 2, the further requirement for vacuum stability sets more stringent limits on the parameter space. In all four cases, the regions of and are favored by perturbativity and vacuum stability. This allowed parameter space in fact can be understood by adopting the boundary matching conditions Eqs. (3.1) and (11), and the running behaviors of ’s. To see this, let us work in the case where . From the matching conditions, we then obtain
(19) Therefore, if were too small, would receive little treelevel threshold effects, and the positive RG running contribution from is also negligible, so the vacuum stability problem remains as serious as that in the SM. On the other hand, if tends to be too large, say, , this threshold contribution would give at low energy boundary scale. Such a large boundary value of would easily cause to run towards the Landau pole. As for the favored region of , it is mainly due to the fact that a large value of would drive to be negative, exactly in the same way as too large topquark Yukawa coupling would lead to a negative value of at highenergy scales.

Lastly and most importantly, we observe that even in the scenario of there is no overlap between the regions allowed by both stability and perturbativity and those favored by the diphoton excess. This indicates that within this simple model one is not able to simultaneously explain the recently observed diphoton excess and cure the electroweak vacuum stability problem. Therefore, additional degrees of freedom are certainly needed. They could be extra vectorlike fermions as discussed in [42], because more vectorlike fermions allow for smaller values of to account for the diphoton excess. Alternatively, one may introduce other scalar fields, which may give additional positive contributions to the quartic coupling of the SM Higgs field [43]. However, the investigations of these possibilities are beyond the scope of current paper.
It is also important to point out that the requirement of the stability of the electroweak vacuum could inversely place restrictive bounds on the new physics models that are intended for the diphoton excess [19].
4 Summary
The recent observations of diphoton excess from both ATLAS and CMS collaborations have inspired a great number of studies on possible interpretations from new physics beyond the SM. Among them, one simple and straightforward explanation is to identify the diphoton excess as the onshell production of a singlet scalar in the channel, as well as a pair of vectorlike fermions participating in both production and decay processes. Considering the possibility that this singlet scalar field may have nonnegligible interactions with the SM Higgs fields, we set out to study the question whether the SM electroweak vacuum instability problem can be cured by this newly introduced scalar field and its associated vectorlike fermions. Unfortunately, we find that simultaneously explaining the observed diphoton excess and solving the vacuum stability problem is not possible within this simple scenario. Moreover, we have identified the reason behind it, namely, accounting for the diphoton production rates in general requires a large value of the Yukawa coupling between the scalar field and the vectorlike fermions, however, such a large value of Yukawa coupling would in turn drive the quartic scalar coupling to be negative, at a scale that is even lower than that present in the SM.
In this respect, among the existing proposals, which explain the diphoton excess but also yield a relatively large value of Yukawa couplings between the introduced scalar fields and fermions, the issue of vacuum instability can be their Achilles’ heel. A simultaneous explanation of both diphoton excess and this vacuum instability issue may call for additional complexities from new physics beyond the SM. We are looking forward to deeper thoughts and more experimental data to sort them out.^{4}^{4}4Recently, both ATLAS [44] and CMS [45] collaborations have updated their analyses of the diphoton data, and found that the signal at 750 GeV becomes slightly stronger. Although this update does not change our results, it is now more encouraging to study new physics scenarios explaining the diphoton excess and take into account both vacuum stability and perturbativity of relevant coupling constants.
Acknowledgements
This work was supported in part by the Innovation Program of the Institute of High Energy Physics under Grant No. Y4515570U1, by the National Youth Thousand Talents Program, and by the CAS Center for Excellence in Particle Physics (CCEPP).
References
 [1] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 710, 49 (2012) [arXiv:1202.1408].
 [2] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 710, 26 (2012) [arXiv:1202.1488].
 [3] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716, 1 (2012) [arXiv:1207.7214].
 [4] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716, 30 (2012) [arXiv:1207.7235].
 [5] S. L. Glashow, Nucl. Phys. 22, 579 (1961); S. Weinberg, Phys. Rev. Lett. 19, 1264 (1967); A. Salam, in Elementary Particle Theory, edited by N. Svartholm (Almqvist and Wiksells, Stockholm, 1968), p. 367; Weak and Electromagnetic Interactions, Conf. Proc. C 680519, 367 (1968).
 [6] P. W. Higgs, Phys. Lett. 12, 132 (1964); F. Englert and R. Brout, Phys. Rev. Lett. 13, 321 (1964); P. W. Higgs, Phys. Rev. Lett. 13, 508 (1964); G. S. Guralnik, C. R. Hagen and T. W. B. Kibble, Phys. Rev. Lett. 13, 585 (1964); P. W. Higgs, Phys. Rev. 145, 1156 (1966); T. W. B. Kibble, Phys. Rev. 155, 1554 (1967).
 [7] G. ’t Hooft, “Naturalness, chiral symmetry, and spontaneous chiral symmetry breaking,” NATO Sci. Ser. B 59, 135 (1980).
 [8] S. R. Coleman and E. J. Weinberg, Phys. Rev. D 7, 1888 (1973).
 [9] H. P. Nilles, Phys. Rept. 110, 1 (1984); H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985); M. F. Sohnius, Phys. Rept. 128, 39 (1985).
 [10] N. Cabibbo, L. Maiani, G. Parisi and R. Petronzio, Nucl. Phys. B 158, 295 (1979); P. Q. Hung, Phys. Rev. Lett. 42, 873 (1979); M. Lindner, Z. Phys. C 31, 295 (1986); M. Lindner, M. Sher and H. W. Zaglauer, Phys. Lett. B 228, 139 (1989); P. B. Arnold, Phys. Rev. D 40, 613 (1989); M. Sher, Phys. Rept. 179, 273 (1989); Phys. Lett. B 317, 159 (1993) [hepph/9307342]; B. Schrempp and M. Wimmer, Prog. Part. Nucl. Phys. 37, 1 (1996) [hepph/9606386].
 [11] G. Altarelli and G. Isidori, Phys. Lett. B 337, 141 (1994); J. A. Casas, J. R. Espinosa and M. Quiros, Phys. Lett. B 342, 171 (1995) [hepph/9409458]; Phys. Lett. B 382, 374 (1996) [hepph/9603227]; T. Hambye and K. Riesselmann, Phys. Rev. D 55, 7255 (1997) [hepph/9610272]; G. Isidori, G. Ridolfi and A. Strumia, Nucl. Phys. B 609, 387 (2001) [hepph/0104016]; G. Isidori, V. S. Rychkov, A. Strumia and N. Tetradis, Phys. Rev. D 77, 025034 (2008) [arXiv:0712.0242]; J. Ellis, J. R. Espinosa, G. F. Giudice, A. Hoecker and A. Riotto, Phys. Lett. B 679, 369 (2009) [arXiv:0906.0954].
 [12] Z. Z. Xing, H. Zhang and S. Zhou, Phys. Rev. D 86, 013013 (2012) [arXiv:1112.3112].
 [13] M. Holthausen, K. S. Lim and M. Lindner, JHEP 1202, 037 (2012) [arXiv:1112.2415]; J. EliasMiro, J. R. Espinosa, G. F. Giudice, G. Isidori, A. Riotto and A. Strumia, Phys. Lett. B 709, 222 (2012) [arXiv:1112.3022];
 [14] G. Degrassi, S. Di Vita, J. EliasMiro, J. R. Espinosa, G. F. Giudice, G. Isidori and A. Strumia, JHEP 1208, 098 (2012) [arXiv:1205.6497]; F. Bezrukov, M. Y. Kalmykov, B. A. Kniehl and M. Shaposhnikov, JHEP 1210, 140 (2012) [arXiv:1205.2893]; S. Alekhin, A. Djouadi and S. Moch, Phys. Lett. B 716, 214 (2012) [arXiv:1207.0980]; I. Masina, Phys. Rev. D 87, no. 5, 053001 (2013) [arXiv:1209.0393]; K. G. Chetyrkin and M. F. Zoller, JHEP 1304, 091 (2013) [arXiv:1303.2890]; D. Buttazzo, G. Degrassi, P. P. Giardino, G. F. Giudice, F. Sala, A. Salvio and A. Strumia, JHEP 1312, 089 (2013) [arXiv:1307.3536]; V. Branchina and E. Messina, Phys. Rev. Lett. 111, 241801 (2013) [arXiv:1307.5193].
 [15] Z. Z. Xing, Int. J. Mod. Phys. A 29, 1430067 (2014) [arXiv:1411.2713].
 [16] Z. Z. Xing and S. Zhou, “Neutrinos in Particle Physics, Astronomy and Cosmology,” SpringerVerlag, Berlin Heidelberg (2011); S. Zhou, arXiv:1511.07255.
 [17] The ATLAS Collaboration, “Search for resonances decaying to photon pairs in 3.2 fb of collisions at = 13 TeV with the ATLAS detector,” ATLASCONF2015081.
 [18] The CMS Collaboration, “Search for new physics in high mass diphoton events in protonproton collisions at 13TeV,” CMSPASEXO15004.
 [19] K. Harigaya and Y. Nomura, Phys. Lett. B 754, 151 (2016) [arXiv:1512.04850]; Y. Mambrini, G. Arcadi and A. Djouadi, arXiv:1512.04913; M. Backovic, A. Mariotti and D. Redigolo, arXiv:1512.04917; A. Angelescu, A. Djouadi and G. Moreau, arXiv:1512.04921; Y. Nakai, R. Sato and K. Tobioka, arXiv:1512.04924; S. Knapen, T. Melia, M. Papucci and K. Zurek, arXiv:1512.04928; D. Buttazzo, A. Greljo and D. Marzocca, arXiv:1512.04929; A. Pilaftsis, Phys. Rev. D 93, 015017 (2016) [arXiv:1512.04931]; S. Di Chiara, L. Marzola and M. Raidal, arXiv:1512.04939; T. Higaki, K. S. Jeong, N. Kitajima and F. Takahashi, arXiv:1512.05295; S. D. McDermott, P. Meade and H. Ramani, arXiv:1512.05326; M. Low, A. Tesi and L. T. Wang, arXiv:1512.05328; B. Bellazzini, R. Franceschini, F. Sala and J. Serra, arXiv:1512.05330; R. S. Gupta, S. Jäger, Y. Kats, G. Perez and E. Stamou, arXiv:1512.05332; C. Petersson and R. Torre, arXiv:1512.05333; E. Molinaro, F. Sannino and N. Vignaroli, arXiv:1512.05334; B. Dutta, Y. Gao, T. Ghosh, I. Gogoladze and T. Li, arXiv:1512.05439; Q. H. Cao, Y. Liu, K. P. Xie, B. Yan and D. M. Zhang, arXiv:1512.05542; S. Matsuzaki and K. Yamawaki, arXiv:1512.05564; A. Kobakhidze, F. Wang, L. Wu, J. M. Yang and M. Zhang, arXiv:1512.05585; R. Martinez, F. Ochoa and C. F. Sierra, arXiv:1512.05617; P. Cox, A. D. Medina, T. S. Ray and A. Spray, arXiv:1512.05618; D. Becirevic, E. Bertuzzo, O. Sumensari and R. Z. Funchal, arXiv:1512.05623; J. M. No, V. Sanz and J. Setford, arXiv:1512.05700; S. V. Demidov and D. S. Gorbunov, arXiv:1512.05723; W. Chao, R. Huo and J. H. Yu, arXiv:1512.05738; S. Fichet, G. von Gersdorff and C. Royon, arXiv:1512.05751; D. Curtin and C. B. Verhaaren, arXiv:1512.05753; L. Bian, N. Chen, D. Liu and J. Shu, arXiv:1512.05759; J. Chakrabortty, A. Choudhury, P. Ghosh, S. Mondal and T. Srivastava, arXiv:1512.05767; A. Ahmed, B. M. Dillon, B. Grzadkowski, J. F. Gunion and Y. Jiang, arXiv:1512.05771; P. Agrawal, J. Fan, B. Heidenreich, M. Reece and M. Strassler, arXiv:1512.05775; C. Csaki, J. Hubisz and J. Terning, arXiv:1512.05776; A. Falkowski, O. Slone and T. Volansky, arXiv:1512.05777; D. Aloni, K. Blum, A. Dery, A. Efrati and Y. Nir, arXiv:1512.05778; Y. Bai, J. Berger and R. Lu, arXiv:1512.05779; E. Gabrielli, K. Kannike, B. Mele, M. Raidal, C. Spethmann and H. Veermäe, arXiv:1512.05961; R. Benbrik, C. H. Chen and T. Nomura, arXiv:1512.06028; J. S. Kim, J. Reuter, K. Rolbiecki and R. R. de Austri, arXiv:1512.06083; A. Alves, A. G. Dias and K. Sinha, arXiv:1512.06091; E. Megias, O. Pujolas and M. Quiros, arXiv:1512.06106; L. M. Carpenter, R. Colburn and J. Goodman, arXiv:1512.06107; J. Bernon and C. Smith, arXiv:1512.06113; W. Chao, arXiv:1512.06297; M. T. Arun and P. Saha, arXiv:1512.06335; C. Han, H. M. Lee, M. Park and V. Sanz, arXiv:1512.06376; S. Chang, arXiv:1512.06426; I. Chakraborty and A. Kundu, arXiv:1512.06508; R. Ding, L. Huang, T. Li and B. Zhu, arXiv:1512.06560; H. Han, S. Wang and S. Zheng, arXiv:1512.06562; X. F. Han and L. Wang, arXiv:1512.06587; M. X. Luo, K. Wang, T. Xu, L. Zhang and G. Zhu, arXiv:1512.06670; J. Chang, K. Cheung and C. T. Lu, arXiv:1512.06671; D. Bardhan, D. Bhatia, A. Chakraborty, U. Maitra, S. Raychaudhuri and T. Samui, arXiv:1512.06674; T. F. Feng, X. Q. Li, H. B. Zhang and S. M. Zhao, arXiv:1512.06696; O. Antipin, M. Mojaza and F. Sannino, arXiv:1512.06708; F. Wang, L. Wu, J. M. Yang and M. Zhang, arXiv:1512.06715; J. Cao, C. Han, L. Shang, W. Su, J. M. Yang and Y. Zhang, arXiv:1512.06728; F. P. Huang, C. S. Li, Z. L. Liu and Y. Wang, arXiv:1512.06732; W. Liao and H. Q. Zheng, arXiv:1512.06741; J. J. Heckman, arXiv:1512.06773; M. Dhuria and G. Goswami, arXiv:1512.06782; X. J. Bi, Q. F. Xiang, P. F. Yin and Z. H. Yu, arXiv:1512.06787; J. S. Kim, K. Rolbiecki and R. R. de Austri, arXiv:1512.06797; L. Berthier, J. M. Cline, W. Shepherd and M. Trott, arXiv:1512.06799; W. S. Cho, D. Kim, K. Kong, S. H. Lim, K. T. Matchev, J. C. Park and M. Park, arXiv:1512.06824; J. M. Cline and Z. Liu, arXiv:1512.06827; M. Bauer and M. Neubert, arXiv:1512.06828; M. Chala, M. Duerr, F. Kahlhoefer and K. SchmidtHoberg, arXiv:1512.06833; D. Barducci, A. Goudelis, S. Kulkarni and D. Sengupta, arXiv:1512.06842; S. M. Boucenna, S. Morisi and A. Vicente, arXiv:1512.06878; C. W. Murphy, arXiv:1512.06976; P. Athron, D. Harries, R. Nevzorov and A. G. Williams, arXiv:1512.07040; A. E. C. Hern¨¢ndez and I. Nisandzic, arXiv:1512.07165; U. K. Dey, S. Mohanty and G. Tomar, arXiv:1512.07212; G. M. Pelaggi, A. Strumia and E. Vigiani, arXiv:1512.07225; J. de Blas, J. Santiago and R. VegaMorales, arXiv:1512.07229; A. Belyaev, G. Cacciapaglia, H. Cai, T. Flacke, A. Parolini and H. Ser ôdio, arXiv:1512.07242; P. S. B. Dev and D. Teresi, arXiv:1512.07243; W. C. Huang, Y. L. S. Tsai and T. C. Yuan, arXiv:1512.07268; S. Moretti and K. Yagyu, arXiv:1512.07462; K. M. Patel and P. Sharma, arXiv:1512.07468; M. Badziak, arXiv:1512.07497; S. Chakraborty, A. Chakraborty and S. Raychaudhuri, arXiv:1512.07527; Q. H. Cao, S. L. Chen and P. H. Gu, arXiv:1512.07541; W. Altmannshofer, J. Galloway, S. Gori, A. L. Kagan, A. Martin and J. Zupan, arXiv:1512.07616; M. Cvetic̆, J. Halverson and P. Langacker, arXiv:1512.07622; J. Gu and Z. Liu, arXiv:1512.07624.
 [20] J. EliasMiro, J. R. Espinosa, G. F. Giudice, H. M. Lee and A. Strumia, JHEP 1206, 031 (2012) [arXiv:1203.0237].
 [21] S. Di Chiara, V. Keus and O. Lebedev, Phys. Lett. B 744, 59 (2015) [arXiv:1412.7036].
 [22] J. McDonald, Phys. Rev. D 50, 3637 (1994) [hepph/0702143]; B. Patt and F. Wilczek, hepph/0605188.
 [23] G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 92, no. 3, 032004 (2015) [arXiv:1504.05511].
 [24] CMS Collaboration, “Search for an Higgs Like resonance in the diphoton mass spectra above 150 GeV with 8 TeV data”, Tech. Rep. CMSPASHIG14006;
 [25] R. Franceschini et al., arXiv:1512.04933.
 [26] S. Knapen, T. Melia, M. Papucci and K. Zurek, arXiv:1512.04928 [hepph].
 [27] J. F. Gunion, H. E. Haber, G. L. Kane and S. Dawson, Front. Phys. 80, 1 (2000).
 [28] S. Heinemeyer et al. [LHC Higgs Cross Section Working Group Collaboration], arXiv:1307.1347.
 [29] M. Bowen, Y. Cui and J. D. Wells, JHEP 0703, 036 (2007) [hepph/0701035].
 [30] F. Bojarski, G. Chalons, D. LopezVal and T. Robens, arXiv:1511.08120.
 [31] A. Falkowski, C. Gross and O. Lebedev, JHEP 1505, 057 (2015) [arXiv:1502.01361].
 [32] G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 91, no. 11, 112011 (2015) [arXiv:1503.05425 [hepex]]; G. Aad et al. [ATLAS Collaboration], JHEP 1508, 105 (2015) [arXiv:1505.04306 [hepex]]; V. Khachatryan et al. [CMS Collaboration], JHEP 1506, 080 (2015) [arXiv:1503.01952 [hepex]]; V. Khachatryan et al. [CMS Collaboration], arXiv:1507.07129 [hepex].
 [33] M. Son and A. Urbano, arXiv:1512.08307; Y. Hamada, T. Noumi, S. Sun and G. Shiu, arXiv:1512.08984.
 [34] M. Kadastik, K. Kannike and M. Raidal, Phys. Rev. D 80, 085020 (2009) [arXiv:0907.1894].
 [35] R. N. Lerner and J. McDonald, Phys. Rev. D 80, 123507 (2009) [arXiv:0909.0520].
 [36] M. Gonderinger, Y. Li, H. Patel and M. J. RamseyMusolf, JHEP 1001, 053 (2010) [arXiv:0910.3167].
 [37] M. Gonderinger, H. Lim and M. J. RamseyMusolf, Phys. Rev. D 86, 043511 (2012) [arXiv:1202.1316].
 [38] M. L. Xiao and J. H. Yu, Phys. Rev. D 90, no. 1, 014007 (2014) [arXiv:1404.0681].
 [39] M. E. Machacek and M. T. Vaughn, Nucl. Phys. B 222, 83 (1983); Nucl. Phys. B 236, 221 (1984); Nucl. Phys. B 249, 70 (1985).
 [40] Q. H. Cao, Y. Q. Gong, X. Wang, B. Yan and L. L. Yang, arXiv:1601.06374; S. F. Ge, H. J. He, J. Ren and Z. Z. Xianyu, arXiv:1602.01801.
 [41] A. V. Bednyakov, B. A. Kniehl, A. F. Pikelner and O. L. Veretin, Phys. Rev. Lett. 115, no. 20, 201802 (2015) [arXiv:1507.08833].
 [42] W. Chao, R. Huo and J. H. Yu, arXiv:1512.05738.
 [43] O. Lebedev, Eur. Phys. J. C 72, 2058 (2012) [arXiv:1203.0156]; P. S. Bhupal Dev, D. K. Ghosh, N. Okada and I. Saha, JHEP 1303, 150 (2013) [arXiv:1301.3453]; E. J. Chun, H. M. Lee and P. Sharma, JHEP 1211, 106 (2012) [arXiv:1209.1303]; C. S. Chen and Y. Tang, JHEP 1204, 019 (2012) [arXiv:1202.5717]; J. Jaeckel, M. Jankowiak and M. Spannowsky, Phys. Dark Univ. 2, 111 (2013) [arXiv:1212.3620]; L. Basso, O. Fischer and J. J. van Der Bij, Phys. Lett. B 730, 326 (2014) [arXiv:1309.6086].
 [44] M. Delmastro, Diphoton searches in ATLAS, 51st Rencontres de Moriond EW 2016. See the slides at https://indico.in2p3.fr/event/12279/session/12/contribution/163/material/slides/1.pdf; The ATLAS Collaboration, “Search for resonances in diphoton events with the ATLAS detector at ”, ATLASCONF2016018.
 [45] P. Musella, Search for high mass diphoton resonances at CMS, 51st Rencontres de Moriond EW 2016, https://indico.in2p3.fr/event/12279/session/12/contribution/218/material/slides/0.pdf; The CMS Collaboration, “Search for new physics in high mass diphoton events in of protonproton collisions at and combined interpretation of searches at 8 TeV and 13 TeV”, CMSPASEXO16018.