Reach and complementarity of µ → e searches

In Effective Field Theory, we describe µ ↔ e flavour changing transitions using an operator basis motivated by experimental observables. In a six-dimensional subspace probed by µ → eγ , µ → e ¯ ee and µ → e conversion on nuclei, we derive constraints on the New Physics scale from past and future experiments, illustrating the complementarity of the processes in an intuitive way


I. INTRODUCTION
Lepton flavor-violating contact interactions -referred to as LFV or CLFV -are excellent probes of physics beyond the standard model (see e.g.[1,2] for reviews).The non-zero neutrino masses and mixing angles imply their existence, and an observation could shed light on the neutrino mass mechanism [3], and even on the matter excess of our universe if generated via leptogenesis [4,5].The rates could be just below the current experimental bounds, in many extensions of the Standard Model that introduce additional CLFV sources [2].
CLFV searches have been conducted in a wide range of reactions, and a subset of current and anticipated experimental constraints are given in Table I.While the multitude of τ channels could contribute to identifying the nature of New Physics (NP), the greater sensitivity in muon sector seems more promising for discovering it.Although the number of muon processes is limited, the corresponding bounds are already quite restrictive, and exceptional improvements are expected in the coming years.
The reach and complementarity of µ → eγ, µ → eēe and µA → eA transitions have been explored from various perspectives.Numerous authors have investigated model preferences and predictions for correlations among CLFV observables [18].An early model-independent analysis was performed by de Gouvea and Vogel [19], using a simple effective Lagrangian to describe NP effects.The presentation of our results follows their intuitive plots.However, their approach was limited to comparing the reach of pairs of processes (e.g.µ → eγ vs. µA → eA) at tree level.More systematic Effective Field Theory (EFT) studies, including more operators and some loop effects, were performed in [20,21]; [20] focused on the sensitivity of the processes to a more complete operator basis, and [21] explored whether the proposed experimental muon program is necessary and sufficient to find µ → e flavour change.
The aim of this work is to graphically illustrate the complementarity and reach of the µ → eγ, µ → eēe and µA → eA processes.We describe the physics of CLFV in an EFT perspective [22][23][24], where the number of operator coefficients is reduced by choosing an operator basis motivated by our observables [21].We quantify complementarity as the degree to which observables probe different operator coefficients, and study the complementarity of observables at the New Physics scale because the aim is to make observations that give distinct information at Λ LF V .The coefficients are translated to Λ LF V using Renormalisation Group Equations.This study extends the analysis described in [21] in several ways: we provide more informative plots of the current and projected experimental reaches, and a more rigorous construction of the basis for the subspace of experimentally accessible operator coefficients.In addition we draw attention to the information loss in matching nucleons to quarks using current theoretical results.Using this formalism to study whether µ → e processes can distinguish among models is an interesting question that we leave for a subsequent publication.
This paper is organized as follows: section II outlines the procedure to take the data parametrized in EFT from the experimental scale to beyond the weak scale.Section III presents constraints from various experimental measurements and projections for future initiatives.The construction of the basis used in this work is discussed in Appendices A and C.An independent issue regarding information loss in relating µA → eA rates to models is finally discussed in Appendix B.

