On neutrinoless double beta decay in the minimal left-right symmetric model

We analyze the general phenomenology of neutrinoless double beta decay in the minimal left-right symmetric model. We study under which conditions a New Physics dominated neutrinoless double beta decay signal can be expected in the future experiments. We show that the correlation among the different contributions to the process, which arises from the neutrino mass generation mechanism, can play a crucial role. We have found that, if no fine tuned cancellation is involved in the light active neutrino contribution, a New Physics signal can be expected mainly from the $W_R-W_R$ channel. An interesting exception is the $W_L-W_R$ channel which can give a dominant contribution to the process if the right-handed neutrino spectrum is hierarchical with $M_1\lesssim$ MeV and $M_2,M_3\gtrsim$ GeV. We also discuss if a New Physics signal in neutrinoless double beta decay experiments is compatible with the existence of a successful Dark Matter candidate in the left-right symmetric models. It turns out that, although it is not a generic feature of the theory, it is still possible to accommodate such a signal with a KeV sterile neutrino as Dark matter.


I. INTRODUCTION
The recent LHC results [1,2] seem to indicate that the Higgs mechanism, with the Higgs mass around 125 GeV, is the responsible for the mass generation of the Standard Model (SM) particles. However, the origin of light neutrino masses, for the existence of which we have compelling evidences from neutrino oscillation experiments, still remains unknown. It is true that the light neutrino masses could also be generated through the Higgs mechanism in a minimally extended SM which includes sterile (right-handed) neutrino fields as SU (2) L singlets and in which the total lepton number is conserved. However, their smallness in comparison with the charged lepton and quark masses calls for a different explanation. In this context, extensions of the SM required to explain the origin of neutrino masses, and compatible with the latest LHC data, arise as quite suggestive models of New Physics (NP).
Among those we find the celebrated seesaw models [3][4][5][6], which can give us, in addition, the key to explain the matter-antimatter asymmetry of the universe through Leptogenesis [7].
Most of those models predict that neutrinos are Majorana particles, something which can be tested in lepton number violating processes such as the neutrinoless double beta (0νββ) decay. The 0νββ decay experiments are the most promising ones in this context but they suffer a serious drawback: the NP contribution to the process is usually short range and thus typically very suppressed compared to that of the light neutrinos. Thanks to the future 0νββ experiments [8][9][10][11][12][13][14], in combination with the complementary information coming from neutrino oscillation experiments and cosmology, we might be able to discover the Majorana nature of neutrinos, but not easily which is the mechanism responsible for the neutrino mass generation [15,16]. In this context, the correlations between the standard light neutrino and NP contribution to the 0νββ decay are crucial, as shown in the case of the type-I [3][4][5][6], type-II [17][18][19][20][21] and type-III [22] seesaw models in Refs. [23,24]. The generation of light neutrino masses in a particular model usually induces important correlations between the different contributions to the 0νββ decay, which should always be considered in a model dependent analysis, helping to understand which type of NP can be feasibly tested in the experiments.
In this work we will focus on the 0νββ decay phenomenology of the minimal left-right symmetric model (MLRSM) [6,21,[25][26][27]. The left-right symmetric models have been widely studied in the literature since, among other features, they provide a natural expla-nation for the smallness of the neutrino masses (some recent analysis in the context of the 0νββ decay can be found in Refs. [28][29][30][31][32][33][34][35]). In our analysis we will assume that no accidental cancellation occurs in the light neutrino mediated W L − W L channel, which involves the exchange of two W L . We will distinguish three regions of the parameter space depending on the mass of the right-handed (RH) neutrinos. First, we will show that if the right-handed (RH) neutrinos are heavier than the 0νββ decay scale (∼ 100 MeV), the 0νββ decay rate is dominated by light neutrino exchange channels with the exception of the channel in which two W R are exchanged (W R − W R channel) mediated by heavy neutrinos 1 . One of the light neutrino mediated channels involves the exchange of one W L and one W R (W L − W R channel); however, it turns out that a NP dominant contribution can come mainly from the W R − W R channel. Secondly, we will study the region of the parameter space where the RH neutrinos are lighter than the 0νββ scale. We have found that in this case the W L − W R contribution cancels out while a NP signal can still be expected from the W R − W R channel.
In this region, the RH neutrinos can give a relevant contribution through the W L − W L channel, as opposed to the type-I seesaw case where the total W L − W L contribution is very suppressed. Finally, we will investigate a mixed scenario with RH neutrinos in both regions below and above the 0νββ decay scale. We have found that this is the only scenario in which the W L − W R channel turns out to be relevant and can be responsible of a future signal (if no cancellation in the W L − W L channel is invoked). In all the cases we will show for which part of the parameter space a NP signal in future 0νββ decay experiments can be expected.
Moreover, we will also analyze if such a signal can be compatible with the existence of a This work is organized as follows. In Section 2 we briefly describe the MLRSM, focusing on the relations among the parameters of the model induced by the neutrino mass generation.
In Section 3 we analyze the neutrinoless double beta decay phenomenology in the MLRSM, studying in particular for which part of the parameter space a 0νββ decay signal coming mainly from NP contributions can be possible. Section 4 is devoted to the analysis of complementary constraints coming mainly from charged LFV experiments and the stability 1 If other contributions coming from different models are not involved.
3 of the light neutrino masses under 1-loop corrections. In Section 5 we study if a successful DM candidate and a NP signal in the future 0νββ decay experiments can be compatible in the MLRSM. Finally, we conclude in Section 6.

II. MINIMAL LEFT-RIGHT SYMMETRIC MODEL AND NEUTRINO MASSES
The Lagrangian of the MLRSM respects an enlarged gauge symmetry SU (3) c ⊗SU (2) L ⊗ SU (2) R ⊗ U (1) B−L plus a discrete left-right symmetry which leads to equal SU (2) L and SU (2) R gauge couplings (g L = g R = g). We are not going into the details of the model since it has been widely studied in the literature (for a recent complete analysis regarding the associated lepton number violating effects, see for instance Ref. [29,33]), but only recall the most relevant features for our analysis. The scalar sector is also augmented by the addition of two scalar triplets (∆ L and ∆ R ) and a bi-doublet scalar under SU (2) L ⊗ SU (2) R , which spontaneously break the electroweak symmetry when they develop vacuum expectation values (vevs).
In this section we will derive the relations which will be used in the phenomenological analysis of the 0νββ decay. Since they come from the neutrino mass generation, let us recall how the complete neutrino mass matrix looks like after the electroweak symmetry breaking: where m i are the light neutrino masses and M i the heavy ones. Notice that in this model the Majorana mass term for the heavy neutrinos is generated dynamically when ∆ R takes a vev (M R = Y ∆ R v R ), while the Majorana mass term M L for the left-handed (LH) neutrinos is generated analogously through the ∆ L vev (M L = Y ∆ L v L ). The neutrino mass matrix is diagonalized as shown above by a 6 × 6 unitary matrix U , through the following rotation between the neutrino flavor and mass eigenstates denoted by α, β = e, µ, τ and i, k = 1, 2, 3, respectively, The diagonalization of the complete neutrino mass matrix presented in Eq.
(1) provides the following useful relations On the other hand, taking into account that the active LH block of U ,Ũ , is unitary to a very good approximation (at least up to the percent level [36]), the complete neutrino mixing matrix can be expanded as, where θ is a 3 × 3 matrix which characterizes the small mixing between the active LH and the heavy RH neutrinos, U pmns is the PMNS matrix and V is a 3 × 3 unitary matrix. From Eqs. (3)(4)(5), we have The discrete (charge conjugation) LR symmetry gives us the following relation between the Yukawa couplings of the triplets: where we have employed the present bounds on v L and v R , namely v L 7 GeV [37] and TeV [38][39][40]). Therefore, the O (M L θ * ) and O M L M R θ 2 can be safely neglected in Eq. (8) and Eq. (9) respectively, and 2 Another option is to consider instead a discrete parity symmetry leading to a similar relation: Of course, θ plays a fundamental role at the phenomenological level since it basically describes the mixing between the active LH neutrinos and the RH ones. It would be very interesting thus to find a useful parametrization of θ as a function of the light neutrino parameters, light neutrino masses and the angles/phases of the PMNS matrix, and the rest of the independent parameters of the model associated with the RH neutrino sector. In principle, an analogous parametrization to the Casas-Ibarra one [41] would be a good candidate [42]. However, the presence of M L in Eq. (7) and the fact that the matrix V is in this case physical, contrary to the type-I seesaw model, makes that parameterization less transparent and more involved than expected. On the other hand, the discrete (charge conjugation) LR symmetry leads to the following constraint and thus, Eq. (11) becomes Plugging this relation and Eq. (9) into Eq. (7), we obtain, and finally with Therefore, θ is completely determined as a function of the light and heavy neutrino masses, m and M , the PMNS matrix, U pmns , v L /v R , and the unitary matrix V [31]. Notice that if this expression is used to obtain θ with the PMNS mixing angles and the solar and atmospheric mass-squared differences as input parameters, we ensure that the model is consistent with the light neutrino mass and mixing pattern measured in neutrino oscillation experiments.

III. NEUTRINOLESS DOUBLE BETA DECAY
In our study of the 0νββ decay in the MLRSM, we will pay special attention to the correlation among all the contributions to the process and, in particular, the connection 6 with the light neutrino masses. We shall see that the correlation between the different contributions and the experimental bounds on the parameters will allow us to safely neglect some of the NP contributions.
As we have already mentioned, we will not analyze the scenario in which a cancellation occurs within the standard light neutrino contribution, which would naively leave the NP channels as the leading contributions [43,44]. Of course, this cancellation can be due to the presence of an extra symmetry added to the model, such as the lepton number which is approximately conserved in the so called inverse or direct seesaw models [45][46][47][48][49]. The problem in this scenario is that, in order for the NP contributions to be measurable, a significant violation of lepton number should be introduced through the NP sector which may not have an impact on the light neutrino masses at tree level but arises naturally at one-loop level, as shown in Ref. [23] in the context of the seesaw models. This makes it very difficult to have a significant contribution from NP channels since the 1-loop correction to the light neutrino masses tends to dominate in the 0νββ decay rate.
We will distinguish three different regions according to the associated 0νββ decay phenomenology: (i) when the RH neutrinos are much heavier than the 0νββ decay scale ( p ≈ 100 MeV), which means heavier than approximately 1 GeV; (ii) when the RH neutrinos are much lighter than the 0νββ decay scale (below 1 MeV); (iii) when the RH neutrinos are in both regions, (i) and (ii).
In the analysis below we have reasonably estimated the NMEs corresponding to some of the channels under study. This is accurate enough for our purposes but, although the associated NMEs errors are still large, in order to be more precise, full calculation of all the NMEs should be considered.

