Lepton flavour violation in composite Higgs models

We discuss in detail the constraints on the partial compositeness coming from flavour and CP violation in the leptonic sector. In the first part we present a formulation of partial compositeness in terms of a flavour symmetry group and a set of spurions, whose background values specify the symmetry breaking pattern. In such a framework we construct the complete set of dimension-six operators describing lepton flavour violation and CP violation. By exploiting the existing bounds, we derive limits on the compositeness scale in different scenarios, characterised by increasing restrictions on the spurion properties. We confirm that in the most general case the compositeness scale should lie well above 10 TeV. However, if in the composite sector the mass parameters and Yukawa couplings are universal, such a bound can be significantly lowered, without necessarily reproducing the case of minimal flavour violation. The most sensitive processes are decays of charged leptons either of radiative type or into three charged leptons, μ→e\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu \rightarrow e$$\end{document} conversion in nuclei and the electric dipole moment of the electron. In the second part we explicitly compute the Wilson coefficients of the relevant dimension-six operators in the so-called two-site model, embodying the symmetry breaking pattern discussed in our first part, and we compare the results with those of the general spurion analysis.


Introduction
If the solution to the gauge hierarchy problem is based on a new symmetry and not on anthropic considerations (for a review see [1]) or on special evolutions of the scalar sector in the early universe [2], new physics at the TeV scale is expected. In most of the existing models new degrees of freedom carrying flavour quantum numbers are present at the a e-mail: feruglio@pd.infn.it b e-mail: paride.paradisi@pd.infn.it c e-mail: pattori@physik.uzh.ch TeV scale, representing potential sources of flavour changing neutral currents (FCNC) and CP violation. So far both direct searches at the LHC and indirect searches in the context of precision tests and flavour physics have brought no conclusive evidence of new physics at the TeV scale. The negative outcome of the search for new physics in the flavour sector is particularly intriguing. Indeed a scale of new physics NP as large as 10 5 TeV [3] is required by an effective operator analysis to preserve the good agreement between observations and theory predictions, unless the new flavour sector is highly non-generic and involves specific mechanisms to suppress FCNC and CP violation at the desired level.
The latter possibility is empirically supported by the huge hierarchies among fermion masses and fermion mixing angles, which can only be explained by some special dynamics. An effective mechanism suppressing FCNC and CP violation can be introduced by observing that in the electroweak theory the symmetry of the flavour sector is broken only by the Yukawa interactions. Minimal flavour violation (MFV) [4] is defined by the assumption that, even including new physics contributions, Yukawa couplings are the only source of such symmetry breaking. In MFV flavour effects from new physics are controlled and damped by the smallness of the fermion masses and mixing angles. In this framework data allow NP to be considerably smaller, close to the TeV scale. MFV provides a useful benchmark for the discussion of the flavour sector, but it does not emerge as a unique framework from the known mechanisms aiming to explain the observed fermion spectrum (for a review see [5]), or from the known models providing a solution to the gauge hierarchy problem.
In this paper we reconsider the possibility that both the origin of fermion masses and the suppression of FCNC and CP violation are due to the mechanism of partial compositeness (PC) [6], as realised in the context of composite Higgs models (for reviews see [7,8]), and we perform a detailed analysis of flavour and CP violations in the leptonic sector. According to PC there are no direct couplings between the elementary fermions and the Higgs doublet. The Higgs dou-blet has potentially strong couplings to a composite sector, including, in the simplest case, a set of vector-like fermions with masses of the order of the compositeness scale. The SM fermions are mostly elementary and get their masses through mixing terms with operators of the composite sector, often modelled by vector-like fermions.
An appealing realisation of this idea involves anarchic Yukawa couplings in the composite sector. In this case the observed hierarchies between SM fermion masses and mixing angles are entirely due to the elementary-composite mixing terms. This is of particular interest, especially for the lepton sector, since the known pattern of neutrino masses and mixing angles as extracted from neutrino oscillation experiments [9] seems to support the idea of an underlying anarchic dynamics [10,11]. Higher-dimensional operators describing low-energy FCNC and CP violations are depleted by both inverse powers of the compositeness scale and by the mixing terms, thus realising an efficient suppression mechanism known as RS-GIM [12]. Quantitative studies in concrete models show that in the anarchic scenario limits from CP violation in the quark sector lead to NP > 10 TeV [8], while the existing bound on the rate of μ → eγ results in NP > 25 TeV [8]. This strongly disfavours the anarchic scenario when the compositeness scale is of the order of 1 TeV. It is also known that PC at the TeV scale can satisfy the bound from flavour physics if Yukawa couplings in the composite sector are non-generic. For instance if we assume that such couplings are universal and, at the same time, that the only irreducible sources of flavour symmetry breaking are proportional to the SM Yukawa couplings, we reproduce exactly the MFV scheme [13].
In this paper we focus on the lepton sector. Several studies of lepton flavour violation (LFV) and CP violation in the lepton sector of composite Higgs models have been realised. Most of the analysis have been performed in the context of five-dimensional (5D) models with a warped spacetime metric, weakly coupled duals to strongly coupled fourdimensional conformal theories, believed to provide a calculable framework for composite Higgs models. Explicit computations with anarchic Yukawa couplings have been carried out in Refs. [14][15][16][17][18][19][20]. They get a lower bound on the masses of the first Kaluza-Klein modes of order 10 TeV. These bounds can be relaxed by requiring discrete [21][22][23][24] or continuous [25][26][27] flavour symmetries. The original motivation for discrete symmetries, tailored to approximately reproduce the tri-bimaximal mixing pattern in the lepton sector, has weaken after the precise measurement of the θ 13 angle and the models considered in [21][22][23][24] require considerable corrections now. Continuous non-abelian symmetries in the composite sector, broken by the elementary-composite mixing terms have been analysed in detail, especially in relation to FCNC and CP violation in the quark sector [28][29][30][31][32]. A review on flavour physics in 5D models with warped space-time metric can be found in Ref. [33]. LFV and CP violation in the lepton sector have been investigated also in SM extensions with extra vector-like heavy leptons [34,35] that can mimic the PC scenario, at least as far as the contributions from heavy fermions is concerned.
In the present work we recast the framework of PC in terms of a generalised flavour symmetry and a suitable set of relevant spurions, much along the lines of Refs. [12,28,32,36]. This involves a certain degree of model dependence, since both the flavour symmetry and the spurions are determined by the specific set of composite leptons that we choose by following the criterion of minimality. In particular we work in the limit of vanishing neutrino masses, turning off the potential effects related to massive neutrinos. It is well known that within PC, LFV and CP violation in the lepton sector are present also in the limit of massless neutrinos. By assuming that the set of adopted spurions are the only irreducible sources of flavour and CP violation, we can construct an exhaustive list of Wilson coefficients related to dimension six operators describing LFV and CP-violating processes. At variance with the previous studies, we include for the first time all Wilson coefficients containing up to four powers of the spurions describing the elementary-composite mixing and we discuss their role in deriving the bounds on the compositeness scale. We also provide a complete list of the LFV Wilson coefficients that can be constructed in the limit of vanishing "wrong" Yukawa couplings. "Wrong" Yukawas in the composite sector are allowed by gauge symmetry, but they do not contribute to SM lepton masses, at the leading order. They directly contribute to the dipole operators describing radiative decays of the charged leptons and setting such Yukawas to zero can relax the bounds on the compositeness scale.
The general scope of our analysis is to check whether there are alternative solutions, beyond MFV, to reconcile PC at the TeV scale with the existing bounds on LFV and CP violation. We also wish to verify if the anarchic scenario is completely ruled out or not. By exploiting the effective Lagrangian of our construction and the existing experimental bounds we estimate the limits on the new physics scale in several scenarios, where our spurions are subjected to a series of increasingly restrictive conditions. We confirm that in general the anarchic scenario is not compatible with new physics at the TeV scale and we provide examples of how PC can be realised at the TeV scale, without necessarily resorting to MFV.
In the second part of our paper we consider as an explicit model realisation of the flavour symmetry and its breaking pattern the so-called two-site model, first introduced in Ref. [37]. Such a realisation contains explicitly vector-like leptons, implementing partial compositness in the lepton sector, as well as a set of spin-one resonances. By integrating out the states at the compositness scale, we evaluate the Wilson coefficients of the relevant LFV and CP-violating dimension six operators and we compare the results with those of the general spurion analysis.
Our paper has the following plan. In Sect. 2 we define the flavour symmetry and the set of spurions of our setup and we characterise the Wilson coefficients of the dimensionsix operators relevant to LFV and CP violation in the lepton sector. In Sect. 3 we perform a phenomenological analysis and we study the bounds on the new physics scale obtained by making different types of assumptions on the available spurions. In Sect. 4 we recall the main aspects of the twosite model, which explicitly incorporates the features of the flavour symmetry breaking defined in Sect. 2. In Sect. 5 we collect our results on the dimension six operators obtained from the model by integrating out heavy fermions and heavy gauge vector bosons. In Sect. 6 we present a phenomenological analysis of LFV in the two-site model. Finally we draw our conclusion. In Appendix B we show the result of our computation of the full one-loop contribution to the electromagnetic dipole operator in the two-site model.

