1 Introduction
###### Abstract

We revisit charged lepton flavour in-flight conversions, in which a beam of electrons or muons is directed onto a fixed target, , and , focusing on elastic interactions with a nucleus . After a general discussion of this observable, we carry a full phenomenological analysis in the framework of minimal Standard Model extensions via sterile neutrinos, with a strong emphasis on the rôle of the increasingly more stringent constraints arising from other (low-energy) charged lepton flavour violation observables. Despite the potential interest of this observable, in particular in the light of certain upcoming facilities with the capability of very intense lepton beams, our study suggests that due to current bounds on three-body decays () and conversion in Nuclei, the expected number of conversions in such a minimal framework is dramatically reduced. An experimental observation of such a conversion would thus signal the presence of another source of flavour violation, possibly at tree-level.

LPT Orsay 16-86

PCCF RI 16-08

IFT-UAM/CSIC-16-138

FTUAM-16-47

July 29, 2019

In-flight cLFV conversion: , and in minimal extensions of the Standard Model with sterile fermions

A. Abada, V. De Romeri, J. Orloff and A.M. Teixeira

Laboratoire de Physique Théorique, CNRS,

Univ. Paris-Sud, Université Paris-Saclay, 91405 Orsay, France Â§

Departamento de Física Teórica and Instituto de Física Teórica, IFT-UAM/CSIC,

Laboratoire de Physique Corpusculaire, CNRS/IN2P3 – UMR 6533,

Campus des Cézeaux, 4 Av. Blaise Pascal, F-63178 Aubière Cedex, France

## 1 Introduction

The quest for a Standard Model (SM) extension capable of addressing its several observational caveats has fuelled intensive experimental searches, encompassing high-energy colliders, high-intensity facilities, as well as numerous astroparticle and cosmological searches.

So far, no direct evidence for the new states has been unveiled in collider searches, and this has in turn intensified the interest for the so-called indirect searches, in which very rare processes, strongly suppressed or even forbidden in the SM, are looked for. Among the many observables that are being studied and explored, those signaling the violation of lepton flavour are powerful probes of New Physics (NP), sensitive to new scales often lying well beyond collider reach. Numerous processes are currently being searched for in high-intensity facilities, and these include charged lepton flavour violating (cLFV) radiative decays, three-body decays and nuclear-assisted transitions; likewise, a vast array of rare transitions and decays is being looked for at high-energy colliders111Other rare processes, such as those violating lepton flavour universality (LFUV), or total lepton number (LNV), can also emerge in relation with cLFV transitions and decays.. The current bounds are already impressive, and many running and/or upcoming experiments (as is the case of MEG II, Mu3e, Mu2e, COMET and LHCb) should improve them in the near future. In Table 1 we summarise the present experimental bounds and future sensitivities for several radiative and 3-body cLFV decays (which will be relevant for our subsequent discussion).

The probing power of cLFV has been at the source of an increasing interest for these processes, leading to further explorations of already existing observables, or to the study of new ones. This was the case of the Coulomb enhanced decays of a muonic atom into two electrons [13, 14], or the lepton flavour and lepton number violation conversion in Nuclei [15, 16, 17, 18].

In the wake of the discovery of oscillations - and of large mixing in the neutral lepton sector - the study of cLFV lepton production in (with denoting a generic nucleon) at high energies [19] was originally proposed. The experimental signature for the cLFV in-flight conversion would be that of a final state composed by a single muon (the tau, despite its large energy, rapidly decaying, ), with a dramatic loss in energy when compared to that of the primary muon beam - the energy loss on target corresponding to the production and subsequent decay of the heavier lepton. First studies also focused on the quasi-elastic in-flight conversion, due to the simpler final state topology and to the associated background. The possibility of having high-intensity (and sufficiently energetic) muon beams (for instance at muon and future neutrino factories) further fuelled the interest for such cLFV observables; as argued in [20], a 50 GeV muon beam, with an expected intensity of muons on target per year could lead to a significant number of events. The original estimation was based on an effective approach, and preceeded the recent stringent bounds on cLFV transitions (many of them collected in Table 1).

Other pioneering studies of in-flight cLFV conversion focused on leptoquark models [19, 21], also highlighting the potential of flavour violating constructions such as -parity ()-violating supersymmetry, or flavour-violating Higgs interactions. Following the model-independent approach of [20], the prospects of supersymmetric extensions of the SM for were discussed [22] in the deep-inelastic scattering (DIS) regime, as were those of low-energy electron-nucleus scattering to probe conversion [24]. In Ref. [25], the impact of massive neutrinos was first considered; contributions arising in the framework of a typical type I seesaw (albeit for low right-handed neutrino masses) were found to enhance those emerging from the presence of three light massive Dirac neutrinos, assuming a CKM-like lepton mixing, by as much as twenty orders of magnitude in the case of dominating photon contributions. Other studies in the DIS regime focused on cLFV conversions induced by “unparticles” [26]. Recent analyses, again based on an effective-Lagrangian approach, included a detailed discussion of the process’ kinematics and hadronic contributions [27]. Associated experimental issues (including a brief overview of backgrounds), and future prospects were discussed in [28].

In view of recent phenomenological and experimental developments, which have led to increasingly severe bounds on the scale of NP mediators and to strong constraints on the strength of possible cLFV couplings, a re-analysis of the in-flight cLFV conversion - and its potential impact on SM extensions - is clearly justified. Expected experimental prospects (such as the capability of high-intensity, high-energy muon beams [22, 23], or a possible electron-ion collider [21]), further motivate revisiting this observable.

Although not necessarily linked to the problem of neutrino masses and mixings (which signal the violation of neutral lepton flavours), cLFV can also emerge in association with SM extensions incorporating a mechanism of neutrino mass generation. Minimal extensions of the SM via additional sterile fermion states are an appealing class of models, in particular those that succeed in explaining oscillation data by the introduction of (not excessively) heavy states. Numerous studies have examined the impact of these models regarding several cLFV observables [29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41], focusing either on specific realisations, or then evaluating the potential contributions of sterile fermions via model-independent, simple constructions (the so-called “3+N” models). Among the many theoretically complete frameworks which simultaneously explain neutrino data, while at the same time having a significant phenomenological impact, one finds several low-scale seesaw models, such as variants of a type I seesaw, the linear seesaw [42, 43], the Inverse Seesaw (ISS) Seesaw [44] or the neutrino minimal SM (MSM) [45].

In the present work we thus revisit cLFV in-flight conversions , carrying a full phenomenological analysis in the framework of SM extensions via sterile neutrinos. In particular, we focus on flavour violating (FV) - and photon-mediated interactions, recomputing their contributions, and comparing our results to previous studies. We consider the three different cLFV channels (, and ), and discuss the corresponding experimental prospects, confronting the latter with other cLFV observables. In this study, we consider quasi-elastic scattering for the in-flight cLFV observables, which offers a first estimate due to the simple final state topology and to the associated background. Although we do discuss the potential of well-motivated low-scale seesaw models (in particular the Inverse Seesaw and the MSM), a first phenomenological approach - and a significant part of the discussion - is done by means of an “ad-hoc” construction, a simple “3+1 toy model”, in which a single massive Majorana state is added to the SM content, with no hypothesis on the underlying mechanism of mass generation.

Our work is organised as follows: after describing the underlying theoretical framework in Section 2, we discuss in-flight conversions, including contributions to the differential cross section and general features of the observables. The experimental prospects, as well as a comparative study with other cLFV observables in minimal SM extensions via sterile neutrinos are collected in Section 4; a brief overview as well as further elements of discussion are summarised in the Conclusions. The relevant expressions of the - and -mediated interactions, together with other relevant form factors, can be found in Appendices A and B.

## 2 Minimal SM extensions via sterile fermions

Motivated by several cosmological and experimental observations, sterile fermions are present as constituent blocks of many SM extensions which encompass a mechanism of neutrino mass generation. If on the one hand sterile neutrinos can indeed provide an explanation to the problem of neutrino masses and mixings, they can also open the door to a rich phenomenology, with potential effects in a large number of observables. This is a direct consequence of their mixings with the light (mostly active) neutrinos, which - if non-negligible - lead to the violation of lepton flavour in both neutral and charged leptonic currents [46, 30].

In the presence of additional sterile (Majorana) neutrinos, the vector and scalar currents are modified as follows222Likewise, the interactions with neutral and charged Goldstone bosons are also modified: ; . (working in the physical basis, i.e., for mass eigenstates):

 LW±=−gw√2W−μ3∑α=13+nS∑j=1Uαj¯ℓαγμPLνj+H.c., LZ0=−gw2cosθwZμ3+nS∑i,j=1¯νiγμ(PLCij−PRC∗ij)νj−gw4cosθwZμ3∑α=1¯ℓαγμ(CV−CAγ5)ℓα, LH0=−gw2MWH03+nS∑i,j=1Cij¯νi(PRmi+PLmj)νj+H.c.. (1)

In the above, denotes the weak coupling constant, , , and are the physical neutrino masses (light and heavy); the indices denote the flavour of the charged leptons, while correspond to the physical (massive) neutrino states. In addition and are the SM coefficients parametrizing the vector and axial-vector -couplings of charged leptons, and . Finally, a rectangular mixing matrix, , parametrizes the mixing in charged current interactions (corresponding to the (unitary) PMNS matrix, in the case of ); the mixing between the left-handed leptons corresponds to a block of , usually denoted . The structure of is at the source of lepton flavour violation in neutral currents, which, as seen from above, is now parametrized by

 Cij=3∑α=1U∗αiUαj. (2)

### 2.1 Constraints on sterile fermions

Due to the presence of the additional sterile states, the modified neutral and charged lepton currents might lead to new contributions to a vast array of observables, possibly in conflict with current data. These SM extensions via sterile fermions must be then confronted to all available constraints arising from high-intensity, high-energy and cosmological observations.

In our subsequent phenomenological analysis, and for the theoretical framework considered, we ensure that compatibility with the following constraints - theoretical (such as perturbativity of the active-sterile couplings) and experimental - is verified at all times.

Sterile states, with a mass above the electroweak (EW) scale, can have sizeable decay widths, a consequence of being sufficiently heavy to decay into a boson and a charged lepton, or into a light (active) neutrino and either a or a Higgs boson. One thus imposes the perturbative unitarity condition [47, 48, 49, 50, 51, 52], . Noticing that the leading contribution to is due to the charged current term, one obtains the following bounds [47, 48, 49, 50, 51, 52]:

 m2νiCii<2M2Wαw(i≥4), (3)

where , and is given in Eq. (2).

Observational constraints on the sterile masses and their mixings with the active states arise from an extensive number of sources. Firstly, and other than requiring compatibility between the left-handed lepton mixing matrix and the corresponding best-fit intervals333We do not impose any constraints on the (yet undetermined) value of the CP violating Dirac phase . defined from -oscillation data [53, 54, 55, 56, 57, 58, 59], we also impose, when relevant, unitarity bounds as arising from non-standard neutrino interactions with matter, on the deviation of from unitarity [60, 61, 62]. Further constraints on the active-sterile mixings (and on the mass regime of new states) arise from electroweak precision observables; these include new contributions to the invisible -decay width (addressed in [63, 65, 64, 66]), which must comply with LEP results on  [67]; moreover, any contribution to cLFV decay modes should not exceed the present uncertainty on the total width [67], . In our study we also take into account current limits on invisible Higgs decays (relevant for ), following the approach derived in [68, 69, 70]. Likewise, negative results from laboratory searches for monochromatic lines in the spectrum of muons from decays are also taken into account [71, 72]. As mentioned in the Introduction, the new states (through the modified currents) induce potentially large contributions to cLFV observables; we evaluate the latter [29, 30, 31, 32, 33, 34, 35, 37, 40] imposing available limits on a wide variety of observables (some of them collected in Table 1). In addition to the cLFV decays and transitions, which can prove instrumental to test and disentangle these extensions of the SM, important constraints arise from rare leptonic and semileptonic decays of pseudoscalar mesons decays (including lepton universality violating, cLFV and lepton number violating modes); we include constraints from numerous , , modes (see [73, 74] for kaon decays, [75, 76] for and decay rates, and [77, 78] for -meson observations), stressing that in the framework of the SM extended by sterile neutrinos particularly severe constraints arise from the violation of lepton universality in leptonic kaon decays (parametrized by the observable [79, 66]. Finally, we also take into account the recent constraints on neutrinoless double beta decay [80]: should the sterile states be Majorana fermions, they can potentially contribute to to the effective mass  [81], which we evaluate following [82, 83].

A number of cosmological observations [84, 71, 85, 86] put severe constraints on sterile neutrinos with a mass below the GeV (in particular below 200 MeV). In our study we will in general explore regimes associated with heavier sterile states ( GeV) so that these constraints are not expected to play a relevant rôle.

### 2.2 Theoretical framework

Several mechanisms of neutrino mass generation, which in addition to accommodating neutrino data, also address in the baryon asymmetry of the Universe and/or put forward a viable dark matter candidate, call upon sterile fermions. Among such models, one encounters appealing SM extensions such as the Inverse Seesaw [44], the MSM [45], or several low-scale type I seesaw variants.

#### 2.2.1 The simple “3+1 model”

As done in previous studies of cLFV in SM extensions via sterile neutrinos, one can use as a first phenomenological approach a minimal “toy model”, consisting in the addition of a single Majorana sterile neutral fermion to the SM field content [38, 40]. This ad-hoc construction makes no assumption on the mechanism of neutrino mass generation; it thus allows to decouple the neutrino mass generation (which could possibly arise at a different, higher scale, or stem from interactions not calling upon the lighter sterile state) from the mechanism at the origin of flavour violation. In such a toy construction, the additional sterile state can also be interpreted as encoding the effects of a larger number of states possibly present in the model.

The simple toy model - which will be adopted in the present study - thus relies on the simple hypothesis that the interaction eigenstates and the physical ones are related via a unitary mixing matrix, . Other than the masses of the three light (mostly active) neutrinos, and their mixing parameters, the simple “3+1 model” can be parametrized via the heavier (mostly sterile) neutrino mass , three active-sterile mixing angles as well as three new CP violating phases (two Dirac and one Majorana). In the numerical analyses we will in general consider a normal ordering for the light neutrino spectra; in what concerns the new degrees of freedom, we will scan over the following range for the mass of the additional heavy state,

 0.5 GeV ≲m4≲106 GeV, (4)

while the active-sterile mixing angles are randomly taken to lie in the interval (as are the different CP violating phases).

#### 2.2.2 Complete theoretical frameworks for neutrino mass generation

Several mechanisms of neutrino mass generation, which in addition to accommodating neutrino data, address in addition the BAU and/or put forward a viable DM candidate, call upon sterile fermions. Inverse seesaw realisations, as well as the MSM, whose main features will be briefly summarised below, are an example of such extensions, known for their rich phenomenological implications.

The (3,3) Inverse Seesaw realisation

The Inverse seesaw mechanism [44] relies in extending the SM via right handed neutrinos and further sterile states. In the present analysis we will consider a realisation of the ISS in which three generations of RH neutrinos as well as three generations of extra singlet fermions are added to the SM, ; both and carry lepton number, . The Lagrangian describing this extension can be cast as

 LISS=LSM−Yνij¯νRi~H†Lj−MRij¯νRiXj−12μXij¯XciXj+H.c., (5)

with and generation indices. The light neutrino spectrum (containing mostly active states) is given by a modified seesaw relation

 mν≈(Yνv)2μXM2R (6)

where is the unique source of lepton number violation in the model. Small values of (which are thus natural in the sense of ’t Hooft) allow to accommodate the smallness of active neutrino masses for sizeable values of , and hence a comparatively low seesaw scale ( lying close to the TeV scale). The spectrum of the (3,3) ISS further contains three nearly degenerate pseudo-Dirac pairs; these heavier, mostly sterile states have masses close to (their degeneracy being lifted by ). The full mass matrix, , can be diagonalised as , with . In the physical charged lepton basis, the leptonic mixings are encoded in the rectangular sub-matrix () defined by the first three columns of , its upper block corresponding to the non-unitarity .

Depending on the specific realisation, and on the regimes for and , the ISS can further account for the observed BAU via leptogenesis [87], as well as provide viable DM candidates whose relic density is in agreement with present observations, and which could also accommodate possible indirect DM detection signals (if confirmed) [89, 88].

The Minimal Standard Model

The MSM minimally extends the SM via the inclusion of three RH neutrinos, aiming at simultaneously addressing the problems of neutrino mass generation, the BAU and providing a viable DM candidate [45, 90, 91, 92]. The new particle content leads to new terms in the leptonic Lagrangian:

 LνMSMmass=−Yνij¯νRi~H†Lj−12¯νRiMMijνcRj+H.c., (7)

where are generation indices, is the lepton doublet and ; denotes the Yukawa couplings, while is a Majorana mass matrix (leading to the violation of the total lepton number, ).

Other than three light (mostly active) neutrinos, the spectrum contains three heavy states (with masses ), whose masses and mixings to the lighter states are strongly constrained in the case in which the MSM is called to successfully address the BAU and the DM problems.

## 3 cLFV in-flight ℓi→ℓj conversion

In what follows we summarise the most relevant points regarding the computation of the observables associated with the in-flight cLFV conversion; due to the underlying process, in which an intense lepton beam hits a fixed target, the observable is also frequently referred to as an “on target” cLFV transition, . As mentioned in the Introduction, there are several possibilities regarding the final state of the nuclei (target) after interaction with the energetic beam: elastic scattering, in which ; quasi-elastic scattering, leading to a final state target composed of several bodies (but conserving the total number of nucleons, with no new hadronic states); inelastic processes (including excited nuclear states), and/or nuclear fragmentation with associated pion or other light hadron production (DIS regime). In the present phenomenological analysis we will focus on the case of elastic scatterings444While inelastic scattering is expected to become dominant at large enough , e.g. above for electron-proton scattering, its description is beyond the purpose of this paper.; quasi-elastic processes (as well as inelastic ones) were also recently addressed in the study of [27]).

The kinematics of the in-flight cLFV conversion requires the beam to have a minimal threshold energy (which depends on the nature of the target and on the mass of final state lepton). Denoting the intervening quadri-momenta as

 ℓi(k)+T(p)→ℓj(k′)+T(p′),%withQ2=−q2=−(k−k′)2=2MTΔEbeam, (8)

with the energy loss of the beam, and the target’s mass, one thus finds that the (threshold) beam energy is555While Eq. (9) leads to an effective lower bound to the beam energy, as previously mentioned we will not enter high-energy regimes leading to DIS phenomena.

 Ebeam>mℓj(1+mℓj2MT), (9)

in which denotes the mass of the heavier lepton in the final state (muon or tau). Moreover, a non-zero momentum transfer to the nuclear system is unavoidable. Depending on the beam’s energy, and the composition of the target, one finds minimal values for the energy transfer - although these do decrease with increasing beam energy and with the (larger) size of the nuclei, non-zero values of are always obtained (see [27] for a comprehensive discussion).

In the framework of NP models in which cLFV occurs via higher-order (loop) transitions (as is the case of R-parity conserving SUSY, seesaw realisations, etc.), the differential cross section for the cLFV conversion of Eq. (8) receives contributions from different processes, depending on the interaction(s) at the source of flavour violation: photon dipole, - and Higgs-penguins, box diagrams, among other contributions. In what follows, we proceed to discuss them.

The differential cross section for the on-target conversion of , exclusively due to photon dipole exchanges (i.e., putting to zero all other contributions), can be written as [27]

 dσi→jdQ2∣∣∣γ=πZ2α2Q4E2beamHγμνLγμνij, (10)

in which denotes the target atomic number. The detailed expression for the hadronic tensor can be found in the Appendix A, while the leptonic tensor can be decomposed as

 Lγμνij=LγijLγμν(k,q), (11)

in which encodes the cLFV (effective) couplings. Important contributions to the on-target cLFV conversion arise from the -mediated interaction. Likewise, and in the limiting case in which only -interactions are present, one can write

 dσi→jdQ2∣∣∣Z=G2F32πE2beamHZμνLZμνij, (12)

with

 LZμνij=LZijLZμν(k,q), (13)

where, and as before, the terms encode the cLFV couplings. Other contributions, such as Higgs mediated interactions (as in the case of SUSY models), box diagrams, etc., might be also present and, depending on the given model (and regime), play a relevant rôle.

The couplings can be interpreted as generic sources of flavour violation at the origin of the cLFV in-flight conversion, in the framework of SM extensions in which cLFV receives important (if not dominant) contributions from penguin loop diagrams; however, and in what follows we focus on a minimal NP model: the SM minimally extended by additional (massive) neutrinos. In such a framework, the most important contributions indeed arise from and photon mediated interactions, mediated box diagrams, and corrections to the lepton propagators, some of them schematically depicted in Fig. 1. In our analysis, we will not take into account the contributions arising from the “handbag” (box) diagrams as in the limiting (unrealistic) case of a real quark, these diagrams would correspond to the usual box contributions common to several observables (such as conversion in Nuclei, or decays). In minimal SM extensions via sterile fermions, and in the large sterile mass regime - which has been shown to be associated with sizeable contributions to the above mentioned decays - the box contributions typically lead to subdominant contributions when compared to the and penguins [37, 38, 40, 39]. Other regimes are known to be associated with important box-diagram contributions [37]. It is worth stressing that the Wilson coefficients for the contribution of boxes, photon and -penguins (cf. Fig. 1) have been evaluated in the SM extended by sterile massive fermions with non-negligible active-sterile mixings [39], with results confirming the above statement.

In this context, can be cast as

 Lγij=α3ws2w64πe2m2ℓjM4W∣∣Gγji∣∣2, (14)

with denoting the photon-lepton dipole coupling, also contributing to other cLFV transitions such as , and which is given in Appendix B; the flavour violating -couplings can be written

 LZij=α4wG2FM4W2(−1/2+sin2w)2+sin4w64∣∣FZji∣∣2, (15)

in which denotes the form factor encoding flavour violating interactions, which is also present in several other cLFV observables (see Appendix B). The full expressions for , as well as that of full leptonic and hadronic tensors are given in Appendix A.

While in low sterile mass regimes the photon penguin does dominate over the , increasing the mass of the sterile neutrinos - which corresponds to regimes typically associated with a significant enhancement of the contributions to many cLFV (in particular to the in-flight differential cross sections under study) - leads to having a -penguin contribution which increasingly dominates over the photon-ones. Although this cannot be straightforwardly inferred by comparing Eq. (10) and Eq. (12) - since the source of cLFV is encoded in the form factors and of of Eqs. (14, 15), respectively - we notice that contrary to diagrams in which a single neutrino and a run in the loop (see Fig. 1, upper-right diagrams), the -penguin further receives contributions from loops where two neutrino states and one boson are present (Fig. 1, upper-left diagram).

As is clear from the above discussion concerning the cLFV couplings, current bounds on many low-energy observables (see Table 1) will play a very constraining rôle on the maximal viable values for the in-flight conversion cross section. Particularly relevant will prove to be the bounds from decays, radiative decays, as well as conversion in Nuclei.

Before entering the study of the prospects for the cLFV on-target conversion in extensions of the SM via sterile fermions, we briefly discuss some issues regarding the nuclear interaction and the beam energy, which can be already understood from the differential cross section, . The nuclear tensors - for both photon and -mediated interactions - can be computed for either spin 0 and spin 1/2 targets. In our phenomenological study, we consider elastic interactions with individual nucleons, that is with spin 1/2 protons and/or neutrons (which corresponds to setting and in the relevant equations). The individual differential cross sections, corresponding to the purely - or -mediated exchanges, for conversion on a neutron target, are displayed on the left panel of Fig. 2 as a function of the momentum transfer, , and for two different beam energies,  GeV. These have been evaluated by simply setting by hand, in a model-blind manner, maximal values for the flavour violating terms , see Eqs. (14, 15). (Leading to the results displayed in this section, no observational bounds have been applied.) Although depending on the actual SM extension under consideration (and in the specific case of additional sterile fermions, on the particular mass regime), -mediated FV conversions often prove to dominate over the photon dipole exchanges (see [37, 38, 40]), the example seen in the left panel of Fig. 2 being typical of heavy sterile masses in the 1-10 TeV range.

Unless otherwise stated, in the following numerical discussion, we will in general consider that -penguins provide the dominant contributions to the observables under study.

On the right panel of Fig. 2 we compare the -mediated contribution for the individual nucleons (proton and neutron). In view of the very similar behaviour for both nucleons, in the following we will for simplicity assume a neutron target (unless otherwise explicitly mentioned). Likewise, and in agreement with the findings of [27], there is only a small difference, typically below , regarding the differential cross section associated with the cLFV conversion of leptons or anti-leptons (cf. Eq. (25), Appendix A); thus in our analysis we will discuss . Even though the results displayed in Fig. 2 correspond to conversion, qualitatively analogous ones have been found for an electron beam (with final state muons or taus).

A second comment concerns the dependency of the differential cross section on the beam’s energy, which was already manifest in the results of Fig. 2. Although both photon and mediated contributions explicitly scale as , the hadronic tensors (see Appendix A) both have non-trivial dependencies (also via - cf. Eq. (8)). The left panel of Fig. 3 generalises the choices of beam energy,  GeV, presented in Fig. 2; for larger values of the beam energy one enters the strong DIS regime - in the latter case, the behaviour of the differential cross section must be interpreted as only illustrative (the results here computed no longer quantitatively hold). For a fixed value of the momentum transfer (which maximises the conversion rate), the dependency of the differential cross section on the beam energy is illustrated on the right panel of Fig. 3. The latter confirms that once the beam energy is sufficiently large to reach the threshold for the in-flight conversion to occur, see Eq. (9), the rate mildly increases until rapidly saturating (in the displayed case at  GeV).

## 4 Experimental prospects

The total expected number of produced leptons for the in-flight conversion can be written as

 Nconver.(ℓi→ℓj)=Nℓi×P(ℓi→ℓj), (16)

where denotes the number of leptons () hitting the target, and the conversion probability. For the case of conversion, the total number of signal events can be directly obtained from the above equation, simply rescaling via parameters associated with the specificity of the target (thickness and density , or equivalently, the target’s mass - expressed in g/cm). For the case of final state tau leptons, their average lifetime implies that they will rapidly decay, and hence one has a further correction factor of BR(), which in the SM is approximately 17.4% [67]. Thus, the final number of expected conversions can be cast as [19]

 Nsignal(ℓi→ℓj)=Nℓi×σ(ℓi→ℓj)×Tm×Np+n[×BR(τ→μνν)], (17)

with the integrated cross section and the total number of nucleons per gramme of target - assuming for simplicity an average value of the contributions from protons and neutrons to the total cLFV conversion cross section. One thus finds

 Nsignal(ℓi→ℓj)=Nℓi×(σ(ℓi→ℓj)fb)×(Tmg cm−2)×6×10−16[×BR(τ→μνν)]. (18)

Recall that in the above two equations, the last term is only present when the final lepton is a . In order to discuss the real expected number of events, one should further take into consideration the detector’s intrinsic efficiency, , as well as the relevant contributions to the background - which we will not address in the present study.

In Table 2 we collect some operating benchmark values (surface density of the target and intensity of the beam), previously considered in former discussions of this cLFV observable.

The simple “3+1 model”

Hereafter focusing on the most minimal “3+1 model”, described in Section 2.2.1, we begin our discussion of the integrated cross section for the several cLFV in-flight conversion modes; as an illustrative case, we present the results obtained for a lepton beam energy of 4 GeV (independent of its nature, electron or muon). Prospects for different (higher) energy beams have already been briefly commented in the previous section, and the qualitative outcome holds for the present discussion. Moreover, and although having carried the numerical computation of both and photon penguin contributions, we only present the contributions of the former, which in our framework are dominant with respect to those of the latter.

The different panels of Fig. 4 display a general survey of the expected contributions to the different cross sections (arising from -mediated cLFV interactions), as a function of the mass of the heavy, mostly sterile state. The left column of Fig. 4 confirms that the cLFV cross section rapidly increases for heavy neutrino masses above the EW scale. Although one could potentially have values for the different observables as large as , current experimental bounds - in particular those arising from the violation of several cLFV bounds - exclude these regimes. In terms of expected number of converted leptons, having at least 10 conversions per year lies beyond realistic prospects for beam intensities: even for the least constrained observable, conversion, very intense muon beams on a dense target cannot account for more than 0.04 converted tau leptons per year (for one would have at best 0.02 converted s, and even lower numbers for conversion).

For final state tau leptons, the strongest cLFV constraint arises from the corresponding 3-body decays (), while for conversion the current bounds on CR(, Au) further add to the already constraining rôle of . The right hand side column of Fig. 4 summarises this discussion, displaying as a function of the flavour violation in -mediated interactions, - see Eq. (15), and Appendix A. Horizontal lines denote the cross sections that would account for a minimum of 10 conversions per year (the different line scheme corresponding to the relevant operating benchmarks of Table 2). Other than the coloured points associated with the leading cLFV constraints, grey points are associated with further exclusions arising from many other observables - as described in Section 2.1.

As extensively discussed in the literature, the interplay of distinct cLFV observables (arising from different sectors, and studied at different energies and experimental setups) is a potentially powerful probe to test flavour violating extensions of the SM. For the case of our minimal framework - extending the SM with one sterile fermion - we illustrate in Fig. 5 the potential synergies between the in-flight conversion rate and other cLFV observables, for which -penguin exchanges are known to provide important (if not dominant) contributions: BR(), BR(), and - in the case of conversion, CR(, N). As could be expected, there is a clear correlation between the in-flight and both high-intensity and high-energy observables.

Should one dispose of an unlimited number of leptons in the beam, the in-flight cLFV conversion could simultaneously probe - or even be complementary to other low - and high-energy cLFV observables. Nevertheless, the small expected number of converted leptons, for what are already optimistic beam configurations, dismisses the latter possibilities.

To finalise the discussion, we briefly comment on the prospects for this cLFV observable in well-motivated mechanisms of mass generation, such as the ISS.

(3,3) Inverse seesaw realisation

The numerical results for the ISS here displayed were obtained relying on a random scan over the neutrino mass matrix (for a detailed discussion of the numerical studies, see for example [38]); we take the following ranges for the and matrices: and , with complex entries for the lepton number violating matrix . In order to accommodate neutrino oscillation data, we use a modified Casas-Ibarra parametrisation [93] for , with complex angles for the matrix which encodes the additional degrees of freedom (these are randomly varied in the interval ), always verifying that the Yukawa couplings are perturbative, i.e. . All bounds referred to in Section 2.1 are taken into account. For the purpose of this section, we consider a NH for the light neutrino spectrum.

We illustrate the synergy between the in-flight conversion and other cLFV observables in the framework of the (3,3) ISS; the distinct panels of Fig. 6 summarise a study similar to that displayed in Fig. 5.

The summary of the ISS prospects, collected in Fig. 6 confirms what had been previously found in studies of other cLFV observables (among them 3-body decays, conversion in Nuclei, or cLFV decays): although such an ISS realisation can in principle account for sizeable values of the in-flight cLFV conversion, experimental bounds preclude the associated regimes. Fot instance, the maximal expected values for the integrated cross section does not exceed (while in the simple “3+1 toy model” one could have regions with above .

While the simple “3+1 toy model” (in which the active-sterile mixing is only constrained from experimental bounds), in the (3,3) ISS realisation the flavour violating structures (i.e., the Yukawa couplings and the LNV matrix) reflect correlations which are a consequence of necessarily accommodating oscillation data. Unlike the simple “3+1 toy model”, which allowed to independently explore different directions in flavour space (and thus, for example, evade sector constraints while enhancing flavour violation), the ISS thus offers a far more constrained scenario. We do not dismis that special textures - i.e., strongly suppressing mixings subject to the most stringent experimental bounds, while enhancing those which play a leading rôle in the observable - could account for higher values [36, 94]. Nevertheless, these are somewhat fine-tuned constructions, which we will not pursue in the present analysis.

MSM

We have also numerically explored the prospects of the MSM; despite the additional degrees of freedom - in particular three new mixing angles (other than the new Dirac and Majorana CP violating phases) - the allowed MSM parameter space [92] leads to very poor results for cLFV observables, with maximal values of the cLFV in-flight cross sections many orders of magnitude below those arising in the framework of the (3,3) ISS realisation above discussed. We notice that due to the very low scale of the new states (typically below 100 GeV), which are accompanied by not excessively large mixings, the general prospects of the MSM for cLFV are not as appealing as those of other low-scale seesaw realisations (see, for example, [38, 40, 95]), a direct consequence of the size of its intrinsic sources of flavour violation.

## 5 Concluding remarks

In the past years, charged lepton flavour violating observables have gained an increasing interest stemming from their potential to probe scenarios of New Physics, even those whose typical scales lie beyond collider reach. In view of upcoming facilities, which are expected to operate with intense lepton beams (for example those dedicated to high-intensity cLFV searches as COMET, NA64, future neutrino factories, or even a Muon Collider), in-flight lepton flavour conversions occurring when the intense beams hit a fixed target, are potentially interesting cLFV observables.

In this study we have thus revisited cLFV in-flight conversion, , and , focusing on elastic interactions with a nucleus (considering moderately energetic beams, with an energy not far from the kinematical threshold). We have studied the different contributions to the differential cross sections, and our findings concerning the derivation of the leptonic and hadronic tensors are in agreement with those of Ref. [27]. Motivated by classes of NP models in which cLFV processes occur at higher order, we have moreover focused on the dipole and -penguin contributions to the in-flight cLFV conversion.

After a general discussion of the observable, we carried a thorough phenomenological analysis in the framework of minimal SM extensions via sterile neutrinos, in which -penguin transitions do indeed dominate over the dipole contributions (and box diagrams as well). Although such minimal frameworks do offer the possibility to have sizeable values for the cross sections, , these values are precluded due to the stringent bounds arising from a number of other cLFV observables. Particularly constraining are those observables in which the -penguin contributions also play a relevant rôle - among them BR(), BR(), BR(), and CR(, Au). Once the latter bounds are taken into account, the distinct cross sections are strongly reduced - at most one expects values of , for the case of conversion (for which the associated low-energy cLFV constraints are less stringent). Even when assuming the possibility of very intense lepton beams, our study suggests that the expected number of conversions lies beyond experimental sensitivity (below ).

Other theoretical frameworks relying on extensions of the SM via several sterile fermions were found to lead to similar (or even worse) prospects: studies of the in-flight cLFV observables in complete models as the (3,3) ISS realisation or the MSM, were carried, and our findings confirmed that such frameworks indeed accounted for smaller predictions to the distinct observables than what is found in the framework of the simple “3+1 model”.

Albeit the results here obtained concern minimal SM extensions via sterile neutrinos, the strong correlation between the in-flight conversion and the cLFV observables which preclude its observability should be common to other NP constructions exhibiting similar features. This is the case of minimal (constrained) SUSY models, where there is typically a strong correlation between radiative decays and the -penguins providing the dominant contributions to 3-body decays; in this sense, our findings confirm those of [24] which pointed out that former bounds on already forbade SUSY contributions to larger than  fb.

It is also worth considering the possibility of having additional sterile states: if on the one hand this might contribute to enhance the cross sections (via a multiplicative factor, thus leading at most to a single order of magnitude enhancement), the additional states would also contribute to the other cLFV observables, so that one does not expect an overall improvement. Likewise, a study in the DIS limit should not qualitatively change the general results here derived.

Should experimental searches for the in-flight cLFV conversion observable be carried in the future, and should an event be observed, then another source of flavour violation, different from - or in addition to - those present in minimal SM extensions via sterile fermions must be necessarily present. Moreover, available (phenomenological) results would suggest that such a NP model would likely exhibit a smaller degree of correlation between different cLFV observables (as is the case of leptoquark models): for example, some transitions occurring at tree-level, while others being mediated via higher order exchanges.

Finally, and as in the case of any cLFV observable, the experimental observation of the in-flight cLFV conversion (as could happen in the near future at NA64 [23, 96]), would clearly signal the presence of New Physics, and allow selecting classes of models (other than those here discussed) which could account for it.

## Acknowledgements

We are grateful to S. Gninenko for interesting remarks and contributions. This work was done within the framework of the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreements No 690575 and No 674896. VDR acknowledges support by the Spanish MINECO through the project FPA2012-31880 (P.I. E. Alvarez Vazquez).

## Appendix A Nuclear and leptonic tensors

We describe the most relevant elements leading to the computation of both the photon- and -mediated interactions.

### a.1 Photonic interaction

The hadronic tensor relevant for the photon mediated on-target conversion, as given in Eq. (10), can be cast as

 Hγμν= −(ημν−qμqν/q2)W1+1M2T(pμ−qμp.q/q2)(pν−qνp.q/q2)W2, (19)

where, and for a spin target, one has

 W1=Q24M2T(F1+F2)2,W2=F21+Q24M2TF22, (20)

with the Dirac and Pauli form factors, which in our analysis refer to the nucleon form factors, . In agreement with [27], one can write the latter as:

 Fp(n)1(Q2)=11+Q2/4M2N[1(0)1+Q2/4M2V+Q24M2Nμp(μn)1+Q2/4M2V], Fp(n)2(Q2)=11+Q2/4M2N[μp−1(μn)1+Q2/4M2V], (21)

where denotes the nucleon mass (), and the relevant scale for the interaction, , and the total magnetic moments, respectively .

Likewise, the leptonic tensor also present in Eq. (10) can be expressed in terms of momenta as666We adopt a similar notation to that of Ref. [27], the results of which we agree with.

 Lγμν= −2[m2ℓ(m2ℓ−q2)(ημν−qμqν/q2)+4q2(kμ−qμk.q/q2)(kν−qνk.q/q2)], (22)

in which denotes the mass of the final state (heavier) lepton.

Bringing all the elements together, the final expression for the photon contribution to the differential cross section is given by

 dσi→jdQ2∣∣∣γ= 2πZ2α2E2Q4Lγij{W1(Q2+m2ℓj)(2m2ℓj−Q2)+ +W2M2T(4Q2(p.k)2+(Q2+m2ℓj)[(p.q)2−4p.qp.k+M2Tm2ℓj])}, (23)

with given in Eq. (14).

### a.2 Z-mediated interaction

Assuming the case of unpolarised lepton beams, the leptonic tensor entering in the -interaction contribution to the differential cross section (see Eq. (12)) can be written as

 LZμν= 16(kμk′ν+kνk′μ−k.k′ημν+iεμνρσkρqσ). (24)

The hadronic tensor can be in general cast in terms of six dimensionless structure functions as

 HZμν= −ημνW1+pμpνM2TW2±iεμνρσpρqσ2M2TW3+qμqνM2TW4+pμqν+pνqμ2M2TW5+ipμqν−pνqμ2M2TW6, (25)

where in the above equation the corresponds to having a lepton (or antilepton) conversion. The different expressions for the structure functions can be found in Ref. [27], with which we agree after an independent derivation.

As above, the contraction of both leptonic and hadronic tensors leads to the following differential cross section for -mediated contribution, which we have used throughout the analysis,

 dσi→jdQ2∣∣∣Z= G2F2πE2LZij{[(Q2+m2ℓj)(W1−12W2)+p.kM2N(2p.k−Q2)W2+ +12(Q2+m2ℓj)m2ℓjM2NW4−(p.k)m2ℓjM2NW5⎤⎦±Q24M2N(4p.k−Q2−m2ℓj)W3⎫⎬⎭, (26)

with has been given in Eq. (15).

## Appendix B cLFV form factors

The relevant form factors for the computation of the diagrams of Fig. 1 are given by [35, 31, 29, 30]:

 Gℓmγ = 3+nS∑j=1UmjU∗ℓjGγ(xj),Fℓmγ=3+nS∑j=1UmjU∗ℓjFγ(xj), FℓmZ = 3+nS∑j,k=1UmjU∗ℓk(δjkFZ(xj)+CjkGZ(xj,xk)+C∗jkHZ(xj,xk)), (27)

where carries the neutrino mass dependency and has been defined in Eq. (2).

The loop functions entering the previous form factors are defined as [35, 31, 29, 30]:

 FZ(x)=−5x2(1−x)−5x22(1−x)2lnx, GZ(x,y)=−12(x−y)[x2(1−y)1−xlnx−y2(1−x)1−ylny], HZ(x,y)=√xy4(x−y)[x2−4x1−xlnx−y2−4y1−ylny], Fγ(x)=x(7x2−x−12)12(1−x)3−x2(x2−10x+12)6(1−x)4lnx, Gγ(x)=−x(2x2+5x−1)4(1−x)3−3x32(1−x)4lnx. (28)

The contributions to different cLFV observables such as radiative and 3-body decays, conversion in Nuclei, or FV decays, which have been evaluated and analysed in the present work (including the relevant loop functions [35, 31, 29, 30]), have been discussed in previous studies (see, for example, [35, 31, 38, 40]), and we will not include them here.

## References

• [1] A. M. Baldini et al. [MEG Collaboration], Eur. Phys. J. C 76 (2016) no.8, 434 [arXiv:1605.05081 [hep-ex]].
• [2] A. M. Baldini, F. Cei, C. Cerri, S. Dussoni, L. Galli, M. Grassi, D. Nicolo and F. Raffaelli et al., arXiv:1301.7225 [physics.ins-det].
• [3] B. Aubert et al. [BaBar Collaboration], Phys. Rev. Lett. 104 (2010) 021802
• [4] T. Aushev, W. Bartel, A. Bondar, J. Brodzicka, T. E. Browder, P. Chang, Y. Chao and K. F. Chen et al., “Physics at Super B Factory,” arXiv:1002.5012 [hep-ex].
• [5] U. Bellgardt et al. [SINDRUM Collaboration], Nucl. Phys. B 299 (1988) 1.
• [6] A. Blondel, A. Bravar, M. Pohl, S. Bachmann, N. Berger, M. Kiehn, A. Schoning and D. Wiedner et al., “Research Proposal for an Experiment to Search for the Decay ,” arXiv:1301.6113 [physics.ins-det].
• [7] K. Hayasaka, K. Inami, Y. Miyazaki, K. Arinstein, V. Aulchenko, T. Aushev, A. M. Bakich and A. Bay et al., Phys. Lett. B 687 (2010) 139 [arXiv:1001.3221 [hep-ex]].
• [8] W. H. Bertl et al. [SINDRUM II Collaboration], Eur. Phys. J. C 47 (2006) 337.
• [9] M. Aoki [DeeMe Collaboration], AIP Conf. Proc. 1441 (2012) 599.
• [10] R. M. Carey et al. [Mu2e Collaboration], FERMILAB-PROPOSAL-0973.
• [11] Y. G. Cui et al. [COMET Collaboration], KEK-2009-10.
• [12] Y. Kuno [COMET Collaboration], PTEP 2013 (2013) 022C01.
• [13] M. Koike, Y. Kuno, J. Sato and M. Yamanaka, Phys. Rev. Lett. 105 (2010) 121601 [arXiv:1003.1578 [hep-ph]].
• [14] Y. Uesaka, Y. Kuno, J. Sato, T. Sato and M. Yamanaka, Phys. Rev. D 93 (2016) no.7, 076006 [arXiv:1603.01522 [hep-ph]].
• [15] L. S. Littenberg and R. Shrock, Phys. Lett. B 491 (2000) 285 [hep-ph/0005285].
• [16] T. Geib, A. Merle and K. Zuber, Phys. Lett. B 764 (2017) 157 [arXiv:1609.09088 [hep-ph]].
• [17] J. M. Berryman, A. de Gouvea, K. J. Kelly and A. Kobach, “On Lepton-Number-Violating Searches for Muon to Positron Conversion,” arXiv:1611.00032 [hep-ph].
• [18] T. Geib and A. Merle, arXiv:1612.00452 [hep-ph].
• [19] S. N. Gninenko, M. M. Kirsanov, N. V. Krasnikov and V. A. Matveev, Mod. Phys. Lett. A 17 (2002) 1407 doi:10.1142/S0217732302007855 [hep-ph/0106302].
• [20] M. Sher and I. Turan, Phys. Rev. D 69 (2004) 017302 [hep-ph/0309183].
• [21] M. Gonderinger and M. J. Ramsey-Musolf, JHEP 1011 (2010) 045 Erratum: [JHEP 1205 (2012) 047] [arXiv:1006.5063 [hep-ph]].
• [22] S. Kanemura, Y. Kuno, M. Kuze and T. Ota, Phys. Lett. B 607 (2005) 165 [hep-ph/0410044].
• [23] S. N. Gninenko, Phys. Rev. D 89 (2014) no.7, 075008 [arXiv:1308.6521 [hep-ph]]; S. Andreas et al., “Proposal for an Experiment to Search for Light Dark Matter at the SPS,” arXiv:1312.3309 [hep-ex]. See also https://na64.web.cern.ch/
• [24] T. Blazek and S. F. King, “Electron to muon conversion in electron-nucleus scattering as a probe of supersymmetry,” hep-ph/0408157.
• [25] K. P. O. Diener, Nucl. Phys. B 697 (2004) 387 [hep-ph/0403251].
• [26] A. Bolanos, A. Fernandez, A. Moyotl and G. Tavares-Velasco, Phys. Rev. D 87 (2013) no.1, 016004 [arXiv:1212.0904 [hep-ph]].
• [27] W. Liao and X. H. Wu, Phys. Rev. D 93 (2016) no.1, 016011 [arXiv:1512.01951 [hep-ph]].
• [28] R. H. Bernstein and P. S. Cooper, Phys. Rept. 532 (2013) 27 [arXiv:1307.5787 [hep-ex]].
• [29] E. Ma and A. Pramudita, Phys. Rev. D 22 (1980) 214.
• [30] M. Gronau, C. N. Leung and J. L. Rosner, Phys. Rev. D 29 (1984) 2539.
• [31] A. Ilakovac and A. Pilaftsis, Nucl. Phys. B 437 (1995) 491 [hep-ph/9403398].
• [32] F. Deppisch and J. W. F. Valle, Phys. Rev. D 72 (2005) 036001 [hep-ph/0406040].
• [33] F. Deppisch, T. S. Kosmas and J. W. F. Valle, Nucl. Phys. B 752 (2006) 80 [hep-ph/0512360].
• [34] D. N. Dinh, A. Ibarra, E. Molinaro and S. T. Petcov, Decays and TeV Scale See-Saw Scenarios of Neutrino Mass Generation,” JHEP 1208 (2012) 125 [Erratum-ibid. 1309 (2013) 023] [arXiv:1205.4671 [hep-ph]].
• [35] R. Alonso, M. Dhen, M. B. Gavela and T. Hambye, JHEP 1301 (2013) 118 [arXiv:1209.2679 [hep-ph]].
• [36] E. Arganda, M. J. Herrero, X. Marcano and C. Weiland, Phys. Rev. D 91 (2015) no.1, 015001 [arXiv:1405.4300 [hep-ph]].
• [37] A. Abada, M. E. Krauss, W. Porod, F. Staub, A. Vicente and C. Weiland, JHEP 1411 (2014) 048 [arXiv:1408.0138 [hep-ph]].
• [38] A. Abada, V. De Romeri, S. Monteil, J. Orloff and A. M. Teixeira, JHEP 1504 (2015) 051 [arXiv:1412.6322 [hep-ph]].
• [39] A. Abada, D. Becirevic, M. Lucente and O. Sumensari, Phys. Rev. D 91 (2015) no.11, 113013 [arXiv:1503.04159 [hep-ph]].
• [40] A. Abada, V. De Romeri and A. M. Teixeira, JHEP 1602 (2016) 083 [arXiv:1510.06657 [hep-ph]].
• [41] V. De Romeri, M. J. Herrero, X. Marcano and F. Scarcella, “Lepton flavor violating Z decays: A promising window to low scale seesaw neutrinos,” arXiv:1607.05257 [hep-ph].
• [42] E. K. Akhmedov, M. Lindner, E. Schnapka and J. W. F. Valle, Phys. Rev. D 53 (1996) 2752 [hep-ph/9509255].
• [43] M. Malinsky, J. C. Romao and J. W. F. Valle, Phys. Rev. Lett. 95 (2005) 161801 [hep-ph/0506296].
• [44] R. N. Mohapatra and J. W. F. Valle, Phys. Rev. D 34 (1986) 1642.
• [45] T. Asaka, S. Blanchet and M. Shaposhnikov, Phys. Lett. B 631 (2005) 151 [hep-ph/0503065].
• [46] J. Schechter and J. W. F. Valle, Phys. Rev. D 22 (1980) 2227.
• [47] M. S. Chanowitz, M. A. Furman and I. Hinchliffe, Nucl. Phys. B 153 (1979) 402.
• [48] L. Durand, J. M. Johnson and J. L. Lopez, Phys. Rev. Lett. 64 (1990) 1215.
• [49] J. G. Korner, A. Pilaftsis and K. Schilcher, Phys. Lett. B 300 (1993) 381 [hep-ph/9301290].
• [50] J. Bernabeu, J. G. Korner, A. Pilaftsis and K. Schilcher, Phys. Rev. Lett. 71 (1993) 2695 [hep-ph/9307295].
• [51] S. Fajfer and A. Ilakovac, Phys. Rev. D 57 (1998) 4219.
• [52] A. Ilakovac, Phys. Rev. D 62 (2000) 036010 [hep-ph/9910213].
• [53] D. V. Forero, M. Tortola and J. W. F. Valle, Phys. Rev. D 86 (2012) 073012 [arXiv:1205.4018 [hep-ph]].
• [54] G. L. Fogli, E. Lisi, A. Marrone, D. Montanino, A. Palazzo and A. M. Rotunno, Phys. Rev. D 86 (2012) 013012 [arXiv:1205.5254 [hep-ph]].
• [55] M. C. Gonzalez-Garcia, M. Maltoni, J. Salvado and T. Schwetz, JHEP 1212 (2012) 123 [arXiv:1209.3023 [hep-ph]].
• [56] D. V. Forero, M. Tortola and J. W. F. Valle, Phys. Rev. D 90 (2014) 093006 [arXiv:1405.7540 [hep-ph]].
• [58] M. C. Gonzalez-Garcia, M. Maltoni and T. Schwetz, JHEP 1411 (2014) 052 [arXiv:1409.5439 [hep-ph]].
• [59] I. Esteban, M. C. Gonzalez-Garcia, M. Maltoni, I. Martinez-Soler and T. Schwetz, “Updated fit to three neutrino mixing: exploring the accelerator-reactor complementarity,” arXiv:1611.01514 [hep-ph].
• [60] S. Antusch, J. P. Baumann and E. Fernandez-Martinez, Nucl. Phys. B 810 (2009) 369 [arXiv:0807.1003 [hep-ph]].
• [61] S. Antusch and O. Fischer, JHEP 1410 (2014) 94 [arXiv:1407.6607 [hep-ph]].
• [62] M. Blennow, P. Coloma, E. Fernandez-Martinez, J. Hernandez-Garcia and J. Lopez-Pavon, “Non-Unitarity, sterile neutrinos, and Non-Standard neutrino Interactions,” arXiv:1609.08637 [hep-ph].
• [63] E. Akhmedov, A. Kartavtsev, M. Lindner, L. Michaels and J. Smirnov, JHEP 1305 (2013) 081 [arXiv:1302.1872 [hep-ph]].
• [64] E. Fernandez-Martinez, J. Hernandez-Garcia, J. Lopez-Pavon and M. Lucente, JHEP 1510 (2015) 130 [arXiv:1508.03051 [hep-ph]].
• [65] L. Basso, O. Fischer and J. J. van der Bij, Europhys. Lett. 105 (2014) 11001 [arXiv:1310.2057 [hep-ph]].
• [66] A. Abada, A. M. Teixeira, A. Vicente and C. Weiland, JHEP 1402 (2014) 091 [arXiv:1311.2830 [hep-ph]].
• [67] C. Patrignani et al. [Particle Data Group Collaboration], Chin. Phys. C 40 (2016) no.10, 100001.
• [68] P. S. Bhupal Dev, R. Franceschini and R. N. Mohapatra, Phys. Rev. D 86 (2012) 093010 [arXiv:1207.2756 [hep-ph]].
• [69] C. G. Cely, A. Ibarra, E. Molinaro and S. T. Petcov, Phys. Lett. B 718 (2013) 957 [arXiv:1208.3654 [hep-ph]].
• [70] P. Bandyopadhyay, E. J. Chun, H. Okada and J. -C. Park, JHEP 1301 (2013) 079 [arXiv:1209.4803 [hep-ph]].
• [71] A. Kusenko, Phys. Rept. 481 (2009) 1 [arXiv:0906.2968 [hep-ph]].
• [72] A. Atre, T. Han, S. Pascoli and B. Zhang, JHEP 0905 (2009) 030 [arXiv:0901.3589 [hep-ph]].
• [73] E. Goudzovski [NA48/2 and NA62 Collaborations], PoS EPS -HEP2011 (2011) 181 [arXiv:1111.2818 [hep-ex]].
• [74] C. Lazzeroni et al. [NA62 Collaboration], Phys. Lett. B 719 (2013) 326 [arXiv:1212.4012 [hep-ex]].
• [75] P. Naik et al. [CLEO Collaboration], Phys. Rev. D 80 (2009) 112004 [arXiv:0910.3602 [hep-ex]].
• [76] H. -B. Li, “Proceedings, 4th International Workshop on Charm Physics (Charm 2010) : Beijing, China, October 21-24, 2010,” Int. J. Mod. Phys. Conf. Ser. 02 (2011).
• [77] B. Aubert et al. [BaBar Collaboration], Phys. Rev. D 77 (2008) 011107.
• [78] I. Adachi et al. [Belle Collaboration], Phys. Rev. Lett. 110 (2013) 131801 [arXiv:1208.4678 [hep-ex]].
• [79] A. Abada, D. Das, A. M. Teixeira, A. Vicente and C. Weiland, JHEP 1302 (2013) 048 [arXiv:1211.3052 [hep-ph]].
• [80] J. B. Albert et al. [EXO-200 Collaboration], Nature 510 (2014) 229 [arXiv:1402.6956 [nucl-ex]].
• [81] P. Benes, A. Faessler, F. Simkovic and S. Kovalenko, Phys. Rev. D 71 (2005) 077901 [hep-ph/0501295].
• [82] M. Blennow, E. Fernandez-Martinez, J. Lopez-Pavon and J. Menendez, JHEP 1007 (2010) 096 [arXiv:1005.3240 [hep-ph]].
• [83] A. Abada, V. De Romeri and A. M. Teixeira, JHEP 1409 (2014) 074 [arXiv:1406.6978 [hep-ph]].
• [84] A. Y. Smirnov and R. Zukanovich Funchal, Phys. Rev. D 74 (2006) 013001 [hep-ph/0603009].
• [85] P. Hernandez, M. Kekic and J. Lopez-Pavon, Phys. Rev. D 90 (2014) no.6, 065033 [arXiv:1406.2961 [hep-ph]].
• [86] A. C. Vincent, E. F. Martinez, P. Hernandez, M. Lattanzi and O. Mena, JCAP 1504 (2015) no.04, 006 [arXiv:1408.1956 [astro-ph.CO]].
• [87] A. Abada, G. Arcadi, V. Domcke and M. Lucente, “Lepton number violation as a key to low-scale leptogenesis,” arXiv:1507.06215 [hep-ph].
• [88] A. Abada, G. Arcadi and M. Lucente, JCAP 1410 (2014) 001 [arXiv:1406.6556 [hep-ph]].
• [89] A. Abada and M. Lucente, Nucl. Phys. B 885 (2014) 651 [arXiv:1401.1507 [hep-ph]].
• [90] T. Asaka and M. Shaposhnikov, Phys. Lett. B 620 (2005) 17 [hep-ph/0505013].
• [91] M. Shaposhnikov, JHEP 0808 (2008) 008 [arXiv:0804.4542 [hep-ph]].
• [92] L. Canetti, M. Drewes, T. Frossard and M. Shaposhnikov, Phys. Rev. D 87 (2013) 093006 [arXiv:1208.4607 [hep-ph]].
• [93] J. A. Casas and A. Ibarra, Nucl. Phys. B 618 (2001) 171 [hep-ph/0103065].
• [94] E. Arganda, M. J. Herrero, X. Marcano and C. Weiland, Phys. Lett. B 752 (2016) 46 doi:10.1016/j.physletb.2015.11.013 [arXiv:1508.05074 [hep-ph]].
• [95] L. Canetti and M. Shaposhnikov, Hyperfine Interact. 214 (2013) no.1-3, 5.
• [96] S. Gninenko, private communication.
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters