The muon g-2 discrepancy: errors or new physics?

The muon g-2 discrepancy: errors or new physics?

M. Passera*, W. J. Marciano^ and A. Sirlin**

*Istituto Nazionale Fisica Nucleare, Sezione di Padova, 1-35131, Padova, Italy
^ Brookhaven National Laboratory, Upton, New York 11973, USA
** Department of Physics, New York University, 10003 New York NY, USA

Abstract. After a brief review of the muon g— 2 status, we discuss hypothetical errors in the Standard Model prediction that
could explain the present discrepancy with the experimental value. None of them looks likely. In particular, an hypothetical
increase of the hadroproduction cross section in low-energy e+e^ collisions could bridge the muon g— 2 discrepancy, but is
shown to be unlikely in view of current experimental error estimates. If, nonetheless, this turns out to be the explanation of
the discrepancy, then the 95% CL upper bound on the Higgs boson mass is reduced to about 130 GeV which, in conjunction
with the experimental 1 14.4 GeV 95% CL lower bound, leaves a narrow window for the mass of this fundamental particle.

Keywords: Muon anomalous magnetic moment. Standard Model Higgs boson
PACS: 13.40.Em, 14.60.Ef, 12.15.Lk, 14.80.Bn

INTRODUCTION

The anomalous magnetic moment of the muon, a^, is one
of the most interesting observables in particle physics.
Indeed, as each sector of the Standard Model (SM) con-
tributes in a significant way to its theoretical prediction,
the precise a^ measurement by the E821 experiment at
Brookhaven [1] allows us to test the entire SM and scru-
tinize viable "new physics" appendages to this theory [2].

The SM prediction of the muon g—2 is conve-
niently split into QED, electroweak (EW) and hadronic
(leading- and higher-order) contributions: a^^ =

af-°+af^+af-°+a™°. The QED prediction, computed
up to four (and estimated at five) loops, currently stands
at flP = 116584718. 10(16)xlO"''[3], while the EW ef-
fects provide a™ = 154(2) x 10""[4]. The latest calcu-
lations of the hadronic leading-order contribution, via the
hadronic e+e^ annihilation data, are in good agreement:
= 6909(44) X 10-ii[5], 6894(46) x 10-"[6, 7],
6921(56) X 10-"[8], and 6944(49) x 10-'^[9]. The
higher-order hadronic term is further divided into
two parts: a™° a™°(vp) + fl™°(lbl). The first

one, -98 (1) X 10"''[6], is the 0{a^) contribution
of diagrams containing hadronic vacuum polarization
insertions [10]. The second term, also of 0{a^), is the
hadronic light-by-light contribution; as it cannot be
determined from data, its evaluation relies on specific
models. Recent determinations of this term vary between
80(40) xlO-'Hll] and 136(25) x 10-^'[12]. The most
recent one, 110 (40) x 1 0^ " [ 1 3] , lies between them. If we
add this result to a^°, for example the value of Ref. [6]
(which also provides the hadronic contribution to the ef-
fective fine-structure constant, later required for our dis-
cussion), and the rest of the SM contributions, we obtain

fl™ = 116591778(61) x lO^'^ The difference with the
experimental value af = 1 16592080(63) x 10"" [1]
is Afl^ = a™ _ ^sM ^ +302(88) x lO^^', i.e., 3.4c7 (all
errors were added in quadrature). Similar discrepan-
cies are found employing the a^° values reported in
Refs. [5, 8, 9]. For recent reviews of aj^ see [7, 14, 15].

The term 0^''° can alternatively be computed in-
corporating hadronic T-decay data, related to those of
hadroproduction in e+e^ collisions via isospin symme-
try [16, 17]. Unfortunately there is a large difference be-
tween the e+e^- and T-based determinations of a^}'°,
even if isospin violation coiTections are taken into ac-
count [18]. The T-based value is significantly higher,
leading to a small (^ Iff) Aa^ difference. As the e+e^
data are more directly related to the calculation than
the T ones, all recent analyses do not include the latter
Also, we note that recently studied additional isospin-
breaking coiTections somewhat reduce the difference be-
tween these two sets of data (lowering the T-based deter-
mination) [19, 20], and a new analysis of the pion form
factor claims that the t and e+e" data are consistent af-
ter isospin violation effects and vector meson mixings
are considered [21].

The 3.4a discrepancy between the theoretical predic-
tion and the experimental value of the muon g—2 can be
explained in several ways. It could be due, at least in part,
to an error in the determination of the hadronic light-by-
light contribution. However, if this were the only cause
of the discrepancy, fl™°(lbl) would have to move up by
many standard deviations to fix it - roughly eight, if we
use the fl™°(lbl) result of Ref. [13] (which includes all
known uncertainties), and more than ten if the estimate
of Ref. [12] is employed instead. Although the errors as-
signed to a™°(lbl) are only educated guesses, this solu-

tion seems unlikely, at least as the dominant one.

Another possibility is to explain the discrepancy Aa^
via the QED, EW and hadronic higher-order vacuum po-
larization contributions; this looks very improbable, as
one can immediately conclude inspecting their values
and uncertainties reported above. If we assume that the
g—2 experiment E821 is correct, we are left with two
options: possible contributions of physics beyond the
SM, or an erroneous determination of the leading-order
hadronic contribution (or both). The first of these
two explanations has been extensively discussed in the
literature; following Ref. [22] we will study whether the
second one is realistic or not, and analyze its implications
for the EW bounds on the mass of the Higgs boson.

ERRORS IN THE HADRONIC CROSS
SECTION?

The hadronic leading-order contribution a^'^ can be
computed via the dispersion integral [23J

1 r

J 4ml

dsK{s)a{s),

(1)

where a{s) is the total cross section for e+e~ annihila-
tion into any hadronic state, with extraneous QED cor-
rections subtracted off, and s is the squared momentum
transfer. The well-known kernel function K{s) (see [24])
is positive definite, decreases monotonically for increas-
ing s and, for large s, behaves as m^/(3s) to a good
approximation. About 90% of the total contribution to
is accumulated at center-of-mass energies ^ be-
low 1.8 GeV and roughly three-fourths of is covered
by the two-pion final state which is dominated by the
p(770) resonance [17]. Exclusive low-energy e+e^ cross
sections were measured at colliders in Frascati, Novosi-
birsk, Orsay, and Stanford, while at higher energies the
total cross section was determined inclusively.

Let's now assume that the discrepancy Aa^ = —
aSM = +302(88) x lO-^\ is due to - and only to - hy-
pothetical errors in g{s), and let us increase this cross
section in order to raise a^°, thus reducing Aa^. This
simple assumption leads to interesting consequences. An
upward shift of the hadronic cross section also induces
an increase of the value of the hadronic contribution to
the effective fine- structure constant at Mz [25],

Mi

" [ ds

JAml

(2)

{P stands for Cauchy's principal value). This integral
is similar to the one we encountered in Eq. (1) for
a^°. There, however, the weight function in the inte-
grand gives a stronger weight to low-energy data. Let

us define a,- = /^'^2 dsfi{s) o{s) (/ = 1,2), where the up-
per limit of integration is Su < M^, and the kernels are
Ms)=Kis)/i47:') and/2(*) = [M^J{M^,-s)]/{4aK^).
The integrals a, with i = 1,2 provide the contributions
to and Aoi^^{Mz), respectively, from 4m| up to Su
(see Eqs. (1,2)). An increase of the cross section g{s)
of the form Aa{s) = ea{s) in the energy range y/s e
[^/sQ — 5/2, -1-5/2], where e is a positive constant
and 2m,i + 5/2 < y/s^ < — 5/2, increases ai by

Afli(V^o,5,e) = ej^^°^gl^2ta{t^)f{t^)dt. If we as-
sume that the muon g—2 discrepancy is entirely due to
this increase in (J{s), so that Aai(v^Q,5,e) = Aa^, the
corresponding increase in Aal^^{Mz) is

Aa2Wso,5)=Aan " ' .

(3)

'Vso-S/2 ■

The shifts Aa2{y/sQ,5) were studied in Ref. [22] for
several bin widths 5 and central values ^/sq.

The present global fit of the LEP Electroweak
Working Group (EWWG) leads to the Higgs bo-
son mass Mh = 84+26 GeV and the 95% confi-
dence level (CL) upper bound M^^ ~ 154 GeV [26].
This result is based on the recent preliminary top
quark mass M, = 172.4(1.2) GeV [27] and the value
A"had(^z) =0-02758(35) [28]. The LEP direct-search
95%CL lower bound is = 1 14.4 GeV [29]. Although
the global EW fit employs a large set of observables, M™
is strongly driven by the comparison of the theoretical
predictions of the W boson mass and the effective EW
mixing angle sin^0^^' with their precisely measured
values. Convenient formulae providing the Mw and
sin^e^^P' SM predictions in terms of Mh, M,, Ao!^l^{Mz),
and av(Mz), the strong coupling constant at the scale
Mz, are given in [30]. Combining these two predic-
tions via a numerical ;|f^-analysis and using the present
world-average values M„ = 80.399(25) GeV [31],
sin20^^P' = O.23153(16) [32], = 172.4(1.2) GeV [27],
as{Mz) = 0.118(2) [33], and the determination
A«had(^z) = 0.02758(35) [28], we get M„ = 89+^^ GeV
and M™ = 156 GeV. We see that indeed the M^ values
obtained from the Mw and sin^O^^' predictions are quite
close to the results of the global analysis.

The Mh dependence of a^*^ is too weak to provide Mh
bounds from the comparison with the measured value.
On the other hand, Aal^^^{Mz) is one of the key inputs of
the EW fits. For example, employing the more recent (and
shghdy higher) value Aa^'^^(Mz) = 0.02768(22) [6] in-
stead of 0.02758(35) [28], the M« prediction shifts down
to Mh = 881^ GeV and M™ = 145 GeV. In [22] we

considered the new values of Aa^^^^{Mz) obtained shift-
ing 0.02768(22) [6] by Aa2{y/sQ, 5) (including their un-

certainties), and computed the corresponding new val-
ues of M™ via the combined j^-analysis based on the

My, and sm^d^X i^^P^ts (for both Aa<'^^(Mz) and aj^^^
we used the values reported in [6]). Our results show
that an increase Ea{s) of the hadronic cross section (in
\/s € [v^Q — 5/2, + adjusted to bridge the

muon g—2 discrepancy Aa^, decreases M™, further re-
stricting the already narrow allowed region for M^. We
concluded that these hypothetical shifts conflict with the
lower Umit M]P when s/sq > 1 .2 GeV, for values of 5 up
to several hundreds of MeV. In [22] we pointed out that
there are more complex scenarios where it is possible to
bridge the Aa^u discrepancy without significantly affect-
ing M™ , but they are considerably more unlikely than
those discussed above.

If T data are incorporated in the calculation of the
dispersive integrals in Eqs. (1,2), 0"'-° significantly in-
creases to 7110(58) X 10-"[17], a™°(vp) slightly de-
creases to — 101(1) X 10~'n6, 15], and the discrepancy
drops to Aa^ = +89(95) x 10"^', i.e. la. While using
T data almost solves the Aa^u discrepancy, it increases
Aa^'^'^(Mz) to 0.02782(16) [34, 17]. In [34] it was shown
that this increase leads to a low prediction which is
suggestive of a near conflict with M\p, leaving a narrow
window for M^. Indeed, with this value of Ao^^|j(Mz)
and the same above-discussed other inputs of the X^'
analysis, we find an M™ value of only 133 GeV.

Recently computed isospin-breaking violations, im-
provements of the long-distance radiative corrections to
the decay n^n^Vx [19], and differentiation of the
neutral and charged p properties [20] reduce to some ex-
tent the difference between t and e+e^ data, lowering the
T-based determination of a^°. Moreover, a recent analy-
sis of the pion form factor below 1 GeV claims that T data
are consistent with the e+e~ ones after isospin violation
effects and vector meson mixings are considered [21]. In
this case one could use the e+e" data below ^ 1 GeV,
confirmed by the T ones, and assume that Aa^ is accom-
modated by hypothetical errors occurring above ~ 1 GeV,
where disagreement persists between these two data sets.
Reference [22] shows that this assumption would lead to
M™ values inconsistent with M]p.

In the above analysis, the hadronic cross section a{s)
was shifted up by amounts Aa(j) = e(y{s) adjusted to
bridge Aa^y. Apart from the impUcations for Mf,, these
shifts may actually be inadmissibly large when compared
with the quoted experimental uncertainties. Consider the
parameter e = A(7(5') /o(s). Clearly, its value depends on
the choice of the energy range [^/sq — 8/2, ^/vq + ^ /2]
where a{s) is increased and, for fixed it decreases
when 5 increases. Its minimum value, ~ 4%, occurs if
c{s) is multiplied by (1 + e) in the whole integration
region, from Inin to infinity. Such a shift would lead to
70 GeV, well below M^. Higher values of e are

obtained for narrower energy bins, particularly if they
do not include the p-o resonance region. For example,
a huge e ~ 52% increase is needed to accommodate
Aa^ with a shift of a{s) in the region from 2m up to
500 MeV (reducing M^^ to 139 GeV), while an increase
in a bin of the same size but centered at the p peak
requires e ~ 8% (lowering M™ to 127 GeV). As the
quoted experimental uncertainty of a{s) below 1 GeV
is of the order of a few per cent (or less, in some specific
energy regions), the possibility to explain Aa^ with these
shifts Aa(s) appears to be unlikely. Lower values of e
are obtained if the shifts occur in energy ranges centered
around the p-co resonances, but also this possibility looks
unUkely, since it requires variations of a{s) of at least
~ 6%. If, however, such shifts A(7(i) indeed turn out
to be the solution of the Aa^^ discrepancy, then M™ is
reduced to about 130 GeV [22].

It is interesting to note that in the scenario where
Aufi is due to hypothetical errors in (y{s), rather than
"new physics", the reduced M)^ 130 GeV induces
some tension with the approximate 95% CL lower bound
Mh 120 GeV required to ensure vacuum stability
under the assumption that the SM is valid up to the
Planck scale [35] (note, however, that this lower bound
somewhat decreases when the vacuum is allowed to be
metastable, provided its Ufetime is longer than the age of
the universe [36]). Thus, one could argue that this tension
is, on its own, suggestive of physics beyond the SM.

We remind the reader that the present values of
sin^0g^' derived from the lep tonic and hadronic ob-
servables are respectively (sin^e^g"), =0.231 13(21) and

(sin^e^^P')^ = 0.23222(27) [32]. In Ref. [22] we pointed
out that the use of either of these values as an input pa-
rameter leads to inconsistencies in the SM framework
that already require the presence of "new physics". For
this reason, we followed the standard practice of employ-
ing as input the world-average value for sin^Q^^' deter-
mined in the SM global analysis. Since M™ also depends
sensitively on , in [22] we provide simple formulae to
obtain the new values derived from different Mf inputs.

CONCLUSIONS

We examined a number of hypothetical errors in the SM
prediction of the muon g—2 that could be responsible
for the present discrepancy Aa^ with the experimental
value. None of them looks likely. In particular, follow-
ing Ref. [22] we showed how an increase Aa{s) =ea{s)
of the hadroproduction cross section in low-energy e+e~
collisions could bridge Aa^. However, such increases
lead to reduced upper bounds (lower than 1 14.4 GeV
- the LEP lower bound - if they occur in energy regions
centered above 1.2 GeV). Moreover, their amounts

are generally very large when compared with the quoted
experimental uncertainties, even if the latter were sig-
nificantly underestimated. The possibility to bridge the
muon g—2 discrepancy with shifts of the hadronic cross
section therefore appears to be unlikely. If, nonetheless,
this turns out to be the solution, then the 95% CL upper
bound M™ drops to about 130 GeV.

If T-decay data are incorporated in the calculation of
a™, the muon g—2 discrepancy decreases to ~ 1 a. While

this almost solves Aa^j, it raises the value of Aa!^l^{Mz)
leading to M™ = 133 GeV, increasing the tension with
the LEP lower bound. One could also consider a scenario,
suggested by recent studies, where the t data confirm the
e+e^ ones below ^ 1 GeV, while a discrepancy between
them persists at higher energies. If, in this case, Aa^u
is fixed by hypothetical errors above 1 GeV, where
the data sets disagree, one also finds values of M™
inconsistent with the LEP lower bound.

It has been suggested [37] that a P-wave electromag-
netic bound state of tt+tt" (pionium) could enter the dis-
persion relations through 1% mixing with the p in a way
that significantly increases the hadronic contribution to
a^. If so, such a state would give little change to the
Higgs mass determination and would seem to refute our
claims. However, Ref. [37] is in error. The required mix-
ing is actually 0.1, not the erroneous 0.01 claim in [37],
and such large mixing is not possible. The actual effect
of pionium on is neghgible [38].

If the Afl^ discrepancy is real, it points to "new
physics", like low-energy supersymmetry where Aa^u is
reconciled by the additional contributions of supersym-
metric partners and one expects Mh ^135 GeV for the
mass of the lightest scalar [39]. If, instead, the deviation
is caused by an incorrect leading-order hadronic contri-
bution, it leads to reduced M™ values. This reduction,
together with the LEP lower bound, leaves a narrow win-
dow for the mass of this fundamental particle. Interest-
ingly, it also raises the tension with the Mh lower bound
derived in the SM from the vacuum stability requirement.

Acknowledgments. We thank A. Vainshtein for discus-
sions and correspondence. The work of W.J.M., M.P.,
and A.S. was supported by U.S. DOE grant DE-AC02-
76CH00016, E.G. contracts MRTN-CT 2004-503369 & 2006-
035505, and U.S. NSF grant PHY-0758032, respectively.

REFERENCES

L G.W. Bennett et al, Phys. Rev. D 73 (2006) 072003;
Phys. Rev. Lett. 92 (2004) 161802; 89 (2002) 101804; 89

(2002) 129903(E).

2. A. Czarnecki and W.J. Marciano, Phys. Rev. D 64
(2001) 013014; D. Stockinger, J. Phys. G 34 (2007) R45;
D. Nomura, these proceedings.

3. T. Kinoshita, M. Nio, Phys. Rev. D 73 (2006) 013003;

70 (2004) 113001; 73 (2006) 053007; T. Aoyama et al,
Phys. Rev. Lett. 99 (2007) 110406; Phys. Rev. D77 (2008)
053012; S. Laporta, E. Remiddi, Phys. Lett. B 301 (1993)
440; 379 (1996) 283; M. Passera, Phys. Rev. D 75 (2007)
013002; A.L. Kataev, Phys. Rev. D 74 (2006) 073011.

4. A. Czarnecki, W.J. Marciano and A. Vainshtein, Phys.
Rev. D 67 (2003) 073006; D 73 (2006) 119901(E);

A. Czarnecki, B. Rrause and W.J. Marciano, Phys. Rev. D
52 (1995) 2619; Phys. Rev. Lett. 76 (1996) 3267.

5. M. Davier, Nucl. Phys. Proc. Suppl. 169 (2007) 288;
S. Eidelman, Acta Phys. Polon. B 38 (2007) 3015.

6. K. Hagiwara et al, Phys. Lett. B 649 (2007) 173.

7. T. Teubncr, these proceedings.

8. F. Jegerlehner, Nucl. Phys. Proc. Suppl. 162 (2006) 22.

9. J.F.de Troconiz,FJ.Yndurain, Phys.Rev.D71(2005)73008.

10. B. Krause, Phys. Lett. B 390 (1997) 392.

11. M. Knecht, A. Nyffeler, Phys. Rev. D 65 (2002) 73034;
M. Knecht etal, Phys. Rev. Lett. 88 (2002) 71802.

12. K.Melmkov, A. Vainshtein, Phys.Rev.D70 (2004) 113006.

13. J. Bijnens, J. Prades, Mod. Phys. Lett. A 22 (2007) 767.

14. F Jegerlehner, Acta Phys. Polon. B 38 (2007) 3021; The
anomalous magnetic moment of the muon. Springer 2007;
J.P Miller, E. de Rafael and B.L. Roberts, Rept. Prog.
Phys. 70 (2007) 795; M. Passera, Nucl. Phys. Proc. Suppl.
169 (2007) 213; 162 (2006) 242; 155 (2006) 365; J. Phys.
G 31 (2005) R75; K. Melnikov and A. Vainshtein, Theory
of the muon anomalous magnetic moment. Springer 2006.

15. M.Davier,W.J.Marciano, Ann.Rev.Nucl.Sci.54(2004)115.

16. R.Alemany,M.Davier,A.H6cker,Eur.Phys.J.C2(1998)123.

17. M.Davier et al, Eur.Phys.J.C27(2003)497;C31(2003)503.

18. W.J. Marciano, A. Sirlin, Phys. Rev. Lett. 61 (1988) 1815;
A. Sirlin, Nucl. Phys. B196 (1982) 83; V. Cirigliano et al,
Phys. Lett. B 513 (2001) 361; JHEP 0208 (2002) 002.

19. F Flores-Baez et al, Phys. Rev. D 74 (2006) 071301.

20. FV. Flores-Baez et al, Phys. Rev. D 76 (2007) 096010.

21. M. Benayoun et al, Eur. Phys. J. C 55 (2008) 199.

22. M. Passera, W.J. Marciano and A. Sirlin, Phys. Rev. D 78
(2008) 013009.

23. C. Bouchiat, L. Michel, J.Phys.Radium 22 (1961) 121;
M. Gourdin, E. de Rafael, Nucl. Phys. B 10 (1969) 667.

24. S. Eidelman and F Jegerlehner, Z. Phys. C 67 (1995) 585.

25. N. Cabibbo and R. Gatto, Phys. Rev. 124 (1961) 1577.

26. LEP EW Working Group, http://lepewwg.web.cern.ch.

27. Tevatron EW Working Group, arXiv:0808.1089 [hep-ex].

28. H. Burkhardt, B. Pietrzyk, Phys. Rev. D72 (2005) 057501.

29. R. Barate et al, Phys. Lett. B 565 (2003) 61.

30. G.Degrassi et al, Phys. Lett. B418(1998)209; G.Degrassi,
P.Gambino, Nucl. Phys. B567(2000)3; A. Ferroglia etal,
Phys.Rev.D65(2002)l 13002; M.Awramik et a/.,Phys.Rev.
D69 (2004) 053006; Phys. Rev. Lett. 93 (2004) 201805;
M.Awramik, M.Czakon, A.Freitas, JHEP0611(2006)048.

31. J. Alcaraz et al, arXiv:hep-ex/0612034; Tevatron EW
Working Group, arXiv:0808.0147 [hep-ex].

32. M. Griinewald et al, Phys. Rept. 427 (2006) 257.

33. C. Amsler et al [PDG], Phys. Lett. B 667 (2008) 1.

34. W.J. Marciano, arXiv:hep-ph/041 1 179.

35. G.Altarelli, G.Isidori, Phys. Lett. B337 (1994) 141; LA.
Casas, et al, Phys.Lett.B342(1995)171; B382(1996)374.

36. G. Isidori et al, Nucl. Phys. B 609 (2001) 387.

37. R.W. Schiel and LP Ralston, Phys. Lett. B 657 (2007) 43.

38. A. Vainshtein, private communication. Prof. Vainshtein
reaches the same conclusions as us.

39. G. Degrassi et al, Eur. Phys. J. C 28 (2003) 133.