Dark Matter in Very Supersymmetric Dark Sectors

If supersymmetry exists at any scale, regardless of whether it is restored around the weak scale, it may be a good symmetry of the dark sector, enforcing a degeneracy between its lowest lying fermions and bosons. We explore the implications of this scenario for the early universe and dark matter, as well as the corresponding signatures. In particular we show that the thermal history of the dark sector results in co-decaying dark matter in much of the parameter space. This implies new phenomenological signatures and presents a new way to discover high scale supersymmetry.


I. INTRODUCTION
The identity of dark matter (DM) has long been a mystery, as has the question of how its observed relic abundance came to be. The lack of observation of weakly interacting massive particle (WIMP) dark matter at operating direct detection experiments has led to a surge in activity in this field, with many proposed new theoretical ideas for dark matter beyond the WIMP (see, e.g., Refs. [1][2][3][4][5][6][7][8][9][10]). Many new experimental ideas to detect such dark matter have similarly been proposed (see, e.g., Refs. [11][12][13][14][15][16][17][18][19]). Supersymmetry (SUSY) has been a focal point for particle physics research since its conception, though null results from the LHC have somewhat shifted the image of what a supersymmetric solution to the hierarchy problem would look like. For the past 35 years, the study of dark matter with supersymmetry has focused on the lightest superpartner (LSP) of the minimal supersymmetric Standard Model (MSSM) as a dark matter candidate, as it can be absolutely stable due to R-parity conservation [20].
In this work, we take a different approach and show that if there is a hidden sector that is approximately supersymmetric, with an interaction portal to the MSSM, then there are profound implications for the thermal history of dark matter. Approximately supersymmetric sectors have been considered in the context of WIMP dark matter in Refs. [21][22][23][24], and the collider physics implications of such sectors have been studied extensively in stealth supersymmetry [25][26][27]. Such sectors can undergo different thermal histories depending on the size of the coupling between the visible and dark sectors; as we will show, this results in codecaying dark matter [9] in much of the parameter space.
In the codecaying dark matter framework, recently proposed in Ref. [9], nearly degenerate particles of a dark sector are in chemical and thermal contact with each other, while the lighter of the states undergoes a slow, out-ofequilibrium, decay into Standard Model (SM) particles. Supersymmetric theories are therefore a natural habitat for the codecaying mechanism, since the near degeneracy of the dark particles is ensured by approximate supersymmetry within the hidden sector.
It is clear that supersymmetry, if it exists at some scale, is badly broken in the visible sector. It may, however, be nearly preserved in a hidden sector. Supersymmetry breaking in a hidden sector can naturally be small if it is mediated from the supersymmetry-breaking sector to the MSSM, but not directly to the dark sector, e.g., via gauge mediation. The dark sector then feels supersymmetry breaking through the portal with the MSSM and through supergravity effects. Provided that the portal is weak, and the gravitino mass is much smaller than the supersymmetric masses in the hidden sector, the dark-sector states will remain approximately degenerate within their supersymmetric multiplets.
Conventional models of gauge-mediated supersymmetry breaking suffer from the gravitino problem, where the gravitino abundance would overclose the Universe unless all the superpartners are close to the tera-electron-volt scale or the scale of reheating is low [28,29]. As we will show, in models undergoing codecay, this tension is alleviated through the entropy dump of the hidden sector into the SM particles. This substantially reduces the predicted relic abundance of the gravitinos, thus opening up the possibility of gauge-mediated split supersymmetry [30] with a higher reheat temperature. This paper is organized as follows. In Sec. II, we discuss some generalities of approximately supersymmetric sectors. Section III reviews the ingredients of codecaying dark matter. An explicit model is presented in Sec. IV, with Sec. V studying the parameter space. We conclude with a discussion in Sec. VI.

