Closing in on minimal dark matter and radiative neutrino masses

We study one-loop radiative neutrino mass models in which one of the beyond-the-standard model fields is either a hypercharge-zero fermion quintet (minimal dark matter) or a hypercharge-zero scalar septet. By systematically classifying all possible one-loop such models we identify various processes that render the neutral component of these representations (dark matter) cosmologically unstable. Thus, our findings show that these scenarios are in general not reconcilable with dark matter stability unless tiny couplings or additional ad hoc symmetries are assumed, in contrast to minimal dark matter models where stability is entirely due to the standard model gauge symmetry. For some variants based on higher-order loops we find that α2 reaches a Landau pole at rather low scales, a couple orders of magnitude from the characteristic scale of the model itself. Thus, we argue that some of these variations although consistent with dark matter stability and phenomenological constraints are hard to reconcile with perturbativity criteria.


JHEP06(2016)108
which the SM and the dark sector transform differently, thus guaranteeing stability. These symmetries can have multiple origins, and can be regarded as remnants of the spontaneous symmetry breaking of a more fundamental symmetry related with e.g. GUTs [13,14] or flavor theories [15].
Another approach that has been explored relates radiative neutrino masses with higherorder electroweak (EW) representations [16][17][18][19][20][21][22][23][24]. In this case the idea is different and relies on the mechanism underlying minimal DM models [25][26][27], namely for higher-order EW representations renormalizable interactions with SM operators cannot be written, and so they are absolutely stable at the renormalizable level (at the renormalizable level the minimal DM Lagrangian exhibits an accidental Z 2 symmetry). Decays are induced by effective operators, thus in these scenarios DM is never absolutely stable but can be cosmologically stable provided the effective operator is sufficiently suppressed. In [25] it was shown that the stability condition, the requirement of EW gauge coupling perturbativity up to M Planck and direct DM searches constraints, singles out two possible representations: hyperchargezero fermion quintet and hypercharge-zero scalar septet, of which the fermion quintet defines the minimal DM model. Recently, in refs. [28] and subsequently in [29] it was shown that for the scalar septet one can write a loop-induced lower dimensional operator that renders the septet a non-viable DM representation, even if the effective scale is assumed to be M Planck . Thus, within the minimal DM context the only viable representation is the fermion quintet, although recently analyses have shown that depending on the DM profile this representation is not consistent with indirect DM searches results [30,31] (see discussion in section 2).
In this paper we show that when the fermion quintet acts as a mediator in one-loop neutrino mass generation, there are always DM decay operators that -under reasonable parameter choices -lead to fast DM decays, regardless of the one-loop neutrino mass model. Some of these operators arise at the tree level but others are loop-induced. The latter being specially important in models where in addition to the fermion quintet, representations beyond quartets are present. We will show that as soon as the minimal DM scheme is extended to include other representations (that allow the construction of one-loop neutrino mass models) there are always Z 2 -breaking couplings, thus fast decay modes are always expected to be present. We consider the case of the scalar septet as well, despite being not consistent with stability even in the minimal DM framework. The reason is that these results allow the identification of DM fast decay modes induced by neutrino physics itself, rather than due to a different type of physics (related with unknown quantum gravity effects).
The rest of this paper is organized as follows. In section 2, we briefly review minimal DM models, introduce our notation and present formulas that will be used in our estimations. In section 2.1, we present the different Z 2 -breaking decay operators and estimate lifetimes for various processes, showing that even for neutrino physics parameters taken to their extreme values the decays are always fast. In section 3, we systematically derive all one-loop realizations of the Weinberg operator, under the condition of the UV completion containing either the fermion quintet or scalar septet. With these results at hand, we then identify the different DM decay operators associated with each of the neutrino mass JHEP06(2016)108 model categories. In section 4, we discuss the main pitfalls of models for Majorana neutrino masses with higher-order EW representations and briefly comment on models with higher-order loops. Finally, in section 5 we summarize and present our conclusions.

Minimal dark matter and possible decay processes
Minimal DM models rely on the observation that for higher-order SU(2) representations, R, renormalizable gauge invariant operators of the form where O SM is an operator involving just SM fields, cannot be written [25][26][27]. Thus, when the SM is endowed with such state, its stability is automatically guaranteed. The stabilization mechanism at work, resembles the SM "mechanism" which assures proton stability, namely the gauge symmetry does not allow for renormalizable operators that might induce decay. Of course, at the non-renormalizable level several effective operators of the form (2.1) can be written, but their renormalizable realizations require degrees of freedom which the minimality criteria do not allow for. 3 The lowest-order representation for which the above argument proves to be true and yields a viable DM candidate is the fermionic R = 5. Let us discuss this in more detail. For R = 3, renormalizable operators can always be written, regardless of whether they are fermions (F ) or scalars (S) and provided Y = 0, 1. For the quartet and beyond it is more subtle. Using the notation R Y X (with X referring to fermion (F ) or scalar (S) and Y to hypercharge, normalized according to: Q = T 3 + Y ), and bearing in mind that one finds that for the quartet non-renormalizable operators are possible depending on its spin, namely However, though involving non-renormalizable operators, R F = 4 F does not provide a viable DM candidate (the R S = 4 S involves renormalizable couplings and so it is of course non-viable either). First of all, even if the cutoff scale is taken to be Λ = M Planck (where arguably quantum gravity effects will generate such operator) its suppression is not sufficiently strong to guarantee cosmological stability. Furthermore, in the context of minimal DM only Y = 0 representations are consistent with direct DM searches: for Y = 0 states,

JHEP06(2016)108
Z-mediated processes induce a tree-level spin-independent DM-nucleon cross section well above the values allowed by current DM direct searches [25][26][27]. We stress that the validity of this statement is subject to the minimality condition. Departures from this assumption, e.g. by introducing another state that mixes with the DM producing a small mass splitting, enable Y = 0 representations (for more detailes see [27]). For R F = 5 F one finds 2 ⊗ 2 ⊗ 2 ⊗ 2 ⊃ 5, which then implies that SU(2) invariance requires four SM fields.
Depending on the spin of the representation, the effective operator in (2.1) is N = 5 (scalar) or N = 6 (fermion), with its full structure determined by the hypercharge of the corresponding state. An upper limit on the size of the representations can be determined by the condition of α 2 being perturbative up to the Planck scale. Ref. [28] has shown that depending on the representation a precise determination of the Landau pole should rely on two-loop RGE. Based on that analysis representations R F > 5 F and R S > 8 S are ruled out, as has been previously pointed out in [25][26][27] using a one-loop RGE calculation. One is left then with R < 5 F,S and 7 S for which the decay lifetimes can be estimated from the effective operator where the Wilson coefficient c N , assumed to be order one for simplicity, can involve flavor indices too. In the limit m DM v (with v = H = ( √ 2G F ) −1/2 246 GeV), the total decay width can be parameterized as follows [28] Γ (N ) where n c refers to the number of Higgs condensate insertions. Of the three potentially viable representations only those for which τ DM 10 26 sec (as required by indirect detection experiments looking for γ-ray, ν, e + or p − signals stemming from DM decay [33][34][35][36]) are consistent. From (2.4) their lifetimes can be calculated. Consistency with τ DM 10 26 sec for 5 0 S requires scales way above M Planck , in contrast 5 0 F and 7 0 S have sufficiently long lifetimes even for scales well below the Planck scale, due to their decays being driven by dim = 6 and dim = 7 operators as can be seen in figure 1.
In the absence of any other degree of freedom, these results hold if no lower dimension operator for these representation can be written. For 5 0 F this proves to be true, for 7 0 S , instead, not. The point is that the previous analysis rely on tree level effective decay modes, but for 7 0 S one-loop radiative dim = 5 operators can be written through the operator 7 0 S (7 0 S 7 0 S )H † HH † H [28] (see discussion below). Thus, this argument rules out the 7 0 S as well. Minimal DM therefore reduces to a single possibility, 5 0 F , for which stringent constraints from indirect DM searches have been found. Of particular relevance are those coming from γ-ray line searches from the galactic center, for which it has been proven that if the Milky Way possesses a Navarro-Frenk-White or Einasto DM profiles this representation is ruled out too [30,31]. 4 It can be however consistently considered in the context of cored profiles such as Burket or Isothermal.

