A Minimal Model for Neutral Naturalness and pseudo-Nambu-Goldstone Dark Matter

We outline a scenario where both the Higgs and a complex scalar dark matter candidate arise as the pseudo-Nambu-Goldstone bosons of breaking a global $SO(7)$ symmetry to $SO(6)$. The novelty of our construction is that the symmetry partners of the Standard Model top-quark are charged under a hidden color group and not the usual $SU(3)_c$. Consequently, the scale of spontaneous symmetry breaking and the masses of the top partners can be significantly lower than those with colored top partners. Taking these scales to be lower at once makes the model more natural and also reduces the induced non-derivative coupling between the Higgs and the dark matter. Indeed, natural realizations of this construction describe simple thermal WIMP dark matter which is stable under a global $U(1)_D$ symmetry. We show how the Large Hadron Collider along with current and next generation dark matter experiments will explore the most natural manifestations of this framework.


Introduction
The Standard Model (SM) of particle physics has great agreement with experiment, however it cannot be the complete theory of nature. One of the most pressing theoretical problems within the SM is the hierarchy between the weak and Planck scales. Both composite Higgs models and constructions which protect the Higgs mass through a new symmetry predict new particles or states with masses at or below the TeV scale.
Beyond this theoretical puzzle, there is overwhelming experimental evidence for dark matter (DM) which also points to new particles and interactions beyond the SM. While there is a vast and varied spectrum of possible DM candidates, weakly interacting massive particles (WIMPs) are perhaps the most theoretically compelling. This is especially the case when viewed through the lens of the hierarchy problem. Then the DM can naturally obtain a weak scale mass and couplings, providing the observed DM density through thermal freeze-out.
However, both symmetry based explanations of Higgs naturalness and thermal WIMPs have become increasingly constrained by experiment. Searches at the Large Hadron Collider (LHC) have pushed the limits on the colored symmetry partners of SM quarks to the TeV scale. At the same time a host of direct detection experiments are driving the limits on WIMP DM cross sections toward the so-called neutrino floor. With the severity of these constraints many new and interesting ideas for both Higgs naturalness and DM have been explored.
In twin Higgs models, the Higgs is a pseudo-Nambu-Goldstone boson (pNGB) of a global SO(8) symmetry breaking to SO (7). The variety in DM candidates typically comes not from the symmetry breaking structure, but by making particular assumptions about the particle content in the twin sector. Other neutral naturalness pNGB constructions [7,11,12,16,17] employ smaller symmetry groups, but this move toward minimality makes it more difficult to accomodate simple DM candidates.
However, as demonstrated in [29,30], the six pNGBs that spring from SO(7)/SO(6) can be associated with the Higgs doublet (respecting the custodial symmetry) along with a complex scalar DM stabilized by a global U (1) D 1 . The mass of the DM and its couplings to the Higgs are determined by the symmetry breaking structure and the low energy fields that transform under the symmetry. This necessarily includes the top quark for the model to address the hierarchy problem. As a consequence, the collider bounds on colored top partners lead to couplings between the Higgs and the DM that are near or beyond the experimental limits [29,32].
In the following section we construct a neutral natural version of the SO(7)/SO(6) symmetry breaking pattern. As expected, the quark symmetry partners are charged under a hidden color group SU (3) c rather than SM color. This mean they can be much lighter, allowing for additional freedom in the Higgs non-derivative coupling to the DM. These SM color neutral top partners are electroweak charged and break the DM shift symmetry, generating the DM potential. Thus, more natural top partner masses can lead to Higgs portal direct detection signals that may not be fully explored until the next generation of dark matter experiments. However, we do find that nearly all natural parameter choices lie above the neutrino floor.
In addition, the new fields related to the top quark exhibit quirky [33] dynamics. These less studied particles can be discovered at the LHC, providing a complementary probe of the model. In Sec. 3 we outline the most promising collider searches, including both prompt and displaced signals. We find that the LHC already bounds the quirks with U (1) D charges. Because these particles determine the coupling between the Higgs and the dark matter, these collider bounds immediately inform the sensitivity of dark matter experiments to the pNGB WIMPs. We also calculate the corrections to the electroweak precision tests (EWPT) due to the presence of the new electroweak charged states.
In Sec. 4 we discuss the DM phenomenology, showing which parameter values lead to the correct thermal relic density and elucidate how direct and indirect searches probe the model. We find collider searches and direct detection experiments provide complementary probes, both delving into the natural parameter space along different directions in parameter space. While current limits allow versions of this framework with ∼ 10% tuning, next generation searches should be able to discover the quirks or DM, often in multipole channels, down to ∼ 1% tuning. Following our conclusions, in Sec. 5, we include two appendices to provide details relating to the work.
2 Neutral Naturalness from SO(7)/SO (6) In this section, we describe a neutral naturalness model which includes the Higgs doublet and a complex scalar DM candidate as pNGBs. This model is related to that of the Refs. [29,30], but crucially has color neutral top partners. The global symmetry structure is SO(7)×SU (7), where SU (7) ⊃ SU (3) c ×SU (3) c ×U (1) X includes the SM color group as well as a hidden sector color denoted SU (3) c . The additional U (1) X ensures SM fields have their measured hypercharges. At some scale f the global SO(7) symmetry is broken to SO Here the SO(4) C SU (2) L × SU (2) R is the familiar custodial symmetry with SU (2) L being the usual SM weak gauge group and SO(2) D = U (1) D is the global symmetry that stabilizes the DM. This construction also breaks the DM's shift symmetry in a new way. In particular, through color neutral vector-like quarks in addition to the color neutral top partners. As we see below, the DM mass and its non-derivative couplings are proportional to the masses of these color neutral vector-like quarks.

