#### Abstract

We investigate the possibility of TeV-scale scalars as low energy remnants arising in the nonsupersymmetric grand unification framework where the field content is minimal. We consider a scenario where the gauge symmetry is broken into the gauge symmetry of the Standard Model (SM) through multiple stages of symmetry breaking, and a colored and hypercharged scalar picks a TeV-scale mass in the process. The last stage of the symmetry breaking occurs at the TeV-scale where the left-right symmetry, that is, , is broken into that of the SM by a singlet scalar field of mass TeV, which is a component of an -triplet scalar field, acquiring a TeV-scale vacuum expectation value. For the LHC phenomenology, we consider a scenario where is produced via gluon-gluon fusion through loop interactions with and also decays to a pair of SM gauge bosons through in the loop. We find that the parameter space is heavily constrained from the latest LHC data. We use a multivariate analysis to estimate the LHC discovery reach of into the diphoton channel.

#### 1. Introduction

After 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 GUT is particularly interesting [9–24] (see [25–31] for analyses of the supersymmetric GUT). Breaking the gauge symmetry into that of the SM can be realized in a single step as well as in multiple steps by various symmetry breaking sequences. The relevant option 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 () [32–39], which is different than the left-right symmetric version since in this case and gauge couplings are different, that is, . Adopting the minimalistic approach and, therefore, keeping the initial field content (the 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 framework [23].

In this paper, we investigate the phenomenology of TeV-scale scalars as low energy remnants of the nonsupersymmetric 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 , 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 of mass around 1 TeV which is assumed to be the excitation of the neutral component of an triplet , denoted as . The field 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 is responsible for the production and decay of through loop interactions.

In our model, we assume two intermediate energy scales between the electroweak scale and the unification scale . At the scale , the is broken into the Pati-Salam group, (). The Pati-Salam group is broken into the group of the left-right model at the first intermediate energy scale , which is followed by the breaking of the left-right model into the SM at the energy scale . In our scenario, is assumed to be in the TeV-scale, while the values of and come out as predictions of the model. Note that the -parity invariance [10, 11, 41], which is a symmetry that maintains the complete equivalence of the left and the right sectors, is broken together with the in the first stage of the symmetry breaking. Therefore, the gauge couplings associated with the and gauge groups, and , evolve under the influence of different particle contents; hence , below the scale . Remember that the -parity is slightly different from the usual Lorentz parity in that the latter does not transform scalars, while the -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 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 and it can successfully be embedded in the minimal nonsupersymmetric 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 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 grand unification framework. We discuss how the two scalars, 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 and 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 , which is assumed to be broken into the SM at the TeV-scale. The breaking is realized by the neutral component ( which we denote as ) of the triplet , which is commonly preferred in the literature. Here, instead of the triplets, the doublet , which originates from the multiplet , 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 which we assume to be produced and decayed through the loop interaction with a color-triplet hypercharged scalar denoted as . originates from the decomposition of component of the multiplet into the SM group as follows:For our purpose, we take the mass of around 1 TeV, while the other components have heavier masses, ~2–5 TeV, and hence their contribution to the production and the decay of 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 acquiring a VEV. The fermion content of the model is the same as the SM. There are seven gauge bosons in the model, , (with ), and , with the gauge couplings , , and , associated with the , , and gauge symmetries, respectively. Using the notation of [21, 23], the symmetry breaking pattern of our model is given bywhere we assume TeV in our analysis.

In choosing the 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 is broken into the Pati-Salam group , is realized by the singlet of . Note that is odd under the -parity [10, 11], and hence it is broken at this stage as well. Therefore, below the scale , we have , 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 , can be accomplished by acquiring a VEV. The breaking of down to the SM gauge group is achieved by the multiplet which belongs to the Pati-Salam multiplet which is a member of the multiplet . In our model, acquires a VEV at around TeV which also set the value of the symmetry breaking scale . Note that is the regular triplet usually used in the literature in order to break the 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 and the EWSB scale , which are and , where the value of is chosen to be 5 TeV.

The TeV-scale left-right model with light colored scalars in the minimal nonsupersymmetric 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 (shown in (1)) into the SM gauge group, which is whose mass is ~1 TeV. Therefore, the renormalization group (RG) running of the gauge couplings at this energy scale is slightly different. The other particle 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 up to with Roman numerals as follows:The boundary/matching conditions we impose on the couplings at the symmetry breaking scales areThe low energy data which we will use as boundary conditions to the RG running are [47, 48]and all are evaluated at GeV, which givesNote that the coupling constants are all required to remain in the perturbative regime during the evolution from down to .

##### 3.2. One-Loop RG Running

For a given particle content; the gauge couplings, in an energy interval , are evolved according to the one-loop RG relationwhere the RG coefficients are given by [49, 50] asHere, the two summations are over irreducible chiral representations of fermions and those of scalars . The coefficient 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 is and is the Dynkin index of each representation. For group, andwhere is the charge, the factor of coming from the traditional normalizations of the hypercharge and charges. The ’s differ depending on the particle content in each energy interval, which changes every time symmetry breaking occurs. We will distinguish the ’s in different intervals with the corresponding roman numeral superscript, cf. (3).

##### 3.3. Results

The scalar content in the energy intervals are It is common in the literature that another scalar Pati-Salam multiplet, , 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 (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])where . Using these equations and the experimentally measured quantities in (8) and demanding TeV, we obtain the following values:The value for the scale 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 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:as . The running of the couplings is given in Figure 1. Similarly, the gauge couplings at TeV are obtained aswhich, together with the values of the symmetry breaking scales in (16), are the main predictions of the model. Notice that the value of is different from the value of , which is expected due to the fact that the -parity invariance is broken together with the symmetry; hence below the unification scale , as mentioned previously. The model also predicts the existence of TeV-scale gauge bosons and whose masses at are given aswhere we choose TeV which, together with (18), yieldsThese are the specific predictions of our model. However, we note that and change significantly with the choice of the symmetry breaking scale . 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 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 to be Majorana in character. Although there is no mechanism that constrains right-handed neutrino mass 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 (which is the case in our model as can be seen in (18)), the 14 TeV LHC searches can probe the range TeV for = 100–700 GeV.

