Leptogenesis from oscillations and dark matter

An extension of the Standard Model with Majorana singlet fermions in the 1–100 GeV range can explain the light neutrino masses and give rise to a baryon asymmetry at freeze-in of the heavy states, via their CP-violating oscillations. In this paper we consider extending this scenario to also explain dark matter. We find that a very weakly coupled B-L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B-L$$\end{document} gauge boson, an invisible QCD axion model, and the singlet majoron model can simultaneously account for dark matter and the baryon asymmetry.


Introduction
The standard model (SM) of particle physics needs to be extended to explain neutrino masses, the missing gravitating matter (DM) and the observed matter-antimatter asymmetry in the universe.
Some of the most minimal extensions of the SM include new fermions, namely two or three sterile Majorana neutrinos (singlets under the full gauge group), which can account for the tiny neutrino masses, through the seesaw mechanism [1][2][3][4], and explain the observed matter-antimatter asymmetry through leptogenesis [5]. The simplest version of leptogenesis establishes the source of the matter-antimatter asymmetry in the CP violating the out-of-equilibrium decay of the heavy neutrinos. This scenario requires however relatively large Majorana masses > 10 8 GeV [6] (or ∼ 10 6 GeV with flavour effects [7] included), which makes these models difficult to test experimentally. For Majorana neutrinos in the 1-100 GeV range, it has been shown by Akhmedov, Rubakov and Smirnov (ARS) [8] and refined by Asaka and Shaposhnikov (AS) [9] that a different mechanism of leptogenesis is at work. In this case the asymmetries are produced at freezein of the sterile states via their CP-violating oscillations. The a e-mail: andrea.caputo@uv.es b e-mail: pilar.hernandez@cern.ch c e-mail: nuria.rius@ific.uv.es original ARS proposal did not include flavour effects and needs at least three Majorana species, while AS have shown that flavour effects can make it work with just two species. In the rest of the paper we will refer indistinctively to both scenarios as baryogenesis from oscillations (BO). In both cases, the lepton asymmetry is reprocessed into a baryonic one by electroweak sphalerons [10]. The extra heavy neutrinos in this case could be produced and searched for in beam dump experiments and colliders (see [11][12][13][14][15][16][17][18][19][20][21][22][23] for an incomplete list of works), possibly giving rise to spectacular signals such as displaced vertices [15,16,18,20,21]. Since only two sterile neutrinos are needed to generate the baryon asymmetry [9,[24][25][26][27][28][29][30][31][32][33][34][35][36][37][38][39], the lightest sterile neutrino in the keV range can be very weakly coupled and play the role of DM [40]. This is the famous νMSM [9]. However the stringent X-ray bounds imply that this scenario can only work in the presence of a leptonic asymmetry [41] significantly larger than the baryonic one, which is quite difficult to achieve. A recent update of astrophysical bounds on this scenario can be found in [42,43].
In this paper our main goal is to consider scenarios compatible with Majorana masses in the 1-100 GeV range and study the conditions under which the models can explain DM without spoiling ARS leptogenesis. 1 In particular, we will focus on models that are minimal extensions of the type I seesaw model with three singlet neutrinos. We will first consider an extension involving a gauged B − L model [45], which includes an extra gauge boson and can explain DM in the form of a non-thermal keV neutrino. We will then consider an extension which includes a CP axion [46] that can solve the strong CP problem and explain DM in the form of cold axions. Finally we consider the majoron singlet model [47,48] which can also explain DM under certain conditions both in the form of a heavy majorana neutrino or a majoron.
The plan of the paper is as follows. We start by briefly reviewing the ARS mechanism and the essential ingredients and conditions that need to be met when the sterile neutrinos have new interactions. In Sect. 3 we discuss the gauged B − L model, in Sect. 4, we study the invisible axion model with sterile neutrinos and in Sect. 5 we reconsider the singlet majoron model. In Sect. 6 we conclude.

Leptogenesis from oscillations
For a recent extensive review of the ARS mechanism see [49]. The model is just the type I seesaw model with three neutrino singlets, N i , i = 1 − 3, which interact with the SM only through their Yukawa couplings. The Lagrangian in the Majorana mass basis is In the early Universe before the electroweak (EW) phase transition, the singlet neutrinos are produced through their Yukawa couplings in flavour states, which are linear combinations of the mass eigenstates. Singlet neutrinos then oscillate, and since CP is not conserved, lepton number L gets unevenly distributed between different flavours. At high enough temperatures T m N i , total lepton number vanishes, in spite of which a surplus of baryons over antibaryons can be produced, because the flavoured lepton asymmetries are stored in the different species and transferred at different rates to the baryons. As long as full equilibration of the sterile states is not reached before the EW phase transition (T EW ∼ 140GeV) , when sphaleron processes freeze-out, a net baryon asymmetry survives. It is essential that at least one of the sterile neutrinos does not equilibrate by t EW . The rate of interactions of these neutrinos at temperatures much higher than their mass can be estimated to be where y α are the eigenvalues of the neutrino Yukawa matrix, T is the temperature and κ = few 10 −3 [50][51][52]. The Hubble expansion rate in the radiation dominated era is where g * is the number of relativistic degrees of freedom (g * ∼ 100 above the EW phase transition). The requirement that no equilibration is reached before t EW is: which implies yukawa couplings of order i.e. not much smaller than the electron yukawa. These yukawa couplings are compatible with the light neutrino masses for Majorana masses in the 1-100 GeV range. Any model that extends the one described above with new fields/interactions should be such that the new interactions do not increase the equilibration rate of the sterile neutrinos for the out-of-equilibrium requirement in the ARS mechanism to be met. We will now consider the implications of this requirement on various extensions of the minimal seesaw model of Eq. (1) that are well motivated by trying to explain also the dark matter and in one case also the strong CP problem.

B-L gauge symmetry
The SM is invariant under an accidental global U (1) B−L symmetry, that couples to baryon minus lepton number. If one promotes this symmetry to a local one [45], the model needs to be extended with three additional right handed neutrinos to avoid anomalies, which interestingly makes the type I seesaw model the minimal particle content compatible with this gauge symmetry. In this case, we have interactions between SM lepton and quark fields with the new gauge boson, V μ , as well as an additional term involving sterile neutrinos where Q f B−L = 1/3, −1 for quarks and leptons respectively. We also assume the presence of a scalar field φ, with charge B − L charge 2: that gets an expectation value φ , breaking B − L spontaneously, 2 and giving a mass to both the gauge boson and the sterile neutrinos: A massive higgs from the B − L breaking, σ , remains in the spectrum with a mass that we can assume to be M σ ∼ φ .
Existing constraints on this model come from direct searches for V in elastic neutrino-electron scattering, V gauge boson production at colliders, Drell-Yan processes and new flavour changing meson decays [53][54][55][56][57][58][59][60]. The status of these searches is summarized in Fig. 1, adapted from [61] (see also [62][63][64]). For masses, 1 GeV ≤ m V ≤ 10 GeV, g B−L is bounded to be smaller than ∼ 10 −4 , while the limit is weaker for larger masses. The improved prospects to search Ref. [61]. The dashed region labeled SHIP is the sensitivity of bremsstrahlung searches in SHIP [66,67]. The solid, dotted and dashed black lines correspond to the correct DM relic abundance in the form of sterile neutrinos of mass 1, 10 and 100 keV respectively. The dotted red line indicates the lower limit on the gauge coupling for which V is in thermal equilibrium. The solid red lines correspond to the upper bounds for successful BO leptogenesis for m N2,3 = 1 and 100 GeV  [68], BBN bounds estimated in [69] and the combined beam dump experiments from [70]. Unshaded regions represent the reach of SHIP from meson decay and bremsstrahlung searches [66]. The curves indicate the values of (g B−L , m V ) where the lightest sterile neutrino N 1 can account for the whole dark matter for three values of the neutrino masses, m N = 1, 10, 100 keV in solid, dotted and dashed lines for right-handed neutrinos exploiting the U (1) B−L interaction have been recently studied in [61], where the authors consider the displaced decay of the N at the LHC and the proposed SHIP beam dump experiment [65].
For m V ≤ 1GeV the strongest constraints come from supernova cooling [68,[71][72][73], beam dump searches [70,[74][75][76][77][78] and big bang nucleosynthesis (BBN) [69,70,[79][80][81]. Recent updates on these bounds are compiled in Fig. 2. The lower mass region labeled BBN is excluded by the effect of ΔN eff on the expansion, while the higher mass BBN region is excluded because the injection of electromagnetic energy from the V decay to charged particles during nucleosynthesis distorts the abundance of light elements. These BBN constraints have been evaluated in detail in Refs. [79,80]. In the relevant region of parameter space, they are seen to depend on the lifetime of the decaying particle and its abundance per baryon prior to decay. The corresponding region in Fig. 2 is a sketch of the excluded region in the latter analysis, which is approximately bounded by the lines corresponding to the lifetimes τ V ∈ [0.1−100]s, the threshold for V hadronic decays, m V ≥ 2m π , and the line corresponding to the fraction of the V decaying to charged particles per baryon, n em V /n b = 1. At least two of the neutrinos will be involved in the BO mechanism and their masses must be in the 1-100 GeV range. We need to ensure that the new B − L interactions do not bring them to thermal equilibrium before t EW , which will set an upper bound on g B−L . It it important to know however if the B − L gauge boson is in thermal equilibrium which also depends on g B−L . For m V 65 GeV, the dominant process is the scattering tt ↔ V H with a rate that can be estimated to be where y t is the Yukawa of the top quark and m t (T ) is its thermal mass. This rate is larger than the Hubble expansion somewhere above the EW transition provided g B−L 10 −6 . For larger masses of the gauge boson (m V 65GeV) the process V ↔ ff kinematically opens up at high temperature and one has to consider the decay and inverse decay with a rate where N C = 3(1) and Q f = 1/3(−1) for quarks(leptons). The sum is over all the standard model fermions whose thermal mass is such that m V (T ) ≥ 2m f (T ). The lower limit on g B−L for the thermalization of the B − L boson is shown as a dotted red line in Fig.1. Provided the B − L boson is in thermal equilibrium we have to consider its interactions with the sterile neutrinos driving leptogenesis. Assuming that the V boson is lighter than 2m N , so that the decay V → N N is kinematically forbidden, the dominant contribution [82] comes from the scattering processes with the fermions and Assuming and using Eq. (8) we get This upper bound is shown by the red solid line in Fig. 1. In the next section we will include these new interactions in the equations for the generation of the baryon asymmetry, and our results confirm the naive estimate in Eq. (13). Finally we consider the production of neutrinos through the decay of σ → N N, which is relevant for T ≥ M σ , since M σ m N . The requirement that this process does not thermalize the sterile neutrinos implies that the decay rate, Γ σ , is slower than the Hubble rate at T ≥ M σ : For m N = 1 − 100 GeV, using Eq. (8), we get On the other hand, for m V ≥ 2m N the dominant production goes via the decay of the gauge boson into two sterile neutrinos V → N N, which, if kinematically allowed, scales with g 2 B−L . In this case the decay rate is Requiring that it is smaller than One may worry if thermal mass corrections can allow the decay V → N N at large temperatures even if m V ≤ 2m N . At high enough temperatures both sterile neutrinos and the gauge boson acquire thermal corrections to the masses of the form The thermal mass of the gauge boson is larger than that of the sterile neutrino, because all fermions charged under B − L will contribute to the former and only the gauge boson loop contributes to the later [83]: We substitute the temperature dependent mass in Eq. (16) and we show in Fig. 4 , for m N = 1 and 100 GeV as a function of g B−L . The upper limit for g B−L are less stringent than those derived from 2 → 2 processes in Eq. (13).We now evaluate in detail the effect on BO induced by the new scatterings of Fig. 3. Leptogenesis in the presence of a new B − L gauge interaction has been recently studied in [82], although not in the context of BO, which as far as we know has not been considered before.

Leptogenesis
The sterile neutrinos relevant for leptogenesis are the heavier ones, m N 2,3 , with masses in the 1-100 GeV, and we focus on the scenario where m V ≤ 2m N 2,3 . We now explain how to include the terms involving the B − L gauge interactions in the quantum kinetic equations for BO leptogenesis as derived in [31]. Following Raffelt-Sigl approach [84], we consider a density matrix, ρ(k), describing the expectation value of number densities of N and a density matrix describing the corresponding anti-particles,ρ(k). 3 These equations are complemented with three equations involving the slow varying chemical potentials, μ B/3−L α . The modification of the kinetic equations induced by the B − L interactions is the addition of new collision terms in the equations for ρ andρ, that have the same flavour structure as the neutral current contribution considered in [84]. As explained above the most relevant contributions come from the scattering processes:f f ↔ NN and V V ↔ NN . The latter is enhanced at high temperatures, but of course will only be relevant if the V s are in thermal equilibrium, which we assume in the following.
The additional collision terms, from the first process, in the equation for the evolution ρ(k) can be writen in the form: where f eq i ≡ f F ( p i ) is the Fermi-Dirac equilibrium distribution function of the particle with momentum p i with p 4 ≡ k; {, } is the anticommutator, and the normalized matrices are: The additional collision terms forρ have the same form with the substitution k ↔ p 3 .
For the second process we have similarly: where f eq i is the equilibrium distribution function of the particle with momentum As usual we are interested in the evolution in an expanding universe, where the density matrices depend on momentum, y ≡ p/T and the scale factor or inverse temperature x ∝ T −1 . We consider the averaged momentum approximation, which assumes that all the momentum dependence factorizes in the Fermi-Dirac distribution and the density r is just a function of the scale factor, ie. ρ(x, y) = f F (y)r (x). In this approximation we can do the integration over momentum and the B − L terms in the equation for r andr become: where H is the Hubble expansion parameter. The averaged rates including the two processes in Eqs. (20) and (22) (assuming in the latter the Z s are in equilibrium) are computed in the appendix with the result: where the two terms inside the brackets correspond respectively to the ff and V V channels, and are valid for Note the different temperature dependence of the two contributions. The growth of the V V ↔ N N at high temperatures originates in the contribution of the longitudinal polarization of the V bosons when the temperature is below the scalar mass, M σ . For higher temperature, T ≥ T max , the contribution of the physical scalar σ has to be included, leading to a rate ∝ T . The new interactions do not modify the chemical potential dependent terms, nor the evolution equation for μ B/3−L α . The equations are therefore those in [31] with the additional B − L terms in Eq. (23).
To illustrate the effect of the B − L gauge interaction, we have considered the test point of Ref. [31] with masses for the heavy steriles m N 2,3 ∼ 0.8 GeV. Within the parameter space of successful leptogenesis, this point was chosen because it leads to charmed meson decays to heavy sterile neutrinos that could be observable in SHIP, and furthermore this measurement, in combination with input from neutrinoless double beta decay and CP violation in neutrino oscillations, could provide a quantitative prediction of the baryon asymmetry. Adding the B − L terms to the equations for r and r of [31], and solving them numerically (for details on the method see [31]) we obtain the curves in Fig. 5. The rates depend on m V so we choose m V = 1GeV. The evolution of

Dark matter
Now we want to discuss possible dark matter candidates in the B − L scenario without spoiling the BO mechanism, which as we have seen imposes a stringent upper bound on the gauge coupling, g B−L . We will be interested in the region where the V boson can decay to the lightest neutrino, ie m V ≥ 2m N 1 . The small value needed for g B−L suggests to consider the possibility of a freeze-in scenario [85,86], where the gauge boson does not reach thermalization, and neither does the lightest sterile neutrino, N 1 . The status of dark matter in a higher mass range through freeze-out has been recently updated in [64].
As it is well known, N 1 in the keV mass range is sufficiently long lived to provide a viable warm DM candidate [40,41]. The B − L model is as we will see a simple extension of the νMSM [9], which avoids the need of huge lepton asymmetries to evade X-ray bounds. In our scenario the keV state is produced from the decay V → N 1 N 1 , while the lifetime of N 1 , relevant in X-ray bounds, is controlled also by mixing, which can be sufficiently small, in a technically natural way, as in the νMSM scenario (provided the lightest neutrino mass is small enough). A similar scenario for DM has been studied in [67]. We now quantify the parameter space for successful DM and leptogenesis in this scenario.
We assume that the abundance of V and N 1 is zero at a temperature below the EW phase transition where all the remaining particles in the model are in thermal equilibrium. All fermions in the model couple to the V and therefore its production is dominated by the inverse decay process: f f → V . The kinetic equation describing the production of V is the following: where f i ( p) are the distribution function of the particle, i, with momentum p, and is the number density, with g i the number of spin degrees of freedom. g f = 4 for a Dirac fermion, g N = 2 for a Majorana fermion and g V = 3 for a massive gauge boson. M is the amplitude for the decay V → f f at tree level.
The sum over f is over all fermions, but we can safely neglect the contribution of the N 1 and also those that are nonrelativistic. We can also neglect the Pauli-blocking and stimulated emission effects ( f i ± 1 ∼ ±1) and approximate the distribution function in equilibrium for fermions and bosons by the Maxwell-Boltzmann, f i ( p i ) = f eq ( p i ) = e −E i /T . Taking into account the relation and the principle of detailed balance the equation can be simplified tȯ As long as f V f eq , the first term on the right-hand side can be neglected and the equation simplifies further to: where K 1 is the first modified Bessel Function of the 2nd kind. The decay width in the V rest frame is given by where N C = 3(1) and Q f = 1/3(−1) for quarks(leptons). As usual we define the yield of particle i as where s is the entropy density and we can assume g * s g * . We also consider the averaged momentum approximation which amounts to assuming that f V has the same momentum dependence as f eq . Changing variable from time to temperature, the final evolution equation for Y V reads: where n eq V = 3 2π 2 m 2 V T K 2 (m V /T ). The production of N 1 is dominated by the decay V → N 1 N 1 . There is also the contribution via mixing with the active neutrinos but this is negligible for mixings that evade present X ray bounds. Neglecting the inverse processes, the evolution equation for n 1 iṡ and in terms of the yield where It is straightforward to solve these equations. In Fig. 6 we show the yields of V and N as function of the inverse temperature for m V = 10 MeV, m N 1 = 10 keV and g B−L = 10 −11.4 . The resulting abundance of N 1 is where s 0 = 2889.2 cm −3 is the entropy today and ρ c = 1.0510 −5 h 2 GeV cm −3 is the critical density. The evolution of Ω N 1 h 2 is shown in Fig. 7 for two values of m V and a fixed value of g B−L . Requiring that Ω N 1 h 2 equals the full DM contribution of Ω DM h 2 0.12 implies a relation between m V and g B−L as shown in the curves of Fig. 2. The values of g B−L corresponding to the right dark matter relic abundance do not affect leptogenesis and lie far below the actual collider limits. Nevertheless some regions of the parameter space are interestingly excluded from supernova and BBN observations.
A final comment concerns the comparison of our calculation of the DM abundance and that in Ref. [67]. In this reference only the evolution of the N 1 is considered, and the collision term corresponds to the scattering process ff → N 1 N 1 , where the narrow width approximation is assumed. We believe this method is only equivalent to ours when all f,f and V distributions are the equilibrium ones, but this is not the case here. In the region they can be compared our results are roughly a factor three smaller than those in [67].

Couplings
According to the previous calculation, the relic DM abundance requires a very small B − L coupling. In order to and therefore φ ∼ 2 · 10 8 GeV.
In order to get the N 1 and N 2,3 in the target range of keV and 1-100 GeV respectively, small and hierarchical h N i couplings are needed: for the heavy sterile neutrinos involved in BO leptogenesis and for the dark matter candidate. Note that the required h N 's couplings are in the same ballpark as the yukawa couplings. The gauged B − L model works nicely to explain neutrino masses, the baryon asymmetry and dark matter. Unfortunately it also requires a very small g B−L which will be very hard to test experimentally. An alternative might be to consider a flavoured U (1), for example L μ − L τ , that might be compatible with a larger g B−L , provided the assignment of charges to the singlet states ensures that not all of them reach thermalization via the flavoured gauge interaction before t EW .

Axion and neutrinos
As a second example we consider an extension of Eq. (1) with a scalar doublet and a scalar singlet. This model is also an extension of the invisible axion model [87] with sterile neutrinos, that was first considered in [46], providing a connection between the Peccei-Quinn (PQ) symmetry breaking scale and the seesaw scale of the neutrino masses. The model contains two scalar doublets, Φ i , and one singlet, φ.
A U (1) P Q global symmetry exists if the two Higgs doublets couple separately to the up and down quarks and leptons so that the Yukawa Lagrangian takes the form: leading naturally to type II two-Higgs-doublet models without FCNC [88,89]. The most general scalar potential of the model compatible with a global U (1) P Q is the following The couplings in this potential can be chosen such that φ gets an expectation value, U (1) P Q is then spontaneously broken and a Nambu-Goldstone boson appears, the QCD axion. Furthermore the Majorana singlets N get a mass. Expanding around the right vacuum, the field can be writen as where σ is a massive field, while a is the axion. Therefore after symmetry breaking we obtain an interaction term between sterile neutrinos and axions The breaking scale f a must be much larger than the vacuum expectation values of the doublets, v 1,2 , so that the axion can evade the stringent bounds from rare meson decays and supernova cooling, which sets a stringent lower bound f a ≥ 4 · 10 8 GeV [90].
The mass of the axion is induced by the QCD anomaly in the sub-eV range: (48) where z = m u /m d . For f a ≥ 4 · 10 8 GeV, we have It is well known that the invisible axion is a viable cold DM candidate, through the misalignment mechanism [87,91,92] (for recent reviews see [93,94]). The DM energy density is given by where θ 0 is the misalignment angle. The constraints on f a depend on whether the breaking of the PQ symmetry happens before or after inflation; in the latter case the misalignment angle can be averaged over many patches so Ω a ≤ Ω DM implies f a 1.2 · 10 11 GeV, (52) with the equality reproducing the observed cold dark matter energy density Ω CDM h 2 ∼ 0.12. This correspond to the solid line in Fig. 8. If the PQ symmetry is broken before inflation, θ 0 is a free parameter and the value of f a to account for DM is inversely proportional to θ 2 0 . The axion can also manifest itself as dark radiation [95], given that it is also thermally produced [96]. This population of hot axions contributes to the effective number of relativistic species, but the size of this contribution is currently well within the observational bounds [97].
In this model the VEV of the scalar singlet gives a Majorana mass to the sterile neutrinos: So, if we want a mass in the electroweak range, O(1 − 10 2 ) GeV and f a ∈ [10 8 , 10 11 ] GeV, we need the coupling h N to be in the range: The hierarchy between f a and the electroweak scale requires that some couplings in the scalar potential in Eq. (44) (k, λ 1φ , λ 2φ ) are very small. Even if not very appealing theoretically, these small numbers are technically natural as already pointed out in [98], where the authors studied the same model with very heavy sterile neutrinos. A relevant question is that of naturalness or fine-tunning of the Higgs mass in this model. In [98], this issue was studied in the context of high-scale thermal leptogenesis, and it was Fig. 8 The light purple region is forbidden by supernova cooling constraints (left band) and the red one by axions over closing the universe (right band). Axions can explain the DM relic density in the vertical line, if inflation happens before PQ symmetry breaking. Within the remaining parameter space, in the region between the two dashed lines successful leptogenesis via oscillations is possible concluded that stability imposes relevant constraints. In particular, a relatively small v 2 30 GeV is necessary to ensure viable leptogenesis for lower m N = 10 5 − 10 6 GeV so that yukawa's are small enough, y ≤ 10 −4 , and do not induce unnaturally large corrections to the Higgs mass. In our case, the yukawa couplings, Eq. (5), are too small to give large corrections to the Higgs mass , so no additional constraint needs to be imposed on v 2 . As a consequence other invisible axion models, such as the KSVZ [99,100], would also work in the context of low-scale m N , but leads to tension with stability bounds in the high-scale version [101,102].

Baryon asymmetry
The possibility to generate the baryon asymmetry in this model a la Fukugita-Yanagida for very heavy neutrinos m N ∼ f a was recognized in the original proposal [46] and further elaborated in [98]. We want to point out here that for much smaller values of h N , the BO mechanism could also work successfully. As explained above the crucial point is whether the new interactions of the sterile states in this model are fast enough to equilibrate all the sterile neutrinos before EW phase transition. The leading order process we have to consider is the decay of the scalar into two sterile neutrinos σ → N N, exactly as we considered in the previous section. The limit of M σ ∼ f a ≥ 2 × 10 5 − 5 × 10 6 GeV derived in Eq. (15) also applies here, which is safely satisfied given the supernova cooling bounds. At second order, we must also consider the new annihilation process of sterile neutrinos to axions N N ↔ aa as shown in Fig. 9. The rate of this process at high tempera- The condition (for m N ∈ [1, 10 2 ] GeV), safely within the targeted range. Fig. 8 shows the region on the ( f a , h N ) plane for which successful baryogenesis through the BO mechanism and DM can work in this model. Even if the necessary condition for BO leptogenesis is met for f a ≥ 10 8 GeV, the presence of the extra degrees of freedom, the axion, the heavy scalar and the second doublet could modify quantitatively the baryon asymmetry. For example, the presence of two scalar doublets could modify the scattering rates of the sterile neutrinos considered in the BO scenario, where the main contributions [35,50,52,103] are: -2 ↔ 2 scatterings on top quarks via higgs exchange -2 ↔ 2 scatterings on gauge bosons -1 ↔ 2 decays or inverse decays including resumed softgauge interactions Sterile neutrinos are coupled to the same Higgs doublet that also couples to the top quarks; in this case nothing changes with respect to the usual calculation, in which the reactions with top quarks are mediated by Φ 1 . However, an alternative model, with a different U (1) P Q charge assignment, is also possible, in which the sterile neutrinos couple to φ 2 and not φ 1 , as done in [98]. In this case top quark scattering does not contribute to sterile neutrino production at tree level, but the baryon asymmetry is not expected to change significantly, since the scattering rate on gauge bosons and the 1 → 2 processes are equally important [50,52]. The process σ → N N is not foreseen to be relevant for M σ ∼ f a , since the scalar is long decoupled when the generation of the asymmetry starts, while the new process N N ↔ aa is expected to be very small according to the above estimates. It could nevertheless be interesting to look for possible corners of parameter space where the differences with respect to the minimal model is not negligible since this could provide a testing ground for the axion sector of the model.

Majoron model
In between the two models described in 3-4, there is the possibility of having a global U (1) spontaneously broken, which is not related to the strong CP problem and we call it lepton number. This is of course the well-known singlet majoron model [47,48,104]. We assume the sterile neutrinos carry lepton number, L N = 1, but Majorana masses are forbidden and replaced by a yukawa interaction as in the B − L model: where Φ is the standard model Higgs doublet, while φ is a complex scalar which carries lepton number L φ = −2.
Then, the complex scalar acquires a VEV and the U (1) L is spontaneously broken giving rise to the right-handed Majorana mass matrix and leading to a Goldstone boson η, the majoron. Consequently the Lagrangian will induce the new scattering processes for neutrinos depicted in Fig. 9. As usual we have to ensure that at least one sterile neutrino does not equilibrate before T EW (see [105,106] for a recent discussion in the standard high-scale or resonant leptogenesis). As in the previous cases we have to consider the decay σ → N i N i and the annihilation into majorons, Fig. 9. The former gives the strongest constraint, as in Eq. (15): These lower bounds for m N = 1 and 100 GeV are shown by the horizontal lines in Fig. 10.

Dark matter
There are two candidates in this model for dark matter that we consider in turn: the Majoron and the lightest sterile neutrino, N 1 .

Majoron
In this model a natural candidate for dark matter is the majoron itself, but it has to acquire a mass, therefore becoming a pseudo Nambu-Goldstone boson (pNGB). One possibility is to appeal to gravitational effects [107,108]. However, the contribution to the mass from gravitational instantons is estimated to be [109][110][111] and therefore extremely tiny, unless f is close to the Planck scale. Another alternative is to consider a flavoured U (1) X and soft symmetry breaking terms in the form of yukawa couplings [112,113]. This possibility has been studied in detail in Ref. [113]. It has been shown that the majoron can be the main component of dark matter for sterile neutrino masses m N ≥ 10 5 GeV, while for masses in the range we are interested in (m N ∼ 1 − −100 GeV) neither thermal production via freeze-out nor via freeze-in works.
The possibility to produce it via vacuum misalignment, analogous to the one which produces the axion relic density has also been discussed in [113]. It was shown to give a negligible contribution compared to the thermal one, because the majoron gets a temperature dependent mass at early times. Even if the mass of the majoron is significantly smaller in our situation, with lighter m N , we find the same result, ie. that only a small fraction of the DM can be produced via misalignmet.
No matter what the production mechanism is if the majoron constitutes the dark matter, there are constraints from the requirement that the majoron be stable on a cosmological timescale and its decay to the light neutrinos should not spoil the CMB anisotropy spectrum [114,115]. This gives constraints on the mass m η and the symmetry breaking scale f , as showed in Fig.10.
As in the axion case there are additional constraints from supernova cooling [116], but they are much weaker and give an upper bound much lower than the range shown in Fig. 10. In the unconstrained region in Fig. 10, ARS leptogenesis and majoron DM could in principle work provided the mechanism to generate the majoron mass does not involve further interactions of the sterile neutrinos.

Sterile neutrino
We want now to consider the sterile neutrino as a dark matter candidate, a possibility already explored in [117,118] in a model with a real scalar field and therefore no Majoron. In our case the presence of the Majoron could make the sterile neutrino unstable, given that it would decay through the channel Thus one has to assume that the Majoron has a larger mass so that this decay is kinematically forbidden. As in the B−L case, BO leptogenesis is driven by the other two heavier neutrino states N 2,3 , while N 1 can be produced through freeze-in from σ → N 1 N 1 decay. Assuming σ is in thermal equilibrium with the bath the Boltzmann equation describing the evolution of the N 1 density: where we have neglected Pauli blocking, and the inverse processes.
Following the standard procedure we end with the contribution to the abundance: Using and requiring that Ω N 1 h 2 matches the observed DM we find The mass of the DM candidate is related to the coupling which regulates the freeze-in process through the VEV of φ therefore if M σ ∼ φ , the coupling needs to be as shown in Fig. 11. If m N 1 is the keV range so that it can satisfy cosmological and astrophysical constraints, the scale of the VEV should be φ ∼ TeV.
As a consequence we see that if one couples the B − L scalar also to the heavier neutrinos, the interactions are too fast and the BO mechanism cannot work since the bound in Eq. (59) is not satisfied. Alternatively, if we assume Eq. (59), then Such a massive neutrino would need extremely small activesterile mixing angle (collectively labeled θ ) to be sufficiently long-lived. The strongest bound from X-rays [119] give while the soft gamma ray bound [120] gives In conclusion, either we consider a global symmetry with different family charges, more concretely a (B − L) 1 , or we need to require extremely tiny yukawa coupling for the DM sterile neutrino, making this last model theoretically unappealing.

Conclusion
The extension of the Standard Model with three heavy majorana singlets at the weak scale can explain neutrino masses and also account for the baryon asymmetry in the Universe via the oscillation mechanism [8,9]. This scenario could be testable in future experiments. Unfortunately the simplest model cannot easily accommodate dark matter. In the νMSM [9], one of the three heavy states is in the keV range and provides a candidate for dark matter, but it requires huge lepton asymmetries that cannot be naturally achieved in the minimal setup.
In this paper we have explored three extensions of the minimal scenario that can accommodate dark matter without spoiling baryogenesis. This is non trivial because new interactions of the heavy singlets can disrupt the necessary out-of-equilibrium condition which is mandatory to generate a lepton asymmetry. We have shown that a extension of the minimal model with a U (1) B−L gauge interaction can achieve this goal. The two heavier majorana fermions take part in the generation of the baryon asymmetry, while the lighest one in the keV range, N 1 , is the dark matter. In contrast with the νMSM the production of the dark matter is not via mixing, but it is dominated by the B − L gauge boson decay. The mixing is however what controls the decay of the N 1 and can be made sufficiently small to avoid the stringent X-ray constraints. The correct DM abundance is achieved for very small B − L gauge couplings, g B−L 10 −8 , which are safely small not to disturb the baryon asymmetry, which remains the same as in the minimal model. Such tiny couplings are far below the reach of colliders. Supernova and BBN provide the most stringent constraints in the relevant region of parameter space, while future searches in SHIP might have a chance to touch on it.
We have also considered an extension involving an invisible axion sector with an extra scalar doublet and a complex singlet. The heavy majorana singlets get their mass from the PQ breaking scale [46]. DM is in the form of cold axions, from the misalignment mechanism and as is well known, the right relic abundance can be achieved for a large value of the PQ breaking scale, f a 10 11 GeV. We have shown that such large scale is compatible with having the heavy neutrinos in the 1-100 GeV scale, and ARS leptogenesis. Finally we have considered the singlet majoron extension of the minimal model, with a global U (1) B−L , that contains two potential DM candidates, the majoron or the lightest heavy neutrino, N 1 . Unperturbed ARS baryogenesis requires a relatively high B − L breaking scale, f 10 6 GeV. Majoron DM requires exotic production scenarios, while neutrino DM works for masses around MeV, which requires extremely small mixings to make it sufficiently long-lived, or alternatively a less theoretically appealing possibity, where the scalar couples to only one sterile neutrino, while the other two have tree level masses or couple to a different scalar with a larger VEV. As a general rule, adding new interactions that affect the heavy Majorana singlets modifies leptogenesis in the minimal model and viable extensions that can explain DM are likely to involve the freeze-in mechanism as in the examples above.
Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: We didn't use real data in this paper, which is mostly a theoretical one. The only data we analyzed were taken from other works and are shown in the constraints plots.] Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .

A Computation of momentum averaged rates in gauged B − L
In this appendix we give some details on the computation of the momentum averaged rates in Eq. (24). The amplitude for f f → N N for vanishing masses is given by while that for V V → N N, in the limit M 2 Defining the Bose-Einstein and Fermi-Dirac distributions and the variables q ± ≡ 1 2 (q 0 ± |q|), where q = p 1 + p 2 . We express all momenta in units of temperature T .
Following the procedure of Ref. [50] the rate ff → N N can be writen as R (1) (k) = A 4(2π) 3 k 0 (r 1 (k) + r 2 (k) + r 3 (k)) , with with and The rate V V → N N can be written as with and A ≡ g 4 The total rate is and the averaged rates γ (0) V and γ (1) V are found to be: where and