Completeness and complementarity for $\mu \rightarrow e\gamma\mu \rightarrow e\bar{e}e$
and $\mu A \rightarrow eA$

S. Davidson

<figure><img src="image_1.png" /></figure>

To cite this version:

S. Davidson. Completeness and complementarity for $\mu \rightarrow e\gamma\mu \rightarrow e\bar{e}e$ and $\mu A \rightarrow eA$. Journal of High Energy Physics, 2021, 02, pp.172. 10.1007/JHEP02(2021)172. hal-02973165

HAL Id: hal-02973165

https://hal.science/hal-02973165

Submitted on 22 Nov 2021

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Completeness and complementarity for $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu A \to eA$

## S. Davidson

LUPM, CNRS, Université Montpellier,
Place Eugene Bataillon, F-34095 Montpellier, Cedex 5, France

E-mail: s.davidson@lupm.in2p3.fr

**ABSTRACT:** Lepton Flavour Violation (LFV) is New Physics that must occur, but is strin-
gently constrained by experiments searching for $\mu \leftrightarrow e$ flavour change, such as $\mu \to e\gamma$,
$\mu \to e\bar{e}e$ or $\mu \to e$ conversion. However, in an Effective Field Theory(EFT) parametri-
sation, there are many more $\mu \leftrightarrow e$ operators than restrictive constraints, so determining
operator coefficients from data is a remote dream. It is nonetheless interesting to learn
about New Physics from data, so this manuscript introduces "observable-vectors" in the
space of operator coefficients, which identify at any scale the combination of coefficients
probed by the observable. These vectors have an overlap $\gtrsim 10^{-3}$ with most of the co-
efficients, and are used to study whether $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu \to e$ conversion give
complementary information about New Physics. The appendix gives updated sensitivities
of these processes, (and a subset of $\tau \to \ell$ decays), to operator coefficients at the weak
scale in the SMEFT and in the EFT below $m_W$.

**KEYWORDS:** Beyond Standard Model, Effective Field Theories

ARXIV ePRINT: 2010.00317
# Contents

<table><tr><td>1</td><td>Introduction</td><td>1</td></tr><tr><td>2</td><td>Completeness</td><td>3</td></tr><tr><td>2.1</td><td>Operators</td><td>3</td></tr><tr><td>2.2</td><td>Experimental bounds</td><td>5</td></tr><tr><td>2.2.1</td><td>Translating μ→e conversion bounds onto quark operators</td><td>6</td></tr><tr><td>2.3</td><td>Accounting for loop effects</td><td>7</td></tr><tr><td>2.4</td><td>Constraints at m<sub>W</sub></td><td>9</td></tr><tr><td>2.4.1</td><td>μ→eγ</td><td>9</td></tr><tr><td>2.4.2</td><td>μ→e<bar>e>e</td><td>11</td></tr><tr><td>2.4.3</td><td>μ→e conversion</td><td>11</td></tr><tr><td>2.5</td><td>What is not probed?</td><td>12</td></tr><tr><td>3</td><td>Complementarity</td><td>13</td></tr><tr><td>4</td><td>Discussion and summary</td><td>17</td></tr><tr><td>A</td><td>Tables of sensitivities</td><td>19</td></tr><tr><td>A.1</td><td>Including a selection of tau decays</td><td>21</td></tr></table>

## 1 Introduction

The observed pattern of neutrino masses implies new particles and interactions (New Physics $\equiv$ NP) in the lepton sector (beyond the Standard Model with left-handed neutrinos), so one of the best-motivated challenges in particle physics is to discover what it is. Knowing where to search would be a useful input to this enterprise. One possibility would to look everywhere where NP is not already excluded; another would be to explore regions suggested by motivated models. However, the guidance from models can be ambiguous, because they are legion, and for a given process, many models may predict rates just beyond the current experimental reach. So this manuscript takes the agnostic view that it is interesting to look for Lepton Flavour Violation (LFV) everywhere it is not excluded. The current aim is to identify where that is; once NP is detected, reconstructing it becomes interesting.

This manuscript focuses on NP that changes lepton flavour $\mu \leftrightarrow e$ (for a review, see e.g. [23]), for simplicity restricted to processes that are otherwise flavour diagonal, such as $\mu \to e\bar{e}e$ or $\mu \to e\gamma\gamma$, but not $K \to \mu^\pm e^\mp$. In order to look for LFV “everywhere that it is not excluded” among such processes, it would be interesting to list observables such that:

1. if $\mu \leftrightarrow e$ flavour-changing NP exists, it would contribute to at least one of the processes. That is, the set of observables is “complete”.
2. the observables are complementary, in the sense that they give independent information about the NP, and cannot be predicted one from the other. In particular, we want to avoid searching for a branching ratio that is already excluded by the upper bound on another observable.

Instead of trying to construct such a list, we ask a more pragmatic question: to what degree are $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu \to e$ conversion on nuclei ($\mu A \to eA$) complete and complementary? These processes are selected because the current experimental bounds on the branching ratios are restrictive ($\lesssim 10^{-12}$) [1–4], and experiments under construction aim to improve the sensitivities to $10^{-14} \to 10^{-16}$ [5–9] or better [10].

Ideally, a complete list of complementary processes should be independent of the theoretical formalism used to establish it — for instance, it should not apply only to some models, or depend on a choice of operator basis in an Effective Field Theory (EFT). Nonetheless, let us start by restricting to “heavy” NP models, where the new particle masses are at a scale $\Lambda_{NP} \gtrsim m_W$. Such models can be parametrised below $\Lambda_{NP}$ in an EFT framework, which allows to separate “known physics” (the data and the Standard Model), from the theoretical speculations above $\Lambda_{NP}$. So LFV is described by operators constructed out of Standard Model(SM) fields and respecting SM gauge symmetries.<fnref n="1" /> In addition, $\Lambda_{NP}$ is assumed large enough to justify retaining only a few terms in the $1/\Lambda_{NP}^n$ expansion, which in the case of LFV, starts at dimension six $\propto 1/\Lambda_{NP}^2$. In this EFT context, a complete list of observables should contain at least as many members as there are operator coefficients — otherwise there can be combinations of coefficients that are not probed (sometimes called “flat directions”). Unfortunately, for $\mu \to e$ flavour change, there are many flat directions: current data on $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu \to e$ conversion impose 12 to 14 [11] bounds on the coefficients of the 80→90 operators listed in section 2.1.

A first step, is nonetheless to explore whether $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu \to e$ conversion are sensitive to the coefficient of each operator in the basis of section 2.1. Such a study is basis-dependent, and corresponds to calculating “one-operator-at-a-time” bounds, or “sensitivities”(the word used in this manuscript): a coefficient smaller than such values could not have been seen. This is consistent with the original aim of identifying where LFV is not excluded, and the results are given in the tables of appendix A. However, the coefficients can be larger than these sensitivities, by sitting along various flat directions, which are discussed in section 2.5.

Loop effects are the backbone of this discussion. This is because the contact interactions which are induced at tree level in models, may not mediate the processes which are stringently constrained by experiment. Consider, for instance, a NP model that induces a tree-level contact interaction $(\bar{b}\gamma^\alpha b)(\bar{e}\gamma_\alpha\mu)$. This mediates the decay $\Upsilon \to \bar{e}\mu$ [12], and contributes at one-loop to $\mu \to e\bar{e}e$. But $\mu \to e\bar{e}e$ will have better sensitivity, both because the bound on the BR is more restrictive, and because the muon lifetime is longer, since it decays via the weak interactions, whereas the $\Upsilon$ decays electromagnetically. SM loop corrections to the contact interactions, are discussed in more detail in section 2.3.

<footnote n="1">¹This overlooks the implications of possible of light NP in other sectors; for instance a light axion/ALP or DM particle.</footnote>
Constraints on LFV contact interactions have been tabulated in many previous publications. This manuscript differs from earlier works, such as [13–17], which mostly listed tree-level bounds, that is, constraints at the experimental scale on the coefficients of operators which contribute at tree level to the process. The sensitivities tabulated here are very similar to those of [18], which apply to coefficients at a short-distance scale ~ 1/$m_W$ and include similar Standard Model loop corrections between $m_W$ and low energy; the results here complete those of [18] (where some operators are missing from the tables), extend them by a subset of $\tau \leftrightarrow \ell$ operators, and also give the sensitivities in the SMEFT basis.

The “observable-vectors” are defined in section 3. When New Physics is parametrised via operators in EFT, these are constructed along theoretical guidelines, and identify the default basis of coefficient space. Observable-vectors could be an alternative basis, corresponding to the directions in coefficient space probed by observables. They appear implicitly in the results of section 2. Their misalignment quantifies the complementarity of the observables: if the vectors are orthogonal, the observables are very complementary, if the vectors are parallel, the observables probe the same thing, so section 3 uses this perspective to plot the complementarity of $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu \to e$ conversion.

## 2 Completeness

Ideally, a list of processes that is “complete” could include processes at any accessible energy scale, and would allow to probe any LFV contact interaction. However in practise, it is difficult to search for processes with more than four legs, and many restrictive constraints come from low-energy processes. So we attempt to construct a list of processes that is sensitive to all LFV contact interactions at “low” energies ~ $m_W$ that have three or four legs. Section 2.1 lists a QCD×QED-invariant basis of operators constructed with three or four lepton, quark photon or gluon fields, that change lepton flavour $\mu \to e$, and are otherwise flavour-diagonal. A quark and gluon basis is used because they appear in loops between $m_W$ and ~ 2 GeV. However, they live inside hadrons; the additional step of matching quark operators onto nucleon operators is discussed with $\mu \to e$ conversion in section 2.2. Section 2.2 gives the Branching Ratios for restrictively-bounded processes in terms of the operator coefficients at the experimental scale, section 2.3 describes how SM loop corrections are included in this analysis, and the branching ratios are given in terms of weak-scale coefficients in section 2.4. Allowing one coefficient at a time to be non-zero in these formulae, gives the “one-operator-at-a-time” bounds, or sensitivities, of each process to each operator, which are collected in the tables in appendix A. The operators can be matched onto the SMEFT at $m_W$, sensitivities are also given in this basis. Finally, since there are more operators than experimental constraints, section 2.5 discusses combinations of coefficients that are not constrained.

### 2.1 Operators

There are ninety operators which are QCD×QED invariant, have three or four legs, and which change flavour $\mu \to e$ (and involve no other flavour change). These are suitable for describing $\mu \to e$ interactions at energies below $m_W$, where the Higgs and SU(2) bosons are not present as external legs. The operators here are constructed with chiral fermions,
because this facilitates matching onto the chiral SMEFT operators at the weak scale, and because in the lepton sector, the electron from muon decays is relativistic, so ≈ chiral, implying that negligible interference between operators involving $e_L$ vs $e_R$. The operators are added to the Lagrangian as

$$ \mathcal{L} = \mathcal{L}_{\text{SM}} + \sum_{\zeta} \sum_{\text{Lor}} \frac{C_{\text{Lor}}^{\zeta}}{v^2} \mathcal{O}_{\text{Lor}}^{\zeta} + \text{h.c.} \quad (2.1) $$

where $1/v^2 = 2\sqrt{2}G_F$ ($v \simeq m_t$), the operator subscript Lor gives the Lorentz structure and chirality of the fermion bilinears, and the superscript $\zeta$ gives the flavour indices. Since this manuscript only considers $\mu \to e$ transitions, all the operators contain $\bar{e}$ and $\mu$ (the $\mu^+ \to e^+$ processes are described by the $+\text{h.c.}$), and the $e\mu$ indices are suppressed from the superscript.

The 22 four-lepton operators are:

$$ \begin{aligned} \mathcal{O}_{V,YY}^{ll} &= (\bar{e}\gamma^\alpha P_Y \mu)(\bar{l}\gamma_\alpha P_Y l), & \mathcal{O}_{V,YX}^{ll} &= (\bar{e}\gamma^\alpha P_Y \mu)(\bar{l}\gamma_\alpha P_X l) \\ \mathcal{O}_{S,YY}^{ll} &= (\bar{e}P_Y \mu)(\bar{l}P_Y l), & \mathcal{O}_{S,YX}^{\tau\tau} &= (\bar{e}P_Y \mu)(\bar{\tau}P_X \tau) \\ \mathcal{O}_{T,YY}^{\tau\tau} &= (\bar{e}\sigma P_Y \mu)(\bar{\tau}\sigma P_Y \tau) & & \end{aligned} \quad (2.2) $$

where $l \in \{e, \mu, \tau\}$, $X, Y \in \{L, R\}$, and $X \neq Y$. Then there are 50 operators with two leptons and a quark bilinear:

$$ \begin{aligned} \mathcal{O}_{V,YY}^{qq} &= (\bar{e}\gamma^\alpha P_Y \mu)(\bar{q}\gamma_\alpha P_Y q), & \mathcal{O}_{V,YX}^{qq} &= (\bar{e}\gamma^\alpha P_Y \mu)(\bar{q}\gamma_\alpha P_X q) \\ \mathcal{O}_{S,YY}^{qq} &= (\bar{e}P_Y \mu)(\bar{q}P_Y q), & \mathcal{O}_{S,YX}^{qq} &= (\bar{e}P_Y \mu)(\bar{q}P_X q) \\ \mathcal{O}_{T,YY}^{qq} &= (\bar{e}\sigma P_Y \mu)(\bar{q}\sigma P_Y q) & & \end{aligned} \quad (2.3) $$

where $q \in \{u, d, s, c, b\}$. And finally, there are 18 operators with two leptons, which include the dipoles and operators with two photons or gluons

$$ \begin{aligned} \mathcal{O}_{D,L} &= m_\mu \overline{e_R} \sigma^{\alpha\beta} \mu_L F_{\alpha\beta} & m_\mu \overline{e_L} \sigma^{\alpha\beta} \mu_R F_{\alpha\beta} \\ \mathcal{O}_{GG,Y} &= \frac{1}{v} (\bar{e} P_Y \mu) G_{\alpha\beta} G^{\alpha\beta}, & \mathcal{O}_{G\tilde{G},Y} &= \frac{1}{v} (\bar{e} P_Y \mu) G_{\alpha\beta} \tilde{G}^{\alpha\beta} \\ \mathcal{O}_{GGV,Y} &= \frac{1}{v^2} (\bar{e} \gamma^\sigma P_Y \mu) G_{\alpha\beta} \partial_\beta G^{\alpha\sigma}, & \mathcal{O}_{G\tilde{G}V,Y} &= \frac{1}{v^2} (\bar{e} \gamma^\sigma P_Y \mu) G_{\alpha\beta} \partial_\beta \tilde{G}^{\alpha\sigma} \\ \mathcal{O}_{FF,Y} &= \frac{1}{v} (\bar{e} P_Y \mu) F_{\alpha\beta} F^{\alpha\beta}, & \mathcal{O}_{F\tilde{F},Y} &= \frac{1}{v} (\bar{e} P_Y \mu) F_{\alpha\beta} \tilde{F}^{\alpha\beta} \\ \mathcal{O}_{FFV,Y} &= \frac{1}{v} (\bar{e} \gamma^\sigma P_Y \mu) F^{\alpha\beta} \partial_\beta F_{\alpha\sigma}, & \mathcal{O}_{F\tilde{F}V,Y} &= \frac{1}{v} (\bar{e} \gamma^\sigma P_Y \mu) F^{\alpha\beta} \partial_\beta \tilde{F}_{\alpha\sigma} \end{aligned} \quad (2.4) $$

where $\sigma^{\alpha\beta} = \frac{i}{2}[\gamma^\alpha, \gamma^b]$, $X, Y \in \{L, R\}$, and $X \neq Y$. The running of the muon mass in the dipole operators is neglected here. The dimension seven operators $O_{GG,L}$, $O_{GG,R}$ were included in $\mu \to e$ conversion in [19], and the other gluon operators will not be further considered here. The $\mu \to e\gamma\gamma$ rate due to the various two-photon operators was calculated in [20], and the Crystal Box experiment [21] set the constraint $BR(\mu \to e\gamma\gamma) \le 7.2 \times 10^{-11}$. The $O_{FF,Y}$ operators also contribute to $\mu \to e$ conversion [22], and the SINDRUMII search for $\mu Au \to eAu$ currently has the best sensitivity to $C_{FF,L}$ and $C_{FF,R}$.
## 2.2 Experimental bounds

We now want to relate the coefficients of these operators to experimental decay rates. We restrict to the bounds on $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu A \to eA$ because the current experimental upper bounds $BR \lesssim 10^{-12}$ [1–4] are restrictive, and will improve by orders of magnitude in coming years [5–9]. Furthermore, the branching ratios compare to a weak decay, so $BR \lesssim 10^{-12}$ probes a new physics scale $\Lambda_{NP} \lesssim 100$ TeV.

The Branching Ratio for $\mu \to e\gamma$ [23], and the current experimental bound [1] are

$$BR(\mu \to e\gamma) = 384\pi^2(|C_{DL}|^2 + |C_{DR}|^2) < 4.2 \times 10^{-13} \quad (2.5)$$

so the dipole coefficients at the experimental scale should be inside the circle in coefficient space given by eq. (2.5), that is, should separately satisfy $C_{D,X} \le 1.05 \times 10^{-8}$.

The branching ratio for $\mu \to e\bar{e}e$ is [23–25]

$$\begin{aligned} BR(\mu \to e\bar{e}e) & \le 10^{-12} \\ &= 2|C_{V,LL} + 4eC_{D,R}|^2 + 2|C_{V,RR} + 4eC_{D,L}|^2 \\ & \quad + |C_{V,RL} + 4eC_{D,L}|^2 + |C_{V,LR} + 4eC_{D,R}|^2 \\ & \quad + \frac{|C_{S,LL}|^2 + |C_{S,RR}|^2}{8} + \left(64 \ln \frac{m_\mu}{m_e} - 136\right) (|eC_{D,R}|^2 + |eC_{D,L}|^2) \end{aligned} \quad (2.6)$$

where $\ln \frac{m_\mu}{m_e} = 5.35$, so $e^2(64 \ln \frac{m_\mu}{m_e} - 136) \simeq 204.8e^2 \simeq 18.78$. Measuring the polarisation of the muon and the angular distribution of the electrons [24, 25], (and even the polarisation of the electrons), could allow to discriminate among these various contributions.

Combining the $\mu \to e\bar{e}e$ and $\mu \to e\gamma$ bounds in a covariance matrix allows to obtain separate constraints on the dipole and vector coefficients (see [11]). Since the current bound on $BR(\mu \to e\gamma)$ is restrictive, this amounts to imposing the bound (2.5) on the dipole coefficients, and then neglecting them in (2.6):

