Lepton g − 2 and W -boson mass anomalies in the DFSZ axion model

: With regard to the leptonic magnetic dipole moment anomaly as well as the W -boson mass excess, we study the DFSZ axion model. Considering theoretical and experimental constraints, we show that the muon and electron g − 2 anomalies can be explained within the parameter space of the model for extra Higgs bosons with mass spectra around the electroweak scale and for an almost equivalent contribution of one-and two-loop diagrams. A negative electron g − 2 could be achieved by introducing heavy neutrinos. Furthermore, the W boson mass excess can be consistently addressed within the mass range of the matter content testable at collider experiments.

QCD axion models are well-motivated solutions to the strong CP problem of the Standard Model (SM) [1][2][3].The axion appears as a pseudo-Nambu Goldstone boson from the spontaneous symmetry breaking of the U (1) PQ Peccei-Quinn (PQ) symmetry at high temperatures and gains a mass at lower scales around the QCD phase transition [4,5].
The axion being neutral and long-lived, it is a viable candidate for cold dark matter.It could have been produced in the early universe by the misalignment mechanism or by decays of topological defects [6][7][8][9].A review of different axion models and their cosmological implications can be found in Refs.[10,11].Although high-energy experiments have not reported so far any significant deviation from the SM, a set of experiments focusing on precision measurements may indicate anomalies in the SM requiring new physics.One of these high-precision experiments in particle physics measures the muon magnetic moment, a µ ≡ (g µ −2)/2, where g µ is the spin g factor of the muon.Fermi National Accelerator Laboratory (FNAL) has recently reported a large deviation from the SM prediction on a µ [12], consistent with the previous measurement of Brookhaven National Laboratory (BNL) in 2006 [13].The result of combining these measurements is ∆a µ ≡ a EXP µ − a SM µ = (251 ± 59) × 10 −11 , ( which shows a 4.2σ discrepancy with the SM prediction .Moreover, improved measurements of the fine-structure constant α allow us to extract a value for the electron magnetic moment, a e , which shows a discrepancy with the SM prediction.Interestingly, there are two experiments whose results on a e are set in opposite directions.Based on the recoil frequency of 133 Cs atoms at Berkeley [35], one has which features a 2.4σ discrepancy with the SM prediction [36,37], while with 87 Rb atoms, the result of the Laboratoire Kastler Brossel (LKB) [38] is with a 1.6σ deviation.We emphasize that even if both experiments show tension with respect to the SM in opposite directions, new physics may be required.Additionally, another precise measured observable is the W -boson mass.Recently, the CDF-II experiment at Fermilab reported a new measurement of the W -boson mass, using data collected from 2002 to 2011 [39] m CDF−II W = 80.4335 ± 0.0094 GeV, (1.4) which is in tension with the SM prediction m SM W = 80.357 ± 0.006 GeV [40] about 7σ and also differs from the previous average of PDG m PDG W = 80.379 ± 0.0012 GeV [40].Once again, if such a deviation is confirmed, new physics is required. 1o far, a large number of different scenarios beyond the SM have been proposed to address the anomalies mentioned above (see, for example, Refs. ).Here, alternatively, we employ the Dine-Fischler-Srednicki-Zhitnitsky (DFSZ) axion model [67,68], which is an attractive framework that not only addresses the strong CP problem but also provides a viable dark matter candidate.The field content of the model, invariant under the global U (1) PQ , includes a Peccei-Quinn complex scalar field and one additional Higgs doublet.We explore possible explanations for anomalous measurements and their phenomenological implications, restricting the parameter space of the model.Considering vacuum stability and perturbativity conditions as well as collider bounds on the mass of extra Higgs bosons [69,70], we show that for Higgs fields with masses at the electroweak (EW) scale, the CDF-II W -boson mass anomaly and the anomalous muon magnetic moment can be explained.Furthermore, the measurement of the anomalous electron magnetic moment can be fitted in the case of a positive contribution.In the case of a negative g e − 2, taking into account the bounds on charged lepton flavor decay processes such as µ → e γ, extra degrees of freedom (d.o.f.) are required; for example, a heavy neutrino contribution can be added, without spoiling the solutions to the W mass or g µ − 2.
As we discussed, QCD axion is motivated as a solution to the strong CP problem and can be a candidate for dark matter.These are the low-energy phenomenological consequences of QCD axion models.Depending on the axion model, it can have different kinds of phenomenology at high energies.Here, our focus is on the DFSZ axion model to verify phenomenological effects imposed by possible explanations of the electron and muon magnetic moment and W -boson mass anomaly on the model.
In the next section, we introduce the model and discuss the physical d.o.f. and the mass spectrum.In Section 4, we discuss the oblique parameters, particularly the T parameter, and its impact on the W mass. Section 5 is devoted to the one-and two-loop-level contributions to the lepton magnetic moment.We conclude in Section 6.

The Model
Regarding the smallness of the experimentally bounded CP-violating θ term in strong interactions, i.e. the strong CP problem, several QCD axion models, all based on a U (1) PQ symmetry, have been proposed.Different QCD axion models correspond to diverse realizations of the PQ symmetry under which the SM is not invariant.In particular, to absorb the θ term, independent chiral transformations of the u and d quarks are required [3].The model can therefore be enlarged with a second Higgs doublet coupled to only the uor to the d-quark types via Yukawa interactions [4,5].
Moreover, the PQ symmetry could be spontaneously broken at an energy scale much higher than the EW scale, leading to the generation of a pseudo-Goldstone boson, the axion.By imposing astrophysical constraints, the symmetry-breaking scale is bounded to 10 8 GeV ≲ f a ≲ 10 17 GeV [11,71,72], where f a stands for the axion decay constant.This decoupling between the electroweak and PQ symmetry-breaking scales naturally occurs in the presence of a SM-singlet scalar field Φ.This kind of QCD axion model, which has such a field content, is called the DFSZ model. 2he Lagrangian of the model is given by where Φ is the PQ scalar field, H u and H d denote the two Higgses, doublets under SU (2) L .The covariant derivative is defined as where W a µ and B µ are the SU (2) L and U (1) Y gauge fields, respectively, whose coupling constants are denoted by g and g ′ , and τ i correspond to the Pauli matrices.The scalar potential is given by Finally, the Yukawa interactions are where The model is invariant under the SM gauge symmetry group SU (3) C ⊗SU (2) L ⊗U (1) Y and under the U (1) PQ , and therefore the PQ charges of the fields in this model can be obtained accordingly [75].In fact, from the κ term of Eq. (2.3) and the orthogonality between PQ and corresponding hypercharge currents [11], we obtain PQ charges of scalar fields as X Hu = cos 2 β, X H d = − sin 2 β, and X Φ = 1.Moreover, supposing that left-handed fermions have no PQ charges, X u R = X Hu , and The Higgs doublets can be expanded as [76] and their vacuum expectation values (VEV) are given by ⟨α [11].
Using the minimization of the potential, three parameters µ d , µ u , and V ϕ , can be fixed (see Appendix A); therefore, in the scalar sector, the free parameters are λ ϕ , tan β = t β ≡ v u /v d , κ, κ d , κ u , λ u , λ, and v ϕ .Finally, we note that λ d is fixed by imposing that one of the Higgs bosons has a mass of 125 GeV, as the SM Higgs.

Physical degrees of freedom in the scalar sector
In this section, we determine the mass spectrum of the scalar sector.From two complex doublets and a complex singlet, there are 10 d.o.f.Three combined Goldstone components are related to gauge bosons and are eliminated by gauge transformations, corresponding to longitudinal parts of the three massive gauge bosons [76,77].After rotating the interaction eigenstates (details of the transformation can be found in Appendix B), the mass eigenstates can be identified.The spectrum counts a pair of charged scalars H ± , a pseudoscalar A, three CP-even scalars h i , and the axion field a.
First, the mass of the pseudoscalar Higgs A is given by From Eq. (2.3), the mass of the charged Higgs fields H ± is Furthermore, there are three neutral scalar states, H, S, and ϱ, which are not mass eigenstates and should be expressed in the mass basis where the mass matrix can be diagonalized by a rotation matrix (Appendix C) where h i are mass eigenstates and the rotation matrix, R, and its components are expressed in Appendix C. The mass spectrum of these scalar fields, up to (2.9) The lighter state h 1 is identified with the SM-like Higgs boson, while h 2 and h 3 are extra heavy neutral CP-even Higgs bosons.We do not consider cases where the mass spectrum is decoupled from the EW energy scale, but cases that can address the aforementioned anomalous problems and lead to rather light masses for the Higgs fields, phenomenologically interesting for collider physics.In this sense, after fixing t β , the mass spectrum of the Higgs fields depends mainly on κ, λ, and λ ϕ .We consider m h 1 as the SM Higgs and also fix λ ϕ = 0.01 that affects m h 3 so that the very massive scalar field can be determined close to v ϕ due to v ϕ ≫ v. 3 The case with λ ϕ → 0 can lead to masses lighter than the SM Higgs mass.SM-like interactions of the Higgs, h 1 , with SM particles can also consistently be obtained, due to the relevant coefficient, a combination of R H1 ∼ 1 and R S1 proportional to v 2 /v 2 ϕ , multiplied to SM Higgs couplings.For more details see Ref. [76].

Collider constraints
Constraints on the mass spectrum of extra Higgs fields in the framework of the Minimal Supersymmetric Standard Model (MSSM) [78] and two-Higgs-doublet models (2HDMs) [77] are extensively studied [40].Regarding the mass of the charged Higgs boson, very light masses, m H ± ≲ 80 GeV have been excluded by LEP [79] and Tevatron [80].Moreover, for m H ± < m t (the top quark mass), the experimental lower bound from direct searches is m H ± > 155 GeV [40].For mass ranges m H ± ∼ 155 − 170 GeV, firm experimental analysis does not exist [81] and a reliable perturbative calculation of the charged Higgs boson production cross section must be performed.For m H ± > m t , ATLAS and CMS have also been excluded regions of the parameter space of MSSM, m H ± < 180(1100) GeV, for t β = 10(60), with a certain luminosity [40].This search is sensitive to the modeling of the top pair production background with extra partons and especially b-quarks.
Also, there have been numerous experimental searches at the LHC for a pseudo-scalar boson, excluding light pseudoscalar Higgs bosons m A ≲ m h 1 /2 [70].In addition, in the context of MSSM, for heavy neutral/pseudo-scalar Higgs bosons in the τ τ final state, the exclusion region is m A ≲ 390(1600) GeV for t β = 10(60) [40].However, in the diphoton channel a slight excess has been observed by CMS at a mass of 95.3 GeV, consistent with the observed excess at LEP, with a local significance of 2.8σ and could not be ruled out by ATLAS at 95% CL [40].
Experimental searches including extra Higgs production crucially depend on accurate theoretical predictions for inclusive cross sections.(For instance see Ref. [82] for precise predictions of pseudo-scalar inclusive cross section specially for the moderate mass range.)According to the matter content of our model and its parameter space, given additional decay channels, for instance H ± → aW ± , A → ah 1 , h 2 → aa, and h 2 → aZ, the bounds may somehow change. 4Here, aside from the exclusion region independent from t β , i.e., m H ± ≲ 155 GeV and m A ≲ m h 1 /2, we consider other mass values, provided that the model parameter t β is in the viable, allowed range.In the following calculations, we show that the anomalies may simultaneously be explained for t β > 100.

W -boson mass
In this section, we study the effects of new physics on the recently reported W -boson mass excess by analyzing the oblique parameters of this model [83,84].Some studies have already considered the impact of having a second Higgs doublet in the W -mass anomaly; see, e.g.Refs.[42, 43, 46, 48, 50, 51, 53-58, 60, 85].The leading correction to the W -boson mass can be approximated as where ∆T is one of the Peskin-Takeuchi parameters, obtained by the following procedure.The one-loop corrected masses of the gauge bosons are given by where Π W W (0) and Π ZZ (0) are the two-point functions of the gauge fields at zero-momentum transfer [76].The deviation from the custodial symmetry can be parameterized as where m 2 V (tree) is the mass at tree level of the gauge boson V , and . (4.5) where R Si is given in Appendix C, f (x, y) ≡ x y ln(x/y)/(x − y) and f (x, x) = x.As a result, a positive value of ∆T can explain the W -boson mass excess.We calculate the W -boson mass excess as a function of m H ± , fixing all other parameters and masses.For a given value of m A and m h 2 , the excess, including the CDF-II measurement, can be justified within m where ∆T = 0. Note that within the considered parameter space the h 3 contribution is negligible.In this case, the low-energy effective theory would be a 2HDM model that has the advantage of avoiding the strong CP problem thanks to the QCD axion.To relate the parameters of the DFSZ potential to those of a 2HDM, see Ref [76].Also, in the whole studied problems, we also subtract the SM-like Higgs contribution in the calculation.As we shall discuss in the following sections, the region of the parameter space that fits the lepton g l − 2 anomalies corresponds to small values of κ, κ d and κ u , therefore, we fix ϕ and κ = κ v 2 /v ϕ with κ < 1. 5 Furthermore, considering vacuum stability conditions (Appendix A), we impose perturbativity criteria as λ, λ d,u ≲ 4π.Therefore, in Fig. 1, we set other parameters as follows.Taking into account m h 1 = 125 GeV, we can find λ d,u , fixing t β and λ d − λ u .In this sense, since we obtain λ d − λ u ≈ λ d , we can fix λ d − λ u < 4π.In our calculations, we take λ d − λ u = 5.Within the considered parameter values, we can approximately express the mass terms of scalar fields, Eq. (2.9), with their first term for each relation.
Additionally, based on the perturbativity constraints from the Yukawa sector Eq. (2.4), we consider a conservative upper bound t β ≤ 300.According to the mass relation of the scalar fields, m h 3 , m h 2 , and m A are also determined by fixing λ ϕ , and then κ.Thus, we can plot m W as a function of m H ± .We consider m H ± as a variable, although its different values can be obtained by fixing λ and κ.As shown in Fig. 1, we can obtain the W -boson mass excess within the parameter space of the theory, showing by three different benchmarks for

Lepton g l − 2 anomalies
In this section, we first compute and discuss the analytical expressions for the lepton g l − 2 anomalies, and then we proceed to calculate the magnetic moment of the muon and electron.On the basis of Lorentz covariance and gauge symmetry, the vertex of the QED interaction of a lepton can be expressed as where F 1 (q 2 ) is related to the electric charge, and the lepton magnetic moment is defined as a l = F 2 (0).Here, we calculate different ∆a l contributions from the model.Regarding the Yukawa couplings, using a unitary transformation on the charged leptons, these Yukawa coupling matrices can be diagonalized, ỹl , so that from Eqs. and (B.4), the Yukawa interactions of the leptonic sector, contributing to a l , at one-loop level, can be expressed as Therefore, from Eqs. (B.7) and (C.10) for mass eigenstates we obtain the one-loop amplitude, Fig. 2, as follows (checked by using Package-X [87]) where x h i = m l /m h i , and x A = m l /m A . 6 The loop functions obtained, F (x) and F H ± (x, z) are given in Appendix D. According to Eq. (2.4), the so-called Barr-Zee two-loop amplitudes can be obtained [89][90][91], where the dominant contribution, Fig. 3, is given by where N f C is the number of colors and Q f is the electric charge.Note that the two-loop contribution is at the α and ỹ2 l level.We consider the dominant contribution from the t and b quarks, and therefore The functions F(ω) and F(ω) are defined in Appendix D.

Muon g µ − 2 anomaly
Having obtained analytical formulas for g l − 2, we investigate possible solutions to the anomalous magnetic moment of the muon and for the recently measured ∆a µ in Eq. (1.1).We show that the observed ∆a µ can be obtained for a range of t β and small values of κ, 6 The contribution of the axion to Eq. (5.3) is proportional to F (xa) and is therefore negligible.
For the contribution of EW gauge bosons, see, e.g.Ref. [88].κ d , and κ u , which is also consistent with the expansion of the scalar field mass spectra in terms of v/v ϕ .In Fig. 4, we fix GeV and obtain ∆a µ for a range of κ and some representative values of t β that simultaneously determine the value of m h 3 , m h 2 and m A , showing in terms of m A .
As it can be seen in Fig. 4, we can find the observed anomaly for a range of t β values.Taking into account the bounds from colliders on m A ≳ m h 1 /2 [70], in addition to the small t β , we can find the observed ∆a µ for t β ≳ 1007 .It should be noted that for 1 ≲ t β ≲ 5 the two-loop Barr-Zee contribution is dominated, whereas for other values of t β the oneand two-loop contributions, which are at the same order in ỹl , are almost equivalent.Additionally, in the small values of t β , the top quark in the two-loop diagram contributes predominantly, while in the large values of t β , the dominant contribution comes from the bottom quark.Furthermore, with λ d − λ u ≲ 0.2, the result of interest is obtained only for large t β .

Electron g e − 2 anomaly
As already mentioned, due to the discrepancy in the measurement of α of Berkeley and LKB experiments, two different values have been reported for ∆a e , with opposite signs.Although results of the two experiments are rather inconsistent, we also study the averaged value case, combing the two results [94] ∆a Avg e = (−2 ± 2.2) × 10 −13 . (5.5) We first try to explore possible explanations for ∆a LBK e and ∆a Avg e .We then study the physics which leads to the explanation of ∆a B e , considering the constraints imposed by charged lepton flavor violation decays.
GeV, and different values of t β , we show that the observed ∆a µ can be obtained, plotting as a function of m A .The gray area shows the excluded region based on the bound m A ≳ m h1 /2 [70].

∆a LBK e
According to the LKB experiment and its result, Eq. (1.3), and because of the muon and electron mass difference, m e /m µ ∼ 5 × 10 −3 , it is expected to obtain positive ∆a e with the similar benchmark chosen for the case ∆a µ .Therefore, as shown in Fig. 5, we can find that ∆a LBK e can be obtained as a function of m A and t β .It is shown that ∆a LBK e can be explained for 1 ≲ t β ≲ 3 where the two-loop Barr-Zee contribution is dominated, while for large t β , this is fulfilled for t β ≳ 500.In these regions of the parameter space, based on our previous discussion, both the W -boson mass and ∆a µ anomalies can be explained.However, as can be seen from Figs. 5, 4, the range of m A values explaining ∆a LBK e and ∆a µ is different and thus both cannot be simultaneously explained.
As shown in Fig. 6, we also find that ∆a Avg e can be obtained within the parameter space of the model.Interestingly, given the averaged value measured for the electron magnetic moment, in addition to the electron magnetic moment, the W -boson mass, and muon magnetic moment anomalies can be simultaneously explained for t β ≳ 120.For instance with the following benchmark, for Figure 7 with red and green color, respectively.

∆a B e
As we discussed in the previous section, the parameter space of the model in which the W -boson mass and ∆a µ anomalies can be explained gives rise to a positive electron g − 2. Therefore, to explain the case of negative electron magnetic moment, Eq. (1.2), the mentioned model should be extended.Here we consider the model supplemented with heavy RH neutrinos, which are PQ charged and SM singlet, and study their contribution to the lepton g l − 2.
The Lagrangian for right-handed neutrinos is given by where N Rα are RH neutrinos.Then, the Yukawa interaction of the charged Higgs and RH neutrinos can be expressed as8 √ 2 t β y l Nα H + l P R N α + H.c. .amplitude should be calculated and (checked by using Package-X [87]) the result would be as follows where and n R is the number of RH neutrinos.In our calculations, we assume n R = 3 and the same properties for all generations.Since this additional contribution, Eq. (5.8), is negative, its effect on the muon magnetic moment should also be noted.For very heavy RH neutrinos, the contribution of this term is negligible.However, as shown in Fig. 9, for y µ N ≳ 0.01 and sufficiently light RH neutrinos, ∆a µ can be negative, and hence we consider it y µ N < 0.01 in the calculations.Before proceeding to calculate the electron magnetic moment for this case, we first discuss another constraint.Charged lepton flavor violations are highly constrained and the bound on the branching ratio of the decay µ → eγ is more severe [95] Br(µ → eγ) < 4.2 × 10 −13 .
(5.9) Thus, we can constrain the model as follows where and the function obtained F H ± is given in Appendix D. 9 According to Eq. (5.9), we should have |a µe | < 2.5 × 10 −14 .To pass the flavor violation decay bound, based on the previous estimation on y µ N < 0.01, with calculating a µe , the hierarchy y µ N ≪ y e N can be imposed, taking y e N < 4π.For example, in Fig. 10, we show that |a µe | < 2.5 × 10 −14 can be obtained for y e N ≤ 0.2 and y µ N = 10 −7 .Considering such bounds on the parameters, we obtain the electron magnetic moment and if ∆a B e is the case, we find that taking y e N < 4π, the measured quantity can be obtained for m N R ≲ 200 TeV.In Fig. 11, we show the result for y e N ≤ 0.2.

Conclusion
Recent measurements of the muon magnetic moment, the W -boson mass, and their deviation from the SM prediction can imply new physics beyond the SM.In this paper, in light of these observations we have studied a QCD axion model, enjoying the PQ symmetry, as an attractive framework that can address other SM shortcomings including the strong CP and the DM problem.
We considered the DFSZ axion model and calculated loop-level diagrams that contribute to the muon and the electron magnetic moment, and W -boson mass.Considering the lower bound on the mass of charged Higgs, we have shown that the CDF-II W boson mass excess can be addressed within the parameter space of the model, except for the case m H ± = m A , m h 2 .We have also shown that the positive muon magnetic moment ∆a µ can  be obtained for a range of t β values with a mass range of extra Higgs bosons at the EW scale, respecting the direct collider bound on the pseudoscalar field and also theoretical constraints.
Regarding the electron g − 2, there are Berkeley and LKB experiments whose measurements are set in opposite directions.We showed that ∆a LBK e can be explained within a different region of the parameter space allowed for ∆a µ .However, we found that considering the averaged value of the two experiments, the three anomalies can be simultaneously explained.The electron magnetic moment measured by another experiment at Berkeley, ∆a B e , is of a different sign and the physics that explains ∆a µ and the W -boson mass excess leads to a positive electron g − 2. To address ∆a µ , ∆a B e , and the W -boson mass anomalies simultaneously, we extended the model with heavy RH neutrinos.Taking into account the flavor violation decay bounds, we found that heavy neutrino-charged Higgs interactions contributing to one-loop diagrams that can be responsible for the negative ∆a e with heavy neutrinos masses up to 200 TeV.under the Supporting TAlent in ReSearch@University of Padova (STARS@UNIPD), Istituto Nazionale di Fisica Nucleare (INFN) through the Theoretical Astroparticle Physics (TAsP) project.FH also thanks the organizers of the "Probing New Physics with Gravitational Waves" Workshop 2022 at the Mainz Institute for Theoretical Physics (MITP) of the Cluster of Excellence PRISMA (Project ID 39083149) in Mainz, for hospitality and partial financial support during the completion of this work.

A Minimization and vacuum stability conditions
We can use the minimization of the potential, Eq. (2.3), and remove three parameters µ d , µ u and V ϕ , expressing in terms of free parameters as follows Assuming that κ d , κ u and κ are sufficiently small and do not affect the vacuum stability conditions obtained from other quartic terms of the potential, we have

B Rotations for the Higgs fields
With the following rotated Higgs doublets, one can find the direction by which unphysical components are removed where c β ≡ cos β and s β ≡ sin β.The upper components are combinations absorbed by W ± and the real part is the component with non-zero vev, thus H ′ d has the structure of SM doublet so that and ⟨H + v⟩ = v.Also, the other rotated doublet would be 3) whose charged components denote the charged Higgs H ± and the real neutral component is a scalar Higgs S such that ⟨S⟩ = ⟨−s β Re [α 0 ] + c β Re [β 0 ]⟩ = 0. Therefore, the gaugeredefined doublets can be expressed as Although Ã remains after gauge redefinition, the field is not yet physical.Moreover, from ϕ = v ϕ + ϱ + i G, G is another non-physical d.o.f.whose combination with Ã identifies the physical one in the pseudoscalar sector [76].In addition to the massless eigenstate, the axion resulting from the PQ current is There is also a massive pseudoscalar eigenstate, evident from the cubic term in the potential, as The rotation matrix can be written as [76]

Figure 2 .
Figure 2. One-loop Feynman diagrams contributing to g l − 2 from Higgs fields.

10 9LKBFigure 5 .
Figure 5. Similar to values chosen for parameters in the Fig. 4, for different values of t β , we show that the observed ∆a LBK e

(5. 7 )Figure 6 .
Figure 6.Similar to values chosen for parameters in the Fig. 4, for different values of t β , we show that ∆a Avg e can be obtained, plotting as a function of m A .The whole range of ∆a Avg e axis in this plot is allowed based on the average experimental bound.Because of the obtained positive results, negative values of Eq. (5.5) in the ∆a Avg e axis are not shown.

Figure 7 .Figure 8 .
Figure 7. Benchmark points, associated with t β = 200 (red) and t β = 300 (green), are represented, satisfying different experimental bounds for M W , ∆a µ and ∆a Avg e discussed in previous sections.Also, different experimental constraints are shown in three panels.

Figure 9 .
Figure 9.For κ d = κ u = v 2 /v 2 ϕ , m h1 = 125 GeV, λ ϕ = 0.01, λ d − λ u = 5, v ϕ = 10 9 GeV, t β = 200, m A = 130 GeV and y µ N = 0.01, we find ∆a µ as a function of m N R and m H ± .The colored regions show different values of the muon magnetic moment which could be even negative for small values of m N R with y µ N = 0.01.For y µ N < 0.01 the negative contribution, Eq. (5.8), to ∆a µ would be negligible.The region on the left of dashed grey line is excluded by the PDG (g − 2) µ bound.