Getting a THUMP from a WIMP

Producing an acceptable thermal relic abundance of dark matter with masses $\gg 10^2$ TeV is a challenge. We propose a novel mechanism where GeV-scale states establish a tiny thermal relic abundance for dark matter, which is later promoted to ultra massive status by a very light scalar. We refer to this dark matter as a THermal Ultra Massive Particle (THUMP). Direct detection of THUMPs can be naturally expected due to large scattering cross sections mediated by low mass states that couple THUMPs to the Standard Model. Our model generically leads to signals for the associated GeV-scale states at accelerator experiments.


INTRODUCTION
The search for clues to the identity of dark matter (DM) continues. So far, no experiment has uncovered conclusive evidence regarding the properties of DM and the range of possibilities for what it could be remains vast. In this situation, one may be compelled to try and look for potential candidates wherever possible. However, launching a search in every possible direction is not feasible and very often one needs some motivating factor to justify the effort. Here, theoretical considerations could point to promising targets.
The weak scale has been, for a long time, considered a likely place for new physics to emerge, in large part due to arguments based on a natural Higgs sector for the Standard Model (SM). This scale is also typical for production of thermal relic DM, a fact that has strongly motivated searches for DM particles of corresponding mass. However, after decades of dedicated searches for both new weak scale particles and DM at accelerators and direct detection experiments, no conclusive evidence of either has been uncovered. While this may point to the elusiveness of such new states or the limitations of current experimental techniques, one could also start to entertain new or less examined possibilities.
Apart from the well-established weak scale, the scale of a possible quantum theory of gravity is another conceptually plausible place where new phenomena may arise. This scale is taken to be set by the Planck mass M Pl ≈ 1.2 × 10 19 GeV, which is also typically assumed to mark the end of short distance physics. Between the weak and Planck scales there could be several other motivated scales where new phenomena may appear. For example, unification of known forces of Nature in a Grand Unified Theory often points to masses 10 15 GeV for new states. Also, the strong CP problem of QCD can be elegantly addressed through the Peccei-Quinn mechanism which typically points to scales f P Q 10 12 GeV, * email: hooman@bnl.gov † email: gmohlabeng@bnl.gov to avoid conflict with cosmology ("overclosing" the Universe) [1][2][3]. However, astrophysical constraints on the associated axion demand f P Q 10 9 GeV [4]. So, 10 9 GeV f P Q 10 12 GeV provides another motivated window at very high scales (the axion could be a good DM candidate at the upper end of the window).
Although the above high scales provide motivation to consider DM masses well above the weak scale, DM relic abundance is generally assumed to be non-thermal for masses above ∼ 100 TeV [5]. In principle, one could construct such models of ultra massive DM, but generically the underlying physics and the signals of this DM may not be accessible, due to the high scales involved [6][7][8][9]. Different mechanisms for setting the thermal relic abundance of ultra massive DM, leading to various astrophysical signatures have been considered in the literature, see for example [10][11][12][13]. It would be interesting to consider possible mechanisms for setting the thermal relic abundance of ultra massive DM that would naturally lead to accessible signals in direct detection experiments or at accelerators. Such a scenario may be motivated by possible multi-scatter [14][15][16][17] or single scatter signals in current or near future detectors.
In this work, we examine whether ultra massive DM particles could be produced initially as a thermal relic, but with a much smaller mass. This would suggest that the thermal relic abundance is set to very small values, which will yield the correct energy density once DM attains a large mass later. We propose a scenario in which the late time slow-roll of an ultra-light modulus scalar field drives the DM mass to values 100 TeV. We note here that similar scenarios were constructed in Refs. [18,19], but with different phenomenology and physics outcomes.
Intuitively, the small relic abundance implies very large annihilation cross sections. We will illustrate through a simple model that this typically points to an initial DM mass GeV and associated dark states that mediate "secluded" annihilations. However, once the DM becomes very heavy, it maintains its strong coupling to light mediators which would naturally lead to detectable signals in the laboratory, such as through direct detection of DM or searches for the mediators in accelerator experiments.
We will refer to our proposed DM candidate as a THermal Ultra Massive Particle (THUMP).
In the following sections, we outline a simple model that realizes the above scenario using a Dirac fermion for DM and a light (sub-GeV) vector particle, which we take to be a "dark photon" coupled to the SM through kinetic mixing. One could invoke a late phase transition (For very recent work in this direction, see Refs. [20][21][22]) or a very light modulus that effects the mass variation of the DM particle from light (GeV scale) to very heavy.
Here, we will focus on the second possibility, which will typically lead to a long range force acting on THUMPs and could possibly offer further astrophysical signatures.