A. Heavy regime
The various contributions to the 0νββ transition rate in this model are described by the Feynman diagrams shown in Fig. 1. We will start describing them one by one in order to show that the contributions of the diagrams in which a heavy fermion (or scalar) is exchanged are subdominant with respect to those of the light neutrino exchange, the only exception being the W R − W R channel.
where M 0νββ (m i ) are the associated nuclear matrix elements (NMEs) following the notation of Ref. [20], where the NMEs were computed as a function of the mass of the neutrino mediating the process for different nuclei. Notice that in this notation the NMEs include the dependence on the propagator. The NMEs corresponding to the light neutrino exchange are independent of the neutrino masses, and then with Eq. (3), the above amplitude can be rewritten as Taking into account that M 0νββ (M i )/M 0νββ (0) 1 [46], the contribution due to the heavy neutrino exchange can be safely neglected. Using again Eq. (3), one obtains, is given by: where ξ is the W L − W R mixing angle. The present experimental bound is given by ξ 10 −2 [38], however, in the minimal left-right symmetric model there is a stronger theoretical upper bound given by M 2 [40]. 3 In the case of the W L − W L contribution, ξ can be safely neglected. M 0νββ are the associated nuclear matrix elements (NMEs) following the notation of Ref. [24], where the NMEs were computed as a function of the mass of the neutrino mediating the process for different nuclei. Notice that in this notation the NMEs include the dependence on the propagator. The NMEs corresponding to the light neutrino exchange are independent 3 Notice that ξ can only saturate this bound if and only if the two vev's of the Higgs doublets are of the same order.
8 of the neutrino masses, and then with Eq. (3), the above amplitude can be rewritten as Taking into account that M 0νββ (M i )/M 0νββ (0) 1 [50], the contribution due to the heavy neutrino exchange can be safely neglected. Using again Eq. (3), one obtains, which is the standard light neutrino contribution.
Using Eq. (6) in the above equation we obtain: Clearly, the second term can be neglected in comparison with the standard contribution due to the double suppression coming from and the active-heavy mixing, at least |θ αi | 2 10 −2 [36,52,53]. The first term, however, can not be neglected, i.e, • W L − W R channel. For the diagram in the top right of Fig. 1, in which W L and W R are exchanged, the amplitude is given by 9 where η ≈ 10 −2 [54][55][56][57][58][59] 4 . Taking into account that U is unitary, we have: and since M 0νββ which implies that the light neutrino mediated contribution of the W L − W R channel is again dominant over the heavy neutrino exchange.
• The amplitude corresponding to the scalar triplet ∆ L exchange (bottom right in Fig. 1 with W L and ∆ L ) is suppressed with the factor with respect to the standard contribution given in Eq. (17). The suppression factor is 1 if no fine tuned cancellation between the two terms in the light neutrino contribution is invoked, i.e., the contribution of this channel is negligible.
For the corresponding "right-handed" version of the diagram the situation is slightly different and the suppression factor now reads It seems that for small enough values of this contribution could be larger than the standard one. However, it is not very easy to achieve a measurable ∆ R contribution, at the reach of the sensitivity of the next-to-next of 0νββ decay experiments (m ββ ∼ 10 −2 eV). Indeed, the corresponding amplitude is given by where we have used The only possibility of having a phenomenologically relevant contribution is to saturate the [40] and M W R ≈ v R g/ √ 2) having at the same time Y ∆ ρ, which is not very feasible since a small value of ρ would render ∆ R too light, contradicting the experimental bound, m ∆ R > 320 GeV [60]. We will thus neglect this contribution.
We have shown that only the contributions coming from the light neutrino exchange can have a significant impact in the 0νββ decay rate, with the exception of the channel mediated by two W R gauge bosons in which the heavy neutrino exchange dominates. In summary, the phenomenologically relevant contributions to the 0νββ decay rate can be recast as where we have made use of Eq. (6) and c LL , c LR and c RR are coefficients which take into account the different chirality of the outgoing electrons. At this point we can make an estimation of the NMEs associated to the heavy neutrino exchange, M 0νββ (M i ), to understand how relevant the remaining NP contributions are. The effective mass becomes, where we have neglected the suppressed interference terms between the different chirality contributions [61].
In the rest of this section, we will first study the bounds that can be extracted from the 0νββ decay experiments if one assumes that the three contributions listed in Eq. (27) are completely independent. After that, we will study the region of the parameter space in which a NP signal in the future 0νββ decay experiments can be expected when the correlations among the different contributions are not ignored. 5 We refer readers to Ref. [21] for more details In Fig. 2 we show the constraints on v R (recall that M W R = gv R / √ 2) and the mixing between ν eL and the lightest heavy neutrino, |θ e1 |, extracted from 0νββ decay experiments when only the contribution from the W L −W R channel (second term in Eq. (27)) is taken into account. In the left panel the mixing ξ saturates the theoretical bound (ξ = M 2 W L /M 2 W R ) while in the right panel ξ is neglected. The shaded region is ruled out by the present constraint, |m ββ | < 0.38 eV [9], while the region between the red dashed lines corresponds to the sensitivity of the next-to-next generation of experiments, 10 −2 eV < |m ββ | < 0.38 eV. In Fig. 4 we show the sensitivity of the next-to-next generation of 0νββ decay experiments (10 −2 eV < |m ββ | < 0.38 eV) to the parameters of the model by including all the relevant 13 contributions and requiring the NP contribution to the 0νββ decay rate (second and third term in Eq. (27)) to be at least 10 times larger than the standard contribution (first term of Eq. (27)). The allowed region is projected onto the v R -(Y ∆ ) ee plane (left panel) and the v R -|θ e1 | plane (right panel). The mixing has been neglected in the upper panels while in the lower panels is fixed to its maximum value ξ = M 2 W L /M 2 W R . The experimental constraints on the W R mass [40] and the active-heavy mixing [36,52,53] have been also included.
We have assumed that the heavy neutrino spectrum is hierarchical ( if the trilinear coupling and the mixing are small enough, 3 · 10 −6 (Y ∆ ) ee 3 · 10 −2 (3 · 10 −6 (Y ∆ ) ee 8 · 10 −2 ) and |θ e1 | 2 · 10 −5 respectively, for ξ = 0 (maximal mixing . This corresponds to a range of heavy neutrino masses from GeV to TeV. Comparing the upper and lower panels we can conclude that including the mixing in the analysis has some impact in the results but it is not very significant.
Comparing Fig. 4 with Fig. 2, where only the W L − W R contribution is included, we see

B. Light regime
If the RH neutrinos are lighter than the 0νββ decay scale, O (100) MeV, the picture significantly changes with respect to the heavy scenario studied above. Eqs. (16), (19) and (21) remain valid but the NMEs associated with the "heavy" neutrino exchange are not suppressed compared to the light neutrino mediated ones. In fact, for M i < 1 MeV, we have in a very good approximation M 0νββ (0) = M 0νββ (M i ). This yields a cancellation within the second term of Eqs. (16) and the corresponding amplitude is then given by while Eq. (19) becomes where we have used Eq. (9). The W L − W R contribution, given by Eq. (21), vanishes due to the unitarity of the 6 × 6 neutrino mixing matrix U . Therefore, from Eqs. (28)-(29), the expression for the effective mass m ββ when the RH neutrinos are lighter than 1 MeV becomes: We conclude that in this regime the 0νββ decay can be completely attributed to the 6 Note that the mass of M W L mainly comes from the SM Higgs vev, not from the ∆ L vev, v L . As a consequence, v L could be very small. Furthermore, a small v L is in better agreement with the ρ(≡ M 2 W L /M 2 Z L cos 2 θ W ) parameter constraints.
In the analysis all the relevant contributions have been simultaneously included. The bounds on M W R [40] have been also included.
Notice that, contrary to the type-I seesaw limit (v L → 0) [24], in this regime A LL does not vanish and the RH neutrinos (second term in the equation above) can contribute to the process. Nevertheless, in this work we will focus on the W R − W R and W L − W R channels.
A dominant NP contribution from W L − W L channel mediated by the RH neutrinos will be investigated elsewhere in more detail.
The future sensitivity of the next-to-next generation of 0νββ decay experiments (10 −2 eV < |m ββ | < 0.54 eV) to the parameters of the model when the W R − W R contribution (second term of Eq. (30)) is at least 10 times larger than that from the W L − W L channel (first term of Eq. (30)) is given in Fig. 5. The allowed region of the parameter space is The bounds on the W R mass [40] have been also included. Finally, comparing the dashed and solid contours we can conclude that the impact of the mixing ξ is not very significant in this region of the parameter space.

C. Mixed Scenario
There is an alternative scenario that has not been studied in the previous sections and consists of the existence of RH neutrinos in both regimes below and above the 0νββ decay scale. In this section we will focus on the particular case in which one of the RH neutrinos is lighter than 1 MeV and the other two are heavier than 1 GeV , i.e., M 1 < 1 MeV and M 2 , M 3 > 1 GeV, but the phenomenology remains similar if two RH neutrinos are lighter than 1 MeV.
As it occurs in the previous section, Eqs. (16), (19) and (21) are also correct in this regime, but only the NMEs associated with the N 2 and N 3 exchange are suppressed compared to the light neutrino mediated ones. The NMEs associated with N 1 satisfy M 0νββ (0) = M 0νββ (M 1 ). As a consequence, in this regime Eqs. (16) and (21) read and Eq. (19) becomes where again we have used Eq. (9) and the fact that M 0νββ (M i )/M 0νββ (0) 1 for i = 2, 3.
Therefore, in this scenario the effective mass m ββ is given by: Contrary to the light regime, in this scenario the W L − W R contribution may be significant.
Notice that if V = I, the W L − W R contribution cancels out, which means that the RH neutrino mixing V is very relevant in this region.
In order to be consistent with the rest of this work, associated with the absence of cancellation in the light active neutrino contribution, we will focus in this section on the limit v L v R → 0, neglecting the first term of the W L − W L contribution. As expected, we have checked that if that term is switched on in the analysis, cancellations between the two terms of the W L − W L contribution can perfectly occur for some part of the parameter space leading to better future sensitivities to the W L − W R and W R − W R contributions. On the other hand, we will not study the possibility of having a NP signal from the W L − W L channel mediated by the RH neutrinos since it would also require some level of fine tuning as demonstrated in Ref. [24] for the type-I seesaw case.
The left and central panel of Fig. 6 show the sensitivity of the next-to-next generation of 0νββ decay experiments (10 −2 eV < |m ββ | < 0.38 eV) to the parameters of the model by including all the relevant contributions and requiring the W L −W R and W R −W R contribution to the 0νββ decay rate (second and third term in Eq. (35)) to be at least 10 times larger than the W L − W L contribution (first term of Eq. (35)) and ξ = M 2 W L /M 2 W R . The allowed region is projected onto the v R -M 1 plane (left panel) and v R -M 2 plane (central panel). We have assumed that U pmns and V are real and fix the light neutrino mass scale to m 1 = 10 −2 eV. As in the previous plots, the bounds on the W R mass [40] and the active-"heavy" mixing have been included [36,52,53]. The X-ray constraints [62], which apply if N 1 is the DM,  [40], the active-"heavy" mixing [53] and the X-ray constraints [62] have been also included.
The PMNS angles and oscillation mass-squared differences have been fixed to the central values given in Ref. [63] and m 1 = 10 −2 eV, while the CP-phases of U pmns and V have been set to zero.
are also shown in Fig. 5: the region between the red dashed lines is ruled out. We have used Eq. (15) in order to be consistent with the light neutrino mass and mixing pattern as mentioned before. However, in order to illustrate better the impact of the mixing θ, we also show in Fig. 6 (right panel) the sensitivity of the next-to-next generation of 0νββ decay experiments to v R and an effective mixing θ ef f defined in the following way. We have  Fig. 6 shows that the v R sensitivity is independent of M 1 is because for a large part of the parameter space the W R − W R channel dominates, where the N 2 and N 3 contribution is important. In summary, a NP 0νββ decay signal can be expected for a quite small light-sterile neutrino mixing (θ ef f 10 −5 ) if one of the RH neutrinos is lighter than 1 MeV while the rest are above the 0νββ decay scale. This is interesting regarding the possibility of accommodating the DM in the left-right symmetric models as we will see in the next section.

IV. COMPLEMENTARY CONSTRAINTS
In this section we will study the impact of the 1-loop corrections on the light neutrino masses and whether the part of the parameter space which can be probed in future 0νββ decay experiments, as described above, is accessible by other experiments.

A. 1-loop corrections
Since, in the scenarios studied here, the light neutrino contribution to the 0νββ decay rate is suppressed with respect to the NP ones, one may expect that this significant NP lepton number violation contribution to the 0νββ decay rate could induce non-negligible 1-loop corrections to the light neutrino masses. Of course, if the 1-loop corrections are larger or similar to the tree-level contribution, they should be included in the analysis, which would modify our previous conclusions. The leading 1-loop correction to the light neutrino masses is given by [64,65]   . This is correct unless some cancellation is at work for the tree-level contribution, which is not the case studied here. Notice that, in this sense, the assumptions made in order to obtain Eq. (37) are quite reasonable.
Therefore, we can conclude that the one-loop corrections to the light neutrino masses are negligible and not relevant in our analysis. The lepton number violation source of the dominant NP contributions studied in the previous sections is the Majorana mass term generated dynamically for the RH neutrinos. Indeed, this source of lepton number violation is related to the light neutrino masses through the seesaw mechanism, and this correlation has been taken into account in the previous analysis. The dominant NP contribution to the 0νββ decay coming from the W R − W R channel (or the W L − W R channel in the mixed scenario) requires the suppression of the standard (and long range) light neutrino one. Since we are not facing the possibility of having any cancellation in the light neutrino contribution, in order to achieve this suppression the Yukawa couplings and v L /v R should be small. The W R − W R contribution can be dominant because the RH mixing is not constrained in contrast with the active-heavy mixing θ, which is necessarily small as the Yukawa couplings. The W L − W R channel can dominate in the mixed scenario (only for large ξ) due to the enhancement coming from the NME and the linear dependence on the active-heavy mixing θ.

B. Other experimental bounds
The charged LFV experiments are also sensitive to the parameters of the model that can be probed in 0νββ decay experiments. Among them, µ → eγ, µ → 3e and µ → e conversion give the stronger bounds. First of all, the small active-heavy mixing required here in order to have a significant NP contribution to 0νββ decay (θ 10 −5 ), renders the type-I seesaw like contribution to this processes completely negligible since the strongest present bound coming from µ → eγ gives the constraint |(θ † θ) eµ | 10 −5 . A complete calculation of the charged LFV branching ratios in the MLRSM can be found in [29]. The most relevant constraint in the context of this work comes from µ → eγ, whose branching ratio, to zeroth order on θ and ξ, is given by: for M ∆ L,R M W R . Applying the present experimental constraint [66] to Eq. (38), the bound on M R reads Saturating the lower bound on M W R , one obtains M R 1 TeV. This is the largest M R that can be probed with 0νββ decay experiments as it can be seen in Fig. 4 (it corresponds 23 to the bottom right corner of the shaded regions in the left panels). This means that the future µ → eγ experiments can be sensitive at least to that corner of the parameter space, which can also give a signal in 0νββ decay experiments. One should, however, keep in mind that the flavour structure of M R plays an important role, being µ → eγ experiments indeed sensitive to (M R M * R ) µe and 0νββ decay mainly to (M R ) ee if the heavy neutrino spectrum is hierarchical. Therefore, a NP signal in 0νββ decay experiments does not necessarily imply also a signal in future µ → eγ experiments. On the other hand, the bound in Eq. (39) has been extracted assuming that M ∆ L,R M W R , but smaller masses of the triplets can clearly enhance the branching ratio [29]. The same applies for the µ → 3e and µ → e conversion case since their branching ratios are inversely proportional to the triple masses. Therefore, we can not extract a bound like Eq. (39) from µ → 3e and µ → e conversion, since for large triplet masses the branching ratios are very suppressed.
So far, in our LFV analysis we have neglected the ξ contribution. If one switches on the left-right mixing ξ, the following constraint from µ → eγ can be extracted [33]: which can be roughly translated into θξ 2 KeV M R , which is basically compatible with most of the parameter space that can give a NP signal in 0νββ decay, since ξ < 10 −3 and θ < 10 −5 .
Only if M R is close to the TeV and ξ saturates the present bound (ξ ≤ M W L /M W R < 10 −3 ), could the mixing ξ have an impact in µ → eγ. Basically, we could probe the same part of the parameter space commented above but for values of ξ close to its present bound, again with the important warning that the flavour structure plays an essential role here.
Finally, in the MLRSM the electric dipole moment (EDM) can be considerably enhanced with respect to the SM result (up to 10 orders of magnitude). This is because the SM contribution to the EDM appears at four loops while the left-right symmetric model can provide a huge enhancement due to the left-right mixing ξ [67]. In fact, the EDM experiments can be sensitive in the future to part of the parameter space studied here [31], mainly through the imaginary part of [(m D ) ee ξ].

V. DARK MATTER
In this section we study the possibility of having a successful DM candidate in the context of the MLRSM when the 0νββ decay rate is dominated by NP contributions. The first question which arises from the results of the previous section is whether N 1 can be DM in the light regime, namely with mass M 1 O(MeV). This reminds us of the Dodelson-Widrow (DW) scenario [68], where a KeV neutrino is produced via neutrino oscillations and can be a viable DM candidate 7 . In the left-right symmetric models, however, a RH KeV neutrino N 1 would be thermally produced via the W R or Z R exchange and decouples from the thermal bath at the freeze-out temperature T f , where g * (T f ) is the number of relativistic degrees of freedom at freeze-out. The rule of thumb to estimate T f is to set the interaction rate equal to the expansion rate of the Universe. Given that we are interested on the region of the parameter space in which the NP dominates the 0νββ decay rate, the W R mass should be in the range M W R ∼ 1-15 TeV (see Figs. 6).
Therefore, N 1 is highly relativistic (M 1 MeV T f ) at freeze-out and the resulting relic density is [70], which, for M 1 ∼ KeV, would be much larger than the observed DM relic density Ω DM = 0.265 [71]. This constraint is much severer than the X-ray constraints shown in Fig. 5 and the Big Bang Nucleosynthesis (BBN) bound, which is indeed still compatible at ∼ 2σ with the existence of one extra relativistic species [72,73]. A possible way out has recently been proposed and studied in detail in Refs. [70,74]. Basically, the idea is to dilute the number density of N 1 by the injection of entropy into the thermal bath after N 1 freezes out. To be more specific, the out-of-equilibrium decay of N 2 and/or N 3 , of mass around GeV, into SM particles can increase the entropy of the Universe, leading to faster Universe expansion and in turn a smaller N 1 density. The set of constraints that should be satisfied if N 1 as 7 A recent study for KeV-neutrino DM on the 0νββ decay in the context of the type-I seesaw can be found in Ref. [69], where the KeV neutrino contribution to the 0νββ rate is subleading due to the X-ray bound.

25
DM was once in thermal equilibrium has been summarized in Ref. [74]. Ref. [70] for the details of the analysis. In any case, as it was already pointed out in Ref. [70], the contribution to the 0νββ decay rate from the W R − W R channel associated with such spectrum can be testable in the future 0νββ decay experiments as we have confirmed in Sec. III C. However, we would like to remark that the W L − W R contribution can also be very relevant in this case, as it was explained in the previous section.
Finally, in the left-right symmetric models, in principle the neutral component of the right-handed triplet ∆ 0 R , which is a singlet under the SM gauge group, could also be a DM candidate. However, it decays at one-loop into two photons via W R exchange [70], i.e., The X-ray constraints on KeV DM resulting from observations on galaxies and clusters of galaxies [75] requires τ ∆ 0 R →γγ = 1/Γ ∆ 0 R →γγ 10 28 sec or Γ ∆ 0 R →γγ 10 −52 GeV. Therefore, a KeV ∆ 0 R would imply a too heavy W R such that the contribution of the NP channels involving W R to the 0νββ decay would be completely negligible. Nevertheless, the X-ray constraints apply only to DM with masses around 1-20 KeV. In principle this leaves another window of ∆ 0 R mass which can be studied. However, other constraints make this possibility quite unfeasible. First, the mass of DM is constrained to be larger than KeV [76] because of the Lyman-α observations. Second, for M ∆ 0 R 20 KeV, τ ∆ 0 R still has to be longer than the age of the Universe, around 10 18 sec, which results again in a very heavy W R that renders any NP contribution to 0νββ decay far beyond the future experimental sensitivity.
In summary, in spite of the existence of various constraints, the left-right symmetric models can accommodate a KeV RH neutrino as a successful DM candidate which can lead to a NP signal in the next-to-next generation of 0νββ decay experiments.

VI. CONCLUSIONS
We have studied the 0νββ decay phenomenology in the MLRSM. In particular, we have analyzed under which conditions a 0νββ decay signal can come mainly from NP contributions associated with this model. Special attention has been paid to the correlation among the different NP contributions and the standard light neutrino one. This correlation emerges from the neutrino mass generation mechanism and should always be considered in the analysis. The scenario in which an accidental cancellation in the W L − W L contribution takes place has not been explored. The role of the W L − W R mixing ξ and the possibility of having a NP dominated 0νββ decay signal compatible with DM is also investigated.
We have distinguished three different regions of the parameter space based on the mass of the RH neutrinos : (i) all the masses heavier than GeV, denoted by heavy regime; (ii) masses lighter than MeV, dubbed light regime; (iii) the lightest mass below the MeV and the rest above GeV, called mixed scenario. Notice that (i) has been extensively studied in the literature, but (ii) and (iii) have not been analyzed before in detail for the left-right symmetric models (at least the 0νββ decay phenomenology).
In the heavy region, we have found that the dominant NP contribution emerges mainly from the W R −W R channel mediated by the heavy neutrinos. To be more precise, it has been shown that this dominant NP contribution could be measured in the next-to-next generation of 0νββ decay experiments for the window M W R ∼ 1 − 15 TeV (M W R ∼ 1 − 20 TeV), if the active-heavy mixing is smaller than ∼ 10 −5 and the right-handed triplet "Yukawa" coupling satisfies 3 · 10 −6 (Y ∆ ) ee 3 · 10 −2 (3 · 10 −6 (Y ∆ ) ee 8 · 10 −2 ), which corresponds to a range of heavy neutrino masses from GeV to TeV for ξ = 0 (ξ = M 2 W L /M 2 W R ). We have also shown that neglecting the present correlation between the various contributions, and that with the light neutrinos in particular, can lead to incorrect results. For instance, it is found that the region of the parameter space which can be experimentally probed when only the W L − W R contribution is included in the analysis shrinks considerably if all the contributions are included at once and their correlations are not neglected.
The results for the light region turn out to be similar to those of the heavy region. We have found that a future 0νββ decay NP signal can come only from the W R − W R channel since the W L − W R contribution cancels out in this regime. In particular, we have shown in which part of the parameter space this is possible and we found a similar but weaker sensitivity Indeed, a signal from the W R − W R channel can be expected in a larger region, even for smaller values of θ since its contribution is independent of the active-heavy neutrino mixing.
In this case we have focused our study on the limit v L /v R → 0, but we have checked that if v L is switched on in the analysis the future sensitivity to v R is much better. However, this can only take place when there is a cancellation between the type-I and type-II seesaw terms in the light neutrino contribution.
In order to study the impact of the W L − W R mixing ξ, we have analyzed the following two extreme limits: ξ = M 2 W L /M 2 W R (maximal) and ξ = 0 (negligible). We have shown that the inclusion of the W L − W R mixing ξ can have some impact on the results but it is not very significant, with the possible exception of the mixed scenario where a large mixing is required in order to have a relevant role of the W L − W R channel. In general, due to the enhancement on the W R − W R contribution for maximal mixing, the sensitivity to M W R is about a factor 1.5 larger for ξ = M 2 W L /M 2 W R than for ξ = 0 in all the regions under study. We have also analyzed the role of the complementary bounds coming from charged LFV processes and the the impact of the 1-loop corrections to the light neutrino masses in the context of this work. It turns out that the light neutrino masses are stable under 1-loop corrections since they might be important only if a cancellation takes places in the light neutrino masses, but not in the case studied here where there is a general suppression of the light masses with small v L /v R and the Yukawa couplings. Future charged LFV experiments might allow us to probe part of the parameter space that can be responsible for a NP signal in 0νββ decay experiments, but only a small region in the heavy regime around M R ∼ 1 TeV (the bottom right corner of the shaded regions in the left panels of Fig. 4). In fact, a more complete study, beyond the scope of this work, including the effect of triplet masses close to their lower bounds, which can enhance the branching ratios, would be required in order to clarify the issue. A large left-right mixing ξ can also be probed in future EDM experiments as it was shown in Ref. [31].
Finally, the following DM-related question has also been addressed. Can a NP dominated 0νββ decay signal be compatible with a successful DM candidate in the left-right symmetric models? We conclude that, regardless of the various strong constraints, it is still possible for the scenario proposed in Ref. [74], where a KeV RH neutrino can be the DM if the scale of the other heavy neutrinos is around 1 − 10 GeV and M W R 10-15 TeV. We have shown that the 0νββ decay signal can be induced by the RH neutrinos through the W L − W R and W R − W R channel. Additionally, Ref. [70] opens a window of M W R ∼ 5 GeV within the horizon of LHC after the QCD phase transition is carefully included.