Sterile neutrino dark matter and low scale leptogenesis from a charged scalar

We show that novel paths to dark matter generation and baryogenesis are open when the standard model is extended with three sterile neutrinos \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_i$$\end{document}Ni and a charged scalar \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta ^+$$\end{document}δ+. Specifically, we propose a new production mechanism for the dark matter particle—a multi-keV sterile neutrino, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_1$$\end{document}N1—that does not depend on the active-sterile mixing angle and does not rely on a large primordial lepton asymmetry. Instead, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_1$$\end{document}N1 is produced, via freeze-in, by the decays of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta ^+$$\end{document}δ+ while it is in equilibrium in the early Universe. In addition, we demonstrate that, thanks to the couplings between the heavier sterile neutrinos \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_{2,3}$$\end{document}N2,3 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta ^+$$\end{document}δ+, baryogenesis via leptogenesis can be realized close to the electroweak scale. The lepton asymmetry is generated either by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_{2,3}$$\end{document}N2,3-decays for masses \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{2,3}\gtrsim $$\end{document}M2,3≳ TeV, or by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N_{2,3}$$\end{document}N2,3-oscillations for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{2,3}\sim $$\end{document}M2,3∼ GeV. Experimental signatures of this scenario include an X-ray line from dark matter decays, and the direct production of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta ^+$$\end{document}δ+ at the LHC. This model thus describes a minimal, testable scenario for neutrino masses, the baryon asymmetry, and dark matter.


Introduction
Sterile neutrinos, that is, fermions singlet under the SU (3) × SU (2) × U (1) gauge symmetry, are a very well-motivated extension of the standard model (SM). On the theoretical side, they are a prediction of left-right symmetric theories, they allow to gauge B − L by removing its anomaly and they are necessary in SO (10) grand unification. On the phenomenological side, they provide a non-vanishing mass to the active SM neutrinos, they allow to realize baryogenesis via leptogenesis, and they are also a viable candidate for dark matter, as long as their mass lies in the keV range. In that case, they are naturally long-lived so that, unlike dark matter candidates at the electroweak scale, no additional symmetries are required to stabilize them. Thus, sterile neutrinos a e-mail: carlos.yaguna@uni-muenster.de can provide a simple solution to the three open problems of the SM: neutrino masses, the baryon asymmetry, and the dark matter.
The minimal model addressing these three issues requires one sterile neutrino N 1 at the keV scale as dark matter candidate [1], and two additional sterile neutrinos N 2,3 for leptogenesis, which is induced either by N -decays, for sterile neutrino masses above the TeV scale [2], or by N -oscillations, for sterile neutrino masses at the GeV scale [3,4]. Light active neutrino masses are easily generated, via the seesaw mechanism, provided the sterile neutrino masses are significantly larger than about 1 eV. Thus, one may argue that no new physics is needed above the electroweak scale to explain neutrino masses, baryogenesis, and dark matter, defining the so-called "ν minimal standard model" (νMSM) [5]. Such a model is undoubtedly economical and very predictive, by reason of the small number of degrees of freedom it contains, but it is also strongly constrained. Present experiments, in fact, already exclude the minimal scenario for sterile neutrino dark matter within the νMSM [6]. In the region of mass and mixing angle compatible with current observations, in particular the X-ray bounds [7][8][9] and the Ly α forest data [10][11][12], one cannot produce a large enough population of sterile neutrinos from standard active-sterile oscillations [1]. The only way out within the νMSM is to generate the dark matter neutrinos via resonant active-sterile oscillations triggered by large initial lepton asymmetries [13,14], but that mechanism requires a high degree of fine-tuning [6,15]. In addition, leptogenesis via N 2,3 -oscillations within the νMSM also requires a significant tuning of parameters, in particular a strong mass degeneracy between N 2 and N 3 [6].
It is important, therefore, to consider alternative ways of realizing leptogenesis and producing sterile neutrinos within extensions of the νMSM. Indeed, most extensions of the SM that are theoretically well-motivated contain new degrees of freedom beside the sterile neutrinos, which may modify significantly the phenomenology of the νMSM. For leptogen-esis, the required tuning of parameters can be released by introducing a second Higgs doublet, as recently studied in [16]. Regarding dark matter, several other possibilities have been considered for the production of the sterile neutrinos. They include the production through inflaton decay [17][18][19], thermal overproduction followed by entropy dilution [20], and the production via the decays of either a neutral scalar in thermal equilibrium [21], or a frozen-in neutral scalar [22].
In this paper we advocate an extension of the SM by three sterile neutrinos, N i , plus a charged scalar, δ + , both of which are naturally present in left-right symmetric or unified extensions of the SM [23]. The charged scalar interacts with SM leptons, both doublets and singlets, and with sterile neutrinos, and its mass can be as low as allowed by current collider bounds-about 200 GeV. We will show that the interplay between N i and δ + modifies the dark matter phenomenology, enables new realizations of baryogenesis, and gives rise to novel experimental signatures at colliders. The decays of δ + , in fact, constitute a new mechanism for the production of sterile neutrino dark matter in the early Universe. An advantage of this mechanism is that the resulting relic density does not depend on the active-sterile mixing angle, allowing to satisfy the strong bounds derived from X-ray observations [7][8][9]. Besides, the charged scalar induces independent sources of leptogenesis that are effective at scales close to the electroweak scale. If the charged scalar mass were accessible at colliders, the present scenario for dark matter and baryogenesis would be testable in a direct way.
The rest of the paper is organized as follows. In Sect. 2 we introduce the model and fix our notation. Section 3 deals with dark matter production via the decays of δ + . We obtain, in particular, the dark matter relic density as a function of the parameters of the model. Baryogenesis via leptogenesis is discussed in Sect. 4. We examine two different ways to generate a lepton asymmetry, and outline the region of the parameter space that successfully realizes baryogenesis in each case. The most interesting experimental signatures of this model are discussed in Sect. 5. Section 6 briefly elaborates on the flavour structure of the masses and Yukawa couplings of the model. We present our conclusions in Sect. 7.

