Novel kinematics from a custodially protected diphoton resonance

Novel kinematics from a custodially protected diphoton resonance

Jack H. Collins    Csaba Csáki    Jeff A. Dror    Salvator Lombardo Department of Physics, LEPP, Cornell University, Ithaca, NY 14853, USA

We study a simple, well-motivated model based on a custodial symmetry which describes the tree-level production of a 750 GeV diphoton resonance from a decay of a singly produced vector-like quark. The model has several novel features. The identification of the resonance as an SU(2) triplet provides a symmetry explanation for suppression of its decays to , , and . Moreover, the ratio of the 13 TeV to 8 TeV cross sections can be larger than single production of a 750 GeV resonance, reaching ratios of up to 7 for TeV scale vector-like quark masses. This eliminates any tension between the results from Run I and Run II diphoton searches. Lastly, we study the kinematics of our signal and conclude that the new production mechanism is consistent with available experimental distributions in large regions of parameter space but, depending on the mass of the new vector-like quarks, can be differentiated from the background with more statistics.

I Introduction

The recent observation of an excess in the diphoton channel around 750 GeV invariant mass by ATLAS and CMS at TeV ATLAS-CONF-2015-081 (); CMS-PAS-EXO-15-004 () has generated much interest in models with a heavy scalar resonance, , that decays to two photons. Most explanations proposed so far are considering loop induced resonance production, typically via heavy vector-like quarks (VLQ) charged under the Standard Model (SM). Otherwise, tree-level decays to SM particles would naturally dominate the branching ratio of , either leading to a diphoton rate too small to explain the excess or a production rate of two SM particles with large invariant mass that is excluded by existing measurements.

In this paper, we propose a novel tree-level production mechanism where arises from the decay of a VLQ. The VLQs can be singly produced due to their mixing with the SM quarks, while the resonance is protected by the SU(2)SU(2) custodial symmetry. In order to have significant mixing between the VLQs and the light quarks without modifying the couplings predicted by the SM, we introduce VLQs in a bidoublet representation of the custodial symmetry, while the resonance is part of a triplet under SU(2). The model has several advantages and new features:

  • It is one of the few viable examples of tree-level production consistent with the excess signal rate, existing experimental constraints, and the kinematic distributions of the diphoton background events.

  • The ratio of production rates between 13 and 8 TeV is different than gluon or quark fusion. Depending on the model the ratio can be as large as about 7 (vs. 4.7 for gluon fusion), eliminating the tension with the 8 TeV diphoton searches.

  • The custodial symmetry protects the resonance from the leading one-loop decays to , , and , while allowing decays to . The suppression of the decay is particularly significant since in most models this coupling will arise at tree-level, making it difficult to reconcile with the expectantly large diphoton branching ratio and unobserved decays. The (as well as and ) decay width is nonvanishing due to the explicit breaking of the custodial symmetry from gauging the U(1) subgroup of SU(2).

  • The custodial symmetry also forbids one-loop gluon fusion production, explaining the dominance of the tree-level production via a decay of a VLQ.

The scenario where is produced primarily from the decay of singly-produced VLQs has not been considered in the diphoton excess literature, although some authors have considered production of through a cascade decay of a heavier parent particle (e.g., Franceschini:2015kwy (); Knapen:2015dap (); Han:2016bus (); Huang:2015evq (); Liu:2015yec (); Berlin:2016hqw (); Bernon:2016dow (); DeRomeri:2016xpb ()). Furthermore, several authors have pointed out the potential to measure VLQ top or bottom partners decaying to  Agrawal:2015dbf (); Kobakhidze:2015ldh (); Dev:2015vjd (). However, to explain the bulk of the excess signal through a decay of pair-produced top or bottom partners would require couplings at their perturbative limits to achieve large enough rate and to explain why the VLQ decay process dominates over gluon fusion Han:2016bus (). Moreover, pair production of two VLQs per event would give several hard jets in the event in addition to the diphoton, which is inconsistent with the kinematic distributions of events in the excess region. The tree-level production mechanism presented in this paper avoids these problems. Interestingly, it has been pointed out in Delaunay:2013pwa () that a top partner need not be a mass eigenstate but rather could be a mixture of top and charm-like mass eigenstates raising the possibility that the vector-like quarks which mix with light quarks could also play a role in solving the little hierarchy problem.

This paper is organized as follows. In Sec. II, we introduce a motivated model consistent with low energy flavor constraints in which a bidoublet of vector-like quarks with mix with the light-SM quarks. In Sec. III we discuss the production and decays of the new particles and find that the model can easily accommodate the current excess in the diphoton data without tension from existing searches, both from the searches sensitive to the 750 GeV resonance as well as searches sensitive to VLQs. Furthermore, we compare the kinematic distributions of our signal and the diphoton background, finding that, depending on the splitting between the VLQ and the scalar, the distributions can be challenging to distinguish without additional data.

Ii Custodial symmetry and light quark mixing

In this section, we present a model in which the resonance is a decay product of an electroweak produced VLQ from the dominant -channel process shown in figure 1. Single production of VLQs requires a large mixing angle between the light quarks and VLQs. One would naively expect such mixing to yield large corrections to the couplings, which are strongly constrained by electroweak precision observables. However, we can protect the couplings via the custodial symmetry by using a bidoublet representation for the VLQs Agashe:2006at (). In addition, we introduce the SU(2) triplet, , whose neutral component (which we denote by ) will play the role of the 750 GeV resonance. This model is an example of a custodially symmetric model commonly considered in composite Higgs and extra dimensional models Atre:2008iu (); Atre:2011ae (); Delaunay:2013pwa (); Redi:2011zi (); Flacke:2013fya (); Redi:2013eaa (); Vignaroli:2012nf (); Vignaroli:2012sf (); Atre:2013ap ().



Figure 1: The dominant production of the diphoton excess from a decaying VLQ (). In addition to the resonance, there are two additional jets. The of the jets (and hence the visibility of the signal) is strongly dependent on the mass of the VLQ.

ii.1 Field content and mixing

We organize the fields into irreducible representations of SU(3)SU(2)SU(2)U(1) in Table 1, where the SU(2) is a global symmetry and the fermions are left-handed Weyl spinors. Hypercharge is embedded in SU(2)U(1) as . The SM quarks are taken to be singlets under SU(2) and we represent the Higgs doublet as a bifundamental, .

Four VLQs form the bidoublet,


The charges of the new quarks are , , and .

Table 1: The representations of relevant fields. The new vector-like quarks are in a single bidoublet , and the 750 GeV resonance is the neutral component of . The Higgs and the light quarks have the usual SM assignments. All fermionic fields are left-handed Weyl spinors.

The Lagrangian for the VLQs (not including the terms coupling to the scalar ) is


where the superscript denotes the quark fields in the SM mass basis of the up-type sector (the basis of diagonal SM Yukawa couplings). We assume the VLQs only mix with a single generation of right-handed up-type quark, however in appendix A.4 we also consider a different U(1) charge for and the case of mixing with one generation of right-handed down-type quark. Note that we have made an important assumption regarding alignment: the bidoublet couples to the up-type quark in the mass basis of the SM. This assumption is to avoid low energy flavor constraints from flavor changing neutral currents but is not a crucial ingredient for the collider phenomenology that follows.

The down-type quark masses are unaffected by the new VLQs. One flavor of the up-type quarks can mix significantly with the new VLQs through the off-diagonal mass matrix:


where is the Yukawa of the up-type quark. From here on, we assume the mixing is with the up or charm quark and neglect the up-type quark mass. The mass eigenstates (, , and ) are related to the gauge eigenstates by:


where and with


where is the vacuum expectation value of the Higgs, and the masses of the VLQs are and . The down-sector does not experience any mixing, i.e. and are mass eigenstates.

The mixing, parameterized by , leads to couplings between a generation of SM quarks, SM gauge bosons, and the VLQs (derived in more detail in appendix A):


where () is the sine (cosine) of the Weinberg angle and () is the SU(2) (QED) coupling constant. If the mixing angle is sufficiently large, these couplings can result in electroweak production of single VLQs (, , or ) which can dominate over the VLQ pair production cross section. Notice that only , which will be responsible for the production of the diphoton resonance, couples to an up-type quark and the . is not produced by electroweak interactions in this model. This is because only the linear combination mixes with the up-type quark, which is a necessary feature for the protection of the coupling in this model.

ii.2 Consequences of a custodial triplet

The 750 GeV diphoton resonance is embedded in an SU(2) triplet scalar as follows.


This allows for a coupling of to the VLQs of the form


where the ellipses refer to terms involving the charged components of . The relative minus sign between the and terms gives rise to the coupling of to and the SM up quark (as opposed to a coupling to ) which is responsible for the production of the resonance.

These interactions will generate couplings of to SM dibosons, such as , , etc. via triangle diagrams with the VLQs. However, in the limit of exact custodial symmetry, these amplitudes are forbidden. For example, the operator has no custodially invariant contraction because the gluon field strength tensor is a custodial singlet. Furthermore, the Higgs coupling to the scalar vanishes since (using ),


In practice, the vanishing of the amplitudes is a consequence of cancellations of the contributions due to different VLQs running in the loop, which contain important relative minus signs as a consequence of the custodial symmetry.

The loop amplitudes for therefore require the insertion of SU(2) violating interactions. The largest such couplings in the SM are the third generation Yukawas, however in the flavor alignment limit these couplings will not directly affect the diphoton resonance sector since we assume mixing is not occurring with the third generation up-type quark. The dominant source of custodial symmetry breaking in this sector will therefore be the embedding of the hypercharge gauge group within the generator of SU(2), and so it is to be expected that the leading loop amplitude will be that coupling to hypercharge gauge bosons, . Indeed, the one-loop contributions to this operator do not cancel among the VLQs due to their differing hypercharges. The other loop amplitudes will be generated at higher order and are suppressed by an additional factor compared to their naive sizes. These two-loop contributions can induce a mixing angle between and the Higgs of order (where denotes the mass of the resonance), however this is much too small to induce sizable decays to . Direct couplings of to the up-type quarks can also arise at two-loops but is suppressed by a Yukawa coupling, making it negligible. We verify the effects of custodial symmetry breaking explicitly in A.3. This custodial protection mechanism generates a natural hierarchy between the decays of the resonance to diphotons and its decays to , , , and also suppresses the gluon fusion production of .

We now briefly compare this to a scenario in which the diphoton resonance is assumed to be a custodial singlet with couplings


In this case, the only -quark-VLQ coupling involves , which does not couple to SM gauge bosons and therefore cannot be produced via VLQ single production. is pair-produced and can decay , however the rate for pair production is subdominant to electroweak production of and insufficient to explain the excess. Furthermore, this singlet does not exhibit custodial protection, which is a consequence of the couplings in eq. 13 adding constructively rather than destructively.

Iii Diphoton cross section

Above we presented a model in which can be produced as a decay product of a singly produced VLQ (we will assume throughout). The dominant production mechanism for the diphoton resonance is depicted in figure 1.111Secondary production modes from production, pair production of VLQs, and direct production mediated by a VLQ make up 10-30% of inclusive diphoton cross section. We explore the size of different contributions in Sec. IV.2, but since the size of the subdominant modes is highly dependent on the detailed parameters of the model, we only include the dominant production when studying the inclusive diphoton rate and kinematics. In this section, we demonstrate that the rate is sufficient to explain the diphoton excess while avoiding constraints from existing VLQ searches and electroweak precision tests. We consider two variations of the model, one where the VLQs mix with the up quark and another with the VLQs mixing with the charm quark.

Since the final state arises from a decay chain, the inclusive cross section into is given in the narrow width approximation by


Each of these contributions has different dependence on the relevant parameters of the model, and . The production cross section of the VLQs, , is proportional to but is independent of .

iii.1 Branching ratios

The complete formulae for the branching ratios of and are given in appendix A.2. We summarize the results here. has two decay channels, and with the dominant decay being . This results in the branching ratio of ranging between 1-10%, proportional to and independent of the mixing angle. has competing decays between a -body tree-level decay and loop-induced 2-body decays. The only tree-level decay of is to through an off-shell with a rate is proportional to , making it highly sensitive to the mixing angle. has additional loop-induced decays into , , and . These decays arise from gauging hypercharge resulting in the relative ratios,


The loop-induced rates are largely independent of the mixing, proportional to .222We have checked that non-zero mixing has at most a 10% effect on the loop-induced rates in the region of parameter space we are interested in. Furthermore, these effects will not have any bearing on the size of the rate and thus we ignore these effects in our analysis.

0.7 0.1
2 0.3
Table 2: Benchmark points for the up quark and charm quark mixing models. For the model, there are more stringent constraints on the mixing angle though a larger production cross section.

[width=7cm]BR_U.pdf \includegraphics[width=7cm]BR_C.pdf

Figure 2: The different branching ratios for . The loop-induced decays to , , and always compete with the -body decay. At large mixing, the -body decay is the preferred decay mode, however for small-mixing, the loop-induced decays (which are roughly independent of the mixing) dominate. Left: The branching ratios for the benchmark point (small-mixing). Right: The branching ratios for the benchmark (large-mixing).

The leading branching ratios of are shown in figure 2 for a benchmark point relevant in the case of up-mixing with (left) and charm-mixing with (right). A couple comments on the choice of benchmark points (displayed in Table 2) are in order. First, regarding the size of the mixing angle, we have provided constraints on the allowed mixing angle for electroweak produced VLQs in appendix B, obtained by reinterpreting the constraints from direct LHC searches on light quark composite partner models Delaunay:2013pwa (). For the up-quark mixing, the mixing angle is experimentally constrained to be , while for charm mixing, the constraints are much weaker with . Larger mixing angles are allowed for the case of charm mixing since the electroweak production cross section of the VLQs are suppressed by the charm parton distribution function.

iii.2 The inclusive cross section

Depending on the size of mixing, the inclusive cross section for diphoton production scales differently with the mixing angle. There are two distinct regimes, large and small mixing angles. If the mixing angle is large, then the dominant decay of is through the -body decay. In this case the branching ratio into diphotons is giving an inclusive cross section


For small mixing, the -body decay is heavily suppressed making the diphoton rate the dominant mode, i.e. , independent of or . In this case, the inclusive cross section scales as


The transition between the two regimes occurs around , and this is the point where the cross section is maximized. Due to the constraints, the up mixing case is always in the small mixing scenario while the charm can be either the large or small mixing regimes. For our chosen benchmark point, the charm scenario corresponds to large mixing.


[width=7cm]XS_U.pdf \includegraphics[width=7cm]XS_C.pdf