#### 4. Phenomenology

In Section 2, we have discussed that the SM-singlet can be as light as ~1 TeV and can potentially be observed at the LHC. Since is a SM-singlet, it can not directly couple to the SM fermions and gauge bosons through any dimension-4 operator due to gauge invariance. Therefore, in order to produce 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 hypercharged scalar that appears naturally in our model and can serve both purposes, production and decay of through loop interactions. Note that the EM charge of 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 to diphoton. We further assume that 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 from gluon-gluon fusion and its decay to two photons through in the loop.

##### 4.1. Production and Decay

The scalar being singlet in nature, there is no tree level couplings of 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 with the SM Higgs doublet, which can lead to a mixing between and the SM Higgs, after EWSB. Consequently, 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 is expected to be small. For simplicity, we consider the - mixing, and therefore the partial widths of to two SM fermions or two Higgs bosons are negligible. Since carries color and hypercharge, it couples to the gluon and the (hypercharge) fields. Note that there is no coupling between and bosons, since is a singlet under . In the effective Lagrangian, we have the following dimension-5 operators for the interactions of with the SM gauge bosons prior to EWSB. where and are the field-strength tensors for and gauge groups, respectively. Effective couplings and are associated with the gluon and the fields, respectively. These couplings can be computed from the knowledge of the trilinear coupling related to the interaction term. In general, for number of colored scalars with hypercharge and for an interaction term , the effective couplings are expressed aswhere is the dimension of the representation (e.g.*, * for triplet and for octet representations) and is the index of the representation (e.g.*, * for triplet and for octet representations). The strong and the electromagnetic couplings are denoted by and , respectively. The cosine of the Weinberg angle is denoted as . The loop function is given byFor only one colored triplet and hypercharged () scalar , , , and . To keep our results as model independent as possible, we assume , where is some new physics scale (this can be chosen as ) for which we choose 5 TeV for all our computations and we keep as a free parameter. The BRs of to , , , and modes are 90.6%, 5.6%, 3.3%, and 0.5% respectively.

Here, we assume that , and therefore cannot decay to a pair. It is important to note that the BR depends only on , not on the other parameters. This is because all the partial widths, and hence the total width, scale as and the loop function (for any values of and ) would be the same for all the partial widths. The scalar has the largest BR in the dijet channel. We expect the BR in the , , and are of similar order, but mode is suppressed due to its phase space factor. The total width is a function of , , and . In Figure 3, we show as functions of for TeV for three different values of assuming TeV. As mentioned previously, scales as and one can easily estimate the total width for other values of 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 FEYNRULES2.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 for all our numerical computations.

