Heavy Z ' bosons in the secluded U ( 1 ) ' model at hadron colliders.

We study Z ' phenomenology at hadron colliders in an U ( 1 ) ' extended MSSM. We choose a U ( 1 ) ' model with a secluded sector, where the tension between the electroweak scale and developing a large enough mass for Z ' is resolved by incorporating three additional singlet superfields into the model. We perform a detailed analysis of the production, followed by decays, including into supersymmetric particles, of a Z ' boson with mass between 4 and 5.2 TeV, with particular emphasis on its possible discovery. We select three different scenarios consistent with the latest available experimental data and relic density constraints, and concentrate on final signals with 2 ℓ + ⧸ E T , 4 ℓ + ⧸ E T and 6 ℓ + ⧸ E T . Including the SM background from processes with two, three or four vector bosons, we show the likelihood of observing a Z ' boson is not promising for the HL-LHC at 14 TeV. While at 27 and 100 TeV, the situation is more optimistic, and we devise specific benchmark scenarios which could be observed.

Introduction

After the discovery of the Higgs boson, the missing piece of the Standard Model (SM), the quest for physics beyond the SM (BSM) has intensified, both from theorists and experimentalists. The searches and analyses are motivated by the fact that while precise theoretical calculations within the SM have been confirmed by a wide range of experimental data establishing SM as a well-tested physics theory, it still lacks explanation for some of the fundamental phenomena, such as matter-antimatter asymmetry of the universe, dark matter, or neutrino oscillations. It is also plagued by theoretical inconsistencies, so it is at best incomplete (for instance, it does not include gravity). In fact, the discovery of the Higgs boson, with mass of the order of the electroweak scale, as expected, points towards a higher structure, because in the SM, it is unclear why the Higgs boson is so much lighter than the Planck mass, as one expects that the large quantum contributions to the square of the Higgs boson mass would inevitably make the mass huge, and comparable to the scale at which new physics appears, unless there is an incredible amount of fine-tuning cancellation between the quadratic radiative corrections and the bare mass. One can explore BSM physics taking a model independent approach, by assuming an effective field theoretic approach [1], which provides a general framework where higher order interactions of independent operators are built and one would be able to match them to explicit ultraviolet complete models in a systematic way. Or, one can enlarge the particle and/or gauge symmetry of the model.

Of the latter, the addition of supersymmetry (SUSY) to the SM is the most popular BSM scenario. It resolves the Higgs mass/gauge hierarchy problem, and provides, in its simplest scenario, the Minimal Supersymmetric Standard Model (MSSM), a natural dark matter (DM) candidate. It does not, however, explain neutrino masses and it provides no resolution for the problem [2-5]. The parameter, so-called higgsino mass term, enters the supersymmetric Lagrangian as , to give masses to the fermionic components of the Higgs bosons, and thus it is expected to be of the order of the SUSY-breaking scale. But the term also enters in the scalar potential, so, for successful electroweak symmetry breaking, its value should be at electroweak scale. Adding an gauge group to the SM/MSSM symmetry group resolves this inconsistency. An additional singlet Higgs field S develops a vacuum expectation value (VEV) which breaks the symmetry and generates the term dynamically, with . In addition, the model contains three singlet right-handed neutrinos that yield masses for the left-handed neutrinos (Dirac or Majorana).

Models with additional groups extend the spectrum of MSSM minimally: in addition to the right-handed neutrino and Higgs field S, they include another neutral gauge boson (as well as theirs supersymmetric partners). This gauge field, a consequence of the additional group, is predicted by many extensions of the SM. String-inspired models [2, 6] and grand-unification (GUT) models usually contain a number of extra U(1) symmetries. The GUT group SO(10) [7] and exceptional group [8-10] are some examples. Here the symmetries are broken at some intermediate energy scales [11], between the GUT and electroweak scales. Phenomenologically, the most interesting option is the breaking around TeV scales, giving rise to extra neutral gauge bosons observable at colliders.

The physics of bosons has been extensively studied in the literature, in models without supersymmetry [12-41], or with [5, 42–53]. The additional neutral gauge bosons have received significant attention from the experimentalists, and have been searched for extensively at the LHC, in dilepton channels [54], dijet channels [55, 56], di-tau [57], or decays [58]. Mass limits of 4 TeV or above, depending of the particular scenario chosen, have hindered extensive analyses of their implications for phenomenology, as the prospects of observing them at the LHC do not appear to be promising.

In supersymmetric models, an additional problem arises. The mass of the boson is, as usual, proportional to the VEV of the singlet Higgs boson S. But this parameter also determines the scale of the chargino/neutralino sector, thus a heavy implies a heavy electroweakino sector, reducing further the interest in such models at the LHC. To avoid this link, we work in a secluded scenario [59, 60], where the scalar sector of the model is augmented by three additional singlet superfields1, whose role is to decouple the mass of the from the scale of chargino and neutralino masses. The disadvantage is extending the particle spectrum of the SM more than the conventional model, while the gain is being able to preserve a large mass while allowing light charginos/neutralinos, and in particular, a light dark matter candidate, which is the lightest supersymmetric particle (LSP).2 In the secluded model, decays into (light) chargino and neutralino pairs, or into sfermions can be significant, and affect the mass limits, albeit slightly.

There have been studies of the boson where the mass constraints were considerably reduced, assuming the model to be leptophobic [61-63] or quark-phobic [64]. For these, couplings to the leptons or quarks are tuned (by assuming family non-universality [65], by using a specific value of the kinetic mixing [63], or by a choice of charges [65]). We propose here take a different point of view. We revisit the model with a secluded sector, and allow for heavy bosons, satisfying mass limit restrictions from LHC. We concentrate on its production, decays and observability at high-luminosity (HL), and/or high-energy (HE) frontier of a future LHC, as well as the hadronic mode of a Future Circular Collider (FCC-hh).

A heavy gauge boson, with mass TeV cannot always be treated as a narrow width, as the ratio could be . In this case, interference effects are important [66-68] and must be included. New physics contributions to the width may significantly decrease the branching ratios into SM particles, and therefore the mass limits quoted by the experiments may have to be revisited. Furthermore, decays into supersymmetric particles represent an excellent tool to investigate the electroweak interactions at the LHC in a phase-space corner that cannot be explored by employing the usual techniques. The secluded model is ideally suited for this analysis, as it allows the electroweakinos to be light. Therefore, the possible discovery of supersymmetry in mediated processes would help to understand the role of in the SUSY breaking and open the path to additional investigations, since one would need to formulate a consistent scenario accommodating both sparticles and heavy gauge bosons.

The scope of this paper is indeed the investigation of the phenomenology of bosons at the LHC, assuming that they are heavy, and that they can decay into both SM and supersymmetric particles. We will analyze the decay channels of the boson, including decays to neutralinos, charginos, sleptons and Higgs bosons, which are normally neglected.