The Gauge Sector
We begin with the interactions amongst the NGBs and the gauge fields. The NGB fields can be parameterized nonlinearly as where Σ 0 = (0, 0, 0, 0, 0, 0, f ) T and Tâ are the broken generators of SO(7) withâ = 1, · · · , 6, see Appendix A for details. We immediately find where |π| ≡ (πâ) 2 . We can then write the leading order NGB Lagrangian as where the covariant derivative is given by Note that the electric charge of fields is defined by Q = T 3 L + T 3 R + X, or the hypercharge is defined as Y = T 3 R + X. The first four NGBs are related to the usual Higgs doublet H = (h + , h 0 ) T by In unitary gauge when h + = 0 and h 0 = h/ √ 2 we have (π 1 , π 2 , π 3 , π 4 , π 5 , π 6 ) = 0, 0, 0, h, where we have defined χ = (π 6 + iπ 5 )/ √ 2 as a complex scalar which is our DM candidate. It is convenient to make the field redefinition [34], We can then write NGB field as in unitary gauge. The NGB Lagrangian has the leading order terms When h gets a vacuum expectation value (VEV) of v ≈ 246 GeV we write 10) to ensure that h is canonically normalized. We also define s v ≡ v/f .

The Quark Sector
The quark fields include particles charged under both SM color and the hidden color group. In terms of SO(7) and SU (6) representations we have the left-and right-handed quarks as Q L = (7, 6) and T R = (1, 6). These can be split up schematically in terms of fields in the 3 of their respective color groups where we have put a hat on fields charged under the hidden color group. More explicitly we write out the low lying left-handed fields as where q L = (t L , b L ) T is the usual SU (2) L quark doublet. This is similar in spirit to Refs. [11,12,29] which use incomplete quark multiplets. One can imagine the other fields lifted out of the low energy spectrum by vector-like masses, or as in extra dimensional models [5,7] that the boundary conditions of the bulk fields are such that the zero modes vanish. In order to obtain the correct hypercharge for t L , t R , and b L , both Q L and T R have a U (1) X charge of 2/3. The Yukawa coupling term Q L ΣT R then implies that the NGBs have zero X charge, which in particular implies that χ has no SM gauge charges.
The top sector couplings follow from where q L = (t L , b L ) T and q L = ( t L , b L ) T are SU (2) L doublets and we have restored the eaten NGBs for the moment into the Higgs doublet H, and defined H = iσ 2 H * . From these interactions we obtain the one-loop diagrams in Fig. 1 relevant to the mass corrections for H and χ. The leading contributions from the top quark are doubled by q L interaction, but this combination is exactly cancelled by T L . Like in [16], the contributions from fields carrying SM color and those carrying hidden color are not equal. Note that the DM shift symmetry is broken by the SM color neutral top partner T and U (1) D charged fermions X L , Y L .
The hidden color fields can be lifted through vector-like mass terms with new heavy states. We can write down the mass terms

The Scalar Potential
We are interested in the obtaining the potentials for both the Higgs and the DM. This is obtained from the Coleman-Weinberg (CW) potential [35] where M 2 is the Dirac fermion mass squared matrices, with masses as functions of h and χ. We note first that there is no quadratic sensitivity to to the cut off because Tr M 2 is independent of the scalar fields. However, we do find logarithmic sensitivity because where we have dropped field independent terms. Any remaining terms in the scalar potential, such as quartic mixing of h and χ or a |χ| 4 term, are independent of Λ UV and so are robustly determined by the low energy physics. Clearly, in order for electroweak symmetry to break we need the Higgs mass parameter to be negative, so we require f λ t > m Q . From Eq. (2.10) we see that Higgs couplings to SM fields will be reduced by c v . As in other pNGB Higgs construction, this implies that f exceeds v by a factor of a few. As in [16] we find there must be a cancellation between independent terms (m Q and λ t f ) to obtain the correct Higgs mass. This motivates defining For simplicity, in this work we take the vector-like masses of the DM sector to be equal This mass scale is related to m Q by the ratio r Q = m 2 V /m 2 Q . In this limit we find one of the dark fermion mass eigenstates is exactly m V , while the others are determined by a cubic equation. We then find the scalar potential, which has the general form of The potential parameters are calculated from the CW potential in Eq. (2.15). We find ln We need the dark U (1) D to remain unbroken so that χ is stable. This means we are interested in vacua with χ = 0 and h = v. With µ 2 h < 0 and µ 2 χ > 0, this is the deepest vacuum as long as λ h λ χ < λ 2 hχ . However, when λ h λ χ > λ 2 hχ the vacuum with χ = 0 becomes a saddle point rather than a minimum, so the deepest stable vacuum still has χ = 0. In this case we find Since we know v 246 GeV and the Higgs mass m h 125 GeV, therefore λ h 0.13 and µ h 89 GeV. The constraints on Higgs couplings (see Sec. 3) imply that f 3v, which means δ m 1. It then makes sense to expand the potential terms to leading order in δ m . We find Here we have taken ln to be order one, as expected for a cutoff of a few TeV.
The Higgs potential has logarithmic dependence on Λ UV . This is similar to both the Twin Higgs [36] and SO(6)/SO(5) constructions [11] where sizable UV contributions lead to the correct Higgs mass. In the limit of small δ m and taking λ MS t = 0.936 we find µ 2 h ≈ −(146) 2 GeV 2 and λ h ≈ 0.13 for m Q = 800 GeV and Λ UV = 3 TeV. These are similar to the SM listed above, so we expect that a suitable UV completion, perhaps composite or holographic, can easily accommodate the measured Higgs mass.
At the same time the quartic couplings that involve χ are determined completely by the low-energy theory. Thus, we can make robust predictions about the DM without knowledge of the UV completion. In Fig. 2 we see that these quartics are order 10 −2 over a wide span of r Q . This gives the value of the DM self-interactions as well as its coupling strength to the Higgs. The value of r Q is constrained by collider production of the hidden color fermions which is taken up in the following section.

Tuning
The Higgs potential obtained above also allows us to determine tuning of the Higgs mass parameter. We use the formula where δµ 2 is the leading one-loop correction to the Higgs mass parameter (2.31) Clearly, this tuning depends sensitively on δ m , and is greatly reduced when δ m 1. It is useful to connect δ m to v/f . This is done by simply minimizing the part of the Higgs potential that depends on ln Λ UV . This leads to the relation similar to what was found in [16]. We rewrite this as to see that δ m roughly tracks the tuning required to misalign the vacuum, as it should, for it is by choosing δ m small that we obtain the correct Higgs mass. This makes clear that taking δ m small is not an additional tuning, but the only tuning required to realize the correct Higgs potential. For instance, when f/v = 3 (10) we find δ m ≈ 0.125 (0.01) which corresponds to ∼10% (1%) tuning.

Collider phenomenology
The collider signals of this model arise from the Higgs sector and the production and decay of the hidden color quirks. To determine both these effects we need the physical mass states of the hidden sector fermions. The relevant mass matrix M F is As noted in the previous section to obtain the correct Higgs mass without introducing additional fine-tuning, we require, where δ m is given in Eq. (2.33). In the following, we fix the vector-like mass for the quirk doublet m Q to the this value. Note that we can use this relation to define f/v in terms of m Q : The physical masses are obtained by diagonalizing the fermion mass matrix by the trans- The mass eigenvalues are given by and the mixing angles are In other words, θ L = −θ R ≡ θ. The unmixed states are described by Dirac fermions T ± with masses m ± , their couplings to SM fields are given in Appendix B.

Scalar Sector
Like other pNGB Higgs models we find the tree level couplings of the Higgs to SM states are reduced. In our case they are reduced by c v , which follows immediately from Eq. (2.10). This leads to the usual bound of f 3v from the LHC measurement of Higgs couplings to gauge bosons. It may also lead to interesting signals at the HL-LHC and future colliders. At the same time, the existence of new fermionic states with electric charge that couple to the Higgs amplifies its coupling to photons. As in the quirky little Higgs model [7], this pushes the rate of h → γγ closer to the SM prediction [37]. Explicitly, we find the Higgs width into diphotons is approximately In Fig. 3 we see how the production of a given Higgs to SM final state rate changes relative to the SM prediction as a function of m Q . The blue curve shows the usual result for tree level Higgs coupling deviations, while the dashed orange curve denotes the decay into two photons. We see that the latter is slightly increased relative to the other rates. However, the deviation is small enough that it would likely require a future lepton collider to measure it [38][39][40]. Current Higgs coupling measurements require this ratio be no less than 0.8, and the HL-LHC is expected to reach a precision corresponding to about 0.9 [41]. We see that these already begin to probe v/f , but do not reach beyond about 10% tuning. The Higgs also develops a loop level coupling to the gluons of the hidden QCD. Similar to coupling to the photon, we find the Higgs coupling to hidden gluons takes the form where α s = g 2 s /(4π) is the hidden sector strong coupling parameter, G a µν is the hidden gluon field strength, and This leads to the Higgs width into hidden gluons which may contribute to a detectable Higgs width at future lepton colliders.
Since the states charged with hidden color carry U (1) X charge, they are electrically charged under the SM. Bounds from LEP imply that such states cannot be lighter than about 100 GeV. Consequently, the lightest hadrons of the hidden confining group are the glueballs. The lightest glueball state is a 0 ++ and has a small mixing with the Higgs. This allows the glueballs to decay with long lifetime to SM states. From [42] we find the glueball partial width into SM states to be where m 0 is the mass of the lightest glueball, is the SM Higgs partial width for a Higgs with mass m 0 . Lattice calculations have determined 4π α s f 0 ++ = 3.1m 3 0 [43]. In addition, the exotic decays of the Higgs into glueballs with displaced decays can lead to striking signals at the LHC [44].
To be more precise we must estimate the mass of the hidden glueball. This is done by estimating the hidden scale Λ QCD using two-loop running. 2 We assume at scales near the cutoff of a few TeV the SM and hidden strong couplings become equal. Thus, we can run the SM strong coupling up to the cutoff and then run the hidden coupling down from the cutoff for a given spectrum. In Fig. 4 we find that the hidden color strong scale varies between about 4.5 to 6.5 GeV for m Q ∈ [800, 1200] GeV. This implies the lightest glueball mass, taken to be about 6.8 Λ QCD , is likely to fall between 30 and 45 GeV. Then using the glueball decay width in Eq. (3.11) we find the glueballs have a decay length of hundreds of meters. The displaced decays from these particles may be quite challenging for the ATLAS and CMS to detect, but may be detected by MATHUSLA-like detectors [46].
There may also be new scalars related to the spontaneous symmetry breaking mechanism. In weakly coupled UV completions there may be a radial mode, a scalar whose mass close to f . As has been detailed for other pNGB realizations of neutral naturalness [47][48][49][50][51][52], this scalar will have order one couplings to all the pNGBs, leading to observable signals at the LHC and future colliders. If the UV completion involves an approximate scale symmetry then a heavy scalar associated with the breaking of scale invariance, the dilaton, can have large coupling to the SM and hidden sector states [53] providing additional collider signals.