$$\begin{aligned} |C_{V,XX}^{ee}| & \le 7.0 \times 10^{-7}, & |C_{V,XY}^{ee}| & \le 10^{-6} \\ |C_{S,XX}^{ee}| & \le 2.8 \times 10^{-6} \end{aligned} \quad (2.7)$$

where the coefficients are evaluated at the experimental scale.

The conversion of a muon to an electron in nuclei is a sensitive probe of $\mu \to e$ flavour change in the presence of quarks. The $\mu^-$ is captured into the 1s state of the nucleus, and can then convert to an electron by interacting with the nucleons or electric field of the nucleus. The SINDRUMII experiment at PSI, with a continuous muon beam, searched for $\mu \to e$ conversion on Titanium and Gold [3, 4], setting bounds $BR(\mu A \to e + A) \lesssim 10^{-12}$. The theoretical rates for (Spin Independent) conversion on many targets are given in [26], and can be written [11]:

$$BR_{SI}(\mu A \to eA) = \frac{\Gamma(\mu A \to eA)}{\Gamma_{cap}(A)} = B_A [|\hat{v}_A \cdot \vec{C}_L|^2 + |\hat{v}_A \cdot \vec{C}_R|^2] \le \begin{cases} 4.3 \times 10^{-12} & \text{Ti} \\ 7 \times 10^{-13} & \text{Au} \end{cases}, \quad (2.8)$$

where the Branching Ratio is normalised to the capture rate $\mu A \to \nu_\mu A'$ on the same nucleus, and is expressed in terms of the coefficients $\{\tilde{C}\}$ of operators constructed with a
nucleons. Several comments:

1. The coefficient subscript gives the Lorentz structure ($V$ or $S$), then the chiral projector of the lepton current, but the nucleon current is not chiral because its more useful in the non-relativistic limit to use a scalar($S$), pseudoscalar($P$), vector ($V$), axial vector ($A$) and tensor ($T$) basis, where the $P, A, T$ components contribute to the Spin Dependent rate [27, 28] so are neglected here. The non-chiral coefficients can be written in terms of chiral coefficients as, eg

$$C_{V,Y}^{ff} = \frac{1}{2} (C_{V,YR}^{ff} + C_{V,YL}^{ff}), \quad C_{A,Y}^{ff} = \frac{1}{2} (C_{V,YR}^{ff} - C_{V,YL}^{ff}) \quad (2.9)$$

2. The formula after the last equality is given in the normalisation of [11], where the coefficients of nucleon operators have been assembled in vectors

$$\vec{C}_L = (\vec{C}_{D,R}, \vec{C}_{S,R}^{pp}, \vec{C}_{V,L}^{pp}, \vec{C}_{S,R}^{nn}, \vec{C}_{V,L}^{nn}) \quad (2.10)$$

(and similarly for $\vec{C}_R$), and the overlap integrals of Kitano, Koike and Okada [26] for target A have been assembled in unit-normalised “target vectors”

$$\vec{v}_A = \left( \frac{D_A}{4}, S_A^{(p)}, V_A^{(p)}, S_A^{(n)}, V_A^{(n)} \right)$$

whose normalisation is absorbed into the $B_A = \frac{32G_F^2 m_\mu^5 |\vec{v}_A|^2}{\Gamma_{cap}(A)}$.

The vectors and normalisation factors for Titanium and Gold are

$$\begin{aligned} \hat{v}_{Ti} &= (0.250, 0.426, 0.458, 0.503, 0.541), & B_{Ti} &= 250 \\ \hat{v}_{Au} &= (0.222, 0.289, 0.458, 0.432, 0.686), & B_{Au} &= 300. \end{aligned} \quad (2.11)$$

Since these vectors are misaligned, Titanium and Gold can measure independent combinations of coefficients [11].

### 2.2.1 Translating $\mu \to e$ conversion bounds onto quark operators

The nucleon coefficients can be written in terms of the quark coefficients (that we wish to constrain) using the expectation values of quark-currents in the nucleus $\{G_\Gamma^{N,q}\}$, defined as, eg $\langle N|\bar{q}(x)q(x)|N \rangle = G_S^{N,q} \bar{u}_N(P_f) u_N(P_i) e^{-i(P_f - P_i)x}$, and given in table 1. For the vector and scalar coefficients:

$$\begin{aligned} \tilde{C}_{V,Y}^{NN} &= \sum_{q \in u,d,s} G_V^{N,q} C_{S,Y}^{qq} \\ \tilde{C}_{S,Y}^{NN} &= \sum_{q \in u,d,s} G_S^{N,q} C_{S,Y}^{qq} + \sum_Q G_S^{N,Q} C_{S,Y}^{QQ} - \frac{8\pi m_N}{9\alpha_s m_t} C_{GG,Y} \end{aligned} \quad (2.12)$$

where $Q \in \{c, b\}$. For the scalars, the first term is the tree contribution of light quarks in the nucleon, the second is the one-loop contribution of heavy quarks to the gluon density (which contributes to the scalar nucleon current), so the explicit gluon contribution given in the last term only contains contributions from $m_W$ or above (such as the $(\bar{t}t)(\bar{e}\mu)$ operator).
<table><caption>**Table 1.** This table is taken from [28]. It gives matching coefficients between nucleon and light-quark-flavour-diagonal operators. The parenthese gives the uncertainty in the last figure(s). The $G_S^{N,q}$ were obtained via EFT methods [33, 34], and an average of lattice results [35] is used for the strange quark. The heavy quark scalar $G_S$ are from [36]. In all cases, the $\overline{\text{MS}}$ quark masses at $\mu = 2$ GeV are taken as $m_u = 2.2$ MeV, $m_d = 4.7$ MeV, and $m_s = 96$ MeV [39]. The nucleon masses are $m_p = 938$ MeV and $m_n = 939.6$ MeV.</caption><thead><tr><th>G<sub>V</sub><sup>p,u</sup> = G<sub>V</sub><sup>n,d</sup> = 2</th><th>G<sub>V</sub><sup>p,d</sup> = G<sub>V</sub><sup>n,u</sup> = 1</th><th>G<sub>V</sub><sup>p,s</sup> = G<sub>V</sub><sup>n,s</sup> = 0</th></tr></thead><tbody><tr><td>G<sub>S</sub><sup>p,u</sup> = <sup>m<sub>p</sub></sup>&frasl;<sub>m<sub>u</sub></sub>0.021(2) = 9.0</td><td>G<sub>S</sub><sup>p,d</sup> = <sup>m<sub>p</sub></sup>&frasl;<sub>m<sub>d</sub></sub>0.041(3) = 8.2</td><td>G<sub>S</sub><sup>p,s</sup> = <sup>m<sub>N</sub></sup>&frasl;<sub>m<sub>s</sub></sub>0.043(11) = 0.42</td></tr><tr><td>G<sub>S</sub><sup>n,u</sup> = <sup>m<sub>n</sub></sup>&frasl;<sub>m<sub>u</sub></sub>0.019(2) = 8.1</td><td>G<sub>S</sub><sup>n,d</sup> = <sup>m<sub>n</sub></sup>&frasl;<sub>m<sub>d</sub></sub>0.045(3) = 9.0</td><td>G<sub>S</sub><sup>n,s</sup> = <sup>m<sub>N</sub></sup>&frasl;<sub>m<sub>s</sub></sub>0.043(11) = 0.42</td></tr><tr><td>G<sub>S</sub><sup>N,c</sup> = <sup>2m<sub>N</sub></sup>&frasl;<sub>27m<sub>c</sub></sub></td><td>G<sub>S</sub><sup>N,b</sup> = <sup>2m<sub>N</sub></sup>&frasl;<sub>27m<sub>b</sub></sub></td><td></td></tr></tbody></table>

This allows to translate the upper bound on the Branching ratio on Gold $BR_{Au}$ to a bound on the quark-level coefficients at a scale $\sim 2$ GeV:

$$4.9 \times 10^{-8} = \sqrt{\frac{BR_{Au}^{exp}}{B_{Au}}} \\ \gtrsim |0.222C_{D,R} + 1.60C_{V,L}^{uu} + 1.83C_{V,L}^{dd} + 6.10C_{S,R}^{uu} + 6.258C_{S,R}^{dd} + 0.303C_{S,R}^{ss} + 0.721\frac{2m_N}{27m_c}C_{S,R}^{cc} + 0.721\frac{2m_N}{27m_b}C_{S,R}^{bb} - 0.721\frac{8\pi m_N}{9\alpha_s v}C_{GG,R}| \quad (2.13)$$

where $v \simeq m_t$, and $C_{S,R}^{bb}$ is at $m_b$ rather than 2 GeV. The same bound applies for $L \leftrightarrow R$.

The SINDRUMII upper bound on the Branching ratio on Titanium, (see eq. (2.8)), is slightly less sensitive to individual operator coefficients, but constrains a different linear combination. In the formulation of [11], a target nucleus of charge $Z$ can be viewed as unit vector $\hat{v}_Z$ in the space of operator coefficients, such that $BR \propto |\hat{v}_Z \cdot \vec{C}|^2$. This allows to write the direction probed by Titanium as the direction probed by Gold plus an orthogonal vector:

$$\hat{v}_{Ti} = \hat{v}_{Au} \cos \phi + \hat{v}_{\perp} \sin \phi$$

where $\cos \phi = \hat{v}_{Ti} \cdot \hat{v}_{Au}$. So one can solve for $\hat{v}_{\perp}$, and with a covariance matrix obtain that the SINDRUM bound on Titanium gives a constraint on $|\hat{v}_{\perp} \cdot \vec{C}|$:

$$\sqrt{\frac{BR_{Ti}^{exp}}{B_{Ti} \sin^2 \phi}} = 6.24 \times 10^{-7} \gtrsim |0.155C_{D,R} - 0.506C_{V,L}^{uu} - 1.168C_{V,L}^{dd} + 9.21C_{S,R}^{uu} + 9.02C_{S,R}^{dd} + 0.446C_{S,R}^{ss} + 1.062\frac{2m_N}{27m_c}C_{S,R}^{cc} + 1.062\frac{2m_N}{27m_b}C_{S,R}^{bb} - 1.062\frac{8\pi m_N}{9\alpha_s v}C_{GG,R}| \quad (2.14)$$

and the same bound applies for $L \leftrightarrow R$.

## 2.3 Accounting for loop effects

In a Wilsonian sense, the coefficients in the branching ratios of the previous section are evaluated at a scale $\Lambda_{\text{low}}$ near the experimental scale. In order to estimate sensitivities to
coefficients at $\Lambda_{NP}$, loop corrections on scales $\Lambda_{low} \to \Lambda_{NP}$ must be “peeled off”. These loops corrections are evaluated here in an EFT perspective, where only SM particles are dynamical below $\Lambda_{NP}$, and new particle propagators are expanded in $p^2/\Lambda_{NP}^2$ around a contact interaction.<fnref n="2" /> These SM loops are then independent of the New Physics, and only need to be calculated once as Renormalisation Group running of operator coefficients, and possibly as matching coefficients as discussed below. As usual with the Renormalisation Group, a one loop calculation generates the $(\log/16\pi^2)^n$ terms for all $n$, a 2 loop calculation generates the $(\log^{n-1}/(16\pi^2)^n)$ terms, and so on. Provided that $\Lambda_{NP}$ is sufficiently large, only a few terms in the $1/\Lambda_{NP}^2$ expansion are required.

The log-enhanced loops included here arise from the one-loop Renormalisation Group Equations (RGEs) of QED and QCD, which can be written in matrix form as

$$ \mu \frac{d}{d\mu} \vec{C} = \frac{\alpha_s(\mu)}{4\pi} \vec{C} \Gamma^s + \frac{\alpha_e}{4\pi} \vec{C} \Gamma, \qquad (2.15) $$

where the operator coefficients are lined up in the row vector $\vec{C}$. The diagonal anomalous dimension matrix $\Gamma^s$ of QCD only renormalises the scalar and tensor operators involving a quark bilinear. Since the QCD coupling $\alpha_s$ is large and runs significantly, its one-loop effects are resummed: $C(\mu_f) = (\frac{\alpha_s(\mu_f)}{\alpha_s(\mu_i)})^{\gamma/2\beta_0}C(\mu_i)$ [40]. Two-loop QCD effects are neglected — although this gives an uncertainty of $\mathcal{O}(10\%)$ — because QCD is less interesting than QED: QED mixes operators, allowing to transform an operator that is difficult to probe experimentally, into one that is tightly constrained. Finally, $\Gamma$ is the well-populated anomalous dimension matrix of QED, augmented by two-loop QED mixing<fnref n="3" /> of vector operators into the dipole [18].

The RGEs are solved analytically as an expansion in $\alpha_e$, which is approximated not to run; we retain the $\mathcal{O}(\alpha_e \log)$ terms, and some of the $\mathcal{O}(\alpha_e^2 \log^2)$ and $\mathcal{O}(\alpha_e^2 \log)$ terms. The aim is to include effects that are $\gtrsim 10^{-3}$. The solution is formally a scale-ordered exponential, which can be expanded in $\alpha_e$ by defining $dt = d\mu/\mu$, and substituting $\frac{\alpha_s}{4\pi}\Gamma^s = -iH_0$, $\frac{\alpha_e}{4\pi}\Gamma^T = -iH_{\text{int}}$, into the solution for the time-translation operator given in chapter 4 of [41]. At $\mathcal{O}(\alpha_e)$, this gives:

$$ \vec{C}(\mu_f) \approx \vec{C}(\mu_i) D_s \left[ I + \frac{\alpha_e}{4\pi} \tilde{\Gamma} \ln \frac{\mu_f}{\mu_i} + \dots \right] \\ \tilde{\Gamma}_{KJ} = \frac{\Gamma_{KJ}^e(\mu_f)}{1 + a_J - a_K - a_d} \frac{1 - \lambda^{a_J - a_K - a_d + 1}}{1 - \lambda} \equiv f_{KJ} \Gamma_{KJ}^e \qquad (2.16) $$

where there is no sum on $KJ$ in $f_{KJ}\Gamma_{KJ}^e$, $D_s = \text{diag}\{1, \dots, \lambda^{-a_{S,T}}\}$ describes the diagonal QCD running of scalar or tensor operators involving quarks with $\lambda = \frac{\alpha_s(\mu_f)}{\alpha_s(\mu_i)}$, $a_T = -4/23$, $a_S = 12/23$, $a_J = 0$ for $J \neq S, T$, and $a_d$ describes the running of $\Gamma_{KJ}^e$ in the case where it includes a running QCD parameter (eg $a_d = a_S$ for anomalous dimensions mixing quark tensor operators to the dipole). At $\mathcal{O}(\alpha_e^2)$, the scale-ordered exponential gives $\lesssim 10\%$ QCD corrections to the mixing of scalars to the dipole (“Barr-Zee”), so these QCD corrections are neglected in the analytic solutions given later.

<footnote n="2">²This allows to implement Wilsonian intuition about Renormalisation Group running of operator coefficients, while computing loops in dimensional regularisation.</footnote>

<footnote n="3">³The two-loop effects are included in the 'Hooft-Veltman scheme for $\gamma_5$, where the one-loop matching contributions vanish, so should give a scheme-independent result.</footnote>
The particle content of the EFT changes when the scale crosses a particle mass, so the operator basis changes too. At each such threshold, the Greens functions in the theory above and below must be identical, which allows to match coefficients from the basis above onto below. (Matching at the weak scale is recently discussed in [43].) Frequently, this matching can be performed at tree level in both theories; in some cases, the leading contribution of an operator arises via loop diagrams of the theory above, in which case these are included (hopefully after checking that they are scheme-independent, which is not always the case [42]). So for instance, the one-loop matching of scalar b, c quark operators onto the gluon operators, at the quark mass scale, is included as given in eq. (2.12).

What level of accuracy is required? The aim is to include the dominant contribution of each operator (at $\Lambda_{NP}$) to each process. Since almost every operator contributes to each observable (see the tables of appendix A), the question is whether the largest contributions have been included. This is difficult to demonstrate, because the dominant contribution may not be the lowest order in every perturbative expansion (loops, couplings, scale ratios); an example is the two-loop Barr-Zee diagrams, which give the dominant contribution of flavour-changing Higgs couplings to radiative decays like $\mu \to e\gamma$. However, the selection of terms included here has been checked against various models [44, 45]. Notice also that including effects $\gtrsim 10^{-3}$ does not mean the calculation is accurate to three figures; rather, the resulting constraints on operator coefficients have uncertainties $\gtrsim 10\%$.

Missing from the results given here, are most strong interaction effects beyond Leading Order, below 2 GeV. For example, quark loop contributions to $\mu \to e\gamma$ and $\mu \to e\bar{e}e$ are included at scales above 2 GeV, but should be replaced by pion loops or resonances from 2 GeV to $m_\mu$. An interesting study in this direction is [46]. The $\mu \to e$ conversion calculation is also at Leading Order in $\chi$PT.

## 2.4 Constraints at $m_W$

This section gives the experimental constraints of eqs. (2.5), (2.7) and (2.8) expressed in terms of operator coefficients at the weak scale.

### 2.4.1 $\mu \to e\gamma$

The MEG bound on coefficients at $m_W$ is

$$ 
\begin{aligned}
1.05 \times 10^{-8} \gtrsim & \left| C_{D,X} \left( 1 - 16 \frac{\alpha_e}{4\pi} \ln \frac{m_W}{m_\mu} \right) \right. \\
& \left. - \frac{\alpha_e}{4\pi e} \left( -8 \frac{m_\tau}{m_\mu} C_{T,XX}^{\tau\tau} \ln \frac{m_W}{m_\tau} + C_{S,XX}^{\mu\mu} \ln \frac{m_W}{m_\mu} + C_{S,XX}^{ee} \frac{m_e}{m_\mu} \ln \frac{m_W}{m_\mu} + C_{2\text{loop}} \ln \frac{m_W}{m_\tau} \right) \right. \\
& \left. + 8 \frac{\alpha_e^2}{e(4\pi)^2} \ln^2 \frac{m_W}{m_\mu} \left( \frac{m_\tau}{m_\mu} C_{S,XX}^{\tau\tau} \right) \right. \\
& \left. - 8 \lambda^{aT} \frac{\alpha_e}{4\pi e} \ln \frac{m_W}{2 \text{GeV}} \left( -\frac{m_s}{m_\mu} C_{T,XX}^{ss} + 2 \frac{m_c}{m_\mu} C_{T,XX}^{cc} - \frac{m_b}{m_\mu} C_{T,XX}^{bb} \right) f_{TD} \right. \\
& \left. + 8 \frac{\alpha_e^2}{3e(4\pi)^2} \ln^2 \frac{m_W}{2 \text{GeV}} \left( \sum_{u,c} 4 \frac{m_q}{m_\mu} C_{S,XX}^{qq} + \sum_{d,s,b} \frac{m_q}{m_\mu} C_{S,XX}^{qq} \right) \right| \quad (2.17)
\end{aligned}
$$
where all the coefficients are to be evaluated at $m_W$, but the masses are on-shell (= $m_ψ(m_ψ)$). Several comments are in order:

1. At one loop, all the tensor operators, and the scalars involving 3 muons or electrons, can mix to the dipole by closing same-flavour legs and attaching a photon. (However, the 3e scalar coefficient can be neglected, because $\mu \to e\bar{e}e$ sets a more stringent bound on this operator.)

(a) The $\tau$-tensor is in the second parenthese, the quark tensors are in the fourth (the light quark tensors are neglected, because they contribute at one loop to SI $\mu \to e$ conversion, unsuppressed by $m_q/m_\mu$). QCD causes $C_{T,XX}^{qq}$ to run, and also the quark mass that appears in the anomalous dimension mixing the tensor to the dipole is taken running. These QCD effects are described by $\lambda = \alpha_s(2\text{GeV})/\alpha_s(m_W) \simeq 2.18$, and

$$f_{TD} = \left( \frac{1}{1 - a_T - a_S} \frac{1 - \lambda^{1-a_T-a_S}}{1 - \lambda} \right) \qquad (2.18)$$

where $a_T = -4/23$, $a_S = 12/23$ are respectively the anomalous dimensions of tensors and scalars/masses in QCD, and $\lambda^{a_T} \simeq 0.873$ and $f_{TD} \simeq 0.862$. In eq. (2.17), the quark masses are at low energy: $m_q(m_q)$.

(b) It is interesting that the coefficients of the dipole and tensor operators can be comparable, allowing for unexpected cancellations. Numerically, the bounds of eq. (2.17) is

$$C_{D,X}(m_\mu) \gtrsim |0.938C_{D,X}(m_W) + 0.981C_{T,XX}^{\tau\tau}(m_W) - 0.75C_{T,XX}^{cc} + \dots| \quad (2.19)$$

2. $C_{2\text{loop}}$ is a combination of vector operator coefficients that mix to the dipole via the 2-loop RGEs [42]

$$C_{2\text{loop}} = -\frac{\alpha_e}{4\pi} \left[ \frac{58}{9} C_{V,YY}^{\tau\tau} + \frac{116}{9} \sum_{l=e,\mu} C_{V,YY}^{ll} + \frac{64}{9} (C_{V,YY}^{uu} + C_{V,YY}^{cc}) + \frac{22}{9} \sum_{q=d,s,b} C_{V,YY}^{qq} - \frac{80}{9} (C_{V,YX}^{uu} + C_{V,YX}^{cc}) - \frac{14}{9} \sum_{q=d,s,b} C_{V,YX}^{qq} - \frac{50}{9} \sum_{l=e,\mu,\tau} C_{V,YX}^{ll} + 4 \sum_{f=b,c,s,\tau} C_{S,YX}^{ff} \frac{Q_f^2 N_f m_f}{m_\mu} \right] \qquad (2.20)$$

where the chirality of the electron outgoing from these operators must be the same as that of the dipole operator into which they mix — so the combination above mix into $C_{D,X}$ for $X \neq Y$.

The $C_{2\text{loop}}$ term in eq. (2.17) is compact but not quite correct: in reality, the lower cutoff of the logarithm should be $m_\tau$ for the $\tau\tau$ operators, $m_b$ for the $b\bar{b}$ operators, 2 GeV for the $cc,ss,uu$ and $dd$ operators, and $m_\mu$ for the $\mu\mu$ and $ee$ operators.

3. the $\tau$ and (heavy) quark scalar operators mix to the tensor at one loop, and therefore the dipole at two-loop. Since the tensor-to-dipole mixing is $O(1)$ for heavy fermions because enhanced by the mass, the scalar-dipole mixing is included for the $\tau$ on the
third parenthese, and for the heavy quarks in the last parenthese. The QCD running is $\lesssim 10\%$, if the heavy quark masses are taken at low energy $m_Q(m_Q)$, so neglected. The lower cutoff of the logarithm for quark operators is approximated as 2 GeV to given a shorter formula, but should be chosen as a function of the quark flavour — eg $m_b$ for the b quark, and 2 GeV for u, d, s.

### 2.4.2 $\mu \to e\bar{e}e$

The decay $\mu \to e\bar{e}e$ constrains the magnitude of several combinations of coefficients, as given in eq. (2.6). At one-loop, there are no operators that mix into the scalar operators $O_{S,...}^{ee}$. All vector operators $O_{V,XY}^{ff}$ and $O_{V,XX}^{ff}$ can mix to $O_{V,XR}^{ee} + O_{V,XL}^{ee}$ by closing the $f$ legs and attaching a photon to the loop, then attaching the photon to $e^+e^-$. These “penguin” diagrams mix the combination

$$Q_e C_{ping1,X} = Q_e \frac{4}{3} \sum_f N_f Q_f (C_{V,XL}^{ff} + C_{V,XR}^{ff}) \quad f \in \{e, \mu, \tau, u, d, s, c, b\} \quad (2.21)$$

into both $C_{V,XR}^{ee}$ and $C_{V,XL}^{ee}$. The 3e and 3$\mu$ vector operators also mix to the 3e vectors via a second penguin diagram which closes a different selection of legs, and the 3e operators are renormalised by photon exchange between the legs. As a result, the constraints from $\mu \to e\bar{e}e$ (see eq. (2.7)), combined with the MEG bounds on $\mu \to e\gamma$, give the following bounds on coefficients evaluated at $m_W$:

$$ \begin{aligned} 7 \times 10^{-7} &> C_{V,XX}^{ee} \left( 1 - 12 \frac{\alpha_e}{4\pi} \ln \frac{m_W}{m_\mu} \right) - \frac{\alpha_e}{4\pi} \left( \frac{4}{3} \left( C_{V,XX}^{\mu\mu} + C_{V,XX}^{ee} \right) \ln \frac{m_W}{m_\mu} + C_{ping1,X} \ln \frac{m_W}{m_\tau} \right) \\ 10^{-6} &> C_{V,XY}^{ee} \left( 1 + 12 \frac{\alpha_e}{4\pi} \ln \frac{m_W}{m_\mu} \right) - \frac{\alpha_e}{4\pi} \left( \frac{4}{3} \left( C_{V,XX}^{\mu\mu} + C_{V,XX}^{ee} \right) \ln \frac{m_W}{m_\mu} + C_{ping1,X} \ln \frac{m_W}{m_\tau} \right) \\ 2.8 \times 10^{-6} &> C_{S,XX}^{ee} \left( 1 + 12 \frac{\alpha_e}{4\pi} \ln \frac{m_W}{m_\mu} \right) \end{aligned} \quad (2.22) $$

where $X \in \{L, R\}$ and $X \neq Y$. The overall logarithm multiplying $C_{ping1,X}$ is an approximation to shorten the formula; more correctly, the $b\bar{b}$ coefficients in $C_{ping1,X}$ should be multiplied by a logarithm cut off at $m_b$, and the logarithm for the $ee$ or $\mu\mu$ coefficients should end at $m_\mu$ (the correct logs are implemented in the numerical bounds in the appendix). (And as usual, there should be hadronic loops to replace the light quarks below 2 GeV.)

### 2.4.3 $\mu \to e$ conversion

The SINDRUMII bound on $\mu$Au$\to e$+Au, given in eq. (2.8) for out-going left-handed electrons, can be expressed as the bound on coefficients at $m_W$:

$$ \begin{aligned} 4.9 \times 10^{-8} \gtrsim & \left| \left( 0.222 C_{D,R}(m_\mu) + 1.602 C_{V,L}^{uu} + 1.830 C_{V,L}^{dd} - 0.721 \frac{8\pi m_N}{9\alpha_s m_t} C_{GG,R} \right) \right. \\ & + \lambda^{aS} \left( 1 + \frac{\alpha}{4\pi} \tilde{\ln} \right) \left( \frac{26}{3} \left[ 6.10 C_{S,R}^{uu} + \frac{1.44 m_N}{27 m_c} C_{S,R}^{cc} \right] \right. \\ & \left. \qquad + \frac{20}{3} \left[ 6.26 C_{S,R}^{dd} + 0.303 C_{S,R}^{ss} + \frac{1.44 m_N}{27 m_b} C_{S,R}^{bb} \right] \right) \\ & - 16 f_{TS} \lambda^{aT} \frac{\alpha}{4\pi} \tilde{\ln} \left( 2[6.100 C_{T,RR}^{uu} + \frac{1.44 m_N}{27 m_c} C_{T,RR}^{cc}] \right. \\ & \left. \qquad - 6.258 C_{T,RR}^{dd} - 0.303 C_{T,RR}^{ss} - \frac{1.44 m_N}{27 m_b} C_{T,RR}^{bb} \right) \\ & \left. - \frac{\alpha}{2\pi} \tilde{\ln} \left( 3.66 C_{A,L}^{dd} + 6.41 C_{A,L}^{uu} - 0.305[C_{V,L}^{\mu\mu} - C_{A,L}^{\mu\mu}] + 0.228 C_{ping,1} \right) \right| \end{aligned} \quad (2.23) $$
where $\tilde{\ln} = \ln(m_W/m_{low})$ should be inside the square brackets because $m_{low} \in \{m_b, 2 \text{ GeV}, m_\mu\}$ depends on the coefficient, the masses are on shell ($m_\psi(m_\psi)$) and the dipole contribution is given at the experimental scale (it can be written in terms of coefficients at $m_W$ with eq. (2.17)). A similar bound applies with $L$ and $R$ interchanged.

The bound from Titanium on an orthogonal combination of coefficients, given in eq. (2.14), can be written with the same conventions as

$$ 
\begin{aligned}
\sqrt{\frac{BR_{Ti}^{exp}}{B_{Ti} \sin^2 \phi}} = 6.24 \times 10^{-7} &\gtrsim \left| (0.155C_{D,R} - 0.506C_{V,L}^{uu} - 1.17C_{V,L}^{dd} - 1.06\frac{8\pi m_N}{9\alpha_s v}C_{GG,R}) \right. \\
&\qquad + \lambda^{as} \left(1 + \frac{\alpha}{4\pi}\tilde{\ln}\right) \left(\frac{26}{3}[9.21C_{S,R}^{uu} + \frac{2.12m_N}{27m_c}C_{S,R}^{cc}] \right. \\
&\qquad \left. + \frac{20}{3}[9.02C_{S,R}^{dd} + 0.446C_{S,R}^{ss} + \frac{2.12m_N}{27m_b}C_{S,R}^{bb}] \right) \\
&\qquad - 16f_{TS}\lambda^{aT}\frac{\alpha}{4\pi}\tilde{\ln}\left(2[9.21C_{T,RR}^{uu} + \frac{2.12m_N}{27m_c}C_{T,RR}^{cc}] \right. \\
&\qquad \left. - 9.02C_{T,RR}^{dd} - 0.446C_{T,RR}^{ss} - \frac{2.12m_N}{27m_b}C_{T,RR}^{bb} \right) \\
&\qquad - \left. \frac{\alpha}{2\pi}\tilde{\ln}(2.576C_{A,L}^{dd} - 2.024C_{A,L}^{uu} - 0.0346[C_{V,L}^{\mu\mu} - C_{A,L}^{\mu\mu}] + 0.32C_{ping,1}) \right|
\end{aligned}
\quad (2.24) $$

where $\tilde{\ln}$ is defined after eq. (2.23) and the masses are on shell ($m_\psi(m_\psi)$).

## 2.5 What is not probed?

This section gives an incomplete list of coefficient combinations that are not probed by $\mu \to e\gamma$, $\mu \to e\bar{e}e$, and $\mu \to e$ conversion. There are 45 $\mu \to e_L$ operator coefficients in section 2.1, upon which the current bounds on $\mu \to e\gamma$, $\mu \to e\bar{e}e$, $\mu\text{Au} \to e_L\text{Au}$ and $\mu\text{Ti} \to e_L\text{Ti}$ set $1+3+1+1=6$ constraints [11]. Including $\mu \to e\gamma\gamma$, which was searched for by Crystal Box [21] would give an additional constraint (on a $\gamma\gamma$ operator) [20, 22], but there remain $\lesssim 40$ unconstrained directions in coefficient space. The case of operators that produce an outgoing $e_R$ is the same and independent, because the Branching Ratios independently constrain both processes.

This section does not give a list of three (or six) dozen flat directions. The subset of flat directions listed here are selected because they are relatively stable under RG running, or “natural”, according to the notion introduced it [11] (for cancellations among operator coefficients). It assumes that model parameters at the New Physics scale $\Lambda_{\text{NP}}$ do not know the experimental scale, so coefficients should not cancel against logarithms of scale ratios — because the scale ratio could be chosen by the observer but the coefficient is determined by the theory.

An elegant way to implement this notion of naturalness, is to only admit cancellations which are RG-invariant. For instance, if one restricts to operators containing a quark bilinear, and only considers their one-loop diagonal QCD running, then cancellations among the vector coefficients, or among the scalars, or the tensors, are RG-invariant. This could be generalised to include the one-loop QCD and QED running, by diagonalising the full anomalous dimension matrix. However, the diagonalisation would be difficult, and the eigenvectors would be curious combinations of different Lorentz structures and external legs — as opposed to the usual “intuitive” basis of section 2.1, where the operators correspond to potentially distinguishable interactions of physical particles. Furthermore, the only
degenerate eigenvalues might be among vector operators involving fermions of the same electric charge (because scalars mix to tensors whose mixing to the dipole is proportional to the fermion mass).

A more pragmatic implementation of this notion of “natural” cancellations is to use the perturbative solution of the RGEs given in eq. (2.16), and allow cancellations among coefficients which 1) run with the same QCD anomalous dimension, and 2) multiply a similar $\alpha_{em} \ln \mu_f / \mu_i$ factor. The formulae for the constraints on coefficients at the weak scale, given in eqs. (2.17), (2.22), (2.23), and (2.24), are arranged to display these “natural” cancellations, or flat directions: each parenthese of the formulae satisfies the two conditions. So any combination of coefficients that causes a parenthese to vanish is a “natural” flat direction.

An example of flatish directions, or operators which are poorly constrained by $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu \to e$ conversion are the axial vector operators of the form

$$ (\bar{e}\gamma_{\alpha} P_X \mu)(\bar{f}\gamma^{\alpha} \gamma_5 f), \quad f \in \{s, c, b, \tau\} $$

which contribute to $\mu \to e\gamma$ at 2-loop with a mass enhancement (see eq. (2.20)). There should be three combinations of these operators, orthogonal to eq. (2.20), which are “flat” to the order calculated here. (The axial vector operators with $f \in \{u, d, \mu\}$ are not in the list, because they contribute at one QED loop to $\mu \to e$ conversion.) A second example would be any combination of vector operators

$$ \sum_f C_f(\bar{e}\gamma_{\alpha} P_X \mu)(\bar{f}\gamma^{\alpha} f), \quad f \in \{s, c, b, \tau\} $$

where the $\{C_f\}$ are chosen orthogonal to $C_{ping,1}$ (see eq. (2.21)), and also to $C_{2loop}$ if one wishes large $C_f$.

Reference [11]’s notion of naturalness precludes cancellations among the combinations of coefficients in parentheses in eqs. (2.17), (2.22), (2.23), and (2.24), so it could be interpreted as transforming the single constraint from $\mu \to e\gamma$ into a constraint on each of the five parentheses, and the 2 bounds from $\mu \to e$ conversion into eight. However, “unnatural” cancellations can occur, see for instance eq. (2.19).

## 3 Complementarity

The aim of this section is to show that the observables considered here ($\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu \to e$ conversion) give independent information about models. To do this, an alternative basis is proposed for coefficient space. This basis is constructed from observables, spans the subspace they probe, and can be defined at all scales.

In an EFT framework, a model $M$ can be represented by the vector $\vec{C}_M(\Lambda)$ of operator coefficients it induces,<fnref n="4" /> which is scale-dependent but relatively simple to calculate at the New Physics scale $\Lambda_{NP}$. An observable can be represented as one or several combinations

<footnote n="4">⁴Since models usually have parameters, there would be a vector for each choice of parameters, and varying the model parameters would scan over the vectors of coefficients that can represent the model.</footnote>
of operators whose coefficients it probes. This also depends on the scale, and is relatively simple to establish at the experimental scale $\Lambda_{\text{expt}}$ — for instance, $\mu \to e_L\gamma$ probes the operator $O_{D,R}(m_\mu)$. However, to obtain the rate for an observable, one must evaluate the matrix element of the operators and integrate phase space. So it is convenient to define observable-vectors $\vec{v}_{\text{obs}}(\Lambda)$ to live in the vector space of the coefficients, such that the rate for an observable can be written

$$ \Gamma \propto \sum |\vec{C}_{\mathcal{M}}(\Lambda) \cdot \vec{v}_{\text{obs}}(\Lambda)|^2 \quad (3.1) $$

at a common scale $\Lambda$. For instance, at the experimental scale, eq. (2.5) implies that the two observable-vectors for $\mu \to e\gamma$ point in the dipole directions, so one can choose

$$ \begin{aligned} \vec{v}_{\mu \to e_L \gamma}(m_\mu) &= \hat{u}_{D,R} \\ \vec{v}_{\mu \to e_R \gamma}(m_\mu) &= \hat{u}_{D,L}, \end{aligned} \quad (3.2) $$

where $\hat{u}_{D,R}$, $\hat{u}_{D,L}$ are unit vectors such that $\vec{C} \cdot \hat{u}_{D,R} = C_{D,R}$. Some observable-vectors for $\mu \to e$ conversion at the experimental scale are already given in eq. (2.11).

