Introduction to Charged Lepton Flavor Violation

: Neutrino masses are evidence of lepton ﬂavor violation, but no violation in the interactions among the charged leptons has been observed yet. Many models of Physics Beyond the Standard Model (BSM) predict Charged Lepton Flavor Violation (CLFV) in a wide spectrum of processes with rates in reach of upcoming experiments. The experimental searches that provide the current best limits on the CLFV searches are reviewed, with a particular emphasis on the muon-based experiments that give the most stringent constraints on the BSM parameter space. The next generation of muon-based experiments (MEG-II, Mu2e, COMET, Mu3e) aims to reach improvements by many orders of magnitude with respect to the current best limits, thanks to several technological advancements. We review popular heavy BSM theories, and we present the calculations of the predicted CLFV branching ratios, focusing on the more sensitive µ → e sector.


Introduction
In the Standard Model (SM) defined with massless left-handed neutrinos, Lepton Flavor (LF) is a conserved quantity.The observation of neutrino oscillations provided clear evidence of non-zero neutrino masses and mixing angles, demonstrating that lepton flavour is not a symmetry of nature.Charged Lepton Flavour Violation (CLFV), defined as a shortrange interaction among the charged leptons, is therefore expected to occur but it is yet to be observed.If neutrinos get their masses through renormalizable Yukawa interactions with the Higgs, the predicted rates for CLFV are typically GIM suppressed G 2 F m 4 ν ∼10 −50 and are practically unobservable.A detection of CLFV would thus be a clear signature of new physics that could shed light on the origin of neutrino masses.Additionally, lepton flavour is an accidental symmetry of the SM that is respected by the most general Lagrangian with gauge invariant renormalizable interactions.Thus, independently from neutrino masses, minimal departures from the SM can easily introduce extra sources of lepton flavour violation and lead to sizeable CLFV rates.
For these reasons, experimental searches of CLFV attract great interest and are a valuable tool in identifying viable Beyond Standard Model (BSM) scenarios.CLFV searches can pinpoint theories at energy scales currently not directly accessible by the collider facilities.Null results from the current experiments significantly constrain the parameter space of new physics models, and the improvements in sensitivity by several orders of magnitude, especially in the µ → e sector, will further probe BSM physics.
In this article, we present an overview of the theoretical and experimental status of CLFV.Excellent complementary reviews on the subject already exist in the literature [1][2][3][4][5].In Section 2, we discuss several SM extensions that could be potentially probed in the upcoming CLFV experimental searches.The focus is on heavy physics models.We briefly review the flavour structure of the SM, and we discuss the LFV phenomenology of models that generate neutrino masses at tree and loop level.In the subsequent sections, we present the CLFV signature of various BSM scenarios, such as Two Higgs Doublets Model and supersymmetric SM.Finally, LFV is studied in the context of effective field theory where model-independent results are retained.
In Section 3, the state of the art and the upcoming experiments looking for CLFV processes are discussed.A particular emphasis is given to those looking for rare muon CLFV decays.Several facilities around the world (Fermilab, PSI and J-PARC) already started building or commissioning new generation experiments with improved sensitivity on the muon CLFV searches (up to four orders of magnitude).This is possible thanks to improvements in the acceleration techniques, necessary to deliver beam with unprecedented intensity ∼10 10 µ/s, and novel detector technologies.The same section also provides an overview of the current best limits achieved on the tau CLFV branching ratios set by general-purpose experiments at e + e − and pp colliders.Also on the tau front, an improved sensitivity on several searches is expected thanks to the unprecedented luminosity of the Large Hadron Collider at CERN and the SuperKEKB collider at KEK laboratory.

CLFV in Standard Model Extensions
The gauge interactions in the Standard Model (SM) are universal among different families, and the Yukawa interactions of leptons and quarks with the Higgs doublet are the only flavor defining couplings: where ℓ i = (ν iL e iL ) T , q i = (u iL d iL ) T are SU(2) left-handed doublets, while u m , d m , e j label right-handed up, down quark and lepton, respectively.H is the Higgs doublet H = (H + H 0 ) T , and H is defined as H ≡ ϵH † , having introduced the SU(2) antisymmetric tensor ϵ 12 = −ϵ 21 = 1.We define the electric charge as Q = T 3 + Y, where Y is the hypercharge and T 3 is the diagonal generator of SU (2).The labels i, j, n, m denote the generation indices, and the Yukawa matrices Y e , Y u , Y d are 3 × 3 complex matrices in flavor space.Any complex matrix can be diagonalized with a bi-unitary transformation where Ŷf are diagonal matrices with non-negative entries.The quark up and down Yukawa matrices cannot be simultaneously diagonalized in the unbroken theory, given that V q d ̸ = V q u .In a basis where the up Yukawa matrix is diagonal, the down Yukawa can be cast in the following form: where V is the Cabibbo-Kobayashi-Masukawa (CKM) unitary matrix, which parameterizes flavour violation in the quark sector.When the Electroweak gauge symmetry is spontaneously broken by the Vacuum Expectation Value (VEV) of the Higgs doublet ⟨H⟩ = 0 v T , with v ≃ 174 GeV, the SM fermions acquire masses through the Yukawa interactions so that the mass eigenstate basis d ′ L for left-handed down quark is related to the gauge interaction basis by a CKM rotation Only the charged currents are affected by a unitary rotation on d quarks, and the flavour-changing interactions with the W bosons are governed by the CKM matrix elements where g = e(sin θ W ) −1 is the SU(2) gauge coupling.In the lepton sector, the Yukawa couplings Y e are the only basis-choosing interaction, and lepton flavour U(1) L e ⊗ U(1) L µ ⊗ U(1) L τ is conserved because mass and gauge eigenstate basis coincide.(The symmetry is anomalous and is not exact at the non-perturbative level.However, B/3 − L e,µ,τ is anomaly free, hence non-baryonic processes like µ → eγ are strictly forbidden in the SM).This holds as long as neutrinos are massless and do not provide additional flavour defining interactions.Since the so-called "solar neutrino problem" of the 1960s, a deficit [6][7][8][9][10][11] in the number of electron neutrinos compared to the prediction of the standard solar model [12][13][14], neutrino oscillations have been confirmed by many observations [15][16][17][18][19] that firmly established non-zero and non-degenerate neutrino masses.Weak eigenstates are superpositions of mass eigenstates with i = e, µ, τ (7) where U is known as Pontecorvo-Maki-Nakagawa-Sakata (PMNS) [20,21] matrix, which is the lepton analogue of CKM.It is a unitary matrix, and it is parametrized by nine real parameters, three angles and six phases.Not all phases are physical, and some can be absorbed in field phase redefinition.The canonical form of PMNS is the following: where we have defined s ij = sin θ ij , c ij = cos θ ij and the matrix P is the identity if neutrinos have Dirac masses, while it contains two extra phases for Majorana neutrinos P = diag 1 e iα 12 e iα 31 .This is because, with self-conjugate left-handed neutrinos, there are less relative field re-definitions that can absorb the matrix phases.Majorana phases are difficult to observe because they contribute as extra sources of CP violation in processes that depend linearly on the masses, hence not in oscillations where the dependence is quadratic.
A recent global fit [22] on neutrino oscillation data gives the following parameters for mass-squared differences and mixing angles: where the 3σ range are displayed and normal ordering m 1 < m 2 < m 3 is assumed.The sign of ∆m 2 21 is determined by solar neutrino observation because relevant matter effects depend on it.On the other hand, atmospheric neutrino data measure only the absolute value of mass squared difference ∆m 2  31 , and they can be consistent with inverted ordering m 3 < m 1 < m 2 , which leads to marginally different best-fit values for the PMNS parameters.

CLFV in Models That Generate Neutrino Mass at Tree Level
Assuming the presence of three right-handed neutrinos ν Ri which are singlets of the SM gauge group, gauge invariance allows for Yukawa couplings between the lepton and the Higgs doublet that generate Dirac masses for neutrinos when electroweak symmetry is spontaneously broken To obtain neutrino masses that are compatible with cosmological constraints m ν ≲ 0.1 eV [23], the neutrino Yukawa couplings must be Y ν ≲ O 10 −12 .Although small Yukawas are technically natural (with technical naturalness we refer to naturalness as defined by t'Hooft.Small values for a parameter c are defined as technically natural when a symmetry is restored in taking the limit c → 0), Dirac masses require a strong hierarchy between the charged and neutral lepton Yukawa sector.
Analogously to CKM, the PMNS matrix is the result of the misalignment between charged lepton and neutrino mass basis, as the neutrino and charged lepton Yukawas cannot be simultaneously diagonalized respecting the electroweak gauge symmetry.Flavor violation is hence parameterized by the presence of the PMNS matrix in the charged lepton currents Charged lepton flavour violation is mediated by neutrino mixing.In Figure 1, we show a representative diagram for µ → eγ decay, illustrating that this is a consequence of the flavour-changing interactions in Equation (11).The amplitude of this process can be generically cast in the following form (see, for instance, Chapter 6, Section 6.2 of [24])): where P R,L = (1 ± γ 5 )/2 are the right-handed and left-handed chiral projectors, while A, B, C are complex numbers.As a consequence of QED gauge invariance, the amplitude satisfies the Ward identity q α M α = 0. On-shell spinors obey the equation of motion ( / p − m)u(p) = 0, and the Ward identity requires which, for on-shell photons q 2 = 0, has the unique solution C R = C L = 0.The only relevant piece is a dipole transition which is chirality-flipping and, thus, proportional to the muon mass (neglecting the electron mass).Equation ( 14) yields the following decay rate [1]: In the diagram of Figure 1, the outgoing electrons are left-handed and only A R is non-zero.The amplitude is proportional to the internal neutrino propagator, which can be expanded for small neutrino masses as We see that the leading term vanishes due to PMNS unitarity, and the amplitude is GIM suppressed by the square of neutrino masses.Indeed, the process is analogous to a flavour-changing neutral current in the quark sector, which features a similar GIM suppression by CKM unitarity [25].The calculation is done in [26] in the R ξ gauge, where additional diagrams replacing W with the charged Goldstones must be included.All diagrams are finite and in the unitary gauge limit ξ → ∞; only the diagram of Figure 1 is non-zero.Dividing by the rate Γ(µ → eν ν) = G 2 F m 5 µ /192π 3 of the dominant LF conserving three-body decay, the resulting branching ratios for µ → eγ is [27][28][29][30][31] Rewriting the sum as and substituting the best-fit values of the mass differences and mixing parameters in Equation ( 9), the predicted branching ratios for the LFV µ → eγ is Br(µ → eγ) = 10 −54 − 10 −55 , which lie beyond any foreseeable experimental reach.In models with Dirac neutrino masses, rates of other LFV processes are similarly GIM suppressed and, thus, too small to be observable.The right-handed neutrinos are sterile, i.e., neutral under the SM gauge group, so SM gauge invariance allows us to add a lepton number violating Majorana mass term having defined ν c = C νT , where C is the Dirac charge conjugation matrix [32].Majorana mass matrices are symmetric because fermion fields are anti-commuting and the charge conjugation matrix C is antisymmetric.Upon electroweak symmetry breaking, the mass Lagrangian can be cast in the following form (suppressing generation indices) where N = ν c L ν R T and M D = vY ν .If we assume that the Majorana masses M R are much larger than the Dirac masses (symbolically M R ≫ M D ), the matrix can be put in block diagonal form that disentangles light left-handed neutrinos and heavy righthanded ones [33] where W is a unitary matrix.At leading M D M −1 R order, the mass matrices are Assuming ∼ O(1) Yukawas, light neutrinos masses can be explained by Majorana masses close to the grand unification scale M R ∼10 15 GeV.This is the celebrated seesaw mechanism [34], specifically known as type I when the SM is extended with singlet righthanded fermions.
A unitary U * diagonalizes the symmetric Majorana matrix M ν with a congruence transformation U * T M ν U * = Mν = diag m 1 m 2 m 3 , but U is not the matrix that appears in the charged currents.Defining U ≡ U ⊗ 1 heavy as acting on the light neutrinos subspace, gauge interaction and mass basis are related by where the matrix W can be expanded at second order as [33] The left-handed weak eigenstates are related to the light mass eigenstates through a non-unitary matrix Lacking unitarity, the GIM suppression no longer operates substituting the U ′ matrix in Equation ( 16), and the rate of µ → eγ becomes [26,35] where ) is a loop function that can be expanded for x i ≪ 1 as F(x i ) ≃ 10/3 − x i .CLFV processes can thus constrain departures from the unitarity of the PMNS matrix [36].Substituting the typical value y 2 ν ≃ m ν M R /v 2 , GIM suppression is replaced by the ratio m ν /M N , which for GUT scale sterile neutrinos predict rates that are nonetheless well below future experimental sensitivity.Seesaw models can predict sizeable CLFV rates if the Majorana right-handed masses are closer to the electroweak scale.In the non-supersymmetric seesaw, this is also desirable to avoid large correction to the Higgs mass [37].However, in a generic setup with TeV scale M R and unsuppressed CLFV rates, fine-tuned cancellations are required in Equation (22) to explain neutrino masses.Fine-tuning is, of course, avoided if a symmetry principle forces the neutrino mass to be small despite having large Yukawa couplings.Observe that neutrino masses are a lepton number violating effect ∼ Y ν M −1 R Y T ν , while the non-unitary matrix that governs CLFV rates is lepton number conserving It is possible to suppress neutrino masses by invoking a small breaking of lepton number conservation while keeping the masses of the sterile neutrinos sufficiently close to the electroweak scale and with no need for small Yukawa couplings.This is, for example, achieved in the inverse seesaw [38,39] and by considering a quasi-degenerate pair of sterile neutrinos [40][41][42].
The seesaw formula can be understood as the result of integrating out the heavy neutrinos.The relevant s and t channel diagrams are shown in Figure 2 and match onto the dimension five Weinberg operator with the coefficient When the Higgs doublet gets a VEV, neutrinos acquire Majorana masses through the Weinberg operator , that agrees with Equation (22).Moreover, the following dimension six operator is generated with a coefficient , which corrects the light neutrinos kinetic terms.The redefinition needed to canonically normalize the fields introduces a non-unitary matrix in the charged currents [43] which again agrees with Equation (25).The advantage of an effective field theory description is that different seesaw scenarios can be described at low energy in a common framework.In Figure 3, we show how extending the Standard Model with particles transforming in different representations of the SM gauge group can generate Majorana neutrino masses via the Weinberg operator.Recent effective field theory analysis of type I and type II seesaw models include the complete one-loop matching onto effective operators [44,45], which are useful resources to study the low-energy CLFV signatures.For more complete reviews on the CLFV phenomenology of seesaw models, we refer the reader to [43,46].

CLFV in Models That Generate Neutrino Masses at Loop Level
New physics not too far from the electroweak scale can account for small neutrino masses if they are generated radiatively via loop diagrams.As a specific example, the so-called scotogenic model [47] adds an additional scalar doublet η (with hypercharge Y η = 1/2) to the Standard Model, together with three generations of sterile neutrinos N. The new particles are assumed to be odd under a discrete Z 2 symmetry, which forbids Yukawa couplings with the SM Higgs between lepton doublets and the sterile neutrinos, as well as constraining the possible interactions in the scalar potential.The Z 2 is also responsible for keeping stable the lightest new particle, which, if neutral, provides a potential dark matter candidate.Omitting the kinetic terms, the scotogenic Lagrangian reads (it is always possible to diagonalize the symmetric Majorana mass matrix of the sterile neutrinos with no loss of generality) where the scalar potential is To preserve the Z 2 symmetry when the electroweak symmetry is spontaneously broken, the potential parameters must be such that η field does not acquire a VEV.We also assume that all parameters in the potential are real and CP is conserved.With this assumption, the real and imaginary parts of the uncharged component η 0 = (η R + iη I )/ √ 2 do not mix.The mass splitting between the two neutral scalars is proportional to λ 5 v 2 , consequently η R,I are approximately degenerate in the limit λ 5 ≪ 1.Note that lepton number is conserved if λ 5 is zero, so small values are technically natural.
The Z 2 symmetry prevents the appearance of tree-level Majorana masses for the lefthanded neutrinos, but it can generate them at the one-loop level through the λ 5 mixing of η with the SM Higgs doublet, as shown in the diagram of Figure 4.The resulting neutrino mass matrix is calculable, and, for λ 5 ≪ 1 (m 0 = m η R ∼ m η I ), can be approximated as [47] [ With respect to the traditional seesaw scenario, the extra suppression ∼ λ 5 /(16π 2 ) can predict small vales of m ν with TeV scale sterile neutrinos and unsuppressed Yukawa couplings.The CLFV signature of the scotogenic model have been studied with particular attention to l i → l j γ processes [48][49][50], while the phenomenology of l i → l j l k l m and µ → e conversion into nuclei has also been discussed [51].In Figure 5, we show a selection of diagrams giving contributions to CLFV processes at the one-loop level.
CLFV processes in the scotogenic model.(a) Diagrams contributing to the l i → l j γ rate.(b) Box diagrams contributing to the l i → l m l n l j rate.(c) Penguin diagrams contributing to the l i → l j l k l k rate and µ → e conversion rate ( f can be a quark or a lepton).
Part of the parameter space of the scotogenic model is excluded by the current experimental LFV searches, while the viable region can give branching ratios within upcoming experimental sensitivities and will be probed in the near future.It is often the case that l i → l j γ is the most constraining LFV channel because the dipole (Figure 5a) contribution to the photon penguin (Figure 5c) can dominate the amplitude of l i → l j f f , leading to the following relation [52] Br(l i → l j γ) (33) However, the box contribution (Figure 5b) can be larger than the photon penguin diagram for mass of the lightest neutrino close to the cosmological upper limit 0.1 eV (Figure 6) so that upcoming µ → 3e searches can constrain the model orthogonally to the MEG bound on µ → eγ.In addition, the penguin diagram of Figure 5c mediates LFV interactions with quarks, contributing to the rate of µ → e conversion in nuclei (we briefly review the µ → e conversion rate calculation in Appendix A).When the dipole dominates the penguin amplitude, µ → e conversion experimental reach is not competitive with µ → eγ searches (although with the future branching ratios sensitivity Br(µA → eA) ∼ 10 −16 , µ → e conversion might be able to probe smaller dipole coefficients than MEG II with Br(µ → eγ) ∼ 6 × 10 −14 ), given that where f (A, Z) is a nucleus dependent factor that is ∼ O(1) for the targets used in experiments [1].As shown in the right plot of Figure 6, the scaling of Equation ( 34) is satisfied for small m N /m η ratios, while the non-dipole penguin amplitude can come to dominate for larger m N /m η values.The upcoming µ → e conversion searches will be a valuable probe for this region of parameter space [53].Another popular model that can generate neutrino masses at loop level is the Zee-Babu Model [54,55], where the SM is extended with two SU(2) singlet and charged scalars k + , k ++ and allow for the Lagrangian terms where ϵ I J is the anti-symmetric SU(2) tensor.Lepton number is not conserved and neutrino masses are generated at the two-loop level, while the interactions also violate lepton flavor.
The CLFV phenomenology of the Zee-Babu model has been studied in [56][57][58].For other models that generate neutrino masses at loop level and their CLFV signatures, we refer the reader to [59].

Two Higgs Doublet Model
Neutrino masses are just one among pieces of evidence and hints of new physics.Dark matter, the baryon asymmetry, B anomalies, the hierarchy and the strong CP problem are some of the observations and theoretical conundrums that call for SM extensions.CLFV is a feature of numerous models that address the above problems and can be employed to constrain or, if CLFV is ever observed, determine the region of parameter space where Beyond Standard Model theories must sit.
One simple extension of the Standard Model features an additional scalar doublet H 2 , which is commonly known as the Two Higgs Doublet Model (2HDM) (for a review see [60]).A second Higgs is strongly motivated by supersymmetry, where one Higgs cannot give masses to all fermions and the second Higgsino, the superpartner of the second doublet, is necessary to cancel the gauge anomalies.Although supersymmetry imposes precise relations among the Higgs masses and couplings, supersymmetry breaking terms lead to modifications so that, at low energy, it is suitable to describe the two Higgs with generic couplings.A general 2HDM (type III) predicts LFV couplings that must be sufficiently suppressed to satisfy the current experimental constraint.Often, additional symmetries are assumed to avoid the appearance of flavour-changing neutral current at tree level.In Type I 2HDM, the SM fermions only couple to one Higgs, while, in the type II, the up quarks couple to a different Higgs with respect to leptons and down quarks, which is the case for supersymmetric SM.
The 2HDM scalar Lagrangian is the following where the potential reads In a region of the potential parameters, the Higgs acquire a VEV that spontaneously breaks the electroweak gauge symmetry, and it is always possible to rotate in a basis where only one has a non-zero expectation value ⟨H 1 ⟩ = 0 v T , ⟨H 2 ⟩ = 0.The doublets are written as Once the Goldstones G are eaten by the gauge bosons, the scalar spectrum contains one ϕ + complex scalar, two CP even neutral scalar ρ 1,2 and one CP odd scalar A. If the potential parameters are real, only the two CP even scalars mix, and we identify two mass eigenstates h, H [61,62] where β − α is the angle that diagonalizes the neutral scalar mass matrix, and h, H have masses m h < m H , respectively.We identify the lighter scalar with the 125−GeV Higgs boson.In a basis where the two doublets Φ i both have VEVs v i , a rotation with angle α diagonalizes the neutral scalar mass matrix, while the angle β such that tan β ≡ v 1 /v 2 allows us to rotate into the H i basis.In the Type III 2HDM, there is no unambiguous way to identify β and α because the two doublets are not distinguishable.On the other hand, β − α is calculable in terms of the potential parameters [62] cos In the fermion mass basis, the Yukawa matrices Y f are diagonalized, while the K f couplings with H 2 are, in general, non-diagonal.The Yukawa interactions between the fermions and the uncharged scalar sector read The LHC measure a h → τ + τ − , µ + µ − rate compatible with the Standard Model prediction [63,64], requiring s βα ∼ 1. Substituting this approximation in Equation ( 41) gives In the decoupling limit small, and in the following, we consider The off-diagonal interaction K e c βα can mediate LFV Higgs boson decay with a rate [66,67] Non-observation of LFV decay modes of the Higgs boson at LHC set the upper limits on the branching fractions reported in Table 1 and directly constrain the size of flavour violating coupling.

Process
Bound Off-diagonal Yukawas are also indirectly bounded by other LFV processes, as they can mediate l i → l j γ through the loop diagrams shown in Figure 7.The two-loop diagrams of Figure 7b,c are relevant and can be numerically larger than one-loop contributions [70] because, in the former, the Higgs line is attached to a heavy particle running in the loop and Yukawa suppression is avoided.In the µ → e sector, µ → eγ is the most sensitive process to LFV Yukawas, which has been extensively studied in the context of 2HDM [62,[71][72][73][74].For the τ ↔ l sector, the bound on the radiative decay Br(τ → lγ) < few × 10 −8 → 10 −9 is less stringent and the Higgs LFV decays are sensitive to smaller off-diagonal Yukawas.In a simplified scenario where only the SM Higgs is present, the author of [75] computed several processes in terms of generic LFV couplings hY ij ēLi e Rj , and in Figure 8, we show the current bounds (sensitivity) set by LFV observables.In 2HDM, there are also contributions from the heavy scalars, which are parametrically of similar size with respect to light Higgs LFV; while they do not suffer from the small mixing angle c βα ∼ v 2 /M 2 , the propagator yields a similar suppression ∼ 1/M 2 .[68]: constraint on LFV Yukawa couplings Y µe , Y eµ from the limits on Br(h → eµ) (observed limit corresponds to the solid blue line, while the expected one is the dashed red line).Shaded regions show the sensitivity of µ → 3e, µ → eγ and µ → e conversion on the LFV Yukawas, from the calculations of [75].Right figure from [69]: same plot for the τ ↔ µ sector.The diagonal line shows the natural limit Y ij Y ji < 2m i m j /v 2 [76].
Contribution to l i → l j l k l k and l i → l j qq also appear at tree level in the 2HDM, but they are suppressed by the small diagonal Yukawas.The same processes receive relevant contributions by attaching a l j l j (qq) current to the photon of the diagrams in Figure 7.

CLFV in Supersymmetry
Supersymmetry (SUSY) is a space-time symmetry that extends Poincare invariance by adding fermionic generators that satisfy the anti-commutation relations of the supersymmetry algebra [77][78][79].SUSY is the largest space-time symmetry that the S−matrix can have given a set of physical assumptions such as unitarity, locality and causality [80,81].Since fermionic operators Q are added to the algebra, irreducible SUSY representation (supermultiplets) contain particles of different spin that are related by the action of Q on one-particle states.In the case of N = 1 SUSY, i.e., with only one pair of conjugate Weyl spinor generators, the spectrum contains boson-fermion degenerate pairs.In the Minimal Supersymmetric SM (MSSM), for every quark q and lepton ℓ, there is a corresponding complex scalar in the same gauge representation, commonly known as the squark q and slepton l.Similarly, the gauginos B, W I , Ga are the fermion superpartner of the gauge bosons, which transforms in the adjoint of the gauge group, while the higgsinos Hu , Hd are the spin-1/2 particles that belong to the supermultiplets of the Higgs doublets.As already discussed in the previous section, a supersymmetric version of the SM requires at least two Higgs doublets.
One of the most attractive features of SUSY is that it provides a solution to the hierarchy problem.The Higgs in the SM is a fundamental scalar, and its mass is quadratically sensitive to any new physics that couples to it.If we assume that a new bosonic or fermionic heavy state with masses ∼ Λ couples to the Higgs, its mass will get corrections that are quadratic in Λ.Therefore, the Higgs tree-level mass must cancel almost precisely with the UV contribution, leaving a small remnant mass at the electroweak scale m h ∼125 GeV.If Λ is taken near the Planck scale M Pl ∼ 10 19 GeV, the fine-tuning required for this accidental cancellation is extreme.If nature is supersymmetric, then the corrections to the Higgs mass cancel precisely among the contributions of superpartners.That is because SUSY predicts the same couplings for all particles in the supermultiplet, and closed fermion loops get minus signs from Grassmann traces.The absence of quadratic divergences can also be understood by observing that fermion-boson pairs are degenerate, and fermion masses are protected by the chiral symmetry and are only logarithmically sensitive to the UV scale.
Degenerate partners of the known SM particles with opposite statistics have never been observed; therefore, SUSY, if realized at all, must be broken at some scale m 0 .To avoid the reappearance of the hierarchy problem, m 0 should not be too far from the Higgs mass and explicit SUSY breaking terms must be soft, i.e., have to contain only terms with strictly positive mass dimension.The solution to the hierarchy problem is also preserved if SUSY is spontaneously broken by the expectation value ∼ m 0 of some scalar field.In models of spontaneous breaking, soft breaking terms appear in the low-energy non-supersymmetric description.Null results from the LHC rule out SUSY breaking scales below few × TeV [82][83][84], although the bounds on superpartners' masses are not completely model-independent.
In the MSSM, the SUSY breaking sector can be a source of lepton flavour violation.The soft breaking terms contain masses for the sleptons and trilinear couplings with the Higgs that introduce LFV if the off-diagonal entries are non-zero in the lepton mass eigenstate basis.For ∼ 100 GeV−TeV soft terms, the current bounds on LFV, and more generally on flavour-changing neutral current in the SM, call for a suppression mechanism of sfermions mass mixing.This is known as the SUSY flavour problem.The spontaneous symmetry breaking of SUSY cannot be triggered by the scalar fields in the MSSM supermultiplets, as this would lead to an unacceptable spectrum.Supersymmetry breaking may occur in a hidden sector that has a very small coupling with the MSSM particle.This is known as the mediation paradigm.If SUSY breaking is communicated via supergravity couplings of the hidden sector to matter, it results in universal and flavour-conserving soft terms at the Planck scale [85,86].Nonetheless, this does not strictly forbid LFV, since mass mixing can still be radiatively generated.In minimal SU(5) Grand Unified Theory (GUT), the matter content of SM is reproduced by three generations of a fermion field in the anti-fundamental 5 of SU(5) ( F), which contains the lepton doublet and the right-handed down-type quark, and one that fills the 10 representation of SU(5) (T), which contains the quark doublet, the right-handed up-type quark and the right-handed charged lepton.They are coupled in the Yukawa sector to two Higgs scalar fields transforming in the 5 and 5 representation.(For the following, we adopt the common convention in SUSY of dropping Dirac notation, and we use chiral Weyl fermions.The product of left-handed Weyl fermion contracted with the anti-symmetric tensor In SUSY GUT, the above equation corresponds to the superpotential W, where T, F, H, H are the superfields that contain the SM particles and the superpartners.Assuming gravity-mediated SUSY breaking, the soft terms at M Pl are flavor blind and characterised by a common mass scale m 0 The top Yukawa is large and loop correction to third generation masses can be sizeable.In a basis where the up Yukawa matrix is diagonalized and neglecting first and second generation couplings, the leading-log correction in the renormalization of the T3 mass is [1,87] ∆ mT,33 ≃ − where M GUT ∼ 10 16 GeV is the GUT scale.The T fields contain the right-handed charged sleptons that have a diagonal but non-universal mass matrix.In the mass eigenstate basis for the charged leptons, the right-handed slepton mass matrix acquires non-diagonal entries where Ŷe is the diagonal lepton Yukawa.In SU(5) GUT, the down and lepton Yukawa are unified Y e = Y T d , and V e correspond to the transpose CKM matrix.In the diagrams of Figure 9, we show how slepton mass mixing can mediate l i → l j γ at loop level, which in most SUSY setups is the largest LFV signal.Box diagrams exist for l i → l j l k l k and l i → l j qq, but the processes are often dominated by the penguin diagrams, where a flavour diagonal current is attached to an off-shell photon in the diagrams of Figure 9.The rate of µ → eγ in minimal SU(5) GUT has been calculated in [87,88].A detectable signal is predicted in upcoming experiments, although the values considered for the sparticles masses are in tension with more recent LHC data [89].LFV can be sizeable in the context of GUT theories with right-handed sterile neutrinos, which has been studied in [88,[90][91][92][93][94][95][96][97][98].In SO (10), right-handed neutrinos naturally appear in the 16 spinor representation that an SM generation fills, and neutrino masses can be explained with a supersymmetric seesaw mechanism.Considering heavy right-handed neutrinos, the superpotential in the lepton sector reads where the notation for the SM superfields is self-explanatory, and N is the superfield that contains sterile neutrinos.As in Equation ( 22), for large Majorana masses M R , the light neutrino mass matrix is The gravity-mediated soft breaking terms regarding sleptons are the following and the left-handed sleptons mass matrix is renormalized in the leading-log approximation as [90] [ The typical size of l i → l j γ branching fraction is [93] where m SUSY is the sparticles mass scale.In general, even knowing neutrino masses and mixing angles, the neutrino Yukawa Y ν is not uniquely defined [99].In a basis where the Majorana masses MR are diagonal, we can use the Casas-Ibarra parametrization [93] for where U is PMNS and R is an unknown orthogonal complex matrix.The matrix that controls the slepton mixing is then and depends on the unknown R matrix.Assuming specific mass hierarchy and degenerate patterns for neutrinos, the free parameters in R are reduced, and the predicted LFV signals are studied when the parameters are varied [100].In SO (10) GUT, the neutrino and up Yukawa are unified at the GUT scale, and different breaking scenarios can lead to lepton flavour change governed by CKM or PMNS mixing with the third generation [94,95].
The PMNS angles are large and lead to an insufficiently suppressed µ → eγ rate, larger than the current upper limit Br(µ → eγ) < 4.2 × 10 −13 [101].In the scenario where LFV amplitudes are proportional to CKM matrix elements, the rate is compatible with the future experimental upper bound Br(µ → eγ) < 6 × 10 −14 [102].This model relates the branching ratios of τ → lγ and µ → eγ via the following relation where V is the CKM matrix, and we have substituted the upcoming MEG branching fraction bound.If a τ → µγ signal is observed by BELLE II with Br(τ → µγ) ∼ 10 −9 [103], it can disfavor the model.In the context of sleptons mixing and LFV, several simplified SUSY scenarios have been more recently studied, complemented with the bounds of the null results of the LHC [104].
The soft-breaking sector is not the only possible source of LFV in supersymmetric SM.Gauge and SUSY invariance allow for the following terms in the superpotential: The λ ′′ term is baryon number violating and can lead to prompt proton decay.To avoid this disastrous outcome, a discrete symmetry known as R−parity is often assumed.The R−parity of a particle is defined as (−1) R ≡ (−1) 3(B−L)+2S , where B, L are the baryon and lepton number, while S is the particle spin.It follows that any SM particle is RP−even and the superpartners are RP−odd.RP−invariance automatically forbids all terms in the superpotential of Equation ( 58), but other discrete symmetries such as baryon parity [105] can allow for lepton flavor violation while conserving baryon number.The first two terms in the superpotential leads to the Lagrangian terms [106] L that can allow for several LFV processes already at tree level.In Figure 10a, we show a diagram for the LFV K 0 decay K 0 L → µe, whose branching fraction is constrained by the current upper limit Br(K 0 L → µe) < 4.7 × 10 −12 [107].Assuming only one non-zero pair of R-parity violating coupling λ ′ * ik1 λ ′ jk2 , the bound implies (adapted from [108]) where we have assumed m ũk ∼ few × TeV.Similarly, Br(µ → 3e) < 10 −12 → 10 −16 [109,110] can set the following constraint on the coupling products λ n21 λ * n11 (λ is anti-symmetric in the first two indices and n ̸ = 1) from the diagrams of Figure 10b: λλ ′ diagram give tree-level contributions to µ → e conversion, and at one loop l i → l j γ is sensitive to λλ, λ ′ λ ′ couplings.For a more complete discussion on LFV in R-parity violating theories, we refer the reader to [111][112][113][114][115][116].10.Examples of LFV tree-level diagrams in the supersymmetric SM with R-parity violation.(a) Diagrams contributing to the LFV meson decay K 0 → l i l j .(b) Diagrams contributing to the l i → l j l k l m rate.

Effective Field Theory for Charged Lepton Flavor Violation
Under the assumption that the new physics responsible for CLFV is heavy, i.e., Λ > TeV, Effective Field Theories (EFT) are powerful model-independent descriptive frameworks.In EFT, the UV physics is hidden in the Wilson Coefficients (WC) of nonrenormalizable interactions among light degrees of freedom that are suppressed by inverse power of the new physics scale For a general introduction on EFT, see [117,118].In a "top-down" perspective, the heavy fields of a particular model can be integrated out, and a tower of higher dimensional operators is generated, whose coefficients are expressed in terms of the couplings and parameters of the model.In this approach, the EFT plays the role of a helpful calculational tool, since observables are more easily computed with contact interactions constructed out of light fields rather than in the full theory.In a "bottom-up" analysis, rates are calculated in the most general EFT that is relevant for the processes considered, and bounds or observations are translated to the operator coefficients, which allow to probe general heavy physics model.Bottom-up calculations should include every contribution to which observables can be sensitive.
At fixed order in the EFT expansion in inverse powers of the heavy scale, EFT renormalization proceeds as in any field theory, and the couplings run with the energy scale.Observables are computed in terms of WC at the experimental energy scale that, for the many LFV channels, lies below the electroweak scale.The EFT permits the separation of energy scales and is most helpful when appropriate for the process considered: it should include the dynamical degrees of freedom and respect the symmetries that concern the probed energy scale.Below the electroweak scale, the top quark, the Higgs boson and the electroweak W, Z boson can be integrated out and the EFT features contact interactions among lighter SM particles, respecting QED and QCD gauge symmetries.We reproduce in Table 2 the operator basis of [119] relevant for µ → e LFV processes that are otherwise flavour diagonal.[119].The experiments can be sensitive to three, four point functions that correspond to operators of dimension from five to eight.X αβ = F αβ , G αβ are the field tensors of photon and gluons, respectively.The chiral projector P Y,Z can be Y, Z ∈ {L, R}, L = R, R = L, while q ∈ {u, d, s, b, c} and l ∈ {e, µ, τ}.
Normalizing the operators of Table 2 where the dipole operator coefficient C D,Y (M) are at M = m µ .If we assume that the dipole coefficient is order one if normalized at the new physics scale, Equation (63) implies Λ NP ∼ 1.6 × 10 4 TeV.The bound can be satisfied with lighter UV physics if the dipole coefficient is loop and/or coupling suppressed.
To use the bound on EFT operator coefficients to constrain generic BSM heavy physics, we should determine the upper limit on coefficients at the new physics scale, where the heavy degrees of freedom are integrated out.This is done by solving the renormalization group equations of the Wilson coefficients, which requires dressing the operator basis with QED and QCD loops.QED amounts to a few percent effect in the renormalization of operator coefficients, but it still plays an important role because it can mix operators.QCD running does not mix operators, but the rescaling of quark scalar and tensors coefficients is numerically relevant (few×10%).Operator mixing allows probing an operator coefficient which is difficult to detect via its mixing to a tightly constrained one.For instance, the tensor operator O ττ T,Y = (ēσ αβ P Y µ)( τσ αβ P Y τ) mixes into the dipole by closing the tau legs in a loop and attaching a photon.The contribution to the dipole coefficient is one-loop suppressed (and log enhanced), but, to close the loop, a chirality flip is necessary, and a τ mass insertion enhances the mixing by m τ /m µ .Complemented with a large anomalous dimension, the mixing is ∼ O(1).The (sensitivity) bound that µ → eγ sets on the tensor coefficient at m W is then [119] We should stress that this is not an exclusion bound, but rather an experimental sensitivity.In coefficient space, µ → eγ constrains one direction that corresponds to the dipole coefficient at the experimental scale m µ .The RGEs can tell us how this direction rotates in the coefficient space at higher energies, but the experiment still imposes a bound in one direction only.In other words, the bound will apply to a single combination of operator coefficients at the high scale; namely, solving the RGEs up to m W , The upper limit in Equation ( 64) corresponds to the case where only C ττ T,Y (m W ) is non-zero and is commonly known as one-operator-at-a-time sensitivity.A sensitivity corresponds to the smallest absolute value which is experimentally detectable, but larger values are possible if cancellations with other contributions occur.Indeed, we can see that if the dipole and the tensor are of similar size and opposite sign, an accidental cancellation can occur in Equation ( 65).This is an example of a flat direction in coefficient space.Flat directions are a general feature of bottom-up EFT analyses of LFV, because the operator basis contains more operators than observables, and the few operators constrained by experiments mix with the rest in the RGEs.Nonetheless, identifying operator coefficients to which observables are most sensitive is a useful guide for model building.The sensitivities of µ → 3e, µ → eγ and µN → eN to Wilson coefficients at m W in the low energy EFT has been extensively studied [119][120][121][122][123]. Spin-dependent µ → e conversion in nuclei [121,122,124,125], although less constraining than the spin-independent searches, allows us to probe different combinations of coefficients and, thus, reduce the number of flat directions.
Leptonic and semi-leptonic rare meson decays such as K 0 L → µe, K + → π + µe are systematically studied in the EFT by adding to the operator basis quark flavor-changing operators, and the sensitivities to Wilson coefficient can be similarly determined [126].Concerning τ ↔ l processes, EFT analysis can be found in [119,127].
Above the electroweak scale, the appropriate EFT is the Standard Model Effective Field Theory (SMEFT), which contains all SM particles, and the operators are SU(3) ⊗ SU(2) ⊗ U(1) invariant.The contact interactions in the low energy EFT for LFV are matched onto SMEFT operators of dimensions between six and eight [128][129][130][131].Running the experimental bound from m W to a higher scale requires solving the RGEs of the SMEFT operators [132][133][134], identifying the contributions to which LFV observables are sensitive.Establishing all the relevant contributions in SMEFT is not trivial, and a complete bottom-up analysis of LFV requires some challenging calculations that are currently missing in the literature [135].

Experimental Review
Searches for evidence of CLFV signals span a broad range of experimental techniques thanks to the large variety of processes one could be looking for, such as rare muon and tau decays (µ + → e + γ , µ ± → e ± e − e + , τ ± → µ ± γ , τ± → e ± γ , τ → 3l ), rare mesons and bosons decays, and direct conversions of a lepton in a nuclear field (µ − N → e − N , µ − N → e + N ′ ).Table 3 summarizes the current best limits on the various channels.The CLFV searches based on muons have been performed with dedicated experiments, usually highly tuned for a specific channel, which took advantage of the facilities capable of delivering a high intensity muon beam (see next section).For all the other cases (tau, mesons and bosons), with the only exception made for the Kaons, it is not possible to deliver a dedicated beam; thus, general-purpose detector systems have been used.
As discussed in the previous sections, the most stringent constraints on various BSM models are set by the direct searches of CLFV decays of muons and taus decays.In the following sections, we describe the most recent and the coming experimental efforts for these two categories: (i) searches using muons, (ii) searches using taus.For each search, a discussion of the peculiarities of the signal topology and of the various experimental challenges is provided.

CLFV Searches Using Muons
In the history of CLFV experiments, muons have been, so far, the most popular.Historically, the first experiment looking for CLFV using muons was performed by Hinks and Pontecorvo using atmospheric muons [158].Since then, the advancements in the muon beam production/acceleration technology at different facilities (PSI, TRIUMPH, LANL, etc.) made available high-intensity muon beams at the level of 10 8 (10 7 )µ + (µ − )/s [3,159], enabling the possibility to search for rare CLFV processes.Facilities under construction at Fermilab (USA) and J-PARC (Japan) [2], or planned, like the High Intensity Muon Beam project at the PSI (Switzerland) [160], have been designed to provide muon beams with an intensity of about 10 10 µ/s.This planned intensity corresponds to 2-3 orders of magnitude improvement with respect to the current state-of-the-art technology.The J-PARC and Fermilab muon experiments will use a novel method for creating the muon beam.A prototype muon beamline, the Muon Science Innovative Channel (MuSIC), was set up at the Research Center for Nuclear Physics (Osaka, Japan) to prove the conceptual idea.The production of an intense muon beam relies on the efficient capture of pions (from protontarget interactions), which subsequently decay to muons, using a novel superconducting solenoid magnet system [161].
The current best limits on the muon-CLFV processes come from experiments that performed dedicated searches for the following processes: µ + → e + γ , µ ± → e ± e − e + , µ − N → e − N and µ − N → e + N ′ .Table 3 summarizes these results.One thing to notice is that all the searches, except the µ − N → e − N and µ − N → e + N ′ , were performed using µ + rather than µ − .This choice is motivated by several advantages: (i) µ + cannot get captured in nuclei, while µ − can undergo nuclear capture events, which produce protons, neutrons, photons and, thus, increase the activity in the detector deteriorating its performance, (ii) the muon beam is obtained from charged pions decay, which are produced in proton-target interactions where π + production is larger; thus, the resulting µ+ beam intensity is higher.
The following sections offer a more detailed description of each of these experimental searches.

µ + → e + γ
In the µ + → e + γ decay, the final state consists of a back-to-back positron and photon with an energy of 52.8 MeV.The background sources for this search can be factorized into two main categories: (i) an intrinsic physics background from the Radiative Muon Decay (RMD) process µ + → e + γν e νµ , where the neutrinos carry off small momenta, and (ii) an "accidental" background where a positron from the Michel decay µ + → e + ν e νµ , together with a photon from an RMD event or an electron-positron annihilation in flight or an e − N nucleus scattering, recreate the topology of the µ + → e + γ decay.While signal and RMD rates are proportional to the muon stopping rate R µ , the accidental background rate is proportional to R 2 µ because both particles come from the beam; the accidental background is, therefore, the dominant enemy of this search.Thus, a continuous muon beam is better suited than a pulsed beam to avoid stripping particles in short bunches, and R µ must be carefully chosen to optimize the sensitivity.The µ + → e + γ searches from the last decades confirmed that the accidental background is dominant, while the intrinsic background accounts for about 10% of the total background budget.
Two different strategies have been applied for designing the experimental apparatus for the µ + → e + γ search: (i) a tracking-only system equipped with a converter to convert the photon in an e + e − pair, or (ii) a tracker combined with a calorimeter for the photon detection.The tracking-only solution has a much better resolution but a cost of a loss in acceptance because converting the photon requires material that spoils the resolution (due to energy loss and multiple scattering) but too little limits the size of the data sample.
One of the first experiments to adopt the calorimetric solution for the photon detection was the Crystal Box experiment at Los Alamos Meson Physics Facility (LAMPF) [162].The experiment, shown in Figure 11, used a surface muon beam at LAMPF with an average intensity of 300 kHz.The detector consists of a cylindrical drift chamber surrounded by 396 NaI(Tl) crystals.A layer of scintillation counters in front of the crystals provided a timing measurement for the electrons and a veto for photons.The energy resolution for electrons and photons was ∼6% (FWHM).The position resolution of the drift chamber was 350 µm, while the time resolution was ∼400 ps for the scintillators and ∼ 1 ns for the crystals.A total of 3 × 10 12 muon were stopped in a thin polystyrene stopping target.A maximum likelihood analysis established a 90% C.L. upper limit of 4.9 × 10 −11 [162].
The next-generation experiment, MEGA [163], was also performed at LAMPF.The MEGA experimental apparatus, shown in Figure 12, used a surface muon beam at the stopped muon channel at LAMPF that was stopped in a 76 µm Mylar foil centered in the 1.5 T magnetic field of a superconducting solenoid.The MEGA detector consisted of a magnetic spectrometer for the positron and three spectrometers for the photon, therefore sacrificing the signal acceptance and efficiency for a better resolution and background rejection.In total, 1.2 × 10 14 muons were stopped during the life of the experiment, and the overall efficiency for the observation of the µ + → e + γ event was ∼3.9 × 10 −3 .The small efficiency was due to the photon conversion probability (∼2.5%) and to the reduced capability of reconstructing the positron tracks in the solenoidal field compared to the design value.For these reasons, the final sensitivity reached by the MEGA experiment, 1.2 × 10 −11 @ 90% C.L. [163], was ∼35 times worse than the design value, proving how challenging it is to deliver progress in this type of search.
The current best limit for the µ + → e + γ branching ratio, 4.2 × 10 −13 @ 90% C.L., comes from the MEG experiment [101].The detector system, shown in Figure 13, covers ∼10% of the solid angle and surrounds a 205 µm-thick polyethylene muon stopping target.The apparatus consists of a positron spectrometer and a liquid-xenon (LXe) calorimeter.
MEG opted for no converter for the photon detection, the opposite of MEGA.This choice avoids the pileup problem in the pattern recognition that limited MEGA but, at the same time, limits the geometrical acceptance.Table 4 summarizes the detector performance measured during the MEG operation [102].A key feature of MEG is the magnetic field design.MEG adopted a graded solenoidal field, set at ∼1.1 T near the center of the apparatus, that sweeps out the positrons emitted at ∼90 deg and provides a constant bending radius for the signal positron essentially independent of the angle of emission.This feature helps in achieving a uniform and efficient signal track reconstruction.Another technological breakthrough from the MEG experiment is the development of the liquid Xe (LXe) calorimeter.The MEG LXe calorimeter is the first application of a large volume of LXe for particle detection and, so far, it proved to have the best performance for the electromagnetic calorimetry detection in the energy range below 100 MeV [164].[102].σ X indicates the resolution of the observable X, ϵ X the detection efficiency for the particle X.For the case of the photon energy resolution σ E γ , the two values refer to the shallow (<2 cm)/deep (>2 cm) events.σ t e + γ is the time resolution on the e + − γ time residual.The reported values for the MEG-II case refer to the updated results from the engineering runs reported in [102].Recently, the MEG collaboration worked on the upgrade of the experiment (MEG II), which aims to reach a sensitivity of 6 × 10 −14 90% C.L. [102].Various improvements on the detector were delivered.The positron spectrometer was replaced with a low-mass singlevolume cylindrical drift chamber with high rate capability.This increased the acceptance of the spectrometer with respect to the MEG configuration by more than a factor of 2. The LXe calorimeter was also upgraded by replacing the MEG photomultiplier tubes (PMTs) with smaller vacuum-ultraviolet sensitive silicon photomultipliers (SiPMs).A novel timing detector for an active suppression of the accidental background was also introduced.The results of the engineering runs showed a fast degradation of the wires of the drift chamber and of the SiPMs [102].Table 4 compares the new detector performance with the previous ones reported for the MEG detector.The MEG-II collaboration plans to build a new chamber to replace the existing one, and they will take advantage of the coming engineering runs to study more carefully the degradation of the SiPMs.Preliminary results show that they can adjust the operation conditions to achieve the desired level of sensitivity [102].

Table 4. Summary of detector performance for the MEG and MEG-II experiments
In the µ ± → e ± e − e + decay, the final state consists of two positrons and one electron emerging from the same vertex with an invariant mass that matches the muon rest mass.In a three-body decay, the energy associated to each product is not a fixed amount.Simple relativistic kinematics consideration show that the maximum energy of one of the decay products is about m µ /2 and that the decay can be described by two independent variables.The energy distribution of each daughter particle depends on the exact dynamics of the underlying unknown physics.In general, the highest energy particle is expected to have a momentum larger than 35 MeV/c, while the distribution of the lowest energy particle peaks near zero and decreases quickly as its energy tends to its upper limit so that only about one half have an energy larger than 15 MeV [3].The background sources for this search can be factorized in two main categories: a physical background coming from the µ + → e + νµ ν e e − e + process, and an uncorrelated component coming from the accidental coincidence of a positron from a Michel decay and a positron-electron pair produced by the interactions of other positrons or muon with the target or the detector material.The accidental background component scales quadratically with the muon beam intensity.As in the µ + → e + γ case, it is more convenient to design an experimental apparatus that uses positive muons.
The current best limit on µ ± → e ± e − e + , 1.0 × 10 −12 [136] at 90% C.L., was set by the SINDRUM experiment at PSI [136] based on ∼10 6 stopped µ + .The SINDRUM apparatus, shown in Figure 14, consisted of a double cone-shaped stopping target in the middle of five concentric multi-wire proportional chambers surrounded by an array of plastic scintillator counters inside a solenoidal magnetic field.For a 50 MeV electron/positron, the detector apparatus had a momentum resolution at the level of ∼1 MeV/c, a timing resolution ≤ 1 ns and a vertex resolution of ∼1 cm.The data reduction was achieved with a multiple stage trigger, taking advantage of track and charge pre-filters that were requiring at least one negatively and two positively charged tracks within a time window of 7 ns.Then, a track-correlator was used to limit the total transverse momentum of the e + e − e + triplet below 17 MeV/c.In the statistical analysis, the event candidates were determined from the two-dimensional distribution of ∑ E i vs. p2 , where p = (p L /σ L ) 2 + (p T /σ T ) 2 (L and T denote the longitudinal and transverse components with respect to the beam axis).This parametrization is particularly convenient because the signal candidates satisfy ∑ E i = m µ and p2 is expected to peak near 0.  A new effort to improve the sensitivity on µ ± → e ± e − e + search is underway at PSI by the Mu3e collaboration [109].The Mu3e experiment aims for a 10 −16 single-event sensitivity, which would correspond to an improvement by four orders of magnitude compared to the limit set by the SINDRUM experiment.Such a leap in sensitivity is enabled by: (i) the availability of high-intensity muon beams, (ii) the use of silicon pixel detectors instead of multi-wire proportional chambers to track the decay products, and (iii) a modern dataacquisition system able to handle the vast amount of data produced by the detector.A first phase of the experiment is currently under construction at the πE5 beamline at PSI, where the intense DC surface muon beam of 10 8 µ + /s will be exploited to achieve a single event sensitivity of 2 × 10 −15 in about 300 days of data taking [110].The Mu3e experimental setup is shown in Figure 15.It is designed to track the two positrons and one electron from the positive muon decaying at rest with a light-weight tracker placed inside a 1 T magnetic field, thereby reconstructing the decay vertex and invariant mass.The muon beam is stopped in a hollow double-cone target placed at the center of the Mu3e solenoid.This allows for the spread out of the decay vertices in z and minimizes the amount of target material traversed by the decay particles.The target is surrounded by the cylindrical central tracker, which consists of an inner silicon pixel detector, a scintillating fiber tracker for time measurements, and an outer silicon pixel detector.A momentum resolution of better than 1 MeV/c at @ 50 MeV/c is achieved by letting the positrons (electrons) re-curl in the magnetic field, either crossing the central tracker again or hitting the outer tracking stations surrounding the upstream and downstream beam pipe.These stations consist of a silicon pixel tracker and a scintillating tile detector mounted on the inside of the pixel tracker.The 5 mm thick tiles enable a time resolution for the tracks reaching these outer stations of better than 100 ps.The material budget, which must be minimized to reduce the multiple scattering and thus deliver the required momentum resolution, was minimized by means of custom High-Voltage Monolithic Active Pixel Sensor (HV-MAPS) [166] based on a commercial 180 nm HV-CMOS process.Together with its support structure, the entire silicon tracking module has a thickness of ∼0.12% radiation lengths, with a single-hit efficiency > 99% and a time resolution of O(10 ns).A gaseous helium cooling system allows the experiment to dissipate 250 mW/cm 2 of power generated by the MAPS modules.A time resolution of about 10 ns is insufficient to determine the direction and thus the charge of the decay particles.A scintillating fiber detector is, therefore, placed between the inner and outer layer of the central silicon-pixel tracker, consisting of a dozen 30 cm long ribbons made from three staggered layers of 250 µm diameter multi-clad round fibers, read out by Silicon Photomultipliers (SiPM) arrays on both sides.Located at the very end of the re-curling particle trajectories hitting the upstream or downstream tracker, where the constraints on the material budget are less stringent, the tile detector provides the needed precise timing information of the particle tracks, in conjunction with the fiber detector, significantly reducing the accidental background.Each tile is read out by a single SiPM.For the tile and fiber detector, a time resolution of <50 ps and <400 ps is achieved, respectively.Mu3e had a successful integration run campaign from May to July 2021 with a reduced detector: 2 pixel layers + fiber detector.
conversion is the process where a muon converts into an electron in the field of a nucleus without producing neutrinos in the final state.This process has the same dynamic of a two-body decay and, therefore, results with a monochromatic electron with an energy E µe : where m µ is the muon mass, E b ∼ Z 2 α 2 m µ /2 is the muonic binding energy and the last term is from nuclear recoil energy up to terms of order 1/m 2 N , neglecting variations of the weak-interaction matrix element with energy [2], where E µ = m µ − E b and m N is the atomic mass.In the case of Al, which is the selected material for the current experiments under construction, E µe ∼ 104.96MeV.In muon conversion experiments, the quantity is measured.The normalization to captures simplifies calculations as many details of the nuclear wavefunction cancel in the ratio [167].The coherent conversion leaves the nucleus intact, and there is only one detectable particle in the final state.As we will see, the resulting electron energy stands out from the background, hence muon-electron conversion does not suffer from accidental background, and extremely high rates can be used.Negative muons stopped in the stopping target can undergo a nuclear capture.Particles generated in the muon capture (n, p and γ) may reach the detector system and create extra activity that can either obscure a conversion electron (CE) track or create spurious hits.As a result, some specific shielding is required to reduce this background.Electrons from the high momentum tail of the muon decay-in-orbit (DIO) represent the intrinsic background source for the µ − N → e − N search.Figure 16 shows the energy spectrum of DIO electrons [168].
The main features of the DIO energy spectrum can be summarized as follows: • the endpoint of the spectrum corresponds to the energy of the electrons from µ − N → e − N conversion (CE); • the overall spectrum is falling as (E µ − E) 5 , where E is the DIO energy; • about 10 −17 of the spectrum is within the last MeV from the endpoint.
Therefore, to reach a high sensitivity at the level O(10 −17 ), the detector resolution is crucial.As the muon beam is generated from charged pions, another relevant background comes from the radiative pion capture (RPC) process π − N → γN * , followed by the electron-positron pair conversion of the photon.Unfortunately, not all pions decay in the transport line, and, consequently, the muon beam is contaminated by pions.This source of background is reduced by taking advantage of the difference between the pion and the muonic atom lifetimes.The pion has a decay constant τ < few tens of ns, while the bound muon has a τ of the order of several hundreds of ns (depending on the Z of the material).Therefore, using a pulsed beam structure, it is possible to set a live gate delayed with respect to the beam arrival, reducing the RPC contribution to the desired level.Other beam-related sources of background are: remnant electrons in the beam that scatter in the stopping target, muon decays in flight and antiprotons interacting in the apparatus.Atmospheric muons can also represent a significant source of background because these particles can interact in the apparatus and eventually generate a signature very similar to the CE.An active shielding is thus required to detect the incoming cosmic muons crossing the apparatus and veto the event candidates on time.Moreover, the detector system has to provide particle identification (PID) capabilities to reject un-vetoed muons that can mimic the CE due to a mis-reconstruction.[169].Figure from [168].
The current best limit on the µ − N → e − N measurement comes from the SINDRUM-II experiment at PSI [138].In SINDRUM-II, a high intensity muon beam was stopped in a target that was surrounded by the detector elements housed in a superconducting solenoid.Figure 17 shows a sketch of the SINDRUM-II apparatus.The detector consisted of two drift chambers, to reconstruct the trajectories of the charged particles, and Cerenkov hodoscopes, to measure the timing of the reconstructed tracks and for providing PID capabilities.With a total of ∼10 14 stopped muons, SINDRUM-II reached a sensitivity at the level of ∼10 −13 on the µ − N → e − N process using different target materials [138].
New experimental concepts have been proposed and are currently under construction at Fermilab (USA) and J-PARC (Japan) to search for µ − N → e − N with unprecedented sensitivity at the level of ∼10 −17 .The Mu2e experiment at Fermilab had its genesis back in the 1980s, behind the Iron Curtain.In a way, Mu2e was born in the Soviet Union.In 1989, the Soviet Journal of Nuclear Physics published a letter to the editor from physicists Vladimir Lobashev and Rashid Djilkibaev, where they proposed an experiment that would perform the most thorough search yet for muon-to-electron flavor violation.In 1992, they proposed the MELC experiment at the Moscow Meson Factory [170], but in 1995, due to the political and economic crisis, the experiment shut down.The same overall scheme was subsequently adopted in the Brookhaven National Laboratory MECO proposal in 1997 [171] and then in the Mu2e and COMET experiments.
The Mu2e apparatus [172], shown in Figure 18, consists of three main superconducting solenoids.The first two, named production and transport solenoid in Figure 18, are used to generate a high-intensity, low-momentum muon beam starting from a 8 GeV proton beam.The third solenoid, named "Detector Solenoid" in Figure 18, contains an Al stopping target, where the muons are stopped to generate the muonic atoms, and downstream to it, we have a low-mass straw-tube tracker [173], followed by a pure-CsI crystal calorimeter [174].Both detectors are left un-instrumented in the inner 38 cm to avoid any interaction with the largest majority (>99%) of the low momenta electrons coming from the muon DIO processes in the stopping target.In Mu2e, the stopping target is not placed in the middle of the tracker as it was done in SINDRUM-II to limit the flux of protons, photons and neutrons (from the muon nuclear captures) in the detector.A graded magnetic field around the stopping target increases the detector geometrical acceptance by reflecting the electrons that initially were emitted in the direction opposite to the detector.The whole detector solenoid and half of the transport solenoid are covered with a cosmic ray veto system designed to detect atmospheric muons with an efficiency ≥99.99%.The design of the COMET experiment at J-PARC, shown in Figure 19, is based on a similar concept.A 8 GeV pulsed proton beam is used to produce pions, which are then captured and transported by a series of superconducting solenoids.The pions decay into muons as they travel along the muon transport channel.The toroidal field of the muon transport channel selects muons with negative charge and momentum less than 75 MeV/c.The major difference with respect to the Mu2e design is that a second transport line is installed between the muon stopping target and the detector regions to select charged particles of momentum centered around 100 MeV/c.The detector system consists of a straw-tube tracker followed by a LYSO crystal calorimeter [175].
COMET plans to operate in two stages: Phase-I and Phase-II.Phase-I will allow the experiment to characterize the beam and the key backgrounds as well as provide enough statistic to reach a 90% C.L. sensitivity at the level of 7 × 10 −15 [175].During Phase-I, COMET will operate with a smaller apparatus that consists of half of the first C-shaped muon transport line directly connected to a solenoid that houses the muon stopping target surrounded by the detector system.For Phase-I, the detector consists of a cylindrical drift chamber and a set of trigger hodoscope counters.
Another experiment, named DeeMe [176], aims to search for the µ − N → e − N process with a single event sensitivity of 1 × 10 −13 using a graphite target.The experiment is conducted at the Materials and Life Science Experimental Facility (MLF) of the J-PARC.Muonic atoms are produced in a primary-proton target itself, which is hit by pulsed proton beams from the Rapid Cycling Synchrotron (RCS) of J-PARC.To detect the electron and measure its momentum, a magnetic spectrometer, consisting of a dipole magnet and four sets of multi-wire proportional chambers (MWPCs) [177], is employed.The spectrometer is expected to reach a resolution of σ p < 0.5 MeV/c at 100 MeV/c.The resolution is needed to reject the DIO background, which is the dominant source of high energy electrons for this search.The number of charged particles hitting the detectors is estimated with Monte Carlo simulation to be approximately 10 8 particles per proton-bunch with an RCS power of 1 MW.The construction of the secondary beamline for DeeMe, the H Line, is now in progress.Meanwhile, the collaboration measured the momentum spectrum of the DIO electrons in the momentum region 48-62 MeV/c at the D2 area at MLF.This measurement will be important for validating the theoretical models used to model the DIO background and characterize the detector performance.Three sets of measurements were performed between the year 2017 and 2019 [176], and the analysis is now underway.3.1.4.µ − N → e + N ′ µ − N → e + N ′ conversion is the process where a muon converts into an positron in the field of a nucleus that undergoes a nuclear transition.This process violates the lepton number (∆L = 2) and the lepton flavor conservation.The experiments looking for the µ − N → e − N process can typically search for the µ − N → e + N ′ as well.The current best limit on the µ − N → e + N ′ process comes from the SINDRUM-II experiment [139] that set a limit at 5.7 × 10 −13 at 90% C.L. The major background source is the radiative muon capture, where the photon can generate (via asymmetric conversion) a positron with an energy close to the signal region.
The search for the µ − N → e + N ′ complements the 0νββ decay searches and is sensitive to potential flavor effects in the neutrino mass-generation mechanism.We refer the reader to [178] for additional information about the current status and future prospects offered by the COMET and Mu2e experiments.
From the experimental point of view, however, a difficulty immediately arises: the tau is an unstable particle, with a very short lifetime (τ = 2.91 × 10 −13 s).As a result, tau beams cannot be realized, and large tau samples must be obtained in intense electron or proton accelerators, operating in an energy range where the tau production cross section is large.At e + e − and pp collider machines, the majority of the tau particles are not produced at rest, which means that, unlike the muon searches discussed before, here we need to deal with decays-in-flight.Thanks to the boost, the decay products could get energy values up to several GeV, which experimentally poses the challenge to deliver wide-range calibrations for the detectors (from a few hundreds of MeV to several GeV).For all these searches, events contain a pair of taus in which one tau undergoes SM decay (tag side), while the signal side is selected on the basis of the appropriate topology of each individual channel.The tagging side accepts the leptonic (τ → lν ν) and 1-prong hadronic decays, while on the signal side, CLFV candidates are selected on the basis of the appropriate topology of each individual channel.For each channel, we report the HFLAV combined limit and the experimental published limits.In some cases, the combined limit is weaker than the limit published by a single experiment.This arises since the CLs method used in the combination can be more conservative compared to other legitimate methods, especially when the number of observed events fluctuates below the expected background [182].
The following paragraphs discuss the current best limits for some of these experimental searches from experiments at B-factories and pp colliders.

τ → lγ
The τ → lγ decay, where l is a light lepton (e, µ), has been one of the most popular CLFV tau channels.The signal is characterized by a l ± − γ pair with an invariant mass and total energy in the center-of-mass (CM) frame (E CM ) close to m τ = 1.777GeV and √ s/2, respectively.The dominant irreducible background comes from τ-pair events containing hard photon radiation and one of the τ leptons decaying to a charged lepton.The remaining backgrounds arise from the relevant radiative processes, e + e − → e + e − γ and e + e − → µ + µ − γ and from hadronic τ decays where a pion is misidentified as an electron or muon.For this decay channel, the current best limits comes from the BaBar and the Belle collaborations.BaBar collected (963 ± 7) × 10 6 τ decays near the Υ(4S), Υ(3S) and Υ(2S) resonances.In the BaBar detector [179], charged particles are reconstructed as tracks with a 5-layer silicon vertex tracker and a 40-layer drift chamber inside a 1.5 T solenoidal magnet.A CsI(Tl) electromagnetic calorimeter is used to identify electrons and photons.A ring-imaging Cherenkov detector is used to identify charged pions and kaons.The flux return of the solenoid, instrumented with resistive plate chambers and limited streamer tubes, is used to identify muons.Signal decays are identified by two kinematic variables: the energy difference ∆E = E CM − √ s/2 and the beam energy constrained τ mass obtained from a kinematic fit after requiring the CM τ energy to be √ s/2 and after assigning the origin of the γ candidate to the point of closest approach of the signal lepton track to the e + e − collision axis (m BC ). Figure 21 shows the distributions of the events for the two decay channels in m BC vs. ∆E.The red dots are experimental points, the black ellipses are the 2σ signal contours and the yellow and green regions contain 90% and 50% of MC signal events.The searches yield no evidence of signals, and the experiment set upper limits on the branching fractions of B(τ ± → e ± γ) < 3.3 × 10 −8 and B(τ ± → µ ± γ) < 4.4 × 10 −8 at 90% confidence level [140].
The Belle experiment [180] reported comparable limits using a data analysis based on 988 fb −1 and a strategy similar to that of BaBar.Kinematical selections on missing momentum and opening angle between particles are used to clean the sample.Figure 22 shows the two-dimensional distribution of∆E/ √ s vs. m BC .The signal events have m BC ∼ m τ and ∆E/ √ s ∼ 0. The most dominant background in the τ ± → µ ± γ (τ± → e ± γ ) search arises from τ + τ − events decaying to τ ± → µ ± ν µ ν τ (τ ± → e ± ν e ν τ ) with a photon coming from initial-state radiation or beam background.The µ + µ − γ and e + e − γ events are subdominant, with their contributions falling below 5%.Other backgrounds such as two-photon and e + e − → q q (q = u, d, s, c) are negligible in the signal region.No significant excess over background predictions from the Standard Model is observed, and the 90% C.L. upper limits on the branching fractions are set at B(τ ± → µ ± γ) ≤ 4.2 × 10 −8 and B(τ ± → e ± γ) ≤ 5.6 × 10 −8 [183].With the full dataset expected for the Belle II experiment [184] (the upgrade of Belle), 50 ab −1 , the upper limit for the branching fraction of LFV decays τ will be reduced by two orders of magnitude.

τ → 3l
The signature for τ → 3l (l = e, µ) is a set of three charged particles, each identified as either an e or a µ, with an invariant mass and energy equal to that of the parent τ lepton.
In the BaBar [185] and Belle [141] analyses, all the six different combinations were explored.Events are preselected requiring four reconstructed tracks and zero net charge, selecting only tracks pointing toward a common region consistent with τ + τ − production and decay.The polar angles of all four tracks in the laboratory frame are required to be within the calorimeter acceptance range, to ensure good particle identification.The search strategy consists of forming all possible triplets of charged leptons with the required total charge and of looking at the distribution of events in the (m BC , ∆E) plane (m BC and ∆E are defined as in the previous section).The backgrounds contaminating the sample can be divided in three broad categories: low multiplicity e + e − → q q (q = u, d, s, c) events, QED events (Bhabha or µ + µ − depending on the specific channel) and SM τ + τ − events.These background classes have distinctive distributions in the (m BC , ∆E) plane.The e + e − → q q (q = u, d, s, c) events tend to populate the plane uniformly, while QED backgrounds fall in a narrow band at positive values of ∆E, and τ + τ − backgrounds are restricted to negative values of both ∆E and m BC due to the presence of at least one undetected neutrino.Figure 20 shows the resulting limit for all the combinations to be at the level of a few 10 −8 for both collaborations.
Even if the results are not yet competitive to those from B-factories, it is interesting to note that experiments at the LHC have also been looking for the τ → 3µ decay.The ATLAS experiment [186] performed a search for the neutrinoless decay τ − → µ − µ + µ − using a sample of W − → τ − ντ decays from a dataset corresponding to an integrated luminosity of 20.3 fb −1 collected in 2012 at a center-of-mass energy of 8 TeV.The LHCb experiment [187] performed the same search using a sample of tau from b and c-hadron decays from a dataset corresponding to an integrated luminosity of 3.0 fb −1 collected by the LHCb detector in 2011 and 2012 at center-of-mass energies of 7 and 8 TeV, respectively.The CMS experiment [188] recently delivered the results for the same search using a sample of τ leptons produced in both W boson and heavy-flavour hadron decays from a dataset corresponding to an integrated luminosity of 33.2 fb −1 recorded by the CMS experiment in 2016 [188].ATLAS, CMS and LHCb reported a 90% C.L. upper limit on the branching ratio of 3.76 × 10 −7 , 8.0 × 10 −8 and 4.6 × 10 −8 , respectively.The Belle-II collaboration studied prospects for the expected sensitivity of this search.This channel has a purely leptonic final state, thus it is expected to be free of background.This allows to scale the experimental uncertainties linearly with the luminosity, which means that at least an improvement of a factor ×50 is expected for Belle-II after accumulating a luminosity of 50 ab −1 [103].

Conclusions
This review intended to provide a general theoretical and experimental overview of CLFV processes.CLFV processes are predicted by a large spectrum of BSM theories.Limits set by the experiments are powerful to rule out or heavily constrain the parameters space of many of these models.In the theory sections, we discussed the CLFV phenomenology of several heavy physics scenarios.We reviewed CLFV signature of BSM theories that account for neutrino masses, discussing tree-level seesaw and models that generate neutrino masses radiatively, such as the scotogenic and Zee-Babu models.Furthermore, we studied CLFV in the 2HDM and supersymmetric SM, and we presented a bottom-up analysis of LFV in effective field theories.
We also highlighted the state-of-the-art experiments involved in direct CLFV searches in the muon sector.Aiming to measure branching ratios below <10 −11 requires very careful optimization of the apparatus to limit the background contamination.The most recent results from the CLFV searches involving taus performed by experiments at the LHC and B-factories were also included.
In the coming decade, we expect to see many results delivered by new muon CLFV direct searches (MEG-II, Mu3e, Mu2e, COMET) and by other experiments at the LHC and B-factories that could potentially improve the sensitivity to levels where many BSM theories expect a signal.This is possible thanks to formidable improvements in the muon beamline, novel detector and accelerator technologies developed to face various experimental challenges.
The state-of-the-art calculations for the conversion rate can be found in [1,167].In their notation, we describe coherent and spin-independent µ → e conversion with LFV contact interactions among leptons and light quark currents q=u,d,s (g LS(q) ēP L µ + g RS(q) ēP R µ) qq (g LV(q) ēγ α P L µ + g RV(q) ēγ α P R µ) qγ α q + h.cwhere F αβ is the photon field tensor.Contributions that depend on the spin-state of the nucleus arise from contact interactions involving axial-vector, pseudo-scalar and tensor quark currents.The spin-independent rate is dominant because it is enhanced by the coherent sum over all nucleons.The effective Lagrangian at the quark-level can be matched onto interactions involving nucleons via the following matrix elements where Γ S = 1, Γ V = γ α and p, n label protons and neutrons, respectively.The vector coefficients are obtained by the quark content of the nucleon = 0 while the scalar charges G S are extrapolated with dispersive relations and lattice results [189,190].The wave function for the muon-bound state is calculated by solving the Dirac equation in the presence of the electric field of the nucleus and is averaged over the proton and neutron densities.Defining g(p) XK = ∑ q G (p,q) K g XK(q) and g(n) XK = ∑ q G (n,q) K g XK(q) , the conversion rate reads where D, S, V are the overlap integrals in Equations ( 19)-( 23) of [167], that involve proton/neutron densities and muon/electron wave functions.All the overlap integral scale with the atomic number ∼ Z, giving the anticipated coherence factor Z 2 that enhances the spin-independent µ → e conversion rate over the muon capture rate.

Figure 2 .
Figure 2. Matching contributions to the Weinberg operator in type I seesaw.

Figure 3 .
Figure 3. Seesaw Majorana neutrino masses generated by integrating out a heavy scalar triplet ∆ (type II) or a heavy fermion triplet Σ (type III).

Figure 4 .
Figure 4. Radiative neutrino mass in the scotogenic model.The loop and λ 5 suppression allows for TeV−scale new physics and small neutrino masses. l

Figure 6 .
Figure 6.The plots show some CLFV branching fractions in the scotogenic model: in the left plot the prediction for Br(µ → eγ) and Br(µ → 3e) for degenerate sterile neutrino masses m N = 4 TeV and m η + = 1 TeV, varying the mass of the lightest neutrino (normal ordering); in the right figure Br(µ → eγ) and Br(µ → e) conversion as a function of (m N /m η + ) 2 .The dashed lines correspond to the current experimental upper limit.Yukawas Y η compatible with neutrino parameters are randomly generated (Figure from [51]).

Figure 7 .
Diagrams for l i → l j γ in the 2HDM, where ϕ = h, H, A. Two-loop Bar-Zee diagrams with a Z exchange also exist.(a) One loop contribution to the l i → l j γ rate in the 2HDM with LFV Yukawa couplings.(b) Two loop Barr-Zee diagram with a top loop contributing to l i → l j γ in the 2HDM with LFV Yukawa couplings.(c) Two loop Barr-Zee diagram with a W loop contributing to l i → l j γ in the 2HDM with LFV Yukawa couplings.

Figure 8 .
Figure 8. Left figure from[68]: constraint on LFV Yukawa couplings Y µe , Y eµ from the limits on Br(h → eµ) (observed limit corresponds to the solid blue line, while the expected one is the dashed red line).Shaded regions show the sensitivity of µ → 3e, µ → eγ and µ → e conversion on the LFV Yukawas, from the calculations of[75].Right figure from[69]: same plot for the τ ↔ µ sector.The diagonal line shows the natural limit Y ij Y ji < 2m i m j /v 2[76].

Figure 14 .
Figure 14.Schematic of the SINDRUM experiment.B, muon beam; S, focussing solenoid; T, target; C, five cylindrical multi-wire proportional chambers; H, hodoscope of 64 scintillators; L, light guides for the hodoscope; P, 128 photomultipliers; A, preamplifiers for the cathode strips and amplifier/discriminators for the anode wires; M, normal conducting coil of the magnet.Figure and caption from [165].
Figure 14.Schematic of the SINDRUM experiment.B, muon beam; S, focussing solenoid; T, target; C, five cylindrical multi-wire proportional chambers; H, hodoscope of 64 scintillators; L, light guides for the hodoscope; P, 128 photomultipliers; A, preamplifiers for the cathode strips and amplifier/discriminators for the anode wires; M, normal conducting coil of the magnet.Figure and caption from [165].
limits on τ LFV decays

Figure 20 .
Figure 20.Tau lepton-flavor-violating branching fraction upper limits combinations summary plot.For each channel, we report the HFLAV combined limit and the experimental published limits.In some cases, the combined limit is weaker than the limit published by a single experiment.This arises since the CLs method used in the combination can be more conservative compared to other legitimate methods, especially when the number of observed events fluctuates below the expected background[182].

Table 2 .
Low-energy QCD⊗QED invariant EFT for µ → e, µ → 3e and µ → e conversion in nuclei, in the notation of