Constraining multi-Higgs flavour models

To study a flavour model with a non-minimal Higgs sector one must first define the symmetries of the fields; then identify what types of vacua exist and how they may break the symmetries; and finally determine whether the remnant symmetries are compatible with the experimental data. Here we address all these issues in the context of flavour models with any number of Higgs doublets. We stress the importance of analysing the Higgs vacuum expectation values that are pseudo-invariant under the generators of all subgroups. It is shown that the only way of obtaining a physical CKM mixing matrix and, simultaneously, non-degenerate and non-zero quark masses is requiring the vacuum expectation values of the Higgs fields to break completely the full flavour group, except possibly for some symmetry belonging to baryon number. The application of this technique to some illustrative examples, such as the flavour groups Delta(27), A4 and S3, is also presented.


Introduction
The origin of the flavour structure of the fermion masses and mixing in the standard model (SM) remains one of the unsolved puzzles in particle physics. Several approaches to this problem have been put forward, most of them based on the use of discrete [1] or continuous [2] flavour symmetries. In particular, mainly motivated by the measurement of several neutrino oscillation parameters, the use of discrete symmetries has recently become more popular (for recent reviews, see e.g. refs. [3]).
In flavour model building, one commonly chooses a flavour symmetry group K and then studies how fermion masses and mixing are constrained by this symmetry. Because of the SM gauge symmetry there is, of course, an additional global hypercharge symmetry U (1) Y in the Lagrangian. Thus, in such studies it is important to distinguish the chosen flavour group K from the full flavour group, denoted henceforth by G = K × U (1) Y . After spontaneous symmetry breaking (SSB) the vacuum expectation values (vev) of the Higgs fields break completely, in many instances, K and G, leaving no residual flavour symmetry in the model. But, in some rather interesting cases, the vevs could fully break K but not G, and thus some residual symmetry remains. Clearly, any accidental symmetry of the Lagrangian, such as U (1) B , corresponding to baryon number, is not relevant for the analysis of the flavour sector neither any residual symmetry of G which is contained in it.
In the literature, there are studies focusing on symmetry groups exclusively applied to the Higgs sector, either with two [4][5][6] or more than two doublets [7,8], as well as studies focusing either on flavour symmetry groups acting in the fermion sector with assumed vev alignments [3,9,10] or on the residual symmetries remaining on the fermion mass matrices [11][12][13]. Our aim is to provide a contribution towards the connection of all these subjects. In particular, we shall study in a systematic fashion the impact of a chosen flavour group on the Higgs potential and its vevs, on the Yukawa coupling matrices and thus on the residual physical properties of the mass matrices and the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix.
The introduction of additional scalars is a common feature to all flavour models, which brings with it a large number of possible model implementations for a given flavour group. Studying all these implementations, even for two or three Higgs doublet models, is not an easy task, specially for large flavour groups. We shall show that most of the unphysical scenarios can be identified just by analyzing how scalars transform under the flavour group and how they break it. Through this analysis we shall prove that under very general conditions some choices of flavour groups are excluded for a given scalar content. Although we shall restrict our discussion to the quark sector of the theory, the results presented here can be easily extended to the lepton sector.
The paper is organized as follows. In section 2, after briefly reviewing our notation and the experimental input, we analyze the invariance of the Yukawa coupling matrices under a full flavour group and then we proceed to the formulation of the main theorem and its proof. This theorem improves that in ref. [11], by including also the constraints arising from the right-handed quark sector. A summary of the possible CKM mixing patterns and quark spectra is presented, according to the existence or not of a residual symmetry in the whole quark sector. Most importantly, we address the problem of knowing whether such a residual symmetry is possible, by introducing the notion of pseudo-invariance of the vev under all subgroups of the flavour group. This is explained in detail in section 3, where we illustrate the technique by applying it to the flavour groups ∆(27), A 4 and S 3 , commonly found in the literature to explain fermion masses and mixing. Finally, our concluding remarks are given in section 4.

Notation and experimental input
Let us consider a model with N Higgs doublets, φ k , and three generations of left-handed quark doublets, Q L = (u L , d L ) T , right-handed up-type quark singlets, u R , and right-handed down-type quark singlets, d R . The Yukawa Lagrangian is where a sum over the Higgs fields (k = 1 . . . N ) is implicit, Γ k and ∆ k are 3 × 3 matrices in the down-and up-type quark family space, respectively. The scalar fields can be combined into Φ = (φ 1 , φ 2 , · · · , φ N ) T . After the spontaneous symmetry breaking, the scalars acquire vevs, Φ = (v 1 , v 2 , · · · , v N ) T , and the Yukawa Lagrangian contains the terms are the down-quark and up-quark mass matrices. Using transformations that keep the SU (2) L doublet structureQ L = ū L ,d L , it is possible to diagonalize M d or M u , but not both matrices. We denote by "left space", the space of vectorsQ L , which has three components in family space, one for each generation of quarks. So, M d and M u cannot be simultaneously diagonalized with left space transformations. To diagonalize them, we must treat u L and d L differently. Indeed, one can always find matrices Often, it is convenient to work with the Hermitian matrices We know from experiment that the quark masses are non-vanishing and non-degenerate. As a result, the matrices H d and H u are invertible and have a non-degenerate spectrum. The CKM mixing matrix is given by It is also known from experiment that V must have three non-trivial angles and that the largest contributor to CP violation in the kaon and B meson systems must be the CKM CP-violating phase, which is proportional to the basis invariant observable Therefore, V cannot be block-diagonal or any permutation of a block-diagonal matrix, it cannot be given by the identity or by some permutation matrix, and J must differ from zero.

Yukawa coupling invariance under a full flavour group G
Let us consider a group G with group elements g. The group can be defined by a set of generators i. For example, we may have two generators, a and b, such that all group elements are obtained by successive products of a and b. Then, one chooses how the various fields transform under the group G. This means that for three quark generations one chooses 3× 3 matrices representing this group, so that the quark fields transform under the generators of the group as where the index i runs from 1 to the number of group generators. Similarly, one chooses N × N matrices, such that the scalar fields transform as These equations hold for matrices G i representing the generators as well as for matrices representing any element of the group, i.e G g . The explicit form of G Li , G u Ri , G d Ri and G i dictates which representations D(Q L ), D(u R ), D(d R ) and D(Φ), respectively, we have chosen for the fields.
We now ask the Yukawa Lagrangian to be invariant under this group before SSB. Using eq. (2.1), this implies the following relations for the couplings: for all group elements g. Equation (2.11) is clearly valid for the identity element g = e. It is also valid for an element g such that G Lg = G Rg = e iθ ½ and G g = ½. The latter relations are obviously always verified since baryon number is conserved at or below the electroweak scale. We call such an element a trivial group element, the group it generates (G g ) a trivial group, and we denote this case by To decide what invariance may remain in the Yukawa Lagrangian after SSB (and, therefore, in the mass matrices), we must know how the vacuum transforms under the group elements. The ensuing analysis may not be straightforward, so we will start with a simple case, and generalize it in steps. For the moment, we make three simplifying assumptions: i) we ignore the up-type quarks; ii) we assume that Φ is in an irreducible representation (irrep) of G; iii) we assume that there is a non-trivial group element, g 1 , leaving the vacuum invariant.
Given assumption ii), we conclude that either no scalar couples to the down-type quarks (which would lead to massless quarks) or else all scalars couple to down-type quarks. This conclusion is a particular case of the following more general assertion: Proposition. For the down-type Yukawa terms in eq. (2.1), each set of Higgs doublets φ k , comprising an irreducible representation of G, either couples to quark fields with linearly independent Γ k , or decouples completely with Γ k = 0.
The same proposition holds for the up-type Yukawa terms and their coefficients ∆ k . The proof is given in Appendix A; its essential idea is that if some Γ k were linearly dependent, then one would be able to identify doublets which do not couple to the down-quark fields and represent an invariant subspace under the action of G, this making the representation reducible. Thus, no matrix Γ k vanishes and for some given group element g 1 . The set of all elements g 1 satisfying eq. (2.12) forms a subgroup of G, which we denote by G q . The Lagrangian after SSB is thus left invariant by a residual symmetry G q . Then, specializing in eq. (2.11) to an element g = g 1 in G q , and using eqs. (2.3) and eq. (2.12), one gets Under these assumptions, M d is symmetric under the group element g 1 and the subgroup generated by it.
If assumption ii) remains valid but we substitute iii) by the assumption that the element g ∈ U (1) B is trivial (in the newly defined sense), then we conclude that M d has no symmetry, i.e., the group G has been completely broken by the vacuum. In this case, we call the residual symmetry trivial, i.e. G q ⊆ U (1) B .
We now turn to the possibility that Φ is in a reducible representation of G, thus negating our assumption ii). First we assume that Φ breaks into the three vectors (2.14) which transform as the irreducible representations ϕ 1 , ϕ 2 , and ϕ 3 of G, respectively. We can still write eq. (2.11). The only novelty is that G g is now a block-diagonal matrix, with one r × r block, one (s − r) × (s − r) block, and one (N − s) × (N − s) block. Said otherwise, we may write where we introduce a label ϕ j denoting the representation. We assume that the ϕ 3 representation is such that the fields (φ s+1 , · · · , φ N ) do not couple to the down-type quarks, i.e., Γ k ϕ 3 = 0. In this case, We further assume that the vev of (φ 1 , · · · , φ r ) is invariant under some subgroup G 1 , while the vev of (φ r+1 , · · · , φ s ) is invariant under some subgroup G 2 . Thus, for g 1 ∈ G 1 and g 2 ∈ G 2 . The common residual symmetry is now G q = G 1 ∩ G 2 . If G q ⊆ U (1) B , then the residual symmetry is trivial. In contrast, if we have a non-trivial G q ⊆ U (1) B , we can pick g ∈ G q and using eqs. (2.15) and (2.16), together with eqs. (2.18) and (2.19), we obtain again eq. (2.13). Notice that, in this example, ϕ 3 is irrelevant for M d but it may matter for M u , should (φ s+1 , · · · , φ N ) couple to the up-type quarks.
We are now ready to generalize the impact of symmetries and to include the up-quark sector. The residual symmetry G q restricted to the down-quark sector is now renamed as G d , while G q is reserved for the residual symmetry of the quark sector as a whole. If Φ is in a reducible representation of G, it may be seen as a sum of irreducible representations Equation (2.11) can now be written as a set of equations for each ϕ j , where no sum is assumed in ϕ j . The matrices Γ k ϕ j and ∆ k ϕ j denote the couplings associated with φ k for a given irreducible representation ϕ j . For example, if Φ is an irreducible quadruplet (φ 1 , φ 2 , φ 3 , φ 4 ) then Γ k , for k = 1 . . . 4, denotes the Yukawa coupling associated with each entry of the scalar representation. If Φ is a reducible quadruplet, composed by two irreducible doublets ϕ 1 = (φ 1 , φ 2 ) and ϕ 2 = (φ 3 , φ 4 ), then the coupling Γ k ϕ i associates two couplings to each irreducible representation, i.e. k = 1, 2. One denotes the vev of the irreducible representations ϕ j as v ϕ j . When the vev breaks the group G, the Lagrangian may still remain invariant under the action of some elements of G.
Let us denote by G d and G u the subgroups of G that are left invariant by the vevs of all Higgs fields coupling to down-type and up-type quarks, respectively. Some fields may interact with both sectors. The relevant residual symmetry of the quark sector is G q = G d ∩ G u . The fields that do not couple to quarks are irrelevant to our discussion. There are several possibilities: there is no non-trivial symmetry left in the quark sector as a whole.
More specifically, , then there is no non-trivial symmetry left in either (up or down) quark sector; (2) If G q ⊆ U (1) B , then there is some non-trivial symmetry left in the quark sector as a whole.
The characteristic in common to the cases in (1) is that Φ breaks the symmetry completely, modulo U (1) B . In contrast, case (2) leaves a common non-trivial symmetry in the whole quark sector. In the latter case, for g ∈ G q , from eqs. (2.21) we find 22) generalizes eq. (2.13) to the case of reducible representations and is the main result of this section. Thus, when using eq. (2.22), as we will do in the theorem to be proved in the next section, it does not matter whether the representation of the group G is reducible or irreducible.

Theorem and proof
Let us consider a model with quarks and scalar fields transforming under some set of representations of a full flavour group G. We denote the representation space where G acts as the flavour space, i.e., the horizontal space of replicated multiplets of quark and scalar fields with the same gauge quantum numbers. The following theorem strongly constrains the viable models.
Theorem (No-Go). Given a group G acting on the flavour space, the only way to obtain a non-block-diagonal CKM mixing matrix and, simultaneously, non-degenerate and non-zero quark masses, is that Φ breaks completely the group G, except possibly for some symmetry belonging to baryon number.
Proof. Suppose that there is a residual symmetry group G q which is left invariant by the vev Φ . Let us denote a generic element of G q as g 1 , and then eq. (2.12) is assumed. 1 We want to show that either g 1 = e, i.e., the trivial element in G, or g 1 acts on the flavour space as a member of the baryon number symmetry.
Equation (2.12) implies that the Lagrangian after SSB is invariant under g 1 . Using the same rationale that led to eq. (2.22), we conclude that Since we shall be interested in the left sector, where the CKM mixing arises, we may turn to the Hermitian combinations in eq. (2.5), which verify The next step is to use the following result (proved in Appendix B): if the generator G Lg 1 given in eq. (2.24) is not proportional to the identity matrix, then the CKM matrix is block-diagonal, in contradiction with experiment. Note that we include any permutation of a block-diagonal matrix, as well as the unit mixing matrix, V = ½, in the category of block-diagonal matrices. We then conclude that Returning now to eqs. (2.23), it follows that because M u and M d are invertible. Applying this to eq. (2.11), we find for l = 1, . . . , N . This implies that  .29) jointly. If Φ is in a faithful representation, then g 1 = e and no residual symmetry is present. In this case, e iθ = 1 in eq. (2.25). If g 1 = e (e iθ = 1), then Φ is unfaithful and some residual symmetry will be present in the final Lagrangian, without constraining the mixing. This residual symmetry should be specifically represented by eqs. (2.25), (2.26) and (2.29), which is just part of the baryon number conservation. This completes the proof of the theorem. 2 At this point two remarks are in order. First, any SM-like Lagrangian, as the one in eq. (2.1), exhibits automatic conservation of the baryon number U (1) B , independently of additional gauge or flavour symmetries. Such accidental symmetry imposes no constraints on masses and mixing, and remains conserved after electroweak symmetry breaking. The proof above shows that a non-block-diagonal V CKM is compatible with some residual symmetry inside G only if the latter is a subgroup belonging to U (1) B . Second, we should stress that eq. (2.29) applies only to the Higgs doublets that couple to quarks. The Higgs doublets that appear solely in the scalar potential are important when analyzing possible vacuum alignments, but they are irrelevant for the statement of the theorem.
Notice that no information of the scalar potential and vevs has been used in this proof. That is, the theorem constrains residual symmetries and can be applied to virtually any model. Nevertheless, our interest in this article lies on NHDM, and, as we will show in section 3, the explicit transformation properties of the Higgs scalars, their vevs, and their relation to the subgroup chains are crucial, in particular, for the application of the theorem to specific symmetries of complete models of scalars and fermions.

Some illustrative examples
The theorem proved in the previous section is completely general; yet it is useful to study in more detail several particular cases that one may come across. We split them into four classes depending on whether the representations D(Q L ) and D(Φ) of G, acting respectively on the left quark space Q L and the scalar space Φ, are faithful or not. We denote the kernel of the representation of D The full symmetry in the left sector is smaller than the one in the scalar sector and, in principle, we could have the proper subgroup ker D(Q L ) of G (or smaller) unbroken by Φ . However, eq. (2.29) should hold and Φ should break G completely.
As an example, let us take This case is automatically excluded unless ker D(Φ) acts like baryon number on quarks. The full symmetry of the potential is G/ ker D(Φ) and Φ can never break ker D(Φ).
As an example, let us take G = A 4 with (φ 1 , φ 2 , φ 3 ) ∼ (1, 1 ′ , 1 ′′ ) and Q L ∼ d R ∼ u R ∼ 3. If all φ k get nonzero vevs, they only break A 4 to Z 2 × Z 2 , which corresponds to the kernel of the representation of Φ and is generated by diag (1, −1, −1) and diag (−1, −1, 1). This subgroup remains in the Q L quark sector and the CKM matrix would be trivial. This case can be discarded given the Yukawa structure in eq. (2.1) and the assumption that at least one of the representations for Q L , Φ, d R , u R is a faithful irrep or contains a faithful irrep (i.e. G, and not a smaller group, is the full flavour symmetry); see also Appendix C.
During the proof and consequent remarks of the theorem, nothing has been said about the reducibility or irreducibility of Φ. If Φ is an irreducible representation of the full flavour group G, then the same field components φ i couple to the up and down sectors. In the case that Φ is a reducible representation of G, the components φ i can be arranged into irreducible multiplets just like in eq. (2.20). In this case, there are many ways of breaking the flavour group G completely in order to allow for models with non-zero and non-degenerate masses with viable V CKM (see end of section 2.2 for the different possibilities).
We shall comment on case (1d) of section 2.2, which is the most predictive and hence one of the most common approaches in model building, since each quark sector has its own non-trivial residual symmetry while the full flavour group G is completely broken. This will lead, in general, to strong constraints to the mixing angles and the CP phase. For instance, in ref. [14], the authors studied which would be the best residual symmetries in order to accommodate the experimental CKM mixing matrix. This was done in the context of the von Dyck flavour groups [13]. It is shown that the flavour groups D N and ∆(6N 2 ), with N an integer multiple of 7, are capable of determining the mixing angle θ 12 close to its experimental value, if we fix θ 23 = θ 13 = 0. Non-zero values for all three mixing angles require, on the other hand, infinite von Dyck groups.
We summarize in Table 1 the different possibilities for the CKM mixing angles and mass spectra depending on the existence (or not) of a residual symmetry in the whole theory. In the presence of a residual symmetry, not contained in U (1) B , we can either have at most three non-zero mixing angles with some massless quarks, as in example (ii) above, or have massive quarks with at most one non-zero mixing angle, as in example (iii).

Residual symmetry
Non-zero mixing angles CP violation Quark spectrum As a final remark, one should notice that until now we have been working within the framework of NHDM. However, most flavour models in the literature keep the number of SU (2) Higgs doublets small, and add additional SU (2) scalar singlets as a way to break flavour. The dynamics of these fields is not taken into account for low energy processes, being solely responsible for the origin of the flavour breaking. These are known as flavon fields. All the analysis done for NHDM can be extended straightforwardly to flavon models. For example, if there is only one Higgs doublet φ carrying no flavour charges, the breaking of flavour symmetries can be performed by flavon fields η = (η 1 , η 2 , ..., η N ) through the field combinations Φ = φ η/Λ, which behave just as flavoured Higgs doublets.

Application of the theorem to flavour models
The no-go theorem previously demonstrated strongly constrains viable flavour models based on a full flavour symmetry G of the Lagrangian. In this section, we shall apply the theorem to some examples of flavour groups commonly found in the literature to explain fermion masses and mixing.
Before proceeding with the examples, it is worth recalling the distinction between what is meant by "full flavour group" and by simply "flavour group". As said before, in model building, we may add symmetries that act only on the flavour space of the model. Those symmetries are usually taken to be subgroups of SU (N ) (in particular SU (3), when we are dealing with three fermion generations). We refer to these symmetries as flavour symmetries K. However, the full flavour symmetry group G is, in general, larger since it contains the global hypercharge transformation, i.e G = K × U (1) Y . 3 The accidental baryon symmetry U (1) B is not included in the full flavour symmetry G because the Higgs fields do not transform under it, but it might happen that G intersects U (1) B . We further assume that there is no additional accidental symmetry in the Yukawa Lagrangian (2.1), apart from G and U (1) B (see refs. [15] for methods for detecting them).
The additional U (1) Y in the full flavour symmetry is of extreme importance. As we shall see, one may break completely the flavour symmetry K but still leave a non-trivial residual symmetry in G. To better understand this scenario, we introduce the notion of pseudo-invariant vacuum. In our discussion we have used the invariance of the vacuum, i.e. eq. (2.12), in order to constrain the possible viable models. This vacuum invariance has been defined in the full flavour group G. Looking only to the flavour group K, we may extend the notion of invariance to pseudo-invariance of an irrep through the condition (3.1) The group elementg 1 changes the vacuum by a global phase transformation. In general, these global phase transformations are not contained in K, but they are contained in G due to the U (1) Y global group. This means that, for a given irrep, any pseudo-invariant vacuum in K can be written as an invariant one in G. Therefore, instead of dealing with the full flavour group G, we may restrict ourselves to the flavour group K and its nontrivial pseudo-invariant vevs. Note that if Φ is composed of more than one irrep ϕ j , as in eq. (2.20), a pseudo-invariant vev where is equivalent to an invariant vev of G only if all vevs in each irrep transform by the same global phase, i.e., α j = α. One consequence immediately follows: the presence of an invariant Higgs doublet, singlet with respect to K, may be relevant. For example, consider down-type quarks interacting only with one n-dimensional irrep ϕ 1 = (φ 1 , · · · , φ n ) whose vev is pseudo-invariant by an elementg 1 = e with α 1 = 0. The mass matrix M d will possess a non-trivial residual symmetry. Now, suppose we add a trivial irrep ϕ 2 = φ n+1 ∼ 1 that also interacts with down-type quarks. Any vev for ϕ 2 will be strictly invariant by K, with α 2 = 0 in eq. (3.2). Hence, M d will no longer possess the previous residual symmetry.
We remark that eq. (3.2) (and consequently eq. (3.1)) is merely an eigenvalue equation. Therefore, from a simple group-theoretical method, and without analyzing the scalar potential, we can extract a set of vacuum alignments that will be automatically excluded by the theorem, namely those corresponding to the eigenvectors of Gg 1 . Notice however that the vev alignments obtained through this procedure may not be global minima of the scalar potential. Yet, this is a straightforward group-theoretical check in the spirit familiar to model building. One could instead follow a geometrical method [8] to find first the global minima of the scalar potential. In the latter case, if all the minima preserve some subgroup of the initial symmetry, the theorem applies directly without the need of solving the eigenvalue equation (3.2). The drawback is that the geometrical approach is not universal, and it is not guaranteed that it could be easily applied to more complicated Higgs sectors, for example to some high symmetry groups in 4HDM.
Next we shall present three examples of flavour groups in the context of 3HDM, two of them with Φ being in a faithful triplet representation, for the ∆(27) and A 4 flavour groups, and one example with Φ in a reducible triplet representation, for the S 3 flavour group.

The flavour group ∆(27) in 3HDM
The flavour group ∆(27) [16] is a subgroup of SU (3) and can be viewed as the group defined by two generators a and b, with presentation We use the specific three-dimensional implementation [7]: where ω = e 2iπ/3 . The group has 27 elements divided into 11 irreducible representations: nine singlets 1 (i,j) , with i, j = 0, 1, 2; and two triplets, 3 (0,1) and 3 (0,2) . There are two elements of particular interest. One is Indeed, since ∆(27) is isomorphic to Z 3 ×Z 3 ⋊Z 3 we can view the first two Z 3 's as generated by a and d (they obviously commute), while the third one is generated by b. The second element of interest is The  Any of the vacuum alignments presented above automatically leads to a non-trivial residual symmetry in the final Lagrangian. Actually, because all elements of the triplet irrep of ∆(27) have eigenvalues 1, ω or ω 2 , any pseudo-invariant vacuum alignment can be seen, with the help of h, as an invariant one of another element of ∆(27). Notice also that h leaves pseudo-invariant any type of vacuum alignment. However, it is not of real interest since h ∈ Z 3 , i.e. is the center of SU (3). As already pointed out, these transformations are included in the U (1) Y part of full flavour group G. Therefore, we just need to look at the ∆(27)/Z 3 group. In this case, some alignments in eq. (3.7) will break the flavour group, but the full flavour group will still be broken to a non-trivial subgroup.
The group theoretical method presented above allow us to exclude several alignments without the need of dealing with the details of the scalar potential. However, this method tell us nothing about the global minima of the potential. Minimizing the general scalar potential for 3HDM is a very hard task, and until now no good method has been developed. When the potential exhibits some symmetry, the geometrical method [8] can give us the answer (specially for large symmetry groups). This method has been used in the study of ∆(27) in 3HDM, where four classes of global minima are found [18], corresponding precisely to the eigenvectors given in eq.

The flavour group A 4 in 3HDM
The flavour group A 4 can be viewed as the group defined by two generators s and t, with presentation The group A 4 has one non-trivial invariant subgroup; Z 2 × Z 2 , generated by s and tst 2 . It also has four subgroups Z 3 , generated by t, tst, st and st 2 , and three subgroups Z 2 generated by s, t 2 st and tst 2 . A summary of the subgroups and their generators is shown in figure 2, based on ref. [17]. A possible matrix representation for the tridimensional irrep is We can use the group theoretical method in order to find the set of alignments forbidden by the theorem. In this case, the eigenvectors (v 1 , v 2 , v 3 ) associated to the generators of the subgroups are where v 1 , v 2 and v 3 are complex parameters. Using the geometrical method, the global minima of the A 4 symmetric scalar potential were found to be [8] (1, 1, 1) , (1, 0, 0) , (±1, e iπ/3 , e −iπ/3 ) , (1, e iα , 0) , (3.11) up to permutations and global phase transformations. All these global minima are included in the list of alignments given in eq. (3.10). Indeed, the first two are obviously there; the third vacuum alignment is (±1, e iπ/3 , e −iπ/3 ) = −(∓1, ω 2 , ω) and is in the list as well. Finally, the last minimum is a particular case when |v 1 | = |v 2 |. In summary, each A 4 global minimum leaves invariant a non-trivial subgroup of the full flavour symmetry G = A 4 × U (1) Y and is thus excluded by the no-go theorem. That is, one cannot construct a viable model based on A 4 where Φ is in a faithful (triplet) irrep [9]. It is interesting to consider an example similar to the one we included in point (ii) of section 2.4, where the representations are This implies that, after up and down quark mass matrix diagonalization, there remain three mixing angles and three non-degenerate quark masses (one of them massless in each sector). This is allowed if we had looked exclusively at the left-handed quark sector, as in ref. [11]. However, it is forbidden once we consider the right-handed quark sector.

