Spin-3/2 dark matter in a simple t-channel model

We consider a spin-3/2 fermionic dark matter (DM) particle interacting with the Standard Model quarks through the exchange of a charged and coloured scalar or vector mediator in a simple t-channel model. It is found that for the vector mediator case, almost the entire parameter space allowed by the observed relic density is already ruled out by the direct detection LUX data. No such bounds exist on the interaction mediated by scalar particles. Monojet + missing energy searches at the Large Hadron Collider provide the most stringent bounds on the parameters of the model for this case. The collider bounds put a lower limit on the allowed DM masses.


Introduction
Many astrophysical and cosmological observations during the last several decades provide strong evidence for the existence of dark matter (DM) in the Universe. The amount of DM has been precisely measured by the Planck satellite mission to be Ω DM h 2 = 0.1188± 0.0010 [1], where the cold dark matter (CDM) content is estimated to comprise roughly 26% of the total energy in the Universe. Investigations into the nature of dark matter particles and their interactions has emerged as an important field of research. Weakly interacting massive particle (WIMP) DM searches constitute an important programme at the Large Hadron Collider (LHC), where the ATLAS and CMS collaborations [2] are looking for DM signatures involving missing energy accompanied by a single or two jet events. It is expected, and there is indeed a real possibility, that the production of DM particles a e-mail: khogali11@gmail.com b e-mail: agoyal45@yahoo.com c e-mail: mukesh.kumar@cern.ch d e-mail: Alan.Cornell@wits.ac.za of any spin at 13 TeV centre-of-mass energy would be detected.
Null results from the direct detection experiments [3][4][5][6], which measure nuclear-recoil in DM-nucleon elastic scattering, have provided the most stringent upper bounds on the spin-independent DM-nucleon elastic scattering crosssection over a wide range of DM masses. This has provided important constraints on the DM models considered in the literature. In addition there are indirect detection experiments whose aim is to detect the signature of annihilating or decaying DM particles into the Standard Model (SM) particles.
Here we consider a spin-3/2 DM particle as an alternative to the conventional scalar, vector or spin-1/2 CDM particles. spin-3/2 CDM has been studied in EFT models and constraints from the relic density, direct and indirect observations obtained [7][8][9][10]. spin-3/2, 7.1 KeV warm dark matter (WDM) has been considered as a means to provide a viable explanation from the anomalous 3.1 KeV X-ray line observed by the XMM Newton [11]. Furthermore spin-3/2 DM with a Higgs portal has also been recently investigated [12]. This idea of using a spin-3/2 DM candidate was also considered in our recent paper [13], where we considered a minimal SM singlet spin-3/2 DM candidate interacting with a spin-1 mediator in a minimal flavour violation (MFV) schannel model. In the case of pure vector interactions, it was found that almost the entire parameter space was allowed by the observed relic density, and is ruled out by the direct detection observation data.
In this paper we consider a minimal spin-3/2 SM singlet DM candidate. The DM particle in this case is a t-channel annihilator. It interacts with the SM particles through the exchange of a scalar (S) or vector (V) particle. A class of such models for scalar and vector mediator couplings with a spin-1/2 DM candidate has been considered in Refs. [14][15][16][17]. The mediators in these t-channel models carry colour or leptonic index. As such we shall describe the model for this study in section 2, and in section 3 we discuss all relevant experimental constraints. The relic density contributed by the DM particles is calculated (taking into account the coannihilation processes), and assuming that the contribution by these spin-3/2 DM particles does not exceed the observed relic density, constraints on the parameters of the model are obtained in section 3.1. With these constraints in place we discuss the compatibility of these constraints from the direct and indirect detection experiments in section 3.2 and section 3.3 respectively. In section 4 we examine the signature of these DM particles at the LHC, where a monojet signal with missing energy is investigated. Section 5 is devoted to the summary and discussion of our main results.

The Model
The model consists of a SM singlet spin-3/2 particle interacting through the mediation of a scalar (S) or a vector (V µ ) which carries a baryonic (colour) or lepton index. In general the mediator couples to right-handed up-type quarks (or leptons), right-handed down-type quarks (or leptons) or left-handed quark (or lepton) doublets. We consider here the right-handed up-type quark case for simplicity, where the other cases are similar. The spin-3/2 free Lagrangian is given by [18]: with Note that χ µ satisfies Λ µν χ ν = 0, and with χ µ being on mass-shell we have The spin sum for spin-3/2 particles are and are given by [18]: In view of the non-renormalisable nature of interacting spin-3/2 theories, we can only write generic interactions which respect to the SM gauge symmetry between the singlet, χ, with SM fermions mediated by a scalar or a vector [19]. We will consider the vector and scalar mediator case separately: (a) Scalar mediator S: For the scalar mediator case we can write the SM gauge invariant interaction as: where i is a generation index and u i R ≡ (u R , c R , t R ). In this case we do not have a dimension-4 interaction term. This is because of the nature of the vector-spinor χ µ , which on the mass-shell satisfies γ µ χ µ = 0, and thus it is not possible to construct a Lorentz-invariant dimension-4 interaction term involving χ µ , S and the Dirac field u R .
(b) Vector mediator V µ : In this case we can write dimension-4 as well as dimension-5 interaction terms, namely and For all calculations we set Λ = 1 TeV. The interaction Lagrangian for the scalar and vector can be written as: The covariant derivative is given by where g s is the QCD strong coupling constant. Unlike the s-channel mediator, where a single mediator is required, in the t-channel model we require a different mediator for each generation.
In general, the interaction given in Lagrangian (7), (8) and (9) induce flavour-changing neutral currents (FCNC), which are strongly constrained by low energy phenomenology. The FCNC constraints can be avoided by imposing a minimal flavour violation (MFV) structure on the Yukawa couplings. The parameter space will consist of the DM candidates mass m χ , the vector (scalar) couplings and the mediator masses m i S (m i V ), for each generation. For simplicity we will set the couplings and mediator masses for all the generations to be equal. If the mediator mass in the kinematically accessible region of the LHC, the decay of the mediator and the ensuing signal will become important. The decay width of the scalar and vector mediators Γ (S i /V i → χū i ), dropping the generation index, are given by: since m i S , m χ m u is true for all quarks, except the top quark, and for dimension-4 and dimension-5 interaction Lagrangians given in Eqs. (8) and (9) respectively.

Constraints
In this section we examine the constraints on the model parameters m χ , m S , m V and the coupling constants from the relic density, direct and indirect observations.

Relic density
In the early Universe the DM relic density is determined by the dominant DM annihilation processes χχ → uū, mediated by the t-channel exchange of scalar/vector mediators.
Since the mediators in this model carry colour and charge, co-annihilation processes like χ S(V ) → u g and S S * (V V * ) → g g (even though exponentially suppressed when mass splitting m S/V − m χ > freeze-out temperature T f ), will play an important role if the DM mass gets close to the mediator mass. The co-annihilation processes χ S(V ) → u g are mediated by t-channel exchanges of mediators, as well as by schannel exchanges of gluons and through the four-point interaction involving the DM, mediator, u-quark and the gluon vertex. These processes will reduce the Yukawa coupling needed to generate the required thermal relic abundance. Self annihilation mediator processes S S * (V V * ) → g g are generated by purely gauge interactions, and are independent of the Yukawa couplings, having the potential to suppress the relic density below the observed value. At freeze-out the DM and mediator particles are nonrelativistic. In the non-degenerate parameter space the channel χχ S/V − − → uū cross-section can be easily evaluated, and in the limit m χ , m S , m V m u are given by and for the scalar-mediator and the vector-mediator dimension-4 and dimension-5 interaction Lagrangians (7), (8) and (9) respectively, with the mass ratio r = m S (m V )/m χ . The thermal relic density of χ's is obtained by solving the Boltzmann equation:  Fig. 1 Contour plots in the allowed DM mass m χ and the mass splitting ratio r − 1 (with r = m S (m V )/m χ ) in the left panels, and in the DM mass m χ and the couplings in the right panels. We have assumed that the DM χ saturates the observed relic density. The top panels are for the dimension-4 interaction term for the vector mediator case. The middle and the bottom panels are for dimension-5 vector and scalar mediator cases respectively. In the left panels the colour gradient corresponds to the Yukawa couplings required to give the observed relic density, which in the right panels, the colour gradient corresponds to the mass splitting consistent with the observed relic density.
where H is the Hubble constant, σ |v is the thermally averaged χ-annihilation cross-section, and η eq χ = 4 The freeze out occurs when the χ's are non-relativistic, and then σ |v can be written as: The Boltzmann equation can be solved numerically to yield [20]: where g * is the number of degrees of freedom at freeze-out temperature T F , X F = m χ /T F is obtained by solving where c is taken to be 1/2. g * (X F ) is the number of degrees of freedom at the freeze-out temperature, and is taken to be 92 for our estimate, and g = 4.
The annihilation cross-section for the co-annihilation processes χ S(V ) → u g in this limit are given by and σ (χ V → u g)|v g V To calculate the relic density we have implemented the t-channel scalar and vector interactions with SM quarks and spin-3/2 DM, including the relevant co-annihilation processes, in micrOMEGAS [21]. Note that this numerically solves the Boltzmann equation by taking the full expressions of the annihilation cross-section. We have checked the relic abundance in the non-degenerate parameter space for some representative values of the parameters, and found them to be in agreement with the numerical calculations done by micrOMEGAS. The necessary model files for micrOMEGAS were built using FeynRules [22]. In Figure 1 we show the contour graphs in the DM mass and the mass splitting ratio r − 1 in the left panels, and between m χ and couplings in the right panels. The colour gradients correspond to the Yukawa couplings in the right panels, and to the mass splitting ratio in the left panels, to confirm to the observed relic density Ω DM h 2 0.12. In the parameter space in which co-annihilation is not important, comparatively large Yukawa couplings are required to obtain the required relic density. In the co-annihilation region on the other hand, we find the couplings to be reduced for almost all DM masses both for the scalar and vector mediator cases. We find that with the increase in DM masses, the co-annihilation channels take over the DM self annihilation processes and the co-annihilation channels involving gauge interactions alone are able to depress the relic density below the observed value. We see from the left-hand panels that there are two regions in the DM mass, one around 80 GeV ≤ m χ ≤ 100 GeV and another one around 300 GeV ≤ m χ ≤ 400 GeV, where the co-annihilation processes result in a sharp drop in the couplings required for the requisite relic density.

Direct detection
Direct detection experiments [3][4][5][6] on elastic nucleon-DM scattering have provided the most stringent bounds on DM mass and interactions in a large number of conventional DM models. In the t-channel spin-3/2 DM model considered here, the cross-sections at zero momentum transfer can be easily calculated [23][24][25]. The dominant contribution to the spinindependent cross-section for the vector mediator case is given by and where µ = m χ m N m χ + m N , f N = 4 for protons and 1 for the neutrons, and we have dropped the terms proportional to the quark mass and momenta in comparison to the leading term. The cross-section given in Eqs. (27) and (28) correspond to the dimension-4 and dimension-5 interaction Lagrangians. The elastic nucleon-DM cross section for the case of scalar mediator is suppressed by terms proportional to quark momenta and has not been considered here.
In Figure 2 we show the predictions for the spin-independent DM-proton scattering cross-sections, σ SI , for the vector mediator case. In the left panels the colour gradient corresponds to the coupling, and in right panels to the mass splitting r − 1. In the left panels we observe that for every DM mass and mass splitting the Yukawa coupling is obtained such that the parameters conform to the observed relic density, whereas in the right panels the required mass splitting is obtained for a given Yukawa coupling. We find that for any DM mass the scattering cross-section generally increases as the degenerate parameter region is approached. This happens because of a resonant enhancement of σ SI near r = 1.  Fig. 2 The spin-independent proton-DM cross-section σ SI . The top and the bottom panels correspond to dimension-4 and dimension-5 vector interactions. In the left and right panel the colour gradients correspond to the Yukawa couplings and mass splittings respectively. All parameters are consistent with the observed relic density. We have also shown the graphs from the observed current upper limits from LUX [6] and PANDAX-II [5] experiments. The projected upper limit for XENON1T [4] has also been shown. Almost the entire parameter space (m χ , m V ) for the vector mediator case considered here is already ruled out from the LUX data.
For the case of dimension-5 vector interactions (bottom panels), we see a drop of several orders of magnitude in the scattering cross-section around the same DM mass regions where the co-annihilation results have a sharp drop in the couplings. In Figure 2 we have also shown the current upper limits from LUX [6], PandaX-II [5] and the projected upper limit for the XENON1T experiment [4].

Indirect detection
The Fermi Large Area Telescope (LAT) collaborations [26] have dedicated detectors to measure cosmic ray fluxes arising from DM annihilation in the Universe. In Figure 3 we show the prediction for the total DM annihilation into uū for the vector/scalar mediated t-channel model. The predictions shown here are for the DM mass, mediator mass and the couplings consistent with the observed relic density. We have also shown the bounds from the 95% CL upper limits on the thermally-averaged cross-section for DM particles annihilating into uū Fermi-LAT observations. As expected in the parameter region where co-annihilation is important (r 1) the χχ annihilation cross-section in the uū channel is greatly suppressed. Even in the region away from resonance (r 1) the Fermi-LAT data does not provide strong bounds on the mass and coupling parameters in the entire range consistent with Ω DM h 2 = 0.12.

Collider bounds
The t-channel mediator model considered here has scalar and vector mediators which carry colour, SU(2) L and U(1) charges. They can thus be singly produced in association with DM particles, or pair produced if they are light enough at the LHC. These processes will contribute to the monojet and dijet signals with missing energy, with distinct signatures that can be searched for in dedicated searches. For the monojet events q g → q χχ are the dominant processes, in comparison to qq → g χχ, because of the large parton distribution probability of the gluon, as compared to quark and antiquark in the proton. The authors of the simplified DM  Fig. 3 The prediction for the DM χχ annihilation rate into uū, as a function of the DM mass m χ . All the parameters are chosen to be consistent with the observed relic density. The top and the middle panels are for dimension-4 and dimension-5 vector interactions respectively. The bottom panels are for the scalar interaction. The colour gradient in the left and right panels correspond to the coupling and mass splitting respectively. Bounds from the Fermi-LAT experiments are also shown [26]. model document [27] have emphasised that the dominance of the associated production channels is a distinct feature of t-channel models. The 8 TeV CMS collaboration data based on an integrated luminosity 19.7 fb −1 [2,28] has been used by the authors of Ref. [16,29] to put bounds on the coupling of fermionic DM as a function of the mediator and DM mass for the case of scalar and vector mediators. In the present study we confine ourselves to constraints arising from the monojet signals using the parameters space (m χ m S/V ) for different values of the couplings g S consistent with the observed relic density. The cross-section for monojet events is obtained by generating parton level events for the process p p → χχ j using MadGraph [30], where the model file for the Lagrangian is obtained from FeynRules, and we use CTEQ611 parton distribution function for five  Figure 4 we find that the collider bounds are much weaker compared to the bounds from the direct detection experiments for the vector mediator case. The scalar mediator case is interesting in this case as the collider bound rules out low mass DM particles.The bounds from the monojet + missing energy cross section puts a lower limit on the DM particle mass, where the limit depends on the coupling, and increases with the coupling.

Summary and discussion
In this paper we have considered a spin-3/2 DM particle interacting with the SM fermions through the exchange of a scalar or a vector mediator in the t-channel. Invoking MFV we restricted ourselves to the coupling of DM candidates with SM singlet right-handed quarks with universal coupling. The thermal relic DM abundance has been computed by taking into account the co-annihilation processes. Coannihilation has the effect of reducing the Yukawa couplings needed to generate the required DM density. The co-annihilation effects are more pronounced in the large m χ regime, where mediator self annihilation into gauge bosons has the potential to suppress the relic density below the observed value. Similar behaviour was observed in t-channel model for spin-1/2 [31] and scalar DM [32] particles. Our main observations are the following: (a) The direct detection experiments, through DM-nucleon elastic scattering data, provide the most stringent bounds for the case of a vector mediator. In this case the entire parameter space allowed by the relic density is already ruled out by the LUX data. This result is consistent with the earlier studies of spin-3/2 DM in the EFT [7] frame work for pure vector couplings, as well as in a simplified s-channel model [13]. The co-annihilation is unable to ameliorate this situation. (b) There are no strong bounds from the the direct detection experiments on the scalar mediated interactions due to the velocity suppression of σ SI . In contrast, in the EFT frame work, both the scalar as well as vector interactions give rise to dominant spin-independent nucleon-DM scattering cross-section and direct detection rules out scalar interaction for spin-3/2 DM particles of mass lying between 10 GeV and 1 TeV [7]. (c) The current constraints from indirect searches like, Fermi-LAT data, are not sensitive enough to put any meaningful constraints. (d) Monojet searches at the LHC do not provide strong bounds at the vector couplings in comparison to the bounds from direct detection experiments. However, for the case of the scalar mediator, where we do not get any strong bounds from the direct detection experiment, collider bounds put a lower limit on the DM mass which is m χ ≥ 300 GeV. This limit rises with the increase in coupling.
Finally, it may be mentioned that bounds from direct detection experiments can, however, be evaded by foregoing the universal coupling between DM mediators and quarks, and letting the DM particles interact with only one generation, say with the third generation quarks (top-philic DM).