A SPECIFIC MODEL
Let us first consider a DM particle χ with mass m i χ in the early Universe, i.e. before Big Bang Nucleosynthesis (BBN). Here, χ is assumed to have charge Q D under a dark U (1) interaction mediated by A µ of mass m A m i χ . We also suppose that the DM χ interacts with a very light scalar φ. The Lagrangian for these interactions is given by In Eq. (1), φ sources a large mass m f χ m i χ for the DM once it starts evolving in its potential. This mass is given by the vacuum expectation value (vev) φ = φ 0 which yields m f χ = m i χ + λφ 0 . The idea here is that dark matter freeze-out is governed by m i χ with φ frozen near φ = 0, leading to an efficient annihilation; we will later discuss how this could be arranged due to thermal effects in our scenario. During this time, the thermal relic abundance of DM is set by the annihilation of χ into pairs of A through a large coupling g D Q D (represented by Fig.  1). The vector A then decays into SM particles before BBN; this is thus a secluded annihilation scenario [23]. The coupling of A to the SM will be assumed to result from kinetic mixing between U (1) and SM hypercharge U (1) Y [24]: where ε 1 is the kinetic mixing parameter, θ W is the weak mixing angle, and F µν refers to the associated field strength tensor for each U (1) interaction. For the rest of this discussion, we set Q D = 1.
After DM has frozen out and the relic abundance has been set, the φ field rolls and oscillates, where its expectation value increases the DM mass to m f χ . Since m A m i χ , the annihilation represented by Fig. 1 is 1: Annihilation of DM χ into "dark photon" pairs, which later decay into SM states (only t-channel shown).
very efficient, especially for g D 1, resulting in a small relic abundance in the early Universe. The DM energy density at late times is given by ρ χ = m f χ n χ and the relic abundance is given by (9) is the Hubble expansion rate scale [4] and ρ c ≈ 1.05×10 −5 h 2 GeV/cm 3 is the critical cosmic energy density. With the rolling of φ after freeze out, raising the mass of χ to m f χ , we obtain the DM relic energy density we expect today. We will next estimate the typical parameters that could give rise to the above qualitatively described scenario.