The model
When gauge singlet chiral fermions N Ri are added to the SM, they have in general a Majorana mass term and a Yukawa coupling to the SM lepton doublets, When a scalar field δ + , singlet under SU (2) and with electromagnetic charge one, is added to the SM, it has in general, besides its gauge and self-interactions, a quartic coupling to the Higgs doublet and an antisymmetric Yukawa coupling to two lepton doublets, When both the sterile neutrinos N Ri and the charged scalar δ + are present, there is an additional Yukawa coupling involving the SM lepton singlets, These three terms, Eqs. (1)-(3), plus the well-known SM ones, constitute the Lagrangian considered in this paper. It is the most general Lagrangian consistent with the SU (3) × SU (2) × U (1) gauge symmetry and with the addition of the singlet fermions and the charged scalar. Without loss of generality, we adopt the basis where M N and the charged lepton mass matrix are real and diagonal. For our following discussion, the parameters of greater relevance are the masses of the singlet fermions (M i , i = 1, 2, 3) and of the charged scalar (M δ ), and the 3 × 3 Yukawa matrices y ν and y R . In this model active neutrino masses are not affected by the existence of δ + and are obtained via the usual seesaw mechanism: m ν y ν M −1 N y T ν v 2 , with v 174 GeV. In Sect. 3 we will identify the sterile neutrino N 1 with a multi-keV dark matter candidate, with very small mixing angles with active neutrinos, corresponding to tiny neutrino Yukawa couplings, (y ν ) α1 3 × 10 −13 . These tiny couplings give rise to a lightest neutrino mass m lightest ν 10 −6 eV. The much larger solar and atmospheric neutrino mass scales are generated instead by the seesaw contribution of the two heavier sterile neutrinos, N 2,3 . The active neutrino mass spectrum is thus hierarchical (either normal or inverted), with one neutrino remaining essentially massless.
The generation of the dark matter energy density and of the baryon asymmetry have been extensively studied in the context of the SM extended with sterile neutrinos only [6]. The SM extension with sterile neutrinos and δ + was considered in [23]. It was shown that, taking the sterile neutrinos N i as light as a few TeVs and heavier than δ + , one can generate the observed baryon asymmetry via leptogenesis, with no need to enhance resonantly the CP asymmetry.
We explore, instead, the possibility that the lightest singlet fermion, N 1 , has a multi-keV mass, M 1 M δ , and accounts for the observed dark matter density. In this framework, the existence of δ + offers an alternative way of producing dark matter that has not been studied before in the literature. Moreover, we will show that this different mass spectrum is still compatible with the leptogenesis mechanism proposed in [23] and, in addition, we will identify a different mechanism for leptogenesis in the presence of δ + based on N 2,3 -oscillations.

Dark matter production from δ + decays
In this model, the only particle that can play the role of dark matter is the lightest singlet fermion, N 1 . It is important, therefore, to determine how it is produced in the early Universe and whether it can account for the observed dark matter density while respecting all other experimental constraints. In this section, after briefly reviewing the standard scenario for sterile neutrino dark matter, we propose a new mechanism for dark matter production via the decays of δ + .
In the absence of δ + , our model is described by the seesaw Lagrangian, Eq. (1), and the region of the parameter space compatible with dark matter is the one of the νMSM. It features a singlet fermion (N 1 ) with a mass at the keV scale, the dark matter particle, and two heavier singlets, N 2,3 . In the νMSM, dark matter is produced at temperatures of order 100 MeV via active-sterile neutrino mixing. The required effective mixing angle, θ 1 , is necessarily small and is related to the Yukawa couplings by θ 2 1 = α=e,μ,τ |(y ν ) α1 | 2 v 2 /M 2 1 . Detailed numerical studies have shown that the observed dark matter density can be reproduced in two different cases [6]. In the non-resonant case (also known as the Dodelson-Widrow mechanism [1]), dark matter neutrinos are thermally and non-resonantly produced with a smooth distribution of momenta. This mechanism yields the minimal amount of dark matter that can be obtained for a given mass and mixing angle. For dark matter masses between 1 keV and 10 keV, the mixing angle required to explain the dark matter density is sin 2 (2θ 1 ) ∼ 10 −8 -10 −9 , as shown e.g. in figure 2 of [6]. These parameters are not consistent with current data [24], which imply M 1 3 − 4 keV from the X-ray line bound [7][8][9] and M 1 8 keV from the Ly α forest observations [10][11][12]. The other case in which the dark matter density can be explained, and the only one that is currently viable in the νMSM, is resonant production (also known as the Shi-Fuller mechanism [13]). In it, the dark matter production rate is resonantly amplified by the presence of a lepton chemical potential in the plasma, which enhances the production for particular momenta as they pass through the resonance, giving rise to a non-thermal momentum distribution that is colder than that obtained in the non-resonant case. To explain the observed dark matter density, a lepton asymmetry |μ α | 8 × 10 −6 at T ∼ 100 MeV is required [14], where μ α = n α /s, s is the entropy density of the Universe, and n α is the total number density of active leptons of flavour α. In principle, this large lepton asymmetry can be generated within the νMSM via CP-violating oscillations of N 2 and N 3 , but only if their mass difference, M 3 − M 2 , is fine-tuned to the order of 10 −11 [15,25].
The existence of δ + allows a new production mechanism for sterile neutrinos in the early Universe. Indeed, thanks to the interaction term in Eq. (3), the decays of the charged scalars can produce sterile neutrinos, δ + → N 1 + . This decay is a thermal process that takes place while δ + is in equilibrium and it is an example of the so-called freeze-in scenario [26] for dark matter production. The defining feature of freeze-in is that the dark matter interactions are so weak that they never reach thermal equilibrium in the early Universe. As a result, the dark matter abundance is negligible at high temperatures and slowly increases as the Universe cools down. The production ceases when the freeze-in temperature is reached, and the dark matter abundance remains constant from then on. The , from δ + decays is obtained by solving the following Boltzmann equation [26] sT where s is the entropy density, H is the expansion rate and γ N 1 (T ) is the thermally-averaged production rate. We have that where K 1 (x) is the Bessel function of the second kind and the sum runs over the different lepton flavours. Since the lepton and N 1 masses are negligible compared to M δ , the δ + decay rates are calculated as and their sum over flavours can be conveniently written as where we have defined y R1 as the combination of couplings that determines the production of N 1 via δ + decays. Using

R1
M δ x max x min where T 0 is the temperature today and we used Y (T M δ ) = 0 as our initial condition. Integrating this equation from x max = ∞ to x min = 0 and setting g s ≈ g ρ ≈ 100 yields The N 1 relic density is then Thus, a keV sterile neutrino can explain the observed dark matter density via δ + decays if M δ ∼ 1 TeV and y R1 ∼ few × 10 −8 . This analytical result is in very good agreement with the numerical calculations we did for similar models of freeze-in dark matter [27,28]. A crucial feature of Eq. (10) is that, unlike the production in the νMSM, it does not depend on the active-sterile mixing angle θ 1 that determines the decay rate of the sterile neutrino. Hence, one can take θ 1 small enough to be in agreement with the X-ray bounds without affecting the predicted dark matter density. In other words, in this model one can decouple the dark matter production (determined by the y R couplings) from the dark matter decay (determined by y ν via mixing angles), opening new viable regions and alleviating the experimental constraints on sterile neutrino dark matter. In particular, M 1 can be larger than the keV scale, as long as the dark matter sterile neutrino is sufficiently long-lived. Since N 1 ∝ M 5 1 sin 2 θ 1 , one cannot exclude, for example, the possibility of vanishing active-sterile mixing, θ 1 = 0, which would allow masses as high as 1 MeV. Above that value, the decay N 1 → νe + e − becomes kinematically available and can proceed through the couplings y R and y L , even for θ 1 = 0. In that case, the strong constraints coming from indirect searches of electrons and positrons come into play, but they compete with the smallness of y R1 , determined by Eq. (10), and of the Yukawa matrix y L , that can be very small too.
Regarding structure formation, the sterile neutrinos produced via δ + decays are colder than those obtained in the νMSM [29,30] because the decays take place at a temperature T ∼ M δ much higher than the QCD scale (∼ 150 MeV). In fact, their free-streaming length and phase space density are identical to those of the so-called chilled sterile neutrinos studied in [31]. Consequently, the Lyman-α bounds on the sterile neutrino mass in our scenario are significantly weaker, allowing for a dark matter mass M 1 smaller than in the νMSM [32].
We assume in the following that the decays of the charged scalars are the dominant source of sterile neutrinos, so that Eq. (10) is satisfied. This will be the case in the absence of a large lepton asymmetry and for small mixing angles, sin 2 2θ 1 < 10 −9 . Next, we show that this assumption is consistent with baryogenesis, and that it leads to new experimental signatures.