By assumption, the dynamics below $\Lambda_{\text{NP}}$ is Standard Model, so the model and observable vectors can be translated in scale by the SM RGEs. Since models are legion and observables are few, it might be more efficient to translate the observable vectors to $\Lambda_{\text{NP}}$, rather than calculating $\vec{C}_{\mathcal{M}}(\Lambda_{\text{expt}})$, as is commonly done. Indeed, if the observable-vectors were known as a function of $\Lambda_{\text{NP}}$, then the predictions of a New Physics model would be simple to obtain from eq. (3.1). Part of the translation $\vec{v}_{\text{obs}}(\Lambda_{\text{expt}}) \to \vec{v}_{\text{obs}}(\Lambda_{\text{NP}})$ appears in section 2, where the combination $\vec{C}_{\mathcal{M}}(m_W) \cdot \vec{v}_{\mu \to e_L \gamma}(m_W)$ is given by eq. (2.17) with $X = R$, $\vec{C}_{\mathcal{M}}(m_W) \cdot \vec{v}_{\mu \to e_L \gamma}(m_W)$ appears in eq. (2.23), and the appropriate inner product for the orthogonal constraint from $\mu \text{Ti} \to e_L \text{Ti}$ is given in eq. (2.24). Similar formulae apply for outgoing $e_R$, but the focus below is mostly on $e_L$. Matching onto the SMEFT, and running up to $\Lambda_{\text{NP}}$ is left for a later work.

The observable vectors corresponding to $\mu \to e\bar{e}e$ can also be constructed. For this, it is convenient to take an idealised theoretical perspective, imagining that the chirality of the four leptons can be observed. Combined with the electron angular distributions [24, 25], this allows to define “observable” vectors for $\mu \to e\bar{e}e$ at the experimental scale

$$ \begin{aligned} \vec{v}_{\mu_L \to e_L \bar{e}_L e_L}(m_\mu) &= \sqrt{2}[\hat{u}_{V,LL}^{ee} + 4e\hat{u}_{D,R}] \\ \vec{v}_{\mu_R \to e_L \bar{e}_R e_L}(m_\mu) &= \frac{1}{2\sqrt{2}}\hat{u}_{S,LL}^{ee} \\ \vec{v}_{\mu_R \to e_R \bar{e}_L e_L}(m_\mu) &= \hat{u}_{V,RL}^{ee}(m_\mu) + 4e\hat{u}_{D,R} \end{aligned} \quad (3.3) $$

plus another three vectors with $L \leftrightarrow R$, and the dipoles. These vectors can be expressed at $m_W$ using eqs. (2.17) and (2.22).

These vectors corresponding to $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu A \to eA$, allow to quantify the complementarity of these processes at any scale. If, at the chosen scale, the vectors remain “relatively orthogonal”, then the observables are complementary. If the vectors are aligned, then the observables become identical probes of New Physics. It is interesting
to study complementarity at $\Lambda_{NP}$, because we are looking for information about the New Physics. However, as an illustrative exercise, we start at the experimental scale.

Notice that the complementarity of observables is a theoretical question, unrelated to the magnitude of current experimental bounds. In this it differs from the correlation matrix that can be constructed from the experimental constraints, which defines the allowed ellipse in coefficient space; the direction and magnitude of the axes of the ellipse depend on the relative magnitude of the constraints. Complementarity is also model-independent, and unrelated to correlations between observables that could arise in some models. Such correlations would occur if the projections of $\vec{C}_M$ onto the observable vectors have similar scaling with the parameters of model M:

$$ \frac{\vec{C}_M \cdot \vec{v}_{obs1}}{\vec{C}_M \cdot \vec{v}_{obs2}} = \text{independent of model parameters.} $$

The notion of “model discriminating power”, present in the literature [19, 47, 48], is related to complementarity, but compares the projection of different observable vectors, onto different model vectors. The discriminating power of observables therefore depends on the models considered, whereas their complementarity is model-independent.

At the experimental scale, $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu \to e$ conversion have different external particles, so naively appear complementary. However, the dipole contributes to all three processes, preventing the observable-vectors from being orthogonal even at the experimental scale. Indeed, at $\Lambda \sim m_\mu$, $C_{D,R}(m_\mu)$ and $C_{V,LL}(m_\mu)$ contribute about equally to $\mu_L \to e_L \bar{e}_L e_L$, so $\vec{v}_{\mu_L \to e_L \bar{e}_L e_L}(m_\mu)$ is misaligned with respect to $\vec{v}_{\mu_R \to e_L \gamma}(m_\mu)$ by $\sim \pi/4$:

$$ \frac{\vec{v}_{\mu_L \to e_L \bar{e}_L e_L} \cdot \vec{v}_{\mu_R \to e_L \gamma}}{|\vec{v}_{\mu_R \to e_L \gamma}| |\vec{v}_{\mu_L \to e_L \bar{e}_L e_L}|} (m_\mu) \sim \cos(\pi/4). \quad (3.4) $$

The dipole also contributes to $\mu \to e$ conversion, but with a smaller weight (see eq. (2.13)), such $\vec{v}_{\mu_R \to e_L \gamma}(2 \text{ GeV})$ and $\vec{v}_{\mu_{Au} \to e_L \text{Au}}(2 \text{ GeV})$ are almost orthogonal (separated by an angle of $\sim 88$ degrees). So at low energy, an approximately orthogonal basis for the observable-vector subspace can be constructed with the two dipoles, and the remaining observable-vectors with the dipole subtracted. For instance,

$$ \begin{aligned} \vec{u}_{\mu \to e_L \gamma}(\Lambda) &= \vec{v}_{\mu \to e_L \gamma}(\Lambda) \\ \vec{u}_{\mu \text{Au} \to e_L \text{Au}}(2 \text{ GeV}) &= \vec{v}_{\mu \text{Au} \to e_L \text{Au}}(2 \text{ GeV}) - 0.222 \vec{v}_{\mu_R \to e_L \gamma}(2 \text{ GeV}) \\ \vec{u}_{\mu_L \to e_L \bar{e}_L e_L}(m_\mu) &= \vec{v}_{\mu_L \to e_L \bar{e}_L e_L}(m_\mu) - 4\sqrt{2} \vec{v}_{\mu_R \to e_L \gamma}(m_\mu) \end{aligned} \quad (3.5) $$

which are called $\vec{u}$ to distinguish them from the observable-vectors $\vec{v}$.

This situation at low energy is illustrated in figure 1, where low energy is taken to be 2 GeV in order to use quark operators for $\mu \to e$ conversion. The plot is in polar coordinates in the subspace probed by $\mu \to e_L \gamma$, $\mu \to e_L \text{Au}$ and $\mu \to e_L \bar{e}_L e_L$, where the vertical z-axis corresponds to the dipole $\vec{u}_{\mu \to e_L \gamma}(2 \text{ GeV})$, and the x and y axes respectively to $\vec{u}_{\mu_L \to e_L \bar{e}_L e_L}(2 \text{ GeV})$ and $\vec{u}_{\mu \text{Au} \to e_L \text{Au}}(2 \text{ GeV})$. This subspace corresponds to a Lagrangian at the experimental scale of

$$ \mathcal{L}_{\text{eff}}(m_\mu) = \frac{1}{\Lambda_{\text{NP}}^2} \left[ c_\theta m_\mu \bar{e}\sigma \cdot FP_R\mu + s_\theta c_\phi (\bar{e}\gamma^\alpha P_L\mu)(\bar{e}\gamma_\alpha P_L e) + s_\theta s_\phi \mathcal{O}_{\text{Au}} \right] \quad (3.6) $$
where $O_{Au}$ is the combination of quark operators probed on Gold, $s_\theta = \sin\theta$, $c_\phi = \cos\phi$ and so on, and the radial coordinate is taken $|\vec{C}| = v^2/\Lambda_{NP}^2$, so that

$$2\sqrt{2}G_F|\vec{C}|^2 = \frac{1}{\Lambda_{NP}^2}.$$

This is the complete low-energy Lagrangian for the considered observables; the only “toy” aspect is that chirality of the electrons is fixed to make the parameter space of plot-table dimension.

The resulting BRs are:

$$ 
\begin{aligned}
BR(\mu \to e_L \gamma) &= 384\pi^2 |\vec{v}_{\mu \to e_L \gamma} \cdot \vec{C}|^2 = 384\pi^2 |\vec{C}|^2 |\cos \theta|^2 \\
BR(\mu \to e_L \overline{e_L} e_L) &= |\vec{v}_{\mu_L \to e_L \overline{e_L} e_L} \cdot \vec{C}|^2 + 18.76 |\vec{v}_{\mu \to e_L \gamma} \cdot \vec{C}|^2 \\
&= |\vec{C}|^2 (2 |\sin \theta \cos \phi + 1.2 |\cos \theta|^2 + 18.76 |\cos \theta|^2) \quad (3.7) \\
BR(\mu \text{Au} \to e \text{Au}) &= 300 |\vec{v}_{\mu A \to e_L A} \cdot \vec{C}|^2 = 300 |\vec{C}|^2 [0.222 \cos \theta + 9.08 \sin \theta \sin \phi]^2.
\end{aligned}
$$

Setting the BRs equal to their experimental limits (given in section 2.2), gives a bound on $|\vec{C}|$ which is plotted in figure 1 as a function of $\theta, \phi$; the coefficient space below the curves is experimentally excluded. The angular dependence of the bounds illustrates that the observables are the complementary at the experimental scale. The $BR(\mu \to e_L \overline{e_L} e_L)$ vanishes at $\phi = \pi/2$ in the second plot (where $\theta = \pi/2$) because both the dipole and four-fermion contributions vanish; the BR does not quite vanish at $\theta = \pi/2$ in the first plot because the subdominant four-fermion contribution is present at $\phi = \pi/4$. One also sees that $\mu \to e$ conversion gives the best bound on $\Lambda_{NP}$, for comparable coefficients of all three operators. This is because the conversion process is coherent on the nucleus (so is enhanced at large atomic number), and because the scalar quark densities in the nucleon are large (~9, see table 1 — the Branching Ratio due to a vector quark current would be ~10 → 100 smaller).

Consider now the complementarity of these observables at $m_W$, $m_W$ being a simple, although inadequate, substitute for $\Lambda_{NP}$. In their evolution from the experimental scale to the weak scale, the basis vectors, like observable-vectors, will rotate in coefficient space and change length (the non-unitary matrix that evolves coefficients in scale can be written as a rotation, a diagonal matrix, and another rotation). The degree of orthogonality between the vectors can be obtained by taking inner products; at the weak scale, the basis vectors for the processes plotted here (= dipoles, plus observable-vectors with their dipole components subtracted) are still orthogonal to three figures despite their significant rotations (this is understandable; the vector space has ~90 dimensions). But their lengths change; in particular the dipole vectors grow by more than a factor two, due to the large loop contributions shown in eq. (2.17).

The observables nonetheless remain complementary at $m_W$, as shown in figure 2. These plots are similiar to those of figure 1, but differ in that the operator coefficients are at $m_W$ (standing in for $\Lambda_{NP}$). The first plot shows $\Lambda_{NP}$ as a function of the angle $\theta$ between the model vector $\vec{C}(m_W)$, and the direction probed by the dipole $\vec{u}_{\mu \to e_L \gamma}(m_W)$ (at $\phi = \pi/4$, that is a model with identical coefficients in the directions $\hat{u}_{\mu_L \to e_L \overline{e_L} e_L}$ and $\hat{u}_{\mu \text{Au} \to e_L \text{Au}}$). The shaded regions are excluded by current constraints. The second plot is the projection
<figure><img src="image_3.png" /><figcaption>**Figure 1.** The bound on the New Physics scale $\Lambda_{NP}$ from current constraints on $\mu \to e_L\gamma$ (thick black), $\mu \to e_L\bar{e}_L e_L$ (dot-dashed blue) and $\mu A \to e_L A$ (thin red), for New Physics parametrised at the experimental scale by the effective Lagrangian given in eq. (3.6). The first figure is as a function of $\theta$ for $\sin\phi = \cos\phi = 1/\sqrt{2}$, and the second is as a function of $\phi$ for $\cos\theta = 0$.</figcaption></figure>

of the model-vector onto the plane perpendicular to the dipole. Were this figure at $\Lambda_{NP}$, rather than $m_W$, and if the three Branching Ratios were observed, it would allow to extract all the information these rates give about New Physics (via dimension six operators and leading order RGEs), at the New Physics scale, so without the blurring effect of SM loops. As expected, $\mu \to e\gamma$ is maximal at $\theta = 0, \pi$, and vanishes for $\theta = \pi/2$, and $\mu \to e$ conversion is larger when the model vector is more orthogonal to the dipole. The $BR(\mu \to e_L \bar{e}_L e_L)$ follows $BR(\mu \to e_L \gamma)$ due to the enhanced (by RG running) contribution of the dipole operator to $\mu \to e_L \bar{e}_L e_L$; it does not vanish with $\vec{C} \cdot \vec{u}_{\mu \to e_L \gamma}$ at $\theta = \pi/2$ due to the four-fermion contribution. Similarly, $BR(\mu \text{Au} \to e_L \text{Au})$ vanishes at $\theta \lesssim \pi$ where there is a cancellation between the negative $\vec{C} \cdot \vec{u}_{\mu \to e_L \gamma}$ and positive $\vec{C} \cdot \vec{u}_{\mu \text{Au} \to e_L \text{Au}}$.

# 4 Discussion and summary

Reconstructing New Physics (NP) from data is a dream for many phenomenologists. If New Physics is heavy, then Effective Field Theory can be a tool in pursuing this dream, because
<figure><img src="image_4.png" /><figcaption>**Figure 2.** The bound on $\Lambda_{\text{NP}}$ from current constraints on $\mu \to e_L\gamma$ (black; vanishes at $\theta = \pi/2$), $\mu \to e_L\bar{e}_L e_L$ (blue; dips at $\theta, \phi \sim \pi/2$) and $\mu A \to e_L A$ (red; vanishes at $\theta, \phi \sim \pi$), as a function of coefficients at $m_W$. Dotted lines indicate possible bounds from future experimental sensitivities. The operator coefficients are parametrised in spherical coordinates in the subspace probed by the experimental processes: the vertical axis is the direction probed by $\mu \to e_L\gamma$ (see eq. (2.17)) and the $xy$ plane is spanned by $\vec{u}_{\mu\text{Au} \to e_L\text{Au}}(m_W)$ and $\vec{u}_{\mu_L \to e_L\bar{e}_L e_L}(m_W)$, see eq. (3.5). The first plot is for $\phi = \pi/4$, and the second for $\theta = \pi/2$ (where the dipole vanishes). See discussion at the end of section 3.</figcaption></figure>

it allows to separate what is known — the Standard Model(SM) and low energy data about NP — from the NP that is pursued. In particular, it allows theoretical travels in energy scale, from the experimental scale towards the NP scale $\Lambda_{\text{NP}}$, because the dynamical degrees of freedom are by assumption in the SM. However, in this manuscript, the experimental constraints are translated only as far as $m_W$; reaching $\Lambda_{\text{NP}}$ is left for later work.

The processes considered here are $\mu \to e\gamma$, $\mu \to e\bar{e}e$, and $\mu \to e$ conversion, because current experimental constraints are stringent ($BR \lesssim 10^{-12}$) and are expected to improve significantly in coming years ($\to 10^{-16}$). These are low-energy processes, occuring at an energy-transfer $\sim m_\mu$, whose experimental Branching Ratios (BRs) are reviewed in section 2.2. Together, they set $\sim$ a dozen constraints on $\mu \leftrightarrow e$ contact interactions at the experimental scale.

We would like to know what these processes can tell us about NP, with as few assumptions about the NP as possible. A first assumption is that the new particles are heavy, with masses at a scale $\Lambda_{\text{NP}} \gtrsim m_W$ (possibly $\gg m_W$). Furthermore, they are assumed to generate some three or four-particle, $\mu \to e$ flavour-changing contact interactions (a list is given in section 2.1).

With these assumptions, section 2 explores whether, if $\mu \leftrightarrow e$ flavour-changing NP exists, we should see $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and/or $\mu \to e\gamma$? In order to reach the short-distance NP interaction that mediates a low-energy $\mu \leftrightarrow e$ transition, the intermediate-scale SM loop corrections must be peeled off, for instance using the Renormalisation Group Equations (RGEs) summarised in section 2.3. This scale evolution can transform one $\mu \leftrightarrow e$ interaction into another, so at short distances/high scales, the experimental constraints apply to lengthly combinations of coefficients. These are given in section 2.4 at the scale $m_W$, and the appendix tabulates the sensitivities to individual operator coefficients (one-at-a-time bounds), both in the basis of section 2.1 and in the SMEFT. These results
include one-loop and some two-loop effects, and should give the leading contribution of each operator to each process (see discussion in section 2.3). The tables show that at this accuracy, $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu \to e$ conversion are sensitive to (almost) all the operators or section 2.1 (the possible exceptions are a few $\mu$-$e$-gluon-gluon operators — see eq. (2.4) and following discussion). Recall that the operator list represents all QED$\times$QCD-invariant LFV three or four-legged contact interactions, in a chiral representation for fermions.

It is unclear if this is reassuring, because cancellations in a matrix element can occur among different coefficients. Such directions in coefficient space are sometimes refered to as flat directions. There are many such flat directions in the coefficient space for $\mu \to e$ flavour change, after imposing the dozen bounds from $\mu \to e\gamma, \mu \to e\bar{e}e$ and $\mu \to e$ conversion, because there are $O(100)$ operators. Indeed, there are so many flat directions, some one could be motivated to introduce a limit on how much cancellation can be “naturally” allowed. The flat directions are discussed in section 2.5.

However, more interesting that the flat directions, are the directions that the observables do constrain, because these are what we can use to discriminate among models. So section 3 introduces vectors in coefficient space, which correspond to observables (these are patterned on the “target vectors” of [11, 28]). At the experimental scale, there can be several observable-vectors for a give process, each selects a combination of coefficients who interfere in the rate:

$$ \Gamma \propto \sum |\vec{C}_{\mathcal{M}}(\Lambda_{\text{expt}}) \cdot \vec{v}_{\text{obs}}(\Lambda_{\text{expt}})|^2 $$

where $\vec{C}_{\mathcal{M}}(\Lambda_{\text{expt}})$ is the vector of coefficients corresponding to model $\mathcal{M}$, evaluated at the experimental scale $\Lambda_{\text{expt}}$. For instance, for $\mu \to e\gamma$, whose BR is given in eq. (2.5), the two observable-vectors at the experimental scale are the unit vectors that select the coefficients of the dipole operators in the vector $\vec{C}_{\mathcal{M}}(\Lambda_{\text{expt}})$. The vectors for the remaining observables are given in section 3. The observable-vectors are interesting, because they are scale-dependent, and can be translated to $\Lambda_{\text{NP}}$ with the Renormalisation Group Equations. Then there would be no need to run operator coefficients down to the experimental scale for every model; rather, rates could be computed from coefficients at $\Lambda_{\text{NP}}$, by dotting them into $\vec{v}_{\text{obs}}(\Lambda_{\text{NP}})$ This is left for future work in the lepton sector.

In section 3, the observable vectors are used to study the complementarity of experimental processes. The vector(s) corresponding to a given observable identifies the subspace of coefficients that the observable probes, so the misalignment angle between the observable-vectors of different processes, quantifies the complementarity of the processes. When the vectors evolve in scale, this misalignment angle can grow or shrink, indicating that the observables become more, or less, complementary at high scales. For the purpose of learning about New Physics, clearly it is desirable for observables to be complementary at $\Lambda_{\text{NP}}$. This manuscript only reaches the weak scale in RGE evolution, so figure 2 illustrates the complementarity of $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu \to e$ conversion at the scale $m_W$.

## A Tables of sensitivities

The tables in this appendix give the sensitivities of various processes (listed in the first row of the tables) to the operator coefficients given in the left column and evaluated at $m_W$.
The “sensitivity” of an (unobserved) process to a coefficient is calculated by allowing only that coefficient to be non-zero at $m_W$. Coefficients smaller than these values are too small to have been observed, but larger coefficients could be allowed, if their contributions cancel against other coefficients.

The tabulated results arise from the experimental bounds given in eqs. (2.5), (2.6), (2.8) for muon decays, and in table 2 for $\tau$ decays.

Tables 3, 4, 8 and 9 contain the QED×QCD-invariant operators relevant below the weak scale, and listed in section 2 for $\mu \to e$ flavour change. The operators are added to the Lagrangian with the normalisation given in eq. (2.1):

$$ \mathcal{L} = \mathcal{L}_{\text{SM}} + 2\sqrt{2}G_F \sum C_{\text{Lor}}^{\zeta} \mathcal{O}_{\text{Lor}}^{\zeta} + \text{h.c.}, \quad 2\sqrt{2}G_F = \frac{1}{v^2}, \quad v \approx m_t $$

where the subscript is the Lorentz contraction, and the superscript $\zeta$ represents the flavour indices. Notice that in the normalisation used here, all operators annihilate muons (or $\tau$s) and create electrons; the reverse process is assured by the $+\text{h.c.}$. So the coefficients of an operator and its conjugate have the same magnitude.

At the weak scale $m_W$, the low energy operator basis can be matched onto the SMEFT, so the sensitivities can be expressed in this basis. Tables 5 to 7 (for $\mu \leftrightarrow e$) and 10 (for $\tau \leftrightarrow e$) apply to operators in the SMEFT basis, added to the Lagrangian as:

$$ \mathcal{L} = \mathcal{L}_{\text{SM}} + 2\sqrt{2}G_F \sum C_J^{\zeta} \mathcal{O}_J^{\zeta} + \text{h.c.} \qquad (A.1) $$

where the sum is over all dimension six operators and all flavour indices. In the usual SMEFT formulation, the $+\text{h.c.}$ is not included for hermitian operators; here there is $+\text{h.c.}$ for all operators but the hermitian operators are defined with a factor 1/2 (which gives the usual normalisation for coefficients, because the operator and its conjugate contribute to the Feynman rule.). The SMEFT convention of summing all flavour indices causes some four-lepton operators to appear several times:

$$ \begin{aligned} \mathcal{O}_{LL}^{e\mu\ell\ell} &= \mathcal{O}_{LL}^{\ell\ell e\mu} = \mathcal{O}_{LL}^{\ell\mu e\ell} = \mathcal{O}_{LL}^{e\ell\ell\mu} \\ \mathcal{O}_{LL}^{e\mu\tau\tau} &= \mathcal{O}_{LL}^{\tau\tau e\mu}, \quad \mathcal{O}_{LL}^{\tau\mu e\tau} = \mathcal{O}_{LL}^{e\tau\tau\mu} \\ \mathcal{O}_{RR}^{e\mu\ell\ell} &= \mathcal{O}_{RR}^{\ell\ell e\mu} = \mathcal{O}_{RR}^{\ell\mu e\ell} = \mathcal{O}_{RR}^{e\ell\mu} \end{aligned} \qquad (A.2) $$

for $\ell \in \{e, \mu\}$ and $l \in \{e, \mu, \tau\}$. So the coefficients of these identical operators are also identical, and the bounds in table 6 apply to the appropriate sum of coefficients:

$$ \begin{aligned} \sum C_{LL}^{e\mu\ell\ell} &= C_{LL}^{e\mu\ell\ell} + C_{LL}^{\ell\ell e\mu} + C_{LL}^{\ell\mu e\ell} + C_{LL}^{e\ell\ell\mu} = 4C_{LL}^{e\mu\ell\ell} \\ \sum C_{LL}^{e\mu\tau\tau} &= C_{LL}^{e\mu\tau\tau} + C_{LL}^{\tau\tau e\mu} = 2C_{LL}^{e\mu\tau\tau} \\ \sum C_{LL}^{\tau\mu e\tau} &= C_{LL}^{\tau\mu e\tau} + C_{LL}^{e\tau\tau\mu} = 2C_{LL}^{e\tau\tau\mu} \\ \sum C_{RR}^{e\mu\ell\ell} &= C_{RR}^{e\mu\ell\ell} + C_{RR}^{\ell\ell e\mu} + C_{RR}^{\ell\mu e\ell} + C_{RR}^{e\ell\mu} \end{aligned} \qquad (A.3) $$

Finally, the CKM matrix was neglected in matching, so eg, low energy operators $\mathcal{O}_{V,RL}^{uu}$ and $\mathcal{O}_{V,RL}^{dd}$ both match onto the SMEFT operator $\mathcal{O}_{EQ}^{dd}$.
Tables 3, 4, and 5 to 7 give the sensitivity of $\mu \to e\gamma$, $\mu \to e\bar{e}e$ and $\mu \to e$ conversion to $\mu \leftrightarrow e$ operators. Since the dipole receives loop contributions from almost all operators, it is sensitive to most coefficients. It also contributes to both $\mu \to e\bar{e}e$ and $\mu \to e$ conversion; when the only contribution of a coefficient to these processes is via the dipole, the sensitivity is in parentheses in tables 3 and 4.

The $\mu \leftrightarrow e$ sensitivities are given to three figures to mitigate rounding uncertainty; modulo misprints and factors of 2, they are expected to have a ~10% uncertainty due to the one-loop QCD running, and a ~50% uncertainty for scalar, tensor and GG operators in $\mu A \to eA$, where the lattice and $\chi$PT determinations of the scalar quark current in the nucleon differ by ~50% [32, 33].

For comparison, tables 8, 9 give estimated sensitivities of a selection of LFV $\tau$ decays to operator coefficients in the QED×QCD-invariant basis. In addition, table 10 gives the sensitivities of some hadronic $\tau$ decays to $2\ell2q$ operators of SMEFT — this illustrates, in the case of SMEFT operators involving quark doublets, the cancellations that can arise between the $u$ and $d$ quark contributions in $\tau$ decays to isotriplet mesons.

**Note added.** During the completion of this manuscript, appeared a comprehensive study of LFV in hadronic $\tau$ decays [49]. The authors calculate decays to a variety of final states using $\chi$PT with resonances, and perform some loop matching calculations in order to include the interesting operators $O_{GG,Y}$ [47]. The results in this manuscript are less precise, include $\tau \to \ell\gamma$ and $\tau \to 3\ell$ but fewer hadronic decays, and include QED loops in the RG running up to the weak scale. (QED effects are interesting, because they can change the external legs and Lorentz structure, transforming difficult-to-detect operators into well-constrained ones.) Husek et al. include the numerically significant QCD running via HEPfit [50], and obtain constraints an SMEFT coefficients, and a correlation matrix.

## A.1 Including a selection of tau decays

A few $\tau$ decays are included here, to illustrate the differences between LFV involving $\tau s$ and the $\mu \leftrightarrow e$ sector. Firstly, the experimental bounds are more restrictive for muon decays: $\widetilde{BR}(\mu \to eX) \lesssim 10^{-12}$, to be compared with $\widetilde{BR}(\tau \to \ell X) \lesssim 10^{-7}$. So there is sensitivity to smaller coefficients in the muon sector. Furthermore, the loop contributions of some operators involving heavy fermion ($\psi \in \{c, \tau, b\}$) bilinears, can be enhanced by the mass ratio $m_\psi/m_{\mu,\tau}$. This, for instance, enhances the sensitivity to $O_{T,LL}^{cc}$ of $\mu \to e\gamma$ with respect to $\tau \to e\gamma$. Both these effects can be seen in the tables.

An advantage of LFV $\tau$ decays is the multitude of different hadronic final states, which each constrain a specific combinations of quark operator coefficients. This differs from $\mu \to e$ conversion, where many quark operators contribute in interference. There are therefore fewer “flat directions” for $\tau \to e$ than for $\mu \to e$. This advantage is exploited in [49], but not here, where only a few decays are considered. Also, only $\tau \to e$ results are listed; the sensitivities to $\tau \to \mu$ operators can be obtained by rescaling, as given in eq. (A.5).

The limits here are a first attempt to include QED loop effects in some LFV $\tau$ decays, allowing to estimate the sensitivity of these processes to operators that contribute via log-enhanced loops. Previous works have included a wider range of $\tau$ decays, but focus mostly
on operators that contribute at tree level. The prospects for discriminating among operator coefficients by studying asymmetries and angular distributions was studied in [51, 52] for leptonic decays (see also [48, 53]) and using hadronic decays was considered in [47]. Black et al. [54] give the sensitivity of various rare decays to a subset of SMEFT coefficients that contribute at tree level; a more complete study in the SMEFT was performed recently in [49]. The intermediate-state contribution of heavy quark mesons to leptonic LFV $\tau$ decays is considered in [17].

Only the $\Delta F = 1$ decays such as $\tau \to e\bar{e}e$ are included here; the operator basis and anomalous dimension matrices for such processes differ from the $\mu \to e$ case only by permutation of lepton flavour indices. $\Delta F = 2$ decays (eg $\tau \to e\bar{\mu}e$) would require additional operators and a dedicated RGE analysis.

Three ingredients are required to calculate the sensitivities tabulated here: the experimental upper bounds on the BRs, the theoretical formulae for the BRs as a function of a complete set of operator coefficients at the experimental scale, and the matrix which accounts for loop contributions by expressing the coefficients at the experimental scale in terms of coefficients at the weak scale. The experimental BRs are in table 2, and the remainder of this appendix gives theoretical formulae for the BRs. The matrix is obtained by solving the RGEs for the LFV operator coefficients (between $m_W$ and $m_\tau$), and since we restrict to $\Delta F = 1$ operators, these RGEs are the same as for $\mu \to e$ decays and conversion [18, 45] (with some index changes). A dedicated analysis would be required to obtain reliable sensitivities, and “observable-vectors” for $\tau$-LFV.

The $\tau$ can decay hadronically, so it is convenient to rescale its LFV Branching Ratios:

$$ \widetilde{BR}(\tau \to \ell X) \equiv \frac{BR(\tau \to \ell X)}{BR(\tau \to \ell \bar{\nu} \nu)}, \qquad (A.4) $$

where $BR(\tau \to \mu \bar{\nu} \nu) = .174$ and $BR(\tau \to e \bar{\nu} \nu) = .178$ [55]. The experimental bounds on the considered $\widetilde{BR}$s are given in table 2.

The masses of the final state leptons are neglected in the rates. This simplification means that the theoretical BRs and RG evolution are identical for the $\tau \to \mu$ and $\tau \to e$ sectors, after interchanging $\mu$ and $e$ indices. Therefore in the tables, only the four-fermion operators describing $\tau \to e$ transitions are listed; for operators with two lepton indices, the sensitivities to $\tau \to \mu$ coefficients can be obtained by rescaling:

$$ C_{\text{Lor}}^{\mu\tau\ldots} \lesssim C_{\text{Lor}}^{e\tau\ldots} \sqrt{\frac{\widetilde{BR}(\tau \to \mu X)}{\widetilde{BR}(\tau \to e X)}} \qquad (A.5) $$

where ... are the indices corresponding to the final state $X$, and the Lorentz structure subscripts should be identical for both coefficients. ICI In the case of four-lepton operators, all the $\mu$ and $e$ indices should be interchanged, e.g. $C_{\text{Lor}}^{\mu\tau ee} \to C_{\text{Lor}}^{e\tau\mu\mu}$, and so on.

The decays $\tau \to \mu e\bar{e}$ and $\tau \to e\bar{\mu}\mu$ are mediated by the $\tau \to e$ and $\tau \to \mu$ operators considered here, but are not included due to temporary discrepancies in the tensor contribution, between my calculation and [51].

The calculation of $\tau$ decays to mesons is pedagogically introduced in [54], and a careful study considering many final states has recently appeared [49]. The decays considered
<table><caption>**Table 2.** Current bounds on selected $\tau$ lepton flavour violating branching ratios, from References [56, 56–61], normalised to leptonic weak decays as in eq. (A.4).</caption><thead><tr><th>BR&#771;</th><th>current bound</th><th>BR&#771;</th><th>current bound</th></tr></thead><tbody><tr><td>&tau; &rarr; &mu;&gamma;</td><td>2.5 &times; 10<sup>-7</sup></td><td>&tau; &rarr; e&gamma;</td><td>1.9 &times; 10<sup>-7</sup></td></tr><tr><td>&tau; &rarr; &mu;&mu;&#773;&mu;</td><td>1.2 &times; 10<sup>-7</sup></td><td>&tau; &rarr; e&epsilon;&#773;e</td><td>1.5 &times; 10<sup>-7</sup></td></tr><tr><td>&tau; &rarr; &mu;&pi;</td><td>6.3 &times; 10<sup>-7</sup></td><td>&tau; &rarr; e&pi;</td><td>4.5 &times; 10<sup>-7</sup></td></tr><tr><td>&tau; &rarr; &mu;&eta;</td><td>3.7 &times; 10<sup>-7</sup></td><td>&tau; &rarr; e&eta;</td><td>5.2 &times; 10<sup>-7</sup></td></tr><tr><td>&tau; &rarr; &mu;&rho;</td><td>6.9 &times; 10<sup>-8</sup></td><td>&tau; &rarr; e&rho;</td><td>1.0 &times; 10<sup>-7</sup></td></tr></tbody></table>

here are $\tau \to e\{\pi_0, \eta, \rho\}$; the results for $\tau \to \mu\{\pi_0, \eta, \rho\}$ can be obtained as in eq. (A.5). These mesons are interesting because they probe complementary combinations of operator coefficients at tree level.

The decays to $\pi_0$ mesons probe axial vector/pseudoscalar operators, in the isospin=1 combination $u-d$. In the notation of [62], where

$$ \langle 0 | \bar{d} \gamma^{\mu} \gamma_5 u | \pi^+ (P) \rangle = i P^{\mu} \sqrt{2} f_{\pi}, \quad \Gamma(\tau \to \pi \nu) = \frac{G_F^2 f_{\pi}^2 m_{\tau}^3}{8 \pi} $$

with $f_{\pi} \simeq 92.2$ MeV, the Branching Ratio in the presence of (axial) vector operators is

$$ \widetilde{BR}(\tau \to \ell \pi_0) = \frac{3\pi^2 f_{\pi}^2}{m_{\tau}^2} |C_{V,XR}^{uu} - C_{V,XL}^{uu} - C_{V,XR}^{dd} + C_{V,XL}^{dd}|^2, \quad (A.6) $$

because

$$ \langle 0 | J_{A-}^{\mu} | \pi^0 (P) \rangle = i P^{\mu} f_{\pi} \quad (A.7) $$

where $J_{A-}^{\mu} = \frac{1}{2}(\bar{u}\gamma^{\mu}\gamma_5 u - \bar{d}\gamma^{\mu}\gamma_5 d)$, and the coefficient of $2\sqrt{2}G_F(\bar{\ell}\gamma_{\mu}P_X\tau)J_{A-}^{\mu}$ in the La-Grangian is$^5 C^{A-} = \frac{1}{2}(C_{V,XR}^{uu} - C_{V,XL}^{uu} - C_{V,XR}^{dd} + C_{V,XL}^{dd})$. RG mixing vanishes for the axial current, but it is renormalised (when attached to a chiral current) so at $m_W$ this becomes the “constraint”

$$ 2.3 \times 10^{-3} \gtrsim |C_{V,XR}^{uu} - C_{V,XL}^{uu} - C_{V,XR}^{dd} + C_{V,XL}^{dd}| - \frac{\alpha_e}{\pi} (2C_{V,XR}^{uu} + 2C_{V,XL}^{uu} + C_{V,XR}^{dd} + C_{V,XL}^{dd}) |\ln \frac{m_W}{m_{\tau}}| \quad (A.8) $$

The correct constraint, and “observable-vector”(s) for this decay, could be obtained from the expression for the BR in terms of all coefficients that can contribute (including e.g. the pseudoscalar operators). However, eq. (A.8) allows to calculate sensitivities, and see cancellations, such as between $C_{V,XL}^{uu}$ and $C_{V,XL}^{dd}$, due to which the SMEFT operators $O_{EQ}$ and $O_{LQ1}$ do not contribute to $\tau \to \ell\pi_0$ at tree level.

It is interesting to also include LFV $\tau$ decays to the isosinglet $\eta$, because there is a contribution from $s$ quarks, and not a cancellation between the $us$ and $ds$. Still in the notation of [62], with the approximation $f_{\eta} \sim f_{\pi}$ [64, 65], one obtains the contribution

$$ \widetilde{BR}(\tau \to \ell \eta) = \frac{\pi^2 f_{\pi}^2}{m_{\tau}^2} |(C_{V,XR}^{uu} - C_{V,XL}^{uu}) + (C_{V,XR}^{dd} - C_{V,XL}^{dd}) - 2(C_{V,XR}^{ss} - C_{V,XL}^{ss})|^2. \quad (A.9) $$

<footnote>$^5$There is a factor 1/2 is missing in [63].</footnote>
At $m_W$, this becomes

$$4.4 \times 10^{-3} \gtrsim (C_{V,XR}^{uu} - C_{V,XL}^{uu}) + (C_{V,XR}^{dd} - C_{V,XL}^{dd}) - 2(C_{V,XR}^{ss} - C_{V,XL}^{ss}) - \frac{\alpha_e}{\pi}(C_{V,XR}^{dd} + C_{V,XL}^{dd} - 2(C_{V,XR}^{uu} + C_{V,XL}^{uu} + C_{V,XR}^{ss} + C_{V,XL}^{ss})\ln\frac{m_W}{m_\tau})$$

