A $\mu$-$\tau$-philic scalar doublet under $Z_n$ flavor symmetry

We propose a minimal model which accommodates the long-standing anomaly of muon magnetic moment based on abelian discrete flavor symmetries. The standard model is extended by scalar doublets charged under a $Z_n$ lepton flavor symmetry. In these models, a large contribution to the muon magnetic moment can be obtained by the chirality enhancement from new scalar mediated diagrams without conflicting with the flavor symmetry. Thanks to the lepton flavor symmetry, these models automatically forbid lepton flavor violation. The minimal model is based on $Z_4$ symmetry with only one extra scalar doublet. In this model, we show that the parameter space favored by the muon $g-2$ can easily be consistent with experimental constraints and theoretical bounds such as the electroweak precision tests, lepton universality, potential stability condition and triviality bound as well as the LHC direct search mass bound. The new contributions to the muon electric dipole moment and the Higgs decay into $\gamma\gamma$ can be indirect signals of the model.


Introduction
The Standard Model (SM) has been established by the discovery of the Higgs boson at the LHC. New particles beyond the SM are also being searched at the LHC. However, there is no signature of new particles until now, and the experimental results are consistent with the SM predictions. Other than the high energy frontier experiment, many of flavor observables are measured very precisely as the luminosity frontier experiment. A striking indication of the beyond the SM would be the muon anomalous magnetic moment (muon g − 2). There is a discrepancy between the measured value and the SM prediction as [1] ∆a µ = a exp µ − a SM µ = 268(63)(43) × 10 −11 , where the numbers in the first and second parentheses represent the statistical and systematic errors, respectively. The total significance of the deviation is 3.5σ far from the SM prediction. 1 Note that there is a non-negligible large theoretical uncertainties in the hadronic contribution due to the light-by-light scattering [3]. Currently, FNAL E989 experiment is ongoing, and will achieve a factor four improvement on its precision at the end of the running [4].
There are many attempts to explain the discrepancy of the muon magnetic moment. For instance, in lepton-specific two Higgs doublet models (THDMs), the new contribution to the muon magnetic moment due to the additional Higgs bosons can be enhanced by a large tan β, which is the ratio of vacuum expectation values (VEVs) of two Higgs doublets [5][6][7][8]. The THDM with tree-level flavor changing neutral currents has also been studied to explain the discrepancy in the light of the h → µτ excess at the LHC [9,10], which has been disappeared. Another way is to consider a light Z ′ gauge boson associated with an extra U(1) Lµ−Lτ symmetry [11,12], or a light hidden photon [13]. In these models, thanks to the new light mediator running in the loop diagram, the muon magnetic moment can be enhanced even with a smaller coupling strength. There is also argument to account for the discrepancy in framework of supersymmetry [14], axion-like particle [15] and fourth generation of leptons [16].
In this paper, we propose a minimal model explaining the discrepancy of the muon anomalous magnetic moment between the SM prediction and the measurement based on abelian discrete symmetries. Models based on the abelian discrete groups easily give a sufficiently large and the correct sign of the contribution to the muon magnetic moment. The model based on a Z 4 symmetry is identified as the minimal model, which is a kind of variant of the inert scalar doublet model based on a Z 2 symmetry. Thanks to the Z 4 symmetry, the lepton flavor violating (LFV) processes such as ℓ → ℓ ′ γ (ℓ, ℓ ′ = e, µ, τ ) are forbidden automatically against severe bounds of their non-observation. As a result, we find a solution to the muon g − 2 anomaly without conflicting with the constraints from the electroweak precision tests and the lepton universality of heavy charged lepton decays Table 1: Particle contents of models based on U(1) Lµ−Lτ and Z n flavor symmetries. The quantum numbers of the SM are also shown in the notation of SU(3) c , SU(2) L U (1) Y . For abelian discrete symmetry Z n , ω is a conjugate of ω, where ω is n-th root of unity. and Z boson leptonic decays. In addition, we examine whether the model is consistent with theoretical bounds of potential stability and triviality. We will formulate analytic expressions of these quantities, and numerically explore the parameter space which can accommodate the discrepancy of the muon anomalous magnetic moment. As further perspective, neutrino mass generation mechanism and a distinctive collider signature, a prediction for muon electric dipole moment induced by new CP phases and influence on the Higgs decay into γγ will also be discussed.

Flavor Charged Scalar Doublets
Let us discuss a simple extension of the SM with a pair of scalar doublets (Φ, Φ ) whose global U(1) Lµ−Lτ flavor charge is (2, −2), where L µ and L τ represent the muon and tau lepton flavor numbers, respectively. Detailed quantum charge assignments are given in Table 1. Under this flavor symmetry, the following new Yukawa interactions are allowed, in addition to the quartic scalar interaction term (H † Φ)(H † Φ). These interactions easily generate sizable contributions to the muon g − 2 by the scalar mediators as shown in Fig. 1. In the ordinary gauged U(1) Lµ−Lτ model, the discrepancy in the muon g − 2 is explained by the new light Z ′ gauge boson [11,12], while in our new proposals a pair of scalar doublets is introduced to give a sizable contribution to the muon g − 2. A similar contribution to the muon g − 2 from the scalar doublets are discussed in the model based on the SU(2) µτ symmetry, which contains the U(1) Lµ−Lτ symmetry as a subgroup [17]. In such cases, a pair of scalar doublets plays the primary role in explaining the muon g −2 anomaly instead of Z ′ bosons. We noted that this new contribution remains even with the unbroken U(1) Lµ−Lτ flavor symmetry limit. From the above consideration in our mind, we begin with a global U(1) Lµ−Lτ symmetry together with a pair of scalar doublets as a simple model for the muon g − 2 anomaly. On the other hand, the U(1) Lµ−Lτ symmetry must be broken in order to realize observed Φ Φ µ τ τ µ γ H H Figure 1: Feynman diagram inducing muon anomalous magnetic moment in U(1) Lµ−Lτ and Z n (n = 2, 3, · · · ) models.
The Yukawa interactions based on Z n flavor symmetries are given by The Z n charge assignment in each model is given in Table 1. For n ≥ 5, an accidental global U(1) Lµ−Lτ symmetry is recovered in the Yukawa interactions taking into account renormalizability. Depending on the chosen abelian discrete flavor symmetry, a specific structure of the Yukawa interaction is predicted. Note that since Φ is identical to Φ in the Z 2 and Z 4 models, the Yukawa interactions of Φ are not shown for these models. In the following, we focus on the models with only one extra scalar doublet Φ, which minimally explain the muon g − 2 anomaly. In the model based on Z 2 or Z 4 , possible large new contributions to the muon g − 2 are retained thanks to the existence of the quartic term (H † Φ) 2 . From the view of exper-imental constraints, the Z 4 model is more favorable because the Z 2 model predicts LFV processes τ → 3µ, eµµ at tree-level, and thus parameter tuning is necessary to suppress these processes. On the other hand, the LFV processes are automatically forbidden in the (unbroken) Z 4 model. From the view of the numbers of parameters in the model, again the Z 4 model is preferable both in the Yukawa sector and the scalar potential. We therefore conclude that the model based on the Z 4 lepton flavor symmetry is the minimal scalar extension of the SM to accommodate the muon g − 2 anomaly.

The Minimal Model for Muon g − 2
Following the argument in the previous section, we introduce a new scalar doublet Φ to the SM, and impose a Z 4 symmetry. The Z 4 charge assignment is shown in Table 1, and all the other fields are trivial under the Z 4 symmetry. The invariant scalar potential is given by This scalar potential is the same as that in the scalar inert doublet model [19], where an exact Z 2 symmetry is preserved in the potential. In general, the quartic coupling λ 5 and the Yukawa couplings y µτ , y τ µ are complex. One of the CP phases can be eliminated by the field redefinition of Φ. Here, we remove the CP phase of λ 5 without loss of generality.
Since we demand a stable vacuum, the potential should be bounded from below. The conditions for these requirements are known as [20] at tree level. The Higgs doublet H develops a VEV as in the SM, and the electroweak symmetry is spontaneously broken. The new doublet scalar Φ is assumed to have a vanishing VEV at leading order. The scalar fields can then be parameterized as A component field h corresponds to the Higgs boson with the mass m h = √ 2λ 1 v = 125 GeV. The electrically neutral component of Φ, φ 0 = (ρ + iη)/ √ 2, splits into the two mass eigenstates ρ and η. The masses of these neutral states and charged component φ + are given by Thus, one can see that the mass splitting between ρ and η is controlled by the quartic In this model, the new contribution to muon anomalous magnetic moment comes from Fig. 2, which is computed as where the loop functions I 1 (a, b) and I 2 (a, b) are defined by Note that the contribution in the first line of Eq. (13) is dominant compared to that in the second line with an enhancement factor m τ /m µ ≈ 17, because of the chirality flipping effect. The numerical value of these loop functions are always positive, and thus the sign of the new contribution is determined by the relative sign of Re (y µτ y τ µ ) and In the numerical analysis, we require that the discrepancy of muon g −2 is improved to be within 2σ range after including the new physics contribution. Thus, ∆a new µ in Eq. (13) should be in the interval [17] 115 4 The Constraints

Electroweak Precision Tests
The new scalar particles ρ, η and φ + affect the electroweak precision observables through vacuum polarization diagrams. These are conveniently parameterized by the socalled S, T, U-parameters [21]. The expression of the S, T, U-parameters in this model is the same as that in the inert doublet model [22] or in the THDM [23,24] with the alignment limit, which are given by where G F is the Fermi constant, α em is the electromagnetic fine structure constant, and the functions F (x, y) and G(x, y) are given by The current experimental bounds on these parameters are summarized as [25] S = 0.05 ± 0.11, T = 0.09 ± 0.13, U = 0.01 ± 0.11, with correlation coefficients 0.90 between S and T , −0.59 between S and U, and −0.83 between T and U, respectively. We impose the requirement that the theoretical prediction on these parameters should be kept in the 2σ range of the experimental values. If relatively light new particles ( m W ) are mediated in a loop, more sophisticated analysis of the electroweak precision tests may be applied as in a lepton-specific THDM [5].

Lepton Universality in Charged Lepton Decays
The new Yukawa couplings y µτ and y τ µ give additional contributions to the decay of charged leptons. First, the new tau decay mode τ → µ ν τ ν µ is induced at tree level. The partial decay width is calculated as Figure 3: Feynman diagrams of the loop corrections to charged lepton currents is the W -boson propagator correction, and r τ γ is the QED radiative correction, which are given by [26] In the numerical evaluation, we use PDG data for the W boson mass, charged lepton masses, the electromagnetic fine structure constant [1]. If we worked out in the mass eigenbasis of neutrinos, one may expect interference effect in τ → µ ν i ν j (i = 1, 2, 3). Such effect is, however, negligible since the chirality flip occurs and it is suppressed by small neutrino masses. Second, one-loop corrections in the charged lepton currents are induced by new Yukawa interactions as shown in Fig. 3. Although each diagram includes a divergence, it cancels out after the sum over all the graphs. We then obtain a finite correction without renormalization. Following the results in Ref. [7] for the Type-X THDM, we define the loop corrections δg W ν ℓ ℓ as g → g 1 + δg W ν ℓ ℓ . The results in the µ-τ -specific scalar doublet model are where the small lepton masses are neglected, and the loop function I L (x, y) is defined by Taking into account above corrections at tree level and the one-loop level, the total leptonic decay widths of muon and tau lepton are summarized as In Eq. (27), the decay widths for the channels τ → µ ν µ ν τ and τ → µ ν τ ν µ are combined since these processes cannot be distinguished in actual measurements. In general, the above leptonic decay widths are conveniently parameterized as The effective weak couplings for leptons g ℓ (ℓ = e, µ, τ ) are severely constrained as [27] g τ g µ = 1.0011 ± 0.0015, g τ g e = 1.0029 ± 0.0015, g µ g e = 1.0018 ± 0.0014, with correlation coefficients 0.53 between g τ /g µ and g τ /g e , −0.49 between g τ /g µ and g µ /g e , and 0.48 between g τ /g e and g µ /g e , respectively. Using this notation, we find analytic expressions for the corresponding quantities, as g τ g e = 1 + δg W ντ τ 1 + m 4 W |y µτ | 2 |y τ µ | 2 g 4 1 + δg W νµµ g µ g e = 1 + δg W νµµ 1 + m 4 W |y µτ | 2 |y τ µ | 2 g 4 1 + δg W νµµ In the numerical analysis, we demand that these quantities should be in the 2σ range of the experimental values.

Lepton Universality in Z Boson Decays
The Z boson leptonic decays are also modified by the Yukawa couplings. In general, the interactions between Z boson and a pair of charged leptons can be written as where g ℓ L and g ℓ R are given by g ℓ L = −1/2 + sin 2 θ W and g ℓ R = sin 2 θ W at tree level, which are universal over lepton flavors. With this convention, the leptonic decay widths are calculated as Figure 4: Feynman diagrams of the loop corrections to Z boson leptonic decays.
where θ W is the Weinberg angle. The couplings g ℓ L and g ℓ R receive the one-loop corrections from the diagrams shown in Fig. 4. The total one-loop correction is finite while each diagram includes a divergence as same as the case of the charged lepton decay vertices. The loop corrections for the neutral current interaction with tau lepton defined by g τ L/R → g τ L/R + δg τ L/R are parameterized as [6] δg Neglecting the small lepton masses, the coefficients a τ L , b τ L , a τ R and b τ R are computed as where ξ a ≡ m 2 a /m 2 Z (a = φ, ρ, η), and the loop functions B Z (ξ),C Z (ξ) and C Z (ξ 1 , ξ 2 ) are defined by [6] Similarly, the loop correction with muon is obtained by replacing y µτ ↔ y τ µ in Eq. (38)- (41), and there is no loop correction for the neutral current interaction with electron at this order.

Collider Limits
The lower bound of the charged scalar mass is given as m φ 93.5 GeV by LEP [29]. There are also LHC bounds which depend on branching ratio of the charged scalar φ + . Since the charged scalar in our model has the same quantum charges with the charged sleptons in supersymmetric models except for the matter parity, the bound for sleptons from the electroweak production can be applied for φ + if the dominant (prompt) decay channels are φ + → τ ν µ , µ ν τ . The slepton mass bound in the massless neutralino limit can be recast to m φ 700 GeV [30,31]. 2 On the other hand, if φ + is heavier than ρ or η the decay channels φ + → W + ρ, W + η open and can be dominant. In such a case, the mass bound for sleptons cannot be simply applied, and m φ can be lighter than 700 GeV. Thus, we choose m φ = 200 and 700 GeV as representative values in the numerical analysis.

Triviality Bound
Even if the couplings in the model are perturbative at electroweak scale, it may become non-perturbative at a high energy scale after including renormalization group running of the couplings. In particular, if the couplings are O(1) at electroweak scale, it can quickly increase, and tends to become non-perturbative around O(10−100) TeV. The β functions for the renormalization group running at one loop level are collected in Appendix A, where the SM Yukawa couplings are neglected except for the top Yukawa coupling y t . We solve the coupled renormalization group equations from the Z boson mass scale to the cut-off scale Λ. In the numerical analysis, we take Λ = 100 TeV. Then, we demand that all the couplings in the model are perturbative until the cut-off scale. Namely, the required conditions are: |λ i | ≤ 4π (i = 1 − 5) and |y t |, |y µτ |, |y τ µ | ≤ √ 4π at the cut-off scale.

Numerical Analysis
We explore parameter space which can explain the discrepancy in the muon g −2 while satisfying the relevant constraints. In Fig. 5, we present the numerical analysis in the (λ 4 , λ 5 ) plane for fixed values of charged scalar mass (m φ ) and Yukawa couplings (y µτ , y τ µ ). In this subsection, we restrict new Yukawa couplings to be real. The upper (lower) two panels show the low (high) mass scenarios with m φ = 200 (700) GeV. In the left panels, we maximize the new physics contributions to the muon g − 2, where y µτ = y τ µ is assumed, while hierarchical Yukawa couplings are taken in the right panels such that the magnitude of the product y µτ y τ µ is retained as same with the left panels so that the parameter space favored by muon g − 2 does not change. The purple region represents the parameter space which can accommodate the muon magnetic moment anomaly at 2σ confidence level (CL). In the top panels in Fig. 5, the left-top and left-bottom region colored by gray is forbidden because the mass of the neutral scalars ρ or η becomes negative. The green region is ruled out by the electroweak precision tests at 2σ CL. This constraint becomes stronger for lighter scalar masses. On the other hand, even if the charged scalar mass is relatively light, the constraint can be evaded if λ 4 ∼ ±λ 5 , which implies that one of ρ and η is nearly degenerate with the charged scalar φ + . The orange region is excluded by the constraint of the lepton universality (Z boson decays). Since the loop corrections to the Z boson decays given by Eq. (38)-(41) are proportional to |y µτ | 2 or |y τ µ | 2 , one can find that the constraint becomes stronger for larger hierarchy between y µτ and y τ µ for a fixed y µτ y τ µ . Note that the loop corrections for the charged lepton currents given by Eq. (25) also have the same dependence on the Yukawa couplings. However, the constraint from muon and tau lepton decays are slightly weaker than the Z boson decays in the above parameter sets. The red region shows the parameter space that the charged scalar φ + becomes the lightest than the neutral scalars ρ and η. In this region, since the charged scalar decays dominantly into a pair of a charged lepton and a neutrino, the LHC mass limit (m φ 700 GeV) is applied. The outside of the dot-dashed curve colored by gray is disfavored by the potential stability conditions given by Eq (8) and the triviality. The negative λ 4 region tends to be excluded by the potential stability conditions while the remaining region is bounded by the triviality of the quartic couplings λ i (i = 1 − 5). Here, we take λ 2 = λ 3 = 0.5 at the Z boson mass scale as an initial condition of the renormalization group equation. Note that if smaller couplings λ 2 and λ 3 are assumed, the bound of the potential stability becomes stronger as we expect from Eq. (8).
In Fig. 6, we show the parameter space in the (y µτ , y τ µ ) plane by fixing the scalar masses m ρ , m η and m φ . The positive Yukawa coupling y µτ is chosen without loss of generality. We here concentrate on the case with negative values of λ 5 , which is favored by the muon g − 2 anomaly together with a positive value of y τ µ . At the same time, we assume a negative λ 4 + λ 5 for m φ = 200 GeV. This parameter choice allows the cascade decay of the charged scalar to other scalars, and therefore we can avoid the strong constraint on the charged scalar mass from the LHC slepton search. The purple region can accommodate the muon g − 2 anomaly at 2σ CL, while the orange and light blue region are excluded by the lepton universality of the Z boson decays and the charged lepton decays, respectively. The constraint of the electroweak precision tests is satisfied in all the plots, which does not depend on the Yukawa couplings. One can see from Fig. 6 that the constraint of the charged lepton decays (light blue) is always stronger than that of the Z boson decay (orange) when the Yukawa couplings are same order (y µτ ∼ y τ µ ). This is due to the the tree level correction given by Eq. (23). In contrast, when the Yukawa couplings are hierarchical, one of the loop corrections for the charged lepton and Z boson decays becomes stronger. The gray region surrounded by the dot-dashed line shows the bounds of the potential stability. In fact, the potential stability bounds are slightly stronger than the triviality bounds. This is because we take the negative quartic couplings λ 4 and λ 5 at the electroweak scale, and the Yukawa couplings involved in the β λ 4 and β λ 5 make λ 4 and λ 5 further negative at the cut-off scale if the Yukawa couplings are O(1). As a result, it conflicts with Eq. (8) at the cut-off scale. Note that this bound is relaxed if a smaller cut-off scale Λ is assumed.

Neutrino Mass Generation Sector
As we mentioned in the beginning, the Z 4 flavor symmetry in our minimal model must be broken in order to fit the observed data of neutrino masses and mixings. As a simple example for neutrino mass generation, we here consider the type-I seesaw mechanism. A SM singlet scalar S with Z 4 charge ω and a three generation of right-handed neutrinos (N 1R , N 2R , N 3R ) with (1, ω, ω) are introduced to the model. The Lagrangian for the neutrino mass generation sector is whereX = iσ 2 X * (X = H, Φ). The singlet S is assumed to have a VEV ǫ S , which breaks the Z 4 symmetry, where ǫ is introduced to count the order of singlet VEVs. At leading order, O(ǫ 0 ), the (symmetric) neutrino mass matrix has non-zero values only in (1, 1) and (2, 3) elements (see also Fig. 7). At this order, due to the vanishing (2, 2) and (3, 3) elements, a large θ 23 mixing is naturally obtained in this model. At the next leading order, O(ǫ 1 ), the matrix takes the two zero minor structure [18,32]. This form of the neutrino mass matrix confronts a severe constraint on the sum of neutrino masses from cosmological observation [33]. In our model, a quartic term, κ S 2 H † Φ, is allowed by the Z 4 flavor symmetry. Through this coupling, a small VEV for Φ, i.e., Φ ∼ κ ǫ 2 ( S 2 /M 2 φ )v is induced from the singlet VEV. As a result, at O(ǫ 2 ) we have additional contributions to the mass matrix. Then, the total structure of the neutrino mass matrix is Therefore, in the present model, the constraints from neutrino data are relatively relaxed as compared with the minimal gauged U(1) Lµ−Lτ model.

Collider Signature
In the previous section, we have taken into account the direct collider search constraint of charged scalars (m φ 700 GeV). This bound will be improved further at the future LHC running by the same search mode. In our model, the neutral scalars (ρ, η) can be lighter than the charged scalar. Such light scalars can be produced at the LHC, and give interesting distinctive signals. Because of the flavor charge conservation in the Z 4 symmetric limit, they are produced in a pair qq (e + e − ) → ρ η at hadron (lepton) colliders and their primary decay modes are µτ pairs. So far no dedicated search has been performed, and it was shown in the Type-X THDM that 2µ2τ final states can be approximately reconstructed even at the hadron collider [34]. Application of this analysis to our model seems to be easy. Firstly, there is no suppression of the signal events by their branching ratio. Secondly, thanks to the collinear approximation of tau leptons, the LFV invariant mass M µτ is fully reconstructable. Then, M µτ is used for a very good discriminant against background events. The study for the discovery potential of (ρ, η) is beyond the scope of this paper, and we leave it for the future.

Muon Electric Dipole Moment
If the new Yukawa couplings y µτ , y τ µ are complex, electric dipole moment (EDM) of muon is induced by the same diagram for muon anomalous magnetic moment in Fig. 2, which is computed as Similar to the case of muon magnetic moment, Eq. (48) has a potentially large contribution from the chirality flipping effect. The current experimental bound for muon EDM is given by the Muon g − 2 Collaboration (BNL) as [35] |d µ | e < 1.9 × 10 −19 cm.
In addition to the current bound, factor 10 improvement is expected by the future FNAL E989 experiment [36], and the future sensitivity of the J-PARC g − 2/EDM Collaboration is roughly |d µ |/e ∼ 10 −21 cm [37].
In the left panel of Fig. 8, we give a contour plot of the muon EDM predictions in the (m ρ , m η ) plane where we assume Im (y µτ y τ µ ) = 1. The yellow region is already excluded by the current muon EDM limit. We see that the current muon EDM limit does not exclude the model without requiring the tuning in the imaginary part of the Yukawa couplings. The solid purple and dashed green lines are the future sensitivities of the FNAL E989 and J-PARC g − 2/EDM, respectively. Although the constraint of the current bound is not so strong, the future experiments can explore parameter space furthermore.

h → γγ
The additional contribution to h → γγ can appear through the charged scalar loop. The decay amplitude including the SM contribution is computed as [38,39] where τ i ≡ 4m 2 i /m 2 h , N c = 1, 3 is color factor, Q f is the electric charge of the SM fermions, ǫ µ/ν is the photon polarization vector, and the loop functions F 1 (τ ), F 1/2 (τ ) and F 0 (τ ) are given by with The last term in Eq. (50) corresponds to the new contribution which is controlled by the quartic coupling λ 3 in the scalar potential. Then, the partial decay width is calculated as The signal strength for h → γγ defined by the ratio of the observed Higgs boson decay to the SM prediction has been reported as µ = 0.99 +0. 15 −0.14 by the ATLAS Collaboration [40], and 1.18 +0. 17 −0.14 by the CMS Collaboration [41]. The signal strength deviates from unity if non-zero value of the quartic coupling λ 3 exists. The constrained parameter space in the (m φ , λ 3 ) plane is shown in the right panel of Fig. 8, where the red region is excluded by the PDG data µ = 1.16 ± 0.18 at 2σ CL [1], and the orange region is excluded by the LEP limit m φ 93.5 GeV. One can see that the parameter space with |λ 3 | = O(1) is ruled out if 100 GeV m φ 200 GeV. There is no substantial constraint if m φ 200 GeV.

Summary and Conclusions
We have studied models based on leptonic flavor symmetries, which can accommodate the long-standing muon g − 2 anomaly. The minimal model is based on a Z 4 lepton flavor symmetry, and includes an inert doublet scalar charged under the flavor symmetry. Large muon anomalous magnetic moment is realized by the chirality enhancement with the factor m τ /m µ ≈ 17 in this model. We have also analytically formulated the constraints from the electroweak precision tests and lepton universality. Taking into account all these constraints, allowed parameter space is explored numerically. For the electroweak precision tests, it has been found that the constraint can easily be evaded if the quartic couplings λ 4 and λ 5 are relatively small or the relation λ 5 ∼ ±λ 4 is satisfied, which corresponds to one of neutral scalars ρ and η is nearly degenerate with the charged scalar φ + . For lepton universality, we have computed tree and one-loop corrections of heavier charged lepton decays, and one-loop correction for Z boson decay. We have found that the tree level correction becomes dominant when the Yukawa couplings are comparable (y µτ ∼ y τ µ ) while the loop correction becomes important for hierarchical Yukawa couplings. In addition, we have numerically examined the potential stability conditions and triviality bounds assuming the cut-off scale of the model, Λ = 100 TeV. We have successfully found that the parameter region where the discrepancy in the muon g − 2 is explained at 2σ level while satisfying all relevant constraints. As further perspective of the minimal Z 4 model, neutrino mass generation with Type-I seesaw mechanism, discriminative collider signatures, indirect signals from muon EDM and Higgs decay width into γγ have also been discussed. We have also found that some parameter space can be explored by the future EDM experiments if rather large CP phase exists in the Yukawa couplings. The signal strength of the Higgs decay width into γγ is influenced by the new contribution if the charged scalar mass is less than 200 GeV.