Electroweak Precision Test
Extensions of the SM are constrained by precision electroweak measurements. The constraints can be expressed in terms of the oblique parameters S, T , and U [54,55]. The contributions to U are typically small, so we only compute the contribution to S and T . These contributions arise from the new electroweak charged fermions inducing important radiative corrections to gauge boson propagators. In addition, the modified coupling of the Higgs boson to gauge bosons leads to an infrared log divergence [56]. We find the leading contributions to be where Λ UV is UV cutoff scale, θ W is the usual weak mixing angle, and the factor of N c comes from the number of dark QCD color. As expected, the contributions from vector-like fermions, X and Y , cancel as well as the power law UV divergences. These contributions are compared to the experimental fits and found to lie within the 68% and 95% allowed regions as provided by the Gfitter collaboration [57]. In

Quirky Signals
The new fermions ( T ± , X, and Y ) can be produced at colliders through Drell-Yan due to their hypercharge of 2/3. We parameterize the couplings of any fermion f to the Z boson and the photon by where c W and s W is the cosine and sine of the weak mixing angle while g and e = gs W are the weak and electric couplings, respectively. As an example, SM fermions have v f = T 3 −2Qs 2 W and a f = T 3 . We then find the partonic cross section for qq → Z, γ → f f to be where α Z ≡ g 2 /(4πc 2 W ). In Fig. 6 we see the fermion cross sections at a 14 TeV protonproton collider. We used MSTW2008 PDFs [58] and a factorization scale of √s /2. All the fermions charged under the hidden color group have masses much above 100 GeV due to LEP bounds on charged particles. The hidden confining scale is of the order of a few GeV, so we expect them to exhibit quirky [33] dynamics, which can give a variety of new signals at colliders [59][60][61][62][63][64]. After production they are connected by a string of hidden color flux which, because there are no light hidden color states, is stable against fragmentation. The quirky pair behaves as though connected by a string with tension σ ∼ 3.6 Λ 2 QCD [65], see also [66].
Much of the subsequent dynamics can be treated semi-classically. Since these quirks carry electric charge the oscillating particles radiate soft-photons, quickly shedding energy until they reach their ground state [67,68]. Annihilation is strongly suppressed in states with nonzero orbital angular momentum, so in nearly every case the quirks do not annihilate until they reach the ground state. Since the quirks are accelerated by the string tension, we can estimate their acceleration as a = σ/m f . Then, using the Larmor formula we can estimate the radiated power as where α = e 2 /(4π). The time it takes the quirky bound state to drop to its ground state is given by K/P, where K is the kinetic energy of the quirks. Taking K ∼ m f we can then estimate the de-excitation time T d as (3.17) Clearly, the de-excitation is very fast, leading to prompt annihilation Depending on the masses of the hidden b quark, the T ± could β-decay by emitting a W . When the mass splitting is small the de-excitation is faster and the states typically annihilate. However, if the splitting is large it is most likely that both top-like states transition to bottom-like states. These would then de-excite and annihilate in the same way, though there would be additional W s in the final state.
If the b quarks are not too heavy, then T ± b combinations can be produced through the W boson. If these states are similar in mass so that β-decay is slow then the bound states can lead to visible signals, like W γ resonances, with appreciable rates. This is because the electric charge of the state prevents its decay into hidden gluons. However, larger splittings allow the heavier state to decay to the lighter promptly, diluting these signals significantly.
Because the quirks are fermions there are four s-wave states, one singlet and three triplet. Following [14] we assume that each of these states is populated equally by production, so we take the total width Γ tot of the bounds state to be where Γ s and Γ t are the widths of the singlet and triplet states respectively.
For the T ± b states which carry weak isospin the dominant quirkonium decays are to W W with a branching fraction of about 75%. This comes from the chiral enhancement in this decay. This signal has been searched for at the LHC by both ATLAS [69,70] and CMS [71,72]. The next largest fractions are into Zh, at the 10% level, which can be compared to ATLAS [73] and CMS [74,75] searches. All other visible final states are suppressed well below the percent level, see Fig. 7. Of these, the most likely LHC signal is a new scalar resonance decaying to W W , though this does depend on the b-quirk mass. As shown in Fig. 8, current searches are not yet sensitive to these signals. Here we assume a production of the T − T − directly, and through production of the T + state which then decays to a soft Z and T − . While the LHC is not yet sensitive to these signals, the high luminosity run (dashed red line) will probe the most natural regions of parameter space [76].  The X and Y particles only couple to visible states through hypercharge, hence there is no rate into Zh and the rate into W W vanishes when the Z mass can be neglected. The largest coupling is to hidden gluons, so this dominates the branching fractions. These gluons shower and hadronize into hidden QCD glueballs, some fraction which may have displaced decays at the LHC [77]. However, they can also annihilate intof f and EW gauge bosons through their hypercharge coupling, see Fig. 7. Of these, dilepton and diphoton channels have the greatest discovery potential because the signal is so clean, which has motivated searches at both ATLAS [78,79] and CMS [80,81]. In the right panel of Fig. 8 we compare the reach of the ATLAS search [79] to the theoretical prediction. We see that quirks below about 550 GeV are in tension with current collider bounds. Seeing that the predicted cross section is near the experimental limit, it is likely that by the end of the LHC run 3, with 300 fb −1 , any quirks of this type below a TeV will be discovered. Further LHC runs can probe even larger m V , but we note that taking this mass larger does not affect the naturalness of the Higgs mass. It does, however, indicate that the DM is heavier, see Eq. (2.27).
When m V > m − + m χ the X, Y quirks will quickly decay, V → T − + χ. In this case the powerful dilepton resonance search will not apply. Instead, the production cross section for T − bound states must include this, in general small, additional mode. A similar story holds if m − > m V + m χ , where now the T − quirk decays promptly to an X or Y and a DM scalar. Then, the dilepton bounds would apply to the T production. For lighter m Q this can strengthen the bound on m V . The red dashed and purple dash-dotted lines on the dilepton bound in Fig. 8 correspond to taking m Q = 800 GeV and m Q = 1000 GeV, respectively, and the DM mass of 100 GeV. By taking the DM heavier these lines would cut off earlier, at m V = m ± − m χ .
In summary, standard collider searches for prompt visible objects do constrain m V > 550 GeV, but the other parameters of the model are less restricted. However, both the displaced searches related to the hidden sector glueballs and dilepton and diboson resonance searches can provide evidence for the hidden QCD sector at the LHC. As we shall see in the next section, this parametric freedom can lead to viable DM, and complementary search strategies from DM experiments.

Dark matter phenomenology
In this section we detail the phenomenology of the DM candidate χ, the complex scalar charge under the global symmetry U (1) D . As mentioned above, this global symmetry stabilizes the DM. All the SM fields and the quirky top partners T ± are U (1) D neutral, whereas the quirky fermions X and Y are charged. The U (1) D global symmetry is exact, so we can associate a discrete dark Z 2 parity under which, but more generally we simply consider particles in this sector as carrying a global dark charge, which prevents their decay. Since the quirky states X and Y have the fractional SM electric charge 2/3 they cannot be the DM. However, the SM neutral complex scalar χ is our DM candidate as long as it is the lightest U (1) D charged particle.
To determine the success of this scalar as explaining the observed DM in the universe, in what follows we calculate the relic abundance and DM-nucleon cross section for the direct detection in our model. We then consider the dark matter annihilation for the indirect detection and impose the collider constraints on our parameter space. We find that much of the natural parameter space of this model has not yet been conclusively probed by experiment, but is expected to be covered next several years.

Relic abundance
The thermal relic density of the scalar χ is obtained using the standard freeze-out analysis. Figures 9 and 10 show the relevant Feynman diagrams for the DM annihilation and semiannihilation/conversion, respectively. The Boltzmann equation for the DM annihilation and semi-annihilation/conversion processes is Figure 9: The Feynman diagrams for the DM annihilation to SM. Figure 10: The Feynman diagrams for the DM semi-annihilation/conversion through the dark quirks.
where φ(φ ) are the SM fields: h, t, W, Z, γ, · · · . Also, H is the Hubble parameter and n i is the number density of species i, whereas then i is its thermal equilibrium value. The quantity σ ijkl v Møl ≡ σ(ij → kl) v Møl is the thermal averaged cross-section of the initial states ij to final states kl with v Møl being the Møller velocity. The last term in the first line of Eq. (4.2) describes the dynamics of the standard DM annihilation to the SM final states as shown in Fig. 9. The second and third lines describe the semi-annihilation and conversion processes shown in Fig. 10.
The dominant DM annihilation channels are to the SM, i.e. χχ * → W W, hh, ZZ, tt, bb, while the semi-annihilation and conversion processes are only relevant if the masses the quirk states ( V , T ± ) are similar to m χ . When the quirk masses are much larger than the DM, their thermal distributions are Boltzmann suppressed, making semi-annihilation or conversion processes very rare as compared to the standard annihilation processes. The relevant Feynman rules to calculate the DM annihilation or semi-annihilation processes are given in Appendix B. The DM relic abundance is computed using the public code micrOMEGAs [82].
Before discussing these results we emphasis some of the features of this model.
• The top partners are SM color neutral, therefore the symmetry breaking scale f may be at or below a TeV. This leads to significant improvements in the fine-tuning while simultaneously allowing a larger window for the pNGB DM masses in comparison to colored top partner models [29,31,32].
• The DM annihilations to SM are dominated by s-channel Higgs exchange. The amplitude for such processes is, where s = 4m 2 χ . The s dependent term originates from the derivative coupling ∂ µ h ∂ µ (χ * χ), while the λ hχ term is a loop induced explicit breaking of the χ shift symmetry, see Eq. (2.19).
• When the standard DM annihilation processes dominate (which we see below is typically the case), the DM relic abundance can be estimated as, where 0.12 is the observed DM relic abundance by the PLANCK satellite [83].
• The thermal averaged annihilation cross section to SM fields via s-channel Higgs exchange is proportional to which implies that in the limit λ hχ → 0, i.e. no explicit shift symmetry breaking, the cross section is proportional to m 2 χ /f 4 . Hence, for a given m χ the relic abundance, Ω χ h 2 , scales as f 4 .
In Fig. 11, we show the relic abundance Ω χ h 2 for two benchmark values of λ hχ = 0.005 and 0.025 as a function of m χ with fixed f/v = 4, 6, 8, 10. Notice that for masses below 50 GeV the DM tends to be overproduced. This is because the thermal averaged crosssection in this region is directly proportional to the portal coupling λ hχ , which direct detection constrains to be relatively small (see below). On the other hand, for m χ ∼ m h /2 the relic abundance drops sharply due to the resonant enhancement of the Higgs portal cross-section. For DM masses m 2 χ ∼ λ hχ f 2 /2 there is cancelation in the cross-section as a result the relic abundance enhances which produces the peaks in Fig. 11. For larger DM masses the cross section is proportional to m 2 χ /f 4 and the relic density drops as DM mass increases.
For the case λ hχ = 0.005 (left-panel), the relic density curves terminate when the DM becomes heavier than the quirk states X, Y . These states are bound by the dark color force into quirky bound states, which then efficiently annihilate due to their electric charge, making them an unsuitable thermal DM candidate. There is also a sharp drop in the relic density at the end of each curve, which is due to an s-channel resonant enhancement of semi-annihilation processes, i.e. χ V → T ± → T ± SM, as shown in Fig. 10. The semiannihilation processes are only significant when m χ ≈ m V ≈ m ± /2 and in most of the parameter space are inefficient as compared to the standard annihilation processes. Since the portal coupling λ hχ is proportional to r Q = m 2 V /m 2 Q it can be reduced for relatively light vector-like quirks V . However, collider searches at the LEP and LHC put a lower bound these vector-like quirks, see Sec. 3.3. We see that the smallest mass that produces the correct DM thermal relic is near the Higgs resonance region, above ∼ 50 GeV. This is fairly independent of f/v and λ hχ . However, the largest DM masses which leads to correct relic abundance does depend on f/v and λ hχ . Since naturalness prefers a smaller f/v and λ hχ is constrained by direct detection (see below), we find that restricting f/v ≤ 10 puts an upper bound of m χ 1 TeV for obtaining the correct relic.

