Signals of new gauge bosons in gauged two higgs doublet model

Recently a gauged two Higgs doublet model, in which the two Higgs doublets are embedded into the fundamental representation of an extra local SU(2)H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$SU(2)_H$$\end{document} group, was constructed. Both the new gauge bosons Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document} and W′(p,m)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W^{\prime (p,m)}$$\end{document} are electrically neutral. While Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document} can be singly produced at colliders, W′(p,m)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W^{\prime (p,m)}$$\end{document}, which is heavier, must be pair produced. We explore the constraints of Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document} using the current Drell–Yan type data from the Large Hadron Collider. Anticipating optimistically that Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document} can be discovered via the clean Drell–Yan type signals at a high luminosity upgrade of the collider, we explore the detectability of extra heavy fermions in the model via the two leptons/jets plus missing transverse energy signals from the exotic decay modes of Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document}. For the W′(p,m)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W^{\prime (p,m)}$$\end{document} pair production in a future 100 TeV proton–proton collider, we demonstrate certain kinematical distributions for the two/four leptons plus missing energy signals have features distinguishable from the Standard Model background. In addition, comparisons of these kinematical distributions between the gauged two Higgs doublet model and the littlest Higgs model with T-parity, the latter of which can give rise to the same signals with competitive if not larger cross sections, are also presented.


Introduction
With the discovery of the 125 GeV scalar boson at the Large Hadron Collider (LHC), the Standard Model (SM) with just one Higgs doublet has now been generally accepted as the standard theory or framework to describe the fundamental strong and electroweak interactions for three generations of elementary particles of quarks and leptons. An extended Higgs sector, however, is often used to address various theoretical puzzles like the neutrino masses, dark matter (DM), a e-mail: timluyu@hotmail.com the matter-antimatter asymmetry, and the hierarchy problem, which remain unexplained in this standard framework. Perhaps the general two Higgs doublet model (2HDM), in particular in the context of supersymmetric theories, is studied most in the literature. Due to its diverse variations, 2HDM has been used as a prototype to address the aforementioned theoretical issues. For reviews of 2HDM and its supersymmetric version, see for example [1][2][3][4]. One of the interesting 2HDM variants is the inert Higgs doublet model [5][6][7][8], in which the neutral component of the second Higgs doublet can be a dark matter candidate due to a discrete Z 2 symmetry imposed on the scalar potential. Detailed analysis for the inert Higgs doublet model can be found for example in [9][10][11][12]. The origin of multiple inert Higgs doublets in the context of grand unification has been addressed in [13].
In a recent work [14], we have proposed a novel model, dubbed the Gauged Two Higgs Doublet Model (G2HDM), where the two Higgs doublets H 1 (2) H , are also introduced to render the model free of gauge anomalies. In this model, an inert Higgs doublet can be naturally realized without imposing the ad hoc Z 2 symmetry mentioned above to accommodate a DM candidate. Flavor changing neutral currents naturally are also absent at tree level. We note that it is widely believed that global symmetry (whether it is discrete or continuous) may be strongly violated by gravitational effects [15,16]. Thus from an effective field theory point of view it is desirable to embed discrete symmetries into local gauge symmetries below the Planck scale [15]. Indeed, besides Z 2 , non-abelian discrete flavor groups like {A 4 , S 4 , A 5 }, {Q 6 , T , O , I }, and  (27), P SL(2, 7)} can be minimally embedded into SO (3), SU (2), and SU (3), respectively [17]. In G2HDM, a distinctive feature is all the SU (2) H gauge bosons Z and W ( p,m) are electrically neutral which is not the case for the Left-Right Symmetry Model (LRSM) [18,19], the Littlest Higgs Model with T-parity (LHT) [20][21][22], and the original Twin Higgs Model (THM) [23,24]. Naturally one might ask: how do we distinguish the G2HDM Z and W ( p,m) from other gauge boson impostors which also arise from the non-abelian group SU (2)? In terms of collider searches, some of the new gauge bosons from the aforementioned models can never be singly produced due to the gauge symmetry involved. For instance, the W ( p,m) in G2HDM and W ± H in LHT may not be singly produced at the LHC. Generally speaking, the Z and W ( p,m) as well as their impostors are short-lived and it is not possible to identify the new gauge bosons using the tracking or displaced vertex techniques designed for long-lived particles. Hence, in addition to their production cross sections, detailed kinematics distributions have to be involved for making the differentiation.
If new gauge bosons can be singly produced, they will be stringently constrained by exotic searches from the LHC [25,26] due to large (resonant) production cross sections. The latest LHC 13 TeV Z resonance searches based on the channels of dilepton [27,28], dijet [25,29], b-quark pair [30], t-quark pair [31], and other bosonic final states [32][33][34] have recently been released. Among these searches, the cleanest dilepton channels yield the most strong constraint on the Z coupling to SM fermions in the light of the small background. Moreover, from the total electric charge of the decay products it is straightforward to tell singly produced charged bosons from neutral ones.
The current bound from the LHC dilepton searches [27,28] for Z R in the minimal LRSM, assuming m Z R = 1.7 m W R and g L = g R , is m Z R > 3.2 TeV [35]. Similarly for Z and W in the Left-Right THM (LRTH) [36] with m W = m Z √ cos 2θ w / cos θ w and g L = g R , one obtains the limit m Z > 3.36 TeV. In some variants of the THM (see e.g. Ref. [37]), where the SU (2) symmetry is doubled and the U (1) symmetry of the twin sector becomes a global symmetry, the exotic W boson can only be doubly produced. In this case, the exotic W might behave similarly to W ( p,m) in G2HDM in terms of collider signatures. To distinguish them, one will need more information such as the new scalar and fermion mass spectra, which are quite different in the two models.
Suppose a new Z is observed via resonance searches at the LHC. We would like to know if this new Z belongs to a new abelian U (1) or a member of a new SU (2). In order to confirm the existence of G2HDM, the next step is to discover the neutral W ( p,m) . In G2HDM, unlike Z , the W ( p,m) do not couple bilinearly to the SM quarks. Thus Z can be singly produced via quark-antiquark annihilations while W ( p,m) must be pair produced via exchange of new heavy fermions or Z . For W ( p,m) produced in pairs, more information like detailed kinematical distributions other than the production cross sections have to be involved so as to make distinguishable signatures from say the W ± H pair in LHT.
Besides the new gauge bosons, in order to give masses to all fermions the scalar sector is also enlarged beyond the two Higgs doublets with one extra doublet H and one extra triplet H of SU (2) H , both of which are singlets under SU (2) L . The particle content of the G2HDM [14] is summarized in Table 1 together with their quantum numbers.
In this work, we will focus on two benchmark mass spectra (Spectrum A and Spectrum B) of the G2HDM for our collider studies. Spectrum A contains heavy and decoupled new quarks, while Spectrum B comprises relatively light new quarks. For all scenarios, new leptons are assumed to be lighter than the additional gauge bosons of interest. Due to the fact Z couples to SM quarks and can be singly produced at the LHC, we first update the bounds on the SU (2) H gauge coupling g H as a function of the Z mass m Z by using the newly released results of the dilepton and dijet searches from the LHC. Next, Z exotic decays into new heavy fermions followed by decays into SM fermions are investigated at the 14 TeV High Luminosity LHC (HL-LHC) and bounds from LHC searches on supersymmetric particles can be applied with simplified assumptions. Then, for the neutral W ( p,m) in G2HDM we propose searching for two channels: two leptons and four leptons with missing transverse energy. We shall demonstrate that the pair production of W ( p,m) can feature quite distinctive kinematical distributions from the W H pair in LHT, which will be chosen as a representative model for comparisons since W can only be pair produced in both models.
This rest of this paper is laid out as follows. In Sect. 2, we briefly review the G2HDM and spell out the relevant gauge interactions for collider searches of interest. In Sect. 3, we discuss the methodology employed in the collider simulations. In Sect. 4, we revisit Z direct search limits from the latest 13TeV LHC data and explore some of its exotic decay channels at the HL-LHC. In Sect. 5, signatures for W at a future 100 TeV proton-proton collider are scrutinized in detail and compared with those from LHT. We summarize our findings and conclude in Sect. 6. For convenience, we also present the scalar potential of G2HDM and the associated scalar mass spectra in two appendices. More details of the scalar sector of G2HDM can be found in [38].

