Pieces of the flavour puzzle

An overview of the flavour problem is presented, with emphasis on the theoretical efforts to find a satisfactory description of the fermion masses and the mixing angles.


Introduction
The origin of the parameters in the flavour sector of the Standard Model (SM), minimally extended to include massive neutrinos, is one of the most enigmatic questions in particle physics. Out of the 22 (20 if B-L is conserved) independent low-energy parameters Y i , which with some abuse of language can be called Yukawa couplings, 18 have been measured. Of the remaining four parameters, the absolute scale of neutrino masses is constrained in a limited range, the leptonic Dirac CP-violating phase starts to be constrained by global fits while the two possible Majorana phases are still unknown. A considerable effort has been devoted to the search for a more economic description, perhaps related to a new principle, such as the gauge principle. Gauge invariance and renormalisability allow one to describe strong and electroweak interactions of three copies of 15 different fermion species in terms of only three parameters. Nothing similar exists so far in the flavour sector and we usually refer to this as the flavour puzzle. Another aspect of flavour is related to the new particle threshold around the TeV scale predicted by all SM extensions addressing the hierarchy problem. Once new TeV particles transforming non-trivially in flavour space are introduced, it is very difficult to maintain the almost perfect agreement between predictions and observations that reigns in the SM. New sources of flavour-changing neutral currents (FCNC) and CP violations appear and the task is to keep them at an acceptable level. This is what we commonly mean by flavour problem, to distinguish it from the flavour puzzle. In this short review the focus will be on the first aspect, I will comment only briefly on the second one. Also, I am not aima e-mail: feruglio@pd.infn.it ing at reviewing all the existing models, but rather at revisiting some main ideas, guided by my own prejudices. Several aspects that have been left aside or only briefly mentioned in this paper are covered in a number of reviews [1][2][3][4]. In particular I will not do justice to the vast literature of fermion masses in grand unified theory [5][6][7], nor to that discussing the flavour puzzle in the framework of string theory [8,9].
There are different approaches to the flavour puzzle, with many intermediate possibilities. We may take a reductionist perspective: the Yukawa couplings Y i should be deduced from first principles. We postulate the existence of a fundamental theory from which Y i can be uniquely determined. Either by proceeding directly from the candidate theory or by appealing to some symmetry or dynamical principle, Y i are then computed in terms of a small set of input parameters. Probably the most striking fact about this program is that nothing approaching a standard theory of Y i exists, despite the decades of experimental progress and theoretical efforts. In another approach a major role is played by chance. There are many variants and practical implementations of this strategy. The Yukawa couplings Y i are typically mapped to a large number of order-one parameters that are considered as irreducible unknowns, like in models with Froggatt-Nielsen abelian flavour symmetries or with fermions living in extra dimensions. Also the simplest version of partial compositeness falls into this class. By scanning the order-one parameters we get probability distributions for masses and mixing angles. Alternatively we start from a fundamental theory, like string theory, which possesses a vast landscape of solutions, with no privileged ground state. The observed Yukawa couplings become environmental quantities and cannot be predicted, like the relative sizes of the solar planetary orbits (For review see [10]). We are allowed to ask much less ambitious questions. For instance, if we have knowledge of the statistical distribution of Y i in an hypothetical multiverse where the laws of physics follow our fundamental theory, we can ask how typical are the Yukawa couplings that we observe. Conversely, barring anthropic selections, we might assume that the observed Y i are typical and try to deduce information on the statistical distribution of Y i in the multiverse [11,12]. Such a variety of open possibilities shows how far we are from the solution of the puzzle and even from identifying the most relevant questions to be addressed.
Most of the parameters Y i are dimensionless and in a dynamical theory of flavour we have essentially no clue about the characteristic scale f . If active neutrinos are Majorana particles and B-L gets violated at a scale , then m 2 atm ≈ 0.05 eV strongly suggests a very large . However, in general and f are independent from each other. Thus there is no clear relation between f and other possible particle physics thresholds such as the TeV scale, relevant to the gauge hierarchy problem, or the grand unified scale. This makes it more difficult to identify unambiguous signatures to confirm or rule out a given model of fermion masses and mixing angles. For instance, the extrapolation of the Yukawa couplings from the scale f down to low-energies where they are measured can involve new particle threshold and/or unknown parameters, thus affecting our ability to test the high-energy theory.

Lessons from the quark sector
A first useful observation is that ratios of charged fermion masses and quark mixing angles can be represented by powers of the Cabibbo angle. Using λ = 0. 22 where all masses have been renormalised at the scale m Z . It is well known that also the elements of the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix can be expressed in terms of powers of λ: For comparison, in the lepton sector, where the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix is denoted by U , we have all |U f i | of order one, except for |U e3 | which is of order λ. The ratio between the solar and the atmospheric neutrino squared-mass differences m 2 21 /| m 2 31 | is of order λ 2 . Focussing on the quark sector, in a pioneering work [13] Froggatt and Nielsen observed that all the small dimensionless parameters of the quark sector such as the quark mass ratios and the CKM mixing angles can be interpreted as powers of the breaking parameter of a flavour symmetry. In this case the flavour symmetry group G f is abelian, G f = U (1) FN . A scalar field ϕ, carrying by convention a negative unit of the abelian charge FN, develops a vacuum expectation value (VEV) that can be parametrised as Quarks carry non-negative U(1) FN charges (the case with charges of both signs can be discussed as well) Under these assumptions the quark Yukawa couplings y u,d are given by where Y u,d are complex matrices with entries of order one, undetermined by the U(1) FN symmetry, while F X are real diagonal matrices, completely specified in terms of λ by the charges FN(X i ): The small quark mass ratios and quark mixing angles originate from the hierarchical structure of the matrices F X . Indeed, by taking FN(q 1 for the matrices V u,d defining the CKM mixing matrix Independently from the specific charge choice, this framework predicts the last equality being correct within a factor of 2. With λ ≈ 0.2, the correct order of magnitudes of the V CKM matrix elements can be reproduced by choosing, for instance, FN(q) = (3, 2, 0). The correct order of magnitudes of the quark mass ratios can be reproduced by choosing, for example r being a non-negative integer. If there is only one Higgs doublet, then we need r to be close to 2 to match the ratio m t /m b . If two Higgs doublets are present, other choices are possible by varying tan β = v u /v d . Several aspects of this class of models have been discussed in Refs. [14][15][16][17][18][19][20][21]. The construction relies on a spontaneously broken abelian flavour symmetry, but the final results (6)- (9) are valid in a more general context, where no symmetry is present to start with. A simple example is provided by a model with an extra spatial dimension, compactified on an orbifold S 1 /Z 2 to allow for 4D chiral fermions. The Lagrangian for a 5D spinor (x, y) reads where the mass m should be odd under the Z 2 parity sending y into −y. A possible choice is M being a real constant and (y) the periodic sign function. The 5D spinor has left (L) and right (R) chiralities in four dimensions with opposite Z 2 parities, such that only the even component developes a massless (zero) mode. Choosing, for instance, L even and R odd, the equation satisfied by the zero mode of L is The solution has an exponential dependence on y where the first factor provides the correct normalisation. The zero mode is localised near y = 0(π R) for M > 0(< 0). In the limit M = 0 the zero mode becomes flat in y. A formally identical solution holds for the zero mode of R , if we choose L odd and R even and we start from a 5D mass term with the opposite sign. If the Higgs field is strictly localised at one of the two branes, for instance the one at y = 0, the Yukawa interactions will be proportional to a Dirac delta δ(y) and we can reproduce the same pattern of Yukawa couplings of Eq. (6) with matrices F X now given by [22] (for early ideas to generate hierarchical Yukawas with extra dimensions see [23]): where μ i and ρ are specified in terms of the bulk quark masses and the geometry of the extra dimension; see Table 1. The suppression factors F X i represent the values that the profiles of the fermion zero modes take at the brane where the Higgs field is localised. On that brane, generic O(1) Yukawa couplings Y u,d with the bulk quark fields are postulated. The role of the Froggatt-Nielsen charges is here played by μ i and ρ, which determine the profiles along the extra dimension of the zero-mode wave functions: Table 1 Parameters μ i and ρ in models with an extra dimension compactified on an interval. In a flat (warped) metric the Higgs field is localised on the brane y = 0 (y = R ) and the ultraviolet cut-off is denoted by (1/R). The fermions are described by five-dimensional spinors, with bulk masses M i . In the warped case, when the framework is applied to the gauge hierarchy problem [24], R, R are length scales of the order of the inverse Planck mass and the inverse TeV scale, respectively There is no flavour symmetry: the hierarchical structure of quark masses and mixing angles is dictated by geometry in the compact space. Similarly, in the partial compositeness scenario [25], light fermions get hierarchical masses from the mixing between an elementary sector and a composite one. As a toy realisation of this idea, consider a model where the composite sector contains, for each SM fermion, a pair of heavy fermions allowing a Dirac mass term of the order of the compositeness scale and a mixing term with the SM fields [26,27] The first line represents the mixing between elementary and composite sector, the second line displays Dirac mass terms for the fermions of the composite sector and the last two lines show the Yukawa interactions that, by assumption, are restricted to the composite sector alone and described by strong couplings Y u,d ,Ỹ u,d ≥ 1. By integrating out the composite sector under the assumption M i v, we get lowenergy Yukawa interactions for the elementary sector whose leading-order (LO) terms have the structure given in Eq. (6) with matrices F X parametrizing the elementary-composite mixing: The same pattern arises when matter chiral multiplets X i of the MSSM are coupled to a superconformal sector in some finite energy range [28][29][30], from an ultraviolet (UV) scale down to a lower scale c . Generic O(1) Yukawa couplings Y i j at the scale 1 .8 × 10 4 3.2 × 10 5 6.9 × 10 −9 2.6 × 10 −11 (c L γ μ u L ) 2 1.2 × 10 3 2.9 × 10 3 5.6 × 10 −7 2.5 × 10 3 3.6 × 10 3 3.9 × 10 −7 1.9 × 10 −7 (b L γ μ s L ) 2 1.4 × 10 2 2.5 × 10 2 5.0 × 10 −5 1.7 × 10 −5 can undergo a sizable renormalisation induced by the corrections to the Kahler potential. In the superconformal window the chiral multiplet X i can have a large positive anomalous dimension γ i and the Kahler potential at the scale c becomes where Moving to a basis of canonical kinetic terms, the Yukawa couplings at the scale c are renormalised and we find again the same pattern of Eq. (6), without imposing any symmetry. In the previous examples the anarchical pattern of Y u,d may result in strong bounds on the scale of new physics NP associated to particles carrying flavour quantum numbers and representing new sources of FCNC and/or CP violation. In the absence of a concrete realisation, it is difficult to estimate reliably the corresponding effects, also because in general the scale of new physics NP and the scale of flavour physics f are independent from each other. A possibility is offered by a spurion analysis [31], analogous to that prescribed by minimal flavour violation (MFV) [32]. To this purpose we assume that the new degrees of freedom have non-trivial flavour properties and that the flavour-violating effects are completely specified by the same spurions that are responsible for fermion masses and mixing angles. Moreover, we assume that the dominant flavour-violating contributions admit an expansion in power series of the spurion fields. We start by noticing that the pattern of Eq.
The full symmetry G f is explicitly broken by both the matrices Y u,d and F X . However, it can be formally restored by treating Y u,d and F X as non-dynamical spurion fields possessing suitable transformation properties. To this aim the Yukawa couplings should transform only under the "hidden" group SU(3) H 3 : The suppression matrices F X are the interface between SU(3) H 3 and SU(3) 3 , and they are assigned appropriate transformations under both factors to guarantee the invariance of the Yukawa interactions described by Eq. (6) under SU(3) 3 ×SU(3) H 3 . The starting point of the spurion analysis is similar to that of MFV. Indeed the maximal flavour symmetry felt by quarks is SU(3) 3 , as in MFV. However, there are more spurions than in MFV, the irreducible ones including now F q , F u c , F d c , Y u and Y d . One of the most dangerous effects originates from the effective operator The contribution of this operator to the CP-violating K parameter is enhanced at the level of both the hadronic matrix element and the QCD corrections and sets one of the most stringent bounds on the scale of new physics NP ; see Table 2.
Assuming a generic O(1) phase for the overall coefficient we need not to spoil the SM prediction for K . This, together with other constraints, suggests that a fully anarchical pattern in Y u,d is probably not tenable if new flavoured physics is present at the TeV scale [35]. When such a spurion analysis is applicable, the estimate of Eq. (26) represents a sort of lower bound on the size of the expected effect and larger contributions are possible [36]. For example in supersymmetric extensions of the SM with a U(1) FN flavour symmetry and gravitymediated supersymmetry breaking the operator considered in Eq. (26) receives contributions from box diagrams with squarks/gluino exchange that are typically larger than the one quoted in Eq. (26). The reason is that in the U(1) FN case the true flavour symmetry is much weaker than SU(3) H 3 × SU(3) 3 and it allows sizable off-diagonal terms in both LL and RR blocks for the first two generations of the down squark mass matrix. For instance, with the charge assignment of Eq. (10), the mass insertions (δ d 12 ) L L and (δ d 12 ) R R are both proportional to λ and the operator (sd c )(s c d) has an overall parametric suppression 1/16π 2 × λ 2 / 2 NP , milder than the one in Eq. (26). This conclusion can be evaded in models with gauge-mediated supersymmetry breaking [37,38].

From quarks to leptons
In the lepton sector we have no evidence for strong hierarchies in mixing angles or in neutrino masses. Hierarchy shows up at the level of charged-lepton masses. In terms of the suppression factors F X i this means For example an acceptable set of charges is Here we focus on Majorana neutrinos. In the context of a type I see-saw mechanism right-handed neutrinos ν c have their own suppression matrices F ν c . Yukawa couplings y ν,e and the mass matrix M of ν c read where Y ν,e,c are complex matrices with unknown entries of order one and M 0 is a mass parameter. At low-energy the active neutrino mass matrix m ν is given by with no dependence on the suppression matrices F ν c . A drastic realisation of this picture is the framework of Anarchy [39][40][41][42][43], which corresponds to the case Table 3 Possible choices of FN charges for the 5 representation in a class of SU(5) grand unified models, from Ref. [46]. The second column shows the value of the FN symmetry breaking parameter optimizing the fit to fermion masses and mixing angles In the anarchic framework the mass matrix for light neutrinos is with undetermined order-one matrix elements. This implies mixing angles and neutrino mass ratios of O(1), in rough agreement with the data. No special values for these quantities is expected. Indeed, before we knew θ 13 from the experiments, Anarchy successfully anticipated values close to the upper bound at the time. Global fits of present data hint at deviations of the atmospheric mixing angle θ 23 from π/4. Today these indications are still weak, as shown by the instability of the best fit value against different fitting procedures. The persistence of these deviations in future tests would further strengthen the case for Anarchy. Anarchy represents an extreme possibility and milder realisation of the relations (32) are possible. For instance, in the context of SU(5) grand unified models, with a Froggatt-Nielsen U(1) FN abelian symmetry, neutrino masses and mixing angles can be reproduced, at the level of order of magnitudes, by several choices of the FN charges for the 5 multiplets hosting the lepton doublets, as shown in Table 3. FN charges for fermions in the 10 representations can be suitably chosen so that, by varying the unknown order-one parameters, reasonable distributions for charged lepton mass ratios, quark mass ratios and quark mixing angles are obtained [44][45][46][47][48]. A naive comparison of the distributions for neutrino masses and mixing angles with data do not appear to favour Anarchy over the other possible charge assignments, as can be seen from Fig. 1. I would personally find more appropriate to use the term Anarchy to denote the approaches giving rise to the results (6) and (30)- (31) where the absence of any special pattern resides in the matrices Y , rather than to indicate the special case defined in Eq. (32). If this framework also comprises new flavoured particles at the TeV scale, severe bounds from lepton flavour violation (LFV) apply, under assumptions analogous to those spelled for the quark sector. The irreducible sources of flavour violation in the lepton sector include the matrices Y e , F e c and  l and LFV can occur even in the limit of vanishing neutrino masses. Notice that, though MFV cannot be extended unambiguously to the lepton sector [49], it predicts no LFV if neutrinos are massless since in this limit the only relevant spurion in the lepton sector is Y e , which can always be chosen diagonal. The dipole operator contributing to LFV is The charged-lepton mass matrix is proportional to (F e c Y e F l ).
In general the combinations (F e c Y e F l ) and (F e c Y e Y † e Y e F l ) are not diagonal in the same basis, not even in the case of universal F l of Eq. (32), and radiative decays of muon and tau are expected. Agreement with the most constraining upper bound, BR(μ → eγ ) < 5.7 × 10 −13 , requires NP well above 10 TeV [50,51]. As in the quark sector, a completely anarchical matrix Y e and flavoured physics at the TeV scale are difficult to reconcile. A sufficient condition for the absence of LFV is that Y e , F e c and F l are diagonal in the same basis, as suggested in some models. Alternatively we can look for special forms of these matrices [52,53], possibly dictated by some symmetry requirements.

Realisations in grand unified theories
A welcome feature of the above description is that it can be adapted to grand unified theories (GUT) where quarks and leptons are hosted in the same multiplet of the gauge group. In SU(5) the gauge symmetry requires Our previous results, Eqs. (10) and (29), come very close to this requirement if we choose r = s and t = 1. If we accept a couple of tunings in the unknown O(1) parameters Y u,d , we can force the equality (35) and still have a decent description of both the quark and lepton mass spectrum. As we have seen F ν c drops from the low-energy quantities. It is instructive to consider also the ansatz F 5 ∝ 1. In this case the hierarchy among fermion masses is entirely due to F 10 . From Eqs. (6) and (30) we see that the mass ratios in the upquark sector are the square of the respective mass ratios in the down-quark and in the charged-lepton sectors, which is correct in first approximation. The large lepton mixing corresponds to a large mixing among d c quarks [54], unobservable in SM weak interactions, but with possible observable effects if transferred from quarks to squarks in SUSY exten-sions of the SM [55]. A minimal model with Higgs bosons in the 5 representation would lead to the unrealistic relation y e = y T d , but the contributions from other Higgs representations or from non-renormalisable operators can solve this problem [56,57] without altering the picture.
At first sight this description does not seem to be compatible with an SO(10) GUT. The most general renormalisable Yukawa interaction of three copies of fermion generations transforming as 16 of SO(10) reads The pattern of Yukawa couplings in Eq. (6) can also be thought to arise from a rescaling of the fermions fields, with the constraint that fermions belonging to a given irreducible representation of the gauge group have to undergo the same renormalisation. By assuming that the matrices Y 10 , Y 120 , Y 126 have complex elements of order one and that the fields 16 undergo a wave function renormalisation we see that all members of a 16 representation are affected in the same way. Even accounting for the Clebsch-Gordan coefficients arising from Eq. (36) and the different overall scales associated with the Higgs VEVs H u,d , we cannot reproduce the observed hierarchies of u, d and e masses. Such a discouraging starting point has been successfully modified in a construction by Kitano and Li [58], recently revisited in Ref. [59]. The model is a SUSY SO(10) GUT realised in a flat fivedimensional space time, the fifth dimension being compactified on an interval [0, π R] whose inverse size is of the order of the GUT scale. The N = 1 5D SUSY corresponds to an N = 2 4D SUSY, which is broken down to N = 1 as a result of appropriate boundary conditions. The model comprises a 5D vector supermultiplet decomposing as a 4D vector 45 V multiplet plus a 4D chiral 45 multiplet. In the bulk there are also three copies of 5D hypermultiplets, equivalent to 4D chiral multiplets 16 and 16 c , with bulk masses M i (i = 1, 2, 3). The boundary conditions allow zero modes only for 45 V and 16. A Yukawa superpotential analogous to Eq. (36) is localised at the brane y = 0. Prior to the SO(10) symmetry breaking, the wave functions of fermion zero modes evaluated at y = 0 effectively drive a rescaling of the Yukawa couplings, as described by (37). The suppression factors in F 16 are (see Eq. (16)): where μ i = M i / and ρ = π R. At this stage the Yukawa interactions are not yet able to reproduce the known hierarchies of u, d and e masses. The key ingredient of the model resides in the gauge interaction of the hypermultiplets. The 5D SUSY gauge interaction contains a 4D Yukawa interaction among 16 i , 16 c i and 45 , controlled by the gauge coupling g 5 , that can be combined with the bulk mass term: The chiral multiplet 45 has no zero mode but can acquire a non-vanishing VEV, 45 = v 3/2 , which breaks SO (10) down to SU(5)×U(1) X . The Yukawa interaction of Eq. (39) gives rise to effective SO(10)-breaking bulk masses: where Q r X is the U(1) X charge of the different SU(5) components inside the 16 multiplet: Q r X = (−1, +3, −5) for r = (10, 5, 1). We are back to the SU(5) case; see Eq. (35): but now the profiles F r i only depend on four free parameters: μ i and k. Neutrinos are described within a type I see-saw mechanism, as in Eq. (31), with masses for heavy Majorana neutrinos originating from the VEV of the SU(5) singlet in the 126 H representation. The model contains many parameters of order one. After rephasing of the relevant fields there are 27 real parameters coming from the matrices Y 10 , Y 120 , Y 126 and eight real parameters describing the embedding of the two light Higgs doublets within 10 H and 120 H . Despite the large number of parameters the agreement with data is not a priori guaranteed, since there are only four profile parameters to describe hierarchical mass ratios and mixing angles. Indeed a fit to an idealised set of 17 observables leads to good agreement only for large values of tan β, for both normal (NO) and inverted (IO) neutrino mass ordering.
A closer inspection reveals that fitting fermion masses and mixing angles in the IO case requires a fine-tuning of the Yukawa parameters. By generating a large sample of random order-one Yukawa parameters, the fit can be repeated by keeping at each iteration only 12 free parameters, four for the profiles and eight for the relevant Higgs combinations. The distributions of the minimum χ 2 over the number of degrees of freedom are shown in Fig. 2 for NO and IO. We see a clear difference between the two cases. While in the IO case we need about 10 5 samples to reach a p-value close to 0.05, in the NO case in about 1 % of the cases we have p > 0.05. The model needs a severe fine-tuning of the "anarchical" parameters in the IO case, while the NO one is realised much more naturally. The most probable values of the profile parameters give F 5 ≈ (0.07, 0.22, 0.63), showing that approximate Anarchy is an output rather than an input of the present construction. Fig. 2 The distributions of minimised χ 2 /ν for NO and IO in neutrino masses and for tan β = 50, from Ref. [59] Focussing on the NO case, there is no preferred value of the leptonic Dirac CP phase. The lightest neutrino mass is predicted below 5 meV, corresponding to a hierarchical neutrino mass spectrum while |m ββ | lies in the range 0.1-5 meV; see Fig. 3. Any positive signal in the current generation of experiments aiming at measuring neutrino masses or |m ββ | in the lab would essentially rule out the model. The hierarchy in the right-handed neutrino spectrum is very pronounced and the corresponding mass distributions are peaked around 10 6 GeV, 10 8 GeV and 10 14 GeV, as shown in Fig. 4.
In summary, fermion masses and mixing angles are well described by the map in Eqs. (6), (30) and (31), in terms of input parameters of order one, the elements of the Y matrices. Such a map can be realised in several different frameworks and does not necessarily need an underlying symmetry. The setup is compatible with both SU(5) and SO(10) grand unification and with the known solution to the gauge hierarchy problem. On the weak side, additional ingredients are probably needed to control the new sources of FCNC and CP-violations arising from new flavoured physics at the TeV scale. Moreover, all entries of the Y matrices are independent free parameters and it is not possible to make absolute Fig. 4 The predictions for the masses of RH neutrinos obtained for χ 2 min /ν < 2.21 in case of NO and tan β = 50, from Ref. [59] predictions, beyond the order-of-magnitude accuracy. This is clearly a major limitation, since we would like to test the theory at the level of the best available experimental precision. Finally the map in Eqs. (6), (30) and (31) might be an oversimplified version of a more accurate description and indeed there are several variants of the frameworks briefly mentioned in Sect. 2 that modify the results of Eqs. (6), (30) and (31) [60].

Flavour symmetries
Given the successful use of symmetries in other contexts of particle physics, it is natural to investigate whether Yukawa couplings can be constrained by a flavour symmetry. The largest possible classical flavour symmetry of a theory with the particle content of the SM is G MFV = U(3) 5 and corresponds to the limit in which the Yukawa couplings are turned off. The observed fermion masses and mixing angles break G MFV almost completely to a residual symmetry that includes the weak hypercharge and the combination B-L. Similarly, in any realistic model based on flavour symme- Fig. 3 The predictions for various observables obtained for χ 2 min /ν < 2.21 in case of NO and tan β = 50, from Ref. [59] tries, the flavour symmetry group G f is broken. In predictive models the breaking is spontaneous and occurs through the VEVs of a set of scalar fields ϕ transforming non-trivially under G f . The VEVs ϕ are either postulated or determined by minimizing a G f -invariant energy functional V (ϕ). The Yukawa couplings become dynamical variables evaluated at the minimum of V (ϕ): Y( ϕ / f ). A huge number of models can be constructed according to this set of rules, depending on the choice of G f (global, local, continuous, discrete, abelian, non-abelian), and on the choice of representations for scalars and fermions. An attempt to start from the full G MFV symmetry is described in Refs. [61][62][63][64][65]. The framework is that of MFV [32]. Under SU (3) 3 quarks transform as in Eq. (24). Yukawa couplings are promoted to spurions transforming as to ensure invariance of the Yukawa interactions under G MFV . By analyzing a general G MFV -invariant potential depending on y u,d , it has been proven the existence of stationary points where pointing to an approximate U(2) 3 symmetry of the quark spectrum and providing a good first-order approximation of quark masses and mixing angles. By extending G MFV to the lepton sector, other stationary points have been identified: Such a solution nicely exhibits maximal θ 23 , but needs m 1 = m to further enforce a large, undetermined θ 12 . Goldstone bosons arising from the breaking of G MFV can be eliminated by gauging the flavour symmetry [66]. Correction terms are needed to promote the leading-order picture into a more realistic theory. If there is new physics close to the TeV scale, the advantage of assuming the largest possible flavour symmetry G MFV is the high degree of protection against FCNC induced by the new flavoured degrees of freedom. If the only sources of flavour symmetry breaking are the SM Yukawa couplings, such a maximal symmetry is minimally violated [32]. The classification of the expected effects is unambiguous in the quark sector, and can be extended in several ways in the lepton sector [49,67].
Continuous flavour symmetry groups such as SO(3) and SU(3) have been proposed [68][69][70], also in the context of grand unified theories [71,72], with the three fermion families assigned to an irreducible triplet representation. Charged fermions of the first two generations are much lighter than those of the third generation and consequently within SU(3) or SO(3) we are forced to introduce large breaking terms. Alternatively we can start from the smaller flavour group U(2) and assign the first two generations to doublets and the third one to singlets. Within the simplest realisation, such an assignment in the quark sector leads to [73] where phases have been omitted. To correctly reproduce quark masses without appealing to cancellations among the matrix elements we need | u,d | | u,d | 1, which corresponds to a sequential breaking of U(2): The following relations can be derived [74,75] by a perturbative diagonalisation of y u,d : up to corrections of relative order , numerically close to few percent. The phase ϕ is arbitrary. These predictions can be translated into constraints on the (ρ, η) plane and compared to the region presently allowed by the data, under the assumption that the SM correctly describes all relevant processes. Using as inputs quark masses, the maximally allowed values for sin 2β and for |V ub /V cb | derived from (49) are too small and the above set of Yukawa matrices are now excluded at the 3σ level [76,77]. 1 Modifications of the ansatz (46), where either the 13 matrix element is non-vanishing or the element 32 is of the same order of the element 33, have been proposed to recover agreement with the data [76,80]. This is a nice example of a model of fermion masses based on symmetry requirements and leading to testable predictions. In this case the predictions were initially supported by data. Later on the precise measurements at the B-factories and the improved theoretical knowledge of both perturbative and non-perturbative QCD effects have considerably reduced the errors on the CKM matrix elements. Also the knowledge of light quark masses has improved and deviations from (49) are now significant. Unfortunately not all present models and ideas can be tested at the same level of accuracy. Independently of their dynamical origin, y u,d in Eq. (46) are an example of textures, a simple-minded but predictive approach to the problem of fermion masses and mixing angles pioneered by Fritzsch [81,82] and Weinberg [1]. In this approach we assume a privileged flavour basis and some special energy scale where some of the entries of the matrices y f vanish. In the quark sector, this can lead to relations between the elements of the CKM mixing matrix and the quark mass ratios which can be precisely tested. Not all zeros give rise to relations among physical quantities. By performing weak basis transformations, we can generate zeros in y f , which have no physical implications [83]. For instance, working with two generations, we can always set to zero the (1, 1) entry of both y d and y u by performing a unitary transformation on the SU(2)-singlet quark fields. This transformation corresponds to a change of basis that has no physical consequences. If, in such a basis, we impose the requirement that y d and y u are symmetric matrices, this represents a non-trivial requirement. By further assuming small off-diagonal matrix elements, we get in first approximation the well-known Gatto-Sartori-Tonin relation [74], Eq. (48). Other unitary transformations, such as permutations, preserves the number of zeros and allow to divide the possible patterns of y d and y u into equivalence classes with identical predictions. It is always possible to enforce texture zeros in arbitrary entries of the fermion mass matrices by means of Abelian symmetries [84].
Special attention has been paid to hermitian textures, y u,d = y † u,d . One such texture is called n-zero if y u and y d have a total of n zeros, the off-diagonal ones counting a half. A typical texture is where D u = D d = 0 in the original proposal by Fritzsch. The 6-zero Fritzsch texture is ruled out since it predicts a too large |V cb |. More general hermitian 6-zero and 5-zero textures have been analyzed under the assumption of hierarchy among the matrix elements [85]. All possible combinations of hermitian 6-zero and 5-zero textures are ruled out by now, 2 because they cannot reproduce |V cb |, with the only exception of y f in (50) with D u = 0 and D d = 0, which is marginally compatible with present data [87]. The 4-zero texture in (50) with both D u and D d non-vanishing is still allowed, provided the hierarchy in the 23 block is mild [88,89]. In particular |B f /C f | should be O(0.1). In this case cancellations among the matrix elements are needed to get m s /m b , m c /m t and |V cb | in the correct range. Notice that such cancellations were excluded by assumption in the case of y u,d of Eq. (46), where the hierarchy between quark masses of second and third generations was attributed to the smallness of . More general textures have been analyzed in Ref. [90]. If no symmetry requirement is imposed, there are viable textures that correctly fit the data. However, they do not show any predictive power with respect to any of the quark masses and mixing parameters. When y u,d are symmetric matrices, there are several independent 4-zero and 5-zero textures that allow to predict some of the light quark masses in terms of the mixing angles and the remaining masses.
In the lepton sector, working in the flavour basis where the charged lepton mass matrix is diagonal, we can study textures of the symmetric matrix for light Majorana neutrinos [91]. Those with three or more zeros are experimentally excluded, while data still allow seven independent 2-zero textures [92,93]. The requirement that zeros of the neutrino mass matrix should be present in the flavour basis can be relaxed and more general textures have been investigated. Both cases of Dirac and Majorana neutrinos have been analyzed and many independent textures in agreement with the existing data have been identified. For a systematic study of all possibilities of texture zeros in the charged-lepton and neutrino mass matrices and for a recent review on the subject we refer the reader to Ref. [94] and Ref. [95], respectively.
Other choices of G f have been considered. For instance Refs. [96][97][98] study the case of G f = U(2) 3 , as an alternative to MFV to suppress FCNC in supersymmetric extensions of the SM. In that context, the advantage of G f = U(2) 3 is that squarks of the first two generations can be taken significantly heavier than the third generation ones. Furthermore the presence of separate U(2) for left and right-handed fields, provides a sufficient protection of flavour-violating effects in the right-handed sector.

Discrete symmetries
The data from neutrino oscillations before 2012 were supporting flavour symmetries, especially through the indication of a vanishing reactor angle θ 13 and a maximal atmospheric mixing angle θ 23 , features that are difficult to attribute to an underlying theory based on pure chance. Today we know with accuracy that θ 13 is neither vanishing nor particularly small, its size being comparable to that of the Cabibbo angle.
Recent global fits [99][100][101] (see Table 4) favour a deviation of the atmospheric mixing angle from the maximal value by several degrees and, at the same time, show a preference for a maximal Dirac CP-violating phase δ CP , though the whole range from 0 to 2π is still allowed at 3σ . These two features are closely related, since they are mainly driven by the direct comparison between P ee measured by reactor experiments, which essentially determines θ 13 , and the conversion probability P μe measured by MINOS and T2K, which in turn is sensitive to a combination of θ 13 , θ 23 and δ CP . There is not yet a clear indication of the octant θ 23 falls in. Furthermore the most precise single experimental determination of the atmospheric angle [102], θ 23 = 0.514 +0.055 −0.056 (0.511 +0.055 −0.055 ) for NO (IO), is still compatible with π/4. Probably it is too early to conclude that a maximal θ 23 is ruled out by data.
There are few known mechanisms to generate a maximal atmospheric angle. We know that θ 23 cannot be made naturally maximal by renormalisation group evolution, barring a fine-tuning of the initial conditions or ad hoc threshold effects [103][104][105][106]. Moreover, in the context of flavour symmetries, θ 23 is completely determined by breaking effects, if we accept that e and μ masses vanish in the limit of exact symmetry [107]. Indeed, if the masses of e, μ and τ are all non-vanishing when the symmetry is exact, then large O(1) breaking terms are needed to explain the relative hierarchy among charged fermion masses. Thus it is more natural to assume that at least m e and m μ are proportional to small symmetry breaking parameters. In this case, when the symmetry is exact, the lepton mixing matrix is determined up to an arbitrary rotation coming from the eμ charged lepton sector and the atmospheric mixing angle can only be determined when the symmetry breaking parameters are turned on. We have no examples of a maximal θ 23 from the breaking of an abelian symmetry. If a nearly maximal atmospheric angle is not due to pure chance, we are left with broken non-abelian symmetries.
Before the measurement of θ 13 a particularly attractive lepton mixing pattern was the tri-bi-maximal one The good agreement between TB mixing and pre-2012 data strongly supported the idea that the true mixing matrix could be described in terms of small corrections to a LO mixing matrix U 0 PMNS , which could be derived from symmetry considerations. The simplest way to reproduce the TB mixing pattern is by exploiting discrete flavour symmetries [108][109][110][111][112]. The theory is invariant under a discrete flavour symmetry G f , broken down in such a way that neutrino and charged-lepton sectors have different residual symmetries G ν and G e , at least in a LO approximation where small effects are neglected. If neutrinos are of Majorana type, the most general group leaving m ν invariant and the individual masses m i unconstrained is Z 2 × Z 2 , a finite group. The subgroup G e can be continuous, but G e discrete is the simplest option. We require a sufficiently large G e to distinguish the three charged leptons. For instance we can choose G e = Z n (n ≥ 3) or G e = Z 2 × Z 2 . Once G e and G ν have been chosen inside G f , the embedding automatically fixes the relative alignment of m † l m l and m ν in flavour space. Lepton masses are unconstrained but U 0 PMNS is determined up to Majorana phases and up to permutations of rows and columns. This freedom apart, this setup predicts the three mixing angles θ 0 i j and the Dirac phase δ 0 CP . In most concrete models, where symmetry breaking is achieved via VEVs of a set of flavons ϕ, the LO results are modified by small corrections of order u = ϕ / f , Before 2012, in the specific case U 0 PMNS = U TB these corrections were expected to be very small, of the order of few percent [113,114], not to spoil the good agreement in the predicted value of the solar mixing angle. On this basis the simplest models reproducing U TB at the LO predicted θ 13 not larger than few degrees, now proven to be wrong by experiments. Discrete flavour symmetries can also be extended to quarks and even incorporated in GUTs, but in the existing constructions the symmetry has to be badly broken in the quark sector. Discrete flavour symmetries are also relevant in the so-called indirect models [110]. In this case the breaking of G f leaves no residual symmetries and the role of the flavour group is mainly to get specific vacuum alignments of the scalar fields that control fermion masses.
Several modifications of the simplest models based on discrete symmetries have been proposed to match the most recent data. If we keep adopting U 0 PMNS = U TB as LO approximation, an economic way to reproduce the actual value of θ 13 is to introduce large correction terms, O(u) ≈ 0.2. This is also viable in other schemes where U 0 PMNS differs substantially from U TB , such as the so-called bi-maximal mixing (BM), To first order in α we have By eliminating (α, δ) we get a relation between sin 2 θ 12 , sin θ 13 , and δ CP , plotted in Fig. 5. This model predicts δ CP close to π in order to reproduce correctly sin 2 θ 12 , as can be seen from Eq. (55).
Another possibility is to relax the symmetry requirements. S 4 is the smallest group reproducing TB mixing through the breaking down to G e = Z 3 and G ν = Z 2 × Z 2 , whose generators are T and (S, U ), respectively 3 . In the basis where T and the charged leptons are diagonal, the element U coincides with the so-called μτ exchange symmetry [118][119][120][121][122][123], directly responsible for the vanishing of θ 13 and for θ 23 being maximal. We can avoid having θ 13 = 0 and θ 23 = π/4 if G ν is a single Z 2 subgroup generated either by the element S or by the element SU . When the preserved parity is S, the mixing pattern, TM 2 , is tri-maximal and corresponds to with both α and δ unconstrained. When the preserved parity is SU, the mixing pattern, TM 1 , is also of tri-maximal type and is given by U TM 1 = U TB U 23 (α, δ), where U 23 (α, δ) is the transformation analogous to U 13 (α, δ), acting in the 23 plane. The mixing angles and the Dirac phase are predicted in terms of (α, δ) and we get two relations among physical quantities, shown in Table 5 [124][125][126][127][128][129][130][131]. The interesting feature of these relations is that the predicted deviations from TB are linear in sin θ 13 for sin 2 θ 23 , and quadratic for sin 2 θ 12 , known with much better precision. One of these relations is plotted in Fig. 5 in the case of TM 2 , from which we see that a substantial improvement in the data is needed to test this possibility. Explicit models based on A 4 realizing the TM 2 breaking pattern were indeed proposed before the measurement of θ 13 [132]. The possibility of reducing the residual symmetry G ν to Z 2 can be systematically investigated [133]. A further possibility is to look for alternative LO approximations where θ 13 is closer to the measured value. Remarkably, an infinite set of groups G f giving rise to LO approximations closer to the data has been found. Of particular interest is the special form of tri-maximal mixing TM 2 in (56), with both α and δ quantised, reproduced by groups of the series (6n 2 ) [134][135][136][137]. For example, choosing n = (4, 8, 10) we have α = (±π/12, ±π/24, ±π/15) and sin 2 θ 0 13 = (0.045, 0.011, 0.029). The Dirac phase is zero (modulo π ). In Fig. 6 the values of |U e3 | are plotted versus n [137]. Other discrete groups have been investigated in Ref. [138]. Very remarkably, a complete classification of all possible mixing matrices |U 0 PMNS | generated from any finite group has been recently carried out in Ref. [139].
Another development consists in combining discrete and CP symmetries [140,141] and exploring the symmetry breaking patterns such a combination can give rise to. A wellknown example is that of the so-called μτ reflection symmetry [142][143][144][145] (not to be confused with the μτ exchange symmetry), which exchanges a muon (tau) neutrino with a tau (muon) antineutrino in the charged-lepton mass basis. If such a symmetry is imposed, the atmospheric mixing angle is predicted to be maximal, while θ 13 is in general non-vanishing for a maximal Dirac phase δ and the Majorana phases vanish. The solar mixing angle remains unconstrained.
A general formalism which combines CP and flavour symmetries [140] can be used to constrain the lepton mixing matrix. A theory symmetric under CP and under a discrete flavour group G f is assumed to have residual symmetries G e , generated by some elements Q i and G ν = Z 2 × CP, generated by a parity Z and a CP transformation X . The action of X in flavour space can be non-trivial [146] and should respect a set of consistency conditions [140,[147][148][149]. The residual symmetries G e and G ν imply the following conditions on m † l m l and m ν : These conditions are strong enough to determine U 0 PMNS completely, up to one real parameter θ , ranging from 0 to π : Mixing angles and phases, both Dirac and Majorana, are then predicted as a function of θ , modulo the ambiguity related to the freedom of permuting rows and columns and to the intrinsic parity of neutrinos. The formalism is completely invariant under any change of basis in field space. The physical results only depend on the initial symmetry and the residual symmetries specified by (Q i , Z , X ). An interesting example is provided by G f = S 4 . An exhaustive analysis has been presented in Ref. [140]. The residual symmetries can be chosen as G e = Z 3 , generated by the element T , and G ν = Z 2 × CP, generated by (Z , X ). The parity transformation Z can be either S (case I) or SU (case IV) and a consistent CP transformation X acting on the lepton doublets coincides with the μ − τ reflection symmetry in the basis where T (and the combination m † l m l ) is diagonal. Thus the predicted mixing pattern has a maximal atmospheric mixing angle, a maximal Dirac phase, vanishing Majorana phases and there is a relation between the solar angle and the reactor angle, shown in Fig. 7.

Hints from empirical relations
Empirical relations among fermion masses and/or mixing angles have been frequently suggested as a clue towards a solution of the flavour puzzle. Here as an example I will comment one of the most striking ones, Koide's formula 4 for charged-lepton masses [172,173]: The experimental values of the pole masses m e , m μ , m τ are [174]: framework. An unsatisfactory feature of the formula is that it requires as inputs pole masses and not running masses, which would be more adequate if masses are believed to originate from a common scale f . This can be seen in Table 6 where the running Yukawa parameters of the charged leptons in the MS scheme are listed for several values of the scale μ. The fourth row shows the value of the τ Yukawa coupling y K τ (μ) as derived from the Koide formula using as inputs y e (μ) and y μ (μ). In the range μ = 0.1 ÷ 10 TeV the mismatch between y K τ (μ) and y τ (μ) is close to 1 %, while the accuracy with which y τ (μ) is estimated in Ref. [175] ranges from 10 −4 to 10 −3 , thus making the discrepancy significant. For example the value of the tau mass at the scale M Z predicted by the Koide relation is smaller by about 20 MeV, while the experimental precision on the same parameter is close to 0.2 MeV. This gap mainly arises as the effect of the pure QED running going from m e to m τ as can be seen from the leading-order relation between pole masses and MS masses: At energies higher than 1 ÷ 10 TeV, the success of the Koide relation depends on unknown physics. If the SM merges into a supersymmetric theory, the evolution of the Yukawa couplings depends on additional parameters, like the superparticle thresholds and tan β. Just above the superpartners mass threshold the Yukawa coupling are known to a much worse precision, around 1 % [175], and the Koide relation might again be compatible with the extrapolated data. At the GUT scale, larger uncertainties are induced on y e,μ,τ (M GUT ) by tan β. More in general, we can ask what is the probability of finding a simple homogeneous relation among charged fermion masses holding to an accuracy similar to the one of the Koide relation at a given scale μ. It is difficult to answer quantitatively this question, but I personally think that such a probability is relatively high. As an example consider the following relation, valid for charged fermion masses rather than for their squared roots: This relation, symmetric under permutations of the flavour labels, produces the outputs in the fifth row of Table 6. We see that the predictions of y τ (μ) are as good as the one from the Koide relation. Probably it is not surprising that scanning thousands of possibilities we can find "simple" relations working at the level of the percent accuracy. It is much more difficult to find, through these relations, a direct link to some unknown fundamental layer of particle physics. Table 6 Values of the running SM Yukawa couplings for the charged fermions in the MS scheme, at selected renormalisation scales μ, from Ref. [175]. The fourth and the fifth rows show the prediction of the mass relations in Eqs. (59) and (62), respectively, using as inputs the first two rows, y e (μ) and y μ (μ). The errors indicated in brackets affect the last digit and are dominated by the uncertainty in y μ (μ)

Conclusion
We are witnessing a continuous experimental progress in flavour physics. In neutrino physics, squared-mass differences and mixing angles are known to an accuracy that approaches the percent level. The reactor angle θ 13 is away from zero by many standard deviations. For the first time global fits hint at a non-trivial Dirac phase. While the new data have been effective in ruling out many models of fermion masses and mixing angles, as a matter of fact no compelling and unique theoretical picture has emerged so far. Present data can still be described within widely different frameworks. Based on our experience with gauge interactions we might hope that the flavour sector becomes simple and symmetric at a high energy scale, with a small number of relevant parameters providing a complete description. It is fair to say that we have not been able to identify a clear symmetry pattern from data so far. Before 2012 discrete symmetries were considered as a good candidate. In particular those predicting a nearly tri-bi-maximal mixing were favoured by data, but the prediction of θ 13 turned out to be wrong. The evidence for discrete symmetries in the quark sector is very poor and in a unified description of all fermions this kind of symmetry is typically badly broken in the quark sector. The whole approach is too much centered on the lepton mixing properties, while a description of the fermion masses seems to need additional ingredients. Several modifications of the simplest schemes to accommodate the present data are still possible and have the advantage of being quantitatively testable. But the real open question is whether a non-trivial implementation of discrete symmetries exists encompassing quark and lepton sectors in a unified picture and providing a quantitative description of both masses and mixing angles.
There are other models where fermion masses and mixing angles are mapped into a large number of irreducible and unconstrained order-one parameters, thus incarnating the Anarchy idea. For their intrinsic nature these models essentially escape experimental tests going beyond the order-ofmagnitude accuracy. However, we cannot fail to be impressed by the fact that they can provide a common description to both fermion masses and mixing angles, that they are compati-ble with grand unified theories and that they can be derived within widely different theoretical frameworks. The fact that this kind of models can be implemented even in a highly constrained setup such as an SO(10) grand unified theory is really remarkable. As a drawback, in these models the bounds on the scale of new flavoured physics is typically pushed above the 10 TeV scale, reducing the possibility of testing these ideas.
Flavour remains a fascinating mystery, still eluding all our attempts to find the rationale underlying our observations. Has this puzzle any solution? Are we misled by the questions we have formulated so far? Will we ever have access to the flavour scale? Man has long been fascinated by the mystery of planet motion. Surprisingly precise measurements and accurate predictions already existed in remote epochs. For a long time the most reliable models were based on the special character of geometrical figures like the sphere or the circle. Attempts to explain the relative sizes of the solar planetary orbits revealed themselves misleading. More accurate observations, perseverance in identifying the correct questions and renunciation of old prejudices opened new perspectives to the scientific thought. Will this be the fate of the flavour puzzle too?