Minimal DM decay in the presence of extra degrees of freedom
Minimal DM instability is induced by UV completions associated with effective operators of the type (2.3). In the minimal model the new degrees of freedom are assumed to be related with quantum gravity effects, something that guarantees stability for 5 0 F . For UV completions associated with "lighter" physics, DM can become cosmologically unstable. For Λ 10 15 GeV, 5 0 F has a lifetime below 10 26 sec. The neutrino mass one gets from one-loop realizations of the Weinberg operator can be estimated as where Y denotes a generic Yukawa coupling. Thus, assuming Y ⊂ [10 −2 , 1] and fixing the neutrino mass to m Atm = 0.05 eV [37][38][39], one can estimate Λ [10 5 , 10 13 ] GeV. Note that in this estimation the one-loop function has been neglected, but in a specific model with an explicit neutrino mass matrix the presence of the loop function will lead to smaller cutoff scales. Thus, with the value derived from (2.5) we are overestimating the DM lifetime. This result then shows that -in principle -any one-loop UV completion of the Weinberg operator that can yield in turn an UV completion for the effective operator responsible for 5 0 F decays will lead to: and therefore to unviable "minimal DM". 5 This conclusion can be understood as follows.
In the absence of further degrees of freedom the minimal DM full Lagrangian reads:

JHEP06(2016)108
This Lagrangian is invariant under the Z 2 transformations 5 F → −5 F and X SM → X SM , that entirely results as a consequence of the gauge symmetry. If the extra representations that allow for the one-loop neutrino mass matrix allow as well for UV completions of the following Lagrangian they will necessarily involve interactions that explicitly break the accidental Z 2 symmetry (as this operator does) and therefore will lead to DM decay processes with typical lifetimes given by (2.6). Instability of the 7 0 S can be understood in the same way. The full Lagrangian is given by where the ellipses refer to additional terms, that as the ones we have explicitly written are In this case the relevant effective Z 2 -breaking term is This operator alone cannot yield DM decays, but when combined with the EW invariant H † H leads to a one-loop dim = 5 decay operator [28]. As in the fermionic DM case, instability of the 7 S can be understood as a consequence of renormalizable Z 2 -violating interactions that in the effective limit reduce to (2.10).
As already pointed out, decay modes depend upon the extra representations and can be sorted according to tree level and one-loop induced decay modes. In what follows whenever calculating lifetimes we will set m DM = 10 4 GeV and the mass of the extra representations, assumed to be universal (m), to 10 9 GeV, which is typically the largest value for the extra degrees of freedom one will get in a specific neutrino mass model. With the relevant masses fixed in this way, the Z 2 -conserving Yukawa couplings (those related with neutrino mass generation: Y ν ) are fixed to one. 6 We start our analysis with the tree level modes: • Tree level decay modes for 5 0 F . These DM decay processes are mediated by scalar or fermion Y = 1/2 quartets and sextets, as required by gauge invariance. In both cases one can distinguish processes with no additional representations (single messenger scenarios) and processes with further representations, e.g. scalar triplets or any other enabled by gauge invariance (multiple messenger scenarios). In single messenger processes, DM decays proceed via off-shell states that in turn directly decay to SM particles, through either renormalizable or non-renormalizable couplings (examples are shown in figure 2, diagrams (a), (c), (d)). In multiple messenger scenarios, instead, DM decay via cascade processes mediated by several off-shell states, see figure 2 diagram (b) for an example. Note that one can as well write multiple messenger decay modes for the sextet by expanding the effective couplings c S 6 and c F 7 in diagrams (c) and (d) in figure 2. In what follows we discuss in more detail the examples shown in figure 2, aiming at showing that in models involving quartets tree level DM decay 6 This parameter choice is a simplification that enables determining the largest lifetimes achievable in each case. In practice, however, these scales are not entirely independent and instead are constrained by phenomenological conditions such as the DM relic abundance. processes always lead to fast DM decay, while for those with sextets these processes do not lead to rapid decay modes, therefore motivating the investigation of DM decay loop-induced processes.