For our phenomenological analysis, we have only three free parameters, namely, , , and (we choose TeV for all our numerical computations). We first derive bounds on the parameters from the latest LHC 13 TeV [62, 63], [64, 65], [66], and [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 byThese values are used in Figure 4(a) where we show the excluded parameter space (colored regions) in plane for TeV. The excluded regions shown in orange, green, blue, and brown are derived from the , , , and resonance search data. We can see that the diphoton data is the most powerful in constraining the parameter space in plane. In Figure 4(b), we present the excluded regions in plane for different with TeV from the latest 13 TeV combined ATLAS and CMS diphoton resonance search data. Cross section ULs () from different experiments and the corresponding uncertainties () are combined statistically using the following relations:where is the combined cross section and is the uncertainty associated with it. In case of asymmetric uncertainties, we get by averaging upper and lower uncertainties. Although uncertainties are used to compute , 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 . If , the decay becomes kinematically allowed and becomes the dominant decay mode of . This will make the diphoton and other branching modes suppressed. Therefore, observing in the , , , and resonance searches becomes much more challenging. One should note that exclusion regions are not very sensitive to the values for a fixed . 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 should remain insensitive for heavier resonances. The reduction in the production cross section as we increase is compensated by the slight change in since quantity is very sensitive to the . 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 -factor of 2 to account for the higher-order effects [69].

**(a)**

**(b)**

##### 4.3. Future Prospects at the LHC

In this subsection, we look at the prospect of discovering 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- algorithm [73] with the clustering parameter, . We use TMVA [74] for the multivariate analysis.

Signal events are generated with up to two jets, that is, jets which are MLM [75] merged at a matching scale GeV. The dominant (roughly 90%) SM background for this signal comes from the process. Similar to the signal, we generate this background by merging jets processes at 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 and 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 once it is properly chosen.

The 13 TeV diphoton data already set an UL on 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 slightly depend on the width of the resonance but we use a fixed width of 1 GeV for all in the following analysis for simplicity. The width of is a function of model parameters, namely, , , and . Instead of choosing a specific benchmark, we use GeV for our MVA. This analysis is insensitive to the actual width choice as long as ; that is, the narrow width approximation is well-valid.

We generate signal and background events with some basic transverse momentum (), pseudorapidity () and separation in plane () cuts as follows:where . We use a strong selection cut on the invariant mass of the photon pair, GeV to 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: , , , , , , and jet multiplicity ( and are the two selected photons ordered according to their ). 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 TeV. We find that the two variables and are very effective in discriminating signal from background. Other variables like and 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 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 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 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 distributions for the signal and the background for TeV at the 13 TeV LHC with . This plot is shown for the significance of where and the numbers of signal and background events that survive after the optimal BDT cut (>0) are with cut efficiency 0.75 and with cut efficiency 0.1, respectively.

**(a)**

**(b)**

In Figures 7(a) and 7(b), we show the discovery contours in plane for different at the 13 TeV LHC for 100 and 300 fb^{−1} integrated luminosities, respectively. As stated earlier, the sky-blue region, that is, , is not considered in our analysis. We observe that the discovery reach for 100 fb^{−1} run in Figure 7(a) is not much improved from the bounds obtained in Figure 4. But for 300 fb^{−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 plane for fb^{−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 significance. As previously mentioned, the limits on 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.

**(a)**

**(b)**

**(c)**

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 TeV. We apply further the following hard cuts on photons, namely, GeV and 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 ~3*σ* 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*σ*. 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 ) 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 grand unification framework. In particular, we investigate the LHC phenomenology of a SM-singlet scalar which interacts with gluons and photons through loop interactions with a color-triplet hypercharged scalar which is remnant from the breaking of the Pati-Salam gauge group, . The part of the model that lies in the TeV-scale is the left-right model, whose gauge group is , augmented with the color-triplet scalar . The scalar is a component of an 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 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 with gluons and photons that lead to a diphoton final state after being produced via gluon fusion. In addition to the decay, can also decay to , , and modes. We present the exclusion region in plane for a benchmark resonance mass 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 at the LHC. As a prospect study, we compute the higher-luminosity LHC discovery reach of by using a state-of-the-art multivariate technique. We present 5*σ* discovery contours for different choices in the plane at the 13 TeV LHC with 100 and 300 fb^{−1} integrated luminosity. From our analysis, we find that for , TeV and TeV can easily be observed with confidence level at the 13 TeV LHC with 300 fb^{−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 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 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.

#### 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.