750 GeV resonance in the gauged U(1)'-extended MSSM

Recently the ATLAS and CMS Collaborations at the LHC announced their observation of a potential 750 GeV di-photon resonance, after analyzing the √ s = 13 TeV LHC data. This observation has signiﬁcant implications for low-energy supersymmetry. Beyond the MSSM and the NMSSM, we study the MSSM- extensions with an extra U ( 1 ) (cid:2) gauge symmetry. The anomaly cancellation and the spontaneous breaking of the non-decoupled U ( 1 ) (cid:2) generally require introducing vector-like supermultiplets (both colored and color-neutral ones) and singlet supermultiplets, respectively. We illustrate that the potential 750 GeV resonance ( Y ) can be accommodated in various mechanisms, as a singlet-like scalar or pseudoscalar. Three benchmark scenarios are presented: (1) vector-like quarks (VLQ) mediated pp → Y → γγ ; (2) scalar VLQ mediated pp → Y → γγ ; (3) heavy scalar (pseudo-scalar) H / A associated production pp → H ∗ / A ∗ → Y H / h . Additionally, we notice that the Z (cid:2) -mediated vector boson fusion production and Z (cid:2) -associated production pp → Yqq (cid:2) , if yielding a signal rate of the observed level, might have been excluded by the searches for Z (cid:2) via Drell–Yan process at the LHC. © 2016 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP 3 .


Introduction
Recently the ATLAS [1] and CMS [2] Collaborations at the Large Hadron Collider (LHC) reported their first results on the di-photon analysis at energy √ s = 13 TeV. Both collaborations have found an excess over the background at the di-photon (γ γ ) invariant mass m γ γ ∼ 750 GeV, suggesting a signal cross section of 10 ± 3 fb at the ATLAS and 6 ± 3 fb at the CMS. The local significance of this signal is modest, about 3.9σ and 2.6σ at the ATLAS and CMS, respectively. Meanwhile, no notable excess was observed in any other channels, including the tt, hh, W W , Z Z and di-jet final states. Also, this resonance was not seen in the 8 TeV Run-1 data, although both the ATLAS and CMS presented a mild upward protuberance in their di-photon analyses. This observation might indicate a new breakthrough achieved at the LHC after the discovery of the 125 GeV standard model (SM)-like Higgs boson, though more measurements and statistics are needed for us to make a conclusive statement. If it turns out to be real, this observation may have deep implications for new physics. So far, many scenarios have been proposed for interpret-ing this di-photon excess either in concrete models or by means of model-independent EFT approach. In either case, new particles such as vector-like fermions [3][4][5][6][7], hypercharged scalars [8] or exotic vector boson [9] are introduced to induce the Y gg and/or Y γ γ loops for the 750 GeV resonance Y , which has null or very small couplings with the SM particles. This idea has also been implemented in the context of non-supersymmetric extensions of the SM with an extra U (1) or enlarged gauge symmetries [10]. On the other hand, recent studies also discuss an axion-like explanation of the di-photon excess [11], in the context of composite models [12]. Beyond the 4D models, Ref. [13] considers the resonance as a radion from the extra dimension at the conformal limit. In addition, one has extensively studied the connection of this potential 750 GeV resonance to dark matter [14], neutrino physics [15] or more fundamental theories [16]. Other than the direct production of di-photon from a heavy resonance at 750 GeV, exotic possibilities that the di-photon signal arises from a three-body decay [17] along with an undetectable particle or from two bunches of photon jets emitted by a pair of pseudoscalars [18] have also been explored.
This naturally raises the question if this resonance can be accommodated in supersymmetric theories because of their popularity in particle physics. This topical issue has been partially addressed in Ref. [19] heavy CP-even Higgs H or the CP-odd Higgs A as the lightest CPeven Higgs h must be at 125 GeV. However, it has been shown in [4,5,20] [21][22][23][24][25][26][27][28][29]). In the next-to-MSSM (NMSSM), the strong tension between the m h and tan β is slightly relaxed. But it is still difficult to accommodate a 750 GeV resonance with the observed properties. The reason is that supersymmetry may impose large corrections to the h couplings, spoiling their SM-like properties, if it plays a non-trivial role in defining the properties of the 750 GeV resonance. As a result, it is wellmotivated to consider the MSSM extensions which can (1) give additional contributions to the mass of the SM-like Higgs boson at tree level [29][30][31][32] either via new F-terms or via new D-terms or via both in the Higgs potential, and (2) introduce new candidates for the 750 GeV resonance, without spoiling the properties of the SM-like Higgs boson.
To this end, we consider the possibility of the MSSM extension with an additional gauge symmetry U (1) , with R-parity conserved and no CP-violation. Such a gauge symmetry extensively exists in supersymmetric theories, e.g., the GUT (for a review, see, e.g., [33]).
In this paper we focus on a gauged Peccei-Quinn (PQ) symmetry, under which the Higgs fields carry non-trivial charges by definition [34]. If its D-term is non-decoupling, the additional gauge symmetry can provide non-trivial contribution to the tree-level mass of the SM-like Higgs boson [35]. Additionally, an effective μ parameter could be dynamically generated in this framework [36,37]. It is known that the bare μ term is forbidden by the U (1) gauge symmetry. Whereas, an effective one can be produced via W H ∼ λSH u H d , yielding μ eff = λv S . Here H d and H u are down-and up-type Higgs supermultiplets, and S is a SM-singlet supermultiplet. In this context, the potentially existing 750 GeV resonance Y can be interpreted as a singlet-like scalar or a singlet-like pseudoscalar.
There are three classes of couplings which are relevant to our studies below. They could be chiral under the U (1) gauge symmetry. It is well-known that the anomaly cancellation for the whole set of the gauge symmetries generically requires the introduction of both colored and color-neutral vector-like supermultiplets, unless some Chern-Simons term is properly incorporated. Because the singlet supermultiplets carry a U (1) charge as well, they may couple with these vector-like supermultiplets, yielding a mass for these exotics. • The resonance couplings with the U (1) gauge boson ∼ Y Z Z , since the singlet supermultiplet carries a U (1) gauge charge.
In the presence of these couplings, Y could be produced with a relatively large boost for the cross section at 13 TeV LHC, compared to that at the 8 TeV LHC. The main production mechanisms include: 1. Vector-like quark (VLQ) mediated gluon fusion. Y is singletlike, with its coupling to the SM fermions and their superpartners suppressed by a small mixing angle. But its coupling with VLQ can yield a large effective coupling with di-gluon (as well as di-photon), resulting in a boost factor as large as 4.7. 2. Scalar VLQ (SLVQ) mediated gluon fusion. 3. H/ A-mediated production. Compared to the first two mechanisms, a relatively large boost factor arises due to the increased energy threshold. 4. Z -mediated vector boson fusion and Z -associated production. Compared to the other mechanisms, these ones are constrained by the current experimental bounds much more tightly, unless the leptonic decay of Z can be greatly suppressed.
With these production mechanisms, the 750 GeV resonance can be explained, given a decay branching ratio Br(Y → γ γ ) of subpercent level or above. This could be achieved with its effective di-photon coupling mediated by VLQ (SVLQ) and VLL (SVLL). We organize the paper as follows. In Section 2, we present a general discussion on the MSSM extension with a PQ-type U (1) gauge symmetry, as well as the properties of the candidate for the 750 GeV resonance. In Section 3, we identify three benchmark scenarios where the experimental data can be properly fit, using an anomaly-free model for illustration. Section 4 contains our conclusions.

