--- title: 'SIMP Dark Matter: Number-Changing Dynamics' url: https://www.emergentmind.com/topics/strongly-interacting-massive-particle-simp type: topic --- # SIMP Dark Matter: Number-Changing Dynamics Strongly Interacting Massive Particle (SIMP) denotes a dark-matter paradigm in which the relic abundance is set primarily by dark-sector number-changing reactions rather than by conventional \(2\to2\) annihilations into Standard Model (SM) states. In the canonical formulation, the dominant process is \(3\to2\), while symmetry-restricted variants can instead realize \(4\to2\). The resulting thermal history is characterized by strong self-interactions, cannibal heating, and a central requirement that the dark sector remain in kinetic equilibrium with a heat sink during freeze-out. In standard radiation domination this often points to MeV–GeV dark matter, but hidden sectors, resonant mediators, and non-standard pre-BBN cosmologies substantially broaden the viable parameter space [2410.02871, 1811.02751]. ## 1. Thermal mechanism and cannibal dynamics The defining kinetic equation of the SIMP framework is the number-changing Boltzmann equation. In the notation used for general \(r\to2\) freeze-out, with \(r=3,4\), \[ \frac{dn_\chi}{dt}+3H n_\chi=-(r-2)\,\langle \sigma_{r\to 2} v^{r-1}\rangle\left(n_\chi^r-n_\chi^2 n_{\chi,\rm eq}^{\,r-2}\right), \] so the relic density is controlled by dark-sector self-annihilations rather than by annihilation into SM particles [2410.02871]. For \(3\to2\), freeze-out occurs when \(n_{\rm eq}^2\langle \sigma_{3\to2}v^2\rangle \simeq H\); for \(4\to2\), when \(n_{\rm eq}^3\langle \sigma_{4\to2}v^3\rangle \simeq H\) [2410.02871]. A distinctive thermodynamic feature is cannibalization. Because a \(3\to2\) or \(4\to2\) reaction converts rest mass into kinetic energy of the remaining dark particles, the dark sector cools more slowly than an ordinary nonrelativistic gas unless it can dump entropy into another bath. In a hidden-sector SIMP with \(T_D\neq T\), the nonrelativistic cannibal phase obeys \[ T_D \simeq \frac{m_D}{3\log(a/\bar a)}, \] and the entropy-conserving yield is \[ Y_{\rm ent}=\frac{T_D}{\xi(m_D+\tfrac52 T_D)}, \] with \(\xi\equiv S/S_D\) the SM-to-dark entropy ratio [1803.07518]. This already shows that SIMP freeze-out is not only a question of cross sections; it is equally a question of thermal contact. The usual contrast with WIMPs is therefore incomplete unless the thermal structure is stated explicitly. WIMP relic density scales inversely with a \(2\to2\) annihilation cross section into SM states, whereas SIMP relic density is controlled by dark-sector number-changing processes and by whether kinetic equilibrium is maintained long enough for the cannibal heat to be removed. A common symmetry-driven variant replaces \(3\to2\) by \(4\to2\): in \(\mathbb Z_2\)-stabilized models, vertices with an odd number of dark fields are forbidden, so \(3\to2\) is absent and \(4\to2\) dominates [1510.08527]. ## 2. Dark pions, chiral dynamics, and vector resonances The best-developed ultraviolet picture of SIMP dark matter is a confining gauge theory in which the dark matter states are pseudo-Nambu–Goldstone bosons. In the dark-pion construction based on \[ SU(3)_L\times SU(3)_R \rightarrow SU(3)_V, \] the low-energy theory is a dark chiral perturbation theory with decay constant \(f_\pi\), and the Wess–Zumino–Witten term generates the five-point interaction responsible for \(3\to2\) [1811.02751]. In the pion-only effective theory, the relevant anomalous operator is \[ - \frac{2N_c}{15\pi^2 f_\pi^5}\epsilon^{\mu\nu\rho\sigma} {\rm Tr}(\pi \partial_\mu \pi\partial_\nu \pi\partial_\rho \pi\partial_\sigma \pi). \] This is the canonical origin of SIMP number-changing dynamics in composite models [1811.02751]. The same framework also makes clear why pure-pion SIMPs are often under pressure. To reproduce the observed relic density using only the WZW interaction, the chiral expansion parameter is driven toward the perturbative boundary. A practical criterion emphasized in the dark-pion literature is \[ \frac{m_\pi}{f_\pi}\sim 2\pi, \] and the WZW-only case violates this bound over the mass range relevant for Bullet Cluster and small-scale structure considerations [1811.02751]. In a related hidden-sector analysis of composite SIMPs, the model-independent parametrization \[ \langle \sigma v^2\rangle_{3\to2}=\frac{\alpha_{\rm eff}^3}{m_D^5},\qquad \frac{\sigma_{\rm scatter}}{m_D}=\frac{a^2\alpha_{\rm eff}^2}{m_D^3} \] makes the same tension manifest: large \(3\to2\) rates and acceptable \(2\to2\) self-interactions are difficult to reconcile if the dark sector temperature tracks the SM exactly [1803.07518]. A controlled resolution is to include dark vector resonances. In the hidden local symmetry construction, the vectors are massive gauge bosons of gauged \(SU(3)_V\), with \[ a \equiv \frac{m_V^2}{(g f_\pi)^2},\qquad m_V=g f_\pi\sqrt a. \] The gauged WZW sector then generates vector-mediated five-point topologies, and near the “3-pion resonance” \[ m_V=3m_\pi\sqrt{1+\epsilon_V}, \] the \(3\to2\) amplitude is resonantly enhanced while the underlying chiral couplings remain moderate [1811.02751]. The narrow-width requirement \[ \frac{m_V\Gamma_V}{9m_\pi^2}\lesssim 0.1 \] implies \(a_{\rm max}\simeq 0.1\) near the \(3\pi\) resonance [1811.02751]. This unitarization mechanism is quantitatively significant. With vectors, “in the mass range of light dark matter \((m_\pi\lesssim 1~{\rm GeV})\),” one can maintain \[ \frac{m_\pi}{f_\pi}<6\ \text{for}\ a=0.1,\qquad \frac{m_\pi}{f_\pi}<4.5\ \text{for}\ a=0.01, \] and “parameter space is wider than the only WZW case” [1811.02751]. The related analysis of dark vector resonances in hidden local symmetry reaches the same conclusion: resonances at \(m_V\simeq2m_\pi\) or \(3m_\pi\) enhance \(3\to2\), reduce the required \(m_\pi/f_\pi\), and bring \(\sigma/m_\pi\) into the \(0.1\)–\(1~{\rm cm^2/g}\) range without leaving the domain of validity of dark ChPT [1801.07726]. A second way to regain perturbative control is to cool the dark sector relative to the SM. In hidden-sector composite SIMPs, a temperature ratio \(R_T\equiv T_D/T<1\) lowers the required \(3\to2\) strength. Explicitly, for \(R_T\simeq 1/3\), a viable mass range \(m_\pi\simeq 40~{\rm MeV}\)–\(1~{\rm GeV}\) overlaps \(\sigma/m\in\{10,1,0.1\}~{\rm cm^2/g}\), while for \(R_T\simeq1/9\) the viable range shifts to \(m_\pi\simeq9\)–\(140~{\rm MeV}\) [1803.07518]. This is one of the cleanest demonstrations that the standard “sub-GeV and strongly coupled” SIMP expectation is not universal, but thermal-history dependent. ## 3. Symmetry classes and representative model realizations Symmetry determines which number-changing channel is even available. Bernal and Chu’s \(\mathbb Z_2\) singlet-scalar construction is the canonical \(4\to2\) example: the potential \[ V=\mu_H^2|H|^2+\lambda_H|H|^4+\mu_S^2S^2+\lambda_S S^4+\lambda_{HS}|H|^2S^2 \] forbids odd-\(S\) interactions, so \(3\to2\) is absent and \[ \langle\sigma v^3\rangle_{4\to2}\sim \frac{27\sqrt{3}}{8\pi}\frac{\lambda_S^4}{m_S^8} \] sets the relic density, while \[ \frac{\sigma_{SS}}{m_S}\simeq \frac{9}{8\pi}\frac{\lambda_S^2}{m_S^3} \] controls self-interactions [1510.08527]. In that scenario, a colder dark sector at freeze-out, with \(T/T'|_F\sim15\)–\(35\), reduces the required \(\lambda_S\) and helps reconcile relic density with cluster bounds [1510.08527]. Other symmetry choices support \(3\to2\) while embedding SIMPs in broader BSM structures. A radiative neutrino-mass model with a \(\mathbb Z_5\) symmetry realizes resonantly enhanced \(3\to2\) through a complex scalar \(X\) and mediator \(S\), while the same Yukawa sector generates neutrino masses at two loops and keeps the dark sector in kinetic equilibrium through \(X\nu\leftrightarrow X\nu\) scattering [1705.00592]. Twin Higgs constructions identify the SIMP states with twin pions in a QCD-like sector, with the WZW term providing \(3\to2\) and a kinetically mixed twin photon maintaining kinetic equilibrium; the preferred regime is \(m_\pi\sim\) few \(\times 10^2~{\rm MeV}\) and \(m_\pi/f_\pi\sim2\pi\) [1805.09345]. A different ultraviolet realization replaces non-Abelian confinement with a strongly interacting \(U(1)\) sector and monopole condensation. In that construction, monopole condensation confines \(U(1)_H\)-charged matter, hidden pions emerge as the low-energy composites, and the radial monopole mode mixes with the Higgs so that no additional singlet portal is required [1606.01628]. The model is notable because baryons would remain \(U(1)_H\)-charged, so there are no low-energy baryon states competing with the pion SIMP. Two-component SIMP sectors introduce additional dynamical structure. In the accidental \(\mathbb Z_4\) model with a complex scalar and a vector-like fermion, \(3\to2\) interactions determine freeze-out, but an unavoidable two-loop induced \(2\to2\) process redistributes the component abundances after chemical freeze-out and “would significantly modify the predictions of the self-interacting cross section of DM compared with other SIMP models” [2201.06856]. A related 2024 analysis considers a pFIMP in the presence of a SIMP and shows how SIMP–pFIMP conversion interpolates between a pure SIMP–FIMP limit and a SIMP–pFIMP limit in which DM–DM conversion reshapes the relic composition [2411.15108]. These models make explicit that “SIMP” can denote a thermal mechanism rather than a single-particle species. ## 4. Portals, kinetic equilibrium, and experimental probes Kinetic equilibrium is structurally central because, as one dark-pion study puts it, “one of the pion masses goes to kinetic energy of dark matter in the \(3\to2\) channel” [1811.02751]. Portal design is therefore not optional. Higgs- and singlet-scalar portals, \(Z'\) portals, axion-like particles, and neutrino-coupled mediators have all been used to transfer entropy from the dark sector to the SM bath. The dark-photon portal is the most extensively developed experimentally. In “SIMP Spectroscopy,” a kinetically mixed vector boson maintains thermal contact while allowing mono-photon missing-mass spectroscopy at \(e^+e^-\) colliders. The recoil photon energy directly measures the dark invariant mass, \[ E_\gamma=\frac{s-M_{\rm inv}^2}{2\sqrt{s}}, \] so the dark resonance structure can, in principle, be reconstructed from the mono-photon spectrum [1512.07917]. The 2023 “SIMPly add a dark photon” analysis pushes this idea further by asking whether one dark photon can simultaneously thermalize the dark sector, resonantly enhance \(3\to2\), and produce velocity-dependent self-interactions. The answer is model-dependent: for \(N_f=3\) the minimal setup is “marginally excluded, as the required kinetic mixing is too small to maintain thermal equilibrium with the SM,” whereas adding an extra dark quark yields viable models with \(m_\pi\sim250\)–\(600~{\rm MeV}\) [2301.04513]. Axion and ALP portals were developed precisely to solve the heat-dump problem without sacrificing the SIMP mechanism. In the ALP-assisted model, semi-annihilation \(\pi\pi\to\pi a\) dominates the relic-setting dynamics, while the ALP maintains SM contact through \[ \mathcal L_{a\gamma}=\frac{1}{4f_{a\gamma}}\,a\,F^{\mu\nu}\tilde F_{\mu\nu}, \qquad \Gamma_a=\frac{m_a^3}{64\pi f_{a\gamma}^2}, \] and the region yielding \(\sigma_{\rm self}/m_{\rm DM}\simeq0.1\)–\(1~{\rm cm^2/g}\) was identified as reachable by SHiP [1704.04505]. A later axion-portal construction based on \(Sp(2N_c)\) dark pions likewise found viable parameter space when the ALP mass is close to the SIMP mass, with strong-scale masses of order a few hundred MeV and near-future coverage by beam-dump and collider experiments [1806.10139]. Higgs-portal thermalization remains viable in singlet-scalar and related models. In the reheating-era singlet-scalar analysis, direct detection, invisible Higgs decays, and future collider sensitivity are all mapped simultaneously: present LZ limits already exclude parts of the allowed SIMP region, while DARWIN/XLZD and precision Higgs measurements at HL-LHC and FCC-ee will probe additional regions [2603.05590]. In the radiative neutrino-mass model, kinetic equilibrium is instead maintained by Yukawa scatterings off neutrinos until freeze-out [1705.00592]. The portal question is thus inseparable from experimental strategy: electron recoils, mono-photons, missing momentum, invisible Higgs width, and displaced low-mass dileptons all arise as portal-specific probes. ## 5. Non-standard cosmologies and the enlargement of mass scales Standard radiation domination leads to stringent mass limits from the combination of relic density and unitarity. For scalar SIMPs freezing out through \(3\to2\) or \(4\to2\), the radiation-dominated upper bounds quoted in the reheating analysis are \[ m_\chi\lesssim1~{\rm GeV}\quad(3\to2),\qquad m_\chi\lesssim7~{\rm MeV}\quad(4\to2), \] assuming s-wave scalar dark matter and nonrelativistic freeze-out [2410.02871]. These values are often taken as indicative of the SIMP scale. Reheating changes that conclusion because the Hubble rate and entropy evolution differ from radiation domination. During a reheating phase driven by a monomial inflaton potential \(V(\phi)\propto\phi^k\), the paper derives different temperature and Hubble scalings for fermionic and bosonic reheating channels, and shows that entropy injection dilutes the final yield. Consequently, “a smaller cross-section” is needed than in the radiation-dominated case, and the unitarity bounds relax dramatically. For quadratic inflaton potential \((k=2)\), the upper mass limits become \[ m_\chi\lesssim10^6~{\rm GeV}\quad(3\to2),\qquad m_\chi\lesssim10^4~{\rm GeV}\quad(4\to2), \] whereas for quartic potential \((k=4)\) the allowed space is much smaller, with \(m_\chi\lesssim300~{\rm GeV}\) for bosonic reheating and \(m_\chi\lesssim68~{\rm GeV}\) for fermionic reheating in the \(3\to2\) case [2410.02871]. A complementary result was obtained in the reheating-era singlet-scalar study of \(4\to2\) SIMP dark matter with \(T_S=T\). In standard radiation domination, this realization “is not viable, as it requires sub-MeV masses and large quartic couplings in tension with bounds on dark matter self-interactions.” During early matter domination, however, the altered Hubble rate and entropy injection make the model viable over an enormous range: \[ 0.1~{\rm GeV}\lesssim m_S\lesssim10^{12}~{\rm GeV}, \] with \[ 5\times10^{-3}\lesssim\lambda_S\lesssim10,\qquad 10^{-5}\lesssim \lambda_{HS}\lesssim O(1), \] subject to BBN, perturbativity, and self-interaction bounds [2603.05590]. This does not erase the standard sub-GeV SIMP regime; it shows instead that it is a statement about cosmological history, not a theorem about number-changing dark matter. ## 6. Conceptual tensions, terminology, and current status Two persistent tensions define the modern SIMP literature. The first is the relic-density versus perturbativity problem: in pure WZW dark-pion models, the \(3\to2\) rate needed for \(\Omega_{\rm DM}h^2\simeq0.12\) frequently drives \(m_\pi/f_\pi\) beyond the reliable regime of chiral perturbation theory [1811.02751]. The second is the heat-dump problem: cannibalization requires kinetic equilibrium with some bath through freeze-out, but that same portal must not reintroduce dominant WIMP-like \(2\to2\) annihilation into the SM. Vector resonances, cooler hidden sectors, ALP portals, and non-standard cosmologies should therefore be regarded not as optional embellishments, but as the principal mechanisms by which the simplest SIMP tension is resolved. This is also why “minimal” constructions can fail. The dark-photon-only setup studied in 2023 is explicitly “marginally excluded” for \(N_f=3\) because the kinetic mixing that preserves the relic density is too small to thermalize the dark pions with the SM, while \(N_f=4\) opens viable space and yields acceptable models for \(m_\pi\sim250\)–\(600~{\rm MeV}\) [2301.04513]. The same study emphasizes that late-time annihilations are non-negligible, making the dark pion “a bit WIMPy” [2301.04513]. That formulation captures a broader lesson: once resonant portals are added, the sharp WIMP–SIMP dichotomy becomes less rigid at the level of phenomenology even when the relic density is still primarily controlled by number-changing interactions. A separate source of confusion is terminological. In the thermal-relic literature, SIMP almost always means a dark matter candidate whose abundance is set by \(3\to2\) or \(4\to2\) reactions inside the dark sector. In collider and direct-detection phenomenology, however, the same acronym is also used for dark matter with large nucleon cross sections. Surface \(\nu\)-cleus data exclude SIMP–nucleon cross sections up to approximately \(10^{-27}\,{\rm cm}^2\) for masses above \(100~{\rm MeV}\) [1708.01484], and trackless-jet searches at the LHC were argued to tentatively exclude roughly \(\sigma_{\chi N}\gtrsim10^{-28}\,{\rm cm}^2\) for \(m_\chi\gtrsim1~{\rm GeV}\) in simplified models [1503.05505]. These studies concern a different interaction regime from the standard thermal SIMP mechanism, even though the acronym is shared. The present status is therefore plural rather than singular. SIMP dark matter encompasses dark-pion WZW freeze-out, \(\mathbb Z_2\) \(4\to2\) singlet scalars, hidden-sector cannibal models with \(T_D