We discuss a minimal realization of the strongly interacting massive particle (SIMP) framework. The model includes a dark copy of QCD with three colors and three light flavors. A massive dark photon, kinetically mixed with the Standard Model hypercharge, maintains kinetic equilibrium between the dark and visible sectors. One of the dark mesons is necessarily unstable but long-lived, with potential impact on CMB observables. We show that an approximate "isospin" symmetry acting on the down-type quarks is an essential ingredient of the model. This symmetry stabilizes the dark matter and allows to split sufficiently the masses of the other states to suppress strongly their relic abundances. We discuss for the first time the SIMP cosmology with sizable mass splittings between all meson multiplets. We demonstrate that the SIMP mechanism remains efficient in setting the dark matter relic density, while CMB constraints on unstable relics can be robustly avoided. We also consider the phenomenological consequences of isospin breaking, including dark matter decay. Cosmological, astrophysical, and terrestrial probes are combined into a global picture of the parameter space. In addition, we outline an ultraviolet completion in the context of neutral naturalness, where confinement at the GeV scale is generic. We emphasize the general applicability of several novel features of the SIMP mechanism that we discuss here.

Dark matter could emerge along with the Higgs as a composite pseudo-Nambu-Goldstone boson $\chi$ with decay constant $f\sim \mathrm{TeV}$. This type of WIMP is especially compelling because its leading interaction with the Standard Model, the derivative Higgs portal, has the correct annihilation strength for thermal freeze-out if $m_\chi \sim O(100)$ GeV, but is negligible in direct detection experiments due to the very small momentum transfer. The explicit breaking of the shift symmetry which radiatively generates $m_\chi$, however, introduces non-derivative DM interactions. In existing realizations a marginal Higgs portal coupling $\lambda$ is generated with size comparable to the Higgs quartic, and thus well within reach of XENON1T. Here, we present and analyze the interesting case where the pattern of explicit symmetry breaking naturally suppresses $\lambda$ beyond the reach of current and future direct detection experiments. If the DM acquires mass from bottom quark loops, the bottom quark also mediates suppressed DM-nucleus scattering with cross sections that will be eventually probed by LZ. Alternatively, the DM can obtain mass from gauging its stabilizing $U(1)$ symmetry. No direct detection signal is expected even at future facilities, but the introduction of a dark photon $\gamma_D$ has a number of phenomenological implications which we study in detail, treating $m_{\gamma_D}$ as a free parameter. Complementary probes of the dark sector include indirect DM detection, DM self-interactions, and extra radiation, as well as collider experiments. We frame our discussion in an effective field theory, motivating our parameter choices with a detailed analysis of an $SO(7)/SO(6)$ composite Higgs model, which can yield either scenario at low energies.

We consider a composite model where both the Higgs and a complex scalar $\chi$, which is the dark matter (DM) candidate, arise as light pseudo Nambu-Goldstone bosons (pNGBs) from a strongly coupled sector with TeV scale confinement. The global symmetry structure is $SO(7)/SO(6)$, and the DM is charged under an exact $U(1)_{\rm DM} \subset SO(6)$ that ensures its stability. Depending on whether the $\chi$ shift symmetry is respected or broken by the coupling of the top quark to the strong sector, the DM can be much lighter than the Higgs or have a weak-scale mass. Here we focus primarily on the latter possibility. We introduce the lowest-lying composite resonances and impose calculability of the scalar potential via generalized Weinberg sum rules. Compared to previous analyses of pNGB DM, the computation of the relic density is improved by fully accounting for the effects of the fermionic top partners. This plays a crucial role in relaxing the tension with the current DM direct detection constraints. The spectrum of resonances contains exotic top partners charged under the $U(1)_{\rm DM}$, whose LHC phenomenology is analyzed. We identify a region of parameters with $f = 1.4\; \mathrm{TeV}$ and $200\;\mathrm{GeV} \lesssim m_\chi \lesssim 400\;\mathrm{GeV}$ that satisfies all existing bounds. This DM candidate will be tested by XENON1T in the near future.