The gauged U (1) -extended MSSM: general discussions
In the gauged U (1) -extended MSSM, a list of chiral supermultiplets in the Higgs sector typically includes (see, e.g., [33]) Here q a and v a denote the U (1) charge carried by H a , and its VEV, respectively.
Then we can write down the tree-level Z − Z mass matrix Here (5) which is constrained to be of or less than O(10 −3 ) by EW precision tests [38].
In this case, the set of {S i } define a so-called secluded U (1) -breaking sector [39]: they mainly serve as an order parameter of the U (1) symmetry breaking, yielding a hierarchy between the Z mass and the EW scale. Below we will be working in this context. For simplicity we further turn off the couplings between the {H d , H u , S} sector and the secluded-U (1) breaking sector in the superpotential and soft SUSY-breaking terms. Then the influence of the U (1) gauge symmetry on the Higgs potential at tree level enters via D-terms only.
The mixing between the {H d , H u , S} sector and the secluded-U (1) breaking sector can be safely neglected. Without loss of generality, we interpret the 750 GeV resonance as an S-like scalar (the analysis is similar in most cases if the S-like pseudoscalar is the candidate for the 750 GeV resonance). With these assumptions, the superpotential and soft supersymmetry-breaking terms of the Higgs sector are given by yielding a tree-level mass matrix for the CP-even Higgs sector in where κ 2 Given the 125 GeV Higgs that has been observed at the LHC is SM-like, we naturally consider the alignment limit with negligibly small tree-level mixing between S and H d , H u . Such a limit can be achieved by an assumption, e.g., Next, using the minimization conditions we replace m 2 Then the tree-level mass eigenvalues are solved to be Here h is the SM-like Higgs, which receives tree-level contributions to its mass from both F-terms and D-terms. H is the MSSM-like CP-even Higgs and Y is the singlet-like one. As discussed above, the secluded U (1) -breaking sector mainly serves as an order parameter of the U (1) symmetry breaking. Indeed, the -dependent terms which involve v S i are gone in Eq. (11). In this context the coupling g Y Hh and g Y hh generically vanish at tree level. We do not expect Y be produced via H → Y h, even if H might be heavier than Y . Interestingly, the coupling could be sizable. The three terms in the first line represent contributions from D-terms, F-terms, and softly SUSY-breaking terms, respectively. To accomplish the second line, the equality relation (10) has been taken and m 2 Y ≈ 8g 2 Z v 2 S has been assumed.
The superpotential and softly supersymmetry-breaking terms for describing the interaction between S and the vector-like fermions and their superpartners are Then the effective couplings between Y and di-gluons and diphotons are yielded by if there is no mixing between F V and F c V . This indicates that κ is a crucial parameter in determining the production and decay patterns of the Y resonance. Additionally, the coupling between Y and di-Z bosons is given by Using these outcomes, it is straightforward to calculate the signal rates in various scenarios.

