Flavored axion in the UV-complete Froggatt-Nielsen models

We propose UV-completions of Froggatt-Nielsen-Peccei-Quinn models of fermion masses and mixings with flavored axions, by incorporating heavy fields. Here, the $U(1)$ Froggatt-Nielsen symmetry is identified with the Peccei-Quinn symmetry to solve the strong CP problem along with the mass hierarchies of the Standard Model fermions. We take into account leading order contributions to the fermion mass matrices giving rise to Nearest-Neighbour-Interaction structure in the quark sector and $A_2$ texture in the neutrino sector. A comprehensive numerical analysis has been performed for the fermion mass matrices. Subsequently, we investigate the resulting axion flavor violating couplings and the axion-photon coupling arising from the model.


INTRODUCTION
The axion [1][2][3] is one of the best motivated elementary particles beyond the Standard Model (BSM). It provides the most cogent solution to the "strong CP problem" 1 of quantum chromodynamics (QCD), where the spontaneous breaking of a global U (1) symmetry, called the Peccei-Quinn (PQ) symmetry U (1) P Q , gives rise to a pseudo-Nambu-Goldstone (pNG) boson [2,3]. It is also a good candidate for cold dark matter within a certain allowed parameter space of the axion scale [6][7][8][9][10]. In the literature, there exist two common approaches to realize the PQ symmetry that give rise to the axion, in one case the Higgs sector is extended (called the DFSZ model) [11,12], whereas in a different scenario heavy quarks are introduced (called the KSVZ model) [13,14]. As per the experimental testability of the axion is concerned, it is the axion couplings to the photon g aγ and the electron g ae that play very important roles. One can establish a relation between these couplings with the PQ charges of fermions [15]. It is the constant E/N , i.e., the ratio of the electromagnetic over the QCD anomaly coefficient which determines the strength of the photon-axion coupling, and it is a free parameter in generic axion models.
Another puzzle in the Standard Model (SM) is the unexplained hierarchies in masses and mixings between elementary particles or the "hierarchical flavor structure" of the Yukawa couplings.
Many flavor models have been studied to address the SM mass hierarchies involving additional symmetries. Among them, the Froggatt-Nielsen (FN) symmetry emerges as a leading candidate to account for the flavor hierarchies [16]. In this mechanism, one introduces a new complex scalar field called flavon, whose vacuum expectation value (vev) generates the hierarchical flavor structure, where a global flavor symmetry U (1) F N is imposed. Now, here one can ask, if it is possible to use the PQ symmetry as a flavor symmetry in order to explain the hierarchical flavor structure while solving the strong CP problem of the SM. Indeed, in [17][18][19] a connection between the PQ symmetry with flavor symmetries has been addressed. Recently, many attempts have been made where the PQ symmetry U (1) P Q is unified with the FN symmetry U (1) F N to address both issues of the SM. The resultant axion that arises from this flavor-dependent framework has been called "flaxion" [20,21], "axiflavon" [22][23][24], or "flavorful axion" [25][26][27].
In this work, we study various low energy phenomenologies of the flavor-dependent axion that arises from the PQ symmetry providing a solution to the strong CP problem. Here, only the third family of quarks is generated at dimension four while other terms are introduced at higher dimensions. We also extended the flavor structure to the leptonic sector. In doing so, we identify the U (1) P Q symmetry with the U (1) F N that gives rise to the fermion hierarchies, where the axion field σ is treated like a flavon. We consider the DFSZ type axion models and show that one can construct a Nearest-Neighbour-Interaction (NNI) structure for both the up-and down-quark mass matrices which was originally proposed in [28]. For the DFSZ models, we affix the PQ charges to all the SM fermions as well as to the Higgs sector in such a way that after the PQ and electroweak (EW) symmetry breaking, one ends up with the NNI structure of the quark mass matrices as given by Eq. 1. Later, we present the UV-complete Froggatt-Nielsen-Peccei-Quinn (FNPQ) models by incorporating heavy fermionic fields 2 . Two examples have been adopted within the framework of the type-I and -II Dirac seesaw mechanisms. In both the scenarios, it is shown that once we integrate out heavy fields, one finds NNI-type quark mass matrices.
Also, to explain the masses and mixings of the leptonic sector, we find the diagonal chargedlepton mass matrix in the FNPQ-framework. Moreover, within the framework of the type-I seesaw mechanism [29][30][31][32][33], we obtain an A 2 -type [34] neutrino mass matrix which is allowed by the latest global analysis of the neutrino oscillation data [35][36][37]. We also perform a numerical analysis of the quarks as well as leptonic mass matrices to find masses and mixings of fermions. Later, we explore flavor violating decays with axions, where we estimate the branching ratios (BR) of K + → π + a and B + → K + a. These further help us to constrain the axion breaking scale. Furthermore, because of the presence of two Higgs doublets in the DFSZ axion model, we calculate the flavor violating process t → hc. Finally, the cosmological consequences of the FNPQ axion are discussed.
We plan the manuscript as below. In Sec. 2 we give a brief overview of the FN symmetry along with the NNI quark mass matrices. A detailed discussion of our theoretical framework that leads to NNI quark textures together with the UV-complete mechanisms are presented in Sec. 3 and its subsequent sub-sections, while the leptonic sector is described in Sec. 3 3.3. The scalar sector and axion couplings are presented in Sec. 3 3.4, and Sec. 3 3.5, respectively. We analyze our numerical results in Sec. 4, and conclusions are summarized in Sec. 5.

