Sequentially loop-generated quark and lepton mass hierarchies in an extended Inert Higgs Doublet model

Extended scalar and fermion sectors offer new opportunities for generating the observed strong hierarchies in the fermion mass and mixing patterns of the Standard Model (SM). In this work, we elaborate on the prospects of a particular extension of the Inert Higgs doublet model where the SM hierarchies are generated sequentially by radiative virtual corrections in a fully renormalisable way, i.e. without adding any non-renormalisable Yukawa terms or soft-breaking operators to the scalar potential. Our model has a potential to explain the recently observed $R_{K}$ and $R_{K^{\ast }}$ anomalies, thanks to the non universal $U_{1X}$ assignments of the fermionic fields that yield non universal $Z^{\prime}$ couplings to fermions. We explicitly demonstrate the power of this model for generating the realistic quark, lepton and neutrino mass spectra. In particular, we show that due to the presence of both continuous and discrete family symmetries in the considered framework, the top quark acquires a tree-level mass, lighter quarks and leptons get their masses at one- and two-loop order, while neutrino masses are generated at three-loop level. The minimal field content, particle spectra and scalar potential of this model are discussed in detail.


I. INTRODUCTION
The origin of various strong hierarchies in the fermion spectra of the Standard Model (SM) still remains a major unsolved problem of contemporary Particle Physics. A symmetry-based understanding of such hierarchies, in the framework of a single universal mechanism, consistent with the current phenomenological bounds on New Physics models, poses a challenging problem for the model-building community. In fact, a number of different potentially realizable mechanisms have been proposed so far, typically with certain limitations and deficiencies, and usually treating the quark, lepton and neutrino hierarchies on a separate footing. The most promising scenarios rely on the existence of horizontal (family) symmetries acting in the space of fermion generations and offering most times a more universal approach to the "fermion hierarchy problem" than other methods. It seems natural and attractive to consider the viable prospect of a high-scale spontaneous discrete symmetry breaking, triggering the radiative generation of new mass operators in the corresponding low-energy effective field theory (EFT) . This way, one arrives at the possible sequential generation (in general, due to a few sequential symmetry breakings at the high-energy scales) of the relevant mass (or SM Higgs Yukawa) terms in the SM, whose values are matched to zero or to a universal non-zero value in the high-scale limit of a more symmetric, and hence more fundamental theory. The search for such ultraviolet (UV) completions, possessing a low-energy SM-like EFT, and explorations of their vast potential for explaining the origin of the SM structure have only began recently [22,23]. A model of radiatively generated fermion masses in the SM via sequential loop suppression has been proposed in Ref. [22]. In this model, the top quark mass is generated at tree level, while the bottom, tau and muon lepton masses arise at one-loop level. Meanwhile, the smaller up, down and strange quark and electron masses are generated at two-loop level, while the light active neutrinos acquire their masses at four-loop level. However, such a natural sequential generation of fermion masses at various orders of the Perturbation Theory comes with a price. Namely, this model, based on the SM gauge symmetry, supplemented with the S 3 × Z 2 discrete group, has a quite low cut-off, since it incorporates non-renormalizable Yukawa terms. In addition, it has another drawback in that the S 3 × Z 2 discrete group is softly broken, yielding an unknown UV completion of the theory. This situation may indicate the need for horizontal continuous symmetries for a consistent description of hierarchies in the matter sectors in a fully renormalizable way. As a follow up to this study, in the later work of Ref. [23], we proposed such a renormalizable model based on the SU 3C × SU 3L × U 1X × U 1L × Z 2 × Z 4 symmetry, which generates the SM fermion mass hierarchies with an emergent sequential loop suppression mechanism. However, that model still retains some drawbacks of the previous formulation. Namely, the Z 2 ×Z 4 symmetry is softly broken, and the scalar and fermion sectors are excessively large, making it very difficult to perform a reliable phenomenological analysis. In addition, that model does not explain the R K and R K * anomalies, recently observed by the LHCb experiment [28][29][30], since it treats the first and second lepton families in the same footing. Furthermore, in the model of Ref. [23], the masses for the light active neutrinos appear at two-loop level, as well as the masses for the electron and for the light up, down and strange quarks. Thus the smallness of the light active neutrino masses with respect to the electron mass does not receive a natural explanation in the model. As a natural step forward, it would be instructive to find a new analogous formulation that enables us to generate the SM charged fermion mass hierarchies via a sequential loop suppression mechanism and to generate three-loop level light active neutrino masses without the inclusion of soft breaking mass terms. It would also be relevant to further explore the potential of such formulations for explaining the LHCb anomalies. In the current work, we propose a first renormalizable model, an extended variant of the Inert Higgs Doublet model (IDM) [24], that enables to generate strong fermion hierarchies via another sequential loop suppression pattern, not yet discussed in the literature, without introducing any soft family symmetry breaking mass terms. Similarly to the previous formulations, in the current model the top quark and exotic fermions do acquire tree level masses, whereas the masses of the remaining SM fermions are radiatively generated. Namely, the masses for the bottom, strange and charm quarks, tau and muon leptons are generated at one-loop level, whereas the masses for the up and down quarks, as well as the electron mass, arise at two-loop level. In variance to the previous version in Ref. [22], the light active neutrinos acquire masses via radiative seesaw mechanisms at three-loop level, whereas in the former the light active neutrino masses were induced at two-loop level. Finally, the minimal field content of the model is not as complicated as in the previous formulations, enabling us to explore several key phenomenological implications of this model, which is the main subject of a follow up study. The current article is organized as follows. In section II we discuss generic conditions for a sequential loop suppression mechanism, providing a motivation for the proposed model. In section III, we set up the formalism for the extended IDM, containing the basic details about the symmetries, particle content, Yukawa interactions and scalar potential crucial for the implementation of the sequential loop suppression mechanism. The scalar mass spectrum is discussed in detail in section IV. In section V we discuss the implications of our model for the radiative generation of the quark mass and mixing hierarchies. The charged lepton and neutrino mass spectra and the corresponding mixing patterns are discussed in section VI. Finally, our conclusions are briefly stated in section VII.

II. SEQUENTIAL LOOP SUPPRESSION MECHANISM
Before describing our model in detail, let us first explain the reasoning behind introducing the additional scalar and fermion degrees of freedom and the symmetries that are required for a consistent implementation of the sequential loop suppression mechanism for generating the SM fermion hierarchies.

A. Quark sector
First of all, it is worth noticing that the top quark mass can be generated at tree level by means of a renormalizable Yukawa operator, with the corresponding coupling of order one, i.e.
where φ 1 is a SU 2L scalar doublet. To generate the charm quark mass at one loop level, it is necessary to forbid the operator: at tree level and to allow other operators instead, which are described in what follows. Obviously, in this case the charm quark mass is not generated in the same way as the top quark mass, i.e., from a renormalizable Yukawa term, since this would imply setting the corresponding tree-level Yukawa coupling unnaturally small. In order to generate a small charm quark mass at one-loop level, we will use the following operators: For this to happen, we need to extend the SM gauge symmetry by adding the U 1X × Z 4 symmetry, where the U 1X is a spontaneously broken gauge symmetry, and Z 4 is a preserved (exact) discrete symmetry (we will explain below why one should use Z 4 instead of Z 2 ) under which the extra φ 2 scalar doublet, the σ 2 electrically neutral scalar singlet, and the q jL , u jR (j = 1, 2, 3) quark fields are nontrivially charged. Furthermore, the SU 2L singlet heavy quarks T L , T R with electric charges equal to 2/3 have to be added to the fermion spectrum in order to implement the one-loop radiative seesaw mechanism that gives rise to the charm quark mass. Besides, an electrically neutral weak-singlet scalar σ 1 is needed in the scalar spectrum to provide a tree-level mass for the exotic T quark and to close the one-loop Feynman diagram that generates the charm quark mass. Let us note that σ 1 is the only scalar neutral under the unbroken Z 4 symmetry, and thus it conveniently acquires a vacuum expectation value (VEV) that breaks the U 1X gauge symmetry.
In addition, we assume that the q nL (n = 1, 2) fields have U 1X charges that are different from the charge of q 3L , while the SU 2L scalar doublets φ 1 and φ 2 have different U 1X charges as well. Thus, the third row of the up-type quark mass matrix is generated at tree level, whereas the first and second row emerge at one-loop level. Note that, since there is only one heavy exotic T quark mediating the one-loop radiative seesaw mechanism that generates the first and second row of the up-type quark mass matrix, the determinant of this matrix is equal to zero. Therefore, the up quark is massless at one-loop level, and in order to generate an up quark mass at two-loop level, the following operators are required: where the extra SU 2L singlet heavy quarks T 1L , T 1R have electric charges equal to 2/3, and ρ 1 is another electrically neutral weak-singlet scalar, charged under the unbroken Z 4 symmetry. Let us explain the reasons for choosing the exact Z 4 discrete symmetry instead of the Z 2 symmetry. Despite the fact that the Z 2 symmetry is sufficient for implementing the one-loop radiative seesaw mechanism that generates the charm quark mass, this symmetry is not enough for generating the up quark mass at two-loop level. Since the σ 2 and ρ 1 scalar fields are charged under an unbroken discrete symmetry, the invariance of the T L ρ 1 T 1R and ρ 2 1 σ 1 σ * 2 operators requires the use of the Z 4 symmetry instead of Z 2 . The operators given above enable to generate the two-loop contributions to the first and second rows of the up-type quark mass matrix, and these contributions yield a nonvanishing determinant for the up-type quark mass matrix, giving rise to a suppressed two-loop up quark mass. Turning now to a possible bottom quark mass generation at one-loop level, the following operators should be forbidden by means of, for example, the U 1X gauge symmetry. The fermion spectrum has to be extended by including the additional SU 2L -singlet heavy quarks B nL , B nR (n = 1, 2) with electric charges equal to −1/3, so that the mass for the bottom quark is generated with the help of the following operators: Note that we have added two SU 2L -singlet heavy quarks B n (n = 1, 2) instead of one, in order to fulfil the anomaly cancellation conditions discussed below. In addition, in order to generate the first and second rows of the down-type quark mass matrix at one-loop level via a radiative seesaw mechanism, we need to incorporate, besides the SU 2L -singlet heavy quarks T L , T R previously introduced in the up-type quark sector, also an electrically charged weak-singlet scalar, namely, ϕ + 2 . Thus, the first and second rows of the down-type quark mass matrix are generated by means of another set of operators given by Furthermore, in order to avoid tree-level masses for the down and strange quarks one has to forbid the terms q nL φ 1 d jR , n = 1, 2, j = 1, 2, 3 .
The latter can be excluded by assigning q nL (n = 1, 2) to be odd, while d jR (j = 1, 2, 3) to be even under the aforementioned exact Z 4 symmetry. Furthermore, since φ 2 and σ 2 have the same Z 4 charge equal to −1, and in order for the operators B nL σ * 1 B mR to be invariant under the Z 4 symmetry, the SU 2L -singlet heavy quarks B nL , B nR (n = 1, 2) should also have Z 4 charge equal to −1 (in non-additive notation). In fact, this condition implies that the q 3L and u 3R quarks should be neutral under Z 4 . Since there is only one fermionic seesaw mediator, i.e., the SU 2L singlet heavy quark T needed to generate the first and second rows of the down-type quark mass matrix at one-loop level, a nonvanishing one-loop strange quark mass emerges, whereas the down quark remains massless at this point. Consequently, similarly to the up-type quark sector, the two-loop contributions to the first and second rows of the down-type quark mass matrix need to be generated in order to give rise to a down-type quark mass at two-loop level. For that purpose, the following operators are required: where an extra electrically charged weak-singlet scalar, ϕ + 1 , as well as an electrically neutral weak-singlet scalar, ρ 2 , have been added to the scalar spectrum. Furthermore, the fermion spectrum has been extended by means of extra SU 2L -singlet heavy quarks T 2L , T 2R with electric charges equal to 2/3. The two-loop contributions to the first and second rows of the down-type quark mass matrix provide its nonvanishing determinant, giving rise to a two-loop down quark mass.

B. Charged lepton sector
Now, consider the sequential loop suppression mechanism capable of explaining the observed hierarchy between the SM charged lepton masses. In what follows, let us discuss a possible way of generating the one-loop tau and muon masses as well as a two-loop electron mass. First of all, the following operators have to be forbidden by using the U 1X gauge symmetry, so that the SM charged lepton mass matrix is generated at one-loop level by means of the terms where weak-singlet charged leptons E jL , E jR (j = 2, 3) have been included in the fermion spectrum. Let us note that these fields mediate the one-loop radiative seesaw mechanism that generates the (1, 1), (2, 2), (3, 3), (1, 3) and (3, 1) entries of the charged lepton mass matrix. Consequently, at this point the determinant of the charged lepton mass matrix is equal to zero and thus only the tau and muon leptons appear to be massive at one-loop level, whereas the electron remains massless. Two-loop corrections to the (1, 1), (3, 3), (1, 3) and (3, 1) entries are needed in order to induce a non-zero electron mass at two-loop level, and an extra weak-singlet charged lepton E 1 would be required for this purpose. To obtain such two-loop corrections, the following operators should be considered where an extra electrically neutral weak-singlet scalar ρ 3 have been introduced to the scalar spectrum. As soon as the two-loop corrections generated by the above operators are included, the determinant of the charged lepton mass matrix becomes nonzero, thus giving rise to a two-loop electron mass, as expected. Let us note that in this setup there is a mixing between the SM charged leptons of the first and third generations, but they do not mix with the second generation, which is a consequence of their U 1X assignments, as will be shown below. The second generation left-handed lepton SU 2L -doublet has a nonvanishing U 1X charge, whereas the first and third generations are not charged under U 1X . In addition, the first and third generations of SM right-handed charged leptons should have the same U 1X charge, which is different from the corresponding charge of the second generation.

C. Light active neutrino sector
Turning to the neutrino sector, in order to generate the SM light active neutrino masses at three-loop level as well as a realistic lepton mixing, a few operators have to be forbidden, namely, while the following operators are required .
where ν iR (i = 1, 2, 3), Ω R and Ψ R are the SM-singlet right-handed Majorana neutrinos, and η, σ 3 and ρ 4 are the extra SM-singlet scalars. The latter fields have to be added in order to ensure three-loop level generation of the SM light active neutrino masses, as well as the lepton mixing parameters sin 2 θ 23 and sin 2 θ 12 . As was mentioned above, the charged lepton mass matrix has a mixing only in the (1,3)-plane such that the lepton mixing parameters sin 2 θ 23 and sin 2 θ 12 emerge from the neutrino sector. Let us note that the U 1X gauge symmetry prevents the light active neutrino masses to be generated at tree level, whereas the Z 2 symmetry, together with the U 1X gauge symmetry, forbid the mixing terms between the right-handed Majorana neutrinos ν kR (k = 1, 3) and ν 2R . This enables us to avoid the appearance of light neutrino masses at one-and two-loop levels. Note also that the Z 2 symmetry, as well as the U 1X gauge symmetry, are crucial to forbid the terms ν 2R σ 1 ν C kR and ν 2R ρ 5 ν C kR (k = 1, 3) that could result in the appearance of SM neutrino masses at one loop.

III. THE EXTENDED IDM MODEL
In this section, we will summarize the main features of the first renormalizable model, an extended variant of the Inert Higgs Doublet model (IDM), that includes a sequential loop suppression mechanism for the generation of the SM fermion mass hierarchies, without the inclusion of soft breaking mass terms and, at the same time, allowing for an explanation of the R K and R K * anomalies, thanks to the non-universal U 1X assignments of the fermionic fields that yield non-universal Z couplings to fermions. In a forthcoming work, we will show in detail how our model can fit the R K and R K * anomalies. As previously stated in the Introduction section, we emphasize that our model, apart from having all the means for explaining the R K and R K * anomalies not previously addressed in the model of Ref. [22], has a more natural explanation for the smallness of the light active neutrino masses than the one provided in Ref. [22], since in the former the masses for the light active neutrinos are generated at three-loop level, whereas in the latter they appear at two loops. Furthermore, unlike the model of Ref. [22], our current model does not include soft breaking mass terms. Let us summarize the structure of a minimal model where the sequential loop suppression mechanism capable of radiative generation of the mass and mixing hierarchies in the SM fermion sectors is realized. The reasons for choosing a particular field content and symmetries have been outlined in the previous section, and will be further substantiated below.
The corresponding Z 2 × Z 4 charges are given by where ω = e 2πi 3 . It is worth noticing that the model does not contain a weak-singlet scalar field with the following simultaneous three features: charged under the Z 2 discrete symmetry, neutral under the unbroken Z 4 symmetry with U 1X charge equal to ±1. Let us note that the SM-type Higgs doublet, i.e., φ 1 , as well as the electroweak-singlet scalars σ 1 and ρ 5 are the only scalar fields neutral under the exact Z 4 discrete symmetry. Since the Z 4 symmetry is preserved, the Higgs doublet φ 1 and the singlets σ 1 and ρ 5 are the only scalar fields that acquire nonvanishing VEVs. The VEV in σ 1 is required to spontaneously break the U 1X local symmetry, whereas the ρ 5 VEV spontaneously breaks the Z 2 discrete symmetry, due to its nontrivial Z 2 charge. Note that the exact Z 4 discrete symmetry guarantees the presence of several stable scalar dark matter candidates in our model. These are represented by the neutral components of the inert SU 2L scalar doublet φ 2 , as well as by the real and imaginary parts of the SM-singlet scalars σ 2 , σ 3 , η and ρ j (j = 1, 2, 3, 4). The set of SU 2L -singlet heavy quarks T L , T R , B kL , B kR (k = 1, 2) represents the minimal amount of exotic quark degrees of freedom needed to implement the one-loop radiative seesaw mechanism that gives rise to the charm, bottom and strange quark masses. Furthermore, in order to ensure the radiative seesaw mechanism responsible for the generation of the up and down quark masses at two-loop level, the SU 2L singlet heavy quarks T kL , T kR (k = 1, 2), as well as the electrically neutral, ρ 1 , ρ 2 , and electrically charged, ϕ + 1 scalar SU 2L -singlets should also be present in the particle spectrum. To summarize, the SM fermion sector of the considered model includes a total of six electrically charged weak-singlet leptons E jL and E jR (j = 1, 2, 3), four right-handed neutrinos ν jR (j = 1, 2, 3), Ω R , and ten SU 2L -singlet heavy quarks T L , T R , T kL , T kR , B kL , B kR (k = 1, 2). It is assumed that the heavy exotic T , T k and B k quarks have electric charges equal to 2/3 and −1/3, respectively. More specifically, the quark sector of the extended IDM under consideration has the following SU 3c ×SU 2L ×U 1Y ×U 1X charges while their Z 2 × Z 4 charge assignments read The radiative seesaw mechanism that generates the charged lepton mass hierarchy is similar to the one that produces the SM down-type quark mass hierarchy. The generation of one-loop tau and muon masses as well as two-loop electron mass is mediated by the electrically charged weak-singlet leptons E jL and E jR (j = 1, 2, 3), by the inert scalar SU 2L -doublet, φ 2 , and by the SU 2L -singlets σ 2 , ρ j (j = 1, 2, 3). Moreover, the three-loop radiative seesaw mechanism responsible for the generation of the light active neutrino masses is mediated by the right-handed neutrinos ν jR (j = 1, 2, 3), Ω R , as well as by the inert scalar SU 2L doublet φ 2 and the SU 2L -singlet σ 2 . To avoid tree-level mixing between the right-handed Majorana neutrinos ν kR (k = 1, 3) and ν 2R triggered by Yukawa interactions with σ 1 , we need to impose a nontrivial Z 2 charge of ν 2R while keeping ν kR (k = 1, 3) Z 2 -neutral.
In particular, the SU 3c × SU 2L × U 1Y × U 1X charges of the leptonic and neutrino fields of the model are defined as follows whereas the corresponding Z 2 × Z 4 charges are given by The U 1X × Z 2 × Z 4 symmetry and the particular assignments listed above are crucial for avoiding the appearance of SM light active neutrino masses at one-and two-loop levels.
Let us note that the left-handed quark SU 2L doublets of the first and second generations are distinguished from the third generation by means of the U 1X charge assignments. In addition, the local U 1X family symmetry distinguishes the second generation left-handed lepton SU 2L doublet from the first and third generation ones. Such nonuniversal U 1X charge assignments in the fermion sector are crucial for implementing the sequential loop suppression mechanism and, hence, for the induced strong hierarchies in the SM fermion mass spectrum. Besides, as will be explicitly demonstrated in a forthcoming paper, such assignments are also relevant for explaining the R K and R K * anomalies. With the above assignments, we have numerically checked that the gauge anomaly cancellation conditions are satisfied in our model.

B. Yukawa interactions
With the above specified particle content and charge assignments, the most general renormalizable Lagrangian of Yukawa interactions and the exotic fermion mass terms, invariant under the symmetry, takes the following form where the Yukawa couplings are O(1) parameters. From the quark Yukawa terms it follows that the top quark mass emerges due to an interaction involving the SM-like Higgs doublet φ 1 only. After spontaneous breaking of the electroweak symmetry, the observed hierarchies of SM fermion masses arise by means of a sequential loop suppression, according to the following pattern: tree-level top quark mass; one-loop bottom, strange, charm, tau and muon masses; two-loop masses for the up, down quarks as well as for the electron. Furthermore, the SM light active neutrinos get their masses by means of a three-loop radiative seesaw mechanism.

C. Scalar potential
The most general renormalizable scalar potential invariant under the gauge and discrete symmetries of our model is given by From the minimization conditions for this potential, we find the following simple relations that will be used below in a discussion of the scalar mass spectrum of the model.

IV. SCALAR MASS SPECTRUM
Considering the scalar potential given above, we find that the squared mass matrices for the CP-even neutral scalar sector are have the following form where M (1) CPeven and M (2) CPeven are the squared mass matrices for the Z 4 -neutral and Z 4 -charged scalars, respectively. The matrix M (1) CPeven in the basis Re φ 0 1 , Re (σ 1 ) , ρ 5 (remind, ρ 5 is a real SM-singlet scalar), takes the form: The second mass form M (2) CPeven in the basis Re ( Since the 126 GeV SM-like Higgs boson is found in the squared mass matrix M CPeven , and considering the fact that the scalar potential has a very large number of parameters, in this first study it is sufficient to diagonalize only M (1) CPeven in the scalar sector. In addition, since this matrix cannot be diagonalized in analytically closed form, and for the sake of simplicity, here we focus on a particular scenario with v σ = v ρ . In this scenario, the matrix M Consequently, the physical scalar states contained in the matrix M CPeven are given by: where h is the 126 GeV SM-like Higgs boson, whereas χ 1 and χ 2 are the physical heavy scalar fields, which acquire masses at the scale of U 1X breaking. The squared masses of these fields are given by Furthermore, we find that the SM-like Higgs boson h has the couplings that are very close to the SM expectation, with small deviations of the order of ∼ v 2 /v 2 σ . Considering the CP-odd neutral scalar sector, we find that the squared mass matrices for the electrically neutral CPodd scalars in the basis, are Im φ 0 are given by CPodd and M CPodd are the squared mass matrices for the CP-odd scalars, neutral and charged under Z 4 , respectively. Note that the squared mass matrix M (1) CPodd (which is written in the basis Im φ 0 1 , Im (σ 1 ) ) is exactly zero, since Im φ 0 1 and Im (σ 1 ) are the Goldstone bosons associated with the longitudinal components of the Z and Z gauge bosons, respectively. Finally, the squared mass matrix for the charged scalar fields in the basis φ such that φ ± 1 are the electrically charged massless scalar states corresponding to the Goldstone bosons associated with the longitudinal components of the W ± gauge bosons.

V. RADIATIVELY GENERATED QUARK MASS AND MIXING HIERARCHIES
The quark Yukawa interactions determined by Eq. (21) give rise to the following up and down mass matrices for the SM quarks, respectively, where the dimensionless parameters ε  Fig. 1.
In what follows, we demonstrate that the mass matrices for SM quarks given above incorporate the observed hierarchies in the SM quark mass spectrum and the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix. To this end, we proceed with a parametrization of the SM quark mass matrices in the following form where l ≈ (1/4π) 2 ≈ 2.0 × λ 4 is the loop suppression factor, and λ = 0.225 is the Wolfenstein parameter. As a consequence, we expect that a We remark that the Feynman diagrams contributing to the entries of the SM fermion mass matrices contain a large number of uncorrelated parameters that belong to the fermion and scalar sectors of our model. Nevertheless, these parameters can be absorbed into a limited number of effective parameters ε ij (n, m = 1, 2 and i, j = 1, 2, 3), which can be used to reproduce the experimental values of the physical observables in the quark sector [33]  Here, m t,u,c,d,s,b are the SM quark masses, θ 12 , θ 23 , θ 13 are the mixing angles, and J is the Jarlskog parameter. While our model does not predict the exact values of the physical observables, it offers a natural explanation of the observed (strong) hierarchies. As was previously mentioned, it only pretends to reproduce the existing pattern of quark masses and mixing caused by a sequential loop suppression predicted by the model. 1.771 .
The above O(1) values exactly reproduce the measured central values of the SM quark masses and CKM parameters given in Eq. (34). Hence, our model is consistent with and successfully reproduces the existing pattern of SM quark masses caused by the sequential loop suppression mechanism, with different quark flavors getting mass at different orders in Perturbation Theory as discussed above.

VI. RADIATIVELY GENERATED LEPTON MASSES AND MIXINGS
The lepton and neutrino Yukawa interactions and mass terms given in Eq. (21) give rise to the characteristic Feynman loop diagrams illustrated in Figs. 2 and 3 that necessarily generate the following SM charged lepton and light active neutrino mass forms: ii , ε (l) 13 , ε (l) 31 are the dimensionless parameters generated at one-loop level, whereas the parameters ε 13 and ε (d) 31 appear at two-loop level. In what follows, we show that the lepton mass matrices given above can accommodate the experimental data on the SM lepton masses and mixing. For this purpose, following the same strategy as for the quark mass forms discussed in the previous section, we parametrize the SM charged lepton mass matrix as follows: In order to fit the measured values of the charged lepton masses, as well as the neutrino mass squared differences and lepton mixing parameters [32], we proceed by solving the eigenvalue problem for the SM lepton and light neutrino   (M ν ) ij O(10 −2 ) (i, j = 1, 2, 3), the experimental values for the physical observables of the lepton sector, i.e., the three charged lepton masses, the two neutrino mass squared splittings, the three leptonic mixing parameters and the leptonic Dirac CP violating phase can be successfully reproduced. Consequently, our model is consistent with and successfully reproduces the existing pattern of SM charged lepton masses generated by the sequential loop suppression mechanism.

VII. DISCUSSIONS AND CONCLUSIONS
We have constructed a first renormalizable extension of the Inert Doublet Model that enables an implementation of a sequential loop suppression mechanism, capable of explaining the observed SM fermion mass hierarchy without invoking soft breaking mass terms. In our model, the SM gauge symmetry is supplemented by the U 1X ×Z 2 ×Z 4 family symmetry, where the gauge U 1X and discrete Z 2 symmetries are spontaneously broken, whereas the Z 4 symmetry is preserved. Our model is consistent with the observed SM fermion mass spectrum and fermionic mixing parameters and allows for an explanation of the recently observed R K and R K * anomalies, thanks to the non-universal Z couplings to fermions.
We focused on an extension of the Inert Higgs Doublet model (IDM) that allows the implementation of the sequential loop suppression mechanism for the generation of SM fermion masses instead of an extension of the inert 3-3-1 model (model based on the SU 3C × SU 3L × U 1X gauge symmetry). As previously mentioned, the extension of the inert 3-3-1 model of Ref [23] does not explain the R K and R K * anomalies and the light active neutrino masses appear at two-loop level like the masses of the light SM charged fermions. Addressing the R K and R K * anomalies in the framework of a 3-3-1 model would require to consider five families of SU (3) L leptonic triplets as done in Ref. [31], in order to have different U (1) X charge assignments for the first and second lepton families, without spoiling the anomaly cancellation conditions. Thus, modifying the inert 3-3-1 model of Ref [23] to account for the R K and R K * anomalies, and to generate the hierarchy of SM fermion masses by sequential loop supression mechanism, with the light active neutrino masses appearing at three-loop level, will require a much larger particle content than the one adopted in the framework of an extended IDM. In our model only the top quark and exotic fermions acquire tree-level masses, whereas the masses of the remaining SM fermions emerge from a radiative seesaw-like mechanism: the masses for the bottom, strange and charm quarks, tau and muon leptons are generated at one-loop level, whereas the masses for the up and down quarks as well as the electron mass appear at two-loop level. Furthermore, light active neutrino acquire masses by means of a radiative seesaw mechanism at three-loop level. Due to an unbroken Z 4 discrete symmetry, our model has several stable scalar dark matter candidates, which can be the neutral components of the inert SU 2L scalar doublet φ 2 as well as the real and imaginary parts of the SM scalar singlets σ 2 , σ 3 , η and ρ j (j = 1, 2, 3, 4). Furthermore, the model can have a fermionic dark matter candidate which is the only SM-singlet Majorana neutrino Ω 1R with a non-trivial Z 4 charge. A study of the phenomenological implications of our model goes beyond the scope of the present paper and will be performed in a forthcoming work.