# LHC Probes of TeV-Scale Scalars in SO(10) Grand Unification. (Research Article).

1. IntroductionAfter the discovery of the Higgs boson at the Large Hadron Collider (LHC) [1,2], the last piece of the triumphant achievement of the high energy physics community, the Standard Model (SM), the great expectations for the observation of some sort of new physics at the LHC, emanated from the paradigms based on the familiar intuitions, some of which have so far lead the community to success, have turned out to be great disappointments as the LHC searches to date have returned empty-handed. Although there have been a couple of noticeable excesses, such as the diphoton [3, 4] (see [5] for a review and the full list of references) and diboson [6-8] anomalies, which caused excitement among the community, these signals have turned out to be statistical fluctuations as more data accumulates in.

While the LHC is still up and running and looking for any hint of trace pointing to physics beyond the SM (BSM), the community has been in an ambitious effort for projecting out the LHC implications of variety of new physics models for a possible future discovery. Among the various search channels, the diphoton resonance search is one of the most important programs at the LHC since this channel provides a comparatively cleaner background. One of the key predictions of many BSM theories is the existence of diphoton resonances around the TeV-scale arising from the decay of TeV-scale scalars present in those models.

One of the most appealing scenarios for a more fundamental picture is the Grand Unified Theory (GUT) framework, in which the SO(10) GUT is particularly interesting [9-24] (see [25-31] for analyses of the supersymmetric SO(10) GUT). Breaking the SO(10) gauge symmetry into that of the SM can be realized in a single step as well as in multiple steps byvarioussymmetrybreakingsequences. Therelevantoption we consider in this paper is the latter, while one possible intermediate phase, which we assume to be in the TeV-scale, is the left-right model whose gauge symmetry is based on [SU(2).sub.L] [cross product] [SU(2).sub.R] [cross product] [U(1).sub.B-L] [cross product] [SU(3).sub.C] ([G.sub.2213]) [32-39], which is different than the left-right symmetric version since in this case [SU(2).sub.L] and [SU(2).sub.R] gauge couplings are different, that is, [g.sub.L] [not equal to] [g.sub.R]. Adopting the minimalistic approach and, therefore, keeping the initial field content (the SO(10) multiplets) minimal, and tempted by the least possible fine-tuning intuition, it seems not possible to obtain a plausible scenario where the left-right model lies in the TeV-scale [23]. For instance, if the Higgs content is determined based on the extended survival hypothesis (ESH) [40], the model does not allow symmetry breaking scale of the left-right model to be in the TeV-scale. Recall that the ESH states that at every step of a symmetry breaking sequence, the only scalars which survive below the corresponding symmetry breaking scale are the ones which acquire vacuum expectation values (VEVs) at the subsequent levels of the symmetry breaking. However, by slightly relaxing the ESH conjecture by allowing one or more colored scalars to become light (at the TeV-scale), it is possible to have a TeV-scale left-right model in the SO(10) framework [23].

In this paper, we investigate the phenomenology of TeV-scale scalars as low energy remnants of the nonsupersymmetric SO(10) GUT. The part of the model that lies in the TeV-scale, as mentioned above, is the left-right model, augmented by a color-triplet scalar [[DELTA].sub.R](1, 3, 2/3, 3), whose one component we assume for our demonstration, has a mass of ~1 TeV, while its other components are heavier in the TeV range. In particular, we explore the phenomenology of a SM-singlet scalar S of mass around 1 TeV which is assumed to be the excitation of the neutral component of an [SU(2).sub.R] triplet [mathematical expression not reproducible] (1, 3, 2, 1), denoted as [mathematical expression not reproducible]. The field [mathematical expression not reproducible] breaks the symmetry of the left-right model into that of the SM by acquiring a VEV presumably at the TeV-scale in our set-up. The scalar [chi] is responsible for the production and decay of S through loop interactions.

In our model, we assume two intermediate energy scales between the electroweak scale [M.sub.z] and the unification scale At the scale the SO(10) is broken into the Pati-Salam group, [SU(2).sub.L] [cross product] [SU(2).sub.R] [cross product] [SU(4).sub.C] ([G.sub.224]). The Pati-Salam group is broken into the group of the left-right model at the first intermediate energy scale [M.sub.c], which is followed by the breaking of the left-right model into the SM at the energy scale [M.sub.R]. In our scenario, [M.sub.R] is assumed to be in the TeV-scale, while the values of [M.sub.U] and [M.sub.c] come out as predictions of the model. Note that the D-parity invariance [10, 11, 41], which is a [Z.sub.2] symmetry that maintains the complete equivalence of the left and the right sectors, is broken together with the SO(10) in the first stage of the symmetry breaking. Therefore, the gauge couplings associated with the [SU(2).sub.L] and [SU(2).sub.R] gauge groups, [g.sub.L] and evolve under the influence of different particle contents; hence [g.sub.R] [not equal to] below the scale [M.sub.U]. Remember that the D-parity is slightly different from the usual Lorentz parity in that the latter does not transform scalars, while the D-parity transforms them nontrivially. Note also that we remain in the minimal picture in terms of the total field content; the model does not have any extra matter field or any scalar SO(10) multiplet other than the ones required to begin with. Thus, the advantage of having a TeV-scale colored scalar is twofold: it is responsible for the production and decay of S and it can successfully be embedded in the minimal nonsupersymmetric SO(10) GUT scheme while maintaining the field content minimal.

In this paper, we identify the region of parameter space of our model constrained from the latest LHC data. By using a multivariate analysis (MVA), we compute the higher-luminosity LHC discovery reach of S into the diphoton channel where, as we will discuss later, the most stringent bounds come from.

The paper is organized as follows. In Section 2, we review the left-right model in the SO(10) grand unification framework. We discuss how the two scalars, S and \ of our interest, arise in our set-up. In Section 3, we discuss the unification of the couplings, derive the values of the intermediate symmetry breaking scales, and present the resulting predictions of the model. In Section 4, we present the phenomenology of S and [chi] including the exclusion limits from the LHC data and future discovery prospects. We summarize our conclusions in Section 5.

2. The Model

We consider a left-right model, whose gauge group is [SU(2).sub.L][cross product] [SU(2).sub.R] [cross product] [U(1).sub.B-L] [cross product] [SU(3).sub.C], which is assumed to be broken into the SM at the TeV-scale. The breaking is realized by the neutral component ([mathematical expression not reproducible] which we denote as S) of the [SU(2).sub.R] triplet [mathematical expression not reproducible] (1, 3, 2, 1), which is commonly preferred in the literature. Here, instead of the SU(2) triplets, the SU(2) doublet (1, 2, 1, 1), which originates from the SO(10) multiplet 16, can also be used. The advantage of the triplet representation is that it can provide a Majorana mass term for the right-handed neutrino and, hence, the seesaw mechanism [42-46] for small neutrino masses.

In this work, we explore the phenomenology of the SM-singlet S which we assume to be produced and decayed through the loop interaction with a color-triplet hyper-charged scalar [mathematical expression not reproducible] (1, 8/3, 3) denoted as % originates from the decomposition of [[DELTA].sub.R](1, 3, 2/3, 3) component of the [G.sub.224] multiplet [[DELTA].sub.R](1, 3,10) into the SM group as follows:

[mathematical expression not reproducible]. (1)

For our purpose, we take the mass of [chi] around 1 TeV, while the other components have heavier masses, ~2-5 TeV, and hence their contribution to the production and the decay of S are relatively suppressed.

