Natural Quasi-Alignment with two Higgs Doublets and RGE Stability

In the context of two Higgs doublet models, we study the conditions required in order to have stable quasi-alignment in flavour space. We show that stability under the RGE imposes strong constraints on the flavour structure of the Yukawa couplings associated to each one of the Higgs doublets. In particular, we find a novel solution, where all Yukawa couplings are proportional to the so-called democratic matrix. This solution is rather unique, since it is the only stable solution which is a good starting point for reproducing the observed pattern of quark masses and mixing. We also show that this stable solution can be obtained by imposing on the Lagrangian a $Z_3 \times Z^\prime_3$ flavour symmetry. Quark masses of the lighter quark generations are generated through the breaking of this discrete symmetry, and, at this stage, scalar mediated FCNC arise, but are naturally suppressed by the smallness of light quark masses. In this way, we relate Higgs alignment to the hierarchy of quark masses through a discrete family symmetry.


Introduction
One of the simplest extensions of the Standard Model (SM) consists of the addition of scalar doublets to the SM spectrum. Multi-Higgs extensions arise in a variety of frameworks, including supersymmetric extensions of the SM, as well as models with family symmetries. A two Higgs doublet model (2HDM) was first introduced by Lee [1], in order to achieve spontaneous breaking of the CP symmetry. If no extra symmetries are introduced, 2HDMs lead to too large tree level scalar mediated flavour-changing-neutral-currents (FCNC) [2]. In order to avoid these potentially dangerous currents, various schemes have been proposed: i) Glashow and Weinberg [3] have pointed out that one can avoid FCNC at tree level by introducing a Z 2 symmetry under which the two Higgs doublets transform differently. The introduction of a Z 2 symmetry in 2HDMs prevents the generation of spontaneous CP breaking [4] unless the symmetry is softly broken [5]. ii) Pich and Tuzon [6] have conjectured the existence of flavour alignment of the two Yukawa matrices, thus avoiding FCNC at tree level. This is an interesting suggestion, but it has the drawback of being an adhoc assumption, not explained by any symmetry. Furthermore, it has been pointed out that in general this scheme is not stable under the renormalization group [7]. There have been attempts at obtaining alignment in various extensions of the SM [8], [9], [10]. iii) Another possibility has been proposed some time ago [11] by Branco, Grimus and Lavoura (BGL) who have pointed out that there is a symmetry which, when imposed on the Lagrangian, constrains the Yukawa couplings in such a way that FCNC do arise at tree level, but are entirely determined by the V CKM matrix, with no other free parameters. In some of the BGL models, one has a strong natural suppression of the most dangerous FCNC, with, for example, the strangeness changing neutral currents, proportional to (V td V * ts ) 2 , which implies a very strong natural suppression of the contribution to the K 0 −K 0 transition. With this suppression, the neutral Higgs masses need not be too large. BGL models have been extended to the leptonic sector [12], their relation to Minimal Flavour Violation models has been studied [13] and their phenomenological implications have been recently analysed [14].
In this paper, we reexamine the question of the stability of flavour alignment under the renormalization group. Assuming that the Yukawa couplings of the two Higgs doublets are aligned, i.e., proportional to each other, we study under what conditions the alignment is maintained by the renormalization group. Apart from the conditions already found in Ref. [7], we find new solutions which can be of great physical interest. One of these solutions, corresponds to having all the Yukawa coupling matrices proportional to the so-called democratic matrix [15], [16], [17]. This solution is rather unique, since on the one hand it is stable under the renormalization group equations (RGE) and on the other hand, it is the only stable solution which provides a good starting point for reproducing the observed pattern of quark masses and mixing. We then point out that this flavour democratic solution can be obtained as a result of a Z 3 × Z 3 flavour symmetry. In the framework that we propose, flavour alignment is exact in the limit where only the third family acquires mass. Once the two light generations acquire a mass, there are small deviations from alignment, which are suppressed by the strong hierarchy of quark masses. As a result, one obtains in this framework, a quasi-alignment of the Yukawa couplings, as a result of the Z 3 × Z 3 symmetry, together with the strong hierarchy of quark masses.
The paper is organised as follows. In the next section, we briefly describe the general flavour structure of the 2HDM, in order to settle our notation. In section 3 we derive all the solutions for the Yukawa couplings, leading to alignment, stable under the renormalization group. In section 4 we show that the flavour democratic solution can be obtained as a result of a Z 3 × Z 3 flavour symmetry and propose an ansatz for the breaking of the Z 3 × Z 3 symmetry. In section 5, we examine the suppression of scalar mediated FCNC in our framework. In section 6, we perform a numerical analysis, showing how the pattern of quark masses and mixing can be obtained in the framework of our ansatz. Finally our conclusions are contained in section 7. In the Appendix we present a full study of the solutions of the alignment conditions.

Yukawa Couplings in the General Two-Higgs-Doublet-Model (2HDM)
For completeness and in order to establish our notation we briefly review the flavour structure of the 2HDM, when no extra symmetries are introduced in the Lagrangian. The Yukawa couplings can be written: where Φ i denote the Higgs doublets and Φ i ≡ i τ 2 Φ * i and Γ i , Ω i , Π i are matrices in flavour space. After spontaneous symmetry breaking, the following quark mass matrices are generated: where v i / √ 2 ≡ | < 0|φ 0 i |0 > | and α denotes the relative phase of the two vacuum expectation values (vevs) of the neutral components φ 0 i of Φ i . The neutral and the charged Higgs interactions with quarks are of the form: where v ≡ v 2 1 + v 2 2 ≈ 246 GeV and H 0 , R are orthogonal combinations of the fields ρ j , arising when one expands [1] the neutral scalar fields around their vevs, φ 0 j = e iα j √ 2 (v j + ρ j + iη j ), choosing H 0 in such a way that it has couplings to the quarks which are proportional to the mass matrices, as can be seen from Eq. (3); similarly, I denotes the linear combination of η j orthogonal to the neutral Goldstone boson. The matrices N 0 d and N 0 u are given by: The quark mass matrices are diagonalized through and the matrices N 0 d and N 0 u in the mass eigenstate basis transform into: There are similar expressions for the leptonic sector. We do not introduce neutrino masses since these are not relevant for our analysis.

Stability of the aligned 2HDM under RGE
The aligned two Higgs doublet model (A2HDM) is defined at tree level by the following relations involving the matrices introduced in Eq. (1): where d, u, e are constants. In this section we analyse the stability of the A2HDM under renormalisation group equations (RGE). The one loop renormalization group equations (RGE) for the Yukawa couplings are [18], [7]: where D ≡ 16π 2 µ (d/dµ) and µ is the renormalization scale. The coefficients a Γ , a Ω and a Π are given by: where g s , g and g are the gauge coupling constants of SU(3) c , SU(2) L and U(1) Y respectively. The alignment relations given by Eq. (7) guarantee the absence of Higgs mediated FCNC at tree level because both matrices M d and N d are proportional to Γ 1 . Similarly both M u and N u are proportional to Ω 1 and M l , N l to Π 1 . In general, these relations are broken at one loop level. From Eqs. (8), (9),(10) one can easily derive: In order to enforce Eq. (7) at one loop level it is easy to realize that it is sufficient to impose: in fact the proportionality constants on the r.h.s. are the running 6 of d, u and e. Therefore Eq. (16) does not impose any constraint: at one loop level the charged lepton sector remains aligned and there are no FCNC in the leptonic sector. This result agrees with the findings of references [19] and [20].
In equations (14), and (15) the pieces that can break the alignment in the quark sector are the terms: Ω 1 Ω † 1 Γ 1 and Γ 1 Γ † 1 Ω 1 respectively 7 . In order to have alignment at one loop level -fulfilling Eqs (17), (18)-there are two types of solutions: With λ Γ and λ Ω complex numbers.
Solutions of type 1 include the usual 2HDM with natural flavour conservation (NFC), where the up and down quarks receive contributions from only one Higgs doublet. It is well known that this can be achieved through the introduction of a Z 2 symmetry. Here, we are not interested in this class of well known solutions. We are interested in the class of solutions of type 2, and in the Appendix, we have studied the complete set of matrices Ω 1 and Γ 1 that obey to the conditions required for this class of solutions. We have shown in the Appendix, that if one requires stability under the RGE and at the same time Yukawa structures which are, in leading order, in agreement with the observed pattern of quark masses and mixing, then one is lead to a unique solution, where the matrices Ω 1 and Γ 1 are of the form: with ∆ the democratic mass matrix: This solution corresponds to the limit where only the top and bottom quarks acquire mass, while the two first generations are massless. The up and down quarks are aligned in flavour space, so the V CKM matrix equals the identity. The other stable solutions of type 2 correspond to non realistic cases like for example having all up or down quarks massless or two up or two down quark masses degenerate or with a V CKM very far from the identity matrix.
It is remarkable to realise that the so called democratic mass matrix is stable under RGE and that precisely this stability also enforces what could be called "quark alignment" in the sense that we also have a proportionality among Γ i and Ω i .

Natural Quasi-Alignment of Yukawa couplings
In this section we search for the minimal symmetry which when imposed on the Lagrangian, leads to the stable solution described in the previous section, corresponding to the democratic Yukawa couplings of Eq. (20). Before describing this symmetry, it is worth analysing another type of alignment which is verified experimentally, the so-called up-down alignment in the quark sector. 7 It can be readily seen that It is worth emphasizing that these structures are precisely the ones obtained in [19] and in [20], which produce FCNC at one loop level.

The up-down alignment in the quark sector
In the quark sector, flavour mixing is small. This means that there is a weak basis (WB) where both M u , and M d are close to the diagonal form. Experiment indicates that not only flavour mixing is small, but there is also up-down flavour alignment in the quark sector in the following sense. We can choose, without loss of generality, a WB where M u = diag(m u , m c , m t ). Of course, this is just a choice of ordering, with no physical meaning. Small mixing implies that in this WB M d is almost diagonal. In principle, since the Yukawa couplings Y u , Y d are not constrained in the SM, there is equal probability of M d being close to M d = diag(m d , m s , m b ). corresponding to up-down alignment, or being close for instance to M d = diag(m b , m s , m d ) in which case there is up-down misalignment. It is clear that in the SM, assuming small mixing and hierarchical quark masses, the probability of obtaining up-down alignment is only 1/6. Given a set of arbitrary quark mass matrices M u , M d , one can derive necessary and sufficient conditions to obtain small mixing and up-down alignment, expressed in terms of WB invariants [21]. Since the experimentally verified up-down alignment is not automatic in the SM, one may wonder whether there is a symmetry which leads to up-down alignment. In the next subsection, we propose a symmetry which leads to up-down alignment in the quark sector and when extended to a 2HDM leads to a natural alignment of the two Higgs doublets in flavour space.

Z 3 × Z 3 symmetry and the two Higgs alignment
We introduce the following Z 3 × Z 3 symmetry under which the quark left-handed doublets Q 0 L i , the righthanded up quarks u 0 L i , and the right-handed down quarks d 0 L i transform in the following way: where Z 3 corresponds to P = 1 + E 1 and Z 3 to P = 1 + E 2 with: and ω = e . This can easily be checked since ∆E 1 = ∆E 2 = 0. We thus conclude that the symmetry of Eqs. (22), (23) leads to the alignment of the two Yukawa coupling matrices, with a democratic flavour structure. Note that this solution also guarantees an up-down alignment in the quark sector, as defined in the previous subsection. In order to give mass to the first two quark generations, the Z 3 × Z 3 symmetry has to be broken. This breaking will also lead to Higgs mediated FCNC, but these couplings will be suppressed by the smallness of the quark masses. In order to illustrate how a realistic pattern of quark masses and mixing can be obtained, we shall assume that the breaking of the Z 3 × Z 3 symmetry occurs in two steps. In the first step the symmetry Z 3 × Z 3 is broken into just one of the Z 3 and the second generation acquires mass and finally in the last step the mass of the quarks u, d is generated. In the first step the symmetry Z 3 × Z 3 is broken to Z 3 generated by P = 1 + E 1 . One can check that: are invariant under this Z 3 symmetry. Note that A E 1 = 0. At this stage the second generation acquires mass. Finally, the lightest quarks, u and d acquire mass through a small perturbation, proportional toδ d,u , which breaks this Z 3 symmetry. We assume that: while with equivalent expressions for the up sector. Here where η is some complex number with modulus of order one. The symmetry is broken, and neither B u nor B d are invariant under the Z 3 × Z 3 symmetry.

Suppression of scalar mediated FCNC
In order to study the suppression of scalar mediated FCNC, it is useful to start by analysing the parameter space in our framework.

The Parameter Space
From Eqs. (2), (25) and (26) it follows that, in leading order Writing we obtain in leading order the following relations: which impose restrictions on the allowed parameter space. A priori, we do not assume any conspiracy between parameters and take t = O(1). It is then clear from Eq. (30), that the c u i are generically of order one, while c d i are smaller, and may assume values of order O( m b m t ). This is an important ingredient which, as we shall see, will play a rôle in the evaluation of the strength of FCNC's and the allowed parameter space for the Higgs masses.
Next we give the structure of the flavour changing neutral Yukawa couplings. For that, it is useful to express the quark mass matrices in Eq. (2) in terms of the perturbations given in Eqs. (25), (26): Then, we derive the expressions for the matrices which couple to the Higgs scalars in Eqs. (3), (4). In the basis where the up and down quark matrices are diagonal, the matrices N d , N u of Eq. ( 6) become where we have used Eqs. (4) and (6) with Eq. (2). Finally, from Eq. (32) combined with Eqs. (25) and Eq.
(31) we find The crucial point is that in our scheme these matrices have an extra suppression factor, proportional to δ d,u .
Using the expressions given in Eqs. (31), (27) and computing the trace, second invariant and determinant for the squared quark mass matrices H u,d ≡ MM † u,d , one can find that in leading order: where λ ≡ 0.2 is of the order of the Cabibbo angle.

From Eq. (31) it follows that in leading order
is the matrix that diagonalizes the exact democratic limit ∆. Thus, taking into account Eq. (27), the matrix contributions from U † d L B d U d R and U † u L B u U u R are both of order one. One can thus conclude that: -for the down sector, with the assumptions made after Eq. (30), we have a total suppression factor of O( m b m t ) · O λ 3 -for the up sector, we have a suppression factor of O λ 5 or smaller depending on the value that we choose to assume for c u 1 , but which, as explained, it is reasonable to take of order one.

Numerical analysis
The matrices of Eq. (31) may be explicitly written as where we have introduced . Although these two coefficients are in general complex, and since the physically meaningful matrices are those defined as H = M M † , both coefficients may be taken as real for our numerical exercise. If one then parametrizes the remaining variables as one is left with twelve real parameters that compose the quark mass matrices in our scheme.
In order to check if this parameter space could accommodate the flavour sector, a numerical survey was made where we looked for one combination that could fit the observed values of the quark masses given at the scale of the Z boson mass [22], It should be noted that the twelve parameters fix not only V CKM and the quark mass spectrum, but also the strength of all the FCNC couplings. In order to evaluate the numerical stability of this reference point, we performed a numerical check, varying the input parameters randomly around the values that produced the reference point above; the new results were then combined in the scatter plots shown in Fig. 1 where the reference point is highlighted. In order to obtain an estimate of the lower bound for the flavour-violating Higgs masses, we consider the contribution to K 0 − K 0 mixing. Apart from the SM box diagram one now has a New Physics contribution arising from the scalar-mediated FCNC tree-level diagrams thus making the total transition amplitude equal to M 12 = M SM 12 + M NP 12 , [14] with: where: and The indices q 1 and q 2 refer to the valence quarks of the meson M, and N is N u or N d depending on the meson system considered. In this framework it is a good approximation to use the matrix F for both U d L and U d R . Using the values we obtained for δ d and taking, as already discussed, t 1 and c d < ∆m K . The region plot that we have obtained is presented in Fig.2. It is clear that in this framework, the masses of the flavour-violating neutral Higgs can be below the TeV scale, so that they could be discovered at the next run of the LHC.

Conclusions
We have studied in detail, in the framework of 2HDM, the question of stability of alignment, under the renormalization group. It was shown that there are new stable solutions, apart from those found in Ref. [7]. Stability under the RGE puts very strict restrictions on the flavour structure of the Higgs Yukawa couplings. If one imposes the stability conditions and at the same time requires that the flavour structure is in agreement with the observed pattern of quark masses and mixing, then one is lead to a unique solution, where all Higgs flavour matrices are proportional to the so-called democratic matrix. We have also shown that these flavour structures leading to stable alignment can be obtained by imposing on the Lagrangian a Z 3 × Z 3 symmetry. In the limit where this symmetry is exact, only the third generation of quarks acquires a mass. Nonvanishing masses for the two first generations are obtained through the breaking of the discrete symmetry which in turn generates scalar mediated FCNC which are suppressed by the smallness of the light quark masses. The scenario presented in this paper establishes a possible intriguing link between stability of alignment in 2HDM and the observed pattern of quark masses and mixing. and we conclude that implying that V CKM = I up to permutations of rows or columns. Denoting the usual bi-unitary diagonalisation procedure by from Eq. (49) we conclude that we can always choose the unitary matrices V Γ L and V Ω L equal to each other and the alignment conditions can be easily reduced to conditions among the diagonal matrices D Γ and D Ω . From Eq. (46), it then follows that: It can be checked that there are only two types of solutions. Those with λ Γ and λ Ω different from zero (solutions 1,2 and 3) and the remaining ones (solutions 4 and 5 ) 1. D Γ = aP 3 and D Ω = αP 3 and changes of P 3 by P 2 or P 1 2. D Γ = a (I − P 1 ) and D Ω = α (I − P 1 ) and changes of P 1 by P 2 or P 3 3. D Γ = aI and D Ω = αI 4. D Γ = 0 and D Ω arbitrary and viceversa.
5. D Γ = aP i and D Ω = α (I − P i ) where P i stand for the projection operators Solutions 2,3 and 4 cannot be good approximations to the actual quark spectra due to the implied degeneracy. Solution 5 gives rise to V CKM matrix very different from the identity matrix. Only solution 1 provides, in leading approximation the correct pattern of quark masses and mixing. In a suitable weak-basis, this solution can be written as a democratic matrix ∆.