Direct detection
The WIMP DM scenario is being thoroughly tested by direct detection experiments. We here highlight the main features of our pNGB DM construction where direction detection null results are explained naturally.
At tree-level the DM-nucleon interaction is only mediated by t-channel Higgs exchange. As discussed above, the DM-Higgs interaction has two sources: (i) the derivative coupling ∼ (∂ µ h)∂ µ (χ * χ)/f 2 , and (ii) the portal coupling ∼ λ hχ hχ * χ. The strength of the derivative interaction in a t-channel process is suppressed by the DM momentum transferred, t/f 2 ∼ (100 MeV) 2 /f 2 1. For all practical purposes we can neglect such interactions. Hence the only relevant interaction for direct detection is the portal coupling λ hχ . 3 In this case, the spin-independent DM-nucleon scattering cross-section σ SI χN can be approximated as (see e.g. [29,31]), where m N is the nucleon mass and f N 0.3 encapsulates the Higgs-nucleon coupling. The current bound on the spin-independent DM-nucleon cross-section for mass range ∼ [50,1000] GeV is by XENON1T with one ton-year of exposure time [84]. For instance, the upper limit on the spin-independent DM-nucleon cross-section for DM mass 300 GeV is χ . Hence to satisfy the direct detection constraints we either need to reduce the portal coupling λ hχ or increase the DM mass.
One feature of this minimal model is that λ hχ is determined by a small number of low-energy parameters: the vector-like masses of the quirks, m V and m Q . However, as noted above in Eq. (3.2), the top partners quirk mass m Q = c v λ t f is fixed in terms of f to obtain the correct Higgs mass. Hence, the free parameters are m χ , f , and r Q ≡ m 2 V /m 2 Q . As discussed above, one can specify f by requiring the correct DM relic abundance and r Q can be exchanged with λ hχ , which is constrained by direct detection.
In Fig. 12 we show the spin-independent DM-nucleon cross section σ SI χN as a function of DM mass m χ . We have performed a random scan of the parameter space for f/v ∈ [3, 10] and m V ∈ [m χ , 4πf ]. The lower value of the f/v = 3 choice is enforced by the SM Higgs coupling measurement and electroweak measurements data, while the upper value of f/v = 10 limits the tuning to ∼1%. The lower value of m V makes sure that χ is the lightest state charged under U (1) D . All the points shown in the plot correspond to the correct relic abundance Ω χ h 2 = Ω obs h 2 ± 5σ, where Ω obs h 2 = 0.12 ± 0.0012 is the observed DM relic density as measured by the Planck satellite [83]. The gray (pentagon) points above the gray line are excluded by the XENON1T [84]. All the colored points (color barcoded with f/v) are allowed by the current XENON1T constraint. The dashed gray line indicates the expected XENONnT bound [84] which covers much of the more natural parameter space. However, there are points allowed below this bound above the so-called neutrino-floor (red dotted), which could be discovered by next generation detectors, e.g. LZ [85] and DARWIN [86].

Indirect detection
We now turn to indirect detection. There are a variety of experiments searching for DM annihilations in the Milky Way galaxy and nearby dwarf galaxies, which are assumed to be dominated by DM. The typical signals of DM annihilation to the SM particles leads to gamma-rays, gamma-lines, and an excess of secondary products like antipositrons and antiprotons in cosmic-rays (CR). In particular, the experimental data can be used to put upper bounds on the various annihilation channels, including W W, ZZ, hh, tt, bb, τ + τ − , · · · . In our model the DM dominantly annihilates into W W, hh, ZZ, tt final states. We calculate the present day DM thermal averaged annihilation to the SM particles at zero velocity by using micrOMEGAs [82]. We find that the DM thermal annihilation cross-section is σv ≈ 2.2 × 10 −26 cm 3 /s for parameter values that produce the correct relic abundance. The fraction of annihilation cross-section to W + W − is ∼ 45% and hh/ZZ ∼ 25% for m χ m h . Whereas the branching fraction is dominantly into In Fig. 13 we show the DM annihilation cross section to W + W − , σv W W , in units of [10 −26 cm 3 /s] as function of m χ . All the data points in this figure produce correct DM relic abundance and satisfy the XENON1T direct detection constraint. Because these points have m χ > m W , the most dominant annihilation channels are the W W, ZZ, hh. In the following we summarize the most sensitive indirect detection probes in the mass range of interest.
Gamma-rays: The most robust indirect detection bounds are due to Fermi-LAT [87] and Fermi-LAT+DES [88] with six years of data from 15 and 45 DM dominated dwarf spheroidal galaxies (dSphs), respectively. Theses constraints are considered robust because the uncertainties associated with propagation of gamma rays are relatively small. The Fermi-LAT results [87] provide upper-limits on the DM thermal annihilation cross section into several SM final states including W W, bb, τ + τ − , whereas, the updated analysis Fermi-LAT+DES [88] only includes the bb and τ + τ − channels. These bounds do not constraint any of the parameter space allowed by the direct detection. However, Fermi-LAT has provided expected 95% C.L. upper-limits for the DM thermal annihilation into bb and τ + τ − channels with 15 years of data and 60 dSPhs [89]. One can interpolate the projected upper-limit from the σv bb to σv W W by a simple rescaling σv W W 1.33 σv bb in the DM mass range of our interest. In Fig. 13 we show the projected 95% C.L. sensitivity on σv W W by Fermi-LAT with 15 years and 60 dSPhs by the solid (gray) curve. This sensitivity sets a lower-limit on the DM mass m χ 150 GeV.
Cosmic-rays: The flux of antipositrons and antiprotons in the cosmic-rays (CR) provides another indirect probe of DM annihilation in the Galaxy. In particular recent precise AMS-02 CR antiproton flux data [90] has led to strong constraints on the DM thermal annihilation. In Refs. [91,92] the AMS-02 antiproton flux data was used to put stringent constraints on DM with masses in range [150, 1000] GeV. The AMS-02 95% C.L. exclusion constraint on σv W W as obtained by CCK [91] is shown in Fig. 13 as dash-dotted (blue) curve. This constraint excludes most of the data points between DM masses m χ ∈ [225, 375] GeV. However, these constraint has large systematic uncertainties, mainly due to CR propagation and diffusion parameters [91]. The updated analysis by (CHKK) [93] reveals a weaker constraint in the W + W − channel, which is also given by a dash-dotted (blue) curve. Even though the updated AMS-02 analysis does not constrain our model, future AMS CR antiproton data are likely to. Another future CR experiment is the Cherenkov Telescope Array (CTA) which is expected to be sensitive to large DM masses [94]. In Fig. 13 we show the projected sensitivity of CTA for DM annihilation to W + W − with Galactic Diffuse Emission (GDE) Gamma model of Ref. [95], as a dashed (red) curve, for two assumptions of systematic error. The most optimistic implies that CTA will probe DM masses above ∼ 300 GeV, though this is quickly weakened when systematic errors are included.

