A complete effective field theory for dark matter

JHEP07(2021)081 Published for SISSA by Springer Received:May 13, 2021 Revised:June 15, 2021 Accepted:June 22, 2021 Published:July 14, 2021 A complete effective field theory for dark matter Juan Carlos Criado,aAbdelhak Djouadi,b,c Manuel Pérez-Victoriaband José Santiagob aInstitute for Particle Physics Phenomenology, Department of Physics, Durham University, South Rd, Durham DH1 3LE, U.K. bFacultad de Ciencias (Edificio Mecenas), Campus Fuentenueva, Granada E-18007, Spain cLaboratory of High Energy and Computational Physics, NICPB, Rävala pst. 10, Tallinn 10143, Estonia E-mail: [email protected],[email protected],[email protected], [email protected] Abstract: We present an effective field theory describing the relevant interactions of the Standard Model with an electrically neutral particle that can account for the dark matter in the Universe. The possible mediators of these interactions are assumed to be heavy. The dark matter candidates that we consider have spin 0, 1/2 or 1, belong to an electroweak multiplet with arbitrary isospin and hypercharge and their stability at cosmological scales is guaranteed by imposing a Z2symmetry. We present the most general framework for describing the interaction of the dark matter with standard particles, and construct a general non-redundant basis of the gauge-invariant operators up to dimension six. The basis includes multiplets with non-vanishing hypercharge, which can also be viable DM candidates. We give two examples illustrating the phenomenological use of such a general effective framework. First, we consider the case of a scalar singlet, provide convenient semi-analytical expressions for the relevant dark matter observables, use present experimental data to set constraints on the Wilson coefficients of the operators, and show how the interplay of different operators can open new allowed windows in the parameter space of the model. Then we study the case of a lepton isodoublet, which involves coannihilation processes, and we discuss the impact of the operators on the particle mass splitting and direct detection cross sections. These examples highlight the importance of the contribution of the various non-renormalizable operators, which can even dominate over the gauge interactions in certain cases. Keywords: Beyond Standard Model, Cosmology of Theories beyond the SM, Effective Field Theories ArXiv ePrint: 2104.14443 Open Access,c The Authors. Article funded by SCOAP3.https://doi.org/10.1007/JHEP07(2021)081 JHEP07(2021)081 Contents 1 Introduction 1 2 Effective field theory for a generic DM multiplet 5 2.1 Symmetries and field content 5 2.2 Suppressing gauge couplings 7 3 Scalar singlet dark matter 9 4 Lepton doublet dark matter 15 5 Conclusions 21 A Operator basis 23 A.1 Operators for a scalar multiplet ϕ24 A.2 Operators for a fermion multiplet χ26 A.3 Operators for a vector multiplet ρ27 B Basis of ϕ†ϕ φ†φD230 1 Introduction The existence of dark matter (DM) in the universe is strongly suggested by several astrophysical and cosmological measurements but its very nature remains enigmatic. Particle physics proposes a plausible and effective solution to this problem in terms of an electrically neutral and weakly interacting massive particle that is stable at cosmological scales [1,2]. DM particles are predicted by many extensions of the Standard Model (SM), including the well motivated ones that address other important theoretical or experimental issues of the model such as supersymmetric [3,4] or extra-dimensional models [5–8].1In fact, basically any extension of the SM with additional neutral particles can accommodate a DM state, provided that a discrete symmetry is imposed to protect it from decaying into SM particles. These DM particles are turning out to be the new holy Grail of contemporary physics, actively searched in many astroparticle [9–14] and collider [15–19] experiments. Because of the large proliferation of DM candidates, it has become customary and quite useful to consider effective field theory (EFT) approaches, which allow to study in a model-independent manner the phenomenology of these particles. It is typically assumed that the new state is either a scalar, a vector or a fermion, although higher spins have 1In these extensions the dark matter candidate is the lightest particle of the new sector, which can be made stable by virtue of a discrete symmetry. This symmetry also helps in relaxing direct and indirect collider limits. – 1 – JHEP07(2021)081 been also considered, see refs. [20,21] for recent accounts. Among the simplest and most economical of these EFTs are the Higgs-portal models, in which a single DM particle is introduced in addition to the SM particles; it interacts in pairs only with the Higgs sector of the theory, which is assumed to be minimal [22] and, hence, involves only the unique Higgs boson observed at the LHC [23,24]; see refs. [25,26] for recent reviews. There are also other possibilities for these simple EFTs and, for instance, Z-portal [27–29] or neutrinoportal [30–32] scenarios have been also extensively discussed. These models are rather predictive as all DM observables can be specified by only a few basic new parameters, for example, the mass of the DM particle and its coupling with the mediator particle. However, it might well be that the relevant new physics extension has also additional new particles, heavy scalars, vectors or fermions that could accompany the DM and/or act as DM mediators. Actually, the simple portal models, with the mediator being exclusively a SM particle, are highly constrained by present data. In this case, another possible approach which also introduces only a handful of free parameters has been put forward: simplified models [33,34] in which an effective Lagrangian is introduced that explicitly includes the mediator particle and its interactions with the SM and DM states. For instance, a model with a singlet DM fermion and a singlet scalar mediator has been studied in ref. [35] and the importance of imposing the full SM gauge symmetry is emphasized; other examples of simplified models can be found in refs. [36–41]. This hybrid approach has many advantages, in particular in the context of collider physics when the mass of the DM and its possible companions are comparable to the collider energy, but in practice introduces, unlike the simple EFT approach with only DM fields, some additional model dependence that makes it rather close to the concrete realizations that it is meant to simplify. We will thus not discuss this approach further here. Instead, in this paper, we follow a genuine EFT approach, including all the operators allowed by the symmetries. We assume that all possible mediators of the DM interactions are heavy and have been integrated out. The degrees of freedom in the EFT are only the SM particles and the DM particle, together with its possible gauge partners, which will be nearly-degenerate in mass. Indeed, whenever the DM field has non-vanishing isospin, it must appear together with other fields, forming complete SU(2)Lmultiplets. In our EFT we consider a DM particle of spin 0, 1 2or 1, with well-defined gauge quantum numbers. That is, the DM field appears in the neutral component of an SU(2)L×U(1)Ymultiplet. For these three different spin assignments and for arbitrary isospin and hypercharge, we consider the most general EFT that describes this scenario, with the additional assumption of a Z2symmetry to stabilize DM. This is nothing but an extension of the SMEFT [42] with new degrees of freedom: an additional scalar, vector or fermionic multiplet containing the DM state.2We perform a systematic classification of the gauge-invariant operators that appear in this extended EFT, and construct a general basis for all relevant operators up to dimension six. These operators capture the effects of the mediators or additional heavy particles, and also of possible non-perturbative UV completions [44]. The Wilson 2Due to the Z2symmetry, the EFT in this paper is complementary to the one considered in ref. [43], in which the extra fields in the SM extension were required to have linear interactions with SM operators. – 2 – JHEP07(2021)081 coefficients of these operators, to be treated as free parameters, can be constrained using present data on the cosmological relic density, direct and indirect DM detection in astroparticle physics experiments and in missing energy searches or invisible Higgs boson decays at high-energy colliders like the LHC. DM multiplets with non-vanishing hypercharge are often discarded because their scattering off nucleons have large cross-sections mediated by a Zboson, in conflict with directdetection bounds. However, it is known that the dominant contribution to the cross-section is avoided if the DM particle is described by a Majorana spinor [45]. This is the case of fermionic DM particles with non-vanishing hypercharge whenever the two components of the Dirac field are non-degenerate. The mass splitting can be produced by mixing with an additional Majorana fermion. More relevant to our EFT setup with one multiplet, it has also been shown in the case of a fermion doublet that it can be induced by nonrenormalizable operators [46]. Here, we argue that this mechanism is very general and works for arbitrary spin whenever there is some breaking of a global U(1) symmetry acting on the DM multiplets. We distinguish the operators in our EFT that break this symmetry and observe that, to dimension 6, they require fields of isospin 1/2. There have been many analyses in the past that have considered EFTs whose field content were that of the SM extended with DM states. In the low-energy regime, nonrelativistic EFTs have been discussed, for example in refs. [47–57], but more focus has been put on the relativistic case. For instance, a complete set of operators of dimension ≤6(written in the broken phase of the electroweak sector) for an EFT containing the SM fields together with a complex scalar DM field in a singlet, doublet or triplet representation is given in ref. [58]. In ref. [59], an EFT consisting of the SM particles together with Majorana or real scalar DM has been considered and the operators of dimension ≤8have been classified in the case where the DM particle is assumed to be a singlet under the SM gauge group and coupling only to fermions. In ref. [60], a basis of operators of dimension ≤6for an EFT with the SM and DM particles has been given, with the non-SM fields being a light right-handed neutrino as well as a singlet scalar, fermion and vector DM fields which were made stable by invoking a Z2symmetry under which the SM fields are uncharged and the DM fields are charged. A basis of dimension-6 operators describing interactions of a singlet-like Majorana DM fermion with SM particles, and its applications for collider and astroparticle searches, has been introduced in refs. [61–63]. Closer to our scope in this paper, a minimal basis of operators of dimension 6 or less describing the interaction of a singlet scalar, a Dirac fermion and a vector DM with quarks and gluons only has been given in ref. [64]. Finally, in ref. [65], a basis of effective operators up to dimension 7 has been given in the case of scalar and fermionic DM embedded in a general SM multiplet. Subsets of effective operators relevant for different applications have been considered in refs. [39,66–75]. In particular, DM interactions with quarks and gluons and their impact on direct detection as well as searches at colliders have been discussed in refs. [76–78] in the case of singlet Dirac and/or Majorana DM and in ref. [78] in the case of scalar DM. The specific case of fermion DM couplings to photons was considered in ref. [79]. An EFT for fermionic and scalar DM interactions with quarks, gluons and photons at low energy scales – 3 – JHEP07(2021)081 has been presented in ref. [80]. In ref. [81], loop effects in an EFT of the SM extended with a singlet Dirac fermion DM particle have been studied. The matching at the electroweak scale between the DM EFTs with and without the standard fermion, gauge and Higgs bosons has been performed in ref. [82]. Finally, the important aspect of co-annihilation has been also discussed and, for instance, an EFT for the DM in this context was constructed in ref. [83] (while a systematic classification of simplified models for co-annihilation has been given in ref. [84]). In this paper, we complete the previous analyses by presenting the most general framework for the description of the interactions between one multiplet that contains the DM particle and SM particles, which we assume to be the only relevant degrees of freedom at sufficiently low energies. The EFT for the DM-SM system in our setting is obtained by extending the SMEFT field content and symmetries, which besides the Lorentz group and the SM gauge invariance include a discrete Z2symmetry in order to stabilize the DM particle. As mentioned before, we study the case of spin-0, 1 2and 1 DM, not constrained to be a SM singlet. In particular, we extend in several ways on the work of ref. [65], which to the best of our knowledge contained the most complete list of operators for generic DM available to this date. First, we discuss the case of vector DM multiplets in addition to the scalar and fermionic cases. Second, we do not assume the dark U(1)Dsymmetry considered in ref. [65], which allows for extra operators for multiplets with hypercharge Y= 1/2.3Third, our non-redundant basis for scalar multiplets with arbitrary hypercharge Y6= 1/2includes two operators that were not considered in ref. [65]. Finally, we present some interesting phenomenological implications of our approach and, in particular, we address the important consequences of the additional operators for multiplets with hypercharge Y= 1/2. As we explain below, these operators are required to produce the mass splitting that helps in avoiding direct-detection limits. On the other hand, we stop at dimension 6 and do not consider dimension-7 operators, which have been included in ref. [65]. The phenomenology of singlet and non-singlet DM particles is quite different and we give an example of each possibility to illustrate the use of our effective theory and to show the relevance of the non-renormalizable operators. We first consider the case of a singlet scalar DM state, where we give semi-analytical expressions for all the relevant DM observables, which are rather simple and convenient for practical use. Indeed, these expressions allow a more efficient exploration of the effect of the various operators and, as an example, we use them to study the possible cancellation of the contributions of two different operators, which allow to ease the stringent experimental constraint from direct DM detection. In a subsequent step, we illustrate the important impact of co-annihilation by considering the case of an iso-doublet of heavy leptons and discuss the simultaneous effect of the various dimension-5 operators on the mass splitting between the neutral DM state and the other charged and neutral leptons, as well as on the annihilation and coannihilation cross sections into SM particles. The paper is organized as follows. In the next section, we discuss the derivation of the complete basis of operators of dimension ≤6which involve two scalar, fermionic or 3We thank the authors of ref. [65] for a discussion on this issue. – 4 – JHEP07(2021)081 vector DM multiplets and SM fermions, gauge and Higgs bosons. In sections 3and 4, we use our EFT to analyze the phenomenology of a scalar singlet and of a fermion doublet, respectively. A short conclusion is given in section 5. In appendix A, we present our complete basis of operators up to dimension 6 and in appendix Bwe comment on the basis of operators of the form ϕ†ϕ φ†φD2. 2 Effective field theory for a generic DM multiplet 2.1 Symmetries and field content Our aim in this study is to set up a general framework for the description of the interactions between the DM and the SM particles in the case in which the possible mediators of the interactions are heavier than these particles and than the energies which are expected to be probed. This DM-EFT can be obtained by extending the SMEFT (the effective field theory constructed with the SM fields only) with an extra field content, a multiplet that we will generically denote X, that involves the DM particle, X0, and it possible companions. The SM gauge group GSM = SU(3)C×SU(2)L×U(1)Ymust be contained in the symmetry group of this EFT, so we shall impose the SM gauge symmetry without losing generality. Our only assumption in this regard is that GSM is linearly realized (as it is the case in the SMEFT) but, of course, we shall also impose Lorentz symmetry. Thus, the extra fields of the DM-EFT can be organized into irreducible representations of GSM and the Lorentz group. As any EFT, our theory is valid only below a cut-off Λ, which should be larger than the electroweak scale v≈246 GeV as well as the DM mass MX. When introducing the various operators that describe the interactions in this EFT, we implement the usual power counting assigning a typical value of Λ4−∆to the Wilson coefficient of an operator of dimension ∆. In order to work with a manageable theory some restrictions on the DM sector need to be imposed. In this paper we make the following assumptions: 1. In order to stabilize the DM particle, we impose a discrete Z2symmetry. 2. The field content of the theory is given by the SM one, including the Higgs doublet φ, and a single extra multiplet Xthat belongs to an irreducible representation of GSM. Under the Lorentz group, Xtransforms either as a scalar, a spinor or a vector. All SM fields are even under Z2, while Xis odd. Due to the restriction to one extra multiplet, our EFT can be regarded as a minimal extension of the SMEFT, much as the minimal dark matter scenario [85–87] is a minimal extension of the SM. We make this assumption for simplicity but, in fact, the same EFT works in the case of several flavors of the same multiplet and also when there is a separate Z2 for each type of multiplet; only an additional labelling of the different fields and couplings would be required. The general case with different types of multiplets would instead involve extra operators, which we do not write here. If the cutoff Λis large, the assumption of a Z2symmetry is not necessary for multiplets in sufficiently large representations of GSM, as the gauge symmetry then forbids operators – 5 – JHEP07(2021)081 linear in Xbelow a given dimension.4The DM particle is then not absolutely stable and could decay, but the corresponding suppression of the decay width could make the lifetime longer than the age of the Universe, as emphasized in ref. [85]. However, we are mostly interested in the case in which Λis not much larger than the TeV scale, so the suppression will not be sufficient except for extremely large representations. Finally, the restriction to spin up to unity is also made for simplicity; DM with higher spins has been recently considered in refs. [20,21]. In order to include a good DM candidate, Xmust be a color singlet and must contain an electrically-neutral component X0. This component is the field associated to the DM particle(s). The existence of X0restricts the possible electroweak quantum numbers of X. Let Tbe its isospin and Yits hypercharge. Then, an electrically neutral component is present if and only if the difference T− |Y|is a natural number. We assume that Yis non-negative without loss of generality: from a multiplet X, one can obtain a new one e X in the same representation of SU(2)Land with opposite hypercharge as e X=Eˆ X,(2.1) where Eis the anti-diagonal matrix with entries (1,−1,1,−1, . . .)and ˆ Xis the conjugate of the field X(ˆϕ=ϕ∗,ˆχ=γ0Cχ∗and ˆ Vµ=V∗ µfor scalar, fermion and vector fields, respectively).5We thus have an infinite but discrete set of possibilities for the DM multiplet X, given by T∈ {0,1/2,1, . . .}, Y ∈ {T, T −1, . . . , T −bTc},(2.2) where bTcdenotes the integer part of the isospin T. The minimal field content of multiplets with vanishing hypercharge and integer isospin corresponds to irreducible representations carried by real scalars, Majorana spinors or real vectors. We note that multiplets with large Tvalues will accelerate the running of the SU(2) gauge coupling g(and the U(1) coupling g0if the hypercharge Yis also large). We do not impose any restriction in this sense because Λis arbitrary and the scale at which g becomes non-perturbative will typically be beyond the regime of validity of the EFT. Our approach covers a large number of the models, both complete in the ultraviolet regime and constructions in the effective theory approach, that have been proposed in the literature and involving a weakly interacting massive DM particle. Some of these scenarios have been mentioned before and, for instance, models in which the DM interacts with the SM particles through the Higgs portal have been reviewed recently in ref. [25] where a rather exhaustive list of references may be found. For most phenomenological purposes, the only relevant operators are those that contain the DM field quadratically. The kinetic and mass terms read L0=ηϕDµϕ†Dµϕ−M2 ϕϕ†ϕ,(2.3) 4For instance, note that no renormalizable operator linear in Xis allowed for representations not included in the list which is given in e.g. ref. [43]. 5This is related to charge conjugation Cby CX C−1=ηC X ˆ X,ηC X=±1, so Aand Bare charge conjugation eigenstates with eigenvalue ηC A,B. – 6 – JHEP07(2021)081 for a scalar multiplet, denoted by ϕ; L0=ηχ¯χi / Dχ −Mχ¯χχ,(2.4) for a fermion multiplet, denoted by χ; and L0=ηρDµρν†Dνρµ−Dµρν†Dµρν+M2 ρρ† µρµ,(2.5) for a vector multiplet, denoted by ρ. The derivatives Dµare covariant with respect to GSM and ηX= 1,1 2for complex and real representations of the multiplet X, respectively. Note that the spin-one extra particles are described in our effective treatment by a Proca vector field. This vector field can represent an extra gauge field of an extended gauge group;6in our minimal setup, this comes along with the assumption that the Higgs fields associated to the spontaneous breaking of the additional gauge invariance are heavier than the scale Λ, and have been integrated out. In order to allow for non-perturbative UV completions, we do not impose here the restrictions on the Wilson coefficients that would arise from the extra gauge invariance. A list of all the operators with dimension D≤6for any electroweak multiplet Xis given in appendix A. This is the main result of this paper. The effective Lagrangian can thus be written as LDMEFT =LSMEFT +L0+X i cH iOH i+X icNH iONH i+h.c.,(2.6) where LSMEFT is the SMEFT Lagrangian; the OH iand ONH ioperators are the hermitian and non-hermitian operators in appendix Arespectively. Also, in the equation above, the cH icoefficients are real, while the cNH iare complex. Not all structures are allowed for all choices of isospin and hypercharge. The restrictions arise because the quantum numbers of the product of two Xmust match the quantum numbers of the combination of SM fields appearing in a given operator. For instance, when D≤6, whenever the latter has non-vanishing hypercharge, Xmust have Y= 1/2for the operator to be gauge invariant. These tables comprise a non-redundant basis of operators: gauge-invariant operators not written here can be expressed as linear combinations of these by the use of algebraic identities, integration by parts and field redefinitions. 2.2 Suppressing gauge couplings The dark matter fields in multiplets with Y= 0 obviously have no coupling to the Zboson. On the other hand, multiplets with Y6= 0 are usually avoided because in general gauge interactions, and in particular the coupling to the Zboson, give too large a contribution to the direct detection cross section, well above the current experimental limits. However, as we discuss in this section, there is a quite generic scenario in which the vector coupling to the Zboson, and therefore its contribution to spin-independent direct detection, vanishes. This leaves only gauge contributions to spin-dependent direct detection processes that are much less constrained experimentally. 6This is actually necessary for perturbativity if the cutoff is not close to the mass of the spin 1 particle. – 7 – JHEP07(2021)081 The electrically neutral component of multiples with Y6= 0 is a complex field X0, which transforms trivially under U(1)Qbut non-trivially under U(1)Y. Hence, it has gauge couplings to the Z boson field with strength eY/(sWcW). This complex field is made of two self-conjugate components: X0=A+iB, with ˆ A=A,ˆ B=B. Consider now the mass term of the fields Aand B. It is a hermitian form over self-conjugate fields, so it is equivalent to a quadratic form and the 2×2mass matrix Mis real and symmetric. This is diagonalized by a real orthogonal transformation and its eigenfields N1and N2are also self-conjugate. If Mhas two degenerate eigenvalues, then Aand Bare mass eigenfields as well, and the full complex X0constitutes the DM candidate. On the other hand, if the eigenvalues are non-degenerate, then only the particle associated to the lighter mass eigenfield, N1, survives and forms the DM today. In this case, only N1is relevant as an initial state in DM-nucleon interactions. Furthermore, if the mass separation is larger than the typical kinetic energy of a DM particle near the Earth, only N1will appear as a final state. Let us consider the non-degenerate case. We are interested here in the trilinear interaction of two N1fields and one Z boson field, which will be of the form jN1 µZµ. The current jN1 µis necessarily invariant under C, since it is bilinear in N1and this field is self-conjugate. Hence, since Aµis Codd, these trilinear interactions cannot conserve C. If the interaction preserves CP, it must be parity-odd. That is, jN1 µmust be a pseudovector current. Vector currents that couple to the Z boson preserving CP(as the ones from minimal coupling) cannot be built with N1alone. These are the currents that enter the spin-independent nuclear-DM cross-sections mediated by a Z, which thus vanish for non-degenerate Ni. In this scenario, small contributions from non-renormalizable operators can dominate direct detection. The relevant question is then: when are the two mass eigenfields non-degenerate? Consider a global U(1)Xsymmetry transformation acting only on the DM multiplet X. This acts as an SO(2) transformation on the two-dimensional real vector space generated by Aand B. In fact, Aand Bspan the vector space of a two-dimensional irreducible representation of this SO(2), since no one-dimensional subspace is left invariant. Indeed, a generic SO(2) transformation transforms any vector into a linearly independent one. The vector space of the irreducible representation is thus the complete space generated by A and B. If the global U(1)Xis a symmetry of the action, then the mass matrix Mwill commute with the SO(2) transformations and all the linear combinations of Aand Bwill be mass eigenfields with the same eigenvalue.7So, if U(1)Xis preserved, we are in the degenerate case. Conversely, any breaking of U(1)Xwill generically produce a splitting of the eigenmasses. This is the case, for instance, of SUSY DM models in which the Higgsino mixes with a Majorana gaugino, which breaks the symmetry. In our effective treatment with only one DM multiplet, the breaking of U(1)Xcan only occur in operators in which the combination of DM fields has non-vanishing hypercharge. The reason is that U(1)Xacts 7An equivalent way to put this is that a 2×2matrix left invariant by arbitrary rotations must be a multiple of the identity. – 8 – JHEP07(2021)081 be trusted away from SM particle production thresholds Mϕ≃mSM (with mSM the mass of any SM particle) and from the Higgs boson resonance Mϕ≃1 2mH. The DM-nucleon scattering elastic cross section is simpler to parametrize. We find the approximate formula σN|s=(Mχ+mN)2= (10−46 cm2) 100 GeV mN+Mχ!2X ij Nij ˜ci˜cj,(3.4) where the dimensionless coefficients Nij are ciand Mχindependent. The sum runs over the φ1,φ4,uφ1and dφ1Wilson coefficients. Again, we obtain the values of Nij from the numerical calculation, and display them in table 2. Similarly, the formula for the branching ratio of the Higgs boson decay into invisible DM particle is given by BR(H→inv) = θ(mH−2Mχ)q1−4M2 ϕ/m2 HPij Hijcicj 1 + q1−4M2 ϕ/m2 HPij Hijcicj ,(3.5) where the Hij coefficients are ciand Mχindependent. The sum runs over the φ1,φ4and φWilson coefficients. The values of the Hij obtained from the numerical calculation are also given in table 2. These expressions greatly facilitate the exploration of the parameter space of the Wilson coefficient. For example, examining the table of Nij coefficients suggests that a cancellation between the contributions of the operators with cφ1and cdφ1coefficients to the DM-nucleon cross section can occur. This opens the possibility of evading the stringent limits from direct detection. We show in figure 3the limits in the [cφ1, cdφ1]plane for two values of the DM particle mass, Mϕ= 50 and 200 GeV, computed numerically using micrOMEGAs. As can be seen there, the semi-analytical treatment has led us to a viable possibility: while masses below the TeV scale are excluded by direct detection limits when only the dimension-4 portal is present, the addition of a dimension-6 operator allows to avoid them while keeping the correct relic DM abundance. For Mϕ>1 2mH, no other limits apply. For Mϕ<1 2mH, the tiny set of values not excluded by direct detection is only ruled out by the limits on the invisible Higgs branching ratio BR(H→inv). Similar results as those shown in figure 3, with in particular a suppression of the direct detection limits, have been obtained in a different context in ref. [41]. 4 Lepton doublet dark matter In this section, we consider the possibility in which the DM particle belongs to a non-singlet SU(2)L×U(1)Ymultiplet and illustrate the use of our approach in the specific case of a DM lepton doublet.10 We start by discussing the contributions of the various operators to the DM annihilation cross sections and to the mass splitting. 10There is a vast literature on DM scenarios with singlet-doublet leptons and, in the context of Higgs portal models, they have been reviewed in e.g. ref. [25]. Doublet leptons alone have been discussed in – 15 – JHEP07(2021)081 Figure 3. Line of points (in black) that give the correct DM abundance, together with the region that gives an overabundance of DM (in gray), the region excluded by direct-detection experiments (in blue) and the region excluded by the upper bound on BR(H→inv)(in orange), in the space of cdφ1vs. cφ1, for different values of Mϕ. Non-singlet multiplets are fundamentally different from singlets in two ways: they have dimension-4 couplings fixed by gauge invariance and more than one particle is present in the spectrum near the DM mass. In order to illustrate this possibility, we consider the case in which the multiplet is a spin-1/2 doublet ∆=(N, E)Twith hypercharge Y=−1/2, so that e ∆corresponds to the Y= 1/2>0multiplet χwe considered in section 2. We choose the opposite hypercharge in order for ∆to have the same quantum numbers as the usual SM lepton doublet. The leading higher-dimensional operators for ∆have dimension 5 and are listed in table 3. Apart from the vector-like mass term in eq. (2.4), the particles contained in the ∆ multiplet receive tree-level corrections to their mass from the four dimension-5 Oφi operators. After the spontaneous breaking of the electroweak symmetry, the mass Lagrangian takes the form −Lmass =Nc LNR v2cφ31 2M−v2 2(cφ1−cφ2) 1 2M−v2 2(cφ1−cφ2)v2cφ4  NL Nc R! + M−v2 2(cφ1+cφ2)!ELER+h.c. .(4.1) The diagonalization of the mass matrix for the states NLand Nc Ris performed by means the context of supersymmetric models with higgsino DM, see for instance refs. [105–107], and in the nonsupersymmetric case, for instance in refs. [46,108]. Non-singlet scalar multiplets have been explored in the context of composite Higgs models in [109,110]. – 16 – JHEP07(2021)081 Operators Constraints Oφ1,Oφ2,Oφ3,Oφ4H→inv, direct detection OB,OW— Table 3. Dimension-5 operators for a DM lepton doublet. All of them contribute to the DM annihilation cross section. of the unitary transformation NL Nc R!= cos θ−eiφ sin θ e−iφ sin θcos θ! e−iα 0 0e−iβ! N1 N2!,(4.2) with mixing angle θgiven by sin2θ=1 2+|cφ3|2−|cφ4|2 2|cφ3|2−|cφ4|22+|cφ3+cφ4| M 2v2−cφ1−cφ2 4−1 .(4.3) After diagonalization, we find that our multiplet contains three particles: two heavy Majorana neutrinos N1and N2, and a heavy charged lepton Ewith masses given by MN1=M"1 + v2 2M−Re cφ1+ Re cφ2−2|cφ3+c∗ φ4|+Oc2 i#,(4.4) MN2=M"1 + v2 2M−Re cφ1+ Re cφ2+ 2|cφ3+c∗ φ4|+Oc2 i#,(4.5) ME=M"1−v2 2M(Re cφ1+ Re cφ2) + Oc2 i#.(4.6) In order for the neutral state N1to be stable, its mass must be smaller than the charged lepton Emass: MN1< ME. This occurs if and only if Re cφ2<|cφ3+c∗ φ4|.(4.7) When both cφ3and cφ4vanish, the masses M1and M2become degenerate, and it is then convenient to keep the neutrino mass matrix anti-diagonal, and view NLand NRand as the leftand right-handed components of a unique Dirac neutrino. The coefficients cφ2,cφ3and cφ4thus provide a general order-1/Λ2parametrization of new physics effects in the mass splitting. Even when no new physics is present, there will be a remaining splitting induced by loops involving the gauge interactions [85–87]. We will assume that the tree-level effects of the dimension-5 coefficients are dominant here, and neglect the loop-induced splitting. Let us briefly summarize the present experimental constraints on the masses and couplings of the vector-like leptons. First, because of the Z2symmetry which has been introduced to stabilize the DM particle, these heavy leptons will not mix with the SM ones and, thus, will not decay into SM particles and cannot be produced in association with them. Hence, the usual stringent constraints on heavy leptons from the anomalous magnetic moments of the muon or electron [111] or from the associated production with charged SM – 17 – JHEP07(2021)081 leptons and neutrinos at LEP2 [112] will not hold.11 In the case of the DM particles, the only relevant constraint will come from the invisible decays of the Zboson into DM pairs (when the particle has non-zero hypercharge) which sets a bound of MN1&45 GeV in the Dirac and MN1&39.5GeV in the Majorana cases [112]; there are also bounds from the invisible decay of the Higgs boson, MN1.1 2mH= 62.5GeV [17,18] if the Higgs-N1N1 coupling is not too small. In the case of the charged lepton, there is a strict lower bound on its mass from the pair production at LEP2: ME&103 GeV for an almost stable E and ME&101 GeV if it decays into a light neutrino and a Wboson [112]. This bounds translates to a bound on MN1when the leptons are close in mass. Searches of these particles have also been conducted at the LHC where they can be produced in the Drell-Yan type processes p→E+E−, NN, E±N(Nstands for a generic neutral lepton) with the latter mode being by far dominant as it is mediated by the exchange of the charged Wboson [113]. However, the constraints are again tight only in the case where the heavy leptons mix with the ordinary ones, leading to prompt electrons or muons and W, Z bosons in the final state [114,115]. In our case, the lightest Nparticle is stable while the next-to-lightest N0and the charged lepton Estates will decay into the DM and Zor Wbosons. If the mass difference between the DM particle and its companions is very small, as it is necessary the case for the co-annihilation mechanism to be relevant, the intermediate gauge bosons would be far off-shell and the available phase space would be tiny, so the processes would correspond to the production of long-lived particles that lead to displaced vertices or particles being eventually stable at the detector level. There are searches by ATLAS and CMS of stable and long-lived charged sleptons and charginos, which constrain these particles to have masses above a few hundred GeV, see e.g. refs. [116,117], and could be relevant in our case. However, as these bounds are model dependent and are very sensitive to the details of the analyses, we will not include them in our study. Dedicated experiments such as MoEDAL, FASER or SHiP for instance could probe longlived particles more efficiently, see e.g. ref. [118] for a review. In our analysis, we will therefore simply use the constraints ME&100 GeV for the charged leptons and MN&45 GeV for the neutral ones, which have a non-zero hypercharge. We will of course also take into account the bounds from the Higgs invisible branching ratios as measured at the LHC, BR(H→NN).20% [17,18] as well as constraints from direct and indirect detection of the DM particle. According to eq. (4.7), there are two classes of Wilson coefficients, those that induce a splitting between the masses MEand MN1, and those that do not, which we denote respectively by cMand c/ M: cM∈ {cφ2, cφ3, cφ4}, c / M∈ {cφ1, cB, cW}.(4.8) An analysis similar to the one in section 3, of the effects of each coefficient individually while turning off all others, is only possible for the cMcoefficients. The relative effects of gauge interactions and dimension-5 operators in the DM annihilation cross section are controlled by two parameters: the relative mass splitting of the 11Some of the dimension-6 operators in table 6would contribute to the lepton anomalous magnetic moments; however, they are not relevant for our discussion here. – 18 – JHEP07(2021)081 charged and lightest neutral leptons ∆≡(ME−MN1)/MN1∼cMv2/MN1, and the effective Yukawa coupling of the neutral lepton YN1∼cMv. Then, there are three different cases to consider: i) ∆ .0.1: in which the co-annihilation processes dominate in the cross sections. In terms of cMand MN1this happens when cM.(0.02 TeV−1)M 10 GeV.(4.9) ii) ∆ &0.1and YN1.gEW: in which DM annihilation processes through the s-channel exchange of the Z-boson dominate. This occurs is when (0.02 TeV−1)M 10 GeV.cM.3 TeV−1.(4.10) iii)YN1&gEW: the annihilation processes through the s-channel exchange of the Higgs boson (with dim-5 operators) dominate, when 3 TeV−1.cM.(4.11) There are two issues with case iii): first, the masses MN1that correspond to the correct abundance are low, about 10 GeV for M, and thus excluded by LEP searches; and second, the perturbative expansion of the EFT is broken if dim-5 interactions have a similar size to dim-4 ones. cMbeing high also means that the masses MN1and MEshould be small, MN1, ME.Λ∼c−1 M. We conclude that in almost all of the allowed parameter space, dimension-5 operators affect the abundance only through the splitting. Not all observables depend on the cMcoefficients through the splitting. The case of direct detection, for instance, can be strikingly different since, as we have discussed in detail in section 2.2, the couplings to the Zboson leading to spin-independent direct detection can be naturally suppressed in this case [119]. Indeed, when N1is a Majorana fermion, its vector current vanishes: N1γµN1=Nc 1γµNc 1=NT 1C†γµCNT 1=−N1γµN1.(4.12) This means that its coupling to the Zboson must be through the axial vector current: LNZ =e cWsW ¯ NγµNZµ ⊃e cWsW ¯ N1(eiα cos θPR−e−i(φ+α)sin θPL)γµ(e−iα cos θPL−ei(φ+α)sin θPR)N1Zµ =−igZ¯ N1γµγ5N1Zµ,(4.13) where we have used eq. (4.2) in passing from the first to the second line, eq. (4.12) in going from the second to the third one, sWand cWare the sine and cosine of the Weinberg angle, and gZ=ecos(2θ) 2sWcW .(4.14) – 19 – JHEP07(2021)081 Thus, this interaction vanishes for the value θ=π/4of the mixing angle, which occurs if cφ3=cφ4. In the case in which cφ3>0,cφ4= 0 and M/cφi v2, the mixing angle and the N1N1Zcoupling are simply given by sin2θ≃1 2+v2|cφ3| M, gZ≃ev2|cφ3| cWsWM≃0.04 cφ3 0.1 TeV−1100 GeV M.(4.15) The coupling of the N1states to the Higgs boson would be given, in this case, by LN1N1H=YN1N1N1H, YN1=√2vcφ3≃0.03 cφ3 0.1 TeV−1.(4.16) Since the N1coupling to the Zboson only contributes to the spin-dependent DM-nucleon cross section while the coupling to the Hboson contributes to the spin-independent one, and the limits on the spin-dependent cross section are at least five orders of magnitude weaker [12], the direct detection limits are essentially controlled by the interaction with the Higgs particle. In figure 4, we show the limits in various regions of parameter space from the relevant experimental constraints: the DM abundance, the direct detection cross section and the searches for heavy leptons at LEP. We also include in each case a red area in which the EFT description is expected to break down, which is the region where ME>Λ = c−1 for each Wilson coefficient c. The two plots at the top correspond to turning on one cMtype coefficient at a time: cφ3and cφ2. We choose negative values for cφ2, as required by eq. (4.7) when cφ3=cφ4= 0. On the bottom plots we turn on more than one coefficient at a time. On the left, we display the limits in [MN1, cφ3]plane when the relation cφ4=cφ3is imposed. On the right, we consider the case in which cφ4=cφ3= 0.1 TeV−1, and explore the [MN1, cB]space. The generation of the correct abundance is dominated by co-annihilations for the black curve on the right of each plot in figure 4. This corresponds to the regime labeled i)above. The curve on the left is dominated by annihilations through the Z, corresponding to case ii). The limit from LEP becomes just the upper bound on MN1for large cM, since this implies a large splitting; and just the upper bound on ME, for small cM. In the three plots in which cφ36= 0,N1is a Majorana particle, and the N1N1Zcoupling does not contribute to the spin-independent DM-nucleon cross section. Thus, the region excluded by direct detection is controlled by the N1N1Hcoupling, which is proportional to the relevant cM coefficient. On the top-right plot, the N1,N2pair is degenerate and can be seen as a single Dirac fermion. Then, there is an spin-independent contribution to the DM-nucleon cross section from the N1N1Zcoupling, leading to strong limits from direct detection. Since the cφ3, cφ4→0limit should lead to the Dirac case, one may wonder how the strong limits from direct detection arise in that case. In this limit, the mass splitting between N1and N2becomes small. Then, the collision between N1and a nucleon may produce an N2particle through the exchange of a Zboson with vector (not axial vector) couplings, thus contributing to the spin-independent cross section. However, in order for this to happen, the mass splitting must be smaller than the typical kinetic energy of a DM particle near the Earth, which is Mv2/2, with v≃1×10−3. This does not happen in the region we consider. – 20 – JHEP07(2021)081 101102103104105 MN1(GeV) 10−2 10−1 100 cφ3(TeV−1) 101102103104105 MN1(GeV) 10−2 10−1 100 −cφ2(TeV−1) 101102103104105 MN1(GeV) 10−2 10−1 100 cφ4(TeV−1) cφ4=cφ3 101102103104105 MN1(GeV) 10−2 10−1 100 cB(TeV−1) cφ3=cφ4= 0.1TeV−1 over abundant direct detection LEP MN1>Λ ME< MN1 Figure 4. Constraints on the coefficients of operators for the DM lepton doublet. The points over the solid black line give the correct DM abundance, while the gray region gives an overabundance of DM. The orange region is constrained by LEP searches of heavy leptons. The region excluded by direct-detection experiments is shown in blue. Indirect detection places not relevant limits. On the bottom-right plot of figure 4, we see that most of the limits on [MN1, cB] space do not depend on cB. This is because its effects are subdominant in each case: the annihilation cross section on the left is dominated by the dimension-4 coupling to the Z, the mass is controlled by Mand cM, and the contribution of cBto direct detection is suppressed by the DM velocity while that of the cMis not. However, for high-enough masses, cBhas a relevant contribution to the annihilation cross section, since then the effective coupling McBbecomes comparable to the dimension-4 gauge couplings and larger than the effective vcMcoupling from the other coefficients. 5 Conclusions Effective field theory has proven to be a useful approach for model independent studies of the phenomenology of the weakly interacting, massive and cosmologically stable parti- – 21 – JHEP07(2021)081 cles that are expected to form the dark matter in the universe. In this paper, we have investigated the possibility of going beyond the simple EFTs which have been discussed in the past, in which the interaction of the DM particles with the SM ones is parameterized in terms of a single or a few dominant operators, renormalisable or not, with the mediators of the interactions assumed to be very heavy and integrated out. Assuming that the DM particles have spin-0, 1 2or 1, and appear as the neutral components of a single SU(2)L×U(1)Ymultiplet with arbitrary isospin and hypercharge, we construct a general and non-redundant basis for all relevant gauge-invariant operators in this theory up to dimension six. This complete EFT can be viewed as being simply the SMEFT that has been widely discussed in recent years, in which an additional scalar, vector or fermionic multiplet containing the DM particle is added to the SM spectrum. We have then illustrated the usefulness of such a general approach with two specific examples. We have first considered a singlet scalar DM state and constrained the Wilson coefficients of all the relevant operators using present data on the DM: its cosmological relic abundance, direct and indirect detection in astroparticle experiments and searches in collider experiments, like through the invisible decays of the SM Higgs boson at the LHC. In order to simplify such an analysis, we have proposed a rather simple and convenient set of semi-analytical expressions for these DM observables, which allows to explore the complex parameter space of the various operators in a very efficient manner. Using our expressions it is, for instance, trivial to find possible cancellations in the contribution to very constrained observables, like direct DM detection, thus opening up the allowed parameter space. In a second illustration, we have considered the example of a vector-like lepton isodoublet, and we have discussed the interplay between the gauge interactions that are present in this case and dimension-5 operators. In particular, we have studied the simultaneous impact of some the higher dimension operators which enter both the mass splitting between the different members of the multiplet and the annihilation and co-annihilation cross sections of these states into SM particles. We have also described how, when the DM particle is Majorana, its gauge coupling to the Zboson leads to a vanishing contribution to the very constrained spin-independent DM-nucleon cross-section relevant for direct detection experiments. In this case the direct detection constraints are dominated by higher-dimensional operators despite the fact that the DM particle has non-zero hypercharge. We have shown that this is a quite generic mechanism, which can be applied to scalar, fermion and vector DM particles with non-vanishing hypercharge, equal to Y= 1/2if we restrict ourselves to operators of mass dimension up to six, as we have done in this work. Acknowledgments Discussions with Giorgio Arcadi, Mikael Chala, Guilherme Guedes, Ennio Salvioni and Jure Zupan are gratefully acknowledged. J.C.C. is supported by the STFC under grant ST/P001246/1. The work of A.D. is supported by the Junta de Andalucia through the Talentia Senior program and, in part, by the ERC Mobilitas Plus grant MOBTT86. This work has been partially supported by the Ministerio de Ciencia e Innovación project PID2019- – 22 – JHEP07(2021)081 106087GB-C22 and by Junta de Andalucia projects FQM-101, A-FQM-211-UGR18, P18FR-4314 and SOMM17/6104/UGR (including ERDF). A Operator basis In this appendix, we present the basis of operators of dimension D≤6containing two multiplets X, for the EFT defined in section 2. The operators for scalar multiplet ϕare shown in tables 4and 5; for a fermion χ, they are in table 6; and for a vector ρ, they are in tables 7,8and 9. We do not list the complex conjugate of non-hermitian operators. We have checked the counting of operators with each field content for the lowest isopins T= 0,1/2,...,4using the code BasisGen [120]. In all the tables use the following notation: •The SU(2) index of X=ϕ, χ, ρ, any SU(2) doublet index and Lorentz spinor indices are implicit. •Lower case letters a, b, . . . are used for SU(2) triplet indices, capital letters A, B, . . . are used for SU(3) octet indices and the letters I, J, . . . are used for SU(2) quadruplet indices. •Tadenotes the SU(2) generators in the Xrepresentation, is the Levi-Civita symbol, CIab denotes the quadruplet-triplet-triplet SU(2) Clebsh-Gordan coefficients. When the SU(2) isospin Tis larger than 1/2, we denote by Qathe unique set of square matrices acting on Xsuch that X†QaXtransforms as a quadruplet. •The tilde symboleis used to denote both the operation defined in eq. (2.1), when used over the DM multiplet and the dual e Fµν =1 2µνρσFρσ when applied to a field-strength tensor F. •GA µν,Wa µν and Bµν are, respectively, the SU(3)C,SU(2)Land U(1)Yfield strengths. The SM matter content is denoted by φfor the Higgs doublet, qand lfor the lefthanded quark and lepton doublets and u,dand efor the charge 2/3 and −1/3quarks and charged lepton singlets. •We define the standard linear combinations of covariant derivatives ↔ Dµ=Dµ−← Dµfor both any multiplet, ↔ Da µ=TaDµ−← DµTawhen applied to Xand ↔ Da µ=σaDµ−← Dµσa when applied to φ. The allowed operators for real representations are those for which both the “SU(2) irrep” and the “Hypercharge” columns of tables 4–8read “any”. In these tables, the operators are listed assuming a complex multiplet. The real case is recovered by identifying X=X†. For spinors this means χL=χc R, so that some operators are duplicated and the copies are to be discarded. Moreover, for vector fields the operators containing ρ† µρν contracted with an anti-symmetric tensor σµν,Fµν or e Fµν will vanish. – 23 – JHEP07(2021)081 A.1 Operators for a scalar multiplet ϕ Name Operator SU(2) irrep Hypercharge Dimension Oφ1(ϕ†ϕ)(φ†φ)any any 4 Oφ2(ϕ†Taϕ)(φ†σaφ)T > 0any 4 Oφ3(e ϕ†Taϕ)(φ†σae φ)T∈Z+ 1/2 1/24 Oφ4(ϕ†ϕ)(φ†φ)2any any 6 Oφ5(ϕ†Taϕ)(φ†φ)(φ†σaφ)T > 0any 6 Oφ6iabc(ϕ†Taϕ)(φ†σbφ)(φ†σcφ)T > 1/2any 6 Oφ7(e ϕ†Taϕ)(φ†σae φ)(φ†φ)T > 0 1/26 Oφ1(ϕ†ϕ)(φ†φ)any any 6 Oφ2(ϕ†Taϕ)D2(φ†σaφ)T > 0any 6 OφD1(ϕ†↔ Dµϕ)(φ†↔ Dµφ)any any 6 OφD2(ϕ†↔ Da µϕ)(φ†↔ Daµφ)T > 0any 6 OφD3(e ϕ†↔ Dµϕ)(e φ†↔ Dµφ)T∈Z+ 1/2 1/26 OφD4(e ϕ†↔ Da µϕ)(e φ†↔ Daµφ)T∈Z+ 1/2 1/26 OB(ϕ†ϕ)(BµνBµν)any any 6 Oe B(ϕ†ϕ)(Bµν e Bµν)any any 6 OW1(ϕ†ϕ)(Wa µνWaµν)any any 6 Oe W1(ϕ†ϕ)(Wa µν f Waµν)any any 6 OW2iabc(ϕ†Taϕ)(Wb µνWcµν)T > 1/2any 6 Oe W2iabc(ϕ†Taϕ)(Wb µν f Wcµν)T > 1/2any 6 OBW (ϕ†Taϕ)(BµνWaµν)T > 0any 6 OBe W(ϕ†Taϕ)(Bµν f Waµν)T > 0any 6 OG(ϕ†ϕ)(GA µνGAµν)any any 6 Oe G(ϕ†ϕ)(GA µν e GAµν)any any 6 Table 4. Basis of operators of dimension ≤6with two scalar DM multiplets and SM bosons only. – 24 – JHEP07(2021)081 [3] J.R. Ellis, J.S. Hagelin, D.V. Nanopoulos, K.A. Olive and M. Srednicki, Supersymmetric Relics from the Big Bang,Nucl. Phys. B 238 (1984) 453 [INSPIRE]. [4] H. Goldberg, Constraint on the Photino Mass from Cosmology,Phys. Rev. Lett. 50 (1983) 1419 [Erratum ibid. 103 (2009) 099905] [INSPIRE]. [5] G. Servant and T.M.P. Tait, Is the lightest Kaluza-Klein particle a viable dark matter candidate?,Nucl. Phys. B 650 (2003) 391 [hep-ph/0206071] [INSPIRE]. [6] H.-C. Cheng, J.L. Feng and K.T. Matchev, Kaluza-Klein dark matter,Phys. Rev. Lett. 89 (2002) 211301 [hep-ph/0207125] [INSPIRE]. [7] K. Agashe, A. Falkowski, I. Low and G. Servant, KK Parity in Warped Extra Dimension, JHEP 04 (2008) 027 [arXiv:0712.2455] [INSPIRE]. [8] G. Panico, E. Ponton, J. Santiago and M. Serone, Dark Matter and Electroweak Symmetry Breaking in Models with Warped Extra Dimensions,Phys. Rev. D 77 (2008) 115012 [arXiv:0801.1645] [INSPIRE]. [9] Planck collaboration, Planck 2018 results. VI. Cosmological parameters,Astron. Astrophys. 641 (2020) A6 [arXiv:1807.06209] [INSPIRE]. [10] XENON collaboration, Dark Matter Search Results from a One Ton-Year Exposure of XENON1T,Phys. Rev. Lett. 121 (2018) 111302 [arXiv:1805.12562] [INSPIRE]. [11] XENON collaboration, Light Dark Matter Search with Ionization Signals in XENON1T, Phys. Rev. Lett. 123 (2019) 251801 [arXiv:1907.11485] [INSPIRE]. [12] XENON collaboration, Constraining the spin-dependent WIMP-nucleon cross sections with XENON1T,Phys. Rev. Lett. 122 (2019) 141301 [arXiv:1902.03234] [INSPIRE]. [13] DarkSide collaboration, Low-Mass Dark Matter Search with the DarkSide-50 Experiment, Phys. Rev. Lett. 121 (2018) 081307 [arXiv:1802.06994] [INSPIRE]. [14] DARWIN collaboration, DARWIN: towards the ultimate dark matter detector,JCAP 11 (2016) 017 [arXiv:1606.07001] [INSPIRE]. [15] ATLAS collaboration, Constraints on new phenomena via Higgs boson couplings and invisible decays with the ATLAS detector,JHEP 11 (2015) 206 [arXiv:1509.00672] [INSPIRE]. [16] CMS collaboration, Searches for invisible decays of the Higgs boson in pp collisions at √s= 7,8, and 13 TeV,JHEP 02 (2017) 135 [arXiv:1610.09218] [INSPIRE]. [17] ATLAS collaboration, Combination of searches for invisible Higgs boson decays with the ATLAS experiment,Phys. Rev. Lett. 122 (2019) 231801 [arXiv:1904.05105] [INSPIRE]. [18] CMS collaboration, Search for invisible decays of a Higgs boson produced through vector boson fusion in proton-proton collisions at √s= 13 TeV,Phys. Lett. B 793 (2019) 520 [arXiv:1809.05937] [INSPIRE]. [19] CMS collaboration, Search for dark matter, extra dimensions, and unparticles in monojet events in proton-proton collisions at √s= 8 TeV,Eur. Phys. J. C 75 (2015) 235 [arXiv:1408.3583] [INSPIRE]. [20] J.C. Criado, N. Koivunen, M. Raidal and H. Veermäe, Dark matter of any spin — an effective field theory and applications,Phys. Rev. D 102 (2020) 125031 [arXiv:2010.02224] [INSPIRE]. – 31 – JHEP07(2021)081 [21] A. Falkowski, G. Isabella and C.S. Machado, On-shell effective theory for higher-spin dark matter,SciPost Phys. 10 (2021) 101 [arXiv:2011.05339] [INSPIRE]. [22] A. Djouadi, The Anatomy of electro-weak symmetry breaking. I: The Higgs boson in the standard model,Phys. Rept. 457 (2008) 1 [hep-ph/0503172] [INSPIRE]. [23] ATLAS collaboration, Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC,Phys. Lett. B 716 (2012) 1 [arXiv:1207.7214] [INSPIRE]. [24] CMS collaboration, Observation of a New Boson at a Mass of 125 GeV with the CMS Experiment at the LHC,Phys. Lett. B 716 (2012) 30 [arXiv:1207.7235] [INSPIRE]. [25] G. Arcadi, A. Djouadi and M. Raidal, Dark Matter through the Higgs portal,Phys. Rept. 842 (2020) 1 [arXiv:1903.03616] [INSPIRE]. [26] G. Arcadi, A. Djouadi and M. Kado, The Higgs-portal for Dark Matter: effective field theories versus concrete realizations,arXiv:2101.02507 [INSPIRE]. [27] R.C. Cotta, J.L. Hewett, M.P. Le and T.G. Rizzo, Bounds on Dark Matter Interactions with Electroweak Gauge Bosons,Phys. Rev. D 88 (2013) 116009 [arXiv:1210.0525] [INSPIRE]. [28] G. Arcadi, Y. Mambrini and F. Richard, Z-portal dark matter,JCAP 03 (2015) 018 [arXiv:1411.2985] [INSPIRE]. [29] J.R. Ellis, A. Fowlie, L. Marzola and M. Raidal, Statistical Analyses of Higgsand Z-Portal Dark Matter Models,Phys. Rev. D 97 (2018) 115014 [arXiv:1711.09912] [INSPIRE]. [30] V. González-Macías, J.I. Illana and J. Wudka, A realistic model for Dark Matter interactions in the neutrino portal paradigm,JHEP 05 (2016) 171 [arXiv:1601.05051] [INSPIRE]. [31] M. Escudero, N. Rius and V. Sanz, Sterile neutrino portal to Dark Matter I: The U(1)B−L case,JHEP 02 (2017) 045 [arXiv:1606.01258] [INSPIRE]. [32] B. Batell, T. Han, D. McKeen and B. Shams Es Haghi, Thermal Dark Matter Through the Dirac Neutrino Portal,Phys. Rev. D 97 (2018) 075016 [arXiv:1709.07001] [INSPIRE]. [33] LHC New Physics Working Group, Simplified Models for LHC New Physics Searches,J. Phys. G 39 (2012) 105005 [arXiv:1105.2838] [INSPIRE]. [34] J. Abdallah et al., Simplified Models for Dark Matter and Missing Energy Searches at the LHC,arXiv:1409.2893 [INSPIRE]. [35] S. Baek, P. Ko, M. Park, W.-I. Park and C. Yu, Beyond the Dark matter effective field theory and a simplified model approach at colliders,Phys. Lett. B 756 (2016) 289 [arXiv:1506.06556] [INSPIRE]. [36] N.F. Bell, Y. Cai and R.K. Leane, Mono-W Dark Matter Signals at the LHC: Simplified Model Analysis,JCAP 01 (2016) 051 [arXiv:1512.00476] [INSPIRE]. [37] F.-Y. Cyr-Racine, K. Sigurdson, J. Zavala, T. Bringmann, M. Vogelsberger and C. Pfrommer, ETHOS — an effective theory of structure formation: From dark particle physics to the matter distribution of the Universe,Phys. Rev. D 93 (2016) 123527 [arXiv:1512.05344] [INSPIRE]. [38] D. Goncalves, P.A.N. Machado and J.M. No, Simplified Models for Dark Matter Face their Consistent Completions,Phys. Rev. D 95 (2017) 055027 [arXiv:1611.04593] [INSPIRE]. – 32 – JHEP07(2021)081 [39] A. De Simone and T. Jacques, Simplified models vs. effective field theory approaches in dark matter searches,Eur. Phys. J. C 76 (2016) 367 [arXiv:1603.08002] [INSPIRE]. [40] T. Alanne and F. Goertz, Extended Dark Matter EFT,Eur. Phys. J. C 80 (2020) 446 [arXiv:1712.07626] [INSPIRE]. [41] T. Alanne, G. Arcadi, F. Goertz, V. Tenorth and S. Vogl, Model-independent constraints with extended dark matter EFT,JHEP 10 (2020) 172 [arXiv:2006.07174] [INSPIRE]. [42] I. Brivio and M. Trott, The Standard Model as an Effective Field Theory,Phys. Rept. 793 (2019) 1 [arXiv:1706.08945] [INSPIRE]. [43] J. de Blas, J.C. Criado, M. Pérez-Victoria and J. Santiago, Effective description of general extensions of the Standard Model: the complete tree-level dictionary,JHEP 03 (2018) 109 [arXiv:1711.10391] [INSPIRE]. [44] S. Bruggisser, F. Riva and A. Urbano, Strongly Interacting Light Dark Matter,SciPost Phys. 3(2017) 017 [arXiv:1607.02474] [INSPIRE]. [45] B. Kayser and R.E. Shrock, Distinguishing Between Dirac and Majorana Neutrinos in Neutral Current Reactions,Phys. Lett. B 112 (1982) 137 [INSPIRE]. [46] A. Dedes, D. Karamitros and V.C. Spanos, Effective Theory for Electroweak Doublet Dark Matter,Phys. Rev. D 94 (2016) 095008 [arXiv:1607.05040] [INSPIRE]. [47] J. Fan, M. Reece and L.-T. Wang, Non-relativistic effective theory of dark matter direct detection,JCAP 11 (2010) 042 [arXiv:1008.1591] [INSPIRE]. [48] A.L. Fitzpatrick, W. Haxton, E. Katz, N. Lubbers and Y. Xu, Model Independent Direct Detection Analyses,arXiv:1211.2818 [INSPIRE]. [49] A.L. Fitzpatrick, W. Haxton, E. Katz, N. Lubbers and Y. Xu, The Effective Field Theory of Dark Matter Direct Detection,JCAP 02 (2013) 004 [arXiv:1203.3542] [INSPIRE]. [50] B. Bellazzini, M. Cliche and P. Tanedo, Effective theory of self-interacting dark matter, Phys. Rev. D 88 (2013) 083506 [arXiv:1307.1129] [INSPIRE]. [51] M. Cirelli, E. Del Nobile and P. Panci, Tools for model-independent bounds in direct dark matter searches,JCAP 10 (2013) 019 [arXiv:1307.5955] [INSPIRE]. [52] R. Catena and P. Gondolo, Global fits of the dark matter-nucleon effective interactions, JCAP 09 (2014) 045 [arXiv:1405.2637] [INSPIRE]. [53] R. Catena, Prospects for direct detection of dark matter in an effective theory approach, JCAP 07 (2014) 055 [arXiv:1406.0524] [INSPIRE]. [54] G. Ovanesyan, T.R. Slatyer and I.W. Stewart, Heavy Dark Matter Annihilation from Effective Field Theory,Phys. Rev. Lett. 114 (2015) 211302 [arXiv:1409.8294] [INSPIRE]. [55] SuperCDMS collaboration, Dark matter effective field theory scattering in direct detection experiments,Phys. Rev. D 91 (2015) 092004 [arXiv:1503.03379] [INSPIRE]. [56] R. Catena, K. Fridell and M.B. Krauss, Non-relativistic Effective Interactions of Spin 1 Dark Matter,JHEP 08 (2019) 030 [arXiv:1907.02910] [INSPIRE]. [57] E. Del Nobile, Appendiciario — A hands-on manual on the theory of direct Dark Matter detection,arXiv:2104.12785 [INSPIRE]. [58] E. Del Nobile and F. Sannino, Dark Matter Effective Theory,Int. J. Mod. Phys. A 27 (2012) 1250065 [arXiv:1102.3116] [INSPIRE]. – 33 – JHEP07(2021)081 [59] A. De Simone, A. Monin, A. Thamm and A. Urbano, On the effective operators for Dark Matter annihilations,JCAP 02 (2013) 039 [arXiv:1301.1486] [INSPIRE]. [60] M. Duch, B. Grzadkowski and J. Wudka, Classification of effective operators for interactions between the Standard Model and dark matter,JHEP 05 (2015) 116 [arXiv:1412.0520] [INSPIRE]. [61] S. Matsumoto, S. Mukhopadhyay and Y.-L.S. Tsai, Effective Theory of WIMP Dark Matter supplemented by Simplified Models: Singlet-like Majorana fermion case,Phys. Rev. D 94 (2016) 065034 [arXiv:1604.02230] [INSPIRE]. [62] S. Matsumoto, S. Mukhopadhyay and Y.-L.S. Tsai, Singlet Majorana fermion dark matter: a comprehensive analysis in effective field theory,JHEP 10 (2014) 155 [arXiv:1407.1859] [INSPIRE]. [63] H. Han, H. Wu and S. Zheng, Effective field theory of the Majorana dark matter,Chin. Phys. C 43 (2019) 043103 [arXiv:1711.10097] [INSPIRE]. [64] A. Belyaev et al., Interplay of the LHC and non-LHC Dark Matter searches in the Effective Field Theory approach,Phys. Rev. D 99 (2019) 015006 [arXiv:1807.03817] [INSPIRE]. [65] J. Brod, A. Gootjes-Dreesbach, M. Tammaro and J. Zupan, Effective Field Theory for Dark Matter Direct Detection up to Dimension Seven,JHEP 10 (2018) 065 [arXiv:1710.10218] [INSPIRE]. [66] R. Harnik and G.D. Kribs, An Effective Theory of Dirac Dark Matter,Phys. Rev. D 79 (2009) 095007 [arXiv:0810.5557] [INSPIRE]. [67] J. Kopp, T. Schwetz and J. Zupan, Global interpretation of direct Dark Matter searches after CDMS-II results,JCAP 02 (2010) 014 [arXiv:0912.4264] [INSPIRE]. [68] J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T.M.P. Tait and H.-B. Yu, Gamma Ray Line Constraints on Effective Theories of Dark Matter,Nucl. Phys. B 844 (2011) 55 [arXiv:1009.0008] [INSPIRE]. [69] K. Cheung, P.-Y. Tseng, Y.-L.S. Tsai and T.-C. Yuan, Global Constraints on Effective Dark Matter Interactions: Relic Density, Direct Detection, Indirect Detection, and Collider, JCAP 05 (2012) 001 [arXiv:1201.3402] [INSPIRE]. [70] M.R. Buckley, Using Effective Operators to Understand CoGeNT and CDMS-Si Signals, Phys. Rev. D 88 (2013) 055028 [arXiv:1308.4146] [INSPIRE]. [71] A. Crivellin and U. Haisch, Dark matter direct detection constraints from gauge bosons loops,Phys. Rev. D 90 (2014) 115011 [arXiv:1408.5046] [INSPIRE]. [72] A. Crivellin, U. Haisch and A. Hibbs, LHC constraints on gauge boson couplings to dark matter,Phys. Rev. D 91 (2015) 074028 [arXiv:1501.00907] [INSPIRE]. [73] M.A. Fedderke, J.-Y. Chen, E.W. Kolb and L.-T. Wang, The Fermionic Dark Matter Higgs Portal: an effective field theory approach,JHEP 08 (2014) 122 [arXiv:1404.2283] [INSPIRE]. [74] J. Hisano, R. Nagai and N. Nagata, Effective Theories for Dark Matter Nucleon Scattering, JHEP 05 (2015) 037 [arXiv:1502.02244] [INSPIRE]. [75] S. Bhattacharya and J. Wudka, Effective Theories with Dark Matter Applications, arXiv:2104.01788 [INSPIRE]. – 34 – JHEP07(2021)081 [76] P.J. Fox, R. Harnik, J. Kopp and Y. Tsai, Missing Energy Signatures of Dark Matter at the LHC,Phys. Rev. D 85 (2012) 056011 [arXiv:1109.4398] [INSPIRE]. [77] J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T.M.P. Tait and H.-B. Yu, Constraints on Light Majorana dark Matter from Colliders,Phys. Lett. B 695 (2011) 185 [arXiv:1005.1286] [INSPIRE]. [78] J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T.M.P. Tait and H.-B. Yu, Constraints on Dark Matter from Colliders,Phys. Rev. D 82 (2010) 116010 [arXiv:1008.1783] [INSPIRE]. [79] C. Arina, A. Cheek, K. Mimasu and L. Pagani, Light and Darkness: consistently coupling dark matter to photons via effective operators,Eur. Phys. J. C 81 (2021) 223 [arXiv:2005.12789] [INSPIRE]. [80] F. Bishara, J. Brod, B. Grinstein and J. Zupan, Chiral Effective Theory of Dark Matter Direct Detection,JCAP 02 (2017) 009 [arXiv:1611.00368] [INSPIRE]. [81] A. Crivellin, F. D’Eramo and M. Procura, New Constraints on Dark Matter Effective Theories from Standard Model Loops,Phys. Rev. Lett. 112 (2014) 191304 [arXiv:1402.1173] [INSPIRE]. [82] R.J. Hill and M.P. Solon, Standard Model anatomy of WIMP dark matter direct detection I: weak-scale matching,Phys. Rev. D 91 (2015) 043504 [arXiv:1401.3339] [INSPIRE]. [83] N.F. Bell, Y. Cai and A.D. Medina, Co-annihilating Dark Matter: Effective Operator Analysis and Collider Phenomenology,Phys. Rev. D 89 (2014) 115001 [arXiv:1311.6169] [INSPIRE]. [84] M.J. Baker et al., The Coannihilation Codex,JHEP 12 (2015) 120 [arXiv:1510.03434] [INSPIRE]. [85] M. Cirelli, N. Fornengo and A. Strumia, Minimal dark matter,Nucl. Phys. B 753 (2006) 178 [hep-ph/0512090] [INSPIRE]. [86] M. Cirelli and A. Strumia, Minimal Dark Matter: Model and results,New J. Phys. 11 (2009) 105005 [arXiv:0903.3381] [INSPIRE]. [87] S. Bottaro, A. Strumia and N. Vignaroli, Minimal Dark Matter bound states at future colliders,JHEP 06 (2021) 143 [arXiv:2103.12766] [INSPIRE]. [88] S. Kanemura, S. Matsumoto, T. Nabeshima and N. Okada, Can WIMP Dark Matter overcome the Nightmare Scenario?,Phys. Rev. D 82 (2010) 055026 [arXiv:1005.5651] [INSPIRE]. [89] A. Djouadi, O. Lebedev, Y. Mambrini and J. Quevillon, Implications of LHC searches for Higgs-portal dark matter,Phys. Lett. B 709 (2012) 65 [arXiv:1112.3299] [INSPIRE]. [90] A. Djouadi, A. Falkowski, Y. Mambrini and J. Quevillon, Direct Detection of Higgs-Portal Dark Matter at the LHC,Eur. Phys. J. C 73 (2013) 2455 [arXiv:1205.3169] [INSPIRE]. [91] V. Silveira and A. Zee, Scalar phantoms,Phys. Lett. B 161 (1985) 136 [INSPIRE]. [92] J. McDonald, Gauge singlet scalars as cold dark matter,Phys. Rev. D 50 (1994) 3637 [hep-ph/0702143] [INSPIRE]. [93] C.P. Burgess, M. Pospelov and T. ter Veldhuis, The Minimal model of nonbaryonic dark matter: A Singlet scalar,Nucl. Phys. B 619 (2001) 709 [hep-ph/0011335] [INSPIRE]. – 35 – JHEP07(2021)081 [94] V. Barger, P. Langacker, M. McCaskey, M. Ramsey-Musolf and G. Shaughnessy, Complex Singlet Extension of the Standard Model,Phys. Rev. D 79 (2009) 015018 [arXiv:0811.0393] [INSPIRE]. [95] S. Andreas, C. Arina, T. Hambye, F.-S. Ling and M.H.G. Tytgat, A light scalar WIMP through the Higgs portal and CoGeNT,Phys. Rev. D 82 (2010) 043522 [arXiv:1003.2595] [INSPIRE]. [96] S. Baek, P. Ko and W.-I. Park, Invisible Higgs Decay Width vs. Dark Matter Direct Detection Cross Section in Higgs Portal Dark Matter Models,Phys. Rev. D 90 (2014) 055014 [arXiv:1405.3530] [INSPIRE]. [97] GAMBIT collaboration, Status of the scalar singlet dark matter model,Eur. Phys. J. C 77 (2017) 568 [arXiv:1705.07931] [INSPIRE]. [98] C. Gross, O. Lebedev and T. Toma, Cancellation Mechanism for Dark-Matter-Nucleon Interaction,Phys. Rev. Lett. 119 (2017) 191801 [arXiv:1708.02253] [INSPIRE]. [99] K.K. Boddy, J. Kumar, A.B. Pace, J. Runburg and L.E. Strigari, Effective J-factors for Milky Way dwarf spheroidal galaxies with velocity-dependent annihilation,Phys. Rev. D 102 (2020) 023029 [arXiv:1909.13197] [INSPIRE]. [100] G. Bélanger, F. Boudjema, A. Pukhov and A. Semenov, MicrOMEGAs 2.0: A Program to calculate the relic density of dark matter in a generic model,Comput. Phys. Commun. 176 (2007) 367 [hep-ph/0607059] [INSPIRE]. [101] G. Bélanger, F. Boudjema, A. Pukhov and A. Semenov, micrOMEGAs_3: A program for calculating dark matter observables,Comput. Phys. Commun. 185 (2014) 960 [arXiv:1305.0237] [INSPIRE]. [102] G. Bélanger, F. Boudjema, A. Goudelis, A. Pukhov and B. Zaldivar, micrOMEGAs5.0: Freeze-in,Comput. Phys. Commun. 231 (2018) 173 [arXiv:1801.03509] [INSPIRE]. [103] M. Ruhdorfer, E. Salvioni and A. Weiler, A Global View of the Off-Shell Higgs Portal, SciPost Phys. 8(2020) 027 [arXiv:1910.04170] [INSPIRE]. [104] M. Frigerio, A. Pomarol, F. Riva and A. Urbano, Composite Scalar Dark Matter,JHEP 07 (2012) 015 [arXiv:1204.2808] [INSPIRE]. [105] M. Drees, M.M. Nojiri, D.P. Roy and Y. Yamada, Light Higgsino dark matter,Phys. Rev. D 56 (1997) 276 [Erratum ibid. 64 (2001) 039901] [hep-ph/9701219] [INSPIRE]. [106] A. Djouadi, The Anatomy of electro-weak symmetry breaking. II. The Higgs bosons in the minimal supersymmetric model,Phys. Rept. 459 (2008) 1 [hep-ph/0503173] [INSPIRE]. [107] H. Baer, V. Barger and D. Mickelson, Direct and indirect detection of higgsino-like WIMPs: concluding the story of electroweak naturalness,Phys. Lett. B 726 (2013) 330 [arXiv:1303.3816] [INSPIRE]. [108] A. Joglekar, P. Schwaller and C.E.M. Wagner, Dark Matter and Enhanced Higgs to Di-photon Rate from Vector-like Leptons,JHEP 12 (2012) 064 [arXiv:1207.4235] [INSPIRE]. [109] A. Carmona and M. Chala, Composite Dark Sectors,JHEP 06 (2015) 105 [arXiv:1504.00332] [INSPIRE]. [110] G. Ballesteros, A. Carmona and M. Chala, Exceptional Composite Dark Matter,Eur. Phys. J. C 77 (2017) 468 [arXiv:1704.07388] [INSPIRE]. – 36 – JHEP07(2021)081 [111] J.C. Criado, A. Djouadi, N. Koivunen, K. Müürsepp, M. Raidal and H. Veermäe, Confronting spin-3/2and other new fermions with the muon g−2measurement, arXiv:2104.03231 [INSPIRE]. [112] Particle Data collaboration, Review of Particle Physics,Phys. Rev. D 98 (2018) 030001 [INSPIRE]. [113] F. del Aguila, J.A. Aguilar-Saavedra and R. Pittau, Heavy neutrino signals at large hadron colliders,JHEP 10 (2007) 047 [hep-ph/0703261] [INSPIRE]. [114] ATLAS collaboration, Search for type-III seesaw heavy leptons in dilepton final states in pp collisions at √s= 13 TeV with the ATLAS detector,Eur. Phys. J. C 81 (2021) 218 [arXiv:2008.07949] [INSPIRE]. [115] CMS collaboration, Search for physics beyond the standard model in multilepton final states in proton-proton collisions at √s= 13 TeV,JHEP 03 (2020) 051 [arXiv:1911.04968] [INSPIRE]. [116] ATLAS collaboration, Search for metastable heavy charged particles with large ionisation energy loss in pp collisions at √s= 8 TeV using the ATLAS experiment,Eur. Phys. J. C 75 (2015) 407 [arXiv:1506.05332] [INSPIRE]. [117] ATLAS collaboration, Search for heavy charged long-lived particles in the ATLAS detector in 36.1fb−1of proton-proton collision data at √s= 13 TeV,Phys. Rev. D 99 (2019) 092007 [arXiv:1902.01636] [INSPIRE]. [118] V.A. Mitsou, MoEDAL, FASER and future experiments targeting dark sector and long-lived particles,PoS LHCP2020 (2021) 112 [INSPIRE]. [119] M. Chala, F. Kahlhoefer, M. McCullough, G. Nardini and K. Schmidt-Hoberg, Constraining Dark Sectors with Monojets and Dijets,JHEP 07 (2015) 089 [arXiv:1503.05916] [INSPIRE]. [120] J.C. Criado, BasisGen: automatic generation of operator bases,Eur. Phys. J. C 79 (2019) 256 [arXiv:1901.03501] [INSPIRE]. – 37 –