Anomalous magnetic and weak magnetic dipole moments of the $\tau$ lepton in the simplest little Higgs model

We obtain analytical expressions, both in terms of parametric integrals and Passarino-Veltman scalar functions, for the one-loop contributions to the anomalous weak magnetic dipole moment (AWMDM) of a charged lepton in the framework of the simplest little Higgs model (SLHM). Our results are general and can be useful to compute the weak properties of a charged lepton in other extensions of the standard model (SM). As a by-product we obtain generic contributions to the anomalous magnetic dipole moment (AMDM), which agree with previous results. We then study numerically the potential contributions from this model to the $\tau$ lepton AMDM and AWMDM for values of the parameter space consistent with current experimental data. It is found that they depend mainly on the energy scale $f$ at which the global symmetry is broken and the $t_\beta$ parameter, whereas there is little sensitivity to a mild change in the values of other parameters of the model. While the $\tau$ AMDM is of the order of $10^{-9}$, the real (imaginary) part of its AWMDM is of the order of $10^{-9}$ ($10^{-10}$). These values seem to be out of the reach of the expected experimental sensitivity of future experiments.


Introduction
The most general dimension-five effective interaction of a neutral V gauge boson (V = γ, Z) to a charged lepton that respects Lorentz invariance can be written in terms of six independent form factors: a marcofis@yahoo.com.mx b tome_gerar@hotmail.com c gtv@fcfm.buap.mx where q is the incoming transfer four-momentum of the gauge boson. The electromagnetic (weak) properties of the lepton are determined by the photon (Z gauge boson) vertex function. The CP-violating terms define the static electric dipole moment (EDM) and the static weak electric dipole moment (WEDM): Since the standard model (SM) predictions for these CPviolating dipole moments are highly suppressed, they may serve to search for new sources of CP violation. In this work we are interested instead in the static anomalous magnetic dipole moment (AMDM) and the static anomalous weak magnetic dipole moment (AWMDM), which are defined in terms of the CP-even form factor as follows: The measurement of the AMDM of a lepton has long been considered a probe for the SM, which considers leptons as point-like objects. The AMDM of the electron a e , which receives its main contributions from quantum electrodynamics (QED), has been calculated up to order of α 5 [1] and the agreement between the theoretical and the experimental values has reached the level of ten significant digits [2], which represents one of the greatest milestones of QED. As far as the muon is concerned, the E821 experiment at Brookhaven National Lab (BNL) measured its AMDM arXiv:1612.09537v3 [hep-ph] 17 Mar 2017 a µ with an unprecedent precision of 0.54 ppm. The current average experimental measurement is [3] a Exp. µ = 11659209.1(5.4)(3.3) × 10 −10 , and further improvement is expected ats future experiments by the (g − 2) µ [4] and J-PARC (g − 2)/EDM [5] Collaborations, which aim to reach a precision of ±0.2 ppm. As for the theoretical prediction of the SM [6] a SM µ = 116591803(1)(42)(26) × 10 −11 , there is still a large uncertainty in the estimate of the hadronic contribution, whereas the QED and electroweak contributions have been determined with a great precision [7]. Thus a more accurate evaluation of the leading order hadronic contribution together with the future experimental measurement are needed to settle down the discrepancy between the SM prediction and the experimental value of a µ , which currently stands at the level of 3.6 standard deviations [6]: µ − a SM µ = 288(63)(49) × 10 −11 .
Since the muon AMDM has become a powerful tool to test the validity of the SM and searching for new physics (NP) effects, a plethora of calculations within the framework of several SM extensions has been reported in the literature in order to explain the ∆ a µ discrepancy [7,8].
As far as the τ lepton is concerned, the SM prediction is a SM τ = 117721(5) × 10 −8 [9,10]. The error of the order of 10 −8 hints that SM extensions predicting values for a τ above this level could be worth studying. Since the SM prediction for a τ is far from the experimental sensitivity, which is one order of magnitude below the leading QED contribution, a more precise determination of the experimental value is necessary. Due to its short lifetime (290.3 ± 0.5 × 10 −15 s) [11], the τ lepton does not allow for a high precision measurement of its AMDM via a spin precession method. The most stringent current bound on a τ is [12] −0.052 < a τ < 0.013, (9) which was obtained using LEP1, SLD and LEP2 data for τ lepton production. It also has been pointed out recently that the τ electromagnetic moments can be probed in γγ and γe collisions at CLIC, which can lead to improved bounds [13,14]. In this regard, it has been pointed out that super B factories could allow for a precise determination of a τ up to the 10 −6 level using unpolarized or polarized electron beams [15][16][17]. Furthermore, due to its large mass, it is expected that the τ AMDM can be very sensitive to NP effects [18] since its electroweak contribution would be ten times larger than the uncertainty of the hadronic contribution [10]. Therefore, it is worth estimating the τ AMDM in any SM extension as future measurements may allow us to search for NP in a rather clean environment.
Unlike the electromagnetic dipole moments of leptons, little attention has been paid to the study of their weak properties. In the experimental arena, the current best limits on the τ AWMDM and WEDM, with 95% C.L., are [19] Re a W τ < 1.14 × 10 −3 , Im a W τ < 2.65 × 10 −3 , Re d W τ < 0.5 × 10 −17 ecm, Im a W τ < 1.1 × 10 −17 ecm, (13) which were extracted from the data collected at the LEP from 1990 to 1995, corresponding to an integrated luminosity of 155 pb −1 . Somewhat weaker bounds were obtained in [20] via a study of the pp → τ + τ − and pp → Zh → τ + τ − h cross sections at the LHC. The current experimental bounds on a W τ are well above the SM theoretical prediction, which was calculated in Ref. [21]: It is thus interesting to analyze whether NP contributions can give a significant enhancement and be at the reach of future experimental detection.
In this work we evaluate the AMDM and AWMDM of a lepton, with special focus on those of the τ lepton, predicted by the simplest little Higgs model (SLHM) [22], which is an appealing SM extension. This model is aimed to deal with the hierarchy problem by conjecturing that the Higgs boson is a pseudo-Goldstone boson arising from a global symmetry broken spontaneously. At the same scale, the local symmetry is also broken by a collective symmetry breaking mechanism. The top quark and the electroweak gauge bosons have heavy partners that give rise to new contributions that exactly cancel the quadratic divergences to the Higgs boson mass at the one-loop level, thereby rendering a mass of about one hundred GeV without the need of fine tuning [23,24].
The rest of this presentation is organized as follows. A brief review on the SLHM is presented in the Section 2, whereas Sec. 3 is devoted to the analytical results for the AWMDM. As a byproduct we will obtain the corresponding expressions for the AMDM. A brief discussion on the current constraints on the parameters of the model, and the numerical analysis of the τ electromagnetic and weak dipole moments is presented in Section 4. Section 5 is devoted to the conclusions, whereas the SLHM Feynman rules as well as explicit expressions for the loop integrals are shown in the Appendices.