The flavour group S 3 in 3HDM
Since the group S 3 is the smallest non-Abelian discrete symmetry, it has been widely used for flavour physics in the quark and lepton sectors (for recent works see e.g. refs. [19,20], and references therein). The flavour group S 3 can be viewed as the group defined by two generators, a and b, with presentation It has the non-trivial invariant subgroup Z 3 , which is generated by a, and three subgroups Z 2 , generated by b, ab and a 2 b. We summarize in figure 3 these subgroups and their generators [17]. The S 3 flavour group has no 3-dimensional irrep. Therefore, any 3HDM with a scalar potential invariant under this group, with a faithful representation, must be built from a 2-dimensional irrep 2 and one of the 1-dimensional irreps, 1 or 1 ′ . The representation of the generators for these irreps can be chosen as  When building an S 3 flavour model we have to decide whether our reducible tridimensional representation is 2 + 1 or 2 + 1 ′ . We may represent them in the compact form with the plus and minus signs assigned to 2 + 1 and 2 + 1 ′ , respectively. Following the group theoretical method, we can find the set of vev alignments forbidden by the theorem. In this case, the eigenvectors associated to the reducible generators of the subgroups will be denoted by ((v 1 , v 2

Conclusions
In this work we have studied the connection between the breaking of flavour symmetries by the Higgs vevs and the existence of residual symmetries in the quark mass matrices. A simple but powerful no-go theorem has been formulated: the vevs of the scalars coupled to the quarks have to break the full flavour symmetry (or break it into a subgroup acting as baryon number) in order to avoid nonphysical quark masses and mixing. In this context, the notion of full flavour group turns out to be important. While, in general, one refers to the flavour group as the group added to the SM acting globally on the flavour space, the SM gauge group already contains the global hypercharge transformation that should also be taken into account.  27) symmetric potentials, this is indeed verified and thus these models are excluded. The same type of analysis can be carried out for other groups and, in particular, for smaller groups which have a more complex scalar potential. While we only treated quarks in our discussion, the extension of the theorem to the whole fermion sector, including leptons, is straightforward. If neutrinos are Dirac-type particles the analogy is direct. There is the experimental possibility of a massless neutrino, however, as shown in our analysis, all cases with CP violation always imply a massless fermion in each sector. The case where neutrinos are Majorana-type particles is more interesting since there is a larger number of possible implementations. Lepton number is no longer conserved and thus the theorem will be slightly modified. either all of them decouple from all up quarks or all down quarks, or none of them can decouple from them. This statement is stronger than just saying that none of the matrices Γ k or ∆ k can be zero in a given Higgs basis. It says that this cannot happen in any Higgs basis. The basis-invariant formulation of this requirement is that Γ k or ∆ k are either all zeros or are linearly independent.
In this appendix, we provide an accurate proof of this statement. We start with the following proposition.
Proposition. Let {Γ k }, k = 1, . . . , N , be a non-zero element of a complex vector space Suppose that a symmetry group G acts in the space of {φ k }, and {φ k } realize an irrep of G. Finally, suppose that there exists an expression L, proportional to Γ k φ k (summation over repeated indices assumed), which is invariant under G. Then, Γ k are linearly independent.
Before proceeding to the proof, a few remarks are in order. First, to make the notation a bit more familiar, we notice that in our case V Γ = C 9 of complex-valued 3 × 3 matrices Γ k , and V φ = C 2 , the space of Higgs doublets. However, the proposition itself is not specific to this particular space. Also, we note that the linear independence among Γ k is a basisinvariant way of saying that no individual Γ k will ever become zero after any basis change in φ k . Finally, {Γ k } being a non-zero element of N V Γ means that not all individual Γ k = 0.
Proof. We shall prove the statement by contradiction: assuming that Γ k are linearly dependent, we will show that none {φ k } can be an irrep.
Suppose that Γ k are linearly dependent, namely, that there exist complex coefficients c k , not all of them being zeros, such that c 1 Γ 1 + · · · + c N Γ N = 0. Denote the number of linearly independent Γ k by p, with 0 < p < N . Then, it is possible to perform a correlated basis transformation in the spaces of Γ's and φ's leaving invariant Γ k φ k , which sets N − p matrices Γ's to zero: for k = 1, . . . , p : Γ k = 0 and they are linearly independent, for k = p + 1, . . . , N : In this basis, we call the corresponding Higgs doublets φ k active for k = 1, . . . , p and passive for k = p + 1, . . . , N . The expression L is thus written only in terms of active fields. Let us now apply a transformation g ∈ G, which leads to The last transformation here is just a relabeling of the indices. Note that we do not require that Γ k φ k itself is invariant under g, because L can contain other overall factors which are also transformed under g (in our particular case, these are left-and right-handed quark fields). These additional factors lead to the appearance of extra transformations (such as G L and G d R ), which are unitary transformations of the space V Γ and transform each individual Γ k , but do not mix different Γ's. In particular, they cannot make a zero Γ k non-zero, and vice versa, nor can they convert a linearly independent set of Γ's into a linearly dependent one.
By construction, the superscript l in the last expression of eq. (A.2) refers to active fields. Since Γ l (G g ) lk is a linear combination of linearly independent Γ's, if k also corresponds to an active field, then this linear combination can stay non-zero, and (G g ) lk can be anything. If k corresponds to a passive field, then this linear combination is zero, which can happen only when (G g ) lk = 0 for all active l's. Since G g is unitary, it also follows that (G g ) lk = 0 for all passive l and active k. In other words, the group transformation G g does not mix active and passive fields, and hence it has a block diagonal form in the space of φ's. Repeating this analysis for all g ∈ G, with the same conclusion, we obtain that active and passive fields represent two invariant subspaces of the space {φ k }. Thus, the representation is reducible. This completes the proof.
The above proposition can be immediately applied to the one given in section 2.2. Suppose that we pick up Higgs doublets which realize an irrep of the flavor group G; this can be either the entire set of doublets of the model or, in the case when the Higgs doublets are in a reducible representation, this can be an irreducible subspace of the Higgs fields. In any of these cases, we apply the above proposition, assuming that N stands for the dimension of the chosen irreducible representation. With this minor modification, we obtain exactly the proposition from the main text.

B Block-diagonal CKM matrix
We prove here the assertion made after eq. (2.24), that if the generator G Lg 1 given in eq. (2.24) is not proportional to the identity matrix, then the CKM matrix is block-diagonal.
The matrix G Lg 1 , being a 3 × 3 unitary matrix, is diagonalizable and we can find three orthonormal eigenvectors. If G Lg 1 is not proportional to the identity, then there is at least one non-degenerate eigenvalue e iθ . Take its associated eigenvector w, for which G Lg 1 w = e iθ w. It follows from eq. (2.24) that i.e., (H d w) = 0 is also an eigenvector of G Lg 1 with eigenvalue e iθ . Since this eigenvalue is non-degenerate, (H d w) should be proportional to w: H d w = m 2 d i w, where obviously m d i is one of the down-quark masses. Hence, w fills one of the columns of the matrix V d which diagonalizes H d . Similar arguments make the same vector w to be in one of the columns of the matrix V u , which diagonalizes H u . The CKM matrix V = V † u V d then contains one row and one column consisting of zeros, except for a unit entry where they cross. Finally, the vector w will remain as a common eigenvector of H d and H u even if H d w = 0 or H u w = 0 (some massless quark). This proves the assertion.

C Yukawa structure
We show here that some combination of representations cannot be assigned to the fields Q L , Φ, d R and u R in order to construct group invariants in the Yukawa interactions of eq. (2.1). For instance, it is not possible that only one ofQ L , Φ, d R or u R is assigned to a faithful irrep or contains a faithful irrep.
We begin by reviewing the following result: given two vectors (fields) x and y that transform under irreducible representations µ and ν of G, the unique invariant present in x × y is x · y = x i y i , only possible when µ = ν * . The proof relies on the Schur lemmas. The invariant in x × y can be written as I = x i y j C ij where C ij are the Clebsch-Gordan coefficients that connect x × y to the trivial representation. Invariance of I requires (D (µ) (g)) T CD (ν) (g) = C , (C.1) for all g in G. We can rewrite it as CD (ν) (g) = (D (µ) (g)) * C. One of the Schur lemmas requires that C = 0 only if µ ≃ ν * . However, if µ = ν * , another Schur lemma ensures that C = ½. Now, this result also applies if x and y are reducible and composed of more than one irreducible piece: complex conjugate irreps need to be contracted to build invariants. Let us now analyzeQ L × Φ × d R . Suppose that Q L and Φ are assigned to unfaithful representations and d R to a faithful irrep. Then the representation ofQ L × Φ is also unfaithful and does not contain any faithful representation that can be contracted to d R .