Conclusion
We have outlined a framework in which the Higgs and a scalar DM candidate arise pNGBs of a broken global symmetry. Because the symmetry partners of the top quark do not carry SM color, the induced scalar potential between the Higgs and the DM, which is UV insensitive, allows for improved fine-tuning and simultaneously explains null results for WIMP DM searches. The quantitative success of this framework is summarized by Fig. 14 in the m χ vs λ hχ plane with the color of scanned points corresponding to values of f/v ∈ [3,10]. This corresponds to fine-tuning in the model of about 10% to 1%, respectively.
The phenomenology can be specified by the DM mass m χ , the global symmetry breaking scale f , and the vector-like mass m V of the quirky fermions, which is the source of breaking the χ shift symmetry. As shown in Sec. 2.3 we can trade m V for λ hχ . Hence, the three free parameters of the model are m χ , f/v, and λ hχ .
The points in Fig. 14 scan in m χ ∼ [50, 1000] GeV and λ hχ ∼ [0.2, 0.0005] while required to produce the correct relic abundance Ω χ h 2 = 0.12 ± 5(0.012). The gray (pentagon) points are excluded at 90% C.L. by the direct detection experiment XENON1T with one year exposure time [84]. Future direct detection XENONnT 90% C.L. reach is overlaid as the dash-dotted (black) curve, which would cover much of the allowed parameter space. Next generation experiments that will descend toward the neutrino floor will fully explore this framework.
The next most stringent constraint is due to the LHC bound on the vector-like mass m V 550 GeV of the quirky fermions X, Y as shown in Fig. 8. This limit from the ATLAS collaboration search for dilepton resonances with 139 fb −1 data is due to the annihilation of quirks V V to + − . We show the bound in Fig. 14 as red (hexagon) points. Since the portal coupling λ hχ is proportional to m 2 V , the lower-bound on m V translates to a DM mass and f/v dependent lower-bound on λ hχ . We have also shown dashed (red) contours of m V = 1 TeV to 10 TeV which shows how future LHC runs may be able to discover quirks in much of the natural parameter space. The complementarity between collider and direct detection could lead to both discovery and confirmation of this construction in the coming years, or its exclusion.
In Fig. 14 we also show how indirect detection gamma-rays 95% C.L. constraints from the Fermi-LAT 15 years with 60 dSphs as blue (star) points. This puts a lower limit on the DM mass m χ 150 GeV. We have not shown in this plot the indirect detection constraints from the cosmic-rays experiments AMS-02 because of their large systematic uncertainties. However, in the future such uncertainties may be reduced, allowing experiments like AMS and CTA) to provide another complementary probe, and hopefully discovery, of this model.
In summary, this framework of WIMP dark addresses the hierarchy problem without colored symmetry partners, and consequently is only tuned at the 10% level while agreeing with all experimental bounds. However, existing experiments will soon be able to discover or exclude these more natural realizations of the model. After the searches of the HL-LHC run and next generation direct detection experiments models with fine tuning at or better than 1% may be thoroughly probed.

A SO(7) Generators
In this appendix we collect all the relevant details. The SO(7) generators in the fundamental representation can be written as, where i, j = 1, . . . , 7. We have chosen the normalization Tr T a T b = δ ab . The unbroken generators T a L,R ij , T ab ij correspond to the SO(6), whereas the broken generators Tâ ij correspond to the SO(7)/SO(6) coset. Note that T a L,R ij correspond to the custodial SO(4) C ∼ = SU (2) L × SU (2) R subgroup of SO(6).

B Feynman rules and Quirk Processes
In this appendix we record formulae for quirk production and decay widths. The relevant Feynman rules are given in Table 1. The decays are typically similar to the results [96,97], using the methods outlined in [98,99]. The couplings of the Z to fermions are taken to be where c W ≡ cos θ W . For convenience we define the following where M is the mass of the relevant bound state. The number of colors in the quirk confining group is N c . We calculate the cross section pp → Z, γ → f f from the quark q initiated partonic cross sectionσ into a quirk Q pair by Table 1: Some of the most relevant Feynman rules of our model are listed in this table, see the text for the corresponding notation.
is defined in terms of the MSTW2008 PDFs [58] f q (x), we take the factorization scale to be √s /2.
Because the quirk states decays from all > 0 states are strongly suppressed [33] we only consider decays of the singlet 1 S 0 and triplet 3 S 1 states. Each of these decay widths depends on the radial wavefunction R(0) of the quirk bound state. This factor is nonperturbative and not exactly known, so we simply give each decay width in units of the unknown factor |R(0)| 2 .
The neutral states are composed of fermionic quirks Q with mass m Q . In this case the Z couplings are labeled v Q and a Q , and the electric charge is denoted Q Q . 4 The mass is denoted m Q and we take the meson mass to be M , which for heavy constituents is approximately 2m Q .
We begin with decays to fermion pairs. These fermions have Z couplings v f and a f as well as electric charge Q f . They also come in N c colors. The decays to f f are, where α W ≡ g 2 /(4π). Next, we turn to decays into Zγ, The decays to ZZ 5 , Next, to Zh, where λ Q is the Yukawa coupling of the quirks to the Higgs. Finally, to hγ, (B.14) One might expect decays to scalar pairs like hh and, in the case of the X and Y quirks, χχ * . However, CP and angular momentum conservation forbid such decays from the s-wave states, though higher angular momentum states do allow these decays. We now turn to decays into W + W − . We label the SU (2) L partner of Q by q, with mass m q etc. The W couplings v W and a W are defined by the interaction We note that this decay depends upon the electric charge of particle that makes up the bound state in a nontrivial way. This is due to the diagrams related to the tor u-channel exchange of the SU (2) L partner of the particle making up the bound state. Mesons made by a quirk with positive charge involve a different diagram than those with negative charge. None of these subtleties affect the singlet case, but we do distinguish the triplet cases as 3 S (+,−) 1 , where the superscript denotes whether the quirk has positive or negative electric charge. The decays to W + W − are We also record the decays involving hidden gluons. These are taken from [100].
where we have denoted the hidden sector strong coupling by α s . Finally, the singlet state can also decay to photons