Baryogenesis via leptogenesis with δ +
In this section we will discuss leptogenesis in the presence of the charged scalar field δ + . We will briefly review the two basic mechanism for leptogenesis in the minimal scenario with sterile neutrinos only: N -decays and N -oscillations. We will show that both mechanisms are still operative when one replaces the Yukawa coupling y ν in Eq. (1) with y R in Eq. (3), that is, when the role of the SM Higgs doublet H is played by δ + , and the role of the SM lepton doublet l L is played by the SM lepton singlet e R .
The baryon asymmetry of the Universe can be generated from a lepton asymmetry, as long as the latter is present before the electroweak phase transition, which occurs at a temperature T EW 150 GeV. Above this temperature the (B + L)violating electroweak sphalerons are in thermal equilibrium, thus converting efficiently leptons into baryons. When the SM is extended by sterile neutrinos, one can generate a lepton asymmetry, provided that the set of couplings (M N ) i j , (y ν ) αi and (y R ) αi violates the C P symmetry, and (some of) these couplings are out-of-equilibrium at some epoch before T EW .

Leptogenesis from N -decays
The traditional leptogenesis mechanism [2] assumes a Majorana mass matrix M N for two (or more) sterile neutrinos, with eigenvalues M i larger than the electroweak scale, and it relies on the out-of-equilibrium decays N i → Hl Lα , at temperatures just below M i . As N 1 plays the role of keV-scale dark matter candidate, the relevant decays are those of N 2 , with T EW < M 2 < M 3 . Since M 2 violates lepton number, the C P-asymmetry that is generated in the decays is also a lepton asymmetry, or equivalently a (B − L)-asymmetry. After all the N 2 particles decayed, the produced (B − L)asymmetry remains constant. In this scenario, the couplings that control leptogenesis are the same Yukawa couplings, y ν , that determine the active neutrino masses. As a consequence, a sufficient amount of baryon asymmetry can be generated only for M 2 10 8 GeV approximately. 1 Thus, no direct test of this scenario can be envisaged.
In the presence of a charged scalar δ + , an analog source of leptogenesis is provided by the decays N i → δ + e Rα , mediated by the Yukawa coupling y R , as shown in [23]. Here we adapt the results of that paper to our new context, where the dark matter particle is the keV-scale sterile neutrino N 1 . We can safely neglect the presence of N 1 for the computation of the lepton asymmetry because its interaction rates (in particular those violating lepton number) are strongly outof-equilibrium above T E W , due to the extreme smallness of M 1 , (y R ) α1 and (y ν ) α1 . Thus, at least two extra sterile neutrinos N 2,3 heavier than δ + are needed to generate the C Pasymmetry. Taking for simplicity M 3 M 2 M δ , the C P-asymmetry in N 2 -decays is given by where we assumed that the charged scalar lepton number is L(δ + ) = −2, since the dominant decay mode is δ + → e + ν, through the Yukawa coupling y L in Eq. (2). 2 To reproduce the observed baryon asymmetry two basic conditions are required [23]: a minimal value for N 2 assuming no washout, and N 2 -decays out-of-equilibrium at T = M 2 to avoid large washout from inverse decays, Thus, one can realize leptogenesis for M 2 as small as a few TeVs, as long as N 2,3 have hierarchical Yukawa couplings: It is worth reminding that, when M 2 approaches the electroweak scale, a sufficiently large CP-asymmetry implies in general strong (B − L)-washout rates, in particular those mediated by N 3 off-shell, unless special conditions are realized. A neat, model-independent discussion of the lower bound on the leptogenesis scale can be found in [36,37]. In the present model with M 2 ∼ a few TeVs, the washout can be Boltzmann suppressed by raising M δ sufficiently close to M 2 , and taking a sufficiently large coupling y L in Eq. (2) to avoid washout from the asymmetry stored in δ + . 3 Note that the asymmetries generated by N 2 -decays through the Yukawa coupling matrix y ν become negligible in the region M 2 10 8 GeV due to the seesaw relation. The y Rentries, on the other hand, are not constrained by the light neutrino masses, and leptogenesis can work close to the electroweak scale as described above. This scenario is, therefore, easier to test directly, by the observation of δ + (and possibly N 2 ) at colliders, as discussed in Sect. 5.

