Lepton flavor violation and scalar dark matter in a radiative model of neutrino masses

We consider a simple extension of the Standard Model that can account for the dark matter and explain the existence of neutrino masses. The model includes a vector-like doublet of SU(2), a singlet fermion, and two scalar singlets, all of them odd under a new Z2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document} symmetry. Neutrino masses are generated radiatively by one-loop processes involving the new fields, while the dark matter candidate is the lightest neutral particle among them. We focus specifically on the case where the dark matter particle is one of the scalars and its relic density is determined by its Yukawa interactions. The phenomenology of this setup, including neutrino masses, dark matter and lepton flavor violation, is analyzed in some detail. We find that the dark matter mass must be below 600 GeV to satisfy the relic density constraint. Lepton flavor violating processes are shown to provide the most promising way to test this scenario. Future μ→3e\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu \rightarrow 3e$$\end{document} and μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu $$\end{document}–e conversion experiments, in particular, have the potential to probe the entire viable parameter space of this model.


Introduction
Neutrino masses and dark matter provide compelling evidence for physics beyond the Standard Model (SM). The gravitational effects of dark matter have been observed, for instance, in galaxies, clusters of galaxies, the large scale structure of the Universe, and in the cosmic microwave background radiation. Recently, the WMAP [1] and Planck [2] collaborations have determined the current dark matter density in the Universe to an unprecedented precision: DM h 2 = 0.1186 ± 0.0031. Oscillations experiments [3][4][5], on the other hand, have demonstrated that neutrinos have non-zero masses and have allowed us to measure the mixing and mass parameters in the neutrino sector [6]. a e-mail: cyaguna@gmail.com Within the SM, neither neutrino masses nor dark matter can be explained, and current data does not tell us how the SM should be extended to account for them. A particularly appealing possibility is that these two problems are not independent, as usually assumed, but arise from the same type of new physics. Moreover, this new physics may appear, as suggested by the WIMP paradigm of dark matter, at the TeV scale -the scale that is currently being probed by the LHC. In [7], many models of this type, featuring neutrino masses at one-loop, were found and classified. They all contain a small number of additional fields and a new discrete symmetry to stabilize the dark matter particle. In this paper, we study one of these models, denoted as T1-3 in [7].
In this model, which had been previously considered in [8,9], the SM is extended with a vector-like SU(2) doublet, a singlet fermion, and two singlet scalars. We focus specifically on the case where the dark matter particle is a scalar and its relic density is determined by the Yukawa interactions (rather than the scalar ones), a possibility that to our knowledge has not been studied before. As we will show, in this setup strong correlations arise between dark matter, neutrino masses and lepton flavor violating (LFV) processes. In our analysis, we will investigate these correlations in some detail, partly relying on a scan over the parameter space of this model. First, we use a large sample of viable models to determine and study the regions that are consistent with the constraints from neutrino masses, μ → eγ , and dark matter. Then, we review the predictions for other LFV processes and the prospects for detection in future experiments. We find, in particular, that future searches for μ → 3e and μ-e conversion in nuclei have the potential to probe the entire parameter space of this model. We also discuss alternative ways to test this scenario, including collider searches as well as direct and indirect dark matter searches.
The rest of the paper is organized as follows. In the next section we describe the model and introduce our notation. Then, in Sect. 3 we analyze semi-quantitatively the most relevant phenomenological aspects of this model, including neutrino masses, LFV processes, and dark matter. Our main results, based on a scan over the parameter space of this model, are presented in Sect. 4. First we examine the viable parameter space and then we demonstrate that the rates for several LFV processes are significant, offering the potential to probe this model thoroughly in the near future. In Sect. 5 we briefly discuss the possibility to test this scenario using collider or dark matter searches. Finally, we summarize and draw our conclusions in Sect. 6.