II. THEORY OVERVIEW
This section gives the Lagrangian and Branching Ratios at the experimental scale, and sketches the transformation of operator coefficients from the experimental scale to Λ LF V (which is described in more detail in [21]).
A challenge of the EFT approach lies in the large number of operators.In the case of µ → e flavour changing processes, about 90 operators [21] are required to parametrize interactions that have ≤ 4 Standard Model legs at low energy and are otherwise flavour-diagonal.The difficulty to constrain and visualize this high-dimensional space is compounded by the fact that there are (only) three processes with excellent sensitivity in the µ → e sector (see Table I), imposing only about a dozen constraints on operator coefficients [25].Improved theoretical calculations and additional µA → eA measurements with different nuclear targets could increase this number to ∼ 20 independent constraints [25].Determining all EFT coefficients appears therefore a daunting task.
This manuscript takes a different perspective, following [21].Since there are three processes with excellent experimental sensitivity, we restrict to the (12-dimensional) subspace of operator coefficients probed by µ → eγ, µ → eēe, and Spin Independent1 µAl → eAl and µAu → eAu.The dimension of the subspace can be further reduced by half since the operator coefficients can be labelled by the helicity (or chirality) of the outgoing relativistic electron, and the results are very similar for either e L or e R , which do not interfere.Restricting the analysis to the subspace corresponding to an outgoing e L , the three muon processes can be described at the experimental scale (∼ m µ ) by the following effective Lagrangian [1]: where Λ LF V is a heavy mass scale, and the dimensionless coefficients {C Z } are lined up in a vector ⃗ C. The first term of this Lagrangian is a dipole operator mediating µ → eγ and contributing to both µ → eēe and µA → eA.The next three contact operators contribute to µ → eēe, while O Alight is a combination of operators probed by light muon conversion targets such as Ti or Al, and O Aheavy⊥ is an orthogonal combination probed by heavy targets such as Au 2 .We take Au and Al as prototypical "heavy" and "light" targets since Au was used by the SINDRUM experiment [10], and Al will be used by the upcoming COMET [11] and Mu2e [12] experiments, in addition to resembling Ti used in the past [10].
The constraint on the dipole operator from µ → eγ is given by: where we introduced unit vectors êA which select coefficients C A in the six-dimensional subspace.The four-lepton operators have negligeable interference in µ → eēe since the electrons are relativistic (≈ chiral), setting the three following constraints: These can conveniently be summarised as where the matrix R µ→e L ēe is proportional to the inverse covariance matrix for µ → eēe, given in eqn (A.3).Similarly, the Conversion Ratios for µA → eA can be written where the R matrices are given in eqn (A.9).These expressions justify a posteriori the basis in eqn (II.1), chosen to be orthogonal, intuitive, and correspond closely to the coefficient combinations probed by observables.The eigenvectors of the covariance matrix could be another basis choice, discussed briefly in Appendix C 3.
From the Lagrangian (II.1), one can easily deduce that the three processes are complementary at the experimental scale: four-fermion interactions with leptons only contribute to µ → eēe, interactions with strongly-interacting particles only contribute to µA → eA (the complementarity between heavy or light targets is discussed at the end of Appendix A), while the dipole contributes to all processes.The complementarity of two processes can be interpreted geometrically as the misalignment between the corresponding vectors in coefficient space, defined as the angle η between the two vectors.In terms of R matrices, the complementarity between processes A and B can be expressed as: which vanishes for perfectly complementary observables, and is equal to one when they contain the same information.The basis in eqn (II.1) was chosen to be perfectly complementary, i.e. orthogonal, at the experimental scale.
In the following, we adopt a different approach to illustrate the complementarity between processes.Instead of using the geometric measure defined above, we show that each process gives independent information about the operator coefficients by plotting the corresponding reach separately.The measured rates can then be combined to identify a point in parameter space.
The degree of complementarity can be evaluated at Λ LF V by translating the coefficients in eqn (II.1) from the experimental scale to Λ LF V .Modifying the scale amounts to changing the separation between lower energy loop effects that are explicitly calculated, and higher energy loops that are implicitly resummed into the coupling constants.At the experimental scale, all the loops via which a New Physics model contributes to an observable are in the operator coefficients, and the rate is straightforward to calculate.On the other hand, the operator coefficients at the heavy LFV scale are easily derived from a New Physics model, but loops must be calculated to evaluate experimental quantities.So the low-energy operator coefficients of eqn ( II.1) can be transformed to the LFV scale Λ LF V via the Renormalisation Group Equations (RGEs) [20,26], which peel off the SM loops in a leading log expansion.Solving the RGEs perturbatively, and modifying the EFT with scale to account for the changing particle content, allows us to write The matrix G is similar to that given in [21].We neglect loop effects in the EFT of nucleons and pions and match at 2 GeV onto a QCD-invariant EFT with gluons and five flavours of quarks (see Appendix A and Table IV).The leading log QED and QCD effects are included up to m W [20,29], where the coefficients are matched (at tree level in the lower-energy EFT) onto dimension six SMEFT operators, augmented by the dimension eight scalar operator [30,31] corresponding to êS , which could be relevant [32].We neglected CKM angles in matching and some other relevant SMEFT operators of dimension eight, and stress that the running from m W → Λ LF V is not included.The Branching Ratios in terms of coefficients at Λ LF V can be expressed as: where the matrix G is not unitary, and does not preserve the orthonormality of the basis since SM loops and matching can change the normalisation and direction of the vectors {ê A }.This is expected since distinct observations at low energy can measure the same high-scale NP coefficients.The changing modulus of the basis vectors is simple to calculate and include, and affects the reach.The changes in direction can affect the complementarity of processes if the vectors become more or less aligned (see Appendix C).