JHEP06(2016)108
Of particular relevance for one-loop neutrino mass models is 4 figure 2). Depending on the different SU(2) contractions this operator chain contains several decay processes, of which l ± W ∓ L Z L Z L is the dominant decay mode. The lifetime can be estimated by inserting c 6 = Y · λ in eq. (2.4) and trading Λ to m. The result reads:  Here h refers to a generic Yukawa coupling (see figure 2, diagram (b)). In single messenger scenarios containing 6 1/2 S,F , tree level processes involve dim = 8 operators (see figure 2, diagrams (c) and (d)). The decay processes will involve two distinct scales, namely m and the scale at which the effective coupling c S 6 or c F 7 is generated. This scale can be related with either the scale at which α 2 = g 2 2 /4π 2 becomes nonperturbative or the Planck scale. Depending on the UV completion leading to c S 6 or c F 7 the former can be rather low, and thus in some instances χ DM decay can be fast. If taken to be M Planck , these processes cannot threaten DM stability. All in all, if one relies only on tree level processes one ends up with the conclusion that scenarios involving representations beyond the quartet (scalar or fermionic) are consistent with DM stability. This conclusion however is misleading since one-loop induced decay processes can always be written and can be fast enough to render χ DM cosmologically unstable.

JHEP06(2016)108
Loop-induced decay processes are of more relevance in models that contain a 6 1/2 S . At the tree level these scenarios involve processes that lead to slow DM decay, and so seem consistent with cosmological stability. In the effective limit, the loop-induced processes are associated with dim = 6 effective operators whose cutoff scale amounts (at most) to m = 10 9 GeV. Thus, an estimation of the DM lifetime can be done from eq. (2.4) by rescaling by the loop factor. For processes such as ν W ± W ∓ Z L the result reads: Which shows that although the tree level decay does not yield fast decays, the loopinduced processes lead to lifetimes that do not amount to 10 26 seconds, under "reasonable" parameter choices. This conclusion is inline with what was first pointed out in ref. [18], showing the infeasibility of the RνMDM [16]. In the case of multiple messenger fields one finds Yukawa couplings of the type R S , which break the Z 2 accidental symmetry. Their presence allow the construction of operators as those shown in figure 3 (diagram to the right), that induce again processes such as χ DM → ν W ± W ∓ Z L . The lifetime for these processes amounts to the value of the single messenger case in (2.13).
We now turn to the 7 0 S DM scenario. As we have already stressed, in the context of minimal DM this representation is not consistent with stability even for Λ = M Planck . Thus, even if a one-loop neutrino mass model does not involve degrees of freedom that can generate at the renormalizable level the operator in (2.10) -arguably -quantum gravity effects will generate it and therefore ϕ DM ⊂ 7 0 S will decay at a fast rate [28]. It is however worth identifying those extra representations that enable writing UV completions of this effective operator, in particular because these results will enable in turn the identification of oneloop neutrino mass models where ϕ DM decays will be due to neutrino physics itself, rather than due to a different type of physics.
As in the fermionic DM scenario, in this case one can identify various Z 2 -breaking terms whose structure depends on the extra representations present. We have found trilinear scalar couplings of the type 7 0 S R 1/2 S (R S ) −1/2 and Yukawa couplings 7 0 S R 1/2 F (R F ) −1/2 , with R = R = 6, 8 and R = 6 and R = 8 (for both, fermions and scalars). With these vertices one can then construct a certain number of loop-induced processes that lead to fast ϕ DM decays. The different operators are shown in figure 4. The different cases can be sorted in single and multiple messenger scenarios, but due to our universal mass simplification they lead to the same lifetime. Possible decay modes for operators to the left in figure 4 are W ± L W ∓ L , while for those to the right l ± l ∓ . In the effective limit, m m DM , the decay lifetime for processes induced by operators as those shown in the diagram to the left in figure 4 has been calculated as [28]:  For processes with leptons modes (induced by the operators to the right in figure 4) the lifetime is larger, as can be readily understood by realizing that the effective operator is dim = 6 rather than 5, as in the pure scalar case. Though larger, the values are always well below what cosmological stability demands.

UV completions of one-loop realizations of the Weinberg operator
In order to be model-independent one has to consider all possible one-loop realizations of the lepton-number-breaking dim = 5 operator. The possibilities have been systematically sorted in [11], and diagrammatically their number reduces to four independent diagrams as shown in figure 5. We start by assuming that one of the loop degrees of freedom is either R = 5 0 F or R = 7 0 S . As soon as this internal representation is fixed the remaining representations are fixed by SU(2) × U(1) Y invariance. It is worth emphasizing that in some instances the resulting representations allow for tree level realizations of higherorder lepton-number-breaking operators (see e.g. ref. [17]). Their contributions to neutrino masses are however negligible provided the BSM fields have masses above ∼ 3 TeV.

The hypercharge-zero quintet case
Results for the possible UV completions that contain a hypercharge-zero fermionic quintet are shown in table 1. The different models can be readily derived with the aid of the SU(2) product decomposition 2 ⊗ R = R ± 1 , In terms of the necessary number of representations, UV completions based on D1 define minimal scenarios. UV completions based on D2 require at least an additional scalar representation, and in contrast to those based on D1 are such that the representation content always allow for two different contributions to the effective neutrino mass matrix: where m D1  F an extra fermionic representation that can be 3 1 F , 5 1 F or 7 1 F . As regards the Z 2 -breaking decay modes, as can be seen in table 1 all possible models always involve representations that enable writing decay modes as those discussed in section 2.1 (see figure 2-4). In table 1 we have collected models according to the category they belong to, the category defined by diagrams D1-D4, and according to the decay mode they lead to (we have specified only loop-induced processes). Representations yielding a Z 2 -violating process can be identified from the following color scheme: yellowish (greenish) cells for representations mean that those representations induce a Z 2 -breaking mode 1 (2).

JHEP06(2016)108
Fermion quintet radiative UV completed models fl-i sl-j Table 1. One-loop UV completions for models with a fermion quintet defined by diagrams D1-D4 in figure 5. For each diagram and for one of the internal fermions fixed to be 5 0 F , the remaining possible representations are listed according to the notation R Y a (with R referring to the SU(2) representation, a = F, S and Y to hypercharge). Representations not allowed by gauge invariance are indicated by a dash, while allowed representations are indicated by checkmarks. For the Z 2breaking decay modes we consider only loop-induced processes. Representations that induce a decay mode are indicated by colors. For example, in the first row 4

The hypercharge-zero septet case
We now turn to the discussion of the scalar septet. As we have already stressed, for this representation there is a loop-induced effective dim = 5 operator which even in the absence of an UV completion (in "pure" minimal DM) leads to fast DM decay. Despite that, here for completeness we identify the different one-loop neutrino mass models for which UV completions of that effective decay operator can be written.
The main results in this case are summarized in table 2. The specification of the models in these cases goes along the same lines that in the quintet DM scenarios, with additional SU(2) product decompositions (e.g. 2 ⊗ 7 = 6 ⊕ 8) given by eq. (3.1). With the aid of these product rules the different UV completions, which we do not list, can be easily derived. Their main features are the following. With one of the loop degrees of freedom fixed according to 7 0 S , UV completions based on D1 are of two types, those involving 6 1 F and those with 8 1 F . Models relying on D2 are of two types: models where the mass matrix is entirely determined by D2-type contributions and scenarios where the neutrino mass matrix has the form: In D3-based scenarios the representations required imply in all cases a neutrino mass matrix of the form Finally, for UV completions involving D4 diagrams the neutrino mass matrix involves only D4 contributions.