Effective field theory for lepton flavour violation
As a first step, we have to choose the flavour symmetry group of our effective theory and its breaking terms. Throughout this paper we will work in the limit of massless neutrinos. The leptons are those of the SM, that is, three copies of SU(2) doublets and singletsẽ. In MFV the flavour symmetry group of the leptonic sector is SU(3) × SU(3)ẽ, corresponding to independent transformations made on andẽ. In our analysis we will instead assume a PC scenario. Charged leptons have no direct coupling to the Higgs doublet and acquire masses via mixing with vector-like heavy fermions. In this framework it is natural to assume as flavour symmetry group (focussing only on the non-Abelian part): under which the lepton fields rotate in generation space in the following way: where V and Vẽ are elements of SU(3) and SU(3)ẽ, respectively. In other words, the SM leptons only transform under the SU(3) ×SU(3)ẽ component of the flavour group G f , and they are invariant under the remaining SU(3) 4 factor. Such a factor will be used to specify the spurions of our effective field theory. We introduce three sets of spurions. We need spurions ( ,˜ ) that mix SM leptons with heavy fermions. Moreover, we allow for spurions (m,m) describing the masses of the heavy fermions in the limit of unbroken electroweak symmetry. Finally the heavy fermion sector can interact with the Higgs doublet and this implies additional spurions (Y * L , Y * R ). Notice that this set of spurions is the most general one compatible with our flavour group G f and with the assumption that the SM Higgs doublet directly couples only to the heavy sector. From these considerations the following transformation properties for our spurions can be deduced: where it is evident from our notation which SU(3) factors of G f are involved in each transformation. An equivalent set of spurions is obtained by replacing the dimensionful quantities ( ,˜ ) with dimensionless combinations In such a case at the leading order (LO) SM charged lepton masses are generated by the operator where ϕ is the Higgs electroweak doublet. It is well known that, if all the spurions (X,X ), (Y * L , Y * R ) and (m,m) are present at the same time and if the Yukawa couplings (Y * L , Y * R ) are anarchic, a severe bound on the new physics scale applies. Indeed in concrete models belonging to the general class we are considering, the one-loop exchange of Higgs and heavy fermions leads to the following electromagnetic dipole operator: If such a contribution is present, and the Yukawa couplings (Y * L , Y * R ) are assumed to be anarchic of O(1), then the electromagnetic dipole operator and the mass operator are not aligned in flavour space and the heavy fermion scale is bounded to be heavier than about 30 TeV, to respect the bound on BR(μ → eγ ). One way to eliminate this dangerous contribution, while maintaining non-vanishing lepton masses, is to assume Y * L = 0. We will come back to this assumption later in this section. For the moment we will adopt it as a working hypothesis. Our purpose is to analyse the flavourviolating contributions surviving in this limit and to estimate the corresponding bounds on the new physics scale.
At variance with MFV, spurions with the dimension of a mass are present in our setup and some additional prescriptions are needed: • First of all we require that our operators are local in the spurions (X,X ). These are mixing parameters that are generically treated as small and provide one set of expansion parameters for our spurion analysis. • We also assume that the operators are local in the Yukawa couplings Y * R , which we restrict in the range 1 ≤ |Y * R | ≤ 4π .
• Each power of Y * R occurs accompanied either by an Higgs electroweak doublet ϕ or by a factor 1/4π . • Masses of the composite sector are described by the spurions (m,m) to which we add an additional parameter M, singlet under the flavour symmetry group, to describe masses of other composite particles, such as for instance a new set of vector boson resonances. These additional states are coupled to the Higgs doublet and to the heavy fermions with a strong coupling constant g * , in the range 1 ≤ g * ≤ 4π . • In our operators masses will always appear in negative powers, to allow for decoupling of all the operators with the scale of new physics.
With this set of assumptions, the LO operator describing charged lepton masses is still O M in Eq. (5). We stress that this operator provides the definition of the spurion Y * R . Any polynomial of the type With no loss of generality we can redefine as Y * R the particular combination of Eq. (7) occurring in O M .
The physical processes we are interested in concern lepton flavour violation in the charged lepton sector as well as the magnetic dipole moments and the electric dipole moments (EDM) of the charged leptons. To this purpose it is convenient to adopt an effective field theory description where the SM Lagrangian is extended by an appropriate set of gauge invariant operators depending on the SM fields [38]: where we have restricted our attention to the lowestdimensional operators relevant to the processes we are interested in, namely those of dimension six. Dots denote higherdimensional operators. We list below a complete set of dimension six operators depending on lepton fields and on the scalar electroweak doublet ϕ [39]. We start with the dipole operators W I μν and B μν are the field strengths for the gauge vector bosons of SU(2) and U (1), respectively. The flavour structure of these two operators is the same and we will focus on the combinations where θ W is the weak mixing angle and Q eW 3 denotes the contribution to Q eW obtained setting to zero W 1,2 μν . The electromagnetic dipole operator Q eγ is the only operator that gives a tree-level contribution to the radiative decays of charged leptons, when i = j. The diagonal elements, i = j, contribute to the anomalous magnetic moments and to the EDM of the charged leptons. We have operators bilinear in the Higgs doublet After the breaking of the electroweak symmetry these operators modify the couplings of the Z boson to leptons, potentially violating both universality, through the diagonal terms, and lepton flavour, through the non-diagonal ones. There is a unique operator trilinear in the Higgs doublet which contributes, with a different weight, to both the masses and the Higgs couplings of the charged leptons. Finally, we have the four-lepton operators which can contribute to muon and tau decays into three charged leptons. We also consider dimension-six operators of the type llqq that can contribute to μ → e conversion in nuclei: 1 1 The operators (Q (u) q ) i j and (Q (d) q ) i j are linear combinations of the operators (Q (1) q ) i j and (Q (3) q ) i j of Ref. [39].
There are other three independent operators of this type [39], but the chosen subset is sufficiently general for the purposes of the present discussion. Notice that each operator carries flavour indices. Hermiticity of the effective Lagrangian is guaranteed either by appropriate symmetry properties of the Wilson coefficients under transposition of the family indices or by addition of the hermitian conjugate operator.
Our aim is to estimate the Wilson coefficients of these operators, by expressing them in terms of the spurions using the set of rules described above. We expand each Wilson coefficient in powers of the mixings (X,X ) and the Yukawa coupling Y * R . Since the spurions (X,X ) control lepton masses, they are expected to be small, of order (0.1/Y * R ) at most. 2 We will stop the expansion in (X,X ) at the fourth order. In the expansion we will go up to the third order in Y * R , since in the anarchic scenario trilinear combinations of Yukawa couplings are in general misaligned with respect to linear ones. Higher orders in Yukawa couplings do not bring any new qualitative feature in our analysis. Finally the correct dimension is provided by negative powers of M or (m,m). Since the operators under study have dimension six, in practice we have two classes of operators, those suppressed by 1/M 2 and those suppressed by the heavy fermion masses. Each of these two classes refers to a specific decoupling limit. When M |m|, |m|, the heavy bosons decouple first and the operators are suppressed by negative powers of (m,m). We call this the heavy boson (HB) case. In the opposite limit, M |m|, |m|, we have a fast heavy fermion decoupling and the operators are suppressed by the smaller scale M. We call this the heavy fermion (HF) case. When the two scales M and (m,m) are comparable, the Wilson coefficient can be a generic function of the ratio of the two scales. For the present discussion the two limiting cases are sufficient to capture the behaviour of the system. For our spurion analysis it is convenient to rewrite the mass matrices m,m in this way where m 0 is a flavour-independent mass parameter, while the flavour dependence is carried by the dimensionless matrices c,c.
Concerning the dependence of the operators on the strong coupling constant g * , we observe that g * is flavour blind; it only affects the overall strength of the operators and can be absorbed in a redefinition of the NP scale. Therefore we will simply set g * = 1 for the rest of this model-independent analysis. In the specific model discussed in Sect. 4, we will focus on the LO contributions to the relevant operators, aris- Table 1 Spurion combination CL R , in a matrix notation, for the lepton bilinear¯ Li (CL R ) i jẽR j . N Y and N X are the orders of the expansion in Y * R and (X,X ), respectively. We restrict the list to N Y ≤ 3 and N X ≤ 4. For convenience we distinguish spurion combinations depending on composite fermion matrices c andc (column HB) from combinations not involving c andc (column HF) ing either at tree level or at one loop, and we will keep explicit track of such a dependence.
To facilitate the identification of the relevant Wilson coefficients, it is useful to identify the combinations of spurions that fit the lepton bilinears In Tables 1, 2 and 3 we collect such combinations, constructed with the rules outlined above.

