The Mystery of Neutrino Mixings

In the last years we have learnt a lot about neutrino masses and mixings. Neutrinos are not all massless but their masses are very small. Probably masses are small because neutrinos are Majorana particles with masses inversely proportional to the large scale M of lepton number (L) violation, which turns out to be compatible with the GUT scale. We have understood that there is no contradiction between large neutrino mixings and small quark mixings, even in the context of GUTs and that neutrino masses fit well in the SUSY GUT picture. Out of equilibrium decays with CP and L violation of heavy RH neutrinos can produce a B-L asymmetry, then converted near the weak scale by instantons into an amount of B asymmetry compatible with observations (baryogenesis via leptogenesis). It appears that active neutrinos are not a significant component of Dark Matter in the Universe. A long list of models have been formulated over the years to understand neutrino masses and mixings. With the continuous improvement of the data most of the models have been discarded by experiment. The surviving models still span a wide range going from a maximum of symmetry, with discrete non-abelian flavour groups, to the opposite extreme of anarchy.


Introduction
Experiments on neutrino oscillations, which measure differences of squared masses as well as mixing angles 1) 2) have established that neutrinos have a mass. Two distinct oscillation frequencies have been first measured in solar and atmospheric neutrino oscillations and later confirmed by experiments on earth, like KamLAND, K2K and MINOS. Two well separated differences need at least three different neutrino mass eigenstates involved in oscillations so that the three known neutrino species can be sufficient. Then at least two ν's must be massive while, in principle, the third one could still be massless. A signal corresponding to a third mass difference was claimed by the LSND experiment (with antineutrinos) but not confirmed by KARMEN. More recently MiniBooNE 3) has reported some possible supporting evidence for the LSND effect in their antineutrino run while no oscillation is observed in the neutrino run. The existence of a third oscillation frequency would imply the need for additional sterile neutrinos (i.e. with no weak interactions, as any new light active neutrino was excluded by LEP) or CPT violation (as, in this case, the masses of neutrinos and antineutrinos can be different).
The main recent developments on the experimental side were the results on θ 13 from T2K 4) and MINOS 5) (very recently also DOUBLE CHOOZ 6) ) and the coming back of sterile neutrinos. As well known, the T2K run was suddenly interrupted by the devastating earthquake that hit Japan on March 11, just minutes away from the scheduled presentation of the first T2K data. Later T2K released the first publication on their data 4) , reporting a 2.5σ signal for sin 2 2θ 13 that indicates a value of θ 13 close to the previous upper bound, of the order of the Cabibbo angle θ C .
On the evidence for sterile neutrinos a number of hints have been recently reported. They do not make yet a clear evidence but certainly pose an experimental problem that needs clarification. First, there is the MiniBooNE experiment 5) that in the antineutrino channel reports an excess of events supporting the LSND oscillation signal (originally observed with antineutrinos). More recently an update of the MiniBooNE data in the antineutrino channel shows less supporting evidence 7) . In the neutrino channel MiniBooNE did not observe a signal in the LSND domain. However, in these data there is a unexplained excess at low energy over the (reliably?) estimated background. In the neutrino data sample, for the search of a LSND-like signal, only the events with neutrino energy above a threshold value E th were used, leaving the issue of an explanation of the low energy excess unanswered. In the antineutrino channel most of the support to the LSND signal appears to arise from an excess above E th but quite close to it, so that there is, in my opinion, some room for perplexity. Then there is the reactor anomaly: a reevaluation of the reactor flux 8) produced an apparent gap between the theoretical expectations and the data taken at small distances from the reactor ( < ∼ 100 m). The discrepancy is of the same order of the quoted systematic error whose estimate, detailed in the paper, should perhaps be reconsidered. Similarly the Gallium anomaly 9) depends on the assumed cross-section which could be questioned. The reactor anomaly and the Gallium anomaly do not really agree on the oscillation parameters that they point to: the ∆m 2 values are compatible but the central values of sin 2 2θ differ by about an order of magnitude, with Gallium favouring the larger angle. Cosmological data allow the existence of one sterile neutrino, while the most stringent bounds arising form nucleosynthesis disfavour two or more sterile neutrinos 10) . Over all, only a small leakage from active to sterile neutrinos is allowed by present neutrino oscillation data, as discussed in refs. 11) If all the indications listed above were confirmed (it looks unlikely) then one sterile neutrino would not be enough and at least two would be needed with sub-eV masses. Establishing the existence of sterile neutrinos would be a great discovery. In fact a sterile neutrino is an exotic particle not predicted by the most popular models of new physics. A sterile neutrino is not a 4th generation neutrino: the latter is coupled to the weak interactions (it is active) and heavier than half the Z mass. A sterile neutrino would probably be a remnant of some hidden sector. The issue is very important so that new and better experimental data are badly needed.
As already mentioned, in neutrino oscillations the leakage from the three active species towards the sterile neutrinos is any case small and, in fact, the best established oscillation phenomena are well described in terms of 3-neutrino models. In this domain the main recent developments have been the T2K and MINOS results on θ 13 . The T2K result 4) , based on the observation of 6 electron events when 1.5 ± 0.3 are expected for θ 13 = 0, is converted into a confidence interval 0.03(0.04) ≤ sin 2 2θ 13 ≤ 0. 28(0.34) at 90% c.l. for sin 2 2θ 23 = 1, |∆m 2 | = 2.4 10 −3 eV 2 , δ CP = 0 and for normal (inverted) neutrino mass hierarchy. Also the MINOS Collaboration released 5) their corresponding 90% c.l. range as 0(0) ≤ sin 2 2θ 13 ≤ 0. 12(0.19), which is displaced towards smaller values with respect to that of T2K. Finally DOUBLE CHOOZ 6) finds (with only the far detector in operation): sin 2 2θ 13 = 0.085 ± 0.051 at 1σ.

Neutrino Masses and Lepton Number Violation
Neutrino oscillations imply non vanishing neutrino masses which in turn demand either the existence of right-handed (RH) neutrinos (Dirac masses) or lepton number L violation (Majorana masses) or both. Given that neutrino masses are certainly extremely small, it is really difficult from the theory point of view to avoid the conclusion that L conservation must be violated. In fact, in terms of lepton number violation the smallness of neutrino masses can be explained as inversely proportional to the very large scale where L is violated, of order M GU T or even M P l .
If L conservation is violated neutrinos can be Majorana fermions. For a Majorana neutrino each mass eigenstate with given helicity coincides with its own antiparticle with the same helicity. As well known, for a charged massive fermion there are four states differing by their charge and helicity (the four components of a Dirac spinor) as required by Lorentz and CPT invariance. For a massive Majorana neutrino, neutrinos and antineutrinos can be identified and only two components are needed to satisfy the Lorentz and CPT invariance constraints. Neutrinos can be Majorana fermions because, among the fundamental fermions (i.e. quarks and leptons), they are the only electrically neutral ones. If, and only if, the lepton number L is not conserved, i.e. it is not a good quantum number, then neutrinos and antineutrinos can be identified. For Majorana neutrinos both Dirac mass terms, that conserve L (ν → ν), and Majorana mass terms, that violate L by two units (ν →ν), are in principle possible. Of course the restrictions from gauge invariance must be respected. So for neutrinos the Dirac mass terms (ν R ν L +h.c.) arise from the couplings with the Higgs field, as for all quarks and leptons. For Majorana masses, a ν T L ν L mass term has weak isospin 1 and needs two Higgs fields to make an invariant. On the contrary a ν T R ν R mass term is a gauge singlet and needs no Higgs. As a consequence, the right-handed neutrino Majorana mass M R is not bound to be of the order of the electroweak symmetry breaking (induced by the Higgs vacuum expectation value) and can be very large (see below).
Some notation: the charge conjugated of ν is ν c , given by ν c = C(ν) T , where C = iγ 2 γ 0 is the charge conjugation matrix acting on the spinor indices (in the following, when dealing with the flavour structure of couplings, the C matrix will be often omitted but understood). In particular (ν c ) L = C(ν R ) T , so that, instead of using ν L and ν R , we can refer to ν L and (ν c ) L , or simply ν and ν c .
Once we accept L non-conservation we gain an elegant explanation for the smallness of neutrino masses. If L is not conserved, even in the absence of heavy RH neutrinos, Majorana masses for neutrinos can be generated by dimension five operators 12) of the form with H being the ordinary Higgs doublet, l i the SU(2) lepton doublets, λ a matrix in flavour space, Λ a large scale of mass, of order M GU T or M P l and a charge conjugation matrix C between the lepton fields is understood. Neutrino masses generated by O 5 are of the order m ν ≈ v 2 /Λ for λ ij ≈ O(1), where v ∼ O(100 GeV) is the vacuum expectation value of the ordinary Higgs. We consider that the existence of RH neutrinos ν c is quite plausible because most GUT groups larger than SU(5) require them. In particular the fact that ν c completes the representation 16 of SO(10): 16=5+10+1, so that all fermions of each family are contained in a single representation of the unifying group, is too impressive not to be significant. At least as a classification group SO(10) must be of some relevance in a more fundamental layer of the theory! Thus in the following we both assume that ν c exist and L is not conserved. With these assumptions the see-saw mechanism 13) is possible. We recall, also to fix notations, that in its simplest form it arises as follows. Consider the SU(3) × SU(2) × U(1) invariant Lagrangian giving rise to Dirac and ν c Majorana masses (for the time being we consider the ν (versus ν c ) Majorana mass terms as comparatively negligible): The Dirac mass matrix m D ≡ y ν v/ √ 2, originating from electroweak symmetry breaking, is, in general, non-hermitian and non-symmetric, while the Majorana mass matrix M is symmetric, M = M T . We expect the eigenvalues of M to be of order M GU T or more because ν c Majorana masses are SU(3)× SU(2)× U(1) invariant, hence unprotected and naturally of the order of the cutoff of the low-energy theory. Since all ν c are very heavy we can integrate them away. For this purpose we write down the equations of motion for ν c in the static limit, i.e. neglecting their kinetic terms: From this, by solving for ν c , we obtain: We now replace in the lagrangian, eq. (2), this expression for ν c and we get the operator O 5 of eq. (1) with and the resulting neutrino mass matrix reads: This is the well known see-saw mechanism result 13) : the light neutrino masses are quadratic in the Dirac masses and inversely proportional to the large Majorana mass. If some ν c are massless or light they would not be integrated away but simply added to the light neutrinos. Notice that the above results hold true for any number n of heavy neutral fermions R coupled to the 3 known neutrinos. In this more general case M is an n by n symmetric matrix and the coupling between heavy and light fields is described by the rectangular n by 3 matrix m D . Note that for m ν ≈ ∆m 2 atm ≈ 0.05 eV (see Table(1)) and m ν ≈ m 2 D /M with m D ≈ v ≈ 200 GeV we find M ≈ 10 15 GeV which indeed is an impressive indication for M GU T .
If additional non-renormalizable contributions to O 5 , eq. (1), are comparatively non-negligible, they should simply be added. For instance in SO (10) or in left-right extensions of the SM, an SU(2) L triplet can couple to two lepton doublets and to two Higgs and may induce a sizeable contribution to neutrino masses. At the level of the low-energy effective theory, such contribution is still described by the operator O 5 of eq. (1), obtained by integrating out the heavy SU(2) L triplet. This contribution is called type II to be distinguished from that obtained by the exchange of RH neutrinos (type I). One can also have the exchange of a fermionic SU(2) L triplet coupled to a lepton doublet and a Higgs (type III). After elimination of the heavy fields, at the level of the effective low-energy theory, the three types of see-saw terms are equivalent. In particular they have identical transformation properties under a chiral change of basis in flavour space. The difference is, however, that in type I see-saw mechanism, the Dirac matrix m D is presumably related to ordinary fermion masses because they are both generated by the Higgs mechanism and both must obey GUTinduced constraints. Thus more constraints are implied if one assumes the see-saw mechanism in its simplest type I version. .

Basic Formulae and Data for Three-Neutrino Mixing
We assume in the following that there are only two distinct neutrino oscillation frequencies, the atmospheric and the solar frequencies. These two can be reproduced with the known three light neutrino species (for more than three neutrinos see, for example, ref. 11) ).
Neutrino oscillations are due to a misalignment between the flavour basis, ν ′ ≡ (ν e , ν µ , ν τ ), where ν e is the partner of the mass and flavour eigenstate e − in a lefthanded (LH) weak isospin SU(2) doublet (similarly for ν µ and ν τ ) and the mass eigenstates ν ≡ (ν 1 , ν 2 , ν 3 ) 14,15) : where U is the unitary 3 by 3 mixing matrix. Given the definition of U and the transformation properties of the effective light neutrino mass matrix m ν in eq. (1): we obtain the general form of m ν (i.e. of the light ν mass matrix in the basis where the charged lepton mass is a diagonal matrix): The matrix U can be parameterized in terms of three mixing angles θ 12 , θ 23 and θ 13 (0 ≤ θ ij ≤ π/2) and one phase ϕ (0 ≤ ϕ ≤ 2π) 16) , exactly as for the quark mixing matrix V CKM . The following definition of mixing angles can be adopted: where s ij ≡ sin θ ij , c ij ≡ cos θ ij . In addition, if ν are Majorana particles, we have the relative phases among the Majorana masses m 1 , m 2 and m 3 . If we choose m 3 real and positive, these phases are carried by m 1,2 ≡ |m 1,2 |e iφ 1,2 17) . Thus, in general, 9 parameters are added to the SM when non-vanishing neutrino masses are included: 3 eigenvalues, 3 mixing angles and 3 CP violating phases. In our notation the two frequencies, ∆m 2 I /4E (I=sun,atm), are parametrized in terms of the ν mass eigenvalues by where ∆m 2 12 = |m 2 | 2 − |m 1 | 2 > 0 and ∆m 2 23 = m 2 3 − |m 2 | 2 . The numbering 1,2,3 corresponds to our definition of the frequencies and in principle may not coincide with the ordering from the lightest to the heaviest state. In fact, the sign of ∆m 2 23 is not known and its determination is one of the existing experimental challenges. A positive (negative) sign corresponds to normal (inverse) hierarchy.
With the above definitions the present data are summarised in Table( Oscillation experiments do not provide information about the absolute neutrino mass scale. Limits on that are obtained 1) from the endpoint of the tritium beta decay spectrum, from cosmology and from neutrinoless double beta decay (0νββ). From tritium we have an absolute upper limit of 2.2 eV (at 95% C.L.) 20) on the antineutrino mass eigenvalues involved in beta decay, which, combined with the observed oscillation frequencies under the assumption of three CPT-invariant light neutrinos, also amounts to an upper bound on the masses of the other active neutrinos. Complementary information on the sum of neutrino masses is also provided by the galaxy power spectrum combined with measurements of the cosmic microwave background anisotropies. According to recent analyses of the most reliable data 21) one obtains i |m i | < 0.60 ÷ 0.75 eV (at 95% C.L.), depending on the retained data and the cosmological model priors assumed. These numbers for the sum have to be divided by 3 in order to obtain a limit on the mass of each light neutrino. The discovery of 0νββ decay would be very important because it would establish lepton number violation and the Majorana nature of ν's, and provide direct information on the absolute scale of neutrino masses. The present limit from 0νββ (with large ambiguities from nuclear matrix elements) is about |m ee | < (0.3 ÷ 0.8) eV 22) (see eq. (13)).
By now, after KamLAND, SNO and the upper limits on the absolute value of neutrino masses, not too much hierarchy in the spectrum of neutrinos is indicated by experiments 18) , 19) : Thus, for a hierarchical spectrum, m 2 /m 3 ∼ √ r ∼ 0.2, which is comparable to λ C ∼ 0.226 (thoughout this article λ C = sin θ C , with θ C being the Cabibbo angle) or m µ /m τ ∼ 0.24. This suggests that the same hierarchy parameter (raised to powers with O(1) exponents) may apply for quark, charged lepton and neutrino mass matrices. In fact, m µ /m τ ∼ 0.06 ∼ λ 2 C and m e /m µ ∼ 0.005 ∼ λ 3−4 C ). For the near future the most important experimental challenges on neutrino oscillation experiments are more precise measurements of the absolute scale of neutrino mass (KATRIN, MARE), the accurate determination of θ 13 (from MINOS, T2K and the reactor experiments DOUBLE CHOOZ, Daya Bay and RENO) and of the shift from maximal of θ 23 , the fixing of the sign of ∆m 2 23 (normal or inverse hierarchy) (e.g. NOνA), the detection of CP violation in ν oscillations. Related to neutrino physics is the issue of the non conservation of the separate e, µ and τ lepton numbers. The recent new limit Br(µ → eγ) < ∼ 2.4 . 10 −12 obtained by the MEG experiment 23) is largely satisfied in the SM but it imposes a strong constraint on SUSY-GUT models.

Importance of Neutrinoless Double Beta Decay
The detection of neutrino-less double beta decay 24) would provide direct evidence of L non conservation and of the Majorana nature of neutrinos. It would also offer a way to possibly disentangle the 3 cases of degenerate, normal or inverse hierachy neutrino spectrum. The quantity which is bound by experiments on 0νββ is the 11 entry of the ν mass matrix, which in general, from m ν = U * m diag U † , is given by : where m 1,2 are complex masses (including Majorana phases) while m 3 can be taken as real and positive and φ is the U P M N S phase measurable from CP violation in oscillation experiments. Starting from this general formula it is simple to derive the bounds for degenerate, inverse hierarchy or normal hierarchy mass patterns shown in Fig.1 25) .
In the next few years a new generation of experiments (CUORE, GERDA, ....) will reach a larger sensitivity on 0νββ by about an order of magnitude. Assuming the standard mechanism through mediation of a light massive Majorana neutrino, if these experiments will observe a signal this would indicate that the inverse hierarchy is realized, if not, then the normal hierarchy case remains a possibility.

Baryogenesis via Leptogenesis from Heavy ν c Decay
In the Universe we observe an apparent excess of baryons over antibaryons. It is appealing that one can explain the observed baryon asymmetry by dynamical evolution (baryogenesis) starting from an initial state of the Universe with zero baryon number. For baryogenesis one needs the three famous Sakharov conditions: B violation, CP violation and no thermal equilibrium. In the history of the Universe these necessary requirements have possibly occurred at different epochs. Note however that the asymmetry generated during one such epoch could be erased in following epochs if not protected by some dynamical reason. In principle these conditions could be fulfilled in the SM at the electroweak phase transition. In fact, when kT is of the order of a few TeV, B conservation is violated by instantons (but B-L is conserved), CP symmetry is violated by the CKM phase and sufficiently marked out-of-equilibrium conditions could be realized during the electroweak phase transition. So the conditions for baryogenesis at the weak scale in the SM superficially appear to be present. However, a more quantitative analysis 26) shows that baryogenesis is not possible in the SM because there is not enough CP violation and the phase transition is not sufficiently strong first order, unless the Higgs mass is below a bound which by now is excluded by LEP. In SUSY extensions of the SM, in particular in the MSSM, there are additional sources of CP violation and the bound on m H is modified but also this possibility has by now become at best marginal after the results from LEP2.
If baryogenesis at the weak scale is excluded by the data it can occur at or just below the GUT scale, after inflation. But only that part with |B − L| > 0 would survive and not be erased at the weak scale by instanton effects. Thus baryogenesis at kT ∼ 10 10 − 10 15 GeV needs B-L violation and this is also needed to allow m ν if neutrinos are Majorana particles. The two effects could be related if baryogenesis arises from leptogenesis then converted into baryogenesis by instantons 27) . The

Models of Neutrino Mixing
Neutrino mixing is important because it could in principle provide new clues for the understanding of the flavour problem. Even more so since neutrino mixing angles show a pattern that is completely different than that of quark mixing: for quarks all mixing angles are small, for neutrinos two angles are large (one is even compatible with the maximal value) and only the third one is small. We first consider the case of models based on discrete flavour groups that have received a lot of attention in recent years 28) . There are a number of special mixing patterns that have been studied in this context. These mixing matrices all have sin 2 θ 23 = 1/2, sin 2 θ 13 = 0 and differ by the value of sin 2 θ 12 (see Fig. 2). The corresponding mass matrices are 2-3 symmetric , i.e. µ − τ symmetric (see, as examples, the early work in ref. 29) and the recent paper ref. 30) ). The observed value of sin 2 θ 12 18) , 19) the best measured mixing angle, is very close, from below, to the so called Tri-Bimaximal (TB) value 31) which is sin 2 θ 12 = 1/3. Alternatively it is also very close, from above, to the Golden Ratio (GR) value 32) , 33) , 34) which is sin is the GR (for a different connection to the GR in this context, see 35) , 36) ). Thus, a possibility is that one or the other of these coincidences is taken seriously and this leads to models where either TB or GR mixing is naturally predicted as a good first approximation. Here I will mainly refer to TB mixing which is the simplest and most studied first approximation to the data.
The TB mixing matrix (in a particular phase convention) is given by: Note that the mixing angles are independent of mass ratios (while, for quark mixings, relations like λ 2 C ∼ m d /m s are typical). In the basis where charged lepton masses are diagonal, the effective neutrino mass matrix in the TB case is given by where: The eigenvalues of m ν are m 1 , m 2 , m 3 with eigenvectors (−2, 1, 1)/ √ 6, (1, 1, 1)/ √ 3 and (0, 1, −1)/ √ 2, respectively. The expression in eq.(15) can be reproduced in models with sequential dominance or with form dominance, discussed by S. King and collaborators 37) .
As we see the most general neutrino mass matrix corresponding to TB mixing, in the basis of diagonal charged leptons, is of the form: This is a symmetric, 2-3 symmetric matrix with a 11 + a 12 = a 22 + a 23 . We now discuss models that naturally produce TB mixing in first approximation. Discrete non-abelian groups naturally emerge as suitable flavour symmetries 28) . In fact the TB mixing matrix immediately suggests rotations by fixed, discrete angles. In a series of papers, started by 38) (for a rather complete list of references see 28) , some recent papers are in ref. 39) ) it has been pointed out that a broken flavour symmetry based on the discrete group A 4 appears to provide a simplest realization of this specific mixing pattern in Leading Order (LO). We recall that A n is the group of even permutations of n objects (n!/2 elements). In the case of GR mixing the simplest choice is the group A 5 34) . Other solutions for TB mixing based on alternative discrete or continuous flavour groups have also been considered 28) , 40) but the A 4 models have a very economical and attractive structure, e.g. in terms of group representations and of field content. We recall that A 4 can be generated by the two elements S and T obeying the relations (a "presentation" of the group): The 12  Note that the squares of the dimensions of all these representations add up to 12, the dimension of A 4 . It is immediate to see that one-dimensional unitary representations are given by: The three-dimensional unitary representation, in a basis where the element T is diagonal, is given by: Note that the generic mass matrix for TB mixing in eq. (17) can be specified as the most general matrix that is invariant under µ − τ symmetry, implemented by the unitary matrix A µτ : and under the S transformation: where S is given in eq. (20). The m mass matrix of eq. (17) is derived in the basis where charged leptons are diagonal. It is useful to consider the product m 2 = m † e m e , where m e is the charged lepton mass matrix (defined as ψ R m e ψ L ), because this product transforms as m ′2 = U † e m 2 U e , with U e the unitary matrix that rotates the left-handed (LH) charged lepton fields. The most general diagonal m 2 is invariant under a diagonal phase matrix with 3 different phase factors: and conversely a matrix m † e m e satisfying the above requirement is diagonal. If T n = 1 the matrix T generates a cyclic group Z n . The simplest case is n = 3, which corresponds to Z 3 (but n > 3 is equally possible) and to the T matrix in eq. (20).
We can now see why A 4 works for TB mixing. It works because S and T are matrices of A 4 (in fact they satisfy eqs. (18)). One could object that the matrix A 23 is not an element of A 4 (because the 2-3 exchange is an odd permutation). But it can be shown that in A 4 models the 2-3 symmetry is maintained by imposing that there are no flavons transforming as 1 ′ or 1 ′′ that break A 4 with two different VEV's: in particular one can assume that there are no flavons in the model transforming as 1 ′ or 1 ′′ 41) .
The group A 4 has two obvious subgroups: G S , which is a reflection subgroup generated by S and G T , which is the group generated by T , which is isomorphic to Z 3 . If the flavour symmetry associated to A 4 is broken by the VEV of a triplet ϕ = (ϕ 1 , ϕ 2 , ϕ 3 ) of scalar fields, there are two interesting breaking pattern. The VEV breaks breaks A 4 down to G T . We have seen that G S and G T are the relevant low-energy symmetries of the neutrino and the charged-lepton sectors, respectively. Indeed we have shown that the TB mass matrix is invariant under G S and, for charged leptons, a diagonal m † e m e is invariant under G T . A crucial part of all serious A4 models is the dynamical generation of this alignment in a natural way. In most of the models A 4 is accompanied by additional flavour symmetries, either discrete like Z N or continuous like U(1), which are necessary to eliminate unwanted couplings, to ensure the needed vacuum alignment (obtained from the minimization of the most general potential compatible with the assumed symmetries) and to reproduce the observed mass hierarchies. Explicit realizations of models for TB mixing based on A 4 can be found, for example, in 41,42,43,44,45) . The possible origin of A 4 from a deeper level of the theory has been discussed in the context of extra dimensions and orbifolding 43) , 46) or as related to the fact that A 4 is a subgroup of the modular group 42) , which plays a role in string theory. In passing, we note that we are mainly interested here in the possibility that the flavour symmetry is valid at the GUT scale and is broken at lower scales.

Applying A 4 to lepton masses and mixings
In the lepton sector a typical A 4 model works as follows 42) . One assigns leptons to the four inequivalent representations of A 4 : LH lepton doublets l transform as a triplet 3, while the RH charged leptons e c , µ c and τ c transform as 1, 1 ′′ and 1 ′ , respectively. Here we consider a see-saw realization, so we also introduce conjugate neutrino fields ν c transforming as a triplet of A 4 . The fact that LH lepton doublets l and, in the see-saw case, also the RH neutrinos ν c , transform as triplets is crucial to realize the fixed ratios of mass matrix elements needed to obtain TB mixing. A drawback is that for the ratio r, defined in eq.(12), one would expect √ r ≈ O(1) to be compared with the experimental value is √ r ≈ 0.2, which implies a moderate fine tuning.
One adopts a supersymmetric (SUSY) context also to make contact with Grand Unification (flavor symmetries are supposed to act near the GUT scale). In fact, as well known, SUSY is important in GUT's for offering a solution to the hierarchy problem, for improving coupling unification and for making the theory compatible with bounds on proton decay. But, in models of lepton mixing, SUSY also helps for obtaining the vacuum alignment, because the SUSY constraints are very strong and limit the form of the superpotential very much. Thus SUSY is not necessary but it is a plausible and useful ingredient. The flavor symmetry is broken by two triplets ϕ S and ϕ T (with the vacuum alignment in eqs. (24,25)) and by one or more singlets ξ. All these fields are invariant under the SM gauge symmetry. Two Higgs doublets h u,d , invariant under A 4 , are also introduced. One can obtain the observed hierarchy among m e , m µ and m τ by introducing an additional U(1) F N flavor symmetry 47) under which only the RH lepton sector is charged (recently some models were proposed with a different VEV alignment such that the charged lepton hierarchies are obtained without introducing a U(1) symmetry 48,45) ). We recall that U(1) F N is a simplest flavor symmetry where particles in different generations are assigned (in general) different values of an Abelian charge. Also Higgs fields may get a non zero charge. When the symmetry is spontaneously broken the entries of mass matrices are suppressed if there is a charge mismatch and more so if the corresponding mismatch is larger. We assign FN-charges 0, q and 2q to τ c , µ c and e c , respectively. There is some freedom in the choice of q. Here we take q = 2. By assuming that a flavon θ, carrying a negative unit of FN charge, acquires a VEV θ /Λ ≡ λ < 1, the Yukawa couplings become field dependent quantities y e,µ,τ = y e,µ,τ (θ) and we have Had we chosen q = 1, we would have needed θ /Λ of order λ 2 , to reproduce the above result. The superpotential term for lepton masses, w l is given by: with dots denoting higher dimensional operators that lead to corrections to the LO approximation. In our notation, the product of 2 triplets (33) transforms as 1, (33) ′ transforms as 1 ′ and (33) ′′ transforms as 1 ′′ . To keep our formulae compact, we omit to write the Higgs and flavon fields h u,d , θ and the cut-off scale Λ. For instance y e e c (ϕ T l) stands for y e e c (ϕ T l)h d θ 4 /Λ 5 . The parameters of the superpotential w l are complex, in particular those responsible for the heavy neutrino Majorana masses, x A,B . Some terms allowed by the A 4 symmetry, such as the terms obtained by the exchange ϕ T ↔ ϕ S , (or the term (ν c ν c )) are missing in w l . Their absence is crucial and, in each version of A 4 models, is motivated by additional symmetries.
As for the neutrino spectrum both normal and inverted hierarchies can be realized. It is interesting that A 4 models with the see-saw mechanism typically lead to a light neutrino spectrum which satisfies the sum rule (among complex masses): 1 A detailed discussion of a spectrum of this type can be found in refs. 42,45) . The above sum rule gives rise to bounds on the lightest neutrino mass. As a consequence, for example, the possible values of |m ee | are restricted. For normal hierarchy we have while for inverted hierarchy In a completely general framework, without the restrictions imposed by the flavor symmetry, |m ee | could vanish in the case of normal hierarchy. In this model |m ee | is always different from zero, though its value for normal hierarchy is probably too small to be detected in the next generation of 0νββ experiments.
In the leading approximation A 4 models lead to exact TB mixing. In these models TB mixing is implied by the symmetry at the leading order approximation which is corrected by non leading effects. Given the set of flavour symmetries and having specified the field content, the non leading corrections to TB mixing, arising from higher dimensional effective operators, can be evaluated in a well defined expansion. In the absence of specific dynamical tricks, in a generic model, all the three mixing angles receive corrections of the same order of magnitude. Since the experimentally allowed departures of θ 12 from the TB value, sin 2 θ 12 = 1/3, are small, numerically not larger than O(λ 2 C ), it follows that both θ 13 and the deviation of θ 23 from the maximal value are expected to also be typically of the same general size. The central values sin θ 13 ∼ 0.16 − 0.11 that can be derived from the experimental results in the two columns of Table(

A 4 , quarks and GUT's
Much attention has been devoted to the question whether models with TB mixing in the neutrino sector can be suitably extended to also successfully describe the observed pattern of quark mixings and masses and whether this more complete framework can be made compatible with (supersymmetric) SU (5) or SO(10) Grand Unification.
The simplest attempts of directly extending models based on A 4 to quarks have not been satisfactory. At first sight the most appealing possibility is to adopt for quarks the same classification scheme under A 4 that one has used for leptons (see, for example, ref. 42) ). Thus one tentatively assumes that LH quark doublets Q transform as a triplet 3, while the antiquarks (u c , d c ), (c c , s c ) and (t c , b c ) transform as 1, 1 ′′ and 1 ′ , respectively. This leads to V u = V d and to the identity matrix for V CKM = V † u V d in the lowest approximation. This at first appears as very promising: a LO approximation where neutrino mixing is TB and V CKM = 1 is a very good starting point. But there are some problems. First, the corrections to V CKM = 1 turn out to be strongly constrained by the leptonic sector, because lepton mixing angles are very close to the TB values, and, in the simplest models, this constraint leads to a too small V us (i.e. the Cabibbo angle is rather large in comparison to the allowed shifts from the TB mixing angles). Also in these models, the quark classification which leads to V CKM = 1 is not compatible with A 4 commuting with SU (5). An additional consequence of the above assignment is that the top quark mass arises from a non-renormalizable dimension-5 operator. In that case, to reproduce the top mass, we need to compensate the cutoff suppression by some extra dynamical mechanism. Alternatively, we have to introduce a separate symmetry breaking parameter for the quark sector, sufficiently close to the cutoff scale.
Due to this, larger discrete groups have been considered for the description of quarks. A particularly appealing set of models is based on the discrete group T ′ , the double covering group of A 4 52) , 53) , 54) . The representations of T ′ are those of A 4 plus three independent doublets 2, 2 ′ and 2 ′′ . The doublets are interesting for the classification of the first two generations of quarks 55) . For example, in ref. 53) a viable description was obtained, i.e. in the leptonic sector the predictions of the A 4 model are maintained, while the T ′ symmetry plays an essential role for reproducing the pattern of quark mixing. But, again, the classification adopted in this model is not compatible with Grand Unification. As a result, the group A 4 was considered by many authors to be too limited to also describe quarks and to lead to a grand unified description. But it has been shown 44) that this negative attitude is not justified and that it is actually possible to construct a viable model based on A 4 which leads to a grand unified theory (GUT) of quarks and leptons with TB mixing for leptons and with quark (and charged lepton) masses and mixings compatible with experiment. At the same time this model offers an example of an extra dimensional SU(5) GUT in which a description of all fermion masses and mixings is accomplished. The formulation of SU (5) in extra dimensions has the usual advantages of avoiding large Higgs representations to break SU(5) and of solving the doublet-triplet splitting problem. The choice of the transformation properties of the two Higgses H 5 and H 5 has a special role in this model. They are chosen to transform as two different A 4 singlets 1 and 1 ′ . As a consequence, mass terms for the Higgs colour triplets are not directly allowed and their masses are introduced by orbifolding,à la Kawamura 56) . In this model, proton decay is dominated by gauge vector boson exchange giving rise to dimension-6 operators, while the usual contribution of dimension-5 operators is forbidden by the selection rules of the model. Given the large M GU T scale of SUSY models and the relatively huge theoretical uncertainties, the decay rate is within the present experimental limits. A see-saw realization in terms of an A 4 triplet of RH neutrinos ν c ensures the correct ratio of light neutrino masses with respect to the GUT scale. In this model extra dimensional effects directly contribute to determine the flavour pattern, in that the two lightest tenplets T 1 and T 2 are in the bulk (with a doubling T i and T ′ i , i = 1, 2 to ensure the correct zero mode spectrum), whereas the pentaplets F and T 3 are on the brane. The hierarchy of quark and charged lepton masses and of quark mixings is determined by a combination of extra dimensional suppression factors and of U(1) F N charges, both of which only apply to the first two generations, while the neutrino mixing angles derive from A 4 in the usual way. If the extra dimensional suppression factors and the U(1) F N charges are switched off, only the third generation masses of quarks and charged leptons survive. Thus the charged fermion mass matrices are nearly empty in this limit (not much of A 4 effects remain) and the quark mixing angles are determined by the small corrections induced by the above effects. The model is natural, since most of the small parameters in the observed pattern of masses and mixings as well as the necessary vacuum alignment are justified by the symmetries of the model. However, in this case, like in all models based on U(1) F N , the number of O(1) parameters is larger than the number of measurable quantities, so that in the quark sector the model can only account for the orders of magnitude (measured in terms of powers of an expansion parameter) and not for the exact values of mass ratios and mixing angles. A moderate fine tuning is only needed to enhance the Cabibbo mixing angle between the first two generations, which would generically be of O(λ 2 C ). The problem of constructing GUT models based on SU (5)  G f = P SL 2 (7) 60) . Clearly the case of SO (10) is even more difficult than that of SU(5) because the neutrino sector is tightly related to that of quarks and charged leptons as all belong to the 16 of SO (10). For a discussion of SO(10) ⊗A 4 models, see 61) . More in general see refs. 62) . In our opinion most of the models are incomplete (for example, the crucial issue of VEV alignment is not really treated in depth as it should) and/or involve a number of unjustified steps and ad-hoc fine tuning of parameters. While A 4 is the minimal flavor group leading to TB mixing, alternative flavor groups have been studied in the literature and can lead to interesting variants with some specific features. Actually, in ref. 63) , the claim was made that, in order to obtain the TB mixing "without fine tuning", the finite group must be S 4 or a larger group containing S 4 . For us this claim is not well grounded being based on an abstract mathematical criterium for a natural model (see also 64) ). For us a physical field theory model is natural if the interesting results are obtained from the most general lagrangian compatible with the stated symmetry and the specified representation content for the flavons. For example, we obtain from A 4 (which is a subgroup of S 4 ) a natural (in our sense) model for the TB mixing by simply not including symmetry breaking flavons transforming like the 1 ′ and the 1 ′′ representations of A 4 . This limitation on the transformation properties of the flavons is not allowed by the rules specified in ref. 63) which demand that the symmetry breaking is induced by all possible kinds of flavons (note that, according to this criterium, the SM of electroweak interactions would not be natural because only Higgs doublets are introduced!). Rather, for naturalness we also require that additional physical properties like the VEV alignment or the hierarchy of charged lepton masses also follow from the assumed symmetry and are not obtained by fine tuning parameters: for this actually A 4 can be more effective than S 4 because it possesses three different singlet representations 1, 1 ′ and 1 ′′ . Models of neutrino mixing based on S 4 have in fact been studied 65) . The group of the permutations of 4 objects S 4 has 24 elements and 5 equivalence classes that correspond to 5 inequivalent irreducible representations, two singlets, one doublet, two triplets: 1 1 , 1 2 , 2, 3 1 and 3 2 . Note that the squares of the dimensions of all these representations add up to 24.

Bimaximal Mixing and S4
The new results showing that probably θ 13 is not far from its former upper bound could alternatively be interpreted as an indication that the agreement with the TB or GR mixing is accidental. Then a scheme where instead the Bimaximal (BM) mixing is the correct first approximation modified by terms of O(λ C ) could be relevant. In BM mixing θ 12 and θ 23 are both maximal while θ 13 = 0 (see. Fig. 2). This is in line with the well known empirical observation that θ 12 + θ C ∼ π/4, a relation known as quark-lepton complementarity 66) , or similarly θ 12 + m µ /m τ ∼ π/4. No compelling model leading, without parameter fixing, to the exact complementarity relation has been produced so far. Probably the exact complementarity relation becomes more plausible if replaced with θ 12 + O(θ C ) ∼ π/4 or θ 12 + O(m µ /m τ ) ∼ π/4 (which we could call "weak" complementarity). One can think of models where, because of a suitable symmetry, BM mixing holds in the neutrino sector at leading order and the necessary, rather large, corrective terms for θ 12 arise from the diagonalization of charged lepton masses 66) . These terms of order O(λ C ) from the charged lepton sector would then generically also affect θ 13 and the resulting value could well be compatible with the present experimental values of θ 13 . A word of caution must be kept in mind: in the presence of these relatively large off diagonal terms in the charged lepton diagonalizing matrix one must arrange that not too large contributions to the decays µ → eγ or τ → µγ are generated 28) .
The BM mixing matrix is given by: Along this line of thought, we have used the expertise acquired with non Abelian finite flavour groups to construct a model 67) based on the permutation group S 4 which naturally leads to the BM mixing at leading order. We have adopted a supersymmetric formulation of the model in 4 space-time dimensions. The complete flavour group is S 4 × Z 4 × U(1) F N . In leading order, the charged leptons are diagonal and hierarchical and the light neutrino mass matrix, after see-saw, leads to the exact BM mixing. The model is built in such a way that the dominant corrections to the BM mixing, from higher dimensional operators in the superpotential, only arise from the charged lepton sector at next-to-the-leading-order and naturally inherit λ C (which fixes the charged lepton mass hierarchies) as the relevant expansion parameter. As a result the mixing angles deviate from the BM values by terms of O(λ C ) (at most), and weak complementarity holds. A crucial feature of the model is that only θ 12 and θ 13 are corrected by terms of O(λ C ) while θ 23 is unchanged at this order (which is essential for a better agreement of the model with the present data). Recently the model was extended to include quarks in a SU(5) Grand Unified version 68) or in a Pati-Salam framework 69) . An SO(10) model is discussed in ref. 70) .

Anarchy versus Symmetry
We now briefly turn to models that do not take seriously any of the coincidences described above (i.e. the proximity of the data to the TB or GR patterns or the quarklepton complementarity: these patterns cannot all be true and it is possible that none of them is true) and are therefore based on a less restrictive flavour symmetry. It is clear that the T2K hint that θ 13 may be large is great news for the most extreme position of this type, which is "anarchy" 71) : no symmetry at all in the lepton sector, only chance. This view predicts generic mixing angles, so the largest angle, θ 23 , should somewhat deviate from maximal and the smallest angle, θ 13 , should be as large as possible within the experimental bounds. Anarchy can be formulated in a SU(5) U(1) context by taking different Froggatt-Nielsen charges 47) only for the SU(5) tenplets (for example 10: (3,2,0) where 3 is the charge of the first generation, 2 of the second, zero of the third) while no charge differences appear in the5:5: (0,0,0). This assignment is in agreement with the empirical fact that the mass hierarchies are more pronounced for up quarks in comparison with down quarks and charged leptons. In a non see-saw model, with neutrino masses dominated by the contribution of the dimension-5 Weinberg operator in eq.(1), the5 vanishing charges directly lead to random neutrino mass and mixing matrices. In anarchical see-saw models also the charges of the SU(5) singlet RH neutrinos must be undifferentiated among the 3 generations: 1: (0,0,0). Anarchy can be mitigated by assuming that it only holds in the 2-3 sector: e.g5: (2,0,0) with the advantage that the first generation masses and the angle θ 13 are then naturally small (see ref. 72) for a recent discussion of this model). In models with see-saw one can alternatively play with the charges for the RH SU(5) singlet neutrinos. If, for example, we take 1: (1, -1, 0), together with5: (2,0,0), it is possible to get a normal hierarchy model with θ 13 small and also with r = ∆m 2 solar /∆m 2 atm naturally small (see, for example, ref . 73) ). In summary anarchy and its variants, all based on chance, offer a rather economical class of models that are among those that are compatible with the recent θ 13 results, with preference with the upper side of the experimentally allowed range.

Conclusion
In the last decade we have learnt a lot about neutrino masses and mixings. A list of important conclusions have been reached. Neutrinos are not all massless but their masses are very small. Probably masses are small because neutrinos are Majorana particles with masses inversely proportional to the large scale M of lepton number violation. It is quite remarkable that M is empirically not far from M GU T , so that neutrino masses fit well in the SUSY GUT picture. Also out of equilibrium decays with CP and L violation of heavy RH neutrinos can produce a B-L asymmetry, then converted near the weak scale by instantons into an amount of B asymmetry compatible with observations (baryogenesis via leptogenesis) 27) . It has been established that most probably active neutrinos are not a significant component of dark matter in the Universe. We have also understood there there is no contradiction between large neutrino mixings and small quark mixings, even in the context of GUTs. This is a very impressive list of achievements. Coming to a detailed analysis of neutrino masses and mixings a long collection of models have been formulated over the years. With continuous improvements of the data and more precise values of the mixing angles most of the models have been discarded by experiment. Still the surviving models span a wide range going from a maximum of symmetry, with discrete non-abelian flavour groups, to the opposite extreme of anarchy. By now, besides the detailed knowledge of the entries of the V CKM matrix we also have a reasonable determination of the neutrino mixing matrix U P M N S . The data appear to suggest some special patterns (recall Fig. 2) like TB or GR or BM mixing to be valid in some leading approximation, corrected by small non leading terms. If one takes these "coincidences" seriously, then non-abelian discrete flavour groups emerge as the main road to an understanding of this mixing pattern. Indeed the entries of e.g. TB mixing matrix are clearly suggestive of "rotations" by simple, very specific angles. It is remarkable that neutrino and quark mixings have such a different qualitative pattern. An obvious question is whether some additional indication for discrete flavour groups can be obtained by considering the extension of the models to the quark sector, perhaps in a Grand Unified context. The answer appears to be that, while the quark masses and mixings can indeed be reproduced in models where TB or BM mixing is realized in the leptonic sector through the action of discrete groups, there are no specific additional hints in favour of discrete groups that come from the quark sector. Further important input could come from µ → eγ and in general from lepton flavour violating processes, from b → sγ and from LHC physics. In fact, new physics at the weak scale could have important feedback on the physics of neutrino masses and mixing.
In conclusion, one could have imagined that neutrinos would bring a decisive boost towards the formulation of a comprehensive understanding of fermion masses and mixings. In reality it is frustrating that no real illumination was sparked on the problem of flavour. We can reproduce in many different ways the observations, in a wide range that goes from anarchy to discrete flavour symmetries) but we have not yet been able to single out a unique and convincing baseline for the understanding of fermion masses and mixings. In spite of many interesting ideas and the formulation of many elegant models the mysteries of the flavour structure of the three generations of fermions have not been much unveiled. I imagine that by the next edition of this by now classic School, we will know the value of θ 13 with a better accuracy, from the continuation of T2K, MINOS and DOUBLE CHOOZ and from the start of Daya Bay and RENO. Some existing models will be eliminated and the surviving ones will be updated to become more quantitative in order to cope with a precisely known mixing matrix. A definitely non vanishing θ 13 value will encourage the planning of long baseline experiments for the detection of CP violation in neutrino oscillations. Along the way the important issue of the existence of sterile neutrinos must be clarified. The on going or in preparation experiments on the absolute value of neutrino masses, on 0νββ, on µ → eγ, on the search for dark matter etc can also lead to extremely important developments in the near future. So this field is very promising and there are all reasons to expect an exciting time ahead of us.

Acknowledgments
I warmly thank Prof. Antonino Zichichi for his kind invitation and for the impres-