III. ILLUSTRATING EXPERIMENTAL CONSTRAINTS
In this section, we illustrate the constraints on New Physics from current and future µ → e searches, and show how these results can be combined to identify the allowed region of coefficient space.We parametrize the coefficient space with spherical coordinates [33] (Table II) assuming that the vector of coefficients ⃗ C is normalised to unity at the experimental scale.The reach of the various experiments in Λ LF V can be calculated as a function of these angles and the branching ratios given in eqn C.3.We stress that we are showing (projected) exclusion curves, as opposed to "one-at-a-time" bounds, since our EFT formulation should account for potential cancellations in the theoretical rate.
In deriving this parametrization, we approximating the operator coefficients as real numbers.This familiar simplification reduces our coefficient space from six complex to six real dimensions, replacing relative phases between interfering coefficients with a relative sign.Furthermore, we focus on a four-dimensional subspace, corresponding approximately to the four processes we examine, by suppressing two of the three four-lepton directions (the four-lepton operators can be distinguished by measuring the angular distribution in µ → eēe [34]).The direction ⃗ e S associated to the scalar four lepton operator interferes with none of the other operators and receives negligible loop corrections, so it is complementary by inspection.We also neglect a linear combination of the vector four-lepton directions ⃗ e V R and ⃗ e V L , since their contributions to µ → eēe have similar form.A judicious choice ensures the approximate orthogonality of the remaining four basis vectors.The full details are given in Appendix C. Modulo these approximations, the parametrisation describes the experimentally constrainable space, so we now plot various slices through the excluded region to illustrate its shape.II.Dimensionless operator coefficients expressed in the angular coordinates.The radial coordinate is 1/Λ 2 LF V , θI : 0..π and ϕ : 0..2π.As discussed in Appendix C 1, the ⃗ eV L × ⃗ eV R plane was projected to a line, deviations from which are measured by θV .In general, the basis vectors {eA} are not unit vectors, and their normalisation is given in Table V and after eqn (C.3) for the primed vectors.• ⃗ e D is negative and µA → eA vanishes when the dipole contribution cancels the remaining contributions.The rate drops abruptly, indicating that the dipole contribution is relatively small and the cancellation only occurs in a narrow region.The valley is broader for µ → eēe, since the contribution of ⃗ C • ⃗ e D is more important, and the rate never vanishes because µ → eēe independently constrains each coefficient contributing to this process, so the rate only vanishes when all the coefficients do (see eqn II.3).

We plot in
Our angular coordinate parametrisation defined a measure on the parameter space that assumes all the coefficients in our subspace are O(1) once the scale Λ LF V is fixed.This might not be the case in some classes of models; for instance ⃗ C • ⃗ e D ≫ four-fermion coefficients can occur (in SUSY [35]), or the dipole could be suppressed, when the four-fermion operators are generated at tree level.To illustrate complementarity when the "natural" size of ⃗ C • ⃗ e D is significantly different from the other coefficients, we also plot the reach in a parametrization similar to that introduced in [19] by defining a variable This non-linear transformation magnifies the regions where the dipole contribution either dominates the four-fermion interactions (θ = 0, π) or is suppressed (θ = π/2).(A similar variable κ V = cotan(θ V ) could be defined to magnify the regions where leptonic four-fermion coefficients are much larger or smaller than those with quarks.)We subtract π/2 in eqn (III.1) in order to have µ → eγ larger at the centre of the plot3 following [19].II.The solid region is currently excluded.
Figure 2 displays the reach as a function of θ V , which is effectively the angle between the µ → eēe and µA → eA four-fermion operators.Results for a vanishing dipole contribution (θ D = π/2) shows that µ → eēe vanishes at θ V = π/2 and µA → eA at θ V = 0, π.Adding a small negative dipole coefficient, µ → eēe doesn't vanish anymore since the dipole contributes independently as well as in interference with the four-fermion contributions, and the rate is reduced when this interference is destructive.The magnitude of the negative dipole coefficient is larger for θ D = 3π/4, exhibiting that µA → eA vanishes when the dipole cancels the four-fermion contributions.
Figure 3 illustrates the complementarity of heavy and light targets for µA → eA, by plotting the conversion ratios as function of ⃗ C • ⃗ e Alight ∝ sin ϕ and ⃗ C • ⃗ e Aheavy⊥ ∝ cos ϕ.Recall that ⃗ C • ⃗ e Aheavy⊥ parametrizes the independent information obtained with Au.This additional contribution to µAu → e L Au causes the rate to vanish at a different value than that of the light targets.The dipole, which also contributes to µA → eA, was taken to either vanish (θ D = π/2), be positive (θ D = 3π/4) or negative (θ D = π/4).This illustrates the impact of ⃗ C • ⃗ e D on the rate: cancellations can occur among the dipole and four-fermion contributions, as well as between the two independent combinations of four-fermion coefficients.
Finally, the dependence of the sensitivity on the angle ϕ and the variable κ D is exhibited in Figure 4.As expected, the µ → eγ and µ → eēe processes are independent of ϕ.The shape of the conversion processes on light and heavy targets are globally similar, although the ridges along which the rates cancel are slightly different.