G2HDM gauge interactions
In this section, we give a brief review of G2HDM, focusing on gauge interactions that are relevant to our study of collider searches. The particle contents summarized in Table 1 have the minimal set of new heavy chiral fermions required for anomaly cancellation and new scalars for facilitating spontaneous electroweak symmetry breaking, as proposed in [14].
As mentioned earlier, the two SU (2) L Higgs doublets H 1 and H 2 are embedded into a doublet H under a non-abelian SU (2) H gauge group. H is also charged under an additional gauged abelian group U (1) X . To provide masses to the additional gauge bosons, we introduce an SU (2) H scalar triplet H and doublet H (both are singlets under the SM gauge group). The vacuum expectation value (vev) of the triplet H not only breaks SU (2) H spontaneously, but it also trig-gers the electroweak symmetry breaking by inducing a vev to the first SU (2) L doublet H 1 , which is identified as the SM Higgs doublet. In contrast, the second Higgs doublet H 2 does not obtain a vev and its neutral component could be the DM candidate. As shown in [14], DM stability is protected by the SU (2) H symmetry and Lorentz invariance. In other words, an inert Higgs doublet H 2 emerges naturally in G2HDM without resorting to the discrete Z 2 symmetry. 1 We specify the most general and renormalizable scalar potential With a non-vanishing 2 , the four Dirac fields d H , u H , e H , and ν H acquire a mass of y d 2 , y u 2 , y e 2 , and y ν 2 , respectively. On the other hand, the SM quarks and leptons obtain their masses from the vev of H 1 via the Yukawa couplings, withH 1,2 = iτ 2 H * 1,2 . Note that in both Yukawa couplings given in (1) and (2) only the SM Higgs doublet H 1 couples bilinearly with the SM fermions. Thus, FCNC interactions for the SM fermions naturally are absent at tree level in G2HDM. It also implies the new heavy fermions can decay into SM fermions plus DM via the Yukawa couplings in (2). For instance, f H R → f L H 0 * 2 where H 0 * 2 is a DM candidate and is manifest as the missing transverse energy. We note that the absence of FCNC interactions in 2HDM by embedding the discrete Z 2 symmetry into an extra U (1) has been studied in [40][41][42][43].
There are SU (2) H gauge bosons, W ( p,m) and W 3 , and the U (1) X gauge boson X , apart from the SM ones. Due to the symmetry breaking pattern, W ( p,m) will not mix with the SM counterparts but W 3 and X mix with the SM SU (2) L W 3 and U (1) Y Y gauge boson. In this setup, besides the SM massless photon corresponding to the unbroken generator Q = T 3 L + Y , there exists a massless dark photon corresponding to the unbroken generator Q D = 4 cos 2 θ w T 3 generator, and θ w is the Weinberg angle. Such a massless dark photon could be cosmologically problematic. To circumvent the problem, one can resort to the Stueckelberg mechanism to give a mass to the U (1) X gauge boson as in [44][45][46][47]. One could take this mass to be large enough so that X is decoupled from the particle spectra. Another way out is to treat U (1) X as a global symmetry as was proposed in [48] and adopted in [14] as well. We will follow the same strategy in what follows.
In this case, after diagonalizing the mass matrix of Y , W 3 and W 3 , one obtains massless γ , massive Z and Z . Furthermore, the Z -Z mixing is constrained to be of order 10 −3 for TeV Z because of the electroweak precision measurements [49]. As a consequence, impacts of the mixing are numerically negligible and will be ignored. The resulting SU (2) H gauge boson mass spectrum is where . Note that W ( p,m) is always heavier than Z in G2HDM.
As the SM right-handed fermions as well as the new fermions are charged under SU (2) H , they couple to the W ( p,m) and Z bosons. The relevant gauge interactions without the Z -Z mixing read Here L(W ) and L(γ ) refer to the charged current mediated by the W boson and the electric current by the photon γ , respectively, where Q f is the corresponding fermion electric charge in units of e. L represents (electrically) neutral current interactions of the massive bosons, Z , Z and W ( p,m) (for demonstration, only the lepton sector is shown, but it is straightforward to include the quark sector): where and The current interactions in L(W ( p,m) ) and L(Z ) will dictate how W ( p,m) and Z decay into SM and heavy fermions, and they determine which final states one should look into for collider searches.

Methodology
To simulate the total cross sections and various distributions for the relevant processes in the colliders, we will follow the standard protocol well established by many collider phenomenologists. We use FeynRules [50] to build up the model files for G2HDM and pass it to Madgraph5 [51] for the matrix element calculation and event generation. We simulate parton showering by using Pythia8.1 [52], and we employ Delphes3 [53] for detector simulations. Finally, the package MadAnalysis5 [54,55] is used to analyze the simulation data.
In the G2HDM, apart from the extra gauge bosons W ( p,m) and Z , additional heavy fermions have to be included to attain gauge invariant Yukawa couplings as explained above. To simplify the analysis, we assume two universal masses for the heavy fermions, one for leptons and the other for quarks. As a result, there are five relevant mass scales in our analysis, namely the masses of the dark matter particle H 0 * 2 , the heavy , and the two heavy gauge bosons W ( p,m) and Z . In addition, the new charged fermions have to be heavier than 100 GeV, a constraint inferred from the combined analysis of the LEP2 run data by the four LEP collaborations [56]. We will study the Table 2 Branching ratios for different decay modes of Z with 1.5 ≤ m Z ≤ 3 TeV Here Q denotes 6 quark flavors (u, d, c, s, t, b) and L (ν) represents 3 lepton flavors (e (ν e ), μ (ν μ ), τ (ν τ )) Table 3 Branching ratios for different decay modes of To achieve m Z m W , one needs v 3v v based on Eq. (3). Note that this setup is different from the previous work [14], where v v v was assumed. Furthermore, for simplicity decays of the heavy gauge bosons into scalar Higgs pairs are presumed to be either kinematically forbidden or negligible. It is justified since all of the new scalars except for DM can be heavier than W and Z as displayed in the last table in Appendix B. Moreover, the coupling between the longitudinal components of W and Z and the DM can in principle be made small by varying the parameters in the scalar potential. In this way, the transverse components (whose coupling to DM is simply g H ) govern decays of W and Z into DM particles but this contribution to the decays is subleading compared to those of the heavy fermions in the final states, given the larger number of the new fermions in the model.
In both scenarios, the new heavy fermions are kinematically allowed to be produced by either Z or W ( p,m) decays. As a result, we propose searches for the new fermions as follows. 2 • For Spectrum A, the heavy charged leptons can be produced via pp → Z → L H L H and pp → W p W m → L H L L L H , and the corresponding final states will be (1) 2l • For Spectrum B, the new quark pairs can also be on-shell produced through pp → Z → Q H Q H , and thus the final states (1) 2 j + E T , (2) 2b + E T , and (3) 2t + E T will be considered. These processes are relevant to the dijet plus missing transverse energy searches for Z . Needless to say, the continuum contributions from QCD to the new quark pair production should be taken into account.
In Tables 2 and 3, we list the branching ratios for the Z and W ( p,m) decays, respectively, in the two scenarios. 3 The Q H Q H final state in the Z decay is kinematically allowed in Spectrum B, resulting in smaller partial decay widths into the SM Q Q and L + L − final states compared to Spectrum A. On the other hand, since W ( p,m) do not decay into SM fermion pairs, the opening of Q H Q + Q Q H final state in Spectrum B affects only the other exotic leptonic channels. These exotic decays can be phenomenologically interesting as we will see in the next section.
Note that Z and W ( p,m) can also decay into scalars, H , H and H , which are charged under SU (2) H . The branching fractions, however, depend on the scalar mixing parameters in the scalar potential [14] which can make the analysis Fig. 1 The Z constraints for G2HDM inferred from the latest ATLAS and CMS 13TeV results. The solid lines denote Spectrum A, while the dashed lines refer to Spectrum B. The main differences between the two scenarios are the branching fractions of Z into the SM quarks and leptons, as shown in Table 2 rather convoluted. As mentioned earlier, to simplify the analysis in this work, we neglect scalar final states and instead focus on the fermion channels at which the corresponding partial decay widths are simply fixed by the SU (2) H gauge coupling as well as the new heavy fermion masses.

Z searches at the LHC
In this section, we first present the Z constraints, derived from the ATLAS and CMS dijet and dilepton searches based on the 13 TeV data. Then we propose potential Z signatures from exotic decay searches which have smaller cross sections than direct Z searches but can be explored at the 14 TeV HL-LHC. In Sect. 5, we will investigate the W ( p,m) searches at a future proton-proton 100 TeV collider. As mentioned before, unlike Z , which can be singly created and probed directly by dilepton and dijet searches at the LHC, the heavier W ( p,m) must be produced pair-wise in the light of the SU (2) H gauge symmetry and therefore are less constrained.

Constraints on Z from current dilepton and dijet searches
In G2HDM, Z is always lighter than W and can be directly probed by resonance searches as mentioned above. Due to the null results of the direct searches, very stringent limits are imposed on any model of Z that directly couples to the SM fermions. For instance, the sequential SM with Z having the same couplings to SM fermions as the SM Z gauge boson is constrained to be heavier than 4 TeV [27,28].
Recently, the ATLAS and CMS collaborations have reported their updated results of Z resonance searches for channels of dilepton [27,28], dijet [25,29], b-quark pair [30], t-quark pair [31], and other bosonic final states [32][33][34] at 13 TeV. In the light of the irreducible QCD background at the LHC, the dilepton channel is the cleanest one to reconstruct the invariant mass of the final state particles and yields the most stringent constraints. In this work, we consider two major type of constraints: dilepton and dijet channels. We calculate the cross sections of pp → Z → l + l − /j j with the help of Madgraph5 [51] and compare them with the latest constraints from the LHC.
In Fig. 1, we present the exclusion regions of g H as a function of m Z . The solid lines correspond to Spectrum A and the dashed lines denote Spectrum B. The major discrepancies between the two scenarios are the branching ratios of the Z decay into SM quarks and leptons, as shown in Table 2. Compared with the previous constraints [14] obtained from the 8 TeV data, the improvement is about a factor of 2 in the region of 1.5 < m Z < 2.25 TeV. In addition, thanks to the higher center-of-mass energy √ s (8 → 13 TeV), the bounds for m Z > 2.25 TeV are significantly improved and become stronger than the LEP limits based on the σ (e + e − → l + l − ) measurements [14].
We note that the constraints on the Z mass in G2HDM is less stringent than the Z in LRSM or LRTHM where a discrete symmetry is imposed to equate the new gauge coupling to the SM SU (2) L one. The price to pay for G2HDM is, of course, a smaller g H .

Z exotic decays into heavy fermions
We now move on to the Z exotic decays which can shed light on the existence of exotic fermions in G2HDM. As noted before, the scalar decay channels of Z depend on details of the complicated scalar potential and hence are ignored in this work. On the other hand, the heavy fermion channels, which are governed by g H and the heavy fermion masses only, can easily be addressed.
Thinking forwardly and optimistically, one can envisage that a Z will be discovered by the direct searches of dilepton and dijet channels in the foreseeable future at HL-LHC. If so, the heavy fermions in G2HDM can also be probed via Z onshell decays if kinematically allowed. In order to perform a more general study for this purpose, we will temporarily relax the mass relations among the heavy fermions, Z , and DM for both Spectrum A and Spectrum B mentioned in Sect. 3. To be specific, in this section and only in this section, we will relax the fixed mass relation in Spectrum A to 2 m L H < m Z < Thus we will choose a nominal value of m DM = 50 GeV in the following analysis. These two modified spectra will be referred to as Spectrum A and Spectrum B in what follows.
Owing to the SU (2) H symmetry, final states of the exotic decay modes are quite similar to those used to search for supersymmetric (SUSY) particles with R-parity conservation in the context of the Minimal Supersymmetric Standard Model (MSSM). For example, for the slepton (l ± ) searches at the LHC the major process is 2l + E T channel, namely Similarly, in G2HDM, decays of Z into a pair of exotic fermions can also lead to the same final states: These two processes with the same final states exhibit analogous event topology, allowing us to apply the same kinematic cuts. 4 We calculate the cross sections for the Z exotic decays and then impose bounds from SUSY searches on these decays. The bounds on the heavy fermion masses can be mitigated in a scenario of the compressed mass spectrum: m (L H ,Q H ) m DM . In this case mono-X (X = γ , g, W , Z , · · · ) + E T , in particular the monojet + E T signal, can be used to search for DM as in the MSSM with the compressed mass spectrum. This scenario, however, will not be considered here as the mass difference m (L H ,Q H ) − m DM is taken to be large.
We now present the constraints on the Z exotic decays. With the modified spectra, we will be able to obtain contours of the production cross sections on the m Z -m f H plane and compare with the LHC limits on the SUSY particle searches.
• Spectrum A : 2 m L H < m Z < 2 m Q H and m DM = 50 GeV.
We concentrate on the following two channels: where l H = (e H , μ H ), and l = (e, μ). In Fig. 2, the contour plots for the cross sections of the processes in (11) and (12)  are shown in terms of σ/g 2 H to factor out the g H dependence. For a specific value of g H , one can simply rescale the contours by g 2 H whose limits for a given value of m Z have been presented in Fig. 1. The black, blue, and red contours correspond to σ/g 2 H = 0.1, 0.05, and 0.01 picobarn (pb) respectively, assuming √ s = 14 TeV.
We employ the recent results of ATLAS SUSY searches for neutralinos and charginos based on the 2l + E T and 2τ + E T channels to constrain the Z exotic decays in G2HDM. The resulting bounds should be regarded as estimated constraints, as the signal regions and efficiency may have some differences between MSSM and G2HDM. For the left panel of Fig. 2 (2l + E T channel), we use the signal region SR2l-A, which refers to a set of event selections listed in Table 1 of Ref. [57]. It gives rise to the constraint σ 95 obs ≤ 1.89 fb at 13 TeV. Assuming that the factor of signal efficiency is of O(1) at 14 TeV, 6 we can infer limits on g H at the 14 TeV LHC in the following way. For instance, to satisfy the SR2l-A bound σ < 1.89 fb, along the black, blue and red contours in the left panel of Fig. 2, the corresponding g H is required to be smaller than 0.137, 0.194, and 0.435, respectively.
Likewise, for the 2τ + E T channel on the right panel of Fig. 2, we utilize the SRC1C1 signal region in Ref. [58]. It yields σ 95 obs ≤ 0.33 fb at 13 TeV, which demands g H to be less than 0.057, 0.081, and 0.182 for the black, blue, and red contours, respectively.
Two comments are in order here, regarding the discrepancies between collider searches discussed here for MSSM and G2HDM.

The signals of 2l + E T and 2τ + E T in Spectrum
A mainly come from Drell-Yan processes for both G2HDM and MSSM. The major difference between the two models is the distribution of the invariant mass of the final charged leptons. In G2HDM, the invariant mass distributions have a cut-off at m Z , i.e., m l + l − , m τ + τ − < m Z , while m l + l − and m τ + τ − are more evenly distributed in MSSM. This is because the underlying Drell-Yan processes are mostly mediated by the on-shell Z in G2HDM, but by the off-shell SM γ and Z for MSSM, as indicated in Eqs. (9) and (10). 2. If the on-shell Z is highly boosted and the mass splitting between Z and the new leptons is large, one will have two collinear outgoing new leptons which result in two collinear SM leptons. In contrast, MSSM will not exhibit such a collinear behavior due to lack of Z , and so one can distinguish G2HDM from MSSM via the event topology of dilepton plus missing transverse energy signals. For the new quarks, they can always be pair produced dominantly by strong processes, like qq, gg → Q H Q H via s-channel gluon exchange or t-channel heavy exotic quark exchange. The cross sections for the strong processes We display the corresponding cross section contours, similar to Fig. 2 in Spectrum A , on the m Z -m j H , m Zm b H , and m Z -m t H planes, respectively, in the right column of Fig. 3. The previous bounds on the new quark masses translated from the SUSY searches are also shown by the dashed lines. Because of the stringent mass limits which push the new quark mass scale beyond TeV, the production cross sections of the new quarks via the Z exotic decays are kinematically suppressed and hence constrained to be small: σ/g 2 H 0.01 pb as can be seen from the right column of Fig. 3.

Future W searches
In the event that Z is discovered at HL-LHC via dijet or dilepton searches, one can ask whether it comes from an additional SU (2) gauge symmetry or simply from an extra U (1) . In this section, we discuss how to look for the electrically neutral W ( p,m) whose existence will help to pin down the SU (2) H as a potential underlying symmetry in nature. In the rest of our analysis, we will switch back to Spectrum A and Spectrum B.

W in different SU (2) models
Before embarking on our detailed analysis, we should point out that there exist, of course, many other wellmotivated models with neutral W gauge bosons in addi-tion to Z , such as THM [23,24], 3-3-1 models [62,63], etc. Certainly, any non-abelian gauge group commutes with the SM SU (2) L ×U (1) Y gauge groups, naturally accommodates neutral W . Here we will not manage to compare the W ( p,m) of G2HDM with all models featuring neutral W . Instead, we focus on collider signals of W ( p,m) in our model and contrast them with LHT with T-parity-a representative composite Higgs model which has neutral Z and charged W [20][21][22].
In the absence of the symmetries in the models, we put a slash in the cells. P T is the T-parity in LHT Both Z and W in LHT can only be produced in pairs due to T-parity and have the same signals, just like W ( p,m) in G2HDM. In other words, LHT is chosen as an illustrative example to underscore differences in the context of collider searches. Since W ( p,m) in G2HDM is always heavier than Z , it might not be easy to produce a pair of W ( p,m) (or even Z ) at the LHC, we will focus on the future 100 TeV proton-proton collider. To identify SU (2) H unambiguously, the discovery of new heavy fermions Q H and L H as well as scalars like H and H will also be necessary on top of W ( p,m) and Z . We note that LRSM [18,19] has the right-handed charged W ± R . It can be singly produced and directly probed by dijet resonance [64,65] or same-sign dilepton plus two jets (l ± l ± j j) searches [66], depending on the right-handed neutrino mass, whereas W ( p,m) in G2HDM must be pair produced. Due to quite different properties between W ( p,m) and W ± R and the stringent bound on Z R : m Z R > 3.2 TeV [35], LRSM will not be considered here.
The LHT model discussed here is based on the coset manifold SU (5)/SO (5), which can be realized as a nonlinear sigma model [20][21][22]. Two different SU (2) × U (1) subgroups of SU (5) are gauged and are broken down to the SM electroweak gauge group SU (2) L ×U (1) Y at a scale f T , which is higher than but not too far away from the electroweak scale so as to provide a possible solution to the fine-tuning problem. In the LHT model, all particles are divided into two classes based on the T-parity (denoted by P T hereafter), which corresponds to the symmetry under the exchange of the two SU (2) × U (1) subgroups. As a result, combinations of different fields of the two subgroups can be formed as having eigenvalue + 1 or − 1 of P T . The lightest T-odd particle is A H , a spin 1 particle, which is ensured to be stable and hence can be a DM candidate. Novel collider signatures like monojet and dijet plus missing transverse energy of A H in LHT was studied in [67]. All of the exotic particles have masses proportional to f T , since the masses are induced from the collective symmetry breaking at the scale f T . Furthermore, the exotic fermions couple to T-odd combinations of the gauge bosons of (SU (2) × U (1)) 2 . Three of the combinations comprise a non-abelian group SU (2) T , which is broken at the scale f T . In the end, the exotic T-odd SU (2) T gauge bosons should couple to one T-even and one T-odd par-ticle so as to conserve T-parity. We will use Z H and W ± H to denote the SU (2) T gauge bosons in the LHT, as opposed to Z and W ( p,m) in G2HDM. The quantum numbers of additional SU (2) L singlet fermions in the G2HDM and LHT are summarized in Table  We will focus on leptonic decay channels for these gauge boson pairs because of the low QCD background as in the Z resonance searches. The final states of two and four leptons plus the missing transverse energy, 2l + E T and 4l + E T , respectively, will be investigated. Take the W ( p,m) pair in G2HDM as an example. The 2l + E T channel comes from the prompt decay of one of the two gauge bosons into one heavy charged lepton plus one SM lepton, while the other boson into one heavy neutrino plus one light neutrino. Each of the two resulting heavy fermions then decays into the DM particle H 0 2 plus the corresponding light SM fermion through the Yukawa couplings from (2). The neutrinos and DM particles in the final state will escape from the detector and become manifest as E T . In Table 5, we summarize the decay chains of the gauge boson pairs into 2l + E T and 4l + E T . In the last column labeled by signal, the particles inside the curly brackets become manifest as E T .

The quantitative study: cross sections
In Fig. 4, we show the cross sections as functions of the corresponding gauge boson mass for the channels 2l + E T and 4l + E T (l = e, μ) with √ s = 100 TeV in both G2HDM (left panel) and LHT (right panel). The solid (dashed) line corresponds to Spectrum A (Spectrum B), while the red (blue) line refers to the channel 2l + E T (4l + E T ). Here we apply the same Spectrum A and Spectrum B for the gauge bosons, heavy exotic fermions, and DM in G2HDM to the corresponding particles in LHT.  The SU (2) T gauge coupling is fixed to be the SM SU (2) L coupling, whereas g H in G2HDM is set to be the maximally allowed value by the Z dilepton searches, presented in Fig. 1 For G2HDM, the upper bound from the Z resonance searches on the gauge coupling g H from Fig. 1 is used. In other words, the region above the line in each case is excluded. That is the reason why the cross sections, which scale as g 4 H , increase when m W (and also m Z ) become larger, since the bound on g H becomes less stringent. By contrast, in LHT the gauge coupling is set equal to the SM electroweak coupling and thus the cross sections decrease as the gauge boson mass increases.
We should point out that in G2HDM the process pp → Z → l H l H is actually the dominant contribution to the final state 2l + E T . If Z is discovered in the resonance searches, as assumed here, one should be able to infer the precise values of g H and m Z . Therefore, the dominant Z contribution can be subtracted from the data so that one can study the contributions from pp → W p W m alone. Alternatively, as we shall see later one can also resort to the 4l + E T final state for W ( p,m) searches. It is a relatively clean channel and is free from pp → Z pollution. The corresponding cross section is only three times smaller than that of the 2l + E T channel. Finally, the Z pair production pp → Z Z will also produce the 2l + E T and 4l + E T signals. The contributions, however, are at least two orders of magnitude smaller than those from pp → W p W m , and therefore will be neglected in our analysis.
In LHT, the cross section of pp → W + H W − H is about one order of magnitude larger than pp → Z H Z H at √ s = 100 TeV. That is because W + H W − H is produced dominantly by the s-channel γ and Z exchange, which is larger than the main contribution from the t-channel q H exchange to the Z H Z H production. Consequently, the cross section for 2l + E T , mostly from the W + H W − H channel, is almost one order of magnitude larger than that of 4l + E T , which arises only from Z H Z H channel. As seen from Fig. 4, however, the cross sections for both 2l + E T and 4l + E T are of the same order in G2HDM. On the other hand, the cross sections for both of these channels in LHT are roughly 1-2 orders of magnitude larger than those in G2HDM, depending on the gauge boson mass. Thus, one can in general distinguish the two models just by measuring the total cross sections of the two channels.

The qualitative study: the kinematical distributions
In this section, we will further investigate the difference between G2HDM and LHT gauge boson decays in terms of three different normalized kinematical distributions. In principle, one should be able to distinguish the electrically charged W ± H from the neutral W ( p,m) by the total charge of the corresponding decay products once they are produced singly. Both W ( p,m) and W ± H , however, have to be pair produced because of the SU (2) H and SU (2) T symmetry, respectively, which lead to the same total charge of the final states. As a consequence, the kinematical distributions of W ( p,m) and W ± H decays are not only interesting but also important to study for further information, even though the production cross sections, as shown in previous Section, are in general much larger in LHT than in G2HDM.
Let X denotes W ( p,m) , W ± H , or Z H , notwithstanding the same symbol has been used for the U (1) X gauge boson which has been assumed to be very heavy and decoupled. In four benchmark points m X = 0.5, 1.5, 3.0, and 4.0 TeV, we will show the normalized kinematical distributions of the spatial separation R e + e − (Fig. 5) and invariant mass M e + e − (Fig. 6) of the electron pair in the 2l + E T channel, and the invariant mass M e + e − μ + μ − (Fig. 7) of four leptons in 4l + E T channel. 7 The muon has the same distributions of R μ + μ − and M μ + μ − as the electron and will not be discussed separately. For comparison, the benchmark point m X = 0.5 TeV is also included because its distribution shape is clearly distinguishable from the SM background. The corresponding coupling g H for m X = 0.5 TeV shall be appropriately small to avoid the current LHC limits as shown in Fig 1. The normalized kinematical distributions, nonetheless, do not depend on the values of g H that we shall keep in mind here and hereafter. In addition, the distributions do not change significantly between the two spectra and we simply choose Spectrum A. The leading order irreducible SM background for each kinematical distribution is also presented for comparison. Further discussions of these three distributions are as follows. First, the spatial separation R e + e − the distribution for W ( p,m) is distinct from that of W ± H . As shown in Table 5, the final lepton pair l + l − is coming from different decay patterns of the gauge bosons W p W m , W + H W − H , or Z H Z H in the two models. In G2HDM, both l + and l − are coming from either W p or W m , while in LHT there are two possible routes-(1) l + from W + H and l − from W − H , and (2) both l + and l − from either one of the Z H in the Z H Z H pair. Overall, since W + H W − H pair production has a much larger cross section than Z H Z H , for the two leptons in the final state, l + mainly comes from W + H and l − from W − H ; it leads to a larger R e + e − in LHT than in G2HDM. Note that the distinction between the two models also depends on to what measure the gauge boson X is boosted. If X is highly boosted (e.g. m X = 0.5 TeV in Fig. 5), the distinction between G2HDM and LHT becomes more visible. For the SM case, the main contributions result from the W + W − pair, each of which decays into a charged lepton and a neutrino, resulting in a large separation R e + e − similar to the LHT case as is clearly reflected in Fig. 5.
Second, for the same reason as in the R e + e − distribution, the invariant mass of the lepton pair M e + e − should be smaller than m W ( p,m) in G2HDM (or m Z H in LHT). From the position of the M e + e − cut-off, one can roughly infer the mass m W ( p,m) . On the other hand, the invariant mass spectrum is distributed more evenly in the X = W ± H case. For the SM, M e + e − is centered around m Z as well as the low-mass region due to W + W − pair, off-shell γ /Z and other non-resonance contributions. The SM contribution is, of course, independent of m X and is shown only in the top-left panel of Fig. 6 to highlight the difference from the new physics. Therefore, one can simply impose appropriate cuts on M e + e − to reduce the SM background and extract either G2HDM or LHT signals for different m X at √ s = 100 TeV. From Fig. 6, one can clearly distinguish among the G2HDM, LHT, and SM by the M e + e − distributions.
Last but not least, the cross section of the 4l + E T channel is about three times smaller than those of 2l + E T for W ( p,m) searches in G2HDM, as can be seen in the left panel of Fig. 4. Due to the facts that only W ( p,m) in G2HDM and Z H in LHT contribute to this channel and both of them are neutral, similar distributions of M e + e − μ + μ − in G2HDM and LHT, as exhibited in Fig. 7, are expected. However, their distributions are clearly distinguishable from those of the SM, which arise from on-shell Z decays and consequently peak toward low-energy regions. Thus, the final state of four leptons plus missing transverse energy can be used to detect physics beyond the SM. However, to distinguish G2HDM from LHT is not easy using the 4l + E T channel unless M X is 0.5 TeV, in which case a much smaller g H is anticipated.
From our studies of the two channels of two/four leptons plus missing transverse energy, one can conclude that in addition to the production cross sections, the kinematical distributions are also indispensable to discriminate the two models, G2HDM and LHT.

Summary and outlook
One of the most interesting features in G2HDM [14], where the two Higgs doublets, H 1 ) can only decay into a new heavy fermion and a SM fermion. The heavy fermion can decay into a SM fermion plus missing energy carried by the DM candidate H 0 2 , whose stability is protected by an emergent Z 2 symmetry.
In this work, we studied collider signals of Z and W ( p,m) which can help us to pinpoint G2HDM. We derived constraints on Z from the LHC 13 TeV data, followed by investigations of the future 14 TeV LH-LHC and 100 TeV protonproton collider searches for Z and W ( p,m) . The main difference between Z and W ( p,m) is that Z can be singly produced and decay into SM fermions, while W has to be pair produced in the light of the SU (2) H gauge symmetry. It leads to stringent bounds on the SU (2) H coupling g H and the mass m Z from the LHC direct searches based on the Drell-Yan type dilepton and dijet final states.
The updated LHC limit shown on the m Z -g H plane in Fig. 1 is roughly a factor of 2 improvement in the low-mass region m Z 2.25 TeV compared to the previous results [14] inferred from the LHC-8 dilepton data. In addition, for the whole region of m Z of interest, the dilepton constraints on g H now become more stringent than those deduced from LEP and the electroweak precision data. Even the dijet bounds, which suffer from the QCD background, are also stronger than the LEP bounds for m Z 2.5 TeV. Moreover, Z can also decay into the new fermions which subsequently decay into SM fermions plus the missing transverse energy. We presented the contours of the rescaled cross sections σ/g 2 H on the plane of the masses of Z and new heavy exotic fermions in Figs. 2 and 3 for Spectrum A and Spectrum B , respectively. For Spectrum A , we considered final states of two leptons (es and μs) and two τ s with missing transverse energy which originate from decays of the exotic heavy leptons. For Spectrum B , two jets, two b-quarks and two t-quarks with missing transverse energy coming from the exotic heavy quarks were considered. Using the LHC constraints on σ 95 obs from SUSY searches for the same final states, one can derive limits for the coupling g H from the contours of σ/g 2 H on the (m Z , m l H ) or m Z -m τ H planes (Fig. 2). While the constraints on g H obtained from exotic decays of Z into heavy exotic leptons are comparable with those from the aforementioned dilepton and dijet searches, it appears that there is no severe bound that can be derived on the masses of heavy leptons from the existing data of SUSY searches. On the other hand, Spectrum B , in which the pair production of the exotic quarks is kinematically allowed, has substantial regions being excluded (Fig. 3). The reason is that the dominant contributions to the exotic quark production arise from the QCD processes which are independent of the gauge coupling g H . Thus, stringent limits from LHC on SUSY squark searches can be directly applied to our case, requiring the new quarks to be heavier than 1 TeV or so.
If Z can be discovered in the future, the neutral W ( p,m) also need to be found in order to identify G2HDM as the underlying theory. Final states of two and four leptons with missing transverse energy from G2HDM and LHT had been studied in detail to underline different signatures between the neutral W ( p,m) in G2HDM and the gauge bosons, Z H and W ± H , in LHT. The total cross sections for 2l + E T and 4l + E T were computed in both models (Fig. 4). In the 2l + E T channel, the two final leptons come from the same W ( p,m) in G2HDM, whereas in LHT they come from different W ± H . Therefore, a smaller spatial separation between the final leptons is expected for W ( p,m) . Indeed as clearly seen in Fig. 5, the spatial separation of e + e − in LHT is completely overlapped with the SM one which has larger R e + e − . Furthermore, by the same reason, the invariant mass of the lepton pair are cut off at the mass of W ( p,m) as opposed to a flatter invariant mass distribution for LHT (Fig. 6). In addition, the 4l + E T channel is also investigated. The invariant mass distributions of four charged leptons behave quite differently between the SM and G2HDM (and LHT), whereas LHT and G2HDM exhibit similar distributions (Fig. 7).
We conclude our study by presenting the search strategies of distinguishing G2HDM from LHT in Table 6. For Z , which can be singly produced, the dilepton final states will be the best search channels. For pair-produced gauge bosons like W ( p,m) in G2HDM, and W ± H and Z H in LHT, apart from the total cross sections, detailed kinematical distributions, such as (i) the spatial separation between the SM lepton pair and (ii) the invariant mass distributions of two and four leptons in the final states, can help us to disentangle these two models.
High luminosity upgrade for the LHC and building a future 100 TeV hadron collider are matters of utmost importance for fully exploring and distinguishing new electroweak scale models, all of which are contrived to address theoretical issues that cannot be answered within the SM. the cubic couplings and dimensionless quartic couplings are necessarily real because every term in the mixed potential V mix (H, H , H ) in (A7) is Hermitian.
The resulting coefficients of the quadratic terms for H 1 and H 2 after spontaneous symmetry breaking of SU (2) H induced by 3 = −v are The gauge coupling g H is fixed to be either 0.08 or 0.04 in order to satisfy the Drell-Yan constraints, shown in Fig. 1 gauge bosons left over after spontaneous symmetry breaking. One of them is naturally identified as the photon and the other one could play the role of either the light dark photon γ D or heavy Z . To render the dark photon γ D or Z massive, one can either resort to the Stueckelberg mechanism or introduce yet another Higgs field X solely charged under U (1) X to break one of the remaining two unbroken generators.
In this work, we confine ourselves to the scenarios with implying y f 0.1. In Table 7, we present four different benchmark points: the second and third columns are for either Spectrum A or Spectrum B, while the last two columns are for either Spectrum A or Spectrum B . Recall that, on the one hand, Spectra A and A are with lighter new quarks, while Spectra B and B are with heavier new quarks. The required Yukawa couplings y f for new fermions in G2HDM are of order 10 −3 for Spectrum A (A ) and order 10 −2 for Spectrum B (B ) which are quite acceptable. On the other hand, the mass spectrum of the Higgses is the same for Spectrum A and Spectrum B and for Spectrum A and Spectrum B . Similarly we have the cases for the mass spectrum of new gauge bosons. We should point out that the mass spectra displayed here are composed of points in the parameter space that satisfy the two theoretical constraints-the scalar potential is bounded from below and all relevant 2 → 2 scattering amplitudes among the scalars are below the unitarity bound. In addition to the theoretical constraints, we also impose the experimental constraints from the 125 GeV Higgs, including its mass and signal strengths decaying into diphoton and τ + τ − . A detailed study of these constraints on the scalar sector of G2HDM is presented in a separate work [38]. One should be aware that we choose a parameter space slightly different from that presented in Ref. [38] where the new gauge sector is much heavier than the Higgs sector. In this analysis, we choose the new gauge bosons Z and W ( p,m) having masses between 1.5 and 3TeV such that they are accessible at the LHC. Nevertheless the parameter space in the two works is chosen to satisfy the same set of aforementioned constraints.