Benchmark scenarios
The PQ-type U (1) with the MSSM matter content only is anomalous. To cancel the mixing anomalies between the SM gauge symmetries and U (1) , vector-like color triplets and EW multiplets carrying U (1) charges are typically required, no matter in the GUT (see, e.g., [40,41]) or in the non-GUT (see, e.g., [42]) frameworks. Such a construction automatically warrants that no new anomalies are introduced for the SM gauge symmetries. As an illustration, we will work on the model built in [42] with a mild modification. The particle spectrum is listed in Table 1, which includes three classes of supermultiplets: the ones defined in the MSSM, the exotic ones charged under the SM gauge symmetries and their signlets. All of these supermultiplets carry the U (1) gauge charges. Note that in   Some relevant terms in the superpotential, W F , for the second class of supermultiplets are given by Here i, j run over 1, 2 and 3, and r, s run over 1 and 2. After the U (1) symmetry is broken, {T i , T c i }, {D r , D c r } and {X r , X c r } obtain vector-like masses by coupling with the supermultiplet S (or S c 1 ). The supermultiplet S 2 may also contribute if its scalar field obtains a non-zero VEV and meanwhile, couples with these vectorlike supermultiplets. An important feature is that these vector-like supermultiplets can decay into the SM particles via proper interactions listed above, avoiding the overproduction in the Universe. Some searches related to these exotics have been pursued at the LHC. The current bounds are dependent on the SM gauge charges carried by these exotics as well as their couplings with the SM particles. {T i , T c i }, in essence, are supermultiplets of leptoquarks.
The LHC searches for leptoquarks have been pursued at both the ATLAS and the CMS [43][44][45]. If they exclusively decay into the SM fermions of third-generation, the data collected at Run-1 set a lower bound for their mass ∼ 600 − 650 GeV [43,45]. As for their fermionic superpartners, if they also exclusively decay into the same final states with one additional lightest superparticle which is neutral, comparable mass bounds are expected. Whereas, {D r , D c r } and {X r , X c r } are colorless supermultiplets. As a matter of fact, the gauge charges carried by {D r , D c r } are exactly the same as the ones carried by {H u , H d }. Thus we do not expect a strong bound for their mass in a general context. In this article, we simply assume a universal lower bound of 650 GeV for the mass of both leptoquarks and their superpartners, as indicated in Fig. 1 and Fig. 2.
In the following we discuss three representative scenarios that could yield the di-photon signal for the 750 GeV resonance. In all cases we assume no mixing between the Higgs doublets and the singlet S, by demanding the relation (10). The pure singlet scalar S accounts for the 750 GeV resonance.  if m H > m Y , the decay width of H → Y H * is too small, compared with that of H → tt, to make this process significant. So we will focus on the former case. We present in Fig. 1 the contours of the di-photon signal rate Remarkably, the desired signal rate σ (pp → Y → γ γ ) of order ∼ O(10) fb can be achieved in both Scenario I and II, and hence fits well both the ATLAS and CMS signals at 13 TeV LHC. The signal rate in Scenario II shows a strong dependence on the mediator mass, simply because the effective di-photon coupling is proportional to a factor of 1/m 2 other than a loop- 1 Unlike the VLQ, the SVLQ may receive contributions to its mass not only from F-terms, but also from D-terms and softly supersymmetry-breaking terms. For sake of brevity, we are not attempting to explicitly write out the full expression. function. Whereas, such a strong dependence on the mediator mass is absent in Scenario I (mediated by a fermion VLQ). Compared to Scenario I, a larger κ 1 is typically required to yield a comparable signal rate in Scenario II. This partially results from the fact that the ratio between a scalar-loop function and a fermionloop function asymptotically approaches a value A f /A s 4 as the mediator mass increases. The di-jet search at the Run-1 have a marginal impact on eliminating the parameter space. The projection of the Run-1 limit to the 13 TeV Run-2 yields a bound σ (gg → Y → gg) < 100 pb [46]. Since the partial width of digluon decay is only ∼ 400 times larger than that of the di-photon decay in our case, the di-jet bound is weak in both Scenario I and Scenario II. The results of Scenario III are displayed in Fig. 3. Clearly, there exists parameter space where a cross section of O(10) fb or above for the Y production in association with a nonstandard Higgs boson H is obtainable. If the di-photon coupling of Y is mediated by VLL or SVLL only, such as the ones in {X 2 , X c 2 }, the branching ra-  [47]. Thus, even if this scenario works, we believe it will be extremely fine-tuning. Finally, we collect in Table 2 three benchmarks points (BMPs), representing the three benchmark scenarios discussed above, respectively. They are all able to produce the di-photon signal rates of the order observed. Should future analyses confirm this signal, it will be of interest to fully test these scenarios using complementary decay channels in experiment. Using the Run-1 data, both the ATLAS and the CMS have analyzed the upper limits on the production cross section of new particle Y at 750 GeV which decay into various final states including Z Z , W W , hh, Zh, Z γ , τ τ and tt. A crude estimation of these upper limits, rescaled by the upper limit of σ (gg → Y → γ γ ) of an order ∼ O(1) fb, has been offered in [48]: R(Z Z, W W , hh, Zh, Z γ , τ τ , tt) (23,46,41,20,7,20,600), with R( X X) = BR(Y → X X)/BR(Y → γ γ ). The Z γ search thus may play a crucial role in this regard. Additionally, these scenarios could be probed by searching for the mediators of the di-photon couplings. We can also probe Scenario III by searching for nonstandard Higgs bosons via its decay into tt. At last, we would point out -the couplings κ 1 and λ of ∼ O(1) may yield a Laudau problem. Particularly in Scenario I and II, κ 1 may hit its Landau pole at a energy scale as low as 10 TeV. We refer to [49] for more discussions in this regard. Turning on more colored mediators, such as the component field in {T 2,3 , T c 2,3 } in this model, in the di-photon and di-gluon loops doesn't help much, since their relatively small hypercharges or electric charges further suppress the branching ration of the di-photon decay. Alternatively, one may turn on colorneutral mediators simultaneously, which may greatly enhances the branching ratio of the di-photon decay. In this case, the kappa value required for generating the signal could be suppressed to be below O(1).