IV. SUMMARY
We use bottom-up EFT to calculate the reach and illustrate the complementarity of experiments searching for NP.This method is particularly well-suited to situations in which the number of observables is much smaller than the number of operators.It provides a complete parametrisation of the rates, without redundancies, and the EFT translation to Λ LF V can be systematically improved.In addition, this formalism allows to explore the complementarity in a self-consistent manner at the same scale at which the theory is defined, and ensure that experiments effectively probe different combinations of NP parameters.This approach is generic and can be applied to many situations.In this manuscript, we use it to study CLFV in the muon sector and derive sensitivity projections for current and future experiments.
At the experimental scale, the Lagrangian given in eqn (II.1) includes all and only the operators contributing at tree level to the observables.The combinations of coefficients constrained experimentally define the operator basis for our subspace, whose dimension is equal to the number of constraints.For µ → e L γ, µ → e L ēe and Spin Independent µAl → e L Al and µAu → e L Au, this subspace is six-dimensional.These coefficients are translated to Λ N P by solving the leading order Renormalisation Group Equations below the weak scale, and matching them to SMEFT at tree level (see eqn (II.7)).Since the number of constraints remains unchanged, the dimension of the subspace cannot grow (but it could decrease, as discussed in appendix B).However, the normalisation and direction of the basis vectors is altered, in order to include, at Λ LF V , the contributions from all the operators to the observables via short-distance effects described in the RGEs.
The ability of different experiments to probe independent operator coefficients -our definition of complementarity -is related to the misalignement between vector of coefficients.While it can be measured in various ways, we observe that a judiciously selected subset of our basis vectors remain approximately orthogonal above the weak scale, and we use various parametrisations (see table II or eqn (III.1)) to plot the experimental exclusion curves using the Branching Ratios given in eqn (C.3).We also display a few projections to illustrate the reach and complementarity of future experiments.
An example of distinct observables probing the same New Physics is recalled in appendix B: µA → eA on various nuclei could distinguish scalar µ → e contact interactions on neutrons from protons, but this may not allow the distinction of LFV scalar operators involving up quarks from those with down quarks.Improving the precision of the scalar qq expectation values in the nucleon would be required to improve the situation.
This work is only a preliminary implementation of bottom-up EFT, relying on theoretical formalism described in [21].In future work, we aim to implement the Renormalisation Group running of our vectors above the weak scale (it was neglected here for simplicity and the lack of knowledge of Λ LF V ), and match models onto the "observable subspace" at Λ LF V .We hope that finding robust distinctions among model predictions could be simplified by the reduced dimension of the subspace.SD thanks A Ibarra for a seminar invitation and discussions.

Appendix A: Operators, Rates and R-matrices at low energy
The operators contributing at the experimental scale ∼ m µ to the tree level amplitudes for µ → e L γ, µ → e L ēe and Spin Independent µA → e L A are given in eqn (II.1),where the µA → eA operators can be expanded on: where N ∈ {p, n}, and the list neglects subdominant operators such as eP X µF F [54].
The Branching Ratios (BRs) for the various processes can be expressed in terms of the operator coefficients ⃗ C and the matrices R: where v ≃ m t is the Higgs vev, and ⃗ C • êD = C D .For µ → e L ēe, the BRs given in the text can be written using: ).For µA → eA, we consider the conversion ratios on prototypical heavy and light targets, taken to be Au and Al 4 .Following Kitano, Koike and Okada (KKO) [42], the Spin-Independent conversion rate can be written where we use the nucleus(A)-and nucleon(N )-dependent "overlap integrals" A , D A given by KKO [42] 5 , and Γ cap is the rate for the muon to transform to a neutrino by capture on the nucleus.Some relevant rates are [55]: KKO observed that the overlap integrals were nucleus-dependent, and measurements of µ → e conversion on different targets could be used to determine the operator coefficients.Reference [25] explored this issue quantitatively, and showed that with current uncertainties, Ti and Au give independent constraints.In this work, we are interested in a slightly different question: whether the observables give different constraints on New Physics heavier than m µ and Λ QCD , instead of understanding if they are independent.Unfortunately, there is "information loss" in matching nucleons onto quarks (see section B), so we match the nucleon onto quark operators before constructing the R-matrices for µAl → e L Al and µAu → e L Au.
Nucleon operators can be matched at 2 GeV onto light quark operators (see [21] for a basis) as x .We also include the two-step matching of scalar b and c operators (ēP R µ)( QQ) [51,56], first onto the gluon operator (ēP R µ)GG, then onto nucleons.As a result, the nucleon and quark coefficients are related as where O ∈ {S, V } and the relevant Gs are given in Table IV.One can then define quark "overlap integrals" for target A as where in this work we use the EFT results [44] for G N,q S , which differ by ∼50% from lattice results [43,47].Assembling the quark overlap integrals for target A into a "target vector" ⃗ u A : allows to write the Conversion Ratio as 4 SINDRUM searched for µ → e conversion on Titanium [10], and we will use these results in our study.However, according to [25], Ti and Al probe the "same" operator coefficients within current uncertainties, so we will apply the Ti bounds along the direction in coefficient space corresponding to Al. 5 The "4/Λ 4 LF V " in eqn (A.4) differs from the "2G 2 F " given in KKO's eqn (14), because −4 C| here = g| KKO in the case of four-fermion operators; and D A is divided by 4 in the above because the dipole normalisation here is identical to KKO.In addition, the KKO overlap integrals are in units of m 3/2 µ , which here sits in front.
where D A is the KKO overlap integrals for the dipole [42].So the R-matrix at a scale of 2 GeV is This translation of nucleon to quark operators neglects higher order QED and "strong interaction" effects between the experimental scale m µ and 2 GeV: the QED running is small for the considered operators, and we did not include recent χPT results [40].
The prospects of distinguishing coefficients by using different targets depend on the misalignment between the target vectors, which can be quantified as an angle [25]: In the quark operator basis at 2 GeV, this angle is given for various potential targets in Table III.The angles are smaller than the misalignmeent angles in the nucleon operator basis [25] (see Figures 6 and 7), because the scalar quark overlap integrals are larger than the vector integrals, and comparable for u and d quarks.In addition, there is currently a large discrepancy, ∼ 50%, between lattice and EFT determinations of the scalar overlap integrals.We assume that this theoretical discrepancy can be solved, so that µA → e L A on Au and Al give independent information.Eqn (A.10) gives the R-matrices for Al and Au, which probe directions in quark coefficient space that are misaligned by ∼ 5 degrees, at a scale of 2 GeV.In the plane spanned by ⃗ u Au and ⃗ u Al , the orthogonal combinations used in eqn (II.1), can be obtained by writing the target vector ⃗ u Au for Au as: where θ AuAl is given in Table III, and ⃗ u ⊥ ≈ .56êuu V + .8êdd V is the direction in coefficient space corresponding to the operator O Aheavy⊥ of eqn (II.1).