II. VERY SUPERSYMMETRIC SECTORS
We begin by summarizing the general features of having a very supersymmetric hidden sector, which we take to be a sector with small boson-fermion splittings within a supersymmetric multiplet. To best characterize the dynamics of such theories, we introduce the relevant scales and schematic shown in the left panel of Fig. 1. The hidden sector has a self-coupling κ and mass m and is connected to the visible sector through a small coupling ε. The SUSY breaking is mediated to both the visible and hidden sectors through M Pl -suppressed operators and to the SM through its gauge interactions at a scale M.
The thermal history of the dark sector and whether there exists a viable dark matter candidate will depend on the couplings. Motivated by the null observation of any baryonor large lepton-number violating processes, we focus on models in which R-parity conservation is assumed. Thus, there is a dark matter candidate in the dark sector if where M SUSY is the mass scale of the superpartners in the visible sector. For the hidden sector to be approximately supersymmetric, the breaking must be sufficiently small, and gravity mediation must be subdominant. This imposes the requirements (note that we use ε schematically here, though we give it an explicit meaning in the context of a model later) where m 3=2 is the gravitino mass. Under these conditions, we call the dark sector very supersymmetric. This leads to the spectrum shown schematically in the right panel of Fig. 1.
The hidden sector will have a lightest R-even state and R-odd state, which we label ϕ andψ, respectively. While R-parity protectsψ from decaying to SM particles, the mass hierarchy of Eq. (2) enables it to decay via couplings to the gravitino, and thereforeψ is only an approximately stable dark matter candidate. These decays ofψ proceed into the gravitino and either an on-shell or off-shell ϕ, depending on the mass hierarchy between the ϕ andψ states: ψ →Gðϕ ðÃÞ → ffÞ. The decay width ofψ is suppressed by either the small mass splitting or by the small coupling to the SM. As a result,ψ is long lived on cosmological scales and is a viable DM candidate.
The R-even state, ϕ, will decay into SM particles via the ε portal to the visible sector, The hidden sector is coupled to the visible sector through an ε-suppressed portal. The hidden sector feels SUSY breaking weakly-indirectly, through the visible sector, and directly, only through M Pl -suppressed operators. Right: Schematic particle spectrum, illustrating the mass relations between particles of the visible sector (SM þ MSSM), the dark-sector lightest multiplet (ϕ;ψ), and the gravitinoG. Note that ϕ is the R-parity even dark-sector particle, which couples directly to the SM, andψ is R-parity odd.
Effectively, the important interactions between the hidden and visible sector can arise from this decay. In essence, the setup we have just described contains all the ingredients of a codecaying dark matter setup: two nearly degenerate particles, interacting with each other and one of them decaying into SM particles. The dynamics in the dark sector and which process sets the relic abundance of the dark matter both strongly depend on the size of these decays.
Regardless of the existence of a dark matter candidate, this setup will generically alleviate the gravitino abundance problem. Since the lightest states in the dark sector are long lived, the dark sector will dominate the energy density of the Universe after they become nonrelativistic. This leads to an early period of matter domination, which has important effects on small scale structure [10]. Once the dark-sector decays, it reheats the SM bath, diluting existing relics. The cosmological gravitino problem-that gravitinos are typically overproduced in the early Universe [28]-is thus addressed.

