REVIEW article

Front. Phys., 04 December 2017

Sec. High-Energy and Astroparticle Physics

Volume 5 - 2017 | https://doi.org/10.3389/fphy.2017.00063

From the Trees to the Forest: A Review of Radiative Neutrino Mass Models

  • 1. School of Physics, Sun Yat-sen University, Guangzhou, China

  • 2. ARC Centre of Excellence for Particle Physics at the Terascale, School of Physics, The University of Melbourne, Melbourne, VIC, Australia

  • 3. ARC Centre of Excellence for Particle Physics at the Terascale, Department of Physics, The University of Adelaide, Adelaide, SA, Australia

  • 4. ARC Centre of Excellence for Particle Physics at the Terascale, School of Physics, The University of Sydney, Sydney, NSW, Australia

  • 5. Instituto de Física Corpuscular (CSIC-Universitat de València), Valencia, Spain

Article metrics

View details

289

Citations

10,9k

Views

4,3k

Downloads

Abstract

A plausible explanation for the lightness of neutrino masses is that neutrinos are massless at tree level, with their mass (typically Majorana) being generated radiatively at one or more loops. The new couplings, together with the suppression coming from the loop factors, imply that the new degrees of freedom cannot be too heavy (they are typically at the TeV scale). Therefore, in these models there are no large mass hierarchies and they can be tested using different searches, making their detailed phenomenological study very appealing. In particular, the new particles can be searched for at colliders and generically induce signals in lepton-flavor and lepton-number violating processes (in the case of Majorana neutrinos), which are not independent from reproducing correctly the neutrino masses and mixings. The main focus of the review is on Majorana neutrinos. We order the allowed theory space from three different perspectives: (i) using an effective operator approach to lepton number violation, (ii) by the number of loops at which the Weinberg operator is generated, (iii) within a given loop order, by the possible irreducible topologies. We also discuss in more detail some popular radiative models which involve qualitatively different features, revisiting their most important phenomenological implications. Finally, we list some promising avenues to pursue.

1. Introduction

The discovery of neutrino oscillations driven by mass mixing is one of the crowning achievements of experimental high-energy physics in recent decades. From its beginnings as the “solar neutrino problem”—a deficit of electron neutrinos from the Sun compared to the prediction of the standard solar model, an anomaly first discovered by the Homestake experiment—through the emergence of the “atmospheric neutrino problem” and its eventual confirmation by SuperKamiokande, to terrestrial verifications by long baseline and reactor neutrino experiments, the existence of non-zero and non-degenerate neutrino masses is now well established [117]. In addition, the existence of oscillations proves that the weak eigenstate neutrinos νe, νμ, and ντ are not states of definite mass themselves, but rather non-trivial, coherent superpositions of mass eigenstate fields called simply ν1, ν2, and ν3, with masses m1, m2, and m3, respectively1. The dynamical origin of neutrino mass is at present unknown, including whether neutrinos are Dirac or Majorana fermions. In the former case, neutrinos and antineutrinos are distinct and have a total of four degrees of freedom, exactly as do the charged leptons and quarks. Majorana fermions, on the other hand, are their own antiparticles, and they have just two degrees of freedom corresponding to left- and right-handed helicity. Dirac neutrinos preserve total lepton number conservation, while Majorana neutrino masses violate lepton number conservation by two units. The purpose of this review is to survey one class of possible models, where neutrino masses arise at loop order and are thus called “radiative.” Almost all of the models we examine are for the Majorana mass case. Before turning to a discussion of possible models, we should summarize the experimental data the models are trying to understand or at least accommodate.

The Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix (Uαi) [18, 19] defines the relationship between the weak and mass eigenstates, through

where α = e, μ, τ and i = 1, 2, 3. The PMNS matrix U is unitary, and may be parameterized by three (Euler) mixing angles θ12, θ23, and θ13, a CP-violating Dirac phase δ that is analogous to the phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix, and two Majorana phases α2,3 if neutrinos are Majorana fermions. The standard parametrization is

where cij ≡ cosθij and sij ≡ sinθij. The neutrino oscillation lengths are set by the ratio of squared-mass differences and energy, while the amplitudes are governed by the PMNS mixing angles and the Dirac phase. The Majorana phases do not contribute to oscillation probabilities. The angles θ12, θ23, and θ13 are sometimes referred to as the solar, atmospheric and reactor angles, respectively, because of how they were originally or primarily measured. The “solar” and “atmospheric” oscillation length parameters are, respectively,

where the distinction between the two atmospheric quantities will be discussed below.

A recent global fit [20] obtains the following 3σ ranges for the mixing angle and Δm2 parameters:

where i = 1, 2 depending on the sign of the atmospheric squared-mass difference (see Forero et al. [21] and Capozzi et al. [22] for earlier fits). The sign of has been measured because the Mikheyev-Smirnov-Wolfenstein or MSW effect [23, 24] in the Sun depends on it. The sign of the atmospheric equivalent is, however, not currently known, and is a major target for future neutrino oscillation experiments. Because of this ambiguity, there are two possible neutrino mass orderings: m1 < m2 < m3 which is called either “normal ordering” or “normal hierarchy”, and m3 < m1 < m2 which is termed “inverted.” The global fit results for the other parameters depend somewhat on which ordering is assumed. In Equations (4) and (5) we quote results that leave the ordering as undetermined. See Esteban et al. [20] for a discussion of these subtleties, but they will not be important for the rest of this review. Note that the convention is i = 1 in Equation (5) for normal ordering and i = 2 for inverted ordering.

At the 3σ level, the CP-violating phase δ can be anything. However, there is a local minimum in χ2 at δ ~ −π/2, which is tantalizing and very interesting. It hints at large CP-violation in the lepton sector, and the specific value of −π/2 is suggestive of a group theoretic origin (but beware that the definition of this phase is convention dependent). As with the mass ordering, the discovery of CP violation in neutrino oscillations is a prime goal for future experiments. One strong motivation for this is the cosmological scenario of baryogenesis via leptogenesis [25], and even if other sources of leptonic CP-violation are involved, it is important to experimentally establish the general phenomenon in the lepton sector. At present, we do not know if neutrinos are Dirac or Majorana fermions, so there is no information about the possible Majorana phases α2,3. Neutrinoless double-beta decay is sensitive to these parameters, as is standard leptogenesis.

The final parameter to discuss is the absolute neutrino mass scale. The square root of the magnitude of the atmospheric Δm2 provides a lower bound of 0.05 eV on at least one of the mass eigenvalues. Laboratory experiments performing precision measurements of the tritium beta-decay end-point spectrum currently place a direct kinematic upper bound of about 2 eV on the absolute mass scale [2628] as quantified by an “effective electron-neutrino mass” , independent of whether the mass is Dirac or Majorana, and the sensitivity of the currently running KATRIN experiment is expected to be about 0.2 eV [29]. With appropriate caution because of model dependence, cosmology now places a strong upper bound on the sum of neutrino masses of about 0.2 eV [30], with the precise number depending on exactly what data are combined. If the neutrino mass sum was much above this figure, then its effect on large-scale structure formation—washing out structure on small scales—would be strong enough to cause disagreement with observations. For Majorana masses, neutrinoless double beta-decay experiments have determined an upper bound on an effective mass defined by

of 0.15−0.33 at 90% C.L., depending on nuclear matrix element uncertainties [31]2. We can thus see that experimentally and observationally, we are closing in on a determination of the absolute mass scale.

The fact that the laboratory and cosmological bounds require the absolute neutrino mass scale to be so low strongly motivates the hypothesis that neutrinos obtain their masses in a different manner from the charged leptons and quarks. A number of approaches have been explored in the literature, with one of them being the main topic of this review: radiative neutrino mass generation. Other approaches will also be briefly commented on, to place radiative models into the overall context of possible explanations for why neutrino masses are so small.

This completes a summary of the neutrino mass and mixing data that any model, including radiative models, must explain or accommodate. As noted above, future experiments and observational programs have excellent prospects to determine the mass ordering, discover leptonic CP violation, observe neutrinoless double beta-decay (0νββ) and hence the violation of lepton number by two units, and measure the absolute neutrino mass scale. In addition, the determination of the θ23 octant—whether or not θ23 is less than or greater than π/4—is an important goal of future experiments. Before turning to a discussion of neutrino mass models, we should review some interesting experimental anomalies that may imply the existence of light sterile neutrinos3 in addition to the active flavors νe,μ,τ (see Gariazzo et al. [32], Kopp et al. [33] for phenomenological fits).

There are three anomalies. The first is > 3σ evidence from the LSND [34, 35] and MiniBooNE [36, 37] experiments of appearance in a beam, with MiniBooNE also reporting a νe signal in a νμ beam. Interpreted through a neutrino oscillation hypothesis, these results indicate an oscillation mode with a Δm2 or order 1 eV2. This cannot be accommodated with just the three known active neutrinos simultaneously with the extremely well-established solar and atmospheric modes that require much smaller Δm2 parameters. This hypothesis thus only works if there are four or more light neutrino flavors, and the additional state or states must be sterile to accord with the measured Z-boson invisible width4. The Icecube neutrino telescope has recently tested the sterile neutrino oscillation explanation of these anomalies through the zenith angle dependence of muon track signals and excludes this hypothesis at about the 99% C.L. [38].

The next two anomalies concern νe and disappearance. Nuclear reactors produce a flux that has been measured by several experiments. When compared to the most recent computation of the expected flux [39, 40], a consistent deficit of a few percent is observed, a set of results known as the “reactor anomaly” [41]. The Gallium anomaly arose from neutrino calibration source measurements by the Gallex and SAGE radiochemical solar neutrino experiments, also indicating a deficit [4245]. Both deficits are consistent with very short baseline transitions driven by eV-scale sterile neutrinos, and a significant number of experiments are underway to test the oscillation explanation. It should be noted that a recent analysis by the Daya Bay collaboration points to the problem being with the computation of the reactor flux rather than being an indication of very short baseline oscillations [46]. The key point is that if a sterile neutrino was responsible, one should observe the same deficit for all neutrinos from the reactor fuel, independent of nuclear species origin, but this was observed to not be the case. There is also a tension between the appearance and disappearance anomalies when trying to fit both with a self-consistent oscillation scheme [32, 33], and there is a cosmological challenge of devising a mechanism to prevent the active-sterile transitions from thermalizing the sterile neutrino in the early universe, as thermalization would violate the ~ 0.2 eV bound on the sum of neutrino masses.

Because the situation with the above anomalies is unclear, and there are challenges to explaining them with oscillations, this review will focus on neutrino mass models that feature just the three known light active neutrinos. If any of the above anomalies is eventually shown to be due to oscillations, then all neutrino mass models will need to be extended to incorporate light sterile neutrinos, including the radiative models that are our subject in this review.

The rest of this review is structured as follows: section 2 provides a general discussion of schemes for neutrino mass generation and attempts a classification. The structure of radiative neutrino mass models is then described in section 3. Section 4 covers phenomenological constraints and search strategies, including for cosmological observables. Detailed descriptions of specific models are then given in section 5, with the examples chosen so as to exemplify some of the different possibilities that the radiative mechanisms permit. We conclude in section 6, where we discuss some research directions for the future. Appendix gives further details on the relative contributions of the different operators to neutrino masses.

2. Schemes for neutrino masses and mixings

In this section, we survey the many different general ways that neutrinos can gain mass, and attempt a classification of at least most of the proposed schemes. As part of this, we place both the tree-level and radiative models in an overarching context—a systematic approach, if you will, or at least as systematic as we can make it. The number of different kinds of models can seem bewildering, so there is some value in understanding the broad structure of the neutrino mass “theory space.”

Under the standard model (SM) gauge group GSMSU(3)c × SU(2)L × U(1)Y, the left-handed neutrinos feature as the upper isospin component of

where on the right-hand (RH) side the first entry denotes the representation with respect to the color group SU(3)c, the second SU(2)L (weak-isospin), and the third hypercharge Y, normalized so that electric charge is given by Q = I3 + Y. In the minimal standard model, there is no way to generate non-zero neutrino masses and mixings at the renormalizable level. Dirac masses are impossible because of the absence of RH neutrinos,

as are Majorana masses because there is no scalar isospin triplet

to which the lepton bilinear could have a Yukawa coupling. Thus, the family-lepton numbers Le, Lμ and Lτ are (perturbatively) conserved because of three accidental global U(1) symmetries. The discovery of neutrino oscillations means that the family-lepton number symmetries must be broken. If they are broken down to the diagonal subgroup generated by total lepton number LLe + Lμ + Lτ, then the neutrinos must be Dirac fermions. If total lepton number is also broken, then the neutrinos are either fully Majorana fermions or pseudo-Dirac5.

The question of whether neutrinos are Dirac or Majorana (or possibly pseudo-Dirac) is one of the great unknowns. The answer is vital for model building, as well as for some aspects of phenomenology. If neutrinos are Majorana, then it is not necessary to add RH neutrinos to the SM particle content. In fact, many of the radiative models we shall review below do not feature them. If RH neutrinos do not exist, then a possible deep justification could be SU(5) grand unification, which is content with a structure per family6. But another logical possibility, motivated by quark-lepton symmetry and SO(10) grand unification, is that RH neutrinos exist but have large (SM gauge invariant) Majorana masses, leading to the extremely well-known type-I seesaw model [4751]. On the other hand, if neutrinos are Dirac, then RH neutrinos that are singlets under the SM gauge group, as per Equation (8), are mandatory and they must not have Majorana masses even though such terms are SM gauge invariant and renormalizable. Thus, at the SM level, something like total lepton-number conservation must be imposed by hand. Most of the radiative models we shall discuss lead to Majorana neutrinos, though we shall also briefly review the few radiative Dirac models that have been proposed.

The choice of Dirac or Majorana is thus a really important step in model building. It is perhaps fair to say that theoretical prejudice, as judged by number of papers, favors the Majorana possibility. There are a couple of reasons for this. One is simply that Majorana fermions are permitted by the Poincaré group, so it might be puzzling if they were never realized in nature, and the fact is that they constitute the simplest spinorial representation. (Recall that a Dirac fermion is equivalent to two CP-conjugate, degenerate Majorana fermions). Another was already discussed above: even if RH neutrinos exist, at the SM level they can have gauge-invariant Majorana masses, leading to Majorana mass eigenstates overall. Yet another reason is a connection between Majorana masses and an approach to understanding electric charge quantization using classical constraints and gauge anomaly cancellation [52, 53]. Nevertheless, theoretical prejudice or popularity in the literature is not necessarily a reliable guide to how nature actually is, so the Dirac possibility should be given due consideration.

2.1. Dirac neutrino schemes

The simplest way to obtain Dirac neutrinos is by copying the way the charged-fermions gain mass. Right-handed neutrinos are added to the SM particle content, producing the gauge-invariant, renormalizable Yukawa term

where the Higgs doublet H transforms as (1, 2, 1/2) with . The Dirac neutrino mass matrix is then

To accommodate the O(0.1) eV neutrino mass scale, one simply takes . The price to pay for this simple and obvious model is a set of tiny dimensionless parameters, some six or seven orders of magnitude smaller than the next smallest Yukawa coupling constant (that for the electron), and smaller even than the value a fine-tuned θQCD needs to be from the upper bound on the neutron electric-dipole moment. This is of course logically possible, and it is also technically natural in the 't Hooft sense [54] because taking yν to zero increases the symmetry of the theory. Nevertheless, it seems unsatisfactory to most people. The really tiny neutrino masses strongly suggest that the generation of neutrino mass proceeds in some different, less obvious manner, one that provides a rationale for why the masses are so small. As well as the Dirac vs. Majorana question, the explanation of the tiny masses has dominated model-building efforts in the literature.

So, how may one produce very light Dirac neutrinos? We highlight three possibilities, but there may be others: (i) a Dirac seesaw mechanism, (ii) radiative models, and (iii) extra-dimensional theories.

2.1.1. Dirac seesaw mechanism

In addition to the νL that resides inside the doublet L, and the standard RH neutrino of Equation (8), we introduce a vector-like heavy neutral fermion NL,R ~ (1, 1, 0) and impose total lepton-number conservation with νL,R and NL,R assigned lepton numbers of 1. In addition, we impose a Z2 discrete symmetry under which νR and a new gauge-singlet real scalar S are odd, with all other fields even. With these imposed symmetries, the most general Yukawa and fermion bare mass terms are

leading to the neutral-fermion mass matrix

where

We now postulate the hierarchy mLmRMN on the justification that the bare mass term has no natural scale so could be very high, and that the symmetry breaking scale of the new, imposed Z2 should be higher than the electroweak scale. The light neutrino mass eigenvalue is thus

and the eigenvector is dominated by the νL admixture so does not violate weak universality bounds. The inverse relationship of the light neutrino mass with the large mass MN is the seesaw effect, with the postulated small parameter mR/MN causing mν to be much smaller than the electroweak-scale mass mL. The above structure is the minimal one necessary to illustrate the Dirac seesaw mechanism (and has a cosmological domain wall problem because of the spontaneously broken Z2), but the most elegant implementation is in the left–right symmetric model [55]. Under the extended electroweak gauge group SU(2)L×SU(2)R×U(1)B−L, the RH neutrino sits in an SU(2)R doublet with B−L = −1, while NL,R remains as gauge singlets. The scalars are a left–right symmetric pair of doublets HL,R with B−L = 1. The usual scalar bidoublet is not introduced. The Z2 symmetry is then a subgroup of SU(2)R, and S is embedded in the RH scalar doublet. The mass and symmetry breaking hierarchy is then . The absence of the bidoublet ensures the zero in the top-left entry of the mass matrix7. Several tree-level Dirac neutrino mass models have been discussed in Ma and Popov [56]: The SM singlet Dirac fermion NL + NR can be obviously replaced by an electroweak triplet. Alternatively a neutrinophilic two Higgs doublet model [57, 58] is an attractive possibility to obtain small Dirac neutrino masses.

2.1.2. Radiative Dirac schemes

A generalization of the symmetry structure of the Z2 Dirac seesaw model discussed above provides us with one perspective on the construction of radiative Dirac neutrino mass models. A basic structural issue with such models is the prevention of the tree-level term generated by the renormalizable Yukawa interaction of Equation (10). Some new symmetry must be imposed that forbids that term, but that symmetry must also be spontaneously or softly broken in such a way that an effective operator is produced. In the case of radiative models, this must be made to happen at loop order. One obvious possibility is to demand that “RH neutrino number” is conserved, meaning that invariance under

with all other SM fields as singlets, is imposed. One may then introduce a complex scalar ρ that transforms, for example, as

whose non-zero expectation value spontaneously breaks the symmetry. The effective operator

produced by integrating out new physics at mass scale Λ, is both SM gauge invariant and invariant under the imposed symmetry 8. It generates a neutrino Dirac mass of order

which will be small compared to the weak scale when . If this operator is “opened up”—derived from an underlying renormalizable or ultraviolet (UV) complete theory—at loop-level, then a radiative neutrino Dirac-mass model is produced. Note that in a loop-level completion, the parameter 1/Λn depends on powers of renormalizable coupling constants and a 1/16π2 per loop as well as the actual masses of new, exotic massive particles. See Ma and Popov [56] for a recent systematic study of 1-loop models based on this kind of idea. Note that the Dirac seesaw model discussed earlier is obtained as a truncated special case: the U(1) symmetry with n = 1 is replaced with its Z2 subgroup, the complex scalar field ρ is replaced with the real scalar field S, and the effective operator is opened up at tree-level.

Obviously, the phase part of ρ will be a massless Nambu-Goldstone boson (NGB), but its phenomenology might be acceptable because it only couples to neutrinos. If one wishes to avoid this long range force, one could find a way to make the new U(1) anomaly-free and then gauge it so that the NGB gets eaten, or one may use a discrete subgroup of the U(1) to forbid Equation (10). See Wang and Han [60] for a discussion of the Z2 case for 1-loop models that also include a dark matter candidate.

The above is simply an example of the kind of thinking that has to go into the development of a radiative Dirac neutrino model—we are not claiming it is the preferred option. To our knowledge, a thorough analysis of symmetries that can prevent a tree-level Dirac mass and thus guide the construction of complete theories has not yet been undertaken in the literature. That is one of the reasons this review will discuss Majorana models at greater length than Dirac models.

2.1.3. Extra-dimensional theories

One way or another, the effective coefficient in front of must be made small. Seesaw models achieve this by exploiting powers of a small parameter given by the ratio of symmetry breaking and/or mass scales. Radiative models augment the seesaw feature with 1/16π2 loop factors and products of perturbative coupling constants. In warped or Randall-Sundrum extra-dimensional theories [61, 62], the geometry of fermion localization in the bulk [63, 64] can lead to the suppression of Dirac neutrino masses through having a tiny overlap integral between the profile functions for the neutrino chiral components and the Higgs boson [63, 6568]. The phenomenological implications of Dirac neutrinos in extra-dimensional set-ups have been studied in De Gouvea et al. [69], where it is shown that these effects can be encoded in specific dimension-six effective operators.

One can also have a “clockwork” mechanism [70, 71] to generate exponentially suppressed Dirac masses. In the same way, it is also useful to have low-scale seesaw [72]. This mechanism can be implemented with a discrete number of new fields or via an extra spatial dimension [73].

2.2. Majorana neutrino schemes

We now come to our main subject: radiative Majorana neutrino mass generation. We also briefly review tree-level seesaw schemes, both for completeness and for the purposes of comparison and contrast to the loop-level scenarios. In the course of the discussion below, an attempt will be made to classify the different kinds of radiative models. This is a multidimensional problem: no single criterion can be singled out as definitely the most useful discriminator between models. Instead, we shall see that several overlapping considerations emerge, including ΔL = 2 effective operators, number of loops, number of Higgs doublets, nature of the massive exotic particles, whether or not there are extended symmetries and gauge bosons, distinctive phenomenology, and whether or not the models address problems or issues beyond just neutrino mass (e.g., dark matter, grand unification, …).

The main distinctive feature of Majorana neutrino mass is, of course, that it violates lepton-number conservation by two units. It is thus extremely useful to view the possibilities for the new physics responsible from a bottom-up perspective, meaning SM gauge-invariant, ΔL = 2 low-energy effective operators that are to be derived from integrating out new physics that is assumed to operate at scales higher than the electroweak. This approach permits the tree-level seesaw [4751, 7480] and radiative models to be seen from a unified perspective.

Taking the particle content of the minimal SM, it is interesting that the simplest and lowest mass-dimension effective operator one can produce is directly related to Majorana neutrino mass generation. This is the famous Weinberg operator [81]

where the SU(2) indices and Lorentz structures are suppressed (one can check that there is only one independent invariant even though there are three different ways to contract the SU(2) indices of the four doublets.). We say the singular “operator” for convenience, but it is to be understood that there are also family indices so we really have a set of operators. This is a mass dimension five operator, so enters the Lagrangian with a 1/Λ coefficient, where Λ is the scale of the new physics that violates lepton number by two units. Replacing the Higgs doublets with their vacuum expectation values (VEVs), one immediately obtains the familiar Majorana seesaw formula,

displaying the required suppression of mν with respect to the weak scale v when ϵ ≡ v/Λ ≪ 1, so that the ΔL = 2 new physics operates at a really high scale.

The Weinberg operator can be immediately generalized to the set

where the number of primes is equal to n. One obtains ever more powerful seesaw suppression,

as n increases.

The task now is to derive, from an underlying renormalizable or UV complete theory, one of the Weinberg-type operators as the leading contribution to neutrino mass. This process has come to be termed “opening up the operator.” The choices one makes about which operator (what value of n) is to dominate and how it is to be opened up determine the type of theory one obtains. Here are some possible choices:

  • Open up O1 at tree-level using only exotic massive fermions and scalars as the new physics.

  • Open up O1 at j-loop level using heavy exotics only.

  • Open up O1 at j-loop level using both light SM particles and heavy exotics.

  • Open up at tree-level using heavy exotics only.

  • Open up at j-loop level using heavy exotics only.

  • Open up at j-loop level using both light SM particles and heavy exotics.

Option 1 leads, in its simplest form, precisely to the familiar type-I [4751], type-II [7479] and type-III [80] seesaw mechanisms, as we review in the next subsection. Option 2 leads to a certain kind of radiative model, to be contrasted with that arising from option 3. The difference between the two can be expressed in terms of the matching conditions used to connect an effective theory below the scale Λ of the ΔL = 2 new physics to the full theory above that scale, as outlined in Figure 1. For scenario 2, the effective Weinberg operator has a non-zero Wilson coefficient at Λ, and for all scales below that. In scenario 3, on the other hand, the Weinberg operator has a coefficient at scale Λ that is loop-suppressed compared to the Wilson coefficients of other, non-Weinberg-type ΔL = 2 operators9 at that scale, where these other operators are obtained by integrating out the heavy fields only. If the matching is performed at tree-level approximation, then the coefficient of the Weinberg operator at Λ in fact vanishes. Under renormalization group mixing, the non-zero ΔL = 2 operators will, however, generate an effective Weinberg operator as the parameters are run to scales below Λ. If the matching is performed at loop-level, then the Weinberg operator will have a non-zero coefficient at scale Λ, but it will be loop-suppressed compared to the coefficients of the relevant non-Weinberg operators. Below Λ, the Weinberg operator coefficient will, once again, receive corrections from the renormalization group running and operator mixing. Option 3 will be a major topic in this review, and it motivates the enumeration of all SM gauge-invariant ΔL = 2 operators, not just those in the Weinberg class, since the non-Weinberg operators describe the dominant ΔL = 2 processes at scale Λ. Opening up the non-Weinberg operators at tree-level then provides a systematic method of constructing a large class of theories that generate neutrino masses at loop order.

Figure 1

Figure 1

Running and matching for (radiative) Majorana neutrino masses. See Appendix for a discussion of the relative contribution of the different operators.

Options 4-6 obviously repeat the exercise, but with two more powers of ϵ which help suppress the neutrino mass. With these options, one needs to ensure that generated from the new physics dominates over O1 and all lower-dimensional operators . Option 6 is similar to 3 in that the effective theory between the weak and new physics scales contains some non-Weinberg type of ΔL = 2 operator(s) that dominate at scale Λ.

2.2.1. Tree-level seesaw mechanisms

The three familiar seesaw models may be derived in a unified way by opening up the Weinberg operator O1 at tree level in the simplest possible way, using as the heavy exotics only scalars or fermions. The available renormalizable interactions are then just of Yukawa and scalar-scalar type. The opening-up process is depicted in Figure 2. The type-I and type-III seesaw models are obtained by Yukawa coupling LH with the two possible choices of (1, 1, 0) and (1, 3, 0) fermions, both of which can have gauge-invariant bare Majorana masses. The type-II model is the unique theory obtained from Yukawa coupling the fermion bilinear to a (1, 3, 1) scalar multiplet, which in turn couples to HH, a cubic interaction term in the scalar potential10. The seesaw effect is obtained in this case by requiring a positive quadratic term for the triplet in the scalar potential, that on its own would cause the triplet's VEV to vanish, but which in combination with the cubic term induces a small VEV for it.

Figure 2

Figure 2

Minimally opening up the Weinberg operator at tree-level using either exotic massive fermions or scalars. (A) Type-I seesaw model. The massive exotic particle integrated out to produce an effective Weinberg operator at low energy is a SM gauge-singlet Majorana fermion, the right-handed neutrino νR. (B) Type-II seesaw model. The massive exotic is a (1, 3, 1) scalar Δ coupling to LL and HH. It gains a small induced VEV from the latter coupling. (C) Type-III seesaw model. The massive exotic is a (1, 3, 0) fermion Σ whose middle component mixes with the left-handed neutrino.

As is clear from Figure 2, there are two interaction vertices for all three cases, and there is only one type of exotic per case. An interesting non-minimal tree-level seesaw model realizing option 4 is obtained by allowing four vertices instead of two, and two exotic multiplets: a (1, 4, −1/2) scalar that couples to HHH and a (1, 5, 0) massive fermion that Yukawa couples to the exotic scalar quadruplet and the SM lepton doublet [8284]. The resulting model produces the generalized Weinberg operator which has mass-dimension nine. This model is a kind of hybrid of the type-II and type-III seesaw mechanisms, because it features both a small induced VEV for the quadruplet and a seesaw suppression from mixing with the fermion quintuplet.

2.2.2. Radiative schemes and their classification

As noted above, there are many different kinds of radiative neutrino mass models and there is probably no single classification scheme that is optimal for all purposes. We thus discuss a few different perspectives, some much more briefly than others. Two will be treated at length: (i) the ΔL = 2 effective operator approach, and (ii) classification by loop-order openings of the Weinberg operator.

A. Standard model ΔL = 2 effective operators. This approach can be considered as stemming from the observations made about options 3 and 6 in section 2.2: when both light SM particles and heavy exotics appear in the neutrino mass loop graph, it is useful to first consider integrating out the heavy exotics at tree-level. This produces effective ΔL = 2 operators that are of non-Weinberg type. They must be of different type, because if they were not, then the heavy exotics would produce the Weinberg operator without participation by light SM particles, leading either to a class 1 model (if O1 is produced at tree-level) or a class 2 model (if O1 is produced at loop level). An exhaustive list of gauge-invariant, non-Weinberg ΔL = 2 operators is thus needed.

Such a list was provided by Babu and Leung (BL) [85], based on the following assumptions: (i) the gauge group is that of the SM only, (ii) no internal global symmetries are imposed apart from baryon number, (iii) the external lines are SM quarks, SM leptons and a single Higgs doublet, and (iv) no operators of mass dimension higher than 11 were considered. We first comment on these assumptions. Clearly, if the gauge symmetry was extended beyond that of the SM, then some combination of effective operators might be restricted to having a single coefficient, and others might be forced to vanish, compared to the SM-gauge-group-only list. Similar observations follow for imposed global symmetries. It is sensible to impose baryon number conservation, because otherwise phenomenological constraints will force the new physics to such high scales that obtaining neutrino masses of the required magnitude (at least one at 0.05 eV) will be impossible. The case of a single Higgs doublet can readily be generalized to multiple Higgs doublets, given that the gauge quantum numbers are the same. This would obviously enrich the phenomenology of the resulting models, and if additional symmetries were also admitted, then it would change the model-building options. The point is simply that HH is invariant under all possible internal symmetries, while is not. (Admitting additional Higgs doublets is also interesting for generalized-Weinberg-operator models, because then a symmetry reason can exist for, say, being generated without also generating what would otherwise be dominant LLH1,2H1,2 operators.) The addition of non-doublet scalar multiplets into the external lines is a more serious complication. Some discussion of the possible roles of additional scalars that gain non-zero VEVs that contribute to neutrino mass generation will be given in later sections. Another restriction worth noting in the BL list is the absence of the gauge-singlet RH neutrinos. In assumption (iv), the point to highlight is the absence of SM gauge fields. Babu and Leung did actually write down the mass-dimension-7 operators containing gauge fields, and Bhattacharya and Wudka [86] further examined them. As far as we know, however, no complete analysis has been undertaken for the dimension-9 and -11 cases. Finally, it is sensible to stop at dimension 11 because at any higher order the contribution to neutrino mass will be insufficiently large. The BL list, as enumerated from O1 to O60, took operators that could be thought of as products of lower-dimension operators with the SM invariants HH and the three dimension-4 charged-fermion Yukawa terms as implicit. de Gouvea and Jenkins [87] extended their list by explicitly including the latter cases, thereby augmenting the operator count to O75.

Operators meeting all of these requirements exist at all odd mass dimensions [85, 88, 89], starting with the Weinberg operator O1 as the unique dimension-5 case (up to family indices). The dimension-7 list is as follows:

We follow the BL numbering scheme, which was based on tracking the number of fermion fields in the operator rather than the mass dimension. The operators are separated in three groups with 2, 4, and 6 fermions. Some comments now need to be made about the schematic notation and what features are suppressed. The field-string defining each operator above completely defines the flavor content of that operator. Thus L ~ (1, 2, −1/2) is the lepton doublet, Q ~ (3, 2, 1/6) is the quark doublet, ec ~ (1, 1, 1) is the isosinglet charged anti-lepton, dc ~ is the isosinglet anti-down, is the isosinglet anti-up, and H ~ (1, 2, 1/2) is the Higgs doublet. The color indices and the different possible Lorentz structures are suppressed. In general, there are a number of independent operators corresponding to each flavor-string. For the dimension-7 list, operators O3 and O4 each have two independent possibilities for the contraction of the isospin indices, as explicitly defined above, but obviously a unique color contraction. Babu and Leung specify the independent internal-index contractions, but only make general remarks on the Lorentz structures, and we shall follow suit. To assist the reader to understand the notation, we write out the above operators more completely in standard 4-component spinor notation, but for scalar and pseudoscalar Lorentz structures only and with isospin indices suppressed:

Of course, these operators feature quark and charged-lepton fields in addition to neutrinos and Higgs bosons, so they do not by themselves produce neutrino masses. The charged fermion fields have to be closed off in a loop or loops to produce a neutrino self-energy graph which then generates a Weinberg-type operator, as per options 3 and 6. In fact, using this procedure and naive dimensional analysis one can estimate their matching contribution to the Weinberg operator, as done in de Gouvea and Jenkins [87]. In addition, every dimension-7 operator in Equation (24) may be multiplied by HH to produce a dimension-9 generalization of that operator, just as is a generalization of O1. At dimension 9, there are many more operators. Six of the flavor strings feature four fermion fields and three Higgs doublets:

Note that the operators O61,66,71 are the products of O1 and the three SM Yukawa operators. Another 12 are six-fermion operators:

Although absent from the BL list another such operator is , which generates the correct neutrino mass scale only for a very low lepton-number violation scale. In case it consists entirely of the first generation SM fermions it is strongly constrained by 0νββ (generated at tree level by this operator). The large number of dimension-11 operators can be found listed in Babu and Leung [85] and de Gouvea and Jenkins [87].

de Gouvea and Jenkins [87] and Angel et al. [90] performed general analyses of diagram topologies for opening up these operators at tree-level using massive exotic scalars and either vector-like or Majorana fermion exotics, and consequently producing neutrino mass at various loop levels. The operators

can give rise to 1-loop neutrino mass models, while

can produce 2-loop models. The set

can form the basis for neutrino mass to be generated at three or more loops.

In each of these cases, one may derive an indicative upper bound on the scale of new physics from the requirement that at least one neutrino mass be at least 0.05 eV in magnitude. For example, for operators involving first generation11 quarks this bound can be estimated as follows: Operator O19, which can be opened up to give a 3-loop neutrino mass contribution, has the lowest upper bound on the new physics scale of about 1 TeV (apart from ). The highest is about 4 × 109 TeV for the 1-loop case of O4a. These estimates come from an examination of the loop contribution to neutrino mass only, and do not take into account other phenomenological constraints that will exist for each complete model. As part of that, any unknown coupling constants, such as Yukawas that involve the exotic fermions and/or scalars were set to unity. In a realistic theory, many of these constants would be expected to be less than one, which would bring the scale of new physics to lower values. In any case, one can see that the required new physics, even for 1-loop models, is typically more testable than the type-I, II, and III seesaw models. Some high loop models, as the O19 case demonstrates, have very low scales of new physics and some may even be ruled out already. At the dimension-11 operator level, so not explicitly discussed here, there are even examples which can at best produce a 5-loop neutrino mass contribution. Those models are definitely already excluded. Examples of full models that are associated with specific operators will be presented in later sections.

B. Number of loops. A complementary perspective on the spectrum of possible radiative neutrino mass models is provided by adopting the number of loops as the primary consideration rather than the type of ΔL = 2 effective operator that dominates the new physics. Equations (28–30) already form the basis for such a classification for type 3 and type 6 scenarios, but a more general analysis will also capture the type 2 and type 5 possibilities.

At j-loop order, neutrino masses are typically given by

for the O1 associated options 2 and 3, and

for the cases of options 5 and 6, where GeV, and Λ is the new-physics scale where lepton number is violated by two units. All coupling constants, and for some models also certain mass-scale ratios, are absorbed in the dimensionless coefficient C. In order to explain the atmospheric mass splitting lower bound of 0.05 eV, we obtain an upper limit on the new physics scale Λ of 105C TeV for 3-loop models and 10 C TeV for 5-loop models corresponding to the O1 cases, and 10 C1/3 TeV for the case at 3-loop order. Constraints from flavor physics severely constrain the scale of new physics and the couplings entering in C. In addition, in models which feature explicit ΔL = 2 lepton-number violation through trilinear scalar interactions, the latter cannot be arbitrarily large because otherwise they have issues with naturalness (see Herrero-García et al. [91] for the case of the Zee model) and charge/color breaking minima (see Frere et al. [92], Alvarez-Gaume et al. [93] and Casas and Dimopoulos [94] for studies in the context of supersymmetry and Herrero-Garcia et al. [95] for the case of the Zee-Babu model). Thus, apart from a few 4-loop models [9698] which compensate the loop suppression by a high multiplicity of particles in the loop, the vast majority of radiative neutrino mass models generate neutrino mass at 1-, 2-, or 3-loop level. We therefore focus on these cases.

1-loop topologies for O 1 = LLHH. The opening up of the Weinberg operator at 1-loop level has been systematically studied in Ma [99] and Bonnet et al. [100]. The authors of Bonnet et al. [100] identified 12 topologies which contribute to neutrino mass. Among all the topologies and possible Lorentz structures, topology T2 cannot be realized in a renormalizable theory. For the other topologies, the expression for neutrino mass and the possible particle content for electroweak singlet, doublet, and triplet representations is listed in the appendix of Bonnet et al. [100]. The divergent ones, T4-1-i, T4-2-ii, T4-3-ii, T5 and T6, need counter-terms to absorb the divergences, which are indeed tree-level realizations of the Weinberg operators. Furthermore, for T4-1-ii, there is no mechanism to forbid or suppress the tree-level contribution from Weinberg operator, such as extra discrete symmetry or U(1). Therefore, there are in total six topologies which generate neutrino mass via a genuine12 1-loop diagram: T1-i, T1-ii, T1-iii, T3, T4-2-i, T4-3-i, which are depicted in Figure 3. Depending on the particle content, the topologies do not rely on any additional symmetry. However, the topologies T4-x-i require a discrete Z2 symmetry in addition to demanding Majorana fermions in the loop with lepton-number conserving couplings. This is difficult to achieve in a field theory, as lepton-number is necessarily broken by neutrino mass. For example, in topology T4-2-i the scalar connected to the two Higgs doublets H is necessarily an electroweak triplet and thus its direct coupling to two lepton doublets L is unavoidable. This coupling induces a type-II seesaw tree-level contribution to neutrino mass. Similar arguments hold for the other topologies T4-x-i.