Pseudoscalar operators can also contribute to the decays $\tau \to \ell\pi_0, \ell\eta$. The operator expectation values

$$\langle 0|\frac{1}{2}(\bar{u}\gamma_5 u - \bar{d}\gamma_5 d)|\pi_0\rangle = \frac{f m_{\pi}^2}{(m_u + m_d)} \\ \langle 0|\frac{1}{2\sqrt{3}}(\bar{u}\gamma_5 u + \bar{d}\gamma_5 d - 2\bar{s}\gamma_5 s)|\eta\rangle = \frac{f m_{\pi}^2}{(m_u + m_d)} = \frac{3f_{\eta}m_{\eta}^2}{(m_u + m_d + 4m_s)} \quad (A.10)$$

give a contribution of pseudoscalar coefficients to the Branching Ratios of

$$\widetilde{BR}(\tau \to e\{\pi_0, \eta\}) = 96\pi^2 \left(\frac{m_{\pi_0}}{m_{\tau}}\right)^4 \left(\frac{f_{\pi}}{m_u + m_d}\right)^2 |C_{\eta,\pi}|^2 \quad (A.11)$$

where in the normalisation of eq. (2.1), the coefficients of the operators of eq. (A.10) are $C_\pi = \frac{1}{2}(C_{S,XR}^{uu} - C_{S,XL}^{uu} - C_{S,XR}^{dd} + C_{S,XL}^{dd})$, and $C_\eta = \frac{1}{2\sqrt{3}}(C_{S,XR}^{uu} - C_{S,XL}^{uu} + C_{S,XR}^{dd} - C_{S,XL}^{dd} - 2C_{S,XR}^{ss} + C_{S,XL}^{ss})$. QED loops can mix tensor operators into (pseudo)scalsars, so this will give some sensitivity to the u, d, s tensor operators.

Finally, decays to the vector $\rho$ meson are normalised to $BR(\tau \to \rho\nu)$, assuming $\rho \to \pi\pi$ (as in [63]; see [49] for a more sophisticated solution) and with the usual factor of 2 for the normalisation of neutral and charged particles:

$$\begin{aligned} \widetilde{BR}(\tau \to \ell\rho_0) & \approx \frac{BR(\tau \to \nu\rho)}{BR(\tau \to \ell\nu\bar{\nu})} \frac{\Gamma(\tau \to \ell\rho_0)}{\Gamma(\tau \to \nu\rho)} \\ & \approx 1.43 \frac{|(C_{V,XR}^{uu} + C_{V,XL}^{uu}) - (C_{V,XR}^{dd} + C_{V,XL}^{dd})|^2}{8|V_{ud}|^2} \end{aligned} \quad (A.12)$$

In the second expression, the contribution of the dipole operator (analogous to the dipole contribution to $\tau \to e\bar{e}e$) is neglected, because the current experimental bounds on $\tau \to e\gamma$ and $\tau \to e\rho_0$ are comparable. QED loops mix vector operators of different quark flavour via penguin diagrams, giving this decay some sensitivity to the coefficients at $m_W$ of vector operators with a heavy quark current:

$$\begin{aligned} 7.5 \times 10^{-4} \gtrsim & (C_{V,XR}^{uu} + C_{V,XL}^{uu}) - (C_{V,XR}^{dd} + C_{V,XL}^{dd}) \\ & + \frac{\alpha_e}{3\pi}(2C_{V,XL}^{uu} - 10C_{V,XR}^{uu} - (C_{V,XR}^{dd} - 5C_{V,XL}^{dd}))\ln\frac{m_W}{m_\tau} \\ & + \frac{2\alpha_e}{3\pi}(2C_{V,XL}^{ee} + C_{V,XR}^{ee} + 2C_{V,XL}^{\tau\tau} + C_{V,XR}^{\tau\tau})\ln\frac{m_W}{m_\tau} \\ & + \frac{2\alpha_e}{3\pi}(C_{V,XR}^{\mu\mu} + C_{V,XL}^{\mu\mu} + C_{V,XR}^{ss} + C_{V,XL}^{ss} - 2(C_{V,XR}^{cc} + C_{V,XL}^{cc}))\ln\frac{m_W}{m_\tau} \\ & + \frac{2\alpha_e}{3\pi}(C_{V,XR}^{bb} + C_{V,XL}^{bb}) \end{aligned} \quad (A.13)$$
<table><caption>**Table 3.** Current sensitivities of the processes in the first row to the coefficients, at $m_W$, of QCD×QED-invariant 2- and 4-lepton operators defined in section 2.1. $X, Y \in \{L, R\}, X \neq Y$.</caption><thead><tr><th>coefficient</th><th>μ → eγ</th><th>μ → e<span style="text-decoration: overline">ee</span></th><th>μA → eA</th></tr></thead><tbody><tr><td>|C<sub>D,X</sub>|</td><td>1.12 × 10<sup>-8</sup></td><td>2.21 × 10<sup>-7</sup></td><td>2.35 × 10<sup>-7</sup></td></tr><tr><td>|C<sub>GG,X</sub>|</td><td></td><td></td><td>5.3 × 10<sup>-7</sup></td></tr><tr><td>|C<sup>ee</sup><sub>V,XX</sub>|</td><td>1.10 × 10<sup>-4</sup></td><td>7.80 × 10<sup>-7</sup></td><td>1.86 × 10<sup>-5</sup></td></tr><tr><td>|C<sup>ee</sup><sub>V,XY</sub>|</td><td>2.55 × 10<sup>-4</sup></td><td>9.31 × 10<sup>-7</sup></td><td>3.77 × 10<sup>-5</sup></td></tr><tr><td>|C<sup>ee</sup><sub>S,XX</sub>|</td><td>1.73 × 10<sup>-4</sup></td><td>2.8 × 10<sup>-6</sup></td><td>(3.64 × 10<sup>-3</sup>)</td></tr><tr><td>|C<sup>μμ</sup><sub>V,XX</sub>|</td><td>1.10 × 10<sup>-4</sup></td><td>5.64 × 10<sup>-5</sup></td><td>1.85 × 10<sup>-5</sup></td></tr><tr><td>|C<sup>μμ</sup><sub>V,XY</sub>|</td><td>2.56 × 10<sup>-4</sup></td><td>1.11 × 10<sup>-4</sup></td><td>3.77 × 10<sup>-5</sup></td></tr><tr><td>|C<sup>μμ</sup><sub>S,XX</sub>|</td><td>8.24 × 10<sup>-7</sup></td><td>(1.63 × 10<sup>-5</sup>)</td><td>(1.73 × 10<sup>-5</sup>)</td></tr><tr><td>|C<sup>ττ</sup><sub>V,XX</sub>|</td><td>3.84 × 10<sup>-4</sup></td><td>1.94 × 10<sup>-4</sup></td><td>3.72 × 10<sup>-5</sup></td></tr><tr><td>|C<sup>ττ</sup><sub>V,XY</sub>|</td><td>4.45 × 10<sup>-4</sup></td><td>1.94 × 10<sup>-4</sup></td><td>3.72 × 10<sup>-5</sup></td></tr><tr><td>|C<sup>ττ</sup><sub>S,XX</sub>|</td><td>5.33 × 10<sup>-6</sup></td><td>(1.05 × 10<sup>-4</sup>)</td><td>(1.12 × 10<sup>-4</sup>)</td></tr><tr><td>|C<sup>ττ</sup><sub>S,XY</sub>|</td><td>3.62 × 10<sup>-5</sup></td><td>(7.28 × 10<sup>-4</sup>)</td><td>(7.75 × 10<sup>-4</sup>)</td></tr><tr><td>|C<sup>ττ</sup><sub>T,XX</sub>|</td><td>1.07 × 10<sup>-8</sup></td><td>(2.11 × 10<sup>-7</sup>)</td><td>(2.25 × 10<sup>-7</sup>)</td></tr></tbody></table>
<table><caption>**Table 4.** Current sensitivities of the processes in the first row to the coefficients, evaluated at $m_W$, of QCD×QED-invariant 2-lepton-2quark operators defined in section 2.1. $X, Y \in \{L, R\}, X \neq Y$.</caption><thead><tr><th>coefficient</th><th>μ → eγ</th><th>μ → e<span style="text-decoration: overline">ee</span></th><th>μA → eA</th></tr></thead><tbody><tr><td>|C<sup>dd</sup><sub>V,XX</sub>|</td><td>1.04 × 10<sup>-3</sup></td><td>2.03 × 10<sup>-4</sup></td><td>5.40 × 10<sup>-8</sup></td></tr><tr><td>|C<sup>dd</sup><sub>V,XY</sub>|</td><td>1.64 × 10<sup>-3</sup></td><td>2.01 × 10<sup>-4</sup></td><td>5.30 × 10<sup>-8</sup></td></tr><tr><td>|C<sup>dd</sup><sub>S,XX</sub>|</td><td>5.79 × 10<sup>-3</sup></td><td>(1.14 × 10<sup>-1</sup>)</td><td>1.03 × 10<sup>-8</sup></td></tr><tr><td>|C<sup>dd</sup><sub>S,XY</sub>|</td><td></td><td></td><td>1.03 × 10<sup>-8</sup></td></tr><tr><td>|C<sup>dd</sup><sub>T,XX</sub>|</td><td>5.57 × 10<sup>-6</sup></td><td>(1.10 × 10<sup>-4</sup>)</td><td>1.90 × 10<sup>-7</sup></td></tr><tr><td>|C<sup>uu</sup><sub>V,XX</sub>|</td><td>3.59 × 10<sup>-4</sup></td><td>1.00 × 10<sup>-4</sup></td><td>6.03 × 10<sup>-8</sup></td></tr><tr><td>|C<sup>uu</sup><sub>V,XY</sub>|</td><td>2.87 × 10<sup>-4</sup></td><td>1.02 × 10<sup>-4</sup></td><td>6.25 × 10<sup>-8</sup></td></tr><tr><td>|C<sup>uu</sup><sub>S,XX</sub>|</td><td>3.09 × 10<sup>-3</sup></td><td>(6.71 × 10<sup>-1</sup>)</td><td>1.03 × 10<sup>-8</sup></td></tr><tr><td>|C<sup>uu</sup><sub>S,XY</sub>|</td><td></td><td></td><td>1.03 × 10<sup>-8</sup></td></tr><tr><td>|C<sup>uu</sup><sub>T,XX</sub>|</td><td>5.95 × 10<sup>-6</sup></td><td>(1.17 × 10<sup>-4</sup>)</td><td>9.65 × 10<sup>-8</sup></td></tr><tr><td>|C<sup>ss</sup><sub>V,XX</sub>|</td><td>1.01 × 10<sup>-3</sup></td><td>2.03 × 10<sup>-4</sup></td><td>3.73 × 10<sup>-5</sup></td></tr><tr><td>|C<sup>ss</sup><sub>V,XY</sub>|</td><td>1.64 × 10<sup>-3</sup></td><td>2.01 × 10<sup>-4</sup></td><td>3.73 × 10<sup>-5</sup></td></tr><tr><td>|C<sup>ss</sup><sub>S,XX</sub>|</td><td>2.92 × 10<sup>-4</sup></td><td>(5.77 × 10<sup>-3</sup>)</td><td>2.13 × 10<sup>-7</sup></td></tr><tr><td>|C<sup>ss</sup><sub>S,XY</sub>|</td><td>1.41 × 10<sup>-2</sup></td><td>(2.78 × 10<sup>-1</sup>)</td><td>2.13 × 10<sup>-7</sup></td></tr><tr><td>|C<sup>ss</sup><sub>T,XX</sub>|</td><td>2.82 × 10<sup>-7</sup></td><td>(5.56 × 10<sup>-6</sup>)</td><td>2.33 × 10<sup>-6</sup></td></tr><tr><td>|C<sup>cc</sup><sub>V,XX</sub>|</td><td>3.59 × 10<sup>-4</sup></td><td>8.99 × 10<sup>-5</sup></td><td>1.68 × 10<sup>-5</sup></td></tr><tr><td>|C<sup>cc</sup><sub>V,XY</sub>|</td><td>2.87 × 10<sup>-4</sup></td><td>9.05 × 10<sup>-5</sup></td><td>1.67 × 10<sup>-5</sup></td></tr><tr><td>|C<sup>cc</sup><sub>S,XX</sub>|</td><td>5.23 × 10<sup>-6</sup></td><td>(1.03 × 10<sup>-4</sup>)</td><td>1.83 × 10<sup>-6</sup></td></tr><tr><td>|C<sup>cc</sup><sub>S,XY</sub>|</td><td>2.37 × 10<sup>-5</sup></td><td>(4.68 × 10<sup>-4</sup>)</td><td>1.80 × 10<sup>-6</sup></td></tr><tr><td>|C<sup>cc</sup><sub>T,XX</sub>|</td><td>1.01 × 10<sup>-8</sup></td><td>(1.99 × 10<sup>-7</sup>)</td><td>2.13 × 10<sup>-7</sup></td></tr><tr><td>|C<sup>bb</sup><sub>V,XX</sub>|</td><td>1.32 × 10<sup>-3</sup></td><td>2.56 × 10<sup>-4</sup></td><td>4.71 × 10<sup>-5</sup></td></tr><tr><td>|C<sup>bb</sup><sub>V,XY</sub>|</td><td>2.05 × 10<sup>-3</sup></td><td>2.54 × 10<sup>-4</sup></td><td>4.71 × 10<sup>-5</sup></td></tr><tr><td>|C<sup>bb</sup><sub>S,XX</sub>|</td><td>1.01 × 10<sup>-5</sup></td><td>(1.98 × 10<sup>-4</sup>)</td><td>7.10 × 10<sup>-6</sup></td></tr><tr><td>|C<sup>bb</sup><sub>S,XY</sub>|</td><td>4.04 × 10<sup>-5</sup></td><td>(7.98 × 10<sup>-4</sup>)</td><td>6.92 × 10<sup>-6</sup></td></tr><tr><td>|C<sup>bb</sup><sub>T,XX</sub>|</td><td>7.81 × 10<sup>-9</sup></td><td>(1.52 × 10<sup>-7</sup>)</td><td>1.64 × 10<sup>-7</sup></td></tr></tbody></table>
<table><caption>**Table 5.** Current sensitivities of the processes in the first row to the coefficients of SMEFT operators at $m_W$, added to the Lagrangian as in eq. (A.1). $C_{e\gamma}^{e\mu} = c_W C_{EB}^{e\mu} - s_W C_{EW}^{e\mu}$.</caption><thead><tr><th>coefficient</th><th>μ → eγ</th><th>μ → e<span style="text-decoration: overline">ee</span></th><th>μA → eA</th></tr></thead><tbody><tr><td>C<sup>eμ</sup><sub>eγ</sub>, C<sup>μe</sup><sub>eγ</sub>*</td><td>1.12 × 10<sup>-8</sup></td><td>2.21 × 10<sup>-7</sup></td><td>2.35 × 10<sup>-7</sup></td></tr><tr><td>C<sup>eμ</sup><sub>HE</sub></td><td>1.18 × 10<sup>-5</sup></td><td>1.20 × 10<sup>-6</sup></td><td>1.42 × 10<sup>-7</sup></td></tr><tr><td>C<sup>eμ</sup><sub>HL1</sub></td><td>1.57 × 10<sup>-5</sup></td><td>1.17 × 10<sup>-6</sup></td><td>1.61 × 10<sup>-7</sup></td></tr><tr><td>C<sup>eμ</sup><sub>HL3</sub></td><td>1.57 × 10<sup>-5</sup></td><td>1.17 × 10<sup>-6</sup></td><td>1.61 × 10<sup>-7</sup></td></tr><tr><td>C<sup>eμ</sup><sub>EH</sub></td><td>7.52 × 10<sup>-7</sup></td><td>1.48 × 10<sup>-5</sup></td><td>5.36 × 10<sup>-5</sup></td></tr></tbody></table>