Appendix B: Are scalar quark currents indistinguishable?
This Appendix is focused on the information loss that occurs in matching nucleons to quarks at 2 GeV.In the data, LFV scalar interactions with neutrons might be distinguishable from those on protons [25].But current theoretical translations of the nucleon results to quarks erase any distinction between LFV interactions with scalar u or d currents (the subdominant s, c, b, t quarks are neglected in this section).
The Spin-Independent Conversion Rate (CR) is given in eqn (A.4) in terms of operator coefficients on nucleons.This result is at "Leading Order" in the low energy theory, and does not include the next order in χPT (parametrically ∼ 10%, two-nucleon effects, pion exchange...) or in the nuclear matrix element.Such effects have been calculated for WIMP scattering [38], and some partial results for µA → eA have been obtained [39][40][41].Improvements of the theoretical calculation of Γ(µA → eA) could change the form of the conversion ratio (this occurs in [40]), or reduce the uncertainties on its parameters, thereby resolving the issue discussed here.
To set the stage, recall from KKO [42] that the overlap integrals for scalar and vector densities of neutrons and protons differ by less than a factor of three.In particular, Al approximately probes (neglecting the dipole which could be constrained/measured elsewhere).In addition, it was pointed out in [25] that measuring µA → eA on another light target with different numbers of n and p would allow the measurement of Then, as noted by KKO, vector overlap integrals dominate over the scalars in heavy targets, so comparing Al to Au could allow to determine However, heavy targets also have more neutrons, so the n − p measurement from comparing two light targets is required to extract the S − V from heavy targets.The one remaining combination of coefficients (an isospin-violating V -S difference) has little impact on the CR given in eqn (A.4), and cannot be extracted with current theoretical uncertainties [25].
TABLE IV.Expectation values of vector and scalar quark currents in the nucleon; the scalar values are sometimes called f N q [43], and for light quarks are related to σπN [44], where G N,q S = m N mq f N q = σπN /mq+ isospin corrections.The first scalar GS induced by u and d quarks were obtained via dispersive techniques [44][45][46] and EFT methods [49] (Some lattice calculations give similar results [50]).The second G N,q S [in brackets] are the lattice results of BMW [47], and an average of lattice results [48] is used for the strange quark.The heavy quark scalar GSs are from the lattice [43], and in parentheses from [51].In all cases, the MS quark masses at µ = 2 GeV are taken as mu = 2.2 MeV, m d = 4.7 MeV, and ms = 96 MeV [52].The running heavy quark masses are [53]  The operator coefficients on nucleons can be transformed to coefficients on quarks according to eqn (A.6), using the matrix G given in Table IV 6 .Focusing on the first generation valence quarks, this can be written as: where the vector coefficients exhibit the expected dominance of u quarks in the proton, and d quarks in the neutron: the quark vector coefficients can be calculated from the nucleon vector nucleon coefficients, and vice versa.In the case of the scalar coefficients, this is almost not the case; for both lattice and EFT determinations of the G N,q S , the determinant of the scalar submatrix in eqn (B.4) is small compared to the product of two )), causing large uncertainties when it is inverted to solve for quark coefficients as a function of nucleon coefficients.In addition, G N,q S ≫ G N,q V increases the sensitivity to scalar coefficients, and reduces the relative contribution of the vector coefficients to about the magnitude of the scalar uncertainties.
The effects of transforming from nucleon operators to quark operators can be seen in Figure 5, where are plotted the quark and nucleon overlap integrals.This illustrates the difficulty to distinguish scalar u vs d, and also that the variation with Z is reduced for the quarks compared to nucleons.
The ability of different targets to distinguish among operator coefficients can be quantified as the angle between the directions they probe in coefficient space [25], where the direction is given by the overlap integrals for the target.Figures 6 and 7 plots this misalignment angle as a function of Z for various past and future targets.The plot on the left is for quark coefficients, and on the right is for nucleon coefficients (notice the difference in the vertical scale).Reference [25] estimated that a misalignment angle > ∼ 0.1 was required to overcome theoretical uncertainties in the nuclear calculation, and obtain independent measurement of distinct nucleon operator coefficients.Even if one neglects the theoretical uncertainty in the G N,q O s, applying this > ∼ 0.1 rule for quark coefficients suggests that the theoretical accuracy needs to be improved.
In summary, upcoming experiments could hope to bound or measure the three combinations of nucleon coefficients given in eqns (B.1-B.3).However, due to the almost-vanishing scalar determinant in eqn (B.4), theoretical progress in calculating the G N q S is required, for these to constrain three independent combinations of quark operator coefficients.From the top-down at Z = 90: vector n (red, dotted), scalar d (blue) and u (green), vector p (black,dotted), scalar n (blue,dotted) and p (green, dotted), and vector d (red) and u (black).The nucleon overlap integrals are defined in Eq. (A.4), and those for quarks in Eq. (A.7); an overlap integral for nucleons (or quarks) in target A is normalised by dividing by the the square root of the sum of the squares of all the nucleon (or quark) integrals for that target.

Appendix C: Including the penguin
This Appendix discusses the orthogonality of the basis used for the experimentally accessible subspace.We need an orthogonal basis to illustrate complementarity in polar coordinates; orthogonality would not be required to evaluate complementarity via eqn (II.6), nor for exploring the model predictions for coefficients in the subspace.In the following, the term "penguins" refer to the operators of eqn (C.1).
The basis vectors at Λ LF V > ∼ m W are G * (m W , m µ )ê A ; we do not give explicit expressions because the translation from the observable-motivated basis to an arbitrary other choice is a technicality more suitable to computers.Expressions for various ⃗ C † G * (m W , m µ )ê A , appropriate for calculating LFV BRs in terms of coefficients at the weak scale, can be found in [21].The norms and inner products among some basis vectors are given in Table V, which shows an overlap of 10 → 20 degrees between the vector operators (e ⊥ has significant quark-vector components).The shrinking of the quark target-vectors is largely due to the shrinking scalar coefficients.
The overlap among the coefficients of (ēγ α P L µ)( f γ α f ) operators arises because the µ → e penguin operators of the SM EFT [57]: (C HL1 + C eµ HL3 )ē Z / P L µ, so that Z exchange gives four-fermion operators in the low energy EFT: where f is any light chiral fermion, we used the SM Z interaction g 2c W (g f L P L +g f R P R ) and c W = cos θ W .As a result, at the weak scale, the basis vectors ⃗ e V L and ⃗ e V R have components g e L and g e R in the two penguin directions, respectively, giving a tree-level contribution to ⃗ e V L • ⃗ e V R of order 2g e L g e R = 4s 2 W (−1 + 2s 2 W ). A more orthogonal basis could be obtained by removing the penguin contribution from the low energy (ēγ α P L µ)( f γ α f ) operators, and adding Z → µe L as an observable at m W .This is analogous to what was already done for the dipole, removing it from the combination of operators contributing to µA → eA, and including µ → eγ as an observable.However, this adds a dimension to the subspace, and tangles the intuitive link between basis vectors and observables.It is pursued in section C 2. Section C 1 describes a simpler approach used to make the plots in the body of the paper.V. Inner products between the basis vectors ⃗ ei(ΛLF V ∼ mW ) = G * (mW , mµ)êi.On the diagonal are the norms (at mW in SMEFT, for unit-normalised êi at low energy), whereas the off-diagonals are cos η, where η is the angle between the vectors.Overlaps smaller than 10 −3 are given as 0.
1.A reduced basis at ΛLF V ∼ mW In this section we outline a method for choosing a penguin-less basis vector corresponding to a linear combination of ⃗ e V L and ⃗ e V R , which is approximately orthogonal to the remaining basis vectors.In addition, the complementarity plots involving ⃗ e V L or ⃗ e V R are similar, so this choice suppresses redundancy.
The dependence of BR(µ → eēe) on C V L and C V R is very similar, as can be seen from eqn (II.3).Since these coefficients can be distinguished via the angular distributions of the final state electrons in µ → eēe [34] rather than by comparing BR(µ → eγ) and BR(µ → eēe), it is sufficient to plot µ → eγ and µ → eēe in terms of a combination of these coefficients.By choosing this combination to suppress penguin contribution, we obtain five approximately orthogonal basis vectors.
As mentioned above, in the SMEFT basis at m W , ⃗ e V L and ⃗ e V R have components ∼ g e L , g e R along the directions of the penguin operators of eqn (C.1).So introducing ⃗ e V ∝ g e R ⃗ e V L − g e L ⃗ e V R , one sees that its penguin component vanishes at tree level, and ⃗ e V will be orthogonal to ⃗ e Aheavy⊥ up to loop effects.We therefore replace the ⃗ e V R ⊗ ⃗ e V L plane by an axis along ⃗ e V and plot the complementarity of the three rates in the resulting 5-dimensional space.
In this reduced basis, the BRs can be written in terms of operator coefficients at where we replace {θ V R , θ V L } → θ V , the un-primed norms of the basis vectors at m W are on the diagonal of Table V and |⃗ e 2. An enlarged basis at ΛLF V ∼ mW This section outlines the approach of adding Z → e ± µ ∓ to the observables, and removing the contribution of the flavour-changing Z penguin from the four-fermion operators.This ensures that the basis vectors are orthogonal at m W to within a degree or two, and highlights the importance of Z → e ± µ ∓ for distinguishing among coefficients and models.
The operators O eµ HL1 and O eµ HL3 of eqn (C.1) mediate flavour-changing Z decays, upon which ATLAS [58] sets the constraint BR(Z → e ± µ ∓ ) < 7.5 × 10 −7 (based on 20 fb −1 of luminosity), implying: The current sensitivity of µA → eA and µ → eēe to these coefficients is three orders of magnitude better, and should improve by another two orders of magnitude with upcoming experiments.Nonetheless, improving the experimental reach in Z → e ± µ ∓ is interesting, because experiments at muon mass scale cannot distinguish these penguin operators from the four-fermion ones 7  We include an additional basis vector in the coefficient subspace above m W : The formulae for the reach, obtained from the Branching Ratios are: where the one-loop matching of penguin operators to the dipole was included.For µ → eēe:

The eigenbasis of the covariance matrix
An alternative basis for the subspace of constrained coefficients, also orthogonal and perhaps more familiar, would be the eigenvectors of the covariance matrix.The inverse covariance matrix V for all the processes can be written as where the coefficients are evaluated at the experimental scale.The inverse eigenvalues give the allowed range of coefficients in the eigenbasis, and the eigenvectors are orthogonal combinations of operators, which correspond to the axes of the allowed ellipse around the origin in coefficient space.This is a convenient basis for plotting, because the limiting values of each parameter are obtained on the axes, so there is no need to do perform scans.However, we prefer the operator basis of eqn (II.1), because it is simple and directly related to the experimental processes.
A covariance matrix for the coefficients at Λ LF V can be obtained by substituting eqn (II.7) into eqn (C.11).This matrix is large, (∼ 90 × 90), so despite that most of the eigenvalues should vanish, finding the eigenvectors of the 12 non-zero eigenvalues would be a numerical exercise which could disconnect the final basis and constraints from the input processes.It has the advantage of giving an orthonormal basis, whose eigenvectors correspond to the axes of the allowed ellipse in coefficient space.

FIG. 1 .
FIG.1.Reach as a function of (left) the angle θD, which parametrizes the relative magnitude of dipole and four-fermion coefficients, and (right) the variable κD = cotan(θD − π/2).The scale Λ is defined in eqn (II.1) with the coefficients normalised according to TableII.The solid region is currently excluded.

7 FIG. 5 .
FIG.5.Normalised quark (solid line) and nucleon (dotted line) overlap integrals as a function of the target charge Z. From the top-down at Z = 90: vector n (red, dotted), scalar d (blue) and u (green), vector p (black,dotted), scalar n (blue,dotted) and p (green, dotted), and vector d (red) and u (black).The nucleon overlap integrals are defined in Eq. (A.4), and those for quarks in Eq. (A.7); an overlap integral for nucleons (or quarks) in target A is normalised by dividing by the the square root of the sum of the squares of all the nucleon (or quark) integrals for that target.

FIG. 7 .
FIG.7.The misalignement angle θ (defined in Eq. (A.11) between quark target vectors; notice the reduced vertical scale and smoother lines as compared to the nucleon plot.The decreasing lines are lead (green,upper) and Au (blue); the rising lines are titanium (blue) and Al (red,dashed, middle) and lithium (black,dashed, upper).
case of an outgoing e L ; for outgoing e R it would be êHE ), and rewrite the operator coefficients in spherical coordinates as θ ping êping + sin θ ping cos θ D êD + sin θ ping sin θ D cos θ S êS,R + sinθ ping sin θ D sin θ S cos θ V R êV,R + sin θ ping sin θ D sin θ S sin θ V R cos θ V L êV,L + sin θ ping sin θ D sin θ S sin θ V R sin θ V L (sin ϕê Al + cos ϕê ⊥ ) , (C.6)where θ J : 0..π and ϕ : 0..2π, and recall that a model would predict the various angles and the scale.The expressions for low energy BRs in this enlarged basis become more complicated, because the low energy fourfermion coefficients are expressed as the component from Z exchange, plus the component from four-fermion operator at the weak scale.The normalisation of some basis vectors changes as well, becoming: |e V L | = 0.944 , |e V R | = 1.041 , |e Al | = 0.768 , |e Aheavy⊥ | = 0.861 (C.7)

TABLE III
. Misalignment angles (in degrees), between target vectors expressed in the quark operator basis.