The SM electroweak symmetry breaking (EWSB) in the left-right model, in general, is achieved by the neutral (diagonal) component of the bidoublet field [phi](2,2,0,1) acquiring a VEV. The fermion content of the model is the same as the SM. There are seven gauge bosons in the model, [W.sup.i.sub.L], [W.sup.i.sub.R] (with i = 1,2,3), and [W.sub.BL], with the gauge couplings and [g.sub.BL], associated with the [SU(2).sub.L], [SU(2).sub.R], and [U(1).sub.B-L] gauge symmetries, respectively. Using the notation of [21, 23], the symmetry breaking pattern of our model is given by

[mathematical expression not reproducible], (2)

where we assume [M.sub.R] = 5 TeV in our analysis.

In choosing the SO(10) multiplets for breaking the symmetries (by acquiring appropriate VEVs), we follow the common tradition in the literature as follows. The first stage of the symmetry breaking, where SO(10) is broken into the Pati-Salam group [G.sub.224], is realized by the singlet [(1,1,1).sub.210] of 210. Note that [(1,1,1).sub.210] is odd under the D-parity [10,11], and hence it is broken at this stage as well. Therefore, below the scale [M.sub.U], we have [g.sub.L] [not equal to] [g.sub.R], since they evolve under the influence of different particle contents below this energy scale according to the ESH and the minimal fine-tuning principle. The second stage, where the Pati-Salam group is broken into the left-right group [G.sub.2213], can be accomplished by [(1,1,15).sub.210] = [summation](1,1,15) acquiring a VEV. The breaking of [G.sub.2213] down to the SM gauge group [G.sub.213] is achieved by the [G.sub.2213] multiplet [mathematical expression not reproducible] which belongs to the Pati-Salam multiplet [(1, 3, 10).sub.126] [equivalent to] [[DELTA].sub.R] (1,3,10) which is a member of the SO(10) multiplet 126. In our model, [mathematical expression not reproducible] acquires a VEV at around 5 TeV which also set the value of the symmetry breaking scale MB. Note that AR is the regular [SU(2).sub.R] triplet usually used in the literature in order to break the [G.sub.2213] symmetry.

3. Unification of the Couplings

In this section, we discuss how the unification of the couplings is achieved and derive the values of the symmetry breaking scales. We have only two intermediate scales in our model in between the unification scale [M.sub.U] and the EWSB scale [M.sub.z], which are Mc and MB, where the value of MR is chosen to be 5 TeV.

The TeV-scale left-right model with light colored scalars in the minimal nonsupersymmetric SO(10) GUT scheme has recently been discussed in [23]. Here, the situation has a slight difference in one of the components in the decomposition of the left-right multiplet [mathematical expression not reproducible] (shown in (1)) into the SM gauge group, which is[mathematical expression not reproducible] whose mass is ~1TeV. Therefore, the renormalization group (RG) running of the gauge couplings at this energy scale is slightly different. The other particle S which, we assume, has a mass also around ~1 TeV, does naturally not contribute to the running since it is a SM-singlet.

3.1. Basics. We label the energy intervals in between symmetry breaking scales starting from [[M.sub.Z], [M.sub.R]] up to [[M.sub.C], [M.sub.U]] with Roman numerals as follows:

[mathematical expression not reproducible]. (3)

The boundary/matching conditions we impose on the couplings at the symmetry breaking scales are

[M.sub.u]: [g.sub.L]([M.sub.u]) = [g.sub.R]([M.sub.u]) = ([M.sub.u]) = [g.sub.4],([M.sub.u]) (4)

[M.sub.c]: [square root of 2/3] [g.sub.BL] ([M.sub.C]) = [g.sub.3] ([M.sub.C]) = [g.sub.4] ([M.sub.C]), (5)

[M.sub.R]: 1/[g.sup.2.sub.1]([M.sub.R]) = 1/[g.sup.2.sub.R]([M.sub.R]) + 1/[g.sup.2.sub.BL]([M.sub.R]), [g.sub.2]([M.sub.R]) = [g.sub.L]([M.sub.R]), (6)

MZ: 1/[e.sup.2]([M.sub.z]) = 1/[g.sup.2.sub.1]([M.sub.z]) + 1/[g.sup.2.sub.1]([M.sub.z]). (7)

The low energy data which we will use as boundary conditions to the RG running are [47,48]

[mathematical expression not reproducible], (8)

and all are evaluated at [M.sub.z] = 91.2 GeV, which gives

[mathematical expression not reproducible]. (9)

Note that the coupling constants are all required to remain in the perturbative regime during the evolution from [M.sub.U] down to [M.sub.Z].

3.2. One-Loop RG Running. For a given particle content; the gauge couplings, in an energy interval [[M.sub.A], [M.sub.B]], are evolved according to the one-loop RG relation

1/[g.sup.2.sub.i]([M.sub.A]) - 1/[g.sup.2.sub.i]([M.sub.B]) = [a.sub.i]/8[[pi].sup.2] ln [M.sub.B]/[M.sub.A], (10)

where the RG coefficients [a.sub.i] are given by [49, 50] as

[a.sub.i] = - 11/3[C.sub.2]([G.sub.i]) + 2/3[summation over[R.sub.f]][T.sub.i]([R.sub.f]) x [d.sub.1]([R.sub.f]) ... [d.sub.n] ([R.sub.f]) + [eta]/3 [summation over[R.sub.s]][T.sub.i]([R.sub.s]) ... [d.sub.n]([R.sub.s]). (11)

Here, the two summations are over irreducible chiral representations of fermions [R.sub.f] and those of scalars [R.sub.s]. The coefficient [eta] is either 1 or 1/2, depending on whether the representation is complex or real, respectively. The quadratic Casimir for the adjoint representation of the group G; is [C.sub.2]([G.sub.i]) and [T.sub.i] is the Dynkin index of each representation. For U(1) group, [C.sub.2](G) = 0 and

[summation over f.s]T = [summation over f.s] [(Y/2).sup.2], (12)

where Y/2 is the U(1) charge, the factor of 1/2 coming from the traditional normalizations of the hypercharge and B - L charges. The a;'s differ depending on the particle content in each energy interval, which changes every time symmetry breaking occurs. We will distinguish the a;'s in different intervals with the corresponding roman numeral superscript, cf. (3).

3.3. Results. The scalar content in the energy intervals are

[mathematical expression not reproducible], (13)

It is common in the literature that another scalar Pati-Salam multiplet, [??](2,2,15), is included in interval III for a rich Yukawa phenomenology [14, 15]. In terms of the RG evolution, which is our main focus here, this extra multiplet would not alter the results noticeably, because its effect in the RG equations would appear as a contribution in the term (-[5a.sub.L] + [3a.sub.R] + [2a.sub.R]) (see (14)), which would be very small compared to the rest of the term. Therefore, for the sake of staying minimal, we do not include this multiplet in our set-up.

The values of the RG coefficients for this Higgs content are listed in Table 1. The relations between symmetry breaking scales, which can be derived by using the one-loop running equations and the boundary/matching conditions, can be obtained as (for derivation see [21, 23])

[mathematical expression not reproducible], (14)

[mathematical expression not reproducible], (15)

where [s.sub.w] = sin [[theta].sub.W]. Using these equations and the experimentally measured quantities in (8) and demanding [M.sub.R]= 5 TeV, we obtain the following values:

[M.sub.c] = [10.sup.15.0] GeV, [M.sub.u] = [10.sup.17.9] GeV. (16)