Figure 3

Figure 3

Feynman diagram topologies for 1-loop radiative neutrino mass generation with the Weinberg operator O1 = LLHH. Dashed lines could be scalars or gauge bosons if allowed.

1-loop topologies for . A similar analysis has been performed for 1-loop topologies that give rise to the dimension-7 generalized Weinberg operator [101]. Of the 48 possible topologies, only the eight displayed in Figure 4 are relevant for genuine 1-loop models. For specific cases, not all of these eight diagrams will be realized. The three-point vertices can be Yukawa, gauge or cubic scalar interactions, while the four-point vertices only contain scalar and gauge bosons.

Figure 4

Figure 4

Topologies that can give rise to genuine 1-loop openings of the dimension-7 Weinberg operator .

2-loop topologies for O 1 = LLHH. A systematic analysis of 2-loop openings of O1 was performed in Aristizabal Sierra et al. [102]. Figure 5 displays the topologies identified in this study as able to contribute to genuine 2-loop models. There are additional 2-loop diagrams – that were termed “class II” – that have the form of one of the 1-loop topologies of Figure 3 with one the vertices expanded into a 1-loop subgraph. They remark the class II topologies may be useful for justifying why a certain vertex has an unusually small magnitude.

Figure 5

Figure 5

Topologies for genuine 2-loop completions of the Weinberg operator O1 = LLHH. Solid dots denote interaction vertices. Crossed lines without a dot at the intersection denote a non-planar configuration.

C. Other considerations. We now briefly survey other perspectives on classifying or discriminating between neutrino mass models.

One suggested criterion is complexity [103]. While recognizing that sometimes nature appears to favor minimal possibilities (in an Occam's razor approach), and at other times not (e.g., the old problem of why there are three families), it does make sense to rank neutrino mass models on some sensible measure of how complex they are. Law and McDonald [103] proposes a hierarchy based on (i) whether or not the model relies on the imposition of ad hoc symmetries, (ii) the number of exotic multiplets required, and (iii) the number of new parameters. Interestingly, they construct radiative models that are even simpler, on the basis of these criteria, than the 1-loop Zee-Wolfenstein model [104, 105]. However, like the Zee-Wolfenstein model, while these models generate non-zero neutrino masses, they fail phenomenologically. Thus, we must conclude that if nature utilizes the radiative mechanism, it will be non-minimal.

Another consideration for Majorana mass models is the important phenomenological connection to 0νββ decay [106108]. Just as Majorana neutrino mass models may be systematically constructed through opening up ΔL = 2 effective operators, models for 0νββ decay can be analysed by opening up the ūūddēē family of operators. The neutrino mass and 0νββ decay considerations are of course connected, but the nature of the relationship is model-dependent. An interesting situation would emerge in a hypothetical future where 0νββ decay is observed, but the standard Majorana neutrino exchange contribution through mββ is contradicted by, for example, cosmological upper bounds on the absolute neutrino mass scale. That would point to a non-minimal framework, which may be connected with radiative neutrino mass generation.

A further interesting aspect is the existence or otherwise of a deep theoretical reason for a given radiative model. At first sight, each such model looks random. However, some of them can be connected with, for example, grand unified theories (GUTs). One simple point to make is that exotics, such as scalar leptoquarks, that often feature in radiative models can be components of higher-dimension multiplets of SU(5) and SO(10). Also, by contributing to renormalization group running, some of them can assist with gauge coupling constant unification [109]. If they are to be light enough to play these roles, while other exotics within the multiplets have, for example, GUT-scale masses, then we face a similar issue to the famous doublet-triplet splitting problem. Nevertheless, this is a starting point for investigating the possible deeper origin of some of the required exotics. Another interesting GUT-related matter was analysed in depth in de Gouvêa et al. [88]. A necessary condition for a ΔL = 2 operator of a certain mass dimension to be consistent with a GUT origin is that it occurs as a term in an effective operator of the same mass dimension derived with grand unified gauge invariance imposed. For example, the dimension-7 operator O3a from Equation (24) does not appear as a component in any SU(5) operator of the same dimension. On the other hand, other SM operators are embedded in the same GUT operator, with only one of them being able of giving the dominant contribution to neutrino masses. In addition to the question of the mere existence of SM-level operators in GUT decompositions, grand unification also imposes relations between SM-level operators, including some that violate baryon number and generate B−L violating nucleon decays and/or neutron-antineutron oscillations, leading to additional constraints. In the end, the authors of de Gouvêa et al. [88] conclude that only a small subset of SM ΔL = 2 operators are consistent with grand unification.

Another strategy for uncovering a deeper origin for a radiative model is by asking if a given model has some close connection with the solution of important particle physics problems beyond just the origin of neutrino mass. One that has been explored at length in the literature is a possible connection to dark matter. Examples of such models will be given in more detail in later sections. Here, we simply mention some systematic analyses of what new symmetries can be imposed in radiative models to stabilize dark matter [110, 111]. Farzan et al. [110] classified the symmetries Gν that can be imposed in order to ensure that the first non-zero contribution to O1 occurs at a given loop order, by forbidding all potential lower-order contributions. All standard model particles are singlets under Gν, implying that the lightest of the exotics that do transform under this symmetry must be stable if the symmetry remains exact, establishing a connection with dark matter. Restrepo et al. [111] performed a systematic analysis of radiative models in a certain class in order to find those that have viable dark matter candidates. The considered models are those that generate mass at 1-loop level using exotics that are at most triplets under weak isospin, and where the stabilizing symmetry is Z2. They found 35 viable models. A similar analysis, but requiring 2-loop neutrino mass generation, can be found in Simoes and Wegman [112].

Besides dark matter, radiative neutrino mass models may also be connected to other physics beyond the SM such as the anomalous magnetic moment of the muon, the strong CP problem, the baryon asymmetry of the Universe or B-physics anomalies, among others. Phenomenology related to radiative neutrino mass models is briefly discussed in section 4 in general and an example of a possible connection to the recent B-physics anomalies is presented in section 5.4.

3. Radiative generation of neutrino masses

We adopt the classification of radiative neutrino mass models according to their Feynman diagram topology13, but refer to the other classification schemes where appropriate. In particular, we indicate the lowest-dimensional non-trivial ΔL = 2 operator which is generated beyond the Weinberg operator LLHH. These ΔL = 2 operators capture light particles which are in the loop to generate neutrino mass and are very useful to identify relevant low-energy phenomenology.

In the subsections 3.1–3.3 we classify Majorana neutrino mass models proposed in the literature according to their topology and specifically discuss models with SM gauge bosons in the loop in section 3.4. In section 3.5 we review Dirac neutrino mass models and briefly comment on models based on the gauge group SU(3)c × SU(3)L × U(1)X in section 3.6.

3.1. 1-loop Majorana neutrino mass models

This section is divided into several parts: (i) 1-loop UV completions of the Weinberg operator, (ii) 1-loop seesaws, (iii) UV completions with additional VEV insertions, (iv) 1-loop UV completions of the higher dimensional operators and (v) other 1-loop models. Notice that the first part includes models with multi-Higgs doublets, while the second part discusses external fields which transform under an extended symmetry. Besides the genuine topologies discussed in section 2, there are models based on the non-genuine 1-loop topologies in Figure 6.

Figure 6

Figure 6

Non-genuine topologies of the Weinberg operator.

3.1.1. Weinberg operator LLHH

We follow the general classification of UV completions of the Weinberg operator at 1-loop [100] discussed in section 2.2.2. The six genuine topologies are shown in Figure 3. Analytic expressions for all 1-loop topologies are listed in the appendix of Bonnet et al. [100].

Here we list the theories falling into respective categories. As the topologies stay the same while incorporating multiple Higgs doublets, theories with more than one Higgs doublet will also be listed here. Models in which the generation of neutrino mass relies on additional VEVs connected to the neutrino mass loop diagram are discussed in section 3.1.3. We first discuss the models based on topology T3, the only one with a quartic scalar interaction, before moving on to the other topologies.

T3

Topology T3 is one of the most well-studied. It was first proposed in Ma [99] and its first realization, the scotogenic model with a second electroweak scalar doublet and sterile fermion singlets (at least two) both odd under a Z2 symmetry, was later proposed in Ma [113]. See section 5.3 for a detailed discussion of the model. Its appeal lies in the simultaneous explanation of dark matter, which is stabilized by a Z2 symmetry. A crucial ingredient is the quartic scalar interaction (Hη)2 (see Equation 95) of the SM Higgs boson H with the electroweak scalar doublet η in the loop. This scalar interaction splits the masses of the neutral scalar and pseudoscalar components of η. Neutrino masses vanish in the limit of degenerate neutral η scalar masses. Several variants of the scotogenic model have been proposed in the literature: with triplet instead of singlet fermions [114116], an extension with an additional singlet scalar [117], one fermionic singlet and two additional electroweak scalar doublets [118], scalar triplets [119], colored scalars and fermions [120, 121], a vector-like fermionic lepton doublet, a triplet scalar, and a neutral [122, 123] or charged [124] singlet scalar, vector-like doublet and singlet fermions and doublet scalar, which contains a doubly charged scalar [125], higher SU(2) representations [126129], an extended discrete symmetry with Z2 × Z2 [130, 131] or Z2 × CP [132], a discrete flavor symmetry based on S3 [133], A4 [134137], Δ(27) [138, 139], which is either softly-broken or via electroweak doublets, and its embedding in (grand) unified theories [137, 140143]. Finally, the authors of Megrelidze and Tavartkiladze [144] proposed the generation of neutrino mass via lepton-number-violating soft supersymmetry-breaking terms. In particular the generation of the dimension-4 term with left-handed sleptons leads to models based on the topology T-3 with supersymmetric particles in the loop. Another variant involves a global continuous dark symmetry [145], Hagedorn, (in prep), termed the generalized scotogenic model.

T1-i

Ma [146] discusses a supersymmetrized version of the scotogenic model, which is based on topology T3 and we discuss in detail in section 5.3. The topology necessarily differs from T3 because the term (Hη)2 is not introduced by D-terms. An embedding of this model in SU(5) is given in Ma [147]. In a non-supersymmetric context, the same topology is discussed in Farzan [117], which introduces one real singlet scalar, in the context of a (dark) left–right symmetric model [148, 149], and in Budhi et al. [150], Kashiwase and Suematsu [151], and Budhi et al. [152], which introduce multiple singlet scalars to connect the two external Higgs fields. The term (Hη)2, which is essential to generate topology T3, is neglected in Budhi et al. [150], Kashiwase and Suematsu [151], and Budhi et al. [152] and thus neutrino mass is generated via topology T1-i. One of the singlet scalars in the neutrino mass model can be the inflaton via a non-minimal coupling with the Ricci-scalar. The term (Hη)2 can be explicitly forbidden by imposing a U(1) symmetry, which is softly broken by the CP-violating mass term χ2 of a complex scalar field χ [153]. Finally the authors of Lu and Gu [154] proposed a model with electroweak singlet and triplet scalars as well as fermions and study the dark matter phenomenology and leptogenesis.

T1-ii

Among the models based on the topology T1-ii, there are four possible operators which models are based on. Besides models with only heavy new particles, there are models with SM charged leptons, down-type quarks, or up-type quarks in the loop, which are based on the operators O2 and O3, respectively. We first discuss the models based on operator O2. The first radiative Majorana neutrino mass model, the Zee model [104], is based on this operator. See section 5.1.1 for a detailed discussion of its phenomenology. Several variants of the Zee model exist in the literature. The minimal Zee-Wolfenstein model [105] with a Z2 symmetry to forbid tree-level FCNCs has been excluded by neutrino oscillation data [155, 156], while the general version with both Higgs doublets coupling to the leptons is allowed [91, 157]. Imposing a Z4 symmetry [158] allows to explain neutrino data and forbid tree-level FCNCs in the quark sector. Previously in Aranda et al. [159] a flavor-dependent Z4 symmetry was used to obtain specific flavor structures in the quark and lepton sector. A supersymmetric version of the Zee model has been proposed in Leontaris and Tamvakis [160], Haba et al. [161], Cheung and Kong [162], and Kanemura et al. [163]. Its embedding into a grand unified theory has been discussed in Zee [104], Tamvakis and Vergados [164] and Fileviez Perez and Murgui [165], and in models with extra dimensions in Chang and Ng [166] and Chang et al. [167].

Other flavor symmetries beyond Z4 have been studied in Babu and Mohapatra [168, 169], Koide and Ghosal [170], Kitabayashi and Yasue, [171], Adhikary et al. [172], Fukuyama et al. [173], Aranda et al. [174, 175], and Okamoto and Yasue [176] studied the Zee model when the third generation transforms under a separate SU(2) × U(1) group. Babu and Mohapatra [168, 169] studied large transition magnetic moments of the electron neutrino, which was an early, now excluded, explanation for the solar neutrino anomaly. General group theoretic considerations about the possible particle content in the loop are discussed in Ma [99].

Models with multiple leptoquarks, which mix among each other, also generate neutrino mass via topology T1-ii. We discuss this possibility in more detail in section 5.4.1. They induce the operator O3 if the leptoquark couples down-type quarks to neutrinos. Well-studied examples of leptoquarks are down-type squarks in R-parity violating SUSY models, which generate neutrino masses, as was first demonstrated in Hall and Suzuki [177]. Specific examples with multiple leptoquarks which mix with each other were discussed in Nieves [178], Chua et al. [179], Mahanta [180], Aristizabal Sierra et al. [181], Helo et al. [182], Päs and Schumacher [183], Cheung et al. [184], Doršner et al. [185]. There are several supersymmetric models [179, 186190] which generate neutrino mass via different down-type quarks or charged leptons in the loop and consequently induce the operators O3 and O2, respectively. Finally, there are models with only heavy particles in the loop such as the inert Zee model [191] or supersymmetric models with R-parity conservation [192, 193].

T1-iii

This topology was first proposed in Ma [99] and it naturally appears in the supersymmetrized version of the scotogenic model [146, 147, 194203] together with topology T1-i. The topology can be used to implement the radiative inverse seesaw [204206], which resembles the structure of the inverse seesaw [207, 208]. This model has been extended by a softly-broken non-Abelian flavor symmetry group [209211] in order to explain the flavor structure in the lepton sector. The SUSY model in Ma and Sarkar [212] generates neutrino mass via sneutrinos and neutralinos in the loop. This mechanism was first pointed out in Hirsch et al. [213]. In the realization of Ma and Sarkar [212], the masses of the real and imaginary parts of the sneutrinos are split by the VEV of a scalar triplet, which only couples to the sneutrinos via a soft-breaking term and thus does not induce the ordinary type-II seesaw. Similarly it has been used in a model with vector-like down-type quarks [214, 215], which requires mixing of the SM quarks with the new vector-like quarks. This model leads to the operator O3.

3.1.2. 1-loop seesaws and soft-breaking terms

For completeness we also include the two possible 1-loop seesaw topologies T4-2-i and T4-3-i which have been identified in Bonnet et al. [100]. Topology T4-2-i always involves a electroweak scalar triplet like in the type-II seesaw mechanism and topology T4-3-i contains an electroweak singlet or triplet fermion like in the type-I or type-III seesaw mechanism, respectively. Based on our knowledge, there are currently no models based on topologies T4-2-i and T4-3-i in the literature.

Finally, although the topology T4-2-ii shown in Figure 6C has been discarded in Bonnet et al. [100], because it is generally accompanied by the tree-level type-II seesaw mechanism, there are three models based on this topology [216218]. They break lepton number softly by a dimension-2 term and thus there is no tree-level contribution by forbidding the “hard-breaking” dimension-4 terms which are required for the type-II seesaw mechanism. Similar constructions may be possible for other topologies and lead to new interesting models.

3.1.3. Additional VEV insertions

The above discussed classification technically does not cover models with additional scalar fields, which contribute to neutrino mass via their vacuum expectation value in contrast to being a propagating degree of freedom in the loop. Inspired by the above classification, we similarly classify these new models according to the topologies in Figure 3 by disregarding the additional VEV insertions.

T1-i

There are several radiative neutrino mass models which are based on a U(1) symmetry, which is commonly broken to a remnant Z2 symmetry: there are models based on a global Peccei-Quinn U(1)PQ symmetry [219, 220], which connects neutrino mass to the strong CP problem, a local U(1)B−L symmetry [221223] and local dark U(1) symmetry [224226]. The authors of Ho et al. [221] systematically study radiative neutrino mass generation at 1-loop (but also 2-, and 3-loop) level based on a gauged U(1)B−L symmetry, which is broken to a ZN symmetry. The models in Chang and Wong [224], Dasgupta et al. [219], Lindner et al. [225], Adhikari et al. [227], Kownacki and Ma [226] also have a contribution to neutrino mass at 2-loop order based on a Cheng-Li-Babu-Zee (CLBZ) topology.

T1-ii

All of the models with additional VEV insertions rely on the breaking of a symmetry: left–right symmetry [228230], a more general SU(2)1×SU(2)2 symmetry [231], a flavor symmetry [232234], U(1)B−L [235], and dilation symmetry [236]. All these models lead to the operator O2. Foot et al. [236] discusses in particular the following two 1-loop models: the scale-invariant Zee model and a scale-invariant model with leptoquarks which induces O3. Finally, there is the inert Zee model with a flavor symmetry [237, 238].

T1-iii

The model in Nomura et al. [239] relies on the VEVs of a scalar triplet and a septuplet which are subject to strong constraints from electroweak precision tests in particular from the T (or ρ) parameter. The minimization of the potential is not discussed, but the VEVs can in principle be introduced via the linear term in the scalar potential, which leads to the operator at 2-loop level, because the linear term for the septuplet is only induced at the 1-loop level. The topology can also be generated by new heavy lepton-like doublets and sterile fermions, which are charged under a new gauged dark U(1) in addition to a Z2 symmetry [240].

T3

There are several variants of the scotogenic model with additional VEV insertions. Most of them are based on an extended symmetry sector, such as a discrete Z3 instead of a Z2 symmetry [241, 242], dilation symmetry [236, 243245], a gauged U(1)B−L [246250], global U(1)B−L [251], a general gauged U(1) [252254], continuous U(1) flavor symmetry [255, 256], a discrete flavor symmetry based on D6 [257], A4 [258262] or S4 [263], and different LR symmetric models without a bidoublet [264]. Apart from additional symmetries, the mixing of the fermionic singlet with a fermionic triplet in the loop requires the VEV of an electroweak triplet with vanishing hypercharge [265267]. Finally, the two models discussed in Okada and Yagyu [268, 269] rely on a similar topology as the scotogenic model, but with triplet VEVs instead of electroweak doublet VEVs.

T4-2-i

Based on our knowledge, there are currently no models based on topology T4-2-i in the literature.

T4-3-i

Wang and Han [270] proposed a model which reduces to topology T4-3-i after breaking of the U(1)B−L symmetry. As the Majorana mass term for the fermionic pure singlet is not introduced, there is no inverse seesaw contribution to neutrino mass after the breaking of the U(1)B−L symmetry and neutrino masses are generated at 1-loop level.

T4-1-i/ii

These types of models contain a triplet scalar which couples to the lepton doublet as per the tree-level type-II seesaw. However, the neutral component of the triplet scalar gets an induced VEV at 1-loop and thus generates neutrino masses effectively at 1-loop. The model in Nomura and Okada [271] is based on topology T4-1-i shown in Figure 6A, which is finite due to additional VEV insertions on the fermion line. The model in Kanemura and Sugiyama [272] is based on topology T4-1-ii shown in Figure 6B. The tree-level contribution is forbidden by a discrete symmetry and renormalizability of the theory. However, at loop-level neutrino mass is generated by a dimension-7 operator with two additional SM singlet fields s1. Note in both cases an extra symmetry such as U(1)B−L or a discrete symmetry and lepton number is needed to forbid the contribution from the tree-level type-II seesaw. Topology T4-1-ii is also induced in the SUSY model in Figueiredo [273] and Franceschini and Mohapatra [201] after the breaking of SUSY and the discrete Z4 symmetry.

3.1.4. Higher-dimensional weinberg-like operators

Apart from UV completions of the Weinberg operator, there are a few models which induce one of the higher dimensional operators with additional Higgs doublets at 1-loop level.

Dimension-7 ()

The first model which induced the dimension-7 operator at 1-loop level in a two Higgs doublet model was proposed in Kanemura and Ota [274]. It was realized using at most adjoint representations and an additional softly-broken Z5 symmetry and an exact Z2 symmetry and thus allows to use the topologies T12 (Figure 4E) and T31 (Figure 7), which would otherwise be accompanied by the dimension-5 operator O1. If the Zee model is extended by a triplet Majoron [275, 276] the operator is induced at tree-level. After closing the loop of charged leptons via topology T3 (Figure 4B), the dimension-7 operator is obtained. Cepedello et al. [101] systematically studies the possible 1-loop topologies of and explicitly shows several models: the only genuine model without representations beyond the adjoint of SU(2) is based on topology T11, while the other models use quadruplets or even larger representations to realize the other genuine topologies.

Figure 7

Figure 7

Non-genuine 1-loop topology T31 for operator .

Dimension-9 ()

In Law and McDonald [277] and Baldes et al. [278] neutrino masses are generated via a radiative inverse seesaw. The mass of the additional SM singlets is induced at tree-level and then first transmitted to the neutral components of new electroweak doublets via a 1-loop diagram, before it induces neutrino mass via the seesaw. It leads to the dimension-9 operator via the four VEV insertions on the scalar line of the 1-loop diagram. There is also a 2-loop contribution, which may dominate neutrino mass depending on the masses of the new particles.

Dimension-11 ()

The model proposed in Aranda and Peinado [279] relies on the VEV of a 7-plet χ, which is induced via a non-renormalizable coupling, linear in χ, to six electroweak Higgs doublets.

As can be seen from the discussion above, in order to generate Weinberg-like effective operators at dimension larger than five, typically extra symmetries (in some cases large discrete symmetries), new large representations, a large number of fields or a combination of all the previous need to be invoked. This makes the model-building of such scenarios much more involved than for the case of the Weinberg operator.

3.1.5. Other 1-loop models

Apart from the models in the general classification [100], it is possible to generate neutrino mass via a radiative inverse seesaw mechanism shown in Figure 8 at 1-loop order, which has been proposed in Ahriche [280]. Tree-level contributions are forbidden by a softly-broken Z4 symmetry. The soft-breaking is indicated by the cross on the scalar line. Note the cross on the fermion line in the loop denotes a Majorana mass term, while the other two denote Dirac mass terms.

Figure 8

Figure 8

Radiative inverse seesaw.

Finally we would like to comment on one further possibility to generate neutrino mass at 1-loop order. If the neutrino masses vanish at tree-level in type-I seesaw model, then 1-loop electroweak corrections give the leading contribution [281]:14 non-zero neutrino masses are induced by finite 1-loop diagrams with either a Z-boson or a Higgs boson. The UV divergent part of the 1-loop corrections to the Weinberg operator cancel due to the absence of a tree-level contribution. This has been explicitly shown in Pilaftsis [281] with a calculation in the mass basis. In terms of the classification of 1-loop topologies, these diagrams correspond to the topologies T3 and T1-iii for the Higgs and Z-boson in the loop, respectively. The vanishing of the tree-level contribution can be achieved using a specific texture in the seesaw model with SM singlet fermions S [285] in addition to the right-handed neutrinos N

in the basis (ν, N, S). In the limit μS → 0 the tree-level contribution to the active neutrinos exactly vanishes and neutrino masses are generated at 1-loop order. This construction has been denoted minimal radiative inverse seesaw [285].

This texture can be obtained by imposing a U(1) symmetry under which S is charged. After it is spontaneously broken by the VEV of a SM singlet scalar η, the Yukawa interaction SNη generates the term MN without generating a Majorana mass term μS for the fermionic singlets S or a coupling of S to the SM lepton doublets L at the renormalizable level.

3.2. 2-loop Majorana neutrino mass models

The possible 2-loop topologies of the Weinberg operator have been discussed in Aristizabal Sierra [102]. We will closely follow this classification. All possible genuine 2-loop topologies are shown in Figure 5. Analytic expressions for the 2-loop diagrams are summarized in the appendix of Aristizabal Sierra [102] and are based on the results in McDonald and McKellar [286] and Angel et al. [287]. Most topologies can be considered as variations of a few 2-loop models discussed in the literature: (i) variations of the Cheng-Li-Babu-Zee (CLBZ) topology [76, 288, 289], (ii) the Petcov-Toshev-Babu-Ma (PTBM) topology [290292], and the so-called rainbow (RB) topology [102]. In the following we will further distinguish between fermion and scalar lines and show in Figures 9, 10B the relevant diagrams of genuine topologies and the internal-scalar-correction (ISC)-type topology which are used in the following discussion. The first two subsections discuss models based on genuine topologies, the third one models based on non-genuine topologies, and the last one models based on multiple topologies.

Figure 9

Figure 9

Relevant genuine 2-loop topologies.

Figure 10

Figure 10

Non-genuine 2-loop topologies.

3.2.1. Genuine 2-loop topologies

The relevant diagrams for the genuine topologies are shown in Figure 9.

CLBZ-1

The topology CLBZ-1 is displayed in Figure 9A. The first model was independently proposed and studied by Zee [288] and Babu [289], and is commonly called Zee-Babu model (See a more detailed discussion in section 5.1.2). It also leads to the operator O9. A scale-invariant version of the model has been proposed in Foot et al. [236]. It has been extended to include a softly-broken continuous LeLμLτ flavor symmetry [293, 294] or discrete flavor symmetry [295], and has been embedded in a SUSY model [296, 297]. The same topology has also been used for models with quarks instead of charged leptons inside the loop. They rely on the introduction of a leptoquark and a diquark [298300] and lead to operator O11. Similarly, there is a version without light fields in the loop [221, 301303]. The models in Ho et al. [221] are part of a systematic study of models based on a gauged U(1)B−L which is broken to a ZN symmetry.

CLBZ-3

Topology CLBZ-3 is depicted in Figure 9B and only differs from topology CLBZ-1 in the way how the Higgs VEVs are attached to the loop diagram: Topology CLBZ-3 has the Higgs VEVs attached to two of the scalar lines, while they are attached to the internal fermion lines for CLBZ-1. Cheng and Li [76] listed several possible neutrino mass models, including the first 2-loop model which was based on topology CLBZ-3 with an effective scalar coupling. A possible UV completion was presented with an electroweak quintuplet scalar. This UV completion leads to the operator (with an additional VEV insertion from an electroweak quintuplet scalar). All models [221, 304307] based on topology CLBZ-3 only contain heavy fields.

CLBZ-8

The topology is shown in Figure 9C. Variants of the Zee-Babu model have also been embedded in grand unified theories [308]. In case of SU(5), there is a 5-plet of matter particles in the loop which leads to the effective operators O9 and O11.

CLBZ-9

Topology CLBZ-9 which is displayed in Figure 9D has been utilized in a model with two diquarks [215].

CLBZ-10

The same paper also introduces another model with two diquarks which is based on topology CLBZ-10, shown in Figure 9E.

PTBM-1

The first model to utilize the topology Figure 9F, although in presence of a tree-level contribution, was presented in Petcov and Toshev [290], Babu and Ma [291], Branco et al. [292] and Babu and Ma [309]. Neutrino mass receives a 2-loop correction via the exchange of two W-bosons as shown in Figure 9F. This idea has been recently revived and experimentally excluded in the context of extra chiral generations [310], but the mechanism can still work in the case of vector-like leptons. Lepton number is violated by the SM singlet Majorana fermion N in the center of the diagram and thus there is a tree-level contribution in addition to the 2-loop contribution to neutrino mass. Lepton number can equally well be broken by the type-III seesaw, when the fermionic singlet is replaced by a fermionic triplet [311]. The model in Babu and Julio [312] has one of the W-bosons replaced by scalar leptoquarks and it is consequently not accompanied by a tree-level contribution. The 1-loop contribution induced by the mixing of the leptoquarks vanishes, because the left-chiral coupling of one of the leptoquarks is switched off [313]. All models with W bosons will lead to operators with derivatives in the classification according to ΔL = 2 operators. Finally, Angel et al. [287] proposed a model with a scalar leptoquark and colored octet fermion.

3.2.2. Genuine topologies with additional VEV insertions

Similar to the 1-loop models, we also categorize the models with additional VEV insertions following the classification of Aristizabal Sierra et al. [102].

CLBZ-1

There are several models based on topology CLBZ-1 (shown in Figure 9A), which all induce the operator O9. Bamba et al. [314] discusses a possible connection of neutrino mass with dark energy. Porto and Zee [315] proposed a model with one electroweak Higgs doublet field per lepton generation, an extension of the so-called private Higgs scenario. Finally, Lindner et al. [316] discusses an extension of the Zee-Babu model by a global U(1)B−L symmetry, which is spontaneously broken to a Z2 subgroup. This implies the existence of a Majoron and a DM candidate.

CLBZ-3

Chang and Keung [317] proposed a variant of the Zee-Babu model with an additional triplet Majoron, which is based on topology CLBZ-3 which is displayed in Figure 9B.

CLBZ-9

The topology CLBZ-9 is depicted in Figure 9D. The model in Guo [318] is based on a dark U(1) symmetry with only heavy fields in the loop.

RB-2

The model proposed in Kajiyama et al. [319] is based on U(1)B−L, which is broken to Z2. Apart from the VEV breaking U(1)B−L, neutrino mass is generated by a diagram with topology RB-2 which is shown in Figure 9G.

3.2.3. Non-genuine topologies

The relevant non-genuine 2-loop topologies are shown in Figure 10.

NG-RB-1

The non-genuine topology NG-RB-1 (Figure 10A) is generated in Nomura and Okada [320]. There are no lower-order contributions due to the U(1) symmetry, which is broken to Z2 as in the above-mentioned models.

Other non-genuine topologies

There are several models which generate vertices or masses of particles at loop level. The models in Aoki et al. [321, 322] realize an ISC-type topology which is shown in Figure 10B by softly breaking lepton number with a dimension-2 scalar mass insertion in the internal scalar loop. Similarly, Ma and Sarkar [323] discusses a supersymmetric model where the scalar-quartic coupling is induced after supersymmetry is softly broken and thus an ISC-type topology is induced for neutrino mass. The models in Kajiyama et al. [324] and Baek et al. [325] have only heavy particles in the loop and can be considered as a 1-loop scotogenic model, where the Majorana mass term for the SM singlet fermions is generated at 1-loop order. Thus, neutrino mass is effectively generated at 2-loop order. It can be considered as an RB-type topology. In contrast to the topology RB-2, the SM Higgs fields are attached to the outer scalar line (the one on the left in Figure 9G). Both models break U(1)B−L to a discrete ZN subgroup. Ghosh et al. [326] proposes another model based on an RB-type topology, where the Higgs fields couple to the fermions in the outer loop. The model features a stable dark matter candidate due to an imposed Z2 symmetry. Moreover, neutrino mass relies on the spontaneous breaking of an extended lepton number symmetry to a discrete Z2 subgroup. The models in Ma [327], Nasri and Moussa [328], Chao [329], Ma and Wudka [330] and Nomura and Okada [331] realize the type-I seesaw by generating the Dirac mass terms at 1-loop order, and the model in Okada and Orikasa [332] generates a radiative type-II seesaw contribution by generating the triplet VEV at 1-loop level, and thus the Weinberg operator at 2-loop level. Finally, Witten [333] and Arbelàez Rodríguez et al. [334] firstly generate the right-handed neutrino mass at 2-loop level in the context of a GUT, which induces the active neutrino mass via the usual seesaw mechanism. Similarly Law and McDonald [335] and Baldes et al. [278] realize a radiative inverse seesaw. The mass of additional singlets is generated at 2-loop order. The model is based on an additional gauged U(1) symmetry (which is spontaneously broken to its Z2 subgroup) to forbid the generation of neutrino mass at tree-level via the seesaw mechanism. The model can explain the matter-antimatter asymmetry of the universe, but not account for the dark matter abundance [278].

3.2.4. Models based on several topologies

Several models in the literature [144, 221, 336351] are based on multiple 2-loop topologies. We highlight three examples. Megrelidze and Tavartkiladze [144] proposed to generate neutrino mass via lepton-number-violating soft supersymmetry-breaking terms using the so-called type-II-B soft seesaw with electroweak triplet superfields. Integrating out the scalar components of the electroweak triplets leads to the dimension-5 lepton-number-violating term . Neutrino mass is generated at 2-loop via a diagram based on topology CLBZ-1 and diagrams which generate the couplings of the scalar component of the electroweak triplet superfield to two lepton doublets L on the one hand and the two electroweak Higgs doublets Hu on the other hand at the 1-loop level. Another interesting class of models are based on internal electroweak gauge bosons, which are based on CLBZ-type topologies and discussed in Chen et al. [345, 346], del Aguila et al. [347], Chen et al. [348], del Aguila et al. [349], King et al. [350] and Geng and Tsai [351]. All of them introduce a doubly-charged scalar and a coupling of the doubly-charged scalar to two W-bosons, which can be achieved via a mixing of the doubly-charged scalar with the doubly-charged scalar in an electroweak triplet scalar. Neutrino mass is typically generated via topologies CLBZ-1 and CLBZ-3 and induces the operator

This possibility is further discussed in section 3.4. Gauge bosons similarly can play an important role in the generation of neutrino mass in extended technicolor (ETC) models as discussed in Appelquist and Shrock [352, 353] and Appelquist et al. [354]. These models contain many SM singlet fermions and only a few elements of the neutral fermion mass matrix are directly generated by condensates, while many elements are generated at 1-loop (or higher loop) level via loop diagrams with ETC gauge bosons. In particular the relevant Dirac mass terms relevant for the active neutrino masses are generated at 1-loop level and thus neutrino mass is effectively generated at 2-loop (or even higher loop) level.

3.3. 3-loop Majorana neutrino mass models

Unlike 1-loop and 2-loop topologies, there is no systematic classification of all 3-loop topologies. Thus, we restrict ourselves to the existing 3-loop models in the literature and do not consider other topologies or different fermion flow for the given topologies. Most of the existing 3-loop models can be categorized in four basic types of diagrams shown in Figure 11 where we do not specify the Higgs insertions. The remaining models are either based on a combination of the listed topologies or the combination of a loop-induced vertex at 1- or 2-loop inside a loop diagram.

Figure 11

Figure 11

3-loop model topologies. Note that we do not specify the Higgs insertions.

3.3.1. The KNT models

The first 3-loop radiative neutrino mass model was proposed in Krauss et al. [355] with the topology shown in Figure 11A by Krauss, Nasri and Trodden (KNT) and it leads to the operator O9. We refer to radiative neutrino mass models sharing the same topology as KNT models and discuss them in more detail in section 5.2. A systematic study with several different variants can be found in Chen et al. [356]. The models of Chen et al. [357], Ahriche and Nasri [358], Ahriche et al. [359], Chen et al. [356], Ahriche et al. [360364] also generate the operator O9, the models of Chen et al. [356], Nomura et al. [365] and Cheung et al. [366] the operator O11 with down-type quarks, while the models in Okada and Okada [367], Chen et al. [356], Okada and Yagyu [368] and Cheung et al. [369] only have new heavy states in the loop.

3.3.2. AKS-type models

Neutrino mass can also arise at 3-loop order from the diagram shown in Figure 11C. The first model of such topology was proposed by Aoki, Kanemura, and Seto (AKS) in Aoki et al. [370] and is based on the operator with two Higgs doublets Hi. We will refer to it as the AKS model and more generally to models based on this topology as AKS-type models. It contains a second Higgs doublet and several SU(2)L singlets. The exotic particles can also be all electroweak singlets [221, 371]. The model in Gu [371] leads to the operator O9. Other variants include colored exotic particles such as leptoquarks [356, 372, 373], which generate the operators O11, 12, or electroweak multiplets [356, 374, 375] generating the operators O1,9. Note cross diagrams may be allowed in specific models.

3.3.3. Cocktail models

The third class of models are based on the two cocktail diagrams shown in Figures 11B,D. The name for the diagram has been coined by Gustafsson et al. [376], which proposed a 3-loop model with two W-bosons based on topology Figure 11B and consequently generated the operator , which are discussed in more detail in section 3.4. The same model has also been studied in Geng et al. [377]. The models in Hatanaka et al. [378] and Alcaide et al.[379] are based on the same topology, but with W bosons replaced by scalars. While Alcaide et al. [379] induces operator , the model of Hatanaka et al. [378] leads to operator O9. Finally, the fermionic cocktail topology Figure 11D is used in the models of Nishiwaki et al. [380], and Kanemura et al. [381], both of which generate operator O9.

Apart from the three classes of models, there are a few models which do not uniquely fit in any of the three classes. The model in Jin et al. [382] is based on topologies Figures 11A,C with two W-bosons and thus generates the operator . Nomura et al. [383] generates the mass of new exotic fermions at 2-loop level via a CLBZ-type diagram, which in turn generate neutrino mass at 1-loop. Geng and Huang [384] studies a 2-loop model based on the operator O8, which itself is generated at 1-loop order.

Most of the 3-loop models need to impose extra discrete symmetries such as Z2 or a continuous U(1) symmetry to forbid lower-loop or tree-level contributions, unless accidental symmetries exist and thus partly require other VEV insertions. One example is to employ higher dimensional representation of SU(2)L [362], e.g., septuplet, in the spirit of minimal dark matter [385, 386] such that undesirable couplings are forbidden by the SM gauge group alone. Due to the existence of the extra imposed or accidental symmetries, 3-loop models serve as a natural playground for DM physics.

3.4. Models with gauge bosons

The first model [290292] with gauge bosons in the loop uses the topology PTBM-1 and leads to operators built from two lepton doublets including covariant derivatives. However, it also has a tree-level contribution, while models based on operators with right-handed charged leptons are genuine radiative neutrino mass models.

In del Aguila [107] two LNV effective operators with gauge bosons, i.e., present in covariant derivatives, were considered, which allowed to have neutrinoless double beta decay rates generated at tree level thanks to new couplings to the SM leptons15. Interestingly, depending on the chirality of the outgoing leptons in 0νββ, there are two new operators (beyond the standard contribution from the Weinberg operator which involves left-handed electrons). For left–right (LR) chiralities of the outgoing electrons, there is a dimension-7 operator:16

For right-right (RR) chiralities, there is a dimension-9 operator as define in Equation (34). After electroweak symmetry breaking, these operators generate the relevant vertices for 0νββ at tree level: and , respectively. The contributions of and to 0νββ are depicted in Figures 12B,C respectively, where the red point denotes the effective operator insertion.

Figure 12

Figure 12

Possible contributions to 0νββ. The red dot indicates the ΔL = 2 effective vertex. Figure reproduced from del Aguila [349].

The lowest order contributions from these operators to neutrino masses occur at 1- and 2-loop orders, respectively, via the diagrams of Figure 13. The dominant contributions come from matching (see also Babu and Leung [85], de Gouvêa et al. [87], Angel et al. [90] and de Gouvêa et al.[88] for estimates of the matching contributions to neutrino masses of LNV operators), which using dimensional analysis can be estimated to be given by del Aguila et al. [107]:

for and by

for . Notice that the appearance of the chirality-flipping charged lepton masses is expected in order to violate lepton number in the LH neutrinos, which naturally generates textures in the neutrino mass matrix.

Figure 13

Figure 13

Lowest order contributions of (left, at 1-loop order) and (right, at 2-loop order) to neutrino masses. The red dot indicates the ΔL = 2 effective vertex. Figure reproduced following del Aguila [107].

Possible tree-level UV completions which have new contributions to 0νββ at tree level were outlined in del Aguila et al. [107]. See also del Aguila [349], which provides a summary of two examples of models generating and , respectively. The UV model of [347] generates 0νββ at tree level, while neutrino masses are generated as expected at 2-loop order. It includes a doubly-charged singlet, a Y = 1 triplet scalar and a real singlet. In order to prevent tree-level neutrino masses as in type-II seesaw via the latter field, a discrete Z2 symmetry, which was spontaneously broken by the VEV of the singlet, was added. Recently a variation has been studied, in which the Z2 symmetry is exact, such that there is a good dark matter candidate, which is a mixture of singlet and triplet [379]. In this case, the contributions to 0νββ and to neutrino masses are further shifted by one extra loop, i.e., they are generated at 1- and 3-loop orders, respectively. Gustafsson et al. [376, 387] studied also a specific model with a dark matter candidate, named the cocktail model, which generated at 1-loop order, i.e., 0νββ at 1-loop order and therefore neutrino masses at 3-loop order. It includes a singly-charged singlet, a doubly-charged singlet and a Y = 1 scalar doublet, together with a discrete symmetry Z2 under which all the new fields except the doubly-charged are odd. Other models generating were presented in Chen et al. [345, 346, 348], King et al. [350], Geng and Tsai [351], and Liu and Gu [388].

3.5. Radiative Dirac neutrino mass models

Although Majorana neutrinos are the main focus of research, Dirac neutrinos are a viable possibility to explain neutrino mass. It is noteworthy that the first radiative neutrino mass model [389] was based on Dirac neutrinos. In recent years, there has been an increased interest in Dirac neutrinos and, in particular, there are a few systematic studies on the generation of Dirac neutrino mass beyond the simple Yukawa interaction, which include both tree-level and loop-level realizations, besides several newly-proposed radiative Dirac neutrino mass models, which we will outline below.

Ma and Popov [56] and Wang and Han [60] performed a study of Dirac neutrino mass according to topology at tree-level and 1-loop level. There are only two possible one-particle-irreducible topologies for the Dirac Yukawa coupling at 1-loop, which are shown in Figure 14. The simplest radiative Dirac neutrino mass models are based on a softly-broken Z2 symmetry, which is required to forbid the tree-level contribution, and generate the topologies in Figure 14. Wang et al. [390] studied scotogenic-type models with a U(1)B−L symmetry at 1- and 2-loop order. Finally, Kanemura et al. [391] takes a model-independent approach and discusses the possible flavor structures of the induced Dirac mass term under a number of constraints: The fermion line only contains leptons and each lepton type can appear at most once.

Figure 14

Figure 14

Generation of the Dirac neutrino Yukawa coupling at 1-loop order.

1-loop models

Many of the proposed 1-loop models are realized in a left–right symmetric context without SU(2) triplet scalars [56, 389, 392397]. Rajpoot [398] attempted the generation of Dirac neutrino masses in the context of a model where hypercharge emerges as diagonal subgroup of U(1)L×U(1)R. To our knowledge the generation of Dirac neutrino mass at 1-loop level with a softly broken Z2 was first suggested in Kanemura et al. [399] based on topology Figure 14B. Farzan and Ma [400] implements the first scotogenic model of Dirac neutrino mass by using a dark Z2 and softly-broken Z2. Both of these possibilities have been studied in more detail in the systematic studies outlined above. Another way to explain the smallness of Dirac neutrino mass is via a small loop-induced VEV [401]. Finally, Borah and Dasgupta [402] discusses a left–right symmetric model with pseudo-Dirac neutrinos. The tree-level Majorana mass terms are not allowed, because the bidoublet is absent and the coupling of the left-handed triplet to leptons is forbidden by a discrete symmetry.

2-loop models

Two explicit models of Dirac neutrino mass have been discussed in Bonilla et al. [403] and Kanemura et al. [404] apart from the general classification [390]. They are both based on a U(1) symmetry, a dark U(1) and lepton number, respectively. The U(1) symmetry is broken to a discrete subgroup and thus both models feature a stable dark matter candidate.

3-loop model

Finally, a Dirac neutrino mass term can also be induced via a global chiral anomaly term [405]. The five-dimensional anomaly term with the pseudo-scalar a and the (dual) field strength tensor Fμν () is induced at 1-loop level and leads to a Dirac mass term at 2-loop order, being effectively a 3-loop contribution.

3.6. 331 models

Another interesting class of models is based on the extended gauge group SU(3)c × SU(3)L × U(1)X. The SM gauge group can be embedded in several different ways and is determined by how the hypercharge generator is related to the generator T8 of SU(3)L and the generator X of U(1)X,

where β is a continuous parameter. In addition to one radiative Dirac neutrino mass model [406], several radiative Majorana neutrino mass models have been proposed at 1-loop level [407423], 2-loop order [424429], 3-loop order [430], and even at 4-loop order [98]. As lepton number violation (LNV) in 331 models and in particular neutrino mass generation has been discussed in a recent review [431], we refer the interested reader to it for a detailed discussion. However, we highlight one model based on gauged lepton number violation [419421], which generates neutrino mass via lepton number violation in the 1-loop diagram shown in Figure 15 with the SU(3)L×U(1)X gauge bosons, where Hi denotes the SM Higgs doublets, 〈1〉χ the VEV in the third component of SU(3)L and Nc the third partner of νL in the triplet of SU(3)L. Note that lepton number is broken by the mixing of the gauge bosons in the vertex at the top of the diagram.

Figure 15

Figure 15

Neutrino mass generation from gauged lepton number violation.

4. Phenomenology

In this section we revisit the most relevant phenomenological implications of radiative neutrino mass models. The possible signals are very model dependent, as each radiative model has its own particularities that should be studied on a case-by-case basis. However, in the following we will try to discuss generic predictions of these models, making use of simplified scenarios and/or of effective operators, and referring to particular examples when necessary.

4.1. Universality violations and non-standard interactions

In the SM, leptonic decays mediated by gauge interactions are universal. Several scenarios of physics beyond the SM have universality violations, that is, decays into different families (up to phase space-factors) are no longer identical17. These may or may not be related to neutrino masses, as lepton number is not violated in these interactions. Indeed, for instance a two Higgs doublet model with general Yukawa interactions breaks universality, irrespective of neutrino masses. In tree-level neutrino mass models, there are also violations of universality, mediated by the (singly) charged scalar boson in the type-II seesaw model, or due to the non-unitarity of the leptonic mixing matrix in type-I and type-III seesaw models when the extra neutral fermions are heavy [432, 433].

In some of the radiative models there can be violations of universality. One illustrative example of this case is due to the presence of a singly-charged singlet h with mass (as in the Zee and Zee-Babu models, see section 5.1). The relevant interaction is , where f is an antisymmetric matrix in flavor space and . Integrating out the singlet, one obtains the following dimension 6 effective operator [434]

One can see that this operator involves left-handed leptons, like charged currents in the SM18. This implies that it interferes constructively with the W boson, modifying among others the muon decay rate [435]. Therefore, the Fermi constant which is extracted from muon decay in the SM, , and that in a model with a singly-charged singlet, , are different, i.e., . Their ratio obeys to leading order in f:

The new Fermi constant is subject to different constraints. For example, from measurements of the unitarity of the CKM matrix, as the Fermi constant extracted from hadronic decays should be equivalent to that from leptonic decays, we can bound f:

Also leptonic decays which in the SM are mediated by charged-current interactions are not universal anymore. The ratio of leptonic decays among the different generations can be tested via the effective couplings given by

All these lead to strong limits on the f couplings depending on the mass on the singlet [95].

Furthermore, the new singly-charged scalar via the effective operator in Equation (39) induces neutrino interactions that cannot be described by W-boson exchange and are termed non-standard neutrino interactions (NSIs). Equation (39) is usually rewritten after a Fierz identity as

where are the NSI parameters given by

These could be in principle probed at neutrino oscillation experiments. However, typically whenever NSIs are induced, lepton flavor violating (LFV) processes are also generated, which are subject to stronger constraints. This is particularly the case for the four-lepton dimension 6 operators, due to gauge invariance. Models with large NSI are difficult to construct, and typically involve light mediators [436, 437]. We refer the reader to Davidson et al. [438], Ibarra et al. [439], Gavela et al. [440], Biggio et al. [441, 442], Antusch et al. [443] and Ohlsson [444] for studies of NSIs and their theoretical constraints.

4.2. Lepton flavor violation

One of the common predictions shared by most neutrino mass models (radiative or not) is the existence of LFV processes involving charged leptons with observable rates in some cases. Indeed, neutrino oscillations imply that lepton flavors are violated in neutrino interactions, and as in the SM neutrinos come in SU(2) doublets together with the charged leptons, also violations of lepton flavors involving the latter are expected. Which is the most constraining LFV observable is, however, a model-dependent question. It is thus convenient to use a parametrization that allows for a model-independent description of these processes. For each of the models one can then compute the relevant coefficients and apply the following formalism. We follow the notation and conventions of Porod et al. [445]19.

The general LFV Lagrangian can be written as

The first term contains the ℓ−ℓ−γ interaction Lagrangian, given by

where e is the electric charge, q is the photon momentum, are the standard chirality projectors and the indices {α, β} denote the lepton flavors. The first term in Equation (48) corresponds to the monopole interaction between a photon and a pair of leptons whereas the second is a dipole interaction term. In this parametrization the form factors vanish when the photon is on-shell, i.e., in the limit of q2 → 0. Similarly, the interaction Lagrangians with the Z and Higgs bosons are given by20

where p is the ℓβ 4-momentum, and

with the SM Higgs h. The general 4-lepton interaction Lagrangian can be written as

where in this case the indices {α, β, γ, δ} denote the lepton flavors and we have defined ΓS = 1, ΓV = γμ and ΓT = σμν. It is clear that the Lagrangian in Equation (51) contains all possible terms allowed by Lorentz invariance. Finally, the general 2ℓ2q 4-fermion interaction Lagrangian (at the quark level) can be split in two pieces

where

Here γ denotes the d-quark flavor and we are neglecting the possibility of quark flavor violation, which is beyond the scope of this review 21.

The parametrization used implies that the operators appearing in Equations (51), (53), and (54) have canonical dimension six. Therefore, the Wilson coefficients , and scale as 1/Λ2, where Λ is the new physics energy scale at which they are generated. Note this scale is unrelated to the scale at which lepton number is violated. The same comment applies to the dipole coefficients in Equation (48). In contrast, the rest of the coefficients discussed in this section, , and SL,R, are dimensionless (although their leading new physics contribution appears at order v22). If we restrict the discussion to flavor violating coefficients, they all vanish in the SM. Therefore, they encode the effects induced by the new degrees of freedom present in specific models.

It should be noted that all operators in the general LFV Lagrangian in Equations (48–54) break gauge invariance. For instance, they contain new charged lepton interactions, but not the analogous new interactions for the neutrinos, their SU(2)L doublet partners which are partly discussed in the previous subsection. This type of parametrization of LFV effects is correct at energies below the electroweak symmetry breaking scale, but it may miss relevant correlations between operators that are connected by gauge invariance in the underlying new physics theory. See for instance Pruna and Signer [450] for a discussion of LFV in terms of gauge-invariant operators.

We now proceed to discuss the LFV processes with the most promising experimental perspectives in the near future. We will provide simple analytical expressions in terms of the coefficients of the general LFV Lagrangian and highlight some radiative neutrino mass models with specific features leading to non-standard expectations for these processes. By no means this will cover all the models constrained by these processes, but will serve as a review of the novel LFV scenarios in radiative neutrino mass models.

Note, however, that there are other processes, which may yield stringent constraints in particular models: for instance in models with leptoquarks, the latter can mediate semi-leptonic τ-decays and leptonic meson decays at tree level. The LFV decays have also been investigated in several radiative models, although they typically have very low rates, see for instance Ghosal et al. [451] and Li et al. [452].

4.2.1. ℓα → ℓβγ

The most popular LFV process is ℓα → ℓβγ. There are basically two reasons for this: (1) for many years, the experiments looking for the radiative process μ → have been leading the experimental developments, with the publication of increasingly tighter bounds, and (2) in many models of interest these are the processes where one expects the highest rates. In fact, many phenomenological studies have completely focused on these decays, neglecting other LFV processes that may also be relevant.

The experimental situation in radiative LFV decays is summarized in Table 1. As one can easily see in this table, muon observables have the best experimental limits. This is due to the existing high-intensity muon beams. The current limit for the μ → branching ratio has been obtained by the MEG experiment, BR(μ → ) < 4.2·10−13 [453], slightly improving the previous bound also obtained by the same collaboration. This bound is expected to be improved by about one order of magnitude in the MEG-II upgrade [454]. The bounds in τ decays are weaker, with the branching ratios bounded to be below ~ 10−8, and some improvements are expected as well in future B-factories [456].

Table 1

LFV process BR Present bound Future sensitivity
μ → 4.2 × 10−13 [453] 6 × 10−14 [454]
τ → 3.3 × 10−8 [455] ~3 × 10−9 [456]
τ → μγ 4.4 × 10−8 [455] ~3 × 10−9 [456]

Current experimental bounds and future sensitivities for ℓα → ℓβγ branching ratios.

The decay width for ℓα → ℓβγ is given by Hisano et al. [457]

where α is the fine structure constant. Only the dipole coefficients , defined in Equation (48), contribute to this process. General expressions for these coefficients can be found in Lavoura [458].

The μ → limit is typically the most constraining one in most radiative neutrino mass models. One can usually evade it by adopting specific Yukawa textures that reduce the μ − e flavor-violating entries (see for example Schmidt et al. [459]) or simply by globally reducing the Yukawa couplings by increasing the new physics scale. However, in some cases this is not possible. A simple example of such situation is the scotogenic model [113] with a fermionic dark matter candidate. The singlet fermions in the scotogenic model only couple to the SM particles via the Yukawa couplings. Therefore, these Yukawa couplings must be sizable in order to thermally produce singlet fermions in the early universe in sufficient amounts so as to reproduce the observed DM relic density. This leads to some tension between the DM relic density requirement and the current bounds on LFV processes, although viable regions of the parameter space still exist [459, 460]. In contrast, in other radiative models the tight connection between neutrino masses and LFV implies suppressed ℓα → ℓβγ rates. This is the case of bilinear R-parity violating models [461463], see section 5.5 for a detailed discussion of this type of supersymmetric neutrino mass models.

4.2.2. ℓα → ℓβδδ

We now consider the ℓα → ℓβδδ 3-body decays. One can distinguish three categories: , (with β ≠ δ) and (also with β ≠ δ). These processes have received less attention even though the experimental limits on their branching ratios are of the same order as for the analogous ℓα → ℓβγ decays. We summarize the current experimental bounds and future sensitivities for the ℓα → ℓβδδ 3-body decays in Table 2. We note that an impressive improvement of four orders of magnitude is expected in the μ → eee branching ratio sensitivity thanks to the Mu3e experiment at PSI [464].

Table 2

LFV process BR Present bound Future sensitivity
μ → eee 1.0 × 10−12 [465] ~10−16 [464]
τ → eee 2.7 × 10−8 [466] ~10−9 [456]
τ → μμμ 2.1 × 10−8 [466] ~10−9 [456]
τeμ+μ 2.7 × 10−8 [466] ~10−9 [456]
τ → μe+e 1.8 × 10−8 [466] ~10−9 [456]
τe+μμ 1.7 × 10−8 [466] ~10−9 [456]
τ → μ+ee 1.5 × 10−8 [466] ~10−9 [456]

Current experimental bounds and future sensitivities for ℓα → ℓβδδ branching ratios.

The ℓα → ℓβδδ decay width receives contributions from several operators of the general LFV Lagrangian. In the case of the first category, , the decay width is given by Porod et al.[445]

in case of the second category, (with β ≠ δ), the expression is given by Abada et al. [467]

whereas for the third category, (with β ≠ δ), the decay width is given by Abada et al. [467]

Here we have defined

The masses of the leptons in the final state have been neglected in Equations (56–58), with the exception of the contributions given by the dipole coefficients , where infrared divergences would otherwise occur.

The dipole coefficients , which contribute to ℓα → ℓβγ, also contribute ℓα → ℓβδδ. It is easy to see how: the Feynman diagram contributing to ℓα → ℓβγ can always be supplemented with a flavor-conserving ℓδ−ℓδ−γ additional vertex resulting in a diagram contributing to ℓα → ℓβδδ22. In fact, such diagrams have been shown to be dominant in many models, the most popular example being the Minimal Supersymmetric Standard Model (MSSM). In this case, known as dipole dominance scenario, a simple proportionality between the decays widths of both LFV decays can be established. For example, in the β = δ case this proportionality leads to

which implies , making the radiative decay the most constraining process.

The dipole dominance assumption is present in many works discussing LFV phenomenology. However, it can be easily broken in many radiative neutrino mass models. This can happen in two ways:23

  • Due to tree-level LFV: In many radiative neutrino mass models the 4-lepton operators receive contributions at tree-level. The most prominent example of such models is the Zee-Babu model, in which the doubly-charged scalar k++ mediates unsuppressed decays. In such case one can easily find regions of parameter space where , see Herrero-Garcia et al. [95] for a recent study.

  • Due to loop-level LFV: Kubo et al. [468], Aristizabal Sierra et al. [469], Suematsu et al. [470], and Adulpravitchai et al. [471] explored the LFV phenomenology of the scotogenic model but only considered μ → . However, this assumption has been shown to be valid only in some regions of the parameter space. In fact, box diagrams contributing to 4-lepton coefficients can actually dominate, dramatically affecting the phenomenology of the scotogenic model [460, 472]. Qualitatively similar results have been found in other variants of the scotogenic model [129, 266]24. In fact, this feature is not specific of the scotogenic model and its variants: one can find other radiative neutrino mass models with loop contributions dominating over the dipole. For instance, Z-penguin contributions have been found to be dominant in the angelic model [90] and RνMDM models [473].

This clearly shows that radiative neutrino mass models typically have a very rich LFV phenomenology with new (sometimes unexpected) patterns and correlations.

4.2.3. μ−e conversion

The most spectacular improvements in the search for LFV are expected in μ − e conversion experiments. Several projects will begin their operation in the near future, with sensitivities that improve the current bounds by several orders of magnitude. The experimental situation is shown in Table 3.

Table 3

LFV process CR Present bound Future sensitivity
μ, Ti → e, Ti 4.3 × 10−12 [474] ~10−18 [475]
μ, Au → e, Au 7 × 10−13 [476]
μ, Al → e, Al 10−15−10−18 [477]
μ, SiC → e, SiC 10−14 [478]

Current experimental bounds and future sensitivities for μ−e conversion in nuclei.

The conversion rate, normalized to the the muon capture rate Γcapt, is given by Kuno and Okada [479] and Arganda et al. [480]

Z and N are the number of protons and neutrons in the nucleus and Zeff is the effective atomic charge [481]. GF is the Fermi constant, α is the electromagnetic fine structure constant, pe and Ee are the momentum and energy of the electron, mμ is the muon mass and Fp is the nuclear matrix element. and (with X = L, R and K = S, V) are effective couplings at the nucleon level. They can be written in terms of effective couplings at the quark level as

The numerical values of the relevant GK factors can be found in Kuno and Okada [479], Kosmas et al. [482] and Porod et al. [445]. For coherent μ − e conversion in nuclei, only scalar (S) and vector (V) couplings contribute and sizable contributions are expected only from the u, d, s quark flavors. The gXK(q) effective couplings can be written in terms of the Wilson coefficients in Equations (48), (53), and (54) as

where Qq is the quark electric charge (Qd = −1/3, Qu = 2/3) and for d-quarks (u-quarks), with X = L, R and K = S, V.

Radiative neutrino mass models can also be probed by looking for μ−e conversion in nuclei. As already pointed out, the search for this LFV process is going to be intensified in the next few years and, in case no observation is made, it will soon become one of the most constraining observables for this type of models. Similarly to the leptonic LFV 3-body decays discussed above, the dipole coefficients also enter the μ − e conversion rate, potentially dominating it. In this case, one can derive a simple relation [483]

where f(Z, N) is a function of the nucleus ranging from 1.1 to 2.2 for the nuclei of interest. The reader is referred to de Gouvea and Vogel [484] and Crivellin et al. [485] for a discussion on the complementarity of μ → and μ − e conversion in nuclei. One can easily depart from this dipole dominance scenario in radiative neutrino mass models due to the existence of sizable contributions to other LFV operators. For instance, non-dipole contributions have been shown to be potentially large in the scotogenic model in Toma and Vicente [472] and Vicente and Yaguna [460]. The dipole coefficients may also be reduced due to partial cancellations in non-minimal models, see for example Ahriche et al. [360, 361] and Rocha-Moran and Vicente [266]. Finally, as already pointed out in the case of ℓα → ℓβδδ decays, some radiative neutrino mass models contain new states that mediate LFV processes at tree level. For instance, in R-parity violating models with trilinear terms (discussed in section 5.5), the superpotential terms induce μ − e conversion at tree level [486]. This easily breaks the expectation in Equation (67).

Finally, we point out that the experiments looking for μ → eee and μ − e conversion in nuclei will soon take the lead in the search for LFV. Therefore, even if dipole contributions turn out to be dominant in a given model, μ → eee and μ − e conversion in nuclei might become the most constraining LFV processes in the near future. Prospects illustrating this point for specific radiative neutrino mass models have been presented in Angel et al. [90], Vicente and Yaguna [460], and Klasen et al. [487].

4.2.4.

In many radiative neutrino mass models, there can also be contributions to lepton-flavor violating Higgs (HLFV) decays, like h → τμ+, τe+ and their CP-conjugates. These same interactions, however, also generate LFV processes such as τ → μ(e)γ, as no symmetry can prevent the latter [488], which are subject to much stronger constraints. In the effective field theory with just the 125 GeV Higgs boson, HLFV decays involving the tau lepton can be sizable, and ATLAS and CMS constraints on its flavor violating couplings (shown in Table 4) are comparable or even stronger than those coming from low-energy observables [489491]. However, in UV models, specially in radiative neutrino mass models, the situation is generally the opposite.

Table 4

HLFV decay BR ATLAS CMS
h → τμ 0.0143 [492] 0.0025 [493]
h → τe 0.0104 [492] 0.0061 [493]

Experimental 95 % C.L. upper bounds on HLFV decays from ATLAS and CMS in the tau sector using the 13 TeV data sets.

The relevant gauge-invariant effective operators that generate HLFV are the Yukawa operator:

and derivative operators like

or

plus their Hermitian conjugates. In Herrero-García et al. [491] all the possible tree-level realizations of these operators were outlined, some of which include particles that are present in radiative neutrino mass models, as we will see below. In Figure 16, we show some possible UV completions of operators , and . The authors concluded that only can have sizable rates, and in particular only for UV completions that involve scalars, like in a type-III two-Higgs doublet model.

Figure 16

Figure 16

Different 1-loop UV completions of the Yukawa operator given in Equation (68), and the derivative operators given in Equation (69) and in Equation (70). F and F1,2 are fermion fields and S1,2 scalar fields. The Zee-Babu and the scotogenic models are examples of radiative models with HLFV generated at 1-loop order. Figure reproduced from Herrero-García et al. [491].

After electroweak symmetry breaking the Yukawa operator gives rise to the interaction Lagrangian in Equation (50). For instance, the SL,R couplings are given by

where Df is the SM flavor-diagonal contribution, not relevant for the present discussion, and CY is the Wilson coefficient of the operator defined in Equation (68). Focusing on the contributions from the Yukawa operator, the branching ratio of the Higgs into a tau and a muon reads:

Most radiative neutrino mass models generate HLFV at 1-loop order [491]25. For instance, the doubly-charged scalar singlet and the singly-charged scalar singlet of the Zee-Babu model (see section 5.1) generate respectively the derivative operators and at 1-loop order. The scotogenic model (see section 5.3) also generates HLFV at 1-loop order ().

We can estimate the loop-induced HLFV in radiative neutrino mass models. Denoting a generic Yukawa coupling of the fermions and scalars with the SM leptons as Y, and a scalar quartic coupling with the Higgs as λih, and taking into account that the amplitude of h → μτ involves a tau mass, one can estimate the dominant contribution to be [491]:

where M is the largest mass in the loop. In all these models, in addition to the loop factor, there are in general limits from charged LFV processes, as usually all radiative neutrino mass models have charged particles that can generate ℓα → ℓβγ. As τ → μγ typically gives the constraint , we get:

well below future experimental sensitivities. Thus, unless cancellations are invoked (which are difficult to achieve in all possible radiative decays), HLFV rates are very suppressed, well below future experimental sensitivities.

One class of models which can have large HLFV are those with another Higgs doublet such that both the SM and the new scalar doublet couple to the lepton doublets [488, 499502]. In such scenarios, both Yukawa couplings cannot be diagonalized simultaneously, which leads to LFV Higgs interactions. One example is the Zee model discussed in section 5.1, which can have BR(h → μτ) up to the percent level [91].

4.3. Anomalous magnetic moments and electric dipole moments

The anomalous magnetic moments (AMMs) and electric dipole moments (EDMs) of the SM leptons receive new contributions in radiative neutrino mass models (see Raidal et al. [503] for a review on the topic). These are contained in the dipole coefficients that also contribute to the radiative ℓα → ℓβγ decays, typically leading to tight correlations between these observables. Using the effective Lagrangian in Equation (48), the anomalous magnetic moment aα and the electric dipole moment dα of the charged lepton ℓα are given by Raidal et al. [503]

The experimental values for the AMMs and EDMs of charged leptons are collected in Table 5. In particular the muon AMM received a lot of intention in recent years due to the discrepancy between the experimentally measured value given in Table 5 and the SM prediction [504]

with the errors due to electroweak, lowest-order, and higher-order hadronic contributions.

Table 5

Lepton AMM a EDM d[ecm]
e (1159.65218091 ± 0.00000026) × 10−6 < 0.87 × 10−28
μ (11659208.9 ± 5.4 ± 3.3) × 10−10 (−1 ± 9) × 10−20
τ [−0.52, 0.013] [−2.20, 4.5] × 10−17+i[−2.50, 8.0] × 10−19

Experimental values for AMMs and EDMs [504].

Both statistical and systematic uncertainties are given for the muon AMM aμ.

There are many examples of radiative neutrino mass models leading to sizable effects in these two observables. For some examples in the case of AMMs see for instance Dicus et al. [234], Babu and Julio [312], Nomura et al. [239], Nomura and Okada [299], Chiang et al. [206], and Lee et al. [505]. In some cases, the new contributions effects can help close the gap between the theory prediction and the experimental measurement of the muon AMM, although in other cases they increase the disagreement, depending on their sign. We refer to the recent review [506] for a guide regarding new physics contributions to the muon AMM.

Regarding lepton EDMs, some examples in radiative neutrino models are given in Borah and Dasgupta [402, 397], and Chiang et al. [206]. In this case one requires CP-violating new physics in the lepton sector, something that is easily accommodated in new Yukawa couplings.

4.4. Neutrinoless double beta decay

One of the main experimental probes to test the Majorana/Dirac nature of neutrinos is neutrinoless double beta decay (0νββ), in which a nucleus (A, Z) decays into another nucleus (A, Z + 2) and two electrons [507]. In order to have sizable 0νββ rates, the nuclei should not have single beta decays. This is achieved with even-even nuclei which, thanks to the nuclear pairing force, are lighter than the odd-odd nucleus, making single beta decays kinematically forbidden. The current strongest experimental limits are obtained using 136Xe by EXO-200 [508] and KamLAND-Zen [509, 510] which yield lower bounds of the lifetime of 1.1·1025 y and 1.9·1025 y at 90 % C.L, respectively. Uncertainties in the nuclear matrix elements translate into uncertainties in the extracted values of |mee| (see Equation 6), whose current strongest upper limits are in the ballpark of ~0.15 eV. For further details regarding the present and future experimental situations see Dell'Oro et al. [511].

The observation of 0νββ decay would imply that lepton number is violated by two units (ΔL = 2), and therefore that neutrinos are Majorana particles [512]. However, quantitatively, this contribution to neutrino masses occurs at 4-loop order and is therefore extremely suppressed, much lighter than the observed neutrino masses (see Duerr et al. [513] for a quantitative study of this statement). So, even if it is true that neutrinos will necessarily be Majorana if 0νββ is observed, the main contribution to their masses may no be necessarily related to 0νββ.

We will mainly focus in this section on radiative models which have new direct contributions to neutrinoless double beta decay beyond the standard ones mediated by the light Majorana neutrinos, which are indirect, as they are generated by the new particles at higher-loop order (via light neutrino masses). For general reviews on the subject the interested reader is referred to Rodejohann [514], Bilenky and Giunti [515] and Dell'Oro et al. [511].

In Päs et al. [516, 517] a general phenomenological formula for the process including both long and short-range interactions was given. The authors considered all possible Lorentz structures for the quarks involved in the process and the outgoing electrons. In del Aguila et al. [107] effective operators that involve gauge bosons were considered, such that there are new effective vertices of the W-boson and the electrons.

In Figure 12 (reproduced from del Aguila et al. [349]) all possible contributions to 0νββ are shown, with the red dot representing the ΔL = 2 vertex. Figure 12A shows the light neutrino contribution, while Figure 12F involves a dimension-9 effective operator. In Bonnet et al. [106] a systematic classification of possible UV models stemming from the dimension 9 operator was performed (Figure 12F). See also Helo et al. [182] for scalar-mediated UV completions and its connection to neutrino masses. Figures 12D,E involve new vertices between quarks, leptons and gauge bosons. Figures 12B,C involve new vertices with just leptons and gauge bosons and no quarks. In del Aguila et al. [107] operators that involve gauge bosons were considered, such that there are new effective vertices of the W-boson and the electrons, as in Figures 12B,C. See section 3.4 for a discussion of the effective operators that generate the latter diagrams and their connection to neutrino masses. A systematic classification of UV models for all the dimension-7 operators was given in Helo et al. [108]. Many (if not all) of these particles can be present in radiative neutrino mass models.

We outline in the following two typical new contributions to 0νββ from radiative neutrino mass models:

  • New particles that couple to quarks. For instance, leptoquarks as in Hirsch et al. [518, 519] and Kohda et al. [298]. In R-parity violating SUSY (see section 5.5) there can be new contributions to 0νββ from new states, see Mohapatra [520], Babu and Mohapatra [521], and Hirsch et al. [522, 523]. Another simple example due to exchange of color octet scalars and fermions that couple to quarks and leptons simultaneously is given in Choubey et al. [524]. A model with two scalar diquarks, a dilepton and a second Higgs doublet is given in Brahmachari and Ma [525]. See other examples in Gu [526] and Helo et al. [182].

  • New particles that open operators that involve gauge bosons del Aguila et al. [107], see discussion in section 3.4.

Let us also mention that, in addition to 0νββ, there are also limits on other lepton number violating elements mαβ of the neutrino mass matrix in flavor basis (where the charged lepton mass matrix is diagonal), different from the mee (which equals mββ) one, stemming from meson decays, tau decays, e+p collider data among other processes [527]. Also indirect bounds using neutrino oscillations and the unitarity of the PMNS matrix can be set [528]. However, both the direct and indirect (even if much stronger than the direct) bounds obtained are typically very weak [527, 528]. μe+ conversion also offers a possibility to test the m element, however typically the rates are not competitive with those of 0νββ, although of course they test a different element and flavor effects could be relevant. A study of the contributions from effective operators was performed in Berryman et al. [529], while a doubly-charged scalar was studied in detail in Geib and Merle [530].

Lepton number violation can also be searched for at colliders. This is specially interesting for channels that do not involve electrons, as it is necessarily the case for 0νββ. Those will be discussed in section 4.5. Also the connection of lepton number violation to the matter-antimatter asymmetry of the universe will be discussed in section 4.6.

4.5. Collider searches

Radiative neutrino mass models generally have a much lower UV scale than the GUT scale, which makes them testable at either current or future colliders. The diversity of exotic particles and their interaction with the SM particles in radiative neutrino mass models leads to an extremely rich phenomenology at colliders. Processes pertaining to the Majorana nature of neutrino masses or LFV couplings between the exotic particles and the SM, i.e., processes violating lepton number and/or lepton flavor, are often chosen as signal regions in collider searches due to the low SM background26. Of course, there are searches for exotic particles in general if they are not too heavy and the couplings are sizable 27. In the following, we sketch different search strategies at colliders, which often utilize the low SM background for LNV and LFV processes. We thus discuss LNV and LFV processes separately before discussing general searches for new particles, which rely on processes without any LNV/LFV.

4.5.1. Lepton number violation

At the LHC, the most sought-after channel of LNV28 are same-sign leptons

where ℓ denotes e or μ, and X can be any number of jets, or other SM objects. The details of the production and the actual content of X are very model-dependent: typically heavy states are produced and decay to final states with same-sign dilepton due to their Majorana nature. We will take a doubly-charged scalar as a simple example to illustrate the basics of this search strategy. A doubly-charged scalar ϕ++ is an SU(2)L singlet with hypercharge Y = 2. They can be pair-produced via Drell-Yan process and subsequently decay to two same-sign dileptons. For large masses the photon-initiated process becomes important and leads to an enhancement [532]. Assuming the branching fraction of ϕ++e+e+ is 100%, the signature for pair-produced doubly-charged scalars is four electrons and thus ZZ production is the main SM background. To reduce the SM background, discriminating variables such as the same-sign dilepton mass, the difference between the opposite sign dilepton mass and the Z boson mass, and the scalar sum of the lepton pT can be utilized. ATLAS [533] has excluded doubly-charged SU(2)L singlet scalar with mass lower than 420 GeV at 95% CL with LHC Run 2 data. The improved limit can be extracted from the CMS search for doubly-charged component of an SU(2)L triplet [534]. In Sugiyama et al. [535], del Aguila and Chala [536, 537] and Kanemura et al. [538] studies of doubly-charged scalars and how to discriminate the multiplet to which they belong were performed.

The sensitivities of 0νββ searches detailed in section 4.4 and the same-sign dilepton searches at the LHC can be compared in any specific model (see for example [539541]). Specifically in Helo et al. [539] and Peng et al. [540] a simplified model with a scalar doublet S ~ (1, 2, 1) and a Majorana fermion F, which has the same matter content as the scotogenic model, is adopted. In this model, the reach of tonne-scale 0νββ generally beats that of the LHC. In the parameter space region where the heavy particle masses are near the TeV scale, however, the two probes are complementary.

4.5.2. Lepton flavor violation

As described in section 4.2, lepton flavor violating processes are commonly predicted in radiative neutrino mass models, which can also be probed at colliders. The actual production topology of the LFV processes varies from model to model. For example, in models with the leptoquark , there are two possible decay channels, or [542]. The dilepton final states are produced from

where X can represent , multiple jets and leptons, and the former contributes dominantly for normal ordering in the minimal model with two leptoquarks. SUSY stop searches in the dilepton final states have the same signatures and their collider bounds can be translated into that of the leptoquark. This has been done for the LHC 8 TeV run [543] and the limit was mS1 ≳ 600 GeV [214]. Note that this limit in LFV channel is stronger than lepton flavor conserving ones (mS1 ≳ 500 GeV) as the SM background is lower. The stop search has been updated for LHC Run 2 [544, 545], though a recast for leptoquarks in LFV dilepton final states still awaits further analysis.

Alternatively, LFV processes can also be studied in an independent manner. In the framework of effective operators with two flavor-diagonal quarks and two flavor-off-diagonal leptons, constraints from LHC searches for LFV final states are interpreted as lower limits on the UV cut-off scale [546]. Compared with the limits derived from low energy precision measurements [449, 546], LHC delivers less stringent limits for light quarks. For heavier quarks, however, competitive limits of ΛUV ≳ 600−800 GeV can already be set for operators with right-handed τ leptons using only LHC Run 1 data.

4.5.3. Searches for new particles

Radiative neutrino mass models may contain exotic particles such as vector-like quarks (VLQs), vector-like leptons (VLLs), scalar leptoquarks, singly- or doubly-charged scalars, colored octet fermions or scalars, and electroweak multiplets. Note that the examples here are far from complete and searches for each individual particle require their own dedicated discussion. In Cai et al. [214], LHC searches for exotic particles in UV complete models based on ΔL = 2 dimension 7 operators are discussed systematically. Here we will only present a simple summary about a handful of new particles.

Vector-like quarks

We refer by VLQs to new SU(3)c triplets which mix with the SM quarks and Higgs via Yukawa couplings [547]. The VLQs include different SU(2) representations: two singlets T and B with hypercharge 2/3 and −1/3; three doublets (T, B), (X, T), and (B, Y) with hypercharge 1/6, 7/6, and −5/6; and two triplets (X, T, B), and (T, B, Y) with hypercharge 2/3 and −1/3. They can be pair produced at the LHC via gluon fusion and quark-antiquark annihilations. Single production is model-dependent and can be dominant for large vector-like quark masses and large mixings [547]. The mass splitting among the components of the fields is suppressed by the mixing angles between the SM quarks and the vector-like quarks, which in turn suppresses the decays between the component fields. Therefore, VLQs will dominantly decay to either a gauge boson or a Higgs plus a SM quark. Both ATLAS and CMS have performed searches for VLQs and have set lower limits on the VLQs masses up to 990 GeV at the 95% confidence level (CL) depending on the representations and the decay branching ratio [548559].

Vector-like leptons

VLLs are the colorless version of VLQs. Similar to VLQs, VLLs mix with the SM leptons via Yukawa couplings with Higgs. Due to the absence of right-handed neutrinos, there are less VLLs: two singlets N and E with hypercharge 0 and 1; two doublets (N, E) and (E, D) with hypercharge 3/2 and 1/2; and triplets (P, N, E) and (N, E, D) with hypercharge 0 and 1, respectively. Detailed studies have been performed in Altmannshofer et al. [560], Falkowski et al. [561], Dermisek et al. [562], and Kumar and Martin [563]. Contrary to the colored VLQs, VLLs are dominantly pair produced at the LHC via Drell-Yan process as the phase space suppression is less significant in the parameter space of interest at the moment. They can also be singly produced in association with W, Z or H, which can be dominant if the pair production channel is phase space suppressed and sizable mixing parameters are assumed. Likewise VLLs decay either to a SM lepton and a boson, W or Z, or Higgs. So far there is no dedicated search for VLLs at colliders, though SUSY searches for sleptons or charginos can be used to derive bounds on VLLs (see Altmannshofer et al. [560] and Hamada et al.[564] for example).

Leptoquarks

Leptoquarks appear frequently in theories beyond the SM such as grand unified theories [565, 566]. As its name suggests, a leptoquark, which can be either a scalar or a vector [542], possesses both non-zero lepton and baryon numbers. Here we will focus on scalar leptoquarks. At hadron colliders, leptoquarks are primarily produced in pairs via gluon fusion and quark-antiquark annihilation. Each leptoquark subsequently decays to one quark and one charged or neutral lepton. Both ATLAS [567, 568] and CMS [569571] have performed searches for leptoquarks in final states with two charged leptons plus multiple jets. Assuming 100% branching fraction of the leptoquark decay into a charged lepton and a quark, current searches at the LHC Run 2 with 13 TeV center of mass energy have excluded leptoquarks with masses less than 1,130 GeV [569], 1,165 GeV [570] and 900 GeV [571] at 95% CL for leptoquark couplings to the first, second and third generations respectively.

Charged scalars singlets

Singly- and doubly-charged scalars are introduced in various radiative neutrino mass models (see Babu and Julio [158], Zee [288], and Babu [289], for instance). As singlets under SU(3)c × SU(2)L, the singly (doubly) charged scalar can only couple to the lepton doublet (right-handed charged lepton) bilinear. So the doubly-charged scalar can only decay to a pair of charged leptons, which leads to LNV signature at colliders (see discussion in section 4.5.1 for details). As for the singly-charged scalar, it decays to a charged lepton and a neutrino whose LNV effects can not be detected at the LHC. Singly-charged scalars are mainly produced in pairs via the Drell-Yan pair process. They are searched for in final states with two leptons plus 29. SUSY searches for sleptons and charginos at the LHC share the same signature as the singly-charged scalars. Thus, we can in principle recast the slepton search in The ATLAS Collaboration [574] and extract the limit for our singly-charged scalars. Note a slepton can also be produced via a W-boson, while singly-charged scalar only via a virtual photon.

Higher-dimensional electroweak multiplet

SU(2)L higher-dimensional representations can also be incorporated in radiative neutrino mass theories [126129, 362, 473, 575, 576]. While the mass splittings among the component fields for scalar multiplets can be generally large due to couplings to the SM Higgs, those for fermion multiplets are only generated radiatively and are typically MeV, with the neutral component being the lightest. This small mass splitting results in lifetimes ns. At the LHC, charged component field can be produced in pair via electroweak interaction and decay to the neutral component plus a very soft pion, which leads to a disappearing track signature. For a triplet with a lifetime of about 0.2 ns, the current LHC searches set the lower mass limit to be 430 GeV at 95% CL [577579].

4.6. Generation of the matter-antimatter asymmetry of the universe

The matter-antimatter asymmetry of the universe has been inferred independently (and consistently) by big bang nucleosynthesis (BBN) predictions of light elements, and by the temperature anisotropies of the cosmic microwave background. In order to generate it, the Sakharov conditions need to be fulfilled [580]. There should be:

  • Processes that involve baryon number violation (BNV).

  • Processes in which both charge conjugation (C) and charge and parity conjugation (CP) are violated.

  • Departure from thermal equilibrium, so that (i) the number densities of particles and antiparticles can be different, and (ii) the generated baryon number is not erased.

In the standard model, it is well-known that due to the chiral nature of weak interactions B+L is violated by sphaleron processes, while B−L is preserved [581]. Also C and CP are violated in the quark sector (in the CKM matrix), although the amount is too small to generate the required CP asymmetry. In the lepton sector (with massive neutrinos) CP can be violated, and there are in fact hints of δ ~ −π/2 [20]. However, the measurement of the Higgs mass at 125 GeV implies that the phase transition is not strongly first-order, with no departure from thermal equilibrium. Therefore, the SM has to be extended to explain the matter-antimatter asymmetry which raises the question whether this new physics is related to neutrino masses or not.

When sphalerons are active and in thermal equilibrium, roughly at temperatures above the electroweak phase transition, B+L can be efficiently violated. Therefore, one natural option in models of Majorana neutrinos is that an asymmetry in lepton number is generated, which is converted by sphalerons into a baryon asymmetry. This is known as leptogenesis [25] (see Davidson et al. [582] for a review on the topic), the most popular example being the case of type-I seesaw, where the out-of-equilibrium decays of the lightest of the heavy right-handed neutrinos into lepton and Higgs doublets and their conjugates, at a temperature equal or smaller than its mass, generate the lepton asymmetry due to CP-violating interactions.

The scotogenic model and its variants, see section 5.3, have been studied in detail regarding the generation of the baryon asymmetry from particle decays with TeV-scale masses. Ma [583] briefly discusses leptogenesis within the scotogenic model. This discussion is extended in Kashiwase and Suematsu [584, 585] and Racker [586] to include resonant leptogenesis. Resonant leptogenesis has also been studied in a gauge extension of the scotogenic model [224226] in Kashiwase and Suematsu [252] and resonant baryogenesis in an extension with new colored states in Dev and Mohapatra [587]. Hambye et al. [588] and Babu and Ma [589] consider extensions of the scotogenic model by an additional charged or neutral scalar to achieve viable non-resonant leptogenesis. The baryon asymmetry can similarly be enhanced by producing the SM singlet fermions in the scotogenic model non-thermally beyond the usual thermal abundance [590]. Leptogenesis via decays of an inert Higgs doublet or a heavy Dirac fermion were studied in Lu and Gu [119, 154] in scotogenic-like models, respectively. In Chen and Law [127] leptogenesis was studied in a scotogenic-like model with fermionic 5-plets and a scalar 6-plet, via the decays of the second-lightest fermionic 5-plet. Baldes et al. [278] demonstrated the feasibility to generate the correct matter-antimatter asymmetry via leptogenesis in the model proposed in Ma [335]. It also showed that any pre-existing baryon asymmetry in the two models proposed in Ma [335] and Law and McDonald [277] is washed out at temperatures above the mass of their heaviest fields.

In radiative models with extra scalars coupled to the Higgs field, the phase transition can generally be stronger, as they contribute positively to the beta function of the Higgs and therefore, they help to stabilize the Higgs potential. Moreover, in these models there are typically extra sources of CP violation. These two ingredients allow the possibility of having electroweak baryogenesis. In particular, the strong first-order phase transition has been discussed using an effective potential in Bertolini et al. [591], and in Aoki et al. [592] for the model of Aoki et al. [370]. Also in the case of a supersymmetric radiative model in Kanemura et al. [202].

However, in general the new states can also destroy a pre-existing asymmetry, irrespective of their production mechanism, as they violate necessarily lepton number by two units [593596]. The new particles typically have gauge interactions, so that they are in thermal equilibrium at lower temperatures than those at which the asymmetry is generated (by high-scale baryogenesis or by leptogenesis, for instance30) potentially washing it out.

Some works have focused on the fact that if LNV is observed at the LHC, one could falsify leptogenesis, as the wash-out processes would be too large [597599]. Similarly, observations of 0νββ rates beyond the one generated by the light neutrinos could impose constraints for the first family [600]. LFV processes could be used to extend it to all families. See Deppisch et al. [531] for further discussions about LNV processes in leptogenesis.

The limits on radiative models due to the requirement of not washing-out any pre-existing asymmetry are model-dependent. A more systematic way to go is to consider the LNV effective operators related to radiative models [81, 85, 87]. These operators lead to wash-out processes if they are in thermal equilibrium above the electroweak phase transition, and therefore their strength can be bounded by this requirement.

4.7. A possible connection to dark matter models

In many radiative neutrino mass models the generation of neutrino masses at tree-level is forbidden by a symmetry, . This symmetry can be global or gauge, continuous or discrete (a typical example is a Z2 parity), imposed or accidental (a by-product of other symmetries in the model). If is preserved after electroweak symmetry breaking, the lightest state transforming non-trivially under it, the so-called lightest charged particle (LCP), is completely stable and, in principle, could constitute the dark matter (DM) of the universe. This opens up an interesting connection between radiative neutrino masses and dark matter. DM may be produced via its coupling to neutrinos and thus the annihilation cross section is closely related to neutrino mass. This has been studied using an effective Lagrangian for light, MeV-scale, scalar DM [601] in a scotogenic-like model and for fermionic DM [460, 468472] in the scotogenic model. A key signature of this close connection is a neutrino line from DM annihilation. The constraints from neutrino mass generation on the detectability of a neutrino line has been recently discussed in El Aisati et al. [602].

Based on the general classification of 1-loop models [100], the authors of Restrepo et al. [111] performed a systematic study for models compatible with DM stabilized by a discrete Z2 symmetry. They focused on the topologies T1-x and T3. The topologies T4-2-i and T4-3-i require an additional symmetry to forbid the tree-level contribution and thus were not studied in Restrepo et al. [111]. A similar classification for 2-loop models has been presented in Simoes and Wegman [112] based on the possible 2-loop topologies discussed in Aristizabal Sierra et al. [102]. Symmetries forbidding tree and lower-order loop diagrams have been discussed in Farzan et al. [110]. In section 5.3 we discuss the prototype example of such models: the scotogenic model.

Besides dark matter being stabilized by a fundamental symmetry, it may be stable due to an accidental symmetry. For example, higher representations of SU(2)L cannot couple to the SM in a renormalizable theory, which leads to an accidental Z2 symmetry at the renormalizable level. This has been dubbed minimal dark matter [385, 386]. After the initial proposal to connect the minimal dark matter paradigm and radiative neutrino mass generation [126], it has been conclusively demonstrated that the minimal dark matter paradigm cannot be realized in 1-loop neutrino mass models [473, 575, 576]. However, there is a viable variant of the KNT model at 3-loop order [362], which realizes the minimal dark matter paradigm without imposing any additional symmetry beyond the SM gauge symmetry.

Finally, the DM abundance in the universe may be explained by a light pseudo-Goldstone boson (pGB) associated with the spontaneous breaking of a global symmetry. It is commonly called Majoron in case the lepton number plays the role of the global symmetry. The possibility of pGB dark matter has been discussed in one of the models in Dasgupta et al. [219] which provides a pGB dark matter candidate after the breaking of a continuous U(1) symmetry to its Z2 subgroup in addition to the LCP. Recently the authors of Ma et al. [603] proposed an extension of the Fileviez-Wise model [120] to incorporate a Majoron DM candidate which simultaneously solves the strong CP problem.

5. Selected examples of models

In the following subsections, we list and discuss different benchmark models for neutrino mass that are qualitatively different. We start with the most well-studied models, which are the Zee model, discussed in section 5.1.1, that is the first 1-loop model for Majorana neutrino masses, and the Zee-Babu model, revisited in section 5.1.2, which is the first 2-loop model. In section 5.2 we discuss the first 3-loop model [355], which was proposed by Krauss, Nasri, and Trodden and is commonly called KNT-model, and its variants. It is also the first model with a stable dark matter candidate. The scotogenic model is discussed in section 5.3. It generates neutrino mass at 1-loop order and similarly to the KNT-model it features a stable dark matter candidate due to the imposed Z2 symmetry. These are the most well-studied models in the literature. However, this preference is mostly due to the historic development (and also simplicity) and we are proposing a few other interesting benchmark models in the following subsections.

5.1. Models with leptophillic particles

There are only three different structures which violate lepton number (LN) by two units that can be constructed with SM fields [76]:

The three different structures can couple respectively to a SU(2) triplet scalar with Y = 1 (we denote it by Δ), a singly-charged SU(2) singlet scalar (we call it h+) and a doubly-charged SU(2) singlet scalar (we call it k++).

In all cases, we could assign LN equal to −2 to the new fields so that such interactions preserve it. However, dimension-3 terms in the scalar potential will softly break LN, as there is no symmetry to prevent them. In the first case, the triplet can have in the potential the lepton-number violating term (with ΔL = 2) with the SM Higgs doublet H

Then, after electroweak symmetry breaking, the triplet gets an induced VEV (strongly bounded by the T parameter to be ≲O(1) GeV), and neutrino masses are generated at tree-level via the type-II seesaw.

If only the singly-charged scalar h+ is present, a ΔL = 2 term can be constructed with two Higgs doublets, the SM Higgs H and an extra Higgs doublet Φ

In this case, however, neutrino masses are not induced by the Higgs VEV at tree-level, but they are generated at 1-loop order. This is known as the Zee model [104, 105].

For the case of the doubly-charged scalar, one can construct the ΔL = 2 term precisely with two singly-charged scalars h+

Notice that no other combination with SM fields exist, given the large electric charge of k++. In this case, neutrino masses are generated at 2-loop order. This is known as the Zee-Babu model [76, 604].

These are the simplest radiative models. By using particles that couple to a lepton and a quark (leptoquarks), one can also have ΔL = 2 interactions and generate neutrino masses at a different number of loops. In the following, we will discuss the Zee and Zee-Babu models.

5.1.1. The zee model

In addition to the SM content with a Higgs scalar doublet H, the Zee model [104, 105] contains an extra Higgs scalar doublet Φ and a singly-charged scalar singlet h+, which is shown in Table 6. It is an example of the operator . Several aspects of the phenomenology of the model have been studied in Petcov [605], Zee [288], Bertolini and Santamaria [606, 607], Yu et al. [608, 609], Frampton and Glashow [610], Jarlskog et al. [611], Ghosal et al. [451], Kanemura et al. [612], Balaji et al. [613], Koide [614], Brahmachari and Choubey [615], Frampton et al. [616], Assamagan et al. [617], He [156], Kanemura et al. [618], and Aristizabal Sierra and Restrepo [619]. While the Zee-Wolfenstein version where just the SM Higgs doublet couples to the leptons has been excluded by neutrino oscillation data [155, 156], the most general version of the Zee model in which both couple remains allowed [157] and has been recently studied in Herrero-García et al. [91] (see also Babu and Julio [158] and Aranda et al. [159] for a variant with a flavor-dependent Z4 symmetry).

Table 6

Field Spin G SM
h + 0 (1, 1, 1)
Φ 0

Quantum numbers for new particles in the Zee model.

The Yukawa Lagrangian is

where and eR are the SU(2) lepton doublets and singlets, respectively, and with τ2 being the second Pauli matrix. Due to Fermi statistics, f is an antisymmetric Yukawa matrix in flavor space, while Y1 and Y2 are completely general complex Yukawa matrices. Furthermore, the charged-lepton mass matrix is given by

where tan β = sβ/cβ = v2/v1 with and and . Without loss of generality, one can work in the basis where mE is diagonal.

Assuming CP-invariance there are two CP-even neutral scalars (one of which is the 125 GeV Higgs boson, with mass mh, and the other is a heavy one with mass mH), one neutral CP-odd scalar with mass mA, and two charged-scalars of masses , whose mixing due to the trilinear term in Equation (81) is given by

Interestingly, μZee cannot be arbitrarily large, as it contributes at 1-loop level to the mass of the light Higgs. Demanding no fine-tuning, we can estimate |μZee| ≲ 4π mh ≃ 1.5 TeV.

The Yukawa couplings of Equation (83), together with the term in the potential given in Equation (81), imply that lepton number is violated by the product mE(Y1v2Y2v1)Zee. Therefore, neutrino masses will be necessarily generated, in particular the lowest order contribution appears at 1-loop order, as shown diagram of Figure 17, where the charged scalars run in the loop. The neutrino mass matrix is given by:

with φ being the mixing angle for the charged scalars given in Equation (85). Therefore, in the Zee model, due to the loop and the chiral suppressions, the new physics scale can be light. From the form of the mass matrix it is clear that if one takes Y2 → 0 (Zee-Wolfenstein model), the diagonal elements vanish, yielding neutrino mixing angles that are not compatible with observations.

Figure 17

Figure 17

1-loop neutrino masses generated in the Zee model in the flavor basis.

Neglecting memμ, mτ and taking f = 0, the following Majorana mass matrix is obtained

Notice that if the term proportional to the muon mass is neglected, one neutrino remains massless. In order to obtain correct mixing angles, we need both and different from zero [91, 491], as they enter in the 1-2 submatrix of Equation (87). This implies that LFV mediated by the scalars will be induced. In fact, in the model large LFV signals are generated, like τ → μγ and μ − e conversion in nuclei. Moreover, also a full numerical scan of the model performed in in Herrero-García et al. [91] showed that large LFV Higgs decays are possible, in particular BR(h → τμ) can reach the percent level. BR(h → τe) is roughly two-orders of magnitude smaller than BR(h → τμ). The singly-charged h also generates violations of universality, as it interferes constructively with the W boson, as well as non-standard interactions, see section 4.1, which however are too small to be observed [91].

In Herrero-García et al. [91] it was also shown that the model is testable in next-generation experiments. While normal mass ordering (NO) provided a good fit, inverted mass ordering (IO) is disfavored, and if θ23 happens to be in the second octant, then IO will be ruled-out. Notice also that the lightest neutrino is required to be massless for IO, as it has also been obtained in He and Majee [157]. Furthermore, future τ → μγ (μ − e conversion) will test most regions of the parameter space in NO (IO). Regarding direct searches at the LHC, the new scalars have to be below ~2 TeV, which implies that they can be searched for similarly as in a two-Higgs doublet model (with an extra charged scalar that could be much heavier). Particularly, the charged scalars are searched for at colliders. See the discussion in section 4.5.

Let us mention that an interesting modification of the Zee model was proposed in Babu and Julio [158] (see also Aranda et al. [159]), where a Z4 symmetry was imposed, being able to reduce significantly the number of parameters. In that case, among the predictions of the model, is that the spectrum should be inverted. Other flavor symmetries beyond Z4 in this framework have been studied in Babu and Mohapatra [168, 169], Koide and Ghosal [170], Kitabayashi and Yasue, [171], Adhikary et al. [172], Fukuyama et al. [173], Aranda et al. [174, 175].

5.1.2. The zee-babu model

The Zee-Babu model contains, in addition to the SM, two SU(2) singlet scalar fields with electric charges one and two, denoted by h+ and k++ [76, 604] as shown in Table 7. It is a UV completion of the operator . Several studies of its phenomenology exist in the literature [95, 435, 620622].

Table 7

Field Spin G SM
h + 0 (1, 1, 1)
k ++ 0 (1, 1, 2)

Quantum numbers for new particles in the Zee-Babu model.

The leptonic Yukawa Lagrangian reads:

where like in the Zee model, due to Fermi statistics, f is an antisymmetric matrix in flavor space. On the other hand, g is symmetric. Charged lepton masses are given by , which be take to be diagonal without loss of generality.

Lepton number is violated by the simultaneous presence of the trilinear term μZB in Equation (82), together with mE, f, g. Note that the trilinear term cannot be arbitrarily large, as it contributes to the charged scalar masses at loop level, and can also lead to charge-breaking minima, if |μZB| is large compared to the charged scalar masses. For naturalness considerations we demand |μZB| ≪ 4π min(mh, mk). See Nebot et al. [435] and Herrero-Garcia et al. [95] for detailed discussions.

As lepton number is not protected, neutrino masses are generated radiatively, in particular at 2-loop order, via the diagram of Figure 18. The mass matrix is approximately given by (see for instance McDonald and McKellar [286], Nebot et al. [435] and Herrero-Garcia et al. [95] for more details)

where M is the heaviest mass of the loop, either that of the singly-charged singlet h+ or of the doubly-charged singlet k++. A prediction of the model is that, since f is a 3 × 3 antisymmetric matrix, det f = 0, and therefore . Thus, at least one of the neutrinos is exactly massless at this order.

Figure 18

Figure 18

2-loop neutrino masses generated in the Zee-Babu model.

In the model, both NO and IO can be accommodated. The phenomenology of the singly-charged scalar is similar to that discussed in the Zee model, apart from the fact that in the Zee model the charged singlet mixes with the charged component of the doublet. Some of the most important predictions of the model are due to the presence of the doubly-charged scalar k++. Firstly, k++ mediates trilepton decays () at tree-level which unlike, in the Zee model, are not suppressed by the small charged lepton masses, as well as radiative decays (ℓi → ℓjγ). Secondly, k++ can be pair-produced at the LHC via Drell-Yan, decaying among other final states into same-sign leptons which yields a clean experimental signature. See the discussion in section 4.5.

5.2. KNT-models

The first radiative neutrino mass model at 3-loop order is the KNT model [355] which has one fermionic singlet N and two singly-charged scalars S1,2 in addition to the SM particles. A discrete Z2 symmetry is imposed, under which only S2 and N are odd. We list the quantum numbers of the exotic particles in Table 8.

Table 8

Field Spin G SM Z 2
S 1 0 (1, 1, −1) +
S 2 0 (1, 1, −1)
N (1, 1, 0)

Quantum numbers for new particles in the original KNT model.

The Z2 symmetry forbids the usual type-I seesaw contribution at tree-level. The relevant Lagrangian is expressed as

where the flavor indices of f and g are all suppressed. With this setup, neutrino masses are generated first at 3-loop order as shown in Figure 19. The neutrino mass matrix is then

where the function F is defined in Ahriche and Nasri [358]. This matrix is, however, only rank one and thus can give exactly one non-zero neutrino mass. Adding more copies of N can increase the rank of the matrix. The phenomenology of this model including flavor physics, dark matter, Higgs decay, electroweak phase transition and collider searches is discussed in detail in Ahriche and Nasri [358].

Figure 19

Figure 19

3-loop neutrino masses generated in the KNT model.

This model is subject to constraints from LFV experiments such as μ → which requires three copies of N for the neutrino mixing to be in agreement with the observations31. Meanwhile in order to be consistent with the measurements of muon anomalous magnetic moment and the 0νββ decay, strong constraints are imposed. For MS1,S2 ⩾ 100 GeV, and , it can satisfy all flavor constraints while reproducing the neutrino mixing data.

Assuming a mass hierarchy MN < MS2, the lightest fermion singlet is stable and serves as a good DM candidate. This is also the first radiative neutrino mass theory with a stable DM candidate running in the loop. If the DM relic density is saturated and all previously discussed constraints are satisfied, the DM mass cannot exceed 225 GeV while the lighter charged scalar S2 cannot be heavier than 245 GeV. If the fermion singlets have very small mass splitting, DM coannihilation effects should be taken into account. With about 5% mass splitting, the DM relic density increases by 50%.

As discussed in section 4.5.3, the singly-charged scalars can be pair-produced at the LHC and subsequently decay to a pair of charged leptons and the fermion singlets which appear as missing transverse energy. This signature is exactly the same as the direct slepton pair production in SUSY theories. ATLAS has performed the search for sleptons in this channel with 36.1 fb−1 data of TeV [574] and has ruled out slepton masses below ~500 GeV in the non-compressed region. The actual constraint on MS2 depends on the decay branching ratio of S2 to different leptons and in principle will be substantially relaxed compared to the ATLAS search.

With the same topology, a lot of variations of the KNT model can be constructed. Chen et al. [356] discusses several possibilities to replace the electron with other SM fermions32 or vector-like fermions. A similar model in which the electron is replaced by a fermion doublet with hypercharge 5/2 and S1,2 with doubly-charged scalar is discussed in Okada and Yagyu [368]. The Z2-odd particles in this model form instead the outer loop.

5.3. The scotogenic model

The most popular model linking dark matter to the radiative generation of neutrino masses is the one proposed by E. Ma in 2006. We will refer to it as scotogenic model [113]33. In the scotogenic model, the SM particle content is extended with three singlet fermions, Ni (i = 1, 2, 3), and one SU(2)L doublet, η, with hypercharge ,

This setup is supplemented with a Z2 parity, under which the new states are odd and all the SM particles are even34. The newly-introduced particles with their respective charges of the scotogenic model are shown in Table 9. The gauge and discrete symmetries of the model allow us to write the Lagrangian terms involving the fermion singlets

Table 9

Field Spin Generations G SM Z 2
η 0 1
N 3 (1, 1, 0)

Quantum numbers of new particles in the scotogenic model.

We do not write the kinetic term for the fermion singlet as it takes the standard canonical form. YN is an arbitrary 3 × 3 complex matrix, whereas the 3 × 3 Majorana mass matrix MN can be taken to be diagonal without loss of generality. We highlight that the usual neutrino Yukawa couplings with the SM Higgs doublet are not allowed due to the Z2 symmetry. This is what prevents the light neutrinos from getting a non-zero mass at tree-level. The scalar potential of the model is given by

Neutrino masses are induced at the 1-loop level via the diagram in Figure 20

where the masses of the scalar ηR and pseudo-scalar part ηI of the neutral scalar are given by

with the electroweak VEV . Neutrino mass vanishes in the limit of λ5 = 0 and thus degenerate masses for the neutral scalars ηR,I, because it is possible to define a generalized lepton number which forbids a Majorana mass term.

Figure 20

Figure 20

1-loop neutrino masses generated in the scotogenic model.

In the scotogenic model, the Z2 parity is assumed to be preserved after electroweak symmetry breaking. This will be so if 〈η〉 = 0. In this case, the lightest Z2-odd state (to be identified with the LCP defined in section 4.7) will be stable and, if neutral, will constitute a potentially good DM candidate. The LCP in the scotogenic model can be either a fermion or a scalar: the lightest singlet fermion N1 or the lightest neutral η scalar (ηR or ηI). As the neutrino Yukawa couplings are generally required to be small to satisfy LFV constraints, the DM phenomenology for a scalar LCP is generally the same as in the inert doublet model [624, 625]. Recently it has been pointed out [626] that late decay of the lightest SM singlet fermion N1 may repopulate the dark matter abundance and thus resurrect the intermediate dark matter mass window between mW, the mass of the W boson, and 550 GeV. In the case of a fermionic LCP, for which the annihilation cross section is governed by the neutrino Yukawa couplings, the connection of the dark matter abundance with neutrino masses leads to a very constrained scenario due to the bounds from lepton flavor violation [460, 468472].

Many scotogenic variations have been proposed since the publication of the minimal model described above. All these models are characterized by neutrino masses being induced by new dark sector particles running in a loop [114135, 138, 139, 141143, 204, 261, 265]. One of them involves a global continuous dark symmetry, instead of a discrete dark symmetry [145], Hagedorn, (in prep). A gauge dark symmetry was considered in Yu [253] and a scale-invariant version presented in Ahriche et al. [245]. The collider [627630] and dark matter [631634] phenomenologies of different scotogenic variants have also been discussed in detail. Finally, we point out that the authors of Merle and Platscher [635] identified a potential problem in this family of models, since some parameter regions lead to the breaking of the Z2 parity at high energies. This problem, how it can be escaped and its phenomenological implications have been explored in Merle et al. [636], Merle and Platscher [267], and Lindner et al. [637].

5.4. Models with leptoquarks

Leptoquarks are common ingredients of radiative neutrino mass models. For example neutrino mass can be generated at loop level by two leptoquarks which mix via a trilinear coupling to the SM Higgs boson [178185]. Neutrino mass generation at 1-loop order with all possible leptoquarks has been systematically studied in in Aristizabal Sierra et al. [181]. At 1-loop order and especially at a higher-loop order, leptoquarks usually appear together with other exotic particles such as vector-like quarks and leptons, charged scalar singlets and electroweak multiplets [214]. We will review two models here, one at 1-loop and one at 2-loop order.

5.4.1. A 1-loop model

Without introducing exotic fermions, the only possible topology that can contribute at 1-loop order to the Weinberg operator is T1-ii shown in Figure 3 as we need the fermion arrow to flip only once. With this topology and leptoquarks as the only exotic particles, the only UV completion we can realize is depicted in Figure 21. The relevant scalar leptoquarks35 are S1, S3 and with quantum numbers detailed in Table 10. The relevant Lagrangian reads

following the convention in in Doršner et al. [185] with all generation indices suppressed. Apparently only the leptoquark component fields with electric charge can contribute. These leptoquarks, in the interaction basis , will mix with each other through the λ1,3 terms in Equation (98)36. We will consider simplified scenarios where either S1 or S3 appears together with . For the model with S1,3, the squared-mass matrix will be diagonalized with angle θ1,3 and the mass eigenvalues are m1 and m2. So the neutrino mass matrix is expressed as [181, 185]

where Md = diag(md, ms, mb) with md,s,b being the down, strange and bottom quark masses. Due to the hierarchy of down-type quark masses, the neutrino mass matrix will be approximately rank-2 with one nearly massless neutrino. Current neutrino oscillation data put lower bounds on the product of Yukawa couplings ranging from 10−12 to 10−7 for leptoquarks with TeV scale masses [181]. On the other hand, low energy precision experiments constrain the Yukawa couplings from above. For example, μ − e conversion in titanium bounds the first generation Yukawa couplings with

for 1 TeV leptoquark masses. Their decay branching fractions are dictated by the same couplings that determine the neutrino masses and mixings, which leads to a specific connection between the decay channels of the leptoquark and the neutrino mixings. Generally LFV decays with similar branching ratios to final states with muon and tau are expected in some leptoquark decays. This neutrino mass model can also be tested at colliders. The leptoquarks running in the loop can be created in pairs and decay to final states containing leptons plus jets with predicted branching ratios. We refer to section 4.5 for further details on searches of leptoquarks at colliders.

Figure 21

Figure 21

1-loop neutrino masses generated via leptoquark mixing.

Table 10

Field Spin G SM
S 1 0
S 3 0
0

Quantum numbers of leptoquarks.

Päs and E. Schumacher [183] explored the possibility to explain the anomalous bsll transitions with S3 and . Different texture of the Yukawa coupling matrices y3 and ỹ2 were considered and leptoquark masses in the the range of 1 to 50 TeV can reproduce the neutrino masses and mixings in addition to RK [639].

5.4.2. A 2-loop model

Based on the gauge-invariant effective operator , which violates lepton number by two units, a UV complete radiative neutrino mass model at 2-loop order containing leptoquark S1 and fermion color octet f can be constructed [287]. We list their quantum numbers in Table 11 for the convenience of the readers.

Table 11

Field Spin Generation G SM B
S 1 0 2 −1
f 1 (8, 1, 0) 0

Quantum numbers of new particles in the Angelic model.

The general gauge invariant Lagrangian for the exotic particles is then expressed as

where generation indices for all parameters and fields are suppressed. We demand baryon number conservation to forbid the terms and which induce proton decay. With this setup, Majorana neutrino mass will be generated at 2-loop order as shown in Figure 22. Generally the contribution to the neutrino mass matrix is proportional to the down-type Yukawa coupling squared which is dominated by the third generation unless strong hierarchy in λLQλdf exists. As a result, we can simplify the formula for the neutrino mass matrix to

with the CKM-matrix element Vtb and Iαβ as a function of mf and mS1 whose exact form can be read from Angel et al. [287]. The indices α and β label the leptoquark copies. This neutrino mass matrix is only rank one if there is only one leptoquark flavor assuming the dominance of the bottom-quark loop37. At least two leptoquarks are needed to fit to the current neutrino oscillation data in this model, where one neutrino mass eigenvalue is nearly vanishing. Among all flavor processes, μ − e conversion in nuclei, μ → and μ → eee give the most stringent constraints.

Figure 22

Figure 22

2-loop neutrino masses generated in the Angelic model.

The leptoquark S1 can explain the recent anomalies observed in semileptonic B decays, i.e., the violation of lepton flavor universality (LFU) of [639] and [640645]. In the parameter space with relatively large , the combination of left- and right-handed couplings induces scalar and tensor operators, which lift the chirality suppression of the semi-leptonic B-decay BD(*)ℓν and produce sizable effects in the LFU observables [313].

5.5. Supersymmetric models with R-parity violation

Supersymmetric models with R-parity violation naturally lead to non-zero neutrino masses and mixings. These models have been regarded as very economical, since no new superfields besides those already present in the MSSM are required. Moreover, their phenomenology clearly departs from the standard phenomenology in the usual SUSY models, typically providing new experimental probes.

With the MSSM particle content, one can write the following superpotential, invariant under supersymmetry, as well as the gauge and Lorentz symmetries,

Here is the MSSM superpotential, whereas

The ϵ coupling has dimensions of mass, {i, j, k} denote flavor indices and gauge indices have been omitted for the sake of clarity. The first three terms in break lepton number (L) whereas the last one breaks baryon number (B). The non-observation of processes violating these symmetries impose strong constraints on these parameters, which are required to be rather small [646]. Also importantly, their simultaneous presence would lead to proton decay, a process that has never been observed and whose rate has been constrained to increasingly small numbers along the years. For this reason, it is common to forbid the couplings in Equation (104) by introducing a discrete symmetry called R-parity. The R-parity of a particle is defined as

where s is the spin of the particle. With this definition, all SM particles have Rp = +1 while their superpartners have Rp = −1, and the four terms in are forbidden. Furthermore, as a side effect, the lightest supersymmetric particle (LSP) becomes stable and can be a dark matter candidate.

However, there is no fundamental reason to forbid all four couplings in . When R-parity is conserved both lepton and baryon numbers are conserved, but in order to prevent proton decay just one these two symmetries suffices. Furthermore, the breaking of R-parity by L-violating couplings generates non-zero neutrino masses, and thus constitutes a well-motivated scenario beyond the standard SUSY models. This scenario (with only L-violating couplings) can be theoretically justified by replacing R-parity by a less restrictive symmetry, such as baryon triality [647].

We can distinguish two types of R-parity violating (RPV) neutrino mass models:

  • Bilinear R-parity violation (b-p): In this case the only RPV term in the superpotential is the bilinear , which breaks lepton number by one unit. This leads to the generation of one mass scale for the light neutrinos at tree-level via a low-scale seesaw mechanism with the neutralinos playing the role of the right-handed neutrinos. The second (necessary) mass scale is induced at the 1-loop level. Therefore, this can be regarded as a hybrid radiative neutrino mass model.

  • Trilinear R-parity violation (t-p): When one allows for the violation of R-parity with the trilinear superpotential terms , lepton number is also broken by one unit and Majorana neutrino masses are generated at the 1-loop level. Therefore, this setup constitutes a pure radiative neutrino mass scenario.

We now proceed to discuss some of the central features of these two types of leptonic RPV models, highlighting the most remarkable experimental predictions. Although in general one can have both types of leptonic RPV simultaneously, we will discuss them separately for the sake of clarity.

Neutrino masses with b-p

Bilinear R-parity violation [177] is arguably the most economical supersymmetric scenario for neutrino masses. The bilinear ϵi = (ϵe, ϵμ, ϵτ) terms in the superpotential come along with new terms in the soft SUSY breaking potential. Therefore, the number of new parameters in b-p with respect to the MSSM is 6, without modifying its particle content, and they suffice to accommodate all neutrino oscillation data. For a comprehensive review on b-p see Hirsch and Valle [648].

The ϵi couplings induce mixing between the neutrinos and the MSSM neutralinos. In the basis , the neutral fermion mass matrix is given by

Here is the standard MSSM neutralino mass matrix and m ∝ ϵ is the matrix containing the neutrino-neutralino mixing. Assuming the hierarchy (naturally fulfilled if ϵ ≪ mW), one can diagonalize the mass matrix in Equation (106) in the seesaw approximation, , obtaining

where Λi = μvi + vdϵi are the so-called alignment parameters. Here M1,2 are the usual gaugino soft mass terms, μ is the Higgsino superpotential mass term, is the VEV and are the sneutrino VEVs (induced by ϵi ≠ 0). The special (projective) form of mν implies that it is a rank 1 matrix, with only one non-zero eigenvalue, identified with the atmospheric mass scale. Furthermore, one can obtain two leptonic mixing angles in terms of the alignment parameters,

The generation of the solar mass scale, which is much smaller (), requires one to go beyond the tree-level approximation. This makes b-p a hybrid radiative neutrino mass model, since loop corrections are necessary in order to reconcile the model with the observations in neutrino oscillation experiments. An example of such loops is shown in Figure 23, where the bottom–sbottom diagrams are displayed. These are found to be the dominant contributions to the solar mass scale generation in most parts of the parameter space of the model. Other relevant contributions are given by the tau-stau and neutrino-sneutrino loops [649651]. In all cases two p projections are required, hence leading to the generation of ΔL = 2 Majorana masses for the light neutrinos.

Figure 23

Figure 23

Bottom–Sbottom diagrams for solar neutrino mass in the b-p model. Open circles correspond to small R-parity violating projections, full circles correspond to R-parity conserving projections and crosses indicate genuine mass insertions which flip chirality. hbYb is the bottom quark Yukawa coupling. Figure reproduced from Diaz et al. [650].

The most important consequence of the breaking of R-parity at the LHC is that the LSP is no longer stable and decays. In fact, this is the only relevant change with respect to the standard MSSM phenomenology. Since the p couplings are constrained to be small, they do not affect the production cross-sections or the intermediate steps of the decay chains, and hence only the LSP decay is altered in an observable way. For instance, the smallness of the p couplings typically imply observable displaced vertices at the LHC, see for instance de Campos et al. [652]. Furthermore, in b-p there is a sharp correlation between the LSP decay and the mixing angles measured in neutrino oscillation experiments [653656]. This connection allows to test the model at colliders. For instance, for a neutralino LSP one finds

A departure from this value would rule out the model completely. Interestingly, these correlations are also found in extended models which effectively lead to bilinear p [657659].

Neutrino masses with t-p

Supersymmetry with trilinear p has many similarities with leptoquark models. Once the trilinear RPV interactions are allowed in the superpotential, the sfermions become scalar fields with lepton and/or baryon number violating interactions, defining properties of a leptoquark. For instance, the right sbottom has the same quantum numbers as the leptoquark S1 discussed in section 5.4.2 and the λ′ coupling in Equation (104) originates a Yukawa interaction exactly like λLQ in Equation (101)38. For this reason, neutrino mass generation takes place in analogous ways, t-p being a pure radiative model.

As already discussed, the breaking of R-parity leads to the decay of the LSP. This is the most distinctive signature of this family of models. However, in contrast to b-p, the large number of free parameters in t-p exclude the possibility of making definite predictions for the LSP decay. Nevertheless, one expects novel signatures at the LHC, typically with many leptons in the final states [661]. Other signatures, already mentioned in section 4.2, include LFV observables, see for instance de Gouvea et al. [486].

6. Conclusions and outlook

The discovery of neutrino oscillations and its explanation in terms of massive neutrinos has been one of the most exciting discoveries in particle physics in recent years and a clear sign of lepton flavor violation and physics beyond the SM. Neutrino masses being the first discovery of physics beyond the SM may be related to the fact that the lowest-order effective operator, the Weinberg operator, generates Majorana neutrino masses. This may point to Majorana neutrinos and consequently lepton number violation introducing a new scale beyond the SM. The magnitude of this scale, and that of lepton flavor violation, are unknown.

The sensitivity to many lepton flavor violating processes will be increased by 2-4 orders of magnitude in the next decade and thus test lepton flavor violation at scales of TeV. In particular the expected improvement of up to 4 orders of magnitude for μ − e conversion and the decay μ → eee, but also other processes, will yield strong constraints on the parameter space of currently allowed models or even more excitingly lead to a discovery. Moreover, the LHC is directly probing the TeV-scale and several possible options for colliders are discussed to probe even higher scales. These exciting experimental prospects, together with the simplicity of the explanation for the smallness of neutrino mass, are the main motivations to study radiative neutrino mass models.

Radiative neutrino mass models explain the lightness of neutrinos without introducing heavy scales. The main idea is that neutrino masses are absent at tree-level, being generated radiatively at 1- or higher-loop orders. This, together with the suppressions due to the possible presence of SM masses and/or extra Yukawa and quartic couplings, implies that the scale of these models may be in the range of TeV. This is also theoretically desirable, because all new particles are light and no hierarchy problem is introduced.

The plethora of neutrino mass models studied in the last decades is overwhelming, reaching the hundreds. We believe that at this point an ordering principle for the theory space is necessary to (i) help scientists outside the field to acquire an overview of the topic, (ii) cover the theory space and spot possible holes, (iii) try to draw generic phenomenological conclusions that can be looked for experimentally, and last but not least (iv) serve as reference for model-builders and phenomenologists.

One can choose to systematically classify the different possibilities and models in different complementary ways: in terms of (i) the effective operators they generate after integrating out the heavy particles at tree-level, (ii) the number of loops at which the Weinberg operator is generated, and (iii) the possible topologies within a particular loop order39. In the first case, the contribution of the matching to the Weinberg operator can be easily estimated, and possible UV completions can be outlined. The second option also sheds light on the scale of the new particles. Finally, the study of possible topologies, which have been analyzed up to 2-loop order, helps to systematically pin down neutrino mass models.

We presented selected examples of radiative neutrino mass models in section 5 which serve as benchmark models and discussed their main phenomenological implications such as lepton flavor-violating processes and direct production of the heavy particles at colliders. The phenomenology is generally very rich and quite model-dependent including extra contributions to neutrinoless double beta decay, electric dipole moments, anomalous magnetic moments, and meson decays. Furthermore, radiative neutrino mass models may solve the dark matter problem with a weakly-interacting massive particle running in the loop generating neutrino mass. Also, the new states can play a crucial role for the matter-antimatter asymmetry, although not necessarily in a positive way, and therefore extra bounds can be set on the lepton number violating interactions.

From our work, we have found that there are several interesting avenues that can be pursued in the future:

  • If anomalies in B-physics [639645], or in the muon anomalous magnetic moment [662], persist, their connection to radiative models should be further pursued.

  • There are only a few studies of the matter anti-matter asymmetry in radiative neutrino mass models and more detailed studies are required.

  • A systematic classification of models generated from effective operators with covariant derivatives40 would help to pin down the possible models involving gauge bosons.

  • Further studies of the symmetries that allow the generation of Dirac masses at loop level.

  • Beyond the LHC, radiative neutrino mass models can be further tested specially if a future collider has initial leptonic states. If those are same sign, one could directly test the neutrino mass matrix by producing for instance the doubly-charged scalar of the Zee-Babu model [622].

To conclude, it is interesting that there are many combinations of what one may call “aesthetically reasonable” particles—those that have SM multiplet assignments and hypercharges that are not too high—that couple to SM particles in such a way as to realize neutrino mass generation at loop level. Radiative mass generation, as well as being a reasonable hypothesis for explaining the smallness of neutrino masses, also provides many phenomenological signatures at relatively low new-physics scales. So, even if nature realizes the seesaw mechanism with heavy right-handed neutrinos, given the difficulty of testing such a paradigm, falsifying radiative models by means of studying in detail their phenomenology and actively searching for their signals seems the only way to strengthen the case of the former by reducing as much as possible the theory space. Not to mention all the useful insights learned on such a journey.

Statements

Author contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Acknowledgments

We acknowledge the use of the TikZ-Feynman package [663]. JHG acknowledges discussions with Arcadi Santamaria and Nuria Rius. This work was supported in part by the Australian Research Council through the ARC Centre of Excellence for Particle Physics at the Terascale (CoEPP) (CE110001104). A.V. acknowledges financial support from the “Juan de la Cierva” program (27-13-463B-731) funded by the Spanish MINECO as well as from the Spanish grants FPA2014-58183-P, Multidark CSD2009-00064, SEV-2014-0398 and PROMETEOII/ 2014/084 (Generalitat Valenciana).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Footnotes

1.^ The possibility of additional neutrino-like states will be discussed below.

2.^ The effective mass mββ depends on the Majorana phases and thus provides a unique probe for them.

3.^ Sterile neutrinos are not charged under the SM gauge group.

4.^ MiniBooNE also has a mysterious excess in their low-energy bins that cannot be explained by any oscillation hypothesis.

5.^ Pseudo-Dirac neutrinos are a special case of Majorana neutrinos where the masses of two Majorana neutrinos are almost degenerate and the breaking of lepton number is small. However, they should not be confused with Dirac neutrinos.

6.^ RH neutrinos could obviously be added as a singlet of SU(5).

7.^ If one does not impose left–right discrete symmetry on the Lagrangian, then there will be no cosmological domain wall problem. The Dirac seesaw mechanism does not require this discrete symmetry.

8.^ This construction resembles the well-known Froggatt-Nielsen mechanism [59].

9.^ The other ΔL = 2 operators also play an important role in the classification of radiative neutrino mass models and will be discussed in detail in section 2.2.2.

10.^ Note that the LL ~ (1, 1, −1) option is irrelevant for tree level mechanisms because it does not produce the required bilinear.

11.^ The bound on the scale of new physics is generally higher for operators involving heavier quarks.

12.^ In a genuine n-loop neutrino mass model, only diagrams starting from n-loop order contribute to neutrino mass. There are no tree level or lower order loop contributions.

13.^ Note that diagrams with scalar or vector bosons are equivalent from a topological point of view.

14.^ The finite 1-loop corrections to the active neutrino mass matrix in the seesaw model were first discussed in Grimus and Neufeld [282] with an arbitrary number of right-handed neutrinos, left-handed lepton doublets, and Higgs doublets. The finite 1-loop corrections are particularly important in case of delicate cancellations in the tree-level neutrino mass terms, which have been studied in Aristizabal Sierra and Yaguna, [283] using the result of Grimus and Lavoura [284].

15.^ In general, 0νββ is generated in these models at a lower order than neutrino masses.

16.^ There are other operators which, however, are simultaneously generated with the Weinberg operator, which dominates as it is dimension 5 [107].

17.^ Higher order effects break universality in a tiny amount due to Higgs interactions, i.e, by the charged lepton Yukawa couplings.

18.^ In models with an extra Higgs doublet coupled to the leptons, other operators can be formed by integrating out the second Higgs doublet. In those cases, the electrons involved are right-handed and therefore there is no interference with the W boson. An example is the Zee model, see Herrero-García et al. [91].

19.^ See Lee et al. [446], Lee and Shrock [447] and Marciano and Sanda [448] for pioneering work on LFV processes.

20.^ Note the different choice of Lorentz structures in Equations (48), (49). The two forms can be related via the Gordon-identity.

21.^ Carpentier and Davidson [449] provides a comprehensive collection of constraints on quark flavor violating operators.

22.^ We clarify that this is only true for the processes and (with β ≠ δ). The process (with β ≠ δ) does not receive contributions from penguin diagrams, but only from boxes.

23.^ In some models, cancellations due to certain Yukawa textures can affect some decays (like μ → ), but it is virtually impossible to cancel all radiative decays simultaneously.

24.^ Interestingly, the authors of Chowdhury and Nasri [129] have shown that in variants of the scotogenic model with higher SU(2) representations the LFV rates become larger due to additive effects from the components of the large multiplets.

25.^ Also in type-I seesaw (and inverse seesaw), and in the MSSM, HLFV is generated at 1-loop order [494498].

26.^ Theoretically there is no SM background. Realistically, however, object misidentification, undetected particles and fake objects can result in similar final states at the detector level.

27.^ Some of the exotic particles may also show up in tree-level neutrino mass models. The interested reader is referred to the recent review [531] for the collider tests of specific tree-level models.

28.^ Strictly speaking the process is not necessarily LNV, because X may carry lepton number as well, for example in form of neutrinos. Currently the searches are limited to electrons and muons. However, τ-leptons may also be used to search for LNV.

29.^ Long-lived charged particles have been searched at the LHC using anomalously high ionization signal [572], also in the context of dark matter [573]. However, charged scalars in radiative neutrino mass models usually have sizable couplings to SM leptons and decay promptly.

30.^ In this last case, of course, the presence of low scale LNV can be regarded as being less motivated, as in principle there would already be an explanation for neutrino masses (at least for one neutrino).

31.^ Less copies of N means less contribution to the neutrino mass matrix, which in turn generally leads to larger Yukawa couplings to generate the same neutrino mass scale and thus more likely to violate constraints from LFV processes.

32.^ The authors of Chen et al. [356] also point out that up-quarks are not feasible due to gauge invariance.

33.^ The scotogenic model has been extensively studied, sometimes referring to it with different names. For instance, some authors prefer the denomination radiative seesaw. In this review we will stick to the more popular name scotogenic model, which comes from the Greek word skotosoτoσ), darkness. scotogenic would then mean created from darkness.