Description of the model
The model we consider belongs, in the generic classification of [10], to the T1-iii models of one-loop neutrino masses. It corresponds, in particular, to the model T1-3-A with α = 0 (where α is related to the hypercharge), as defined in [7]. In this model the SM particle content is extended with a vector-like fermion (or two chiral fermion) doublet under SU(2), D, D , one left-handed singlet fermion, S, and two real scalar singlets φ i (i = 1, 2). To guarantee the stability of the dark matter candidate and to prevent tree-level neutrino masses, the SM gauge group is extended with a Z 2 discrete symmetry, under which all the new fields are odd, while the SM fields are even. Thus, the charges of the new particles under the SU(2) ⊗ U(1)⊗ Z 2 symmetry are given by for the fermion doublets and for the fermion and scalar singlets, respectively. The SM electroweak sector under the same symmetry transforms instead as for the doublets (i = 1 · · · 3) and for the SU(2) singlets. The most general Lagrangian consistent with the symmetry and the particle content of our model can be written as where L SM is the SM Lagrangian, V (H, φ i ) is the scalar potential, andH = iσ 2 H * . This Lagrangian includes a Dirac mass term for the vector-like doublet fermion, a Majorana mass term for the singlet fermion, Yukawa terms between the new singlet and doublet fermions, and a new Yukawa interaction for the SM lepton doublets. It is precisely this new interaction that allows to obtain non-zero neutrino masses. At tree-level, neutrinos remain massless because the scalar fields φ i do not acquire a vev -the Z 2 symmetry must remain exact to explain the dark matter. It is only at one-loop that neutrinos will get a Majorana mass. In fact, it is not difficult to verify that this Lagrangian generically violates lepton number. Being real scalar fields, φ i should have zero lepton number and likewise the singlet Majorana fermion S. But then the terms with coefficients α i j , β 1 and β 2 cannot simultaneously conserve lepton number. Lepton number would still be a good symmetry in the limits α i j = 0 or β 1 = β 2 = 0. Neutrino masses, which are lepton-number violating, must therefore vanish in such cases.
Once expanded, the Lagrangian includes the following mass terms and the following couplings with the SM particles Defining := (ψ, ψ , S), the Majorana mass matrix for the odd neutral fermions can be written, from Eq. (6), as This matrix can be diagonalized via the transformation where ξ is the mixing matrix and χ i are the physical mass eigenstates with masses m i such that m i < m j for i < j. In addition to these three Majorana fermions, the spectrum contains also a charged Dirac fermion with mass μ ≡ m E -see Eq. (6) -constrained by collider searches [11] to be larger than about 103.5 GeV, and two neutral scalars with masses m φ 1,2 .

Phenomenology
In this section, the main phenomenological features of this model are discussed, in particular regarding neutrino masses, dark matter, and lepton flavor violating processes.