In our study, we will allow the parameters to run within suitable ranges, taking into account the recent experimental limits. Throughout this work, we will focus especially on the decay of the into slepton, chargino and neutralino pairs, eventually leading to multilepton final states or or . To test the observability of such signals at the LHC, we devise benchmark scenarios and test their features at high integrated luminosity and at the center-of-mass and 100 TeV. We include calculations of he SM background from VV, VVV, and VVVV processes and we present a simulation analysis for the HL-LHC and future hadron colliders, indicating the significance of each scenario, and most promising observable for each signal.

Our paper is organized as follows. In Sect. 2 we describe the secluded model with particular emphasis on its neutral gauge (in Sect. 2.1) and neutralino (in Sect. 2.2) sectors. We then proceed to analyze the implications of the model at colliders in Sect. 3, focusing first in choosing three benchmarks, which obey experimental constraints, and which are able to reproduce the correct relic density, while maximizing decays into supersymmetric particles, in Sect. 3.1. Then we proceed with the analysis of production and decays. We concentrate our analysis on multilepton signals Sect. 3.2, looking at (in Sect. 3.2.1), (in Sect. 3.2.2), and (in Sect. 3.2.3) signals. We summarize our findings and conclude in Sect. 4.

The secluded model

We present here the main ingredients of the secluded model, with particular emphasis on the boson. The model is based on the gauge group , which breaks to the SM/MSSM . The additional Abelian group introduces, in addition to the MSSM superfields, three right-handed neutrino superfields (one for each generation), four singlet superfields , , and and an additional neutral gauge boson and gaugino, corresponding to the gauge sector of . While only one scalar field S is needed to break the symmetry, three additional singlets , and (the secluded sector) are introduced to split the mass scale of the additional gauge boson from that of electroweakinos.

Unfortunately, anomaly cancelation requires the presence of additional superfields (namely, the exotics and ), with exotic quantum numbers, which are assumed to be heavy and decoupled form the rest of the spectrum. We list the superfields in the model, together with the number of generations and charge assignments under the gauge group in Table 1. The superpotential in this model including the exotic fields is given bywhere the fields are the exotics, is a large mass scale and is the Yukawa coupling responsible for generating neutrino masses. This non-renormailzable term is added in the original formulation of the model to account for Dirac neutrino masses (see, for example, [69] for the origin of the term). Thus, in this form, neutrinos are Dirac particles, whose masses imply, for the Yukawa coupling [70], .

Table 1

Superfield (SF) configuration in the secluded model, including notations for fermionic and bosonic states, number of generations, and charges under

