#
, Higgses and heavy neutrinos in models:

from the LHC to the GUT scale

###### Abstract

We study a class of non-exotic minimal extensions of the Standard Model, which includes all scenarios that are anomaly-free with the ordinary fermion content augmented by one Right-Handed neutrino per generation, wherein the new Abelian gauge group is spontaneously broken by the non-zero Vacuum Expectation Value of an additional Higgs singlet field, in turn providing mass to a state. By adopting the example, whose results can be recast into those pertaining to the whole aforementioned class, and allowing for both scalar and gauge mixing, we first extract the surviving parameter space in presence of up-to-date theoretical and experimental constraints. Over the corresponding parameter configurations, we then delineate the high energy behaviour of such constructs in terms of their stability and perturbativity. Finally, we highlight key production and decay channels of the new states entering the spectra of this class of models, i.e., heavy neutrinos, a second Higgs state and the , which are amenable to experimental investigation at the Large Hadron Collider. We therefore set the stage to establish a direct link between measurements obtainable at the Electro-Weak scale and the dynamics of the underlying model up to those where a Grand Unification Theory embedding a can be realised.

## 1 Introduction

The evolution of the three Standard Model (SM) gauge couplings through the Renormalisation Group (RG) equations shows a remarkable convergence, although only approximate, around GeV. This feature is even more evident in a supersymmetric context and represents one of the most solid hints in favour of a Grand Unification Theory (GUT). One of the main predictions of such theories is the appearance of an extra gauge symmetry which can be broken at energies accessible at the Large Hadron Collider (LHC). If this new Abelian group is broken around the TeV scale, the associated neutral gauge boson could be observed as a resonance in the di-lepton channel or elsewhere.
The search for extra neutral currents at the LHC is for sure one of the most
important topics of current interest and is expected to be so also in the future. In fact, stronger and stronger exclusion limits on the mass of an hypothetical
extra neutral gauge boson, which have been and are being generated by the large amount of data collected at Run 1 and Run 2 by the ATLAS and CMS collaborations, do not actually rule out the possibility of discovering such neutral currents altogether. This is due to the structure of the parameter space of the class of models embedding these states which, even in their most economical form, can easily comply with experimental constraints.

There are several realisations of GUTs that may predict bosons (see, for instance, Langacker:1980js (); Hewett:1988xc ()), such as string theory motivated , and Left-Right (LR) symmetric models. In fact, bosons may also appear in effective theories from string compactifications Faraggi:1990ita (); Faraggi:2015iaa (); Faraggi:2016xnm () or from Kaluza-Klein constructions Randall:1999ee (). Each of these is characterised by distinct mass spectra, scalar potential and flavour structures.
Even in simpler bottom-up approaches, where the issue of the unification of the underlying gauge symmetry into groups of higher rank is not addressed, several possibilities remain wide open and allow the search to continue. In more general scenarios based on GUTs, however, experimental searches must confront with a larger number of allowed gauge bosons, with the appearance of several scales of sequential breakings of the gauge symmetry and, last but not least, of new matter belonging to higher representations.

In this work we address a class of minimal models which are the most economical Abelian and renormalisable extension of the SM with only few additional free parameters Appelquist:2002mw (); Chankowski:2006jk (); Langacker:2008yv (); Erler:2009jh (); Salvioni:2009mt (). At the Electro-Weak (EW) scale, these Abelian extensions of the SM, in which the gauge symmetry is described by the group, are also characterised by a new complex scalar field, heavier than the SM Higgs, whose Vacuum Expectation Value (VEV) can lie in the TeV range and provides the mass for the . Moreover, an enlarged flavour sector is also present. Indeed, in this class of models, the cancellation of the gauge and gravitational anomalies naturally predicts right-handed neutrinos Appelquist:2002mw (). The same heavy scalar develops, thorough its spontaneous symmetry breaking, a Majorana mass at the TeV scale for the SM singlet fermions, dynamically realising a low-scale seesaw mechanism. We consider three right-handed neutrinos, one for each family, as predicted by an anomaly free and flavour universal Abelian extension. See Martinez:2013qya (); Martinez:2014ova (); Martinez:2014rea (); Martinez:2015wrp () for some non-universal examples.

The special case of a “pure” extension, where the conserved charge of the extra Abelian symmetry is the number, with and the baryon and lepton numbers, respectively, has been extensively scrutinised in the literature, see Khalil:2007dr (); Basso:2008iv (); Basso:2010hk (); Basso:2010pe (); Basso:2010yz (); Basso:2010jm (); Accomando:2010fz (); Basso:2011na (); Basso:2012sz (); Basso:2012ux (); Accomando:2013sfa (); Accomando:2015cfa (); Accomando:2015ava (); Okada:2016gsh (). Here, we surpass these previous studies in several directions.

Firstly, the model setup therein is characterised by the condition of vanishing mixing between the two and gauge groups at the EW scale, an assumption that is relaxed in our analysis. Indeed, in the framework discussed in this work, an additional factor can always be described by a linear combination of the hypercharge and of the quantum number, with coefficients parameterised by the new Abelian gauge coupling and the mixing . Therefore, choosing the as a reference gauge group extension, we can explore an entire class of minimal Abelian models through the ratio of the gauge couplings . Furthermore, the addition of the mixing parameter allows to explore the phenomenology of different scenarios which are characterised by distinctive features with respect to the “pure” case. For instance, new decay channels for the will be opened which, albeit being interesting on their own, will increase the size of the total width and, therefore, the effect of the interference term in the Drell-Yan (DY) process Accomando:2013sfa (), ultimately resulting in significantly altered experimental constraints.

Secondly, we will establish a direct connection between the yield of a variety of possible LHC measuments and the high scale structure of our construction, thereby directing experimental investigation towards key analyses enabling one to make an assessment of the model stability based on their low energy spectra. Both the enlarged boson and fermion content of our setup afford one with a variety of specific signatures. For example, the extended flavour sector allows for the possibility of a decaying into long-lived heavy neutrinos with very clear multi-leptonic signatures Basso:2008iv (); Kang:2015uoc (), which highlight the striking features of these minimal Abelian extensions providing an optimal channel, in addition to the di-lepton one, for their identification and characterisation. Furthermore, the presence of an heavy scalar, which can mix with the SM Higgs boson, offers the possibility Basso:2010yz () to search for its decay into a pair of heavy neutrinos, which is a completely novel signature with respect to many beyond the SM scenarios. Then, the non-zero scalar mixing angle also provides the decay of a heavier Higgs boson state into two light ones, which represents a unique way to probe the scalar sector and the mechanism of spontaneous EW Symmetry Breaking (EWSB).

In essence, with extended scalar, gauge and flavour sectors, it is natural to ask whether the vacuum of the this class of models is in a stable configuration both at the classical and at the quantum level and what is the impact of the new physics on the SM EW ground state Basso:2010yz (); Basso:2010jm (); Basso:2013vla (); Datta:2013mta (); Chakrabortty:2013zja (); Coriano:2014mpa (); DiChiara:2014wha (); Oda:2015gna (); Das:2015nwk (); Coriano:2015sea (); Das:2016zue (). Indeed, the extrapolation of the SM to high energy scales, through the RG equations, exhibits a scalar potential with a non-trivial structure: its minimum does not correspond to the EW vacuum which is found to be in a metastable configuration, very close to a phase transition Bezrukov:2012sa (); Buttazzo:2013uya (). This scenario, dubbed near-criticality, has been clarified by the Higgs discovery and by the measurement of its mass at the LHC Run 1. For the measured value of the Higgs mass, the scalar quartic coupling is driven, during the RG evolution, to a negative value by the Yukawa coupling of the top quark. This implies that the EW vacuum is only a local minimum and can possibly decay into the true ground state through quantum tunnelling. Fortunately, this dangerous situation is safely avoided because the lifetime of our vacuum is found to be much larger than the age of the Universe. This result, however, critically depends on the pole mass of the top quark and on its extraction from experimental measurements Masina:2012tz (); Alekhin:2012py (); Bezrukov:2014ina (). This picture strictly holds under the big desert condition, namely, the absence of significant new physics effects between the EW and Planck scales and it is therefore natural to ask which are the variations induced by new degrees of freedom, such as those embedded in our model. Clearly, an extra Abelian symmetry with an augmented scalar and flavour sectors allows one to conceive the possibility of new scenarios in which the issue of stability is resolved and the criticality due to the top mass fades away.

In this analysis we will enforce the requirements of stability of the EW vacuum and perturbativity of the couplings through the RG evolution in order to identify the allowed regions in the parameter space of these minimal Abelian extensions. For this purpose we will adopt two-loop RG equations equipped with one-loop matching conditions at the EW scale Coriano:2015sea (). These results will be combined with bounds from di-lepton analyses at the LHC at 8 TeV and with exclusion limits from Higgs searches, therefore delineating the viable parameter space from both a phenomenological perspective and its theoretical consistency. Ultimately, we will assess which of the aforementioned LHC signatures ought to be accessed and when, so as to establish a direct connection between accessible EW scale spectra and a potential underlying GUT structure.

The work is organised as follows. In section 2 we briefly introduce the model and comment on its main aspects. Then, in sections 3 and 4 we look at current experimental bounds and how they impinge on the RG analysis, respectively. In section 5 we discuss the most relevant phenomenological features of our scenario at the LHC. Finally, we draw our conclusions in section 6.

## 2 The model

The model under study is described by the SM gauge group augmented by a single Abelian gauge factor which we choose to be . We investigate a minimal renormalisable Abelian extension of the SM with the only minimal matter content necessary to ensure the consistency of the theory. In practise, in order to satisfy the cancellation of the gauge and the gravitational anomalies, we enlarge the fermion spectrum by a Right-Handed (RH) neutrino, one for each generation (we assume universality between the flavour families), which has charge and is singlet under the SM gauge group. Concerning the scalar sector, in addition to the SM-like Higgs doublet , we introduce a complex scalar field to achieve the spontaneous breaking of the extra Abelian symmetry. The new scalar field is a SM singlet with charge and its vacuum expectation value , which we choose in the TeV range, provides the mass to the gauge boson and to the RH neutrinos. The latter acquire a Majorana mass through the Yukawa interactions and thus naturally implement the type-I seesaw mechanism. More complex matter configurations, usually characterised by extra fermionic degrees of freedom or by non-universal charge assignments, can also be realised but require, consequently, a different charge assignment for the extra complex scalar.

The details of the model can be found in Basso:2008iv (); Coriano:2015sea (), here we restrict the discussion on its salient features. The most general kinetic Lagrangian for the two Abelian gauge fields, allowed by gauge invariance, admits a kinetic mixing between the corresponding field strengths

(1) |

By a suitable redefinition of the gauge fields it is possible to remove the kinetic mixing , thus recasting the kinetic Lagrangian into a canonical form. In this basis the Abelian part of the covariant derivative acting on the charged field is non-diagonal and reads as

(2) |

where and are the gauge fields of the and gauge groups, respectively, while , and , are the corresponding couplings and charges. The mixing between the two Abelian groups is described by the new coupling . For the sake of simplicity we choose to work with a canonically normalised kinetic Lagrangian and a non-diagonal covariant derivative but other parameterisations are nevertheless equivalent.

Usually, an effective coupling and charge are introduced as and some specific benchmark models can be recovered by particular choices of the two gauge couplings at a given scale (usually the EW one). For instance, the pure model is obtained enforcing () which implies the absence of mixing at the EW scale. Moreover, the Sequential SM (SSM) is reproduced by requiring (), the extension is realised by the condition while the arising from unification is described by .

Therefore, a continuous variation of the mixing coupling allows to span over the entire class of anomaly-free Abelian extensions of the SM with three RH neutrinos. We stress that there is no loss of generality in choosing the gauge group to parameterise this class of minimal models because the charges of an arbitrary symmetry can always be written as a linear combination of and , as a result of the anomaly cancellation conditions.

The scalar sector is described by the potential

(3) |

which is the most general renormalisable scalar potential of a doublet and a complex scalar . The stability of the vacuum, which ensures that the potential is bounded from below, is achieved by the following conditions

(4) |

obtained by requiring the corresponding Hessian matrix to be definite positive at large field values.

After spontaneous EWSB, the mass eigenstates of the two scalars are defined by the orthogonal transformation

(5) |

where the scalar mixing angle is given by

(6) |

and the masses of the physical scalars are

(7) |

with and identified with the 125.09 GeV Higgs boson.

Eqs. (6)–(7) can easily be inverted as

(8) |

relations which can be used to define the initial conditions on the quartic couplings through the physical masses , the VEVs and the mixing angle . Notice that the light (heavy) Higgs boson couples to SM particles proportionally to (), while the interaction with the and heavy neutrinos is provided by the complementary angle ().

When the two scalars acquire non-vanishing VEVs, the neutral component of the gauge sector becomes massive and with mass eigenstates determined by two mixing angles: the usual Weinberg angle and a new mixing angle , for which

(9) |

with values in the interval .
In contrast, the charged gauge bosons are unaffected by the presence of the extra Abelian factor and their masses remain as in the SM.

The mixing angle is completely defined in terms of the mass of the , through the VEV of the singlet scalar, and of its gauge couplings. In general it is always non-vanishing unless which corresponds to the pure model. The EW Precision Tests (EWPTs) considerably constrain the mixing angle to small values, namely, Cacciapaglia:2006pk (); Salvioni:2009mt (), in which case

(10) |

This relation can be satisfied provided either or , the latter condition allowing a generous range of values for .

Finally, the Yukawa Lagrangian is

(11) |

where is the SM contribution. Notice that a Majorana mass term for the RH neutrinos, , is dynamically generated by the VEV of the complex scalar and, therefore, the type-I seesaw mechanism is automatically implemented through spontaneous symmetry breaking. The light physical neutrinos emerge as a combination of the left-handed SM neutrinos and a highly damped sterile RH component, while the heavier ones are mostly RH. The damping term in such combinations is proportional to the ratio of the Dirac and Majorana masses.

## 3 Constraints from EWPTs and LHC searches

The parameter space is subjected to well established bounds coming from EWPTs extracted from LEP2 data. These bounds can be recast into constraints for a well-defined set of higher-dimensional operators Cacciapaglia:2006pk () which describe the effects of new physics. For the Abelian extension under study, these operators have been computed in Salvioni:2009mt () in terms of the mass and gauge couplings , neglecting, however, the impact of the heavy neutrinos and of the extended scalar sector.
To these constraints we add the one drawn from the more recent data of the first Run of LHC at 8 TeV and fb, based on a signal-to-background analysis for the di-lepton (electrons and muons) channel. Next-to-Next-to-Leading-Order (NNLO) Quantum Chromo-Dynamics (QCD) effects are taken into account through a -factor correction.
We show in fig. 1 the exclusion limit at 95% Confidence Level (CL) from both EWPTs and DY studies for three values of the , namely and 3 TeV.
For the masses of the under our investigation, the LHC studies represent a strong improvement with respect to the EW related ones, with the only comparable case being the one with TeV. Consequently, we will employ all such tight bounds in the following sections.

The sequential () and pure () models are strongly constrained while the leptophobic direction in which the coupling to leptons is minimal obviously represents the least bounded charge assignment. Moreover, we have explicitly verified that the bounds from the DY analysis are not considerably modified by the values of the heavy neutrino mass and the parameters of the extra scalar sector, such dependence entering only in the total width of the boson.

The extra scalar sector is strongly constrained by Higgs searches at LEP, Tevatron and LHC experiments. The present exclusion limits are enforced using HiggsBounds arXiv:0811.4169 (); arXiv:1102.1898 (); arXiv:1301.2345 (); arXiv:1311.0055 (); arXiv:1507.06706 () and the agreement of the model with the signal strength measurements of the discovered 125.09 GeV Higgs scalar is taken into account via HiggsSignals Bechtle:2013xfa (). The results in the plane are reported in fig. 2. The most sensitive exclusion channels are depicted with different colours depending on the mass (fig. 2(a)). The most effective exclusion search, covering almost all the mass interval from 150 GeV to 450 GeV, is of a Higgs boson decaying into a pair of and bosons Khachatryan:2015cwa () (blue region). In particular, the fully leptonic and semileptonic decay channels have been considered for while for the final states containing four charged leptons, two charged leptons and two quarks or two neutrinos have been studied. Finally, in fig. 2(b), we show a compatibility fit with the Higgs signal measurements in the plane. We have chosen a fixed reference value for the mass, namely GeV, and for the heavy neutrino mass, GeV, so that only SM-like decay channels are open for the lightest scalar . The requirement of a compatibility at results anyway into a weaker bound with respect to the exclusion limits that we have taken into account in all the following analyses.

## 4 The RG analysis

The parameter space available at the EW scale can be further constrained through the requirements of perturbativity of the couplings and stability of the vacuum through the RG evolution. In this analysis we use RG equations at two-loop order and matching relations at one-loop level in perturbation theory. The latter supply the initial conditions of the RG running at a particular scale (we choose the top pole mass as a starting point). These relations provide the initial value of the running couplings computed in the renormalisation scheme as a function of the physical on-shell parameters, namely, the pole masses and the physical couplings such as the Fermi constant. The details of the implementation of the matching conditions for various Abelian extensions of the SM can be found in Coriano:2015sea (). Before moving into the details of the analysis, we stress that the role of the two-loop corrections on the functions and of the one-loop matchings on the initial conditions is crucial in order to reach definitive conclusions on the dynamics of these extensions of the SM. Indeed, it has been shown in Coriano:2015sea () that they can lead to sizeable effects on the RG evolution and usually improve the stability of the scalar potential above and beyond the SM conditions.

We identify the regions in which the class of the Abelian extensions under study possesses a stable vacuum (described by the conditions given in Eq. (4)) and is characterised by a weakly coupled regime (couplings are required to be less than ^{1}^{1}1The parameters upon which the perturbative expansion is performed are usually in the form of rather than . Different and less conservative choices are also viable as the one that can be inferred applying Naive Dimensional Analisys (NDA). This latter case would shape a perturbative regime with couplings less then with related to the degrees of freedom inside the loops. Both options have been explored in our analysis with very rare occasions of noticeable discrepancies among the two choices.) along the RG evolution, up to some given scales which will be specified below. These regions are defined in the space of the new parameters, evaluated at the EW scale, introduced by these minimal Abelian extensions. We will focus, in particular, on the impact of the Abelian gauge couplings , , of the scalar mixing angle and of the masses of the heavy Higgs and boson.

For ease of reference, the legend of the stability and perturbativity regions, according to the maximum scale up to which the vacuum is stable and the model remains perturbative, is depicted in fig. 3. In the cyan region the new parameters of these Abelian extensions are such that the stability and/or the perturbativity is lost at a scale lower than the instabilitiy scale of the SM. A model with gauge couplings lying in this region of the parameter space worsen the high energy behaviour of the SM and clearly calls, with more urgency, for an embedding into a complex scenario, such as GUT unification, already appearing below the GeV. In the blue region the extensions behave, from the stability point of view, as the SM, whose instability scale GeV. In contrast, the green and yellow regions delineate portions of the parameter space in which this class of models is more stable than the SM, up to the GUT scale and above, thus identifying them as compelling extensions of the EW theory.

In fig. 4 we show the regions of stability and perturbativity, up to some given scale , as a function of the two Abelian couplings , and for different values of the scalar mixing angle .
These results have been obtained for TeV and GeV which corresponds to a Yukawa coupling of order . This value is too small to affect the RG evolution of the quartic scalar couplings, therefore the destabilising effect of new fermionic degrees of freedom is completely suppressed. Indeed, a is, at least, required to appreciate the impact of the heavy RH neutrinos in the running of the scalar sector Coriano:2015sea (). This roughly corresponds to for .

The constraints coming from di-lepton searches at the LHC with TeV and TeV strongly restrict the allowed parameter space in the plane completely leaving out the cyan regions and therefore only selecting the configurations in which the model is at least as stable as the SM.
The dashed lines correspond to three particular and very common extensions which can be described, in our conventions, by straight lines in the plane. These are, in anti-clockwise direction, the pure , the and the extensions, while the sequential SM lies on the axis. The black dots represent the benchmark models usually addressed in the literature in which the Abelian gauge couplings are fixed to specific values. Notice also that these reference points, although allowed by EWPTs, are excluded by LHC data if TeV.

As one can easily notice from fig. 4, the effect of the mixing angle of the two scalars is crucial for identifying the regions in which the vacuum is stable. Indeed, scalar degrees of freedom, contrary to the fermionic ones, usually drive the instability scale towards higher values improving the stability of the potential. In the case (which corresponds to ), the extra scalar sector is decoupled from the SM Higgs doublet and the RG evolution of the extension shares the same behaviour of the SM if the new Abelian gauge couplings are sufficiently small. If moves away from zero, the two scalar sectors begin to communicate and the stability effect of the complex scalar becomes quickly significant, preventing the decay of the vacuum up to the GUT scale and above.

In fig. 5 we show the same study for TeV. The regions defined by the RG analysis are unchanged with respect to the TeV case but the LHC bounds become looser. This allows to explore bigger values of the Abelian gauge couplings which can even fall in a region in which the perturbativity is spoilt (cyan), although only for a small slice of the parameter space. For these values, a bigger is ineffective to increase because the poor behaviour of the model is due to the loss of perturbativity in the Abelian sector and not to the instability of the vacuum. For heavier s, the constraints from di-lepton searches at the LHC are overtaken by EWPTs which still enclose this Abelian extension in a configuration almost as stable as the SM, provided .

A similar study is carried out in the plane in order to emphasise the impact of the extended scalar sector. The results are presented in fig. 6(a) where the hatched area has been excluded by LHC data using the HiggsBounds tool. The Abelian gauge couplings used for this particular analysis are and , which have been chosen on the contour line. We have explicitly verified that different values on the same curve do not lead to any qualitative change. Interestingly, the cyan region, in which the RG behaviour of these models worsens with respect to the SM case, is completely disallowed for GeV.
Notice also that, for GeV, both stability and perturbativity are satisfied, up to the GUT scale and above, mainly for a highly-mixed scalar sector while, for heavier , the mixing angle is bounded from above and the same regions extend horizontally. These regions will eventually shrink at bigger values of the heavy Higgs mass due to a loss of perturbativity of the quartic coupling. ^{2}^{2}2The parameter space studied in fig. 6(a) is the unique case showing sensitivity to our different definitions of the perturbativity condition. Using the NDA prescription with we found a broadening of the stability and perturbativity regions at high values of , excluded by Higgs searches, with the cyan region completely covering the white space. The portion of the parameter space allowed by the HiggsBounds analysis remains unaffected by the perturbativity condition being dominated by the stability requirement.

A similar analysis, focusing on the role of an extra real scalar field, has been presented in Falkowski:2015iwa () in which analogue results have been obtained concerning the one-loop stability of the scalar potential up to the Planck scale in the plane.
The scenario presented in Falkowski:2015iwa () is not embedded into an extended gauge sector and therefore the VEV of the extra scalar is unconstrained, while in our case its value is intimately related to the mass^{3}^{3}3Moreover, we require a complex scalar since it has to provide the longitudinal degree of freedom to the new massive gauge boson.. Nevertheless, the source of the main differences is in the perturbative order of the RG evolution. Indeed,
as emphasised in Coriano:2015sea (), the role of a NLO RG study, with two-loop functions and one-loop matching equations, is to improve the stability of the potential and therefore it is found to be necessary to draw conclusive statements on the high-energy behaviour of the vacuum. To highlight the impact of a NLO analysis we show in fig. 6(b) the region of stability and perturbativity up to the GUT scale at NLO (yellow region) in comparison to a LO only (region enclosed in the dashed curve) study in which only one-loop functions and tree-level matching conditions are employed. It is evident that, in a NLO analysis, the parameter space providing a well-behaved theory up to high energies broadens towards smaller values of the scalar mixing angle , which are, quite interestingly, in the region allowed by Higgs searches at the LHC.

In fig. 7 we show regions of stability and perturbativity as a function of the top pole mass and of (fig. 7(a) with GeV) or (fig. 7(b) with ). The bold numbers on the boundaries of the different coloured regions represent the maximum scale of stability and perturbativity GeV. The dashed line corresponds to the central value of the top mass GeV ATLAS:2014wva () which is an average from the combined analysis of ATLAS, CMS, CDF and D0, extracted through Monte Carlo (MC) modelling of production and decay of the top quark in hadronic collisions. Due to its origin, the measurement leads to the so-called MC mass which does represent neither the pole mass nor the mass. Usually, one assumes that the MC mass is sufficiently close to the pole mass with differences estimated of the order of 1 GeV Buckley:2011ms (); Juste:2013dsa (); Frixione:2014ala () and then extracts its value using matching conditions at the EW scale. The corresponding Yukawa coupling is then determined and fed to the RG equations. Unfortunately, this procedure is plagued by many sources of uncertainty and therefore it would be much better, due to its critical role Bezrukov:2012sa (); Buttazzo:2013uya (); Masina:2012tz (); Alekhin:2012py (); Bezrukov:2014ina (), if the MC event generators were defined directly in terms of the Yukawa parameter. The analysis presented in fig. 7, far from being exhaustive, has the only purpose of showing how the impact of the top mass, investigated in a window of 1 and 3 according to ATLAS:2014wva (), is affected by the parameters of the enlarged scalar potential. As one naively expects, the mixing angle weakens the destabilising effect of the top (fig. 7(a)) and eventually completely overcomes it for . The restoration of the vacuum stability, for a fixed value of the top mass, also appears as one increases the mass of the heavy Higgs (fig. 7(b)). Contrary to , the effect of is softened and, in the range 150 GeV 500 GeV with , only shifts the instability induced by the top quark to higher values of its mass.

## 5 LHC Phenomenology

In this section we explore the possible experimental signatures that characterise the class of models encompassed by our general analysis.

### 5.1 production and decay

In fig. 8 the different branching ratios of a decaying in fermions are displayed for the values of and TeV and for different (dashed regions are excluded according to fig. 1). The partial decay width of the decaying into leptons and quarks is

(12) |

where is the colour factor while are the left/right-handed couplings of the fermion to the boson. These are given by

(13) |

where we have used the short-hand notation , , and , with being the third component of the weak isospin, the electric charge in unit of and . The hypercharge is normalised as while is the charge which we have identified with the quantum number. The decay width of the into heavy neutrinos is

(14) |

It is clear how the favourite channel for the pure is in two charged leptons Basso:2008iv (). This decay mode provides nearly 40% of the width, the remainder being almost equally shared by the decays into light quarks, heavy and light neutrinos (note that we considered the branching for charged leptons, light and heavy neutrino states summed over generations). When we turn our attention to the gauge mixing, the decay mode hierarchy is drastically changed. In the limit of a sequential , which is recovered for , a preference for light quark decays reaches the highest value of 60% for the Branching Ratio (BR). The leptonic decay mode is sub-dominant in this range but starts becoming sizeable with increasing . This is to be expected given the restoration of the pure case in the limit of large (equivalent to ). The BR in heavy neutrinos shows a variability with . Moving from zero in the sequential limit it reaches the highest contribution at of the BR in the pure case. Indeed, the partial width is independent of and it is solely controlled by the Abelian coupling .

The possibility to explore different model configuations is enabled by gauge mixing, which opens new decay channels into SM bosons, which are absent in the pure , namely, , and . The corresponding partial decay widths are given by

where and are the momentum and the energy of the boson in the rest frame. Moreover, in the previous equations, we have defined , , , where and are, respectively, the charges of the SM doublet and singlet scalar which, in our case, are and .

The interplay between the mixing in the Abelian and scalar sectors is visible in the corresponding BRs as given in fig. 9. The decays into charged gauge bosons and represent the main patterns
regardless of the value of the scalar mixing angle in the range ( is a very conservative choice, larger values are possibile depending on the mass, see fig. 2) with kinematics accounting for the main differences.
The non-zero scalar mixing also clears the way for a channel but with a highly dumped BR not exceeding the 0.1% value. To understand these features, it is instructive to study the partial widths in eq. (5.1) in the regime taking into account the smallness of the gauge mixing angle through Eq. (10). In this limit we obtain

(16) |

which clearly describes the behaviour depicted in fig. 9.

We concentrate now on the on-shell production of a gauge boson through DY mode to accomodate the discovery/exclusion opportunities of our model in LHC Run 2. The computation has been performed using CalcHep Belyaev:2012qa () and the corresponding model file implementation Basso:2010jm (); Basso:2011na () on the
High Energy Physics Model Data-Base (HEPMDB) hepmdb ().
From this perspective, we present in fig. 10 the corresponding cross section at 13 TeV as a function of and for different values of
and mass. We consider the bounds coming from the previous significance analysis from DY at LHC Run 1 and highlight the excluded with dotted lines.
The of the pure model, which is strongly constrained in terms of , is characterised by a cross section up to fb for
TeV and up to fb for TeV. Increasing may increase the coupling to quarks and also allow higher values of and consequently more sizeable cross sections but
without exceeding the fb, a value approached at TeV and .

#### 5.1.1 Hallmark LHC signatures from a

The production of heavy neutrinos from decay is a smoking-gun signal of the particular minimal class of models considered, where an extended fermion sector is required to cure the anomalies of the new gauge boson. The successive decays of the heavy neutrino may result in distinctive multi-lepton signatures which have been under recent investigation (see, for instance, Emam:2007dy () for the 2-lepton, Basso:2008iv () for the 3-lepton and Huitu:2008gf (); Khalil:2015naa () for the 4-lepton channel). We explore here the role played, in this process, by the new Abelian couplings and different assignments of and masses computing the cross section for the production of heavy neutrinos from a decaying . The results are plotted in fig. 11 with contour plots computed for a Centre-of-Mass (CM) energy of 13 TeV.

### 5.2 Higgs production and decay

In this section we address the collider perspectives for a scalar signal of origin at the LHC. In our setup the parameter space defining the new scalar sector consists of the mass of the physical heavy scalar and the related scalar mixing angle . The mixing angle plays, as expected, a central role scaling all the interactions with SM-like particles by () when involving a () and with the complementary angle when involving particles in the peculiar spectrum of the model (as and heavy neutrinos). Also, the case with offers the chance to investigate new decay channels with multi-scalars, an important hallmark of the mechanism responsible for our extended spontaneous EWSB. Here, we build on the results presented in Basso:2010yz () for a pure scenario taking into account the new exclusion data from Higgs searches and the impact of the gauge mixing coupling .

#### 5.2.1 Standard production mechanisms

The most important set of mechanisms exploited to reveal the SM-like GeV Higgs boson at LHC involve gluon-gluon fusion, vector-boson fusion, associated production and Higgs-strahlung. The cross sections for these standard production channels of the light scalar , which we assume to coincide with the 125.09 GeV Higgs boson, can be simply obtained from the SM results by a rescaling with a factor. Here, instead, we present in fig. 12 the cross sections related to such processes for the case of the heavy scalar Higgs () production as function of its mass and with the benchmark value of the mixing angle for and TeV as CM energy at the LHC. The hierarchy of the cross sections is the same as for the SM Higgs case, the couplings to SM particles being rescaled by a factor of .

#### 5.2.2 Non-standard production mechanisms

The connection of the extended scalar sector with the remaining particles allows for new mechanisms for heavy Higgs production. Among these, the associated production with the boson is of great importance, opening a window towards the -specific spectrum. In fig. 13 we plot the variation of the cross section for the process with respect to the scalar mixing angle. A fixed value of the heavy scalar mass has been taken and different benchmarks of mass and couplings have been considered. Notice that, due the mixing, the same final state can be obtained with a exchanged in the -channel. (We have verified that this contribution and its interference with the diagram are non-negligible). The influence of the gauge sector in this production mechanism is translated in the enhancing effect from the Abelian gauge couplings and leads to a maximum value of fb. Despite the small cross section, this is the only accessible production channel for when . The ensuing couplings have been chosen appropriately within the 95% CL area of fig. 1 and compensate for the dumping effect in the cross section due to the increasing mass.

#### 5.2.3 BRs and widths of the Higgs bosons

We now move to the investigation of the various decay modes of in two particle final states and the role played by the related unknown parameter space. We begin by studying the variation of the branchings of for a change of its mass in the range 150–500 GeV. Two benchmark points have been considered with two assignments of the scalar mixing angle, consistent with the bounds extracted from Higgs searches, and a common value for the heavy neutrino and masses, as for the Abelian gauge couplings set at and . The resulting BRs are shown in fig. 14. With respect to the SM case, new decay channels are accessible, namely, the and , the former almost ubiquitous in many extensions of the scalar sector, the latter being a hallmark of scenarios. For both values of the mixing considered, and 0.28, the main channel is represented by the decay into charged gauge boson, a predominance which is weakly challenged only by the decay in two s and, when overcomes the threshold at GeV, by the one in two light scalars. Indeed, the hierarchy of the different decay modes in SM final states is the same as that of the SM Higgs, the partial decay widths being rescaled by a factor of . The scalar mixing enters critically in the BRs into heavy neutrinos. When the corresponding kinematical region is allowed, it is evident that a heavy Higgs mainly projecting onto the SM scalar singlet (for smaller values of ) has, in our model, a stronger interaction with the heavy neutrinos, and, at the same time, a weaker coupling to SM particles. The corresponding BR endures a one order of magnitude suppression when is raised to 0.28.

In fig. 15 we show the BRs as a function of the scalar mixing angle for two values of its mass in order to explore different kinematical regions. Indeed, moving from the GeV to the GeV case, the decays in a top quark pair and in two become accessible. As mentioned before, the role of , for the interaction structure of our model, is clarified by the interplay between the decay in heavy neutrino and the other modes. In both cases shown in fig. 15 the increase in causes the dropping of the heavy neutrino decay mode and a growth of the SM-like decay channels. Notice also that the mode does not have a trivial dependence on .

In fig. 16 the dependence on and of the heavy Higgs total width is illustrated. In fig. 16(a) the heavy scalar masses were allowed to span in the range 150 GeV 500 GeV while three different assignments have been considered. The case with zero mixing singles out in showing a recognisable threshold due to the heavy neutrino decay being the only allowed channel. The values of the width rapidly grow when such threshold is exceeded reaching the MeV order. Further, with the increase of the scalar mixing the width experiences another sizeable growth due to the now open SM decay channels. Also for such cases the channel is available resulting in a mild threshold in the width plot. We can appreciate how the non-zero mixing causes a large increment in the width allowing values of order GeV to be reached for high values. The critical role of the scalar mixing angle is more visible in fig. 16(b) where we considered the variation of the width respect to in the range . For the given choices of , the constraints coming from Higgs searches at the LHC have been taken into account excluding a large sector (dashed lines) of the values of in the plot.

Let us now turn to the decay patterns of the SM-like Higgs state, . When a new interesting channel become accessible to it, (into heavy neutrinos), otherwise it behaves as the SM Higgs boson, with the same BRs and a total width rescaled by a factor of . We show in fig. 17(a) the light Higgs decay mode into a pair of heavy neutrinos for GeV and for three different benchmark points. For comparison we also show the BRs of some decay channel of the SM Higgs boson. Quite interestingly the neutrino BR spans from 0.1% to 1% becoming comparable to, or even exceeding, the mode of the SM Higgs. The behaviour of the depicted curves can be understood by scrutinising the structure of the vertex. This is proportional to and therefore, for fixed , can be increased by growing the ratio . Taking into account the LHC limits on the Abelian gauge couplings discussed in section 3, which are obviously more constraining for lower masses, we find a bigger ratio for TeV, in which case is allowed to vary up to 0.6. For completeness, we depict in fig. 17(b) the BR values for the process at the LHC with TeV CM energy, which can reach 100 fb. Notice that the production cross section scales with a factor of with respect to the SM case, which is reproduced by a vanishing scalar mixing angle. In such case pb LHCHXSWG () which has been used to normalise our cross section.

#### 5.2.4 Hallmark LHC signatures from Higgs states

The production cross sections and decay BRs of can be combined with the recent limits, coming from LHC search on the extended Higgs sector, to probe realistic
discovery opportunities. Our phenomenological scenario calls for a TeV CM energy and an integrated luminosity of 100 fb, as expected
to be collected at LHC.
From what has been illustrated in the previous analysis, the heavy scalar decay can reveal its presence and that of the remaining beyond-SM spectrum
through peculiar decay channels. Such distinctive signatures involve heavy neutrinos and light scalars. Considering production from gluon-gluon fusion we
project in the plane, fig. 18(a), the contour of equal value for the cross section times BR of the process
.
We kept the heavy neutrinos at a common degenerate mass of 95 GeV, summing the final state over generations, and considered the benchmark point in the
extended gauge sector with 2 TeV, and . The values of ,0.2 fb and 0.5 fb illustrate the magnitude involved and
the number of neutrino events that can be expected. We crossed the results with the stability/perturbativity implications of a given choice of the parameter space.
We notice how the request to exceed 50 events selects a restricted area of the heavy scalar mass,
roughly 200 GeV 250 GeV, with values of the scalar mixing not excluded (hatched area) by LHC data. The same area covers a region with a scale of
stability/perturbativity breaking greater than the SM case.
A more generous response is obtained when the gluon-gluon cross section is multiplied for the branching of .
In fig. 18(b) are drawn, for the latter process and the same setting of masses and gauge parameters of the previous figure,
the contours with 0.1, 1, 100 and 200 fb. Above the threshold GeV,
the scalar mixing angle can critically raise the value of leading potentially to 100 events.
The LHC limits severely intervene to exclude large value of with the resulting effect of an upper bound of 200 events in the space investigated.

The decay in light scalars or heavy neutrinos states represents a peculiar feature of our minimal class of models, nevertheless
a search aimed to a heavy scalar discovery would favour different channels.
From the previous analysis of the BRs (see figs. 14-15), decays in , and are
the main candidates as search channels. Consequently, we proceed by testing the gluon-gluon induced cross section of such channels against the LHC exclusion limits in fig. 19.
The corresponding contours of equal value for the cross section of and are illustrated
in figs. 19a-b. The two cases share the absence of a threshold in the interval of considered
and a cross section increasing with the scalar mixing. At the highest values of mixing allowed the decay is more capable to get close to 1 pb while the decay has a
weaker growth as can be read off from the path of the line pb. The process completes our survey. The threshold is
sufficiently high to concern only a small section of the plane. The values of the cross section times BR depicted are for 10,25,50 fb.

## 6 Conclusions and outlook

In summary, we have shown how production and decay patterns peculiar to a class of models (of which we have taken the case as an example) involving the entirity of their additional particle spectrum, i.e., heavy neutrinos, a second Higgs state and a , at times interplaying with each other in experimental signatures accessible at the second stage of the LHC, can be linked to the high scale behaviour of such scenarios. This has been made possible by combining the description of their low and high energy dynamics through an advanced RG analysis which specifically used as boundary conditions only those potentially accessible by experiment at present and in the near future.

Our study has so far been confined to the inclusive level only, as no dedicated MC analysis has been attempted yet to uncover the interesting signatures at the LHC and extract from these the values of the fundamental parameters defining these scenarios. Indeed, building on the results obtained here, we intend to carry out this task in a forthcoming publication which, once accomplished, will translate into the ability to use LHC data in order to open a window on high-scale physics, possibly providing circumstantial evidence of its ultimate structure. In fact, our initial study has already made clear that several signatures potentially accessible at the CERN machine during Run 2 cover a sizeable region of the parameter space which points towards a well-behaved scenario, stable and perturbative up to the GUT scale, thus supporting these models as viable extensions of the SM.

L.D.R. thanks L. Basso for useful discussions during the development of this manuscript. E.A., J.F. and S.M. are supported in part through the NExT Institute. C.C. thanks The Leverhulme Trust for a Visiting Professorship to the University of Southampton, where part of his work was carried out. The work of L.D.R. has been supported by the "Angelo Della Riccia” foundation.

## References

- (1) P. Langacker, Grand Unified Theories and Proton Decay, Phys. Rept. 72 (1981) 185.
- (2) J. L. Hewett and T. G. Rizzo, Low-Energy Phenomenology of Superstring Inspired E(6) Models, Phys. Rept. 183 (1989) 193.
- (3) A. E. Faraggi and D. V. Nanopoulos, A SUPERSTRING AT O (1-TeV) ?, Mod. Phys. Lett. A6 (1991) 61–68.
- (4) A. E. Faraggi and M. Guzzi, Extra s and s in heterotic-string derived models, Eur. Phys. J. C75 (2015), no. 11 537, [arXiv:1507.0740].
- (5) A. E. Faraggi and J. Rizos, The 750 GeV di-photon LHC excess and extra s in heterotic-string derived models, Eur. Phys. J. C76 (2016), no. 3 170, [arXiv:1601.0360].
- (6) L. Randall and R. Sundrum, A Large mass hierarchy from a small extra dimension, Phys. Rev. Lett. 83 (1999) 3370–3373, [hep-ph/9905221].
- (7) T. Appelquist, B. A. Dobrescu, and A. R. Hopper, Nonexotic neutral gauge bosons, Phys.Rev. D68 (2003) 035012, [hep-ph/0212073].
- (8) P. H. Chankowski, S. Pokorski, and J. Wagner, Z-prime and the Appelquist-Carrazzone decoupling, Eur.Phys.J. C47 (2006) 187–205, [hep-ph/0601097].
- (9) P. Langacker, The Physics of Heavy Gauge Bosons, Rev. Mod. Phys. 81 (2009) 1199–1228, [arXiv:0801.1345].
- (10) J. Erler, P. Langacker, S. Munir, and E. Rojas, Improved Constraints on Z-prime Bosons from Electroweak Precision Data, JHEP 08 (2009) 017, [arXiv:0906.2435].
- (11) E. Salvioni, G. Villadoro, and F. Zwirner, Minimal Z-prime models: Present bounds and early LHC reach, JHEP 11 (2009) 068, [arXiv:0909.1320].
- (12) R. Martinez, F. Ochoa, and J. P. Rubio, Some phenomenological aspects of a new model, Phys. Rev. D89 (2014), no. 5 056008, [arXiv:1303.2734].
- (13) R. MartÌnez, J. Nisperuza, F. Ochoa, and J. P. Rubio, Scalar dark matter with CERN-LEP data and search at the LHC in an model, Phys. Rev. D90 (2014), no. 9 095004, [arXiv:1408.5153].
- (14) R. Martinez, J. Nisperuza, F. Ochoa, J. P. Rubio, and C. F. Sierra, Scalar coupling limits and diphoton Higgs decay from LHC in an U(1)? model with scalar dark matter, Phys. Rev. D92 (2015), no. 3 035016, [arXiv:1411.1641].
- (15) R. Martinez and F. Ochoa, Spin-independent interferences and spin-dependent interactions with scalar dark matter, JHEP 05 (2016) 113, [arXiv:1512.0412].
- (16) S. Khalil and A. Masiero, Radiative B-L symmetry breaking in supersymmetric models, Phys.Lett. B665 (2008) 374–377, [arXiv:0710.3525].
- (17) L. Basso, A. Belyaev, S. Moretti, and C. H. Shepherd-Themistocleous, Phenomenology of the minimal B-L extension of the Standard model: Z’ and neutrinos, Phys. Rev. D80 (2009) 055030, [arXiv:0812.4313].
- (18) L. Basso, S. Moretti, and G. M. Pruna, Constraining the coupling in the minimal Model, J.Phys. G39 (2012) 025004, [arXiv:1009.4164].
- (19) L. Basso, A. Belyaev, S. Moretti, G. M. Pruna, and C. H. Shepherd-Themistocleous, discovery potential at the LHC in the minimal extension of the Standard Model, Eur. Phys. J. C71 (2011) 1613, [arXiv:1002.3586].
- (20) L. Basso, S. Moretti, and G. M. Pruna, Phenomenology of the minimal extension of the Standard Model: the Higgs sector, Phys.Rev. D83 (2011) 055014, [arXiv:1011.2612].
- (21) L. Basso, S. Moretti, and G. M. Pruna, A Renormalisation Group Equation Study of the Scalar Sector of the Minimal B-L Extension of the Standard Model, Phys.Rev. D82 (2010) 055018, [arXiv:1004.3039].
- (22) E. Accomando, A. Belyaev, L. Fedeli, S. F. King, and C. Shepherd-Themistocleous, Z’ physics with early LHC data, Phys. Rev. D83 (2011) 075012, [arXiv:1010.6058].
- (23) L. Basso, S. Moretti, and G. M. Pruna, Theoretical constraints on the couplings of non-exotic minimal bosons, JHEP 08 (2011) 122, [arXiv:1106.4762].
- (24) L. Basso, K. Mimasu, and S. Moretti, Z’ signals in polarised top-antitop final states, JHEP 09 (2012) 024, [arXiv:1203.2542].
- (25) L. Basso, K. Mimasu, and S. Moretti, Non-exotic signals in , and final states at the LHC, JHEP 11 (2012) 060, [arXiv:1208.0019].
- (26) E. Accomando, D. Becciolini, A. Belyaev, S. Moretti, and C. Shepherd-Themistocleous, Z’ at the LHC: Interference and Finite Width Effects in Drell-Yan, JHEP 10 (2013) 153, [arXiv:1304.6700].
- (27) E. Accomando, A. Belyaev, J. Fiaschi, K. Mimasu, S. Moretti, and C. Shepherd-Themistocleous, Forward-backward asymmetry as a discovery tool for Z? bosons at the LHC, JHEP 01 (2016) 127, [arXiv:1503.0267].
- (28) E. Accomando, A. Belyaev, J. Fiaschi, K. Mimasu, S. Moretti, and C. Shepherd-Themistocleous, as a discovery tool for bosons at the LHC, arXiv:1504.0316.
- (29) N. Okada and S. Okada, portal dark matter and LHC Run-2 results, Phys. Rev. D93 (2016), no. 7 075003, [arXiv:1601.0752].
- (30) Z. Kang, P. Ko, and J. Li, New Avenues to Heavy Right-handed Neutrinos with Pair Production at Hadronic Colliders, Phys. Rev. D93 (2016), no. 7 075037, [arXiv:1512.0837].
- (31) L. Basso, Minimal Z’ models and the 125 GeV Higgs boson, Phys. Lett. B725 (2013) 322–326, [arXiv:1303.1084].
- (32) A. Datta, A. Elsayed, S. Khalil, and A. Moursy, Higgs vacuum stability in the extended standard model, Phys.Rev. D88 (2013), no. 5 053011, [arXiv:1308.0816].
- (33) J. Chakrabortty, P. Konar, and T. Mondal, Constraining a class of B?L extended models from vacuum stability and perturbativity, Phys. Rev. D89 (2014), no. 5 056014, [arXiv:1308.1291].
- (34) C. Coriano, L. Delle Rose, and C. Marzo, Vacuum Stability in U(1)-Prime Extensions of the Standard Model with TeV Scale Right Handed Neutrinos, Phys. Lett. B738 (2014) 13–19, [arXiv:1407.8539].
- (35) S. Di Chiara, V. Keus, and O. Lebedev, Stabilizing the Higgs potential with a Z, Phys. Lett. B744 (2015) 59–66, [arXiv:1412.7036].
- (36) S. Oda, N. Okada, and D.-s. Takahashi, Classically conformal U(1)? extended standard model and Higgs vacuum stability, Phys. Rev. D92 (2015), no. 1 015026, [arXiv:1504.0629].
- (37) A. Das, N. Okada, and N. Papapietro, Electroweak vacuum stability in classically conformal B-L extension of the Standard Model, arXiv:1509.0146.
- (38) C. Coriano, L. Delle Rose, and C. Marzo, Constraints on abelian extensions of the Standard Model from two-loop vacuum stability and , JHEP 02 (2016) 135, [arXiv:1510.0237].
- (39) A. Das, S. Oda, N. Okada, and D.-s. Takahashi, Classically conformal U(1)’ extended standard model, electroweak vacuum stability, and LHC Run-2 bounds, arXiv:1605.0115.
- (40) F. Bezrukov, M. Y. Kalmykov, B. A. Kniehl, and M. Shaposhnikov, Higgs Boson Mass and New Physics, JHEP 1210 (2012) 140, [arXiv:1205.2893].
- (41) D. Buttazzo, G. Degrassi, P. P. Giardino, G. F. Giudice, F. Sala, et. al., Investigating the near-criticality of the Higgs boson, JHEP 1312 (2013) 089, [arXiv:1307.3536].
- (42) I. Masina, Higgs boson and top quark masses as tests of electroweak vacuum stability, Phys.Rev. D87 (2013), no. 5 053001, [arXiv:1209.0393].
- (43) S. Alekhin, A. Djouadi, and S. Moch, The top quark and Higgs boson masses and the stability of the electroweak vacuum, Phys. Lett. B716 (2012) 214–219, [arXiv:1207.0980].
- (44) F. Bezrukov and M. Shaposhnikov, Why should we care about the top quark Yukawa coupling?, J. Exp. Theor. Phys. 120 (2015) 335–343, [arXiv:1411.1923]. [Zh. Eksp. Teor. Fiz.147,389(2015)].
- (45) G. Cacciapaglia, C. Csaki, G. Marandella, and A. Strumia, The Minimal Set of Electroweak Precision Parameters, Phys.Rev. D74 (2006) 033011, [hep-ph/0604111].
- (46) P. Bechtle, O. Brein, S. Heinemeyer, G. Weiglein, and K. E. Williams, HiggsBounds: Confronting Arbitrary Higgs Sectors with Exclusion Bounds from LEP and the Tevatron, Comput. Phys. Commun. 181 (2010) 138–167, [arXiv:0811.4169].
- (47) P. Bechtle, O. Brein, S. Heinemeyer, G. Weiglein, and K. E. Williams, HiggsBounds 2.0.0: Confronting Neutral and Charged Higgs Sector Predictions with Exclusion Bounds from LEP and the Tevatron, Comput. Phys. Commun. 182 (2011) 2605–2631, [arXiv:1102.1898].
- (48) P. Bechtle et. al., Recent Developments in HiggsBounds and a Preview of HiggsSignals, PoS CHARGED2012 (2012) 024, [arXiv:1301.2345].
- (49) P. Bechtle et. al., HiggsBounds-4: Improved Tests of Extended Higgs Sectors against Exclusion Bounds from LEP, the Tevatron and the LHC, Eur. Phys. J. C74 (2014) 2693, [arXiv:1311.0055].
- (50) P. Bechtle, S. Heinemeyer, O. Stal, T. Stefaniak, and G. Weiglein, Applying Exclusion Likelihoods from LHC Searches to Extended Higgs Sectors, arXiv:1507.0670.
- (51) P. Bechtle, S. Heinemeyer, O. StÂl, T. Stefaniak, and G. Weiglein, : Confronting arbitrary Higgs sectors with measurements at the Tevatron and the LHC, Eur. Phys. J. C74 (2014), no. 2 2711, [arXiv:1305.1933].
- (52) CMS Collaboration, V. Khachatryan et. al., Search for a Higgs Boson in the Mass Range from 145 to 1000 GeV Decaying to a Pair of W or Z Bosons, JHEP 10 (2015) 144, [arXiv:1504.0093].
- (53) CMS Collaboration, Properties of the Higgs-like boson in the decay H to ZZ to 4l in pp collisions at sqrt s =7 and 8 TeV, CMS-PAS-HIG-13-002.
- (54) CMS Collaboration, Update on the search for the standard model Higgs boson in pp collisions at the LHC decaying to W + W in the fully leptonic final state, CMS-PAS-HIG-13-003.
- (55) CMS Collaboration, Combination of standard model Higgs boson searches and measurements of the properties of the new boson with a mass near 125 GeV, CMS-PAS-HIG-12-045.
- (56) A. Falkowski, C. Gross, and O. Lebedev, A second Higgs from the Higgs portal, JHEP 05 (2015) 057, [arXiv:1502.0136].
- (57) ATLAS, CDF, CMS, D0 Collaboration, First combination of Tevatron and LHC measurements of the top-quark mass, arXiv:1403.4427.
- (58) A. Buckley et. al., General-purpose event generators for LHC physics, Phys. Rept. 504 (2011) 145–233, [arXiv:1101.2599].
- (59) A. Juste, S. Mantry, A. Mitov, A. Penin, P. Skands, E. Varnes, M. Vos, and S. Wimpenny, Determination of the top quark mass circa 2013: methods, subtleties, perspectives, Eur. Phys. J. C74 (2014) 3119, [arXiv:1310.0799].
- (60) S. Frixione and A. Mitov, Determination of the top quark mass from leptonic observables, JHEP 09 (2014) 012, [arXiv:1407.2763].
- (61) A. Belyaev, N. D. Christensen, and A. Pukhov, CalcHEP 3.4 for collider physics within and beyond the Standard Model, Comput. Phys. Commun. 184 (2013) 1729–1769, [arXiv:1207.6082].
- (62) M. Bondarenko, A. Belyaev, L. Basso, E. Boos, V. Bunichev, et. al., High Energy Physics Model Database : Towards decoding of the underlying theory (within Les Houches 2011: Physics at TeV Colliders New Physics Working Group Report), arXiv:1203.1488.
- (63) W. Emam and S. Khalil, Higgs and Z-prime phenomenology in B-L extension of the standard model at LHC, Eur. Phys. J. C52 (2007) 625–633, [arXiv:0704.1395].
- (64) K. Huitu, S. Khalil, H. Okada, and S. K. Rai, Signatures for right-handed neutrinos at the Large Hadron Collider, Phys. Rev. Lett. 101 (2008) 181802, [arXiv:0803.2799].
- (65) S. Khalil and S. Moretti, The Supersymmetric Standard Model with Inverse Seesaw at the Large Hadron Collider, arXiv:1503.0816.
- (66) LHC Higgs Cross Section Working Group, https://twiki.cern.ch/twiki/bin/view/LHCPhysics/CERNYellowReportPageAt13TeV, .