<table><caption>**Table 6.** Current sensitivities of the processes in the first row to the coefficients of SMEFT operators at $m_W$, added to the Lagrangian as in eq. (A.1). $\sum C_{LL}^{\zeta}$ and $\sum C_{RR}^{\zeta}$ are defined in eq. (A.3).</caption><thead><tr><th>coefficient</th><th>μ → eγ</th><th>μ → e<span style="text-decoration: overline">ee</span></th><th>μA → eA</th></tr></thead><tbody><tr><td>ΣC<sup>eμee</sup><sub>EE</sub>, ΣC<sup>eμee</sup><sub>LL</sub></td><td>1.10 × 10<sup>-4</sup></td><td>7.87 × 10<sup>-7</sup></td><td>1.85 × 10<sup>-5</sup></td></tr><tr><td>C<sup>eeμ</sup><sub>LE</sub>, C<sup>eμee</sup><sub>LE</sub></td><td>2.55 × 10<sup>-4</sup></td><td>9.31 × 10<sup>-7</sup></td><td>3.77 × 10<sup>-5</sup></td></tr><tr><td>ΣC<sup>eμμμ</sup><sub>EE</sub>, ΣC<sup>eμμμ</sup><sub>LL</sub></td><td>1.10 × 10<sup>-4</sup></td><td>5.67 × 10<sup>-5</sup></td><td>1.85 × 10<sup>-5</sup></td></tr><tr><td>C<sup>μμeμ</sup><sub>LE</sub>, C<sup>eμμμ</sup><sub>LE</sub></td><td>2.55 × 10<sup>-4</sup></td><td>1.11 × 10<sup>-4</sup></td><td>3.77 × 10<sup>-5</sup></td></tr><tr><td>ΣC<sup>eμττ</sup><sub>EE</sub>, ΣC<sup>eμττ</sup><sub>LL</sub></td><td>3.84 × 10<sup>-4</sup></td><td>1.95 × 10<sup>-4</sup></td><td>3.72 × 10<sup>-5</sup></td></tr><tr><td>ΣC<sup>eττμ</sup><sub>LL</sub></td><td>3.84 × 10<sup>-4</sup></td><td>1.95 × 10<sup>-4</sup></td><td>3.72 × 10<sup>-5</sup></td></tr><tr><td>C<sup>ττeμ</sup><sub>LE</sub>, C<sup>eμττ</sup><sub>LE</sub></td><td>4.40 × 10<sup>-4</sup></td><td>1.91 × 10<sup>-4</sup></td><td>3.75 × 10<sup>-5</sup></td></tr><tr><td>C<sup>eττμ</sup><sub>LE</sub>, C<sup>τμeτ</sup><sub>LE</sub></td><td>1.80 × 10<sup>-5</sup></td><td>3.64 × 10<sup>-4</sup></td><td>3.88 × 10<sup>-4</sup></td></tr></tbody></table>
<table><caption>**Table 7.** Current sensitivities of $\mu \to e\gamma$, $\mu \to e\bar{e}e$, and $\mu A \to eA$ to the coefficients of SMEFT operators at $m_W$.</caption><thead><tr><th>coefficient</th><th>μ → eγ</th><th>μ → e<span style="text-decoration: overline">e</span>e</th><th>μA → eA</th></tr></thead><tbody><tr><td>C<sup>eμuu</sup><sub>EU</sub></td><td>3.59 × 10<sup>-4</sup></td><td>1.00 × 10<sup>-4</sup></td><td>6.03 × 10<sup>-8</sup></td></tr><tr><td>C<sup>eμuu</sup><sub>LU</sub></td><td>2.87 × 10<sup>-4</sup></td><td>1.02 × 10<sup>-4</sup></td><td>6.25 × 10<sup>-8</sup></td></tr><tr><td>C<sup>eμuu</sup><sub>LEQU</sub></td><td>3.09 × 10<sup>-3</sup></td><td>6.71 × 10<sup>-1</sup></td><td>1.03 × 10<sup>-8</sup></td></tr><tr><td>C<sup>eμuu</sup><sub>T,LEQU</sub></td><td>5.95 × 10<sup>-6</sup></td><td>1.17 × 10<sup>-4</sup></td><td>9.65 × 10<sup>-8</sup></td></tr><tr><td>C<sup>eμdd</sup><sub>LQ1</sub></td><td>2.67 × 10<sup>-4</sup></td><td>1.98 × 10<sup>-4</sup></td><td>2.85 × 10<sup>-8</sup></td></tr><tr><td>C<sup>eμdd</sup><sub>LQ3</sub></td><td>5.47 × 10<sup>-4</sup></td><td>6.71 × 10<sup>-5</sup></td><td>5.09 × 10<sup>-7</sup></td></tr><tr><td>C<sup>eμdd</sup><sub>EQ</sub></td><td>2.44 × 10<sup>-4</sup></td><td>2.05 × 10<sup>-4</sup></td><td>2.87 × 10<sup>-8</sup></td></tr><tr><td>C<sup>eμdd</sup><sub>ED</sub></td><td>1.04 × 10<sup>-3</sup></td><td>2.03 × 10<sup>-4</sup></td><td>5.40 × 10<sup>-8</sup></td></tr><tr><td>C<sup>eμdd</sup><sub>LD</sub></td><td>1.64 × 10<sup>-3</sup></td><td>2.01 × 10<sup>-4</sup></td><td>5.30 × 10<sup>-8</sup></td></tr><tr><td>C<sup>eμdd</sup><sub>LEDQ</sub></td><td></td><td></td><td>1.01 × 10<sup>-8</sup></td></tr><tr><td>C<sup>eμcc</sup><sub>EU</sub></td><td>3.59 × 10<sup>-4</sup></td><td>8.99 × 10<sup>-5</sup></td><td>1.68 × 10<sup>-5</sup></td></tr><tr><td>C<sup>eμcc</sup><sub>LU</sub></td><td>2.87 × 10<sup>-4</sup></td><td>9.10 × 10<sup>-5</sup></td><td>1.67 × 10<sup>-5</sup></td></tr><tr><td>C<sup>eμcc</sup><sub>LEQU</sub></td><td>5.23 × 10<sup>-5</sup></td><td>1.03 × 10<sup>-4</sup></td><td>1.83 × 10<sup>-6</sup></td></tr><tr><td>C<sup>eμcc</sup><sub>T,LEQU</sub></td><td>1.01 × 10<sup>-8</sup></td><td>1.99 × 10<sup>-7</sup></td><td>2.09 × 10<sup>-7</sup></td></tr><tr><td>C<sup>eμss</sup><sub>LQ1</sub></td><td>2.67 × 10<sup>-4</sup></td><td>1.60 × 10<sup>-4</sup></td><td>3.02 × 10<sup>-5</sup></td></tr><tr><td>C<sup>eμss</sup><sub>LQ3</sub></td><td>5.47 × 10<sup>-4</sup></td><td>6.20 × 10<sup>-5</sup></td><td>1.15 × 10<sup>-5</sup></td></tr><tr><td>C<sup>eμss</sup><sub>EQ</sub></td><td>2.44 × 10<sup>-4</sup></td><td>1.64 × 10<sup>-4</sup></td><td>2.98 × 10<sup>-5</sup></td></tr><tr><td>C<sup>eμss</sup><sub>ED</sub></td><td>1.04 × 10<sup>-3</sup></td><td>2.03 × 10<sup>-4</sup></td><td>3.73 × 10<sup>-5</sup></td></tr><tr><td>C<sup>eμss</sup><sub>LD</sub></td><td>1.64 × 10<sup>-3</sup></td><td>2.01 × 10<sup>-4</sup></td><td>3.73 × 10<sup>-5</sup></td></tr><tr><td>C<sup>eμss</sup><sub>LEDQ</sub></td><td>1.41 × 10<sup>-2</sup></td><td>2.78 × 10<sup>-2</sup></td><td>2.09 × 10<sup>-7</sup></td></tr><tr><td>C<sup>eμbb</sup><sub>ED</sub>, C<sup>eμbb</sup><sub>LQ1</sub></td><td>1.32 × 10<sup>-3</sup></td><td>2.56 × 10<sup>-4</sup></td><td>4.71 × 10<sup>-5</sup></td></tr><tr><td>C<sup>eμbb</sup><sub>LQ3</sub></td><td>1.32 × 10<sup>-3</sup></td><td>2.56 × 10<sup>-4</sup></td><td>4.71 × 10<sup>-5</sup></td></tr><tr><td>C<sup>eμbb</sup><sub>LD</sub>, C<sup>eμbb</sup><sub>EQ</sub></td><td>2.07 × 10<sup>-3</sup></td><td>2.51 × 10<sup>-4</sup></td><td>4.72 × 10<sup>-5</sup></td></tr><tr><td>C<sup>eμbb</sup><sub>LEDQ</sub></td><td>4.04 × 10<sup>-5</sup></td><td>7.98 × 10<sup>-4</sup></td><td>6.92 × 10<sup>-6</sup></td></tr></tbody></table>
<table><caption>**Table 8.** Current sensitivities of the processes in the first row, to the coefficients $C_{\text{Lor}}^{e\tau\dots}$, evaluated at $m_W$, of $\tau \leftrightarrow e$ flavour-changing, QCD$\times$QED-invariant 2- and 4-lepton operators, defined as in section 2.1 with $\mu \to \tau$. $X, Y \in \{L, R\}, X \neq Y$.</caption><thead><tr><th>coefficient</th><th>&tau; &rarr; e&gamma;</th><th>&tau; &rarr; e&macr;e&macr;</th><th>&tau; &rarr; e&rho;</th></tr></thead><tbody><tr><td>|C<sub>D,X</sub>|</td><td>7.35 &times; 10<sup>-6</sup></td><td>6.36 &times; 10<sup>-5</sup></td><td>...</td></tr><tr><td>|C<sup>ee</sup><sub>V,XX</sub>|</td><td>1.29 &times; 10<sup>-1</sup></td><td>2.83 &times; 10<sup>-4</sup></td><td>6.55 &times; 10<sup>-2</sup></td></tr><tr><td>|C<sup>ee</sup><sub>V,XY</sub>|</td><td>3.00 &times; 10<sup>-1</sup></td><td>3.78 &times; 10<sup>-4</sup></td><td>1.31 &times; 10<sup>-1</sup></td></tr><tr><td>|C<sup>ee</sup><sub>S,XX</sub>|</td><td>3.37</td><td>(1.06 &times; 10<sup>-3</sup>)</td><td></td></tr><tr><td>|C<sup>&mu;&mu;</sup><sub>V,XX</sub>|</td><td>2.59 &times; 10<sup>-1</sup></td><td>7.65 &times; 10<sup>-2</sup></td><td>1.31 &times; 10<sup>-1</sup></td></tr><tr><td>|C<sup>&mu;&mu;</sup><sub>V,XY</sub>|</td><td>3.00 &times; 10<sup>-1</sup></td><td>7.49 &times; 10<sup>-2</sup></td><td>1.31 &times; 10<sup>-1</sup></td></tr><tr><td>|C<sup>&mu;&mu;</sup><sub>S,XX</sub>|</td><td>9.76 &times; 10<sup>-1</sup></td><td>(8.46)</td><td></td></tr><tr><td>|C<sup>&mu;&mu;</sup><sub>S,XY</sub>|</td><td>4.17 &times; 10<sup>-1</sup></td><td>(3.61)</td><td></td></tr><tr><td>|C<sup>&mu;&mu;</sup><sub>T,XX</sub>|</td><td>2.03 &times; 10<sup>-3</sup></td><td>(1.76 &times; 10<sup>-2</sup>)</td><td></td></tr><tr><td>|C<sup>&tau;&tau;</sup><sub>V,XX</sub>|</td><td>1.29 &times; 10<sup>-1</sup></td><td>3.82 &times; 10<sup>-2</sup></td><td>6.55 &times; 10<sup>-2</sup></td></tr><tr><td>|C<sup>&tau;&tau;</sup><sub>V,XY</sub>|</td><td>3.00 &times; 10<sup>-1</sup></td><td>7.49 &times; 10<sup>-2</sup></td><td>1.31 &times; 10<sup>-1</sup></td></tr><tr><td>|C<sup>&tau;&tau;</sup><sub>S,XX</sub>|</td><td>9.68 &times; 10<sup>-4</sup></td><td>(8.39 &times; 10<sup>-3</sup>)</td><td></td></tr></tbody></table>
<table><caption>Table 9. Similar to table 8, but for 2-lepton-2-quark operators.</caption><thead><tr><th>coefficient</th><th>τ → eγ</th><th>τ → e<span style="text-decoration: overline">e</span>e</th><th>τ → eρ</th><th>τ → eπ</th><th>τ → eη</th></tr></thead><tbody><tr><td>|C<sup>dd</sup><sub>V,XX</sub>|</td><td>7.03 × 10<sup>-1</sup></td><td>7.84 × 10<sup>-2</sup></td><td>7.61 × 10<sup>-4</sup></td><td>2.32 × 10<sup>-3</sup></td><td>4.44 × 10<sup>-3</sup></td></tr><tr><td>|C<sup>dd</sup><sub>V,XY</sub>|</td><td>1.10</td><td>7.79 × 10<sup>-2</sup></td><td>7.48 × 10<sup>-4</sup></td><td>2.28 × 10<sup>-3</sup></td><td>4.36 × 10<sup>-3</sup></td></tr><tr><td>|C<sup>dd</sup><sub>S,XX</sub>|</td><td>66.0</td><td>(572)</td><td></td><td>3.74 × 10<sup>-4</sup></td><td>6.96 × 10<sup>-4</sup></td></tr><tr><td>|C<sup>dd</sup><sub>S,XY</sub>|</td><td></td><td></td><td></td><td>3.74 × 10<sup>-4</sup></td><td>6.96 × 10<sup>-4</sup></td></tr><tr><td>|C<sup>dd</sup><sub>T,XX</sub>|</td><td>6.28 × 10<sup>-2</sup></td><td>(5.44 × 10<sup>-1</sup>)</td><td></td><td>6.92 × 10<sup>-3</sup></td><td>1.29 × 10<sup>-2</sup></td></tr><tr><td>|C<sup>uu</sup><sub>V,XX</sub>|</td><td>2.42 × 10<sup>-1</sup></td><td>3.88 × 10<sup>-2</sup></td><td>7.46 × 10<sup>-4</sup></td><td>2.26 × 10<sup>-3</sup></td><td>4.33 × 10<sup>-3</sup></td></tr><tr><td>|C<sup>uu</sup><sub>V,XY</sub>|</td><td>1.94 × 10<sup>-1</sup></td><td>3.93 × 10<sup>-2</sup></td><td>7.72 × 10<sup>-4</sup></td><td>2.34 × 10<sup>-3</sup></td><td>4.48 × 10<sup>-3</sup></td></tr><tr><td>|C<sup>uu</sup><sub>S,XX</sub>|</td><td>35.3</td><td>(305)</td><td></td><td>3.73 × 10<sup>-4</sup></td><td>6.93 × 10<sup>-4</sup></td></tr><tr><td>|C<sup>uu</sup><sub>S,XY</sub>|</td><td></td><td></td><td></td><td>3.73 × 10<sup>-4</sup></td><td>6.93 × 10<sup>-4</sup></td></tr><tr><td>|C<sup>uu</sup><sub>T,XX</sub>|</td><td>6.71 × 10<sup>-2</sup></td><td>(5.81 × 10<sup>-1</sup>)</td><td></td><td>3.46 × 10<sup>-3</sup></td><td>6.44 × 10<sup>-3</sup></td></tr><tr><td>|C<sup>ss</sup><sub>V,XX</sub>|</td><td>7.03 × 10<sup>-1</sup></td><td>7.84 × 10<sup>-2</sup></td><td>1.31 × 10<sup>-1</sup></td><td></td><td>2.22 × 10<sup>-3</sup></td></tr><tr><td>|C<sup>ss</sup><sub>V,XY</sub>|</td><td>1.11</td><td>7.79 × 10<sup>-2</sup></td><td>1.31 × 10<sup>-1</sup></td><td></td><td>2.18 × 10<sup>-3</sup></td></tr><tr><td>|C<sup>ss</sup><sub>S,XX</sub>|</td><td>3.23</td><td>(28.0)</td><td></td><td></td><td>3.48 × 10<sup>-4</sup></td></tr><tr><td>|C<sup>ss</sup><sub>S,XY</sub>|</td><td>15.9</td><td>(138)</td><td></td><td></td><td>3.48 × 10<sup>-4</sup></td></tr><tr><td>|C<sup>ss</sup><sub>T,XX</sub>|</td><td>3.08 × 10<sup>-3</sup></td><td>(2.66 × 10<sup>-2</sup>)</td><td></td><td></td><td>6.44 × 10<sup>-3</sup></td></tr><tr><td>|C<sup>cc</sup><sub>V,XX</sub>|</td><td>2.42 × 10<sup>-1</sup></td><td>3.46 × 10<sup>-2</sup></td><td>5.84 × 10<sup>-2</sup></td><td></td><td></td></tr><tr><td>|C<sup>cc</sup><sub>V,XY</sub>|</td><td>1.94 × 10<sup>-1</sup></td><td>3.50 × 10<sup>-2</sup></td><td>5.84 × 10<sup>-2</sup></td><td></td><td></td></tr><tr><td>|C<sup>cc</sup><sub>S,XX</sub>|</td><td>6.11 × 10<sup>-2</sup></td><td>(5.29 × 10<sup>-1</sup>)</td><td></td><td></td><td></td></tr><tr><td>|C<sup>cc</sup><sub>S,XY</sub>|</td><td>2.68 × 10<sup>-1</sup></td><td>(2.32)</td><td></td><td></td><td></td></tr><tr><td>|C<sup>cc</sup><sub>T,XX</sub>|</td><td>1.06 × 10<sup>-4</sup></td><td>(1.01 × 10<sup>-3</sup>)</td><td></td><td></td><td></td></tr><tr><td>|C<sup>bb</sup><sub>V,XX</sub>|</td><td>8.80 × 10<sup>-1</sup></td><td>9.80 × 10<sup>-2</sup></td><td>1.64 × 10<sup>-1</sup></td><td></td><td></td></tr><tr><td>|C<sup>bb</sup><sub>V,XY</sub>|</td><td>1.38</td><td>9.73 × 10<sup>-2</sup></td><td>1.64 × 10<sup>-1</sup></td><td></td><td></td></tr><tr><td>|C<sup>bb</sup><sub>S,XX</sub>|</td><td>(1.16 × 10<sup>-1</sup>)</td><td>(1.00)</td><td></td><td></td><td></td></tr><tr><td>|C<sup>bb</sup><sub>S,XY</sub>|</td><td>4.57 × 10<sup>-1</sup></td><td>3.96</td><td></td><td></td><td></td></tr><tr><td>|C<sup>bb</sup><sub>T,XX</sub>|</td><td>8.82 × 10<sup>-5</sup></td><td>(7.65)</td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>**Table 10.** Current sensitivities of selected hadronic $\tau$ decays to the coefficients evaluated at $m_W$ of $2\ell2q$ SMEFT operators.</caption><thead><tr><th>coefficient</th><th>τ → eγ</th><th>τ → eρ</th><th>τ → eπ</th><th>τ → eη</th></tr></thead><tbody><tr><td>C<sup>erdd</sup><sub>LQ1</sub></td><td>1.80 × 10<sup>-1</sup></td><td>3.74 × 10<sup>-2</sup></td><td>8.93 × 10<sup>-2</sup></td><td>2.19 × 10<sup>-3</sup></td></tr><tr><td>C<sup>erdd</sup><sub>LQ3</sub></td><td>3.69 × 10<sup>-1</sup></td><td>3.77 × 10<sup>-4</sup></td><td>1.14 × 10<sup>-3</sup></td><td>1.71 × 10<sup>-1</sup></td></tr><tr><td>C<sup>erdd</sup><sub>LD</sub></td><td>1.10</td><td>7.48 × 10<sup>-4</sup></td><td>2.28 × 10<sup>-3</sup></td><td>4.36 × 10<sup>-3</sup></td></tr><tr><td>C<sup>erdd</sup><sub>EQ</sub></td><td>1.65 × 10<sup>-1</sup></td><td>2.38 × 10<sup>-2</sup></td><td>8.93 × 10<sup>-2</sup></td><td>2.21 × 10<sup>-3</sup></td></tr><tr><td>C<sup>erdd</sup><sub>ED</sub></td><td>7.04 × 10<sup>-1</sup></td><td>7.61 × 10<sup>-4</sup></td><td>2.32 × 10<sup>-3</sup></td><td>4.44 × 10<sup>-3</sup></td></tr><tr><td>C<sup>erdd</sup><sub>LEDQ</sub></td><td></td><td></td><td>3.74 × 10<sup>-4</sup></td><td>6.96 × 10<sup>-4</sup></td></tr><tr><td>C<sup>eruu</sup><sub>EU</sub></td><td>2.42 × 10<sup>-1</sup></td><td>7.46 × 10<sup>-4</sup></td><td>2.26 × 10<sup>-3</sup></td><td>4.33 × 10<sup>-3</sup></td></tr><tr><td>C<sup>eruu</sup><sub>LU</sub></td><td>1.94 × 10<sup>-1</sup></td><td>7.72 × 10<sup>-4</sup></td><td>2.34 × 10<sup>-3</sup></td><td>4.48 × 10<sup>-3</sup></td></tr><tr><td>C<sup>eruu</sup><sub>LEQU</sub></td><td>35.3</td><td></td><td>3.73 × 10<sup>-4</sup></td><td>6.93 × 10<sup>-4</sup></td></tr><tr><td>C<sup>eruu</sup><sub>T,LEQU</sub></td><td>6.71 × 10<sup>-2</sup></td><td></td><td>3.46 × 10<sup>-3</sup></td><td>6.44 × 10<sup>-3</sup></td></tr><tr><td>C<sup>erss</sup><sub>LQ1</sub></td><td>1.80 × 10<sup>-1</sup></td><td>1.05 × 10<sup>-1</sup></td><td></td><td>2.22 × 10<sup>-3</sup></td></tr><tr><td>C<sup>erss</sup><sub>LQ3</sub></td><td>3.69 × 10<sup>-1</sup></td><td>4.04 × 10<sup>-2</sup></td><td></td><td>2.22 × 10<sup>-3</sup></td></tr><tr><td>C<sup>erss</sup><sub>LD</sub></td><td>1.06</td><td>1.31 × 10<sup>-1</sup></td><td></td><td>2.18 × 10<sup>-3</sup></td></tr><tr><td>C<sup>erss</sup><sub>EQ</sub></td><td>1.65 × 10<sup>-1</sup></td><td>1.05 × 10<sup>-1</sup></td><td></td><td>2.18 × 10<sup>-3</sup></td></tr><tr><td>C<sup>erss</sup><sub>ED</sub></td><td>7.04 × 10<sup>-1</sup></td><td>1.31 × 10<sup>-1</sup></td><td></td><td>2.22 × 10<sup>-3</sup></td></tr><tr><td>C<sup>erss</sup><sub>LEDQ</sub></td><td>16</td><td></td><td></td><td>3.48 × 10<sup>-4</sup></td></tr><tr><td>C<sup>ercc</sup><sub>EU</sub></td><td>2.42 × 10<sup>-1</sup></td><td>5.84 × 10<sup>-2</sup></td><td></td><td></td></tr><tr><td>C<sup>ercc</sup><sub>LU</sub></td><td>1.94 × 10<sup>-1</sup></td><td>5.84 × 10<sup>-2</sup></td><td></td><td></td></tr><tr><td>C<sup>ercc</sup><sub>LEQU</sub></td><td>6.11 × 10<sup>-2</sup></td><td></td><td></td><td></td></tr><tr><td>C<sup>ercc</sup><sub>T,LEQU</sub></td><td>1.16 × 10<sup>-4</sup></td><td></td><td></td><td></td></tr><tr><td>C<sup>erbb</sup><sub>LQ1</sub></td><td>8.79 × 10<sup>-1</sup></td><td>1.64 × 10<sup>-1</sup></td><td></td><td></td></tr><tr><td>C<sup>erbb</sup><sub>LQ3</sub></td><td>8.79 × 10<sup>-1</sup></td><td>1.32 × 10<sup>-3</sup></td><td></td><td></td></tr><tr><td>C<sup>erbb</sup><sub>LD</sub></td><td>1.38</td><td>1.64 × 10<sup>-1</sup></td><td></td><td></td></tr><tr><td>C<sup>erbb</sup><sub>EQ</sub></td><td>8.79 × 10<sup>-1</sup></td><td>1.64 × 10<sup>-1</sup></td><td></td><td></td></tr><tr><td>C<sup>erbb</sup><sub>ED</sub></td><td>8.79 × 10<sup>-1</sup></td><td>1.64 × 10<sup>-1</sup></td><td></td><td></td></tr><tr><td>C<sup>erbb</sup><sub>LEDQ</sub></td><td>4.57 × 10<sup>-1</sup></td><td></td><td></td><td></td></tr></tbody></table>
Acknowledgments