The simplest Little Higgs Model
We now present an overview of the SLHM focusing only on the details relevant for our calculation. For a detailed account of this model and the study of its phenomenology we refer the reader to Refs. [22,[25][26][27][28], which we will follow closely in our discussion below. The SLHM is the most economic version of simple-group little Higgs models, which have the feature that the SM gauge group is embedded into a larger simple gauge group instead of a product gauge group. The SLHM has a [SU(3) ×U(1)] 2 global symmetry and a SU(3) L ×U(1) X gauge symmetry, which requires the introduction of nine gauge bosons. At the TeV scale, the global symmetry is broken down spontaneously to [SU(2) ×U(1)] 2 via the vacuum expectation values (VEVs) f 1 and f 2 of two sigma fields Φ 1 and Φ 2 , giving rise to ten Goldstone bosons. At the same scale the gauge group breaks to the SM gauge group SU(2) L ×U(1) Y and five Goldstone bosons are eaten by the heavy fields: a charged gauge boson X ± , a no-self conjugate neutral boson Y 0 = Y † 0 , and an extra neutral gauge boson Z , which thus acquire masses of the order of the scale f 1 ∼ f 2 . The remaining Goldstone bosons are accommodated in a complex doublet (the SM one) and a real singlet of SU (2). The Goldstone bosons can be parametrized by the triplets where the pion matrix is here f = f 2 1 + f 2 2 , h is the SU(2) complex doublet of the SM, and η is a real scalar field. The normalization is chosen to produce canonical kinetic terms. The dynamics of the Goldstone bosons is described by a non-linear sigma model with the SU(3) L ×U(1) X covariant derivative where T a (a = 1 . . . 8) are the SU(3) L generators in the fundamental representation, A a µ are the SU(3) L gauge fields, B X µ is the U(1) X gauge field, and Q X = −1/3 for Φ 1 and Φ 2 .
The new gauge bosons accommodate in a complex SU(2) × U(1) doublet (X + , Y 0 ) with hypercharge 1 2 and a neutral singlet Z 0 . The matching of the gauge coupling constants yields with t W = s W /c W the tangent of the Weinberg angle θ W . As mentioned above, after the first stage of symmetry breaking, there emerge the pair of charged gauge bosons X ± , the no-self conjugate neutral gauge boson Y 0 , and the neutral gauge boson Z 0 (following Ref. [25] we will denote the gauge eigenstates with the subindex 0) with masses of the order of f . Four gauge fields remain massless at this stage. While A 1 , A 2 , A 3 identify with the SU(2) L gauge bosons W a , the charged gauge bosons W ± and the hypercharge gauge boson are After the electroweak symmetry breaking (EWSB), the weak gauge bosons W ± and Z 0 acquire mass and the heavy gauge boson masses get corrected. Up to order (v/ f ) 2 W , X and Y 0 coincide with the mass eigenstates and their masses are [26] with t β = tan β = f 1 / f 2 . If higher order terms are considered W and X need to be rotated to obtain the physical states [26]. On the other hand, the photon and the light neutral Z 0 gauge boson are given by Finally, Z 0 and Z 0 need to be rotated to obtain the mass eigenstates Z and Z , which are given by The respective masses, up to order (v/ f ) 2 , are [26] 4 The kinetic Lagrangian of the gauge bosons gives rise to the trilinear gauge boson couplings necessary for our calculation. It can be written as with the Abelian and non-Abelian gauge strength tensors and with f abc the structure constants of the SU(3) group. From the relations between gauge eigenstates and mass eigenstates (20)- (24) and (28)- (31) we can obtain after some lengthy algebra the Feynman rules listed in Appendix A for the VV ± j V ∓ j vertices, namely AW ± W ∓ , AX ± X ∓ , ZW ± W ∓ , and ZX ± X ∓ .
In the lepton sector of the SLHM, for each generation there is a left-handed triplet L T m = (ν L m , L m , iN Lm ), which is completed with a new neutral lepton N Lm , and two righthanded singlets Rm and N Rm . The Yukawa Lagrangian can be written, in the basis where flavor and mass N m eigenstates coincide, as (37) where Λ = 4π f is the cut-off of the effective theory. Here m and n are generation indices, whereas i, j, and k are SU(3) indices. After EWSB this Lagrangian yields the lepton masses and the heavy neutrino masses up to order (v/ f ) 2 [26] L mass = − f s β λ m where represents the mixing between a heavy neutrino and a SM neutrino of the same generation. Notice also that the rotation that diagonalizes λ N does not necessarily diagonalizes λ so there is mixing between the charged leptons and the heavy neutrinos mediated by the charged gauge bosons. The charged lepton mass eigenstates Lm are thus related to the flavor eigenstates Lm0 by the rotation where V mi is a CKM-like mixing matrix. Also, in each generation, the SM and heavy neutrino mass eigenstates are obtained through where again the 0 subindex stands for flavor eigenstates. The lepton masses up to order (v/ f ) 2 are [26] where y is the eigenvalue of the λ matrix and The vertices of a gauge boson to a lepton pair are obtained from the lepton kinetic Lagrangian, which can be written as where the covariant derivative was given in Eq. (18), with Q X = −1/3, 0 and 1 for L m , N m and m . We need to introduce the mass eigenstates to obtain the interactions of the physical gauge bosons Z, W , X, and Z to a lepton pair, which are necessary for our calculation of the AMDM and AWMDM of a lepton. They are given by and Notice that there is lepton flavor violation mediated by the charged gauge bosons. Finally, the interactions of the photon with a charged lepton pair are dictated by QED The scalar Higgs bosons H and η also contribute to a lepton AMDM and AWMDM. From the Lagrangian (17) we can obtain the respective interactions with the Z and Z gauge bosons. After some algebra one can extract the vertices ZZ H and ZZH, which are given as follows to the leading order in (v/ f ) as together with the ZZH interaction From the Yukawa Lagrangian (37) we can also obtain the interactions of the scalar Higgs bosons H and η to leptons, which we need for our calculation. They are diagonal and are given to leading order in (v/ f ) 2 by [25] and The ηZZ vertex vanish since η is CP-odd scalar, but the HηZ coupling does arise, though its contribution to the AMDM and AWMDM of a lepton vanishes. A similar result was found in Ref. [29], where the contributions of two-Higgs doublet models (THDMs) to the AWMDM of a fermion were calculated.
Other details of this model are irrelevant for our calculation, so we refrain from presenting a discussion of the quark sector and the Coleman-Weinberg scalar potential.

Anomalous magnetic and weak magnetic dipole moments in the SLHM
We now turn to present our results. All the Feynman rules necessary for our calculation follow straightforwardly from the above interaction Lagrangians and are presented in Appendix A. Since we are interested in the¯ V µ vertex with all the particles on their mass shell, the loop amplitudes will be gauge independent. We used the unitary gauge as it is best suited for our calculation method. In order to solve the loop integrals, we used both Feynman parametrization and the Passarino-Veltman reduction scheme. We will first present the results for the AWMDM, from which the results for the AMDM will follow easily. As far as the CP-violating properties are concerned, we will assume that there is no new sources of CP-violation in this model's version, so both the EDM and the WEDM will vanish. In the most general scenario, CP-violation could arise from new additional phases in the extended Yukawa sector of the SLHM.

Anomalous weak magnetic dipole moment
Feynman diagrams that contribute to the AWMDM of charged lepton i at the one-loop level in the gauge sector of the SLHM. Here V can be either the W gauge boson or the new charged X gauge boson, ν i and N j stand for a SM neutrino and a new heavy one predicted by the SLHM, respectively. Notice that diagram (4) involves the non-diagonal vertex Zν i N j .
In the SLHM, in addition to the pure SM contributions, the AWMDM receives new physics contributions arising from the loops carrying only new particles, but also from loops involving only SM particles. The latter are due to corrections to the SM vertices and appear as a series of powers of v/ f , so it is enough to consider the leading order terms. The AWMDM of lepton i can thus by written as where a W −SM i stands for the SM contributions and a W −NP i for the new physics ones, which can be written as with a W −Gauge i a W −Scalar i the contributions arising from the gauge (scalar) sector of the SLHM.
In the gauge sector, the NP contributions to the AWMDM arise from the Feynman diagrams of Fig. 1, where V stands for the charged gauge bosons W and X, ν i is a SM neutrino, and N k is a heavy neutrino. According to our discussion above, for the loops involving only W gauge bosons and SM neutrinos we consider the leading order v/ f contributions arising from the corrections to the W¯ ν vertex. The contributions to the AWMDM of each Feynman diagram will be written as follows where the three-letter superscript stands for the particles circulating in each loop diagram (A 2 and A 3 are the particles attached to the external Z gauge boson, whereas A 1 is the particle attached to the external leptons). Here f are coefficients involving all the couplings appearing in each amplitude, whereas I stand for the loop integrals, which depend on the masses of the virtual particles. We present in Appendix B both the f coefficients and the loop integrals in terms of parametric integrals and Passarino-Veltman scalar functions. We have verified that the contribution of each diagram to the AWMDM is free of ultraviolet divergences. The full contribution of the gauge sector can be written as where the index j runs over the three lepton families.
In the scalar sector of the SLHM there are contributions to the AWMDM of a lepton arising from both the SM scalar boson H and the new pseudoscalar boson η via the Feynman diagrams of Fig. 2. The contributions of the SM Higgs boson arise from corrections of the order of (v/ f ) 2 to the SM vertices H and HZZ. The respective contributions to the AWMDM can also be written as in Eq. (56), where the f functions are presented in Appendix B. Again we have verified that the contribution of each diagram to the AWMDM is ultraviolet finite.
The full scalar contribution is thus

Anomalous magnetic dipole moment
In the gauge sector of the SLHM, the AMDM of a lepton arises from the Feynman diagrams (1) and (2)  similar to diagram (1) of Fig. 2. The corresponding contributions to the AMDM can be obtained straightforwardly from those to the AWMDM by considering the limit m Z → 0 and substituting the Z coupling constants by those of the photon. We can write the contributions to a arising from each diagram as where again the three-letter superscript corresponds to the three particles circulating in the loop. Explicit expressions for the f functions are presented in Appendix B in terms of parametric integral and Passarino-Veltman scalar functions.
The overall NP contribution of the SLHM to a i is thus where the index j runs over the three lepton families. It is worth mentioning that our results for the τ AMDM in terms of parametric integrals, obtained by a limiting procedure from our results for the AWMDM, agree with previous calculations presented in the literature [7,30]. This serves as a cross-check for our calculation.

Numerical analysis
We now present our numerical results for the AMDM and AWMDM of the τ lepton in the context of the SLHM. We will briefly review the existing bounds on the free parameters of the model and afterwards analyze the potential contributions to the AMDM and the AWMDM of the τ lepton for parameter values consistent with these bounds.

Bounds on the parameter space of the SLHM
The SLHM parameters involved in our calculation are f , t β , m N k , m η , δ ν , and the matrix elements V mi l . We will discuss  [31] Electroweak precision data 5.6 (95% C.L.) [32] the current bounds on these parameters obtained from the study of experimental data of several observables as reported in the literature.
Symmetry breaking scale f : bounds on this parameter arise from several observables. We list the most relevant in Table 1. We can observe that the most stringent bound f ≥ 5.6 TeV arises from electroweak precision data (EWPD), whereas the weakest limit f ≥ 1.7 TeV arises from parity violation in Cesium.
f 1 to f 2 ratio t β : a fit on 21 electroweak precision observables from LEP, SLC, Tevatron, and the Higgs boson data reported by the LHC collaborations ATLAS and CMS, allowed the authors of Ref. [33,34] to find out the allowed region in the t β -f plane, which we will take into account for our numerical analysis. For the strongest bound f ≥ 5.6 TeV, the allowed interval of t β values is 1-9. We will analyze below the dependence on t β of the τ AMDM and AWMDM in the allowed interval.
Mixing between light and heavy neutrinos δ ν : this parameter is experimentally constrained to be small [26], with the corresponding bound being flavor dependent: δ ν e ≤ 0.03, δ ν µ ≤ 0.05, and δ ν τ ≤ 0.09 with 95% C.L. Since we are interested in the study of the τ lepton, we need to make sure that we use values of f and t β consistent with the bound 09, which in turn translates into the bound f t β 1932 GeV. Such a constraint is fulfilled for the values of f and t β chosen in our analysis.
Pseudoscalar mass m η : this parameter is basically dependent on the µ parameter (m η ∼ µ) appearing in the scalar potential via the term −µ 2 (Φ † 1 Φ 2 + H.c.). In our analysis we will explore values consistent with the lower bound m η ≤ 7 GeV, which arises from the non-observation of the ϒ → ηγ decay [35]. Although the dominant contribution to the AMDM can arise from a very light pseudoscalar, with mass of the order of 10 GeV, this requires relatively large values of t β , which are already excluded according to the above discussion [33,34].
Mixing matrix elements V mi : previous studies on LFV within the SLHM [26-28, 36, 37] have parametrized the mixing matrix V and considered bounds on the mixing angles from experimental data on LFV processes. A simple approach was taken by the authors of Ref. [26] in which an scenario with mixing between the first and second families only was considered. It was found that the respective angle is tightly constrained: the limit on e− µ conversion yields an upper bound on sin 2θ 12 of the order of 0.005. As we are interested in the τ AMDM and AWMDM, we will assume an scenario with mixing between the second and third families only, namely we will consider the following mixing matrix In the loop circulate the charged W and X gauge bosons and a heavy neutrino N k .
We will analyze whether current experimental bounds on processes such as the muon AMDM and the τ → µγ decay can be helpful to find a bound on the mixing angle θ . We already have presented the results for the AMDM and AWMDM of a lepton, we will now present the SLHM contribution to the τ → µγ decay in terms of both parametric integrals and Passarino-Veltman scalar functions. The Feynman diagrams inducing this decay at the one-loop level are shown in Fig. 3. The corresponding amplitude can be written, in the limit of massless j , as with q = p i − p j and F L given by for the X gauge boson. The f (x k ) function is presented in Appendix C. We have verified that our results are in agreement with [26]. The i → j γ decay width is given by The current experimental limit BR(τ → µγ) ≤ 4.4 × 10 −8 [6] can translate into a bound on sin 2θ . For this purpose, we introduce the mass splitting δ 23 ≡ z 3 −z 2 with z k = (m N k /m X ) 2 and show in Fig. 4 the contours of the branching ratio of the τ → µγ decay in the δ 23 vs t β plane for z 2 = 1 and f = 2000 GeV. We conclude that even for a large splitting δ 23 , the branching ratio BR(τ → µγ) would hardly reach a level above 10 −8 for sin 2θ of the order of unity, thereby yielding a very weak constraint on this parameter. On the other hand, the muon AMDM also does not yields a useful bound on θ as the heavy neutrino contribution is negative and cannot account for the muon AMDM discrepancy of Eq. (8). As we will see below, the only relevant contribution to the AMDM and WAMDM involving the mixing matrix el-ements is the contribution of the heavy neutrino, which has the generic form It turn out that the second term is subdominant for a small mixing angle θ and a small splitting δ 23 . Following the authors of Ref. [27] in their study of LFV hadronic τ decays, we will consider values of sin 2θ of the order of 10 −1 . Under this assumption, the term proportional to sin 2 θ becomes subdominant, which is equivalent to consider an approximately diagonal mixing matrix. In addition, we have found out that there is little sensitivity of our results to a small change in the value of sin 2θ . A larger value of this parameter would increase slightly the τ AMDM and AWMDM, but an enhancement larger than one order of magnitude would hardly be attained. Heavy neutrino mass: we will follow the approach of Ref. [26], in which m N k is parametrized through the ratio z k = (m N k /m X ) 2 . As observed in Fig. 5, for the most stringent limit f = 5.6 TeV, m N k ≥ 0.836, which corresponds to the value z k = 0.1, whereas m N k ≥ 2.644 TeV for z k = 1. Again, the τ AMDM and AWMDM contributions arising from the heavy neutrinos show little sensitivity to moderate changes in the value of z 3 and the mass splitting δ 23 , so we will use as reference values z 3 1 and δ 23 ≤ 0.5.  Table 1).
In conclusion, we will use the set of values shown in Table 2 for the SLHM parameters involved in our numerical analysis. We have found that except for f and t β there is little sensitivity of the τ AMDM and AWMDM to a mild change in the values of the remaining parameters as far as they lie between the allowed intervals. Table 2 Values used in our numerical analysis for the parameters involved in the AMDM and AWMDM of the τ lepton.

Parameter
Value In order to estimate the τ AMDM we used the Mathematica numerical routines to evaluate the parametric integrals involved in our calculation. A cross-check was done by evaluating the results expressed in terms of Passarino-Veltman scalar functions via the numerical FF/LoopTools routines [38,39]. As already mentioned, it is convenient to analyze the behavior of the AMDM and AWMDM as functions of the symmetry breaking scale f since the mass of the new particles and the corrections to the SM couplings and particle masses depend on it. Also, since the mixing angle δ ν depends on t β , it is worth examining the dependence on this parameter in the allowed interval.

Anomalous magnetic dipole moment of the τ lepton
In the left plot of Fig. 6 we show the absolute values of the main partial contributions to a NP τ along with the total sum as a function of f for t β = 9, whereas in the right plot we set f = 4 TeV and show the dependence of a NP τ on t β . For the remaining parameters we use the values shown in Table  2. We have refrained from showing the curves for the most suppressed contributions.
We first discuss the behavior observed in Fig. 6(a). Notice that the magnitude of each contribution depends highly on the respective f A 1 A 2 A 3 γ coefficient and to a lesser extent on the magnitude of the loop integral, which in turn dictates its behavior. Therefore, the ν τ XX, N τ WW , and ν τ WW contributions, which are not shown in the plot, are the most suppressed ones, with values below the 10 −10 level. This stems from the fact that the f A 1 A 2 A 3 γ coefficients associated with these contributions include two powers of the coupling constants g V n L , which are of the order of δ ν ∼ v/ f , thereby being considerably suppressed for large f . Although the Hττ and Z ττ contributions are less suppressed, they are below the 10 −9 level, whereas the ηττ and N τ XX contributions are the largest ones and can reach values up to the order of 10 −8 for f around 2 TeV, which is a result of the fact that the respective f coefficients have no (v/ f ) 2 suppression factor. Another point worth to mention is that all the partial contributions are negative except for the Hττ and N τ WW ones. Since these contributions are relatively small, they will not interfere with the dominant contributions, which will add up constructively. In conclusion both the ηττ and N τ XX contributions will represent the bulk of the total contribution to a NP τ , which is of the order of 10 −8 for f = 2 TeV, but has a decrease of about one order of magnitude as f increases up to 6 TeV, as observed in the plot.
We now turn to discuss the dependence of a NP τ on the t β parameter as depicted in Fig. 6(b). We observe that the N τ XX contribution, which has a very slight dependence on t β indeed, is the dominant one, with marginal contributions arising from other diagrams. In the allowed t β interval, the Hττ contribution is negligible and is not shown in the plot. For low t β , the ν τ WW and N τ WW contributions can be as large as the N τ XX one, but the ηττ contribution is the one that becomes important when t β increases. Since these contributions are directly proportional to the square of the mixing parameter δ ν = v/( √ 2t β f ), they get suppressed by two orders of magnitude as t β goes from 1 to 9. We observe that the total contribution of the SLHM to a NP τ remains almost constant in this t β interval as it is dominated by the N τ XX contribution.
The behavior of the SLHM contribution to a NP τ is best illustrated in Fig. 7, where we plot the contours of |a NP τ | in the f vs t β plane for the parameter values of Table 2. We observe that the SLHM contribution to a NP τ is of the order of 10 −8 for f between 2 TeV and 5 TeV, but it is below 10 −9 for f above 5 TeV and decreases rapidly as f increases. So, if we consider the most stringent constraint on f , namely 5.6 TeV, we can expect values of a NP τ of the order of 10 −9 . We also observe that there is little dependence of a NP τ on the value of t β , but such dependence is more pronounced for large f . It is interesting to make a comparison with the typical predictions of some popular extension models as reported in the literature. In this respect, several extension models predict values for a NP τ lying in the interval between 10 −9 and 10 −6 [40][41][42][43]. We note that although the SLHM contribution is of the same order than the potential contribution of leptoquark models (LQM) [40], it is disfavored with respect to the contributions of THDMs [29,44], the minimal supersymmetric standard model (MSSM) [45,46], and unparticles (UP) [41], which can reach values as high as 10 −6 .

Anomalous weak magnetic dipole moment of the τ lepton
We now present the analysis of the NP contribution of the SLHM to the AWMDM of the τ lepton. There are some differences with respect to the behavior of the AMDM: apart  Table 2. that the AWMDM receives extra contributions, it can develop an imaginary part, which arises from the diagrams where the external Z gauge boson is attached to a couple of particles whose total mass is lower than m Z . This occurs when there are two internal SM neutrinos or τ leptons in the loop. Again, the magnitude of each partial contribution to the AWMDM is highly dependent on the corresponding f A 1 A 2 A 3 Z coefficient, while its behavior is dictated by the loop integral. We show in Fig. 8  We first discuss the behavior of Re[a W −NP τ ] as a function of f as shown Fig. 8(a). Numerical evaluation shows that for f around 2 TeV the magnitude of the partial contributions ranges from 10 −14 to 10 −9 , with an additional suppression of at least one order of magnitude for f around 5 TeV. The most suppressed contributions are the ν τ XX, Xν τ ν τ and W N τ N τ ones, which are proportional to the (v/ f ) 2 factor and are below the 10 −13 level. Other contributions such as τHZ , ν τ WW , W N τ ν τ , N τ WW , ηττ, Hττ, and XXN τ are less suppressed but are also below the 10 −10 level. In fact, only the contributions shown in the plot, τHZ, N τ XX, and W ν τ ν τ , are relevant for the total sum. While the τHZ contribution is the dominant one, the N τ XX contribution plays a subdominant role, which again is due to the fact that the f N τ XX Z coefficient is not suppressed by the (v/ f ) 2 factor. We also observe that the W ν τ ν τ contribution, which together with the τHZ contribution are absent in the AMDM, is below the 10 −9 level. Very interestingly, the Hττ and ηττ contributions, which were not very suppressed in the AMDM case, are now negligible as they are proportional to the small g Z V coupling. Note also that all the contributions shown in the plot are positive except for the τHZ one. As a result the total contribution will have an additional suppression as the main contributions will have a large cancelation due to their opposite signs. This becomes evident in the curve for the total contribution, which appears below the curve for the τHZ contribution.
As far as the behavior of the real part of a W −NP τ as a function of t β is concerned , we observe in Fig. 8(b) that the τHZ and W ν τ ν τ contributions are highly dependent on t β : in the interval where this parameter increases from 1 to 10, the τHZ contribution is negative and increases by one order of magnitude, whereas the W ν τ ν τ contribution is positive and decreases by one order of magnitude. On the other hand, the N τ XX contribution is positive and remains almost constant throughout this t β interval. It is interesting that due to the opposite signs of the partial contributions, there is a flip of sign of the total contribution around t β = 6.8. For low t β , the total contribution is positive as it arises mainly from the N τ XX and W ν τ ν τ contributions, with the τHZ contribution being subdominant. As t β increases, the magnitude of the τHZ contribution increases, whereas that of the W ν τ ν τ contribution decreases. Therefore, the total sum cancels out around t β = 6.8 and becomes negative above this value since the τHZ contribution becomes the dominant one. This effect is evident in the large dip of the total contribution curve, which is due to the flip of sign of the AWMDM. This behavior of the partial and total contributions around t β = 6.8 is best illustrated in the zoomed region displayed at the bottom-right corner of Fig. 8(b), where we show the τ AWMDM contributions without taking their absolute values.
The behavior of the real part of a W −NP τ as a function of f and t β is best illustrated in Fig. 9, where we show the contours of the real part of the total SLHM contribution to a W −NP τ in the f vs t β plane. It can be observed that the real part of a W −NP τ reaches its largest values, of the order of 10 −9 , for low and high t β , irrespective of the value of f . There is a band centered around t β ∼ 7 where the lowest values of real part of a W −NP τ are reached. Such a band (darkest region) widens as f increases . We observe that for t β around 10, there is a slow decrease of Re[a W −NP W ] as f increases but in general its magnitude is of the order of 10 −10 − 10 −9 . A comparison of the SLHM contribution with typical predictions of some SM extensions allow us to conclude that for f up to 4 TeV, the SLHM contribution can be above the 10 −9 level and can be larger than the contributions predicted by THDMs, the MSSM and UP, which are of the order of 10 −10 . For f 4 TeV, the SLHM contribution decreases and it is expected to be smaller than the values predicted by other extension models.
We now turn to examine the behavior of the imaginary part of the partial contributions of the SLHM to a W −NP τ . We will follow the same approach as that used for the analysis of the real part. There are only four contributions to a W −NP τ that can develop an imaginary part. In Fig. 10(a) we show the behavior of such contributions as a function of f for t β = 9. Contrary to what happens with the real part, the imaginary part of the W ν τ ν τ contribution is the dominant one by far, at the level of 10 −10 , with the imaginary parts of the ηττ, Hττ, and Z ττ contributions suppressed by more than one order of magnitude. Therefore, the imaginary part of a W −NP τ will be completely dominated by the W ν τ ν τ contribution. In fact both the curve for the W ν τ ν τ contribution and the total contribution overlap. In this region of the parameter space of the SLHM, the imaginary part of a W −NP τ is positive, with a magnitude of the order of 10 −10 , which slightly decreases as f increases. As far as the behavior of the imaginary part of a W −NP τ as a function of t β is concerned [ Fig. 10(b)], there is no considerable change in the analysis as in the allowed t β interval the imaginary part of the W ν τ ν τ contribution is dominant, whereas the remaining contributions are negligibly small. For low t β , the imaginary part of the total contribution to a W −NP τ can be of the order of 10 −9 , but it decreases by almost one order of magnitude as t β goes up to 10.
Finally we present the contours of the imaginary part of a W −NP τ in the f vs t β plane in Fig. 11. We observe that the imaginary part of a W −NP τ can be of the order of 10 −9 for low values of t β , irrespective of the value of f . This is also true for t β ∼ 10 and f ≤ 3 TeV, but there is a pronounced decrease of about one order of magnitude for larger values of f . We also can observe that Im[a W −NP τ ] decreases mildly as f increases. As far as the values predicted by other extension models, although the imaginary part of the total contribution of the SLHM is larger than that predicted by type-I and type-II LQMs (of the order of 10 −10 ) it is well below the contributions predicted by THDMs and the MSSM (of the order of 10 −7 ).

Conclusions
In this work we have calculated analytical expressions, both in terms of parametric integrals and Passarino-Veltman scalar functions, for the one-loop contributions to the static anomalous magnetic and weak magnetic dipole moments of a charged lepton in the context of the SLHM. We have considered the scenario in which there is no CP violation in the model and thereby there are no electric nor weak electric dipole moments. The expressions presented for the weak properties are very general and can be useful to compute the weak properties of a charged lepton in other extension models. For the numerical analysis we have focused on the case of the τ lepton since their electromagnetic and weak properties are the least studied in the literature and also because they have great potential to be experimentally tested in the future. For values of the parameters of the model allowed by current experimental data we find that the respective contribution to the τ AMDM is of the order of 10 −9 , whereas the real (imaginary) part of the τ AWMDM is of the order of 10 −9 (10 −10 ). The SLHM contribution to the τ AMDM could have some enhancement in the scenario in which there is a very light pseudoscalar boson η, with a mass of the order of about 10 GeV, and t β is of the order of 20. However such a value of t β are already excluded according to the bounds obtained from experimental data. Proposed future experiments are expected to reach a sensitivity to the τ AMDM and the real part  Table 2. of the AWMDM of the order of 10 −6 [15][16][17] and 10 −4 [21]. Therefore, the values predicted for these observables by the SLHM would be out of the reach of experimental detection.
One further remark is in order here. Little Higgs models are effective theories valid up to the cut-off scale Λ = 4π f . Below this scale a good prediction is provided by the effective theory, but at higher energies the physics would become strongly coupled and the effective theory must be replaced by its ultraviolet (UV) completion, which would be a QCDlike gauge theory with a confinement scale around 10 TeV (see for instance [47]). This gives rise to the possibility that the EWSB is driven by strong dynamics such as occurs in technicolor theories. It is thus possible that the τ AMDM and AWMDM can receive some enhancement from the UV completion, but an analysis along these lines is beyond the scope of the present work.  The same as in Fig. 9 but for the imaginary part of a W −NP τ . Table 3 Feynman rules for the trilinear gauge boson couplings V µ (p 1 )V + ν (p 2 )V − ρ (p 3 ). The Feynman rule for all these couplings have the generic form where the g V i V j V j constants are shown in the right column. All the four-momenta are taken incoming.
and 2g

Appendix B: Loop integrals
The AWMDM and AWMDM of a lepton are given by Eqs. (56) and (59). The f coefficients are shown in Table 6 and the loop functions I A 1 A 2 A−3 V (V = Z, γ) will be presented below. The loop integration was performed via both Feynman parametrization and the Passarino-Veltman method [48]. For the Dirac algebra and the Passarino-Veltman reduction we used the Feyncalc routines [49], and a further simplification was done with the help of the Mathematica symbolic algebra routines. We will first present the results in terms of parametric integrals. Table 5 Feynman rules for the vertices involving the scalar bosons H and η at the leading order in v/ f in the SLHM.  (59), respectively, to compute the AWMDM and AMDM of a charged lepton arising from loops carrying the A 1 A 2 A 3 particles. Here n is the SM neutrino ν or a heavy one N. The coupling constants g A 1 A 2 A 3 are shown in Appendix A.

: Parametric integrals
After introducing Feynman parameters and integrating over the 4-momentum space, the I loop integrals can be cast in the following form after one Feynman parameter is integrated out: We will first present the F A 1 A 2 A 3 Z (x) functions necessary to calculate the AWMDM of a lepton.
Appendix B.1.1: Anomalous weak magnetic dipole moment We start by introducing the functions where we have introduced the dimensionless variable x a = (m a /m V ) 2 .
As for the contribution of diagram (2), the F Z Z function can be written as with y a = (m a /m Z ) 2 and We now present the contribution of diagram (3): Finally, the loop function arising from diagram (4) is given by where the functions f respectively. In addition We now turn to the parametric integrals for the diagrams of Fig. 2. For diagram (1) we obtain with w a = (m a /m H ) 2 and Also, for the contribution of the pseudoscalar η we obtain As for the contribution of diagram (2) and the one obtained by exchanging the internal gauge boson with the scalar boson, it is as follows

.2: Anomalous magnetic dipole moment
For completeness we present the contributions to the AMDM, which can be obtained from the AWMDM results after taking the limit m Z → 0 and substituting the Z coupling constants by the photon ones. The f coefficients of Eq. (59) are presented in Table 6, whereas the respective loop integrals are of the form of (B.1).
As far as Fig. 1 is concerned, there are only contributions from diagrams (1) and (2), but with the external Z gauge boson replaced by the photon. The contribution of diagram (1) can be written as As for the contribution of diagram (2), the F Z γ function is given by an analogous expression to Eq. (B.8) but now and In the limit of small x the integration of the above functions is straightforward and one obtains Finally, there are also contribution of the scalar bosons H and η arising from the diagram (1) of Fig. 2, with the Z gauge boson replaced by the photon. The respective F All of the above results agree with previous calculations of the AMDM of a lepton (see for instance [7,30]).

Appendix B.2: Passarino-Veltman scalar functions
We now present the results for the AWMDM and AMDM of a lepton in terms of Passarino-Veltman scalar functions.

Appendix B.2.1: Anomalous weak magnetic dipole moment
We first introduce the following ultraviolet finite functions given in terms of two-point Passarino-Veltman scalar integrals: and use a shorthand notation for the following three-point scalar functions The I loop functions arising from the diagrams of Fig. 1 are given as follows. For diagram (1) we obtain whereas the loop function arising from diagram (2) is given by a similar expression to Eq. (B.8): and As for diagram (3), the respective loop integral is x n (x Z + 6) − 5x Z − 6) + 2x 2 ∆ 5 + 2 x 2 (2x n + x Z + 12) − x 2x n (x Z − 5) + x 2 n + 8x Z + 17 + x 2 n (x Z + 6) − x n (7x Z + 12) − x 3 + 2 x 2 Z + 4x Z + 3 C 3 . (B.58) Finally, the loop function arising from diagram (4) obeys We now present the loop functions for the diagrams of Fig. 2. We will use the following additional Passarino-Veltman scalar functions and for the contribution of the new pseudoscalar η where ∆ 11 , ∆ 12 , and C 6 are obtained from ∆ 9 , ∆ 10 and C 5 , respectively, after the replacement m H → m η . As for diagram (2) of Fig. 2, it yields x n x (5x + 12) − 7x 2 n x − 9x n + 3x 3 n − x ((x − 12) x + 17) + 6 C 1 + 2 (x n − 1) (x n + x + 2) ∆ 1 + (15 − 4x n ) x + 3 x n + x 2 n − 2 + x 2 ∆ 2 − 2x (x n + x + 2) , (B.71) where the primed scalar functions ∆ i and C i are obtained from the unprimed ones by setting m Z = 0. We note that all the three-point functions C i appearing in the AMDM are of the generic type C 0 (m 2 A , m 2 A , 0, m 2 B , m 2 C , m 2 B ), which can be written in terms of two-point scalar functions as follows [50] with λ (x, y, z) = (x − y − z) 2 − 4yz and ∆ X = B 0 (0, m 2 X , m 2 X ) − B 0 (m 2 A , m 2 B , m 2 C ). As for diagram (2) we obtain a similar expression to Eq. Again the primed scalar functions ∆ i and C i are obtained from the unprimed ones by setting m Z = 0.
Appendix C: Decay i → j γ In this appendix we present the amplitude for the i → j γ decay both in terms of parametric integrals and Passarino-Veltman scalar functions. The contributions arise from the Feynman diagrams of Fig. 3. The decay width is given in Eq. (64). We have obtained the f L (x k ) function [x k = (m N k /m V ) 2 ] appearing in Eq. (63) via the Feynman parameters technique using the approximation of massless final lepton m j = 0. We first define the function , whereas the f L function is given as f L (x k ) = For the sake of completeness we also include the amplitude in terms of Passarino-Veltman scalar functions. The result reads as where we have defined