34.^ The Z2 symmetry can obtained from the spontaneous breaking of an Abelian U(1) factor, see for instance Aristizabal Sierra et al. [623].

35.^ We follow the nomenclature in Doršner et al. [638] and Buchmuller et al. [542] for the names of the leptoquarks, where subscripts indicate dimension of the SU(2)L representations.

36.^ Aristizabal Sierra et al. [181] considered the most general interactions with all possible leptoquarks and found in total four mass matrices for leptoquarks with electric charges , , and .

37.^ The contributions of the strange and down quarks are suppressed by and thus have been neglected in the discussion of Angel et al. [287].

38.^ There are, however, additional couplings that supersymmetry forbids but would be allowed for general leptoquarks. Therefore, t-p can then be regarded as a constrained leptoquark scenario. See Deshpande and He [660] for a paper on t-p as a possible explanation for the B-meson anomalies that highlights the similarities between this setup and leptoquark models.

39.^ A fourth complementary classification in terms of particles can be done, which will appear in a future publication Cai, (in prep).

40.^ All possible dimension-7 operators with SM fields and right-handed neutrinos have been listed in Bhattacharya and Wudka [86].

References

  • 1.

    Hirata KS Inoue K Kajita T Kifune T Kihara K Nakahata M et al . Results from one thousand days of real time directional solar neutrino data. Phys Rev Lett. (1990) 65:1297300. 10.1103/PhysRevLett.65.1297

  • 2.

    Hirata KS Inoue K Ishida T Kajita T Kihara K Nakahata M et al . Observation of a small atmospheric νμe ratio in Kamiokande. Phys Lett. (1992) B280:14652. 10.1016/0370-2693(92)90788-6

  • 3.

    Fukuda Y Hayakawa T Ichihara E Inoue K Ishihara K Ishino H et al . Evidence for oscillation of atmospheric neutrinos. Phys Rev Lett. (1998) 81:15627. 10.1103/PhysRevLett.81.1562

  • 4.

    Cleveland BT Daily T Davis R Jr Distel JR Lande K Lee CK et al . Measurement of the solar electron neutrino flux with the Homestake chlorine detector. Astrophys J. (1998) 496:50526. 10.1086/305343

  • 5.

    Hampela W Handta J Heussera G Kikoa J Kirstena T Laubensteina M et al . GALLEX solar neutrino observations: Results for GALLEX IV. Phys Lett. (1999) B447:12733.

  • 6.

    Abdurashitov JN Gavrin VN Gorbachev VV Gurkina PP Ibragimova TV Kalikhov AV et al . Measurement of the solar neutrino capture rate with gallium metal. Phys Rev. (1999) C60:055801.

  • 7.

    Fukuda S Fukuda Y Ishitsuka M Itow Y Kajita T Kameda J et al . Solar B-8 and hep neutrino measurements from 1258 days of Super-Kamiokande data. Phys Rev Lett. (2001) 86:56515. 10.1103/PhysRevLett.86.5651

  • 8.

    Ahmad QR Allen RC Andersen TC D Anglin J Barton JC Beier EW et al . Direct evidence for neutrino flavor transformation from neutral current interactions in the Sudbury Neutrino Observatory. Phys Rev Lett. (2002) 89:011301. 10.1103/PhysRevLett.89.011301

  • 9.

    Altmann M Balata M Belli P Bellotti E Bernabei R Burkert E et al . Complete results for five years of GNO solar neutrino observations. Phys Lett. (2005) B616:17490. 10.1016/j.physletb.2005.04.068

  • 10.

    Ahn MH Aliu E Andringa S Aoki S Aoyama Y Argyriades J et al . Measurement of neutrino oscillation by the K2K experiment. Phys Rev. (2006) D74:072003. 10.1103/PhysRevD.74.072003

  • 11.

    Michael DG Adamson P Alexopoulos T Allison WW Alner GJ Anderson K et al . Observation of muon neutrino disappearance with the MINOS detectors and the NuMI neutrino beam. Phys Rev Lett. (2006) 97:191801. 10.1103/PhysRevLett.97.191801

  • 12.

    Abe S Ebihara T Enomoto S Furuno K Gando Y Ichimura K et al . Precision measurement of neutrino oscillation parameters with KamLAND. Phys Rev Lett. (2008) 100:221803. 10.1103/PhysRevLett.100.221803

  • 13.

    Abe K Abgrall N Ajima Y Aihara H Albert JB Andreopoulos C et al . Indication of electron neutrino appearance from an accelerator-produced off-axis muon neutrino beam. Phys Rev Lett. (2011) 107:041801. 10.1103/PhysRevLett.107.041801

  • 14.

    Abe Y Aberle C Akiri T dos Anjos JC Ardellier F Barbosa AF et al . Indication of reactor disappearance in the double chooz experiment. Phys Rev Lett. (2012) 108:131801. 10.1103/PhysRevLett.108.131801

  • 15.

    An FP Bai JZ Balantekin AB Band HR Beavis D Beriguete W et al . Observation of electron-antineutrino disappearance at Daya Bay. Phys Rev Lett. (2012) 108:171803. 10.1103/PhysRevLett.108.171803

  • 16.

    Ahn JK Chebotaryov S Choi JH Choi S Choi W Choi Y et al . Observation of reactor electron antineutrino disappearance in the RENO experiment. Phys Rev Lett. (2012) 108:191802. 10.1103/PhysRevLett.108.191802

  • 17.

    Abe K Adam J Aihara H Akiri T Andreopoulos C Aoki S . Precise measurement of the neutrino mixing parameter θ23 from Muon neutrino disappearance in an off-axis beam. Phys Rev Lett. (2014) 112:181801. 10.1103/PhysRevLett.112.181801

  • 18.

    Pontecorvo B . Inverse beta processes and nonconservation of lepton charge. Sov Phys JETP. (1958) 7:1723.

  • 19.

    Maki Z Nakagawa M Sakata S . Remarks on the unified model of elementary particles. Prog Theor Phys. (1962) 28:87080. 10.1143/PTP.28.870

  • 20.

    Esteban I Gonzalez-Garcia MC Maltoni M Martinez-Soler I Schwetz T . Updated fit to three neutrino mixing: exploring the accelerator-reactor complementarity. J High Energy Phys. (2017) 1:087. 10.1007/JHEP01(2017)087

  • 21.

    Forero DV Tortola M Valle JWF . Neutrino oscillations refitted. Phys Rev. (2014) D90:093006. 10.1103/PhysRevD.90.093006

  • 22.

    Capozzi F Fogli GL Lisi E Marrone A Montanino D Palazzo A . Status of three-neutrino oscillation parameters, circa 2013. Phys Rev. (2014) D89:093018. 10.1103/PhysRevD.89.093018

  • 23.

    Mikheev SP Smirnov AYu . Resonant amplification of neutrino oscillations in matter and solar neutrino spectroscopy. Nuovo Cim. (1986) C9:1726. 10.1007/BF02508049

  • 24.

    Wolfenstein L . Neutrino oscillations in matter. Phys Rev. (1978) D17:236974.

  • 25.

    Fukugita M Yanagida T . Baryogenesis without grand unification. Phys Lett. (1986) B174:457.

  • 26.

    Lobashev VM . The search for the neutrino mass by direct method in the tritium beta-decay and perspectives of study it in the project KATRIN. Nucl Phys. (2003) A719:15360. 10.1016/S0375-9474(03)00985-0

  • 27.

    Kraus C Bornschein B Bornschein L Bonn J Flatt B Kovalik A et al . Final results from phase II of the Mainz neutrino mass search in tritium beta decay. Eur Phys J. (2005) C40:44768. 10.1140/epjc/s2005-02139-7

  • 28.

    Aseev VN Belesev AI Berlev AI Geraskin EV Golubev AA Likhovid NA et al . An upper limit on electron antineutrino mass from Troitsk experiment. Phys Rev. (2011) D84:112003. 10.1103/PhysRevD.84.112003

  • 29.

    Osipowicz A Blümer H Drexlin G Eitel K Meisel G Plischke P et al . KATRIN: a next generation tritium beta decay experiment with sub-eV sensitivity for the electron neutrino mass. Lett. Intent (2001). Available online at: https://arxiv.org/pdf/hep-ex/0109033.pdf

  • 30.

    Ade PAR Aghanim N Arnaud M Ashdown M Aumont J Baccigalupi C et al . Planck 2015 results. XIII. Cosmological parameters. Astron Astrophys. (2016) 594:A13. 10.1051/0004-6361/201525830

  • 31.

    Agostini M Allardt M Bakalyarov AM Balata M Barabanov I Baudis L et al . Background free search for neutrinoless double beta decay with GERDA Phase II. (2017) Nature544:47. 10.1038/nature21717

  • 32.

    Gariazzo S Giunti C Laveder M Li YF . Updated global 3+1 analysis of short-baseline neutrino oscillations. J High Energy Phys. (2017) 6:135. 10.1007/JHEP06(2017)135

  • 33.

    Kopp J Machado PAN Maltoni M Schwetz T . Sterile neutrino oscillations: the global picture. J High Energy Phys. (2013) 5:050. 10.1007/JHEP05(2013)050

  • 34.

    Athanassopoulos C Auerbach LB Burman RL Caldwell DO Church ED Cohen I et al . Evidence for νμ → νe neutrino oscillations from LSND. Phys Rev Lett. (1998) 81:17747. 10.1103/PhysRevLett.81.1774

  • 35.

    Aguilar-Arevalo A Auerbach LB Burman RL Caldwell DO Church ED Cochran AK et al . Evidence for neutrino oscillations from the observation of anti-neutrino(electron) appearance in a anti-neutrino(muon) beam. Phys Rev. (2001) D64:112007. 10.1103/PhysRevD.64.112007

  • 36.

    Aguilar-Arevalo AA Bazarko AO Brice SJ Brown BC Bugel L Cao J et al . A Search for electron neutrino appearance at the Δm2 ~ 1eV2 scale. Phys Rev Lett. (2007) 98:231801. 10.1103/PhysRevLett.98.231801

  • 37.

    Aguilar-Arevalo AA Brown BC Bugel L Cheng G Church ED Conrad JM et al . Improved search for oscillations in the miniBooNE experiment. Phys Rev Lett. (2013) 110:161801. 10.1103/PhysRevLett.110.161801

  • 38.

    Aartsen MG Ackermann M Adams J Aguilar JA Ahlers M Ahrens I et al . Search for sterile neutrino mixing using three years of IceCube DeepCore data. Phys Rev. (2017) D95:112002. 10.1103/PhysRevD.95.112002

  • 39.

    Mueller TA Lhuillier D Fallot M Letourneau A Cormon S Fechner M et al . Improved predictions of reactor antineutrino spectra. Phys Rev. (2011) C83:054615. 10.1103/PhysRevC.83.054615

  • 40.

    Huber P . On the determination of anti-neutrino spectra from nuclear reactors. Phys Rev. (2011) C84:024617. 10.1103/PhysRevC.84.024617

  • 41.

    Mention G Fechner M Lasserre T Mueller TA Lhuillier D Cribier M et al . The reactor antineutrino anomaly. Phys Rev. (2011) D83:073006. 10.1103/PhysRevD.83.073006

  • 42.

    Hampel W Heusser G Kiko J Kirsten TA Laubenstein M Pernicka E et al . Final results of the Cr-51 neutrino source experiments in GALLEX. Phys Lett. (1998) B420:11426.

  • 43.

    Abdurashitov JN Gavrin VN Girin SV Gorbachev VV Ibragimova TV Kalikhov AV et al . Measurement of the response of the Russian-American gallium experiment to neutrinos from a Cr-51 source. Phys Rev. (1999) C59:224663.

  • 44.

    Abdurashitov JN Gavrin VN Gorbachev VV Gurkina PP Ibragimova TV Kalikhov AV et al . Measurement of the solar neutrino capture rate with gallium metal. III: results for the 2002–2007 data-taking period. Phys Rev. (2009) C80:015807. 10.1103/PhysRevC.80.015807

  • 45.

    Kaether F Hampel W Heusser G Kiko J Kirsten T . Reanalysis of the GALLEX solar neutrino flux and source experiments. Phys Lett. (2010) B685:4754. 10.1016/j.physletb.2010.01.030

  • 46.

    An FP Balantekin AB Band HR Bishai M Blyth S Cao D et al . Evolution of the reactor antineutrino flux and spectrum at Daya bay. Phys Rev Lett. (2017) 118:251801. 10.1103/PhysRevLett.118.251801

  • 47.

    Minkowski P . μ → eγ at a Rate of One Out of 109 Muon Decays?Phys Lett. (1977) B67:4218. 10.1016/0370-2693(77)90435-X

  • 48.

    Yanagida T . Horizontal symmetry and masses of neutrinos. Conf Proc. (1979) C7902131:959.

  • 49.

    Gell-Mann M Ramond P Slansky R . Complex spinors and unified theories. Conf Proc. (1979) C790927:31521.

  • 50.

    Mohapatra RN Senjanovic G . Neutrino mass and spontaneous parity violation. Phys Rev Lett. (1980) 44:912. 10.1103/PhysRevLett.44.912

  • 51.

    Glashow SL . The future of elementary particle physics. NATO Sci Ser B (1980) 61:687. 10.1007/978-1-4684-7197-7_15

  • 52.

    Foot R Joshi GC Lew H Volkas RR . Charge quantization in the standard model and some of its extensions. Mod Phys Lett. (1990) A5:272132. 10.1142/S0217732390003176

  • 53.

    Babu KS Mohapatra RN . Is there a connection between quantization of electric charge and a majorana neutrino?Phys Rev Lett. (1989) 63:938. 10.1103/PhysRevLett.63.938

  • 54.

    't Hooft G . Naturalness, chiral symmetry, and spontaneous chiral symmetry breaking. NATO Sci Ser B (1980) 59:13557.

  • 55.

    Roy P Shanker OU . Observable neutrino dirac mass and supergrand unification. Phys Rev Lett. (1984) 52:7136. 10.1103/PhysRevLett.52.713

  • 56.

    Ma E Popov O . Pathways to naturally small dirac neutrino masses. Phys Lett. (2017) B764:1424. 10.1016/j.physletb.2016.11.027

  • 57.

    Wang F Wang W Yang JM . Split two-higgs-doublet model and neutrino condensation. EuroPhys Lett. (2006) 76:38894. 10.1209/epl/i2006-10293-3

  • 58.

    Gabriel S Nandi S . A new two Higgs doublet model. Phys Lett. (2007) B655:1417. 10.1016/j.physletb.2007.04.062

  • 59.

    Froggatt CD Nielsen HB . Hierarchy of quark masses, cabibbo angles and CP violation. Nucl Phys. (1979) B147:27798.

  • 60.

    Wang W Han ZL . Naturally small dirac neutrino mass with intermediate SU(2)L multiplet fields. J High Energy Phys. (2016) 4:166. 10.1007/JHEP04(2017)166

  • 61.

    Randall L Sundrum R . A Large mass hierarchy from a small extra dimension. Phys Rev Lett. (1999) 83:33703. 10.1103/PhysRevLett.83.3370

  • 62.

    Randall L Sundrum R . An alternative to compactification. Phys Rev Lett. (1999) 83:46903. 10.1103/PhysRevLett.83.4690

  • 63.

    Grossman Y Neubert M . Neutrino masses and mixings in nonfactorizable geometry. Phys Lett. (2000) B474:36171. 10.1016/S0370-2693(00)00054-X

  • 64.

    Gherghetta T Pomarol A . Bulk fields and supersymmetry in a slice of AdS. Nucl Phys. (2000) B586:14162. 10.1016/S0550-3213(00)00392-8

  • 65.

    Huber SJ Shafi Q . Neutrino oscillations and rare processes in models with a small extra dimension. Phys Lett. (2001) B512:36572. 10.1016/S0370-2693(01)00736-5

  • 66.

    Moreau G Silva-Marcos JI . Neutrinos in warped extra dimensions. J High Energy Phys. (2006) 1:048. 10.1088/1126-6708/2006/01/048

  • 67.

    Agashe K Okui T Sundrum R . A common origin for neutrino anarchy and charged hierarchies. Phys Rev Lett. (2009) 102:101801. 10.1103/PhysRevLett.102.101801

  • 68.

    Chang WF Ng JN Wu JMS . Discrete dirac neutrino in warped extra dimensions. Phys Rev. (2009) D80:113013. 10.1103/PhysRevD.80.113013

  • 69.

    De Gouvea A Giudice GF Strumia A Tobe K . Phenomenological implications of neutrinos in extra dimensions. Nucl Phys. (2002) B623:395420. 10.1016/S0550-3213(01)00621-6

  • 70.

    Choi K Im SH . Realizing the relaxion from multiple axions and its UV completion with high scale supersymmetry. J High Energy Phys. (2016) 1:149. 10.1007/JHEP01(2016)149

  • 71.

    Kaplan DE Rattazzi R . Large field excursions and approximate discrete symmetries from a clockwork axion. Phys Rev. (2016) D93:085007. 10.1103/PhysRevD.93.085007

  • 72.

    Hambye T Teresi D Tytgat MHG . A clockwork WIMP. J High Energy Phys. (2017) 7:47. 10.1007/JHEP07(2017)047

  • 73.

    Giudice GF McCullough M . A clockwork theory. J High Energy Phys. (2017) 2:036. 10.1007/JHEP02(2017)036

  • 74.

    Magg M Wetterich C . Neutrino mass problem and gauge hierarchy. Phys Lett. (1980) B94:614. 10.1016/0370-2693(80)90825-4

  • 75.

    Schechter J Valle JWF . Neutrino masses in SU(2) × U(1) theories. Phys Rev. (1980) D22:2227.

  • 76.

    Cheng TP Li LF . Neutrino masses, mixings and oscillations in SU(2) × U(1) models of electroweak interactions. Phys Rev. (1980) D22:2860.

  • 77.

    Lazarides G Shafi Q Wetterich C . Proton lifetime and Fermion masses in an SO(10) model. Nucl Phys. (1981) B181:287300.

  • 78.

    Wetterich C . Neutrino masses and the scale of B-L violation. Nucl Phys. (1981) B187:34375.

  • 79.

    Mohapatra RN Senjanovic G . Neutrino masses and mixings in Gauge models with spontaneous parity violation. Phys Rev. (1981) D23:165.

  • 80.

    Foot R Lew H He XG Joshi GC . Seesaw neutrino masses induced by a triplet of leptons. Z Phys. (1989) C44:441.

  • 81.

    Weinberg S . Baryon and Lepton nonconserving processes. Phys Rev Lett. (1979) 43:156670. 10.1103/PhysRevLett.43.1566

  • 82.

    Kumericki K Picek I Radovcic B . TeV-scale seesaw with quintuplet fermions. Phys Rev. (2012) D86:013006. 10.1103/PhysRevD.86.013006

  • 83.

    Picek I Radovcic B . Enhancement of h → γγ by seesaw-motivated exotic scalars. Phys Lett. (2013) B719:4048. 10.1016/j.physletb.2013.01.056

  • 84.

    Liao Y . Cascade seesaw for tiny neutrino mass. J High Energy Phys. (2011) 6:98. 10.1007/JHEP06(2011)098

  • 85.

    Babu KS Leung CN . Classification of effective neutrino mass operators. Nucl Phys. (2001) B619:66789. 10.1016/S0550-3213(01)00504-1

  • 86.

    Bhattacharya S Wudka J . Dimension-seven operators in the standard model with right handed neutrinos. Phys Rev. (2016) D94:055022. 10.1103/PhysRevD.94.055022

  • 87.

    de Gouvea A Jenkins J . A survey of lepton number violation via effective operators. Phys Rev. (2008) D77:013008. 10.1103/PhysRevD.77.013008

  • 88.

    de Gouvêa A Herrero-García J Kobach A . Neutrino masses, grand unification, and baryon number violation. Phys Rev. (2014) D90:016011. 10.1103/PhysRevD.90.016011

  • 89.

    Kobach A . Baryon number, lepton number, and operator dimension in the standard model. Phys Lett. (2016) B758:4557. 10.1016/j.physletb.2016.05.050

  • 90.

    Angel PW Rodd NL Volkas RR . Origin of neutrino masses at the LHC: ΔL = 2 effective operators and their ultraviolet completions. Phys Rev. (2013) D87:073007. 10.1103/PhysRevD.87.073007

  • 91.

    Herrero-García J Ohlsson T Riad S Wirén J . Full parameter scan of the Zee model: exploring Higgs lepton flavor violation. J High Energy Phys. (2017) 4:130. 10.1007/JHEP04(2017)130

  • 92.

    Frere JM Jones DRT Raby S . Fermion masses and induction of the weak scale by supergravity. Nucl Phys. (1983) B222:119.

  • 93.

    Alvarez-Gaume L Polchinski J Wise MB . Minimal low-energy supergravity. Nucl Phys. (1983) B221:495.

  • 94.

    Casas JA Dimopoulos S . Stability bounds on flavor violating trilinear soft terms in the MSSM. Phys Lett. (1996) B387:10712.

  • 95.

    Herrero-Garcia J Nebot M Rius N Santamaria A . The zee–babu model revisited in the light of new data. Nucl Phys. (2014) B885:54270. 10.1016/j.nuclphysb.2014.06.001

  • 96.

    Nomura T Okada H . A four-loop Radiative Seesaw Model. Phys Lett. (2017) B770:30713. 10.1016/j.physletb.2017.05.002

  • 97.

    Nomura T Okada H . Four-loop neutrino model inspired by diphoton excess at 750 GeV. Phys Lett. (2016) B755:30611. 10.1016/j.physletb.2016.02.022

  • 98.

    Cárcamo Hernández AE Kovalenko S Schmidt I . Radiatively generated hierarchy of lepton and quark masses. J High Energy Phys. (2017) 2:125. 10.1007/JHEP02(2017)125

  • 99.

    Ma E . Pathways to naturally small neutrino masses. Phys Rev Lett. (1998) 81:11714. 10.1103/PhysRevLett.81.1171

  • 100.

    Bonnet F Hirsch M Ota T Winter W . Systematic study of the d=5 Weinberg operator at one-loop order. J High Energy Phys. (2012) 7:153. 10.1007/JHEP07(2012)153

  • 101.

    Cepedello R Hirsch M Helo JC . Loop neutrino masses from d = 7 operator. J High Energy Phys. (2017) 7:079. 10.1007/JHEP07(2017)079

  • 102.

    Aristizabal Sierra D Degee A Dorame L Hirsch M . Systematic classification of two-loop realizations of the Weinberg operator. J High Energy Phys. (2015) 3:40. 10.1007/JHEP03(2015)040

  • 103.

    Law SSC McDonald KL . The simplest models of radiative neutrino mass. Int J Mod Phys. (2014) A29:1450064. 10.1142/S0217751X1450064X

  • 104.

    Zee A . A theory of lepton number violation, neutrino majorana mass, and oscillation. Phys Lett. (1980) B93:389. 10.1016/0370-2693(80)90349-4

  • 105.

    Wolfenstein L . A theoretical pattern for neutrino oscillations. Nucl Phys. (1980) B175:936. 10.1016/0550-3213(80)90004-8

  • 106.

    Bonnet F Hirsch M Ota T Winter W . Systematic decomposition of the neutrinoless double beta decay operator. J High Energy Phys. (2013) 3:55. 10.1007/JHEP03(2013)055

  • 107.

    del Aguila F Aparici A Bhattacharya S Santamaria A Wudka J . Effective Lagrangian approach to neutrinoless double beta decay and neutrino masses. J High Energy Phys. (2012) 6:146. 10.1007/JHEP06(2012)146

  • 108.

    Helo JC Hirsch M Ota T . Long-range contributions to double beta decay revisited. J High Energy Phys. (2016) 6:6. 10.1007/JHEP06(2016)006

  • 109.

    Hagedorn C Ohlsson T Riad S Schmidt MA . Unification of gauge couplings in radiative neutrino mass models. J High Energy Phys. (2016) 9:111. 10.1007/JHEP09(2016)111

  • 110.

    Farzan Y Pascoli S Schmidt MA . Recipes and ingredients for neutrino mass at loop level. J High Energy Phys. (2013) 3:107. 10.1007/JHEP03(2013)107

  • 111.

    Restrepo D Zapata O Yaguna CE . Models with radiative neutrino masses and viable dark matter candidates. J High Energy Phys. (2013) 11:11. 10.1007/JHEP11(2013)011

  • 112.

    Simoes C Wegman D . Radiative two-loop neutrino masses with dark matter. J High Energy Phys. (2017) 4:148. 10.1007/JHEP04(2017)148

  • 113.

    Ma E . Verifiable radiative seesaw mechanism of neutrino mass and dark matter. Phys Rev. (2006) D73:077301. 10.1103/PhysRevD.73.077301

  • 114.

    Ma E Suematsu D . Fermion triplet dark matter and radiative neutrino mass. Mod Phys Lett. (2009) A24:5839. 10.1142/S021773230903059X

  • 115.

    Chao W . Dark matter, LFV and neutrino magnetic moment in the radiative seesaw model with fermion triplet. Int J Mod Phys. (2015) A30:1550007. 10.1142/S0217751X15500074

  • 116.

    von der Pahlen F Palacio G Restrepo D Zapata O . Radiative type III seesaw model and its collider phenomenology. Phys Rev. (2016) D94:033005. 10.1103/PhysRevD.94.033005

  • 117.

    Farzan Y . A minimal model linking two great mysteries: neutrino mass and dark matter. Phys Rev. (2009) D80:073009. 10.1103/PhysRevD.80.073009

  • 118.

    Hehn D Ibarra A . A radiative model with a naturally mild neutrino mass hierarchy. Phys Lett. (2013) B718:98891. 10.1016/j.physletb.2012.11.034

  • 119.

    Lu W-B Gu P-H . Mixed inert scalar triplet dark matter, radiative neutrino masses and leptogenesis. Nucl Phys. (2017) B924:279311. 10.1016/j.nuclphysb.2017.09.005

  • 120.

    Fileviez Perez P Wise MB . On the origin of neutrino masses. Phys Rev. (2009) D80:053006. 10.1103/PhysRevD.80.053006

  • 121.

    Liao Y Liu JY . Radiative and flavor-violating transitions of leptons from interactions with color-octet particles. Phys Rev. (2010) D81:013004. 10.1103/PhysRevD.81.013004

  • 122.

    Farzan Y Pascoli S Schmidt MA . AMEND: a model explaining neutrino masses and dark matter testable at the LHC and MEG. J High Energy Phys. (2010) 10:111. 10.1007/JHEP10(2010)111

  • 123.

    Okada H Orikasa Y . Radiative neutrino model with an inert triplet scalar. Phys Rev. (2016) D94:055002. 10.1103/PhysRevD.94.055002

  • 124.

    Brdar V Picek I Radovcic B . Radiative neutrino mass with scotogenic scalar triplet. Phys Lett. (2014) B728:198201. 10.1016/j.physletb.2013.11.045

  • 125.

    Aoki M Kanemura S Yagyu K . Doubly-charged scalar bosons from the doublet. Phys Lett. (2011) B702:3558. 10.1016/j.physletb.2011.07.017

  • 126.

    Cai Y He XG Ramsey-Musolf M Tsai LH . RνMDM and lepton flavor violation. J High Energy Phys. (2011) 12:054. 10.1007/JHEP12(2011)054

  • 127.

    Chen CH Law SSC . Exotic fermion multiplets as a solution to baryon asymmetry, dark matter and neutrino masses. Phys Rev. (2012) D85:055012. 10.1103/PhysRevD.85.055012

  • 128.

    Law SSC McDonald KL . A class of inert N-tuplet models with radiative neutrino mass and dark matter. J High Energy Phys. (2013) 9:92. 10.1007/JHEP09(2013)092

  • 129.

    Chowdhury TA Nasri S . Lepton flavor violation in the inert scalar model with higher representations. J High Energy Phys. (2015) 12:40. 10.1007/JHEP12(2015)040

  • 130.

    Chen CH Geng CQ Zhuridov DV . Neutrino masses, leptogenesis and decaying dark matter. J Cosmol Astropart Phys. (2009) 0910:001. 10.1088/1475-7516/2009/10/001

  • 131.

    Patra S Sahoo N Sahu N . Dipolar dark matter in light of the 3.5 keV x-ray line, neutrino mass, and LUX data. Phys Rev. (2015) D91:115013. 10.1103/PhysRevD.91.115013

  • 132.

    Ferreira PM Grimus W Jurciukonis D Lavoura L . Scotogenic model for co-bimaximal mixing. J High Energy Phys. (2016) 7:10. 10.1007/JHEP07(2016)010

  • 133.

    Fortes ECFS Machado ACB Montaño J Pleitez V . Lepton masses and mixing in a scotogenic model. (2017) arXiv:1705.09414.

  • 134.

    Ma E . Dark scalar doublets and neutrino tribimaximal mixing from A4 symmetry. Phys Lett. (2009) B671:3668. 10.1016/j.physletb.2008.12.038

  • 135.

    Bhattacharya S Ma E Natale A Rashed A . Radiative scaling neutrino mass with A4 symmetry. Phys Rev. (2013) D87:097301. 10.1103/PhysRevD.87.097301

  • 136.

    Carcamo Hernandez AE de Medeiros Varzielas I Kovalenko SG Päs H Schmidt I . Lepton masses and mixings in an A4 multi-Higgs model with a radiative seesaw mechanism. Phys Rev. (2013) D88:076014. 10.1103/PhysRevD.88.076014

  • 137.

    Campos MD Cárcamo Hernández AE Kovalenko S Schmidt I Schumacher E . Fermion masses and mixings in an SU(5) grand unified model with an extra flavor symmetry. Phys Rev. (2014) D90:016006. 10.1103/PhysRevD.90.016006

  • 138.

    Ma E . Neutrino mixing and geometric CP violation with delta(27) symmetry. Phys Lett. (2013) B723:1613. 10.1016/j.physletb.2013.05.011

  • 139.

    Ma E Natale A . Scotogenic Z2 or U(1)D model of neutrino mass with Δ(27) symmetry. Phys Lett. (2014) B734:4035. 10.1016/j.physletb.2014.05.070

  • 140.

    Gu PH Sarkar U . Radiative seesaw in left-right symmetric model. Phys Rev. (2008) D78:073012. 10.1103/PhysRevD.78.073012

  • 141.

    Adulpravitchai A Lindner M Merle A Mohapatra RN . Radiative transmission of lepton flavor hierarchies. Phys Lett. (2009) B680:4769. 10.1016/j.physletb.2009.09.042

  • 142.

    Parida MK . Radiative seesaw in SO(10) with dark matter. Phys Lett. (2011) B704:20610. 10.1016/j.physletb.2011.09.016

  • 143.

    Ma E . Unified framework for matter, dark matter, and radiative neutrino mass. Phys Rev. (2013) D88:117702. 10.1103/PhysRevD.88.117702

  • 144.

    Megrelidze L Tavartkiladze Z . Soft see-saw: radiative origin of neutrino masses in SUSY theories. Nucl Phys. (2017) B914:55376. 10.1016/j.nuclphysb.2016.11.019

  • 145.

    Ma E Picek I Radovčić B . New scotogenic model of neutrino mass with U(1)D gauge interaction. Phys Lett. (2013) B726:7446. 10.1016/j.physletb.2013.09.049

  • 146.

    Ma E . Supersymmetric model of radiative seesaw majorana neutrino masses. Ann Fond Broglie (2006) 31:285. Available online at: https://arxiv.org/abs/hep-ph/0607142

  • 147.

    Ma E . SU(5) completion of the dark scalar doublet model of radiative neutrino mass. Phys Lett. (2008) B659:8857. 10.1016/j.physletb.2007.12.013

  • 148.

    Khalil S Lee HS Ma E . Generalized lepton number and dark left-right gauge model. Phys Rev. (2009) D79:041701. 10.1103/PhysRevD.79.041701

  • 149.

    Ma E . Variants of the dark left-right gauge model: neutrinos and scotinos. Phys Rev. (2009) D79:117701. 10.1103/PhysRevD.79.117701

  • 150.

    Budhi RHS Kashiwase S Suematsu D . Inflation in a modified radiative seesaw model. Phys Rev. (2014) D90:113013. 10.1103/PhysRevD.90.113013

  • 151.

    Kashiwase S Suematsu D . Lepton number asymmetry via inflaton decay in a modified radiative seesaw model. Phys Lett. (2015) B749:60312. 10.1016/j.physletb.2015.08.062

  • 152.

    Budhi RHS Kashiwase S Suematsu D . Inflation due to a nonminimal coupling of singlet scalars in the radiative seesaw model. Phys Rev. (2016) D93:013022. 10.1103/PhysRevD.93.013022

  • 153.

    Arhrib A Bœhm C Ma E Yuan TC . Radiative model of neutrino mass with neutrino interacting MeV dark matter. J Cosmol Astropart Phys. (2016) 1604:49. 10.1088/1475-7516/2016/04/049

  • 154.

    Lu WB Gu PH . Leptogenesis, radiative neutrino masses and inert Higgs triplet dark matter. J Cosmol Astropart Phys. (2016) 1605:040. 10.1088/1475-7516/2016/05/040

  • 155.

    Koide Y . Can the zee model explain the observed neutrino data?Phys Rev. (2001) D64:077301. 10.1103/PhysRevD.64.077301

  • 156.

    He XG . Is the zee model neutrino mass matrix ruled out?Eur Phys J. (2004) C34:3716. 10.1140/epjc/s2004-01669-8

  • 157.

    He XG Majee SK . Implications of recent data on neutrino mixing and lepton flavour violating decays for the zee model. J High Energy Phys. (2012) 3:23. 10.1007/JHEP03(2012)023

  • 158.

    Babu KS Julio J . Predictive model of radiative neutrino masses. Phys Rev. (2014) D89:053004. 10.1103/PhysRevD.89.053004

  • 159.

    Aranda A Bonilla C Rojas AD . Neutrino masses generation in a Z4 model. Phys Rev. (2012) D85:036004. 10.1103/PhysRevD.85.036004

  • 160.

    Leontaris GK Tamvakis K . The supersymmetric singlet. Phys Lett. (1987) B191:4216. 10.1016/0370-2693(87)90633-2

  • 161.

    Haba N Matsuda M Tanimoto M . Large neutrino flavor mixings and gauge mediated supersymmetry breaking scenario. Phys Lett. (2000) B478:3517. 10.1016/S0370-2693(00)00275-6

  • 162.

    Cheung Km Kong OCW . Zee neutrino mass model in SUSY framework. Phys Rev. (2000) D61:113012. 10.1103/PhysRevD.61.113012

  • 163.

    Kanemura S Shindou T Sugiyama H . R-parity conserving supersymmetric extension of the zee model. Phys Rev. (2015) D92:115001. 10.1103/PhysRevD.92.115001

  • 164.

    Tamvakis K Vergados JD . Neutrino masses and oscillations in an unconventional model of lepton number violation. Phys Lett. (1985) B155:37380.

  • 165.

    Fileviez Perez P Murgui C . Renormalizable SU(5) unification. Phys Rev. (2016) D94:075014. 10.1103/PhysRevD.94.075014

  • 166.

    Chang WF Ng JN . Neutrino masses in 5-D orbifold SU(5) unification models without right handed singlets. J High Energy Phys. (2003) 10:036. 10.1088/1126-6708/2003/10/036

  • 167.

    Chang WF Chen IT Liou SC . Neutrino masses via the zee mechanism in 5D split fermions model. Phys Rev. (2011) D83:025017. 10.1103/PhysRevD.83.025017

  • 168.

    Babu KS Mohapatra RN . Model for large transition magnetic moment of the νe. Phys Rev Lett. (1989) 63:228. 10.1103/PhysRevLett.63.228

  • 169.

    Babu KS Mohapatra RN . Large transition magnetic moment of the neutrino from horizontal symmetry. Phys Rev. (1990) D42:377893. 10.1103/PhysRevD.42.3778

  • 170.

    Koide Y Ghosal A . Bimaximal neutrino mixing in a Zee type model with badly broken flavor symmetry. Phys Rev. (2001) D63:037301. 10.1103/PhysRevD.63.037301

  • 171.

    Kitabayashi T Yasue M . Large solar neutrino mixing and radiative neutrino mechanism. Phys Lett. (2002) B524:30818. 10.1016/S0370-2693(01)01368-5

  • 172.

    Adhikary B Brahmachari B Ghosal A Ma E Parida MK . A4 symmetry and prediction of Ue3 in a modified Altarelli-Feruglio model. Phys Lett. (2006) B638:3459. 10.1016/j.physletb.2006.05.051

  • 173.

    Fukuyama T Sugiyama H Tsumura K . Phenomenology in the zee model with the A4 symmetry. Phys Rev. (2011) D83:056016. 10.1103/PhysRevD.83.056016

  • 174.

    Aranda A Bonilla C Ramos R Rojas AD . A renormalizable fermion mass model with the double tetrahedral group. (2010) arXiv:1011.6470.

  • 175.

    Aranda A Bonilla C Ramos R Rojas AD . Model of flavor with quaternion symmetry. Phys Rev. (2011) D84:016009. 10.1103/PhysRevD.84.016009

  • 176.

    Okamoto Y Yasue M . Comment on radiative neutrino mass matrix with a sterile neutrino. Prog Theor Phys. (1999) 101:111927. 10.1143/PTP.101.1119

  • 177.

    Hall LJ Suzuki M . Explicit R-parity breaking in supersymmetric models. Nucl Phys. (1984) B231:41944. 10.1016/0550-3213(84)90513-3

  • 178.

    Nieves JF . Baryon and lepton number nonconserving processes and intermediate mass scales. Nucl Phys. (1981) B189:182204.

  • 179.

    Chua CK He XG Hwang WYP . Neutrino mass induced radiatively by supersymmetric leptoquarks. Phys Lett. (2000) B479:2249. 10.1016/S0370-2693(00)00325-7

  • 180.

    Mahanta U . Neutrino masses and mixing angles from leptoquark interactions. Phys Rev. (2000) D62:073009. 10.1103/PhysRevD.62.073009

  • 181.

    Aristizabal Sierra D Hirsch M Kovalenko SG . Leptoquarks: neutrino masses and accelerator phenomenology. Phys Rev. (2008) D77:055011. 10.1103/PhysRevD.77.055011

  • 182.

    Helo JC Hirsch M Ota T Pereira dos Santos FA . Double beta decay and neutrino mass models. J High Energy Phys. (2015) 5:92. 10.1007/JHEP05(2015)092

  • 183.

    Päs H Schumacher E . Common origin of RK and neutrino masses. Phys Rev. (2015) D92:114025. 10.1103/PhysRevD.92.114025

  • 184.

    Cheung K Nomura T Okada H . Testable radiative neutrino mass model without additional symmetries and explanation for the bs+ anomaly. Phys Rev. (2016) D94:115024. 10.1103/PhysRevD.94.115024

  • 185.

    Doršner I Fajfer S Košnik N . Leptoquark mechanism of neutrino masses within the grand unification framework. Eur Phys J. (2017) C77:417. 10.1140/epjc/s10052-017-4987-2

  • 186.

    Kong OCW . LR scalar mixings and oneloop neutrino masses. J High Energy Phys. (2000) 9:37. 10.1088/1126-6708/2000/09/037

  • 187.

    Davidson S Losada M . Neutrino masses in the Rp violating MSSM. J High Energy Phys. (2000) 5:21. 10.1088/1126-6708/2000/05/021

  • 188.

    Koide Y Sato J . R-parity violation in a SUSY GUT model and radiative neutrino masses. Phys Rev. (2003) D68:056004. 10.1103/PhysRevD.68.056004

  • 189.

    Koide Y . Radiative neutrino masses in a SUSY GUT model. Phys Lett. (2003) B574:828. 10.1016/j.physletb.2003.09.002

  • 190.

    Koide Y . Neutrino masses without seesaw mechanism in a SUSY SU(5) model with additional . Phys Lett. (2004) B595:46975. 10.1016/j.physletb.2004.05.072

  • 191.

    Longas R Portillo D Restrepo D Zapata O . The inert zee model. J High Energy Phys. (2016) 3:162. 10.1007/JHEP03(2016)162

  • 192.

    Ma E Raidal M Sarkar U . Neutrino masses in supersymmetry: R-parity and leptogenesis. Phys Lett. (1999) B460:35964. 10.1016/S0370-2693(99)00755-8

  • 193.

    Suematsu D Toma T Yoshida T . Neutrino masses and μ terms in a supersymmetric extra U(1) model. Int J Mod Phys. (2010) A25:403353. 10.1142/S0217751X1005024X

  • 194.

    Ma E . Supersymmetric U(1) gauge realization of the dark scalar doublet model of radiative neutrino mass. Mod Phys Lett. (2008) A23:7215. 10.1142/S0217732308026753

  • 195.

    Ma E . Axionic extensions of the supersymmetric standard model. Mod Phys Lett. (2009) A24:133542. 10.1142/S0217732309030837

  • 196.

    Fukuoka H Kubo J Suematsu D . Anomaly induced dark matter decay and PAMELA/ATIC experiments. Phys Lett. (2009) B678:4016. 10.1016/j.physletb.2009.06.048

  • 197.

    Suematsu D Toma T . Dark matter in the supersymmetric radiative seesaw model with an anomalous U(1) symmetry. Nucl Phys. (2011) B847:56789. 10.1016/j.nuclphysb.2011.02.007

  • 198.

    Fukuoka H Suematsu D Toma T . Signals of dark matter in a supersymmetric two dark matter model. J Cosmol Astropart Phys. (2011) 1107:001. 10.1088/1475-7516/2011/07/001

  • 199.

    Bhattacharya S Ma E Wegman D . Supersymmetric left-right model with radiative neutrino mass and multipartite dark matter. Eur Phys J. (2014) C74:2902. 10.1140/epjc/s10052-014-2902-7

  • 200.

    Kanemura S Machida N Shindou T Yamada T . A UV complete model for radiative seesaw scenarios and electroweak baryogenesis based on the supersymmetric gauge theory. Phys Rev. (2014) D89:013005. 10.1103/PhysRevD.89.013005

  • 201.

    R Figueiredo AJ . Neutrino masses from SUSY breaking in radiative seesaw models. Eur Phys J. (2015) C75:99. 10.1140/epjc/s10052-014-3244-1

  • 202.

    Kanemura S Machida N Shindou T . Radiative neutrino mass, dark matter and electroweak baryogenesis from the supersymmetric gauge theory with confinement. Phys Lett. (2014) B738:17886. 10.1016/j.physletb.2014.09.013

  • 203.

    Hundi RS . μ → eγ in a supersymmetric radiative neutrino mass model. Phys Rev. (2016) D93:015008. 10.1103/PhysRevD.93.015008

  • 204.

    Fraser S Ma E Popov O . Scotogenic inverse seesaw model of neutrino mass. Phys Lett. (2014) B737:2802. 10.1016/j.physletb.2014.08.069

  • 205.

    Restrepo D Rivera A Sánchez-Peláez M Zapata O Tangarife W . Radiative neutrino masses in the singlet-doublet fermion dark matter model with scalar singlets. Phys Rev. (2015) D92:013005. 10.1103/PhysRevD.92.013005

  • 206.

    Chiang C-W Okada H Senaha E . Dark matter, muon g−2, electric dipole moments, and in a one-loop induced neutrino model. Phys Rev. (2017) D96:015002. 10.1103/PhysRevD.96.015002

  • 207.

    Mohapatra RN . Mechanism for understanding small neutrino mass in superstring theories. Phys Rev Lett. (1986) 56:5613. 10.1103/PhysRevLett.56.561

  • 208.

    Mohapatra RN Valle JWF . Neutrino mass and baryon number nonconservation in superstring models. Phys Rev. (1986) D34:1642. 10.1103/PhysRevD.34.1642

  • 209.

    Ma E . Transformative A4 mixing of neutrinos with CP violation. Phys Rev. (2015) D92:051301. 10.1103/PhysRevD.92.051301

  • 210.

    Natale A . A radiative model of quark masses with binary tetrahedral symmetry. Nucl Phys. (2017) B914:20119. 10.1016/j.nuclphysb.2016.11.006

  • 211.

    Ma E . Soft A4Z3 symmetry breaking and cobimaximal neutrino mixing. Phys Lett. (2016) B755:34850. 10.1016/j.physletb.2016.02.032

  • 212.

    Ma E Sarkar U . Scalar neutrino as asymmetric dark matter: radiative neutrino mass and leptogenesis. Phys Rev. (2012) D85:075015. 10.1016/S0375-9601(96)00888-2

  • 213.

    Hirsch M Klapdor-Kleingrothaus HV Kovalenko SG . B-L violating masses in softly broken supersymmetry. Phys Lett. (1997) B398:3114. 10.1016/S0370-2693(97)00234-7

  • 214.

    Cai Y Clarke JD Schmidt MA Volkas RR . Testing radiative neutrino mass models at the LHC. J High Energy Phys. (2015) 2:161. 10.1007/JHEP02(2015)161

  • 215.

    Popov O White GA . One Leptoquark to unify them? Neutrino masses and unification in the light of (g−2)μ, and RK anomalies. Nucl Phys. (2017) B923:32438. 10.1016/j.nuclphysb.2017.08.007

  • 216.

    Ma E . Derivation of dark matter parity from lepton parity. Phys Rev Lett. (2015) 115:011801. 10.1103/PhysRevLett.115.011801

  • 217.

    Fraser S Kownacki C Ma E Popov O . Type II radiative seesaw model of neutrino mass with dark matter. Phys Rev. (2016) D93:013021. 10.1103/PhysRevD.93.013021

  • 218.

    Guo SY Han ZL Liao Y . Testing the type II radiative seesaw model: from dark matter detection to LHC signatures. Phys Rev. (2016) D94:115014. 10.1103/PhysRevD.94.115014

  • 219.

    Dasgupta B Ma E Tsumura K . Weakly interacting massive particle dark matter and radiative neutrino mass from Peccei-Quinn symmetry. Phys Rev. (2014) D89:041702. 10.1103/PhysRevD.89.041702

  • 220.

    Ma E . Syndetic model of fundamental interactions. Phys Lett. (2015) B741:2024. 10.1016/j.physletb.2014.12.045

  • 221.

    Ho SY Toma T Tsumura K . Systematic U(1)BL extensions of loop-induced neutrino mass models with dark matter. Phys Rev. (2016) D94:033007. 10.1103/PhysRevD.94.033007

  • 222.

    Ma E Pollard N Popov O Zakeri M . Gauge BL model of radiative neutrino mass with multipartite dark matter. Mod Phys Lett. (2016) A31:1650163. 10.1142/S0217732316501637

  • 223.

    Nomura T Okada H . Radiative neutrino mass in an alternative U(1)BL gauge symmetry. (2017) arXiv:1705.08309.

  • 224.

    Chang WF Wong CF . A model for neutrino masses and dark matter with the discrete gauge symmetry. Phys Rev. (2012) D85:013018. 10.1103/PhysRevD.85.013018

  • 225.

    Lindner M Schmidt D Watanabe A . Dark matter and U(1)' symmetry for the right-handed neutrinos. Phys Rev. (2014) D89:013007. 10.1103/PhysRevD.89.013007

  • 226.

    Kownacki C Ma E . Gauge U(1) dark symmetry and radiative light fermion masses. Phys Lett. (2016) B760:5962. 10.1016/j.physletb.2016.06.024

  • 227.

    Adhikari R Borah D Ma E . New U(1) gauge model of radiative lepton masses with sterile neutrino and dark matter. Phys Lett. (2016) B755:4147. 10.1016/j.physletb.2016.02.039

  • 228.

    Babu KS Mathur VS . Radiatively induced seesaw mechanism for neutrino masses. Phys Rev. (1988) D38:3550. 10.1103/PhysRevD.38.3550

  • 229.

    Fileviez Perez P Murgui C . Lepton flavour violation in left-right theory. Phys Rev. (2017) D95:075010. 10.1103/PhysRevD.95.075010

  • 230.

    Fileviez Perez P Murgui C Ohmer S . Simple left-right theory: lepton number violation at the LHC. Phys Rev. (2016) D94:051701. 10.1103/PhysRevD.94.051701

  • 231.

    Hue LT Arbuzov AB Ngan NTK Long HN . Probing neutrino and Higgs sectors in SU(2)1 × SU(2)2 × U(1)Y model with lepton-flavor non-universality. Eur Phys J. (2017) C77:346. 10.1140/epjc/s10052-017-4866-x

  • 232.

    Ma E . Gauged B−3Lτ and radiative neutrino masses. Phys Lett. (1998) B433:7481. 10.1016/S0370-2693(98)00599-1

  • 233.

    Ma E Roy DP . Phenomenology of the B−3Lτ gauge boson. Phys Rev. (1998) D58:095005.

  • 234.

    Dicus DA He HJ Ng JN . Neutrino - lepton masses, Zee scalars and muon g-2. Phys Rev Lett. (2001) 87:111803. 10.1103/PhysRevLett.87.111803

  • 235.

    Sahu N Sarkar U . Extended zee model for neutrino mass, leptogenesis and sterile neutrino like dark matter. Phys Rev. (2008) D78:115013. 10.1103/PhysRevD.78.115013

  • 236.

    Foot R Kobakhidze A McDonald KL Volkas RR . Neutrino mass in radiatively-broken scale-invariant models. Phys Rev. (2007) D76:075014. 10.1103/PhysRevD.76.075014

  • 237.

    Chen CS Chou CH . Neutrino masses, muon g-2, dark matter, lithium problem, and leptogenesis at TeV-scale. Phys Lett. (2011) B699:6873. 10.1016/j.physletb.2011.03.050

  • 238.

    Arbeláez C Cárcamo Hernández AE Kovalenko S Schmidt I . Radiative seesaw-type mechanism of fermion masses and non-trivial quark mixing. Eur Phys J. (2017) C77:422. 10.1140/epjc/s10052-017-4948-9

  • 239.

    Nomura T Okada H Orikasa Y . SU(2)L septet scalar linking to a radiative neutrino model. Phys Rev. (2016) D94:055012. 10.1103/PhysRevD.94.055012

  • 240.

    Nomura T Okada H . A model with isospin doublet U(1)D gauge symmetry. (2017) arXiv:1706.05268.

  • 241.

    Haba N Tsumura K . ν-two higgs doublet model and its collider phenomenology. J High Energy Phys. (2011) 6:68. 10.1007/JHEP06(2011)068

  • 242.

    Okada H Toma T . 3.55 keV X-ray line signal from excited dark matter in radiative neutrino model. Phys Lett. (2014) B737:1626. 10.1016/j.physletb.2014.08.046

  • 243.

    Okada H Orikasa Y . Classically conformal radiative neutrino model with gauged BL symmetry. Phys Lett. (2016) B760:55864. 10.1016/j.physletb.2016.07.039

  • 244.

    Ahriche A Manning A McDonald KL Nasri S . Scale-invariant models with one-loop neutrino mass and dark matter candidates. Phys Rev. (2016) D94:053005. 10.1103/PhysRevD.94.053005

  • 245.

    Ahriche A McDonald KL Nasri S . The scale-invariant scotogenic model. J High Energy Phys. (2016) 6:182. 10.1007/JHEP06(2016)182

  • 246.

    Kanemura S Seto O Shimomura T . Masses of dark matter and neutrino from TeV scale spontaneous U(1)BL breaking. Phys Rev. (2011) D84:016004. 10.1103/PhysRevD.84.016004

  • 247.

    Kajiyama Y Okada H Toma T . Light dark matter candidate in B-L gauged radiative inverse seesaw. Eur Phys J. (2013) C73:2381. 10.1140/epjc/s10052-013-2381-2

  • 248.

    Okada H Toma T . Fermionic dark matter in Radiative Inverse Seesaw Model with U(1)BL. Phys Rev. (2012) D86:033011. 10.1103/PhysRevD.86.033011

  • 249.

    Dasgupta A Hati C Patra S Sarkar U . A minimal model of TeV scale WIMPy leptogenesis. (2016) arXiv:1605.01292.

  • 250.

    Seto O Shimomura T . Atomki anomaly and dark matter in a radiative seesaw model with gauged BL symmetry. Phys Rev. (2017) D95:095032. 10.1103/PhysRevD.95.095032

  • 251.

    Machado ACB Montano J Pleite V Rodriguez MC . Scotogenic model with BL symmetry and exotic neutrinos. (2017) arXiv:1705.01009.

  • 252.

    Kashiwase S Suematsu D . Radiative neutrino mass model with degenerate right-handed neutrinos. Eur Phys J. (2016) C76:117. 10.1140/epjc/s10052-016-3964-5

  • 253.

    Yu JH . Hidden gauged U(1) model: unifying scotogenic neutrino and flavor dark matter. Phys Rev. (2016) D93:113007. 10.1103/PhysRevD.93.113007

  • 254.

    Ko P Nomura T Okada H . A flavor dependent gauge symmetry, predictive radiative seesaw and LHCb anomalies. Phys Lett. (2017) B772:54752. 10.1016/j.physletb.2017.07.021

  • 255.

    Baek S Okada H Yagyu K . Flavour dependent gauged radiative neutrino mass model. J High Energy Phys. (2015) 4:49. 10.1007/JHEP04(2015)049

  • 256.

    Baek S . Dark matter and muon (g−2) in local U(1)LμLτ-extended Ma Model. Phys Lett. (2016) B756:15. 10.1016/j.physletb.2016.02.062

  • 257.

    Kajiyama Y Kubo J Okada H . D(6) family symmetry and cold dark matter at LHC. Phys Rev. (2007) D75:033001. 10.1103/PhysRevD.75.033001

  • 258.

    Ahn YH Chen CS . Non-zero Ue3 and TeV-leptogenesis through A4 symmetry breaking. Phys Rev. (2010) D81:105013. 10.1103/PhysRevD.81.105013

  • 259.

    Ahn YH . Probing lepton flavor violation with nonzero Ue3 and leptogenesis through A4 symmetry breaking. (2010) arXiv:1006.2953.

  • 260.

    Ahn YH Okada H . Non-zero θ13 linking to dark matter from non-abelian discrete flavor model in radiative seesaw. Phys Rev. (2012) D85:073010. 10.1103/PhysRevD.85.073010

  • 261.

    Ma E Natale A Rashed A . Scotogenic A4 neutrino model for nonzero θ13 and large δCP. Int J Mod Phys. (2012) A27:1250134. 10.1142/S0217751X12501345

  • 262.

    Holthausen M Lindner M Schmidt MA . Lepton flavor at the electroweak scale: a complete A4 model. Phys Rev. (2013) D87:033006. 10.1103/PhysRevD.87.033006

  • 263.

    Mukherjee A Borah D Das MK . Common origin of non-zero θ13 and dark matter in an S4 flavour symmetric model with inverse seesaw. Phys Rev. (2017) D96:015014. 10.1103/PhysRevD.96.015014

  • 264.

    Borah D . Light sterile neutrino and dark matter in left-right symmetric models without a Higgs bidoublet. Phys Rev. (2016) D94:075024. 10.1103/PhysRevD.94.075024

  • 265.

    Hirsch M Lineros RA Morisi S Palacio J Rojas N Valle JWF . WIMP dark matter as radiative neutrino mass messenger. J High Energy Phys. (2013) 10:149. 10.1007/JHEP10(2013)149

  • 266.

    Rocha-Moran P Vicente A . Lepton flavor violation in the singlet-triplet scotogenic model. J High Energy Phys. (2016) 7:78. 10.1007/JHEP07(2016)078

  • 267.

    Merle A Platscher M Rojas N Valle JWF Vicente A . Consistency of WIMP dark matter as radiative neutrino mass messenger. J High Energy Phys. (2016) 7:13. 10.1007/JHEP07(2016)013

  • 268.

    Okada H Yagyu K . Radiative generation of lepton masses. Phys Rev. (2014) D89:053008. 10.1103/PhysRevD.89.053008

  • 269.

    Okada H Yagyu K . Radiative generation of lepton masses with the U(1)′ gauge symmetry. Phys Rev. (2014) D90:035019. 10.1103/PhysRevD.90.035019

  • 270.

    Wang W Han ZL . Radiative linear seesaw model, dark matter, and U(1)BL. Phys Rev. (2015) D92:095001. 10.1016/j.physletb.2011.02.014

  • 271.

    Nomura T Okada H . Loop induced type-II seesaw model and GeV dark matter with U(1)BL gauge symmetry. Phys Lett. (2017) B774:57581. 10.1016/j.physletb.2017.10.033

  • 272.

    Kanemura S Sugiyama H . Dark matter and a suppression mechanism for neutrino masses in the Higgs triplet model. Phys Rev. (2012) D86:073006. 10.1103/PhysRevD.86.073006

  • 273.

    Franceschini R Mohapatra RN . Radiatively induced type II seesaw models and vectorlike 5/3 charge quarks. Phys Rev. (2014) D89:055013. 10.1103/PhysRevD.89.055013

  • 274.

    Kanemura S Ota T . Neutrino masses from loop-induced d ≥ 7 operators. Phys Lett. (2011) B694:2337. 10.1016/j.physletb.2010.09.064

  • 275.

    Chang D Keung WY Pal PB . Spontaneous lepton number breaking at electroweak scale. Phys Rev Lett. (1988) 61:24203. 10.1103/PhysRevLett.61.2420

  • 276.

    Santamaria A . The hyperchargeless triplet majoron model. Phys Rev. (1989) D39:2715. 10.1103/PhysRevD.39.2715

  • 277.

    Law SSC McDonald KL . Inverse seesaw and dark matter in models with exotic lepton triplets. Phys Lett. (2012) B713:4904. 10.1016/j.physletb.2012.06.044

  • 278.

    Baldes I Bell NF Petraki K Volkas RR . Two radiative inverse seesaw models, dark matter, and baryogenesis. J Cosmol Astropart Phys. (2013) 1307:29. 10.1088/1475-7516/2013/07/029

  • 279.

    Aranda A Peinado E . A new radiative neutrino mass generation mechanism with higher dimensional scalar representations and custodial symmetry. Phys Lett. (2016) B754:113. 10.1016/j.physletb.2016.01.007

  • 280.

    Ahriche A Boucenna SM Nasri S . Dark radiative inverse seesaw mechanism. Phys Rev. (2016) D93:075036. 10.1103/PhysRevD.93.075036

  • 281.

    Pilaftsis A . Radiatively induced neutrino masses and large Higgs neutrino couplings in the standard model with Majorana fields. Z Phys. (1992) C55:27582. 10.1007/BF01482590

  • 282.

    Grimus W Neufeld H . Radiative neutrino masses in an SU(2) X U(1) model. Nucl Phys. (1989) B325:1832. 10.1016/0550-3213(89)90370-2

  • 283.

    Aristizabal Sierra D Yaguna CE . On the importance of the 1-loop finite corrections to seesaw neutrino masses. J High Energy Phys. (2011) 8:13. 10.1007/JHEP08(2011)013

  • 284.

    Grimus W Lavoura L . One-loop corrections to the seesaw mechanism in the multi-Higgs-doublet standard model. Phys Lett. (2002) B546:8695. 10.1016/S0370-2693(02)02672-2

  • 285.

    Dev PSB Pilaftsis A . Minimal radiative neutrino mass mechanism for inverse seesaw models. Phys Rev. (2012) D86:113001. 10.1103/PhysRevD.86.113001

  • 286.

    McDonald KL McKellar BHJ . Evaluating the two loop diagram responsible for neutrino mass in Babu's model. (2003) e-Print: hep-ph/0309270.

  • 287.

    Angel PW Cai Y Rodd NL Schmidt MA Volkas RR . Testable two-loop radiative neutrino mass model based on an LLQdcQdc effective operator. J High Energy Phys. (2013) 10:118. 10.1007/JHEP10(2013)118

  • 288.

    Zee A . Quantum numbers of majorana neutrino masses. Nucl Phys. (1986) B264:99110. 10.1016/0550-3213(86)90475-X

  • 289.

    Babu KS . Model of ‘calculable’ majorana neutrino masses. Phys Lett. (1988) B203:1326. 10.1016/0370-2693(88)91584-5

  • 290.

    Petcov ST Toshev ST . Conservation of lepton charges, massive majorana and massless neutrinos. Phys Lett. (1984) B143:1758. 10.1016/0370-2693(84)90829-3

  • 291.

    Babu KS Ma E . Natural hierarchy of radiatively induced majorana neutrino masses. Phys Rev Lett. (1988) 61:674. 10.1103/PhysRevLett.61.674

  • 292.

    Branco GC Grimus W Lavoura L . The seesaw mechanism in the presence of a conserved lepton number. Nucl Phys. (1989) B312:492508. 10.1016/0550-3213(89)90304-0

  • 293.

    Lavoura L . New model for the neutrino mass matrix. Phys Rev. (2000) D62:093011. 10.1103/PhysRevD.62.093011

  • 294.

    Kitabayashi T Yasue M . Neutrino oscillations induced by two loop radiative mechanism. Phys Lett. (2000) B490:23641. 10.1016/S0370-2693(00)00988-6

  • 295.

    Araki T Geng CQ . μ − τ symmetry in Zee-Babu model. Phys Lett. (2011) B694:1138. 10.1016/j.physletb.2010.09.046

  • 296.

    Aoki M Kanemura S Shindou T Yagyu K . An R-parity conserving radiative neutrino mass model without right-handed neutrinos. J High Energy Phys. (2010) 7:84. 10.1007/JHEP07(2010)084

  • 297.

    Haba N Shindou T . Tiny neutrino mass from SUSY and lepton number breaking sector. Phys Lett. (2011) B701:22933. 10.1016/j.physletb.2011.05.050

  • 298.

    Kohda M Sugiyama H Tsumura K . Lepton number violation at the LHC with leptoquark and diquark. Phys Lett. (2013) B718:143640. 10.1016/j.physletb.2012.12.048

  • 299.

    Nomura T Okada H . An extended colored zee-babu model. Phys Rev. (2016) D94:075021. 10.1103/PhysRevD.94.075021

  • 300.

    Chang WF Liou SC Wong CF Xu F . Charged lepton flavor violating processes and scalar leptoquark decay branching ratios in the colored Zee-Babu model. J High Energy Phys. (2016) 10:106. 10.1007/JHEP10(2016)106

  • 301.

    Ma E . Z3 dark matter and two-loop neutrino mass. Phys Lett. (2008) B662:4952. 10.1016/j.physletb.2008.02.053

  • 302.

    Ding R Han ZL Liao Y Xie WP . Radiative neutrino mass with ℤ3 dark matter: from relic density to LHC signatures. J High Energy Phys. (2016) 5:30. 10.1007/JHEP05(2016)030

  • 303.

    Nomura T Okada H . Generalized Zee–Babu model with 750 GeV diphoton resonance. Phys Lett. (2016) B756:295302. 10.1016/j.physletb.2016.03.034

  • 304.

    Okada H Toma T Yagyu K . Inert extension of the Zee-Babu model. Phys Rev. (2014) D90:095005. 10.1103/PhysRevD.90.095005

  • 305.

    Aoki M Toma T . Impact of semi-annihilation of ℤ3 symmetric dark matter with radiative neutrino masses. J Cosmol Astropart Phys. (2014) 1409:016. 10.1088/1475-7516/2014/09/016

  • 306.

    Ho S-Y Toma T Tsumura K . A radiative neutrino mass model with SIMP dark matter. J High Energy Phys. (2017) 7:101. 10.1007/JHEP07(2017)101

  • 307.

    Baek S Okada H Orikasa Y . A radiative neutrino model linking to a monochromatic gamma-ray line. (2017) arXiv:1703.00685.

  • 308.

    Wu Dd . The masses of the neutrinos in the unification gauge theories and the neutrino - anti-neutrino oscillations. Phys Rev. (1981) D23:2038. 10.1103/PhysRevD.23.2038

  • 309.

    Babu KS Ma E Pantaleone JT . Model of radiative neutrino masses: mixing and a possible fourth generation. Phys Lett. (1989) B218:2337. 10.1016/0370-2693(89)91425-1

  • 310.

    Aparici A Herrero-García J Rius N Santamaria A . Neutrino masses from new generations. J High Energy Phys. (2011) 7:122. 10.1007/JHEP07(2011)122

  • 311.

    Liao Y Liu JY Ning GZ . Radiative neutrino mass in Type III seesaw model. Phys Rev. (2009) D79:073003. 10.1103/PhysRevD.79.073003

  • 312.

    Babu KS Julio J . Two-loop neutrino mass generation through leptoquarks. Nucl Phys. (2010) B841:13056. 10.1016/j.nuclphysb.2010.07.022

  • 313.

    Cai Y Gargalionis J Schmidt MA Volkas RR . Reconsidering the one leptoquark solution: flavor anomalies and neutrino mass. J High Energy Phys. (2017) 10:047. 10.1007/JHEP10(2017)047

  • 314.

    Bamba K Geng CQ Ho SH . Radiative neutrino mass generation and dark energy. J Cosmol Astropart Phys. (2008) 0809:001. 10.1088/1475-7516/2008/09/001

  • 315.

    Porto RA Zee A . Neutrino mixing and the private higgs. Phys Rev. (2009) D79:013003. 10.1103/PhysRevD.79.013003

  • 316.

    Lindner M Schmidt D Schwetz T . Dark matter and neutrino masses from global U(1)BL symmetry breaking. Phys Lett. (2011) B705:32430. 10.1016/j.physletb.2011.10.022

  • 317.

    Chang D Keung WY . Higgs mediated neutrinoless double beta decay and neutrino mass in a majoron model. Phys Rev. (1989) D39:1386.

  • 318.

    Guo G He XG Li GN . Radiative two loop inverse seesaw and dark matter. J High Energy Phys. (2012) 10:044. 10.1007/JHEP10(2012)044

  • 319.

    Kajiyama Y Okada H Yagyu K . Two loop radiative seesaw model with inert triplet scalar field. Nucl Phys. (2013) B874:198216. 10.1016/j.nuclphysb.2013.05.020

  • 320.

    Nomura T Okada H . Two-loop induced majorana neutrino mass in a radiatively induced quark and lepton mass model. Phys Rev. (2016) D94:093006. 10.1103/PhysRevD.94.093006

  • 321.

    Aoki M Kubo J Takano H . Two-loop radiative seesaw mechanism with multicomponent dark matter explaining the possible γ excess in the Higgs boson decay and at the Fermi LAT. Phys Rev. (2013) D87:116001. 10.1103/PhysRevD.87.116001

  • 322.

    Aoki M Kubo J Takano H . Multicomponent dark matter in radiative seesaw model and monochromatic neutrino flux. Phys Rev. (2014) D90:076011. 10.1103/PhysRevD.90.076011

  • 323.

    Ma E Sarkar U . Revelations of the E(6)/U(1)(N) model: two-loop neutrino mass and dark matter. Phys Lett. (2007) B653:28891. 10.1016/j.physletb.2007.08.019

  • 324.

    Kajiyama Y Okada H Toma T . Multicomponent dark matter particles in a two-loop neutrino model. Phys Rev. (2013) D88:015029. 10.1103/PhysRevD.88.015029

  • 325.

    Baek S Okada H Toma T . Two loop neutrino model and dark matter particles with global B-L symmetry. J Cosmol Astropart Phys. (2014) 1406:027. 10.1088/1475-7516/2014/06/027

  • 326.

    Ghosh A Mondal T Mukhopadhyaya B . Heavy stable charged tracks as signatures of non-thermal dark matter at the LHC : a study in some non-supersymmetric scenarios. (2017) arXiv:1706.06815.

  • 327.

    Ma E . Radiative quark and lepton masses in a left-right gauge model. Phys Rev Lett. (1989) 63:1042. 10.1103/PhysRevLett.63.1042

  • 328.

    Nasri S Moussa S . Model for small neutrino masses at the TeV scale. Mod Phys Lett. (2002) A17:7718. 10.1142/S0217732302007119

  • 329.

    Chao W . Neutrino masses and A TeV scale seesaw mechanism. Phys Rev. (2010) D82:016008. 10.1103/PhysRevD.82.016008

  • 330.

    Ma E Wudka J . Vector-boson-induced neutrino mass. Phys Lett. (2012) B712:3915. 10.1016/j.physletb.2012.05.008

  • 331.

    Nomura T Okada H . Loop suppressed light fermion masses with U(1)R gauge symmetry. Phys. Rev. (2017) D96:015016. 10.1103/PhysRevD.96.015016

  • 332.

    Okada H Orikasa Y . Two-loop neutrino model with exotic leptons. Phys Rev. (2016) D93:013008. 10.1103/PhysRevD.93.013008

  • 333.

    Witten E . Neutrino masses in the minimal O(10) theory. Phys Lett. (1980) B91:814. 10.1016/0370-2693(80)90666-8

  • 334.

    Arbeláez Rodríguez C Kolešová H Malinský M . Witten's mechanism in the flipped SU(5) unification. Phys Rev. (2014) D89:055003. 10.1103/PhysRevD.89.055003

  • 335.

    Ma E . Radiative inverse seesaw mechanism for nonzero neutrino mass. Phys Rev. (2009) D80:013013.

  • 336.

    Babu KS Ma E . Radiative hierarchy of majorana neutrino masses. Phys Lett. (1989) B228:50812. 10.1016/0370-2693(89)90983-0

  • 337.

    Babu KS Julio J . Radiative neutrino mass generation through vector-like quarks. Phys Rev. (2012) D85:073005. 10.1103/PhysRevD.85.073005

  • 338.

    Kanemura S Nabeshima T Sugiyama H . TeV-scale seesaw with loop-induced dirac mass term and dark matter from U(1)BL gauge symmetry breaking. Phys Rev. (2012) D85:033004. 10.1103/PhysRevD.85.033004

  • 339.

    Ma E . Radiative origin of all quark and lepton masses through dark matter with flavor symmetry. Phys Rev Lett. (2014) 112:091801. 10.1103/PhysRevLett.112.091801

  • 340.

    Baek S Okada H Toma T . Radiative lepton model and dark matter with global U(1)′ symmetry. Phys Lett. (2014) B732:8590. 10.1016/j.physletb.2014.03.021

  • 341.

    Kanemura S Matsui T Sugiyama H . Neutrino mass and dark matter from gauged U(1)BL breaking. Phys Rev. (2014) D90:013001. 10.1103/PhysRevD.90.013001

  • 342.

    Kashiwase S Okada H Orikasa Y Toma T . Two loop neutrino model with dark matter and leptogenesis. Int J Mod Phys. (2016) A31:1650121. 10.1142/S0217751X16501219

  • 343.

    Nomura T Okada H Orikasa Y . Radiative neutrino mass in alternative left–right model. Eur Phys J. (2017) C77:103. 10.1140/epjc/s10052-017-4657-4

  • 344.

    Nomura T Okada H Orikasa Y . Radiative neutrino model with SU(2)L triplet fields. Phys Rev. (2016) D94:115018. 10.1103/PhysRevD.94.115018

  • 345.

    Chen CS Geng CQ Ng JN . Unconventional neutrino mass generation, neutrinoless double beta decays, and collider phenomenology. Phys Rev. (2007) D75:053004. 10.1103/PhysRevD.75.053004

  • 346.

    Chen CS Geng CQ Ng JN Wu JMS . Testing radiative neutrino mass generation at the LHC. J High Energy Phys. (2007) 8:22. 10.1088/1126-6708/2007/08/022

  • 347.

    del Aguila F Aparici A Bhattacharya S Santamaria A Wudka J . A realistic model of neutrino masses with a large neutrinoless double beta decay rate. J High Energy Phys. (2012) 5:133. 10.1007/JHEP05(2012)133

  • 348.

    Chen CS Geng CQ Huang D Tsai LH . Many high-charged scalars in LHC searches and Majorana neutrino mass generations. Phys Rev. (2013) D87:077702. 10.1103/PhysRevD.87.077702

  • 349.

    del Águila F Aparici A Bhattacharya S Santamaria A Wudka J . Neutrinoless double β decay with small neutrino masses. PoS (2013) Corfu2012:028. Available online at: https://arxiv.org/abs/1305.4900

  • 350.

    King SF Merle A Panizzi L . Effective theory of a doubly charged singlet scalar: complementarity of neutrino physics and the LHC. J High Energy Phys. (2014) 11:124. 10.1007/JHEP11(2014)124

  • 351.

    Geng CQ Tsai LH . Study of two-loop neutrino mass generation models. Ann Phys. (2016) 365:21022. 10.1016/j.aop.2015.11.010

  • 352.

    Appelquist T Shrock R . Neutrino masses in theories with dynamical electroweak symmetry breaking. Phys Lett. (2002) B548:20414. 10.1016/S0370-2693(02)02854-X

  • 353.

    Appelquist T Shrock R . Dynamical symmetry breaking of extended gauge symmetries. Phys Rev Lett. (2003) 90:201801. 10.1103/PhysRevLett.90.201801

  • 354.

    Appelquist T Piai M Shrock R . Fermion masses and mixing in extended technicolor models. Phys Rev. (2004) D69:015002. 10.1103/PhysRevD.69.015002

  • 355.

    Krauss LM Nasri S Trodden M . A Model for neutrino masses and dark matter. Phys Rev. (2003) D67:085002. 10.1103/PhysRevD.67.085002

  • 356.

    Chen CS McDonald KL Nasri S . A class of three-loop models with neutrino mass and dark matter. Phys Lett. (2014) B734:38893. 10.1016/j.physletb.2014.05.082

  • 357.

    Cheung K Seto O . Phenomenology of TeV right-handed neutrino and the dark matter model. Phys Rev. (2004) D69:113009. 10.1103/PhysRevD.69.113009

  • 358.

    Ahriche A Nasri S . Dark matter and strong electroweak phase transition in a radiative neutrino mass model. J Cosmol Astropart Phys. (2013) 1307:035. 10.1088/1475-7516/2013/07/035

  • 359.

    Ahriche A Nasri S Soualah R . Radiative neutrino mass model at the ee+ linear collider. Phys Rev. (2014) D89:095010. 10.1103/PhysRevD.89.095010

  • 360.

    Ahriche A Chen CS McDonald KL Nasri S . Three-loop model of neutrino mass with dark matter. Phys Rev. (2014) D90:015024. 10.1103/PhysRevD.90.015024

  • 361.

    Ahriche A McDonald KL Nasri S . A model of radiative neutrino mass: with or without dark matter. J High Energy Phys. (2014) 10:167. 10.1007/JHEP10(2014)167

  • 362.

    Ahriche A McDonald KL Nasri S Toma T . A model of neutrino mass and dark matter with an accidental symmetry. Phys Lett. (2015) B746:4305. 10.1016/j.physletb.2015.05.031

  • 363.

    Ahriche A McDonald KL Nasri S . A radiative model for the weak scale and neutrino mass via dark matter. J High Energy Phys. (2016) 2:38. 10.1007/JHEP02(2016)038

  • 364.

    Ahriche A McDonald KL Nasri S . Scalar sector phenomenology of three-loop radiative neutrino mass models. Phys Rev. (2015) D92:095020. 10.1103/PhysRevD.92.095020

  • 365.

    Nomura T Okada H Okada N . A colored KNT neutrino model. Phys Lett. (2016) B762:40914. 10.1016/j.physletb.2016.09.038

  • 366.

    Cheung K Nomura T Okada H . Three-loop neutrino mass model with a colored triplet scalar. Phys Rev. (2017) D95:015026. 10.1103/PhysRevD.95.015026

  • 367.

    Okada H Orikasa Y . X-ray line in radiative neutrino model with global U(1) symmetry. Phys Rev. (2014) D90:075023. 10.1103/PhysRevD.90.075023

  • 368.

    Okada H Yagyu K . Renormalizable model for neutrino mass, dark matter, muon g−2 and 750 GeV diphoton excess. Phys Lett. (2016) B756:33744. 10.1016/j.physletb.2016.03.040

  • 369.

    Cheung K Ishida H Okada H . Accommodation of the dirac phase in the Krauss-Nasri-Trodden model. (2016) arXiv:1609.06231.

  • 370.

    Aoki M Kanemura S Seto O . Neutrino mass, dark matter and baryon asymmetry via TeV-scale physics without fine-tuning. Phys Rev Lett. (2009) 102:051805. 10.1103/PhysRevLett.102.051805

  • 371.

    Gu PH . High-scale leptogenesis with three-loop neutrino mass generation and dark matter. J High Energy Phys. (2017) 4:159. 10.1007/JHEP04(2017)159

  • 372.

    Cheung K Nomura T Okada H . A three-loop neutrino model with leptoquark triplet scalars. Phys Lett. (2017) B768:35964. 10.1016/j.physletb.2017.03.021

  • 373.

    Culjak P Kumericki K Picek I . Scotogenic RνMDM at three-loop level. Phys Lett. (2015) B744:23743. 10.1016/j.physletb.2015.03.062

  • 374.

    Okada H Yagyu K . Three-loop neutrino mass model with doubly charged particles from isodoublets. Phys Rev. (2016) D93:013004. 10.1103/PhysRevD.93.013004

  • 375.

    Ko P Nomura T Okada H Orikasa Y . Confronting a new three-loop seesaw model with the 750 GeV diphoton excess. Phys Rev. (2016) D94:013009. 10.1103/PhysRevD.94.013009

  • 376.

    Gustafsson M No JM Rivera MA . Predictive model for radiatively induced neutrino masses and mixings with dark matter. Phys Rev Lett. (2013) 110:211802. 10.1103/PhysRevLett.110.211802

  • 377.

    Geng CQ Huang D Tsai LH . Loop-induced neutrino masses: a case study. Phys Rev. (2014) D90:113005. 10.1103/PhysRevD.90.113005

  • 378.

    Hatanaka H Nishiwaki K Okada H Orikasa Y . A three-loop neutrino model with global U(1) symmetry. Nucl Phys. (2015) B894:26883. 10.1016/j.nuclphysb.2015.03.006

  • 379.

    Alcaide J Das D Santamaria A . A model of neutrino mass and dark matter with large neutrinoless double beta decay. J High Energy Phys. (2017) 4:49. 10.1007/JHEP04(2017)049

  • 380.

    Nishiwaki K Okada H Orikasa Y . Three loop neutrino model with isolated k±±. Phys Rev. (2015) D92:093013. 10.1103/PhysRevD.92.093013

  • 381.

    Kanemura S Nishiwaki K Okada H Orikasa Y Park SC Watanabe R . LHC 750 GeV diphoton excess in a radiative seesaw model. PTEP (2016) 2016:123B04. 10.1093/ptep/ptw164

  • 382.

    Jin LG Tang R Zhang F . A three-loop radiative neutrino mass model with dark matter. Phys Lett. (2015) B741:1637. 10.1016/j.physletb.2014.12.034

  • 383.

    Nomura T Okada H Orikasa Y . Radiative seesaw model with degenerate majorana dark matter. Phys Rev. (2016) D93:113008. 10.1103/PhysRevD.93.113008

  • 384.

    Geng CQ Huang D . Large ν- oscillations from high-dimensional lepton number violating operator. J High Energy Phys. (2017) 3:103. 10.1007/JHEP03(2017)103

  • 385.

    Cirelli M Fornengo N Strumia A . Minimal dark matter. Nucl Phys. (2006) B753:17894. 10.1016/j.nuclphysb.2006.07.012

  • 386.

    Cirelli M Strumia A . Minimal dark matter: model and results. New J Phys. (2009) 11:105005. 10.1088/1367-2630/11/10/105005

  • 387.

    Gustafsson M No JM Rivera MA . Radiative neutrino mass generation linked to neutrino mixing and 0νββ-decay predictions. Phys Rev. (2014) D90:013012. 10.1103/PhysRevD.90.013012

  • 388.

    Liu Z Gu PH . Extending two Higgs doublet models for two-loop neutrino mass generation and one-loop neutrinoless double beta decay. Nucl Phys. (2017) B915:20623. 10.1016/j.nuclphysb.2016.12.001

  • 389.

    Cheng TP Li LF . On weak interaction induced neutrino oscillations. Phys Rev. (1978) D17:2375.

  • 390.

    Wang W Wang R Han ZL Han JZ . The BL scotogenic models for dirac neutrino masses. (2017) arXiv:1705.00414.

  • 391.

    Kanemura S Sakurai K Sugiyama H . Probing models of dirac neutrino masses via the flavor structure of the mass matrix. Phys Lett. (2016) B758:46572. 10.1016/j.physletb.2016.05.046

  • 392.

    Mohapatra RN . A model for dirac neutrino masses and mixings. Phys Lett. (1987) B198:6972. 10.1016/0370-2693(87)90161-4

  • 393.

    Mohapatra RN . Left-right symmetry and finite one loop dirac neutrino mass. Phys Lett. (1988) B201:51724.

  • 394.

    Balakrishna BS Mohapatra RN . Radiative fermion masses from new physics at Tev scale. Phys Lett. (1989) B216:34952.

  • 395.

    Gu PH Sarkar U . Radiative neutrino mass, dark matter and leptogenesis. Phys Rev. (2008) D77:105031. 10.1103/PhysRevD.77.105031

  • 396.

    Borah D Dasgupta A . Observable lepton number violation with predominantly dirac nature of active neutrinos. J High Energy Phys. (2017) 1:72. 10.1007/JHEP01(2017)072

  • 397.

    Borah D Dasgupta A . Naturally light dirac neutrino in left-right symmetric model. J Cosmol Astropart Phys. (2017) 1706:003. 10.1088/1475-7516/2017/06/003

  • 398.

    Rajpoot S . Model with calculable Dirac neutrino masses. Phys Rev. (1991) D44:2258. 10.1103/PhysRevD.44.225

  • 399.

    Kanemura S Nabeshima T Sugiyama H . Neutrino masses from loop-induced dirac yukawa couplings. Phys Lett. (2011) B703:6670. 10.1016/j.physletb.2011.07.047

  • 400.

    Farzan Y Ma E . Dirac neutrino mass generation from dark matter. Phys Rev. (2012) D86:033007. 10.1103/PhysRevD.86.033007

  • 401.

    Kanemura S Matsui T Sugiyama H . Loop suppression of dirac neutrino mass in the neutrinophilic two higgs doublet model. Phys Lett. (2013) B727:1516. 10.1016/j.physletb.2013.09.061

  • 402.

    Borah D Dasgupta A . Common origin of neutrino mass, dark matter and dirac leptogenesis. J Cosmol Astropart Phys. (2016) 1612:034. 10.1088/1475-7516/2016/12/034

  • 403.

    Bonilla C Ma E Peinado E Valle JWF . Two-loop dirac neutrino mass and WIMP dark matter. Phys Lett. (2016) B762:2148. 10.1016/j.physletb.2016.09.027

  • 404.

    Kanemura S Sakurai K Sugiyama H . Neutrino mass, dark matter and baryon asymmetry without lepton number violation. (2017) arXiv:1705.07040.

  • 405.

    Pilaftsis A . Anomalous fermion mass generation at three loops. Mod Phys Lett. (2013) A28:13500831. 10.1142/S0217732313500831

  • 406.

    Gutierrez DA Ponce WA Sanchez LA . Study of the SU(3)c × SU(3)L × U(1)X model with the minimal scalar sector. Int J Mod Phys. (2006) A21:221735. 10.1142/S0217751X06029442

  • 407.

    Okamoto Y Yasue M . Radiatively generated neutrino masses in SU(3)L × U(1)N gauge models. Phys Lett. (1999) B466:26773. 10.1016/S0370-2693(99)01134-X

  • 408.

    Kitabayashi T Yasue M . Nearly bimaximal neutrino mixing in an SU(3)L × U(1)N gauge model with radiative neutrino masses. (2000) e-Print: hep-ph/0006040.

  • 409.

    Kitabayashi T Yasue M . Radiatively induced neutrino masses and oscillations in an SU(3)L × U(1)N gauge model. Phys Rev. (2001) D63:095002. 10.1103/PhysRevD.63.095002

  • 410.

    Kitabayashi T Yasue M . Two loop radiative neutrino mechanism in an SU(3)L × U(1)N gauge model. Phys Rev. (2001) D63:095006. 10.1103/PhysRevD.63.095006

  • 411.

    Kitabayashi T . Comment on neutrino masses and oscillations in an SU(3)L × U(1)N model with radiative mechanism. Phys Rev. (2001) D64:057301. 10.1103/PhysRevD.64.057301

  • 412.

    Kitabayashi T Yasue M . The interplay between neutrinos and charged leptons in the minimal SU(3)L × U(1)N gauge model. Nucl Phys. (2001) B609:6182. 10.1016/S0550-3213(01)00276-0

  • 413.

    Kitabayashi T Yasue M . S2L permutation symmetry for left-handed μ and τ families and neutrino oscillations in an SU(3)L × U(1)N gauge model. Phys Rev. (2003) D67:015006. 10.1103/PhysRevD.67.015006

  • 414.

    Chang CHV Chang WF Ng JN . Neutrino masses in a 5-D SU(3)W TeV unification model. Phys Lett. (2003) B558:92102. 10.1016/S0370-2693(03)00219-3

  • 415.

    Chang WF Ng JN . Phenomenology of a 5-D orbifold SU(3)W unification model. Phys Rev. (2004) D69:056005. 10.1103/PhysRevD.69.056005

  • 416.

    Ponce WA Zapata O . Lepton masses and mixing without Yukawa hierarchies. Phys Rev. (2006) D74:093007. 10.1103/PhysRevD.74.093007

  • 417.

    Salazar JC Ponce WA Gutierrez DA . Phenomenology of the SU(3)c × SU(3)L × U(1)X model with exotic charged leptons. Phys Rev. (2007) D75:075016. 10.1103/PhysRevD.75.075016

  • 418.

    Dong PV Long HN . The economical SU(3)c × SU(3)L × U(1)Xmodel. Adv High Energy Phys. (2008) 2008:739492. 10.1155/2008/739492

  • 419.

    Boucenna SM Morisi S Valle JWF . Radiative neutrino mass in 3-3-1 scheme. Phys Rev. (2014) D90:013005. 10.1103/PhysRevD.90.013005

  • 420.

    Boucenna SM Fonseca RM Gonzalez-Canales F Valle JWF . Small neutrino masses and gauge coupling unification. Phys Rev. (2015) D91:031702. 10.1103/PhysRevD.91.031702

  • 421.

    Deppisch FF Hati C Patra S Sarkar U Valle JWF . 331 models and grand unification: from minimal SU(5) to minimal SU(6). Phys Lett. (2016) B762:43240. 10.1016/j.physletb.2016.10.002

  • 422.

    Okada H Okada N Orikasa Y . Radiative seesaw mechanism in a minimal 3-3-1 model. Phys Rev. (2016) D93:073006. 10.1103/PhysRevD.93.073006

  • 423.

    Cárcamo Hernández AE Long HN . A highly predictive A4 flavour 3-3-1 model with radiative inverse seesaw mechanism. (2017) arXiv:1705.05246.

  • 424.

    Kitabayashi T Yasue M . Large mixing angle MSW solution in an SU(3)L × U(1)N gauge model with two loop radiative mechanism. Phys Lett. (2001) B508:8595. 10.1016/S0370-2693(01)00397-5

  • 425.

    Aizawa I Ishiguro M Kitabayashi T Yasue M . Bilarge neutrino mixing and mu - tau permutation symmetry for two-loop radiative mechanism. Phys Rev. (2004) D70:015011. 10.1103/PhysRevD.70.015011

  • 426.

    Sanchez LA Ponce WA Mira JM . Hierarchical radiative masses in a supersymmetric three-family model without Higgsinos. Eur Phys J. (2005) C42:20515. 10.1140/epjc/s2005-02286-9

  • 427.

    Chang D Long HN . Interesting radiative patterns of neutrino mass in an SU(3)c × SU(3)L × U(1)X model with right-handed neutrinos. Phys Rev. (2006) D73:053006. 10.1103/PhysRevD.73.053006

  • 428.

    Benavides RH Ponce WA Giraldo Y . SU(3)cSU(3)LU(1)X models with four families. Phys Rev. (2010) D82:013004. 10.1103/PhysRevD.82.013004

  • 429.

    Cárcamo Hernández AE Kovalenko S Long HN Schmidt I . A novel 3-3-1 model for the generation of the SM fermion mass and mixing pattern. (2017) arXiv:1705.09169.

  • 430.

    Cárcamo Hernández AE . A novel and economical explanation for SM fermion masses and mixings. Eur Phys J. (2016) C76:503. 10.1140/epjc/s10052-016-4351-y

  • 431.

    Fonseca RM Hirsch M . Lepton number violation in 331 models. Phys Rev. (2016) D94:115003. 10.1103/PhysRevD.94.115003

  • 432.

    Antusch S Biggio C Fernandez-Martinez E Gavela MB Lopez-Pavon J . Unitarity of the leptonic mixing matrix. J High Energy Phys. (2006) 10:084. 10.1088/1126-6708/2006/10/084

  • 433.

    Abada A Biggio C Bonnet F Gavela MB Hambye T . Low energy effects of neutrino masses. J High Energy Phys. (2007) 12:061. 10.1088/1126-6708/2007/12/061

  • 434.

    Bilenky MS Santamaria A . One loop effective Lagrangian for a standard model with a heavy charged scalar singlet. Nucl Phys. (1994) B420:4793. 10.1016/0550-3213(94)90375-1

  • 435.

    Nebot M Oliver JF Palao D Santamaria A . Prospects for the Zee-Babu Model at the CERN LHC and low energy experiments. Phys Rev. (2008) D77:093013. 10.1103/PhysRevD.77.093013

  • 436.

    Farzan Y . A model for large non-standard interactions of neutrinos leading to the LMA-Dark solution. Phys Lett. (2015) B748:3115. 10.1016/j.physletb.2015.07.015

  • 437.

    Farzan Y Shoemaker IM . Lepton flavor violating non-standard interactions via light mediators. J High Energy Phys. (2016) 7:33. 10.1007/JHEP07(2016)033

  • 438.

    Davidson S Pena-Garay C Rius N Santamaria A . Present and future bounds on nonstandard neutrino interactions. J High Energy Phys. (2003) 3:11. 10.1088/1126-6708/2003/03/011

  • 439.

    Ibarra A Masso E Redondo J . Systematic approach to gauge-invariant relations between lepton flavor violating processes. Nucl Phys. (2005) B715:52335. 10.1016/j.nuclphysb.2005.03.017

  • 440.

    Gavela MB Hernandez D Ota T Winter W . Large gauge invariant non-standard neutrino interactions. Phys Rev. (2009) D79:013007. 10.1103/PhysRevD.79.013007

  • 441.

    Biggio C Blennow M Fernandez-Martinez E . General bounds on non-standard neutrino interactions. J High Energy Phys. (2009) 8:90. 10.1088/1126-6708/2009/08/090

  • 442.

    Biggio C Blennow M Fernandez-Martinez E . Loop bounds on non-standard neutrino interactions. J High Energy Phys. (2009) 3:139. 10.1088/1126-6708/2009/03/139

  • 443.

    Antusch S Baumann JP Fernandez-Martinez E . Non-standard neutrino interactions with matter from physics beyond the standard model. Nucl Phys. (2009) B810:36988. 10.1016/j.nuclphysb.2008.11.018

  • 444.

    Ohlsson T . Status of non-standard neutrino interactions. Rept Prog Phys. (2013) 76:044201. 10.1088/0034-4885/76/4/044201

  • 445.

    Porod W Staub F Vicente A . A flavor kit for BSM models. Eur Phys J. (2014) C74:2992. 10.1140/epjc/s10052-014-2992-2

  • 446.

    Lee BW Pakvasa S Shrock RE Sugawara H . Muon and electron number nonconservation in a v-a six quark model. Phys Rev Lett. (1977) 38:937. 10.1103/PhysRevLett.38.937

  • 447.

    Lee BW Shrock RE . Natural suppression of symmetry violation in gauge theories: Muon - Lepton and electron lepton number nonconservation. Phys Rev. (1977) D16:1444. 10.1103/PhysRevD.16.1444

  • 448.

    Marciano WJ Sanda AI . Exotic decays of the muon and heavy leptons in gauge theories. Phys Lett. (1977) 67B:3035.

  • 449.

    Carpentier M Davidson S . Constraints on two-lepton, two quark operators. Eur Phys J. (2010) C70:107190. 10.1140/epjc/s10052-010-1482-4

  • 450.

    Pruna GM Signer A . The μ → decay in a systematic effective field theory approach with dimension 6 operators. J High Energy Phys. (2014) 10:014. 10.1007/JHEP10(2014)014

  • 451.

    Ghosal A Koide Y Fusaoka H . Lepton flavor violating Z decays in the Zee model. Phys Rev. (2001) D64:053012. 10.1103/PhysRevD.64.053012

  • 452.

    Li B Liao Y Ma XD . Charged Lepton flavor-violating transitions in color octet model. Eur Phys J. (2016) C76:615. 10.1140/epjc/s10052-016-4462-5

  • 453.

    Baldini AM Bao Y Baracchini E Bemporad C Berg F Biasotti M et al . Search for the lepton flavour violating decay μ+e+γ with the full dataset of the MEG experiment. Eur Phys J. (2016) C76:434. 10.1140/epjc/s10052-016-4271-x

  • 454.

    Baldini AM Cei F Cerri C Dussoni S Galli L Grassi M et al . MEG upgrade proposal. (2013) arXiv:1301.7225.

  • 455.

    Aubert B Karyotakis Y Lees JP Poireau V Prencipe E Prudent X et al . Searches for lepton flavor violation in the decays τ±e±γ and τ± → μ±γ. Phys Rev Lett. (2010) 104:021802. 10.1103/PhysRevLett.104.021802

  • 456.

    Aushev T Bartel W Bondar A Brodzicka J Browder TE Chang P et al . Physics at super B factory. (2010) arXiv:1002.5012.

  • 457.

    Hisano J Moroi T Tobe K Yamaguchi M . Lepton flavor violation via right-handed neutrino Yukawa couplings in supersymmetric standard model. Phys Rev. (1996) D53:244259. 10.1103/PhysRevD.53.2442

  • 458.

    Lavoura L . General formulae for f(1) → f(2)γ. Eur Phys J. (2003) C29:1915. 10.1140/epjc/s2003-01212-7

  • 459.

    Schmidt D Schwetz T Toma T . Direct detection of leptophilic dark matter in a model with radiative neutrino masses. Phys Rev. (2012) D85:073009. 10.1103/PhysRevD.85.073009

  • 460.

    Vicente A Yaguna CE . Probing the scotogenic model with lepton flavor violating processes. J High Energy Phys. (2015) 2:144. 10.1007/JHEP02(2015)144

  • 461.

    Cheung Km Kong OCW . μ → from supersymmetry without R-parity. Phys Rev. (2001) D64:095007.

  • 462.

    Abada A Davidson S Losada M . Neutrino masses and mixings in the MSSM with soft bilinear R(p) violation. Phys Rev. (2002) D65:075010. 10.1103/PhysRevD.65.075010

  • 463.

    Carvalho DF Gomez ME Romao JC . Charged lepton flavor violation in supersymmetry with bilinear R-parity violation. Phys Rev. (2002) D65:093013. 10.1103/PhysRevD.65.093013

  • 464.

    Blondel A Bravar A Pohl M Bachmann S Berger N Kiehn M et al . Research proposal for an experiment to search for the decay μ → eee. (2013) arXiv:1301.6113.

  • 465.

    Bellgardt U Otter G Eichler R Felawka L Niebuhr C Walter HK et al . Search for the decay μ+e+e+e. Nucl Phys. (1988) B299:16.

  • 466.

    Hayasaka K Inami K Miyazaki Y Arinstein K Aulchenko VM Aushev T et al . Search for lepton flavor violating τ decays into three leptons with 719 million produced τ+τ pairs. Phys Lett. (2010) B687:13943. 10.1016/j.physletb.2010.03.037

  • 467.

    Abada A Krauss ME Porod W Staub F Vicente A Weiland C . Lepton flavor violation in low-scale seesaw models: SUSY and non-SUSY contributions. J High Energy Phys. (2014) 11:48. 10.1007/JHEP11(2014)048

  • 468.

    Kubo J Ma E Suematsu D . Cold dark matter, radiative neutrino mass, μ → , and neutrinoless double beta decay. Phys Lett. (2006) B642:1823. 10.1016/j.physletb.2006.08.085

  • 469.

    Aristizabal Sierra D Kubo J Restrepo D Suematsu D Zapata O . Radiative seesaw: warm dark matter, collider and lepton flavour violating signals. Phys Rev. (2009) D79:013011. 10.1103/PhysRevD.79.013011

  • 470.

    Suematsu D Toma T Yoshida T . Reconciliation of CDM abundance and μ → in a radiative seesaw model. Phys Rev. (2009) D79:093004. 10.1103/PhysRevD.79.093004

  • 471.

    Adulpravitchai A Lindner M Merle A . Confronting flavour symmetries and extended scalar sectors with Lepton flavour violation bounds. Phys Rev. (2009) D80:055031. 10.1103/PhysRevD.80.055031

  • 472.

    Toma T Vicente A . Lepton flavor violation in the scotogenic model. J High Energy Phys. (2014) 1:160. 10.1007/JHEP01(2014)160

  • 473.

    Cai Y Schmidt MA . Revisiting the RνMDM models. J High Energy Phys. (2016) 5:28. 10.1007/JHEP05(2016)028

  • 474.

    Dohmen C Groth K-D Heer B Honecker W Otter G Steinrücken B et al . Test of lepton flavor conservation in μ → e conversion on titanium. Phys Lett. (1993) B317:6316. 10.1016/0370-2693(93)91383-X

  • 475.

    The PRIME Working Group Collaboration. Search for the μ → e Conversion Process at an Ultimate Sensitivity of the Order of 10−18 with PRISM. LOI to J-PARC 50-GeV PS, LOI-25. Available online at: http://www-ps.kek.jp/jhf-np/LOIlist/pdf/L25.pdf

  • 476.

    Bertl W Engfer R Hermes EA Kurz G Kozlowski T Kuth J et al . A search for muon to electron conversion in muonic gold. Eur Phys J. (2006) C47:33746. 10.1140/epjc/s2006-02582-x

  • 477.

    Pezzullo G . The Mu2e experiment at Fermilab: a search for lepton flavor violation. Nucl Part Phys Proc. (2017) 285–286:37. 10.1016/j.nuclphysbps.2017.03.002

  • 478.

    Natori H . DeeMe experiment - an experimental search for a mu-e conversion reaction at J-PARC MLF. Nucl Phys Proc Suppl. (2014) 248–250:527. 10.1016/j.nuclphysbps.2014.02.010

  • 479.

    Kuno Y Okada Y . Muon decay and physics beyond the standard model. Rev Mod Phys. (2001) 73:151202. 10.1103/RevModPhys.73.151

  • 480.

    Arganda E Herrero MJ Teixeira AM . mu-e conversion in nuclei within the CMSSM seesaw: universality versus non-universality. J High Energy Phys. (2007) 10:104. 10.1088/1126-6708/2007/10/104

  • 481.

    Chiang HC Oset E Kosmas TS Faessler A Vergados JD . Coherent and incoherent (μ, e) conversion in nuclei. Nucl Phys. (1993) A559:52642.

  • 482.

    Kosmas TS Kovalenko S Schmidt I . Nuclear muon- e- conversion in strange quark sea. Phys Lett. (2001) B511:203. 10.1016/S0370-2693(01)00657-8

  • 483.

    Albrecht J Artuso M Babu K Bernstein RH Blum T Brown DN et al . Working Group Report: Charged Leptons. In: Proceedings, 2013 Community Summer Study on the Future of U.S. Particle Physics: Snowmass on the Mississippi (CSS2013), July 29-August 6, 2013. Minneapolis, MN (2013). Available online at: http://inspirehep.net/record/1265506/files/arXiv:1311.5278.pdf

  • 484.

    de Gouvea A Vogel P . Lepton flavor and number conservation, and physics beyond the standard model. Prog Part Nucl Phys. (2013) 71:7592. 10.1016/j.ppnp.2013.03.006

  • 485.

    Crivellin A Davidson S Pruna GM Signer A . Renormalisation-group improved analysis of μ → e processes in a systematic effective-field-theory approach. J High Energy Phys. (2017) 5:117. 10.1007/JHEP05(2017)117

  • 486.

    de Gouvea A Lola S Tobe K . Lepton flavor violation in supersymmetric models with trilinear R-parity violation. Phys Rev. (2001) D63:035004. 10.1103/PhysRevD.63.035004

  • 487.

    Klasen M Lamprea DR Yaguna CE . Lepton flavor violation and scalar dark matter in a radiative model of neutrino masses. (2016) arXiv:1602.05137.

  • 488.

    Doršner I Fajfer S Greljo A Kamenik JF Košnik N Nišandžic I . New physics models facing lepton flavor violating higgs decays at the percent level. J High Energy Phys. (2015) 6:108. 10.1007/JHEP06(2015)108

  • 489.

    Blankenburg G Ellis J Isidori G . Flavour-changing decays of a 125 GeV higgs-like particle. Phys Lett. (2012) B712:38690. 10.1016/j.physletb.2012.05.007

  • 490.

    Harnik R Kopp J Zupan J . Flavor violating higgs decays. J High Energy Phys. (2013) 3:26. 10.1007/JHEP03(2013)026

  • 491.

    Herrero-García J Rius N Santamaria A . Higgs lepton flavour violation: UV completions and connection to neutrino masses. J High Energy Phys. (2016) 11:84. 10.1007/JHEP11(2016)084

  • 492.

    Aad G Abbott B Abdallah J Abdinov O Abeloos B Aben R et al . Search for lepton-flavour-violating decays of the Higgs and Z bosons with the ATLAS detector. Eur Phys J. (2017) C77:70. 10.1140/epjc/s10052-017-4624-0

  • 493.

    CMS Collaboration . Search for lepton flavour violating decays of the Higgs boson to μτ and eτ in proton-proton collisions at . CMS-PAS-HIG-17-001 (2017).

  • 494.

    Pilaftsis A . Lepton flavor nonconservation in H0 decays. Phys Lett. (1992) B285:6874.

  • 495.

    Arganda E Curiel AM Herrero MJ Temes D . Lepton flavor violating Higgs boson decays from massive seesaw neutrinos. Phys Rev. (2005) D71:035011. 10.1103/PhysRevD.71.035011

  • 496.

    Arganda E Herrero MJ Marcano X Weiland C . Imprints of massive inverse seesaw model neutrinos in lepton flavor violating Higgs boson decays. Phys Rev. (2015) D91:015001. 10.1103/PhysRevD.91.015001

  • 497.

    Thao NH Hue LT Hung HT Xuan NT . Lepton flavor violating Higgs boson decays in seesaw models: new discussions. Nucl Phys. (2017) B921:15980. 10.1016/j.nuclphysb.2017.05.014

  • 498.

    Arganda E Herrero MJ Marcano X Morales R Szynkman A . Effective lepton flavor violating HâŞiâṢj vertex from right-handed neutrinos within the mass insertion approximation. Phys Rev. (2017) D95:095029. 10.1103/PhysRevD.95.095029

  • 499.

    Davidson S Grenier GJ . Lepton flavour violating Higgs and tau to mu gamma. Phys Rev. (2010) D81:095016. 10.1103/PhysRevD.81.095016

  • 500.

    Aristizabal Sierra D Vicente A . Explaining the CMS higgs flavor violating decay excess. Phys Rev. (2014) D90:115004. 10.1103/PhysRevD.90.115004

  • 501.

    Omura Y Senaha E Tobe K . Lepton-flavor-violating Higgs decay h → μτ and muon anomalous magnetic moment in a general two Higgs doublet model. J High Energy Phys. (2015) 5:28. 10.1007/JHEP05(2015)028

  • 502.

    Botella FJ Branco GC Nebot M Rebelo MN . Flavour changing higgs couplings in a class of two higgs doublet models. Eur Phys J. (2016) C76:161. 10.1140/epjc/s10052-016-3993-0

  • 503.

    Raidal M van der Schaaf A Bigi I Mangano ML Semertzidis Y Abel S et al . Flavour physics of leptons and dipole moments. Eur Phys J. (2008) C57:13182. 10.1140/epjc/s10052-008-0715-2

  • 504.

    Patrignani C Agashe K Aielli G Amsler C Antonelli M Asner DM et al . Review of particle physics. Chin Phys. (2016) C40:100001. 10.1088/1674-1137/40/10/100001

  • 505.

    Lee S Nomura T Okada H . Radiatively induced neutrino mass model with flavor dependent gauge symmetry. (2017) arXiv:1702.03733.

  • 506.

    Lindner M Platscher M Queiroz FS . A call for new physics : the muon anomalous magnetic moment and lepton flavor violation. (2016) arXiv:1610.06587.

  • 507.

    Furry WH . On transition probabilities in double beta-disintegration. Phys Rev. (1939) 56:118493. 10.1103/PhysRev.56.1184

  • 508.

    Albert JB Auty DJ Barbeau PS Beauchamp E Beck D Belov V et al . Search for Majorana neutrinos with the first two years of EXO-200 data. Nature (2014) 510:22934. 10.1038/nature13432

  • 509.

    Aad G Abajyan T Abbott B Abdallah J Abdel Khalek S Abdinov O et al . Limit on neutrinoless ββ decay of 136Xe from the first phase of KamLAND-Zen and comparison with the positive claim in 76Ge. Phys Rev Lett. (2013) 110:062502. 10.1103/PhysRevLett.110.062502

  • 510.

    Asakura K Gando A Gando Y Hachiya T Hayashida S Ikeda H et al . Results from KamLAND-Zen. AIP Conf Proc. (2015) 1666:170003. 10.1063/1.4915593

  • 511.

    Dell'Oro S Marcocci S Viel M Vissani F . Neutrinoless double beta decay: 2015 review. Adv High Energy Phys. (2016) 2016:2162659. 10.1155/2016/2162659

  • 512.

    Schechter J Valle JWF . Neutrinoless double beta decay in SU(2) × U(1) theories. Phys Rev. (1982) D25:2951.

  • 513.

    Duerr M Lindner M Merle A . On the quantitative impact of the schechter-valle theorem. J High Energy Phys. (2011) 6:91. 10.1007/JHEP06(2011)091

  • 514.

    Rodejohann W . Neutrinoless double beta decay and neutrino physics. J Phys. (2012) G39:124008. 10.1088/0954-3899/39/12/124008

  • 515.

    Bilenky SM Giunti C . Neutrinoless double-beta decay: a brief review. Mod Phys Lett. (2012) A27:1230015. 10.1142/S0217732312300157

  • 516.

    Päs H Hirsch M Klapdor-Kleingrothaus HV Kovalenko SG . Towards a superformula for neutrinoless double beta decay. Phys Lett. (1999) B453:19498. 10.1016/S0370-2693(99)00330-5

  • 517.

    Päs H Hirsch M Klapdor-Kleingrothaus HV Kovalenko SG . A Superformula for neutrinoless double beta decay. 2. The short range part. Phys Lett. (2001) B498:359. 10.1016/S0370-2693(00)01359-9

  • 518.

    Hirsch M Klapdor-Kleingrothaus HV Kovalenko SG . New leptoquark mechanism of neutrinoless double beta decay. Phys Rev. (1996) D54:R4207210.

  • 519.

    Hirsch M Klapdor-Kleingrothaus HV Kovalenko SG . New low-energy leptoquark interactions. Phys Lett. (1996) B378:1722.

  • 520.

    Mohapatra RN . New contributions to neutrinoless double beta decay in supersymmetric theories. Phys Rev. (1986) D34:345761. 10.1103/PhysRevD.34.3457

  • 521.

    Babu KS Mohapatra RN . New vector - scalar contributions to neutrinoless double beta decay and constraints on R-parity violation. Phys Rev Lett. (1995) 75:227679. 10.1103/PhysRevLett.75.2276

  • 522.

    Hirsch M Klapdor-Kleingrothaus HV Kovalenko SG . Supersymmetry and neutrinoless double beta decay. Phys Rev. (1996) D53:132948.

  • 523.

    Hirsch M Klapdor-Kleingrothaus HV Kovalenko SG . New supersymmetric contributions to neutrinoless double beta decay. Phys Lett. (1995) B352:17. 10.1016/0370-2693(95)00460-3

  • 524.

    Choubey S Duerr M Mitra M Rodejohann W . Lepton Number and Lepton Flavor Violation through Color Octet States. J High Energy Phys. (2012) 5:17. 10.1007/JHEP05(2012)017

  • 525.

    Brahmachari B Ma E . Neutrinoless double beta decay with negligible neutrino mass. Phys Lett. (2002) B536:25962. 10.1016/S0370-2693(02)01869-5

  • 526.

    Gu PH . Significant neutrinoless double beta decay with quasi-Dirac neutrinos. Phys Rev. (2012) D85:093016. 10.1103/PhysRevD.85.093016

  • 527.

    Quintero N . Constraints on lepton number violating short-range interactions from |ΔL| = 2 processes. Phys Lett. (2017) B764:605. 10.1016/j.physletb.2016.10.056

  • 528.

    Rodejohann W . Neutrino oscillation experiments and limits on lepton number and lepton flavor violating processes. Phys Rev. (2000) D62:013011. 10.1103/PhysRevD.62.013011

  • 529.

    Berryman JM de Gouvêa A Kelly KJ Kobach A . Lepton-number-violating searches for muon to positron conversion. Phys Rev. (2017) D95:115010. 10.1103/PhysRevD.95.115010

  • 530.

    Geib T Merle A . μe+ conversion from short-range operators. Phys Rev. (2017) D95:055009. 10.1103/PhysRevD.95.055009

  • 531.

    Deppisch FF Bhupal Dev PS Pilaftsis A . Neutrinos and collider physics. New J Phys. (2015) 17:075019. 10.1088/1367-2630/17/7/075019

  • 532.

    Babu KS Jana S . Probing doubly charged higgs bosons at the LHC through photon initiated processes. Phys Rev. (2017) D95:055020. 10.1103/PhysRevD.95.055020

  • 533.

    The ATLAS Collaboration . Search for doubly-charged Higgs bosons in same-charge electron pair final states using proton-proton collisions at with the ATLAS detector. ATLAS-CONF-2016-051. (2016).

  • 534.

    Collaboration C . A search for doubly-charged Higgs boson production in three and four lepton final states at . CMS-PAS-HIG-16-036. (2017).

  • 535.

    Sugiyama H Tsumura K Yokoya H . Discrimination of models including doubly charged scalar bosons by using tau lepton decay distributions. Phys Lett. (2012) B717:22934. 10.1016/j.physletb.2012.09.044

  • 536.

    del Águila F Chala M . LHC bounds on Lepton Number Violation mediated by doubly and singly-charged scalars. J High Energy Phys. (2014) 3:27. 10.1007/JHEP03(2014)027

  • 537.

    del Aguila F Chala M Santamaria A Wudka J . Discriminating between lepton number violating scalars using events with four and three charged leptons at the LHC. Phys Lett. (2013) B725:3105. 10.1016/j.physletb.2013.07.014

  • 538.

    Kanemura S Yagyu K Yokoya H . First constraint on the mass of doubly-charged Higgs bosons in the same-sign diboson decay scenario at the LHC. Phys Lett. (2013) B726:3169. 10.1016/j.physletb.2013.08.054

  • 539.

    Helo JC Hirsch M Päs H Kovalenko SG . Short-range mechanisms of neutrinoless double beta decay at the LHC. Phys Rev. (2013) D88:073011. 10.1103/PhysRevD.88.073011

  • 540.

    Peng T Ramsey-Musolf MJ Winslow P . TeV lepton number violation: from neutrinoless double-β decay to the LHC. Phys Rev. (2016) D93:093002. 10.1103/PhysRevD.93.093002

  • 541.

    Gonzales L Helo JC Hirsch M Kovalenko SG . Scalar-mediated double beta decay and LHC. J High Energy Phys. (2016) 12:130. 10.1007/JHEP12(2016)130

  • 542.

    Buchmuller W Ruckl R Wyler D . Leptoquarks in Lepton - Quark collisions. Phys Lett. (1987) B191:4428. 10.1016/0370-2693(87)90637-X

  • 543.

    Aad G Abbott B Abdallah J Abdinov O Aben R Abolins M et al . Search for direct top-squark pair production in final states with two leptons in pp collisions at 8TeV with the ATLAS detector. J High Energy Phys. (2014) 6:124. 10.1007/JHEP06(2014)124

  • 544.

    CMS Collaboration . Search for direct stop pair production in the dilepton final state at 13 TeV. CMS-PAS-SUS-17-001. (2017).

  • 545.

    Carrà S . Search for top squark with two leptons in the final state at LHC Run 2 with the ATLAS detector. Nuovo Cim. (2017) C40:26. 10.1393/ncc/i2017-17026-5

  • 546.

    Cai Y Schmidt MA . A case study of the sensitivity to LFV operators with precision measurements and the LHC. J High Energy Phys. (2016) 2:176. 10.1007/JHEP02(2016)176

  • 547.

    Aguilar-Saavedra JA . Identifying top partners at LHC. J High Energy Phys. (2009) 11:030. 10.1088/1126-6708/2009/11/030

  • 548.

    Aad G Abbott B Abdallah J Abdinov O Aben R Abolins M et al . Analysis of events with b-jets and a pair of leptons of the same charge in pp collisions at TeV with the ATLAS detector. J High Energy Phys. (2015) 10:150. 10.1007/JHEP10(2015)150

  • 549.

    Aad G Abbott B Abdallah J Abdel KS Abdinov O Aben R et al . Search for pair and single production of new heavy quarks that decay to a Z boson and a third-generation quark in pp collisions at TeV with the ATLAS detector. J High Energy Phys. (2014) 11:104. 10.1007/JHEP11(2014)104

  • 550.

    Aad G Abbott B Abdallah J Abdinov O Aben R Abolins M et al . Search for production of vector-like quark pairs and of four top quarks in the lepton-plus-jets final state in pp collisions at TeV with the ATLAS detector. J High Energy Phys. (2015) 8:105. 10.1007/JHEP08(2015)105

  • 551.

    Khachatryan V Besancon M Couderc F Dejardin M Denegri D Fabbro B et al . Search for vector-like charge 2/3 T quarks in proton-proton collisions at sqrt(s) = 8 TeV. Phys Rev. (2016) D93:012003. 10.1103/PhysRevD.93.012003

  • 552.

    Aaboud M Aad G Abbott B Abdinov O Abeloos B Abidi SH et al . Search for pair production of vector-like top quarks in events with one lepton, jets, and missing transverse momentum in TeV pp collisions with the ATLAS detector. J High Energy Phys. (2017) 8:052. 10.1007/JHEP08(2017)052

  • 553.

    Aaboud M Aad G Abbott B Abdallah J Abdinov O Abeloos B et al . Search for top squarks in final states with one isolated lepton, jets, and missing transverse momentum in TeV pp collisions with the ATLAS detector. Phys Rev. (2016) D94:052009. 10.1103/PhysRevD.94.052009

  • 554.

    The ATLAS collaboration . Search for new physics using events with b-jets and a pair of same charge leptons in 3.2 fb−1 of pp collisions at TeV with the ATLAS detector. ATLAS-CONF-2016-032 (2016).

  • 555.

    Sirunyan AM Tumasyan A Adam W Ambrogi F Asilar E Bergauer T et al . Search for pair production of vector-like T and B quarks in single-lepton final states using boosted jet substructure in proton-proton collisions at TeV. J High Energy Phys. (2017) 11:085. 10.1007/JHEP11(2017)085

  • 556.

    Chatrchyan S Besancon M Couderc F Dejardin M Denegri D Fabbro B et al . Inclusive search for a vector-like T quark with charge in pp collisions at = 8 TeV. Phys Lett. (2014) B729:14971. 10.1016/j.physletb.2014.01.006

  • 557.

    Aad G Abbott B Abdallah J Abdinov O Aben R Abolins M et al . Search for pair production of a new heavy quark that decays into a W boson and a light quark in pp collisions at TeV with the ATLAS detector. Phys Rev. (2015) D92:112007. 10.1103/PhysRevD.92.112007

  • 558.

    Khachatryan V Sirunyan AM Tumasyan A Adam W Asilar E Bergauer T et al . Search for pair-produced vectorlike B quarks in proton-proton collisions at =8 TeV. Phys Rev. (2016) D93:112009. 10.1103/PhysRevD.93.112009

  • 559.

    Aad G Abbott B Abdallah J Abdinov O Aben R Abolins M et al . Search for vector-like B quarks in events with one isolated lepton, missing transverse momentum and jets at 8 TeV with the ATLAS detector. Phys Rev. (2015) D91:112011.

  • 560.

    Altmannshofer W Bauer M Carena M . Exotic leptons: higgs, flavor and collider phenomenology. J High Energy Phys. (2014) 1:60. 10.1007/JHEP01(2014)060

  • 561.

    Falkowski A Straub DM Vicente A . Vector-like leptons: higgs decays and collider phenomenology. J High Energy Phys. (2014) 5:92. 10.1007/JHEP05(2014)092

  • 562.

    Dermisek R Hall JP Lunghi E Shin S . Limits on vectorlike leptons from searches for anomalous production of multi-lepton events. J High Energy Phys. (2014) 12:13. 10.1007/JHEP12(2014)013

  • 563.

    Kumar N Martin SP . Vectorlike leptons at the Large Hadron Collider. Phys Rev. (2015) D92:115018. 10.1103/PhysRevD.92.115018

  • 564.

    Hamada Y Kawai H Kawana K Tsumura K . Models of the LHC diphoton excesses valid up to the Planck scale. Phys Rev. (2016) D94:014007. 10.1103/PhysRevD.94.014007

  • 565.

    Pati JC Salam A . Lepton number as the fourth color. Phys Rev. (1974) D10:27589. 10.1103/PhysRevD.10.275

  • 566.

    Georgi H Glashow SL . Unity of all elementary particle forces. Phys Rev Lett. (1974) 32:43841. 10.1103/PhysRevLett.32.438

  • 567.

    Aaboud M Aad G Abbott B Abdallah J Abdinov O Abeloos B et al . Search for scalar leptoquarks in pp collisions at = 13 TeV with the ATLAS experiment. New J Phys. (2016) 18:093016. 10.1088/1367-2630/18/9/093016

  • 568.

    Aad G Abbott B Abdallah J Abdinov O Aben R Abolins M et al . Searches for scalar leptoquarks in pp collisions at = 8 TeV with the ATLAS detector. Eur Phys J. (2016) C76:5. 10.1140/epjc/s10052-015-3823-9

  • 569.

    Khachatryan V Sirunyan AM Tumasyan A Adam W Bergauer T Dragicevic M et al . CMS Collaboration. Search for pair-production of first generation scalar leptoquarks in pp collisions at with 2.6 fb−1. CMS-PAS-EXO-16-043 (2016).

  • 570.

    Khachatryan V Sirunyan AM Tumasyan A Adam W Bergauer T Dragicevic M et al . CMS Collaboration. Search for pair-production of second-generation scalar leptoquarks in pp collisions at with the CMS detector. CMS-PAS-EXO-16-007 (2016).

  • 571.

    CMS Collaboration . Search for the third-generation scalar leptoquarks and heavy right-handed neutrinos in ττhjj final states in pp collisions at 13 TeV. CMS-PAS-EXO-16-023 (2016).

  • 572.

    Aad G Abbott B Abdallah J Abdinov O Aben R Abolins M et al . Search for heavy long-lived multi-charged particles in pp collisions at TeV using the ATLAS detector. Eur Phys J. (2015) C75:362. 10.1140/epjc/s10052-015-3534-2

  • 573.

    Khoze VV Plascencia AD Sakurai K . Simplified models of dark matter with a long-lived co-annihilation partner. J High Energy Phys. (2017) 6:41. 10.1007/JHEP06(2017)041

  • 574.

    The ATLAS Collaboration . Search for electroweak production of supersymmetric particles in the two and three lepton final state at TeV with the ATLAS detector. ATLAS-CONF-2017-039 (2017).

  • 575.

    Ahriche A McDonald KL Nasri S Picek I . A critical analysis of one-loop neutrino mass models with minimal dark matter. Phys Lett. (2016) B757:399404. 10.1016/j.physletb.2016.04.022

  • 576.

    Aristizabal Sierra D Simoes C Wegman D . Closing in on minimal dark matter and radiative neutrino masses. J High Energy Phys. (2016) 6:108. 10.1007/JHEP06(2016)108

  • 577.

    Aad G Abajyan T Abbott B Abdallah J Khalek SA Abdinov O ., et al. Search for charginos nearly mass degenerate with the lightest neutralino based on a disappearing-track signature in pp collisions at TeV with the ATLAS detector. Phys Rev. (2013) D88:112006.

  • 578.

    Khachatryan V Sirunyan AM Tumasyan A Adam W Bergauer T Dragicevic M et al . Search for disappearing tracks in proton-proton collisions at TeV. J High Energy Phys. (2015) 1:96. 10.1007/JHEP01(2015)096

  • 579.

    The ATLAS collaboration . Search for Long-lived Charginos based on a Disappearing-Track Signature in pp Collisions at TeV with the ATLAS Detector. (2017) Report number: ATLAS-CONF-2017-017.

  • 580.

    Sakharov AD . Violation of CP invariance, c asymmetry, and baryon asymmetry of the universe. Pisma Zh Eksp Teor Fiz. (1967) 5:325. 10.1070/PU1991v034n05ABEH002497

  • 581.

    Kuzmin VA Rubakov VA Shaposhnikov ME . On the anomalous electroweak baryon number nonconservation in the early universe. Phys Lett. (1985) B155:36. 10.1016/0370-2693(85)91028-7

  • 582.

    Davidson S Nardi E Nir Y . Leptogenesis. Phys Rept. (2008) 466:10577. 10.1016/j.physrep.2008.06.002

  • 583.

    Ma E . Common origin of neutrino mass, dark matter, and baryogenesis. Mod Phys Lett. (2006) A21:177782. 10.1142/S0217732306021141

  • 584.

    Kashiwase S Suematsu D . Baryon number asymmetry and dark matter in the neutrino mass model with an inert doublet. Phys Rev. (2012) D86:053001. 10.1103/PhysRevD.86.053001

  • 585.

    Kashiwase S Suematsu D . Leptogenesis and dark matter detection in a TeV scale neutrino mass model with inverted mass hierarchy. Eur Phys J. (2013) C73:2484. 10.1140/epjc/s10052-013-2484-9

  • 586.

    Racker J . Mass bounds for baryogenesis from particle decays and the inert doublet model. J Cosmol Astropart Phys. (2014) 1403:025. 10.1088/1475-7516/2014/03/025

  • 587.

    Dev PSB Mohapatra RN . TeV scale model for baryon and lepton number violation and resonant baryogenesis. Phys Rev. (2015) D92:016007. 10.1103/PhysRevD.92.016007

  • 588.

    Hambye T Kannike K Ma E Raidal M . Emanations of dark matter: muon anomalous magnetic moment, radiative neutrino mass, and novel leptogenesis at the TeV scale. Phys Rev. (2007) D75:095003. 10.1103/PhysRevD.75.095003

  • 589.

    Babu KS Ma E . Singlet fermion dark matter and electroweak baryogenesis with radiative neutrino mass. Int J Mod Phys. (2008) A23:181319. 10.1142/S0217751X08040299

  • 590.

    Suematsu D . Leptogenesis in a neutrino mass model coupled with inflaton. Phys Lett. (2016) B760:53843. 10.1016/j.physletb.2016.07.048

  • 591.

    Bertolini S Di Luzio L Kolešová H Malinský M . Massive neutrinos and invisible axion minimally connected. Phys Rev. (2015) D91:055014.

  • 592.

    Aoki M Kanemura S Seto O . A model of TeV scale physics for neutrino mass, dark matter and baryon asymmetry and its phenomenology. Phys Rev. (2009) D80:033007. 10.1103/PhysRevD.80.033007

  • 593.

    Fischler W Giudice GF Leigh RG Paban S . Constraints on the baryogenesis scale from neutrino masses. Phys Lett. (1991) B258:458.

  • 594.

    Campbell BA Davidson S Ellis JR Olive KA . Cosmological baryon asymmetry constraints on extensions of the standard model. Phys Lett. (1991) B256:48490.

  • 595.

    Nelson AE Barr SM . Upper bound on baryogenesis scale from neutrino masses. Phys Lett. (1990) B246:1413.

  • 596.

    Harvey JA Turner MS . Cosmological baryon and lepton number in the presence of electroweak fermion number violation. Phys Rev. (1990) D42:334449.

  • 597.

    Antaramian A Hall LJ Rasin A . Hypercharge and the cosmological baryon asymmetry. Phys Rev. (1994) D49:38815.

  • 598.

    Frere JM Hambye T Vertongen G . Is leptogenesis falsifiable at LHC?J High Energy Phys. (2009) 1:51. 10.1088/1126-6708/2009/01/051

  • 599.

    Deppisch FF Harz J Hirsch M . Falsifying high-scale leptogenesis at the LHC. Phys Rev Lett. (2014) 112:221601. 10.1103/PhysRevLett.112.221601

  • 600.

    Deppisch FF Harz J Hirsch M Huang WC Päs H . Falsifying high-scale baryogenesis with neutrinoless double beta decay and lepton flavor violation. Phys Rev. (2015) D92:036005. 10.1103/PhysRevD.92.036005

  • 601.

    Boehm C Farzan Y Hambye T Palomares-Ruiz S Pascoli S . Is it possible to explain neutrino masses with scalar dark matter?Phys Rev. (2008) D77:043516. 10.1103/PhysRevD.77.043516

  • 602.

    El Aisati C Garcia-Cely C Hambye T Vanderheyden L . Prospects for discovering a neutrino line induced by dark matter annihilation. J Cosmol Astropart Phys. (2017) 1710:021. 10.1088/1475-7516/2017/10/021

  • 603.

    Ma E Ohata T Tsumura K . Majoron as the QCD axion in a radiative seesaw model. Phys Rev. (2017) D96:075039. 10.1103/PhysRevD.96.075039

  • 604.

    Babu KS Macesanu C . Two loop neutrino mass generation and its experimental consequences. Phys Rev. (2003) D67:073010. 10.1103/PhysRevD.67.073010

  • 605.

    Petcov ST . Remarks on the zee model of neutrino mixing (μ → , heavy neutrino → light neutrino gamma, etc.). Phys Lett. (1982) B115:4016.

  • 606.

    Bertolini S Santamaria A . The doublet majoron model and solar neutrino oscillations. Nucl Phys. (1988) B310:71442.

  • 607.

    Bertolini S Santamaria A . The strong CP problem and the solar neutrino puzzle: are they related?Nucl Phys. (1991) B357:22240.

  • 608.

    Smirnov AYu Tao Zj . Neutrinos with Zee mass matrix in vacuum and matter. Nucl Phys. (1994) B426:41533. 10.1016/0550-3213(94)90318-2

  • 609.

    Smirnov AYu Tanimoto M . Is zee model the model of neutrino masses? In: Neutrino '96. Proceedings, 17th International Conference on Neutrino Physics and Astrophysics, June 13-19, 1996Helsinki (1996). p. 30914.

  • 610.

    Frampton PH Glashow SL . Can the Zee ansatz for neutrino masses be correct?Phys Lett. (1999) B461:958. 10.1016/S0370-2693(99)00824-2

  • 611.

    Jarlskog C Matsuda M Skadhauge S Tanimoto M . Zee mass matrix and bimaximal neutrino mixing. Phys Lett. (1999) B449:240.

  • 612.

    Kanemura S Kasai T Lin GL Okada Y Tseng JJ Yuan CP . Phenomenology of Higgs bosons in the Zee model. Phys Rev. (2001) D64:053007. 10.1103/PhysRevD.64.053007

  • 613.

    Balaji KRS Grimus W Schwetz T . The Solar LMA neutrino oscillation solution in the Zee model. Phys Lett. (2001) B508:30110. 10.1016/S0370-2693(01)00532-9

  • 614.

    Koide Y . Prospect of the zee model. Nucl Phys Proc Suppl. (2002) 111:2946. 10.1016/S0920-5632(02)01726-7

  • 615.

    Brahmachari B Choubey S . Viability of bimaximal solution of the Zee mass matrix. Phys Lett. (2002) B531:99104. 10.1016/S0370-2693(02)01366-7

  • 616.

    Frampton PH Oh MC Yoshikawa T . Zee model confronts SNO data. Phys Rev. (2002) D65:073014. 10.1103/PhysRevD.65.073014

  • 617.

    Assamagan KA Deandrea A Delsart PA . Search for the lepton flavor violating decay A0/H0 → τ±μ at hadron colliders. Phys Rev. (2003) D67:035001. 10.1103/PhysRevD.67.035001

  • 618.

    Kanemura S Ota T Tsumura K . Lepton flavor violation in Higgs boson decays under the rare tau decay results. Phys Rev. (2006) D73:016006. 10.1103/PhysRevD.73.016006

  • 619.

    Aristizabal Sierra D Restrepo D . Leptonic charged higgs decays in the zee model. J High Energy Phys. (2006) 8:36. 10.1088/1126-6708/2006/08/036

  • 620.

    Aristizabal Sierra D Hirsch M . Experimental tests for the Babu-Zee two-loop model of Majorana neutrino masses. J High Energy Phys. (2006) 12:052. 10.1088/1126-6708/2006/12/052

  • 621.

    Ohlsson T Schwetz T Zhang H . Non-standard neutrino interactions in the Zee-Babu model. Phys Lett. (2009) B681:26975. 10.1016/j.physletb.2009.10.025

  • 622.

    Schmidt D Schwetz T Zhang H . Status of the Zee–Babu model for neutrino mass and possible tests at a like-sign linear collider. Nucl Phys. (2014) B885:52441. 10.1016/j.nuclphysb.2014.05.024

  • 623.

    Aristizabal Sierra D Dhen M Fong CS Vicente A . Dynamical flavor origin of ℤN symmetries. Phys Rev. (2015) D91:096004. 10.1103/PhysRevD.91.096004

  • 624.

    Deshpande NG Ma E . Pattern of symmetry breaking with two higgs doublets. Phys Rev. (1978) D18:2574.

  • 625.

    Lopez Honorez L Nezri E Oliver JF Tytgat MHG . The inert doublet model: an archetype for dark matter. J Cosmol Astropart Phys. (2007) 702:28. 10.1088/1475-7516/2007/02/028

  • 626.

    Borah D Gupta A . A new viable region of inert higgs doublet dark matter model with scotogenic extension. (2017) arXiv:1706.05034.

  • 627.

    Ho SY Tandean J . Probing scotogenic effects in higgs boson decays. Phys Rev. (2013) D87:095015. 10.1103/PhysRevD.87.095015

  • 628.

    Ho SY Tandean J . Probing scotogenic effects in e+e colliders. Phys Rev. (2014) D89:114025. 10.1103/PhysRevD.89.114025

  • 629.

    Hessler AG Ibarra A Molinaro E Vogl S . Probing the scotogenic FIMP at the LHC. J High Energy Phys. (2017) 1:100. 10.1007/JHEP01(2017)100

  • 630.

    Díaz MA Rojas N Urrutia-Quiroga S Valle JWF . Heavy higgs boson production at colliders in the singlet-triplet scotogenic dark matter model. J High Energy Phys. (2017) 8:017. 10.1007/JHEP08(2017)017

  • 631.

    Molinaro E Yaguna CE Zapata O . FIMP realization of the scotogenic model. J Cosmol Astropart Phys. (2014) 1407:15. 10.1088/1475-7516/2014/07/015

  • 632.

    Faisel G Ho SY Tandean J . Exploring X-ray lines as scotogenic signals. Phys Lett. (2014) B738:3805. 10.1016/j.physletb.2014.09.063

  • 633.

    Ho SY Faisel G Tandean J . Exploring X-ray lines as scotogenic signals. PoS (2016) DSU2015:68.

  • 634.

    Chowdhury TA Nasri S . The sommerfeld enhancement in the scotogenic model with large electroweak scalar multiplets. J Cosmol Astropart Phys. (2017) 1701:041. 10.1088/1475-7516/2017/01/041

  • 635.

    Merle A Platscher M . Parity problem of the scotogenic neutrino model. Phys Rev. (2015) D92:095002. 10.1103/PhysRevD.92.095002

  • 636.

    Merle A Platscher M . Running of radiative neutrino masses: the scotogenic model - revisited. J High Energy Phys. (2015) 11:148. 10.1007/JHEP11(2015)148

  • 637.

    Lindner M Platscher M Yaguna CE Merle A . Fermionic WIMPs and vacuum stability in the scotogenic model. Phys Rev. (2016) D94:115027. 10.1103/PhysRevD.94.115027

  • 638.

    Doršner I Fajfer S Greljo A Kamenik JF Košnik N . Physics of leptoquarks in precision experiments and at particle colliders. Phys Rept. (2016) 641:168. 10.1016/j.physrep.2016.06.001

  • 639.

    Aaij R Adeva B Adinolfi M Affolder A Ajaltouni Z Akar S et al . Test of lepton universality using B+K++ decays. Phys Rev Lett. (2014) 113:151601. 10.1103/PhysRevLett.113.151601

  • 640.

    Lees JP Poireau V Tisserand V Garra Tico J Grauges E Palano A et al . Evidence for an excess of decays. Phys Rev Lett. (2012) 109:101802. 10.1103/PhysRevLett.109.101802

  • 641.

    Lees JP Poireau V Tisserand V Grauges E Palano A Eigen G et al . Measurement of an excess of decays and implications for charged higgs bosons. Phys Rev. (2013) D88:072012. 10.1103/PhysRevD.88.072012

  • 642.

    Huschle M Kuhr T Heck M Goldenzweig P Abdesselam A Adachi I et al . Measurement of the branching ratio of relative to decays with hadronic tagging at Belle. Phys Rev. (2015) D92:072014. 10.1103/PhysRevD.92.072014

  • 643.

    Sato Y Iijima T Adamczyk K Aihara H Asner DM Atmacan H et al . Measurement of the branching ratio of relative to decays with a semileptonic tagging method. Phys Rev. (2016) D94:072007. 10.1103/PhysRevD.94.072007

  • 644.

    Hirose S Iijima T Adachi I Adamczyk K Aihara H Al Said S et al . Measurement of the τ lepton polarization and R(D*) in the decay . Phys Rev Lett. (2017) 118:211801. 10.1103/PhysRevLett.118.211801

  • 645.

    Aaij R Adeva B Adinolfi M Affolder A Ajaltouni Z Akar S et al . Measurement of the ratio of branching fractions . Phys Rev Lett. (2015) 115:111803. 10.1103/PhysRevLett.115.111803

  • 646.

    Barbier R Berat C Besancon M Chemtob M Deandrea A Dudas E et al . R-parity violating supersymmetry. Phys Rept. (2005) 420:1202. 10.1016/j.physrep.2005.08.006

  • 647.

    Ibanez LE Ross GG . Discrete gauge symmetries and the origin of baryon and lepton number conservation in supersymmetric versions of the standard model. Nucl Phys. (1992) B368:337. 10.1016/0550-3213(92)90195-H

  • 648.

    Hirsch M Valle JWF . Supersymmetric origin of neutrino mass. New J Phys. (2004) 6:76. 10.1088/1367-2630/6/1/076

  • 649.

    Hirsch M Diaz MA Porod W Romao JC Valle JWF . Neutrino masses and mixings from supersymmetry with bilinear R parity violation: a Theory for solar and atmospheric neutrino oscillations. Phys Rev. (2000) D62:113008. 10.1103/PhysRevD.62.113008

  • 650.

    Diaz MA Hirsch M Porod W Romao JC Valle JWF . Solar neutrino masses and mixing from bilinear R parity broken supersymmetry: analytical versus numerical results. Phys Rev. (2003) D68:013009. 10.1103/PhysRevD.71.059904

  • 651.

    Grossman Y Rakshit S . Neutrino masses in R-parity violating supersymmetric models. Phys Rev. (2004) D69:093002. 10.1103/PhysRevD.69.093002

  • 652.

    de Campos F Eboli OJP Magro MB Porod W Restrepo D Das SP et al . Probing neutralino properties in minimal supergravity with bilinear R-parity violation. Phys Rev. (2012) D86:075001. 10.1103/PhysRevD.86.075001

  • 653.

    Mukhopadhyaya B Roy S Vissani F . Correlation between neutrino oscillations and collider signals of supersymmetry in an R-parity violating model. Phys Lett. (1998) B443:1915.

  • 654.

    Choi SY Chun EJ Kang SK Lee JS . Neutrino oscillations and R-parity violating collider signals. Phys Rev. (1999) D60:075002.

  • 655.

    Romao JC Diaz MA Hirsch M Porod W Valle JWF . A Supersymmetric solution to the solar and atmospheric neutrino problems. Phys Rev. (2000) D61:071703. 10.1103/PhysRevD.61.071703

  • 656.

    Porod W Hirsch M Romao J Valle JWF . Testing neutrino mixing at future collider experiments. Phys Rev. (2001) D63:115004. 10.1103/PhysRevD.63.115004

  • 657.

    Hirsch M Vicente A Porod W . Spontaneous R-parity violation: lightest neutralino decays and neutrino mixing angles at future colliders. Phys Rev. (2008) D77:075005. 10.1103/PhysRevD.77.075005

  • 658.

    Bartl A Hirsch M Vicente A Liebler S Porod W . LHC phenomenology of the μνSSM. J High Energy Phys. (2009) 5:120. 10.1088/1126-6708/2009/05/120

  • 659.

    Liebler S Porod W . On-shell renormalization of neutralino and chargino mass matrices in R-parity violating models - Correlation between LSP decays and neutrino mixing angles revisited. Nucl Phys. (2012) B855:774800. 10.1016/j.nuclphysb.2011.10.027

  • 660.

    Deshpande NG He XG . Consequences of R-parity violating interactions for anomalies in and b+μ. Eur Phys J. (2017) C77:134. 10.1140/epjc/s10052-017-4707-y

  • 661.

    Dreiner HK Staub F Vicente A Porod W . General MSSM signatures at the LHC with and without R-parity. Phys Rev. (2012) D86:035021. 10.1103/PhysRevD.86.035021

  • 662.

    Bennett GW Bousquet B Brown HN Bunce G Carey RM Cushman P et al . Final report of the muon E821 anomalous magnetic moment measurement at BNL. Phys Rev. (2006) D73:072003. 10.1103/PhysRevD.73.072003

  • 663.

    Ellis J . TikZ-Feynman: Feynman diagrams with TikZ. Comput Phys Commun. (2017) 210:10323. 10.1016/j.cpc.2016.08.019