I thank B. Echenard for encouraging interest and useful plots.

**Open Access.** This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.

References

[1] MEG collaboration, *Search for the lepton flavour violating decay $\mu^+ \to e^+\gamma$ with the full dataset of the MEG experiment*, Eur. Phys. J. C 76 (2016) 434 [arXiv:1605.05081] [INSPIRE].

[2] SINDRUM collaboration, *Search for the Decay $\mu^+ \to e^+e^+e^-$*, Nucl. Phys. B 299 (1988) 1 [INSPIRE].

[3] SINDRUM II collaboration, *A search for muon to electron conversion in muonic gold*, Eur. Phys. J. C 47 (2006) 337 [INSPIRE].

[4] SINDRUM II collaboration, *Test of lepton flavor conservation in $\mu \to e$ conversion on titanium*, Phys. Lett. B 317 (1993) 631 [INSPIRE].

[5] A. Blondel et al., *Research Proposal for an Experiment to Search for the Decay $\mu \to eee$*, arXiv:1301.6113 [INSPIRE].

[6] COMET collaboration, *Conceptual design report for experimental search for lepton flavor violating $\mu^- e^-$ conversion at sensitivity of $10^{-16}$ with a slow-extracted bunched proton beam (COMET)*, KEK-2009-10 [INSPIRE].

[7] COMET collaboration, *Overview of the COMET Phase-I experiment*, PoS FPCP2015 (2015) 059 [INSPIRE].

[8] MU2E collaboration, *Proposal to search for $\mu^- N \to e^- N$ with a single event sensitivity below $10^{-16}$*, FERMILAB-PROPOSAL-0973, [INSPIRE].

[9] MEG II collaboration, *The design of the MEG II experiment*, Eur. Phys. J. C 78 (2018) 380 [arXiv:1801.04688] [INSPIRE].

[10] PRISM collaboration, *An Experimental Search for a $\mu N \to eN$ Conversion at Sensitivity of the Order of $10^{-18}$ with a Highly Intense Muon Source: PRISM*, unpublished, J-PARC LOI, (2006).

[11] S. Davidson, Y. Kuno and M. Yamanaka, *Selecting $\mu \to e$ conversion targets to distinguish lepton flavour-changing operators*, Phys. Lett. B 790 (2019) 380 [arXiv:1810.01884] [INSPIRE].

[12] CLEO collaboration, *Search for Lepton Flavor Violation in Upsilon Decays*, Phys. Rev. Lett. 101 (2008) 201601 [arXiv:0807.2695] [INSPIRE].

[13] S. Davidson, D.C. Bailey and B.A. Campbell, *Model independent constraints on leptoquarks from rare processes*, Z. Phys. C 61 (1994) 613 [hep-ph/9309310] [INSPIRE].

[14] F. Cuypers and S. Davidson, *Bileptons: Present limits and future prospects*, Eur. Phys. J. C 2 (1998) 503 [hep-ph/9609487] [INSPIRE].
[15] M. Carpentier and S. Davidson, *Constraints on two-lepton, two quark operators*, Eur. Phys. J. C 70 (2010) 1071 [arXiv:1008.0280] [INSPIRE].

[16] A. Crivellin, S. Najjari and J. Rosiek, *Lepton Flavor Violation in the Standard Model with general Dimension-Six Operators*, JHEP 04 (2014) 167 [arXiv:1312.0634] [INSPIRE].

[17] C.O. Dib, T. Gutsche, S.G. Kovalenko, V.E. Lyubovitskij and I. Schmidt, *Bounds on lepton flavor violating physics and decays of neutral mesons from $\tau(\mu) \to 3\ell$, $\ell\gamma\gamma$-decays*, Phys. Rev. D 99 (2019) 035020 [arXiv:1812.02638] [INSPIRE].

[18] A. Crivellin, S. Davidson, G.M. Pruna and A. Signer, *Renormalisation-group improved analysis of $\mu \to e$ processes in a systematic effective-field-theory approach*, JHEP 05 (2017) 117 [arXiv:1702.03020] [INSPIRE].

[19] V. Cirigliano, R. Kitano, Y. Okada and P. Tuzon, *On the model discriminating power of $\mu \to e$ conversion in nuclei*, Phys. Rev. D 80 (2009) 013002 [arXiv:0904.0957] [INSPIRE].

[20] J.D. Bowman, T.P. Cheng, L.-F. Li and H.S. Matis, *New Upper Limit for $\mu \to e\gamma\gamma$, Phys. Rev. Lett. 41 (1978) 442 [INSPIRE].

[21] R.D. Bolton et al., *Search for Rare Muon Decays with the Crystal Box Detector*, Phys. Rev. D 38 (1988) 2077 [INSPIRE].

[22] S. Davidson, Y. Kuno, Y. Uesaka and M. Yamanaka, *Probing $\mu e\gamma\gamma$ contact interactions with $\mu \to e$ conversion*, Phys. Rev. D 102 (2020) 115043 [arXiv:2007.09612] [INSPIRE].

[23] Y. Kuno and Y. Okada, *Muon decay and physics beyond the standard model*, Rev. Mod. Phys. 73 (2001) 151 [hep-ph/9909265] [INSPIRE].

[24] Y. Okada, K.-i. Okumura and Y. Shimizu, $\mu \to e\gamma$ and $\mu \to 3e$ processes with polarized muons and supersymmetric grand unified theories, Phys. Rev. D 61 (2000) 094001 [hep-ph/9906446] [INSPIRE].

[25] Y. Okada, K.-i. Okumura and Y. Shimizu, *CP violation in the $\mu \to 3e$ process and supersymmetric grand unified theory*, Phys. Rev. D 58 (1998) 051901 [hep-ph/9708446] [INSPIRE].

[26] R. Kitano, M. Koike and Y. Okada, *Detailed calculation of lepton flavor violating muon electron conversion rate for various nuclei*, Phys. Rev. D 66 (2002) 096002 [Erratum ibid. 76 (2007) 059902] [hep-ph/0203110] [INSPIRE].

[27] V. Cirigliano, S. Davidson and Y. Kuno, *Spin-dependent $\mu \to e$ conversion*, Phys. Lett. B 771 (2017) 242 [arXiv:1703.02057] [INSPIRE].

[28] S. Davidson, Y. Kuno and A. Saporta, *“Spin-dependent” $\mu \to e$ conversion on light nuclei*, Eur. Phys. J. C 78 (2018) 109 [arXiv:1710.06787] [INSPIRE].

[29] G. Bélanger, F. Boudjema, A. Pukhov and A. Semenov, *Dark matter direct detection rate in a generic model with MicrOMEGAs 2.2*, Comput. Phys. Commun. 180 (2009) 747 [arXiv:0803.2360] [INSPIRE].

[30] HERMES collaboration, *Precise determination of the spin structure function $g_1$ of the proton, deuteron and neutron*, Phys. Rev. D 75 (2007) 012007 [hep-ex/0609039] [INSPIRE].

[31] J. Green et al., *Up, down, and strange nucleon axial form factors from lattice QCD*, Phys. Rev. D 95 (2017) 114502 [arXiv:1703.06703] [INSPIRE].

[32] S. Dürr et al., *Lattice computation of the nucleon scalar quark contents at the physical point*, Phys. Rev. Lett. 116 (2016) 172001 [arXiv:1510.08013] [INSPIRE].
[33] M. Hoferichter, J. Ruiz de Elvira, B. Kubis and U.-G. Meißner, *High-Precision Determination of the Pion-Nucleon σ Term from Roy-Steiner Equations*, Phys. Rev. Lett. 115 (2015) 092301 [arXiv:1506.04142] [INSPIRE].

[34] J.M. Alarcon, J. Martin Camalich and J.A. Oller, *The chiral representation of the πN scattering amplitude and the pion-nucleon sigma term*, Phys. Rev. D 85 (2012) 051503 [arXiv:1110.3797] [INSPIRE].

[35] P. Junnarkar and A. Walker-Loud, *Scalar strange content of the nucleon from lattice QCD*, Phys. Rev. D 87 (2013) 114510 [arXiv:1301.1114] [INSPIRE].

[36] M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, *Remarks on Higgs Boson Interactions with Nucleons*, Phys. Lett. B 78 (1978) 443 [INSPIRE].

[37] H.-Y. Cheng, *Low-energy Interactions of Scalar and Pseudoscalar Higgs Bosons With Baryons*, Phys. Lett. B 219 (1989) 347 [INSPIRE].

[38] T. Bhattacharya, V. Cirigliano, R. Gupta, H.-W. Lin and B. Yoon, *Neutron Electric Dipole Moment and Tensor Charges from Lattice QCD*, Phys. Rev. Lett. 115 (2015) 212002 [arXiv:1506.04196] [INSPIRE].

[39] PARTICLE DATA GROUP collaboration, *Review of Particle Physics*, Chin. Phys. C 38 (2014) 090001 [INSPIRE].

[40] A.J. Buras, *Weak Hamiltonian, CP-violation and rare decays*, hep-ph/9806471 [INSPIRE].

[41] M.E. Peskin and D.V. Schroeder, *An introduction to quantum field theory*, Addison-Wesley, Reading, U.S.A. (1995).

[42] M. Ciuchini, E. Franco, L. Reina and L. Silvestrini, *Leading order QCD corrections to $b \to s\gamma$ and $b \to sg$ decays in three regularization schemes*, Nucl. Phys. B 421 (1994) 41 [hep-ph/9311357] [INSPIRE].

[43] W. Dekens and P. Stoffer, *Low-energy effective field theory below the electroweak scale: matching at one loop*, JHEP 10 (2019) 197 [arXiv:1908.05295] [INSPIRE].

[44] S. Davidson, $\mu \to e\gamma$ in the 2HDM: an exercise in EFT, Eur. Phys. J. C 76 (2016) 258 [arXiv:1601.01949] [INSPIRE].

[45] S. Davidson, $\mu \to e\gamma$ and matching at $m_W$, Eur. Phys. J. C 76 (2016) 370 [arXiv:1601.07166] [INSPIRE].

[46] W. Dekens, E.E. Jenkins, A.V. Manohar and P. Stoffer, *Non-perturbative effects in $\mu \to e\gamma$*, JHEP 01 (2019) 088 [arXiv:1810.05675] [INSPIRE].

[47] A. Celis, V. Cirigliano and E. Passemar, *Model-discriminating power of lepton flavor violating $\tau$ decays*, Phys. Rev. D 89 (2014) 095014 [arXiv:1403.5781] [INSPIRE].

[48] M. Giffels, J. Kallarackal, M. Krämer, B. O'Leary and A. Stahl, *The lepton-flavour violating decay $\tau \to \mu\bar{\mu}$ at the CERN LHC*, Phys. Rev. D 77 (2008) 073010 [arXiv:0802.0049] [INSPIRE].

[49] T. Husek, K. Monsalvez-Pozo and J. Portoles, *Lepton-flavour violation in hadronic tau decays and $\mu - \tau$ conversion in nuclei*, JHEP 01 (2021) 059 [arXiv:2009.10428] [INSPIRE].

[50] J. De Blas et al., *HEPfit: a code for the combination of indirect and direct constraints on high energy physics models*, Eur. Phys. J. C 80 (2020) 456 [arXiv:1910.14012] [INSPIRE].

[51] R. Kitano and Y. Okada, *P odd and T odd asymmetries in lepton flavor violating tau decays*, Phys. Rev. D 63 (2001) 113003 [hep-ph/0012040] [INSPIRE].
[52] A. Matsuzaki and A.I. Sanda, *Analysis of lepton flavor violating $\tau^{\pm} \to \mu^{\pm}\mu^{\pm}\mu^{\mp}$ decays*, Phys. Rev. D **77** (2008) 073003 [arXiv:0711.0792] [INSPIRE].

[53] B.M. Dassinger, T. Feldmann, T. Mannel and S. Turczyk, *Model-independent analysis of lepton flavour violating tau decays*, JHEP **10** (2007) 039 [arXiv:0707.0988] [INSPIRE].

[54] D. Black, T. Han, H.-J. He and M. Sher, $\tau - \mu$ flavor violation as a probe of the scale of new physics, Phys. Rev. D **66** (2002) 053002 [hep-ph/0206056] [INSPIRE].

[55] PARTICLE DATA GROUP collaboration, *Review of Particle Physics*, PTEP **2020** (2020) 083C01 [INSPIRE].

[56] BABAR collaboration, *Searches for Lepton Flavor Violation in the Decays $\tau^{\pm} \to e^{\pm}\gamma$ and $\tau^{\pm} \to \mu^{\pm}\gamma$*, Phys. Rev. Lett. **104** (2010) 021802 [arXiv:0908.2381] [INSPIRE].

[57] BELLE collaboration, *New search for $\tau \to e\gamma$ and $\tau \to \mu\gamma$*, decays at Belle, Phys. Lett. B **666** (2008) 16 [arXiv:0705.0650] [INSPIRE].

[58] K. Hayasaka et al., *Search for Lepton Flavor Violating Tau Decays into Three Leptons with $\gamma$19 Million Produced Tau+Tau- Pairs*, Phys. Lett. B **687** (2010) 139 [arXiv:1001.3221] [INSPIRE].

[59] BABAR collaboration, *Search for Lepton Flavor Violating Decays $\tau^{\pm} \to \ell^{\pm}\pi^0$, $\ell^{\pm}\eta$, $\ell^{\pm}\eta'$*, Phys. Rev. Lett. **98** (2007) 061803 [hep-ex/0610067] [INSPIRE].

[60] BELLE collaboration, *Search for lepton flavor violating $\tau^{-}$ decays into $\ell^{-}\eta$, $\ell^{-}\eta'$ and $\ell^{-}\pi^{0}$*, Phys. Lett. B **648** (2007) 341 [hep-ex/0703009] [INSPIRE].

[61] BELLE collaboration, *Search for Lepton-Flavor-Violating tau Decays into a Lepton and a Vector Meson*, Phys. Lett. B **699** (2011) 251 [arXiv:1101.0755] [INSPIRE].

[62] A. Pich, *Effective Field Theory with Nambu-Goldstone Modes*, Les Houches Lect. Notes **108** (2020) [arXiv:1804.05664] [INSPIRE].

[63] S. Davidson and M. Gorbahn, *Charged lepton flavor change and nonstandard neutrino interactions*, Phys. Rev. D **101** (2020) 015010 [arXiv:1909.07406] [INSPIRE].

[64] T. Feldmann, *Quark structure of pseudoscalar mesons*, Int. J. Mod. Phys. A **15** (2000) 159 [hep-ph/9907491] [INSPIRE].

[65] H. Leutwyler, *On the $1/N$ expansion in chiral perturbation theory*, Nucl. Phys. B Proc. Suppl. **64** (1998) 223 [hep-ph/9709408] [INSPIRE].