SF	Spin 0	Spin 12	Generations	Charges
Q^	Q~	Q	3	(16,2,3,QQ′)
L^	L~	L	3	(-12,2,1,QL′)
H^_d	H_d	H~_d	1	(-12,2,1,QH_d′)
H^_u	H_u	H~_u	1	(12,2,1,QH_u′)
D^	D~R∗	DR∗	3	(13,1,3¯,Qd′
U^	U~R∗	UR∗	3	(-23,1,3¯,Qu′)
E^	E~R∗	ER∗	3	(1,1,1,Qe′)
N^	N~^∗	N^∗	3	(0,1,1,Qv′)
S^	S	S~	1	(0,1,1,Qs′)
S^_1	S_1	S~_1	1	(0,1,1,Qs_1′)
S^_2	S_2	S~_2	1	(0,1,1,Qs_2′)
S^_3	S_3	S~_3	1	(0,1,1,Qs_3′)
Q^	Q~	Q	3	(Y_Q,1,1,QQ′)
Q¯^	Q¯~	Q¯	3	(Y_Q¯,1,1,Q_Q¯^′)
L^	L~	L	2	(Y_L,1,1,QL′)
L¯^	L¯~	L¯	2	(Y_L¯,1,1,Q_L¯^′)

The effective term is generated dynamically as . The scalar potential includes the F-term, given bywhile the D-term scalar potential iswhere , and are the coupling constants for the , and gauge groups while is the charge of the field . Finally, the potential includes the SUSY-breaking soft terms, expressed in terms of soft-SUSY breaking mass parameters and triple scalar couplings asThe symmetry-breaking sector of the model is very complex, and finding an acceptable minimum of the Higgs potential, even at the tree level, is non-trivial [60]. Once a minimum is found, the mass of the lightest Higgs boson can be fine tuned to 125 GeV by small variations in the parameter . Setting masses for the additional scalars in the TeV range insures that the mixing with the lightest Higgs boson is small, and thus this Higgs will obey mass [71] and signal bounds [72] consistent with the SM-like Higgs found at the LHC. Additional Higgs states, in particular the lightest pseudoscalar, being heavy, will also satisfy constraints from branching ratio [73].

The charges of the fields satisfy conditions arising the requirement of cancellation of gauge and gravitational anomalies. For instance, the charges for Higgs fields in the model are chosen so thatThe charge of the quark doublet is kept as a free parameter after the normalizationA complete list of conditions for anomaly cancellations in the model, and a choice of charge assignments of the SM and exotic quarks and leptons in the model can be found in [70].

Gauge boson masses and mixing

Through spontaneous breakdown of the group to the Higgs acquire the VEVsHere the first two VEVs are required to break the gauge symmetries of the SM, and the third to break . After symmetry breaking, one massless state (the photon) and two massive states (the and bosons which are not yet the physical eigenstates due to a non-zero mass mixing term, to be introduced below) arise as orthonormal combinations of , and gauge bosons. The and combine to form , the charged vector bosons in the model. Unlike in the MSSM, the boson is not a physical state by itself but mixes with the boson. This mass mixing term arises from the fact that the Higgs doublets are charged under each factor of , and the associated mass-squared matrix is given byin the basis, where the matrix elements arewhere . The physical neutral vector bosons, , are obtained by diagonalizing :whereis their mass mixing angle, andare their squared masses of the corresponding mass eigenstates. The collider searches plus various indirect observations require the – mixing angle to be at most a few times [15], where unavoidable model dependence arises from couplings. This bound requires either to be large enough (well in the range) or to be sufficiently suppressed by the vacuum configuration, that is, . Which of these options is realized depends on the charge assignments and the soft-breaking mass parameters in the Higgs sector. Having large term in Eq. 6 insures a small mixing angle.

We expand more on the reason for introducing the extra scalars and . In their absence, term is equal to the one given in Eq. 7 without the summation term at the end. Hence, the mass of boson will be determined by the charges under , and the vacuum expectation value of the singlet scalar and thus effectively proportional to the VEV . At the same time, the masses of charginos and neutralinos in the model would be proportional to the effective , and the two scales, barring serious fine-tuning of the coupling, are connected. Introducing the extra scalars, their VEVs will be added to the contribution to the term (eventually determining the mass of boson), while , determining the mass scale of charginos and neutralinos, remains unchanged. As strong constraints are imposed on the boson mass at LHC, while the electroweakinos can remain relatively light, decoupling these two scales is desirable.

Thus, in our investigations, the mass of the boson will be expected to be heavy, TeV. In this case, the mixing angle between and becomes rather small,3 and the mass eigenstates are almost identical to the the original gauge states . However, despite the smallness of the mixing angle , we keep it in our numerical analysis.

Neutralinos in the secluded model

While the chargino sector is the same as in MSSM, the neutralino content is significantly enlarged by the additional gaugino, and the four additional singlinos and . In the basis where and are the neutral gauge fermions of and , the neutralino mass matrix is given asand is diagonalized by =diag, The parameters introduced in the neutralino mass matrix elements in Eq. 11 are defined asThe gaugino mass parameters and are free parameters of the model, and we introduce the ratiosThese parameters, representing mixing of and gauginos, and mass parameter of the gaugino, measured relative to the gaugino mass parameter, will be seen to be important in scanning over the parameter space, as the underlying physics will be sensitive to their variation.

boson in the model at the current and future hadron colliders

We now proceed to the main analysis in this work, looking at the consequences of a heavy neutral gauge boson at the collider. As we shall see, and as found before, bosons satisfying all collider, cosmological and low energy constraints, do not offer promising prospects for observability at the present LHC, even operating at 3 ab. Thus, we will also analyze the prospects of observing a signal at the HE-LHC operating at 27 TeV as well as at the FCC-hh. As the parameter space is large, choosing realistic benchmarks is a more transparent method to show physics results than a scan. Some previous analyses of at present and future colliders exist, e.g., in [74], but in that case, the authors have considered an -inspired leptophobic model, constrained to yield light masses. While we share with that analysis a light chargino-neutralino sector (insured in our case, by the existence of a secluded model), ours is an analysis for signals from a heavy boson scenario.

benchmark points and relic density

In order to give definite predictions for the production and decay rates of the boson, we scan the parameter space for benchmark scenarios to showcase the salient points of the model.

The benchmark points chosen must obey five important conditions:

The parameters chosen had to insure the stability of the vacuum;

The points had to satisfy relic density constraints from WMAP of cold dark matter [75] for the LSP, assumed here to be the lightest neutralino;

The mass of the boson has to satisfy mass constraints from ATLAS and CMS, as discussed in the next subsection;

Of the parameter points satisfying the above two conditions, benchmarks were chosen to enhance the supersymmetric decay signals of the boson; and

In each scenario, the lightest Higgs boson is SM-like and has GeV.

Table 2

The parameters characterizing benchmarks BP1, BP2 and BP3 for the secluded model. The values of dimensionful parameters are given in GeV

Parameters	BP1	BP2	BP3
g1′	0.2	0.12	0.15
tanβ	1.345	1.198	1.175
QQ′	0.6	0.1	- 0.81
μ_eff	260	280	250
(h_ν,h_s,h¯_s)	(1.0, 0.739, 0.1)	(1.0, 0.7235, 0.1)	(1.0, 0.724, 0.1)
(A_s,A_s¯)	557.7	(557.7, 2200)	(557.7, 1200)
(v_s_1,v_s_2,v_s_3)	(8675, 8650, 8675)	(6675, 15,600, 14,675)	(12,100, 14,550, 14,500)
(M_Y~,M_W~,M_g~)	(- 200, 2000, 2500)	(- 760, 750, 2500)	(- 260, 300, 2500)
(R_Y^′,R_YY^′)	(5.0,4.8)	(1.0, 0.01)	(1.0, 0.01)
(M_ν~_e_R,M_ν~_μ_R,M_ν~_τ_R)	500	3000	500
(M_L_1,M_L_2,M_L_3)	520	450	200
(M_E_1,M_E_2,M_E_3)	450	2125	1700
(M_Q_1,M_Q_2,M_Q_3)	(2200, 2200, 2400)	(2200, 2200, 2400)	(2200, 2200, 2400)
(M_U_1,M_U_2,M_U_3)	(2200, 2200, 2500)	(2200, 2200, 2500)	(2200, 2200, 2500)
(M_D_1,M_D_2,M_D_3)	(2300, 2300, 2500)	(2300, 2300, 2500)	(2300, 2300, 2500)
(MSS_12,MSS_22,MS_1S_22)	(-9×10^6,-9×10^6,0)	(-9×10^6,-9×10^6,0)	(-9×10^6,-9×10^6,0)
(A_t,A_b)	(- 697.75, - 959.66)	(- 697.75, - 959.66)	(- 697.75, - 959.66)

To analyze the model, we used CalcHEP [77], micrOMEGAs [78], PYTHIA8 [79, 80], Delphes [81-83], and MadAnalysis [84] to prepare the model, calculate the mass spectrum and branching ratios, calculate the relic density, generate events and eventually carry out the simulation. Our goal was to find benchmarks that satisfy cosmological constraints on dark matter and satisfy collider constraints (the invisible width of the Z boson, limits on charged sparticle masses, charginos mass, first and second-generation squark masses, lightest Higgs boson mass, Br), Br), and various others), as outlined in [85, 86], and choose those exhibiting distinct decay features, while offering some promise for collider observability.

In general, the boson in this model can decay into all SM fermions, into supersymmetric particles: squark, slepton, sneutrino, neutralino, chargino, in addition to Higgs-boson pairs, W-boson pairs and ZH.

The three benchmark points, and all the parameters associated with them, are given in Table 2. We give VEVs, Yukawa couplings, trilinear couplings, mass ratios and mixings for the gauginos and soft scalar fermion mass parameters. The low value of is consistent with constraints from branching ratio [73]. For each benchmark scenario, the mass spectra for the supersymmetric partners obtained are given in Table 3. The mass of the additional boson is TeV, and consistent with the ATLAS [54] analyses on dilepton decays. As seen in Table 2, the VEVs of the additional scalars ( and ) are mostly taken above the TeV scale so that the mass bound is satisfied independent of the value of the chosen VEV of the scalar field S. For convenience, the parameters and are taken as free parameters and the VEV of S is determined using the relationThe differences between the benchmarks are the following. In BP1, the gaugino mass parameters , is large, while the light (left-handed) sneutrinos and sleptons have mass GeV and are approximately degenerate. The right-handed sneutrinos are slightly heavier. In BP2, the gauginos have intermediate mass parameters , is very small, while the masses of the light (left-handed) sneutrinos and light sleptons are degenerate and around 500 GeV. The heavy ones split from the light sector significantly and can have much larger masses (up to GeV). In BP3, the gaugino mass parameters are both light , is very small, while, unlike the other two scenarios, the right-handed sneutrino masses are much lighter than those of the left-handed sneutrinos. The masses of the sleptons run in this range. The neutralino parameters affect the LSP and its composition, while the slepton and sneutrino masses affect branching ratios of into sfermions.

Table 3

The mass spectra (in GeV) for the supersymmetric sector and the relic density values of the benchmark points given in Table 2 for the Secluded

Masses	BP1	BP2	BP3	Bounds
m_Z^′	4250	4069	5195	3900 [76]
m_Hi0,i=1,…,6	(125.9, 543, 671, 1077, 4237, 17719)	(125.0, 557, 1148, 2418, 4171, 19151)	(125.3, 524, 1045, 1611, 5210, 22171)	m_H10=125.2 [76]
m_Ai0,i=1,…,4	(550, 719, 1012, 17718)	(563, 592, 20769, 19151)	(531, 572, 1882, 22170)	m_A10>93.4 [76]
m_χ~i0,i=1,…,5	(51, 167, 262, 312, 613)	(48, 269, 328, 762, 763)	(52, 195, 264, 303, 360)	m_χ~10>50 [76]
m_χ~i0,i=6,…,9	(1226, 2004, 4222, 4638)	(1170, 1740, 4047, 4237)	(1036, 1939, 4908, 5551)	–
(m_χ~1±,m_χ~2±)	(256, 2004)	(267, 763)	(192, 359)	m_χ~1±>103.5 [76]
m_H^±	540.9	554.9	522.4	m_H^±>80 [76]
(m_e~_L,m_μ~_L,m_τ~_1)	(503, 503, 457)	503	(1412, 1412, 473)	m_ℓ~>(82-107) [76]
(m_e~_R,m_μ~_R,m_τ~_2)	(457, 457, 503)	1850	(473, 473, 1412)	m_ℓ~>(82-107) [76]
(m_ν~_e,m_ν~_μ,m_ν~_τ)	501	501	1412	m_ν~_ℓ>41 [76]
(m_ν~_eR,m_ν~_μR,m_ν~_τR)	553	3472	645	–
Ω_DMh^2	0.117	0.121	0.119	0.111 [75]

The calculation of the relic density is performed importing the model files from CalcHEP [77] into the MicrOmegas package [78]. All the numbers obtained are within the range of the WMAP result obtained from the Sloan Digital Sky Survey [75]The relic density of dark matter is very sensitive to the parameter . The value of the relic density is shown in Table 3, where we also give explicit values for masses of the physical eigenstates in the Higgs and sparticle sectors. In addition, we checked that the benchmarks satisfy low energy data. For this, various flavor observables are calculated with the help of the packages SARAH [95, 96] and SPheno version 4.0.4 [97, 98]. The results, normalized with the corresponding SM values, are listed in Table 4. The values are all consistent with the current available data. In the same table, we give the corrections of the secluded to the SM values for the anomalous magnetic moments of electron, muon and tau, . The measured values for the first two indicate a departure from the SM, in opposite directions for the electron [87-89] and muon [90, 91]:In our benchmarks, the contributions are too small to saturate these differences, as these were chosen to instead yield interesting phenomenology. If we would aim to satisfy constraints on , we would choose values for the chargino, neutralino, slepton and sneutrino masses consistent with anomalous magnetic moment constraints. Since the chargino and neutralino sector masses are not directly connected to the boson mass, they would affect the phenomenology through decays. As the chargino-sneutrino loop is expected to be similar to the one in MSSM, the difference will arise from the lightest neutralino contribution. Its mass and composition are important, and so is the mass of the slepton. In addition, in both MSSM and models, the anomalous magnetic moment depends almost linearly on [99, 100]. A larger value of may increase the contributions to , but for our benchmarks, this is in conflict with flavor constraints.

Table 4

Leptonic anomalous moments corrections and flavor observables for each benchmark scenario considered in this study

Observable	BP1	BP2	BP3	Bounds
Δa_e	5.68×10^-16	1.14×10^-15	3.37×10^-16	-(8.7±3.6)×10^-13 [87–89]
Δa_μ	2.43×10^-11	4.86×10^-11	1.44×10^-11	(2.7±0.9)×10^-9 [90, 91]
Δa_τ	7.82×10^-9	1.12×10^-8	-4.78×10^-10	|Δa_τ|<1.75×10^-5 [92]
Br(B→X_sγ)Br(B→X_sγ)_SM	1.18	1.17	1.15	1.05± 0.11 [76]
Br(Bs,d0→μ^+μ^-)Br(Bs,d0→μ^+μ^-)_SM	1.09	1.11	1.10	0.83± 0.25 [76]
Br(B^+→τ^+ν_τ)Br(B^+→τ^+ν_τ)_SM	0.991	0.991	0.991	1.04± 0.34  [76]
ΔM_B_(s,d)/ΔMB_(s,d)SM	(1.10, 1.04)	(1.12, 1.04)	(1.12, 1.04)	(1.00±0.15,0.86±0.28) [76]
R_K/RKSM	1.00	1.00	1.00	1.00±0.17 [76]
ϵ_K/ϵKSM	1.00	1.00	1.00	0.99±0.18 [76, 93]
ΔM_K/ΔMKSM	1.00	1.00	1.00	1.24±0.16 [76, 94] ^a

In the theoretical calculation of the quoted value for ratio assumes contributions from both the so-called short distance and long distance physics. However, the latter part is not very reliable and needs improvement which might drive the value in either direction. The ratio becomes if only the short distance contribution is kept

Both ATLAS and CMS collaborations have searched for bosons. The assumption is that they are produced in pp, then decay into SM particles. The decay channels explored are jj [101-103], [104], [105, 106], [54, 107], [54, 107], [57, 108], [109, 110], and [110], within a variety of models with extended and SU(2) gauge groups. Of these channels, the leptonic decays and impose the most stringent constraints on the mass, normally TeV. However, all these analyses assumed non-supersymmetric scenarios. It has been shown that, including supersymmetry, these bounds can be reduced by GeV [63]. Hence, as we wish to explore the largest parameter space possible, we shall assume that TeV.

We calculated the branching ratios of the decaying into various final states for the three selected benchmark points, and show the results for the dominant ones in Table 5. As expected, branching ratios for decays into quarks (BP2 and BP3) or neutrinos (BP1) dominate over those into supersymmetric particles. The benchmarks were chosen for non-negligible decays into SUSY particles pairs, and are dominated by decays into sneutrinos and chargino pairs (BP1), chargino and neutralino pairs (BP2), and into slepton pairs (BP3). In Table 5 we also test the width/mass ratio. For all benchmarks considered, remains safely under 10%, justifying treating as a narrow resonance. The branching ratios of are very sensitive to variations in , the charge for the left-handed quark doublet. To highlight decays into supersymmetric channels, we fixed all parameters for the chosen benchmarks BP1, BP2 and BP3 as in Table 2, except , which is allowed to vary freely. In Fig. 1 we plot the branching ratios of the as a function of for each benchmark scenario. The particular choice for for each benchmark as given in Table 2, is obtained by requiring that some branching ratios into supersymmetric particles be maximal. We indicated these choices in each panel of Fig. 1 as a vertical grey line.

Table 5

Decay width (in GeV), width over mass ratios, and dominant branching ratios (in %) of boson decay channels for the three scenarios considered. The total branching ratios for decay modes with are also shown separately

Width[GeV]andbranchingratios[%]	BP1
Γ_Z^′	386
Γ_Z^′/M_Z^′[%]	9.0
BR(Z^′→∑_ℓν~_ℓ_Rν~_ℓ_R)	15.69
BR(Z^′→χ~1±χ~1∓)	2.93
BR(Z^′→χ~30χ~40)	2.09
BR(Z^′→∑_ℓν_ℓν¯_ℓ)	38.70
BR(Z^′→∑_qq_dq¯_d)	15.39
BR(Z^′→∑_qq_uq¯_u)	12.33
BR(Z^′→∑_ℓℓℓ¯)	4.08
∑_i[BR_i(Z^′→others)<1%]	8.79

Fig. 1

decay branching ratios as a function of for BP1 (left), BP2 (middle) and BP3 (right). For each benchmark we fix all other parameters as in Table 2 except for , which is allowed to vary. The choice of for each benchmark, chosen to maximize decays into supersymmetric particles, is indicated in each panel as a vertical grey line. Note that the branching ratios of into the SM fermions and their superpartners as well as into the right-handed scalar neutrinos are all depicted after summing over three generations

Typically, the SUSY decay modes include (i) or , (ii) , (iii) , (iv) , or , or etc. Such pure leptonic modes give rise to a signature consisting of charged-leptons and large missing energies, which are particularly well suited for observation at the LHC. To determine and classify all possible signals for the three scenarios we look into the decay topology of these particles. We classify signals according to the final number of leptons present in the signal events. The generic Feynman diagrams contributing dominantly to channels leading to signals with leptons and missing energy are shown in Fig. 2 for BP1 (a), BP2 (b) and BP3 (c).

Fig. 2

The generic Feynman diagrams for the decay channels of the in the secluded model for BP1, BP2, and BP3 in a–c, respectively

The events are generated at the partonic level with CalcHEP [77] and they are subsequently passed to PYTHIA8 [79, 80] for decay, showering and hadronization. Events which are saved in HepMC format [111] are then passed to MadAnalysis [84] for applying cuts and further data analysis. Delphes [81-83] is used for fast detector simulations. We simulated events for the , and signals at the LHC with center-of-mass energy, the HE-LHC with , and as well as the FCC-hh with . We used different PDF sets to model parton distributions for the colliders at different center of mass energies. While CTEQ6l1 PDF set [112] was used for the LHC case, the PDF set from the PDF4LHC15 collaboration [113] was used for both HE-LHC and FCC-hh. In the numerical study, for the calculation of signal significance, we have taken the integrated luminosities ab, 15 ab, and 30 ab for the HL-LHC, HE-LHC, and FCC-hh, respectively.

For each benchmark, we use two different formulas for the significance of the signals [114-116], where the first expression corresponds to the usual definition, and the second is more useful for smaller number of background events:where S is the number of signal events, B the number of background events, and the standard deviation for background events. We generated the events for the signal in each scenario, and we also simulated the SM background for the three benchmarks, separately for 14, 27 and 100 TeV.

We concentrate on leptonic final states, considered as golden channels in experimental searches at LHC. To exploit these features, this study will be focused on the decays of the boson into supersymmetric particles, leading to final states with leptons and missing energy, due to the presence of neutralinos or neutrinos. In the following, we present a study of decays into multileptonic final states for a given set of the secluded sector model parameters (BP1, BP2 and BP3), dividing our analysis into , and signals.

Multilepton analysis

In this analysis, for each final state, we impose cuts on the kinematical observables to suppress the SM background, as given in Table 6. Given the event topologies, stricter cuts on the leading lepton transverse momentum favor events with and . While the cuts on the angular variables and lepton separation remain the same, the kinematic cuts increase (in general), as expected going from 14 to 27 TeV and eventually to 100 TeV. Very stringent cuts are needed for the case of signal, and this is valid for the three center of mass energies but especially at 100 TeV. The final set of cuts are obtained by requiring to maximize the signal significance. We proceed in turn to analyze each of the final states, , and signals and discuss their potential for observability.

Table 6

The set of kinematical cuts used to isolate signal events from background

@14TeV	@27TeV	@100TeV
2ℓ+⧸E_T
|η|<2.5	|η|<2.5	|η|<2.5
ΔR_ℓℓ≥0.5	ΔR_ℓℓ≥0.5	ΔR_ℓℓ≥0.5
p_T(ℓ_1)>475 GeV	p_T(ℓ_1)>500 GeV	p_T(ℓ_1)>2000 GeV
p_T(ℓ_2)>50 GeV	p_T(ℓ_2)>300 GeV	p_T(ℓ_2)>1000 GeV
⧸E_T>50 GeV	⧸E_T>400 GeV	⧸E_T>2300 GeV
4ℓ+⧸E_T
|η|<2.5	|η|<2.5	|η|<2.5
ΔR_ℓℓ≥0.5	ΔR_ℓℓ≥0.5	ΔR_ℓℓ≥0.5
p_T(ℓ_1)>100 GeV	p_T(ℓ_1)>100 GeV	p_T(ℓ_1)>100 GeV
p_T(ℓ_2)>50 GeV	p_T(ℓ_2)>50 GeV	p_T(ℓ_2)>50 GeV
p_T(ℓ_3)>25 GeV	p_T(ℓ_3)>25 GeV	p_T(ℓ_3)>25 GeV
p_T(ℓ_4)>15 GeV	p_T(ℓ_4)>15 GeV	p_T(ℓ_4)>15 GeV
⧸E_T> 400 GeV	⧸E_T>350 GeV	⧸E_T>800 GeV
6ℓ+⧸E_T
|η|<2.5	|η|<2.5	|η|<2.5
ΔR_ℓℓ≥0.2	ΔR_ℓℓ≥0.2	ΔR_ℓℓ≥0.2
p_T(ℓ_1)>50 GeV	p_T(ℓ_1)>50 GeV	p_T(ℓ_1)>100 GeV
p_T(ℓ_2)>20 GeV	p_T(ℓ_2)>20 GeV	p_T(ℓ_2)>50 GeV
p_T(ℓ_3)>20 GeV	p_T(ℓ_3)>20 GeV	p_T(ℓ_3)>20 GeV
p_T(ℓ_4)>20 GeV	p_T(ℓ_4)>20 GeV	p_T(ℓ_4)>20 GeV
p_T(ℓ_5)>10 GeV	p_T(ℓ_5)>10 GeV	p_T(ℓ_5)>15 GeV
p_T(ℓ_6)>5 GeV	p_T(ℓ_6)>5 GeV	p_T(ℓ_6)>5 GeV
⧸E_T>100 GeV	⧸E_T>100 GeV	⧸E_T>100 GeV

Two lepton signal:

The main decay modes of giving rise to dilepton final states are:In the following figures, we first show the relevant kinematic variables for signals and background at 14, 27 and 100 TeV, before any cuts were imposed. We plot the differential cross-section, normalized to unity, with individual bin contents divided by the sum of all the data in the available bins. This way, the uncertainties are correlated across the bins, such that the uncertainties on the total integrated luminosity cancel. The resulting normalized differential fiducial cross-section is plotted as a function of various representative kinematic variables [117, 118]. Let us define the variable transverse mass [119]for a system composed of a lepton and the invisible transverse momentum available in each event. Here . We show, in Fig. 3, of the leading lepton (top panels) and next-to-leading lepton (bottom panels), and in Fig. 4, the leading lepton transverse momentum (top panels) and next-to-leading lepton transverse momentum (bottom panels). The left-hand side panels in both figures correspond to signals and backgrounds at 14 TeV, the middle at 27 TeV, and the right-hand side for 100 TeV.

Fig. 3

The transverse mass , for the leading lepton (top panels) and next-to-leading lepton (bottom panels), for the signal and background in the scenario. (Left-hand) signals and backgrounds at 14 TeV, (middle) at 27 TeV, and (right-hand) at 100 TeV. The main backgrounds (di-bosons) are indicated in solid lines while the signals are plotted in dotted lines: green for BP1, black for BP2 and orange for BP3

Fig. 4

The transverse momentum , for the leading lepton (top panels) and next-to-leading lepton (bottom panels), for the signal and background in the scenario. (Left-hand) signals and backgrounds at 14 TeV, (middle) at 27 TeV, and (right-hand) at 100 TeV. The main backgrounds (di-bosons) are indicated in solid lines while the signals are plotted in dotted lines: green for BP1, black for BP2 and orange for BP3

The main backgrounds (di-bosons for ) are indicated in solid lines while the signals are plotted in dotted lines, color coded: green for BP1, black for BP2 and orange for BP3. For our analysis we have checked cross sections of other potential backgrounds (tri-bosons) for , such as WWZ and ZZZ and found that their cross sections at 14 TeV are about three orders of magnitude smaller than those of the di-bosons. Such tri-boson background suppressions hold also at 27 TeV and 100 TeV center of mass energies. Thus they are subdominant and ignored in our analysis. We repeated the background calculations at each energy for the other signals (namely, and ) and only the dominant backgrounds, tri-bosons for and ZZZZ for , are kept in this analysis. As our signal is leptonic only, we selected certain number of isolated leptons + MET and vetoed any jets to avoid consistently signals and backgrounds due to QCD interactions. As an extensive background calculation is beyond the scope of our work, we restricted ourselves with backgrounds without jets, thus neglecting any contribution from soft jets, or jets which will not survive cuts.

Both of these graphs show clearly that at large and , the signal dominates the background, and the graphs justify our choice of kinematic cuts. Table 7 gives values for signal and background cross-sections after each cut. We also show the signal significance, for both and , for each benchmark, at proposed total integrated luminosity: for 14 TeV at ab, for 27 TeV at ab and for 100 TeV at ab. While BP3 appears to be most promising, the significance for all benchmarks at 14 TeV is very low, dispelling any hope for observing the boson in the final state. This is not surprising, and in complete agreement with other findings [63, 120]. However, at 27 TeV the cuts imposed are especially effective for BP3, and we obtain large significances for both and . At 100 TeV, the cuts imposed to reduce the background wiped out the signal, and we were not able to gain any predictable features. We had here two available options: either loose all the signal events due to the stringent cuts applied, or allow the signals to be overwhelmed by large background events.

Table 7

Signal selection strategy and cuts imposed in the scenario at 14, 27 and 100 TeV. We give the cross-section for background and benchmark scenarios in fb. Statistical significances and of signal are given for each energy

2ℓ+⧸E_T@14TeV	Background [fb]	BP1[fb]	BP2[fb]	BP3[fb]
Nocut	7.15×10^2	2.26×10^-2	8.25×10^-3	9.11×10^-3
|η_i|<2.5,ΔR_12≥0.5	6.99×10^2	1.14×10^-2	4.54×10^-3	9.06×10^-3
p_T(ℓ_1)>475 GeV	7.13×10^-2	1.00×10^-3	1.61×10^-3	6.29×10^-3
p_T(ℓ_2)>50 GeV	7.13×10^-2	1.52×10^-4	1.28×10^-3	6.01×10^-3
⧸E_T>50GeV	7.10×10^-2	1.5×10^-4	1.3×10^-3	5.9×10^-3
Significance:σ_A	L=3ab^-1	0.031σ	0.26σ	1.2σ
σ_B		0.022σ	0.18σ	0.85σ
2ℓ+⧸E_T@27TeV
Nocut	1.59×10^3	3.54×10^-1	1.25×10^-1	1.94×10^-1
|η_i|<2.5,ΔR_12≥0.5	1.55×10^3	1.52×10^-1	6.43×10^-2	1.94×10^-1
p_T(ℓ_1)>500 GeV	2.14	1.54×10^-2	1.18×10^-2	1.6×10^-1
p_T(ℓ_2)>300 GeV	9.22×10^-1	3.35×10^-4	9.76×10^-3	1.12×10^-1
⧸E_T>400 GeV	1.8×10^-1	1.0×10^-4	4.4×10^-3	6.0×10^-2
Significance:σ_A	L=15ab^-1	0.029σ	1.2σ	14.9σ
σ_B		0.021σ	0.88σ	11.3σ
2ℓ+⧸E_T@100TeV
Nocut	1.89×10^4	9.11	3.23	6.52
|η_i|<2.5,ΔR_12≥0.5	8.41×10^3	3.36	6.62×10^-1	5.75
p_T(ℓ_1)>2000 GeV	2.01×10^1	0.0	1.62×10^-3	5.15×10^-2
p_T(ℓ_2)>1000 GeV	8.94	0.0	1.2×10^-4	1.5×10^-2
⧸E_T>2300GeV	8.1×10^-1	0.0	0.0	3.5×10^-4
Significance:σ_A	L=30ab^-1	0.0σ	0.0σ	0.067σ
σ_B		0.0σ	0.0σ	0.048σ

Four lepton signal:

The main decay modes of the boson yielding signals areIn the following figures, we first plot the relevant kinematic variables signals and background at 14, 27 and 100 TeV, before any cuts were imposed. We show, in Fig. 5, the total missing transverse energy . For each figure we plot, in the left-hand columns, signals and backgrounds at 14 TeV, the middle columns at 27 TeV, and right-hand columns, at 100 TeV. Both these plots indicate that the signal over background, before cuts, seems most promising at higher energies, and that is overall a better variable to differentiate between signal and background. We also note that, as for the case, the scenario BP3 is the most promising.

Fig. 5

The total missing energy for the signal and background for the signal. (Left-hand) signals and backgrounds at 14 TeV, (middle) at 27 TeV, and (right-hand) at 100 TeV. The main backgrounds (three-bosons) are indicated in solid lines while the signals are plotted in dotted lines: green for BP1, and orange for BP3

In Fig. 6 we plot the transverse mass for the leading lepton (top row panels) and next-to-leading lepton (bottom row panels) for the signal and background for the signal. In Fig. 7, we show the transverse momenta of the leading lepton (top row) and next-to-leading lepton (bottom), for the signal and background. The transverse momenta of the third and fourth leptons have similar distributions to the next-to leading lepton, and thus we do not plot them. The left-hand side panels in both figures correspond to signals and backgrounds at 14 TeV, the middle at 27 TeV, and the right-hand side for 100 TeV. The main backgrounds (three-bosons) are indicated in solid lines while the signals are plotted in dotted lines: green for BP1, and orange for BP3.4 These graphs show clearly that at large and , the signal dominates the background, and the graphs justify our choice of kinematic cuts. The distribution is rather similar for the leading leptons, and the observable is better than the at distinguishing signals from the backgrounds. In all the plots, the signal from BP3 is most promising, particularly at high energy/momenta.

Fig. 6

The transverse mass , for the leading lepton (top row) and next-to-leading lepton (bottom row), for the signal and background for the signal. (Left-hand) signals and backgrounds at 14 TeV, (middle) at 27 TeV, and (right-hand) at 100 TeV. The main backgrounds (three-bosons) are indicated in solid lines while the signals are plotted in dotted lines: green for BP1, and orange for BP3

Fig. 7

The transverse momentum , for the leading lepton (top panels) and next-to-leading lepton (bottom panels), for the signal and background for the signal. (Left-hand) signals and backgrounds at 14 TeV, (middle) at 27 TeV, and (right-hand) at 100 TeV. The main backgrounds (three-bosons) are indicated in solid lines while the signals are plotted in dotted lines: green for BP1, and orange for BP3

Table 8 gives values for signal and background cross sections after each cut. We also show the signal significance, for both and , for each benchmark, at proposed total integrated luminosity: for 14 TeV at ab, for 27 TeV at ab and for 100 TeV at ab. Unlike the case of , we keep the cuts constant for different centre-of-mass energies. Again, the significance for all benchmarks at 14 TeV for observing the boson in the final state is low. However, at 27 TeV both BP1 and BP3 show some promise, and we obtain large significances of or more for . At 100 TeV, though there are many uncertainties and unknowns, and our results should be interpreted as estimates only, both BP1 and BP3 show significant promise for observability.

Table 8

Signal selection strategy and cuts imposed in the scenario at 14, 27 and 100 TeV. We give the cross-section for background and benchmark scenarios in fb. Statistical significances and of signal are given for each energy

4ℓ+⧸E_T@14TeV	Background [fb]	BP1[fb]	BP3[fb]
NoCut	3.45×10^-2	1.99×10^-3	1.86×10^-3
|η_i|<2.5,ΔR_12≥0.5	2.89×10^-2	5.31×10^-4	4.63×10^-4
p_T(ℓ_1)>100 GeV	1.7×10^-2	5.1×10^-4	4.61×10^-4
p_T(ℓ_2)>50 GeV	1.58×10^-2	5.02×10^-4	4.61×10^-4
p_T(ℓ_3)>25 GeV	1.52×10^-2	4.87×10^-4	4.59×10^-4
p_T(ℓ_4)>15 GeV	1.4×10^-2	4.44×10^-4	4.31×10^-4
⧸E_T>400GeV	1.0×10^-4	5.5×10^-5	1.3×10^-4
Significance:σ_A	L=3ab^-1	0.24σ	0.46σ
σ_B		0.19σ	0.38σ
4ℓ+⧸E_T@27TeV
NoCut	9.59×10^-2	3.03×10^-2	3.86×10^-2
|η_i|<2.5,ΔR_12≥0.5	7.52×10^-2	5.9×10^-3	4.14×10^-3
p_T(ℓ_1)>100 GeV	4.7×10^-2	5.74×10^-3	4.13×10^-3
p_T(ℓ_2)>50 GeV	4.39×10^-2	5.68×10^-3	4.13×10^-3
p_T(ℓ_3)>25 GeV	4.2×10^-2	5.52×10^-3	4.12×10^-3
p_T(ℓ_4)>15 GeV	3.86×10^-2	5.08×10^-3	3.9×10^-3
⧸E_T>350GeV	1.1×10^-3	1.2×10^-3	2.2×10^-3
Significance:σ_A	L=15ab^-1	3.1σ	4.6σ
σ_B		2.6σ	4.1σ
4ℓ+⧸E_T@100TeV
NoCut	2.48	9.38×10^-1	1.64
|η_i|<2.5,ΔR_12≥0.5	6.56×10^-1	1.33×10^-1	1.16×10^-1
p_T(ℓ_1)>100 GeV	4.46×10^-1	1.3×10^-1	1.15×10^-1
p_T(ℓ_2)>50 GeV	4.1×10^-1	1.29×10^-1	1.14×10^-1
p_T(ℓ_3)>25 GeV	3.83×10^-1	1.22×10^-1	1.01×10^-1
p_T(ℓ_4)>15 GeV	3.1×10^-1	9.9×10^-2	8.44×10^-2
⧸E_T>800GeV	7.1×10^-3	8.0×10^-3	2.3×10^-2
Significance:σ_A	L=30ab^-1	11.0σ	23.0σ
σ_B		9.4σ	22.0σ

Six lepton signal:

The dominant decay mode of the gauge boson, yielding a signal isQuite clearly, the signal requires that this decay have a non-negligible branching ratio, which occurs for the BP3 scenario, as seen from Table 5, where . Thus as expected, this will be the only signal of interest for the signal. In Fig. 8 we plot the missing transverse energy , for the signal and background for the signal: (left-hand) signals and backgrounds at 14 TeV, (middle) at 27 TeV, and (right-hand) at 100 TeV. We show the main (four-bosons) backgrounds in solid lines while the the signal BP3 is given by a dotted green line.

Fig. 8

The missing transverse energy , for the signal and background for the signal. (Left-hand) signals and backgrounds at 14 TeV, (middle) at 27 TeV, and (right-hand) at 100 TeV. The main backgrounds (four-bosons) are indicated in solid lines while the dotted green line represents the signal BP3

In Fig. 9 we plot transverse momentum for the leading lepton (top panels), the next-to-leading lepton (second panels), and the next-to-next to leading lepton (bottom panels) for the signal and background for the signal. The transverse momentum plots for the fourth, fifth, and sixth lepton are similar and we do not show them here. As before the left side panels indicate signals and backgrounds at 14 TeV, the middle at 27 TeV, and the right side at 100 TeV. The main backgrounds (four-bosons) are indicated in solid lines while the dotted green line represent the signal BP3. As expected, the leading lepton distribution is most promising in distinguishing this signal from background, with other lepton distributions slightly less so.

Fig. 9

The transverse momentum , for the leading lepton (top panels), next-to-leading lepton (second panels), and next-to-next to leading lepton (bottom panels) for the signal and background for the signal. (Left-hand) signals and backgrounds at 14 TeV, (middle) at 27 TeV, and (right-hand) at 100 TeV. The main backgrounds (four-bosons) are indicated in solid lines while the dotted green line represents the signal BP3

Similar to the and signals, the effects of various cuts on cross sections are listed in Table 9. We also show the signal significance, for both and , for BP3, at proposed total integrated luminosity: for 14 TeV at ab, for 27 TeV at ab and for 100 TeV at ab. The signal significance can be around at 27 TeV and even greater than at 100 TeV.

Table 9

Signal selection strategy and cuts imposed in the scenario at 14, 27 and 100 TeV. We give the cross-section for background and benchmark scenarios in fb. Statistical significances and of signal are given for each energy

6ℓ+⧸E_T@14 TeV	Background [fb]	BP3 [fb]
No Cut	3.12×10^-7	8.1×10^-5
|η_i|<2.5,ΔR_12≥0.2	2.91×10^-7	4.37×10^-5
	2.91×10^-7	4.37×10^-5
⧸E_T>100GeV	1.2×10^-7	3.8×10^-5
Significance:σ_A	L=3ab^-1	0.34σ
σ_B		0.40σ
6ℓ+⧸E_T@27TeV
No Cut	8.4×10^-7	1.71×10^-3
|η_i|<2.5,ΔR_12≥0.2	7.87×10^-7	7.03×10^-4
	7.80×10^-7	7.03×10^-4
⧸E_T>100GeV	3.9×10^-7	6.5×10^-4
Significance:σ_A	L=15ab^-1	3.1σ
σ_B		3.7σ
6ℓ+⧸E_T@100TeV
No Cut	2.27×10^-5	9.52×10^-2
|η_i|<2.5,ΔR_12≥0.2	1.04×10^-5	4.92×10^-3
p_T(ℓ_1)>100 GeV	8.27×10^-6	4.92×10^-3
p_T(ℓ_2)>50 GeV , p_T(ℓ_3)>20 GeV	8.17×10^-6	4.92×10^-3
p_T(ℓ_4)>20 GeV	8.07×10^-6	4.66×10^-3
p_T(ℓ_5)>15 GeV	7.5×10^-6	3.65×10^-3
p_T(ℓ_6)>5 GeV	6.77×10^-6	2.88×10^-3
⧸E_T>100GeV	3.9×10^-6	2.3×10^-3
Significance:σ_A	L=30ab^-1	8.3σ
σ_B		9.7σ

Summary and conclusion

In this work we have analyzed the LHC, HE-LHC and FCC-hh discovery prospects of a new neutral gauge boson () through its supersymmetric decay modes, as a promising signal for supersymmetry in an extended gauge structure. We have assumed that the originates from an additional symmetry, and we decoupled its mass scale from that of supersymmetry breaking, assumed to be the same as the scale of breaking . This allows the boson to be heavy, as indicated by lower limits TeV from the experimental searches, while the electroweakinos remain light. For this, we relied on the secluded model, where three additional singlet superfields are added to the model. Unlike the VEV of the singlet scalar field which breaks and affects the mass of supersymmetric particles, the VEVs of the additional scalars enter only in the expression for the mass.

This scenario provides a fertile ground for analyzing decays into chargino, neutralino, slepton and sneutrino pairs. As LHC is particularly sensitive to events containing one or more leptons, we looked for production of followed by decays into multileptons plus missing energy. For this, we devised three benchmarks (BP1, BP2 and BP3) where branching ratios into some supersymmetric particles are enhanced. For instance in BP1, the decay into sneutrinos and chargino is enhanced, in BP2 it is the decay into chargino and neutralino pairs, while for BP3, the decay into right sleptons and lightest staus is important. The benchmarks are chosen also so they satisfy collider and relic density constraints.

We proceed by analyzing the observability of the signals at = 14, 27 and 100 TeV, looking separately at , and signals. Throughout our benchmarks, the ratio remains under 10%, so we can treat as a narrow resonance. For each signal, we perform a Monte Carlo simulation of the signal and background, and devise cuts to disproportionately suppress the latter. We present the results before and after the cuts, and calculate the significance in two different ways.

Overall, our findings indicate that the probability of observing through supersymmetric decays at 14 TeV is not good, even at high total integrated luminosity ab. This occurs across all , and signals and for all benchmarks. This confirms past analyses for [60, 74], which indicated that, unless is leptophobic, and thus much lighter, the signal significance is small. However, that is not necessarily so at 27 or 100 TeV and across all signals. At 27 TeV, benchmark BP3 gives a significance well above 5 in signal. While for the signal, we obtain significances of 3–4 for both BP1 and BP3, and much higher at 100 TeV. For the signal, only BP3, were the decay into right sleptons is important, gives any significant contributions. The significance at 27 TeV with total integrated luminosity ab is 3–4, and can reach 8–9 at 100 TeV with total integrated luminosity ab.

Of course, analyses at 27 TeV are plagued by uncertainties, and those at 100 TeV can be interpreted as merely estimates. However, our analysis shows that HE/HL-LHC and FCC-hh can be promising grounds for observing consequences of both supersymmetry and extended gauge symmetry, of which an additional neutral gauge boson is one of the simplest examples. A heavy boson, accompanied by a light electroweakino sector would also be indicative of a model with a secluded sector, as this set-up facilitates the splitting of the two scales.