Low-scale Leptogenesis with Minimal Lepton Flavour Violation

We analyse the feasibility of low-scale leptogenesis where the inverse seesaw (ISS) and linear seesaw (LSS) terms are not simultaneously present. In order to generate the necessary mass splittings, we adopt a Minimal Lepton Flavour Violation (MLFV) hypothesis where a sterile neutrino mass degeneracy is broken by flavour effects. We find that resonant leptogenesis is feasible in both scenarios. However, because of a flavour alignment issue, MLFV-ISS leptogenesis succeeds only with a highly tuned choice of Majorana masses. For MLFV-LSS, on the other hand, a large portion of parameter space is able to generate sufficient asymmetry. In both scenarios we find that the lightest neutrino mass must be of order $10^{-2}\text{ eV}$ or below for successful leptogenesis. We briefly explore implications for low-energy flavour violation experiments, in particular $\mu \rightarrow e\,\gamma$. We find that the future MEG-II experiment, while sensitive to MLFV in our setup, will not be sensitive to the specific regions required for resonant leptogenesis.

The type-I seesaw model [1][2][3][4][5] is the simplest known extension of the Standard Model (SM) that simultaneously addresses the origin of neutrino mass and the generation of the cosmological matter-antimatter asymmetry, the latter through thermal leptogenesis [6] (for reviews see [7][8][9]). One extra heavy Majorana sterile fermion is added per generation, where their masses are required to be higher than the critical temperature of the electroweak phase transition.
This scenario has been extensively studied and it has long been known that in the standard hierarchical 1 type-I scenario of thermal leptogenesis, generation of the necessary asymmetry places a lower-bound [10,11] on the heavy sterile neutrinos (SNs) of order 10 9 GeV.
This prevents low-scale realisations of this model in its simplest form. However, an upperbound [12,13] can be placed on the SN neutrino mass of order 10 7 GeV by requiring that radiative corrections to the Higgs mass parameter µ 2 driven by the high seesaw scale does not exceed 1 TeV 2 . A clear tension exists as there is no overlap between these two bounds.
This can be resolved in a number of ways; by invoking supersymmetry to stabilise the Higgs mass, permitting a departure from strictly hierarchical heavy sterile neutrino masses [14,15], adding a second Higgs doublet [16,17] or adopting a modified lepton sector which allows for low-scale leptogenesis 2 from much lighter SN masses. Note that the first solution, while preventing large corrections, may generate other problems, e.g. cosmological gravitino overproduction [18][19][20].
Lowering the thermal leptogenesis scale is an attractive and obvious way to resolve the above tension, and forms the basis of the present analysis. The scale of leptogenesis can be significantly lowered if a quasi-degenerate (meaning not exactly degenerate) spectrum of masses for the SNs is assumed [21,22]. In the minimal type-I scenario this has the added consequence of suppressing the Yukawa couplings and therefore reducing discovery signals such as charged lepton flavour violation (cLFV). Additionally, a quasi-degeneracy should be motivated by some symmetry if such models are to be taken seriously. It is natural to expect that the origin of such a theory would tie in with the flavour problem of the SM.
A possible method of lowering the scale of thermal leptogenesis with potential cLFV signals is by introducing two SN states (with opposite lepton number assignments) per light neutrino and simultaneously promoting lepton number to being a "good symmetry" of the theory. For certain regimes of the couplings a double suppression of the active neutrino masses occurs, allowing significantly larger Yukawa couplings for SNs at TeV-scale masses compared to the minimal scenario. Additionally, due to the weakly broken lepton number symmetry, the heavy SNs form pseudo-Dirac states with mass splittings proportional to the small lepton number breaking terms. For specific choices of parameters, two popular limiting cases constitute the "inverse seesaw (ISS)" [23][24][25][26][27][28][29][30] and the "linear seesaw (LSS)" [31][32][33][34] models. The two cases can be linked by a rotation [35] only in the case where no additional 1 Hierarchical is typically taken to mean 200m N1 < ∼ m N2 < ∼ m N3 . 2 We define low-scale leptogenesis to occur at O(10 − 100 TeV) or below. symmetry precludes such a rotation (e.g. left-right or Pati-Salam symmetry [4,[36][37][38][39]). We will therefore be treating the two as independent scenarios in what follows.
In the specific case of the LSS above the electroweak phase transition, the heavy SNs are degenerate in mass and therefore self-energy diagrams vanish, leaving only the highlysuppressed vertex contributions to produce the CP asymmetry [40]. In the case of the ISS, a mass splitting does exist allowing for a resonant enhancement in the self-energy contribution for a specific choice of the mass splitting. However, decreasing the mass splitting also increases the efficiency of washout making asymmetry generation difficult in its simplest realisations [41][42][43][44][45] for ISS at the TeV scale.
In addition to the intra-family SN degeneracies, there is also the possibility of a quasidegeneracy between SNs from different families, a flavour degeneracy in other words. To that end we explore a scenario in which an SU (3) 2 × SO(3) 2 flavour symmetry, to be defined below, is utilised in the lepton sector. Due to the SO(3) symmetry, mass terms for the heavy SNs are proportional to the identity matrix and lead to a flavour degeneracy amongst the SN states in both the ISS and LSS scenarios. However, radiative effects induce SO(3)breaking terms which in turn break these degeneracies, enabling a resonant enhancement in the asymmetry generation for both the LSS and ISS scenarios.
We work in the framework of an extended version of the Minimal Lepton Flavour Violation (MLFV) [46,47] scenario, which is itself an extension of the well-established idea of Minimal Flavour Violation (MFV) within the quark sector [48][49][50]. The MFV scheme is an ansatz designed to address the very stringent bounds placed on new physics from various rare hadronic decays and neutral meson oscillations. It is motivated by the idea that taking the Yukawa sector couplings to zero recovers five separate 3 flavour space rotations which leave the Lagrangian invariant. Therefore an MFV or MLFV Ansatz is based on the core principle that these are "good" flavour symmetries which are broken solely by the SM Yukawa fields.
Under this ansatz, Yukawa couplings are said to arise from dynamical fields typically known as 'spurions' which transform under the flavour symmetries such that the Lagrangian is flavour invariant. These spurion fields are taken to have non-zero vacuum expectation values (VEVs) which spontaneously break the flavour symmetries and thereby induce the masses and flavour mixing we observe. As a result, in its most minimal realisation, any flavour changing process can be predicted from the well-measured Cabibbo-Kobayashi-Maskawa 3 In the absence of neutrino mass. If neutrino mass is included more rotations are present.
(CKM) mixing matrix and the fermion masses.
Due to the multiple distinct ways in which neutrino masses can be incorporated into the SM there is no unique way of defining an MLFV Ansatz, in contrast with MFV in the quark sector. The potential connection between an MLFV Ansatz and resonant leptogenesis 4 has been identified previously [52][53][54] in the type-I seesaw context. Since then it has been discovered that, due to developments within flavoured leptogenesis, there is a significant reduction in the allowed parameter space 5 [55][56][57] specific to MLFV leptogenesis. Additionally, in order to produce large leptonic flavour violation whilst suppressing neutrino masses in such models, a separation between the scale of lepton-number violation (LNV) and lepton-flavour violation (LFV) should exist which is not present in the type-I seesaw alone [53,58]. Neither of these issues are present for MLFV in the ISS and LSS, allowing the possibility of MLFVinduced resonant leptogenesis in these scenarios. We assume an ansatz such that the heavy sterile Majorana neutrinos introduced have exactly degenerate masses in the Lagrangian, motivated by the flavour symmetry of the theory. Radiative effects will break this down to a quasi-degeneracy between the SNs from which a resonant enhancement in the asymmetry generated per decay of SN will occur. For some choices of parameters which we will identify, this enhancement is able generate the necessary baryon asymmetry observed in the universe.
We will explore the viability of leptogenesis both from PMNS phases (particularly the Dirac phase) alone [43,[59][60][61] as well as from CPV parameters related to the heavy SNs. This paper is organised as follows. Section II reviews the basic idea behind the MFV and MLFV Ansätze including justification from hadronic observables for the MFV hypothesis.
Section III describes our technique for calculation of the baryon asymmetry in the ISS and LSS models with MLFV, while Sec. IV provides a description of our choice of scans and a discussion of the results. Here we also briefly assess the discovery potential of these models at current and future cLFV experiments. We conclude in Sec. V. 4 See [51] for an example of an MLFV Ansatz which incorporates leptogensis without resonance effects. 5 Recently it was found that flavour effects can still be relevant at much higher temperatures than thought previously, potentially reducing the allowed parameters space even further [15]. The Minimal Flavour Violation ansatz [48][49][50] recognises that, in the massless limit of The Abelian transformations are flavour-blind and we have the freedom to identify them with baryon number and two axial rotations [62]. The non-Abelian transformations do act in flavour space; they govern interactions between the different flavours and hence are responsible for the flavour violation in the theory. Table I  Unique transformations under the flavour symmetries are assigned to the spurions in order to make the would-be flavour-breaking terms invariant. This is a hypothesis motivated by the lack of experimental evidence for new flavour-violating physics; it should be treated as 6 In this work, promoting a Yukawa coupling to be a dynamical field will be represented by Y i → Y i .
V CKM is the quark mixing matrix measured by experiment and here, by the usual convention, we have chosen the VEVs of the spurion fields such that the down-type matrix is diagonal and the up-type is not.
The MFV ansatz is based on two assumptions about the high-scale, renormalisable theory: • The flavour symmetry present in the SM quark kinetic terms is a 'good' symmetry which the high-scale sector respects. The only sources of flavour symmetry breaking within the theory at the high scale are from the VEVs of the Yukawa coupling spurions. 7 The dynamics by which these VEVs are generated lie in the unspecified high-scale theory.
• The SM is an effective theory for which all renormalisable and non-renormalisable operators must respect both the gauge and flavour symmetries of the theory. All operators which are not formally invariant under the flavour symmetry are made so by insertions of the appropriate combinations of Yukawa fields.
A number of consequences follow from this ansatz.
SM effective field theory operators which describe flavour changing processes require insertions of spurion combinations in order to become flavour invariant. For example, the four fermion operator (Q L γ µ Q L )(Q L γ µ Q L ) is forced to be flavour preserving by the flavour symmetry. In order to get a related ∆F = 1 operator, a spurion insertion must be made:   [50].
Of particular importance when adapted to resonant leptogenesis, which will be discussed in more detail later, all terms in the lowest-order MFV Lagrangian will receive corrections to their couplings from spurion terms which transform in the same way as the basic spurion under the flavour symmetry. For example, in the higher-order term Y d Y † d Y d transforms in the same way as Y d and therefore under the ansatz has to be included. As Yukawa spurion VEVs are small in magnitude these corrections are generally not significant.

B. MFV vs Experiment
The motivation behind adopting the MFV ansatz is due to the strong limits placed on the mass scale of new physics from flavour-violating hadronic processes involving the first generation, as illustrated in table III. If new physics appears at scales lower than these bounds it would indicate that either flavour violation arises solely, or dominantly, due to the SM CKM matrix indicating an MFV-like theory, or that it could couple dominantly to the heavier generations which have less stringent bounds.
A key advantage of the former assumption is that strict relations between different flavour changing processes arise due to precise measurements of the CKM parameters. Therefore Hadronic Process Bound on Λ NP observables [64][65][66], some of which are summarised in table IV. They can be compared to their corresponding experimental measurement or limit [67][68][69][70][71] showing strong correlation.
Most channels have measurements (or upper limits) in agreement with MFV predictions, however a discrepancy exists for B d → µ + µ − between the CMS/LHCb [68] and ATLAS [71] collaborations. CMS/LHCb measure the rate for this process at the limit of the MFV prediction, whereas ATLAS place an upper-bound consistent with MFV. More precise measurements are required and could potentially suggest a deviation from an exact MFV theory.
Certainly, however, we have strong evidence that flavour violation is dominantly generated by the Yukawa couplings.

C. MLFV
While the MFV ansatz for the quark sector can be uniquely implemented (with strong agreement with measurement) there is an issue when extending this concept to leptons. Currently we do not know the origin of LFV and most importantly whether the light neutrinos are predominantly Dirac or Majorana. One may expand the MFV ansatz to the lepton sector, which is termed MLFV, but a dependence exists on the specific SM extension adopted for neutrino mass generation and a freedom exists in choosing which couplings are flavour vi-Hadronic Process Measured value or upper limit MFV prediction olating [72]. The simplest approach is to extend the SM by three right-handed neutrinos N R in the usual way where lepton number violation occurs through gauge-invariant Majorana neutrino mass terms, where L k ss is the usual kinetic term and The numbers e,D and µ N are to be distinguished from the matrices (in flavour space) Y e,D,N which will ultimately act as sources for the breaking of the Abelian and non-Abelian flavour symmetries respectively once promoted to spurions having non-zero VEVs.
Similar to the quarks and MFV, the gauge sector of the leptons obey a flavourtransformation invariance which can be defined analogously to eq. (1), G L A G L

NA
In the absence of L m ss , the Lagrangian is invariant under the following SU (3) rotations implying the following representations for the fields: Once again, flavour invariance can be recovered (once L m ss is reintroduced) by promoting the relevant couplings to spurions with fixed transformation properties under the flavour symmetries, implying the following representations for the fields Lepton masses and mixings appear once these spurionic fields acquire non-zero VEVs along with the SM Higgs. We are free to choose a basis such that It is clear from eq. (12) that it is not possible to assign unique background values to the spurions Y D and Y N in terms of physically measurable parameters. Rather, it is the combination in the seesaw formula that is fixed by the physical observables. As a consequence each individual spurion cannot be uniquely written in terms of physical masses and mixings.
In the quark sector, under the MFV ansatz, higher-dimensional operators which control rare processes are made flavour invariant from the necessary combination of spurion fields. As their background values are only dependent on quark sector masses and mixings, this prevents any sources of new physics from contributing in a way not aligned with the SM flavour violation. For the lepton sector, however, rare flavour-violating processes will not in general be made invariant by the exact combination Y D Y −1 N Y D , but through different combinations of these two spurions (see [47] for a list of examples). Therefore specific flavour-changing processes cannot be uniquely written in terms of physically measurable parameters in the same way they can be for MFV.
Predictivity can be recovered if only one of the spurion fields is taken to have non-trivial flavour transformations, with the other not acting as a source of flavour-symmetry breaking.
Usually, it is assumed that the Majorana mass term does not act as a source of lepton-flavour breaking and, and as in the quark sector, only the Yukawa couplings are responsible. If it is assumed that the Majorana mass term is lepton-flavour blind, then Y N must necessarily be equal to the identity matrix (in flavour space). Under this assumption the non-Abelian flavour symmetry of the theory 8 reduces to The lepton-flavour symmetry is broken by Y D and Y e which can now be written uniquely in terms of lepton masses and mixings The measured masses and mixings in the lepton sector fix the combination Y D Y T D in this setup and not the spurion Y D itself. However, often the combination Y D Y † D (which transforms as a (8,1,1) under G L ) is required to be inserted for many operators. Without the imposition of CP invariance, additional and yet unmeasured phases will appear in the predictions for flavour-changing processes. However, in general we expect CP violation (CPV) to be present in the lepton sector.
Allowing for leptonic CPV, as necessary for leptogenesis, a slightly less minimal flavour ansatz is required, namely Lifting the assumption of CP invariance introduces new non-SM phases, as well as allowing non-trivial phases in the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix. While they spoil the absolute predictivity of the theory, they generally lead to only small deviations in cLFV observables [52][53][54]. This is not in conflict with experiment (as it might be in the quark sector) as currently no cLFV processes have been measured. Additionally, although not definitive, there is interesting evidence for CPV in neutrino oscillations [74] motivating the less-minimal definition of the MLFV principle where phases are allowed and their effect on measurable processes are taken into account. redefinition. Under this assumption no degeneracy is implied amongst the right-handed neutrinos and CP invariance is not necessary. [73] As already emphasised, MFV and MLFV should be understood as hypotheses motivated by experiment and do not provide a complete UV description of the flavour sector. Attempts have been made in moving from this ansatz to a renormalisable model. In particular in the case of the simplest type-I MLFV Ansatz, the spurion scalar potentials have been studied [75][76][77] with suggestive conclusions. For hierarchical masses amongst the charged leptons, large mixing angles appear when Majorana neutrinos are considered, whereas in the quark sector small mixing angles are predicted. It has been suggested that the disparity in the mixing between the lepton and quark sector is therefore due to the Majorana properties of the right-handed neutrinos. These results suggest that UV completions of the MFV and MLFV Ansätze could ultimately explain the origin of the flavour structure of the SM.

D. MLFV and the inverse and linear seesaw models
It was noted in [46,53,58] that in order to generate measurable LFV effects whilst explaining the smallness of the active neutrino masses, a decoupling between the LFV scale and LNV scale is required such that Λ LFV Λ LN . This observation, coupled with any future measurement of cLFV, motivates the incorporation of such a scale separation into the MLFV hypothesis. Note that LNV arises from the breaking of Abelian symmetries in eq. (7), whereas LFV arises from the breaking of the non-Abelian symmetries.
Scale separation, however, does not occur in the minimal version of the type-I and type-III seesaw mechanisms where only one new scale is introduced such that Λ LFV = Λ LNV = µ N .
Under an MLFV Ansatz the flavour spurions decouple when the LNV scale is taken to infinity [58], preventing significant LFV. In contrast it was argued that the type-II seesaw does exhibit such a behaviour as where M ∆ is the mass of the triplet and µ ∆ is the dimensionful cubic coupling constant between the triplet and the SM Higgs doublet. This therefore defines the only minimal seesaw model for which an MLFV Ansatz may lead to significant LFV.
The simplest 9 fermionic completion which achieves a natural scale separation features in the ISS and LSS mechanisms (and a combination of the two) where small LNV parameters are introduced in order to ensure an approximate U (1) L symmetry. Here additional sterile neutrino degrees of freedom S L are introduced alongside the right-handed neutrinos N R , where it is conventional to choose L( L ) = L(N R ) = 1 and L(S L ) = −1.
As before, Y i correspond to dimensionless 3 × 3 matrices in flavour space associated with the breaking of G L NA . We impose an SO(3) rather than SU (3) flavour invariance for the heavy singlets to ensure that only the Yukawa interactions with the Higgs doublet act as sources of flavour-symmetry breaking. Therefore all bare mass terms related to the heavy singlets are proportional to the identity matrix. The flavour symmetry of the theory is defined analogously to eq. (7), The mass matrix has the form where for the ISS ( L = 0) and for the LSS (µ N , µ S = 0) The imposition of small lepton number violation means that µ S , µ N and L are small. This provides the desired separation of scales.
As with the quark sector and the example of the minimal type-I seesaw model, the kinetic terms in eq. (17) are invariant under the following flavour rotations where U i are 3 × 3 special unitary matrices, O i are 3 × 3 real-orthogonal matrices and for example L = (e L , µ L , τ L ). The leptons transform as fundamentals under their respective non-Abelian rotations It follows that eq. (18) can be made invariant by promoting the relevant couplings to dynamical spurions with the flavour transformation properties, implying that Lepton masses and mixings result as usual from non-zero VEVs for the spurion fields where, similarly to eq. (12), a basis can be chosen such that in the case of the inverse seesaw and in the case of the linear seesaw.

III. LEPTOGENESIS
The leading-order flavour degeneracy amongst the heavy SNs is built into the ansatz, with higher order spurion effects breaking the degeneracy. This allows for a resonant enhancement in asymmetry generation.
The case of the simplest MLFV scenario has been explored in detail [52][53][54]78]. Recent developments related to flavoured leptogenesis [55][56][57] have, however, highlighted the importance of a term not previously included, such that the flavoured CP-asymmetry is proportional to Equation 29 requires the existence of real, off-diagonal entries in Y † D Y D in order for non-zero asymmetry generation to appear at lowest order in the couplings. It turns out that this presents a problem, as we now review, and provides an independent motivation for the ISS and LSS extended seesaw schemes in the MLFV and leptogenesis context.
In the case of the MLFV Ansatz implemented in the type-I seesaw scenario, the radiative corrections to the SN Majorana masses are incorporated through where µ N is defined in eq. (6) and the SO(3) breaking source δµ N arises from all spurion insertions that transform such that the combination N c R δµ N N R starts off flavour invariant. From here on, Majorana masses with tildes will denote that corrections have been included whereas without will denote the SO(3) conserving bare mass. At lowest order they are where n 1 is a Wilson coefficient which at the effective level is a free parameter. The two terms appearing in eq. (31) can be checked to transform in the appropriate way by applying the rotations defined in eqs. (24) and (25) to the terms of eq. (31). For the MLFV Ansatz these terms must be included, but we remain agnostic as to how they are generated.
In order to calculate the CP-asymmetry one must first rotate the SN into their mass basis. The unitary matrix responsible for diagonalising the heavy SNs in this case is real (and therefore orthogonal): Due to the degeneracy required in Y N and the allowed flavour invariant terms of eq. (31) the rotation matrix O is exactly the matrix which diagonalizes Re Y † D Y D . Therefore the CP asymmetry for flavoured leptogenesis shifts to which has no off-diagonal terms and thus the CP asymmetry vanishes exactly. A nonzero asymmetry could potentially be recovered by including higher-order terms of O(Y 6 D ) of appropriate form. However, due to the suppression from the additional powers of couplings in the assumed perturbative regime, a substantial asymmetry will be difficult to generate even if resonance effects are present.
Alternatively, the consideration of higher-order corrections to the Majorana mass which are not all aligned in flavour space will cause non-zero real entries in the relevant term of eq. (29). Higher-order terms which transform in the same way as in eq. (31) can be constructed from the charged-lepton spurion, as per where, for clarity, only the terms not flavour aligned have been written explicitly. (The additional terms in eq. (34) will be fully written out later on.) Due to the inclusion of Y e , Changing to the mass basis now allows off-diagonal real terms and therefore the generation of a non-zero asymmetry. The necessary off-diagonal entries, however, are generated by the misalignment between the two spurions, and thus they will be suppressed compared to if they had been generated at leading order.
Due to the extra fields introduced, the above flavour alignment problem does not occur in general for the extended seesaw models. Consider the full ISS + LSS scenario where the equivalent radiative corrections to the relevant Majorana masses are included, Hereμ i = diag(. . . ) and the rotation matrices O and V diagonalise δµ N and δµ S respectively.
The spoiling of the flavour alignment should be clear: in the process of diagonalising the Majorana mass terms µ i , the Dirac mass between the two heavy SNs m R 1 3 has been rotated as well. Diagonalising the lower right 2 × 2 block will now no longer rotate the Dirac mass matrices m D and m L by the same orthogonal rotation that diagonalises Y † D Y D . The same alignment effect can occur only for very specific choices between the entries of eq. (35) and is not a general feature. We will provide an illustrative example below.

A. Baryon Asymmetry
In this section we briefly describe the procedure and formalism we employ for estimating the efficiency of asymmetry generation. A more detailed explanation of our conventions, definitions and numerical calculations can be found in [43]. These, in turn, are based on the Boltzmann equations (BEs) derived in [21,22] specifically for flavoured leptogenesis in the resonant regime.
As is conventional, we work in the heavy sterile neutrino mass basis, i.e. the bottom-right sub-matrix of Eq. (20) is real, positive and diagonal. The block diagonalisation of the lower right 2 × 2 block defines the rotation, leading to where in every second line we have expanded up to O( Once the spurionic radiative corrections are included in the Majorana self-energies, the O(3) 2 invariance is lifted and the matrices no longer commute. In the regime where these corrections are small these expressions will remain approximately valid. However, once the corrections are included, we perform the rotations numerically during our scan.

The combination
, is now the relevant combination for flavoured leptogenesis. The entries in the diagonal blocks of h † h control the CP asymmetry generated between two different generations of heavy SNs with same sign mass splitting The off-diagonal blocks control the CP asymmetry generated between heavy SNs of opposite mass splitting (same generation or otherwise) Clearly, in the absence of radiative effects where we have µ N , µ S ∝ 1 3 only ∆m OS i,j = 0, with ∆m SS i,j = 0. Once the corrections are included, however, the SO(3) 2 invariance is broken allowing for both ∆m OS i,j = 0 and ∆m SS i,j = 0. From eqs. (38) and (39) it can be seen that for example if m D = m L 10 along with µ N = µ S then m D = 0 but m L = 0. No real, offdiagonal entries would exist and the same alignment problem of the type-I seesaw would occur. This is not a general feature in the extended MLFV seesaw (as it is in the minimal type-I MLFV seesaw) and would require an additional symmetry to enforce the necessary relations between the couplings and masses such that this flavour alignment would occur.
We consider a low-scale scenario in which the sterile masses are set to O(TeV). This implies a low-temperature regime of lepton asymmetry generation where individual lepton flavours are in equilibrium with the thermal bath and distinguishable. We employ flavourdependent Boltzmann equations [22,79] for which we define the CP asymmetry generated from the decays of an SN, N i , to a specific lepton flavour α , These have been calculated previously [21,22] for the general case relevant for both hierarchical and degenerate scenarios with potential resonance effects [80] included. The result, separated into vertex ε V and self-energy ε S contributions, is and the loop function for the vertex-diagram contribution f (x ij ), is given by [6] f A resonant enhancement in the CP asymmetry will occur when which leads to the simplified condition to move from eq. (45) to (46).
We choose the regulator of ε i α to be of the form m i Γ j , which as discussed in Appendix A of [55] has consistent behaviour for models with approximate lepton-number conservation.
For other choices of regulators unphysical behaviour may be encountered in such a scenario.
For example choosing a regulator of the form (m i Γ i −m j Γ j ) diverges in the scenario m i → m j combined with Γ i → Γ j which occurs as the LNV parameters are taken to zero.
Specifically for the LSS scenario, it is important to note that the second term in eq. (43), due to the one-loop self-energy correction ε S , identically goes to zero when m N i → m N j .
By contrast, the first term arises from the one-loop vertex correction ε V . In the limit m N i → m N j the asymmetry it generates is non-zero. However, A α ij = −A α ji and therefore for m N i = m N j . Therefore the asymmetry produced by the decay of N i to flavour α is almost equal and opposite to that of N j and in the limit h † h ii → h † h jj the asymmetry will vanish identically once the decays of all SNs are included. If the two SNs have different decay widths a non-zero asymmetry is possible, albeit highly suppressed.
We will work under the condition that, due to the strong washout nature of the temperature regime we favour [43], asymmetry generation occurs predominately well before the electroweak phase transition crossover and so we need not consider the changing sphaleron rate as the temperature approaches and crosses its critical value. The baryon asymmetry is then expressed as a fraction of the lepton asymmetry generated through where the factor of 28/51 arises from the fraction of lepton asymmetry reprocessed into a baryon asymmetry by electroweak sphalerons [81], while the dilution factor 1/27 arises from photon production until the recombination epoch [82].

IV. NUMERICAL RESULTS
In full analogy to the type-I scenario, all appropriate flavour invariant operators must be included. This leads to corrections to the Majorana mass terms from spurion combinations transforming the appropriate way such that when coupled with the heavy SN fields, the term is flavour invariant at the high scale. Now in the case of the extended seesaw both Majorana masses receive corrections, as per As can be checked explicitly by applying the transformations defined in eq. (25), the combinations of spurions above are flavour invariant when coupled to N c R N R and S L S c L . All terms at next-to-leading order in the flavout-invariant spurion insertions will be explicitly written below.
The coefficients n i and s i are dimensionless Wilson coefficients which are treated as free parameters in the absence of an explicit high-scale, renormalisable theory. They are conventionally either taken to be O(1) numbers [46,47,52,83], or arising from radiative effects [53,54] so that n 1 s 1 1 16π 2 and s We consider the general scenario where one copy of N R and S L is added for each generation of active neutrino. The tree-level light neutrino mass matrix can be expressed in powers of the LNV parameters [84] which in complete generality (without assuming MLFV) can be expanded to where for completeness we have written all leading and next-to-leading terms valid for ||µ S ||, ||µ N ||, ||m L || ||M R ||. The above equation is valid for the general case of noncommuting sub-matrices, but for our MLFV setup M R = m R 1 3 , m L is diagonal and In the limiting case of the ISS, only the first line of eq. (50) contributes at leading order to the active neutrino masses. This mass vanishes in the limit µ S → 0 even if µ N = 0. Even for µ S = 0 the contribution from µ N is highly suppressed, including for a large regime where At one-loop order, additional terms are generated for the active neutrino masses which can be important [84], where is a one-loop function and eq. (51) is only valid for M R = m R 1 3 . Different terms in eq. (51) can contribute significantly for different hierarchies amongst µ N , µ S and m L which are free parameters in our scan.
Combining the tree-level and one-loop contribution at leading order leads to a general light-neutrino mass matrix in our MLFV Ansatz of the form one-loop For both the ISS and LSS we fix m R = 1 TeV. For the ISS specifically, the LNV parameters are randomly (and independently) scanned over the ranges where we allow for larger values of µ N over µ S as its contribution to the active neutrino mass is suppressed by a loop factor. For the LSS we vary where, in order to reduce the number of degrees of freedom, we take it to be diagonal but not proportional to the identity matrix in order for its inclusion to break the flavour symmetry.
We fit active-neutrino data by fixing the Dirac mass matrix m D with an approximate 11 Casas-Ibarra parameterisation [85] which for the ISS is where we define an effective mass For the LSS, with an equivalently defined effective mass The matrix R is a complex-orthogonal matrix R R T = 1 3 which can be parameterised with three complex mixing angles and we have defined cos(x) ≡ c x and sin(x) ≡ s x in the usual way.
In order to reduce the number of free parameters we fix the real components of the mixing angles to where there is no signficance to the values chosen. The asymmetry is not sensitive to the values of the real angles [43] and therefore should apply for any choice of their values. We fix θ c 1 = θ c 3 = 0 for simplicity and scan over 11 The parameterisation is only approximate as we ignore the corrections to the mass terms e.g.μ i = µ i +δµ i generated through spurion insertions when solving for m D . This approximation is valid in the regime where δµ i < µ i . We additionally check that once these corrections are included the active-neutrino mass differences are not spoiled.
By contrast to the complex-orthogonal R, the matrix C instead satisfies C + C T = 1 3 and therefore must be an anti-symmetric matrix of the specific form, with a i = a r i + ia c i where we choose to fix and scan over The active-neutrino masseŝ are inputs for eqs. (56) and (58) where for simplicity we assume normal ordering which is currently favoured over inverted ordering [74], implying and unless stated otherwise we fix the lightest neutrino mass m ν 1 to 0.01 eV. Other active neutrino parameters are fixed to their current best fit values [74] and listed in table V, while the Dirac phase has been fixed to be maximally CP-violating.
We work in the perturbative regime of the Yukawa couplings generated by the Casas-Ibarra parameterisations such that |Y ij | 2 < ∼ 4π, which prevents very large choices for the complex parameters. We fix the unknown Wilson coefficients s i and n i to 1 16π 2 for definiteness unless stated otherwise.
To check our numerical solutions to the Boltzmann equations, we compare the results to a known approximate analytic expression of the baryon asymmetry that is valid in the strong washout regime (K α eff ≥ 5) [ [74].
with m N i = (M ± ) ii . The asymmetry is approximated by η approx. .
The numerical versus analytic approximation comparison illustrated in table VI provides strong evidence that our numerical routines are accurate. where for these example parameters we have set all Wilson coefficients at lowest order to 1/16π 2 .

A. Inverse seesaw (ISS)
For the ISS (Y L = 0,) the only LNV parameters present are the two small Majorana From eq. (49) only µ N receives radiative corrections from SO(3) N R breaking terms proportional to the Yukawa spurion Y D .
We separately consider three scenarios for the ISS: first where no corrections are intro-duced, second where corrections are introduced at the next lowest order, labelled N 1 and, third, corrections up to next-to-leading order, labelled N 2 : where and we fix n (i) 2 = (n 1 ) 2 for computational convenience. Due to the perturbative regime in which we operate, it will always be the case that ||N i || 1. All terms are formally flavour invariant and serve to break the flavour degeneracy amongst the heavy SNs when these spurions acquire non-zero VEVs.
For the ISS as m L = 0 the combination in eq. (39) simplifies to where we are now dealing with unprimed matrices. Due to the fact that M R ∝ µ S ∝ 1 3 , in the case where only N 1 is included the corrected Majorana mass µ N can be diagonalised without affecting the other entries in the 2 × 2 sub-block, with under the rotation N R → ON R and S L → OS L . Similarly to the minimal type-I scenario, if only N 1 is included then Re m † D m D ∝ N 1 = δµ N and the combinations of m † D m D appearing in eq. (75) will have no real off-diagonal components. We emphasize that this only occurs for the ISS because m L = 0. Differently to the type-I scenario however the off-diagonal blocks of eq. (75) have an additional complex phase such that the imaginary This therefore allows for non-zero generation of lepton asymmetry unlike the type-I scenario.
As only the off-diagonal blocks of h † h contain the necessary terms, only diagrams between two SNs of different mass splittings, i.e. from ∆m OS i,j , will contribute to asymmetry. The real components of the (1, 1) and (2, 2) blocks are diagonal, on the other hand, meaning that SNs with the same sign mass splitting ∆m SS i,j will not contribute to ε i α . Once N 2 is included however, the flavour mis-alignment between δµ N and m † D m D will generate real off-diagonal components in m † D m D and therefore all four blocks of eq. (75) will contribute to ε i α . In fig. 1 we plot the baryon asymmetry numerically calculated as a function of both µ N and µ S . All three scenarios are simultaneously plotted. Both in the case where radiative effects are ignored as well as when only N 1 is included, similar behaviour occurs. The inclusion of N 1 breaks the mass degeneracy between all six SNs, as opposed to without its inclusion where two groups of identical mass SNs form. However, as discussed above, due to flavour alignment only diagrams involving opposite mass split SNs contribute to the asymmetry. Therefore the mass difference is generated exclusively from with While the spurion corrections are present we have µ N N 1 and therefore they are not the dominant source for the mass splittings in this situation. In this regime resonant leptogenesis is not feasible [43][44][45]. The inclusion of N 2 generates non-zero, off-diagonal entries in the (1, 1) and (2, 2) entries of eq. (75) such that new contributions turn on arising from the mass splitting We emphasis that even when only N 1 is included, the mass degeneracy between the same sign mass split SNs is still broken ∆m SS i,j = 0. No resonant enhancement occurs between two such SNs in this case not because of a mass degeneracy between them (as is the case when no corrections are included) but because of the flavour alignment issue described above.   it is also scaled by the Majorana mass µ N . For small LNV the mass splitting will be much less than the decay width but for large values forcing a resonant enhancement to occur. In fig. 2 we distinguish between two different hierarchical cases for the Majorana masses: µ N > µ S in purple and µ S > µ N in green.
The asymmetry in the resonant regime is independent of µ S as in our setup this mass does not receive any radiative corrections. If a stricter setup was considered such that it did  washout of any generated asymmetry. Figure 5 plots the washout as a function of the two Majorana masses for the hierarchy µ N > µ S . Both the naïve washout K α ∝ Γ i m i /H(m i ) and the effective washout K eff α defined in eq. (68) (and relevant for scenarios with small LN violating parameters [82,86,87] due to the sensitivity to 2 ↔ 2 scatterings) are plotted together.
While the naïve washout grossly overestimates the efficiency of washout, it is clear that for GeV induced by ∆m OS i,j independent of the radiative spurion contributions, in agreement with [43]. However, this region is accompanied by a much larger effective washout and cannot accommodate sufficient asymmetry generation. In contrast only when both N 1 and there is a slight dependence on the size of θ c 2 which can vary the asymmetry generated by one or two orders of magnitude. For small values of θ c 2 the asymmetry is generated by the CP-violating phase δ CP within the low-energy mixing matrix. Similar behaviour occurs when N 2 is included, with an overall preference for the range 0.1 < θ c 2 < 1. In the case where only N 1 is included, we confirm its strong dependence on the angle θ c 2 similarly to the type-I scenario [52]. However, we note that in our case the asymmetry is not zero exactly unless θ c 2 = 0. Asymmetry generation can only occur for non-zero values of θ c 2 independent of any low-energy CP violating phases as phases from U PMNS cancel in the term m † D m D in combination with the flavour alignment described above. Leptogenesis is viable both through the Dirac phase as well as high scale CPV but we note that within MLFV-ISS the region in which it is possible is very narrow.

B.
i → j γ and MLFV-ISS Here we briefly discuss the consequences of MLFV-ISS on low-energy cLFV processes, specifically the impact that the introduction of CPV has on predictions for i → j γ. We also Here we have fixed µ N 950 GeV and µ S 10 GeV to be on resonance. In blue no radiative effects are included, in red only to leading order and in green next-to-leading order effects are included. The behaviour of the asymmetry as a function of the complex angle is clear. At small values it freezes out to that provided by the Dirac phase δ CP however if only the leading radiative effects are included (red), the asymmetry tends to zero with decreasing θ c 2 even for non-zero δ CP . This effect is similar to what occurs in the type-I MLFV scenario [52,54,88]. Asymmetry generation is possible both with low-energy phases as well as complex entries of the R matrix where we loosely find the criterion 0.1 < ∼ θ c 2 < ∼ 1 for sufficient asymmetry generation when δ CP = 0.
suppressed. These operators are not invariant under the flavour symmetry defined for MLFV in eq. (19). Insertion of spurion combinations transforming as (3, 3, 1, 1) make these terms formally invariant. At lowest order this is simply the spurion Y e . However, as we work in a basis where this matrix is diagonal, it does not contribute to flavour-violating processes.
The lowest order spurion combination that makes the relevant operators flavour invariant and contains non-diagonal entries is ∆Y e , where The combination transforms as ∆Y e → U L ∆Y e U † e R as required. A simple expression for the branching ratio is obtained [46] in the limit m j m i , where C accounts for a combination of Wilson coefficients of the relevant operators. As is conventional when considering cLFV in MLFV, we set |C| = 1. Note that in a specific UVcomplete model based on MLFV these parameters will be fixed by the high-scale dynamics and may have a different magnitude.
In eq. (81) the only parameter which carries family indices is the spurion combination ∆ ij . Therefore useful predictive parameters for MLFV are ratios of branching ratios for different lepton flavour decays, In these ratios, the unknown scale of LFV (Λ LFV ) cancels out. While it can be identified with M R in the ISS model, it is also possible that it could arise as the result of some unknown dynamics not directly related to the seesaw mechanism.
Many detailed explorations of cLFV in MLFV have been made for the minimal type-I seesaw scenario [46,47,53,83], the results of which should also hold in the ISS. Here we briefly consider the impact of the complex angle θ c 2 and the Dirac phase δ CP on the predictions for specific ratios of cLFV observables. In particular we compare the CP-conserving scenario (assumed in the simplest version of MLFV) with the CP-violating scenario and the size of deviation their inclusion introduces. As we have demonstrated that leptogenesis is viable with both the low-energy phase δ CP and with the inclusion of the CP-violating angle θ c 2 , we analyse the two limiting cases where CPV arises solely from one of these angles. Figure 9 plots three ratios of cLFV observables as a function of the lightest neutrino mass m ν 1 for θ c 2 = 0 while varying δ CP . The red and orange lines correspond to the CP-conserving scenarios δ CP = 0 or π, which allows for full reconstruction of Y D from low energy observables. The blue shaded region corresponds to all other values of δ CP where the green line Here the complex angle θ c 2 has been switched off but the low-energy phase δ CP is varied. Lines in red and orange correspond to the CP-conserving cases δ CP = 0 and δ CP = π respectively. In green is the maximally violating case of δ CP = π/2 or δ CP = 3π/2. The shaded blue region corresponds to δ CP = (0, 2π)\{π/2, π, 3π/2}. All results are in full agreement with [47] for the CP-conserving cases. specifically corresponds to the maximally CP-violating cases δ CP = π/2 or 3π/2.
From eq. (49) as long as Y L ∝ 1 3 there is a radiative contribution to both the (2, 2) and (3, 3) entry which will break the SO(3) 2 invariance, leading tõ where As before, we consider separately the cases where no corrections are included, corrections are included at leading order and corrections are included at next-to-leading order. Here the dimensionful parameters µ N and µ S cannot be identified with bare Majorana mass terms. Rather, they arise from the unknown UV complete dynamics. As we remain agnostic about these dynamics, yet under the MLFV Ansatz we must include such terms, we choose to fix the values of µ N and µ S in some plausible way. For example, it seems reasonable to take µ i v 2 /M X , which would be true if µ i arose from some effective coupling with the SM Higgs doublet mediated by a heavy new particle X, generating an operator similar to the Weinberg operator.
We fix these parameters to µ N = µ S = 1 GeV which, in the example given above, would arise if m X O(10 − 1000) TeV depending on the strength of the relevant couplings. We also make the reasonable assumption that this mass scale is independent from m L . the self-energy contribution to the CP-asymmetry is switched off but a non-zero asymmetry is generated from the vertex contribution and differences in the decay widths of each SN.
The most asymmetry (albeit much too small) is generated in the region where the washout is minimised, as shown fig. 13, in agreement with [43].
Points in orange correspond to spurion insertions at lowest order and points in cyan include all relevant terms. Now, due to the inclusion of the radiative Majorana masses, mass splittings occur between the six heavy SNs. The self-energy component of eq. (43) turns on and increases the overall CP-asymmetry (due to a resonance between the SN masses) by several orders of magnitude. Here, including the next-to-leading order contributions does not change the size of the asymmetry generated. By contrast with with MLFV-ISS, the flavour mis-alignment between m D and m L generates non-diagonal real entries in eq. (39) at lowest order which also allows for larger values of ε α i . Figure 14 plots the CP-asymmetry into a specific lepton flavour α in the three scenarios, where once again due to the anarchic nature each flavour has the same behaviour and overall size. It is clear that, for all values of m L , the CP-asymmetry is orders of magnitude larger when the radiative Majorana masses are included. In other words, for the entire region of the parameter space, roughly the same resonant enhancement is occurring. This forcedresonance is occurring as the mass splitting between the heavy SNs is related to the small Majorana parameters, Contrary to the case of MLFV-ISS, corrections at leading order are sufficient to generate the necessary asymmetry. The case without any radiative corrections is highly suppressed as the self-energy contributions to ε α i are identically zero. However a non-zero asymmetry is generated due to the differences in the decay widths generating slight deviations in the vertex contribution for each SN, c.f. eq. (47). leading to where we take the parametersμ N andμ S to be independent of m L . Therefore the same level of resonant enhancement occurs for the entire parameter space. Unlike with MLFV-ISS, both ∆m OS i,j = 0 and ∆m SS i,j = 0 occurs at lowest order in the corrections and they both contribute to the asymmetry generated.
To illustrate this, fig. 15 plots the mass splitting ∆m ij and the decay width Γ i as a function of m L . The parameter which most significantly impacts the level of resonance (where the mass splitting overlaps with the decay width) is the combination of µ N n i and µ S s i . The resonant enhancement is maximised when Points in green correspond to when N i and S i are not included and points in purple correspond to when they are included (at lowest order and next-to-leading order the points are identical). Clearly a resonance due to their inclusion increases the overall CP-asymmetry by many orders of magnitude. The mass splittings induced in this scenario are directly proportional to the decay widths and therefore for any choice of m L the mass splittings will always be in a resonant regime.
three LNV parameters based on the two scenarios here will allow for minimised washout with mass splittings related to the heavy SN decay widths for the necessary resonance to occur. However, as resonant leptogenesis is already feasible [43] in this scenario without the need for radiative resonant leptogenesis, we do not consider this scenario further.
Finally, in fig. 17 we plot the behaviour of the asymmetry as a function of the complex component of the C matrix when low-energy CPV is included and when it is not. Similarly to the ISS case, asymmetry generation can be predominately due to either the Dirac phase δ CP or the complex component a c 2 of the C matrix. We find consistent behaviour for the baryon asymmetry as these CPV parameters are taken to zero. It is clear from both figs. 12 and 17 that a larger portion of the parameter space provides the necessary asymmetry generation allowing for smaller sizes of the CPV parameters, decreasing their contribution to flavour-violating processes.
GeV. This is plotted against the decay width Γ i,j in brown. As the combination N 1 + S 1 ∝ Γ i,j the mass splitting induced will always be on resonance independent of the value of m L . Here there is no distinction between ∆m SS i,j and ∆m OS i,j and they both behave in a similar way.

D.
i → j γ and MLFV-LSS Once again we briefly consider the prospects of detection of MLFV-LSS through cLFV processes. Unlike in the MLFV-ISS scenario a much less tuned region of parameter space is required in order to generate the necessary baryon asymmetry which may lead to improved detection prospects.
As before we require insertions of spurions transforming as a (3, 3, 1, 1) in order to make the necessary dimension six effective operators invariant. Off-diagonal terms are required in order for LFV processes to occur. The lowest order combination which satisfies this is once While the combination Y L Y † L may transform in the correct way, it does not contain the necessary off-diagonal terms in our scan and therefore cLFV will once again be controlled all radiative spurion effects are included. In this scan we fixed a c 2 = 0.7 and in the left plot set n i , s i = 1/16π 2 whereas we fixed m ν 1 = 0.01 eV in the right plot. In both plots similar behaviour compared to MLFV-ISS is found. Figure 18 plots the three ratios of cLFV observables as in the ISS case for the scenario where a c 2 = 0 while δ CP is varying. Similar predictions on the ratio R (µ,e) [τ,µ] are made compared to the ISS scenario and therefore the branching ratio BR(τ → µγ) should be larger by one or two orders of magnitude compared to BR(µ → eγ). Here there is a cancellation in τ → eγ for specific values of the lightest neutrino m ν 1 and the phase δ CP . In these regions a strong suppression of this channel is predicted. Outside the regions of strong cancellation the LSS similarly predicts R (τ,e)[τ,µ] < 1 but always predicts R (µ,e)[τ,e] > 1 unlike the type-I and ISS case where this varied depending on the value of δ CP . As before the introduction of CPV appears to bring the ratios closer together. Figure 19 plots the ratios of cLFV observables in the scenario where δ CP = 0 and a c 2 is varied. We find that overall the presence of the CPV parameter becomes more impactful for larger values of the lightest neutrino mass, while for a hierarchical spectrum its impact is less significant. Overall similar predictions to the case where δ CP was varied are obtained and a significant portion of the parameter space is not strongly sensitive to the presence of CPV. This implies the inclusion of CPV necessary for leptogenesis does not significantly modify the predictions given by the CP-conserving case.
Finally, fig. 20 plots the predictions for the branching ratio BR(µ → eγ) as a function of the LNV parameter m L for Λ LFV = 1 TeV. As before, we plot cases with no CPV violation Here the complex parameter a c 2 has been switched off but the low-energy phase δ CP is varied. Lines in red and orange correspond to the CP-conserving cases δ CP = 0 and δ CP = π respectively. In green is the maximally violating case of δ CP = π/2 or δ CP = 3π/2. The shaded blue region corresponds to δ CP = (0, 2π)\{π/2, π, 3π/2}. scale for a given LFV scale to be probed. A future measurement of the process µ → eγ may rule out LSS as a leptogenesis candidate in the absence of additional physics introduced to lower the overall strength of the washout present for smaller values of m L .

V. CONCLUSION
We have studied a well-motivated way in which small mass splittings between heavy SNs from different families may arise, within both the ISS and LSS frameworks, such that leptogenesis is possible despite the strong washout present in the theory. Previously it was found that while a mass splitting naturally exists for the ISS it is not sufficient in order for resonant leptogenesis to occur. For the LSS the degeneracy amongst the SNs at the high scale prevents significant asymmetry generation. While leptogenesis is feasible when all LNV terms are switched on (the ISS+LSS case) we explore the potential for ISS or LSS leptogenesis to occur independently, where additional symmetries may prevent both terms from existing.
In the context of broken flavour symmetries and the MLFV hypothesis, a degeneracy amongst the heavy SNs is naturally produced, for the purposes of having a predictive theory. The degeneracy is then broken by higher-order spurion VEV contributions, leading to a parameter region consistent with resonant asymmetry generation. In order for the desired splitting the occur in the intended way during cosmological evolution, the critical tempera-  [89,90] respectively for this decay mode.
Currently very small values of the LNV parameter are probed for MLFV-LSS which correspond to a very large seperation of the LNV and LFV scales. In order to allow for necessary asymmetry generation we require approximately 10 −5 < ∼ m L /GeV < ∼ 10 −3 which will not be probed by MEG-II in the near future. The inclusion of CPV of any type does not significantly impact the predictions. ture at which the spurions acquire their non-zero VEVs must be assumed to be above the scale of thermal leptogenesis.
We found that for MLFV-ISS only a small region of very large Majorana masses is able to generate the required asymmetry. Here asymmetry generation requires next-to-leading order corrections to be included, therefore suppressing the overall size of the CP-asymmetry generated per decay of SN. We briefly discussed the impact of CPV on potential low-energy observables, in particular how various ratios of cLFV decays are impacted compared to the CP conserving MLFV scenario. The region compatible with leptogenesis will not, however, be probed by current and near future cLFV experiments.
For MLFV-LSS a large region of parameter space is capable of satisfying the resonance condition simultaneously with the minimised washout required for successful asymmetry generation. Here corrections at lowest order are not flavour aligned and therefore much larger values of the CP-asymmetry are possible compared to the ISS. Similarly, we studied the impact of CPV on cLFV observables and find that small deviations occur due to their inclusion. In this case, relatively small values of the CPV parameters allow for sufficient asymmetry generation allowing for even smaller deviations as compared to the ISS case.
In both cases we estimated the impact of the lightest neutrino mass m ν 1 on the asymmetry generated. We find a clear preference for small values where the light neutrinos are hierarchical and estimate that m ν 1 < ∼ 10 −2 eV is required. Unsurprisingly we find that MLFV leptogenesis favours larger values for the Wilson coefficients related to the Majorana mass corrections in both scenarios.

ACKNOWLEDGMENTS
This work was supported in part by the Australian Research Council. TPD thanks Vincenzo Cirigliano and Luca Merlo for clarifications related to MLFV. TPD is grateful to Uli Felzmann for access to computing hardware on which our scans were conducted.