Neutrino masses
Radiative mechanisms are a very attractive way of explaining neutrino masses and have been implemented within many different scenarios -see e.g. [9,[12][13][14][15][16][17][18][19][20][21][22][23]. In the model we consider, neutrino masses are obtained through a radiative mechanism at one-loop. The relevant diagram, shown in the left panel of Fig. 1, contains the neutral odd particles (scalars and fermions) running in the loop. Since one of these neutral particles will be the dark matter candidate, some connection between dark matter and neutrino masses is necessarily present. The evaluation of the diagram yields, for a given neutral fermion (k) and a given singlet scalar (r ), the following contribution to the neutrino mass matrix where r lm = α rl ξ 2m is the coupling between φ r , χ m , and ν c L . In this equation, ∝ 1 ( = 4 − d) whereas μ R is an arbitrary renormalization constant. The neutrino mass matrix element m ν i j can therefore be computed as a sum over all the fermions and scalars that run in the loop: where we have used the fact that Notice that, as expected, the final result is finite and thus independent of the unphysical renormalization factor μ R . As anticipated, the neutrino mass matrix vanishes in the limit α i j = 0, where lepton number is conserved. To see that it also vanishes in the limit β 1 = β 2 = 0, we can diagonalize Eq. (8) analytically at leading order in β 1 , β 2 1 to obtain Hence, the neutrino masses in Eq. (11) can be expressed as which indeed goes to zero in the limit β 1,2 → 0. Since the experimentally determined neutrino masses are tiny, either α i j or β 1,2 must be suppressed. To estimate the typical values of the parameters that give rise to neutrino masses consistent with the data, we can evaluate Eq. (14) for masses around the TeV scale to obtain Notice that neutrino masses cannot distinguish whether it is the α i j or the β 1,2 that are suppressed but, as we will see later, the dark matter constraint only allows the latter option.
Since α i j and β 1,2 enter both quadratically in m ν , the required coupling suppression is not as large as in the wellknown scotogenic model [16], where λ 5 can be of order 10 −10 [24]. In our model the coupling suppression is only of order 10 −5 .
It is easy to verify that, in our scenario, the neutrino mass matrix has only two non-zero eigenvalues. This result is a direct consequence of the two real scalar fields that interact with the lepton doublets. To obtain three non-zero neutrino masses, one more scalar field would be needed. But, since current neutrino data is consistent with one massless neutrino, we will stick to this minimal framework, as outlined in the previous section. The neutrino spectrum, therefore, is necessarily hierarchical and we assume it henceforth to be of normal type (NH).
The resulting neutrino mass matrix, Eq. (11), can be written in a way reminiscent of the seesaw mechanism as with This form of the neutrino mass matrix allows us to use a slightly modified version of the Casas-Ibarra parametrization [25]. Explicitly, we have that where For definiteness we take all the parameters entering into Eq. (18) to be real. That is, the Dirac and Majorana phases present in U as well as the phase of θ are assumed to be equal to zero. It is important to stress that these assumptions are not expected to affect our results in a significant way. Regarding lepton flavor violating processes, it is known that CP phases can modify only the angular distribution of the final state particles [27], but not the rates of these processes, which are the quantities we impose constraints on. Regarding the relic density, the CP phases in this model can manifest only via loop diagrams, and are therefore expected to have a negligible effect on the calculation of the dark matter density, which is the other constraint we rely on.
Equation (18) provides, for a given set of scalar and fermion mass parameters, the general form of the α i j couplings that is consistent with the observed neutrino data, which we use as input parameters. In other words, the constraints on neutrino masses and mixing angles are automatically incorporated into the structure of α as given by Eq. (18), simplifying enormously the analysis of the viable parameter space.

Lepton flavor violating processes
The existence of non-zero neutrino masses and mixing angles also imply the violation of the lepton flavor. In consequence, processes involving the charged leptons and where the lepton flavor is not conserved -lepton flavor violating (LFV) processes -such as μ → eγ and μ → 3e should occur at some level in this model. The relevant diagrams appear first at one-loop, as illustrated in the right panel of Fig. 1 for μ → eγ . Notice that the fermion running in the loop is now the Dirac charged fermion rather than the neutral ones.
In this model, the μ → eγ branching ratio can be evaluated to be where An analogous expression holds for the related processes τ → μγ and τ → eγ . For other processes, the analytical expressions are more complicated and are therefore better studied numerically (with the help of FlavorKit [28]), as we will do in the next section. Typically, μ → eγ provides the most stringent constraint among LFV processes, so many works have focused on this process. But, as emphasized in [24] this situation might drastically change in the near future thanks to the significant improvements that will likely be achieved for other LFV processes, as shown in Table 1. Particularly relevant will be the  [37,38] μ − , Au → e − , Au 7 × 10 −13 [39] μ − , Al → e − , Al 10 −15 − 10 −18 [40] μ − , SiC → e − , SiC 10 −14 [41] expected limits on μ → 3e and μ-e conversion in nuclei, which may improve, respectively, by up to four and six orders of magnitude [42]. It is important, therefore, not to limit the discussion on LFV processes to μ → eγ but to consider also these other promising processes when analyzing future detection prospects. A crucial difference between neutrino masses and LFV processes, both of which appear at one-loop in this model, is that the former violate lepton number, while the latter do not. As a result, LFV processes depend only on the α i j couplings but not on β 1,2 or m S , as shown explicitly above for μ → eγ .
To obtain observable rates for LFV processes, two conditions must generically be satisfied: the particles in the loop must no be that heavy and the relevant couplings, α i j in this case, should be of order 1-0.1. As we saw in the previous subsection, this second requirement is compatible with neutrino masses as long as the beta parameters are tiny, β 1,2 1. Next, we show that, under certain circumstances, these requirements are actually enforced by the dark matter constraint.

Dark matter phenomenology
The dark matter candidate in this model is the lightest neutral odd (under the Z 2 ) particle in the spectrum, which can be one of the three fermions, or one of the two scalars. Since the fermion and scalar masses are free parameters, both of these possibilities can be realized in this model. When the dark matter is a fermion, the resulting phenomenology is very similar to that of the so-called singlet-doublet fermion model [43][44][45], which has been extensively studied in the recent literature -see e.g. [46][47][48][49][50]. It has been shown, in particular, that by adjusting the mixing between the singlet and the doublet one can obtain the correct relic density, via freeze-out in the early Universe, for dark matter masses below 1 TeV or so. The additional scalars present in this model may slightly modify this picture due to coannihilation effects, as recently demonstrated in [9].
When the dark matter is the scalar, two different scenarios can be distinguished depending on its dominant interactions. Scalars, in fact, not only have the Yukawa interactions explicitly shown in Eq. (7) but also scalar interactions with the SM Higgs boson, implicitly included in V (H, φ i ). These scalar interactions give rise to the so-called singlet scalar or Higgs-portal model [51][52][53], where dark matter annihilations are mediated by the Higgs boson. These models are quite predictive and have been the subject of many previous analyses [54][55][56][57][58][59][60][61][62][63]. In this paper we assume instead that the scalar dark matter interacts dominantly via the Yukawa interactions in Eq. (7). Thus, it will annihilate into leptons via t-channel fermion-mediated diagrams with a cross section proportional to α 4 i j -see Fig. 2. Consequently, non-trivial correlations between neutrino masses, LFV processes and Fig. 2 General diagram contributing to the relic density computation dark matter are expected. In fact, we can already state that to obtain a relic density in agreement with the observations, α i j should be of order one and the mediators cannot be that heavy, which are essentially the same conditions that ensure observable rates for LFV processes, as we saw in the previous subsection. Throughout the rest of the paper, we will be working on this specific framework where the dark matter particle is a scalar, denoted by φ 1 , that interacts dominantly via the Yukawa terms with the SM lepton doublets. Notice that as a result, the dark matter turns out to be leptophilic, with important implications for the direct detection prospects, as will be discussed in Sect. 5. In this setup, dark matter annihilates into charged leptons with a cross section given, in the non-relativistic limit, by [64] σ Our expression accounts for the possibility of annihilation into lepton pairs that are not flavor-diagonal. In our numerical analysis we will see that, due to the constraints from μ → eγ , it is the τ + τ − final state that actually dominates the cross section. The cross section gets suppressed as the dark matter mass or the charged fermion mass increases. Notice also that it has a strong velocity dependence (∝ v 4 ). During freeze-out, this velocity suppression is not that important but today, when v ∼ 10 −3 in our Galaxy, it pretty much prevents the annihilation of dark matter particles into two charged leptons, significantly affecting the dark matter indirect detection signatures in this model. For neutrinos, the corresponding annihilation cross section receives contributions from t-and u-channel exchanges of all three neutral fermions χ l or χ m , The leading non-relativistic contribution of O(v 2 ) is generated by mixing in the neutral fermion sector. As a result of this mixing, all odd neutral fermions are Majorana particles, and t-and u-channel diagrams interfere destructively.
In the limit of negligible mixing, the terms proportional to v 2 vanish, and one is left with O(v 4 ) terms similar to Eq. (22) pertinent to Dirac fermions.
To obtain an accurate prediction of the dark matter relic density, we implemented this model into DarkSUSY [65] (micrOmegas [66] gives numerical errors in the evaluation of the annihilation cross section). Specifically, we use the σ v above to compute the invariant annihilation rate as defined in the DarkSUSY manual, and use the function dsrdens() to obtain the relic density. We also verified that the computations done in this way match the approximate analytical results known in the literature.

Numerical results
In this model non-trivial correlations between dark matter, neutrino masses and lepton flavor violating processes are expected. The reason is that, on the one hand, the dark matter constraint requires sizable α i j couplings between the dark matter and the SM leptons. On the other hand, these couplings cannot be flavor diagonal because they determine the structure of the neutrino mass matrix. Consequently, these couplings induce significant rates for LFV processes. In this section, we use a scan over the parameter space of this model to numerically study these correlations. To that end, we first obtain a large sample of viable models that we use to analyze the regions in the parameter space of the model that are consistent with current bounds. Then, we analyze the prediction for LFV processes and demonstrate that future LFV experiments have the potential to probe most of the viable models.

The viable parameter space
The relevant free parameters of this model are just seven: m φ 1,2 , μ, m S , β 1,2 and θ . Our sample of viable models is obtained after scanning over the allowed range of these parameters and imposing the following constraints: the dark matter relic density ( DM h 2 ∼ 0.12), the μ → eγ upper limit -BR(μ → eγ ) < 5.7 × 10 −13 -, and perturbativity (somewhat arbitrarily, we require all dimensionless couplings to be smaller than three). Regarding the dark matter relic density, we assume that it is obtained via thermal freezeout in the early Universe and that coannihilation effects play no role in its determination, as is generically the case. To enforce this latter condition, we require all other masses to be larger than 1.2m φ 1 . The upper bound on the mass parameters is taken to be 10 TeV. By construction, all our models are Fig. 3 The viable models projected onto the plane (m φ ,α i1 ). The α i1 couplings determine the dark matter annihilation rate and must therefore be non-negligible consistent with neutrino data, for we use as input the experimental data on neutrino masses and mixing angles according to Eq. (18). In our analysis, we made use of SARAH and the BSM Tool Box scripts [67], which facilitate the whole process.
Let us begin by examining the viable regions of the parameter space. Figure 3 displays our set of viable models in the plane (m φ ,α i1 ), with i = e, μ, τ . Notice that, as a result of the strong bounds from μ → eγ , α τ 1 α e1 , α μ1 . Thus, dark matter annihilates dominantly into the τ + τ − final state. In fact, α τ 1 increases with the dark matter mass, going from about one for m φ ∼ 50 GeV to three -the perturbativity limit we imposed -for m φ ∼ 600 GeV. Consequently, the dark matter particle in this scenario must be light, lying below 600 GeV.
The condition α 1τ α 1e , α 1μ is not satisfied for generic values of the α i j couplings. According to Eq. (18), we have that If we now require α 1τ α 1e , α 1μ we get that, for typical values of the neutrino parameters, θ ∼ 0.35. The left panel of Fig. 4 shows the value of θ for our set of viable models. We see that, indeed, θ varies only within a narrow range around  (m φ ,θ). Notice that the experimental constraints essentially select a single value of θ. Right: the viable models projected onto the plane (β 1 , β 2 ). Neutrino masses require β 1,2 to be small  (m φ ,m χ1 ). The relic density constraint favors a χ 1 not much heavier than the dark matter (φ). Right: the viable models projected onto the plane (m φ ,m χ3 ). Notice that the heaviest odd fermion, χ 3 , rarely has a mass above 1 TeV 0.35. Thus, the dark matter constraint and the limits on LFV processes select a rather specific value for θ .
As we saw in the previous section, the parameters β 1 and β 2 must be tiny so that neutrino masses are sufficiently suppressed. The right panel of Fig. 4 displays the viable models in the plane (β 1 , β 2 ) and confirms that this is really the case. These parameters are never larger than about 10 −4 . β 1 could be much smaller than that, whereas β 2 varies only between 10 −4 and about 10 −6 . Notice that the dark matter constraint plays also a role in this case, as it enforces a sizable value of the α i j couplings.
The relic density depends strongly on the masses of the dark matter particle and of the Dirac fermion that mediates the annihilation processes. The left panel of Fig. 5 displays the ratio m E /m φ 1 as a function of the dark matter mass for our sample of viable points. The lower limit on m E /m φ 1 is 1.2 (to prevent coannihilations) and we see from the figure that it tends to that value at the highest dark matter masses.
In any case, m E is never much larger than the dark matter mass. The dark matter constraint thus requires the existence of a charged fermion with a mass below 720 GeV.
Two of the three neutral fermions have a mass very close to m E ∼ μ as a result of the small mixing induced by the parameters β 1,2 1 (the mixing among the neutral fermions is very small). The other neutral fermion, with mass ∼ m S , can have a much larger mass. The right panel of Fig. 5 displays the dark matter mass versus the mass of the heaviest neutral fermion, M χ 3 . Even though this mass can reach its upper limit (10 TeV), in most models it tends to be close to the dark matter mass. The spectrum in this setup is thus rather compressed.

Implications for LFV processes
In this section we study in detail the predicted rates for the most relevant LFV processes and compare them against cur- Fig. 6 The branching ratios τ → μγ (blue squares) and τ → eγ (red stars) as a function of the dark matter mass for our sample of viable models. The solid lines show the current bound for each process while the dashed line corresponds to the expected future sensitivity Fig. 7 The rates of the most relevant μ-e LFV processes as a function of the dark matter mass for our sample of viable models rent limits and expected sensitivities in planned experiments. As we will see, this scenario predicts sizable rates for several LFV processes, particularly μ → 3e and CR(μ − e). These processes, in fact, provide the most promising way of probing this model in the near future.
Let us begin our analysis with LFV τ decays. Figure 6 displays the branching ratios for the decay processes τ → μγ (blue squares) and τ → eγ (red stars) as a function of the dark matter mass. Both processes feature the same behaviour -its maximum value decreasing with the dark matter mass -and similar values at a given mass. These branching ratios vary over a wide range at low masses (from 10 −5 to 10 −13 or so) but tend to concentrate between 10 −8 and 10 −10 at higher dark matter masses. For comparison, current limits are also shown in the figure as solid lines, as well as the expected future sensitivity for both processes (dashed line). Interestingly, we see that current limits can be violated for dark matter masses below 250 GeV or so. In addition, future experiments will be able to probe a significant fraction of models over the entire range of dark matter masses.
Let us now switch to μ-e LFV processes, which are typically more relevant. Figure 7 shows, as a function of the dark matter mass, the μ → eγ (blue squares) and μ → 3e (red stars) branching ratios as well as the μ-e conversion rate in Titanium (green circles). Notice that whereas the current limit on BR(μ → 3e) can be violated only at low masses, the one on CR(μ-e, Ti) can be exceeded over the entire viable range of dark matter masses. It is also important to stress that the rates of all these μ-e LFV processes do not extend to arbitrarily low values, and tend instead to lie up to few orders of magnitude below present bounds. Figure 8 explores the correlations between the different μ-e LFV processes in this model. The left panel shows the viable points in the plane BR(μ → eγ ) vs. BR(μ → 3e), whereas the right panel shows BR(μ → eγ ) vs CR(μ-e, Ti).  We also displayed, for each process, its current limit (solid lines) and its expected future sensitivity in planned experiments (dashed lines). As we have already observed in the previous figure, current limits on μ → 3e and μ-e conversion in nuclei are not necessarily fulfilled and can be violated by more than two orders of magnitude. Thus, satisfying the μ → eγ limit does not guarantee compatibility with other μ-e LFV experiments. In other words, contrary to the naive expectation, μ → 3e and μ-e in nuclei may give a more stringent constraint than μ → eγ for some models. Regarding future prospects, we see from the figure that the improvement in μ → eγ though important will not be decisive as many models feature much smaller branching ratios. Completely different is the situation for μ → 3e and μ-e in nuclei. There, the improvements will be more significant, and their impact will be crucial for this scenario. In fact, very few models lie beyond the expected sensitivity of future μ → 3e or μ-e conversion experiments.
To emphasize this point, we plot instead BR(μ → 3e) versus CR(μ-e, Ti) in Fig. 9. As noticed in the previous figure, only few points lie beyond the expected future sensitivity for either process. But now we can also observe that not a single point in our scan lies beyond the expected sensitivity for both processes. Future experiments searching for μ → 3e and μ-e conversion in nuclei therefore have the potential to probe most, if not all, of the viable parameter space of this model.

Discussion
As we have seen, this model could be probed in future experiments via LFV processes. Other ways to test it include the direct [68] and indirect detection [69] of dark matter as well as searches at the LHC. In this section we give a brief look at these alternatives.
In our setup, the direct detection prospects are not very encouraging because the dark matter particle is leptophilic. Thus, it has a vanishing tree-level scattering cross section with nuclei. Radiative corrections will generate a non-zero cross section but it will still be highly suppressed. At oneloop, both the spin-independent interaction (mediated by the Higgs) and the anapole-moment (mediated by the photon) can be generated in this model. Recently, these one-loop effects were calculated explicitly within a similar radiative model [70] and it was found that the predicted signal lies well below the expected sensitivity of future experiments for the anapole-moment and just barely within it for the spin-independent one. In our scenario the situation is worse because the spin-independent cross section will be further suppressed by lepton masses or the β 1,2 parameters. We can safely conclude, therefore, that direct detection experiments cannot probe the viable parameter space of this model.
In spite of the velocity-suppressed annihilation cross section, the indirect detection of dark matter in this model is not entirely out of question. Notice that the velocity suppression in this setup is stronger than for Majorana fermions, featuring a v 4 rather than a v 2 behavior. Given that v ∼ 10 −3 for dark matter particles in the galactic halo, this v 4 suppression implies an annihilation rate today about twelve orders of magnitude below the thermal one (σ v ∼ 3 × 10 −26 cm 3 s −1 . This suppression can however be avoided with the emission of an additional photon, the so-called internal bremsstrahlung process [71], which gives rise to a gamma-ray feature at high energy. In [64,72], where this effect was studied for scalar dark matter, it was found that, depending on the specific parameters of the model, this feature could actually be observed in current and planned gamma-ray telescopes. The other indirect detection channels -neutrinos, positrons and antiprotons -are not expected to play any role in constraining or testing this model.
As we have seen, our model predicts the existence of several particles with masses below the TeV scale. In particular, the charged fermion should have a mass below 720 GeV while the dark matter particle must be lighter than about 600 GeV. Naively, one would think that such particles should be easily produced and detected at the LHC, but that is not really the case, as we now explain.
The LHC collaborations ATLAS and CMS have so far not performed any specific searches for models with radiative neutrino masses and dark matter candidates. Possible LHC constraints can therefore only be derived from analyses with leptonic signatures that are either model-independent or have been performed for models with similar additional physical states.
At first sight, our model bears some similarity with the Minimal Supersymmetric Standard Model (MSSM), where the heavy charged and neutral Z 2 -odd fermions are the charginosχ ± 1 and the (often mass-degenerate) neutralinos χ 0 2 . However, contrary to the MSSM our fermions do not decay via cascades to a "golden" charged three-lepton final state, a neutrino and two light neutralinosχ 0 1 carrying away missing transverse energy E T , but directly to two charged leptons, a charged lepton and a neutrino or two neutrinos plus two invisible scalar dark matter particles. The trilepton ATLAS [73] and CMS [74] analyses do therefore not apply.
Two-lepton final states with E T have also been analysed by ATLAS [75] and CMS [76], but under specific assumptions for the masses of the intermediate sleptons. While these assumptions are absent in the dilepton searches for Z bosons [77,78], the latter lack the E T criterion. Single-lepton final states have been analysed in searches for W bosons [79,80] which rely, however, on the partial reconstruction of an schannel resonance from the observed transverse mass not applicable here.
The ATLAS and CMS analyses cited above thus currently do not impose any direct limits on our model. The LHC data would have to be reanalysed with full signal, background and detector simulations. And even if we were to do so, which is beyond the scope of the present paper, the LHC data are not expected to significantly constrain the parameter space of this model.

Conclusions
We analyzed a TeV-scale extension of the SM that can simultaneously explain neutrino masses and the dark matter. The SM particle content was enlarged with a vector-like doublet of SU(2), a singlet fermion, and two singlet scalars, all of them assumed to be odd under a Z 2 symmetry that guarantees the stability of the dark matter particle. In this scenario, neutrino masses are generated radiatively via one-loop processes mediated by the new fields. We examined the case where the dark matter particle -the lightest odd field -is a scalar and its relic density is determined by its Yukawa interactions. In this case, interesting correlations appear between dark matter, neutrino masses, and lepton flavor violating processes. We studied analytically and numerically the phenomenology of this scenario. We found that the dark matter constraint can only be satisfied for a dark matter mass below 600 GeV and a charged fermion mass below 720 GeV. We argued that neither existing collider searches at the LHC nor the dark matter direct or indirect detection searches can significantly constrain this scenario. A generic prediction of this setup is instead the existence of sizable rates for several LFV processes. In fact, they provide the most promising way of probing this model in the near future. Future searches for μ → 3e and μ-e conversion in nuclei, in particular, have the potential to probe the entire parameter space of this model.