Figure 3: The inclusive cross section into as a function of VLQ mass and varying the values of . In gray we show the rough cross section necessary to explain the excess with a narrow width (Kamenik:2016tuv (). Left: The cross section for the benchmark point. Right: The cross section for the benchmark point.

To reproduce the excess, we simulate the production of at leading order using a custom FeynRules model Alloul:2013bka () with MadGraph5 Alwall:2014hca (). To roughly estimate the size of next-to-leading order (NLO) effects, we compute the production cross section with an additional jet, finding that it makes up about 50% of the leading order cross section. This suggests an NLO -factor of about 1.5, and we use this correction throughout. We compute the diphoton rate using equation 15 and the branching ratios given in appendix A.2. Figure 3 shows the inclusive diphoton production cross section for different VLQ masses and values of for the two benchmark points.

We see in the left of figure 3 that we need to get enough cross section to explain the excess in the up quark variation of the model. Larger Yukawa couplings are required in the charm-mixing benchmark point, requiring to achieve the minimum cross section needed to explain the excess.

As with many models explaining the diphoton excess, such Yukawa couplings can lead to non-perturbativity of the model before the GUT scale. The up-mixing benchmark becomes non-perturbative at around 100 TeV, although there is some parameter space where the coupling remains perturbative beyond the GUT scale. The charm-mixing model is more problematic given the Yukawa coupling runs to its perturbative limit at a few TeV, putting into question the validity of our analysis. However, this problem can be easily overcome by adding additional flavors of bidoublets (these may or may not mix with the SM quarks) which feed into the running, but also boost the diphoton decay rate as the square of the number of flavors allowing for much smaller couplings, and a much higher scale of strong coupling.



Figure 4: The ratio of to TeV cross section of the up-mixing signal (blue) and charm-mixing signal (red). We have included the scaling properties of the down-type version of this model where the down (green), strange (pink), or bottom (light blue) quark mix with VLQs. The scaling of other proposed production processes are shown as dashed lines and were taken from Franceschini:2015kwy ().

iii.3 Eliminating tension with TeV data

One of the puzzling features of the 13 TeV diphoton excess is its seemingly large cross section compared to cross section limits from 8 TeV searches. A 750 GeV resonance is not ruled out by Run I searches but, depending on the production mechanism and width, may be in tension with Run I limits Falkowski:2015swt (); Kamenik:2016tuv (). Thus to reproduce the excess, it is important to have sufficiently large scaling, , defined as the ratio of the cross section at 13 TeV to that at 8 TeV. We compute the scaling for our model as a function of the VLQ mass (the scaling is independent of the couplings) and assuming the -factor is constant from to TeV. The results are shown in figure 4 alongside the scaling of other proposed models, including gluon fusion, production Gao:2015igz (), and photon fusion Csaki:2015vek (); Fichet:2015vvy (); Csaki:2016raa (). The up-mixing model inherits the scaling from the production at GeV but grows with the mass of the VLQ due to the higher center-of-mass energy. For heavy VLQ masses near 1500 GeV, the scaling is comparable to gluon fusion. The charm variation of the model, however, scales much better due to the parton distribution function of the charm in the initial state. For VLQ masses nearly degenerate with ( GeV GeV) for which the extra jet from the decay is relatively soft, the scaling is as large as for production. For larger masses, the charm-mixing scenario has , but in this region of parameter space, the would be accompanied by a high- jet in the final state (we explore the plausibility of this scenario in Sec. IV). We conclude that, depending on the mass of the VLQ, our signal can achieve larger cross section scaling from 8 to 13 TeV than any proposed model of single resonance production.

For simplicity, we have only considered mixing with up-type quarks. It is possible to construct a similar model in which new VLQs mix with the down-type quarks. The terms in the Lagrangian responsible for production of a down-type variation of this model are presented in appendix A.4. We include the scaling properties of production from down, strange, and bottom quark mixing in figure 4. The bottom quark mixing scenario has the largest ratio of 13 to 8 TeV cross section, while the down and strange quark scenarios interpolate between the scaling of the up and charm scenarios.

Iv Kinematics

iv.1 Comparing with ATLAS


[height=4cm]pTaa_withgg.pdf \includegraphics[height=4cm]Njets_withgg.pdf

Figure 5: The kinematic distributions of the sum of the signal and background for vector-like quark mass of (blue) and (red) compared to the distributions observed by ATLAS (black) with 3.2 fb of data. We also provide gluon fusion kinematics (green) for comparison. The and distributions for the background and observed events are obtained from the slides presented by ATLAS Moriond_ATLAS ().

[height=4cm]pTj1_withgg.pdf \includegraphics[height=4cm]etaFJ_withgg.pdf

Figure 6: The leading-jet and forward jet pseudorapidity distributions of the signal for vector-like quark mass of (blue) and (red) along with a gluon fusion signal (green) for comparison. For VLQs almost degenerate with the resonance, the signal is difficult to differentiate from the QCD background or a resonance produced via gluon fusion since these events also contain soft, forward jets from initial state radiation.

In addition to the diphoton resonance signature at , our signal has two additional jets, with one of the jets typically in the forward direction. ATLAS and CMS have remarked that events in the excess region are consistent with the background kinematics. Furthermore, ATLAS has recently provided kinematic distributions of the excess events Moriond_ATLAS () for the number of jets, (ATLAS defined a jet using for and for ), and the transverse momentum of the resonance, . Furthermore, ATLAS provided estimates for the expected SM diphoton background from simulations. The distributions are provided for the region and with the requirements on the leading and subleading photon energies and . In this bin, ATLAS found a total of events, about 10 of which are diphoton excess candidates.

To compare the compatibility of the kinematics of the excess with our signal we simulate our signal using a combination of MadGraph5, Pythia 8.2 Sjostrand:2007gs (), and Delphes 3 deFavereau:2013fsa () making use of the NNPDF2.3LO parton distribution functions Ball:2012cx (). To compare with the distributions observed by ATLAS, we perform a weighted sum of our signal and the background estimates provided by ATLAS


This is done for each bin and we take . We also simulate gluon fusion at leading order and perform this procedure in order to compare our signal with the kinematics of single production. By comparing the distributions in figure 5, we conclude that although our signal has two extra jets in the final state, the distributions for and are consistent with the data provided by ATLAS (as is also the case for the gluon fusion signal). Furthermore, the mass of the VLQ has only a mild effect on these distributions since the number of jets from the hard process is independent of . The signal does have distinctive features in other distributions, however, and we explore these features in the next section.

iv.2 Additional signatures

Our signal predicts observable jet signatures that can be used to discern this process from background events or from other resonance production mechanisms. In particular, we expect a forward jet as well as one central jet with higher , depending on . The distributions for the of the leading jet () and the absolute value of the pseudorapidity of the most forward jet () for the signals with and 1000 GeV are shown in figure 6 along with the gluon fusion signal for comparison. Note that ATLAS and CMS did not provide the background or observed kinematic distributions for these observables, so we did not combine the background and the signal in these plots.

The distributions have some distinctive characteristics. Firstly, we see that for small splitting is peaked around zero since, at these splittings, the central jet, which will typically be the leading jet, has low . However, for larger splitting the distribution has a kinematic edge with the end-point at the splitting between the VLQ and the resonance. This is prototypical of a jet arising from a heavy particle decay to a second heavy particle. Interestingly, such a distribution can suggest the mass of the VLQ. The forward jet in the event is most easily probed using its pseudorapidity. The distribution has a dip at as a consequence of the cuts used by ATLAS for the jet definition (see above). As expected, the signal has a jet with large , however this is also true of the dominant background, . In this background, jets are emitted from the initial state and hence tend to be in the forward region. This can result in the feature being well hidden inside the SM background of the searches.

In addition to the features of the dominant production mode, there are secondary production modes of the excess. In principle, one may expect that any of the VLQs could decay into the resonance, but this is not the case. Due to the custodial structure of the model, only couples to the resonance and a SM quark (see equation 10), and hence its the only single VLQ production mode. However, there are three subdominant modes which can contribute significantly to the cross section, single VLQ production through , QCD pair production of VLQs through , and direct production of the resonance through a -channel VLQ, . The cross section composition depends strongly on the choice of mixing angle and mass of the VLQ. The various contributions to the total cross section as a function of the mass for our benchmark points are shown at leading order in figure 7 (for simplicity we do not apply -factors when comparing between these different channels). We conclude that the additional production modes make up 10-30% of the inclusive diphoton cross section for reasonable choices of parameters. With more statistics, excesses in these subleading channels could be used to differentiate our signal.


[width=7cm]XS_U_comparison.pdf \includegraphics[width=7cm]XS_C_comparison.pdf

Figure 7: The cross sections for the different production modes at our benchmark points. We see that production through single VLQ dominates with the secondary production modes providing up to 10-30% corrections on the inclusive diphoton cross section.

Lastly we note that in principle the charged scalars, which are almost degenerate with , can also be observed as they are singly produced by a similar mechanism as . However, the loop-induced decays to and both vanish in the custodial limit, rendering their -body decays dominant in almost all of parameter space. These -body decays could be probed, but such studies are likely less sensitive then other searches.

V Conclusions

We have presented a model describing a 750 GeV diphoton resonance arising from a custodial triplet which is produced as a decay product of a singly-produced VLQ. Our model has novel kinematics compared to other proposed production mechanisms and eliminates the tension from the 8 TeV diphoton searches while maintaining consistency with the kinematic distributions in the excess region. With additional statistics, our signal could be confirmed by the presence of a forward jet in the diphoton events or as a kinematic edge in the leading-jet distributions if the VLQ mass is significantly heavier than 800 GeV. The scalar resonance enjoys custodial protection, explaining the dominance of the decay rate over , , and dijet decays.

Additional signatures of the model include a corresponding excess in the and channels with fixed rates with respect to the rate. Furthermore, searches for single production of VLQs in Run II will probe deep into the viable parameter space of this model.

We now note some interesting model building possibilities which we leave for future studies. First, in this work we focus on the case where the new scalar arises from an triplet but is uncharged under . However, many of the benefits enjoyed by this model are present in similar representation choices, in particular if is a or under . These models also forbid gluon fusion and the tree-level production mechanism can dominate. Another interesting possibility is if the VLQs are related to the the top sector, as one might expect in a composite Higgs model, and the custodial symmetry is broken explicitly by the top Yukawa. Such a breaking can induce production through gluon fusion, perhaps in a controlled manner such that the decays to as well as decays to are still suppressed. Finally, we comment on some ways to further reduce the size of the Yukawas. In this work we focused on a single flavor of VLQ for simplicity. However, if there are additional flavors (with or without mixing to the other SM quarks), this can greatly enhance the diphoton decay rate, reducing the size of Yukawas necessary to reproduce the excess. An additional possibility is if SU(2) is gauged. In this case, the additional gauge bosons will propagate in the diphoton loop giving a significant enhancement to the diphoton rate.

We thank Eric Kuflik for collaboration in the early stages of this project. We are grateful to Kiel Howe and Wee Hao Ng for useful discussions. This work was supported in part by the NSF through grant PHY-1316222. JD is supported in part by the NSERC Grant PGSD3-438393-2013.

Appendix A Model details

a.1 Couplings

In this section, we derive the couplings relevant for the model between the quarks and the vector-like quarks, beginning with the couplings. The boson interactions with the up-type quarks in the interaction basis are given by (we define ) 333Note that is the upper component of an SU(2) doublet, while is the lower component of a second doublet.


Notice that we can split the coupling matrix into two pieces,


The first matrix is diagonal and commutes with the rotation to the mass basis while the second matrix yields new couplings between the VLQs and the quarks upon moving to the mass basis. Performing the rotation (the rotation matrices are given in eq. 4) gives,


Notice that the rotation to the mass basis has left the top-left entry of the coupling matrix unchanged. This is very important as it means the mixing with the VLQs does not affect the coupling which is tightly constrained experimentally. We see that we have a new coupling between the quark and the VLQ:


Now consider the -boson couplings. The right-handed up quark does not couple to the in the gauge basis, so the relevant couplings are simply:


where we moved to the mass basis in the last line.

a.2 Decay rates

In this section, we present formulae for the different decay rates used in the text for both the scalar resonance and the quark.

a.2.1 decays

We begin by considering the tree-level decays of . The dominant contribution is


where is the number of colors, , and


The other conceivable -body decays, , , and vanish identically due to the custodial production.

In addition has several loop induced decays to vector bosons as well as the Higgs. All the loop induced decays decays violate custodial symmetry. This can easily be seen at the operator level, where the terms


(where and represent the hypercharge, , and QCD field strength tensors respectively) all violate custodial symmetry and vanishes identically. There is a large breaking of this symmetry from gauging hypercharge, which induces decays into , , and . Since gauging hypercharge only breaks , the interactions are suppressed by powers of the Weinberg angle, resulting in these being generically subdominant to the photon decays. The general computation of these decay rates is made complicated due to the mixing of the VLQs with the up quark, however since these contributions are suppressed by powers of the mixing angle they are generically small. We have checked the size of these corrections by computing the rates numerically using FeynArts3, FormCalc8, and LoopTools2 Hahn:1998yk (); Hahn:2000kx () and we find that the effect is at most in the interesting region of parameter space (though often much smaller), and we neglect these effects for simplicity.

The decay rate of a scalar into two photons mediated by VLQs with mass is Gunion:1989we ()


where and (for ) . The sum runs over all VLQs and for a bidoublet the sum is


where the only non-zero contribution arises from the and quarks.

The decay to two gluons mediated by VLQs is


where is the strong coupling constant. For a bidoublet of VLQs with a triplet scalar, the sum is equal to zero showing that gluon fusion is custodially protected as expected.

The decay to mediated by VLQs is


The sum for the bidoublet is:


The decay to is


where and


and and . For the bidoublet,


The decays obey the expected relationship when they all arise from :


The decay of to is


For a bidoublet the sum vanishes identically as expected.

Lastly, the operator vanishes at tree-level and at one-loop by custodial symmetry but will be generated at two-loops by custodial symmetry breaking.

a.2.2 decays

The vector-like quarks can decay in a couple ways. We will assume such that the VLQs can decay to the scalar. Furthermore, we will focus on since that’s the only VLQ that will play a role in the phenomenology.

The decays rates are


Notice that the decay is enhanced by due to the longitudinal polarization of the . Thus in order for the decay to be substantial one needs larger Yukawas. The branching ratio into is shown in figure 8 for different Yukawas.



Figure 8: Branching ratio of the decay for different value of the couplings. The fraction is independent of the mixing angle.

a.3 Custodial symmetry breaking

In this work we have assumed that the couplings and masses of the VLQs and the triplet preserve the custodial symmetry, which enforces a cancellation in the loop amplitudes corresponding to the gluon fusion production of as well as the decays to , , and . Assuming no cancellations or large mass hierarchies, a generic scalar coupling to VLQs with coupling would acquire an effective coupling to gluons of the form


with similar expressions for the other amplitudes. The explicit breaking of custodial symmetry due to the gauging of hypercharge means that these amplitudes will still be generated, but with an additional suppression of compared to the above estimate.



Figure 9: Prototypical loop contributions to the custodial symmetry breaking amplitudes. Such two-loop contributions can induce gluon fusion production and decays two gluons. Similar diagrams can give rise to decays to and .

In particular, the VLQ mass renormalization and the renormalization of the couplings due to hypercharge gauge boson loops, illustrated in Figure 9, will contribute operators of the form . Since the mass and vertex renormalizations are logarithmically divergent, they require counterterms which are not calculable in the effective theory. Instead, we calculate the size of the IR contributions and take this as an estimate of the overall size of the irreducible contributions.

The mass renormalization of the VLQs introduces a mass splitting between the different states of size


Similarly, the vertex and wavefunction renormalization provide a contribution to the operator of size


where is the coupling constant. Now, custodial symmetry violating amplitudes of the kind in Figure 9 can be generated either with an insertion of instead of , or with the coupling . Therefore, the amplitude is suppressed by a factor


In the same spirit, one can also generate a mixing between the Higgs and the new scalar . Such a mixing is induced at two loops from the operator of the form, . The coefficient of this operator is of order


which results in a mixing angle between the Higgs and of order


where we have substituted for the mixing angle. This mixing will induce decays of to , but due to the smallness of the coupling, we do not expect this decay to be observable in the near future.

a.4 Down-type model

We now present the down-type model which can have mixing between the SM down-type quarks and the VLQs. The model is identical to the up-type model but assigning the bidoublet a U(1) charge of (as opposed to ). This gives the following fields


where , , and . As in the up-mixing case, a mixing is generated between a SM quark and a VLQ through:


where denotes the down quark in the SM mass basis of the down-type sector. The mixing produces a coupling resulting in electroweak production of , which can decay into the diphoton resonance.

Appendix B Experimental constraints



Figure 10: The constraints on the bidoublet model (reinterpreted from the work of Delaunay:2013pwa ()) arising from an ATLAS 7 TeV dedicated search for single production of VLQs ATLAS:2012apa (), a CMS 8 TeV search for /-tagged dijet resonances CMS-PAS-EXO-12-024 (), and electroweak precision (EWP). Areas above the lines are excluded. Here we neglected effects due to additional decay channels of the vector-like-quarks into the scalars in our model.

The LHC has performed searches with significant sensitivity to models with light-quark mixing. The constraints were studied in detail in Delaunay:2013pwa () for both up-quark mixing and charm-quark mixing in the context of a composite model and in Atre:2013ap () for the up-type mixing model. The dominant constraints arise from charged current production of and quarks. There are additional constraints from production of the charged quarks, but since they are always subdominant, we omit these. Instead of recasting the constraints ourselves we make use of the recast performed in Ref. Delaunay:2013pwa (). The authors recast two searches: a 7 TeV search by ATLAS searching for the bidoublet model (without the additional scalar triplet) ATLAS:2012apa () and an TeV search for excited quarks CMS:2013fea () with a similar final state which is not optimized for the single production of vector-like-quarks but shares a similar final state. The two searches have competitive limits. Additionally there are constraints on pair production of VLQs, however these are subdominant in the mass ranges we are interested in. In particular ATLAS has performed a search for VLQs decaying to finding a limit around for a single VLQ Aad:2015tba (). With two copies of such VLQs the limits strengthen but do not extend past . Furthermore, electroweak precision places an additional constraint from additional contributions to the parameter Atre:2011ae (). One might worry that the additional scalar would complicate the limits, in particular the VLQs can decay to the scalar weakening the constraints. In general these branching ratios are and we ignore such effects in our discussion.

Our goal is now to convert the single production limits quoted Ref. Delaunay:2013pwa () into our (closely related) framework. In Ref. Delaunay:2013pwa () the authors study a bidoublet model but with an additional VLQ singlet which they denote as . We can decouple the particle to match with our framework. Multiplying their cross sections by the correction factor,


with , , and gives the cross sections in our case. Employing this procedure we obtain the limits shown in figure 10.


  • (1) ATLAS collaboration, Search for resonances decaying to photon pairs in 3.2 fb of collisions at = 13 TeV with the ATLAS detector, Tech. Rep. ATLAS-CONF-2015-081, CERN, Geneva, Dec, 2015.
  • (2) CMS collaboration, Search for new physics in high mass diphoton events in proton-proton collisions at TeV, Tech. Rep. CMS-PAS-EXO-15-004, CERN, Geneva, 2015.
  • (3) R. Franceschini, G. F. Giudice, J. F. Kamenik, M. McCullough, A. Pomarol, R. Rattazzi et al., What is the gamma gamma resonance at 750 GeV?, 1512.04933.
  • (4) S. Knapen, T. Melia, M. Papucci and K. Zurek, Rays of light from the LHC, 1512.04928.
  • (5) X.-F. Han, L. Wang, L. Wu, J. M. Yang and M. Zhang, Explaining 750 GeV diphoton excess from top/bottom partner cascade decay in two-Higgs-doublet model extension, 1601.00534.
  • (6) F. P. Huang, C. S. Li, Z. L. Liu and Y. Wang, 750 GeV Diphoton Excess from Cascade Decay, 1512.06732.
  • (7) J. Liu, X.-P. Wang and W. Xue, LHC diphoton excess from colorful resonances, 1512.07885.
  • (8) A. Berlin, Diphoton and diboson excesses in a left-right symmetric theory of dark matter, Phys. Rev. D93 (2016) 055015, [1601.01381].
  • (9) J. Bernon, A. Goudelis, S. Kraml, K. Mawatari and D. Sengupta, Characterising the 750 GeV diphoton excess, 1603.03421.
  • (10) V. De Romeri, J. S. Kim, V. Martin-Lozano, K. Rolbiecki and R. R. de Austri, Confronting dark matter with the diphoton excess from a parent resonance decay, 1603.04479.
  • (11) P. Agrawal, J. Fan, B. Heidenreich, M. Reece and M. Strassler, Experimental Considerations Motivated by the Diphoton Excess at the LHC, 1512.05775.
  • (12) A. Kobakhidze, F. Wang, L. Wu, J. M. Yang and M. Zhang, 750 GeV diphoton resonance explained as a heavy scalar in top/bottom seesaw model, 1512.05585.
  • (13) P. S. B. Dev, R. N. Mohapatra and Y. Zhang, Quark Seesaw, Vectorlike Fermions and Diphoton Excess, 1512.08507.
  • (14) C. Delaunay, T. Flacke, J. Gonzalez-Fraile, S. J. Lee, G. Panico and G. Perez, Light Non-degenerate Composite Partners at the LHC, JHEP 02 (2014) 055, [1311.2072].
  • (15) K. Agashe, R. Contino, L. Da Rold and A. Pomarol, A Custodial symmetry for , Phys. Lett. B641 (2006) 62–66, [hep-ph/0605341].
  • (16) A. Atre, M. Carena, T. Han and J. Santiago, Heavy Quarks Above the Top at the Tevatron, Phys. Rev. D79 (2009) 054018, [0806.3966].
  • (17) A. Atre, G. Azuelos, M. Carena, T. Han, E. Ozcan, J. Santiago et al., Model-Independent Searches for New Quarks at the LHC, JHEP 08 (2011) 080, [1102.1987].
  • (18) M. Redi and A. Weiler, Flavor and CP Invariant Composite Higgs Models, JHEP 11 (2011) 108, [1106.6357].
  • (19) T. Flacke, J. H. Kim, S. J. Lee and S. H. Lim, Constraints on composite quark partners from Higgs searches, JHEP 05 (2014) 123, [1312.5316].
  • (20) M. Redi, V. Sanz, M. de Vries and A. Weiler, Strong Signatures of Right-Handed Compositeness, JHEP 08 (2013) 008, [1305.3818].
  • (21) N. Vignaroli, Early discovery of top partners and test of the Higgs nature, Phys. Rev. D86 (2012) 075017, [1207.0830].
  • (22) N. Vignaroli, Discovering the composite Higgs through the decay of a heavy fermion, JHEP 07 (2012) 158, [1204.0468].
  • (23) A. Atre, M. Chala and J. Santiago, Searches for New Vector Like Quarks: Higgs Channels, JHEP 05 (2013) 099, [1302.0270].
  • (24) J. F. Kamenik, B. R. Safdi, Y. Soreq and J. Zupan, ’Comments on the diphoton excess: critical reappraisal of effective field theory interpretations, 1603.06566.
  • (25) A. Alloul, N. D. Christensen, C. Degrande, C. Duhr and B. Fuks, FeynRules 2.0 - A complete toolbox for tree-level phenomenology, Comput. Phys. Commun. 185 (2014) 2250–2300, [1310.1921].
  • (26) J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer et al., The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations, JHEP 07 (2014) 079, [1405.0301].
  • (27) A. Falkowski, O. Slone and T. Volansky, Phenomenology of a 750 GeV Singlet, JHEP 02 (2016) 152, [1512.05777].
  • (28) J. Gao, H. Zhang and H. X. Zhu, Diphoton excess at 750 GeV: gluon-gluon fusion or quark-antiquark annihilation?, 1512.08478.
  • (29) C. Csáki, J. Hubisz and J. Terning, Minimal model of a diphoton resonance: Production without gluon couplings, Phys. Rev. D93 (2016) 035002, [1512.05776].
  • (30) S. Fichet, G. von Gersdorff and C. Royon, Scattering Light by Light at 750 GeV at the LHC, 1512.05751.
  • (31) C. Csaki, J. Hubisz, S. Lombardo and J. Terning, Gluon vs. Photon Production of a 750 GeV Diphoton Resonance, 1601.00638.
  • (32) ATLAS collaboration, “Diphoton searches in ATLAS.” Rencontres de Moriond, 2016.
  • (33) T. Sjostrand, S. Mrenna and P. Z. Skands, A Brief Introduction to PYTHIA 8.1, Comput. Phys. Commun. 178 (2008) 852–867, [0710.3820].
  • (34) DELPHES 3 collaboration, J. de Favereau, C. Delaere, P. Demin, A. Giammanco, V. Lemaître, A. Mertens et al., DELPHES 3, A modular framework for fast simulation of a generic collider experiment, JHEP 02 (2014) 057, [1307.6346].
  • (35) R. D. Ball et al., Parton distributions with LHC data, Nucl. Phys. B867 (2013) 244–289, [1207.1303].
  • (36) T. Hahn and M. Perez-Victoria, Automatized one loop calculations in four-dimensions and D-dimensions, Comput. Phys. Commun. 118 (1999) 153–165, [hep-ph/9807565].
  • (37) T. Hahn, Generating Feynman diagrams and amplitudes with FeynArts 3, Comput. Phys. Commun. 140 (2001) 418–431, [hep-ph/0012260].
  • (38) J. F. Gunion, H. E. Haber, G. L. Kane and S. Dawson, The Higgs Hunter’s Guide, Front. Phys. 80 (2000) 1–448.
  • (39) ATLAS collaboration, Search for Single Production of Vector-like Quarks Coupling to Light Generations in  fb of Data at TeV, .
  • (40) CMS collaboration, Search for heavy resonances in the W/Z-tagged dijet mass spectrum in pp collisions at 8 TeV, Tech. Rep. CMS-PAS-EXO-12-024, CERN, Geneva, 2013.
  • (41) CMS collaboration, CMS, Search for heavy resonances in the W/Z-tagged dijet mass spectrum in pp collisions at 8 TeV, .
  • (42) ATLAS collaboration, G. Aad et al., Search for pair production of a new heavy quark that decays into a boson and a light quark in collisions at TeV with the ATLAS detector, Phys. Rev. D92 (2015) 112007, [1509.04261].
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description