Family Symmetries and Multi Higgs Doublet Models

Imposing a family symmetry on the Standard Model in order to reduce the number of its free parameters, due to the Schur’s Lemma, requires an explicit breaking of this symmetry. To avoid the need for this symmetry to break, additional Higgs doublets can be introduced. In such an extension of the Standard Model, we investigate family symmetries of the Yukawa Lagrangian. We find that adding a second Higgs doublet (2HDM) does not help, at least for finite subgroups of the U(3) group up to the order of 1025.


Introduction
At currently achievable energies, the Standard Model (SM) of fundamental interactions is a very good working theory. However, it is commonly believed that it is only an effective theory which at higher energies need to be modified. One of the signs of this state of things is a large number of free parameters (more than 20) which now need to be fitted from experiments. The main parameters are: masses, mixing angles and CP violating phases for quarks and leptons. The SM does not explain those parameters but introduces the mechanism by means of which all particles acquire masses by the so called Higgs mechanism. One of several proposals how to restrict number of a free parameters in the SM is to introduce symmetry between Yukawa constants in Yukawa SM interaction in such a way that after spontaneous symmetry breaking get masses and mixing matrix parameters for quarks and leptons which are consistent with experience. Such symmetry is known in the literature as a flavor symmetry [1] (but also family or horizontal symmetry). Lepton sector and especially neutrino physics is an attractive area to search for such a symmetry due to the so called lepton mixing matrix [2,3].
There exist direct links between the mixing and lepton masses. Charged lepton and neutrino mass matrices M l(ν) are diagonalized by biunitary transformations (for Majorana neutrinos (U ν ) R = (U ν ) * L ): The lepton mixing matrix U PMNS is composed from the charged lepton (U l ) L and neutrino (U ν ) L diagonalizing matrices: where: are 3 dimensional representations matrices for the left handed lepton doublets for some N-order flavour symmetry group G. As a direct consequence of the Schur's Lemma-since M l M † l and M ν M † ν are proportional to the identity matrix, the lepton mixing matrix U PMNS becomes trivial.
In the literature there are some ideas about how to escape from the trivialisation of a matrix (2). One approach is to break the family symmetry group by scalar singlet-so called "flavons" (e.g., [6,7]). Non trivial mixing can be also achieved by extending the Higgs sector by additional multiplets (e.g., [8,9]). A proposal for the two Higgs doublet model (2HDM) [10] was widely discussed in [11]. In the most general situation, the mass generation mechanism allows couplings with various Higgs fields. In this context, extensions of the SM assuming the existence of different numbers of doublets and Higgs triplets are allowed. Theoretical proposals assuming the existence of two Higgs fields are not the only possible and potentially experimentally verifiable space for applying the symmetry implementation. In this paper, the methodology proposed previously for 2HDM only [11] is extended to any number of additional Higgs fields.
Additionally, new forms of results, equivalent to [11], for 2HDM are given. The obtained results depend on many phases and we present here a more detailed discussion concerning relations between them. We hope that it may help to determine the analytical origin of these solutions.

Multi Higgs Doublet Description
Discussion below stands for Dirac neutrinos. To describe the coupling between lepton fields and the Higgs field we take the n-Higgs doublet Yukawa interaction term of the form: where i = 1, 2, . . . , n and α, β = e, µ, τ. The charged lepton states l βR and neutrinos ν βR are right-handed SU(2) singlets and then the gauge doublets for the left-handed lepton and Higgs fields are: where ϕ 0 i and ϕ − i are complex scalar fields in spacetime for i = 1, 2, . . . , n and σ 2 = The 3 × 3 Yukawa matrices h l i and h ν i each define the couplings of left-handed doublets with right-handed singlets via the i-th Higgs doublet. Due to the form of the Higgs potentials, ground states occur at non-zero ϕ, with the vacuum expectation values: for some complex-valued v i , where: Mass matrices for charged leptons and neutrinos read as follows: For some finite flavour group G, the family symmetry means that after fields transformations (A L , A R l and A R ν are 3 dimensional representations): and (A Φ is a n dimensional representation): the Lagrangian does not change: Symmetry conditions can be written as an eigenproblem of a direct product of unitary group representations to the eigenvalue 1. For any group elements we have: It is sufficient to check the above equations for group generators only as then they will automatically be satisfied for all group elements.
In such a model, the invariance equations for the mass matrices are not trivial, so we avoid the consequences of Schur's Lemma: The same conclusion is valid if one assumes that neutrinos have Majorana nature. In such a frame, the Yukawa interaction Lagrangian has to be rewritten in an appropriate way (see [11]) producing family symmetry condition in the form:

Two Higgs Doublet Model (2HDM) Results
As a potential flavour symmetry group G, we chose finite, non-abelian subgroups of U(3), up to the order of 1025. This class of groups is very important in practice [12], even though there exist models in which the flavour symmetry group cannot be conceived as a subgroup of U(3) (the upper limit on the group order was of course due to the calculation time). Using the GAP [13] system for computational discrete algebra, with the included SMALL GROUPS LIBRARY [14] and the REPSN [15] package for constructing representations of finite groups, we have found groups which fulfil the requirements of our model and impose flavour symmetry on the Yukawa Lagrangian. Next we have calculated the Yukawa matrices and created mass matrices and mixing matrices. The last step was to check agreement with experimental data. In total we have found 10862 groups with at least one 2 dimensional and at least one 3 dimensional irreducible representation. Only 413 of these groups are subgroups of U(3). Either a group has at least one faithful 3 dimensional irreducible representation (there are 173 such groups) or it has at least one faithful 1+2 reducible representation (there are 240 such groups).

Results for Dirac Combinations
All obtained solutions for Yukawa matrices for Dirac neutrinos and charged leptons can be expressed through seven base forms. Putting ω = e 2π i k / 3 and allowing any integer k value (note that ω 3 = 1 and ω 2 = ω * ), the first three forms are: The next three forms can be obtained from the above ones through a simple interchange h 1 h 2 . The last, seventh form, valid only for k = . . . , −2, 1, 4, . . . , is (note the diagonality of these matrices): Dirac neutrinos are always defined by "ordered" pairs: where φ are some real phases (which depend on the actual group and its representations' combinations).
For each of such Dirac neutrinos' solutions, there always exist two different solutions for charged leptons, defined by the two corresponding "ordered" pairs: where δ l = 0, π. Assuming complex c x , real v 1 , φ 1 , v 2 , φ 2 and putting: we get exactly the same set of three eigenvalues of M M † for any of the seven above Yukawa matrices forms (note that this set of eigenvalues is invariant with respect to v 1 v 2 and/or φ 1 φ 2 , hence we can safely use them for all cases h 1 h 2 , and that if v 1 v 2 ≥ 0 then m 2 The mass squared differences are (note that if v 1 v 2 ≥ 0 then ∆m 2 31 ≥ ∆m 2 32 ≥ ∆m 2 21 ≥ 0 when 7π/6 ≤ φ 2 − φ 1 ≤ 4π/3 and ∆m 2 31 ≤ ∆m 2 32 ≤ ∆m 2 21 ≤ 0 when π/6 ≤ φ 2 − φ 1 ≤ π/3): (so, no neutrino mixing is possible at all), while for the first six forms it is:

In all cases
and the unitary matrix that diagonalizes it (so that U † M M † U = diag m 2 1 , m 2 2 , m 2 3 ) is: Note that the U matrix does not depend on the phase difference φ 2 − φ 1 at all so, it will be exactly the same for neutrinos (U ν ) L and charged leptons (U l ) L . Hence, the U PMNS = (U l ) † L (U ν ) L = I (so again, no neutrino mixing is possible at all).
In order to directly apply the above equations to Dirac neutrinos, one should put where we assume that the vacuum expectation values are v 1 e i φ 1 and v 2 e i φ 2 (the same for neutrinos and charged leptons, of course). When moving between Dirac neutrinos and charged leptons, we can easily notice that all equations are invariant with respect to v 1 v 2 and the phase difference φ 2 − φ 1 is simply shifted by δ l = 0, π. That means that, for δ l = π, all mass squared differences ∆m 2 ij will change signs, so that their mass ordering schemes will be reversed, while for δ l = 0 there will be no change at all.
The ratios of experimental mass squared differences for neutrinos and charged leptons are: and so they cannot be reproduced in this theoretical frame (they would need to be exactly equal while they differ by a factor of about 8 to 9).
Moreover, Equation (27) are not able to reproduce experimental masses of charged leptons at all. They return, in all cases, 1 ≤ m τ /m µ ≤ 2 and requiring m µ /m e ≈ 206.8 or m τ /m e ≈ 3477 returns m τ ≈ m µ .

Results for Majorana Combinations
All obtained solutions for Majorana neutrinos, which do not produce a scalar mass matrix M ν , are defined by "ordered" quadruples: where h 1 and h 2 are defined by Equation (23), φ are some real phases (which depend on the actual group and its representations' combinations), δ = 0, π and v 0 e i φ 0 is a free complex parameter. For each of such Majorana neutrinos' solutions, there always exist two different solutions for charged leptons, defined by the two corresponding "ordered" pairs { h l 1 , h l 2 } = { h 2 , e −i (δ l +φ) h 1 }, where δ l = 0, π (note that δ l is independent of δ so, there always exist four different combinations for each φ).
Assuming complex c ν (usually denoted as g/(2M)), real v 0 , φ 0 , v 1 , φ 1 , v 2 , φ 2 and putting (where we assume that the vacuum expectation values are v 1 e i φ 1 and v 2 e i φ 2 ): we find that the neutrino mass matrix M ν is a diagonal matrix and The mass squared differences are: For δ = π (e i δ = −1), we can also reuse Equations (27) and (28), when we put c ). Note that the mass ordering schemes will be reversed between δ = 0 and δ = π, which can easily be seen if one puts v 0 = 0 in Equations (34) and (35). For the corresponding charged leptons, we can also simply reuse Equations (27) and (28), where we put . As already mentioned, Equation (27) is not able to reproduce experimental masses of charged leptons at all. Moreover, as all relevant matrices are diagonal (for both neutrinos and charged leptons), no neutrino mixing is possible at all.

Conclusions
Multiple Higgs doublet models are in general promising that, in order to get a non trivial lepton mixing matrix, one will not need to explicitly break the family symmetry. However, our results for the 2HDM are utterly negative. The big open question is why. First, we select finite, non-abelian subgroups of U(3) which are provided by the SMALL GROUPS LIBRARY [14] in GAP [13]. Is it possible that we miss some vital groups (due to the constrains that we impose when selecting groups or due to the used library itself)? Then, is it possible that non trivial solutions can only be obtained for models with more than two Higgs doublets? Finally, is it possible that the Equations (18), (19) and (21) will always lead to only trivial solutions, even though these models seem to be free from the consequences of the Schur's Lemma (20)? In order to address, at least partially, the last two questions, we are currently working on processing groups which are not subgroups of U(3) for 2HDM and on a family symmetry approach with three Higgs doublets, hoping that it will give some positive outcome.