Dipole operator constraints on composite Higgs models

Flavour- and CP-violating electromagnetic or chromomagnetic dipole operators in the quark sector are generated in a large class of new physics models and are strongly constrained by measurements of the neutron electric dipole moment and observables sensitive to flavour-changing neutral currents, such as the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B{\,\rightarrow \,} X_s\gamma $$\end{document}B→Xsγ branching ratio and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon '/\epsilon $$\end{document}ϵ′/ϵ. After a model-independent discussion of the relevant constraints, we analyze these effects in models with partial compositeness, where the quarks get their masses by mixing with vector-like composite fermions. These scenarios can be seen as the low-energy limit of composite Higgs or warped extra dimensional models. We study different choices for the electroweak representations of the composite fermions motivated by electroweak precision tests as well as different flavour structures, including flavour anarchy and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U(3)^3$$\end{document}U(3)3 or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U(2)^3$$\end{document}U(2)3 flavour symmetries in the strong sector. In models with “wrong-chirality” Yukawa couplings, we find a strong bound from the neutron electric dipole moment, irrespective of the flavour structure. In the case of flavour anarchy, we also find strong bounds from flavour-violating dipoles, while these constraints are mild in the flavour-symmetric models.


Introduction
The discovery of the Higgs boson [1,2] has made the question what stabilises the electroweak scale more acute and has reduced the implementation of weak-scale naturalness to two possibilities: supersymmetry and Higgs compositeness. In this paper, we study the latter possibility, which arguably has received less attention in the literature. This is in part due to the difficulty in treating the strong interactions that are responsible for the Higgs bound state. Much progress has been made in recent years by warped compactifications of a e-mail: koenigma@thep.physik.uni-mainz.de b e-mail: neubertm@uni-mainz.de c e-mail: david.straub@tum.de higher-dimensional space-times [3][4][5], providing a weakly coupled dual description of the strong interactions in four dimensions. On the other hand, purely four-dimensional (4D) models have been constructed as well [6][7][8][9][10], with a particularly well-motivated example being models in which the Higgs arises as a pseudo Nambu-Goldstone boson, explaining its lightness with respect to the other, as yet unobserved, resonances. 1 In all these models, indirect constraints from low-energy precision observables play a crucial role. Generating fermion masses without excessive flavour violation singles out the mechanism of partial compositeness where, from a 4D effective theory point of view, the elementary Standard Model (SM) fermions obtain masses by mixing linearly with composite vector-like fermion resonances [12]. Since the degree of compositeness is required to be smaller for light quarks, tree-level flavour-changing neutral currents (FCNCs) mediated by composite resonances are automatically suppressed [13][14][15]. However, unless one is willing to accept a fine tuning of a few per cent, for TeV-scale resonances this suppression mechanism is not quite strong enough to suppress CP violation in kaon mixing, if no additional flavour symmetry is assumed [16][17][18][19]. In addition to flavour constraints, the mixing of the SM fermions with composite states with different electroweak quantum numbers leads to potentially large corrections to electroweak precision observables. In particular, custodial symmetry should be imposed on the strong sector to avoid a tree-level correction to the T parameter, and the representations of the composite fermions have to be chosen to avoid large tree-level corrections to the Zb L b L vertex [20,21].
In addition to electroweak precision observables and treelevel flavour-changing processes, an important class of constraints on these models is given by loop-induced dipole operators that contribute to radiative FCNC decays or fermionic dipole moments. The presence of heavy vector-like fermions charged under the electroweak gauge group implies a potential enhancement of these chirality-violating operators compared to the SM. A number of studies of dipole operators have been presented in the literature, either for Randall-Sundrum models in the KK basis [16,[22][23][24][25], for purely 4D models [26], or genuine five-dimensional (5D) calculations [27][28][29]. The aim of this work is to exploit the computational simplicity of the 4D models to study the impact of different choices for the fermion representations and of different flavour symmetries on the constraints from observables sensitive to dipole operators. To this end, we will use a generalisation of the framework of Ref. [6], considering one set of composite partners for each SM field, plus the additional states required for custodial protection of T and Z → b LbL . Our work can be seen as a complement to similar studies of electroweak and tree-level flavour constraints [30] and Zmediated rare decays [31] in the same models.
The remainder of this paper is organised as follows. In Sect. 2 we define our model setup. Section 3 is devoted to a model-independent discussion of dipole operators, the observables probing them and the constraints obtained from existing measurements. In Sect. 4 we derive approximate analytical expressions for the leading contributions to the Wilson coefficients of the dipole operators within our setup and estimate the size of subleading contributions. These results are then used in Sect. 5 to obtain numerical bounds on the masses of composite resonances derived from the various observables sensitive to dipole operators. While depending weakly on the choice of fermion representations, these bounds will depend strongly on the presence or not of a flavour symmetry in the strong sector. We will discuss the explicit examples of a U (3) 3 or a U (2) 3 flavour symmetry, only broken by the left-or right-handed composite-elementary mixings [32][33][34][35]. Since some results in the literature partially overlap with our results, it is mandatory that we compare our findings to them; we do this in Sect. 6. Section 7 contains our conclusions.

Setup
We consider simple 4D models, in which partial compositeness is implemented along the lines of Ref. [ 1) X and a set of fermion resonances that fill complete representations of G c . SM fermions and gauge bosons obtain masses from linear mixing terms with the composite states. Since G c is larger than the SM gauge group, there is some freedom in the choice of the fermion representations. In addition to the simplest case, involving just one SU (2) L doublet and one SU (2) R doublet ("doublet model"), we consider two cases ("triplet model", "bidoublet model"), which are motivated by the custodial protection of the Zbb coupling.
Using a notation where lower-case letters refer to elementary fields, while upper-case letters denote composite states, the part of the Lagrangian involving fermions reads • In the doublet model, where transforms as a (1, 2) 1/6 . Here and in the following, i, j are flavour indices and a, b, c are SU (2) L or SU (2) R indices. Fermion fields without a chirality index are to be understood as ψ = ψ L + ψ R , so that the mass terms for these fields are of the where L is a bidoublet transforming as a (2, 2) 2/3 , and R and R are SU (2) L and SU (2) R triplets, transforming as (3, 1) 2/3 and (1, 3) 2/3 , respectively. In component notation, the multiplets are given by In the Yukawa couplings, we have also used the triplets rewritten as 2 × 2 matrices, where L U transforms as a (2, 2) 2/3 and L D transforms as a (2, 2) −1/3 under the composite gauge group. U and D are singlets with the U (1) X charge 2/3 and −1/3, respectively. The components of the multiplets are After rotating to the mass basis, the light and mostly elementary SM fermions couple to the Higgs through their mixings λ with the composite states. For example, in the doublet model, the mass matrix of light quarks, after removing the mixing with the heavy fermions but before rotating to the mass basis, can be written as where v = 246 GeV is the Higgs vacuum expectation value, and similar expressions hold in the bidoublet and triplet models. At leading order in v/m Q,R , only the Yukawa couplings Y (and notỸ ) enter the mass matrix, which is why the latter are sometimes called "wrong-chirality" Yukawa couplings. Although they are not necessary for the generation of quark masses, they are present in many models, and we will see that they play a crucial role in the generation of dipole operators, so we keep them in our Lagrangians.

Effective Hamiltonian
We are interested in the electromagnetic and chromomagnetic dipole operators involving quarks, both flavour violating and flavour conserving. The relevant effective Hamiltonian can be written as where q = u, d and V = γ, g. We define the dipole operators as In the flavour-conserving case, one has C qqV = C * qqV , so in total there are 18 magnetic and 18 chromomagnetic quark dipole operators. Among those, the most phenomenologically relevant ones are the first-generation flavour-conserving operators contributing to the neutron EDM, the flavourviolating ones with down-type quarks contributing to FCNCs with B and K mesons, as well as Q ( ) cuV relevant for charm FCNCs. Before discussing the observables probing these operators in turn, we briefly summarise the QCD evolution that is necessary to relate the operators generated at a high new physics scale to the low-energy observables.

QCD corrections
The operators O ( ) q i q j γ and O ( ) q i q j g are subject to QCD renormalisation and undergo mixing. They evolve according to (omitting flavour indices) and equivalently for the primed coefficients. For the running from some high new physics matching scale μ h down to the top mass m t , one has at leading logarithmic order [36] For the evolution from m t down to some low scale μ l , the number of active quark flavours change and quark mass thresholds have to be taken into account. We list numerical values of the evolution coefficients η for the evolution from some exemplary high scale values to m t , as well as from m t to phenomenologically relevant low scales, in Table 1.
In writing Eq. (13), we have neglected the mixing of neutral or charged current-current (four-quark) operators into the dipole operators [37]. Although such operators are generated at tree level in our setup, we expect that their contributions to the dipole operators are small, since they are suppressed by additional powers of the composite-elementary mixing angles.
where μ l is a hadronic scale of order 1 GeV. The calculation of the contributions of the quark (C)EDMs to the neutron EDM is plagued by considerable hadronic uncertainties. An estimate obtained using QCD sum rules [38] yields Experimentally, the neutron EDM is already strongly constrained [39], |d n | < 2.9 × 10 −26 e cm at 90 % C.L.
Several experiments are in construction that plan to improve this bound by up to two orders of magnitude [40]. Indirectly, the neutron EDM is also sensitive to the CEDMs of second-and third-generation quarks. In the QCD evolution of the CEDMs to low energies, when integrating out a heavy quark, a finite threshold correction is generated to the three-gluon Weinberg operator, which directly contributes to the neutron EDM and mixes under renormalisation with the first-generation quark (C)EDMs [41]. Taking these effects into account, the bound (17) can be translated into bounds on the charm, bottom and top CEDMs, which read [42][43][44] |d c | < 1.0 × 10 −22 cm,

Down-type FCNCs
The most well-measured flavour-changing dipole transitions are the b → sγ /g processes probed in the inclusive decay B → X s γ . The corresponding decay probing the b → dγ /g transitions is even rarer in the SM due to the stronger CKM suppression and consequently is measured less precisely. Normalizing the current experimental measurements to the SM expectations for the branching ratios, one has at present [45][46][47][48] R bsγ = 1.13 ± 0.11, Beyond the SM, these quantities are modified as [37,49] where For the numerical bounds on the Wilson coefficients in the next sections, we imposed the constraints (20) at 2σ .
The s → dγ /g transitions are less constrained experimentally, since the long-distance dominance in K decay processes makes it difficult to relate experimental observables to the short-distance contributions. Nevertheless, a meaningful bound on the Wilson coefficients C ( ) sdg can be obtained from the measurement of the parameter / . With the conservative assumption that the new physics contribution to / should not exceed its experimental central value, one obtains the bound [50]

Charm FCNCs
Recent experimental hints that the direct CP asymmetry difference A CP between D → K K and D → ππ decays is larger than the SM expectation have attracted a lot of interest as a possible sign of new physics, also in the context of models with partial compositeness [25,51]. But even if the observed effect is not due to new physics, the upper bound on A CP can be used to put a constraint on the charm chromomagnetic dipole operator Q cug [52]. Following [44], we impose in the numerical analysis that the new physics contribution

Model-independent bounds
Given all the experimental constraints discussed above, we can derive model-independent bounds on the Wilson coefficients of the dipole operators. We list them in Table 2 at a renormalisation scale of 1 TeV, considering one purely real or purely imaginary Wilson coefficient at a time. The only operators in the effective Hamiltonian (10) we have not considered are the flavour-changing ones involving top quarks. Although they are not yet strongly constrained, they will be probed at LHC in the future through the decays t → qγ and t → qg, where q = u, c.

Analytical results for the Wilson coefficients
In this section, we derive approximate analytical expressions for the Wilson coefficients of the dipole operators for the three different choices of fermion representations. We denote by M ∼ m Q,R a generic composite mass, by λ a generic composite-elementary mixing parameter, by g an elementary gauge coupling and by g ρ the coupling of the composite vector resonances. Our goal is to obtain expressions for the Wilson coefficient to a given order in the small ratios v/M, λ/M, and g/g ρ . To this end, we first consider the case of a single generation of elementary and composite fermions. The relevant mass matrices arising in the three models are collected in Appendix 8. We diagonalise these matrices at a given order in the small ratios, rotate all couplings to the mass eigenstate basis and compute the Wilson coefficients. The resulting oneloop expressions for the contributions to the Wilson coefficients involving scalars or vectors are listed along with the relevant loop functions in Appendix 9. It turns out that the dominant contributions typically arise from diagrams with a heavy fermion-lifting the chirality suppression-together with a W , Z or Higgs in the loop. In Sect. 4.1, we first discuss these contributions in detail, before qualitatively discussing the additional contributions in Sect. 4.2.

Leading contributions
For a single generation of fermions, to leading order in the small parameters v/M, λ/M and in the limit of heavy vector resonances, we find that the Wilson coefficients can be written in the form 3 where q = u, d and V = γ, g. In this limit, the only relevant diagrams feature a Higgs, W or Z boson as well as a heavy fermion in the loop. We have computed all coefficients a qV in the doublet, triplet and bidoublet models and list them in Table 3. To illustrate our procedure, we give a detailed account of our calculation of a dγ in the bidoublet model in Appendix 10. We note that the value a uγ = 0 in the doublet model is not due to a symmetry but rather due to an accidental cancellation between the W , Z and Higgs contributions. An important result of our calculation is that, at leading order, there is no quadratic term in Y , as was also emphasised in Ref. [25] in the context of the c → ug dipole transition 4 . This means that in models in which the "wrong-chirality" Yukawa couplings are absent or suppressed, the dipole operators will be suppressed as well. We will discuss other contributions, which become the leading ones in the limitỸ → 0, in Sect. 4.2.
The result in Eq. (25) is only valid in the unrealistic case of a single generation of fermions. Taking into account all three generations and an arbitrary flavour structure, the full analytic diagonalisation of the mass matrices is clearly not feasible. Still, it is possible to obtain an approximate analytical expression valid for three generations of elementary and composite quarks by promoting Eq. (25) to a matrix equation in flavour space. Concretely, for i ≥ j, one has and analogously for up-type quarks, 5 where U Lq,Rq are the matrices diagonalizing the quark mass terms (9). We checked numerically that Eqs. (26)-(28) indeed give a very good approximation to the exact results obtained by numerically diagonalizing the mass matrices. 4 The proportionality to the wrong-chirality Yukawa was also found for the μ → eγ dipole in Ref. [53]. 5 In the triplet model, there is an additional factor of 1 √ 2 in front of u compared to Eqs. (26)- (28), cf. the mass matrices in Appendix 8.

Subleading contributions
Going beyond the leading order in the expansion of composite-elementary mixings and v/M and beyond the limit of heavy vector resonances, there are several classes of contributions that can become relevant in some cases, in particular in models without wrong-chirality Yukawa couplings. This can be the case e.g. in Randall-Sundrum models with the Higgs field strictly localised on the IR brane, where the presence of wrong-chirality Yukawa couplings is not required by 5D gauge invariance (see e.g. [54]). In general, no simple analytical expressions can be given for these subleading contributions, so our discussion will remain qualitative. We start by discussing the contributions of type (a). In the case of b → sγ /g and b → dγ /g, there are two contributions to the W -top loop that are only suppressed by the degree of compositeness of the right-handed top quark compared to the leading contribution. They read

Higher orders in elementary-composite mixing
where q = d, s, x t = m 2 t /m 2 W and the W coupling of the SM quarks is written as In general, one has where the coefficients a and b are given in Table 4.
The contribution in Eq. (29) due to the right-handed W tb coupling was first discussed in Ref. [26]. We emphasise that the contribution in Eq. (30) can be equally important in specific models, although it depends on the flavour structure of δg R W , while the contribution in Eq. (29) is present even for a flavour diagonal δg R W . Concerning the contributions of type (b), which involve heavy fermions in the loop, we merely note that they are parameterically of the same order as the ones of type (a), but can have a different flavour structure and are therefore more model dependent. Contributions involving the degree of compositeness of the left-handed bottom quark are suppressed by m b /m t , and we will neglect them.

Higher orders in inverse powers of the composite mass scale
We now consider terms that do not involve additional composite-elementary mixings, but are present even in the limitỸ → 0. These contributions are relevant for Wilson coefficients not involving the third generation in models where the "wrong-chirality" YukawasỸ are absent or suppressed. Such contributions arise for example from an expansion of the loop functions of the diagrams with a W , Z or Higgs and a heavy fermion to higher order in the ratios where m ψ is the heavy fermion mass. In fact, the analytic cancellation of the contributions proportional to Y 2 works only at the leading order of the expansion of the loop functions. As an example, we discuss the Higgs contribution to the down-type quark dipole operator for a single generation in the bidoublet model forỸ = 0. We find There are similar contributions suppressed by m 2 W /m 2 ψ and m 2 Z /m 2 ψ . In the case of the W contribution, there is the special feature that, in addition to the quadratic term in the mass ratio, there is also a logarithm that becomes dominant for large fermion resonance mass. In view of this complicated dependence, we refrain from giving full analytical expressions for this type of subleading contributions in all models, but simply keep in mind that, in the absence of wrong-chirality Yukawa couplings and sizable compositeelementary mixings, the Higgs, W and Z contributions to   Fig. 1, which is familiar from the two-loop chargino contribution to the EDM in split supersymmetry [55]. We estimate the contribution of this diagram to the Wilson coefficient in the limit Y = 0 as

Higher loop orders
up to an O(1) factor. We see that it can be safely neglected with respect to the leading contribution (25) even forỸ ∼ Y , but it can dominate compared to the contribution (33) if m ψ ∼ m Q ∼ m R is in the multi-TeV regime.

Diagrams with heavy vector resonances
Until now, we have only considered one-loop diagrams with a heavy fermion and a W , Z or Higgs in the loop, but there are also diagrams with a heavy vector resonance and a heavy fermion. These contributions are always parametrically suppressed by a factor g 2 ρ /m 2 ρ , where g ρ is the coupling and m ρ the mass of the vector resonance. In general, the analytical expressions for these contributions are complicated, since, in contrast to the W or Z contributions considered above, one has to keep the full dependence of the loop functions if the fermion and vector resonance masses are comparable. However, it is important to notice that in the limit where all the fermion resonances are degenerate, the contribution to the dipole operators from these diagrams is real and diagonal in the mass basis and thus does not contribute to any of the observables we consider, which always feature either flavour or CP violation. 6 In general, we find that in the case of non-zeroỸ , these contributions are always suppressed by v 2 g 2 ρ /m 2 ρ and/or m 2 ψ /m 2 ρ (which is preferred to be smaller than 1 since naturalness prefers light fermion resonances and electroweak precision tests require heavy vector resonances) with respect to the leading contribution (25), and we confirmed with a numerical scan that they are typically small. We will not consider this class of contributions in the following, but one should keep in mind that, in particular corners of the parameter space, they might be relevant in specific models and would tighten the bounds considered below.

Higher-dimensional operators
Finally, in a more complete theory like a composite Higgs model, there can be additional contributions that are not captured by our Lagrangians defined at the beginning of Sect. 2 and are therefore not calculable in our setup. This means that the bounds we obtain below can be viewed as conservative estimates. It is possible that there are additional contributions that make the bounds more severe; but, on general grounds, there is no reason to expect that these additional effects conspire with the calculable ones to eliminate the constraints.

Phenomenological analysis
We now proceed to a numerical analysis of the bounds on partial compositeness from observables sensitive to dipole operators. Since all these observables probe either flavour or CP violation (or both), the bounds crucially depend on the assumptions made on the flavour structure of the model. We start with the most popular assumption of flavour anarchy, which arises in models aiming at a geometrical explanation for the quark mass and mixing hierarchies but is known to have a problem (assuming TeV-scale resonance masses) with excessive CP violation in K mixing, unless one is willing to accept an O(10 −2 ) fine tuning of the relevant CP-violating phase. We then also consider models with a global flavour symmetry in the strong sector-either U (3) 3

or U (2) 3 -only broken by the composite-elementary mixings.
Our aim in this section is not to perform a full numerical analysis of these models and the contributions to dipole operators. Rather, we aim to provide analytical expressions for the dominant contributions to the dipole operators and use them to extract approximate lower bounds on the resonance masses from the experimental measurements. These results can then be used to judge how severely a model with a given choice of fermion representations and with a given flavour structure is constrained by the observables sensitive to dipole operators.

Models with flavour anarchy
We first consider the case of flavour anarchy [13][14][15][16][17][18][19], where all the couplings in the strong sector are assumed to have O(1) off-diagonal elements and phases. In general, all coefficients then depend on complicated functions of the anarchic Yukawa and mass matrix elements. To give simplified approximate expressions one can use the fact that, up to O(1) factors, the quark Yukawa couplings and the CKM matrix elements can be written in terms of the degrees of compositeness s L ,R (see Appendix 8 for their definition) as where Y can be understood as an "average" Yukawa coupling. In the following, we provide simplified expressions for the Wilson coefficients in terms of "average" parameters Y ,Ỹ and m i that keep track of how the quantities scale with the parameters, but we neglect O(1) factors coming from the flavour structure. We do, however, take into account the numerical factors derived in Sect. 4.
The leading contributions to the Wilson coefficients at the matching scale then read where V = γ, g. Arbitrary phases and O(1) factors are understood in all cases. Concerning the relative importance of the primed and unprimed flavour-changing Wilson coefficients, it is interesting to note that in b → s and b → d transitions the flavour prefactor is an order of magnitude larger for the primed coefficients, so observables in B decays sensitive to the primed Wilson coefficients, i.e. to right-handed flavour-changing neutral currents, are particularly promising in the anarchic model (see [56] for an overview of promising observables). For the s → d transition, the prefactors of primed and unprimed coefficients are comparable, and for the c → u transition the unprimed coefficient has a prefactor that is about a factor 30 larger than the unprimed one. For the b → s and b → d transitions, there is an additional important contribution that is only suppressed by the degree of compositeness of the right-handed top quark, as discussed in Sect. 4.2. Here, we give only a crude parametric estimate of this contribution, For all Wilson coefficients, there is in addition a subleading contribution not involvingỸ that is parametrically suppressed by m 2 h,W,Z /(m Q m R ) compared to the leading one (forỸ ∼ Y ), as discussed in Sect. 4.2.2.
Having fixed the parametric dependencies of the Wilson coefficients up to O(1) factors, we can proceed to put numerical bounds on the combination YỸ /(m Q m R ) and the corresponding quantities for the subleading contributions. These bounds are listed in Table 5. We make the following observations: • The strongest bounds come from the down quark (C)EDM and constrain the fermion resonance masses to be above 4-5 TeV for Y ∼Ỹ ∼ 1. • For Y ∼Ỹ ∼ 1, there is a multitude of bounds in the ballpark of 1-2 TeV. Since these refer to operators with different phases and flavour structures, we conclude that it will be hard to avoid all of them by fortuitous cancellations, even if the bounds listed here are subject to O(1) uncertainties. Consequently, if Y ∼Ỹ , dipole operators alone imply that sub-TeV fermion resonances are borderline and require a Y not much larger than 1. 7 • In models withỸ = 0, the bounds turn out to be quite mild and an anarchic flavour and CP structure is compatible with sub-TeV fermion resonances for Y 3, if only constraints from dipole operators are considered.

Flavour-symmetric models
Since the flavour anarchic model is not only plagued by strong constraints from dipole operators but also from meson-antimeson mixing induced at tree level, it has been suggested that the strong sector is invariant under a flavour symmetry that is only broken by the composite-elementary mixings of one chirality. The simplest case is a U (3) 3 symmetry broken by the composite-elementary mixings of right-handed quarks ("left-handed compositeness") or of left-handed quarks ("right-handed compositeness") [32][33][34]. Among the three models considered here, right-handed compositeness can only be realised in the bidoublet model, as it requires different mixings for left-handed up-and downtype quarks. While the U (3) 3 models successfully suppress FCNCs, they are strongly constrained by electroweak and dijet constraints, since they predict a significant degree of compositeness for one chirality of light quarks [30,34]. This problem is avoided in models with a U (2) 3 flavour symmetry in the strong sector, again broken only by one chirality of composite-elementary mixings [30,35].

EDM constraints in U (2) 3 and U (3) 3 models
In U (3) 3 flavour models with left-or right-handed compositeness, the parameters in the strong Lagrangian are generation invariant, e.g. for the triplet model, and analogously for the other models. It can be shown that in all models the only physical phases apart from the CKM phase reside in the wrong-chirality Yukawa couplingsỸ [34]. In flavour models based on a U (2) 3 symmetry, one has instead etc. As a result, there is an additional phase in the compositeelementary mixings related to the flavour symmetry-breaking spurions, but in the strong sector it is true as well that the only physical phases can be chosen to be the ones of theỸ couplings, which can be different for the third and the first two generations. Below, we will adopt a phase convention where Y is real. Consequently, in both U (3) 3 and U (2) 3 models, there is a clear-cut prediction for the flavour-conserving firstgeneration Wilson coefficients relevant for (C)EDMs, where in U (2) 3 , the masses and Yukawa couplings refer to those of the first two generations of composite fermions. Note that, in contrast to the anarchic model above, we have used a "=" sign, since there is no further O(1) factor in front. This leads to the bounds on the combination Y ImỸ m Q m R shown in Table 6. We conclude that sub-TeV fermion resonances in U (3) 3 models, or sub-TeV fermion resonances of the first two generations in U (2) 3 models, require Y ImỸ 0.05. As discussed above, in the limitỸ → 0 the strong sector carries no new phase both in U (3) 3 and U (2) 3 models. The remaining contributions to the EDMs involving the phases in the composite-elementary mixings are tiny, and hence there is no relevant bound.
In the case of the U (2) 3 model, ifỸ = 0 or the firstgeneration fermion partners are decoupled, the leading contribution to the up-and down quark (C)EDMs is absent. But also the third-generation wrong-chirality YukawaỸ 3 can contribute to the neutron EDM, if it is complex. On the one hand, it will generate a contribution to the top CEDM via the Wilson coefficient which leads to the bound shown in the last row of Table 6. On the other hand, a two-loop contribution to the first-generation EDMs proportional to Im(Ỹ 3 ) can arise, e.g. from the diagram in Fig. 1. Estimating this contribution naively as C qqV ∼ g 2 YỸ /(16π 2 m 2 ψ ), one would obtain a similar bound on m ψ of the order of 0.4 TeV for Y 3 ∼ Im(Ỹ 3 ) ∼ 1.

Flavour violation in U (3) 3 models
The leading contributions to the flavour-changing dipole operators in Eq. (25) vanish in models with U (3) 3 flavour symmetry and left-or right-handed compositeness. Subleading contributions to the unprimed Wilson coefficients arise, as discussed in Sect. 4.2. The strongest bound is on the coefficient C bsV , for which a crude estimate yields leading to the bound

Flavour violation in U (2) 3 models
In U (2) 3 flavour models with left-handed compositeness, the leading contributions to the Wilson coefficients read and all other coefficients are negligible. Again, there are no additional O(1) factors. Since the coefficients relevant for b → d and b → s transitions are correlated in these models, it is sufficient to quote the (stronger) bound derived from the B → X s γ branching ratio. It is shown in Table 7.
In U (2) 3 models with right-handed compositeness, the Wilson coefficients vanish at leading order in the compositeelementary mixings. Beyond the leading order, there are contributions both in left-and right-handed compositeness analogous to the ones in U (3) 3 models. They give rise to a bound similar to Eq. (50).

Comparison with the literature
Since some of the dipole operators have been considered in the literature in various models similar to the ones we studied here, we present below a detailed comparison of our findings with those of existing analyses. We find mostly agreement, but also some important differences.
• In Ref. [23], the B → X s γ branching ratio has been calculated in a model similar to our doublet model, the difference being that the right-handed quarks do not mix with an SU (2) R doublet, but with two singlets, such that the strong Yukawa couplings explicitly break custodial symmetry. For the leading-order contribution to the Wilson coefficient from W , Z or Higgs loops, this difference is, however, irrelevant. Up to an overall sign, we agree with the result for the charged Goldstone (W ) contribution [(a dγ ) W = 5/12 in our language], but disagree with the result for the neutral contribution [we find (a dγ ) h = −1/8 and (a dγ ) Z = −1/24]. • In Ref. [24], the B → X s γ branching ratio and the observable / have been estimated in a Randall-Sundrum framework. In the anarchic doublet model, which most closely resembles their setup, the bounds we obtain from these processes are consistent with the ones found in that reference. • Ref. [34] has given an estimate of the leading contribution to the B → X s γ branching ratio from loops with a Higgs boson or a charged or neutral Goldstone boson (W or Z ), corresponding to our Eq. (25), in the anarchic bidoublet model. We disagree with the fact that the leading-order contribution does not involve the wrong-chirality Yukawa couplings. In the same reference, EDMs in U (3) 3 models with left-handed compositeness were discussed, and it was claimed that the new CP-violating phase does not enter the EDM at leading order, since it can be shifted tõ Y . Our analysis shows that the converse is true: the leading contribution to the EDM is proportional toỸ . The bounds we obtain are shown in Table 6. • The authors of Ref. [28] have performed a 5D calculation of b → qγ processes in a Randall-Sundrum setup. The choice of fermion representations is similar to our triplet model, but the right-handed up-type quarks couple to a singlet. Furthermore, the model effectively has Y = Y . Our numerical estimates for the bound from the B → X s,d γ branching fractions are compatible with the numerical analysis presented in that work. • Ref. [26] contains a thorough analysis of the B → X s γ branching ratio and the observable / , closely following [23], in the triplet and bidoublet models (denoted TS5 and TS10, respectively) with flavour anarchy, settingỸ = Y . While we agree on the overall form of the results, we have several differences in the coefficients a qV . We present the details of our calculation in the bidoublet model in Appendix 10. • In Ref. [25], the c → ug dipole transition was considered in Randall-Sundrum models in the context of A CP in D → K K, ππ decays. In particular, the authors emphasise the dependence of the leading contribution on the wrong-chirality Yukawa coupling, and we confirm their findings in our 4D setup.

Conclusions
Dipole operators with quarks and an on-shell photon or gluon are generated at the one-loop level in theories based on the mechanism of partial compositeness, where the quarks get their masses by mixing with heavy vector-like "composite" fermions. Paradigm examples in this class of models are composite Higgs models or warped extra dimensions. The dipole operators contribute to numerous observables, like EDMs or FCNC decays, which can then be used to constrain these models. In this paper, we have performed an analysis of all dipole operators in the quark sector that are constrained by experiment within a simple 4D setup with a single set of vector-like fermions. We have chosen this simple framework so as to be able to study the effects of choosing different representations for the composite fermion fields and of imposing different flavour structures in the strong sector. Our main findings can be summarised as follows: •  (1), this implies that the average fermion resonance mass scale should be above 4 TeV or so. For larger Yukawa couplings, the bounds become even stronger. Apart from the neutron EDM bound, there is a multitude of bounds in the 1-2 TeV ballpark, as summarised in Table 5. • In models in which the wrong-chirality Yukawa couplings are absent, the bounds from dipole operators are mild, even for an anarchic flavour and CP structures. • In models featuring a U (3) 3 flavour symmetry broken only by left-or right-handed composite-elementary mixings, there is a bound from the neutron EDM that is as strong as in the anarchic case. It can be avoided by assuming the wrong-chirality Yukawa coupling to be real (or absent). • In models with a U (2) 3 flavour symmetry broken only by left-or right-handed composite-elementary mixings, the EDM bound can be avoided alternatively by raising the mass of the composite fermions of the first two generations. • In U (2) 3 flavour models with left-handed compositeness, there are bounds from flavour-violating dipoles, which are, however, very mild.
While our results have been obtained in the simple framework defined in Sect. 2, in more complete models the details of the analysis can be different. Here we only briefly comment on models where the Higgs is a pseudo Nambu-Goldstone boson (PNGB) [4,5], which are particularly well motivated in view of the lightness of the Higgs boson. Strictly speaking, these models are not a special case of the Lagrangian defined in Sect. 2. Given that the dominant contributions to the dipole operators come from diagrams with a heavy fermion and a Higgs, W or Z and does not depend on the details of the composite spin-1 sector, it is instructive to compare the fermion mass matrices in the two cases. For example, in the minimal composite Higgs model referred to as MCHM5 [57], the fermion mass matrix can be written in a form (see e.g. [58]) that, to leading order in the expansion in the Higgs vacuum expectation value, corresponds to the mass matrix in our bidoublet model (see Appendix 8) withỸ = Y . Thus we expect that with the identificationỸ = Y our results for the leading contributions to the Wilson coefficients also hold in composite PNGB models, up to O(1) factors.
There are several ways to extend our analysis. Also in the charged lepton sector, dipole operators arise and contribute to i → j γ decays, to the electron EDM or to the muon anomalous magnetic moment [53,59,60] (see also [61][62][63]). Finally, a global numerical analysis of all contributions to F = 1 and F = 2 processes, taking into account electroweak constraints, would be interesting. We leave this to a future publication.

Appendix A: Fermion mass matrices
In this section we list the mass matrices of the heavy fermion resonances in all three models. In the doublet model, they are given by In the triplet model, they are and In the bidoublet model, they are In the discussion of our results we have switched from the mixing parameters λ i to the sines of the mixing angles determining the degree of compositeness. The compositeelementary mixings are in general given by s i ≡ sin ϕ i ≈ tan ϕ i = λ i /m i ; more specifically, for the doublet and triplet models whereas in the bidoublet model

Appendix B: Model-independent formulae for the Wilson coefficients
Here we give the exact analytical expressions for the oneloop Wilson coefficients, which were used to obtain the approximate expressions given in the text. The Wilson coefficients of the q i → q j γ dipole operators, as defined in Eqs. (10)- (12) with i > j, can be written as where ψ denotes the fermion entering the loop and X can be either vector, scalar or a heavy gluon resonance. The parameter x is given by x = m 2 ψ /m 2 X . The expression for the primed Wilson coefficient C q i q j γ,g can be obtained from Eq. (59) by interchanging L ↔ R. The loop functions are defined as Note that the charge parameters Q ψ,V,S in the loop functions are not necessarily the electric charges of the corresponding particles but can be colour factors in the cases of either the external gauge field being a gluon or the loop involving heavy gluon resonances. A complete reference is given in (64)

Appendix C: Calculation of leading contribution in the bidoublet model
Here we illustrate the calculation for the leading-order correction to C qqγ in the bidoublet model for one generation.
This contribution is governed by diagrams with a heavy fermion and a W , Z or Higgs in the loop, as shown in Fig. 2. We obtain where m ψ and m X are the masses of the fermion and the boson in the loop. The V L ,R bψ X are the fermion-gauge couplings in the mass eigenbasis, as defined in Appendix 9. For the loop functions, we use the approximations where we only kept the first non-vanishing order for x → ∞. The gauge couplings up to quadratic order of the compositeelementary mixings are given in Table 9. Table 10 lists the mass eigenstates to order O(v). We have followed the convention of setting m Qu = m Qd → m Q and m U = m D → m R everywhere. Putting all the pieces together, we find the contributions listed in Table 11. Summing up these contributions, we have obtained the result from Sect. 4, i.e.