A minimal axion model for mass matrices with five texture-zeros

A model with fermion and scalar fields charged under a Peccei–Queen (PQ) symmetry is proposed. The PQ charges are chosen in such a way that they can reproduce mass matrices with five texture zeros, which can generate the fermion masses, the CKM matrix, and the PMNS matrix of the Standard Model (SM). To obtain this result, at least 4 Higgs doublets are needed. As we will see in the manuscript this is a highly non-trivial result since the texture zeros of the mass matrices impose a large number of restrictions. This model shows a route to understand the different scales of the SM by extending it with a multi-Higgs sector and an additional PQ symmetry. Since the PQ charges are not universal, the model predicts flavor-changing neutral currents (FCNC) at the tree level, a feature that constitutes the main source of restrictions on the parameter space. We report the allowed regions by lepton decays and compare them with those coming from the semileptonic decays K±⟶πν¯ν\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^{\pm }\longrightarrow \pi {\bar{\nu }}\nu $$\end{document}. We also show the excluded regions and the projected bounds of future experiments for the axion–photon coupling as a function of the axion mass and compare it with the parameter space of our model.


I. INTRODUCTION
The discovery of the Higgs with a mass of 125 GeV, by the ATLAS [1] and CMS [2] collaborations, is very important because it provides experimental support for spontaneous symmetry breaking, which is the mechanism that explains the origin of the masses of fermions and gauge bosons.Additionally, it opens up the possibility of new physics in the scalar sector, such as the two Higgs doublet model [3][4][5][6][7][8][9][10][11][12][13][14][15], models with additional singlet scalar fields [16], or scalar fields that could be candidates for Dark Matter [15,[17][18][19][20]. On the other hand, in the Standard Model (SM) [21][22][23], symmetry breaking generates a coupling of the Higgs to fermions, proportional to their masses, which is consistent with experimental data.However, there are several orders of magnitude between the fermion mass hierarchies that cannot be explained within the context of the SM.Six masses must be defined for the up and down quarks, three Cabibbo-Kobayashi-Maskawa (CKM) mixing angles, and a complex phase that involves CP violation.On the other hand, in the lepton sector, there are three masses for charged leptons, two squared mass differences for neutrinos, three mixing angles, and a complex phase that involves CP violation in the lepton sector.In this case, it is necessary to determine the mass of the lightest neutrino and the character of neutrinos, whether they are Dirac or Majorana fermions.
In the Davis experiment [24], which was designed to detect solar neutrinos, a deficiency in the solar neutrino flux was first observed.According to the results of Bahcall, only one-third of solar neutrinos would reach the Earth [25].Neutrino oscillation was first proposed by Pontecorvo [26], and the precise mechanism of solar neutrino oscillations was proposed by Mikheyev, Smirnov, and Wolfenstein, involving a resonant enhancement of neutrino oscillations due to matter effects [27,28].These observations have been confirmed by many experiments from four different sources: solar neutrinos as in Homestake [24], SAGE [24], GALLEX & GNO [29,30], SNO [31], Borexino [32,33] and Super-Kamiokande [34,35] experiments, atmospheric neutrinos as in IceCube [36], neutrinos from reactors as KamLAND [37], CHOOZ [38], Palo Verde [39], Daya Bay [40], RENO [41] and SBL [42], and from accelerators as in MINOS [43], T2K [44] and NOνA [45].Neutrino oscillations depend on squared mass differences.On the other hand, the lightest neutrino mass has not been determined yet, but from cosmological considerations, none of the neutrino masses can exceed 0.3 eV, which implies that the neutrino masses are much smaller than the charged fermion masses.However, unlike quarks and charged leptons, in the SM the neutrinos are massless, which is explained by assuming that neutrinos are left-handed.Therefore, the discovery of neutrino masses implies new physics beyond the SM.By adding righthanded neutrinos, the Higgs mechanism of the SM can give neutrinos the same type of mass acquired by charged leptons and quarks.It is possible to add right-handed neutrinos ν R to the SM, as long as they do not participate in weak interactions.With the presence of right-handed neutrinos, it would be possible to generate Dirac masses m D , similar to those of charged leptons and quarks.In principle, it is also possible to give Majorana masses to left-handed neutrinos, and similarly, right-handed neutrinos can have Majorana masses M R .For a very large M R , it would give effective Majorana masses for left-handed neutrinos as m eff ≈ m 2 D /M R .The presence of large Majorana masses allows to explain the tiny neutrino masses compared to the charged fermion masses [38].To explain the smallness of neutrino masses, there are three types of seesaw mechanisms in the literature: type I with three electroweak neutrinos and three heavy right-handed neutrinos, type II [46,47], type III [48], and inverse seesaw [49,50].One way to explain the fermion mass hierarchies and the CKM and PMNS mixing angles is through zeros in the Yukawa couplings of fermions (this is known as texture-zeros or simply textures of the mass matrices, and these zeros are usually chosen by hand).It is common in the literature to consider Fritzsch-type textures [51,52], or similar [53][54][55][56][57][58], for the neutrino and charged lepton mass matrices.
Due to the fact that there are three up-type quarks and three down-type quarks, the mass operators are 3×3 complex matrices with 36 degrees of freedom.If we consider these operators to be Hermitian [117][118][119], the number of free parameters reduces to 18, which cannot be fully determined from the 10 available physical quantities, namely masses and mixing angles [120].This provides freedom to reduce the number of free parameters in the matrices and search for matrix structures with zeros that provide eigenvalues and mixing angles consistent with the masses and mixing matrices of the fermions.One way to find zeros in the mass matrices that is automatically consistent with experimental data is based on weak basis transformations (WBT) for quarks and leptons [118,119,121,122].Fritzsch proposed an ansatz with six zeros [123,124,[124][125][126][127][128], but the value of |V ub /V cb | ≈ 0.06 is too small compared to the experimental value |V ub /V cb | exp ≈ 0.09 [129].For this reason, the use of 4 and 5 zero-textures was proposed [117,118,128,[130][131][132][133].Reference [117,119] showed that matrices with five zero-textures could reproduce the mass hierarchy and mixing angles of the CKM matrix.
The strong CP problem arises from the fact that the QCD Lagrangian has a non-perturbative term ("θ-term") that explicitly violates CP in strong interactions.Peccei and Quinn proposed a solution to the strong CP problem [134,135], where it is assumed that the SM has an additional global chiral symmetry U (1), which is spontaneously broken at a large energy scale f a .One consequence of this breaking is the existence of a particle called the axion, which is the Goldstone boson of the broken U (1) P Q symmetry [136,137].Due to the fact that the PQ symmetry is not exact at the quantum level, as a result of a chiral anomaly, the axion is massive and its mass (see Appendix D) is given by: where z = 0.56 is assumed for the up and down quark mass ratio, while f π ≈ 92 MeV and m π = 135 MeV are the pion decay constant and mass, respectively.The effective couplings of axions to ordinary particles are inversely proportional to f a , and also depend on the model.It was originally thought that the PQ symmetry breaking occurred at the electroweak scale, but experiments have ruled this out.The mass of the axion and its coupling to matter and radiation scale as 1/f a , making its direct detection extremely difficult.The combined limits from unsuccessful searches in nuclear and particle physics experiments and from stellar evolution imply that f a ≥ 3 × 10 9 GeV [138].Furthermore, there is an upper limit of f a ≤ 10 12 GeV that comes from cosmology, since light axions are produced in abundance during the QCD phase transition [139][140][141][142][143]. Hence, these models are generically referred to as "invisible" axion models and remain phenomenologically viable.There are two classes of invisible axion models in the literature: KSVZ (Kim, Shifman, Vainshtein, and Zakharov) [138,144] and DFSZ (Dine, Fischler, Srednicki, and Zhitnitsky) [145,146].The main difference between KSVZ-type and DFSZ-type axions is that the former do not couple to ordinary quarks and leptons at tree level, but instead require an exotic quark that ensures a nonzero QCD anomaly to generate CP violation.Depending on the assumed value of f a , the existence of axions could have interesting consequences in astrophysics and cosmology.The emission of axions produced in stellar plasma through their coupling to photons, electrons, and nucleons would provide a new mechanism for energy loss in stars.This could accelerate the evolutionary process of stars and, therefore, shorten their lifespan.Axions can also exist as primordial cosmic relics produced copiously in early times and could be candidates for dark matter.Of course, if axions are candidates for dark matter [147][148][149][150][151][152][153], they should be very light, with m a ≤ 10µeV .However, data from numerous laboratory experiments and astrophysical observations, together with the cosmological requirement that the contribution to the mass density of the Universe from relic axions does not saturate the Universe, restricts the allowed values of the axion mass to a range of 10 −5 eV < m a < 10 −2 eV.One source of axions would be the Sun, which, coupled to two photons, could be produced through the Primakoff conversion of thermal photons in the electric and magnetic fields of the solar plasma.The limits are primarily useful for complementing the arguments of stellar energy loss [154] and the searches for solar axions by CAST at CERN [155] and the Tokyo axion helioscope [156].
The axion-photon coupling (see Appendix D) can be calculated in chiral perturbation theory as [134,135].
This coupling and the axion mass are related to each other through the relation E/N , which depends on the model and can be tested in experiments.
The strongest limits on the axion-electron coupling are derived from observations of stars with a dense core, where bremsstrahlung is very effective.These conditions are realized in White Dwarfs and Red Giant Stars, where the evolution of a White Dwarf is a cooling process by photon radiation and neutrino emission, with the possible addition of new energy loss channels such as axions.Current numerical analysis suggest a limit of g ae ≤ 2.8 × 10 −13 [157].In particular, using data from the Sloan Digital Sky Survey (SDSS) and SuperCOSMOS Sky Survey (SCSS) [158], they showed that the axion-electron coupling is approximately 1.4 × 10 −13 .More recent analysis of the data in Ref. [158] found g ae = 1.5 +0.6 −0.9 × 10 −13 (95% CL) [159,160].The two groups studying the axion-electron coupling are M5 [154] and M3 [161].Their combination yields the limit g ae = 1.6 +0. 29 −0.34 × 10 −13 .For a recent and comprehensive review of axion physics, see [157].
This document is organized as follows: In Section 2, we review the textures for the quark and lepton mass matrices that will be used in this work.We also write the real parameters of these matrices in terms of the masses of the SM fermions and two free parameters.In Section 3, we present the particle content of our model and the necessary PQ charges to generate the mass matrix textures presented in Section 2. In Section 4, we adjust the Yukawa couplings to obtain the masses of the charged leptons and neutrinos.It is important to note that we cannot use the VEVs to adjust the lepton masses, as these were already adjusted to reproduce the quark masses.It is also important to note that by using a seesaw mechanism, we can avoid adjusting the Yukawas, however, that is not our purpose in the present work.In Section 5, we show the Lagrangian of our model.In Section 6, we present some constraints in the parameter space, as well as projected constraints for upcoming experimental results, both for experiments under construction and in the data-taking phase.

II. THE FIVE TEXTURE-ZERO MASS MATRICES
The reason for dealing with texture zeros in the Standard Model (SM) and its extensions is to simplify as much as possible the number of free parameters that allow us to see relationships between masses and mixings present in these models.The Yukawa Lagrangian is responsible for giving mass to SM fermions after spontaneous symmetry breaking.A first simplification, without losing generality, is to consider that the fermion mass matrices are Hermitian, so the number of free parameters for each sector of quarks and leptons is reduced to 18, but there is still an excess of parameters to reproduce the experimental data provided in the literature.Due to the lack of a model to make predictions, discrete symmetries can be used to prohibit some components in the Yukawa matrix, generating the so-called texture zeros for the mass matrices.In many works, instead of proposing discrete symmetries, texture zeros are proposed as practical and direct alternatives.The advantage of this approach is that it is possible to choose each mass matrix in an optimal way for the analytical treatment of the problem, and at the same time adjust the mixing angles and the masses of the fermions.

A. Quark Sector
We should keep in mind that six-zero textures in the SM have already been discarded because their predictions are outside the experimental ranges allowed; but, five-zero textures for quark mass matrices is a viable possibility [61,121,[162][163][164][165].Specifically, we chose the following five-zero textures because they fit well with experimental quark masses and mixing parameters [117,119,166]: In addition, the phases in M D can be removed by a weak basis transformation (WBT) [117,118,122], so that they are absorbed by the off-diagonal terms in M U .In this way, the mass matrices (3) can be rewritten as: By applying the trace and the determinant to the mass matrices (4), before and after the diagonalization process, the free real parameters of M U and M D can be written in terms of their masses: A possibility that works very well is to consider the second generation of quark masses to be negative, i.e., with eigenvalues −m c and −m s .And A u is a free parameter, whose value, determined by the quark mass hierarchy, must be in the following range: The exact analytical procedure for diagonalizing the mass matrices ( 4) is indicated in Appendix C.

B. Lepton sector
In this work, we will consider Dirac neutrinos.This is achieved, in part, by extending the SM with right-handed neutrinos.In this way, we can carry out a treatment similar to that of the quark sector, that is, the mass matrices of the lepton sector can be considered Hermitian and the weak basis transformation (WBT) can be applied [117,118].In the literature, work has been done considering various texture-zeros for the Dirac mass matrices of the lepton sector [59,[167][168][169][170][171][172][173][174][175][176][177][178].In our treatment, we are going to consider the following five-zero texture model studied in the paper [132], which can accurately reproduce the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix V PMNS (mixing angles and the CP violating phase), the charged lepton masses, and the squared mass differences in the normal mass ordering.
Without loss of generality, by using a WBT, the phases of the charged lepton mass matrix, M E , can be absorbed into the entries C ν and B ν of the neutrino mass matrix, M N .Similarly, as was done in the case of the quark sector, the parameters present in the mass matrices of the lepton sector ( 7) can be expressed in terms of the masses of the charged leptons m e , m µ and m τ and the masses of the neutrinos m 1 , m 2 and m 3 , in the normal ordering (m 1 < m 2 < m 3 ): where the values of the masses and the parameter A ν are given in Table V.Furthermore, for the adjustment of the mass matrices (7) it is very convenient to assume that the eigenvalues associated with the masses of the second family, −m 2 and −m µ , are negative quantities.The exact diagonalizing matrices of the mass matrices (7) are shown in appendix C, equations (C6) and (C5).

III. PQ SYMMETRY AND THE MINIMAL PARTICLE CONTENT A. Yukawa Lagrangian and the PQ symmetry
The texture-zeros of the mass matrices defined in the equations ( 4) and ( 7) can be generated by imposing a Peccei-Queen symmetry U (1) P Q on the Lagrangian model, Eq. ( 9) [65,179,180].As will be explained below, the minimal Lagrangian that allows us to implement this symmetry is given by [64,181] As it was shown in reference [64], at least four Higgs doublets are required to generate the quark mass textures, therefore α = 1, 2, 3, 4. In (9) i, j are family indices (there is an implicit sum over repeated indices).The superscripts U , D, E, N refer to up-type quarks, down-type quarks, electron-like and neutrino-like fermions, respectively; and D µ = ∂ µ + iΓ µ is the covariant derivative in the SM.The scalar potential V (Φ, S 1 , S 2 ) is shown in appendix A (for further details, see reference [64]).In Eq. ( 9) ψ stands for the SM fermion fields plus the heavy quark Q (see Tables I  and II).As it is shown in Table II the PQ charges of the heavy quark can be chosen in such a way that only the interaction with the scalar singlet S 2 is allowed.We assign Q PQ charges for the left-handed quark doublets (q L ): x qi , right-handed up-type quark singlets (u R ): x ui , right-handed down-type quark singlets (d R ): x di , left-handed lepton doublets ( L ): x i , right-handed charged leptons (e R ): x ei and right-handed Dirac neutrinos (ν R ): x νi for each family (i = 1, 2, 3).We follow a similar notation for the scalar doublets, x φα (α = 1, 2, 3, 4), and the scalar singlets In this work, the PQ charges assigned to the quark sector and the scalar sector, as well as the VEVs assigned to the scalar doublets, will be the same as those assigned in [64] (Tables (I) and (II)), and we will adjust the PQ charges of the lepton sector to reproduce the texture-zeros given in Eq. (7).To forbid a given entry in the lepton mass matrices, the corresponding sum of PQ charges must be equal to zero, so that we can obtain texture-zeros by imposing the following conditions: where Since the PQ charges of the Higgs doublets (α = 1, 2, 3, 4) are already given, the possible solutions of ( 10) and ( 11) are strongly constrained.Table I provides a solution for the PQ charges of the lepton sector.
The subindex i = 1, 2, 3 stand for the family number in the interaction basis.Columns 6-8 are the Peccei-Quinn charges, QP Q, for each family of quarks and leptons in the SM.s1, s2 and α are real parameters, with s1 = s2.In our model we include two scalar singlets S 1 and S 2 that break the global symmetry U (1) P Q .The QCD anomaly of the PQ charges is where A Q = x QL − x QR is the contribution to the anomaly of the heavy quark Q, which is a singlet under the electroweak gauge group, with left (right) Peccei-Quinn charges x QL,R , respectively.We can write the charges as a function of N (since N must be different from zero), such that where ŝ1 and are arbitrary real numbers.To solve the strong CP problem with N = 0 and simultaneously generate the texture-zeros in the mass matrices, it is necessary to maintain = 9(s1−s2) N = 0.With these definitions for Flavor-Changing Neutral Currents (FCNC) observables, the relevant parameters are ŝ1 and .This parameterization is quite convenient (for those cases where the parameters α q and α are not relevant) because by fixing N and f a , we can vary ŝ1 and for a fixed Λ PQ = f a N in such a way that the parameter space naturally reduces to two dimensions.

IV. NATURALNESS OF YUKAWA COUPLINGS
A. The mass matrices in the quark sector In reference [64], it was shown that to generate five texture zeros in the quark mass matrices (3), as a consequence of a PQ symmetry, it is necessary to include at least four scalar doublets in the model.After spontaneous symmetry breaking, the quark sector mass matrices take on the following form: where the vi are defined in terms of the vacuum expectation values, vi = v i / √ 2. In [64] it was shown that the five-texture zeros ( 4 As we can see, the hermiticity of the mass matrices is not fully achieved, but it is good to impose it for several reasons: (i) In the SM and its extensions, in which the right chirality fields are singlets under SU (2), the mass matrices can be assumed Hermitian without losing generality, (ii) the previous fact allows us to consider Hermitian mass matrices, even after imposing an additional P Q symmetry in the model, (iii) we can implement the WBT method [117], and (iv) there is an extensive literature on physically viable Hermitian mass matrices.It is important to noticing that the mass matrices in Eq. ( 14) are Hermitian.

B. The mass matrices in the lepton sector
We can obtain the lepton mass matrices by starting from the Yukawa Lagrangian (9), which is invariant under the Peccei-Quinn U (1) P Q symmetry, and taking into account the Yukawa parameters and expectation values (15).After the spontaneous symmetry breaking, the mass matrices for neutral and charged leptons are given respectively by [132]: As we previously mentioned, at least four Higgs doublets are needed to obtain the five texture-zeros for the chosen quark mass matrices.Our goal in this work is to keep the same number of Higgs doublets and their respective P Q charges to generate the mass matrices and texture zeros for the lepton sector, Eq. ( 7).To get an Hermitian mass matrix M N , it is necessary to impose y N On the other hand, to obtain a symmetric mass matrix, M E , for the charged leptons, it is sufficient to assume that the Yukawa couplings are Hermitian.Through these choices it is possible to avoid additional Higgs doublets.
Based on the results of Table V, Appendix C, and the relationships established in (8), we find the following values for the Yukawa couplings of the lepton sector: To reproduce the neutrino masses quoted in [132], in the SM is required a Yukawa coupling around 10 −14 .In our case, the smallest Yukawa coupling is 10 −7 , which significantly reduces the fine-tuning in comparison to that given by the SM.

V. THE EFFECTIVE LAGRANGIAN
The strongest constraints on non-universal PQ charges come from the FCNC.To determine these constraints, we start by writing the most general next-to-leading order (NLO) effective Lagrangian as [182,183]: c aΦ α and c 1,2,3 are Wilson coefficients; α 1,2,3 = 4π , where g 1,2,3 are the coupling strengths of the electroweak interaction in the interaction basis; and the Wilson operators are: where B, W a and G a correspond to the gauge fields associated with the SM gauge groups U (1) Y , SU (2) L and SU (3) C , respectively.a is the axion field which corresponds to the CP odd component of S 1 .It is possible to redefine the fields by multiplying by a phase In this expression, x ψ corresponds to the PQ charges of the SM fermions, i.e., {x ψ L,R } = {x qi , x ui , x di , x li , x ei , x νi } and {x Φ α } are the PQ charges of the Higgs doublets {Φ α }.Replacing these definitions in the kinetic terms of Eq. ( 9), we obtain new contributions to the effective Lagrangian Eq. ( 18) (the NLO contributions in the non-derivative terms cancel out).The leading order (LO) terms in Λ −1 PQ can be written as [181,182]: where with The field redefinitions (20) induce a modification in the measure of the functional path integral whose effects can be obtained from the divergence of the axial-vector current: where the hypercharge is normalized by Q = T 3L + Y .The relation ( 24) is an on-shell relation, which is consistent with the momentum of an on-shell axion.Substituting this result into we obtain new contributions to the leading-order Wilson coefficients [185] where Σq ≡ x q1 + x q2 + x q3 .The corresponding NLO Lagrangian is It is convenient to define c eff 3 = c 3 − 2Σq + Σu + Σd − A Q = −N .In our case, c i = 0 and the only contributions to c eff i come from the anomaly.It is customary to define Λ PQ = f a |c eff 3 | to include the factor c eff 3 in the normalization of the PQ charges.From now on, we will assume that all the PQ charges are normalized in this way, so that x ψ corresponds to x ψ /c eff 3 .For normalized charges, c eff 3 = 1, therefore, we still maintain the general form despite writing all the expressions in terms of the effective scale f a .
The scalar fields and their PQ charges are the same as in the reference [64], so the scalar potential V (Φ, S) is identical to that of the mentioned reference.With the VEVs and couplings given in [64], the model reproduces the mass of the SM Higgs, while the masses of the exotic scalars are above the TeV scale.This potential has the appropriate number of Goldstone bosons to give masses to the SM gauge bosons Z 0 , W ± and has an extra field that can be identified with the axion a.

A. Flavor changing neutral currents
Due to the non-universal PQ charges in our model, a tree-level analysis of flavor-changing neutral currents is necessary.As mentioned in reference [157], the strongest limits on the axion-quark FCNC couplings come from meson decays in light mesons and missing energy.The decays K ± → π ± a provide the tightest limits (NA62 Collaboration [186]) for the axion mass [157].Currently the most restrictive limits come from the semileptonic decays of kaons K ± → π ± νν and leptons 1 → 2 +missing energy.From the term ∆L K ψ , we obtain the vector and axial couplings for a multi-Higgs sector model, as shown in references [64,157] where Collaboration Upper bound N62 Collaboration [186] B K + → π + a < (10.6 +4.0 −3.4 stat ± 0.9syst) × 10 −11 TRIUMF [189] B µ + → e + a < 2.6 × 10 −6 Crystal Box [190] B µ + → e + γa < 1.1 × 10 −9 where . In these expressions, F stands for U, D, N or E and the U F L,R are de diagonalizing matrices (see Appendix C).In Eq. ( 28), we normalize the charges using c eff 3 , as explained in the last paragraph of section V (in other references |c eff 3 | = |N | is considered, corresponding to the SU (3) × U (1) P Q anomaly).The branching ratio for lepton decays i → j a is given by [187] Br in this expression, the vector and axial couplings contribute in the same way where m is the mass of the leptons and Γ( i ) is the total decay width of the particle j .
For the lepton decay i → j a γ, we can relate this branching ratio to the branching ratio of the process without the photon in the final state, according to the expression: where α is the fine structure constant, and the function, f (x, y) = (1−x)(2−y−xy) y 2 (x+y−1) , depends on the mass and the energies x = 2E 2 /m 1 and y = 2E γ /m 1 .For the lepton decay µ → e a γ, the constraints come from the Crystal Box experiment [188], with cut energies E γ , E e > 30 MeV, θ eγ > 140 • , where cos θ eγ = 1 + 2(1−x−y) xy , so that dx dyf (x, y) ≈ 0.011.In our model, there is a natural alignment between the Φ 3 (which is quite similar to H 1 in the Georgi basis [192]) and the standard model Higgs boson as a consequence of the large suppression of the VEVs of the scalar doublets v i , with i = 1, 2, 4, respect to v 3 , the VEV of Φ 3 .To some extent, this alignment avoids FCNC involving the SM Higgs boson [192]; however, after alignment, there are other sources of FCNC associated with the additional scalar doublets, which cannot be avoided by any means; however, as argued in Ref. [64] they are suppressed by a factor 1/M 4 (where M > 1TeV is the mass of the exotic scalar doublets), and therefore, our model avoids these potential sources of FCNC in agreement with the general argument presented in [192].
From astrophysical considerations we have: bounds from black holes superradiance and the SN 1987A upper limit on the neutron electric dipole moment, which, when combined, impose a constraint on the axion decay constant in the range [157] (see Figure 1) : 0.8 × 10 6 GeV ≤ f a ≤ 2.8 × 10 17 GeV.

B. Constraints on the axion-photon coupling
There are several experiments designed to look for exotic particles.The sources studied in the search for axions are: the solar axion flux (helioscopes experiments), dark matter halo (haloscopes experiments), and axions produced in the laboratory.

VII. SUMMARY AND CONCLUSIONS
We have presented a model in which the fermion and scalar fields are charged under a U (1) P Q Peccei-Quinn symmetry.A recent work [64] showed that at least four Higgs doublets are required to generate Hermitian mass matrices in the quark sector with five texture-zeros, reproducing the quark masses, the mixing angles, and the CPviolating phase of the CKM mixing matrix.In this work, we show that using the same number of Higgs doublets, without changing the PQ charges in the quark and Higgs sectors, it is possible to generate Hermitian mass matrices in the lepton sector that reproduce the neutrino mass-squared differences in the normal mass ordering, the mixing angles, and the CP-violating phase of the PMNS mixing matrix.This result is quite non-trivial as we maintain the same four Higgs doublets required in the quark sector to generate a different texture pattern in the lepton sector.When compared to the SM, our model has almost all Yukawa couplings close to 1 in the quark sector.In the neutrino sector, the smallest Yukawa coupling is of the order of 1.6 × 10 −7 , which is seven orders of magnitude larger than the corresponding Yukawa coupling in the SM, so it requires less fine-tuning than the SM.
The polar decomposition theorem [216,217] allows any matrix to be written as the product of a Hermitian matrix and a unitary matrix.In the SM and in theories where the right-handed fermion fields are singlets under the gauge group, it is possible to absorb the unitary matrix into the right-handed fields by redefining them; from this procedure, we can write any mass matrix as a Hermitian matrix.In our work, we assume that the mass matrices are Hermitian in the interaction space, this hypothesis has been used in previous studies on textures [61,117,121,[218][219][220], and it is quite useful for studying the flavor problem.In our work, we have normalized the PQ charges with the QCD anomaly −N in such a way that by keeping the parameter = 0, we obtain the textures of the mass matrices, addressing the flavor and strong CP problems simultaneously.
If nature is not fine-tuned in a more fundamental high-energy theory, we expect that, eventually, it will be possible to find a texture that allows us to obtain all the scales of the SM from the VEVs of a Higgs sector with a minimal scalar content without the need to adjust the Yukawa couplings.
In our analysis, we report the constraints from lepton decays and compare them with the constraints from the search for neutrino pairs in charged Kaon decays K ± −→ π ± νν.The results are shown in figure 1, where the allowed region in the parameter space generated by and the axion decay constant f a is displayed.This figure shows that the strongest constraints come from the semileptonic meson decay K ± → πν ν.It is important to note that the lepton decays do not further constrain the parameter space of our model (compared to the region excluded by the meson decay).We also show the excluded regions for the axion-photon coupling as a function of the axion mass; these results are summarized in Figure 2; the gray region corresponds to the parameter space of our model in the interval −1 < < 1.
In this article, we have demonstrated that with four Higgs doublets, it is possible to fit the textures of the mass matrices, both in the lepton and quark sectors.These matrices generate the masses and the mixing matrices for quarks and leptons within the experimentally reported values in the literature.The introduction these doublets improves the fine-tuning problem of the Yukawa couplings and shows that this approach is a viable way to tackle the flavor problem.We hope to improve our results in future work by using the Seesaw mechanism in the lepton sector.
where a sum is assumed on repeated indices.Here i, j run over 1, 2, 3 and α over 1, 2, 3, 4. The Higgs boson doublet fields are parameterized as follows: Similar to the two Higgs doublet model [221] we rotate the Higgs fields to the (generalized) Georgi basis, that is, where the orthogonal matrices where tan , tan and tan β 3 = v4 v3 , and T .This basis is chosen in such a way that only the neutral component of H 1 acquires a vacuum expectation value and F = U, D, N, E; where we have defined With these definitions, equation (A1) becomes It is necessary to rotate to the fermion mass eigenstates, i.e., where the diagonalization matrices U L,R are defined below, in section C. From the Lagrangian for the charged currents it is possible to obtain the CKM (V CKM = U U L U D † L ) and PMNS (V PMNS = U E L U ν † L ) mixing matrices by rotating to the fermion mass eigenstates.In particular, we are interested in the coupling of the axial neutral current to the axion in the mass eigenstates. where . In these expressions the mass functions in the interaction basis are: where v = H 0 1 is the Higgs vacuum expectation value.

Appendix B: Scalar potential
As studied in [64], the scalar sector requires four scalar doublets φ α to reproduce the mass textures of the fermion sector correctly, and two scalar singlets S 1 and S 2 that break the PQ symmetry while generating a phenomenologically viable scalar mass spectrum.The S 2 singlet also gives mass to the heavy quark.The most general potential allowed by the PQ symmetry according to the charges established in Table II is: where the terms proportional to F i are allowed by the particular choice of PQ charges and these couplings F i have units of mass.After spontaneous symmetry breaking (SSB), the four Higgs doublets acquire VEVs that give mass to all the SM particles.The scalar doublets and singlets are written as follows: where the VEVs satisfy the following hierarchy: The scalar singlets S 1 and S 2 break the PQ symmetry at the high energy scale given by v s1 ≈ v s2 .The last two terms in equation (B1) correspond to the softbreaking masses of the imaginary and the real parts of S 2 , which are generated at one loop in the Coleman-Weinberg potential from the interaction term λ Q S 2 QR Q L + h.c.Additionally, we choose numerical values for the parameters of the potential (B1) in order to obtain a scalar sector mass spectrum consistent with the existing phenomenology.The values of these parameters are: In particular, the value of λ 3 adjusts the SM Higgs mass.The v i are determined from the SM fermion masses and the quark mass matrix textures, Eq. ( 15).The VEV v s1 remains a free parameter; however, this parameter is important for the axion physics due to the relationship [159], In our calculations we took v s1 ≈ v s2 ≈ 10 6 GeV.It is important to emphasize that in our model, f a can take arbitrary values; nevertheless, a small f a restricts (Eq.13) to values close to zero.Taking into account all these considerations, including Eq. (B3), the scalar mass spetrum (in GeV) is: CP even = {1.73× 10 The mass spectrum of the scalar fields is above the TeVs scale, except for the SM Higgs, which is at 125 GeV.The pseudoscalar sector (CP odd fields) have two massless eigenstates, the axion field and the Goldstone boson which is absorbed by the longitudinal component of the SM Z boson.A similar result is obtained in the charged sector, where it is possible to identify the two Goldstone bosons required to give mass to the SM W ± fields.

Appendix C: diagonalization matrices
To compare with physical quantities, it is necessary to rotate fields to the mass eigenstates, i.e., f L,R = U F L,R f L,R , where the prime symbol stands for the interaction basis.In our formalism the quark mass matrices are Hermitian, so the right-and left-handed diagonalizing matrices are identical; additionally, we establish that the eigenvalues of the second family of quarks are negative in order to generate texture-zeros in some diagonal terms of the mass matrices, as indicated in [118].This sign is taken into account by introducing the identity matrix written as I 2 I 2 = 1 with where the U F diagonalization matrices are defined below.It is important to stress that the texture-zeros pattern in the matrix Y F 1 ij are identical to those in the original Yukawa couplings y F α ij , since the sum over α does not mix the i, j indices.In fact, according to equations ( 14) and (16) The diagonalization matrices are: Au(mc+mu)(mt−mu) −e i(φ Bu +θ2u) mc(mt−Au)(Au−mu) Au(mc+mt)(mc+mu) e i(φ Bu +θ3u) (Au+mc)mt(Au−mu) where θ 1u , θ 2u , θ 3u , θ 1d and θ 2d are arbitrary phases (a third phase for the diagonalization matrix (C4) can be absorbed by the remaining phases) that are useful for conforming to the V CKM = U U L U D † L matrix convention.Taking as input the SM parameters at the Z pole, the best fit values are given in Table IV Similarly, in the lepton sector, the diagonalization matrices of the mass matrices (7) are: where θ 1 , θ 2 , θ 1ν , θ 2ν , θ 3ν are necessary phases in order to adjust to the established convention for the PMNS mixing matrix [222] 1 ; and c ν and b ν are the phases of C ν and B ν in the neutral mass matrix M N in Eq. (7).The best fit values for these quantities are shown in Table V where Σf ≡ f 1 + f 2 + f 3 is the sum of the PQ charges of the three families.There are similar definitions for the interaction of the axion with the gluons where c eff 3 = c 3 − 2Σq + Σu + Σd − A Q , in our particular case c i = 0.In axion phenomenology, it is usual to define The decay widths of an axion decaying in two photons and a Z decaying in an axion and a photon are [184] Γ(a → γγ) = 4πα The axion-photon interaction is given by where α = e 2 4π .Due to the gluon-axion interaction, the axion gets a mass term, which is described at low energies as an axion-pion interaction [224] m a = 5.7 (7) ) are flexible enough to set the quark Yukawa couplings close to 1 for most of them (except for y U 2 23 , y D3 23 and y U 1 13 ), in this way we obtain: v1 = 1.71GeV, v2 = 2.91 GeV, v3 = 174.085GeV, v4 = 13.3MeV.

FIG. 2 :
FIG.2:The excluded parameter space by various experiments corresponds to the colored regions, the dashed-lines correspond to the projected bounds of coming experiments looking for axion signals.The gray region corresponds to the parameter space scanned by our model.

TABLE II :
Beyond the SM fields and their respective PQ charges.The parameters s1, s2 are reals, with s1 = s2 and xQ R = xQ L .

TABLE III :
These inequalities come from the window for new physics in the branching ratio uncertainty of the meson decay in a pair νν.

TABLE IV :
. MeV 1.2684 MeV 633.197MeV 171268 MeV 3.14751 MeV 56.1169 MeV 2910.01MeV Best-fit point of the mass matrix parameters with respect to experimental data for the masses and mixing angles of the quark sector at the Z pole.

TABLE V :
. Best fit values.