III. CODECAYING DARK MATTER
In this section, we briefly review the dynamics of decaying dark sectors with emphasis on the codecaying regime. A detailed study of the thermal history is contained in Ref. [9]. Codecay takes place when the lowest lying states in a hidden sector are degenerate and in thermal contact, and part of that sector is unstable. The density of the stable states will track that of the unstable states as they decay away, until the interactions between the two freeze out, and the stable states remain as relics. Codecaying DM will be a generic feature of nearly supersymmetric dark sectors, in which R-parity stabilizes one or more of the light states, while the others decay via a weak portal to the MSSM sector.
The relic density for a codecaying sector can be given in terms of the standard sudden freeze-out approximation, where g ⋆ is the number of relativistic degrees of freedom, x ≡ m=T, with m the mass of the degenerate states and T the SM bath temperature, and hσvi f is the thermally averaged annihilation cross section between the stable and unstable states.
Here, subscripts f refer to quantities evaluated at the time of freeze-out, T ¼ T f . For s-wave scattering of degenerate states, the cross section approaches a constant, and therefore the thermally averaged cross section is proportional to the average velocity. In terms of the on-threshold cross section σ, where x 0 ¼ m=T 0 and T 0 is the temperature of the dark sector. Using this, the relic abundance of Eq. (4) is If the decay Γ of the unstable state is faster than the Hubble rate at T 0 ¼ m, Γ ≳ H m , then the dark sector will decay and freeze out while in equilibrium with the SM, in which case the thermal history resembles that of a WIMPlike thermal relic, and The phenomenology of such a scenario was studied in detail in Ref. [31].
Another possibility is that in the early Universe the dark sector decouples from the SM bath. The heavier states freeze out and decouple, leaving just the lightest states in the sector. Assuming the decays of the unstable light states occur out of equilibrium, Γ ≪ H m , the comoving number density does not deplete when the nearly degenerate light states become nonrelativistic. Instead, the exponential suppression of the number density is delayed until the unstable states decay. After decays begin (at a later time), the stable population continues to follow the unstable population until the system can no longer maintain chemical equilibrium. At this point, the stable population freezes out, and the unstable particles decay away. The relic density of the stable particles will depend on the cross section, σ, between the stable and unstable states, as well as the lifetime of the decaying states, Γ −1 .
The evolution of the dark matter density of the stable stateψ and the long-lived state ϕ are given by the Boltzmann equations, where hΓ ϕ i TðT 0 Þ is the thermally averaged decay rate at the given temperature. Here, n, ρ, and P are the number density, energy density, and pressure, respectively. These equations can be solved, along with the Friedman equations, to obtain the relic abundance ofψ, along with the entropy dump of the out-of-equilibrium ϕ decays, for given mass and interactions. We solve these equations numerically for the explicit model given below. Useful analytical approximations that roughly capture the behavior of the relevant quantities can be found as follows. For an estimate of the relic abundance in terms of the parameters, one must determine the temperatures of the dark sector and SM at the time of freeze-out. The dark sector undergoes several stages of evolution: (1) Cannibalization.-When the dark sector becomes nonrelativistic, it undergoes cannibalization as the result of 3 → 2 processes. (2) Nonrelativistic cooling.-After 3 → 2 processes freeze out, the dark sector will redshift as a noninteracting nonrelativistic particle. (3) Decay.-When the age of the Universe reaches the average lifetime of the decaying states, i.e., Γ ¼ H, the dark matter particles will begin to decay away. (4) Freeze-out.-The decays will stop depleting the DM abundance when the dark-sector interactions freeze out. The temperatures at freeze-out can be obtained by evolving the dark-sector density and temperatures through these different stages in the dark matter evolution.
Using conservation of entropy and the second law of thermodynamics, one can show that the number density of the dark sector after it has become nonrelativistic up until the onset of decays (m < T 0 < T 0 Γ ) is where ξ is the ratio of dark to SM degrees of freedom at decoupling and s is the entropy density of the photon bath [9]. During this period, the dark matter can come to dominate the energy density of the Universe, and therefore the Hubble scale has important contributions from both the energy density of the radiation and the nonrelativistic dark matter: We refer to the temperatures at which cannibalization ends as T c and T 0 c . These temperatures can be found from the instantaneous freeze-out approximation, and from Eq. (8), using μ c ¼ 0. Here, subscripts c refer to quantities evaluated at T ¼ T c . Next, the dark sector redshifts as a decoupled nonrelativistic species until decays begin, roughly when The temperatures and chemical potential are simply given by Equations (12) and (13) can be solved for x Γ and x 0 Γ .
Finally, the temperature at freeze-out can be determined using the instantaneous freeze-out approximation for the 2 → 2 annihilations and the number density for a decaying particle. Further details can be found in the original work [9]. We find that and the dark sector approximately continues to redshift until freeze-out, at Equations (6)- (15) can be solved for the relic abundance when the decays occur out of equilibrium. As stated above, however, in presenting our results in the next sections, we use the full numerical computation.