The value for the scale [M.sub.c] is sufficiently high to ensure that the effects induced by the presence of scalar and vector-leptoquarks are suppressed adequately enough to remain consistent with the experimental constraints [51]. Besides, the unification scale [M.sub.u] is high enough to escape the bound on the proton decay induced by gauge boson exchanging operators. We should also note that we have light color-triplets in our model, and as well known they lead to scalar-induced dimension-6 operators that contribute to the proton decay amplitude. Although these contributions are typically suppressed by small Yukawa couplings, the color-triplets being as light as the TeV-scale can cause a potentially dangerous situation [52]. In such a case, a mechanism is required to adequately suppress these interactions, such as the ones proposed in [53, 54].

The value of the unified gauge coupling can be found via the following equation:

2[pi]/[[alpha].sub.s] - 2[pi]/[[alpha].sub.u] = [a.sup.III.sub.4] ln [M.sub.U]/[M.sub.C] + [a.sup.II.sub.3] ln [M.sub.C]/[M.sub.R] + [a.sup.I.sub.3] ln [M.sub.R]/[M.sub.Z] (17)

as [[alpha].sup.-1.sub.u] [equivalent] 47.2. The running of the couplings is given in Figure 1. Similarly, the gauge couplings at [M.sub.R] = 5 TeV are obtained as

[mathematical expression not reproducible], (18)

which, together with the values of the symmetry breaking scales in (16), are the main predictions of the model. Notice that the value of [g.sub.R](5 TeV) is different from the value of [g.sub.2](5 TeV) = [g.sub.L](5 TeV) [equivalent] 0.63, which is expected due to the fact that the D-parity invariance is broken together with the SO(10) symmetry; hence [g.sub.R] [not equal to] below the unification scale [M.sub.U], as mentioned previously. The model also predicts the existence of TeV-scale gauge bosons [W.sub.R] and [Z.sub.R] whose masses at [M.sub.R] are given as

[mathematical expression not reproducible], (19)

where we choose [mathematical expression not reproducible] which, together with (18), yields

[mathematical expression not reproducible]. (20)

These are the specific predictions of our model. However, we note that [mathematical expression not reproducible] change significantly with the choice of the symmetry breaking scale [M.sub.R]. Therefore, these mass values are not very distinctive predictions of the model. The more reliable and robust prediction is rather the values of the gauge couplings in the TeV-scale, given in (18), which do not change noticeably with the choice of the value of MR due to their logarithmic dependence on the energy scale.

Recall that our model is just the left-right model augmented by a colored scalar at the TeV-scale. Therefore, similar to the usual left-right model it allows the right-handed neutrino Nr to be Majorana in character. Although there is no mechanism that constrains right-handed neutrino mass MNr in the left-right models, there exist bounds obtained from various low energy processes [55]. The LHC implications of TeV-scale left-right models regarding a heavy Majorana right-handed neutrino for variety of mass ranges have been studied in the literature [56, 57]. As for the future runs of the LHC; as recently studied in [58], for [g.sub.R]/[g.sub.L] ~ 0.79 (which is the case in our model as can be seen in (18)), the 14 TeV LHC searches can probe the range [mathematical expression not reproducible].

4. Phenomenology

In Section 2, we have discussed that the SM-singlet S can be as light as ~1 TeV and can potentially be observed at the LHC. Since S is a SM-singlet, it can not directly couple to the SM fermions and gauge bosons through any dimension4 operator due to gauge invariance. Therefore, in order to produce S at the LHC, it is necessary to introduce extra colored particles that present in the loop. Similarly, for its decay to pair of EW gauge bosons, we need particles in the loop with nonzero hypercharge. These particles can be scalar, vector, or fermionic in nature. As mentioned previously, we would like to keep the matter and gauge sectors minimal and want to do a simplistic phenomenological study of that scenario. We, therefore, choose only one colored and hyper-charged scalar [chi] that appears naturally in our model and can serve both purposes, production and decay of S through loop interactions. Note that the EM charge of [chi] is 4/3 which is the largest among the TeV-scale colored scalars in our model. Therefore, it couples to photon with a relatively greater strength which implies large BR of S to diphoton. We further assume that [chi] is the lightest among all the colored and EM charged scalars of our model and contributes most in our analysis. We neglect any small contamination from other particles in the loop assuming that they are heavier and thus their effects are relatively suppressed. In Figure 2, we present the Feynman diagram of the production of S from gluon-gluon fusion and its decay to two photons through [chi] in the loop.

4.1. Production and Decay. The scalar S being singlet in nature, there is no tree level couplings of S to the SM fermions and gauge bosons. It can decay to a pair of SM gauge bosons only through nonrenormalizable dimension-5 operators. In the potential of the model, there could be some interaction terms which connect S with the SM Higgs doublet, which can lead to a mixing between S and the SM Higgs, after EWSB. Consequently, S can decay to a pair of the SM particles at the tree level. We know from experiments that the 125 GeV scalar observed at the LHC is very much the SM-like Higgs and therefore its mixing with S is expected to be small. For simplicity, we consider the S-h mixing, and therefore the partial widths of S to two SM fermions or two Higgs bosons are negligible. Since [chi] carries color and hypercharge, it couples to the gluon and the [B.sub.[mu]] (hypercharge) fields. Note that there is no coupling between [chi] and W bosons, since [chi] is a singlet under [SU(2).sub.L]. In the effective Lagrangian, we have the following dimension-5 operators for the interactions of S' with the SM gauge bosons prior to EWSB.

L [contains] - 1/4 [[kappa].sub.g]S[G.sup.a.sub.[mu]v][G.sup.a.[mu]v] - 1/4 [[kappa].sub.B]S[B.sub.[mu]v][B.sub.[mu]v][B.sup.[mu]v], (21)

where [G.sup.a.sub.[mu]v] and [B.sub.[mu]v], are the field-strength tensors for [SU(3).sub.C] and [U(l).sub.y] gauge groups, respectively. Effective couplings [[kappa].sub.g] and [[kappa].sub.B] are associated with the gluon and the B fields, respectively. These couplings can be computed from the knowledge of the trilinear coupling related to the S[absolute value of [chi].sub.i].sup.2] interaction term. In general, for [N.sub.j] number of colored scalars [[chi].sub.i] with hypercharge [Y.sub.i] and for an interaction term [y.sup.i.sub.S][absolute value of [chi].sup.2], the effective couplings are expressed as

[mathematical expression not reproducible], (22)

where [d.sup.i.sub.R] is the dimension of the SU(3) representation (e.g., [d.sub.R] = 3 for triplet and [d.sub.R] = 8 for octet representations) and [C.sup.i.sub.R] is the index of the SU(3) representation (e.g., [C.sub.R] = 1/2 for triplet and [C.sub.R] = 3 for octet representations). The strong and the electromagnetic couplings are denoted by as and a, respectively. The cosine of the Weinberg angle is denoted as [c.sub.w]. The loop function 70 is given by

[I.sub.0] ([tau]) = -3[tau] [1 - [tau] [[[sin.sup.-1] (1/[square root of [tau]).sup.2]] . (23)

For only one colored triplet and hypercharged (Y = 8/3) scalar [chi], [N.sub.f] = 1, [d.sub.R] = 3, and [C.sub.R] = 1/2. To keep our results as model independent as possible, we assume [y.sub.S] = [kappa][LAMBDA], where [lambda] is some new physics scale (this can be chosen as [M.sub.R]) for which we choose 5 TeV for all our computations and we keep [kappa] as a free parameter. The BRs of S to gg, [gamma] [gamma], Z[gamma], and ZZ modes are 90.6%, 5.6%, 3.3%, and 0.5% respectively.

Here, we assume that [M.sub.[chi]] > [M.sub.S]/2, and therefore S cannot decay to a [chi] pair. It is important to note that the BR depends only on Y, not on the other parameters. This is because all the partial widths, and hence the total width, scale as [[kappa].sup.2][[LAMBDA].sup.2] and the loop function 70 (for any values of [M.sub.S] and [M.sub.[chi]]) would be the same for all the partial widths. The scalar S has the largest BR in the dijet channel. We expect the BR in the [gamma][gamma], Z[gamma], and ZZ are of similar order, but ZZ mode is suppressed due to its phase space factor. The total width [[GAMMA].sub.S] is a function of [M.sub.S], [M.sub.[chi]], and [kappa][lambda]. In Figure 3, we show [[GAMMA].sub.S] as functions of [M.sub.[chi]], for [M.sub.s] = 1 TeV for three different values of [kappa] assuming [LAMBDA] = 5 TeV. As mentioned previously, [[LAMBDA].sub.S] scales as [[kappa].sup.2][[LAMBDA].sup.2] and one can easily estimate the total width for other values of [kappa][LAMBDA] from this plot.

4.2. Exclusions from LHC Data. To derive bounds on the model parameters from the LHC data and related numerical analysis, we implement the Lagrangian given in (21) in FeynRules 2.0 [59] to generate the model files for the MadGraph5 [60] event generator. We use CTEQ6L1 [61] parton distribution functions (PDF) to compute cross sections. We fix the factorization and renormalization scales at Ms for all our numerical computations.

For our phenomenological analysis, we have only three free parameters, namely, [M.sub.S], [M.sub.[chi]], and [kappa] (we choose [LAMBDA] = 5 TeV for all our numerical computations). We first derive bounds on the parameters from the latest LHC 13 TeV yy [62, 63], Z[gamma] [64, 65], ZZ [66], and jj [67, 68] resonance search data. The observed upper limit (UL) at 95% confidence level (CL) on the cross sections for the resonance mass of 1 TeV of four type of resonances is given by

[mathematical expression not reproducible]. (24)

These values are used in Figure 4(a) where we show the excluded parameter space (colored regions) in [M.sub.[chi]]-[kappa] plane for [M.sub.S] = 1 TeV. The excluded regions shown in orange, green, blue, and brown are derived from the [gamma][gamma], Z[gamma], ZZ, and jj resonance search data. We can see that the diphoton data is the most powerful in constraining the parameter space in [M.sub.[chi]]-[kappa] plane. In Figure 4(b), we present the excluded regions in [M.sub.S]-[M.sub.[chi]] plane for different [kappa] with [LAMBDA] = 5 TeV from the latest 13 TeV combined ATLAS and CMS diphoton resonance search data. Cross section ULs ([[sigma].sub.i]) from different experiments and the corresponding uncertainties (Aa;) are combined statistically using the following relations:

1/[([DELTA][[sigma].sub.c]).sup.2] = [summation over i]1/[([DELTA][[sigma].sub.i]).sup.2], [[sigma].sub.c])/[([DELTA][[sigma].sub.i]).sup.2] = [summation over i][[sigma].sub.i]/[([DELTA][[sigma].sub.i]).sup.2], (25)

where [[sigma].sub.c] is the combined cross section and [DELTA][[sigma].sub.c] is the uncertainty associated with it. In case of asymmetric uncertainties, we get [DELTA][[sigma].sub.i] by averaging upper and lower uncertainties. Although uncertainties are used to compute [[sigma].sub.c], we have not shown the uncertainty bands in the exclusion plots for simplicity. The sky-blue regions in these plots cannot be probed in our set-up as we always assume [M.sub.[chi]] > [M.sub.S]/2. If [M.sub.[chi]] < [M.sub.[chi]]/2, the S [right arrow] [chi][chi] decay becomes kinematically allowed and becomes the dominant decay mode of S. This will make the diphoton and other branching modes suppressed. Therefore, observing S in the [gamma][gamma], Z[gamma], ZZ, and jj resonance searches becomes much more challenging. One should note that exclusion regions are not very sensitive to the [M.sub[chi]] values for a fixed k. This is because, for heavier resonances, the cross section ULs are not very sensitive to the resonance mass due to lack of statistics and therefore the quantity [sigma]([M.sub.S], [M.sub.[chi]]) x BR should remain insensitive for heavier resonances. The reduction in the production cross section as we increase [M.sub.S] is compensated by the slight change in [M.sub.chi]] since [sigma]([M.sub.S], [M.sub.[chi]] quantity is very sensitive to the [M.sub.chi]]. Note that these bounds are derived from the observed 95% CL UL on the cross sections. Consideration of uncertainties on the cross sections limits would relax the derived bound somewhat. In all our computations, we have considered a next-to-leading order Kfactor of 2 to account for the higher-order effects [69].

4.3. Future Prospects at the LHC. In this subsection, we look at the prospect of discovering S at the 13 TeV LHC runs with high integrated luminosities. In previous subsection, we find that the most stringent bounds come from the diphoton data. Therefore, we only focus on the diphoton final state for the present prospect study. After event generation, we use Pythia6 [70] for parton shower and hadronization. The subsequent detector simulation is done using Delphes3 [71] package. Jets are clustered with FastJet [72] using the anti[[kappa].sub.T]- algorithm [73] with the clustering parameter, R = 0.4. We use TMVA [74] for the multivariate analysis.

Signal events are generated with up to two jets, that is, pp [right arrow] S([right arrow] [gamma][gamma]) + 0,1,2 jets which are MLM [75] merged at a matching scale [Q.sub.cut] = 50 GeV. The dominant (roughly 90%) SM background for this signal comes from the qq - [gamma][gamma] process. Similar to the signal, we generate this background by merging pp - [gamma][gamma] + 0,1,2 jets processes at [Q.sub.cut] = 15 GeV. We only consider this dominant background in our analysis. Appropriate matching scales for signal and background are determined by assuring smooth transition in the differential jet-rate distributions between events with N and N +1 jets and matched cross sections are within ~10% of the zero jet contribution. We also check the stability of the matched cross section with the variation of [Q.sub.cu]t once it is properly chosen.

The 13 TeV diphoton data already set an UL on [sigma] x BR ~ 1 fb for the resonance mass of around 1 TeV. Therefore, it is very challenging to observe such a signal over the large SM background. ATLAS and CMS collaborations use cut-based technique in their diphoton resonance searches at the 13 TeV LHC. In this paper, to obtain better sensitivity, we use a MVA to discriminate tiny signal from the large SM background. ATLAS and CMS ULs on a x BR slightly depend on the width of the resonance but we use a fixed width of 1 GeV for all [M.sub.S] in the following analysis for simplicity. The width of S is a function of model parameters, namely, [M.sub.S], [M.sub.[chi]], and [kappa][LAMBDA]. Instead of choosing a specific benchmark, we use [[GAMMA].sub.S] = 1 GeV for our MVA. This analysis is insensitive to the actual width choice as long as [[GAMMA].sub.S] [much less than] [M.sub.S]; that is, the narrow width approximation is well-valid.

We generate signal and background events with some basic transverse momentum ([p.sub.T]), pseudorapidity ([eta]) and separation in [eta]-[phi] plane ([DELTA]R) cuts as follows:

[mathematical expression not reproducible], (26)