Conclusions
Motivated by the recent observation of a potential 750 GeV di-photon resonance at the √ s = 13 TeV LHC, we studied its implications for low-energy supersymmetric theories. Though it is generically difficult to explain the observed signal rate in the MSSM and in the NMSSM, due to the lack of an efficient mechanism to generate the observed signal rate, we pointed out that such a resonance can be accommodated in the MSSM-extensions with a PQ-type U (1) gauge symmetry. For illustration, we explore several benchmark scenarios in which the di-photon resonance is produced either via (scalar) VLQ-mediated gluon fusion process or in association with a heavy Higgs scalar. Additionally, we studied the possibility of Z -mediated vector boson fusion production and Z -associated production: pp → Y qq . Compared to the other scenarios, this possibility is constrained by the current experimental bounds much more strongly. Though our studies were focused on a CP-even candidate resonance, the discussions can be extended to a CP-odd candidate resonance, with a couple of exceptions. First, the Z -mediated mechanism can not be applied to the latter case, unless CP-violation is allowed. Additionally, in this context a CP-odd singlet prefers to couple with a MSSM-like CP-odd Higgs A and a SM-like Higgs h, instead of the A plus a heavier Higgs H . Therefore, we may well expect a significant production of the resonance in association with a SM-like Higgs h, mediated by an off-shell A.
As a candidate theory to explain the potential 750 GeV diphoton resonance, this class of supersymmetric models have three main advantages. First, the singlet supermultiplet (the one providing the candidate of the 750 GeV resonance) is generically required for dynamically generating an effective μ parameter. Second, the vector-like supermultiplets are generically required by anomaly cancellation of the whole set of gauge symmetries, unless some Chern-Simons term is properly incorporated. Third, several mechanisms work well in this context, enabling the 750 GeV resonance to be accommodated in multiple ways, if it does exist. All of them do not require a mixing between the singlet and the Higgs doublets in the 750 GeV resonance. This certainly broadens our understanding on both the potential 750 GeV resonance and the low-energy scale supersymmetric theories.
In spite of these advantages, the resonance exclusively decays into gg, γ γ , Z γ and Z Z in Scenarios I and II and into γ γ , Z γ and Z Z in Scenario III. This yields a partial width of Y → γ γ typically narrow, in comparison to ∼ 45 GeV which is indicated by the ATLAS measurement [1]. This effect may not cause a decay length longer than sub-millimeter, e.g., by raising the coupling value between Y and {X 2 , X c 2 } in Scenario III. Nonetheless, it is necessary to think about how to reconcile this discrepancy. Given the LHC may release more precise information on the properties of this resonance soon, we leave a further study in this regard to future work.