RELIC ABUNDANCE IN THE EARLY UNIVERSE
In this work we consider only χ as the DM and not φ. However in principle φ could contribute to the DM budget in the Universe, depending on its mass and when it starts oscillating. As indicated in Eq. 1, the potential for φ is given by At early times when H > m φ , the φ field is overdamped and does not oscillate up until H m φ . Since we are interested in DM number densities that are far smaller than the usual thermal relic values, the freeze out parameter x f 20. Here x f ≡ m χ /T f , where T f is the freeze out temperature. Given that the large annihilation cross sections required will lead us to m χ GeV, we expect that T f is O(10 MeV), larger than ∼ MeV, to avoid disrupting BBN. Therefore, a typical expectation in our scenario is that the freeze out Hubble scale In what follows, we will consider final masses for DM up to m f χ ∼ 10 9 GeV. We find that pushing beyond this limit could require dealing with non-perturbative effects that could complicate the analysis and render our estimates unreliable. However, there should be no obstacle, in principle, to push beyond the above mass. We note that the final value of the scalar field φ 0 could be much larger than ∼ 10 9 GeV, as long as λ 1. Let us consider φ 0 ∼ 10 15 GeV, as may be expected in an ultraviolet (UV) framework, corresponding to λ ∼ 10 −6 . If m φ 10 −14 eV, then φ will start to roll after the freeze out era discussed earlier, corresponding to typical temperatures ∼ 10 MeV. Since we are assuming that the initial value of φ is given by φ i ≈ 0, the energy density initially stored in φ will be given by m 2 φ φ 2 0 10 20 eV 4 . Absent any significant interactions, this energy density would redshift by (T eq /T f ) 3 ∼ 10 −21 , where T eq ∼ 1 eV marks the era of matter-radiation equality, with energy density O(T 4 eq ). Hence, with the above sample values of parameters, it seems that if φ attains its ground state instantaneously, it does not have enough energy to account for the final energy stored in DM. Below, we will show that this apparent mismatch can be accounted for once the interactions of φ with the relic χ particles is properly considered 1 . We also note that the above initial energy density stored in φ is much less than the radiation energy density of ∼ (10 MeV) 4 , so the initial potential energy of φ does not lead to inflation.
For the above scenario to yield initially small χ mass, it is required that φ i ≈ 0 to begin with. One can just assume that this is the case, but it would be more interesting if the underlying physics could lead to such an initial condition. Here, we argue that this is indeed the case in our model (see also Ref. [18]). Let us declare upfront that we will ignore possible quantum corrections to the potential for φ. This question may need to be addressed in a fuller theory describing the above effective interactions and could require new dynamics beyond what has been considered here. While we do not claim a rigorous connection to the Higgs hierarchy problem in the SM, we simply note that the seeming insular nature of the SM could be hinting at a more subtle effect for quantum contributions to scalar masses. We do, however, consider the effect of the ambient medium on scalar masses.
For the assumed interactions, the thermal effects of the DM on the scalar φ can be parametrized by δm 2 φ ∼ λ 2 T 2 , [25] where T 10 MeV and λ ∼ 10 −6 . Hence, the mass of φ, considering thermal corrections will be given by where we have left out O(1) factors that would not affect our main conclusions. With the above contribution, the temperature dependent value of the scalar vev will be given by For typical values of parameters, m φ ∼ 10 −14 eV, λ ∼ 10 −6 , and T ∼ O(10 MeV), we see that φ(T ) ∼ 10 −30 φ 0 , which is a severe suppression and will ensure that the field is pegged near φ = 0 for T m i χ . 1 We thank Yue Zhao for discussions on this point.
The above formula is valid when the DM is relativistic and in thermal equilibrium with the SM. However, we would like to keep φ φ 0 until freeze-out. Since for the largest values of the final DM mass m f χ ∼ 10 9 GeV considered here we have m i χ /m f χ ∼ 10 −10 , we would like to have φ(T f )/φ 0 m i χ /m f χ , so that DM mass remains low enough throughout the freeze-out process. One can easily check that the discussion presented for this choice of parameters can be applied to other large values of m f χ considered in this work.
We note that m f χ ∼ 10 9 GeV for the THUMP mass corresponds to a relic number density of n χ ∼ 10 −18 T 3 f at freeze out, compared to ∼ T 3 in the relativistic regime. The in-medium "tadpole" term λφχχ → λφ n χ ∼ 10 −3 eV 3 φ would then dominate the vacuum term m 2 φ φ 0 φ ∼ (10 −4 eV 3 )φ. The in-medium "plasma frequency," given by ω 2 φ ∼ λ 2 n χ /m i χ ∼ 10 −17 eV 2 , is much larger than the vacuum parameter m 2 φ ∼ 10 −28 eV 2 . Given the dominance of the medium effects at T f ∼ 10 MeV, we note that the value of φ at this temperature is roughly given by χ /m f χ , which implies that the value of m χ is near its initial value all the way down to freeze out temperature and the relic number density calculations are governed by m i χ , as assumed in our treatment. The above discussion implicitly leads to the conclusion that due to the non-negligible coupling to DM, φ remains in thermal contact with the SM until temperatures near freeze-out. In our analysis the smallest freezeout temperature is ∼10 MeV below which φ contributes to radiation in a similar way as neutrinos. Below temperatures associated with electron-positron annihilation (T ∼ 1 MeV), the contribution of φ to the effective number of neutrino species is given by ∆N ef f = 4/7 ∼ 0.6. For larger freeze-out temperature, this number would be somewhat smaller and hence a typical prediction of our scenario is ∆N ef f 0.6. Interestingly, these values could potentially lead to easing of the tension in the Hubble parameter implied by local and early universe measurements [26][27][28][29][30].
The preceding analysis also implies that the value of φ remains frozen until the "tadpole" term becomes dominated by its vacuum value. For this to happen, n χ must be diluted by the expansion of the Universe corresponding to a temperature T * at which (T * /T f ) 3 ∼ 1/10, so that the vacuum tadpole starts to dominate. We note that this also suggests that the evolution of φ is not instantaneous and the aforementioned mismatch between the initial energy density m 2 φ φ 2 0 and the eventual DM energy content is remedied. This is achieved by the dilution of the relic density n χ for T ∼ T * at which φ starts to move away from φ(T f ). The above treatment only considers in-medium terms of linear and quadratic order in φ. This truncation is justified since φ starts out very small due to thermal effects prior to freeze out, as in-dicated by Eq. (4). Having illustrated that the mass of χ remains close to m i χ with our assumptions about the model parameters, we will next examine the annihilation processes that set the relic abundance of χ.
The DM particle χ could give the correct relic abundance through the annihilation process given in Fig.1  (plus u-channel). The annihilation cross-section can be expanded in powers of velocity and written as Using the results of Refs. [11,31,32], we find Here, r = m A /m i χ , and the a and b coefficients represent the s-wave and p-wave contributions respectively. For the values considered in this work, a bv 2 , making this annihilation process s-wave dominated. Expanding Eq. (6) around r and in the limit that r 1, we get Though this is a hidden sector annihilation, the dark sector and SM are in thermal equilibrium in the early Universe. DM χ freezes out when it is non-relativistic at a temperature corresponding to x = x f , which is given by (8) Here we take c ∼ 0.5, g χ counts the internal degrees of freedom of χ, which in our case is a Dirac fermion, and hence g χ = 4 [25]. The quantity g * counts the relativistic degrees of freedom and is evaluated at freeze-out. After freeze-out, the relic energy density of DM is given by where s 0 = 2891.2 cm −3 is the late time entropy density of the Universe, ρ c is the critical density defined above and m f χ is the final mass of the DM after φ → φ 0 . The quantity Y ∞ represents the late time comoving number density of DM which is given by The observed value of the DM energy density is given by Ω ob h 2 ≈ 0.12 [28]. For any set of input parameters {g D , m i χ , m A }, we can determine the a and b terms in Eq. (6), x f using Eq. (8), and Y ∞ in Eq. (10), which together yield Ω χ h 2 in Eq. (9).
On the left panel of Fig. 2, we plot the DM relic abundance in the m f χ vs m i χ parameter space and in the m i χ vs x f parameter space, on the right panel. We assume m i χ = 2m A (in general, it is sufficient to assume m i χ m A ), for two values of the dark sector coupling α D = g 2 D /4π. The solid (dashed) diagonal line corresponds to the observed relic abundance for α D = 1.0 (α D = 0.5). In this work, we limit the dark sector coupling to α D = 1.0 in order to avoid our theory being non-perturbative, as we discussed above. In order to explain the observed relic abundance, we show that the maximum m f χ we can attain in this model is ∼ 10 9 GeV (∼ 10 8 GeV) for α D = 1.0 (α D = 0.5). Furthermore, we note that if we were to choose α D ∼ 4π the maximum final mass we would obtain is ∼ 10 11 GeV, though we do not consider this estimate to be reliable. In Fig. 2, we see that near m i χ ∼ 10 GeV the relic abundance lines change shape; this is due to the fact g * (T ) in our scenario is calculated at freeze-out, and the change in shape corresponds to the QCD phase transition in the early Universe. DM could annihilate during BBN, thus disrupting elemental abundances during this era. To be conservative, we assume that DM freezes out before BBN. The shaded blue region shows this BBN bound, which is independent of m f χ , but dependent on α D and m i χ . To set this constraint, we assumed that T BBN ∼ 5 MeV [33] and used Eq. (8) to find the initial DM mass m i χ corresponding to freeze-out during BBN. For the benchmark values we have chosen here the boundary occurs near m i χ ∼ 200 MeV, as shown by the solid (dashed) line corresponding to α D = 1.0 (α D = 0.5). The figure on the right panel can be translated to obtain the temperature at which DM freezes-out for a given m i χ , since x f = m i χ /T f .

CONSTRAINTS AND PHENOMENOLOGY
In this section we include the various experimental bounds that apply to our model space. In our Galactic neighborhood the virial DM velocity is v ∼ 10 −3 . If a DM particle of mass m f χ scatters with a nucleon of mass m n ∼ 1 GeV, it will transfer momentum q ∼ m n v 1 MeV. Throughout this study, we assume m A ≤ 0.5 m i χ and in Fig. 2, the largest allowed value of m f χ ∼ 10 9 GeV corresponds to m i χ ∼ 300 MeV. Hence for m A ∼ 150 MeV, m A q and the DM scattering cross-section is given by [34]  where µ nχ is the reduced DM-nucleon mass. For the DM masses considered here, the reduced mass to very good approximation is given by µ nχ ≈ m n .
In Fig. 3 we show the current model constraints, in the ε vs m A parameter space, for benchmark α D = 1.0 and m f χ = 5 × 10 7 GeV. For the purposes of this study, we consider a mediator mass up to 10 GeV. Given our assumption that m i χ = 2m A , this corresponds to m f χ = 5×10 7 GeV, which we choose as our benchmark in Fig. 3. The grey shaded areas are current bounds on the visibly decaying dark photon from various fixed target and collider experiments 2 [36][37][38][39][40][41][42][43][44][45][46][47][48], and from Supernova 1987A [49][50][51]. The most stringent direct-detection bounds on weak scale DM come from the XENON1T experiment and for m χ ∼ 1 TeV, the limit lies at σ χn ∼ 10 −45 cm 2 [52]. The lack of a DM signal in XENON1T allows us to place a bound on our parameter space, which we obtain by rescaling the XENON1T limit with our mass of m f χ = 5 × 10 7 GeV. This gives us a cross-section upper bound of ∼ 5 × 10 −41 cm 2 , which is represented by the red line in Fig. 3. The red shaded region is ruled out by the XENON1T experiment. In a similar fashion, the blue dot-dashed line corresponds to σ χn ∼ 1.5 × 10 −43 cm 2 and is the projected sensitivity of the DARWIN experiment with a 200 ton-year exposure [53]. The orange dashed line is the boundary of the neutrino background, below which neutrinos start dominating the nuclear recoil spectrum [54]. 2 For the constraints not depicted here, see Ref. [35].
We also show projected sensitivity of future fixed target experiments. The region bounded by the purple solid line is the expected sensitivity reach for SHiP [48,55], while the green solid line bounds the expected reach for SeaQuest/DarkQuest with 10 20 protons on target [48,56,57]. For simplicity, here we only show projections from SeaQuest/DarkQuest and SHiP; for other experiments with similar reach, such as LongQuest or NA62, see Ref. [48]. To obtain the bounds shown in grey, we used the publicly available DarkCast code [58] from Ref. [35]. In Fig. 3, we show that future experiments such as SHiP or SeaQuest/DarkQuest would be able to probe the ∼ GeV-scale mediators of THUMP DM, corresponding to DM-nucleon cross-sections of 10 −45 − 10 −40 cm 2 . Some of this parameter space may be probed by next generation direct detection experiments, such as DAR-WIN. However the GeV-scale mediator signals available to accelerator experiments may provide hints of the dark sector below the neutrino floor.
We note that our model contains a long range force that acts on THUMPs with range typically of the order of 1/m φ ∼ 10 4 km. This could affect the distribution of DM on macroscopic scales with potentially observable effects. An examination of such effects is beyond the scope of this work. We also point out that a two-loop process mediated by DM and A couples φ to ordinary charged matter, however this effect is too suppressed to have any observable effects.

DISCUSSION AND CONCLUSIONS
In this work, we proposed a new mechanism for generating the thermal relic abundance of DM with masses above the so-called unitarity limit of ∼ 10 2 TeV. We showed that DM would be composed of GeV-scale WIMPs in the early Universe, before decoupling from the thermal bath. If these WIMPs are coupled significantly to other lighter dark sector mediators, they could annihilate efficiently enough that the DM relic abundance after freeze-out is too small, given its initial mass, to constitute the observed DM energy density today. We postulate that the DM is coupled to an ultra-light scalar φ which is held near φ = 0 in its potential before the DM freezes out. After freeze-out, φ starts rolling to some large minimum, sourcing a large mass for the DM, thereby raising its energy density to the observed levels. We call the DM in our scenario a THermal Ultra Massive Particle or "THUMP." In a specific model, we assumed that DM, which is initially at GeV scale can couple to a lighter dark photon mediator with strengths of up to α D = 1.0. The strength of this coupling determines how low the DM number density can be after freeze-out, setting the maximum final THUMP mass that would be sourced by the scalar φ in order to get the correct relic abundance.
This model provides a connection between ultra heavy DM and GeV-scale mediators. Hence, our proposal motivates DM-nucleon cross-sections that are within reach of the next generation direct detection experiments as well as GeV-scale mediators that can be probed at accelerators. Our scenario takes advantage of the various complementary methods of searching for DM, at both the cosmic and intensity frontiers. We further point out that the model discussed here may be extended to other examples such as a scalar mediator, which would give different phenomenology, however, we leave this for future work.