{'fraction_non_alphanumeric': 0.037233865878367155, 'fraction_numerical': 0.011791277775932507, 'mean_word_length': 4.296042216358839, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 7, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 2, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Baryogenesis from L-violating Higgs-doublet decay in the density-matrix formalism

We compute in the density-matrix formalism the baryon asymmetry generated by the decay of the Higgs doublet into a right-handed (RH) neutrino and a Standard Model lepton. The emphasis is put on the baryon asymmetry produced by the total lepton-number violating decay. From the derivation of the corresponding evolution equations, and from their integration, we find that this contribution is fully relevant for large parts of the parameter space. This confirms the results found recently in the CP-violating decay formalism with thermal corrections and shows in particular that the lepton-number violating processes are important not only for high-scale leptogenesis but also when the RH-neutrino masses are in the GeV range. For large values of the Yukawa couplings, we also find that the strong washout is generically much milder for this total lepton-number violating part than for the usual RH-neutrino oscillation flavour part.


I. INTRODUCTION
In the framework of the type-I seesaw model of neutrino masses with right-handed (RH) neutrinos below the electroweak scale, there exists a well-known mechanism to account for the baryon asymmetry of the Universe, through oscillations of right-handed neutrinos [1][2][3][4][5][6][7][8][9][10][11][12][13]. This Akhmedov-Rubakov-Smironv (ARS) scenario is based on the generation of particle-antiparticle asymmetries for the various lepton flavours. These asymmetries cancel each other in the total Standard-Model (SM) lepton-number asymmetry but, thanks to washout effects, which do not affect the different flavours in the same way, a net lepton asymmetry remains. In this framework, since the relevant processes do not involve a RH neutrino Majorana mass insertion, a lepton number can be assigned to the two helicities of the RH neutrinos. Thus, in the ARS scenario, the total lepton number L, i.e. the sum of the SM and the RH-neutrino ones, is conserved, but not both components separately, due to flavour effects. The SM lepton-number asymmetry component which is produced in this way before the sphalerons decouple is reprocessed in part into a baryon asymmetry, unlike the other component. The evolution of the lepton asymmetries as a function of the temperature of the thermal bath can be calculated in the densitymatrix formalism, which properly takes into account the coherences between various RH neutrinos and their associated oscillations.
In the different CP-violating decay formalism usually used for leptogenesis, it has been shown recently [14] that the total lepton-number violating decay of the Higgs doublet into a RH neutrino and a SM lepton, i.e. the decays which do involve a Majorana mass insertion, could also account for the baryon asymmetry. This is possible thanks to thermal effects which induce a non-zero absorptive part for the self-energy of the RH neutrino in * thambye@ulb.ac.be † daniele.teresi@ulb.ac.be the final state of this decay, see Fig. 1 (as this absorptive part vanishes at zero temperature). In Ref. [14] this thermal cut of the self-energy has been computed in the Kadanoff-Baym formalism. The production takes advantage of the fact that, for RH neutrino masses below the sphaleron cut, the Yukawa interactions do not thermalize the RH neutrinos so easily as for higher masses. This results in a large departure from equilibrium for the RHneutrino number densities in the final state, boosting the asymmetry production. This total (SM + RH neutrino) lepton number violating scenario is in many respects different from the ARS total lepton-number conserving scenario. One difference is that it gives a non-vanishing asymmetry already for one lepton flavour. Another is that, as a result of the fact that it involves a Majorana mass insertion, it gives rates for the relevant processes that are proportional to m 2 N /T 2 , relatively to the total lepton-number conserving ARS piece. As a result, the asymmetry is typically produced at lower temperatures than in the ARS case, basically not long before the sphaleron decoupling. Also, since this mechanism does not require an asymmetric washout for the different flavours, its contribution in the weak-washout regime is proportional to 4 powers of the Yukawa couplings, rather than to 6 for the ARS one.
The purpose of this paper is twofold. The first goal is to determine how this L-violating contribution can develop itself in the framework of the density-matrix formalism, rather than in the usual leptogenesis CPviolating decay formalism considered in [14], and to com- pare the results obtained in both formalisms in the single lepton-flavour case, where the ARS mechanism is not operative. We will find below a good agreement between the results of both approaches showing that, also in the density-matrix formalism, the L-violating contribution is fully relevant. Moreover, as will be discussed in detail below, for large Yukawa couplings, the density-matrix formalism gives even a larger asymmetry, because the actual washout for the total lepton number is actually much smaller than the total one, considered in [14]. The second purpose is to compare, within the same formalism, which of the L-violating (LV) and L-conserving (LC) contributions dominates the production of the baryon asymmetry as a function of the parameters. In fact, evolution equations including the L-violating contribution have been derived long ago in the density-matrix formalism [3,7], and also very recently [15,16], but this was considered as numerically irrelevant for the generation of the baryon asymmetry. This stems from the fact that L-violating processes have been thought to lead to a negligible contribution due to the m 2 N /T 2 suppression factor discussed above (and thus approximated as zero in the solution of the evolution equations for baryogenesis [7]) or, very recently, because the corresponding obtained rates were considered as numerically "not significant" [15]. Numerically we will find below that the L-violating contribution is clearly significant, as it dominates over the flavour violating piece in large portions of the parameter space, which we determine.
We will do all that considering only the contribution of the Higgs doublet decay/inverse decay process. Other processes, such as the top-quark and gauge scattering [15,17,18] and the infrared-enhanced gauge corrections to Higgs decay [19], are also expected to be fully relevant. But they are expected to lead to an asymmetry in a similar way to the Higgs-doublet decay because they involve as a subprocess the same Higgs-doublet transition into a RH neutrino and a SM lepton. Thus we anticipate that the relative importance of the two contributions in the various regions of parameter space is well captured by the analysis presented in this paper at the level of this decay basic building block. A fully quantitative calculation of the asymmetry produced, including these processes, is left for a further work. All the numerical results below will be obtained considering 2 RH neutrinos, which simplifies the discussion.

II. DENSITY-MATRIX EQUATIONS INCLUDING LEPTON-NUMBER VIOLATION
We start by deriving in this Section the density-matrix equations which allow to determine the evolution of the number density of the RH neutrinos N α (which have + helicity), of the conjugated statesN α (which have -helicity) and of the lepton asymmetries for the various lepton flavours. The Lagrangian that we consider is the one of the usual type-I seesaw neutrino mass model, with noth-ing else, where H ≡ iσ 2 H with H = (H + , H 0 ) T and L l = (ν l L , l − L ) T . In the same way as in Ref. [14], we will make the approximation that the 4 scalar states in the doublet remain on the same footing at all relevant temperatures, meaning in particular that they have same masses at all relevant temperatures.
The matrix number densities for the two helicities of the RH neutrinos are and analogously for n L l ≡ n L ll ,n L l ≡n L ll (having neglected flavour coherences between the SM leptons) and n H ≈ n H ≈ n H eq (having neglected the asymmetry stored in all the SM species but the lepton doublets).
We start from the Markovian master equation [20] (see [21] for a detailed derivation): where H N 0 is the free Hamiltonian for the RH neutrinos and where the dispersive terms, giving rise to thermal corrections to the masses, are not written explicitly. Proceeding by a set of standard manipulations, we obtain where E N ≡ E N (k). We have left the t-dependence implicit and the production and absorption rates Γ ≶ are given by where iS ≶ and i∆ ≷ denote, respectively, the SM-lepton and Higgs-doublet absolutely-ordered Wightman propagators, with the Higgs-doublet one approximated to be in thermal equilibrium. Their explicit form is well-known and can be found, for instance, in [22]. Alternatively, the transport equation (5) can be obtained directly from the field-theoretic Kadanoff-Baym formalism [12,23,24], once the so-called quasi-particle approximation is performed, and generalized to other processes by taking into account the appropriate self-energies. By assuming kinetic equilibrium for all relevant species, we find the rate equation for the (momentumintegrated) number-density matrix n N αβ : where the thermally-averaged energy matrix is given by The equilibration reaction densities are where the phase-space integral is As for the washout reaction densities, they are γ LC W C,l,αβ ≡ dΠ PS n N eq (k)n L eq (p) γ LV W C,l,αβ ≡ dΠ PS n N eq (k)n L eq (p) The "Q" label refers to terms which are present only due to the quantum character of the particle statistics, unlike the "C" labelled terms which are also present for "classical" Maxwell-Boltzmann statistics. The equation forn N can be obtained easily from (5) and (6) by replacing the helicity-plus spinors with the minus ones in the traces, and by making use of the identities leading to The equation for the asymmetry δn L l can be obtained by noticing that the combinations l δn L l + α (n N αα − n N αα ) and l δn L l − α (n N αα −n N αα ) are not changed, respectively, by LC and LV processes. Therefore where the l index is not summed. The reaction densities γ LC l and γ LV l are defined as in (10) and (11), taking in them only the term with l indexed Yukawa couplings.
The general structure of the evolution equations (8), (19) and (20) agrees with the ones presented for the LC part in [7,12,13,15] and for the LV part in [7,15], within the various approximations performed. More precisely, for the LC part these equations agree with [13], and also with [7,12], up to the fact that we do not make the extra approximation to take all washout terms as given by a same term whose form is similar to the source term, and up to the fact that [7] involves an extra 1−2n N eq (k) factor in the equilibration rates of Eqs. (10) and (11) above. For the LV part we also agree with [7] up to the same 1 − 2n N eq (k) factor difference. Finally with respect to the very recent Ref. [15], the approximations made for both the LC and LV parts are not exactly the same, resulting in different equations, but the effect of these different approximations is expected to be small.

III. REACTION DENSITIES
Next in this Section we compute the "γ" reaction densities entering in the density-matrix evolution equations above. In the ultra-relativistic regime, i.e. for m N |k|, E N (k) |k| + m 2 N /(2|k|) and from (9) we have where in the last equality we have dropped the flavour independent first term because it drops in the commutators in the evolution equations. The thermal averaged T /|k| is equal to 0.46 for the Fermi-Dirac statistics.
The spinor traces appearing in the lepton-number conserving and violating reaction densities are found to be where θ is the angle between k and p. One recognizes that for negative helicity, i.e. for LV transitions, a m 2 N factor appears, as it should be, since LV requires a Majorana mass insertion. For the washout reaction densities, denoting by a tilde their value once the Yukawa couplings have been factorized out, e.g. γ LC l,αβ ≡ γ LC h * lα h lβ , we find where (neglecting m N compared to the other masses) and we have included g iso = 2 to account for the two isospins of the processes (e.g. H 0 → Nν and H + → N e + ). For the thermal masses of H and L we have used the high-temperature asymptotic masses (i.e. the ones in the large 3-momentum limit)

IV. ANALYTIC SOLUTIONS IN THE WEAK-WASHOUT AND MANY-OSCILLATIONS REGIME
Before presenting the numerical results obtained by integrating the density-matrix evolution equations, it is instructive to compare the LC and LV analytic solutions which can be obtained in the weak-washout regime for the case of many oscillations taking place, i.e. when oscillations start long before sphaleron decoupling. The derivation of these analytic solutions is given in Appendix A. They are determined at first non-vanishing order in powers of the Yukawa couplings. We give them for the lepton asymmetry yield for each flavour l, defined as the lepton asymmetry obtained for this flavour over the entropy density, i.e. Y l ≡ n L l /s. Notice that in the weak-washout regime considered in this section the total asymmetry is simply given by the sum of the ARS and the LV H-decay one. For the LC contribution, summing over the lepton flavour, Y ≡ l Y l , at the sphaleron decoupling temperature, T = T s 131.7 GeV where the baryon asymmetry freezes, we get the following total lepton asymmetry, which involves 6 powers of the Yukawa couplings, and M 0 M P L /(1.66 g 1/2 * ) 7 × 10 17 GeV. This result involves three rates whose value are given in the previous Section. For the LV contribution, instead, one gets an asymmetry which involves 4 powers of the Yukawa couplings and 2 rates with and δ LV ≡ l δ LV l . Comparing Eqs. (35) and (33) one observes a series of differences. First of all, both contributions involve different combinations of the Yukawa coupling, δ LV versus l δ LC l (hh † ) ll for the LV and LC contributions respectively. As already mentioned in the Introduction, the LC combination vanishes in the case of a single flavour, unlike the LV one. Generically, the combination of Yukawa couplings in the LV case is much larger because it involves 2 Yukawa couplings less. Moreover, unless one chooses precisely specific phases in the Yukawa entries, the ratio of l δ LC l (hh † ) ll to δ LV turns out to be smaller than the naive estimate m ν m N /v 2 by one order of magnitude or more. Next, these 2 contributions display different powers of the Majorana mass scale, of the mass splitting and of the effective Planck scale, M 0 m 2 N /∆m 2 N versus M 7/3 0 /(∆m 2 N ) 2/3 . As discussed above, the m 2 N factor results from the Majorana-mass insertion necessary for the LV contribution. The different powers of the mass splitting result from the different T -dependence of the LC and LV rates, and in particular from the relative size of the commutator term and the "γ" rates as functions of the temperature. As for the different powers in the effective Planck mass M 0 , they result from the integration over t of the different temperature dependence of the contributions. Such a scale ratio is clearly in favour of the LC contribution, especially for small m N and large mass splittings. For instance, with m N = 10 GeV and ∆m N /m N = 10 −10 or with m N = 1 GeV and ∆m N /m N = 10 −5 one gets a relative factor equal to 5.9 × 10 −20 and 5.9 × 10 −23 respectively. Finally, beside a different numerical factor in favour of the LC contribution, -18.5 versus 7.9, the rate factors are very different, α LC α LV versus (α LC ) 2 α LC W . As said above, the LV piece involves 2 rates whereas the LC piece involves 3 of them, as a result of the fact that it requires the asymmetric washout to play a role. Hence, the rate factor is much larger for the LV contribution. Numerically, from Eqs. To compare the results obtained above for the LV contribution in the density-matrix formalism with the ones obtained in the CP-violating decay formalism [14], it is convenient to consider only one lepton flavour and 2 RH neutrinos. Considering only one flavour makes sense because in this case the LC ARS contribution identically vanishes and does not interplay with the LV contribution in a non-linear way (as it does, instead, in the strong washout regime with several flavours).
With one lepton flavour and 2 RH neutrinos there is only a very limited number of parameters, the RH neutrino scale m N , the mass splitting between the two RH neutrinos, ∆m N , the Yukawa coupling of N 1 which can be taken to be real, h 1 , and the Yukawa coupling of N 2 which then has a phase, h 2 = |h 2 |e iθ . The last 3 parameters can be traded for the effective massm ≡ v 2 h 2 1 /m N , which parametrizes the speed of the decay of N 1 , for the ratio of decay widths Γ 22 /Γ 11 and for the CP-violating parameter Starting from a situation where there is no RH neutrinos at T >>> T s , in Figs. 2 and 3 we plot the results for the baryon asymmetry obtained by integrating numerically Eqs. (8), (19) and (20). This is shown as a function of m N andm for the same value of the other parameters as considered in [14]. Also shown on these figures are the results obtained from the CP-violating decay formalism in [14]. These figures nicely show a same qualitative pattern. Quantitatively, for not too small m N , ∆m N /m N and not too large Yukawa couplings, the asymmetries obtained are in general of the same order of magnitude. For small m N and ∆m N /m N the difference is due to the simplified treatment of RH-neutrino thermal masses in the density-matrix equations above, which would necessitate a T -dependent diagonalization of the N mass matrix, not taken into account here for simplicity. For largem the fact that the produced asymmetry obtained can be much larger in the density matrix formalism is not due to a real discrepancy between both methods. This difference stems from the fact that in [14] the washout suppression has been overestimated. As showed above, in the densitymatrix formalism one gets 2 different types of washout. First there is the usual ARS washout term γ LC W ∝ T 4 , Eq. (26), which violates SM lepton number but does not violate total lepton number (as it does nothing more than transferring the asymmetry between the SM-lepton and RH-neutrino sectors). This piece washes out the ARS LC asymmetry produced, but does not wash out the total lepton-number asymmetry produced. Second, there is the γ LV W ∝ m 2 N T 2 washout part, Eq. (27), which violates total lepton number but is not as effective, because suppressed by a factor of m 2 N /T 2 with respect to the other washout. In the CP-violating decay formalism of [14], the fact that a lepton-number asymmetry can be stored in the RH-neutrino sector is not taken into account (but could, in principle, in a more involved decay formalism). Thus, the part of the washout ∝ T 4 incorrectly washes out the final lepton lepton asymmetry too, resulting in an overestimation of the washout suppression. In the density-matrix formalism this is instead correctly taken into account.

B. Results with 3 lepton flavours: comparison of the LC results with the total LC + LV results
We now consider the realistic situation with 3 SM lepton flavours, with seesaw parameters fulfilling all neutrino-mass constraints (with, here too, 2 RH neutrinos for simplicity). This allows us to identify the regions of parameter space where the LV Higgs-decay contribution is important as compared to the LC ARS one.
We fit the light-neutrino oscillation data by using the so-called Casas-Ibarra parametrization [25] where v = 246 GeV, with m diag ν and m diag N denoting, in the mass eigenbasis, the light and heavy neutrino mass matrices, respectively. For the PMNS matrix we use the standard parametrization of [26], fixing all parameters to their best-fit values [26] and δ = −π/2. The Majorana phase α 1 plays no role in what follows, unlike the α 2 phase. Thus we will present our results for different values of α 2 . For definiteness we choose a normal hierarchy of masses. In the 2 RH-neutrino case considered here, the orthogonal matrix R can be expressed in terms of a single complex angle ω = w + iγ. We fix w to the value which maximizes δ LV , w = π/4. For w fixed in this way the Casas-Ibarra parameter γ has a one-to-one relation with the effective neutrino mass m ≡ v 2 (h † h) 11 /m N .
Starting here too from no RH neutrinos at T >>> T s , we have integrated numerically the evolution equations. The final n B /n γ baryon asymmetry obtained is plotted as a function of m N and γ in the left panels of Fig. 4-9,      for various values of ∆m N /m N and α 2 1 . These plots also show the relevant existing constraints and future prospects of the proposed SHiP [27] and FCC-ee [28] experiments. We have also calculated the ARS leptonnumber conserving contribution, obtained by setting to zero the LV rates in the evolution equations, as done previously in the existing literature. The results for the ratio of the full (LV Higgs decay + LC ARS) result to the ARS one are shown in the right panels of Fig. 4-9. The regions left in very light grey on these plots are the ones where this ratio lies between 0.8 and 1.2, i.e. the regions where the incorporation of the LV part does not change sizeably the final result.
Let us start with a mass splitting ∆m N /m N = 10 −10 . Fig. 4-7 show the results for this mass splitting and α 2 = 0, π/2, π, 3π/2. Although the precise quantitative values (and sign) of the ratio between the full and ARS results depend strongly on the phase α 2 , the overall picture is quite stable. For small γ, i.e. for values of the Yukawa couplings of order of the ones given by the basic seesaw scaling, the two contributions have similar size, with the LV part dominating in a large region. This behaviour is already captured by the analytic solutions given in Sec. IV, and has been argued already in [14].
In addition, for large values of γ, the LV Higgsdecay contribution becomes generically dominant, even by many orders of magnitude 2 with respect to the ARS one. These large values of γ (or equivalently m) can be obtained generically if lepton number is a symmetry which is approximately conserved. Such a symmetry gives at the same time a quasi-degenerate spectrum for 2 RH neutrinos, as necessary here. This large-γ region covers a large fraction of the region of parameters accessible at SHiP and FCC-ee, i.e. where these experiments could directly produce the RH neutrino(s). The origin of this LV dominance is clear: for large γ we are in the strong washout regime. In this regime, as already explained at the end of Section V A, the γ LC W washout part will effectively wash out the ARS LC part, whereas the γ LV W washout part, which is suppressed by an extra m 2 N /T 2 factor (see (26) and (27)), will wash out less the LV part (and only at later times due to this factor), resulting in a dominant LV part.
In Fig. 8 and 9 we plot the results for ∆m N /m N equal to 10 −8 and 10 −6 , respectively. As already emphasized in Sec. IV, for the weak-washout regime m ≈ m ν , the larger is ∆m N /m N , the less is important the LV contribution. However, the LV effects are still important in the strong-washout region and still dominate by many orders of magnitude. Nevertheless, in this case the parameter space leading to successful leptogenesis shrinks, and the region of successful leptogenesis where the LV contribution dominates becomes smaller, even if still significant. For ∆m N /m N = 10 −6 , from the right panel of Fig. 9 we see that the LV contribution is still largely dominant for m N ≈ 3 GeV, γ ≈ 4 − 5. For larger m N this dominance is expected to increase further but in this case the numerical solution of the density-matrix equations becomes computationally challenging and a more sophisticated numerical code would be needed (especially for even larger values of ∆m N /m N ).
We conclude this section by discussing the dependence on the initial "reheating" temperature T in , where one has assumed that there are no RH neutrinos to start with. As already argued in [14], the LV Higgs-decay mechanism is a low-scale one, in the sense that the asymmetry is generated mainly at temperatures close to the sphaleron decoupling one. Thus, it does not rely on UV physics well above the electroweak scale. The LC flavoured asymmetries in the ARS scenario, instead, are mainly produced at T ≈ T osc = (M 0 ∆m 2 N ) 1/3 (at least in the weakwashout regime), which easily lies well above T s . We illustrate this in Fig. 10, where we plot the LC and LV contributions as a function of the initial temperature T in in the weak-washout regime for ∆m N /m N = 10 −7 . The former is obtained, again, by setting to zero the LV rates, whereas the later is obtained as the total one minus the LC one, since in the weak-washout regime the two contributions are additive. For the parameters chosen in Fig. 10, even though the LV contribution is subdominant when T in > T osc , it becomes dominant, even by 2 orders of magnitude, if the initial temperature of the evolution is smaller.

VI. CONCLUSIONS
In summary, we have computed, in the density-matrix formalism, the contribution to the baryon asymmetry of the total lepton-number violating Higgs-doublet decay and inverse decay. This contribution has been obtained by deriving the corresponding density-matrix evolution equations, and by integrating them. The results obtained confirm the ones derived in the CP-violating decay formalism [14]: this contribution is fully relevant for large parts of the parameter space. This is the case even if the L-violating processes are suppressed by a factor of m 2 N /T 2 with respect to the L-conserving ones, by virtue of a different parametric dependence, in particular in the Yukawa couplings and rates.
In the one lepton-flavour case the results for the LV contribution of both formalisms can be easily compared and show a nice qualitative agreement, or even quantitative agreement for not too large values of the Yukawa couplings. For larger Yukawa couplings, leading to strong washout, the density-matrix formalism gives a larger asymmetry because the washout of the LV contribution turns out to be much less efficient than what we could have expected at first sight.
In the realistic 3 lepton-flavour case, and for relatively small RH-neutrino mass splitting, the LV contribution dominates over the ARS LC contribution for Yukawa couplings of the order of the ones given by the usual neutrinomass seesaw formula without cancellations, and also for Yukawa couplings much larger than these. This latter dominance is, here too, due to the fact that the washout of the LV contribution is not easily efficient, and is actually much less efficient than the washout of the ARS LC contribution. For larger mass splitting this dominance of the LV contribution remains generic for large Yukawa couplings. Also, as a result of the fact that the LV part involves the mass insertion m 2 N factor, the larger m N the more the LV part dominates over the ARS LC part. Generic dominance of the LV contribution is also obtained for low reheating temperature.
All the above takes into account exclusively the baryon asymmetry produced by the Higgs-doublet decay and inverse decay. It is known that other processes can be fully relevant, in particular the top-quark and gauge scattering and the infrared-enhanced gauge corrections to Higgs decay [15,17,18]. These processes are expected to generate a baryon asymmetry in the same way because they involve as a subprocess the same Higgs-doublet transition into a RH neutrino and a SM lepton. One could therefore anticipate that, including them, the rates will be enhanced by a factor of 3 to 6 in the same way as for the ARS LC contribution [17], thus enhancing the asymmetry produced too. A fully quantitative calculation of the baryon asymmetry produced would nevertheless re-quire a proper inclusion of all these relevant LC and LV processes.
number conserving CP violating Yukawa coupling combination 22 (A4) In the regime ∆m 2 N µ 2 osc , Imf (z) performs many oscillations (of increasing frequency) before the sphalerons decouple, and its integral in (A2) saturates at early times to the value (A7) In order to obtain a non-vanishing total asymmetry one has to perform one iteration more in the solution of the Boltzmann equations, in particular taking into account at leading order the effect of a flavour-asymmetric washout. By exploiting the fact that, neglecting lepton-number violating processes, Y LC ≡ l Y l = tr(ρ −ρ), one finds at O(h 6 ): which gives approximately Eq. (33).

L-violating Higgs-decay contribution
In the same regime as before we can now calculate analytically the contribution of the lepton-number violating processes. Proceeding iteratively, there are two contributions: the first one comes from γ LV in the equations for ρ andρ and γ LC in the equation for Y l . We find this contribution to be and where δ LV l is the total lepton-number violating and CP-violating combination of the Yukawa couplings 22 . (A11) Notice that we have already summed over SM-lepton flavours, since δ LV ≡ l δ LV l = 0, i.e. the LV Higgsdecay mechanism gives rise to an unflavoured asymmetry already at O(h 4 ), differently from the ARS one. In the regime ∆m 2 N µ 2 osc the integral of Im g(z ) approaches the behaviour up to a residual oscillating term (which we neglect) whose frequency increases and amplitude decreases with time. Thus, The second contribution comes from the iteration with γ LC in the equations for ρ andρ and γ LV in the equation for µ. This is given by The integral is again found to approach the behaviour up to an oscillating term which with time increases its frequency but does not change its amplitude (equal to (µ 2 osc /∆m 2 N ) 4/3 ). Neglecting this term we find that this second contribution is the same as the first one (A13). The sum of both terms leads to Eq. (35).