U (1) F N AND NEAREST-NEIGHBOUR-INTERACTION
Here, we start with a brief overview of the Froggatt-Nielsen symmetry and Nearest-Neighbour-Interaction structure for the fermion mass matrices. In order to address the hierarchical structure of the fermion masses Froggatt and Nielsen developed a symmetry called FN-symmetry [16], where they introduced a SM singlet "flavon" field σ. In their framework, one can express the necessary Yukawa terms for the quark sector as where n u/d ij are complex numbers, σ is the flavon field and Λ is the scale of flavor dynamics. Also, Q i , u Ri , d Ri for i, j = 1, 2, 3 represent the left-handed quark doublet, right-handed up-type quark, and right-handed down-type quark, respectively. Notice that one can write similar terms for the leptonic sector. Once, the singlet acquires its vev and breaks the symmetry, it leads to hierarchy in the masses of quarks and leptons. Thus, one can notice from Eq. 2 that the hierarchy of the Yukawa couplings can be correlated to the smallness of σ /Λ, and hence to the fermion mass matrices m ij = y ij v( σ /Λ) n , with v being the SM Higgs vacuum expectation value.
Furthermore, this global U (1) symmetry can be used to obtain NNI-textures in the quark sector. The NNI-texture for a N × N mass matrix consists of a general matrix m with m ij = 0 only for i = j ± 1 or i = N = j. Notice that one can connect the NNI-texture with the Fritzsch texture [38,39] with an additional constraint of Hermiticity on both the quark mass matrices.
However, the Fritzsch texture was ruled out as it failed to simultaneously generate the small value of the CKM matrix element V cb together with the large value of top-quark mass m t [40,41]. The NNI texture is phenomenologically compatible with the quark sector, but contains more degrees of freedom in contrast to the experimental constraints, i.e. 6 masses, 3 angles and a CP-violating phase of the quark sector. In the SM, the NNI texture may be obtained through Weak-Basis (WB) transformations from an arbitrary Yukawa matrix [28], and thus imposes no physical constraints.
However, in models with more than one scalar a global symmetry that imposes the NNI on the quark mass matrix imposes a particular structure of scalar interactions [42], potentially leading to physical observables. It is well known that the NNI structure can reproduce the CKM mixing pattern, while demanding a hierarchy between its elements to obtain that of the quark masses [43].
One particular feature is that the (3, 3) element is of order m t (m b ) in the up (down) sector, while the other entries are heavily suppressed. The proposed framework for the flavored PQ symmetry accommodates the DFSZ style of axion models, while leading to the NNI structure of the quark mass matrices, explaining naturally the largeness of the (3, 3) element, and the A 2 structure of the neutrino mass matrix in the diagonal charged lepton basis.

FERMION MASSES IN FLAVORED DFSZ AXION MODELS
Now, we propose UV-complete models of quark masses, where a U (1) global symmetry enforces NNI texture in the up and down quark mass matrices. The (3, 3) entry of the matrix is generated by a dimension-4 operator, while the other entries are obtained by dimension-5 operators with the addition of a symmetry-breaking flavon, thus obtaining mass hierarchies in a Froggatt-Nielsen style. Additionally, through the color anomalies of the global U (1) symmetry we obtain a QCD axion from the flavon fields, which solves the strong CP problem and is a Dark Matter candidate.
We present two realizations within the framework of the DFSZ style of QCD axions.
In the DFSZ axion model, two Higgs doublets are needed in addition to the scalar singlet which introduces the PQ symmetry breaking scale. We list the relevant field contents with their corresponding SU (2) L × U (1) Y and U (1) P Q charges in Table I. We have derived the PQ charges by imposing the appearance of the NNI structure at dimension-5, while choosing opposite charges for left and right handed fields in the spirit of simplicity. Additionally we have followed the convention followed in DFSZ models where the charges of the Higgs doublets are normalized to unity.
U (1) P Q (9/2, -5/2, 1/2) (-9/2, 5/2, -1/2) (-9/2, 5/2, -1/2) 1 1 1 The effective Lagrangian that describes the up-quark sector in the model with the SM × U (1) P Q charges of Table I can be written as where C u ij represents coupling constant and Λ is the cut-off scale of the model. Similarly, one can write the Lagrangian for the down-quark sector with a similar meaning for C d ij and Λ as the up-quark sector. Notice that these Higgs doublets couple to both the quark sectors. After the symmetry breaking, and counting terms up to dimension-7, we find NNI-type quark-mass matrices as 3 where, ε = σ /Λ or σ * /Λ. For typical values of ε ∼ 0.2, one can safely neglect terms proportional to ε 4 and ε 8 as has been pointed in Eqs. 3, 4. In this framework the mass hierarchy between the third family and the first two is naturally explained by the fact that the (3, 3) entry is generated at dimension-4 and the rest of the entries are generated at dimension-5 . Additionally, one can explain the hierarchy between m t and m b by a hierarchy between v u and v d , as in two Higgs doublet models with natural flavor conservation.

The UV-completion: DFSZ type-I Seesaw
First, we consider a type-I Dirac seesaw model, where the heavy mediators are vector-like quarks.
Furthermore, this restriction will exclude them from contributing to the anomalous couplings of the axion. The UV-complete operator for the up-quark sector, as given by Eq. 3, within the DFSZ type-I seesaw formalism can be achieved as follows where F ij qC are the vector like fields and their PQ-charges are given in Table II. Similarly, the UV-completion of the down-quark sector as given by Eq. 4 can be achieved by using F ij qC (see Table II for their charges) and the corresponding Lagrangian can be written as In Fig. 1 we show the Feynman diagram corresponding to Eqs. 6 and 7.
3 Notice that the higher dimensional operators are either suppressed by a high enough energy scale compared to dimension-4 or 5, or one can forbid them by using an additional symmetry, and hence one can safely deal with the NNI-type quark texture. After the PQ and EW symmetry breaking by the vevs v σ , H u and H d we obtain the following quark mass matrices in the ( where the 4 submatrices of M u/d are given by Now, the light quark mass matrix in this Dirac type-I seesaw scenario can be written as at leading order. The resulting form of this matrix is given by where A, B, C, D, and E are complex entries.

The UV-completion: DFSZ type-II Seesaw
As an alternative, we consider a type-II Dirac seesaw framework of a DFSZ model. Within this model, the UV-completion of the quark-sectors (see Eqs. 3, 4) can be achieved by introducing Here, no additional heavy quark states are added. We express the UV-complete Lagrangian for the up-quark sector as follows: The SM × U (1) P Q charges for the remaining fields are given by Table I. Similarly, the down-quark sector can be written as The resulting quark mass matrices after the PQ and EW symmetry breaking can be read as where the vevs v Φ u/d of the additional doublets are determined by the scalar potential 4 . We can write at leading order both vevs as where M Φ u/d are heavy scalar masses. Finally, we like to remind here that the quark mass matrices in Eqs. 14 and 17 have a NNI structure.
UV complete diagram within the DFSZ type-II seesaw framework as apparent from Eqs. 15, and 16.

Lepton Sector
In our formalism, the Yukawa Lagrangian invariant under SM × U (1) P Q for charged-leptons and neutrinos is given by The charge assignments for the leptonic fields are given in Table III. This leads to diagonal charged-lepton as well as Dirac neutrino mass matrices, respectively. On the other hand, the Lagrangian involving right-handed neutrinos is given by 5 Therefore, the Majorana neutrino mass matrix takes the form Now, in the type-I seesaw formalism, a light neutrino mass matrix is given by 5 Notice that the same leptonic Lagrangian was obtained in the context of a gauged U (1) symmetry in [45]. and it corresponds to the type A 2 neutrino mass matrix.
It is to be noted here that the inclusion of σ allows a non-zero dimension-6 terms in the quark sector, which are zero at dimension-5 as given by Eq. 5. The Yukawa Lagrangian at dimension-6 for the up-quark sector can be written as However, we find that the UV-completion of the first term of Eq. 23 can only be achieved using the vector like fermions as given by Table II, whereas one needs new fermions to do the same for second and third terms, respectively. We explicitly show the UV-completion of the first term and the corresponding Lagrangian can be written as Similarly, the down-quark sector can be written as Here, we also find that one can achieve the UV-completion only for the first term of Eq. 25 and the UV-complete Lagrangian can be written as Therefore, in presence of dimension-6 terms, we generate a non-zero (1, 3) element of the quark mass matrix as given by Eq. 5. However, in our analysis, we focus on the leading dimension-4 and 5 contributions to the mass matrices, noting that given the matter content proposed, this is the minimal set of operators needed to obtain the fermion masses and mixings compatible with observations. We assume higher order operators are either suppressed by a high enough energy scale to be subdominant, or forbidden by an additional symmetry.

Scalar sectors
In minimal DFSZ models the SU (2)  We define the scalar fields present in the first DFSZ model as whereas the second DFSZ model contains two additional SU (2) L doublets Φ u and Φ d , which we define as The scalar potentials for both models have been provided in Appendix B. For each model we may write the Goldstone mode eaten by the Z-boson as where Y i is the hypercharge of each scalar field, v i its vev, for the type-I seesaw model A i ∈ {A, A , A u , A d } and A i ∈ {A, A , A u , A d , A u , A d } for the type-II seesaw. The Goldstone boson related to the PQ symmetry is likewise given by where X i is the PQ charge of the scalar field to which A i belongs. The physical axion, however, is orthogonal to the Z-boson Goldstone, so to obtain it from these two equations, we must perform the following subtraction [15] The possible values of v i for the doublet fields are bounded by the SM, i Y 2 i v 2 i = 246 GeV. On the other hand, the possible values of i X 2 i v 2 i are bounded from below by its relationship to the axion decay constant f a . This means that there must be a strong hierarchy between the vevs of the σ and σ fields and the vevs of the doublet fields. From Eqs. 29 -31 we can see that this results in the physical axion being composed to a high degree of only the σ and σ fields. Additionally, we choose to adopt the hierarchy v σ << v σ as it leads to negligible mixing between the A and A fields, resulting in a stronger coupling of the axion to the quark fields, as in our framework quark masses are derived from the couplings of σ at leading order. With these considerations in mind we extract the axion-quark couplings by setting a ∼ A. We would also like to mention that the details about the energy scale of f a and related phenomenology are analyzed in subsequent sections.

Axion couplings and mass
The couplings of the axion to the gluon and the photon are determined by the electromagnetic anomaly factor E and the QCD anomaly factor N. The effective Lagrangian that relates the axion couplings with anomalies (E and N ) can be read as where where this sum runs over all scalars with PQ charge X i and vev v i . For simplicity we will consider the dominant contribution to f a be from σ, so we may write f a ∼ 2v σ /N . The axion mass in terms of the decay constant is [47] m a = 5.70 µeV Finally, the coefficients of the anomalous axion-gluon coupling and axion-photon coupling are given by [47] In the Georgi-Glashow model [48] of gauge unification one can show that E/N = 8/3, using this in Eq. 36, the axion-photon coupling can be expressed as where the Fine-structure constant α EM = 1/137 has been used. Here, for the DFSZ models the SM fermions are the only contributors to the anomalies, resulting in N = 5 and E = 28/3. Therefore, compared to this benchmark the model presented here has the suppression For completeness, we briefly summarize here that the axion-photon interactions in the (g aγ − m a ) plane for m a ≥ 1 eV is mostly constrained from cosmology and astrophysics as can be seen from figure-1 of Ref. [49]. Besides this, for m a ≤ 1 eV various haloscope detectors put the tightest constraints on the axion mass and axion-photon coupling as has been outlined in figure We now proceed to discuss various low energy phenomenologies within the FNPQ-model in subsequent sections. As per the quark sector goes, we consider the quark mass matrix as given by Eq. 14 with a vanishing (1, 3) element and we adopt this as it can be seen from Eq. 5 that the (1, 3) element arises from dimension-6 which is 1/Λ 2 suppressed.

Masses and mixings of fermions
Here we perform a χ 2 analysis to find masses and mixing parameters of both quarks and leptons sectors. First, a global χ 2 fit is conducted to find the values of the parameters of the up and down quark matrices in our framework. The χ 2 function takes the following form where the sum runs over all observables. Also, µ f it represent the masses and mixings calculated using fitting parameters, as given by Eq. 40, we list them in Table IV. Here µ exp and σ exp are the observables and their standard deviation [5,52] as given by  and down-quark mass matrices that are used in the fit are given by Also, as pointed out in Appendix A, it is the difference in the up-and down-quark phase matrices that are relevant for the quark mixing matrix, hence in our numerical analysis, phases α d and β d have been fixed to zero. There are 12 parameters that need to be fitted, 10 amplitudes (5 for each matrix) and 2 phases. We fit these 12 parameters to account for the 10 physical observables related to them, the 6 quark masses, the 3 CKM angles and the CP-violating phase.
The observables are obtained from these matrices using the MPT package [53]. The fit is done at the energy scale M Z [52]. The initial values for the fitting procedure are randomized. The results from the best fit are given in Table IV, where χ 2 q = 1.0901. Following a similar approach, we also find the masses and mixing angles of the leptonic sector.
In this case only the elements of the neutrino mass matrix are fitted, since the charged lepton  best-fit as well as their 1σ error [5,52] for the various observables are given in the second and third columns of the lower table. Also, the best-fit values of the various observables are listed in the last column of the lower table.
matrix is already diagonal. The neutrino mass matrix has the A 2 -texture and takes the form The χ 2 function is identical as in Eq. 39, the only differences being the observables and their standard deviations. The leptonic masses and mixings obtained from the fit, which are compatible with the latest global fit data, can be seen in Table V at χ 2 l = 2.0053.

Flavor Violating decays with axions
The models described here contain flavor violating Yukawa couplings. Of particular interest are the axion flavor violating couplings to quarks as they lead to the decays with final state axions q i → q j a. These processes can be probed by meson decays with final state axions. These flavor violating couplings are common to both models presented. From the effective Lagrangian of Eqs.
3 and 4 we can write the Yukawa couplings to σ in the interaction basis as where the Yukawa coupling matrices y u( ) ij and y d( ) ij are In the quark mass basis the couplings to sigma are obtained by transforming the fermion fields to the physical basis where u /d denote quark mass eigenstates, U L and U R diagonalize the up quark mass matrix and V L and V R diagonalize the down quark mass matrix, respectively. By defining λ u( Finally, by neglecting axion mixing with other scalars, we write the quark couplings to the axion as where u,d ij = (λ ij − λ † ij )/2 and u,d ij = (λ ij + λ † ij )/2 . From these couplings we may calculate the branching ratio of flavor violating decays with axions. The most sensitive tests of neutral flavor violation with a final state axion is the K + → π + a process. The decay ratio for the Kaon decay to axion and pion is given by where B S is a non-perturbative parameter B S ∼ 4.6 [54]. To estimate the Kaon decay we use the best-fit point, as tabulated in Table IV, in Eqs. 44 -45 to find | d 21 | 2 , this leads to Using the latest constraint of the BR of Kaon decay from the E949 Collaboration [55] i.e., we may constrain v σ 2.5 × 10 10 GeV. For the models described above we have v σ ≈ √ 2f a N with N = 5, and hence using this approximation one can translate bounds of v σ to the axion decay constant f a , where f a 7 × 10 9 GeV .
Analogously, we may consider the B + → K + a decay, where the bottom to strange quark transition is probed. The decay width of this process is given by with f K 0 (0) ∼ 0.33 [56]. The experiment Belle-II may constrain the branching ratio of this process to BR(B + → K + a) < 10 −6 − 10 −8 [22], which would lead to v σ 1.8 × (10 7 − 10 8 ) GeV .
The axion decay constant f a > 7 × 10 9 GeV bound translates in our models to the equivalent limits m a < 0.7 × 10 −3 eV, and |g aγ (GeV −1 )| < 0.8 × 10 −14 . These bounds are two orders of magnitude stronger than the limits from astrophysics, see figure-1 of [49], for completeness. We end this section by noting that these bounds are obtained by neglecting σ − σ mixing. When this mixing is sizable these bounds on f a are relaxed.

Flavor Violating Higgs couplings
The where the matrices C u/d i are given by FIG. 4. Exclusion region plot in the (tan α − tan β) plane obtained from the non-observation of the t → hc flavor violating decay. The gray colored region is excluded by ATLAS data [58], the purple colored region is expected to be probed in the future by the HL-LHC experiment [59] and the orange region will be further probed by ILC or CLIC [60]. The uncolored region (see white thin band) predicts a branching ratio beyond the sensitivity of these experiments. The dashed line indicates limit of no flavor violation in light Higgs Yukawa couplings.
In terms of the physical scalar and quark fields we have where the matrices C u/d i are defined as while future experiments HL-LHC [59], ILC and CLIC [60] project the following sensitivities to this process Thus, the branching ratio can be calculated from Eq. 56 as [61] Γ where the coupling C tc is defined as m t is the top quark mass, R h is the higgs to top mass ratio R h = m h /m t , and R c is the charm to top mass ratio R c = m c /m t . Using the experimental value for the total width of the top quark [5] we derive the constraints on the free parameters tan α and tan β as shown in Fig. 4. We use the best fit point given in Table IV, to obtain a numerical value of Eq. 61 and we derive the approximate constraint cos α sin β for a given experimental input of the decay width Γ Exp t→hc . We see from Fig. 4 that, as expected from Eq. 63, small values for tan β allow only large values of tan α and vice versa. It can also be seen that ATLAS data has already ruled out a large portion of the parameter space (gray-region) and HL-LHC (purple-region) and CLIC (orange-region) will leave only a small region around the tan β = cot α limit unprobed (see Fig. 4 caption for details). Finally, we would like to mention that the t → hu and h → cu decays can also place constraints on α and β, however we find these to be numerically much weaker than the constraints from t → hc, given that the hierarchy in the C u i (see Eq. 55) matrices is preserved in the physical C u i matrix (see Eq. 57). Moreover, currently there are no experimental constraints on h → bs decays, although they are phenomenologically interesting in the context of two Higgs Doublet models [62].

Flavored Axion as Dark Matter candidate
The axion can be a good dark matter candidate, provided a sufficient amount of them was produced in the early universe. There are several production mechanisms for axions. The relic density produced by the misalignment mechanism is [63,64] Ω a h 2 ≈ 2 × 10 4 f a 10 16 GeV where θ a,i is the initial misalignment angle of the cosmological axion field and it lies in the range (0, 2π). Now, we notice that for the axion breaking scale 5 × 10 10 < f a < 1 × 10 15 (GeV), one can match the axion relic density to the observed dark matter relic abundance Ω DM h 2 ∼ 0.12 for 0 < θ a,i < 2π. It is worthwhile to mention that the N > 1 prediction of DFSZ models induce the formation of stable domain walls in the universe. The existence of these stable domain walls is incompatible with the standard cosmology [65]. One way to avoid the effect of domain walls on the observed universe is to embed this type of models in a cosmological model where inflation happens after the formation of these walls, thereby inflating away the density of the walls. Another possible resolution of the domain wall problem is to destabilize the walls with a dynamical breaking of the PQ symmetry [66,67].

CONCLUSION
In this work, to address the hierarchal flavor structure of fermion masses as well as the strong CP problem of QCD, we adopt a formalism based on the Froggatt-Nielsen symmetry. It is well known that the FN symmetry is one of the leading mechanisms to explain fermion masses and their mixings. On the other hand, the Peccei-Quinn symmetry provides an elegant solution to the strong CP problem. Here, an attempt has been made by identifying the PQ symmetry U (1) P Q with the FN symmetry U (1) F N to address both the drawbacks of the SM. In doing so, we assign PQ charges to the SM and BSM fields (see Table I), within the framework of the DFSZ style axion models, in such a manner that the (3, 3) entry of the quark mass matrices are generated by a dimension-4 operator, whereas remaining entries are obtained at the dimension-5 level. We end up with the Nearest-Neighbor-Interaction structure for the quark mass matrices in two different UV-completions, namely the type-I, and -II Dirac seesaw mechanisms for quark masses. Moreover, upon the breaking of the assigned symmetry, a flavored-axion from the flavon fields has been obtained, which solves the strong CP problem and also is a Dark Matter candidate. For the lepton sector, we also extended the model by considering the type-I seesaw mechanism. We obtained the A 2 -type Majorana neutrino mass matrix in the diagonal charged lepton mass basis.
A comprehensive numerical analysis has been performed to find low energy fermion masses and mixings, which are summarized in Table IV, V and are consistent with the latest global-fit data. Besides this, based on the ratio between the electromagnetic to the QCD anomaly factor, we have determined the axion-photon coupling, which is suppressed by a factor of 14 (as given by Eq. 38) compared to the SU (5) GUT model. The models discussed here also contain axion flavor violating couplings to quarks. Using these couplings, we have calculated the branching ratios of K + → π + a and B + → K + a processes that involve the flavored-axion. These indicate that the axion decay constant is f a > 7 × 10 9 GeV, whereas the axion mass limit is m a < 0.7 × 10 −3 eV, and the axion-photon coupling is |g aγ (GeV −1 )| < 0.8 × 10 −14 . We also have the possibility of flavor violating neutral currents mediated by scalar particles due to the presence of two Higgs doublets.
The branching ratio for the t → hc decay channel has been calculated for the model and the results are summarized in Fig. 4, where the latest bounds of ATLAS as well as the future experiments limits of HL-LHC, ILC and CLIC have been incorporated. Finally, we have pointed out that when the axion breaking scale lies in the range f a ∈ (5 × 10 10 , 1 × 10 15 ) GeV, then the axion relic density can be matched to the observed dark matter relic abundance. The up-and down-quark mass matrices as given by Eq. 14 can be written as with A, B, C, D, E, α and β as real parameters. The above mass matrix m u/d can be diagonalized by bi-unitary transformation of the form where L and R depict the left-and right-chiral fields, respectively. Also, V L = P L O and V R = P R O are the unitary matrices that diagonalize m † m and mm † and P L = diag(1, e iα , e iβ ), P R = diag(e iρ 1 , e iρ 2 , e iρ 3 ) are the diagonal phase matrices, respectively. Notice that we drop subscript (u/d) to demonstrate phase redefinition as the following formalism is same for both the up and down sector. Now, given the form of P L and P R , one can construct a real solution of the quark mass matrix (Eq. A1) following the transformation P † L mP R as mentioned by Eq. A2. Thus, the quark mass matrix in terms of real parameters can be written as Notice that for most of the application V R is irrelevant, it is the V L that enters in the CKM parameterization. Therefore, one can write quark mixing matrix as We see here that it is the difference in the diagonal-phase matrix that enters in the V CKM . Thus, in our numerical an analysis we have used the quark mass matrix of the form (A5)