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 |Vus|≃|Vcd|\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|V^{}_{us}| \simeq |V^{}_{cd}|$$\end{document}, |Vcb|≃|Vts|\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|V^{}_{cb}| \simeq |V^{}_{ts}|$$\end{document} and |Vub|/|Vcb|<|Vtd|/|Vts|\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|V^{}_{ub}|/|V^{}_{cb}| < |V^{}_{td}|/|V^{}_{ts}|$$\end{document}, 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σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\sigma $$\end{document} level and call upon a further investigation of the realization of these structures from a top–down prospective.


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 V [7,8], and the moduli of its nine matrix elements |V i j | 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 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 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 established [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 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 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 a 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 cannot 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 Sect. 2, 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). Section 3 is devoted to an approximate and analytical approach to diagonalize hierarchical quark mass matrices. The expressions of CKM matrix elements have been derived. In Sect. 4, 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 Sect. 5.

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 righthanded 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 Moreover, the mismatch between those diagonalizations results in the CKM matrix V = V † u V d , which enters into the flavor-changing chargedcurrent interactions and thus can be precisely measured in experiments.
It is evident that a weak-basis transformation, defined as unitary matrix, keeps the CKM matrix V unchanged. 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 (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 onto 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.

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 ith 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 have 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 parametriza- is the rotation matrix in the i j-plane with an angle θ q i j (for i j = 12, 13, 23), we have O T q M q O q =M q , whereM 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 i j obtained in the following discussion.
Due to the hierarchical structure of M q , we expect that the rotation angles θ q i j 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 i j ≡ sin θ q i j and c q i j ≡ cos θ q i j 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 i j 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 Note that the rotation angles have been expressed in terms of quark masses and four free parameters {ε u , d u } and {ε d , d d }. 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 q = 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, which 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) [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).

Numerical analysis
Now we proceed to explore the parameter space of the nonparallel four-zero textures of quark mass matrices in Eqs. (4) and (5) 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 V ud V * ub + V cd V * cb + V td V * tb = 0, e.g.
etc., which can be determined independently from the measurements of CP violation in B decays [9].

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  where the off-diagonal elements of the real and symmetric matrices M u and M d have been omitted. A numerical investigation 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 • has been extracted 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 we observe agreement with hierarchical structures in Eq. (8), e.g.
and the following allowed ranges for the two phases are obtained:

Case II
Now we turn to the second case in Eq. (5) with d u = 0, ε u = 0, d d = 0 and ε d = 0. Hence, there are also four free parameters {d u , ε d } and {φ 1 , φ 2 }. As in the 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: 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

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 non-parallel 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 a possible extension of the Fritzsch ansatz has been carried out under the assumption of factorizable phases.
We are able to extract vital clues toward 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 the 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.