Leptogenesis from N -oscillations
A complementary mechanism for leptogenesis [3] relies on N -oscillations rather than on N -decays. It also requires two (or more) sterile neutrinos, coupled to the SM through the Yukawa matrix y ν . In the early Universe, the thermal population of lepton doublets l Lα produces, through small, out-ofequilibrium couplings (y ν ) αi , a coherent superposition N α of the sterile neutrino mass eigenstates N i . These sterile neutrinos, at temperatures much larger than their masses, coherently oscillate among the different flavours α. Such oscillations conserve lepton number (the N α conserve their helicity), but violate lepton flavour numbers. If the C P-symmetry is also violated, one generates non-zero flavour asymmetries between the opposite helicities of the N α . The asymmetry in the flavour α is transferred efficiently to l Lα , as long as the coupling (y ν ) αi goes into equilibrium for some i. Since the total lepton asymmetry, that is, the sum over α of the flavour asymmetries vanishes, one needs that some but not all flavours go into equilibrium before T EW , so that a net lepton asymmetry remains stored in the sterile neutrino sector, and an opposite one is available in the SM sector to be transferred to the baryons by electroweak sphalerons. The Yukawa (out-of-)equilibrium condition at T E W reads |(y ν ) αi | 10 −7 for some α and i, |(y ν ) βi | 10 −7 for some β = α, ∀i.
A number of additional constraints, on the size of y ν -entries and on the values of M i , must be satisfied for this leptogenesis mechanism to work. To facilitate the comparison between the scenarios with and without the charged scalar, we will describe them in some detail. Due to the dark matter constraint, N 1 plays no role in leptogenesis and the asymmetry must be generated by N 2 and N 3 only. 4 Successful leptogenesis implies a few upper bounds on the y ν -entries, beside the flavour-dependent one in Eq. (14). First of all, note that the Majorana masses M i violate lepton number. Therefore, they play the role of washout in this scenario, as they transform opposite helicities of the N α into one another. The condition to keep the lepton number violation rate out-of-equilibrium down to T E W is approximately |(y ν ) αi | 10 −5 GeV M i , ∀α and ∀i, which combined with Eq. (14) implies sterile neutrino masses below the electroweak scale, M i 100 GeV. 5 In addition, for the case of two sterile neutrinos N 2,3 , analytic and numerical studies [4,16,39] show that a sufficient baryon asymmetry requires a strong degeneracy between M 2 and M 3 , with M/M 10 −5 . This comes from the interplay of a number of subtle effects: the oscillation time increases as M decreases, and this allows for larger asymmetries because the Yukawa interaction rates are closer to equilibrium at later times; note also that larger Yukawa couplings (for some but not all the flavours) enhance the flavour asymmetries, but also tend to spoil the coherence of the sterile neutrinos, as active-sterile transitions may become faster than the oscillation time. We extrapolate the resulting constraint from figure 7 of [6], that can be written as The masses M 2,3 are below collider energies, but it is difficult to produce them directly since they have small couplings. Still, the mixing with active neutrinos can be sufficiently large to have an observable effect in various neutrino experiments [40]. The resulting upper bound on active-sterile mixing is also shown in figure 7 of [6]. An experimental proposal to improve significantly the present bound can be found in [41]. In less minimal models, e.g. with three GeV scale sterile neutrinos, there are good detection perspectives in meson decay experiments [42].
On the other hand, several observables other than the baryon asymmetry put lower bounds on the y ν -entries, that add to the one in Eq. (14). The lower bound on active neutrino masses, |(m ν ) αβ | 0.025 eV for some α and β, implies a lower bound on the neutrino Yukawa couplings through the seesaw formula, for some α and i.
Indeed, since N 1 gives a negligible contribution to m ν , one active neutrino is approximately massless, therefore |(m ν ) αβ | 0.05 eV for all α and β. Then, the right-hand 5 We remark that most numerical studies in the literature assume that the lepton number violating rates mediated by M i are negligible; in view of the size of the Yukawa couplings relevant for leptogenesis, this assumption seems to be justified only for M i GeV. For larger M i , one should include the lepton number violating rates in the Boltzmann equations. In this regime, it remains conceivable that the C P-asymmetries are large enough to compensate the relatively strong washout. This possibility was recently entertained to make this scenario work even for M i larger than T EW [38]. side of Eq. (17) provides the natural value of the Yukawa couplings for i = 2, 3: the largest (y ν ) αi are, the strongest is the cancellation needed among the contributions of N 2 and N 3 to m ν . These GeV-scale sterile neutrinos should decay (e.g. in 3ν's or νe + e − ) before T 1 MeV, not to spoil nucleosynthesis (see [43] for a detailed analysis). We roughly estimate this constraint as , for some α, ∀i.
Comparing with Eq. (16), this implies M i 0.1 GeV. All the constraints above confirm that N 1 plays no role in leptogenesis.
Let us show that the same mechanism of leptogenesis through N -oscillations is operative in the presence of the charge scalar δ + , by replacing the role of y ν with y R . Analogously to the previous case, the couplings (y R ) αi must be small to remain out-of-equilibrium while they slowly produce coherent sterile neutrino states. The lepton flavour asymmetries generated by N -oscillations are (partially) transferred to the SM lepton singlets e Rα . The latter are in equilibrium with l Lα through the charged lepton Yukawa couplings (at least for α = μ, τ ), thus electroweak sphalerons transfer the asymmetries to baryons as usual.
One needs that the y R -entries satisfy the same inequalities as the y ν -entries in Eqs. (14) and (15), so that some flavour asymmetries are transferred to baryons, while the others remain stored in the sterile sector, and the washout from lepton number violating scattering is small. The only differences amount to (i) order one factors to account for the singlet (doublet) nature of e R (l L ) in scattering rates; (ii) the range of temperatures where M δ (M H ) can be neglected: note that the asymmetry is generated at some scale well above T EW through scattering processes mediated by δ + (H ).
We also expect that, when both y ν and y R take values relevant for leptogenesis, there is more freedom to generate large asymmetries, e.g. because of the presence of extra C P-violating phases, and the requirement of a strong degeneracy between M 2 and M 3 could be relaxed. As in the case of y ν , one cannot raise too much the value of y R -entries, as the coherence of sterile neutrino oscillations requires the scattering rate with the SM leptons to be out-of-equilibrium at the time of oscillations; this should translate in a bound similar to the one in Eq. (16). A numerical study is needed to establish more precisely the allowed parameter space, and the differences with respect to the νMSM one. The effect of the two sets of Yukawa couplings could be dramatic, as is the effect of a third GeV scale sterile neutrino [42,44], or of a second Higgs doublet [16].
There is no constraint on y R coming from the seesaw relation, so one can take couplings smaller than in Eq. (17), slightly enlarging the region of parameters of leptogenesis.
Perhaps more importantly, one can take (y R ) αi significantly larger than the right-hand side of Eq. (17), with no need of fine-tuning to keep m ν small. Note that the coupling y R does not induce any mixing with active neutrinos. Indeed, direct searches of active-sterile mixing are presently sensitive to y ν -entries much larger than in Eq. (17): if the seesaw parameters take their natural values, no direct signal of active-sterile mixing is expected.
Finally, the three-body decay rate of N i through the coupling y R is proportional to |(y R ) αi (y L ) βγ | 2 /M 4 δ , to be compared with |(y ν ) αi | 2 /(M 2 i v 2 ) for a decay through the mixing with active neutrinos. Since the latter is typically much faster, the nucleosynthesis bound applies to y ν only, in the form of Eq. (18). Then, leptogenesis through y R could work even for M i < 0.1 GeV, but in this region nucleosynthesis demands large values for the y ν -entries, that require strong cancellations in the seesaw.

Experimental signatures for N 1 and δ +
Two important differences between our model and the νMSM are the additional region of the parameter space where the dark matter constraint can be satisfied and the presence of the extra charged scalar δ + . They both give rise to new experimental signatures that may allow to distinguish one model from the other.

Dark matter indirect detection
A multi-keV dark matter neutrino is inherently unstable and decays into three light neutrinos at tree-level (N 1 → 3ν α ) and radiatively into a light neutrino and a photon (N 1 → ν α γ ). This radiative decay produces an X-ray line at E γ M 1 /2 that can be searched for and used to constrain the model or to help determine its parameter space [45]. The presence of the additional scalar δ + does not affect the decay modes of the dark matter neutrino, which are still determined by its mass and its mixing with the active neutrinos, just as in the νMSM, but it modifies the regions that are consistent with the dark matter constraint, allowing, in particular, for smaller mixing angles. This fact has important implications, as we show next.
Recently, the detection of an unidentified spectral line at about 3.5 keV has been reported from two independent data sets [46,47]. Arguments in favour or against the dark matter decay interpretation of the signal can be found in [48][49][50][51]. If confirmed, that signal would provide compelling evidence for keV-scale dark matter and, in particular, for dark matter in the form of sterile neutrinos. Within that framework, the signal can be explained if M 1 7 keV and sin 2 2θ 1 5×10 −11 . In the minimal scenario for sterile neutrino dark matter, where they are produced non-resonantly, such parameters lead to a relic density way below the range determined by cosmological observations. Thus, an additional source of sterile neutrinos is required. A simple possibility for that new source are the decays of the δ + particle, as explained in the previous section. These decays allow to decouple the dark matter production (determined by the y R couplings) from the dark matter decay (determined by y ν via mixing angles), opening new viable regions consistent with all bounds. Specifically, the freeze-in production of sterile neutrinos that we have examined in this paper enables to explain the tentative 3.5 keV line and to simultaneously account for the observed relic density, even in the absence of a large, primordial lepton asymmetry. If that line signal turns out to be spurious, X-ray observations will continue to be the main way in which the dark matter sector of this model can be tested in the foreseeable future.

Charged scalar searches at colliders
The singly-charged isosinglet scalar δ + could be directly produced at colliders, if its mass were within their energy range. The dominant production channel is the Drell-Yan process ψψ → γ /Z → δ + δ − , with ψ = e at LEP and ψ = q at Tevatron and LHC, with a partonic cross-section given e.g. in [52]. The δ + decays into one anti-lepton and one anti-neutrino, either through the Yukawa coupling (y L ) αβ , or (y R ) αi when M i < M δ . 6 The standard experimental searches for charged scalars usually assume production and decay modes different from those above. For a singly-charged scalar H + , the decay channel into lepton plus missing energy has been analyzed, but only for an isodoublet produced from a top-bottom vertex, as in type II two-Higgs doublet models, allowing for an effective background reduction [53,54]. Indeed, the signal over background ratio is significantly smaller in the case of δ + ; an analysis of the cuts required to maximize the signal has been presented in section 5.2 of [55].
Presently, it appears that the best way to constrain M δ is to use the experimental searches for supersymmetric particles. Since δ + and a right-handed slepton have the same gauge quantum numbers, when the latter is directly pair-produced and decays into a lepton and a light neutralino, it behaves very much as the former, and the same bounds apply to both particles. Recently, that setup was studied in the context of simplified supersymmetric models, for the case of a righthanded s-electron degenerate with a right-handed s-muon, decaying with branching ratio 100 % into electron and muon, respectively, plus missing energy [56,57]. Assuming a similar efficiency for electron and muon detection and reducing the signal by a factor two, to take into account that the δ + production cross-section is the same as for a single right-handed slepton, we can extract a 95 % C.L. lower bound M δ 200 GeV at ATLAS [56] (see also the auxiliary figures 11a and 12a available online 7 ) and M δ 160 GeV at CMS [57].
Note, however, that these bounds apply as long as the δ + branching ratio into taus is negligible. The pair-production of s-taus decaying into taus is presently unconstrained at the LHC (see figure 9 in [58]). As a consequence, if δ + decays mostly into taus, its mass can be significantly smaller, say down to the LEP bound ∼ 100 GeV. Indeed, in minimal flavour models the δ + Yukawa couplings are hierarchical as the SM Yukawa couplings (see Sect. 6), and a large branching ratio into taus is expected.

Flavour structure of the model
We accomplished our goal of demonstrating that this model can account for dark matter, neutrino masses, and baryogenesis. That is, we have determined, in the previous sections, that there exist regions in the parameter space of the model where the masses of the new particles and their Yukawa couplings are such that these three issues are simultaneously explained. One may wonder, nonetheless, whether those values of the masses and Yukawa couplings have a generic flavour structure, or if they require a special tuning, indicating that some flavour symmetries are operative. Previous works along this line include [59][60][61][62][63], which focused on the generation of the keV scale in neutrino models, and on the flavour structure of the νMSM. In this section, we do not aim to construct complete flavour models, but rather to recap the order of magnitude of the parameters needed phenomenologically, and to suggest some rationale to explain them.
In view of the hierarchical values of masses and couplings that are required for dark matter and leptogenesis, it is useful to describe the required flavour structure in terms of a U (1) F family symmetry, with different charges assigned to the various fermions. In this framework, the coefficient of each fermion bilinear ψ Li ψ L j is suppressed by a power q i + q j of a small parameter (we take the fermion U (1) Fcharges to be positive). One may assume some underlying flavour dynamics, that generates = φ / F 1, that is, the vev of a spurion field with q φ = −1, over the cutoff of the flavour theory.
We need M 1 M 2 ∼ M 3 . Note that in the νMSM a strong degeneracy of M 2 and M 3 is required, while in our scenario this constraint is relaxed. The straightforward way to realize this pattern is to take q N 1 > q N 2 = q N 3 , so that M 1 /M 2 ∼ 2n , with n ≡ q N 1 − q N 2 . For N -oscillation 7 https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/ SUSY-2013-11/.
(N -decay) leptogenesis, one needs 2n ∼ 10 keV/GeV = 10 −5 ( 2n 10 keV/TeV = 10 −8 ). This charge assignment implies automatically a hierarchy n : 1 : 1 among the three columns of the matrix y R (as well as of y ν ), in other words, (y R ) 2 αi /(y R ) 2 α j ∼ M i /M j . Since this relation is not satisfied in certain regions of parameters relevant for dark matter, leptogenesis and neutrino masses, an additional flavour symmetry must be operative. The Yukawa couplings can be further suppressed by introducing parity symmetries Z (i) 2 :N i → −N i , that allow for M i but forbid (y ν ) αi and (y R ) αi ; these couplings must then be proportional to a small symmetry-breaking parameter i . To suppress (y ν ) αi and not (y R ) αi (or vice versa), one may argue that under the same parity the combination e Rα δ + (or l Lα H ) is also odd. 8 Coming to the U (1) F charges of the three families of lepton doublets and singlets, q Lα and q Rα , first of all they determine the hierarchy of the charged lepton masses, m e : m μ : m τ ∼ q Le +q Re : q Lμ +q Rμ : q Lτ +q Rτ . In turn, the hierarchy among the rows of y R (y ν ) is determined by the charges q Rα (q Lα ). Note that the charges of l Lα and e Rα are important for charged lepton and active neutrino masses, as well as for leptogenesis from N -oscillations; however they are not very relevant for dark matter production nor for leptogenesis from N -decays.
Let us confront these simple flavour symmetries with the values of the parameters needed for dark matter, neutrino masses and leptogenesis. The freeze-in of the desired amount of N 1 from δ + decays requires α |(y R ) α1 | 2 2.5× 10 −24 M δ /M 1 (see Eq. (10)). The N 1 -production from activesterile oscillations is negligible for |(y ν ) α1 | 10 −13 (M 1 /keV) (barring a large primordial lepton asymmetry). To generate large enough active neutrino masses, Eq. (17) must be satisfied. Clearly N 1 does not contribute significantly to m ν because of the dark matter constraint, therefore N 2 and N 3 are responsible to generate the atmospheric and solar mass scales. Since m 2 sol / m 2 atm 0.2, Eq. (17) must hold for both i = 2 and 3, up to a factor of a few. Coming to leptogenesis, the two mechanisms discussed in Sect. 4 correspond to two very different regions of parameters that we discuss in turn.
In the case of leptogenesis from N -decays through the y R -interaction, one requires M 2,3 > TeV, as well as |(y R ) α3 | > 10 −3 √ M 3 /M 2 for at least one α, and |(y R ) α2 | < 10 −4 M 2 /(10 8 GeV) for all α. For definiteness, let us take M 1 ∼ keV, M δ ∼ TeV and M 2,3 ∼ 100 TeV (leptogenesis through y ν is irrelevant at such low scales, barring resonances). Taking = 0.1 and the U (1) F charges q N 1 = 7, q N 2 = q N 3 = 2, q Rτ = 0 ≥ q Rμ , q Re , q Lτ = 2 ≥ q Lμ , q Le , we can reproduce the correct size of the M N , y R and y ν entries, except for the second column of y R , that must be further suppressed by R2 10 −3 , and the first column of y ν , to be further suppressed by ν1 10 −4 . As discussed above, R2 ( ν1 ) can be associated to the breaking of N 2 -parity (N 1parity). Note that here M 2,3 ∼ 4 L , with a lepton number violation scale L ∼ 10 9 GeV.
In the case of leptogenesis from N -oscillations coupled to the SM through y R , one requires for i = 2, 3, M i < 100 GeV, |(y R ) αi |M i 10 −5 GeV, as well as |(y R ) αi | > 10 −7 for some flavour α and < 10 −7 for some different flavour β. For definiteness, let us take M 1 ∼ 10 keV, M 2,3 ∼ 10 GeV and M δ ∼ TeV. Then, one can reproduce the correct size of the M N and y R entries taking = 0.1 and U (1) F charges q N 1 = 8, q N 2 = q N 3 = 5, q Rτ = 1, q Rμ = 2 and q Re = 3, and charged lepton masses require q Lτ = 1 ≥ q Lμ , q Le . Coming to the structure of y ν , one needs an extra source of suppression for (y ν ) α1 , by a factor ν1 10 −3 , to avoid N 1 -overproduction from active-sterile mixing; as before, ν1 can be understood as the breaking parameter of N 1 -parity. We note that, for the chosen values of parameters, both y R and y ν couplings to N 2,3 are relevant for leptogenesis. Here lepton number is broken at the scale L ∼ M 2,3 / 10 ∼ 10 11 GeV.
In summary, the regions of parameters where our scenario is effective require a well-defined hierarchy of masses and couplings, that point to specific but relatively natural flavour structures, with no need of extreme tuning of parameters.

Conclusions
We have shown that, when the SM is extended with three sterile neutrinos N i and a charged scalar δ + , baryogenesis via leptogenesis and the production of dark matter can both be mediated by δ + .
The dark matter candidate in this scenario is the lightest sterile neutrino N 1 , with a mass M 1 larger than about a few keVs. Such sterile neutrinos are produced by the decays of δ + while it is in thermal equilibrium in the early Universe. We showed that this novel mechanism, which is an example of freeze-in dark matter production, can account for the observed dark matter density. Moreover, since the resulting relic density does not depend on the active-sterile mixing angle θ 1 , the constraints from X-ray data are alleviated. The recent hint of a signal at 3.5 keV can be easily explained by N 1 -decays, for θ 1 ∼ 3 × 10 −6 . Alternatively, for a vanishingly small θ 1 , the dark matter mass M 1 can be raised to the MeV scale, or beyond.
Leptogenesis can be realized through the coupling y R between δ + and N 2,3 in two different ways, either via Ndecays or via N -oscillations. Both mechanisms are well-known in the context of the SM plus sterile neutrinos only, as they can both proceed through the neutrino Yukawa coupling y ν . The latter is constrained by light neutrino masses, restricting significantly the allowed parameter space and the testability of leptogenesis. In contrast, in our scenario y R is not involved in the neutrino mass generation. As a consequence, the N -decay mechanism is successful for M 2,3 as small as a few TeVs, with no need of a resonant enhancement. The parameter window for N -oscillation leptogenesis is similar when one employs the coupling matrix y R instead of y ν : one needs roughly M 2,3 ∼ GeV and y R -entries smaller than 10 −6 . However, y R introduces additional sources of C Pviolation, and moreover it is not constrained neither by light neutrino masses nor by nucleosynthesis. Thus one can largely relieve the fine-tuning of parameters necessary in the case with y ν only.
Neutrino masses are generated by a low scale seesaw mechanism involving the two heavier singlets, N 2,3 , the N 1contribution being negligible. In this way one can accommodate the current data on neutrino masses and mixing angles, with the prediction of an almost massless lightest neutrino. At colliders such as the LHC, this scenario can be probed and constrained via the direct production of δ + . This minimal extension of the SM thus provides a simple and testable scenario to explain neutrino masses, the dark matter, and the baryon asymmetry of the Universe.