JHEP06(2016)108
Scalar septet radiative UV completed models Table 2. One-loop UV completions for models with a scalar septet defined by diagrams D1-D4 in figure 5. For each diagram and for one of the internal scalars fixed to be 7 0 S , the remaining possible representations are listed according to the notation R Y a (with R referring to the SU(2) representation, a = F, S and Y to hypercharge). Representations not allowed by gauge invariance are indicated by a dash, while allowed representations are indicated by checkmarks. For the Z 2breaking decay modes we consider only loop-induced processes. Representations that induce a decay mode are indicated by colors. For example, in the first row 6 Regarding Z 2 -breaking modes, as shown in table 2 all UV completions of one-loop realizations of the Weinberg operator always involve a Z 2 -violating decay mode. As we have pointed out in section 2.1, compared with the fermion quintet these processes lead to faster DM decay since they are related with dim = 5 decay operators. The color scheme follows the same conventions of table 1: yellowish (greenish) cells for representations mean that those representations induce a Z 2 -breaking mode 1 (2), with the different modes referring to figure 4.

Pitfalls of DM from higher-order representations and neutrino masses
The different decay modes we have identified depend upon neutrino-related parameters and Z 2 -breaking couplings. In the estimations of lifetimes we have taken the former to their extreme values, while keeping perturbative Yukawa couplings (|Y ν | < 1). Thus, if one takes the extra representations to be heavier the largest neutrino mass will be below m Exp Atm 0.05 eV. The only way that enables the increasing of the lifetime is through couplings not related with neutrino physics that can be either Z 2 -breaking or Z 2 -conserving. For neutrino-related parameters taken to their extreme values, couplings of order 10 −11 or so (see e.g. eq. (2.13)) will lead to lifetimes of order 10 26 seconds. Stability, however, is not assured anymore by an accidental symmetry but from the smallness of a parameter. Whether such small coupling is justified depends on the UV completed model, e.g. it could be that in some UV realizations the coupling is radiatively induced or generated from the spontaneous symmetry breaking of a symmetry present in the UV completed theory (it has a dynamical origin 7 ), case in which a small value would be justified. One could as well JHEP06(2016)108 argue that for Z 2 -breaking couplings their smallness is technically "natural", in 't Hooft's sense [41], since tiny values increase the symmetry of the Lagrangian. This statement, however, depends on their RGE stability and therefore on their nature. For Z 2 -breaking Yukawas, diagrams contributing to RGE running are proportional to the Yukawa itself (or certain power) and so one expects their values to be RGE stable. For scalar quartic couplings, instead, this is not necessarily the case and strongly depends on the dimension of the representations involved [42]. Thus, in the absence of an explicit UV completed model DM slow decays are probably better justified if they arise from small Yukawas.
Note that even smaller couplings are required if one lowers the mass scale of the messenger fields. The point is that with m ∼ 10 9 GeV any possible low-energy indirect probe (e.g. lepton-flavor-violating or collider-related observables) is well below current and future experimental sensitivities. Thus, testability of these scenarios requires m 10 9 GeV. For scales comparable to those of DM the lifetimes are reduced by several order of magnitudes and so smaller couplings, order 10 −21 or so, are required to achieve consistent lifetimes. Thus, in summary, in neutrino mass models that contain extra representations as those we have identified in the previous section (in addition to 5 0 F or 7 0 S , or both) one expects DM to decay fast, with the decays determined by Z 2 -violating operators (see figures 2-4). The explicit breaking of the accidental Z 2 symmetry renders DM cosmologically unstable. Tuning of couplings is possible, but in that case cosmological stability is unrelated with the mechanism underlying minimal DM models.
A more plausible way of relating Majorana neutrino masses with accidentally stabilized DM is through higher-order loops. In this case, SU(2) × U(1) Y invariance is less restrictive due to the topologies of the possible diagrams. In the two-loop case, the topologies allow for overall shifts in the hypercharge of the fields flowing in the loops, while keeping Y = 0 for the DM representation (for the different topologies see ref. [12]). This "hypercharge freedom" might enable the construction of consistent two-loop models. For three-loop models, as far as we are aware, two examples have been pointed out [22,23]. Ref. [22] considers a model with three 7 0 F , one 7 1 S and one 1 1 S . The neutral component of the lightest 7 0 F is cosmologically stable due to an accidental Z 2 symmetry. Ref. [23], instead, considers a model with three copies of the 5 0 F , one copy of the 7 0 S and one copy of the 5 1 S , in addition to a two-Higgs-doublet model. In this case, the neutral component of the lightest 5 0 F is cosmologically stable. Although consistent with DM stability, neutrino physics data and low-energy lepton-flavor-violating constraints, we have found that in both cases α 2 = g 2 2 /4π reaches a Landau pole at relatively low energy scales: Λ (1) Landau-pole 10 m X , Λ (2) Landau-pole 10 2 m X (see figure 6), where m X refers to the masses of the heavy representations which we have taken to be universal, for simplicity. This result has been derived by integrating two-loop RGEs and neglecting Yukawa couplings contributions. Note that since in threeloop models Yukawas are order one, their inclusion will change the location of the Landau pole, with the precise value increasing or decreasing depending on numerical details.
This result then shows that in these scenarios perturbativity is lost at scales a couple orders of magnitude higher from the characteristic energy scale of the models (taken to be m X ), with such behavior being more pronounced for the model in [22]. In the absence of a strongly coupled SU(2) × U(1) Y sector at that scale, ideally one should demand pertur- . Energy scale at which α 2 reaches a Landau pole as a function of m X /m DM in the three-loop neutrino models of [22,23] (m X refers to the masses of the heavier fermion and scalars). The plot to the left is the result for the model in [22] with the DM mass fixed to 20 TeV and the singlet scalar (1 1 S ) mass to 500 GeV. The plot to the right is for [23] with the DM mass fixed to 10 TeV and the second doublet mass to 1 TeV. These results have been derived from two-loop RGEs and neglecting Yukawa contributions [43] (see text for more details).
bativity to hold at least up to the GUT scale, where one could expect that a new gauge sector assures α 2 < 1. Thus, arguably, this may be a drawback of these approaches. The presence of multiple higher-order representations can lead to non-perturbative scenarios at relatively low energy scales. There is of course something one should bear in mind. One could in principle understand these scenarios as effective limits of a larger UV theory, in which it could be that by taken into account all the degrees of freedom and intercations perturbativity could be restored. Furthermore, in that theory it could be as well that the Z 2 -violating couplings that we have identified as responsible for DM fast decay are either absent or dynamically suppressed, case in which DM will become cosmologically stable.

