Radiative two-loop neutrino masses with dark matter

Using the Weinberg operator, we present a full collection of genuine two-loop models for neutrino mass generation, which contain a dark matter particle as one of the internal messengers. These models can be constructed simply by adding new fields that are singlets or doublets of SU(2)L. We ensure the stability of the dark matter candidate by the addition of a Z2 symmetry that will also be used to forbid tree level or one-loop diagrams. Thus we only present models where the main contribution for neutrinos masses is generated from the corresponding two-loop diagram. We also discuss a short outline corresponding to some phenomenological characteristics of these models.


Introduction
The existence of neutrino masses and dark matter (DM) have corroborated the need for physics beyond the standard model (BSM). While neutrino experiments [1] have continuously supplied us with better information about parameters of the neutrino sector, direct DM detection experiments [2][3][4], given the lack in the direct detection signals, have only provided with bounds on the DM mass. Nevertheless, cosmological evidence [5] has shown arguably enough proof of its existence. With the addition of each new field (or any source of new physics), there will be a growing number of parameters added to the Standard Model (SM). Therefore, connecting this expanding number of parameters via physical arguments (that should arise from concrete models) helps reduce the freedom of these variables by tracing them to a single origin. With this idea in mind, many proposed models have attempted to connect neutrino masses to JHEP04(2017)148 DM [6][7][8][9][10][11][12][13][14], one of the most appealing of which is the Ma-Scotogenic model [15]. Radiative neutrino models as the above mentioned are very attractive, they not only give a natural explanation to the smallness of neutrino masses, but also given the new fields added, a DM candidate may be used as a messenger, ensuring a relation between neutrinos and DM.
The exact nature of neutrinos is still unclear; they may behave similarly to the charged leptons, i.e., Dirac particles, or entirely different by being their own antiparticles, referred to as Majorana particles. The existence of heavy Majorana neutrinos is a very appealing scenario to explain the smallnesses of SM neutrino masses [16].
If one wants to write an effective theory for Majorana neutrino masses, one can start with the SM Lagrangian and add higher order non-normalizable operators [17], where Λ is the energy scale, C n is a constant and O n is an operator of order n.
It can be shown [18] that there is only one possibility for n = 5, commonly referred in literature as the Weinberg dim=5 operator, where H is the Higgs field, L the SM lepton field, and the SU(2) L contractions have not been written for simplicity. This operator violates lepton number by two units and after the Higgs gets a vacuum expectation value (VEV), H = v, one can use this operator to calculate Majorana neutrino masses. Given the smallness of neutrino masses, assuming a mass of the order m ν ∼ O(10 −1 ) eV [1], it is possible to conclude that for the dim=5 operator, the constant C on eq. (1.1) is very small and proportional to the mass scale C ∼ 10 −12 × (Λ/TeV). This equation shows that is possible to have small neutrino masses at tree level with a very high scale, Λ ∼ (10 15 -10 16 ) GeV, where new physics is required. Conversely, one can reproduce small neutrino masses by using a radiative method, with C ∼ (y 2 /4π) m and m ≥ 2, where m = 2 corresponds to one-loop, m = 3 to two-loop and so on. An approximate calculation shows that for m = 5 (i.e., four-loop diagrams) the value of the Λ scale is too low (∼ eV), in other words, within the dim=5 Weinberg operator, it is only possible to reproduce the scale of neutrino masses up to models with three-loops. Also higher dimensional operators are possible [19][20][21][22], although the Weinberg operator remains a more elegant and simple way to generate the masses.
At tree level Majorana neutrinos can acquire mass via three mechanisms [23], called seesaw type-I [24][25][26][27], -II [28][29][30][31][32] and -III [33,34]. To be able to recreate neutrino experimental data type-I needs the addition of (at least) two fermion singlets, type-II requires one scalar triplet and type-III needs (at least) two fermion triplets. At one-loop there are four possible diagrams for the Weinberg operator [35], a study of these one-loop radiative neutrino mass models with viable dark matter candidates was done in ref. [36]. At the two-loop level the Weinberg operator has twenty genuine different diagrams (by "genuine" two-loop models we mean the ones for which, given a specific set of fields, there are no JHEP04(2017)148 contributions to the neutrino masses either from tree level (seesaw) or from one-loop.); a systematic classification of these two-loop realizations was presented in ref. [37].
As mentioned above, another hint for the need of new physics is the cosmological evidence of a new type of matter. If DM is a particle, present constrains dictate that it should be color neutral, without electric charge (or a very small one [38][39][40][41]), and fit experimental constraints, as the ones coming from relic abundance and direct detection experiments [42]. Also the particle must be stable in the sense that must be long-lived (experimental bounds set this limit around seven orders of magnitude higher than the lifetime of the universe [43][44][45][46]). But given the lacking evidence of DM via direct detection, indirect detection or colliders, its basic nature has not been determined yet [47], in that sense DM could be scalar [48][49][50][51][52][53], fermionic [54][55][56], vectorial [57][58][59][60] or even have multiple components [61][62][63][64][65][66].
The electric charge of a particle inside a multiplet is calculated using its isospin I 3 and hypercharge Y . If we require a DM candidate to fit the restriction of non-existing electrical charge, Q = I 3 + Y /2 = 0, we arrive to the condition that only a particle where Y = −2 I 3 can be DM. If one wants a particle to be stable, the most common mechanism is to include an extra unbroken symmetry that will forbid Lagrangian terms that would induce decay. Particularly, the cyclic Abelian group Z n is most commonly used, Z 2 being the simplest of these symmetries. Ten years ago, a new method of stabilizing DM was introduced [67][68][69]: in this minimal dark matter scenarios the high SU(2) L representation of the new fields does not permit the construction of renormalizable operators that allow the decay of the neutral component of the new field. Nevertheless it was proven that this scenario is not compatible with neutrino masses at one-loop, although higher loops might still be possible [70][71][72].
Electroweak (EW) multiplets that contain a DM candidate will have tree level interactions with quarks via Z bosons, giving a direct detection cross section via nucleon recoil that are proportional to Y 2 (the hypercharge squared) [73]. Experimental results have eliminated these particles as DM, if their mass is bigger than 10 GeV, they would produce scattering cross sections big enough, that would have been seen given current limits. This condition eliminates multiplets with even number of fields (i.e., doublets, quartets etc.), except those with Y = 0 [67][68][69]. In models with multiple added fields it is possible to have other kind of candidates if they mix with the above-mentioned particles. Also, there is an exception in the case of the scalar doublet with hypercharge Y = 1. For these models, that have two scalar doublets (SM Higgs plus the one added) a mass splitting can be enforced between the scalar and pseudo-scalar [74], this can eliminate the coupling with the Z at tree level, leaving only loop induced detection, allowing the neutral component of the doublet to have a detection rate that is still lower than experimental bounds.
In this article we will present multiple models with two-loop neutrino masses and a DM candidate that helps to generate the mass of the neutrinos, all of the candidates being either fermionic or scalar singlets, or a scalar doublet (or a mixture of the latter two). All of these models are unique and present interesting phenomenological features consistent with current experimental bounds. We will follow closely the results presented in [37], although we have two main differences. The first difference is that while they have shown that there JHEP04(2017)148 are twenty independent ways to create genuine two-loop diagrams, we only present models for nineteen diagrams, the reason for this is that one of the diagrams does not give any viable model that fit our conditions. The second difference is that restrictions for genuine two-loop models presented in [37] get relaxed with the implementation of the Z 2 symmetry; one can see that some fields and interactions that where forbidden since they allow tree and one-loop contributions, are now allowed in our models.
In section 2 we present the models, including for all of them the corresponding Lagrangian interactions and the neutrino mass equations, while in section 3 we present a small review on phenomenological applications of some of the fields, that can be used to discriminate some models. In appendix A we show all the integrals that are needed to calculate the neutrino masses., and in appendix B we show two examples of the models to clarify some of the aspects of their construction.

The models
In this section we present all genuine two-loop neutrino models that are possible to construct by only adding to the SM colorless SU(2) L singlets and doublets. We do not have a specific rule to explicitly forbid the one-loop models from appearing, instead, we construct all possible models that generate a two-loop diagram and eliminate the ones with particle content that can be used to build any tree level or one-loop diagram. It is important to keep in mind, that for all of the models that will be presented, all BSM fields contain only colorless particles, we will not discuss particles that have SU(3) c charges as possible messenger inside the loops, doing so would increase the number of possible models greatly. Nevertheless, the extension of the models to include colored particles is quite simple: one just needs to assign the same SU(3) c charge to all particles around one of the loops, while simultaneously avoiding charging the particle that would be the DM candidate. This is always possible for all models that will be presented. These two-loop setups do allow for a portion of the particles to be charged under color but still have other colorless particles, including the DM candidate. Although it is not the focus of this article, we want to remark that this is one of the more interesting aspects about these models and it embodies one of the important differences between one-loop and two-loop scenarios.
The results will be divided into two main sections according to the number of fields needed to complete the diagrams: those that require seven particles (section 2.1) and those with six particles (section 2.2). The theoretical difference between models in sections 2.1 and 2.2 is that models in the former have solely Yukawa and cubic interactions, while in the latter a scalar quartic coupling is always present.
While some of the fields used to complete the diagrams can be SM, there is an exception: the Higgs, a scalar with quantum numbers (2, 1) under SU(2) L × U(1) Y , cannot be present in the internal lines since it will fracture the two-loop diagram and generate a one-loop mass contribution after EW symmetry breaking, this effect has been carefully explained in section 3.3 of ref. [74]. Furthermore, all new scalars are taken to be inert, i.e., vanishing VEV, to avoid as well one-loop diagrams from appearing after SU(2) L breaking. The SM fields that might be used to complete the loop are then the fermions L = (2, 1)

JHEP04(2017)148
and e R = (1, 2). The field L appears explicitly in the Weinberg operator and therefore is used as two of the external legs in the diagrams, it has a left-handed chiral structure encoded in the projector P L = (1 − γ 5 )/2. BSM fermions needed to construct the models presented lack this chiral structure and are vector-like in nature.
There are multiple ways of drawing the two-loop diagrams. We will choose a particular way for which the figure does not have any crossing lines in a plane (i.e., any vertex in the drawing will represent a term in the Lagrangian). Sketching them in this fashion will make the two-loops clear and visually explicit. Based on the drawings, we will divide each section into two categories and we will label the models in each of them by the category's number and a lowercase letter, e.g., for Category 1, with eight possible diagrams, the diagrams will be numbered from 1.a to 1.h. The diagrams are traced in the standard way, the dashed lines refer to scalars and the solid lines to fermions. The particles in the loops are generically labeled 1 as X i with i = 1 to 7 in section 2.1 and i = 1 to 6 in section 2.2.
We assume the presence of a Z 2 symmetry in all models. This symmetry, introduced with the purpose of ensuring the stability of the DM candidate, will forbid (when necessary) Lagrangian terms that can be used to construct type-I seesaw masses and, in most cases, forbids one-loop diagrams. Since we are considering only singlets and doublets of SU(2) L , we do not need to be concerned about scenarios with seesaw type-II and type-III. Additionally, since we do not allow fermion singlets with zero hypercharge and uncharged under Z 2 (i.e., right-handed neutrinos) it is not possible to write neutrino Dirac masses.
Following what was explained in the introduction, we are only interested in models that explicitly have an DM particle functioning as an internal messenger in the loops. Therefore we have excluded all of the models that, while theoretically allow for a genuine two-loop neutrino mass and have a field that contains a DM candidate, do not have an electrically neutral particle in the loop. This can happen, for example, when one decomposes the fields into their corresponding components and the only contribution to the neutrino mass comes from diagrams where all of the particles inside the loops have an electric charge.
In order to build the two-loop diagrams we need to specify the SU(2) L × U(1) Y × Z 2 quantum numbers of each particle X i . Models of Categories 1, 2 and 3 have four different SU(2) L possibilities while models of Category 4 have just two, we give these possibilities labeled by roman numerals (i -iv ). The hypercharge Y (with normalization Q = I 3 + Y /2) is encoded in the parameters α and β, and has integer values −4 < α, β < 4. While models can be built using higher hypercharges, we will not be using those since they would involve particles with "exotic" electric charge such as Q = 3 or particles with higher charges that spoil the perturbativity of α EW .
In what Z 2 is concerned, the SM particles transform as even (+) while the BSM particles transform as either even (+) or odd (−). This ambivalence results from multiple possibilities to assign the Z 2 charges while allowing for all the interactions in the loop. Ignoring the trivial assignments we have a total of 3 different combinations for the Z 2 charges, identifiable by capital letters (A, B, C).

JHEP04(2017)148
In each section we present the relevant Lagrangian terms for the vertices of the diagram as well as the neutrino mass matrix that can be calculated after integration. We will not be concerned either with the exact composition of the Yukawa matrices or with the possibility to recreate neutrino mixing since one can always use the Casas-Ibarra parametrization [75] to compute it.
The models are presented in the form of tables; for each Class (diagram) we give the hypercharge parameters α and β, the DM candidate appearing in the model, the minimum number of BSM fields needed to complete the diagram (considering that fields with the same quantum numbers are the same and taking into account possible SM fields inside the loop). We also present the number of BSM scalar doublets, with hypercharge Y = 1, 3 (2 1 S , 2 3 S ), and finally we mark if the models contain a doubly charged fermion and/or scalar without distinguishing their nature, i.e., singlets with hypercharge Y = 4 or doublets with Y = 3.

Seven-particle models
In this section we discuss all two-loop models for which seven particles are needed. The Lagrangian for these models will only have Yukawa and scalar cubic interactions (although some of these will only have the first type).
The diagrams can be split in two ways: those for which all external legs couple to outside particles in the loop -the planar diagrams or Category 1 (figure 1) and those in which two external legs couple to an internal line in the loop -the non-planar diagrams or Category 2 (figure 2). One has in total 15 different diagrams built with seven particles, 8 in Category 1 and 7 in Category 2.
There are three different Z 2 charge assignments for each category which are given in the upper right corner of tables 1 and 9.

Category 1
Category 1 is composed of 8 diagrams, shown in figure 1. Each diagram defines a Class; there are two classes with six Yukawas and no cubic interactions (1.d and 1.g), two diagrams for which there are five Yukawas and one cubic interaction (1.a and 1.e), three diagrams with four Yukawas and two cubic interactions (1.b, 1.f and 1.h) and just one diagram that presents three Yukawas and three cubic interactions (1.c).
In table 1 we show the different assignments (i -iv ) for the SU(2) L (upper left), the hypercharge (bottom) and the Z 2 charge (upper right) assignments. Notice that the hypercharge assignment is diagram dependent, i.e., depends on the geometry of the diagram, and hence one has three different assignments (one for each line of diagrams in figure 1) Class 1.a There are four models in Class 1.a for which it is possible to generate neutrino masses and have stable DM candidates. In all the cases the DM candidate is a scalar: 1 0 S in one model and 2 1 S in the remaining ones. There is an interesting model for which only four BSM particles are required however this model does not have any doubly charged particle possible testing at LHC. The results are presented in table 2. X 1 X 2 X 3 X 4 X 5 X 6 X 7 i 1 1 1 2 1 1 2 ii 1 2 2 1 2 1 2 iii 2 1 2 2 1 2 1 iv 2 2 1 1 2 2 1 The most generic Lagrangian is given by,   The mass matrix is then given as,

JHEP04(2017)148
The masses m AB and the angles θ AB refer to the mass eigenstate of a particle, X AB , that got mixed between particles X A and X B and its mixing angle, respectively. The indices α and β indicate the different components of the mass eigenstates. Hereinafter we will use this terminology for all the classes. The integrals of I 1, q 2 , k 2 , (q+k) 2 χ type are shown in the appendix and are computed explicitly in [37]. We should consider for this class χ ≡ {S 6 , F 17 α , S 3 , S 5 , F 24β }. Class 1.b This class has five models for which one can generate neutrino masses. All models have scalar DM candidates stabilized by Z 2 symmetry whose assignment is given in table 1. The most appealing model needs five BSM fields to be drawn and presents a DM candidate that is a mixture between 1 0 S and 2 1 S , referred as * in the table of results,  table 3. The Lagrangian for this class is for which the neutrino mass matrix is given by  Class 1.c There are four models generating neutrino masses via diagram 1.c: two of them lead to scalar DM and the other two to fermion singlet DM candidates. The models in the first and last column of table 4 are similar in all respects other than the Z 2 assignment. The Lagrangian for this class is given by

JHEP04(2017)148
Class 1.d This class presents fourteen realizations (see results in table 5) for the two-loop neutrino mass generation. The three Z 2 charge assignments as well as the four possible SU(2) L assignments are present. There is no fermion DM candidate coming from this class JHEP04(2017)148 and the mass matrix Class 1.e There are three models that can generate neutrino masses via Class 1.e diagrams, solutions in table 6 below. In all the models, the DM candidate is a fermionic singlet. The Lagrangian is written as follows 9) and the neutrino mass matrix is given by The neutrino mass matrix is given by Class 1.g Diagram 1.g leads to four models divided by two Z 2 charge assignments (B and C) and the SU(2) L assignment iv. The DM can be both fermion or scalar singlet. The solutions are shown in table 8. The Lagrangian is given by while the neutrino mass matrix by Class 1.h This class has no models that allow for a DM messenger inside the two-loop diagram with the discrete symmetry considered. We give the Lagrangian and the mass matrix just for completeness:

Category 2
This category is composed of 7 diagrams, shown in figure 2. As in the previous category, the diagrams define a class and they have the following number of Yukawa/cubic scalar interactions: six/zero (diagrams 2.b, 2.d and 2.f), five/one (diagram 2.g), four/two (diagram 2.e) and three/three (diagrams 2.a and 2.c). The main difference between Category 1 and Category 2 is the geometry: in Category 2 the two external legs are placed inside the loops making these non-planar diagrams.
The SU(2) L , Z 2 and U(1) Y assignments are given, respectively, in the upper left, upper right and bottom of table 9. This category has four different hypercharge assignments which are related with the way how the SM fields are organized in the loop.

Class 2.b
For this class, we found seven models three of them with 1 0 S as a DM candidate. There are two similar sets of solutions for SU(2) L assignment iii but with different Z 2 charges. See the results on table 11. The Lagrangian of this class is given by the neutrino mass matrix is then with χ ≡ {S 2 , F 4 , F 345α , F 1 , S 7 }.

JHEP04(2017)148
for which one extracts  this diagram is that leads to the neutrino mass matrix given by with ρ ≡ {F 24 α , S 3 , S 56 β , S 7 , F 1 }.
Class 2.e There are fourteen models leading to neutrino mass via diagram 2.e whose the Lagrangian is The matrix M ν is given by for ρ ≡ {S 24 α , S 6 , F 35 β , S 7 , F 1 }, Among the solutions, presented in table 14, there are models with a scalar singlet as the DM candidate and solutions where the DM is the mixture between 1 0 S and 2 0 S . All the models apart from two have at least one doubly-charged particle.  figure 2.

JHEP04(2017)148
where the neutrino mass M ν is given by Class 2.g The last class of this section contains five models with either scalar singlets or fermion singlets as DM candidates. All the models have doubly-charged fermions and scalars. The solutions given in table 16 have the following Lagrangian,   and the neutrino mass matrix is given by for η ≡ {S 24 α , F 3 , S 5 , F 6 , F 17 β }.

Six-particle models
In this section we separate the possible diagrams into two categories based on the quartic coupling they have: diagrams that have a coupling of H S i S j S k type (section 2.2.1) and those with the coupling H H S i S j (section 2.2.2), where S i,j,k refer to BSM scalars that can be singlets or doublets.

Category 3
This category is defined by the presence of a term of the form H S i S j S k in the scalar potential. The diagrams within this category are shown in figure 3 while the assignments for SU(2) L , Z 2 and U(1) Y are given, respectively, in the upper left, upper right and bottom charts of table 17 (Class 3.a and 3.b) and table 18 (Class 3.c). Table 17. SU(2) L (upper left), Z 2 (upper right) and U(1) Y (bottom) assignments for particles X 1 to X 6 of diagrams 3.a and 3.b in figure 3. Class 3.a We compile in table 19 the fifteen models found for this class. Among the solutions seven have no doubly-charged particle. The DM candidate is, in all models, a scalar. The Lagrangian is given by

JHEP04(2017)148
leading to the neutrino mass matrix where χ ≡ {S 5 , S 24 α , S 3 , F 1 , S 6 }.    figure 3. For each set of SU(2) L quantum numbers (roman numerals) and Z 2 assignments (capital letters A, B, C) we give the hypercharge in terms of α and β parameters, the DM candidates, the number of beyond SM particles (#), the number of 2 1 S and 2 3 S present as well as the existence (or not) of doubly-charged fermions (F ++ ) and doubly-charged scalars (S ++ ) particles. The symbol (×) means present (absent) while the symbol * , if present in the DM row, means mixing between 1 0 S and 2 1 S . The SU(2) L , Z 2 and U(1) Y assignments are given in table 17. In red we marked models with at least one even-charged Z 2 scalar doublet with hypercharge Y = 1 that generically can create FCNCs.  which have the following neutrino mass matrix,

JHEP04(2017)148
Class 3.c From the seven models of Class 3.c only one has a pure scalar singlet as DM candidate being the remaining ones a mixture between 1 0 S and 2 1 S . The result of this class with the Lagrangian given by   X 1 X 2 X 3 X 4 X 5 X 6 ii 1 2 2 1 2 2 iii 2 1 2 2 2 1 Table 22. SU(2) L (left) and U(1) Y (right) assignments for particles X 1 to X 6 of diagrams in figure 4.

Category 4
This category contains two different classes for which a scalar term H H S i S j is present. The SU(2) L and U(1) Y assignments are given in  Table 23. Z 2 assignments for particles X 1 to X 6 of diagrams 4.a (left) and 4.b (right) of figure 4.  figure 4. For each set of SU(2) L quantum numbers (roman numerals) and Z 2 assignments (capital letters A, B, C) we give the hypercharge in terms of α and β parameters, the DM candidates, the number of beyond SM particles (#), the number of 2 1 S and 2 3 S present as well as the existence (or not) of doubly-charged fermions (F ++ ) and doubly-charged scalars (S ++ ) particles. The symbol (×) means present (absent) while the symbol * , if present in the DM row, means mixing between 1 0 S and 2 1 S . The SU(2) L , Z 2 and U(1) Y assignments are given in table 22. In red we marked models with at least one even-charged Z 2 scalar doublet with hypercharge Y = 1 that generically can create FCNCs. requires only three BSM fields to be complete. The interactions in this class are given by which lead to the following neutrino mass matrix   The mass matrix is given by

Some remarks on the models phenomenology
In this section we briefly comment on some phenomenological aspects of the models presented above.

The scalar sector
Depending on the specific diagram, two-loop radiative neutrino masses require six or seven internal messengers to close the loop. More concretely, the number of additional scalars that need to be added are controlled by the number of Yukawa interactions and cubic scalar terms of each diagram.  1.b, 2.a and 2.c) employ five additional scalars. By inspection of the models one can see that the simplest scalar sector is for a model for which one has the SM Higgs plus two singlets, while the more complicated scalar sector will be for a model with the SM Higgs, plus three doublets and two singlets. Although each model has a specific scalar potential, it will be a combination of the SM Higgs potential plus the contribution coming from singlets with hypercharges Y = 0, 2, 4 and/or doublets with hypercharges Y = 1, 3. The scalar potential for the pure SM is given by

JHEP04(2017)148
while the generic scalar potential for SM plus scalar doublets is given as where a and b represent different fields. φ refers to a BSM scalar SU(2) L doublet with Y = 1 and χ is a SU(2) L doublet with Y = 3. For a or b = 1, φ 1 = H, is the SM Higgs.
In the case of doublets plus singlets (S i ), the generic scalar potential has the extra terms, where, in addition to the labels in the last equation, we have the indices i and j that represent different singlets with the same hypercharge. Finally, the potential with the cubic terms is given by Note that the scalar potential given above is very generic and each specific model contains different terms. This is not only because of the difference between the particle content of each model but also due to the presence of the Z 2 symmetry that forbids certain terms from appearing. The expressions given above in eqs. (3.1)-(3.4) are generic and need to be taken together with the Z 2 charge assignments.
When writing the potential we have considered H as the only scalar that will acquire a VEV, i.e., the BSM scalars are all inert. 2 Although in most cases H is, in fact, the SM Higgs, this is not the case for all models, the reason for this has to do with the fact that in some models it mixes with other scalars. It is very common in the literature to add an ad hoc Z 2 symmetry which forbids the terms that mix the SM Higgs and other doublets (λ Hi 4 and m Hi terms). In our case the Z 2 charges are fixed by the two-loop diagrams, this allows some models to have a scalar doublet φ with the same Z 2 charge as the SM Higgs and hence mixing them. This implies that there could be cases where the SM Higgs is the lightest of a two-eigenstate system {H 0 , φ 0 }.
Note that one could also think of models where µ 1 mixes a hypercharge Y =1 scalar doublet, φ, with a neutral scalar singlet, S 0 . In this generic case, the SM Higgs is now the lightest of the system {H 0 , φ 0 , S 0 }. Given the Z 2 assignments, none of the models in this article have these three eigenstates at once and we do not need to be concerned about these cases.
In section 2 we showed for each model the number of scalar doublets needed. In twenty four of the models the scalar sector only consists of the SM Higgs plus singlets. Also, it can be seen that some of the models that were presented have a scalar potential that has been studied in multi-Higgs-doublet-model (MHDM), more exactly there are fifteen 2HDM, nine 3HDM and five 4HDM. These models are a simple extension of the SM JHEP04(2017)148 obtained by adding scalar doublets and they have been thoroughly studied in the literature, see e.g. [76]. The most common MHDM contains only SU(2) L doublets with hypercharge Y = 1, nevertheless, we can see that most models exhibit more of a modified structure since they include also doublets with hypercharge Y = 3, that have been studied a little in literature [77][78][79][80][81] and still require more analysis. The addition of inert singlets does not affect the phenomenology [82] of MHDM models.

Flavor changing neutral currents
In the simplest of these cases, when a second Higgs scalar is added to the SM, then the two Higgs doublets will allow SM fermions to couple to up-type, down-type fermions and leptons. Therefore, generic Yukawa couplings are not allowed without inducing dangerous contributions to flavor-changing neutral currents (FCNCs). In other words, in the NHDM generic case the are terms in the Yukawa Lagrangian that will create unwanted FCNCs [76,[83][84][85].
This will be translated into our models, that unless some other physics is involved, there should be a restriction for the models to contain an even-charged Z 2 scalar doublets with hypercharge Y = 1 (other than the SM Higgs) [86]. For this reason we kept all such models marked in red in the results shown before. Luckily models in Category 1 and 2 never have this field. The reason is that whenever this field appears it is always accompanied by fields that allow for one-loop realizations. Nevertheless, this is not the case for Category 3 and 4 where only a limited number of such models appear.
There are many possibilities of how to avoid FCNCs in these multi-Higgs models, the most common is to chose the Z 2 of the scalars to forbid in the Yukawa terms the dangerous contribution [87].
However, the models in this paper lack this freedom, since the Z 2 charges are assign uniquely based on what the two loops diagrams permit. Other option is to chose the VEV alignment of the scalars, but once again in all the models presented only the SM Higgs gets a VEV, and therefore this is not a solution to eliminate FCNCs. Nevertheless, there are options of how to eliminate FCNCs, for example, hierachical, Yukawa-alinged, VEV-aligned, minimal flavor violation models [88][89][90][91][92][93][94].

Addition on triplets
Interestingly, none of the models presented have both fermion and scalar as possible DM candidates at the same time. This statement changes if we relax the condition of having just singlets and doublets and introduce SU(2) L triplets. In fact, all the diagrams presented can be extended to include SU(2) L triplets, to do so, one needs to trade in each vertex of the type doublet-doublet-singlet, the singlet for a triplet and in the singlet-singlet-singlet vertices, two singlets for two triplets or all singles for triplets. Such an extension allows us to find solutions that were previously not viable due to our restrictions, particularly those that can create neutrino masses but have no DM particle inside the loop or those for which the external legs were not all neutral. However, the introduction of triplets escalates the number of allowed models. As an example, we use diagram 1.d; without triplets there are 14 possible models, with them we found 162. Even if other diagrams do not behave as JHEP04(2017)148 dramatically as 1.d, it is clear that the number of models, if one wishes to include triplets, would be too great to be presented in a concise manner. One needs to be aware that the addition of SU(2) L triplets to the models needs to be made carefully since they might induce neutrino masses at tree level via type-II and/or type-III seesaw mechanism.

Doubly-charged particles
In section 2 we showed all models that can generate genuine two-loop neutrino masses when considering the condition of adding just singlets and doublets to the SM. This condition was implemented, mainly, to keep the amount of presented solutions to a manageable number but also to restrain how "exotic" the new fields are. Nevertheless, most of these models do contain at least one peculiar particle characterized by having an electric charge Q = 2.
The total number of models presented is 141, out of those 118 have the presence of at least one doubly-charged particle: 15 have doubly-charged fermions, 33 doubly-charged scalars and 70 have both. Given our construction these particles can only be singlets with hypercharge Y = 4 or the component of a doublet with Y = 3.
The connection between neutrino masses and doubly-charged particles is not new and it has been studied since the Babu-Zee model [95] which includes a doubly-charged scalar singlet.
The possibility of detecting such a particle at colliders has been discussed in the literature [95][96][97][98][99][100] and several searches are being done at the CMS [101] and ATLAS experiments [102]. Because they have not been found, these experiments have set bounds on their masses. For example, a thorough study has been done in ref. [103] of possible LHC discovery of doubly-charged particles, considering spins 0, 1/2 and 1 together with different SU(2) L assignments, singlets, doublets and triplets. One of the interesting characteristic for these searches is that the SM background is suppressed making these signatures easy to analyze. Any positive signature would be excellent to discriminate between many of the models.
Doubly-charged particles coming from singlets or doublets will have different decays and therefore different signatures. There are multiple channels for production of doublycharged particles, but the production via s-channel gauge-boson exchange in the one with the biggest branching ratio: In eq. (3.5a) X and Y represent two different particles from the same doublet while in eq. (3.5b) X represents either a singlet or a doublet. One of the possible phenomenological differences between singlets and doublets is that singlets can only decay through the second channel, eq. (3.5b), given they do not have W -boson interactions.
Once a doubly-charged particle is produced its subsequent decay depends on the specificity of the model, namely on the nature of the particle, fermion or scalar. Therefore the final collider signatures and their corresponding branching ratios would be able to distinguish between singlets and doublets, and fermions or scalars [77][78][79]. In particular, the shape of the invariant mass distribution m for the final states will allow to distinguish JHEP04(2017)148 between doubly-charged fermion or scalars in the initial state, while the scalar will exhibit a sharp peak around its original mass, the fermion will not display this peak [80].

Parameter space and scale
We will show that the atmospheric neutrino scale m atm = 0.05 eV can be obtained for a benchmark point where we keep the numerical value of the Yukawas perturbative, Y ν < √ 4 π. Furthermore the value of the cubic parameters, µ, is restricted by the one-loop induced quartic coupling (λ eff < 1) which is about the mass of the Higgs times a constant that depends on the masses hierarchy. For benchmark purposes we will use µ ∼ 4m h [95] where we consider the mass of all fields heavier than m h . In this exercise we will determine the approximate scale of the messenger particles inside the loops.
After the Higgs acquires a VEV, v = 174 GeV, the internal propagators connected to the Higgs leg will mix and that will reduce the number of internal messengers in the loops. All diagrams will have five messengers. Diagrams 1.g, 1.h, 2.a and 2.b are special in the sense that the mixture creates three-state eigenvalues. Diagrams of Category 3 have an H S i S j S k vertex that does not create mixing, leaving only one vertex that does so (and reducing the number of messengers from 6 to 5). Diagrams of Category 4 have an H H S i S j vertex that creates the mixing but since it is the only one that does so (the same as in Category 3) the internal messengers only get reduced by one unit. This has the consequence that most of these diagrams have the same generic integrals, with the exception of 1.g, 1.h, 2.g and 4.b (one can see that this is due to the position and composition of the external legs, compared to the other diagrams of the same category), although the integrals of these diagrams are still similar to the aforementioned.
We will use the results from ref. [37], in it, all of the two-loop integrals needed for the masses were calculated in a model-independent fashion.
When calculating the neutrino masses we will always use the mass-eigenstate basis, in this basis these models will always have particles that get mixed by the Higgs in the external legs. Depending on the model there are two types of mixing angles, those that are proportional to a Yukawa, Y (these angles mix two fermions), and the angles that mix two scalars (these are proportional to a cubic coupling µ), where the indices 1, 2 refer to the two different mass eigenstates. It is convenient to introduce a variable that parametrizes the degeneracy between two mass parameters, M 2 = deg × M 1 . Also, the scale of the new fields added is such that 2Y v, generically speaking the value of the mixing angles will be very small, unless we consider that the mass eigenstates are quasidegenerate.
To illustrate the scale of the new physics required in two-loop neutrino mass models, we use class 1.e as an example. To keep the analysis manageable, we will use a simplified parameter space, where all Yukawas are set equal to each other, the cubic couplings are all equal and have a value of µ = 500 GeV (this being a limit case, as explained above). The mass of the particle in the middle of the loop, i.e., the one with no interactions with JHEP04(2017)148 charged leptons (which in most cases is the DM candidate) will be set to M DM = 5000 GeV, all other masses will be scaled to this one. Also, we will set both parameters deg equal to each other and use values away from the degenerate and quasidegenerate cases. Figure 5 shows the value of the generalized Yukawa for different values of the messenger masses, parametrized by the ratio squared of the mass divided by the DM mass. Values above the dashed line will require Yukawas that are non perturbative, making this region of the parameter space unacceptable. All of the lines have peaks that are caused in the calculation by accidental cancellations of parameters. One can see from the general trend of the plots that one or two of the masses can have a very high scale (higher than PeV), as long as all other are kept at a TeV scale.

Number of generations
After symmetry breaking there are two types of models: the ones where  For the first case, if there is one generation per BSM fields then Y ab is a constant C and hence the neutrino matrix is given by , in this scenario there are two possibilities, either α = γ, which implies that M ν only has one nonzero eigenvalue, or α = γ which means that M ν has two nonzero eigenvalues. In the second case the Yukawas Y ab and Y bc always connect three BSM fields and they must have the same dimensions, which translates into the fields having the same number of generations.
If it is only one, then Y ab Y bc is a constant, and similarly to the previous case the neutrino JHEP04(2017)148 mass is given by ( , which has one or two non zero eigenvalues depending if α = δ or not. This can be resumed by the fact that if both Yukawa interactions for the external legs connecting the SM charge lepton with BSM fields are the same, then one needs at least two generations of each BSM field to be able to reproduce neutrino data, or three generations if one wishes to reproduce three nonzero neutrino masses. If these Yukawa interactions are different, then one BSM generation suffices (but two generations are necessary to have all three nonzero neutrino masses).
Higher loop diagrams could be use to generate the smallest of neutrino masses and therefore reduce the number of BSM generations, but this calculation is beyond the scope of this paper.

Conclusions
In this work we provided a list of genuine two-loop models for neutrino mass generation where one of the internal messengers is a dark matter candidate. We obtained these twoloop realizations of the Weinberg operator by adding only singlets or doublets of SU(2) L to the SM. Given the bounds set by direct detection experiments, DM will be either fermionic or scalar singlets with null hypercharge, or a scalar doublet with Y = 1 (or a mixture of the latter two). Moreover, we ensured the stability of the DM candidate by the introduction of a Z 2 symmetry.
As it was shown, there are 20 different diagrams that generate genuine two-loop neutrino masses via a Weinberg operator, one of these diagrams can not accommodate simultaneously the restrictions on the fields and the existence of a DM candidate as a messenger leaving 19 different diagrams with results. Therefore we found in total 141 different sets which fill our constrains, these models are generated via a total of 19 different diagrams which we classify into two types, those with seven particles and those with six-particles, with two subsections each, called categories. The results are compiled in tables 2-8 for Category 1, tables 10-16 for Category 2, tables 19-21 for Category 3 and finally in tables 24 and 25 for Category 4.
A short discussion on some phenomenological aspects of the models has been presented in section 3, it was mentioned the possibility of these models to have FCNCs if a scalar (2, 1, +), under SU(2) L × U(1) Y × Z 2 , is present. We also mentioned how the models in this paper can be discriminated by the use of doubly-charged particle signatures. Furthermore, we included a short comment regarding the possibility of extending these models with SU(2) L triplets; this is simple to do but increases the number of models presented to an unpractical amount. We noted the fact that most models only need one generation of BSM fields to recreate neutrino data, but two generations might be necessary in some specific cases or if one wants to have three non-zero eigenvalues for the neutrino masses.
Finally, we showed that for a particular set of parameters all the models fit the atmospheric neutrino mass whose latest experimental result is placed around m atm = 0.05 eV [1]. It was not our intention to show all the parameter space, nevertheless, given the multiple free parameters of the models as well as the solutions we did find, we expect a reasonable parameter space that allows for realistic neutrino mass spectra and a viable DM candidate.

A Integrals
In this appendix we explicitly write the integrals appearing in section 2. The integrals are specified by the indices χ, ρ, η and γ.
For diagrams 1.a-1.f, 2.a-2.b, 3.a-3.b and 4.a the integral is of the form with χ ≡ {x, y, z, r, s}. For diagrams 2.c-2.f the integral is given by with ρ ≡ {x, y, z, r, s} while for diagram 2.g the integral reads as with η ≡ {x, y, z, r, s}. Finally for diagram 1.g, 1.h and 4.b it is given by with γ ≡ {x, y, z, r, s}, where ζ refers to 1, k 2 , q 2 and (q + k) 2 . The solution to these integrals is computed in a model-independent fashion explicitly in [37].

B Examples
B.1 Z 2 symmetry to forbid one-loop contributions Diagrams 1.a, 1.b and 4.b have bilateral symmetry, i.e., if one draws a vertical line right through the middle of the diagram there will be two mirror images on both sides of the lines. In these three cases it is not necessary to consider the Z 2 assignments A and B of table 1, and only one suffices. Diagrams 1.c, 1.e and 1.f also have a bilateral symmetry over an inclined line that also reduce the number of possible models. Nevertheless, that is not the case for all other models, and the specific Z 2 assignment is not arbitrary and have real physical consequences for the models.
To explain this in more detail we use an example, with no particular reason, we choose a model from diagram 1.d, with α = −2 and β = −1, as well as ii under SU(2) L assignment as given by table 1. With these quantum numbers we look at two different models: in one we assign Z 2 charge to the fields in the left loop (Z 2 assignment A), and for a second model the Z 2 charges are for the right loop (Z 2 assignment B). In table 26 we show the particle content for these models.
In figure 6 we show the Feynman diagrams for the Majorana neutrino masses. In the upper panel we show the main contribution given by a two-loop diagram 1.d, while in the JHEP04(2017)148 SU(2) L U(1) Y Z 2 (model 1) Z 2 (model 2)  bottom panel we show the two main contributions, in red we marked the interactions in the two-loop contribution that can be used to generate a one-loop. One can clearly see that if one wishes to search for all possible models then the Z 2 assignment for the fields is important and one can not just assign charges to one-loop. Hence, as can be seen by this example the Z 2 not only stabilizes a DM candidate but also helps to forbid contributions from one-loop or tree level.

B.2 A model with DM that does not contribute to neutrino mass
In section 2 it was mentioned that: "[. . . ] we have excluded all of the models that, while theoretically allow for a genuine two-loop neutrino mass and have a field that contains a DM candidate, they do not have an electric neutral particle in the loop [. . . ]". Since this idea is difficult to grasp, we will explain it by using an example. In figure 7 we show two models both generate a Majorana neutrino mass at a genuine two-loop level with a diagram of the type 1.d, the difference between them is that for the second one we have rescaled the hypercharge of one of the loops by 2. On the left we show the diagrams constructed with the fields and their SU(2) L × U(1) Y × Z 2 numbers. On the right we show the same model and diagram but with the explicit field components for each vertex of the diagram. As can be seen in both models there is a DM candidate given by a doublet scalar with hypercharge Y = 1, nevertheless in the first model the neutral component of the scalar doublet is the one that participates in the loop, while in the second model is the charged component of this doublet. Therefore, the second model does not have an explicit DM candidate contributing to the neutrino mass.