Hypercharged Naturalness

We present an exceptional twin-Higgs model with the minimal symmetry structure for an exact implementation of twin parity along with custodial symmetry. Twin particles are mirrors of the Standard Model yet they carry hypercharge, while the photon is identified with its twin. We thoroughly explore the phenomenological signatures of hypercharged naturalness: long-lived charged particles, a colorless twin top with electric charge $2/3$ that once pair-produced, bounds via twin-color interactions and can annihilate to dileptons or a Higgs plus a photon or a $Z$, and glueballs produced from Higgs decays and twin-quarkonium annihilation that either decay displaced, or are stable on collider scales and eventually decay to diphotons. Prospects for detection of these signatures are also discussed.


Introduction
In an era of active experimental research but no direct signs of physics beyond the Standard Model, especially of the dynamics behind the insensitivity of the scalar Higgs field to higher mass scales, solutions to the electroweak hierarchy problem with non-standard low-energy signatures are called-for.
Recent progress in this direction ranges from theories where the gauge quantum numbers of the BSM states cutting off the UV contributions to the Higgs potential are unconventional, the prime realization of this idea being the so-called twin Higgs [1][2][3][4], to scenarios relying on a non-trivial cosmological history to explain the preference for a light Higgs sector [5,6].The former alternative, which we focus on in this work, is easy to understand and to motivate: the prevailing exploration of the high-energy territory is due to the LHC, a hadron machine where colored BSM particles can be produced with large cross sections.This collider has already shown its constraining power on such type of new physics, e.g.colored top partners in composite-Higgs models (see [7][8][9] for reviews), whose mass should now exceed a TeV [10,11], and similar bounds apply to colored supersymmetric particles.While this might indicate that the characteristic scale of the new dynamics is higher than expected from naturalness considerations, there is also the possibility that nature has just chosen none of the simpler low-energy manifestations of the TeV solutions to the hierarchy problem, specifically those with light colored states.
Indeed, in composite-twin-Higgs models, the Higgs arises as a Nambu-Goldstone boson (NGB) from the spontaneous breaking of a global symmetry at a scale f ≈ TeV, as in standard theories of Higgs compositeness, yet crucially the leading radiative contributions to the Higgs potential, from loops of SM fields, are cut off by uncolored states, their twins, owing to a Z 2 symmetry that relates them.In particular, the large explicit breaking of the Higgs shift-symmetry from the top Yukawa coupling no longer results in a squared Higgs mass proportional to (y t /4π) 2 m 2 * , with m * = g * f the typical mass of the colored composite resonances, because a color-neutral twin top introduces a complementary breaking of the shift-symmetry, enforced by twin parity and resulting in a (y t /g * ) 2 reduction.A twin top mass m t m t / √ ξ, where ξ = v 2 /f 2 , is predicted, while the resonances can be pushed to m * ∼ 5 TeV, beyond LHC reach and without compromising the fine-tuning of the EW scale.Consequently, the collider phenomenology of twin-Higgs models is substantially different from that of a standard composite-NGB Higgs, see e.g.[12][13][14][15][16].
In this paper, we present a new realization of the twin-Higgs paradigm, in fact the last missing construction with minimal NGB content and a bona-fide twin-Higgs mechanism, i.e. with a twin parity that could be enforced exactly.The model is based on a global SO (7) symmetry broken to the exceptional group G 2 , and its main phenomenological characteristic is that the twin fields carry hypercharge, and in fact they have the same electric charge as the SM states.This feature gives rise to new and exciting signals at colliders, which include long-lived charged twin scalars, twin quarks behaving as microscopic quirks, that after pair production form bound states that annihilate into SM final states, dileptons and γh, Zh, or twin gluons, and some twin glueballs which, while likely stable on collider scales, eventually decay almost exclusively to diphotons.Albeit some of these signatures have been already experimentally explored, here they are linked to the naturalness of the EW scale, a strong motivation that warrants further exploration at the LHC, in particular in its high-luminosity phase: since production cross sections are of electromagnetic size, large data samples mean excellent prospects to probe this scenario.
The same structural reason that requires the SM particles and their twins to have identical electric charges, also explains the absence of a twin photon.Because the global symmetry of the composite sector is a rank three SO (7), there is no place to embed both the SM and twin hypercharge U (1), after the SM weak SU (2) L and its twin SU (2) L have been accommodated; by construction the photon and its twin must be identified.The twin-Higgs model presented here is therefore smaller, in terms of its global symmetries, than previous constructions based on SO(8)/SO (7) [17][18][19], while retaining custodial symmetry.There is however a tradeoff: if the SU (2) L is gauged by twin weak gauge bosons, a Z coupled to the SM fermions is predicted, for which there exist strong experimental bounds from the LHC, rendering the model less plausible.This issue is resolved in the same way as in standard twin-Higgs constructions [13,19], where the twin photon is avoided by decoupling the twin gauge bosons at energies above the confinement scale.Such a departure from a complete, or fraternal [13], implementation of twin parity has the advantage of introducing the required source of explicit Z 2 symmetry to misalign the vacuum with respect to the twin symmetric one, i.e. ξ 1/2, as Higgs couplings measurements and EW precision constraints require.This the minimal scenario is therefore preferred from both phenomenological and theoretical considerations; yet we find the study of the fraternal version of the model illustrative at the very least.
Finally, let us briefly comment on how our exceptional twin Higgs compares with previous models of neutral naturalness.Because the twin quarks are uncolored but carry hypercharge as well as twin color, their phenomenological signatures resemble those of folded supersymmetry [20] or the quirky little Higgs [21].However, the scenario presented in this paper is not supersymmetric, and it enjoys a custodial symmetry, in contrast to the quirky little Higgs.Furthermore, because the twin sector is hypercharged but SU (2) L -neutral, its phenomenology shows a novel predilection for photons and Z's.Particularly, the twin top, which is the minimal light degree of freedom for a successful implementation of the twin-Higgs mechanism, can be produced at colliders with γ, Z-mediated cross sections that, although small, could be eventually probed by the HL-LHC, for instance in γh, Zh final states.
The rest of the paper is organized as follows.In the next section we present the symmetry structure and NGBs of the composite sector giving rise to the exceptional twin Higgs.In Section 3 we show how the elementary gauge and fermion fields are coupled to the strong dynamics, both for the SM particles and their twins.We study two different options regarding the latter: in the first (fraternal), all the twin global symmetries are gauged, thus it includes, beyond a twin top, a twin bottom and twin leptons, while in the second (minimal), no twin gauge bosons are considered and the twin NGBs remain uneaten.The scalar potential, for the Higgs and for the extra scalars (in the minimal realization) is analyzed in Section 4. Lastly, the exceptional phenomenology of the model is presented in Section 5, where both indirect and directs constraints are discussed, as well as the future prospects for discovery.

Higgs sector
The twin-Higgs mechanism relies on the Higgs arising as a NGB from the spontaneous breaking of a global symmetry.In this work we assume that such a breaking is driven by a strongly interacting sector that confines at a scale m * close to the TeV.The characteristic mass and coupling of the composite resonances is set by m * and g * respectively, related by the symmetry-breaking order parameter, f , as m * = g * f .
The novel features of our scenario stem from a different global symmetry breaking pattern, SO(7) → G 2 , with respect to previous twin-Higgs models.This breaking has the peculiarity of being the smallest that gives rise to seven NGBs while leaving an unbroken custodial symmetry, SU (2) L × SU (2) R ⊂ G 2 . 1 The twin-Higgs mechanism can still be operative in this coset, with the Higgs complex doublet H and its twin H embedded in the spinorial 8 representation of SO (7).Interestingly, SO (7) has rank three and it contains as a subgroup the product SU (2) L × SU (2) L × SU (2) R , out of which the unbroken SU (2) R is the diagonal combination of the SU (2) L and SU (2) R .As it will become clear in the following, a discrete Z 2 symmetry, a.k.a.twin parity, requires the SM particles and their twins to share the same quantum numbers under SU (2) R , which in turn implies they have identical electric charges, with the SM photon being its own twin.
We parametrize the SO(7)/G 2 coset with the Σ vector associated with the vacuum expectation value of a spinor of SO (7), that is Σ ∼ 8. 2 Given the decomposition 8 = (2, 1, 2) + (1, 2, 2) under SU (2) L × SU (2) L × SU (2) R , the Higgs and its twin are identified as where h ≡ π 4 .It then follows that under the custodial SU (2) L × SU (2) R symmetry, the Higgs has the proper quantum numbers, i.e.H ∼ (2, 2), while the twin Higgs decomposes as H ∼ (1, 1) ⊕ (1, 3), which correspond to the singlet radial component, σ, and a triplet of SU (2) R .Their kinetic terms take the same form as in other spherical cosets, Twin parity in the strong sector is realized as a discrete SO(7) transformation of the form whose action on Σ enforces H ↔ H.This parity also interchanges the SU (2) L and SU (2) L , while leaving the SU (2) R intact.This explains why the twin NGBs, ω± ≡ f (π 6 ± iπ 5 )/ √ 2 and ω0 ≡ f π 7 have the same electric charge as those eventually eaten by the W ± and the Z.
Let us finally recall that in composite-twin-Higgs models the strong sector is also symmetric under SU (3) C ×SU (3) C ×Z 2 transformations, with the Z 2 exchanging C ↔ C.Moreover, to reproduce the hypercharges of the SM fermions an extra U (1) X abelian symmetry is introduced, along with its counterpart U (1) X and with twin parity enforcing X ↔ X. 3

Gauge and fermion sectors
To specify the couplings of the external (i.e.weakly coupled in comparison to g * , elementary) SM gauge and fermion fields and their twins to the strong sector, we need to identify the operator content of the latter.Besides the scalar operator responsible for the spontaneous breaking of SO(7) to G 2 , there is a vector current associated to each global symmetry of the strong dynamics: J A ∼ 21 of SO (7), which under G 2 decomposes in unbroken J a ∼ 14 and broken J â ∼ 7 components, J X, X µ , and the color and twin color currents J C, C µ .The weak gauging of a given global symmetry then gives rise to the coupling of the elementary gauge fields with the corresponding composite currents, i.e. g i A i µ J µi .Besides, in order to implement the mechanism for fermion mass generation known as partial compositeness [24], we assume the strong sector contains fermionic operators Ψ in non-trivial representations of SO (7) The interaction of the elementary fermion fields ψ with the strong dynamics proceeds via linear couplings of the form y ψΨ.The actual set of fermionic operators depends on which twin fermions are present at energies around m * , that in turn depends on which of the global symmetries are gauged.We discuss the most interesting possibilities in the following.

Fraternal
A fully fledged analog of the standard twin-Higgs models, which are based on a global SO (8) symmetry and where twin partners are introduced for all the SM gauge bosons and for at least one complete generation of SM fermions [13], can also be constructed upon the symmetry group SO (7).In this case however, because of the smaller rank of SO(7) compared to SO (8), the gauge content of the twin sector is necessarily reduced.In particular, elementary massless vectors gauge SU (2) L × SU (2) L × U (1) Y , with where the three gauge bosons associated to the broken directions become massive, while the unbroken hypercharge group is identified with the diagonal combination of U (1 Only after the Higgs field develops a VEV, H ∼ v, does the SU (2) L ×U (1) Y EW symmetry gets broken to U (1) Q and the W and Z become massive.From such a pattern of symmetry breaking it follows that electric charge is given by x i 5 g T + w P n 8 A S J 7 j k w = < / l a t e x i t > h < l a t e x i t s h a 1 _ b a s e 6 4 = " A v q j 6 D R n e I U 3 5 9 F 5 c d 6 d j 0 V r w c l n j u E P n M 8 f z 0 e M 8 A = = < / l a t e x i t > W < l a t e x i t s h a 1 _ b a s e 6 4 = " y x U g N F p p S s e b k U q o 0 H , , ⇢, . . .

< l a t e x i t s h a 1 _ b a s e 6 4 = " H / I
y E e X q I p u U Q 3 V E U W P 6 B m 9 o j f n y X l x 3 p 2 P 2 e i S k 3 s O 0 J 9 y P n 8 A i b 6 W P Q = = < / l a t e x i t > ⌧ < l a t e x i t s h a 1 _ b a s e 6 4 = " V r Y G E U W w n j e F e M V j k b < l a t e x i t s h a 1 _ b a s e 6 4 = " P F A X q 8 y o y Q P thus twin parity implies that the SM particles and their twins have the same electric charge.
The masses of the gauge bosons follow from Eq.( 3) upon introducing covariant derivatives In the unitary gauge, where h = (v + ĥ)/f , with v ≈ 246 GeV and ĥ the Higgs boson.This equation clearly shows the upshot of twin parity, which recall exchanges SU (2) L and SU (2 As a result, we can consistently impose a discrete Z 2 symmetry acting on the elementary gauge fields as W α µ ↔ W α µ and trivially on B µ , which enforces g = g while leaving ĝ free.This constitutes a novel implementation of the twin-Higgs mechanism in the gauge sector, leading to important phenomenological departures from the SM, but no (massless) twin photon.To see this right away, let us work in the limit ξ ≡ v 2 /f 2 → 0, that is neglecting EW symmetrybreaking (EWSB) effects.The SM hypercharge gauge field is given by B = c θ B + s θ W 3 with tan θ = s θ/c θ = ĝ/g and gauge coupling g = gs θ.The orthogonal combination, and, by virtue of the mixing s θ with B, it couples to the SM fermions, which carry non-zero B-charges (the same as their twins), given by 1.As we will show in Section 5.3, such a Z contributes to the EW precision tests (specifically Table 1: EW gauge quantum numbers of the SM fields and their twins.In the minimal scenario, with no twins for the SM W and Z, only SU (2) L × U (1) Y are gauged and the twin low-energy content reduces to tL,R and ω±,0 ⊂ H.
to the Y -parameter) and, most importantly, can be produced at the LHC with significant cross sections. 5The only other extra gauge boson, the W ± , gets a mass m W = gf √ 1 − ξ/2 and it has electric charge ±1, being the twin of the SM W ± .The contribution of the twin gauge bosons to the Higgs potential will be discussed in Section 4.
Regarding the elementary fermions, anomaly cancellation requires a twin partner for every SM fermion, at least for one complete generation: in practice this means for the third generation, since the twin top is the crucial player leading to a successful cancellation of the Higgs potential induced by the top.The EW gauge quantum numbers of the twin particles simply follow from twin parity and are given in Table 1.Of particular importance for the phenomenology of this scenario is the fact that the twin fermions have the same electric charge as their SM partners (since twin fermions are SU (2) L neutral).
Besides, as in previous twin-Higgs models we assume the twin quarks are not colored but instead carry twin color, i.e. they are fundamentals 3 of SU (3) C , which is gauged by the twin gluons.Twin parity exchanges them with the SM gluons, from where it follows g 3 = g3 , up to terms that explicit break the twin parity. 6he SM fermions and their twins get mass via their interactions with the strong sector, that we assume to be of partial compositeness type, at least for the top sector.The top quark fields t R , q L and their twins tR , qL then linearly couple, with strengths y L,R and ỹL,R , to composite operators.For a proper implementation of the twin-Higgs mechanism, the latter should transform under SO (7) ) .We recall that the spinorial representation of SO (7) At low energies the interactions of q L and qL can be simply obtained via the embeddings while those of t R and tR are trivial, being SO( 7) singlets.For instance, the top Yukawa couplings, generated at the scale m * with y t ∼ y L y R /g * and likewise for ỹt , are given by from where the twin-top mass follows, It is evident that a Z 2 symmetry acting on the elementary fields as q L , t R ↔ qL , tR leads to y t = ỹt .The rest of SM fermions and their twin partners can get masses in a similar manner, with twin parity enforcing the equality of their Yukawa couplings and thus m ψ = m ψ ξ/(1 − ξ).However, since as long as ỹψ y t the associated contribution to the Higgs potential will be negligible (see Section 4), the approximate equality y ψ = ỹψ need not be enforced.

Minimal
The construction presented in the previous section demonstrates that it is feasible to build a twin-Higgs model based on SO (7) as global symmetry and where twin parity could be exact.However, we advance here that the experimental constraints on the twin Z are quite severe, disfavouring such a twin-symmetric model in view of fine-tuning considerations.Besides, it is well known that the twin parity cannot be an exact symmetry at m * if the Higgs VEV is to be misaligned from the Z 2 -symmetric vacuum H 2 + H 2 = f 2 .For these reasons we consider in this section a realization that is more minimal in terms of the (elementary) twin particle content.
We assume that the gauged SU (2) L symmetry has been spontaneously broken at a scale Λ / Z 2 significantly above the strong sector confinement scale, such that at m * only SU (2) L × U (1) Y is gauged, in a similar fashion as the recent twin-Higgs constructions of [27,28].This has several important consequences: the formerly eaten twin NGBs ω±,0 are now physical and become massive (much in the same way as the QCD pions because of EWSB), and the twin fermionic content can be reduced to the essential, that is a (vector-like) twin top.
R n e I U 3 5 9 F 5 c d 6 d j 0 V r w c l n j u E P n M 8 f z 0 e M 8 A = = < / l a t e x i t > W < l a t e x i t s h a 1 _ b a s e 6 4 = " y x U g N F p p S s e b k U q o 0 H b < l a t e x i t s h a 1 _ b a s e 6 4 = " P F A X q 8 y o y Q P z X q y X q x 3 6 2 P e u m I V M 4 f o D 6 z P H y c 4 k s E = < / l a t e x i t > , , ⇢, . . .

< l a t e x i t s h a 1 _ b a s e 6 4 = " H / I
x i 5 g T + w P n 8 A S J 7 j k w = < / l a t e x i t > The kinetic terms for the twin Goldstones, along with derivative Higgs and self-interactions, follow from Eq. (3), where D µ ω+ = ∂ µ ω+ − ig B µ ω+ (ω + has electric charge +1, while ω0 is neutral).It is clear that in this scenario there can be no twin cancellation of the gauge radiative correction to the Higgs potential.However this can be kept small, i.e. without compromise for fine-tuning, as long as EW composite resonances are sufficiently light (see Section 4).
The top and its twin couple to the strong sector as in Eq. (10), where however we can consider the twin bottom is no longer present in the low-energy spectrum but decoupled above m * ,7 or that, even if it is light, its couplings differ from those of t.As a result the Yukawa couplings are now given by Given the explicit breaking of the twin parity from the absence of the twin W 's and the twin bottom (as well as the twins of the light flavors), it is no longer the case that y L,R = ỹL,R (thus neither that y t = ỹt ) at m * .As we will show in Section 4, this indicates there is an upper bound on the scales (like Λ / Z 2 ) where the sources of Z 2 breaking originate, such that the twin cancellation of the top radiative correction to the Higgs potential is still effective.

Scalar potentials
In this section we discuss the radiative generation of the Higgs potential and examine the conditions as well as the amount of fine-tuning required to achieve EWSB.We also compute the potential of the extra pseudo-NGBs in the case that no twin symmetries are gauged.
Fraternal: Let us discuss first the gauge contributions to the Higgs potential in the fraternal model, where the SU (2) L × SU (2) L × U (1) R symmetries are gauged.The leading order (LO) terms in the gauge couplings arise at O(g 2 i ) from 1-loop diagrams with either the W 's, the twin W 's or the B, where we estimated the coefficient C g based on naive dimensional analysis (NDA), with an O(1) uncertainty encoded in the parameter c g .This radiative correction is quadratically sensitive to the compositeness mass scale, which we parametrized by m ρ ∼ g ρ f /2 m * .In the twin-symmetric limit g = g, Eq. ( 14) does not depend on the Higgs, which is the primary feature of twin-Higgs models.Note also that in our realization the twin cancellation of the O(ĝ 2 ) correction is automatic, since B is its own twin or, in other words, H and H have identical U (1) R charges.In consequence, the usual hypercharge contribution to the Higgs potential is absent.
In contrast with most composite-NGB-Higgs models, ( 14) is in fact the most relevant correction to the Higgs potential in our fraternal SO (7) twin-Higgs model.This is because of the experimental constraints on the Z , the mass eigenstate of the W 3 -B system up to O(ξ) terms, see Section 5.3 for the details.Since the couplings of the Z to SM fields scale with s θ = g /g and its mass Eq. ( 9) is well approximated by gf /2, the bounds can only be satisfied either if the twin gauge coupling g is substantially above g, thus breaking twin parity and spoiling the cancellation of Eq. ( 14), or if f is pushed to several TeVs.As we show in Section 5.3, either case implies significant fine-tuning of the Higgs potential.
Minimal: In view of these considerations, we now turn to the scalar potential in the minimal model, where only SU (2) L ×U (1) Y are gauged.The LO gauge contribution formally reads as the l.h.s. of Eq. ( 14), but now evaluates to where the last term is from a B loop, since the twin ω± carries hypercharge.The contribution to the Higgs potential is similar to the LO gauge correction in standard composite-NGB-Higgs models, where in fact it is usually considered subleading and relatively unimportant compared to the corrections arising from the top (at least in partial compositeness).In contrast, Eq. ( 15) is certainly important in the present twin-Higgs model, since it can introduce the required amount of Z 2 -breaking to accomplish v/f 1, as explained below.EW gauge corrections at O(g 4 i ) are generically subleading and we neglect them in the following.The other important contributions to the scalar potential arise from the explicit SO(7) symmetry-breaking terms sourced by the top and its twin, of which two different types are relevant: those generated at m * , which we denote as UV (one such type, for gauge fields, is Eq. ( 15)), and those generated from IR loops, independent of the details of the strong dynamics.Let us discuss them in turn.
The LO corrections from 1-loop UV diagrams with q L and tL arise at O(y 2 L ) and O(ỹ 2 L ) and are given by8 where we used NDA to estimate the size of C y = C y , with the equality between the coefficients following from twin parity in the strong sector.Note that while the overall sign of this contribution (i.e. the sign of c y = O(1)) cannot be predicted without explicit information on the strong dynamics, a prediction is obtained instead for the relative sign of the mass terms of the Higgs and of the charged scalar.The latter arises because no bL loop has been included, since e.g. it acquired a vector-like mass (i.e.independent of f ) m b > m * [29] or because its couplings to the strong sector depart from those of tL , these two states no longer related by a gauged SU (2) L symmetry.In other words, if at m * a twin bottom were present (i.e. in the fraternal model or if m b m * ) and its coupling to the strong sector were still fixed by ỹL , then the |ω + | 2 term would be absent.In contrast, a potential for the neutral scalar ω0 is automatically absent because neither y L nor ỹL break the corresponding U (1) L− R shift symmetry.Finally, an exact twin parity would enforce ∆y 2 L ≡ y 2 L − ỹ2 L = 0, thus exactly cancelling the |H| 2 term.This cancellation is the raison d'être of twin-Higgs models, and in the present realization does indeed take place at LO in the elementary (weak) couplings.However, since the Z 2 symmetry is not exact, ∆y 2 L will generically be non-vanishing at the relevant scale, m * , due to renormalization group evolution (RGE) from the Z 2 -breaking couplings, i.e. at NLO.In our case the latter are the EW gauge couplings g, g , which contribute as The coefficients A and A parametrize our ignorance on the strong dynamics at scales above m * and are a priori O(1) in size.Another potential source of Z 2 breaking percolating to the top sector at one loop depends on the vector-like mass of the twin b.If m b m * , then loops of q L = (t L b L ) cannot be matched by those of tL , inducing (∆y 2 L ) b ∼ (y 4 L /16π 2 ) log m b /m * .We should note that this correction is model dependent, in particular there is no experimental reason to completely decouple the twin bottom.Besides, if twin colored fermions are decoupled, a differential running of the SU (3) C and SU (3) C gauge couplings is induced, which eventually adds to ∆y 2 L (formally a 2-loop effect, i.e.NNLO, but could easily be numerically important).
Other potentially relevant UV radiative corrections to the potential arise from 1-loop diagrams with 4 insertions of either y L or ỹL , where D y = D y from twin parity.Even if Z 2 symmetric, (18) contains both a Higgs mass and quartic terms.The latter could be particularly important in order to reproduce the physical Higgs mass, depending on the size of the IR contributions to the Higgs potential, which we discuss in the following.Below m * , loops of the top and its twin give further corrections to the scalar potential via the Yukawa couplings in (13).The leading logarithmic (LL) term reads as in standard twin-Higgs models, where log t ≡ log(m 2 * /m 2 t ) and likewise for log t.This contribution to the Higgs potential is generically identified as the leading one in twin-Higgs models, being logarithmically enhanced in comparison to UV terms such as (18).In our scenario the same would be true for the charged twin potential if it was not for Eq.(18), which is quadratically sensitive to m * (there identified with m Ψ ).
Finally, we should note that NLL corrections to Eq. ( 19) from RGE due to the top/twintop Yukawas and color/twin-color interactions have been shown in [27,30,31] to be numerically important, in particular for extracting the physical Higgs mass.We expect a similar analysis could be performed in our scenario, leading to similar results, at least qualitatively if not quantitatively.In this work we simply bear in mind such corrections when presenting our O(1) estimates in the next section.

EWSB and Higgs mass
The set of contributions to the Higgs potential presented above can be simply parametrized as (focusing only on the relevant component h) [19], where and we have taken ỹt = y t and y L = ỹL in ( 19) and ( 18) respectively, which is a good approximation at the order we are working.The Z 2 -breaking term α depends on whether the model is fraternal or minimal, Minimal: where we have taken m ρ = g ρ f /2 and m Ψ = g Ψ f in (14, 15) and ( 16), respectively, and assumed A = A in Eq. ( 17).
The first point to note is that in our scenario the Higgs quartic is approximately the same as in standard twin-Higgs models.In particular the IR contribution to the physical Higgs mass provides a significant fraction of the observed value, where we have evaluated the top Yukawa at high scales, y t (1 TeV) ≈ 0.85, in order to roughly include NLL effects, and taken m * = 5 TeV and m 2 t m 2 t /ξ with ξ = v 2 /f 2 = 0.1.It is then clear that UV corrections, proportional to d y in (21), can easily and naturally provide the missing fraction of m h .
Therefore, to a good approximation the degree of fine-tuning required in our construction is determined by how unlikely it is to achieve a realistic EWSB, something that is directly controlled by the size of the h 2 terms in Eq. (20).As in most twin-Higgs models, there is a minimum amount of tuning (∆) as a direct consequence of twin parity: the Z 2 -symmetric term in Eq. ( 20) leads to ∆ min = 2ξ, while ξ is bounded from above from direct and indirect measurements of the Higgs couplings.Another way to see this is that twin parity implies that the minimum of the potential is at h 2 = ξ = 1/2, thus some other finely-tuned Z 2 -breaking contribution is needed to misalign the vacuum at ξ 1.This is in fact the reason why sources of explicit Z 2 -breaking are needed in twin-Higgs models, whose leading effect in the Higgs potential we have encoded in the term proportional to α (α > 0 in order to accomplish such a misalignment).Note that in the case that α is controlled by a single Z 2 -breaking source of the right size to reproduce a given value of ξ, then the associated tuning is simply given by ∆ min .However, in more complicated situations with several Z 2 -breaking sources, it is possible that a tuning between them is needed to achieve a given ξ, thus increasing the overall tuning.A better measure of fine-tuning, applicable to either case, is ∆ i = m 2 h /4α i f 2 .Fraternal: The fraternal scenario illustrates both of these possibilities.The twin gauge coupling g could be larger than g such that, provided c g < 0, the first term in (22) can be tuned to the Z 2 -symmetric piece (2β log a), how much tuned determined by ∆ min ; the hierarchy between f and v is however no longer determined by the bounds on the Higgs couplings, but by the bounds on the Z mass (Section 5.3), leading to ∆ min 1% for g ρ 4 (thus m ρ 6 TeV).Otherwise, f can be kept relatively low consistently with the Z bounds if there is another source of Z 2 -breaking, e.g. the second term in (22), and g is sufficiently large.In such a case the tuning is no longer tied to v/f , but instead is well approximated by m 2 h /4αf 2 , clearly worse than ∆ min for fixed m ρ .Minimal: The situation in the minimal model is certainly better, since that is, none of the Z 2 -breaking terms in (23)

Twin pseudo-NGBs
Let us focus first on the mass of the charged twin scalar present in the minimal model.This arises from both gauge and top/twin-top radiative corrections, Eq. ( 15) and Eqs.(16,18,19) respectively.The largest of these is the UV O(ỹ 2 L ) contribution in (16), which is quadratically sensitive to the cutoff, there parametrized by m Ψ .As long as the coefficient c y is negative, 9then the mass-squared is positive and of size The other corrections give rise to smaller contributions: the gauge contribution is suppressed by O(g 2 /y 2 L ) and yields ∆m 2 ω+ ∼ (100 GeV) 2 (m ρ /3 TeV) 2 , while the IR twin-top contribution (which is negative) is suppressed by O(m 2 t /m 2 Ψ ).After EWSB there also small corrections from the quartic coupling |H| 2 |ω + | 2 , suppressed by ξ.Other contributions to the ω+ mass could also be present, e.g. from loops of twin taus, if these are in the IR spectrum.
The attentive reader will have already noticed that none of the radiative corrections discussed above gave rise to a potential for the neutral twin scalar.This was to be expected, since neither gauging only U (1) Y nor considering a left-handed twin bottom with couplings different than those of the twin top (or just a decoupled b), breaks the U (1) L− R shift symmetry protecting the ω0 .However, an exact global U (1) is not at all guaranteed, on the contrary quite generically a source of explicit symmetry breaking will be present, lifting the ω0 .In fact, it is easy to imagine examples for such a source.One instance, perhaps not the best but certainly simple, is to consider a f -independent mass for the twin top m t tR tL , generated below Λ / Z 2 where no gauge symmetry forbids it, but naturally small m t m * since still protected by a chiral symmetry [29].Such a mass in fact breaks explicitly U (1) L− R , and gives rise to an extra radiative correction to the scalar potential where we used NDA to estimate the size of C m , ỹt ∼ ỹL ỹR /g * y t , and included a possible phase difference γ between ỹt and m t .Besides the tadpole term (which vanishes for sin γ = 0), ( 29) contains a ω0 mass of size (taking cos γ = 1 and c m > 0) The f -independent mass for the twin top is just one example of how the singlet twin could be lifted.Yet one should be aware that regardless of how ω0 gets a potential, it should be such that the Higgs potential is not significantly altered, not to raise the fine-tuning of the EW scale.In our example this would happen if m t 100 GeV, which explains why the reference value used in Eq. ( 29).Of course m t could be much smaller and the neutral twin much lighter accordingly.Let us note in this regards that ω0 does not linearly couple to SM fermions, but if kinematically allowed it opens a new Higgs-decay channel, which would force f 1.2 TeV to ensure consistency with LHC data.

Phenomenology
In this section we discuss the collider signals of the exceptional twin-Higgs model.While we focus mostly on its minimal incarnation, the main constraints on the fraternal scenario will also be discussed at the end of this section to understand why such a case is less attractive.
While sharing some of standard signatures of twin-Higgs models, such as the modification of the Higgs couplings and non-standard Higgs decays, these models present novel phenomenological features stemming from the fact that the twin particles carry hypercharge.We will first examine the indirect effects of the twin states on the well-measured properties of the SM particles.These are primarily the Y -parameter (which encodes new physics contributions to the hypercharge propagator), the running of the hypercharge gauge coupling g , and the Higgs decay rate to photons.These departures from the SM arise at the 1-loop level (with the exception of those mediated by the twin Z in the fraternal model).Other indirect effects, common in constructions where the Higgs is a (custodially protected) composite pNGB, such as a the S-parameter or deviations in other Higgs couplings, will not be discussed here further and we will merely recall when necessary their implications on the parameters of the strong sector, f and m * ∼ m ρ , m Ψ [31]. 10irect production of the twin particles at colliders gives rise to the most interesting, though in some cases challenging, signatures of the exceptional twin Higgs.The twin quarks, being hypercharged, can be pair-produced via Drell-Yan (DY) and, because of the existence of twin strong interactions, exhibit the typical characteristics of the so-called quirks [21,35,36].Such quirks are microscopic and lead to the formation of bound states that predominantly decay either to twin glueballs or back to the SM, while the glueballs in turn decay via an offshell Higgs or, whenever this channel is forbidden, to photons.Another attractive signature arises from the charged twin pNGB ω± or the twin leptons: they are pair-produced in DY and likely stable on detector scales, thus subject to current long-lived charged particle (LLCP) searches.Finally, we also present the stringent constraints on the twin Z .These are the main reason why the minimal model is preferred, since such bounds are absent and one is left with a light but almost inert ω0 , whose phenomenology is that of a twin axion-like particle that only couples to SM hypercharge, the study of which we defer to a future study. 11

Indirect effects
Since the twin particles carry hypercharge, they induce at the loop level a non-standard selfenergy for the hypercharge field. 12At leading order in a momentum expansion, such type of corrections are customarily parametrized by the dimension-6 operator −(Y /4m 2 W )(∂ ρ B µν ) 2 .Indeed, N ψ heavy fermions of mass m ψ and hypercharge q ψ generate a Y -parameter of size For e.g. the twin top, N t = 3, q t = 2 3 and m t m t / √ ξ, leading to a very small contribution, Y t ∼ 1×10 −5 for ξ = 1/4, in comparison to the per mille constraints from LEP [37] and more recently from the LHC [38,39].In fact, for any type of twin fermion we find that Y ψ is below LEP sensitivity for masses m ψ 100 GeV.For an even lighter twin fermion, encoding its effects in the Y -parameter is no longer adequate, nor it is for the twin top at the LHC [39].In these cases however, constraints could a priori still be placed by considering the contribution to the running of the hypercharge gauge coupling, encoded in the beta function coefficient ∆b Y in (30).However, the present accuracy on neutral DY processes is too low to see any such effect [40].The charged pNGB ω± also gives a small contribution to the Y -parameter at one loop, which reads as in Eq. ( 30) with m 2 ψ → 2m 2 ω+ , see Eq. ( 27), and ∆b Y = 1/3.More important are the 1-loop contributions of the twins to the Higgs decay to diphotons.In composite-Higgs models deviations of such a rate are generically O(ξ) and a consequence of the modified couplings of the Higgs to the top and the W .The latter are also present in twin-Higgs models, since they are intrinsic of the NGB nature of the Higgs [41].On top of these, because the twin particles are electrically charged in our scenario, a direct O(ξ) contribution is generated as well, specifically from a twin top loop.Since m t m h , such an effect can be parametrized by the dimension-6 operator (c γ g Note that in contrast to standard composite-Higgs models [41], this operator is not suppressed by y 2 t /g 2 * . 13The associated contribution to the hγγ vertex is (c γγ ) t/(c γγ ) t −ξ/(1 − ξ), where (c γγ ) t is the top loop contribution in the SM; a result that matches the expectation from twin parity once ξ = 1/2.The contribution from the twin top therefore, being opposite in sign to the top's, increases by O(ξ) the Higgs coupling to photons.Since the standard corrections lead to a reduction of hγγ, c γγ (c γγ ) SM √ 1 − ξ, the twin-top loop alleviates, rather than aggravates, the departure from the SM prediction, a fact that can become relevant as the precision on Higgs couplings measurements improves.Besides, since the twin top is uncolored, it does not affect the Higgs coupling to gluons, thus breaking the correlation between hγγ and hgg usually found in composite-Higgs models.

Direct production
All our production cross sections have been computed using the MSTW2008NNLO parton distribution functions [42].

Long-lived charged particles
Some of the hypercharged twin particles are potentially stable on detector scales.This is the case of the charged twin scalar in the minimal model, since its leading interactions with the SM respect an accidental Z 2 symmetry ω± → −ω ± , while we expect the decay channels ω → t + b, τ + ντ to be kinematically closed, the latter because the experimental bound on a collider-stable twin tau is stronger than for ω, as we show in the following (thus the decay τ → ω + ντ is open).Besides, we recall that minimality does not require neither the twin bottom nor the twin tau to be present in the IR, while in the fraternal model both of them are required for anomaly cancellation.
Searches at the LHC for such type of particles have been carried out at √ s = 13 TeV, the latest analysis from ATLAS with 36.1 fb −1 of data [43].In Fig. 3 we show the corresponding constraints on the production cross sections for a |Q| = 1 charged scalar (left panel), identified with a stau in the ATLAS analysis, and a |Q| = 1 DY charged fermion (right panel), identified with a chargino.These constraints readily apply to our twin ω+ and τ , respectively.The comparison of the predicted cross section with the experimental bound for the former shows that masses below 420 GeV are excluded.This should be compared with our expectation for the twin scalar mass m ω+ ∼ 550 GeV, Eq. (27).Even though the Figure 3: Upper 95% CL limits on the cross section for DY pair-production of |Q| = 1 heavy stable scalars (left) and fermions (right) from an ATLAS search in 36.1 fb −1 of data at √ s = 13 TeV [43] (red).The theoretical predictions for the twin pNGB ω± and the twin tau τ are also shown (blue).current bound does not quite reach such values, a factor of a few improvement in experimental sensitivity would start probing the relevant parameter space.For what regards the twin tau, the situation is certainly different, since twin parity would imply m τ m τ / √ ξ, while data indicates that only masses above 700 GeV are allowed.Of course, explicitly breaking twin parity in the lepton sector is a priori an option.However, the experimental bound is so stringent that the twin-tau Yukawa should be larger than the twin-top's to avoid it, thus a large contribution to the Higgs potential should be expected.This extra complication is yet another reason why the fraternal implementation of our scenario is disfavoured versus the minimal model, where the bounds on the charged scalar are less consequential and the τ could be lifted by giving it a vector-like mass.
While LLCPs could be regarded as a differential phenomenological feature of the exceptional twin Higgs, it is important to note that neither the stability of τ nor of ω+ are due to super-selection rules.Indeed, at low energies these states only carry electric charge and therefore decays such as τ → τ + γ, Z or ω → + ν, q + q, W + γ, Z are allowed.Let us discuss in some detail the decay of the twin tau (a similar discussion holds for ω+ ).The first point to note is that, since τ and τ have different X (and X) charges, the interaction mediating the decay cannot be generated by the strong dynamics alone, thus it should have its origin at some UV scale Λ.One instance is an interaction of the form τR σ µν H ˜ L B µν , which could be generated at low energies with a coefficient ĝ/m * Λ, leading to a τ decay rate Therefore, model-dependent UV considerations could render the twin tau (and/or ω+ ) shortlived, displaced, or long-lived.It would be interesting to study these types of signatures at the LHC, e.g.similar to those of an excited tau if the twin tau decays promptly.Finally, we restate that given the constraints in Fig. 3, the decay τ → ω + ντ would also proceed in the minimal model, and dominate over (32).

Hypercharged quirks
The lightest twin quark of the exceptional twin Higgs behaves as a quirk [35], a heavy stable particle that interacts via a new unbroken non-abelian gauge group and carries SM charges, in our case twin color and hypercharge, respectively.Once pair-produced, via DY in our scenario, quirks do not hadronize but instead form (meta-)stable strings.This is because any of the twin quarks are heavy in comparison with the scale where twin-color interactions become strong, i.e. m q ΛQCD : in such a case the breaking of the string by pair production takes an exponentially large time, t break ∼ (4π 3 /m q) exp(c m q/ ΛQCD ) 2 , with c an O(1) factor that depends on the precise definition of ΛQCD .The absence of twin quarks lighter than ΛQCD follows from several considerations.On the one hand, naturalness of the Higgs potential indicates that the SM and twin-color gauge couplings are not substantially different at m * , gs ≈ g s , thus ΛQCD will not be far apart from the scale Λ QCD ≈ 250 MeV where the SM color becomes strong.The difference arises mainly from Z 2 -breaking radiative effects that originate from the different colored content of SU (3) C and SU (3) C below m * , since we assumed a single generation of twin quarks (and m q = m q ).Such an effect leads to ΛQCD ∼ 2−10 GeV, depending on m * , on the exact values of gs and g s at m * , and on the twin quark masses [13].On the other hand, experimental constraints on the twin quarks force them to be considerably heavier than ΛQCD .In particular, the precise measurement at LEP1 of the Z decay width, ∆Γ Z /Γ Z ≈ 9 • 10 −4 [44], rules out twin quarks (with Y q = 2/3 or −1/3) for which the decay channel Z → q q is kinematically open, that is m q m Z /2 ≈ 45 GeV.This constraint is only relevant for the twin bottom, since it implies m b > m b / √ ξ, away from the Z 2 -symmetric relation, while the twin top is always heavier.In summary, we conclude that m q/ ΛQCD 4 and that t break is very large, much longer, as we show next, than the time it takes for the string to annihilate, which is the other possible fate of the string. 14 The typical annihilation time depends a priori on the annihilation rates of the different states of energy and angular momentum the string can be in.In practice a good proxy is to consider only the low-lying low angular momentum ( = 0) bound states, whose lifetime is proportional to the classical crossing time of the quirks (or length of the string), 14 Heavier twin quarks, up to m q ≈ 103.5 GeV, should have also been pair produced at LEP2 through an off-shell photon or Z.However, being the cross sections of O(g 2 ) and the quirks losing a significant fraction of their energy through electromagnetic radiation (see below), we expect LEP did not have the required sensitivity.Besides, even though the twin quarks modify at one loop SM processes like dilepton pair production (this being the same effect we encoded in the Y -parameter for twin masses beyond LEP reach, see Section 5.1), we checked the effect is below experimental uncertainties [45].We note that twin-QCD dynamics could also affect such processes, e.g.via Z-glueball mixing, however we expect the effects to be small; also, the decay Z → ggg is below current constraints on non-standard Z decays, BR(Z → X BSM ) 10 −4 .
i.e. our twin strings are microscopic, L ∼ 0.1 fm.This approach might naively seem at odds with the fact that most twin quarks will be produced relativistic, √ ŝ − 2m q ∼ m q, therefore in highly excited states.However excited bound states have annihilation rates suppressed by their large angular momentum [35].Therefore the low-lying states are eventually reached after the quirks radiate away most of their energy into (relatively) soft twin glueballs and photons [46].For the latter the typical radiation time can be estimated as t QED rad ∼ (3/8πq 2 q α)m 3 q / Λ4 QCD ≈ (2/q 2 q ) × 10 −22 sec for the same parameters as in (33) and q q the twin electric charge.This means that the de-excitation process is fast, in particular no displaced vertices in quirk pair-production and annihilation should be expected.Besides, twin glueball emission is also expected to contribute to the process, at least for energies where the kinematical suppression from the non-zero glueball mass is irrelevant. 15et us focus therefore on the lightest twin quarkonium bound states with = 0, which are a pseudoscalar η −+ and a vector Υ −− , both electrically neutral (see also [47,48]). 16Their mass is given to a good approximation by m η,Υ ≈ 2m q, while their decay rates can be found e.g. in [49][50][51].Of special significance is the fact that the decay of Υ −− to a pair of twin gluons is forbidden, thus enhancing its branching ratio to SM final states.Neglecting the masses of the daughter particles as well as m Z , where Q Υf V = q q(q f + q V Zf /c 2 θ W ) and Q Υf A = q q q A Zf /c 2 θ W with q i the electric charge of i and q V,A Zf the vector and axial Z-charges of the SM fermions, e.g. for f = e, q V Ze = 1 2 (− 1 2 + 2s 2 θ W ) and q A Ze = − 1 4 .The twin-QCD structure constant αs is to be evaluated at m Υ , 17 and ψ(0) is the radial wave function of the bound state at the origin.For a Υ −− made of twin bottoms (q q = − 1 3 ), the decay to three twin gluons is the largest for all twin quarks masses of interest, m q ∈ (0.1, 1) TeV, while for a bound state of twin tops (q q = 2 3 ), decays to γh and Zh dominate, owing to a large twin Yukawa coupling.Importantly, decays to SM fermions are always non-negligible, e.g.BR(Υ −− → e + e − ) ≈ 4−10%, where the lower end is found for a twin top, quite independently of its mass.In contrast, the η −+ decays almost exclusively to two twin gluons, while its other allowed decay channels, to a pair of SM neutral gauge Figure 4: (Left) Pair production cross sections of twin quarks with hypercharge Y q = 2 3 (solid) or Y q = − 1 3 (dashed).(Right) Upper 95% CL limit on Z section times branching ratio to dileptons from ATLAS searches in 36.1 fb −1 [52] (black solid) and in 139 fb −1 [53] (red solid) of data at √ s = 13 TeV, where the Z is identified with the bound state Υ −− ∼ q q of mass m Υ = 2m q, with theory predictions for Y q = 2 3 (blue solid) and Y q = − 1 3 (blue dashed).Also shown the projected bound after 3 ab −1 of integrated luminosity (see text for details).bosons (i.e.γγ, ZZ, Zγ), are suppressed by the small electromagnetic coupling constant, e.g.
For the relevant range of quirk masses BR(η −+ → γγ) ≈ q 4 q (1 − 4)%, which means that even for a twin-top bound state the branching ratio is always below the per cent.
The substantial branching ratio of the Υ −− to dileptons make its production at the LHC one of the most promising avenues for detection of the twin quarks.To assess the sensitivity of current LHC searches for dilepton resonances, we plot in Fig. 4 (left panel) the production cross section of pairs of twin quarks with hypercharge either as that of the twin top, Y q = 2 3 , or as the twin bottom, Y q = − 1 3 (recall q q = Y q).We keep the mass m q a free parameter, but we recall that in our twin-Higgs model m t m t / √ ξ while m b only needs to be large enough for Z → bb to be kinematically forbiden, but otherwise it is not bounded by naturalness considerations; we effectively vary m q by changing the corresponding vector-like mass m q (i.e.keeping the Yukawa coupling ỹq to its Z 2 -symmetric value).The fraction of events that go through the formation and decay of Υ −− bound states rather than of η +− is expected to depend on their total decay rates, Γ Υ,η , as r = 3Γ Υ /(3Γ Υ + Γ η ) [54], where the factor of three accounts for the number of degrees of freedom in the vector bound state versus in the pseudoscalar.Because of the large decay rate of the η +− to gg, we find r is always below 25% for a twin-top bound state and below 5% for the twin bottom.With these results we have computed the cross section for production and decay to e + e − (and µ + µ − ) of the Υ −− , which we show in Fig. 4 (right panel) for a bound state made of twin quarks with either top-like of bottom-like hypercharges.Comparing them with the latest ATLAS bounds on a Z decaying to dileptons [53], we find that the current 139 fb −1 of data at √ s = 13 TeV are not quite enough to probe the small cross sections associated with the hypercharged quirks.However, the order of magnitude increase in sensitivity needed to probe the relevant twin-top cross sections could potentially be achieved at the LHC with 3 ab −1 of integrated luminosity, as we also shown in Fig. 4.Such a expected a limit has been obtained from a naive luminosity rescaling, after contrasting the improvement of the expected bounds from 36 fb −1 [52] to 139 fb −1 [53] of data (note that the HL-LHC projections in [55] do not go below 2 TeV resonance masses).In addition, searches for γh resonances are also relevant in view of the large branching ratios of the twin-top Υ −− to this final state.Current analyses [56] set a constraint σ • BR 10 fb, roughly one order of magnitude above the cross section predicted in our scenario for m Υ ∼ 1 TeV, which makes this a very interesting signature for the high-luminosity phase of the LHC.
Similarly, the cross section for production and decay to diphotons of the pseudoscalar, σq q • (1 − r) • BR(η −+ → γγ) ≈ 1.6 × 10 −2 fb for m η = 1 TeV and Y q = 2 3 , is for the most part down by one/two orders of magnitude compared to current LHC bounds on diphoton resonances [57].Once again the interesting conclusion is that the HL-LHC should have enough sensitivity to probe a twin-top bound state decaying to γγ of mass as expected from twin parity, therefore reaching the relevant parameter space.We note that the situation is certainly different if twin quarks had larger hypercharges, for instance if Y q = 2 current diphoton searches would already exclude m q 1 TeV.

Twin glueballs
In our exceptional twin Higgs, twin glueballs are the lightest states of the twin sector.They are produced from Higgs decays, as in other scenarios of neutral naturalness [13,58], as well as from the decays of the twin-quark bound states, as discussed in the previous section. 18hile the former process dominates for heavy twin quarks, it is interesting that the decays of twin-bottom bound states lighter than approximately 300 GeV dominate the inclusive glueball production at √ s = 13 TeV: σ (h) g ≈ 67 fb (ξ/0.1) 2 while σ (Υ+η) g is very near the total quirk production cross section, shown on the left panel of Fig. 4, since bound-state annihilation to glueballs always dominates.
The mass, lifetime and decay modes of a given glueball are determined by its quantum numbers, see e.g.[59,60].Therefore, the actual signatures associated to glueball production depend on which and how many glueballs are produced from the aforementioned decays.These questions however cannot be reliably answered without a proper understanding of the non-perturbative twin-color dynamics (see however [61]).We will therefore focus our attention on the lightest glueball 0 ++ , whose dominant decays are mediated by the Higgs, via the dimension-6 operator |H| 2 G µν G µν , as well as on the lightest glueball for which such decays are forbidden, i.e. the 0 −+ , which decays almost exclusively to two photons.Both of Figure 5: (Left) Branching ratio of the lightest glueball 0 ++ to diphotons (orange) and lab decay length (blue) from production via either Higgs decays (dotted) or η −+ decays (solid), for ξ = 0.1.(Right) Lab decay length of the lightest parity-odd glueball 0 −+ (red), of mass m 0 −+ ≈ 1.5 m 0 , from production via either Higgs decays (dotted) or η −+ decays (solid).The decay lengths have been computed under the assumption that the parent particle produces two glueballs only, thus they should be understood as upper bounds.
these glueballs, of mass m 0 ++ ≡ m 0 ≈ 6.9 ΛQCD and m 0 −+ ≈ 1.5m 0 , are expected to represent a significant low-energy population of the energetic twin gluons originally produced, with the 0 ++ expected to dominate if the kinetic energy of the twin gluons is low, being the lightest.Other glueballs could also yield interesting and somewhat distinct phenomenology, however the 0 ++ and 0 −+ are good representatives of the particularities of our exceptional twin Higgs with respect to previous models.For instance, we checked that the 2 ++ glueball decays predominantly to diphotons and is stable on collider scales, as the 0 −+ , see Eq. (37).
Significant branching ratios to diphotons is in fact the identifying feature of the glueballs in our scenario.These arise via dimension-8 operators of the form c F µν F µν G ρσ G ρσ generated by a twin-quark loop, of size which is therefore enhanced at low twin-quark masses.For this reason in the following we will consider the twin bottom as a key player in glueball phenomenology, and comment when relevant on the differences that would arise if the b is decoupled and only the twin top contributes, e.g.Eq.( 36) would be enhanced due to the larger hypercharge of t but suppressed by its larger mass.In the left panel of Fig. 5 we show the branching fraction to diphotons of the lightest glueball in the (m 0 , m b) plane, in the range of glueball mass expected from our estimates of ΛQCD .As anticipated, only for light twin bottoms this channel dominates over the standard Higgs-mediated decay to a pair of SM bottoms.The decay length of the 0 ++ can be approximated, for sufficiently large m b and/or ξ by cτ 0 ++ ≈ (0.4 m)(15 GeV/m 0 ) 7 (0.1/ξ) 2 as in standard twin-Higgs constructions.Therefore, such glueballs are relativity long-lived, giving rise to displaced vertices in a significant fraction of the parameter space.This decay length, but in the lab frame, is also shown in the left panel of Fig. 5, where the boost factor is either from the decay of a Higgs (blue dotted) or the decay of the twin-bottom bound states (blue solid), to two, and only two, glueballs.More realistically, such decays will give rise to more glueballs, meaning these lengths should be considered as upper bounds.One of the new and exciting features of our glueballs is found in the right panel of Fig. 5, where we show the lab decay length of the 0 −+ , which almost exclusively decays to γγ and it is very long-lived, in fact stable on collider scales in a large portion of parameter space.This can be easily understood from the fact that its rest-frame decay length scales as which in the lab frame translates to up to 250 m.Since the dependence on either m 0 or m b is very strong, these glueballs could give rise either to displaced photons or to missing energy signatures in the detector.In the latter case, the glueballs would eventually decay to a pair of photons, which could potentially be detected by one of the proposed detectors dedicated to long-lived particle searches [62][63][64].Besides, notice that if the twin bottom is decoupled, cτ 0 −+ would be substantially larger, since the decay would only be mediated by the heavier twin top.Before discussing the experimental prospects for such collider-stable glueballs, let us understand the sensitivity of current and future LHC data to the 0 ++ glueballs.Since these are long-lived but always decay within the detector, they give rise to displaced vertices.This kind of signals have been studied in detail in several works, e.g.[66,[68][69][70][71], and we reinterpret here some of their results to account for the effects of our hypercharged twin quarks.In the left panel of Fig. 6 we show the constraints on m 0 and m b from several searches for displaced vertices.These in general focus on decays either in the inner detector, cτ γ 25 cm (γ the boost factor), therefore sensitive to heavier 0 ++ glueballs, or in the hadronic calorimeter and/or muon spectrometer, that is 2 m cτ γ 10 m and thus of relevance for lighter glueballs.The current most sensitive experimental analysis that is well suited to the features of our twin glueballs is from ATLAS at 8 TeV and with 20 fb −1 of data [65] (small orange region in Fig. 6).While this search barely reaches the interesting region, the theory projections of [66] for 13 TeV and 20 fb −1 (dotted red line) show that with current data one could already probe twin-bottom masses up to 150 GeV, at least for glueballs that decay within the inner detector. 19We should note however that these projections are somewhat uncertain, in particular to date there is no experimental analysis demonstrating the sensitivity to such heavy glueballs, while the region at small m 0 , which overlaps with Figure 6: (Left) Current and projected LHC constraints on the 0 ++ glueball mass m 0 and the twin-bottom mass m b from searches for displaced vertices, for ξ = 0.1.The bounds shown correspond to 0 ++ production from annihilation of twin-bottom bound states at 8 TeV in 20 fb −1 of ATLAS data [65] (orange), and projected searches from [66] at 13 TeV with 20 fb −1 (dotted red) and 3 ab −1 (dotted purple), as well as from the HL-LHC [67] (purple).Shown also the HL-LHC bounds on 0 ++ production from Higgs decays (green).(Right) Expected number of glueballs decaying to diphotons inside the MATHUSLA detector, from twin gluons produced at the 13 TeV LHC after 3 ab −1 of integrated luminosity, via decays of the Higgs and twin-bottom bound-states (black) or Higgs and twin-top bound-states with m t m t / √ ξ (blue), for ξ = 0.1.
the ATLAS 8 TeV analysis [65], has not been corroborated by a recent 13 TeV analysis that studies glueballs decaying in the hadronic or the outer edge of the electromagnetic calorimeters [72].We have checked that this experimental analysis, based on 10.8 fb −1 of data, does not yield any constraint on out parameter space.Besides, the analysis of [66] is not sensitive to glueballs originating from Higgs decays (for ξ 0.1), but only to those from the annihilation of twin-bottom bound states.In Fig. 6 we also show the region covered by a naive (by 1/ √ L) rescaling of the projections of [66] to L = 3 ab −1 (dotted purple).This simple extrapolation shows that the high-luminosity phase of the LHC could considerably extend the coverage up to m b ∼ 300 GeV.Moreover, we have also recasted the HL-LHC projections from [67], which shows sensitivity to light 0 ++ glueballs produced from twinbottomonium annihilation up to m b ∼ 400 GeV (purple region), as well as glueballs from Higgs decays (green).The latter region is independent of m b since the twin bottom does not contribute significantly to the Higgs branching ratio into twin gluons, being instead solely determined by ξ.Indeed, for ξ 0.01 the green region in Fig. 6 disappears.In addition, both the constraints associated to Higgs and twin quarkonium annihilation extend to higher m 0 values for smaller ξ, since the decay rate of the 0 ++ , dominated by Higgs exchange, decreases.Let us also comment on the dependence of these constraints on the electric charge of the twin quarks, in particular for Y t = 2/3 as it corresponds to the twin top.A larger Y q implies a larger production cross section of the twin quarks for the same mass, thus the main effect is that the constraints from searches for displaced vertices extend to heavier m q.We find that 3 ab −1 could probe up to m t ∼ 520 GeV for light glueballs or m t ∼ 450 GeV for heavy ones, close to the twin-top mass expected from twin parity; we could forecast that a dedicated search would achieve the required sensitivity.Let us add that in the above we have always assumed that two, and only two, 0 ++ glueballs are produced during twin-color glueballization.
The very slow decay of twin glueballs such as the 0 −+ , see Eq. ( 37), could potentially be detected by a surface detector such as MATHUSLA [62], granted its efficiency to identify photons is not negligible.In the right panel of Fig. 6 we show a (very) rough estimate of the number of photons that would decay inside the detector: L • σ 0 −+ • P (d), with L = 3 ab −1 , the 13 TeV cross section σ 0 −+ includes glueballs from twin gluons produced in both Higgs (ξ = 0.1) and twin-bottomonium decays, and P (d) = e −d i /d (1 − e −(d i +∆d)/d )Ω M /4π is our naive estimate of the probability for the glueball to decay inside MATHUSLA, where d i ≈ 225 m (distance from the interaction point to center of the closer horizontal edge of the detector), ∆d ≈ 45 m (the corresponding maximum distance to the farther edge of the detector), d is to be identified with the lab-frame decay length of the glueball, (cτ γ) 0 −+ , and Ω M ≈ 0.3 the solid angle covered by the detector. 20Once again we have assumed for simplicity that at least two, but only two, 0 −+ glueballs arise from the hard-scattered twin gluons.We therefore conclude that large numbers of glueballs, up to O(10 3 ) after 3 ab −1 of integrated luminosity, could potentially decay to diphotons within MATHUSLA.Interestingly, this conclusion does not depend on the presence of light twin bottoms in the spectrum, as shown by the blue contours in Fig. 6, corresponding to glueballs from the decays of Higgs and twin-top bound states at their twin-symmetric mass (ξ = 0.1).
We conclude this section by noticing that another potentially relevant signature of our glueballs is given by the process pp → gg → 0 ++ 0 −+ + X, where the 0 ++ decays displaced to γγ (even if this is a subleading decay channel, see Fig. 5), while 0 −+ escapes the detector and thus shows up as missing energy [73,74].

Fraternal Z
In the the fraternal model, the spontaneous breaking of SU (2) L × U (1) Y to U (1) Y gives rise to a tree-level contribution to Y -parameter, which decouples with both the twin gauge coupling g and the symmetry-breaking scale f as Figure 7: (Left) Upper 95% CL limits on Z cross section times branching ratio to dileptons as a function of the Z mass from a published search in 36.1 fb −1 of data at √ s = 13 TeV [52] (expected bound as black solid and observed as dashed) and from a more recent search with 139 fb −1 [53] (red), as well as the prediction for the twin Z for different SU (2) L gauge couplings: g = g, √ 2g, 2g (blue dotted, dashed, solid respectively).(Right) In the plane of Higgs decay constant f and g/g, excluded regions from Z bounds (red shaded), the Yparameter (purple shaded) and ξ < 0.1 (green shaded).Also shown contours of constant twin Z mass (blue solid) and Higgs tuning (grey dashed), the latter computed with g ρ = 6.
where on the r.h.s.we neglected subleading terms in ξ and g /g.The constraint from LEP Y 1 • 10 −3 then leads to an important upper bound on the parameter combination g4 f 2 .Nevertheless, except for very large g or unless some ad hoc assumption is made on the decays of the Z , this constraint is milder than the one arising from direct Z searches at the LHC, as shown in the right panel of Fig. 7.
We recall that the twin Z is the axial combination of the B and W 3 gauge bosons, with mass m Z ∼ gf /2.Its couplings to the SM and twin fermions can be written as where s θ = g /g and J Q, Y ,Y are the U (1) Q, Y ,Y currents made of either SM fermions ψ or twin fermions ψ, with twin parity enforcing Y ( ψi ) = Y (ψ i ) and Q( ψi ) = Q(ψ i ).The O(ξ) term stands for corrections induced after EWSB, which we safely neglected in the following.
From the last term in brackets one can immediately see that the Z can be produced from a SM fermion pair, e.g.q q at the LHC, with a cross section that scales as σ Z ∼ g 2 (g /g) 2 for fixed Z mass and g g .In this limit the width of the Z is dominated by decays to twin fermions, with branching ratios to SM fermions scaling as BR( ψψ) ∼ (g /g) 2 .Therefore, for large values of g both production and SM decays are suppressed.
We find the strongest direct constraints on the Z arise from ATLAS searches for a resonant peak in dilepton (e + e − , µ + µ − ) invariant mass distributions [52,53].The bounds on σ • BR are shown in the left panel of Fig. 7, along with the theoretical prediction for three different values of the twin SU (2) L gauge coupling and assuming one full generation of twin fermions the Z can decay to. 21In the twin-symmetric case g = g, we find m Z 4.1 TeV, a stringent bound that, according to Eq. ( 9), implies f 10 TeV, way beyond natural values.This motivates larger values of g, for which the constraints on m Z are milder and the corresponding amount of fine-tuning is smaller.This is shown in the right panel of Fig. 7, where we plot the excluded regions from LHC direct searches and from the Y -parameter Eq. (38) in the (f , g/g) plane, along with contours of fixed m Z and tuning ∆ g associated to the gauge contributions to the Higgs potential for g ρ = 6; values of ∆ g 0.5% can be achieved for g 2g.
In view of these results, the question arises of how g could happen to be larger than g at low energies.The naive possibility of decoupling two full twin generations at a high-energy scale above m * where g ≈ g, is at odds with the requirement of approximately equal color and twin-color gauge couplings at m * .Therefore, other ideas seem to be required (e.g.large threshold corrections or extra SU (2) L -charged matter).These difficulties nevertheless suggest that the easiest realization of our twin-Higgs scenario is the minimal one, in which the twin SU (2) L gauge symmetry has been broken above m * .

Conclusions
Twin Higgs models enjoy phenomenological signatures that, while linked to the naturalness of the electroweak scale, significantly differ from those of the standard TeV solutions of the hierarchy problem.Instead of collider signals associated with production of colored particles, e.g.top partners, stops and gluinos, that subsequently decay to top/bottom and Higgs/EW gauge bosons or missing energy, in the twin-Higgs scenario the phenomenology at low energies is dominated by states that predominantly couple to the SM via the Higgs, thus giving rise for instance to exotic Higgs decays or heavy Higgs-like signatures.This kind of signals are more difficult to detect at colliders such as the LHC, which is the reason why current constrains are mild.More important is however the fact that large luminosities, expected to be delivered by the LHC in the forthcoming years, could potentially set the scales even.
In this paper, we have presented a new twin-Higgs construction that nicely exemplifies the exciting prospects for discovery of these unorthodox low-energy manifestations of theories solving the hierarchy problem, where the states that account for the little hierarchy between the Higgs mass and the cutoff m * of the Higgs sector are colorless but hypercharged.This is a consequence of the symmetry structure of the strong sector, whose dynamics is ultimately responsible for shielding the EW scale from UV thresholds while giving rise to a seemingly elementary Higgs in the IR.Besides a Z 2 parity that exchanges color and its twin, a global SO(7) symmetry, spontaneously broken to the G 2 , delivers a Higgs and a twin Higgs that 21 The constraints would be mildly weaker if three twin generations instead of one were to be considered.carries hypercharge.These couple to the SM fermions and their twins, in particular to a twin top with 2/3 of electric charge, eliminating the leading sensitivity of the Higgs potential to m * .While this exceptional twin-Higgs model admits a full mirroring of the SM, we find that only when the twin gauge bosons are decoupled above m * , the Higgs VEV and mass can be reproduced with a mild level of fine-tuning, to a good approximation given by 2ξ ≈ 10%.The misalignment ξ = v 2 /f 2 between the EW scale and the Higgs decay constant is a common requisite in all the constructions with a composite-NGB Higgs, since its couplings to the SM depart from those of an elementary Higgs proportionally to ξ, and such departures are constrained by both EW precision tests and Higgs couplings measurements.This characteristic of twin-Higgs models becomes ever more important to test, something that could well be achieved at a future Higgs factory.
The exceptional twin Higgs brings about a different set of thrilling probes characterizing hypercharged naturalness.The leading indirect effect is found in the loop-induced Higgs coupling to photons, which receives a contribution from the twin top which makes it more SM-like in comparison to other composite-Higgs models.The fondness of the twins for hypercharge, that is for photons and Z's, shows up in a striking way in processes where the twins are directly produced.In the fraternal scenario, this affinity quickly drives the model into fine-tuned territory, due to the relatively light Z that couples to the SM quarks and leptons via their mixing with the B. Much more compelling is the case where we strip our construction to the bare minimum: the twin pNGBs, the twin top and possibly the twin bottom.While the neutral pseudoscalar ω 0 is almost inert, the electrically charged scalar ω ± is long-lived or stable on collider scales, giving rise to charged tracks that the LHC could soon detect.The twin top/bottom exhibits quirky phenomenology: once pair-produced they lose energy up until they are brought back together by their twin-QCD interactions, eventually forming low-lying twin-quarkonium states, either a pseudoscalar η −+ or a vector Υ −− .The former annihilates mostly to twin gluons, thus becoming another source of twin glueballs besides those from Higgs decays as in other twin-Higgs models.The latter either annihilates to SM fermions, thus it could be detected as a Z -like resonance (although with much smaller cross sections than the fraternal Z ), or to γh and Zh, which can be considered a smoking gun of our exceptional twin Higgs.The production cross sections are small but kinematically accessible at the LHC, and prospects for a detection at the HL-LHC are promising.Finally, the twin glueballs, in particular the lightest 0 ++ , exhibit similar phenomenology as those in other twin-Higgs models: they are long-lived if the twin-QCD scale is not considerably above Λ QCD .Amusingly, some of the glueballs such as the 0 −+ decay almost certainly to diphotons, and they live very long, the longer the heavier the twin bottom, a reasonable possibility being this state merely an spectator for what concerns the stability of the electroweak scale.In such a case, these glueballs would leave the LHC detectors as missing energy, and eventually reappear as a pair of photons, to be picked up by a future detector such as MATHUSLA.
These are experimental opportunities that we believe should not be missed, given the high stakes in the search for the dynamics, even if eerie, that protect us from the unknown at the highest energy scales. 22 u y p X 6 t e l 6 l 0 W R x 7 O 4 B w u w Y M b q M I D 1 K A B D B C e 4 R X e n C f n x X l 3 P p a t O S e b O Y U / c D 5 / A O R / j P 4 = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " e y t h I 2 o p g x t N i U b g r f 8 8 i p p 1 6 r e R b X W v K z U b / I 4 i n A C p 3 A O H l x B H e 6 g A S 1 g w O E Z X u H N e X R e n H f n Y 9 F a c P K Z Y / g D 5 / M H 4 X e M / A = = < / l a t e x i t > b < l a t e x i t s h a 1 _ b a s e 6 4 = " Y J j h R 7 R Y 5 h y N t V L B H / M e r r m O Q 7 I = " > A A A B 6 H i c b V B N S 8 N A E J 3 U r 1 q / q h 6 9 L B b B U 0 m q o M e i F 4 8 t 2 F p o Q 9 l s J + 3 a z S b s b o Q S + g u 8 e F D E q z / J m / / G b Z u D t j 4 Y e L w 3 w 8 y 8 I B F c G 9 f 9 d g p r 6 x u b W 8 X t 0 s 7 u 3 v 5 B + f C o r e N U M W y

4 8 t 2 F
p o Q 9 l s J + 3 a z S b s b o Q S + g u 8 e F D E q z / J m / / G b Z u D t j 4 Y e L w 3 w 8 y 8 I B F c G 9 f 9 d g p r 6 x u b W 8 X t 0 s 7 u 3 v 5 B + f C o r e N U M W y x W M S q E 1 C N g k t s G W 4 E d h K F N A o E P g T j 2 5 n / 8 I R K 8 1 j e m 0 m C f k S H k o e c U W O l 5 q h f r r h V d w 6 y S r y c V C B H o 1 / + 6 g 1 i l k Y o D R N U 6 6 7 n J s b P q D

5 9 3 5 W
L S u O f n M C f y B 8 / k D t Y O M 3 w = = < / l a t e x i t > Z < l a t e x i t s h a 1 _ b a s e 6 4 = " F O 5 N R j T 0 D a p z d a Q Z Y F w 0 P t q J G 0 U = " > A A A B 6 H i c b V D L T g J B E O z F F + I L 9 e h l I j H x R H b R R I 9 E L x 4 h k U e E D Z k d e m F k d n Y z M 2 t C C F / g x Y P G e P W T v P k 3 D r A H B S v p p F L V n e 6 u I B F c G 9 f 9 d n J r 6 x u b W / n t w s 7 u 3 v 5 B 8 f C o q e N U M W y w W s s j j y c w C m c g w d X U I U 7 q E E D G C A 8 w y u 8 O Y / O i / P u f C x a c 0 4 2 c w x / 4 H z + A L o P j O I = < / l a t e x i t > t < l a t e x i t s h a 1 _ b a s e 6 4 = " v w q O i O G F r 5 i b Q x M v l 0 o r V e 5 6 g S c = " > A A A B 7 3 i c b V B N S 8 N A E J 3 U r 1 q / q h 6 9 L B b B U 0 m q o M e i F 4 8 V 7 A e 0 o W w 2 m 3 b p Z h N 3 J 0 I p / R N e P C j i 1 b / j z X / j t s 1 B W x 8 M P N 6 b Y W Z e k E p h 0 H W / n c L a + s b m V n G 7 t L O 7 t 3 9 Q P j x q m S T T j D d Z I h P d C a j h U i j e R I G S d 1 L N a R x I 3 g 5

Figure 1 :
Figure 1: Sketch of the fraternal spectrum of the exceptional twin Higgs.
4 u a s 3 r I o 4 y O k Y n 6 A y 5 6 B I 1 0 S 1 q o T a i K E P P 6 B W 9 W U / W i / V u f S x a S 1 Y x c 4 T + w P r 8 A S 7 M k s Y = < / l a t e x i t > !+ < l a t e x i t s h a 1 _ b a s e 6 4 = " I M r J L O M p Y t E f P Z g r c h I E O a 6 r e 0 s = " > A A A B + H i c b V B N S 8 N A E N 3 4 W e t H o x 6 9 L B Z B E E p S B T 0 W v X i s Y D + g i W W z m b R L N 5 u w u x F q 6 C / x 4 k E R r / 4 U b / 4 b t 2 0 O 2 v p g 4 P H e D D P z g

Figure 2 :
Figure 2: Sketch of the minimal spectrum of the exceptional twin Higgs.