Quark Flavor Mixings from Hierarchical Mass Matrices

In this paper, we extend the Fritzsch ansatz of quark mass matrices while retaining their hierarchical structures and show that the main features of the Cabibbo-Kobayashi-Maskawa (CKM) matrix $V$, including $|V^{}_{us}| \simeq |V^{}_{cd}|$, $|V^{}_{cb}| \simeq |V^{}_{ts}|$ and $|V^{}_{ub}|/|V^{}_{cb}|<|V^{}_{td}|/|V^{}_{ts}|$, can be well understood. This agreement is observed especially when the mass matrices have non-vanishing $(1,3)$ and $(3,1)$ off-diagonal elements. The phenomenological consequences of these for the allowed texture content and gross structural features of `hierarchical' quark mass matrices are addressed from a model independent prospective under the assumption of factorizable phases in these. The approximate and analytical expressions of the CKM matrix elements are derived, and a detailed analysis reveals that such structures are in good agreement with the observed quark flavor mixing angles and the CP-violating phase at the $1\sigma$ level and call upon a further investigation of the realization of these structures from a top-down prospective.


I. INTRODUCTION
The ATLAS [1] and CMS [2] experiments at the CERN Large Hadron Collider (LHC) have recently discovered a scalar particle of mass m h ≈ 125 GeV, which is likely to be the Higgs boson in the standard model (SM), witnessing a great success of the SM in describing the electroweak interactions. However, the flavor structures of Yukawa interactions between the Higgs boson and the fermions are poorly understood. For instance, the strong hierarchy among quark masses, the smallness of quark flavor mixing angles and the origin of CP violation remain to be explained. As the quark masses and flavor mixing parameters have been precisely measured, it is an important and urgent task to explore the flavor dynamics of quarks [3][4][5].
To be explicit, we evaluate the running quark masses at the energy scale of M Z = 91.2 GeV. The mass ratios are found to be [6]  On the other hand, the quark flavor mixing is described by the CKM matrix * Electronic address: rohitverma@ihep.ac.cn † Electronic address: zhoush@ihep.ac.cn V [7,8], and the moduli of its nine matrix elements |V ij | for i = u, c, t and j = d, s, b have been determined from current experimental data with high precisions [9]. A strongly hierarchical structure of CKM matrix elements is also observed, namely, |V tb | > |V ud | > |V cs | |V us | > |V cd | |V cb | > |V ts | |V td | > |V ub |. These strong hierarchies in quark mass spectra and flavor mixings are probably correlated with each other.
Unfortunately, we have not yet had a convincing flavor theory to realize such a correlation. Therefore, phenomenological approaches on specific textures of quark mass matrices have been widely adopted to obtain vital clues on the underlying flavor dynamics [10][11][12][13][14][15][16][17]. However, experience shows that there exist a vast number of texture zero structures, [18][19][20][21][22][23][24][25][26][27][28][29][30][31] both for the quark and lepton sectors, which in principle, may not be ruled out by experiments. This suggests that the texture approach alone is perhaps not sufficient enough to account for the flavor dynamics and one may require to use additional assumptions (based on observations) to extract some useful information on the gross possible structures of these matrices from a model-independent prospective, particularly, for the quark sector.
A natural explanation for quark flavor dynamics requires a translation of the observed hierarchies in the quark masses and mixing angles onto the corresponding matrix elements in these mass matrices. One typical example is the Fritzsch mass matrices with six texture zeros [32,33] where q = u or d denotes the up-type or down-type quark mass matrix, and a pair of off-diagonal texture zeros in the Hermitian or symmetric matrix are usually counted as one zero. As a consequence of texture zeros [32][33][34][35][36][37][38][39][40][41][42], a few simple but instructive relations among flavor mixing angles and quark mass ratios can be estab-lished [34,36,[43][44][45][46][47][48][49] and tested in experiments. Unfortunately, the Fritzsch ansatz failed to reconcile a large top-quark mass with the observed value of |V cb | [50,51]. Possible extensions of the Fritzsch ansatz have been considered in the literature [3]. For instance, a nonzero element in the (2, 2) position of M q , i.e., has been proved to be important to accommodate the experimental data on quark flavor mixings [52]. A remarkable difference between the Fritzsch ansatz in Eq.
(2) and the parallel four-zero texture in Eq. (3) is that a hierarchical structure is not allowed in the latter case [39,40,53,54]. Even more, it has recently been found in Ref. [42] that the relation |D q | |B q | |C q | is valid. By contrast, the Fritzsch ansatz with |A q | |B q | |C q | seems to reflect the strong hierarchies of quark mass spectra and flavor mixings in a more transparent and natural way. However, one has not been able to find [55] any strongly hierarchical structure for quark mass matrices which explains the observed flavor mixings.
The purpose of this paper is to extend the Fritzsch mass matrices of Eq. (2) assuming hierarchical structures for the corresponding mass matrices following a natural mixing among the three quark generations, see Eq.(8). However, a texture zero at the diagonal (3, 3) position disturbs such a natural hierarchy among the mass matrix elements and is therefore ruled out in this context. Particularly, we find that the non-vanishing matrix elements at the off-diagonal (1, 3) and (3, 1) positions are crucial in retaining a natural mass matrix hierarchy. For example, one may consider nonzero elements in both (1,3) and (2, 2) positions of M u , while M d takes the form of the Fritzsch ansatz: or vice versa. Another possibility can be obtained from the first type by just moving the nonzero element in the (1, 3) position of M u to the same place of M d , i.e., or vice versa. More examples of such such viable texture structures were recently [24] obtained for symmetric mass matrices. However, the study was mainly limited to a numerical investigation only while a desirable phenomenological explanation of these was not explored and requires an urgent investigation. We undertake a phenomenological study of such texture structures, when M q (3, 3) = 0, and find that these have a common unique feature of displaying strong hierarchy among the mass matrix elements and that these features can still be retained for Hermitian mass matrices within the framework of the Standard Model. It is worth stressing that the matrix elements D q of M q can not be vanishing simultaneously, for texture specific hierarchical mass matrices, due to the constraint from the precise measurement of |V cb |, as we shall show later. In addition, the nonzero element F u or F d is indispensable to fit the experimental value of |V ub | while preserving the hierarchical structures of quark mass matrices. In Sec. II, based on the weak-basis transformation, we argue that the non-parallel four-zero textures in Eqs. (4) and (5) are as general as the parallel four-zero texture in Eq. (3). Sec. III is devoted to an approximate and analytical approach to diagonalize hierarchical quark mass matrices. The expressions of CKM matrix elements have been derived. In Sec. IV, we explore the allowed regions of model parameters by performing a complete numerical analysis and find an excellent agreement with the experimental data. Finally, we summarize in Sec. V.

II. WEAK BASIS TRANSFORMATIONS
In the SM, the Yukawa interactions between the Higgs boson and the fermions play a very important role. First, for a light Higgs boson of mass m h ≈ 125 GeV, the dominant decay channel is h → bb and the decay width is governed by the Yukawa coupling of bottom quark. Second, the quark mass spectra, flavor mixing angles and CP violation are determined by the complex Yukawa coupling matrices. As the right-handed quark fields are gauge singlets, it is always possible to make the quark mass matrices M u and M d Hermitian through a unitary transformation of right-handed quark fields in the flavor space, if there is no flavor-changing right-handed current as in the SM. The diagonalizations of M u and M d lead to quark mass eigenvalues

which enters into the flavor-changing charged-current interactions and thus can be precisely measured in experiments.
It is evident that a weak-basis transformation, defined However, it indeed changes the flavor structure of quark mass matrices. Through the weak-basis transformations, one can prove that the quark mass matrices can be recast into the following forms with three texture zeros: as demonstrated in Ref. [56,57] or as pointed out in Ref. [58]. Notice that Eqs. (6) and (7) are valid as well when the texture structures of M u and M d are exchanged. Note also that the crosses in the mass matrices in Eqs. (6) and (7) stand for nonzero elements, and the Hermiticity is implied. The texture zeros stemming from weak-basis transformations carry no physical information, but additional ones lead to correlations among quark masses and flavor mixing angles. Two comments are in order: • In Eq. (7), the element in the (1, 1) position of M d has been found to be redundant [59][60][61] in the sense that the quark mass matrices are compatible with experimental data even in the limit of (M d ) 11 = 0. In this limit, we obtain the Fritzsch-like mass matrices with four texture zeros in Eq. (3). Therefore, the parallel four-zero textures can be derived from the most general ones in Eq. (7) under the minimal additional assumption of a vanishing (1, 1) diagonal element. Furthermore, most of the recent investigations involving texture specific quark mass matrices have been limited to parallel texture zeros [55,60,61] or based on numerical investigations only [24].
• A comparison of mass matrices in Eqs. (4), (5) and Eq. (6) reveals that one matrix element in the diagonal position of M u or M d in Eq. (6) has been taken to be zero. Hence, the non-parallel four-zero textures are as general as the parallel ones in Eq. (3) and deserve a dedicated study.
In the following section, we validate an approximate and analytical diagonalization of such non-parallel quark mass matrices and consequently establish simple relations between the CKM matrix elements and quark masses. With the help of these relations, we can translate the observed strong hierarchies among the quark masses and mixing angles into the hierarchical structures of quark mass matrices. Such a translation is not possible for the parallel four-zero textures in Eq. (3) where |V ub | is obtained from a next-to-leading order mixing of the neighboring quark generations owing to a suppression of a leading-order mixing among the first and third quark generations through the position of texture zeros in these.

III. ANALYTICAL RESULTS
In view of the strong hierarchy of quark masses, we expect naturally the hierarchical structures of quark mass matrices, which in the most general case can be formulated as follows where m i denotes the mass of i-th generation up-type (or down-type) quark for q = u (or d), and each matrix element is expressed in terms of its modulus and argument. In this formulation, the hierarchical structure of a mass matrix should be understood as the rightmost one in Eq. (8) up to a factor of O(1) for each nonzero matrix element. The hierarchical structures have been extensively discussed in the literature [3][4][5], for example, we In order to achieve an analytical diagonalization of the mass matrices and realize simple empirical relations among the quark mixing angles and the corresponding mass ratios, we assume that the phases are factorizable in M q i.e. γ q = α q + β q . This also reduces the number of free parameters in the associated quark mass matrices, to be compared with ten physical observables. As a result, one can rewrite the mass matrix as M q = P q M q P † q , where P q ≡ diag{e iα q , 1, e −iβ q } and M q is a real and symmetric matrix. The diagonalization of M q can now be achieved by an orthogonal matrix O q , which can be parametrized in terms of three rotation angles {θ q 12 , θ q 13 , θ q 23 }. Adopting a convenient parametrization O q ≡ R 13 (θ q 13 )R 12 (θ q 12 )R 23 (θ q 23 ), where R ij (θ q ij ) is the rotation matrix in the ij-plane with an angle θ q ij (for ij = 12, 13, 23), we have O T q M q O q =M q , wherẽ M q = diag{k q m 1 , −k q m 2 , m 3 } is a diagonal matrix with quark mass eigenvalues and k q = ±1. The diagonalization of the mass matrices using O T q M q O q =M q for the above hierarchical mass matrices, with factorizable phases, requires one negative mass eigenvalue to ensure real values for trigonometric sines of the associated rotation angles θ q ij obtained in the following discussion. Due to the hierarchical structure of M q , we expect that the rotation angles θ q ij are associated with the quark mass ratios and thus in general very small. To the leading order of small rotation angles, the orthogonal matrix is given by where s q ij ≡ sin θ q ij and c q ij ≡ cos θ q ij have been defined.
Then, the CKM matrix is determined by where φ 1 ≡ α u − α d and φ 2 ≡ β u − β d . The final task is to calculate the rotation angles θ q ij from the quark mass matrices. and in the leading-order approximation. Taking e q = 0, as for the most general mass matrices in Eq. (6), one obtains They can be substituted into Eq. (11) to obtain the other matrix elements ofM q . According to Eqs. (9) and (10), it is easy to derive the following simple relations for the off-diagonal elements of the CKM matrix: Some remarks on the above approximate and analytical results are helpful. First of all, the relation |V us | ≈ |V cd | ≈ m d /m s ≈ λ holds as in the case of Fritzsch mass matrices. Second, if d u = 0 is assumed, we can derive a lower bound on |V cb | from Eq. (14), namely, where the quark mass ratios in Eq. (1) have been used. The lower bound in Eq. (15) is in contradiction with the experimental value |V cb | = 0.0414 ± 0.0012 that has been precisely measured in the exclusive and inclusive semileptonic decays of B mesons into charm [9]. Hence, it is evident why the Fritzsch ansatz with a hierarchical structure and d u = d d = 0 is excluded. A nonzero d q is of crucial importance to accommodate the experimental data. Third, we observe from the last two lines in Eq. (14) that where the quark mass ratios are used and ε u = ε d = 0 is assumed. The second relation in Eq. (16) is compatible with |V td |/|V ts | = 0.216±0.011, which has been extracted from the mass difference for the neutral B mesons in the B 0 d -B 0 d system and that in the B 0 s -B 0 s system [9]. But, the estimated value of |V ub |/|V cb | ≈ 0.0469 in Eq. (16) is smaller than the observed one 0.0857 [9] almost by a factor of two.
Based on the above observations, we conclude that ε u and ε d cannot be vanishing simultaneously, and likewise for d u and d d . To emphasize this further, we carry out a detailed investigation of mass matrix structures in Eqs. (4) and (5) as the minimal extension of the Fritzsch ansatz where hierarchical structure is retained. The analytical results in Eq. (14) not only hold for the general quark mass matrices in Eqs. (6), (8) but are also applicable to all viable Hermitian texture structures obtained in Ref. [24] with M q (3, 3) = 0 under the assumption of factorizable phases in these.
Furthermore, from the expression for s q 23 in Eq. (13), it is evident that d q = 0 is not allowed for k q = −1 corresponding to the negative m 1 eigenvalue. As a result, one is compelled to consider m 2 negative for the corresponding mass matrices with d d = 0 in Eq. (4) and d u = 0 in Eq. (5). However, a negative sign associated with either of the two mass eigenvalues is consistent with the data for mass matrices with d q = 0 and it is trivial to check that the allowed parameter space for d q is not very different in such cases. In particular, the modular d q values are of the same order but marginally larger for k q = −1 as compared to the k q = +1 which is also clearly evident through the expression for V cb ≈ s d 23 − s u 23 e iφ2 when d d = 0 (or d u = 0). As a result, we restrict our analysis with k q = +1 for obtaining greater agreement with the notion of 'hierarchical' mass matrices assumed in Eq. (8).

IV. NUMERICAL ANALYSIS
Now we proceed to explore the parameter space of the non-parallel four-zero textures of quark mass matrices in Eqs. (4) and (5). We use the running quark masses at M Z = 91.2 GeV as inputs where the central values with 1σ errors are taken from Ref. [6]. The mass matrices are claimed to be consistent with experimental data if the absolute values of three CKM matrix elements and one inner angle of the unitarity triangle can be reproduced [9], i.e., where the 1σ errors are attached to the best-fit values. At the same time, the allowed regions of free model parameters can be obtained. As usual, three inner angles of the unitarity triangle correspond to the orthogonal condition etc., which can be determined independently from the measurements of CP violation in B decays [9].

A. Case I
We consider the quark mass matrices in Eq. (4). In this case, we have f d ∝ ε d = 0 and d d = 0, but f u ∝ ε u = 0 along with d u = 0. Hence, the free model parameters include {d u , ε u } and two phases {φ 1 , φ 2 }. In our numerical calculations, we exactly diagonalize the quark mass matrices, and freely vary the physical ranges of free model parameters {d u , ε u } in their allowed parameter space consistent with Eqs. (8) and (14). Both φ 1 and φ 2 are allowed to freely vary in [0 • , 360 • ] interval. Given quark masses and the values of four model parameters, we are able to calculate the CKM matrix elements. Only when the experimental constraints in Eq. (18) are satisfied, we output the values of model parameters and reconstruct the quark mass matrices.
The allowed regions of f u and d u are given in Fig. 1, where the 1σ ranges of |V ub | and |V cb | are also shown as vestigation of the above mass matrices, without using the condition of factorizable phases, has been more recently carried out in [62,63] pointing to possible hierarchical structures for these matrices. In Fig. 2, the predictions for the other two inner angles α and γ are given, since the free parameters ε u , d u , φ 1 and φ 2 are restrictively constrained by the experimental data. Currently, the angle α = 85.4 •+3.9 • −3.8 • has been extracted from the CP violation in the B → ππ, ρπ and ρρ modes, while γ = 68 •+8.0 • −8.5 • from that in the B ± → DK ± and B 0 → D * ± π ∓ decay modes. Thus, the precise determinations of α and γ in B factories and future collider experiments (e.g., LHCb) will test the unitarity of the CKM matrix. Additionally, we also find the Jarlskog invariant [64][65][66] J CP ≡ Im [V ts V * tb V * us V ub ] = (2.8 · · · 3.4) × 10 −5 , which is in a perfect agreement with the data J CP = 3.06 +0. 21 −0.20 × 10 −5 [9]. We also reconstruct below the viable quark mass matrices corresponding to an interchange of textures for the above case, and observe an agreement with hierarchical structures in Eq. (8), e.g.
M u = 0 0.02 · · · 0.03 0 0 9.8 · · · 10.8 171 · · · 173 GeV , 0.01 · · · 0.02 0.003 · · · 0.01 −0.07 · · · − 0.04 0.07 · · · 0.21 2.8 · · · 3.0 GeV , and the following allowed ranges for the two phases are obtained i.e. previous case, these model parameters are allowed free variations within the corresponding allowed parameter spaces governed by Eqs. (8) and (14). With the help of experimental constraints in Eq. (18), we obtain the following quark mass matrices that are consistent with the data at the 1σ level: GeV , M d = 0 0.010 · · · 0.014 0.01 · · · 0.02 0 0.3 · · · 0.4 2.8 · · · 3.0 GeV . (22) In addition, the allowed regions of two phases turn out to be φ 1 ∈ [250 • , 300 • ] and φ 2 ∈ [13 • , 18 • ]. In order to understand the dependence of the CKM matrix elements on the nonzero elements d u and f d , we show the values of (|V ub |, f d ) and (|V cb |, d u ) in the upper and lower panels of Fig. 3, respectively. Here we emphasize the importance of a non-vanishing f d for |V ub | and observe that f d ∼ √ m d m s ∼ a d accounts for the observed |V ub | within the 1σ interval. This is understandable from Eq. (14), since the down-type quark masses are less hierarchical than the corresponding masses of the up-type quarks. The allowed regions of the three quark flavor mixing angles and the three inner angles of the unitarity triangle remain the same as in the previous case. The Jarlskog invariant J CP = (2.82 · · · 3.44) × 10 −5 is also well consistent with the experimental observations. For the purpose of completion, we present below the viable structures for the quark mass matrices when the texture zeros are interchanged in the above case, e.g.

V. CONCLUSIONS
In the present paper, based on the weak-basis transformations [56], we argue that the non-parallel four-zero textures of quark mass matrices are as general as the parallel ones, which have attracted a lot of attention in recent times [55]. Moreover, one advantage of the nonparallel textures is that these offer a natural translation of the observed strong hierarchies in the quark masses and mixing angles onto the corresponding mass matrix elements. A phenomenological analysis of generic quark mass matrices with non-parallel textures incorporating possible extension of Fritzsch ansatz has been carried out under the assumption of factorizable phases. We are able to extract vital clues towards the formulation of these matrices using the available data on quark masses and CKM parameters in a model-independent way. Two interesting observations are made. First, the experimental data from |V ub | measurements requires a nonzero element in the (1, 3) position of quark mass matrices, namely, f u ∝ √ m u m t or f d ∝ √ m d m s ∼ a d , which is also crucial in realizing such hierarchical structures. Second, the experimental constraints from |V cb | measurements forbid d u = d d = 0 in these quark mass matrices. We also derive general empirical relations for the CKM matrix elements, which offer a natural description of |V us | |V cd |, |V cb | |V ts | and |V ub |/|V cb | < |V td |/|V ts |.
The study indicates that there should exist, not one [55], but several viable sets of hierarchical and Hermitian quark mass matrices, which may provide further possible clues for model building in top-down approaches [67][68][69][70][71][72][73][74]. In this context, it has been recently shown [42] that the texture zeros in the quark sector are essentially stable against the one-loop renormalization-group-evolution of energy scale required in the top-down investigations of such mass matrices. It should be interesting to see if those specific structures of quark mass matrices can be verified in the future B-physics experiments.