IV. EXPLICIT MODEL
We now present an explicit realization of a very supersymmetric hidden-sector model and explore all the codecaying parameter space. The model forms a useful benchmark to demonstrate in detail how codecay occurs within supersymmetric theories, including both in-and outof-equilibrium decays.
Consider the MSSM coupled via a small coupling ε, to a chiral superfield Φ with self-coupling, κ, and mass parameter, m. The superpotential is where We follow the notation conventions of Ref. [32]. Note that a similar Lagrangian is considered in the NMSSM [33][34][35] and in stealth-supersymmetry models [25][26][27], but here the Φ field is not responsible for generating the μ-term, and the collider signatures are not stealthy. Contained in the Φ superfield are the R-even complex scalar ϕ ¼ 1 ffiffi 2 p ðϕ R þ iϕ I Þ and the R-odd fermionψ. The fermion,ψ, will be taken to be the lightest R-odd particle (up to the gravitino). After the scalar Higgses acquire their vacuum expectation values (VEVs), the portal interaction will induce mixing between theψ and the MSSM neutralinos and between the ϕ and the MSSM Higgses. Since we consider the portal to be very weak, ε ≪ 1, the hiddensector mass eigenstates will contain a small amount of MSSM fields, so we will continue to call the mass eigenstates by ϕ R , ϕ I , andψ. Full details of the mixing can be found in Appendix A.
All the ingredients for the codecay mechanism are present in the superpotential (16). The supersymmetric nature of the dark sector imposes approximate degeneracy between the fermion and scalar fields; R-parity stabilizes the fermion; the self-coupling κ induces interactions within the dark sector that keep it in equilibrium; and the small portal coupling, ε, will result in the scalars decaying to the SM. We review these ingredients here and leave the detailed calculations for the Appendix.
Mass splitting between the scalars and the fermion is generated by the mixing with the Higgs sector and via supergravity. A correction of OðεÞ for the fermion mass is also generated, coming from the VEV of ϕ and the κΦ 3 interaction. The masses of the hidden-sector particles at this order are given by Note that in general the gravity-mediated effects are not precisely calculable without a complete model of the SUSY-breaking sector; above, we simply show the order of magnitude expected in a generic model. We provide more details above Eq. (25) in Sec. V B. For the results in Sec. III to apply, the mass splitting between the fermion and scalars Δm must be smaller than the temperature at freeze-out, jΔmj < T 0 f . If the mass splitting is larger, the model can still exhibit a viable dark matter candidate, but the parameters differ from those previously described.
For m < m h ; m H ; m A , the scalars decay to SM pairs via mixing with the MSSM Higgs bosons, The trilinear superpotential coupling, κΦ 3 , generates self-interactions within the dark sector. Theψψ → ϕϕ annihilation cross section is given by and is responsible for maintaining chemical equilibrium within the dark sector until freeze-out. The interaction also generates 3 → 2 annihilations that can cannibalize the dark sector when nonrelativistic. There are many different possible 3 → 2 processes, involving different numbers of scalars and fermions, some of which are shown in Fig. 2. The temperature when cannibalization ends can be estimated from the total 3 → 2 thermally averaged rate, averaged over all fermions and scalars that can appear in the initial and final states. We find this to be The temperature of freeze-out of the 3 → 2 processes can then be determined by the instantaneous freeze-out condition, n 2 ψ;ϕ hσ 3→2 v 2 i ¼ H. Assuming the mass splitting is small, the relic abundance depends on m, κ, and ε, with mild dependence on the standard MSSM parameters. 1 In presenting our results throughout this paper, we choose a representative set of values for these, but note that the qualitative features are unchanged for variations in these values. In Fig. 3, we show the κ vs ε relation that gives the observed dark matter relic The relic abundance depends on the MSSM parameters tan β , m A , and μ through the dark scalar decay rate, Γ ϕ . We find that varying tan β from 5-50 and m A from 700-10; 000 GeV results in maximal deviations of Oð2%Þ of Γ ϕ . The sensitivity of Γ ϕ to μ is more significant: Γ ϕ ∝ μ 2 for large μ [see Eqs. (B1) and (B2)], but the qualitative features remain unchanged. abundance, for several values of m, as indicated by the solid colored curves. The gray shaded region indicates where decays occur after neutrino decoupling and is thus excluded by big bang nucleosynthesis (BBN); see the discussion below.

V. CONSTRAINTS
Having presented the model, we move to discuss the constraints on the parameter space. The important parameters are m 3=2 , m, ε, and κ, with the supersymmetric parameters taken illustratively as the benchmark values in Eq. (22). To obtain the correct relic abundance, 2 we takẽ ψ to make up the entirety of dark matter, which fixes the value of κ in terms of ε and m, as demonstrated in Fig. 3. The resulting parameter space is presented in Fig. 4 in terms of the ðε; mÞ parameter space. Below, we describe the various constraints of each region.

A. Lifetime of ϕ
The decays of the scalars can occur either in or out of equilibrium. If either scalar has a sufficiently large decay rate [see Eq. (19)], it will bring the entire sector into equilibrium with the SM. The differentiation between these two regions can roughly be estimated by comparing the decay width to the Hubble parameter when the temperature is of order the particle's mass, Γ ϕ ≃ 0.1HðmÞ. This distinction is represented by the blue dashed curve in Fig. 4: above (below) this curve, ϕ decays roughly occur in (out) of equilibrium.
If the decays are active during BBN, this will dump entropy into the SM as well as induce photodisintegration of light nuclei. Both these processes are heavily constrained by measurements of the ratios of abundances of light elements today [36]. To this end, we require that both scalars decay before neutrino decoupling, T ≃ MeV. The resulting constraint is indicated by the shaded gray region labeled "BBN" in Fig. 4 (and by the shaded gray region in Fig. 3).

B. Lifetime ofψ
Requiring the mass splitting between the fermion and the scalars of the dark sector is small imposes m 3=2 ≪ m. Thus, the gravitino is the LSP, and R-parity can only render the dark matter candidateψ approximately stable. The light gravitino induces the decays of the dark matter: three-body decaysψ →Gðϕ Ã → ffÞ and two-body decaysψ →Gϕ (if mψ > m 3=2 þ m ϕ ). The decay rates are given by These decay rates are observable both cosmologically and today. Decays of particles between recombination and reionization will dump energy into the SM bath, which can induce distortions to the cosmic microwave  4. Constraints on the parameter space. The gray shaded region indicates where decays occur after neutrino decoupling and is excluded by BBN. The blue dashed curve indicates the transition between in-equilibrium decays (above this curve) and out-of-equilibrium decays (below the curve). The green shaded region shows where the mass splitting (ignoring gravity contributions) is larger than the freeze-out temperature, in which case the relic abundance analysis is altered from that presented in this work. The red (purple) shaded region indicates where telescope experiments searches would have observed DM decay (annihilations).
background [37]. Furthermore, unstable dark matter decay in DM-rich regions can be observed using telescope experiments, placing a constraint on its decay rate [38,39]. In general, both these probes depend on the final state products and the mass of the decaying particle. Which probe is most stringent depends on the mass ofψ and the particular decay mode; however, we find that, across the range of interest, constraints from the isotropic gamma ray background measured by FERMI-LAT [40] are the most powerful [38].
The dominant fermion decay rate depends on the mass splitting in the dark sector. From Eq. (18), the splitting between masses squared of the hidden-sector fermion and scalars is with the last term above indicating the gravity-mediated effects, which are not precisely calculable without a complete model of the SUSY-breaking sector. However, there are calculable unavoidable contributions from anomaly mediation. Given the dimensionful coupling in the superpotential m and including the conformal compen- Additional terms of the same order will arise from the supergravity Lagrangian [26]. Generically, the two-body decay is fast if kinematically allowed; it is shut off if The irreducible contributions are of order this size, but we require that the sum of these contributions adds to satisfy Eq. (26) such that the two-body decay is turned off and the three-body decay dominates. The dominant constraints in our mass region of interest are then from decays toψ →Gbb and related fermionic channels (we sum all relevant channels when estimating the constraints). Such decays have similar kinematics to those of DM →bb decays, which have been well studied in the literature, and we use the results of Ref. [38] to estimate the constraints. Noting that the three-body decay rate (denoted here simply by Γψ ) is inversely proportional to the square of the gravitino mass-namely that Γψ m 2 3=2 is constant-the lifetime constraint can be written as a lower bound on the gravitino mass m 3=2 in terms of the other parameters, where Γ lim is the observed limit, and is a function of the dark matter mass. If m 3=2 is too large, a large mass splitting between the hidden-sector particles will be induced, which renders the results of Sec. III inapplicable. For simplicity, we therefore choose to focus on the region of parameter space where the splitting is smaller than the dark-sector temperature at freeze-out, Δm < T 0 f . To achieve this without too much tuning, we require Combined, to satisfy both Eqs. (27) and (28), we impose that As long as this constraint is satisfied, there exists a value of m 3=2 where the mass splitting is sufficiently small, and the three-body decay rate is sufficiently slow. The red shaded region in Fig. 4 is ruled out by a failure to meet this constraint. The nonsmooth behavior of the constraint for small ε is from rapid changes in T 0 f as cannibalization stays on the entire time or not depending on m.

C.ψψ annihilations
In addition to slow dark matter decays, telescope experiments also have the ability to detect DM annihilations which cascade into SM states throughψψ → ϕ R=I ϕ R=I . The produced final SM state depends on the mass, with the dominate state beingbb for the range of interest. We recast previous limits on direct DM annihilations [38,41], while rescaling the rates and kinematics to match a 2 → 2 process to a 2 → 4 process. Due to the inevitable velocity suppression in the annihilation cross section, we use the constraints from diffuse emission instead of dwarf galaxies or the cosmic microwave background, which are typically stronger but have a larger velocity suppression in rates.
The constraints are shown in the purple shaded region in Fig. 4. Unfortunately, the recast constraints for the final states of interest are only available for dark matter masses above 10 GeV (which corresponds to 20 GeV for our 2 → 4 processes). Nevertheless, since the experiments detect emission well below this threshold, it should be possible to extend these constraints in future work.

D. Gravitino abundance
A necessary requirement of very supersymmetric sectors is the lightness of the gravitino; otherwise, there will be inevitable large splitting in the hidden sector. Light gravitinos are a form of hot dark matter and cannot make up a substantial fraction of dark matter without being ruled out by the observation of structure at small scales. Unfortunately, gravitinos are typically overproduced in what is known as the cosmological gravitino problem [28] (see, e.g., Ref. [29] for a recent review).
Postinflationary gravitino production occurs through two main processes: freeze-in production from superpartners dropping out of the SM bath and decays of the (second) lightest R-odd states in the spectrum into the gravitino. The gravitino problem for very supersymmetric sectors is alleviated in two different ways. First, by assumption that the DM is roughly cosmologically stable, there may not be significant gravitino production from decays. As discussed in Sec. V B, this decay rate can be heavily suppressed if ϕ is heavier thanψ. Second, if freeze-out occurs while the hidden sector is out of equilibrium with the SM [Γ ≪ HðmÞ], this can result in an early period of (dark) matter domination. Once the decays become efficient, they reheat the SM bath, diluting existing relics, Y f ≃ Y Γ S Γ =S RH , where Y ≡ n=s is the comoving number density.
The ratio of entropy after reheating to the onset of the decay can be estimated using energy conservation as (see, e.g., Ref. [42]) Entropy of the dark sector is conserved from the end of cannibalization (T c 0 ) until the decays begin. Using conservation of entropy then allows us to relate the ratio n Γ =s Γ to T c 0 leading to where T 0 c is the temperature in the dark sector when cannibalization stops (namely 3 → 2 process freeze-out). Taking the extreme case that T 0 c ∼ m and the decay happens at T ∼ MeV, one can see the entropy dump can be as large as approximately m=MeV. For a generic codecaying model, this can in principle be as large 10 7 [9], significantly diluting any preexisting gravitino abundance.
The initial gravitino abundance can be estimated using standard techniques [29]. In Fig. 5, we show the constraints on the parameter space for different dilution factors assuming negligible production of gravitinos from the epoch of reheating after inflation. Note that codecaying dark matter greatly reduces the gravitino abundance, allowing for relatively high scale inflation and superpartners well above the tera-electron-volt scale. This opens up the possibility of gauge-mediated split supersymmetry [30]. Note that the entropy dump solution to the gravitino problem presented here does not depend on reproducing the relic abundance. In terms of the concrete model used here, the entropy dump has only a mild dependence on the parameters of the theory through the log-dependence of T 0 c .

VI. CONCLUSIONS
We have shown that very supersymmetric dark sectors, in which there is near degeneracy between bosons and fermions, can have distinct thermal histories which differ significantly from an ordinary WIMP. In much of the parameter space, the dark matter relic abundance is set by a codecay process, where the decays of the scalar occur either in or out of equilibrium. We presented a toy model which exemplifies all the features of such sectors and studied the resulting constraints on the parameter space. Furthermore, we have demonstrated how very supersymmetric dark sectors can alleviate the gravitino problem via large entropy dumps, which are natural in such thermal histories.
These supersymmetric dark sectors have distinct phenomenology and offer a new way to discover SUSY through the energy, lifetime, indirect detection, and cosmological frontiers.
At colliders, SUSY dark sectors will generically yield long-lived particles. If kinematically accessible, colored particles introduced to solve the hierarchy problem will be produced at the LHC. These will decay into the hidden sector through the small coupling ε between the dark and visible sectors, rendering these decays long lived. A distinct feature of such sectors is decays both into the hidden sector and back again, which results in two displaced vertices; see Fig. 6. These vertices will have correlated decay lengths as they are both suppressed by the small coupling ε. For the parameter space considered here, the lifetimes of ϕ vary from a millimeter to effectively stable on collider scales. The lifetime for the lightest observable SUSY particle is model dependent, but a similar range of lifetimes is expected. With the proposal of new experiments designed to discover longlived particles such as MATHUSLA [43], FASER [44], CODEX-b [45], and AL3X [46], the LHC could be primed to see one displaced vertex, while the future detector could see the second, longer, displaced vertex.
Telescope experiments are also uniquely poised to observe very supersymmetric dark sectors via two distinct signatures. First, dark matter is predicted to inevitably decay to gravitinos and SM particles, with rates predicted by Eq. (23). These decay rates depend on the couplings to the SM, the mass of the hidden sector, and the gravitino mass. In addition, one can also observe the annihilation process which governed the thermal freeze-out as in Eq. (20). The annihilation rate depends only on the coupling within the hidden sector and its mass, which offers a complementary probe of the supersymmetric nature of the sector. An observation of both these signals could in principle be differentiated by the energy dependence of the final state.
The above signatures persist for both the in-and out-ofequilibrium regions of parameter space. However, a unique prediction of the out-of-equilibrium codecaying regime is the possibility of a period of matter domination before the onset of BBN. This will produce enhanced power at small scales in the matter-power spectrum. As a result, compact dark objects can form, which could be seen through indirect detection [47], microlensing [48,49], astrometry [50], and pulsar timing arrays [51]. Unlike the collider signatures described above, enhanced power at small scales could be present even if SUSY is present well above the weak scale; this forms a unique opportunity to probe high scale supersymmetry.