Conclusions
In this paper, we have presented a model-independent study of one-loop neutrino mass models in which one of the loop messengers is a hypercharge-zero fermion quintet or a hypercharge-zero scalar septet. In the absence of additional higher-order EW representations the neutral component of the quintet is absolutely stable at the renormalizable level, and cosmologically stable when the leading-order effective operator responsible for its decay (dim = 6) is included (assuming a cutoff above 10 15 GeV, see figure 1). We have shown that such feature no longer holds in the presence of the extra degrees of freedom that define the one-loop neutrino mass model, regardless of the model. For completeness we have extended our analysis to include the case of a hypercharge-zero scalar septet, although this representation cannot be reconciled with cosmological stability even in minimal DM scenarios.
We have systematically classified all such models and the corresponding DM decay operators in each case. Our main findings are summarized in tables 1 and 2, combined with the different DM decay operators given in figures 3 and 4. From these results our main JHEP06(2016)108 conclusion is that one-loop neutrino mass models that use minimal DM representations are not consistent with DM stability, unless tiny couplings -O ∼ 10 −21 -10 −11 (with the value depending on the highest scale of the corresponding neutrino mass model) -or ad hoc symmetries are assumed. The former being hard to reconcile in the absence of a dynamical mechanism assuring such extreme small values, while the latter missing the whole point behind higher-order EW representations. 8 We have commented on neutrino mass models variations based on higher-order loops. We argued that in the two-loop case, models consistent with DM stability might be possible writing due to the "hypercharge freedom" that the two-loop topologies offer. We have stressed that these models and three-loop scenarios might be hard to reconcile with perturbativity criteria. We have illustrated that for two specific three-loop models for which we have found that α 2 reaches a Landau pole at relatively low energy scales (figure 6), Λ Landau-pole ∼ (10-10 2 )m X , with m X the characteristic scale of the model.