where x, y = {y, j}. We use a strong selection cut on the invariant mass of the photon pair, [absolute value of M([gamma][gamma])- [M.sub.S]] < 100 GeVto reduce the huge diphoton background before passing events to TMVA. For MVA, we use the Boosted Decision Tree (BDT) algorithm where we feed the following seven kinematic variables: [P.sub.T]([[gamma].sub.1]), [P.sub.T]([[gamma].sub.2])[absolute value of [eta]([[gamma].sub.1])], [absolute value of [eta]([[gamma].sub.2])],[DELTA]R([[gamma].sub.1], [[gamma].sub.2], M([[gamma].sub.1], [[gamma].sub.2]), and jet multiplicity ([[gamma].sub.1] and [[gamma].sub.2], are the two selected photons ordered according to their [p.sub.T]). In Figure 5, we show the signal (blue) and background (red) distributions of these variables used in MVA. We choose these simple variables which are less correlated and have sufficiently good discriminating power. In Table 2, we show the relative importance (RI) of these variables for the benchmark mass [M.sub.S] = 1 TeV. We find that the two variables M([[gamma].sub.1], [[gamma].sub.2]) and [DELTA]R([[gamma].sub.1], [[gamma].sub.2]) are very effective in discriminating signal from background. Other variables like [p.sub.T] and [eta] of photons also have reasonably good discriminating power. We obtain the cut efficiency of almost 75% for the signal but as small as 10% for the background for the whole range of MS we considered. It is important to mention that this set of seven variables used might not be the optimal one. There is always a scope to improve the analysis with cleverer choices of variables.

The BDT algorithm is prone to overtraining and therefore one should always be careful while using it in MVA. Overtraining of the signal and background test samples can usually happen due to the improper choices of BDT tuning parameters. Whether a test sample is overtrained or not can be checked by using the Kolmogorov-Smirnov (KS) statistics. Generally, a test sample is not overtrained if the corresponding KS probability lies within the range 0.1 to 0.9. In our analysis, we use two statistically independent samples for each [M.sub.S] choice, one for training and the other for testing the BDT. In Figure 6(a), we show the BDT response of the signal and background for the benchmark mass [M.sub.S] = 1 TeV. From the BDT response, one can see that a BDT cut around ~0 can effectively separate the signal from the background and lead to best significance. In Figure 6(b), we show the M([[gamma].sub.1], [[gamma].sub.2]) distributions for the signal and the background for [M.sub.S] = 1 TeV at the 13 TeV LHC with L = 300 [fb.sup.-1]. This plot is shown for the significance of 5[sigma] where [sigma] = [N.sub.s]/[square root of [N.sub.s] + [N.sub.B]] and the numbers of signal and background events that survive after the optimal BDT cut (>0) are [N.sub.s] = 69 with cut efficiency 0.75 and [N.sub.B] = 120 with cut efficiency 0.1, respectively.

In Figures 7(a) and 7(b), we show the 5[sigma] discovery contours in [M.sub.S]-[M.sub.[chi]] plane for different [kappa] at the 13 TeV LHC for 100 and 300 [fb.sub.-1] integrated luminosities, respectively. As stated earlier, the sky-blue region, that is, [M.sub.[chi]] < [M.sub.S]/2, is not considered in our analysis. We observe that the discovery reach for 100 [fb.sub.-1] run in Figure 7(a) is not much improved from the bounds obtained in Figure 4. But for 300 [fb.sub.-1] run, a substantially bigger region of parameter space can be probed. In Figure 7(c), we show the expected 95% CL exclusion plot in [M.sub.S]-[M.sub.[chi]] plane for L = 300 [fb.sub.-1]. It is obvious that the parameter space which can be excluded with 95% CL is much bigger than the parameter space which can be discovered with 5[sigma] significance. As previously mentioned, the limits on [sigma] x BR for scalar decay to diphoton are already very strict. Therefore, to observe such a scalar at the LHC is very challenging and we need a more dedicated analysis for that.

In this paper, we choose to use a MVA for the LHC prospect study to achieve better sensitivity to the parameter space compared to a cut-based analysis. To give the readers a rough idea of gain in sensitivity, we wish to present here a quantitative comparison between the two types of analyses for the benchmark mass [M.sub.S] = 1 TeV. We apply further the following hard cuts on photons, namely, [p.sub.T]([[gamma].sub.1]), [p.sub.T]([[gamma].sub.2]) > 200 GeV and [absolute value of M([[gamma].sub.1], [[gamma].sub.2]) - [M.sub.S]] < 50 GeV on the events that are used for the BDT analysis. In context of Figure 6(b), we have discussed previously that the numbers of signal and background events which survive after the optimal BDT cut (around ~0) are 69 and 120, respectively. The corresponding signal and background events that survive after the cut-based analysis are 65 and 432, respectively, which leads to a ~3a significance. One can see, therefore, that BDT analysis is very effective in terms of background reduction compared to a cut-based analysis. Note that this set of cuts is not fully optimized (but fairly good) and one can vary these cuts to find the optimized set of cuts to improve the significance from ~3[sigma]. But an optimized BDT analysis is always expected to perform better than an optimized cut-based analysis as long as a clever set of variables are used. A BDT analysis is usually more effective than a cut-based analysis especially in the low mass (here low [M.sub.S]) region. For heavier masses, where the SM background is expected to be very small compared to the signal, an optimized cut-based analysis can compete to an optimized BDT analysis.

5. Summary

In this paper, we explore the phenomenology of TeV-scale scalars in the nonsupersymmetric SO(10) grand unification framework. In particular, we investigate the LHC phenomenology of a SM-singlet scalar S which interacts with gluons and photons through loop interactions with a colortriplet hypercharged scalar [chi] which is remnant from the breaking of the Pati-Salam gauge group, [SU(2).sub.L] [cross product] [SU(2)sub.R] [cross product] [SU(4).sub.C]. The part of the model that lies in the TeV-scale is the left-right model, whose gauge group is [SU(2).sub.L] [cross product] [SU(2)sub.R] [cross product] [U(1).sub.B-L] [cross product] [SU(3).sub.C], augmented with the color-triplet scalar The scalar S is a component of an [SU(2)sub.R] triplet scalar which is responsible for the breaking of the left-right model into the SM. Note that we have stayed in the minimal picture in terms of the total field content; the model does not have any extra matter fields or any SO(10) multiplets in the scalar content other than the ones required to begin with.

The colored scalar in our set-up effectively induces the interaction terms of S with gluons and photons that lead to a diphoton final state after being produced via gluon fusion. In addition to the [gamma][gamma] decay, S can also decay to jj, [gamma]Z, and ZZ modes. We present the exclusion region in [M.sub.[chi]]-[kappa] plane for a benchmark resonance mass [M.sub.S] = 1 TeV using the latest LHC data. We find that the most stringent bounds on the parameter space of our model come from the diphoton resonance search data. Therefore, we consider the diphoton channel as the most promising channel for the discovery of S at the LHC. As a prospect study, we compute the higher-luminosity LHC discovery reach of S by using a state-of-the-art multivariate technique. We present 5[sigma] discovery contours for different [kappa] choices in the [M.sub.S]-[M.sub.[chi]] plane at the 13 TeV LHC with 100 and 300 [fb.sup.-1] integrated luminosity. From our analysis, we find that for [kappa] ~1, [M.sub.S] ~ 0.5-2 TeV and [M.sub.[chi]] ~ 1 TeV can easily be observed with 5[sigma] confidence level at the 13 TeV LHC with 300 [fb.sup.-1] integrated luminosity. Note that the role of various systematic uncertainties is always important to consider in an analysis for robust and accurate prediction. But in the current scope, we do not consider systematic uncertainties for simplicity.

The unification of the couplings in the model is successfully realized, where the TeV-scale colored triplet plays an important role. As discussed in [23], it is very difficult to achieve a successful SO(10) grand unification set-up with a TeV-scale left-right model. Slightly modifying the low energy scalar content by allowing a number of colored scalars, originated from various Pati-Salam multiplets, to become light generates the possibility of accommodating a TeV-scale left-right model in the SO(10) grand unification framework. Among a number of low energy scalar configurations, the ones with the very color-triplet selected in our model appear to particularly stand out [23]. We also note that the values obtained for the intermediate scale (where the Pati-Salam is broken) and the unification scale are sufficiently high to remain compatible with the experimental constraints regarding the leptoquark induced effects and the proton decay.

https://doi.org/10.1155/2017/7498795

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors would like to thank Ilia Gogoladze for stimulating conversations. This work is supported by the Swedish Research Council under Contract 621-2011-5107. Tanumoy Mandal is supported by the Carl Trygger Foundation under Contract CTS-14:206.

References

[1] S. Chatrchyan et al., "Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC," Physics Letters B, vol. 716, pp. 30-61, 2012.

[2] G. Aad et al., "Observation of a new particle in the search for the standard model higgs boson with the ATLAS detector at the LHC," Physics Letters B, vol. 716, pp. 1-29, 2012.

[3] M. Aaboud et al., "Search for resonances in diphoton events at Vs = 13 TeV with the ATLAS detector," Journal of High Energy Physics, vol. 9, no. 1, 2016.

[4] V. Khachatryan et al., "Search for resonant production of high-mass photon pairs in proton-proton collisions at [square root of s] = 8 and 13 TeV," Physical Review Letters, vol. 117, no. 5, Article ID 051802, 2016.

[5] A. Strumia, "Interpreting the 750 GeV digamma excess: a review," 2016, arXiv:1605.09401.

[6] G. Aad et al., "Search for high-mass diboson resonances with boson-taggedjetsinproton-proton collisions at [square root of s] = 8 TeV with the ATLAS detector," Journal of High Energy Physics, vol. 12, no. 55, 2015.

[7] G. Aad et al., "Search for resonant diboson production in the ll[bar.qq.] final state in pp collisions at [square root of s] 8 TeV with the ATLAS detector," The European Physical Journal C, vol. 75, no. 69,2015.

[8] G. Aad et al., "Search for production of WW/WZ resonances decaying to a lepton, neutrino and jets in pp collisions at [square root of 8] TeV with the ATLAS detector," The European Physical Journal C, vol. 75, no. 5, p. 209, 2015.

[9] D. Chang, R. N. Mohapatra, and M. K. Parida, "Decoupling of parity- and SU(2)E-breaking scales: a new approach to left-right symmetric models," Physical Review Letters, vol. 52, no. 13, pp. 1072-1075, 1984.

[10] D. Chang, R. N. Mohapatra, and M. K. Parida, "New approach to left-right-symmetry breaking in unified gauge theories," Physical Review D, vol. 30, no. 5, article 1052,1984.

[11] D. Chang, R. N. Mohapatra, J. M. Gipson, R. E. Marshak, and M. K. Parida, "Experimental tests of new SO(10) grand unification," Physical Review D, vol. 31, no. 7, pp. 1718-1732,1985.

[12] M. K. Parida and C. C. Hazra, "Superheavy-Higgs-scalar effects in effective gauge theories from SO(10) grand unification with low-mass right-handed gauge bosons," Physical Review D, vol. 40, no. 9, pp. 3074-3085, 1989.

[13] N. G. Deshpande, E. Keith, and P. B. Pal, "Implications of the CERN LEP results for SO(10) grand unification," Physical Review D, vol. 46, no. 5, pp. 2261-2264,1992.

[14] B. Bajc, A. Melfo, G. Senjanovic, and F. Vissani, "Yukawa sector in nonsupersymmetric renormalizable SO(10)," Physical Review D--Particles, Fields, Gravitation and Cosmology, vol. 73, no. 5, Article ID 055001, 2006.

[15] S. Bertolini, L. Di Luzio, and M. Malinsky, "Intermediate mass scales in the nonsupersymmetric SO(10) grand unification: A reappraisal," Physical Review D--Particles, Fields, Gravitation and Cosmology, vol. 80, no. 1, Article ID 015013, 2009.

[16] K. S. Babu and R. N. Mohapatra, "Coupling unification, GUT-scale baryogenesis and neutron-antineutron oscillation in SO(10)," Physics Letters, 2012.

[17] R. L. Awasthi, M. K. Parida, and S. Patra, "Neutrino masses, dominant neutrinoless double beta decay, and observable lepton flavor violation in left-right models and SO(10) grand unification with low mass WR, ZR bosons," Journal of High Energy Physics, vol. 2013, no. 8, article 122, 2013.

[18] B. P. Nayak and M. K. Parida, "New mechanism for Type-II seesaw dominance in SO(10) with low-mass Z/,RH neutrinos, and verifiable LFV, LNV and proton decay," European Physical Journal C, vol. 75, article no. 183, 2015.

[19] M. K. Parida, R. L. Awasthi, and P. K. Sahu, "Proton decay and new contribution to 0v2ss decay in SO(10) with low-mass Z' boson, observable n-n oscillation, lepton flavor violation, and rare kaon decay," Journal of High Energy Physics, vol. 1, no. 45, 2015.

[20] T. D. Brennan, "Two loop unification of non-SUSY SO(10) GUT with TeV scalars," Physical Review D, vol. 95, no. 6, Article ID 065008, 2017.

[21] U. Aydemir, D. Minic, C. Sun, and T. Takeuchi, "Pati-salam unification from noncommutative geometry and the TeV-scale WR boson," International Journal of Modern Physics A, vol. 31, no. 1, Article ID 1550223, 2016.

[22] T. Bandyopadhyay, B. Brahmachari, and A. Raychaudhuri, "Implications of the CMS search for WR on grand unification," Journal of High Energy Physics, vol. 2, no. 23, 2016.

[23] U. Aydemir, "SO(10) grand unification in light of recent LHC searches and colored scalars at the TeV-scale," International Journal of Modern Physics A, vol. 31, no. 8, Article ID 1650034, 2016.

[24] M. K. Parida, B. P. Nayak, R. Satpathy, and R. L. Awasthi, "Standard coupling unification in SO(10), hybrid seesaw neutrino mass and leptogenesis, dark matter, and proton lifetime predictions," Journal of High Energy Physics, vol. 2017, no. 75,58 pages, 2017, arXiv:1608.03956.

[25] N. G. Deshpande, E. Keith, and T. G. Rizzo, "SO(10) grand unification with a low-energy SU(2)E-symmetry-breaking scale Me," Physical Review Letters, vol. 70, no. 21, pp. 3189-3192,1993.

[26] T. Fukuyama, A. Ilakovac, T. Kikuchi, S. Meljanac, and N. Okada, "General formulation for proton decay rate in minimal supersymmetric SO(10) GUT," The European Physical Journal C, vol. 42, pp. 191-203, 2005.

[27] T. Fukuyama, A. Ilakovac, T. Kikuchi, S. Meljanac, and N. Okada, "SO(10) group theory for the unified model building," Journal of Mathematical Physics, vol. 46, no. 3, Article ID 033505, 2005.

[28] S. K. Majee, M. K. Parida, A. Raychaudhuri, and U. Sarkar, "Low intermediate scales for leptogenesis in supersymmetric SO(10) grand unified theories," Physical Review D, vol. 75, no. 7, Article ID 075003, 2007.

[29] M. K. Parida, "Intermediate left-right gauge symmetry, unification of couplings and fermion masses in SUSY SO(10)XS(4)," Physical Review D, vol. 78, no. 5, Article ID 053004, 2008.

[30] P. S. B. Dev and R. N. Mohapatra, "TeV scale inverse seesaw model in SO(10) and leptonic nonunitarity effects," Physical Review D, vol. 81, Article ID 013001, 2010.

[31] M. K. Parida and A. Raychaudhuri, "Inverse seesaw mechanism, leptogenesis, observable proton decay, and [[DELTA].sup.[+ or-][+ or-].sub.R] in supersymmetric SO(10) with heavy [W.sub.R]," Physical Review D, vol. 82, no. 9, Article ID 093017, 2010.

[32] J. C. Pati and A. Salam, "Erratum: Lepton number as the fourth color"' Physical Review D, vol. 10, pp. 275-289,1975.

[33] R. N. Mohapatra and J. C. Pati, "'Natural' left-right symmetry," Physical Review D, vol. 11, no. 9, pp. 2558-2561,1975.

[34] G. Senjanovic and R. N. Mohapatra, "Exact left-right symmetry and spontaneous violation of parity," Physical Review D, vol. 12, no. 5, pp. 1502-1505,1975.

[35] R. N. Mohapatra and G. Senjanovic, "Neutrino mass and spontaneous parity nonconservation," Physical Review Letters, vol. 44, p. 912,1980.

[36] R. N. Mohapatra and G. Senjanovic, "Neutrino masses and mixings in gauge models with spontaneous parity violation," Physical Review D, vol. 23, article 165,1981.

[37] P. Duka, J. Gluza, and M. Zralek, "Quantization and renormalization of the manifest left-right symmetric model of electroweak interactions," Annals of Physics, vol. 280, no. 2, pp. 336-408, 2000.

[38] U. Aydemir, D. Minic, and T. Takeuchi, "The Higgs mass and the emergence of new physics," Physics Letters, Section B: Nuclear, Elementary Particle and High-Energy Physics, vol. 724, no. 4-5, pp. 301-305, 2013.

[39] U. Aydemir, D. Minic, C. Sun, and T. Takeuchi, "Higgs mass, superconnections, and the TeV-scale left-right symmetric model," Physical Review D--Particles, Fields, Gravitation and Cosmology, vol. 91, no. 4, Article ID 045020, 2015.

[40] F. Del Aguila and L. E. Ibanez, "Higgs bosons in SO(10) and partial unification," Nuclear Physics, Section B, vol. 177, no. 1, pp. 60-86,1981.

[41] A. Maiezza, M. Nemevsek, F. Nesti, and G. Senjanovic, "Left-right symmetry at LHC," Physical Review D--Particles, Fields, Gravitation and Cosmology, vol. 82, no. 5, Article ID 055022, 2010.

[42] P. Minkowski, "[mu] [right arrow] e[gamma] at a rate of one out of [10.sub.9] muon decays?" Physics Letters B, vol. 67, pp. 421-428,1977.

[43] M. Gell-Mann, P. Ramond, and S. Richard, "Complex spinors and unified theories," Conf. Proc, Article ID C790927, pp. 315-321,1979, https://arxiv.org/abs/1306.4669.

[44] T. Yanagida, "Horizontal symmetry and masses of neutrinos," Conf. Proc, Article ID C7902131, pp. 95-99,1979.

[45] J. Schechter and J. W. F. Valle, "Neutrino masses in SU(2) [cross product] U(1) theories," Physical Review D, vol. 22, p. 2227,1980.

[46] J. Schechter and J. W. F. Valle, "Neutrino decay and spontaneous violation of lepton number," Physical Review D, vol. 25, p. 774, 1982.

[47] C. Patrignani et al., "Review of particle physics," Chinese Physics C, vol. 40, no. 10, Article ID 100001, 2016.

[48] S. Schael et al., "Precision electroweak measurements on the Z resonance," Physics Reports, vol. 427, no. 5-6, pp. 254-457,2006.

[49] D. R. T. Jones, "The two loop beta function for a G(1) x G(2) gauge theory," Physical Review D, vol. 25, no. 581,1982.

[50] M. Lindner and M. Weiser, "Gauge coupling unification in left-right symmetric models," Physics Letters. B. Particle Physics, Nuclear Physics and Cosmology, vol. 383, no. 4, pp. 405-414, 1996.

[51] J. L. Evans, N. Nagata, K. A. Olive, and J. Zheng, "The ATLAS diboson resonance in non-supersymmetric SO(10)," Journal of High Energy Physics, vol. 2016, no. 2, article no. 120, pp. 1-25, 2016.

[52] S. Bertolini, L. Di Luzio, and M. Malinsky, "Seesaw scale in the minimal renormalizable SO(10) grand unification," Physical Review D--Particles, Fields, Gravitation and Cosmology, vol. 85, no. 9, Article ID 095014, 2012.

[53] G. R. Dvali, "Light colour-triplet Higgs is compatible with proton stability: an alternative approach to the doublet-triplet splitting problem," Physics Letters B, vol. 372, no. 1-2, pp. 113-120, 1996.

[54] S. Rakshit, G. Raz, S. Roy, and Y. Shadmi, "Light grand unified theory triplets and Yukawa splitting," Physical Review D, vol. 69, no. 9, Article ID 095006, 2004.

[55] V. Tello, M. Nemevek, F. Nesti, G. Senjanovi C, and F. Vissani, "Left-right symmetry: from the LHC to neutrinoless double beta decay," Physical Review Letters, vol. 106, no. 15, Article ID 151801, 2011.

[56] R. N. Mohapatra, "Neutrino mass as a signal of TeV scale physics," Nuclear Physics B, vol. 908, pp. 423-435, 2016.

[57] J. C. Helo, M. Hirsch, and S. Kovalenko, "Erratum: Heavy neutrino searches at the LHC with displaced vertices [Phys. Rev. D 89, 073005 (2014)]," Physical Review D, vol. 93, no. 9, Article ID 099902, 2016.

[58] R. Ruiz, Lepton number violation at colliders from kinematically inaccessible gauge bosons, arXiv:1703.04669 hep-ph.

[59] A. Alloul, N. D. Christensen, C. Degrande, C. Duhr, and B. Fuks, "FeynRules 2.0--a complete toolbox for tree-level phenomenology," Computer Physics Communications, vol. 185, no. 8, pp. 2250-2300, 2014.

[60] J. Alwall, R. Frederix, and S. Frixione, "The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations," Journal of High Energy Physics, vol. 2014, no. 7, article 079,2014.

[61] J. Pumplin, D. R. Stump, J. Huston, H.-L. Lai, P. Nadolsky, and W.-K. Tung, "New generation of Parton distributions with uncertainties from global QCD analysis," Journal of High Energy Physics, vol. 7, no. 12, 2002.

[62] The ATLAS collaboration, "Search for scalar diphoton resonances with 15.4 [fb.sup.-1] of data collected at [square root of s] = 13 TeV in 2015 and 2016 with the ATLAS detector," ATLAS-CONF-2016-059, 2016.

[63] V. Khachatryan et al., "Search for high-mass diphoton resonances in proton-proton collisions at 13 TeV and combination with 8 TeV search," Physics Letters B, vol. 767, pp. 147-170, 2017.

[64] CMS Collaboration, "Search for high-mass resonances in Z[gamma] [right arrow] [e.sup.+][e.sup.-][gamma]/[[mu].sup.+][[mu].sup.-][gamma] final states in proton-proton collisions at [square root of s ] = 13 TeV"' CMS-PAS-EX0-16-034.

[65] CMS Collaboration, "Search for high-mass resonances in Z(qq)[gamma] final state in pp collisions at [square root of s] = 13 TeV with 12.9 [fb.sup.-1]," CMS-PAS-EX0-16-035.

[66] ATLAS collaboration, "Searches for heavy ZZ and ZW resonances in the llqq and wqq final states in pp collisions at [square root of s] = 13 TeV with the ATLAS detector," ATLAS-CONF-2016-082, 2016.

[67] The ATLAS collaboration, "Search for New Phenomena in Dijet Events with the ATLAS Detector at [square root of s] =13 TeV with 2015 and 2016 data," ATLAS-CONF-2016-069, 2016.

[68] CMS Collaboration, "Searches for narrow resonances decaying to dijets in proton-proton collisions at 13 TeV using 12.9 inverse femtobarns," CMS-PAS-EX0-16-032.

[69] D. de Florian et al., "Handbook of lhc higgs cross sections: 4. deciphering the nature of the higgs sector," Tech. Rep., 2016, arXiv:1610.07922 hep-ph.

[70] T. Sjostrand, S. Mrenna, and P. Skands, "PYTHIA 6.4 physics and manual," Journal of High Energy Physics, vol. 5, article 026, 2006.

[71] J. de Favereau, C. Delaere, P. Demin et al., "DELPHES 3: a modular framework for fast simulation of a generic collider experiment," Journal of High Energy Physics, vol. 2014, article 57, 2014.

[72] M. Cacciari, G. P. Salam, and G. Soyez, "FastJet user manual: (For version 3.0.2)," European Physical Journal C, vol. 72, no. 3, article 1896, pp. 1-54, 2012.

[73] M. Cacciari, G. P. Salam, and G. Soyez, "The Anti-k(t) jet clustering algorithm," Journal of High Energy Physics, vol. 2008, article 063, 2008.

[74] A. Hocker et al., "TMVA--toolkit for multivariate data analysis," in Proceedings of PoS ACAT, vol. 40, 2007.

[75] M. L. Mangano, M. Moretti, F. Piccinini, and M. Treccani, "Matching matrix elements and shower evolution for top-pair production in hadronic collisions," Journal of High Energy Physics, vol. 2007, no. 1, article no. 013, 2007.

Ufuk Aydemir and Tanumoy Mandal

Department of Physics and Astronomy, Uppsala University, P.O. Box 516, 75120 Uppsala, Sweden

Correspondence should be addressed to Tanumoy Mandal; tanumoy.mandal@physics.uu.se

Received 14 February 2017; Accepted 24 April 2017; Published 30 May 2017

Academic Editor: Anna Cimmino

Caption: Figure 1: Running of the gauge couplings for the model. The vertical dotted lines from left to right correspond to the symmetry breaking scales [M.sub.z], [M.sub.R], and [M.sub.c], which also indicate the beginning of the energy intervals I, II, and III, respectively. For [[alpha].sup.-1.sub.1] and [[alpha].sup.-1.sub.BL] we plot the redefined quantities [[??].sup.-1.sub.1] [equivalent to] (3/5) [[alpha].sup.-1.sub.1] and [[??].sup.-1.sub.BL] = (3/2) [[alpha].sup.-1.sub.BL]. Note that the discontinuity on the [[??].sup.-1.sub.BL] plot at the energy scale [M.sub.R] occurs due to the boundary condition given in (6).

Caption: Figure 2: The Feynman diagram of the production and decay of S at the LHC through [chi] in the loop.

Caption: Figure 3: Total width of S for [M.sub.S] = 1 TeV as functions of [M.sub.[chi]] for [kappa] = 2,4,6 assuming [LAMBDA] = 5 TeV.

Caption: Figure 4: (a) The excluded region in [M.sub.[chi]]-[kappa] plane for [M.sub.S] = 1 TeV from the 13 TeV LHC data. The orange, green, blue, and brown regions represent the region ruled out by the [gamma][gamma], Z[gamma], ZZ, and jj resonance search data, respectively. (b) The excluded region in [M.sub.S]-[M.sub.[chi]] plane for different [kappa] with [LAMBDA] = 5 TeV using the 13 TeV combined ATLAS and CMS diphoton resonance search. ATLAS and CMS data are combined statistically using (25). The widths of the resonance assumed by ATLAS and CMS are 4MeV [62] and 140 MeV [63]. The sky-blue region cannot be probed in our analysis as we assume [M.sub.[chi]] > [M.sub.S]/2.

Caption: Figure 5: Signal (blue) and background (red) distributions of the input variables used for MVA. These distributions are generated for [M.sub.S] = 1 TeV assuming [[GAMMA].sub.s] = 1 GeV. We do not present the jet multiplicity distribution here since it has the smallest RI in the MVA as shown in Table 2 and therefore the [N.sub.jet] distribution would not differ much for the signal and the background.

Caption: Figure 6: (a) BDT response of the signal (blue) and the background (red) for [M.sub.S] = 1 TeV. (b) M([[gamma].sub.1], [[gamma].sub.2]) distributions for the signal (blue) and the background (red) for [M.sub.S] = 1 TeV after applying the optimal BDT cut at ~0 to obtain 5[sigma] significance at the 13 TeV LHC with L = 300 [fb.sup.-1].

Caption: Figure 7: The 5[sigma] discovery contours in [M.sub.S]-[M.sub.chi]] plane for different k at 13 TeV LHC for (a) 100 [fb.sup.-1] and (b) 300 [fb.sup.-1] integrated luminosity. (c) The expected 2[sigma] exclusion plot for L = 300 [fb.sup.-1]. The sky-blue region cannot be probed in our analysis as we assume [M.sub.[chi]] > [M.sub.S]/2.

Table 1: The Higgs content and the RG coefficients in the energy intervals for our model. Interval Higgs content III [phi](2,2,1), [[DELTA].sub.R](1,3,10), [summation](1,1,15) II [mathematical expression not reproducible] I [phi](2,1,1), S(1,1,1), [chi](1,8/3,3) Interval RG coefficients III [([a.sub.L],[a.sub.R],[a.sub.4]).sup.III] = (-3,11/3,-7)) II [([a.sub.L],[a.sub.R],[a.sub.BL],[a.sub.3]).sup.II] = (-3,-1/3,4,-13/2)) I [([a.sub.1],[a.sub.2],[a.sub.3]).sup.II] = (155,18,-19/6, -41/6)) Table 2: Input variables used for MVA to separate the signal from the background and their relative importance (RI). These numbers are shown for MS _ 1 TeV. Variables RI x [10.sup.-1] [p.sub.T]([[gamma].sub.1]) 1.22 [p.sub.T]([[gamma].sub.2]) 1.31 [eta]([[gamma].sub.2]) 1.02 [eta]([[gamma].sub.1]) 1.11 [DELTA]R([[gamma].sub.1], [[gamma].sub.2]) 1.30 M([[gamma].sub.1], [[gamma].sub.2]) 3.25 [N.sub.jet] 0.78

Printer friendly Cite/link Email Feedback | |

Author: | Aydemir, Ufuk; Mandal, Tanumoy |
---|---|

Publication: | Advances in High Energy Physics |

Date: | Jan 1, 2017 |

Words: | 10141 |

Previous Article: | Entropy Product Formula for Gravitational Instanton. |

Next Article: | Cosmic Microwave Background as a Thermal Gas of SU(2) Photons: Implications for the High-z Cosmological Model and the Value of [H.sub.0]. |