Appendix

On the relative contribution of operators

Oftentimes the effective ΔL = 2 operators are discussed using a cutoff regularization scheme. In the following, however, we outline the relative contribution of the different ΔL = 2 operators to neutrino mass using dimensional regularization with a momentum-independent renormalization scheme such as renormalization. Power counting in the SM effective theory establishes that the dominant contributions to neutrino mass are given by (i) the lowest-dimensional Weinberg-like operator which is induced via matching at the new physics scale Λ and (ii) the contributions induced by mixing via renormalization group running of the operator into the Weinberg operator or other lower-dimensional Weinberg-like operators.

Using naive dimensional analysis we discuss in more detail the relative contribution to neutrino mass from each operator in the SM effective field theory. Note that here we follow the matching and running from low energy scale to high energy scale. Below the electroweak scale effective operators that can contribute to neutrino masses should contain two neutrinos and possibly additional fields. Those additional fields have to be closed off and their contribution to neutrino masses vary: for photons and gluons, the contribution from the tadpole diagram vanishes; for fermions f, the contribution is proportional to a factor m3/16π2Λ3 per fermion loop. Thus, the contribution of operators with additional fields to neutrino mass either vanishes or is generally suppressed. Matching at the electroweak scale may similarly include loops with electroweak gauge bosons or the top quark and lead to a suppression of the respective operator. Additional Higgs fields yield a factor v/Λ each. Above the electroweak scale the operators generally mix. Higher-dimensional operators also mix into lower dimensional ones. For example although the operator mixes into the operator O1 via renormalization group running and thus it is an operator of lower dimension, its contribution to the Wilson coefficient is suppressed by a factor of order and therefore it is of the same order as the operator . At the new physics scale the relative size of the Wilson coefficients is determined by the couplings and the loop level at which they are generated. The Wilson coefficient of the Weinberg-like operators at the new physics scale may be suppressed by a loop factor compared to other operators, but the other operators receive a further loop-factor suppression when matching onto the effective interactions at the electroweak scale or finally onto the neutrino mass term at a lower scale. The contributions of all operators to neutrino mass has at least the same loop-factor suppression as the leading Weinberg-like operator which is induced by matching at the new physics scale. Higher-dimensional Weinberg-like operators will induce the lower-dimensional ones via mixing when running the Wilson coefficients to the low scale, but the contribution of the induced operator is still of the same order as the original higher-dimensional operator. In summary, an order of magnitude estimate of neutrino mass can be obtained from the leading Weinberg-like operator which is induced from matching at the new physics scale keeping in mind that its contribution to lower-dimensional Weinberg-like operators will be of a similar order of magnitude.

Summary

Keywords

neutrino masses, lepton flavor violation, lepton number violation, beyond the standard model, effective field theory, model building, LHC, dark matter

Citation

Cai Y, Herrero García J, Schmidt MA, Vicente A and Volkas RR (2017) From the Trees to the Forest: A Review of Radiative Neutrino Mass Models. Front. Phys. 5:63. doi: 10.3389/fphy.2017.00063

Received

18 July 2017

Accepted

14 November 2017

Published

04 December 2017

Volume

5 - 2017

Edited by

Frank Franz Deppisch, University College London, United Kingdom

Reviewed by

Mayumi Aoki, Kanazawa University, Japan; Koji Tsumura, Kyoto University, Japan; Hiroaki Sugiyama, Toyama Prefectural University, Japan

Updates

Copyright

*Correspondence: Juan Herrero García

This article was submitted to High-Energy and Astroparticle Physics, a section of the journal Frontiers in Physics

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics