Exploring the quark flavor puzzle within the three-Higgs double model

# Exploring the quark flavor puzzle within the three-Higgs double model

David Emmanuel-Costa    J.I. Silva-Marcos Centro de Física Teórica de Partículas, CFTP, and Departamento de Física, Instituto Superior Técnico, Universidade de Lisboa, Avenida Rovisco Pais nr. 1, 1049-001 Lisboa, Portugal    Nuno Rosa Agostinho Departament d’Estructura i Constituents de la Matèria and Institut de Ciencies del Cosmos, Universitat de Barcelona, Diagonal 647, E-08028 Barcelona, Spain
###### Abstract

We extend the standard model with two extra Higgs doublets. Making use of a symmetry principle, we present flavour symmetries based on cycle groups that oblige each Higgs doublet to contribute to the mass of only one generation. The Higgs doublets couple to the fermions with different strengths and in this way accommodate the quark mass hierarchy. We systematically search for all charge configurations that naturally lead to the alignment in flavour space of the quark sectors, resulting in a CKM matrix near to the identity, determined by the quark mass hierarchy, and with the correct overall phenomenological features. The minimal realization is by the group . We show that only a limited number of solutions exist, and that any accidental global symmetry that may occur together with the discrete symmetry is necessarily anomalous. A phenomenological study of each class of solutions concerning predictions to the flavour changing neutral current (FCNC) phenomena is also performed: for some solutions, it is possible to obtain realistic quark masses and mixing, while the flavour violating neutral Higgs are light enough to be accessible at LHC.

## I Introduction

The discovery of the Higgs boson in 2012 at the LHC has attested the success of the standard model (SM) in describing the observed fermions and their interactions. However, there exist many theoretical issues or open questions that have no satisfactory answer. In particular, the observed flavour pattern lacks of a definitive explanation, i.e., the quark Yukawa coupling matrices and , which in the SM reproduce the six quark masses, three mixings angles and a complex phase to account for CP violation phenomena, are general complex matrices, not constrained by any gauge symmetry.

Experimentally the flavour puzzle is very intricate. First, there is the quark mass hierarchy in both sectors. Secondly, the mixings in the SM, encoded in the Cabibbo-Kobayashi-Maskawa (CKM) unitary matrix, turns out to be close to the identity matrix. If one takes also the lepton sector into account, the hierarchy there is even more puzzling Emmanuel-Costa:2015tca (). On the other hand, in the SM there is in general no connection between the quark masses hierarchy and the CKM mixing pattern. In fact, if one considers the Extreme Chiral Limit, where the quark masses of the first two generations are set to zero, the mixing does not necessarily vanish Botella:2016krk (), and one concludes that the CKM matrix  being close to the identity matrix has nothing to do with the fact that the quark masses are hierarchical. Indeed, in order to have , one must have a definite alignment of the quark mass matrices in the flavour space, and to explain this alignment, a flavour symmetry or some other mechanism is required Botella:2016krk ().

Among many attempts made in the literature to address the flavour puzzle, extensions of the SM with new Higgs doublet are particularly motivating. This is due to fact that the number of Higgs doublets is not constrained by the SM symmetry. Moreover, the addition of scalar doublets gives rise to new Yukawa interactions and as a result it provides a richer framework in approaching the theory of flavour. On the other hand, any new extension of the Higgs sector must be very much constrained, since it naturally leads to flavour changing neutral currents. At tree level, in the SM, all the flavour changing transitions are mediated through charged weak currents and the flavour mixing is controlled by the CKM matrix Cabibbo:1963yz (); Kobayashi:1973fv (). If new Higgs doublets are added, one expects large FCNC effects already present at tree level. Such effects have not been experimentally observed and they constrain severely any model with extra Higgs doublets, unless a flavour symmetry suppresses or avoids large FCNC Branco:2011iw ().

Minimal flavour violating models Joshipura:1990pi (); Antaramian:1992ya (); Hall:1993ca (); Mantilla2017 (); Buras:2001 (); Dambrosio:2002 () are examples of a multiHiggs extension where FCNC are present at tree-level but their contributions to FCNC phenomena involve only off-diagonal elements of the CKM matrix or their products. The first consistent models of this kind were proposed by Branco, Grimus and Lavoura (BGL) Branco:1996bq (), and consisted of the SM with two Higgs doublets together with the requirement of an additional discrete symmetry. BGL models are compatible with lower neutral Higgs masses and FCNC’s occur at tree level, with the new interactions entirely determined in terms of the CKM matrix elements.

The goal of this paper is to generalize the previous BGL models and to, systematically, search for patterns where a discrete flavour symmetry naturally leads to the alignment of the flavour space of both the quark sectors. Although the quark mass hierarchy does not arise from the symmetry, the effect of both is such that the CKM matrix is near to the identity and has the correct overall phenomenological features, determined by the quark mass hierarchy, Branco:2011aa (). To do this we extend the SM with two extra Higgs doublets to a total of three Higgs . The choice for discrete symmetries is to avoid the presence of Goldstone bosons that appear in the context of any global continuous symmetry, when the spontaneous electroweak symmetry breaking occurs. For the sake of simplicity, we restrict our search to the family group , and demand that the resulting up-quark mass matrix is diagonal. This is to say that, due to the expected strong up-quark mass hierarchy, we only consider those cases where the contribution of the up-quark mass matrix to quark mixing is negligible.

If one assumes that all Higgs doublets acquire vacuum expectation values with the same order of magnitude, then each Higgs doublet must couple to the fermions with different strengths. Possibly one could obtain similar results assuming that the vacuum expectation values (VEVs) of the Higgs have a definite hierarchy instead of the couplings, but this is not considered here. Combining this assumption with the symmetry, we obtain the correct ordered hierarchical pattern if the coupling with gives the strength of the third generation, the coupling with gives the strength of the second generation and the coupling with gives the strength of the first generation. Therefore, from our point of view, the three Higgs doublets are necessary to ensure that there exists three different coupling strengths, one for each generation, to guarantee simultaneously an hierarchical mass spectrum and a CKM matrix that has the correct overall phenomenological features e.g. ,  and denoted here by

Indeed, our approach is within the BGL models, and such that the FCNC flavour structure is entirely determined by CKM. Through the symmetry, the suppression of the most dangerous FCNC’s, by combinations of the CKM matrix elements and light quark masses, is entirely natural.

The paper is organised as follows. In the next section, we present our model and classify the patterns allowed by the discrete symmetry in combination with our assumptions. In Sec. III, we give a brief numerical analysis of the phenomenological output of our solutions. In Sec. IV, we examine the suppression of scalar mediated FCNC in our framework for each pattern. Finally, in Sec. V, we present our conclusions.

## Ii The Model

We extend the Higgs sector of the SM with two extra new scalar doublets, yielding a total of three scalar doublets, as , , . As it was mentioned in the introduction, the main idea for having three Higgs doublets is to implement a discrete flavour symmetry, that leads to the alignment of the flavour space of the quark sectors. The quark mass hierarchy does not arise from the symmetry, but together with the symmetry the effect of both is such that the CKM matrix is near to the identity and has the correct overall phenomenological features, determined by the quark mass hierarchy.

Let us start by considering the most general quark Yukawa coupling Lagrangian invariant in our setup

 −LY=(Ωa)ij¯¯¯¯QLi ˜ϕa uRj+(Γa)ij¯¯¯¯QLi ϕa dRj+h.c., (1)

with the Higgs labeling and are just the usual flavour indexes identifying the generations of fermions. In the above Lagrangian, one has three Yukawa coupling matrices , , for the up-quark sector and three Yukawa coupling matrices , , for the down sector, corresponding to each of the Higgs doublets , , . Assuming that only the neutral components of the three Higgs doublets acquire vacuum expectation value (VEV), the quark mass and are then easily generated as

 Mu =Ω1⟨ϕ1⟩∗+Ω2⟨ϕ2⟩∗+Ω3⟨ϕ3⟩∗, (2a) Md (2b) where VEVs ⟨ϕi⟩ are parametrised as
 ⟨ϕ1⟩=v1√2,⟨ϕ2⟩=v2eiα2√2,⟨ϕ3⟩=v3eiα3√2, (3)

with , and being the VEV moduli and , just complex phases. We have chosen the VEV of to be real and positive, since this is always possible through a proper gauge transformation. As stated, we assume that the moduli of VEVs are of the same order of magnitude, i.e.,

 v1∼v2∼v3. (4)

Each of the couples to the quarks with a coupling which we take be of the same order of magnitude, unless some element vanishes by imposition of the flavour symmmetry. In this sence, each and will generate it’s own respective generation: i.e., our model is such that by imposition of the flavour symmmetry, , , will generate respectivelly , that , , will generate respectivelly , and that , , will generate respectivelly . Generically, we have

 v1∣∣(Ω1)ij∣∣ ∼mu,v2∣∣(Ω2)ij∣∣∼mc,v3∣∣(Ω3)ij∣∣∼mt, (5a) v1∣∣(Γ1)ij∣∣ ∼md,v2∣∣(Γ2)ij∣∣∼ms,v3∣∣(Γ3)ij∣∣∼mb, (5b) which together with Eq. 4 implies a definite hierarchy amongst the non-vanishing Yukawa coupling matrix elements:
 ∣∣(Ω1)ij∣∣ ≪∣∣(Ω2)ij∣∣≪∣∣(Ω3)ij∣∣, (6a) ∣∣(Γ1)ij∣∣ <∣∣(Γ2)ij∣∣≪∣∣(Γ3)ij∣∣. (6b) Next, we focus on the required textures for the Yukawa coupling matrices Ωa and Γa that naturally lead to an hierarchical mass quark spectrum and at the same time to a realistic CKM mixing matrix. These textures must be reproduced by our choice of the flavour symmetry. As referred in the introduction, we search for quark mass patterns where the mass matrix Mu is diagonal. Therefore, one derives from Eqs. (2a), (6a) the following textures for Ωa
 Ω1=⎛⎜⎝x00000000⎞⎟⎠,Ω2=⎛⎜⎝0000x0000⎞⎟⎠,Ω3=⎛⎜⎝00000000x⎞⎟⎠. (7)

The entry means a non zero element. In this case, the up-quark masses are given by , and .

Generically, the down-quark Yukawa coupling matrices must have the following indicative textures

 (8)

We distinguish rows with bold in order to indicate that it is mandatory that at least one of matrix elements within that row must be nonvanishing. Rows denoted with may be set to zero, without modifying the mass matrix hierarchy. These textures ensure that not only is the mass spectrum hierarchy respected but it also leads to the alignment of the flavor space of both the quark sectors Branco:2011aa () and to a CKM matrix . For instance, if one would not have a vanishing, or comparatively very small, entry in the , this would not necessarily spoil the scale of , but it would dramatically change the predictions for the CKM mixing matrix.

In order to force the Yukawa coupling matrices and to have the indicative forms outlined in Eqs. (7) and (8), we introduce a global flavour symmetry. Since any global continuous symmetry leads to the presence of massless Goldstone bosons after the spontaneous electroweak breaking, one should instead consider a discrete symmetry. Among many possible discrete symmetry constructions, we restrict our searches to the case of cycle groups . Thus, we demand that any quark or boson multiplet transforms according to as

 χ→χ′=eiQ(χ)2πNχ, (9)

where is the -charge attributed for the multiplet .

We have chosen the up-quark mass matrix to be diagonal. This restricts the flavour symmetry . We have found that, in order to ensure that all Higgs doublet charges are different, and to have appropriate charges for fields and , we must have . We simplify our analysis by fixing and choose:

 Q(QLi) =(0,1,−2), (10a) Q(uRi) =(0,2,−4), (10b) In addition, we may also fix
 Q(QLi)=Q(ϕi) (11)

It turns out that these choices do not restrict the results, i.e. the possible textures that one can have for the matrices. Other choices would only imply that we reshuffle the charges of the multiplets.

With the purpose of enumerating the different possible textures for the matrices implementable in , we write down the charges of the trilinears corresponding to each as

 Q(¯¯¯¯QLiϕ1dRj)=⎛⎜⎝d1d2d3d1−1d2−1d3−1d1+2d2+2d3+2⎞⎟⎠, (12a) Q(¯¯¯¯QLiϕ2dRj)=⎛⎜⎝d1+1d2+1d3+1d1d2d3d1+3d2+3d3+3⎞⎟⎠, (12b) Q(¯¯¯¯QLiϕ3dRj)=⎛⎜⎝d1−2d2−2d3−2d1−3d2−3d3−3d1d2d3⎞⎟⎠, (12c) where di≡Q(dRi). One can check that, in order to have viable solutions, one must vary the values of di∈{0,1,−2,−3}. We summarise in Table 1 all the allowed textures for the Γa matrices and the resulting Md mass matrix texture, excluding all cases which are irrelevant, e.g. matrices that have too much texture zeros and are singular, or matrices that do not accommodate CP violation. It must be stressed, that these are the textures obtained by the different charge configurations that one can possibly choose. However, if one assumes a definite charge configuration, then the entire texture, Md and Mu and the respective phenomenology are fixed. As stated, the list of textures in Table 1 remains unchanged even if one chooses any other set than in Eqs. (10), (11). As stated, that all patterns presented here are of the Minimal Flavour Violation (MFV) type Joshipura:1990pi (); Antaramian:1992ya (); Hall:1993ca (); Mantilla2017 (); Buras:2001 (); Dambrosio:2002 (). Pattern I in the table was already considered in Ref. Botella:2009pq () in the context of Z8. We discard Patterns IV, VII and X, because contrary to our starting point, at least one of three non-zero couplings with ϕ1 will turn out be of the same order as the larger coupling with ϕ2 in order to meet the phenomenological requirements of the CKM matrix. Notice also, that the structure of other Md’s cannot be trivially obtained, e.g. from Pattern I, by a transformation of the right-handed down quark fields. Our symmetry model may be extended to the charged leptons and neutrinos, e.g. in the context of type one see-saw. Choosing for the lepton doublets Li the charges Q(Li)=(0,−1,2), opposite to the Higgs doublets in Eq. (11), and e.g. for the charges Q(eRi)=(0,−2,4) of the right-handed fields eRi, we force the charged  lepton mass matrix to be diagonal. Then for the right-handed neutrinos νRi, choosing Q(νRi)=(0,0,0), we obtain for the neutrino Dirac mass matrix a pattern similar to pattern I. Of course, for this case, the heavy right-handed neutrino Majorana mass matrix is totally arbitrary. In other cases, i.e. for other patterns and charges, in particular for the right-handed neutrinos, we could introduce scalar singlets with suitable charges, which would then lead to certain heavy right-handed neutrino Majorana mass matrices. Next, we address an important issue of the model, namely, whether accidental U(1) symmetries may appear in the Yukawa sector or in the potential. One may wonder whether a continuous accidental U(1) symmetry could arise, once the Z7 is imposed at the Lagrangian level in Eq. (1). This is indeed the case, i.e., for all realizations of Z7, one has the appearance of a global U(1)X. However, any consistent global U(1)X must obey to the anomaly-free conditions of global symmetries Babu:1989ex (), which read for the anomalies SU(3)2×U(1)X, SU(2)2×U(1)X and U(1)2Y×U(1)X as
 A3≡123∑i=1(2X(QLi)−X(uRi)−X(dRi))=0, (13a) A2≡123∑i=1(3X(QLi)+X(ℓLi))=0, (13b) A1≡163∑i=1(X(QLi)+3X(ℓLi)−8X(uRi)−2X(dRi)−6X(eRi))=0, (13c) where X(χ) is the U(1)X charge of the fermion multiplet χ. We have properly shifted the Z7-charges in Eq.(10) and in Table 1 so that X(χ)=Q(χ), apart of an overall U(1)X convention. In general to test those conditions, one needs to specify the transformation laws for all fermionic fields. Looking at the Table 1, we derive that all the cases, except the first case corresponding to di=(0,0,0), violate the condition given in Eq. (13a) that depends only on coloured fermion multiplets. In the case di=(0,0,0), if one assigns the charged lepton charges as X(ℓLi)=X(QLi) one concludes that the condition given in Eq. (13b) is violated. One then concludes that the global U(1)X symmetry is anomalous and therefore only the discrete symmetry Z7 persists. We also comment on the scalar potential of our model. The most general scalar potential with three scalars invariant under Z7 reads as
 V(ϕ)=∑i[−μ2iϕ†iϕi+λi(ϕ†iϕi)2]+∑i

where the constants , , and are taken real for . Analysing the potential above, one sees that it gives rise to the accidental global continuous symmetry , for arbitrary , which upon spontaneous symmetry breaking leads to a massless neutral scalar, at tree level. Introducing soft-breaking terms like can erase the problem. Another possibility without spoiling the symmetry is to add new scalar singlets, so that the coefficients are effectively obtained once the scalar singlets acquire VEVs.

## Iii Numerical analysis

In this section, we give the phenomenological predictions obtained by the patterns listed in Table 1. Note that, although these patterns arrize directly from the chosen discrete charge configuration of the quark fields, one may further preform a residual flavour transformation of the right-handed down quark fields, resulting in an extra zero entry in . Taking this into account, all the parameters in each pattern may be uniquely expressed in terms of down quark masses and the CKM matrix elements . This follows directly from the diagonalization equation of

 V †Md W=diag(md,ms,mb)⟹Md=V diag(md,ms,mb) W† (15)

with being the CKM mixing matrix, since is diagonal. Because of the zero entries in , it is easy to extract the right-handed diagonalization matrix , completely in terms of the down quark masses and the . Thus, all parameters, modulo the residual transformation of the right-handed down quark fields, are fixed, i.e., all parameters in each pattern may be uniquely expressed in terms of down quark masses and the CKM matrix elements , including the right-handed diagonalization matrix of . More precisely, all matrix elements of are written in terms of Wolfenstein real parameters , , and , defined in terms of rephasing invariant quantities as

 λ≡|Vus|√|Vus|2+|Vud|2,A≡1λ∣∣∣VcbVus∣∣∣, (16a) ¯¯¯ρ+i¯¯¯η≡−VudV∗ubVcdV∗cb (16b) and diag(md,ms,mb) in Eq. (15)
 √mdms=√kdks λmsmb=ks λ2⟹md=kd λ4mbms=ks λ2mb (17)

with phenomenologically, and being factores of order one. Writing in Eq. (15) as , with the vectors formed by the -th column of , we find e.g. for pattern II,

 v3=1n3⎛⎜ ⎜ ⎜⎝mdmbV11msmbV12V13⎞⎟ ⎟ ⎟⎠×⎛⎜ ⎜ ⎜⎝mdmbV31msmbV32V33⎞⎟ ⎟ ⎟⎠ (18)

where is the norm of the vector obtained from the external product of the two vectors. Taking into account the extra freedom of transformation of the right-handed fields, we may choose , corresponding to in Table 1, and we conclude that

 v1=1n1⎛⎜ ⎜ ⎜⎝mdmbV31msmbV32V33⎞⎟ ⎟ ⎟⎠×v∗3 (19)

Obviously, then . This process is replicated for all patterns. Thus, and , are entirely expressed in terms of Wolfenstein parameters and and of Eq. (17). These two matrices will be later used to compute the patterns of the FCNC’s in Table 3. Indeed, in this way, we find e.g. for pattern II,  in leading order order,

 Md=mb ⎛⎜⎝−kd λ3(¯¯¯ρ−i¯¯¯η) A λ30−kd λ2A λ2−ks λ3010⎞⎟⎠ (20)

which corresponds to the expected power series where the couplings in to the first Higgs are comparatively smaller than then couplings in , and these smaller to the couplings in . Similar results are obtained for all patterns in Table 1, except for patterns IV, VII and X, where e.g. for pattern IV, we find that the coupling in is proportional to, which is too large and contradicts our initial assumption that all couplings in to the first Higgs must be smaller than the couplings in to the second Higgs . Therefore, we exclude Patterns IV, VII and X.

We give in Table 2 a numerical example of a Yukawa coupling configuration for each pattern. We use the following quark running masses at the electroweak scale :

 mu =1.3+0.4−0.2MeV,md=2.7±0.3MeV,ms=55+5−3MeV, (21a) mc =0.63±0.03GeV,mb=2.86+0.05−0.04%GeV,mt=172.6±1.5GeV. (21b) which were obtained from a renormalisation group equation evolution at four-loop level 1674-1137-38-9-090001 (), which, taking into account all experimental constrains Charles:2015gya (), implies:
 λ =0.2255±0.0006,A=0.818±0.015, (22a) ¯¯¯ρ =0.124±0.024,¯¯¯η=0.354±0.015. (22b)
 −LNeutral Yukawa= H0v(¯¯¯¯¯¯dLDddR+¯¯¯¯¯¯uLDuuR)+1v′¯¯¯¯¯¯dLNd1(R1+iI1)dR (23) +1v′¯¯¯¯¯¯uLNu1(R1−iI1)uR+1v′′¯¯¯¯¯¯dLNd2(R2+iI2)dR +1v′′¯¯¯¯¯¯uLNu2(R2−iI2)uR+h.c.

where the are the matrices which give the strength and the flavour structure of the FCNC,

 Nd1=1√2V†(v2Γ1−v1eiα2Γ2)W, (24a) Nd2=1√2V†(v1Γ1+v2eiα2Γ2−v21+v22v3eiα3Γ3)W, (24b) Nu1=1√2(v2Ω1−v1e−iα2Ω2), (24c) Nu2=1√2(v1Ω1+v2e−iα2Ω2−v21+v22v3e−iα3Ω3). (24d) Since in our case the Nui are diagonal, there are no flavour violating terms in the up-sector. Therefore, the analysis of the FCNC resumes only to the down-quark sector. One can use the equations of the mass matrices presented in Eq. (2) to simplify the Higgs mediated FCNC matrices for the down-sector:
 Nd1 =v2v1Dd−v2√2(v2v1+v1v2)eiα2V†Γ2W−v2v3v1√2eiα3V†Γ3W (25a) Nd2 =Dd−v2v3√2eiα3V†Γ3W (25b) In order to satisfy experimental constraints arising from K0−¯¯¯¯¯¯¯K0, B0−¯¯¯¯¯¯B0 and D0−¯¯¯¯¯¯¯D0, the off-diagonal elements of the Yukawa interactions Nd1 and Nd2 must be highly suppressed Botella:2014ska () AndreasCrivellin2013 (). For each of our 10 solutions in Table 1, we summarize in Table 3 all FCNC patterns, for each solution, and for v1=v2=v3 and α2=α3=0. These patterns are of the BGL type, since in Eq. (25) all matrices can be expressed in terms of the CKM mixing matrix elements and the down quark masses. As explained, to obtain these patterns, we express the CKM matrix V and the matrix W in terms of Wolfenstein parameters. The tree level Higgs mediated ΔS=2 amplitude must be suppressed. This may allways be achieved if one chooses the masses of the flavour violating neutral Higgs scalars sufficiently heavy. However, from the experimental point of view, it would be interesting to have these masses as low as possible. Therefore, we also estimate the lower bound of these masses, by considering the contribution to B0−¯¯¯¯¯¯B0 mixing. We choose this mixing, since for our patterns, the (3,1) entry of the matrix Nd1 is the less suppressed in certain cases and would require very heavy flavour violating neutral Higgses. The relevant quantity is the off-diagonal matrix element M12, which connects the B meson with the corresponding antimeson. This matrix element, MNP12, receives contributions Botella:2014ska () both from a SM box diagram and a tree-level diagram involving the FCNC:
 M12=MSM12+MNP12, (26)

where the New Physics (NP) short distance tree level contribution to the meson-antimeson contribution is:

 MNP12= 2∑if2BmB96v2m2Ri{[(1+(mBmd+mb)2)(aRi)12]−[(1+11(mBmd+mb)2])(bRi)12} (27) +2∑if2BmB96v2m2Ii{[(1+(mBmd+mb)2)(aIi)12]−[(1+11(mBmd+mb)2])(bIi)12}

with and

 (aRi)12=[(Ndi)∗31+(Ndi)13]2(aIi)12=−[(Ndi)∗31−(Ndi)13]2 ,(bRi)12=[(Ndi)∗31−(Ndi)13]2(bIi)12=−[(Ndi)∗31+(Ndi)13]2 ,i=1,2 (28)

In order to obtain a conservative measure, we have tentatively expanded the original expression in Botella:2014ska () and, for the three Higgs case, included all neutral Higgs mass eigenstates.

Adopting as input values the PDG experimental determinations of , and and considering a common VEV for all Higgs doublets, we impose the inequality . The following plots show an estimate of the lower bound for the flavour-violating Higgs masses for two different patterns. We plot two masses chosen from the set , while the other two are varied over a wide range. In Fig. 1, we illustrate these lower bounds for Pattern III, which are restricted by the entry of matrix and suppressed by a factor of . For Pattern VIII, in Fig. 2 we find the flavour violating neutral Higgs to be much lighter and possibly accessible at LHC.

## V Conclusions

We have presented a model based on the SM with 3 Higgs and an additional flavour discrete symmetry. We have shown that there exist flavour discrete symmetry configurations which lead to the alignment of the quark sectors. By allowing each scalar field to couple to each quark generation with a distinctive scale, one obtains the quark mass hierarchy, and although this hierarchy does not arise from the symmetry, the effect of both is such that the CKM matrix is near to the identity and has the correct overall phenomenological features. In this context, we have obtained 7 solutions fulfilling these requirements, with the additional constraint of the up quark mass matrix being diagonal and real.

We have also verified if accidental symmetries may appear in the Yukawa sector or in the potential, particularly the case where a continuous accidental symmetry could arise, once the is imposed at the Lagrangian level. This was indeed the case, however we shown that the anomaly-free conditions of global symmetries are violated. Thus, the global symmetry is anomalous and therefore only the discrete symmetry persists.

As in this model new Higgs doublets are added, one expects large FCNC effects, already present at tree level. However, such effects have not been experimentally observed. We show that, for certain of our specific implementations of the flavour symmetry, it is possible to suppress the FCNC effects and to ensure that the flavour violating neutral Higgs are light enough to be detectable at LHC. Indeed, in this respect, our model is a generalization of the BGL models for 3HDM, since the FCNC flavour structure is entirely determined by CKM.

###### Acknowledgements.
This work is partially supported by Fundação para a Ciência e a Tecnologia (FCT, Portugal) through the projects CERN/FP/123580/2011, PTDC/FIS-NUC/0548/2012, EXPL/FIS-NUC/0460/2013, and CFTP-FCT Unit 777 (PEst-OE/FIS/UI0777/2013) which are partially funded through POCTI (FEDER), COMPETE, QREN and EU. The work of D.E.C. is also supported by Associação do Instituto Superior Técnico para a Investigação e Desenvolvimento (IST-ID). N.R.A is supported by European Union Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 674896. N.R.A is grateful to CFTP for the hospitality during his stay in Lisbon.

## References

• (1) D. Emmanuel-Costa, N. R. Agostinho, J. I. Silva-Marcos, and D. Wegman, “Novel parametrization for the leptonic mixing matrix and violation,” Phys. Rev. D, vol. 92, no. 1, p. 013012, 2015.
• (2) M. N. R. F.J. Botella, G. C. Branco and J. I. Silva-Marcos, “What if the masses of the first two quark families are not generated by the standard higgs?,” Phys. Rev. D 94 115031, 2016.
• (3) N. Cabibbo, “Unitary Symmetry and Leptonic Decays,” Phys. Rev. Lett., vol. 10, pp. 531–533, 1963. [,648(1963)].
• (4) M. Kobayashi and T. Maskawa, “CP Violation in the Renormalizable Theory of Weak Interaction,” Prog. Theor. Phys., vol. 49, pp. 652–657, 1973.
• (5) G. C. Branco, P. M. Ferreira, L. Lavoura, M. N. Rebelo, M. Sher, and J. P. Silva, “Theory and phenomenology of two-Higgs-doublet models,” Phys. Rept., vol. 516, pp. 1–102, 2012.
• (6) A. S. Joshipura and S. D. Rindani, “Naturally suppressed flavor violations in two Higgs doublet models,” Phys. Lett. B, vol. 260, pp. 149–153, 1991.
• (7) A. Antaramian, L. J. Hall, and A. Rasin, “Flavor changing interactions mediated by scalars at the weak scale,” Phys. Rev. Lett., vol. 69, pp. 1871–1873, 1992.
• (8) L. J. Hall and S. Weinberg, “Flavor changing scalar interactions,” Phys. Rev. D, vol. 48, pp. R979–R983, 1993.
• (9) S. Mantilla and R. Martinez, “3hd+1hs u(1)’ family non-universal model,” arXiv 1704.04869, 2017.
• (10) M. G. S. J. A.J. Buras, P. Gambino and L. Silvestrini, “Universal unitarity triangle and physics beyond the standard model,” Phys.Lett. B500 (2001) 161-167, 2001.
• (11) G. I. G. D’Ambrosio, G.F. Giudice and A. Strumia, “Minimal flavor violation: An effective field theory approach,” Nucl.Phys. B645 (2002) 155-187, 2002.
• (12) G. C. Branco, W. Grimus, and L. Lavoura, “Relating the scalar flavor changing neutral couplings to the CKM matrix,” Phys. Lett. B, vol. 380, pp. 119–126, 1996.
• (13) G. C. Branco and J. I. Silva-Marcos, “Invariants, Alignment and the Pattern of Fermion Masses and Mixing,” Phys. Lett., vol. B715, pp. 315–321, 2012.
• (14) F. J. Botella, G. C. Branco, and M. N. Rebelo, “Minimal Flavour Violation and Multi-Higgs Models,” Phys. Lett. B, vol. 687, pp. 194–200, 2010.
• (15) K. S. Babu and R. N. Mohapatra, “Quantization of Electric Charge From Anomaly Constraints and a Majorana Neutrino,” Phys. Rev. D, vol. 41, p. 271, 1990.
• (16) K. Olive and P. D. Group, “Review of particle physics,” Chinese Physics C, vol. 38, no. 9, p. 090001, 2014.
• (17) J. Charles et al., “Current status of the Standard Model CKM fit and constraints on New Physics,” Phys. Rev., vol. D91, no. 7, p. 073007, 2015.
• (18) F. J. Botella, G. C. Branco, A. Carmona, M. Nebot, L. Pedro, and M. N. Rebelo, “Physical Constraints on a Class of Two-Higgs Doublet Models with FCNC at tree level,” JHEP, vol. 07, p. 078, 2014.
• (19) A. A. Crivellin and C. Greub, “Flavor-phenomenology of two-higgs-doublet models with generic yukawa structure,” Phys.Rev. D87 (2013) no.9, 094031, 2013.
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters