Multi-Higgs-doublet models and singular alignment

We consider a 4-Higgs-doublet model in which each Higgs doublet gives mass to one of the fermion sets {mt}, {mb, mτ, mc}, {mμ, ms}, and {md, mu, me}. The sets have the feature that within each of them the masses are similar. Our model explains the mass hierarchies of the sets by hierarchies of the vacuum expectation values of the Higgs doublets associated to them. All Yukawa couplings are therefore of order one. Neutrino masses are generated by a type-I seesaw mechanism with PeV-scale singlet neutrinos. To avoid the appearance of tree-level flavour changing neutral currents, we assume that all Yukawa matrices are singularly aligned in flavour space. We mean by this that the Yukawa matrices are given as linear combinations of the rank 1 matrices that appear in the singular value decomposition of the mass matrix. In general, singular alignment allows to avoid flavour changing neutral currents in models with multiple Higgs doublets.


Introduction
An understanding of fermion masses and mixing is still lacking. In particular, the mass values display unexplained patterns and hierarchies; this is the case when one considers the three generations as well as the species: 1 (2) ? m ν2 (1) ? m ν1(3) 1 That is, any of the four masses within the same generation.

JHEP07(2019)036
We can summarize the situation by asking the following questions: • Why is the top quark mass the only fermion mass of the order of the electroweak (EW) scale, m t ≈ v EW with v EW 174 GeV?
• Why is the top quark mass so much heavier than the rest of fermion masses, m t m f ?
• Why do all charged fermions satisfy the hierarchy, m 3 m 2 m 1 ?
• Why are for the first generation the masses (except for neutrinos) closer to each other than for the other two generations, m d ∼ m u ∼ m e versus m c m s ∼ m µ and m t m b ∼ m τ ?
• What could the interspecies hierarchy, e.g. m t m b > m τ m ν3 , be telling us?
• Why are neutrino masses much smaller than the charged fermions, 2 m ν ∼ 10 −7 m e ?
This is commonly referred to as the problem of mass [1]. Part of the mystery lies in the contrast of expecting Yukawa couplings to be order one, y f = O(1), whereas the observed values with a single Higgs doublet are much smaller than 1, except for the top quark, y f 1 (f = t). In the following, we assume Yukawa couplings to be order one, y f = O(1), and try to understand the fermion mass patterns through a theory with multiple Higgs doublets. The most extreme approach along this line would be the "private Higgs" scenario, in which among other things, for each fermion a Higgs doublet is introduced [2,3], see also [4][5][6]. The mass hierarchies are explained by hierarchies of vacuum expectation values of the individual Higgs doublets: m f v f , where v f is the vacuum expectation value of the Higgs that is responsible for the fermion f = u, d, c, s, t, b, e, µ, τ .
In general, in a model with N Higgs doublets, Φ i (i = 1, 2, . . . , N ), where each of their neutral components acquires a vacuum expectation value (vev), Φ 0 j = v j e iθ j , a relation among these vacua is satisfied: Here v EW 174 GeV, v i ≥ 0, and all doublets share the same hypercharge Y = 1 2 . Now, if we consider that each single Higgs is fully responsible for the mass of one single fermion (where N should equal the number of fermions in the theory), then the previous relation is modified to

JHEP07(2019)036
Furthermore, if we consider that Yukawa couplings should be order one numbers, y f = O(1), we could approximately say that, to good approximation In the case of the Standard Model (SM), with N = 12 fermions, the previous equation is fulfilled. We will call this relation the mass-vacuum relation. An amusing possibility from this relation is that if all N doublets have the same vev, one would have N fermions with mass of about 174/ √ N GeV, which would be about 50 GeV for 12 fermions. If two doublets have vev v EW / √ 2 and the rest a vanishing vev, then there would be two fermions with mass v EW / √ 2 123 GeV. In turn, if only one doublet has a vev, there is only one fermion with mass v EW . Forcing the mass-vacuum relation to be fulfilled and assuming that only one Higgs acquires a vev leaves hardly any mass for the other fermions and explains the top quark's dominance. Moreover, this same argument could help us to understand why neutrinos are so light when assumed as Dirac fermions. 3 The particle content in the main scenario discussed in this paper is smaller than that for a private Higgs-like scenario. Our observation is that the fermion masses can be grouped into four different sets: In each set the masses are quite similar and can in fact be explained by similar O(1) Yukawa couplings to an individual Higgs doublet Φ t,b,µ,d . Such a 4-Higgs-Doublet Model has to the best of our knowledge not been considered before. We find several attractive and testable features of the model, and demonstrate that it is not in conflict with measured Higgs couplings and other tests. Our model traces the hierarchy of the mass values of the different fermion sets to hierarchies of vevs of their respective Higgs doublets. We show that the smaller vevs can be induced by the larger vevs, and the hierarchy among them arises because the four vevs are protected by different symmetries.
The main problem in multi-Higgs doublet models is of course the presence of flavour changing neutral currents (FCNC). Theories which through the use of symmetries naturally avoid those FCNC are said to possess Natural Flavour Conservation (NFC). Options to evade FCNC include, next to arranging the additional scalar particles to be very heavy, suppressing dangerous Yukawa couplings [9], separating the Yukawa matrices such that only one scalar doublet couples to a given right-handed fermion field [10,11], or Yukawa alignment [12,13], in which the different Yukawa matrices are proportional to each other. As a proof of principle that FCNC can be entirely avoided in our setup, we assume here another solution. We note that if the Yukawa matrices are proportional to any of the rank-one matrices that appear in the singular value decomposition of the fermion mass matrices, FCNC are absent. We denote this as "singular alignment".
The paper is organized as follows: in section 2 we present singular alignment and discuss some of its features. The model with four Higgs doublets to explain the masses of the individual sets {m t }, {m b , m τ , m c }, {m µ , m s }, and {m d , m u , m e } is presented and JHEP07(2019)036 analyzed in section 3. Conclusions are presented in section 4, and some technical details are delegated to appendices.

Singular alignment
In general, having multiple Higgs doublets coupling to fermions with the same electric charge will produce tree-level FCNC, which are experimentally strongly constrained. Three main possibilities to overcome this problem have typically been studied: (i) assume "dangerous" Yukawa couplings to be sufficiently suppressed at tree-level [9]; (ii) assume the corresponding Yukawa matrices of each type of fermion (up-type quarks, down-type quarks and charged leptons) to be proportional to the mass matrix [12,13]; (iii) impose an adequate symmetry such that each fermion type couples exactly to one of the doublets [10,11]. In the following, we comment only on the last two possibilities and introduce singular alignment.
Let us start from the most general case for a Yukawa Lagrangian in a N HDM, where F L and f R are three dimensional vectors in family space and transform as a doublet and as a singlet under SU(2) L , respectively. The N Higgs doublets acquire a vev, v a = Φ 0 a . In general, Yukawa couplings will couple all fermions to all Higgses. Therefore, the most general form of a mass matrix is Each Yukawa matrix, Y i , is a 3 × 3, arbitrary, and complex matrix with rank 3. The appearance of tree-level FCNC is automatic within this setup as diagonalization of the mass matrices does not mean, in general, simultaneous diagonalization of the individual Yukawa matrices. However, to avoid introducing dangerous tree-level FCNC the following can be done: NFC theories. Adequate symmetries are imposed in such a way that each of the three charged fermions will only couple to a single Higgs [10,11], i.e. for each fermion type holds where no sum over k is intended. In this case diagonalization of the l.h.s. means diagonalization of the r.h.s. . For N Higgs doublets, the easiest way to achieve this is via a symmetry of the form where in order for this symmetry to be realizable = N − 1 should hold. Realizable symmetries are a set of allowed discrete symmetries of the scalar potential which have no accidental larger groups that could give rise, for example, to massless Goldstone bosons [14].

JHEP07(2019)036
Now, before turning to the next possiblity, let us comment on the Singular Value Decomposition (SVD) of a mass matrix: (2.5) Here L and R are unitary matrices which rotate independently the left-and right-handed fermion fields and Σ = diag(m 1 , m 2 , m 3 ) with m i > 0. Realize that the SVD may also be written as a sum of three rank 1 matrices, where P i are three projector operators, P 2 i = P i and i P i = 1 3×3 , which have the form In the following, we will denote each rank 1 matrix appearing in the SVD by and call it singular matrix.
Yukawa alignment. As each Yukawa term in eq. (2.2) is a rank 3 matrix, a second possibility to avoid FCNC, is to assume that each of them is proportional to the full SVD [12,13]: Here the ζ i are real and we have that Diagonalization of the l.h.s. means diagonalization of the r.h.s. . This is understandable as each Yukawa matrix is rank 3 and thus if related to the singular matrices should be composed of the three independent singular matrices. Furthermore, one has the constraint N j=1 ζ j = 1 . (2.11) Singular alignment. A more general scenario is that in which each Yukawa matrix is given by a linear combination of the singular matrices, i.e.

JHEP07(2019)036
Appendix A gives a straightforward proof of the absence of FCNC in case the Yukawa matrices take this form. Comparing with the full mass matrix, which can be written as Hence, all fermion masses are independent linear combinations of the different vevs and all Higgs doublets can be responsible for giving mass to all fermions. In practice, models may also lead to Yukawa matrices with ranks less than 3. In this case the singular alignment can still hold and the only new difference would be to have some of the constants η i , Ω i , Λ i appearing in eq. (2.12) equal to zero. In short, singular alignment is the very strong Ansatz of choosing Yukawa matrices to be related to the rank 1 matrices appearing in the SVD. Through this alignment, no treelevel FCNC appear for any number of Higgs doublets. Let us consider now some explicit examples.

The two-fermion family case
We assume N Higgs doublets for two generations of charged fermions. In this case, the mass matrix is If no symmetry is imposed all Higgs doublets are allowed to couple to our 6 fermions. Therefore, all the Yukawa matrices are rank 2. Let us implement the singular alignment. For this purpose, the SVD of the mass matrix is written as where, without any loss of generality, we have chosen to work in the basis where the righthanded fermions have been already transformed and we have explicitly written the most general expression for a unitary matrix in two dimensions. The two singular matrices are Singularly aligning our Yukawa matrices in flavour space means

JHEP07(2019)036
We identify the masses as Regarding FCNC, note that in the mass basis we have Hence, no FCNC are introduced in this model of N > 1 Higgs doublets that couple to all fermions. Also, if matrices of lower rank are obtained through the use of convenient symmetries, then our general expressions in eq. (2.20) will still hold but with some of the parameters η i or Ω i vanishing.

The three-fermion family Case
The next example deals with three generations and three Higgs doublets. To be singularly aligned, each rank 1 Yukawa matrix should be seen as a column vector satisfying unitarity conditions (recall eq. (2.12)): Here we have denoted (a 1 , a 2 , a 3 ) T ≡ |a and similarly for the other columns. Notice that the Yukawa couplings should not enter into these expressions. A practical way to implement all these conditions is to make use of an explicit parametrization of a unitary matrix. Then, a singularly aligned mass matrix could take the form where T and Q are diagonal phase matrices with two phases each and we have used the shorthand notation for the sine and cosine functions. Recall that a 3 × 3 unitary matrix possesses 6 complex phases (one of which is global) and 3 real parameters. Each column is proportional to a given singular matrix. At last, realize that masses and mixing get completely decoupled when singularly aligning the Yukawa matrices. Recall in this context that any set of singular vectors corresponding to a set of non-degenerate singular values is always orthonormal.

Hierarchical fermion masses
A shared feature among all the charged fermions is that their masses are hierarchical, To theoretically understand this in a N HDM with singular alignment, see eq. (2.14), one must understand under what conditions this property gets always realized. We are not interested in any fine-tuned scenario where through adequate values for the set of parameters {η, Ω, Λ} we generate hierarchical masses, we are assuming that Furthermore, we are actually interested in the minimal number of scalar doublets necessary to explain all the observed patterns in the fermion masses. For the moment, notice that one possibility is to couple a single Higgs to each different flavour with the same electric charge.
In this case we have It is obvious then that the only way to achieve hierarchical masses with This fact is connected to the mass-vacuum relation. The maximal setup, if neutrinos are assumed as Dirac particles, would require 12 Higgs doublets. However, this large number of scalars can be significantly reduced if one notices that among the different masses there are majorly 4 (5) mass scales, where the (5) corresponds to Dirac neutrino masses. This is what we will deal with in section 3. In case of Majorana neutrinos there are four possibilities depending on from which Higgs doublet the Dirac mass matrix of the type-I seesaw mechanism stems. We will come back to this point later. Of course, neutrino mass could also be independent of the Higgs doublets.

Beyond singular alignment
If a small amount of flavour violation via neutral mediators is permitted, then a less restrictive venue can be obtained through the following conditions: (i) the third Yukawa matrix for all fermion species is the only rank 1 matrix and proportional to the third singular matrix, (2.28) (ii) the first and second Yukawa matrices are no longer proportional to the singular matrices, so they may in general produce FCNC; (iii) however, to produce a hierarchy between the first and second generation, the second Yukawa matrix should be at most rank 2 and have no contributions to the first family masses; (iv) the first Yukawa matrix can be rank 3, 2 or 1. In other words, the three Yukawa matrices should imply the sequential symmetry breaking chain

JHEP07(2019)036
where F might either be baryon or lepton number. The introduction of flavour violation as allowed by the two lightest families means no risk as this set of flavour transitions will be sequentially suppressed by the approximately conserved symmetries at each step.

Radiative stability
In the absence of a specific symmetry protection, one-loop quantum corrections may induce misalignment in the different singularly aligned Yukawa matrices and bring about FCNC's at the loop level. It is important to know if this effect is small and compatible with current experimental constraints. The study of this issue can be directly related to the work of ref. [13] wherein the issue of radiative stability was investigated for the most generalized Yukawa aligned-like form given by where Ξ i is a complex 3 × 3 matrix subject to the condition This generalized Yukawa-alignment means breaking flavour universality. Notice that the normal Yukawa-alignment, eq. (2.9), can be recovered when all diagonal elements in the r.h.s. of eq. (2.31) are equal (flavour universal). In ref. [13], it was shown that the induced misalignment is a quite small effect, as the initial alignment in the multi-Higgs Lagrangian has some residual flavour symmetries, which tightly limit the type of FCNC operators that can be generated at higher orders. This can be easily understood as the Yukawa alignment is a linear realization of the minimal flavour violation hypothesis [15] and could be derived from it [16]. This hypothesis states that the only source of flavour breaking should come from the Yukawa matrices, even in the presence of new particles and interactions [17][18][19][20].
The previous discussion also applies to the Singular Alignment as it is possible to show that it is equivalent to the generalized Yukawa-alignment via substitution in eq. (2.30) of the relations Therefore, the ansatz of singularly aligning Yukawa matrices in flavour space, in order to avoid FCNC's at tree level, has a sufficiently small misalignment, induced by one-loop quantum corrections, consistent with all known phenomenological tests.

The minimal setup: a 4HDM
Now we discuss a 4HDM which takes into account that among the measured fermion masses four different sets can be identified: for the masses in their respective sets. 4 The corresponding mass-vacuum-like relation in analogy to eq. (1.3) would take the form The model can be constructed by first imposing fields to transform under the symmetry Z 2 × Z 2 × Z 2 , as shown in table 1. The Yukawa Lagrangian implied by the charge assignment is The way in which we have employed the charge assignment to couple fermions with Higgs doublets has given us a model where all Yukawa matrices for the charged fermions are rank Similar expressions can be given for the other fermion species. Notice we are employing a conventional notation for the Yukawa couplings, y f i , in order to distinguish at this point generic Yukawa matrices from those which have been singularly aligned. Now, to singularly align these matrices, we demand that each column should be given by a single singular matrix (in order to have a hierarchy of masses with order one Yukawa couplings, cf. section 2.3), in our up-type example this means: The explicit form of these singular matrices was given in section 2.2, they correspond to one of the three columns in eq. (2.24). We can also write them as ∆ u,c,t = L † P 1,2,3 R, see the discussion around eq. (2.8).
The model presented here arranges that a certain Higgs doublet will couple to a given set of fermions, even if they possess different electric charge. All corresponding Yukawa matrices will already be rank 1. Through the special requirement that Yukawa matrices should be singularly aligned in flavour space, as discussed in section 2, it is possible to avoid flavour violation at tree-level.

JHEP07(2019)036
The model allows to reproduce fermion mixing, as shown in appendix B. Neutrino masses are generated via the type-I seesaw mechanism. We have associated the three right-handed neutrinos to the Higgs doublet Φ d . This implies that via the type-I seesaw mechanism the heavy neutrino mass scale M should be around PeV, where we have assumed that m D Φ 0 d 10 MeV and m ν m 2 D /M 0.1 eV. The contributions to some lepton-flavor-violation processes coming from the admixture of the heavy righthanded neutrinos with the left-handed ones can already be estimated via the standard formulae of type-I seesaw models [21]. This calculation is greatly simplified in the limit M M W , which is our case. 5 The following upper bound to various processes of interest may be obtained: B th r (µ → eγ) < 10 −14 , B th r (µ → 3e) < 10 −18 , B th r (τ → 3µ) < 10 −7 , and B th r (τ → µγ) < 10 −4 . Notice how, in general, these numbers will still get suppressions by small mixing-like angles of the order of m D /M ∼ 10 −8 times order one numbers (at most) arising from corresponding form factors. The present experimental upper limits on these decays at 90% C.L. are given by: B exp r (µ → eγ) < 4.2×10 −13 [22], B exp r (µ → 3e) < 10 −12 [23], B exp r (τ → 3µ) < 4.6×10 −8 [24] and B exp r (τ → µγ) < 3.3×10 −8 [25]. The smallness of the estimated branching ratios is of no surprise, as the high-scale type-I seesaw is known for giving very suppressed rates, see for example [26] and references therein. On the other hand, possible contributions coming from the scalar mediators at the loop level can also be expected to be sufficiently small and consistent with phenomenological tests as suggested by analyses of the minimal lepton-flavour violation hypothesis [27,28] and as discussed in section 2.5. Notice, however, that this set of flavour-violating processes have a strong dependence in the ratio between the two scales (Λ LN /Λ LFV ) 4 , where the first and second one correspond to the scale where lepton number (LN) is broken and lepton-flavour-violation (LFV) is produced. Therefore, if a large hierarchy exists between these two scales one may obtain observable effects. In our case, we may estimate this ratio as (v d M/v 2 τ ) 2 ∼ 10 8 which is still sufficiently small. For example, after substitution in B th r (µ → eγ) = 1.6 × 10 24 (Λ LN /Λ LFV ) 4 we obtain B th r (µ → eγ) ∼ 10 −16 , where the previous relation was taken from ref. [27]. Hence, we see that for the particular purposes of this work the rates for LFV processes are expected to be in agreement with the current upper bounds.

The scalar potential
The most general, renormalizable and gauge invariant scalar potential of the model is (3.5) Here a, b = t, b, µ, d, and for the sake of simplicity we are assuming all couplings to be real.

JHEP07(2019)036
In order to generate a hierarchy among the vevs we choose the particular case where µ 2 t < 0 and µ 2 b,µ,d > 0 , (3.6) such that the only Higgs acquiring a vev is Φ t : We are following the convention Φ 0 . As Φ t has no charge under any of the three Abelian symmetries, see table 1, its vev preserves the symmetry. Equivalently, the symmetries are protecting the other scalars from acquiring a vev. Thereafter, through the following subset of soft-breaking terms, where µ 2 ab v 2 t , µ 2 a , we induce vevs for the other three Higgs doublets. To be more specific, the particular choice of soft-breaking terms is motivated by the fact that each of them will only break a particular piece of the whole symmetry. That is, only break Z 2 , Z 2 , and Z 2 , correspondingly. Therefore, once the EW symmetry is spontaneously broken, the first soft-breaking term will induce a vev to Φ b which in return will induce a vev to Φ µ until finally reaching Φ d . It is possible to show that within this limit the minimization conditions are satisfied if the vevs are given as together with eq. (3.7) and where (XY Z) ab = X ab + Y ab + Z ab and µ 2 ab < 0. By virtue of this choice, the vevs are naturally small and obey the desired hierarchy

Fermionic couplings to the SM-like Higgs
The introduction of the soft breaking terms (eq. (3.8)) in the Higgs potential will produce a small mixing among the four Higgs doublets. For the moment, let us focus on the neutral scalars. We assume all parameters in the scalar potential to be real. Through this choice we consider it to be CP -symmetric. Hence, no admixture between the real and imaginary components of the neutral fields is allowed as they have definite CP quantum numbers. To compute their couplings to all fermions we start from the Yukawa Lagrangian in the mass basis which is written as where we have changed our notation {η, Ω, Λ} to the conventional one, y f . We can bring the CP -even scalar sector to its mass basis via where R is an orthogonal matrix, R T R = RR T = 1 4×4 , and h 0 is the lightest state with a mass of m h 0 125 GeV. Now, in order to find out how fermions couple to the SM-like Higgs, h 0 , we substitute φ k = R 1k h 0 in eq. (3.11) to obtain We can define the following four classes of fermion-scalar couplings: sin α 1 sin α 2 sin α 3 . (3.14) The angles α i in these relations are We note an attractive and testable feature of the model, namely that the couplings between fermions and the SM-like Higgs are modified in the same way for each set. That is, the couplings of the sets {m t }, {m b , m τ , m c }, {m µ , m s }, and {m d , m u , m e } are changed with respect to the SM-case by the same amount for each set, see figure 2. The coupling to the top quark is always essentially SM-like, ξ t h = 1. This is understood because R 11 and cos α 1 are both very close to 1, which is caused by the vev hierarchy v t v b,µ,s . Notice that, even though the mixing R ik in all cases is proportional to the soft-breaking parameters, the implied smallness in |R ik | may be compensated by α i 1, and therefore, in general, ξ f h should not be expected to be small. In fact, within this scenario we can have four different possibilities: (i) hyper-couplings with ξ f h > 1, aligned couplings with JHEP07(2019)036 ξ f h = 1, hypo-couplings with ξ f h < 1, and a mixture of any of these (the ξ f h can even be negative). One has to confront the couplings in this model with present measurements of Higgs couplings. We adopt the following numbers from combined fits of data taken at √ s = 13 TeV [29]: No useful information about the couplings to first and second generation fermions exist, except for the muon, where the uncertainties are nevertheless very large. In our case κ Z,W can be reproduced as in any multi-Higgs doublet model. The values of κ t,τ,b need to be compared with our ξ f h , which is what the plots in figure 2 do for the four benchmark scenarios to be discussed next. 6

Numerical examples
A thorough analysis of the Higgs potential is beyond the scope of this work, nevertheless, we will present four numerical benchmark scenarios. They obey the following conditions and constraints:
• Vacuum stability: where (a = b, µ, d). This set of conditions was computed from the requirement that the squared mass matrices for the charged scalars and pseudo-scalars should be positive definite, for further details see appendix C.
• Contributions to the ρ parameter: that is, it should be consistent with the maximum allowed deviation from the SMexpectation [34]. The first and second uncertainty originates whether the oblique parameter U is fixed to zero or not within the multi-parameter fit. For our calculations, we employ the one-loop contribution coming from a generic N -Higgs doublet model obtained in ref. [35], for further details see appendix D. Our analysis is consistent with ref. [36] where the interplay between the maximum number of N -Higgs doublets and their allowed masses in the oblique parameters is discussed.
• Charged Higgs masses above the lower bound [37]: • Recently, a search for a Higgs-like particle, φ, decaying into a pair of bottom quarks with at least one additional bottom in proton-proton collision was reported [38]. The following mass range was excluded with 95% confidence level: 100 GeV < m φ < 300 GeV . While not directly comparable with our scenario, our benchmark points nevertheless obey this constraint.  Table 3. Outputs for each of the four numerical benchmark scenarios. Scalar masses are given in GeV.

Conclusions
Within the SM the huge hierarchy of Yukawa couplings remains a puzzle. In this regard we used the fact that the observed fermion masses indicate that the following sets have similar Yukawa couplings: {m t }, {m b , m τ , m c }, {m µ , m s } and {m d , m u , m e }. We have shown that a 4HDM can be constructed that explains this feature. Each set of fermions has its own Higgs doublet. Their vevs are hierarchical which explains the mass hierarchy of the sets. In the model a flavour symmetry was introduced to generate rank 1 Yukawa matrices. Soft breaking was included in the potential, which makes it possible to induce the smaller vevs by the larger ones, where each smaller vev corresponds to a different broken symmetry, and is thus protected by it. All Yukawa couplings take on "natural" values of order 1.
In the model neutrino masses are generated via a type-I seesaw mechanism with a Dirac neutrino mass matrix of order of the down-quark mass scale, hence the right-handed singlet Majorana neutrinos are of PeV-scale. We have demonstrated that fermions of a given set couple to the SM-like Higgs with the same modified factor. In this regard, the clearest signal for this kind of models is to investigate their coupling to the SM-like Higgs and determine if they are grouped. The top quark couples to the SM-like Higgs essentially with the same strength as in the SM. Benchmark scenarios with definite predictions for those couplings as well as for scalar masses were provided.
Multi-Higgs doublet models face of course problems with FCNC. By singularly aligning the Yukawa matrices we have shown explicitly that those can be evaded. This alignment JHEP07(2019)036 assumes that the Yukawa matrices are related to the rank 1 matrices that appear in the singular value decomposition of the mass matrices. In this manner, it is in general, not only in our model, possible to avoid FCNC while simultaneously coupling several Higgs doublets to an individual given fermion. Moreover, its equivalence to the most general Yukawa alignment also allows us to state that its misalignment at the one-loop level is sufficiently small and is consistent with all known phenomenological tests.
The model as well as aspects of singular alignment allow for several follow-up studies regarding both model building and phenomenology.
A Proof of FCNC disappearance in the singular basis Consider a given fermion type coupled to N different scalar doublets. Its mass matrix would be given by where we have assumed that each scalar doublet acquires a vev. On the other hand, the SVD of the mass matrix is where L and R are unitary transformations acting independently on the left-and righthanded fields. Using Dirac notation, the SVD can be rewritten as Singular alignment requires assuming each Yukawa matrix to be related to the rank 1 matrices, | i r i |, of the SVD. In general, we can express the Yukawa matrices as a linear combination of the rank one singular matrices where the parameters {η, Ω, Λ} are real.
In the mass basis, each Yukawa matrix would take the form, Therefore, through singularly aligning we have avoided the appearance of dangerous treelevel FCNC. For last, notice that after substitution of the previous relation in eq. (A.1) we obtain

JHEP07(2019)036
B Numerical example for quark mixing The following singular matrices allow us to reproduce exactly the observed mixing in the quark sector as recently reported in the PDG 2018 [34]: By virtue of those conditions one may easily derive from the charged scalar matrix that: For the pseudo-scalar matrix: , (a = b, µ, d) . (C.8) We have not neglected here the off-diagonal contributions to the CP -even scalar matrix, as even though they are very small, they can still influence the Higgs-fermion couplings as already previously discussed.

D Constraints from the ρ parameter
The one-loop level contribution to the ρ parameter from a theory with N Higgs doublets has been calculated in ref. [35] and is expressed as where F (x, y) ≡  an interesting property, as it grows linearly with max(x, y), that is, quadratically with the heaviest-scalar mass, when that mass becomes very large. As long as the difference in the scalar masses is small, δ 200 GeV, the maximum value of this function lies within the 3σ deviation in ∆ρ, as shown in figure 3, even if the masses become very heavy, M > 300 GeV. This can be seen from Taylor expanding the function, δ x, Moreover, realize that these contributions will get further suppressed by the factors coming from the off-diagonal matrix elements in the product of the orthogonal matrices. In the limit in which we are working, mass matrices can be considered to a very good degree of accuracy to be diagonal, therefore, the maximum amount of contributions in this model will take the form