Dipole operators
The dipole operators in Eq. (9) involve the lepton bilinear Li (CL R ) i jẽR j and their possible Wilson coefficients are the (CL R ) i j combinations listed in Table 1. We expect that the dipole operators are loop generated in perturbation theory. The naive loop suppression factor 1/(16π 2 ) is not present in Table 1 and should be included in the Wilson coefficient, The new physics scale associated to dipole operators is M in the HF case (M |m|, |m|) and m 0 in the HB one (M |m|, |m|). The first important outcome of our analysis is that there is a large set of potentially lepton flavour-violating combi- Table 2 Spurion combination CL L , in a matrix notation, for the lepton bilinear¯ Li (CL L ) i j L j . N Y and N X are the orders of the expansion in Y * R and (X,X ), respectively. We restrict the list to N Y ≤ 2 and N X ≤ 4. For convenience we distinguish spurion combinations depending on composite fermion matrices c andc (column HB) from combinations not involving c andc (column HF)

Table 3
Spurion combination CR R , in a matrix notation, for the lepton bilinearẽ Ri (CR R ) i jẽR j . N Y and N X are the orders of the expansion in Y * R and (X,X ), respectively. We restrict the list to N Y ≤ 2 and N X ≤ 4. For convenience we distinguish spurion combinations depending on composite fermion matrices c andc (column HB) from combinations not involving c andc (column HF) nations, beyond that of Eq. (6). The actual appearance of these combinations in concrete models containing our set of spurions will depend on the specific dynamics of the model under consideration. We will discuss the phenomenological implications of the new structures for the Wilson coefficients (C eγ,Z ) i j in Sect. 3.

Scalar operator
The scalar operator (Q eϕ ) i j = (ϕ † ϕ)(¯ LiẽR j ϕ) has exactly the same flavour structure of the dipole operators and we can directly read from Table 1 the list of possible Wilson coefficients, provided we have at least three Yukawa couplings, N Y ≥ 3, since this operator contains three Higgs doublets. Up to an overall flavour-independent coefficient, which is expected to be of order one, we have Notice that there is no loop suppression in this case, since this operator can be generated at tree level. As for the dipole operators, the scale of new physics is M in the HF case and m 0 in the HB one. This operator, together with the mass operator of Eq. (5), contributes to both masses and Yukawa couplings of charged leptons. If the Wilson coefficients are not exactly aligned in flavour space, we have flavour-violating decays of the Higgs. Even in the case where there is a perfect alignment among all Wilson coefficients, the overall strength of Yukawa couplings is altered compared to the case of the SM, where it is completely fixed by the fermion mass and by the electroweak VEV. We will discuss the impact of these modifications in Sect. 3.

Vector operators
In the operators (Q (1,3) ϕl ) i j , (Q ϕe ) i j of Eq. (11) new flavour structures arise. The corresponding Wilson coefficients can be read from Table 2, for (Q (1,3) ϕl ) i j and Table 3, for (Q ϕe ) i j , but excluding the case N Y = 0, since the vector operators are bilinear in the Higgs doublets and this requires at least two powers of the Yukawa couplings. We have There is no loop suppression, in general, and the new physics scale is identified as for the previous operators. We see that, in general, the Wilson coefficients for the vector operators are not diagonal in the mass basis. This leads to violation of flavour in the couplings of the Z boson, with consequences that we discuss in Sect. 3.

Contact operators
Finally we consider the contact operators (Q ll ) i jmn , (Q ee ) i jmn and (Q le ) i jmn given in Eq. (13). In these operators we recognise combinations of the flavour structures already discussed. We expect the following possible factorisations: Since the contact terms contain no Higgs doublets, we have no restrictions on N Y . For the same reason, each bilinear in the Yukawa coupling should be accompanied by a loop suppression factor 1/(16π 2 ). 3 Concerning the new physics scale , when the two factors C come from the HF column we have = M. When one the factors comes from the HF column and the other from the HB column, we have = m 0 . Similarly, for the llqq operators we have with no restrictions on N Y .

Stability of the solution Y
Since we are interested in the scenario where Y * L is nearly vanishing, a legitimate question is whether and to which extent such a limit is stable under quantum corrections. To this purpose it is better to distinguish the two regimes v When v m 0 M, we can consider the following spurion combinations: that behave as effective Yukawas of type Y * L since they have the same transformation properties as Y * L . We expect that terms like those in (25) arise, in perturbation theory, through threshold corrections induced by the one-loop exchange of heavy gauge bosons and heavy fermions, from which we estimate where g * is a coupling constant of the strong sector and k is either an additional loop factor 1/16π 2 stemming from a Higgs loop, or a factor coming from the electroweak VEV, v 2 /M 2 , if the two Yukawas are attached to external Higgs legs. The contributions of Eq. (26) are small if there is an hierarchy between heavy fermions and heavy gauge boson masses. They are also suppressed in the semi-perturbative regime g * , Y * R ≈ 1. Similarly, exchange of ordinary gauge bosons and heavy fermions lead to effective Y * L couplings of the type 3 If contact terms originate from dimension eight operators, the loop factor can be effectively replaced by v 2 / 2 .
These contributions are naturally suppressed by the ratio m 2 W,Z /m 2 0 and can easily be kept at the percent level. In the other regime, v M m 0 , beyond the combinations of Eq. (27), we can also consider which also transform as the spurion Y * L . These combinations formally decouple in the limit of infinitely large m 0 . There are, however, also combinations that do not decouple, such as We have checked, through the one-loop computation in the two-site model that we will consider in Sect. 4, that this kind of contributions are generated. They cannot be parametrically suppressed by mass ratios and the corresponding Y * L can be depleted only in the semi-perturbative regime g * , Y * R ≈ 1. Notice that for m 0 ≈ M, the effective Yukawa couplings arising from Eqs. (26) and (28), are insensitive to the overall mass scale. In summary, in the regime of heavy gauge bosons much heavier than heavy fermions, the effective Y * L can remain close to the percent level even in a strongly coupled regime, while in the opposite regime the solution Y * L ≈ 0 is typically stable under quantum corrections (up to the % level) only if we assume the semi-perturbative regime g * , Y * R ≈ 1.

Model-independent bounds on the scale of new physics
Model-independent studies of LFV processes have already been done in the literature [53][54][55]. In this section we collect the bounds on the Wilson coefficients of the Lagrangian of Eq. (8) and discuss their impact on our spurion analysis. The main bounds on the coefficients of the dipole operators O eγ and O eZ are given in Table 5, showing the results obtained in Refs. [54,55] and derived from Ref. [56] in the case of μ − Au → e − Au. They come from the present limits reported in Table 4. The amplitudes for these processes get a tree-level contribution from the off-diagonal elements of the electromagnetic dipole operator O i j eγ , while O i j eZ contributes at one loop. Each bound has been derived by assuming a single non-vanishing Wilson coefficient at the time. This also applies to all the bounds discussed in this section. The bounds on the coef-  [41] μ → 3e 1.0 × 10 −12 [42] ≈10 −16 [43] μ − Au → e − Au 7.0 × 10 −13 [44] ?   [1-loop] ficients C ji eγ,Z are equal to the bounds on the coefficients C i j eγ,Z . The diagonal elements of the dipole operators contribute to electric and magnetic dipole moments of the charged leptons. From the present bounds reported in Table 4 we have [54] Im(C ee eγ ) 1 TeV 2 < 3.9 × 10 −12 , Given the current deviation a μ = a EXP μ − a SM μ in the muon anomalous magnetic moment a μ = (g − 2) μ /2 [57,58] a μ = (29 ± 9) × 10 −10 , we would need to account for the central value of the discrepancy. Also the scalar operator O eϕ is mostly bounded by the limits on radiative lepton decays [54,55,[59][60][61]. The scalar operator contributes to lepton masses and to higgs couplings with a different weight: where and y SM i j are the Standard Model Yukawa couplings. Radiative charged lepton decays constrain the off-diagonal elements of Y i j in the basis where the mass matrix M i j is diagonal. To convert these constraints in bounds on the Wilson coefficients (C eϕ ) i j it is convenient to work in the basis where the SM couplings y SM i j are diagonal and expand the unitary matrices that diagonalise M in powers of v 2 / 2 . We found that, to first order in this parameter, the off-diagonal elements Y i j (i = j) in the lepton mass basis are given by: By translating the bounds on Y i j (i = j) given in Refs. [59][60][61] into bounds on (C eϕ ) i j , we get the results shown in Table  6. The bounds on the coefficients (C eϕ ) ji are equal those on the coefficients (C eϕ ) i j . These bounds are dominated by twoloop contributions of the corresponding operator to the radiative lepton decay, through Barr-Zee type diagrams, assuming that the top Yukawa coupling is as in the SM. One-loop contributions to charged lepton radiative decays and tree-level contributions to → 3 decays lead to less severe bounds than the ones given in Table 6.  [2-loop] Coming to the vector operators (Q (1,3) ϕl ) i j , (Q ϕe ) i j , they lead to lepton flavour-violating Z decays, but the corresponding limits on the Wilson coefficients, assuming = 1 TeV, are of order 10 % [54]. Through one-loop diagrams they also contribute to radiative decays of the charged leptons [54,55]. It turns out that the most restrictive bounds come from the processes μ − Au → e − Au and → 3 whose branching ratios satisfy the experimental limits of Table 4.
We collect the corresponding bounds in Table 7. Also the contact operators can contribute to μ − Au → e − Au, → 3 and, through one-loop diagrams, to the radiative decays of the charged leptons. The most significant bounds are given in Table 8.
We can translate the bounds collected in Tables 5, 6, 7 and 8 into limits on the masses M and m 0 of our spurion analysis. To do this we should determine or make some assumptions on the parameters c,c, X ,X , Y * R and Y * L , through which we can express all the Wilson coefficients, as explained in Sect. 2. By exploiting the flavour symmetry of our setup, we see that it is not restrictive to work in the basis where the mixing matrices X andX are diagonal, real and non-negative and we will adopt this choice, unless otherwise stated. The LO mass matrix of the charged fermions is diagonalised by a bi-unitary transformation Assuming that Y * R is anarchic, it is straightforward to find that where factors of order one have been omitted from the matrix elements. Similarly the lepton masses are approximately given by The most favorable scenario to minimise FCNC effects is realised when X i =X i , which we assume, for the time being.
We have which will be used in our estimates. Table 7 The bounds on off-diagonal Wilson coefficients (C (1,3) ϕl ) i j and C i j ϕe from Refs. [54,55]. The bounds from μ − Au → e − Au have been derived from Ref. [56]. In the second column we list the upper bound on the Wilson coefficients assuming = 1 TeV, while in the third column we set to unity the coefficients and we list the corresponding lower bound on , in TeV If c,c, Y * R and Y * L are all generic matrices, the dipole operator, Eq. (6), leads to the well-known limit where Y denotes a suitable average of the Y * L ,R couplings and similarly for c . Such a bound can be evaded only if special features of the flavour parameters c,c, X ,X , Y * R and Y * L are adopted. One possibility is to recover the framework of minimal flavour violation (MFV), as proposed in Ref. [13]. This can be done by assuming c,c, X ,X , Y * R and Y * L all proportional to the identity matrix, except either X (right-handed compositeness) orX (left-handed compositeness). In this case, when neutrino masses are neglected, the only spurion that breaks the flavour symmetry becomes proportional to the charged lepton Yukawa couplings, exactly as in MFV. A choice of basis where such a spurion is diagonal is always possible and LFV is only present when neutrino masses are turned on. In the latter case one can reconcile LFV with a scale of new physics close to the TeV scale, however, we think it is interesting to explore other options allowing for a TeV scale .
We start by taking Y * L = 0. To guarantee the stability of this condition we should also assume either a perturbative regime, where |Y * R | ≤ 1, or the hierarchy M m 0 , as explained in Sect. 2.5. In this case the dipole operator is dominated by the contribution (N Y , N X ) = (1, 2) (see Table  1) and the bound on μ → eγ leads to m 0 c > 33 TeV. As a next step we consider the case of Y * L = 0 and universal heavy fermion masses, namely c andc proportional to the unit matrix. With this assumption the Wilson coefficients (N Y , N X ) = (1, 2) of the dipole operator are diagonal in the mass basis, to LO. This, however, is not the case, in general, for the coefficients (N Y , N X ) = (3, 2), which are comparable in size to those of Eq. (6). Therefore, even working in the limit of vanishing "wrong" Yukawa coupling Y * L , and universal heavy fermion masses, we expect that with anarchic Y * From the discussion in Sect. 2.5, we know that a nonvanishing Y * L can be generated. Since such a Y * L will only depend on c,c and Y * R , it will be universal too. Admitting a universal Y * L would not change our conclusions and, for simplicity, we keep the assumption Y * L = 0. Moreover, from now on we set to unity the matrices c andc, absorbing their effect in the overall scale m 0 of the heavy fermion masses. The remaining coupling Y * R is described by a single parameter y: At variance with MFV, in IFV the sources of flavour breaking are both X andX . In the basis defining IFV the mixing matrices X andX are not diagonal, in general. They are generic complex matrices. At the LO, the lepton mass matrix is proportional to the product XX † . When we move to the charged lepton mass basis, in general, only the product XX † becomes diagonal, not X andX individually and they both can lead to LFV.
By exploiting the freedom related to our flavour symmetry, we can provide an alternative, but equivalent, description of the IFV scenario. We can choose a basis where X andX are diagonal, real and non-negative. In this case by means of the symmetry transformations of Eqs. (3) and (4) we can still maintain c andc proportional to the identity, but the matrix Y * R becomes a generic unitary matrix. In such a basis the lepton mass matrix XY * RX † is non-diagonal and LFV is now ascribed to the interplay between Y * R and X ,X . Within IFV the special case X =X , which was previously assumed to minimise FCNC effects, forbids any LFV effect. Indeed, if X =X , there exist a basis where both X and the mass matrix XY * RX † are diagonal at the same time. As long as neutrino masses are neglected there is no source of LFV and all the Wilson coefficients are diagonal in flavour space. The only difference with respect to MFV is that the spurion X is not proportional to the charged lepton Yukawa couplings, but to their square root. Therefore in discussing IFV we consider the general case where X andX are not equal.
In IFV no LFV is generated from the Wilson coefficients of the dipole and scalar operators with N X = 2, to LO. These coefficients are of the type X AX † , where A is any combination of c,c and Y * R , and they are automatically aligned to the charged lepton mass matrix, to LO. To estimate the other Wilson coefficients we should specify the choice of X andX . Rather than scanning for the most general possibility, here we provide an example. We assume that, in the basis where X andX are diagonal, their elements are given, up to coefficients of order one, by where λ ≈ 0.22 is a small parameter of the order of the Cabibbo angle. The choice of elements of nearly the same order of magnitude for X is motivated by a possible role that the matrix X can play in describing large lepton mixing angles, once neutrino masses are turned on [5]. The value λ 3 ÷ λ 4 is chosen here for convenience, to adequately suppress LFV. The choice for the elements ofX is fixed, at the level of orders of magnitude, by the relation (39). The transformations needed to diagonalise the lepton mass matrix are given in Eq. (38). Now we can complete our discussion concerning the Wilson coefficients that are bilinear in the spurions X ,X . Coefficients (N X , N Y ) = (2, 2) of the vector operators are of the type X AX † orXÃX † where A depends on c,c and Y * R . From Tables 2 and 3, we see that Y * R enters A only in the combina- Therefore, up to an overall coefficient of order one, in the basis where X andX are diagonal and real we have When we move to the charged lepton mass basis, we get The most stringent limits of Table 7 arise from μ − Au → e − Au and require which translate, respectively, into > √ y 2.2 TeV.
Here and in the rest of this section stands for either M or m 0 . The limit on the decay τ → 3μ gives rise to a similar bound: resulting in > √ y 1.9 TeV. Coefficients of the type 0) can arise for the contact operators of the type llqq and we expect, in the mass basis: The limits of Table 8 from μ − Au → e − Au require and we get > 5.0/ √ y TeV, respectively.
We are left with the coefficients that are quadrilinear in X ,X . In the IFV scenario the off-diagonal elements of the dipole operator (Q eγ ) i j are dominated by the coefficients which, in the lepton mass basis are approximately given by The limit (C eγ ) eμ,μe 2 < 2.5 × 10 −10 TeV −2  Tables 5, 6, 7 and 8 do not lead to more restrictive bounds on . For instance, the bound on the contact operators, produces the bound > (0.3/y) TeV. The bounds from the vector operators scale as √ y, while those from the dipole or contact operators scale as 1/y or 1/ √ y. Therefore a new physics scale around few TeV is still acceptable provided the Yukawa coupling y is close to one. A new physics scale ≈ 1 TeV is too large to allow, in the present framework, for an explanation of the central value of the a μ anomaly. Indeed, in the lepton mass basis we have and to match the required value, Eq. (32), we need = 0.56 TeV. Since the contribution to a μ scales with the inverse square of , choosing = 1 TeV gives a μ = 9 × 10 −10 , less than one third of the central value of the current anomaly. Concerning the electron EDM, we assume as in MFV that the sources of CP violation and LFV are the same. Since the Wilson coefficients of the dipole operators with N X = 2 are aligned in flavour space with the mass operator, we identify in (51) the dominant coefficients that can contain non-trivial phases. We estimate which, for > (0.16/ √ y) TeV, respects the bound of Eq. (30). In summary, even in the limit of "wrong" Yukawa coupling negligibly small, an anarchic Y * R requires a scale of new physics well above 10 TeV. One way to lower this bound consists in mimicking the case of MFV, where there is a single non-universal spurion, either X orX . In this case, when neutrino masses are neglected, LFV is absent. By discussing a large class of possible LFV effects, we have shown that MFV is not the only possibility to reconcile LFV bounds with a new physics scale around few TeV. In the example analyzed here, both X andX are non-universal and represent potential sources of LFV, which is allowed also in the limit of vanishing neutrino masses. The corresponding bounds on approach the TeV scale, provided the Yukawa couplings Y * R are close to one.

Two-site model
In the previous sections, we focussed on PC scenarios for charged leptons from a general perspective, exploiting a spurionic analysis. Now, instead, we consider a specific simplified composite Higgs model, the so-called two-site model [37]. Its relevant features are: In Table 9 is summarised our notation for generators, boson fields and coupling constants associated with each simple subgroup. The group for G comp has been chosen in order to provide a custodial symmetry. PC is assumed to arise from an unknown dynamical spontaneous symmetry breaking mechanism, which takes place at higher energies and which breaks G el × G comp into the diagonal group (which can be recognised as the electroweak gauge group). This spontaneous symmetry breaking mechanism, which in turns triggers the PC scenario, will be effectively described by linear couplings between elementary and composite bosons. charged under [SU(2) L × SU(2) R ] comp , which will be identified with the composite Higgs field, interacting only with the composite fermions. Table 10 summarises the quantum numbers for fermions and Higgs doublet.
We are ready now to introduce the Lagrangian of the twosite model: with where lower (upper) case letters denote elementary (composite) fields and the 'tilde' denotes SU(2) L singlets. Moreover, among the heavy vector bosons, we distinguish between those mixing with the SM gauge bosons, ρ * μ = {W * μ , B * μ }, and those that do not,ρ μ = {W μ ,B μ }.
In L el we can recognise a SM-like Lagrangian (with F a μν collectively denoting the field strength tensors for the elementary gauge bosons) without the Higgs doublet, which is now a composite particle. L comp contains mass terms for composite leptons, kinetic terms for composite leptons and bosons (we have collectively denoted the field strength tensors with ρ a μν ), their interactions, the Higgs sector and Yukawa interactions. We have allowed the Yukawa couplings to explicitly break the SU(2) R symmetry, which applies to the bosonic part of the composite sector and protects the model against large contributions to the ρ parameter. Such more general Yukawa interactions are technically acceptable, since the cutoff-dependent, loop-induced contributions to the ρ parameter from heavy fermion exchange are negligible, at least in the perturbative regime, as explicitly shown in [37]. Finally, L mix contains linear mass-mixing terms among elementary and composite particles. We have also included in L mix a mass terms for composite vector bosons. L mix explicitly breaks G el × G comp down to the diagonal subgroup, and it effectively reproduces the mechanism of PC.
From the above Lagrangian we can easily recognise the flavour symmetry group of Eq. (1) and convince ourselves that promoting Y * L , Y * R , m,m, and˜ to spurions with the transformation properties of Eq. (3) we actually restore G f .
The Lagrangians of Eqs. (60)- (62) are expressed in the elementary/composite basis. For the bosons, in order to switch to the mass basis (before EWSB), we diagonalise the mass mixings in L mix by the following field transformations: At this stage the fields A μ are massless, while ρ * μ andρ μ have masses of order M * . Taking into account EWSB, the diagonalisation of the mass terms becomes rather involved. To a good approximation, Eq. (63) still diagonalises the boson fields, while the mass matrices for the fermions read With this notation, one can then perform a suitable rotation that transforms to the mass basis: where V L ,R , U L ,R are defined in order to have Using the rotations (63) and (65) one can recast the Lagrangian (59) in terms of lepton mass eigenstates. Using an approximate expression for the rotation matrices V L ,R , U L ,R , assuming universal masses for the heavy leptons, (m i = m,m i =m), and retaining the leading terms relevant for FC neutral-current (FCNC), one gets From Lagrangian (66) one can easily read some of the most relevant features of FCNC in this class of models: • The first line of Eq. (66) accounts for FCNC interactions among light and heavy leptons and the Higgs. They contribute to low-energy processes such as i → j γ at loop level.
As a result, after EW symmetry breaking, we generate the operatorsē Li / Ze L j andē Ri / Zẽ R j , which are suppressed by v 2 /m 2 and v 2 /m 2 factors, respectively. The leading effects induced by these terms are the tree-level FCNC decay modes Z → i j as well as i → 3 j and μN → eN . Concerning the flavour structure of the various interaction terms, we remember that each Higgs electroweak doublet ϕ or its vacuum expectation value v occur accompanied by a power of Y * R or Y * L , while the lepton fields e L andẽ R come with X andX , respectively.
Before discussing LFV in this model we comment on the issues of renormalizability, gauge-invariance and UV sensitivity. The two-site model is a non-renormalizable effective theory and therefore it is valid only up to energies of the order of an UV cut-off UV . Such a cut-off has been estimated in Ref. [37], where it has been found that UV ≈ 8π M * /g * , by analysing self-interactions in the composite sector. The non-renormalizability is a consequence of the explicit breaking of the gauge symmetry G el × G comp , by the heavy vector boson masses, by the Yukawa couplings Y * L and Y * R in the composite sector and by the mixing between elementary and composite fermions and˜ .
This explicit breaking raises the question of the reliability and consistency of our results. We can always promote our effective Lagrangian to a gauge-invariant theory by interpreting the sources of gauge symmetry breaking as spurions. Providing the spurions with suitable transformation properties under the gauge group we can recover the gauge invariance under the full local group G el × G comp . Once we treat the spurions as dynamical fields, this procedure defines a (non-unique) embedding of the effective theory in a possible UV completion. The spurions include two types of degrees of freedom: the would-be Goldstone bosons, eaten by the heavy gauge vector bosons of the composite sector through the Higgs mechanism, and physical scalars i . By working in the unitary gauge, which we adopt in our computations, a generic amplitude comprises two separate contributions: one coming only from the exchange of the physical polarisations of the gauge vector bosons and one including the exchange of some physical scalar degrees of freedom i . Clearly our computation retains only the first one, while the second is missing and makes our results sensitive to the details of the UV completion. If the masses of the extra scalars i are close the cut-off scale UV , we expect that the contributions we are neglecting in our estimates of the Wilson coefficients are generically of order (M/ UV ) 2 and/or (m,m/ UV ) 2 .
The embedding of the model in a UV completion also shows that it is not consistent to deal with the HF case, within the effective theory. Indeed to keep the masses of the elementary leptons non-vanishing in the HF limit, we have to consider at the same time large ( ,˜ ), such that the ratios (X,X ) remain constant. In a UV completion ,˜ are proportional to VEVs that break the gauge symmetry G el × G comp and contribute to the masses of the heavy gauge vector bosons of the composite sector. Therefore, barring tuning of the parameters, we cannot make m,m M.

Lepton flavour violation in the two-site model
In this section we will present our results for the Wilson coefficients of the various LFV operators considered in Sect. 2, in the context of the two-site model introduced above. We work at the leading order in the loop expansion and we pay particular attention to the spurionic structure of the coefficients. Throughout this section we assume universal masses for the heavy leptons: m i = m andm i =m.

Dipole operators
Dipole operators are among the most important operators in the leptonic sector as they generate LFV radiative decays such as μ → eγ , leptonic EDMs and (g − 2)s. We will consider the one-loop contribution to dipole operators in the twosite model, assuming the validity of the perturbative expansion. The result of our computation of the full one-loop contribution to the electromagnetic dipole operator is shown in Appendix B. Given the fact that the composite Higgs scenario is naturally characterised by a strongly interacting regime, the use of perturbation theory can be questioned. Here we will restrict ourselves to the portion of the parameter space that is compatible with perturbation theory. This allows us to con-  Fig. 1. In particular, summing over the h, Z and W amplitudes, we find that the dominant effects for class (i) are given by in agreement with the result of Refs. [35,62]. As we can see, the amplitude of Eq. (67) has precisely the same spurionic structure as the one of Eq. (6). In particular, according to our spurionic classification of Sect. 2, it turns out that (C eγ ) h+Z +W is of order (N Y , N X ) = (3, 2), where N Y and N X are the orders of the expansion in (Y * R , Y * L ) and (X,X ), respectively. Interestingly, we observe that contributions of order (N Y , N X ) = (3, 2) containing only the Yukawa Y * R , which would be allowed by the flavour symmetries of our model, are absent, to one-loop order.
The leading effects for class (ii) start from the order (N Y , N X ) = (1, 2) and read where M stands for a common heavy boson mass, y = m 2 /M 2 , z =m 2 /M 2 and the loop functions are defined in the appendix. Notice that the first term of (C eγ )B +B * and (C eγ ) W * are aligned with the mass operator as long as m,m ∝ 1, at least. Therefore, as already pointed out in the literature [12], they cannot induce neither flavour nor CPviolating effects after switching to the mass basis for the SM leptons. Yet, we stress that they contribute to the leptonic g − 2. However, we point out here that the second term of (C eγ )B +B * , not discussed in the literature to our knowledge, is not aligned with the mass operator and therefore generates flavour and CP-violating effects. From Eq. (67) one can easily find that, in the anarchic scenario, the bound from BR(μ → eγ ) imposes that √ mm/ Y 10 TeV. In order to relax such a strong bound while keeping Y * R anarchic, we analyze here in great detail the solution with Y * L = 0, as already done in the modelindependent analysis of Sect. 2.
Setting Y * L = 0, the NLO effects stemming from class (i) are given by: where x = m 2 /m 2 . The first (second) terms of (C eγ ) h , (C eγ ) Z and (C eγ ) W arise from diagrams such as those shown in the left (right) plot of Fig. 2   last term of (C eγ ) h as well as the third terms of (C eγ ) Z stem from dimension eight operators such as and the corresponding diagrams are shown in the left plot of Fig. 3. Such contributions are of order (N Y , N X ) = (5, 2) and turn out to be suppressed by a factor of v 2 /m 2 compared to dimension six contributions. As a result, the parametric ratio between dimension six and dimension eight operators is of order X 2m2 /Y *  2) provided v 2 /m 2 1/16π 2 . The last term of (C eγ ) Z stems from dimension eight operators of the form 4) with loop exchange of heavy gauge bosons and heavy fermions Y * L = 0, we find that: where in the mass basis for the charged leptons we have with Therefore, the coefficient C eϕ of the dimension six oper- Notice that for Y * L = 0 the corrections to y SM are proportional to y SM itself, as in scenarios where the Higgs is a pseudo-Goldstone boson [63][64][65].

Vector operators
The Z boson interactions with charged leptons are also modified after we integrate out at tree level the heavy fermions. In particular, we find that where g L and g R are defined as follows: Switching to the coefficients C (1) ϕl , C ϕl and C ϕe we find that:

Contact operators
Contact operators induced by the exchange of heavy gauge bosons can be easily derived starting from the coefficients C (ρ) R R entering the interaction Lagrangian of heavy gauge bosons and SM leptons which we report in the appendix. Integrating out the heavy gauge bosons at tree level, we find that and where higher-order effects suppressed by additional factors of (g SM /g * ) 2 , v 2 /m 2 (m 2 ) and X 2 (X 2 ) have been neglected.

Phenomenological analysis
In this section, we evaluate the most relevant low-energy processes in the charged lepton sector in the context of the twosite model by making use of the Wilson coefficients derived in the previous section. One of the most interesting features of the two-site model is the absence of the Wilson coefficients (N Y , N X ) = (3, 2) from dimension six dipole operators, at least in the one-loop approximation. Such coefficients are particularly dangerous for LFV. As we saw in Sect. 3, their presence is not compatible with anarchic Yukawas Y * R if the compositeness scale is close to 1 TeV. We do not know if this feature of the model is an accident of the one-loop approximation or if it persists at higher loop orders. In the following we will neglect the (N Y , N X ) = (3, 2) contribution to dimension six dipole operators, assuming that higher loops provide a sufficient suppression. Our purpose is to check if under these conditions anarchic Yukawas are still viable or not when we consider a compositeness scale around 1 TeV.
Among the most interesting LFV channels are μ → eγ , μ → 3e, μ → e conversion in nuclei as well as τ LFV processes. However, hereafter, we focus on processes with an underlying μ → e transition since they are the best probes of composite Higgs models with anarchic Y * R . Concerning flavour-conserving processes, we are interested in the electron EDM and the muon g − 2. The current status and future experimental sensitivities for the above processes are collected in Table 4. As recalled in Sect. 3, there is a ∼3.5σ discrepancy between the SM prediction and the experimental value of the muon g − 2.
The dipole transition i → j γ is responsible for both LFV radiative decays (when i = j) like μ → eγ and flavourconserving processes like the electron EDM and the muon g −2 when i = j = e or μ, respectively. The branching ratio for the process μ → eγ can be written as Radiative LFV transitions are tightly related to the magnetic and electric leptonic dipole moments, which are also extremely sensitive probes of new physics. In particular, one can find that In concrete scenarios as our two-site model, a , d and BR( → γ ) are expected to be correlated. However, such correlations depend on the flavour and CP structure of the couplings which are unknown. In our discussion, we assume order one CP-violating phases and anarchic Y * R . Here we provide the dominant contribution to μ → eγ which arises from dimension eight operators. Focussing on the HB scenario (M m,m) with anarchic Y * R and X =X , it turns out that where we made use of Eq. (40) to eliminate X,X and where Y * R now stands for an average element of the anarchic matrix Y * R i j . Notice that the above expression is valid only in the mass insertion approximation which requires that m,m, |m− m| vY * R . As an example, if |m 2 −m 2 | ≈ m 2 , one can find implying that μ → eγ saturates its current experimental bound for m ≈ 1.5 TeV and Y * R ≈ 1. This is confirmed by our numerical results shown in Fig. 5 where we have reported the predictions for BR(μ → eγ ) as a function of the heavy fermion mass m. In Fig. 5, as well as in all other plots, we have assumed anarchic Y * R and Y * L and X =X so that the relevant flavour mixing angles entering μ → e transitions are X e / X μ =X e /X μ ∼ m e /m μ . The lower (upper) red line in the left plot refers to the case where . The most prominent features emerging by this plot are: and therefore the electron EDM poses a severe constraint on the heavy fermion scale at the same level of μ → eγ . Setting Y * L = 0, the EDM bound is satisfied with much lighter masses. In particular, we find that where in the last equality we have assumed that |m 2 −m 2 | ≈ m 2 .
In Fig. 6 we show our numerical results for d e . Indeed, Fig. 6 represents the analogous of Fig. 5 for the electron EDM case and similar conclusions drawn for μ → eγ apply here too.

μ → 3e
The process μ → 3e receives more contributions than μ → eγ as it is sensitive also to four-fermion operators. In particular, starting from the general expression of Refs. [53,66] and neglecting scalar and non-dipole photonic contributions which are negligible in our model, we find the following expression for its branching ratio: where and the explicit expression for the Wilson coefficients C ab can be found in Sect. 5. In the above equations, we have neglected scalar four-fermion operators, since they are very suppressed in our model. We remind the reader that, whenever the dipole operator is dominant in μ → 3e, there exists a model-independent correlation between the branching ratio of μ → 3e and μ → eγ given by As a result, the current MEG bound BR(μ → eγ ) ≤ 5.7×10 −13 already implies that BR(μ → 3e) 3.8×10 −15 . However, as shown in our numerical analysis and illustrated in Fig. 8, μ → 3e turns out to be dominated by non-dipole operators and the correlation of Eq. (95) is significantly violated. This is a relevant result, as within composite Higgs models with Y * L = 0, as well as in supersymmetric scenarios, Eq. (95) holds to an excellent approximation. Therefore, the two-site model with Y * L = 0 can be disentangled among other models if both μ → 3e and μ → eγ will be observed.
In particular, in the HB scenario (M m,m) with anarchic Y * R , m =m and X =X , we find implying that μ → 3e saturates its current experimental bound for m ≈ 1 TeV and Y * R ≈ 1. This result is well reproduced by our numerical analysis shown in the left plot of Fig. 7 where the behaviour of BR(μ → 3e) as a function of the heavy fermion mass m is shown for different values of Y * R . Notice that BR(μ → 3e) ∼ Y * R 2 /m 4 and therefore it has a much milder growth with Y * R and a slower decoupling with m compared to BR(μ → eγ )

μN → eN
As the last process governed by an underlying μ → e transition we consider μ − e conversion in nuclei. As in the μ → 3e case, this process receives effects from both dipole and four-fermion operators. In particular, the μ − e conversion branching ratio is defined as where ω capt is the muon capture rate while ω conv is given by the expression The flavour-conserving interactions of the heavy vector resonances with the light quarks have been derived from the respective couplings with the light leptons simply rescaling them according to the different charges of the quarks under Y and T 3 . Finally, the quantities D and V p,n refer to the overlap integrals between the wave functions and the nucleon densities [56] which we report for few relevant nuclei in Table 11.
An inspection of Eqs. (97), (98) and Table 11 shows that at present μ − Au → e − Au is the most sensitive probe of new physics among the various μ − e conversion in nuclei processes. In scenarios with dipole dominance, the following model-independent relation holds: and therefore BR(μ − Au → e − Au) 2.2 × 10 −15 after imposing the bound BR(μ → eγ ) ≤ 5.7 × 10 −13 . However, as in the case of μ → 3e, in turns out that in our model BR(μ − Au → e − Au) is dominated by non-dipole operators and the correlation of Eq. (105) is not at work. This is shown in Fig. 8. In particular, in the HB scenario (M m,m) with anarchic Y * R , m =m and X =X , we find and therefore the current experimental bound is saturated for m ≈ 3 TeV and Y * R ≈ 1. These expectations are fully confirmed numerically as shown by the right plot of Fig. 7 where we report BR(μ − Au → e − Au) as a function of the heavy fermion mass m for different values of Y * R . Since both BR(μ − Au → e − Au) and BR(μ → 3e) are proportional to Y * R 2 /m 4 , the consideration done for the case of μ → 3e hold here too. However, we point out that in our model BR(μ − Au → e − Au) is much enhanced with respect to BR(μ → 3e) as opposite to scenarios with dipoledominance where BR(μ − Au → e − Au) ≈ 0.6 × BR(μ → 3e). As a result, the simultaneous observation of these processes would enable us to disentangle among the underlying theory at work.

Conclusion
We have reanalysed the bounds on the compositeness scale derived from LFV and CP violation in the lepton sector, in the framework of PC. In the generic case of anarchic Yukawa couplings, such bounds are known to be quite severe, stronger than the analogous bounds from the quark sector and requiring a compositeness scale above 30 TeV. In this work we focussed on the case of vanishing neutrino masses and vanishing "wrong" Yukawa couplings, in the hope of minimising FCNC and CP-violating effects. We performed a general effective operator analysis, where the Wilson coefficients of the relevant dimension-six operators are determined by a flavour symmetry group and by a set of spurions. We have considered all lowest-dimensional operators leading to LFV violation and CP violation that can be constructed with the Higgs doublet and the SM leptons. We have formally expanded the Wilson coefficients in powers of the Yukawa couplings of the composite sector and in powers of the elementary-composite mixing terms. Our expansion includes, for the first time, terms quadrilinear in the elementary-composite mixing (X,X ). By exploiting the known limits on the Wilson coefficients we have shown that, even in the case of vanishing "wrong" Yukawa couplings, the anarchic scenario is not compatible with a compositeness scale of 1 TeV, barring accidental cancellations not incorporated in our general spurion analysis. We have also shown that there are schemes where a TeV scale can be accommodated, without necessarily reproducing the case of MFV. In the example we have provided this is achieved in a semiperturbative regime, where the universal Yukawa coupling y of the composite sector is close to unity. In this respect it is interesting to note that while the bound on NP coming from operators bilinear in (X,X ) scales with √ y, the limit from the operators quadrilinear in (X,X ) scales with 1/y, thus making essentially impossible to lower NP below the TeV. This also shows the importance of keeping terms quadrilinear in the mixing in the expansion of the Wilson coefficients. We have also derived the low-energy effective Lagrangian relevant to LFV and CP violation in the lepton sector in the two-site model, that includes heavy vector-like fermions as well as heavy spin-one particles. We focussed on the case Y * L = 0, finding a strong suppression of LFV effects compared to a generic composite Higgs scenario where Y * L = 0. In a perturbative analysis where we evaluated the relevant amplitudes at the LO, we found that the μ → eγ transition is dominated by dimension eight operators, at variance with our spurion analysis where the most important contribution comes from operators of dimension six. As a result the best probes of our scenario are μ → e conversion in nuclei and the electron EDM. Al least in that portion of the parameter space where perturbation theory is applicable, LFV allows a compositeness scale close to the TeV, even in the case of anarchic Yukawas. In this regime there are interesting relations among the various LFV transitions, which allow one to disentangle the model from other possibilities.

B One-loop contributions to the dipole operator
In the following, all lepton fields are assumed to be in the mass basis. We introduce a compact notation to collectively denote such fields, using an upper case Latin index to distinguish among them: The flavour index, when necessary, will be made explicit with a lower case Latin letter. For most of the expressions, however, this index will be omitted to simplify the notation. In those cases one should interpret the expressions as 3x3 matrices in flavour space.
In Fig. 9 we see the different kinds of loop diagrams contributing to C i j eγ . From the Lagrangian (59), after performing the rotations (63) and (64) to the mass basis, one can read the Feynman rules of interest for these computations.
We summarise them in Fig. 10 where g * = g cot θ 1 and g * = g cot θ 2 . Remember that this is a compact notation in which the flavour indices are omitted. Then each element of the above matrices is really a 3x3 matrix in flavour space (such as Y * R , Y * L ). We can now summarise the leading results for the different loop diagrams. In these expressions, terms of order (3) Fig. 9 One-loop diagrams contributing to the dipole operator in our model. In diagrams (1) and (2) we have A = 1, 2, 3, in diagram (3) we have B = 1, 2 L ,R can be found in the appendix   Notice that the result is expressed in a block matrix notation. For example, taking Eq. (123) and recalling that ϒ is a 9 × 9 matrix, ϒ 12 should be read as the 3 × 3 block of ϒ in position (1, 2). In this notation, the expression (ϒ 12 1 m ϒ 21 ) i j stands for the matrix element: