TwoLoop Effects in LowEnergy Electroweak Measurements
Abstract
We outline the recent results on the twoloop electroweak contributions to the electronelectron scattering cross sections and asymmetries. Although the twoloop corrections are strongly suppressed relative to the oneloop corrections, they still contribute a few percent to the polarization asymmetry, and even this small contribution cannot be ignored at for ultraprecision experiments such as MOLLER planned at JLab. The NNLO calculation techniques we developed for the electronelectron scattering can be adapted for electronproton processes, electronpositron collisions, and other lowenergy experiments involving leptons.
I Motivation
There are three major ways to look for new physics beyond the Standard Model (SM): the energy frontier (highenergy colliders), the intensity/precision frontier (intense beams) and the cosmic frontier (underground experiments, ground and spacebased telescopes). At the precision frontier, one of the most promising processes is polarized electronelectron (Møller) scattering with parity violation, potentially allowing the indirect detection of hypothetical new physics particles coupling to the SM sector through the kinetic mixing. The first measurement of the Møller scattering cross section was done in 1932 M1932 (), but only the recent improvements in precision allowing to measure the parityviolating leftright asymmetry made Møller scattering a candidate for the newphysics sector detection.
The first observation of parity violation in Møller scattering was made by the E158 experiment at SLAC 2 (), which studied Møller scattering of 45 to 48GeV polarized electrons on the unpolarized electrons in a hydrogen target. Its result at low , E158 () allowed one of the most important parameters in the Standard Model  the sine of the Weinberg angle (the weak mixing angle)  to be determined with an accuracy of 0.5% ( in the scheme).
A recentlycompleted JLab experiment measuring the electronproton scattering asymmetry, QWeak (), aims to determine with relative precision of 0.3%. The results of commissioning run, constituting about 4% of the data collected QWeak2013 (), give the leftright asymmetry of , the smallest and most precise asymmetry ever measured in ep scattering, and lead to the first determination of the weak charge of the proton, , in agreement with the SM prediction of . From the theory aspect, the precision can be improved by the better control of hadronic corrections and accounting for the NNLO contributions to the electron line discussed in this work. With that, and with the full set of data analyzed, has potential to place tight constraints on the possible SM extensions. Another PV ep experiment, P2 proposed for the newlyfunded MESA facility at Mainz, aims to determine even more precisely, to 2%.
The nextgeneration experiment to study ee scattering, MOLLER (Measurement Of a Lepton Lepton Electroweak Reaction) Moller2014 (), planned at JLab following the 11 GeV upgrade, will offer a new level of sensitivity and measure the PV asymmetry in the scattering of longitudinally polarized electrons off an unpolarized target to a precision of 0.73 ppb. That would allow a determination of the weak mixing angle with an uncertainty of JLab12 (), or about 0.1%, an improvement of a factor of five in fractional precision when compared with the E158 measurement. At such precision, any inconsistency with the Standard Model will signal new physics, so the MOLLER experiment, building on the concept of indirect probes, can provide access to physics at multiTeV scales. The experiment will undoubtedly be more challenging than previous parityviolating electron scattering experiments, but the MOLLER collaboration has extensive experience in the similar experiments such as MIT Bates, SLAC E158, JLab G0, HAPPEX, PREX and . The major advantage of the Møller scattering is that ee scattering asymmetry is much less affected by the uncertainties in the hadronic corrections then ep asymmetry, and the gammaZ box radiative correction to PV elastic ep scattering calculated at 11 GeV in N42 () has an accuracy sufficient to keep the uncertainty from this background within the limits of the MOLLER experiment. The rest of the electroweak radiative corrections (EWC), although extensive, can in principle be controlled at sub1% level, with the SM predictions carried out with full treatment of oneloop radiative corrections and at least leading twoloop corrections.
It was repeatedly shown in the literature that even oneloop radiative corrections modify the treelevel prediction for the asymmetry quite significantly (Czar1996 (), Petr2003 (), EM2005 (), ABIZprd (), arx2 ()), so it is essential to have them under a very firm control. (Please see KK2013 () for a review of the lowenergy measurements of the weak mixing angle and additional references.) In ABIZprd (), we found the total correction calculated specifically for 11 GeV ee scattering to be close to , with no significant theoretical uncertainties. A much larger theoretical uncertainty in the prediction for the asymmetry will come from the twoloop corrections, so, for the newgeneration precision measurements, predictions for its scattering asymmetry must include not only a full treatment of oneloop radiative corrections (NLO) but also leading twoloop corrections (NNLO).
We approach the NNLO EWC in stages, by dividing the corrections to the Born () cross section into two classes: the part induced by quadratic oneloop amplitudes (), and the part corresponding to the interference of the Born and twoloop diagrams (). The details of our calculations for the quadratic oneloop amplitudes, the part, are shown in Qpart (), where, following the same approach we used for NLO EWC, we performed a tuned stepbystep comparison between different calculation approaches verifying the results obtained by a semiautomatic approach based on FeynArts, FormCalc, LoopTools and Form with the results from the equations derived by hand. As we found in Qpart (), for the MOLLER kinematic conditions, the part of the NNLO EWC can increase the asymmetry by up to 4%, and depends quite significantly on the energy and scattering angles.
In this paper, we discuss a set of contributions corresponding to the interference of the Born and twoloop diagrams (the Tpart), including the gauge invariant set of boson self energies and vertices of twoloop amplitude , and discuss work still to be done in the future.
Ii Cross Section and Asymmetry
The cross section of polarized Møller scattering with the Born kinematics:
(1) 
can be expressed as:
(2)  
where and is the scattering angle of the detected electron with 4momentum in the centerofmass system of the initial electrons. The 4momenta of initial ( and ) and final ( and ) electrons generate a standard set of Mandelstam variables:
(3) 
is the Born () amplitude shown in Fig.1. The amplitudes (Fig.2) and (Figs.48) correspond to oneloop () and twoloop () matrix elements, respectively.
The oneloop amplitude consists of the boson selfenergy (BSE) (Fig.2a), vertex (Ver) (Fig.2b,c) and box diagrams (Fig.2d,e). We use the onshell renormalization scheme from BSH86 (); Denner (), so there are no contributions from the electron selfenergies.
We present the oneloop complex amplitude as the sum of IR and IRfinite parts . The IRfinite part can be found in Qpart () and for the IR part we have:
(4) 
where is the photon mass and the complex value can be presented in the following form (see, for example, KuFa ()):
(5) 
Analogously, the twoloop amplitude is the sum , where
(6) 
Note that the structure of first term in (6) is the same as in (4) in terms of the soft photon factorization.
To cancel the infrared divergences, we split the differential cross sections into dependent (IRDterms) and independent (infraredfinite) parts:
(7) 
The oneloop cross section is already carefully evaluated with full control of the uncertainties in ABIZprd (). The simplest form for IRDterms are:
(8)  
The imaginary part of the total cross section could be removed in the sum  and parts due to following properties: , and . Thus, in the following sections we can ignore the imaginary part, i. e. and .
Iii Bremsstrahlung for NLO and NNLO
Bremsstrahlung for both NLO (Fig.3a) and NNLO (Fig.3b,c) is needed to cancel the infrared divergences in the oneloop and twoloop amplitudes, correspondingly. (Radiation from only one lepton line is shown in Fig.3, but all four lepton lines are accounted for in our calculations, of course.) To evaluate the cross section induced by the emission of one soft photon with energy less then , we follow the methods of HooftVeltman () (see also KT1 ()). Then, this cross section can be expressed as: , where have the similar factorized structure based on the factorization of the softphoton bremsstrahlung:
(9) 
where
(10) 
The first part of the softphoton cross section, , cancels the IRD at the oneloop order, while the second part, , cancels the IRD at the twoloop order, with half of going to the cancellation of the IRD in the part and the other half going to treat IRD in the part. At last, the cross section induced by the emission of two soft photons with a total energy less then is calculated in Qpart () as:
(11) 
where is a statistical factor and .
Bringing all terms together, we arrive at the result that is free from infrared divergences. For one loop, the logarithms will cancel out:
(12) 
For the second loop, the cancellation proceeds in a more involved way, that is
(13) 
where
(14) 
Iv Numerical Results
For the numerical calculations at the central kinematic point of MOLLER (=11 GeV, ) we use , , and lepton masses as input parameters in accordance with PDG08 (). The effective quark masses which we use for the vector boson selfenergy loop contributions are extracted from shifts in the fine structure constant due to hadronic vacuum polarization =0.02757 jeger (). For the mass of the Higgs boson, we take and for the maximum soft photon energy we use , according to ABIZprd () and 5DePo ().
Let us define the relative corrections to the Born cross section due to a specific type of contributions (labeled by ) as
In the text below the term "part" corresponds to the contributions of a gauge invariant set of the BSE and vertices only. The parityviolating asymmetry is defined in a traditional way:
(15) 
and the relative corrections to the Born asymmetry due to contribution are defined as
In general, corrections from different diagrams are not additive. Their total contribution is given by
(16) 
where summation is performed over the index . There are reasons to believe that the correction is small, but we can not say the same about .
In the table below we bring together all contributions derived for relative corrections to the unpolarized cross section to the asymmetry, including contributions that stem from the gaugeinvariant set of twoloop vertex and bosonselfenergy diagrams. The three dots mean the contribution from the line above, so we progressively add new contributions as we have them ready. So far, as one can see from the table, the part induced by quadratic oneloop amplitudes (), and the contributions to part corresponding to the interference of the Born and twoloop diagrams () considered here shift the result in the same direction. Under the kinematic conditions of the MOLLER experiment, the asymmetry that takes into account the concerted effect of one and twoloop contributions decreases by about 62.7%. For comparison, the oneloop contribution yield a value of 69.5%. Clearly, there is still a lot to be done, and no definite conclusion can be made until all contributions are accounted for, but it looks like if no major cancellations are introduced by the remaining twoloop contributions, the NNLO effect on the PV asymmetry may be more significant that previously believed. Thus, it is safe to say even now that the research program for the MOLLER experiment must include evaluation of the full set of twoloop corrections. Although our numerical calculations are done for the MOLLER the kinematic conditions, the analytics will be directly applicable for the collider experiments, so we assume the NLO and NNLO contributions will be affecting their cross section asymmetry significantly as well, and will need to be evaluated for these measurements in the future.
Type of  Published  
contribution  
NLO  0.1145  0.6953  ABIZprd () 
…+(O+B)/2+Q  0.1125  0.6536  Qpart () 
…+(O+B)/2+BBSE  0.1201  0.6420  YadFiz2013 () 
+VVer+VerBSE  
…+ double boxes  0.1201  0.6534  EurPh2012 () 
…+NNLO QED  0.1152  0.6500  
…+SE and  0.1152  0.6503  
Ver in boxes  
…+NNLO EW Ver  under way 
V Conclusion
As one can see from our numerical data, at the MOLLER kinematic conditions, the part of the NNLO EWC we considered in this work can increase the asymmetry by up to 7%. The  and parts do not cancel each other but, on the contrary, they are adding up to increase the physical PV effect. Clearly, the large size of the investigated parts demands a detailed and consistent consideration of the rest of the part, which will be the next task of our group. Since the problem of EWC for the Møller scattering asymmetry is rather involved, a tuned stepbystep comparison between different calculation approaches is essential. To make sure that our calculations are errorfree, we control our results by comparing the data obtained from the equations derived by hand with the numerical data obtained with a semiautomatic approach based on FeynArts, FormCalc, LoopTools and Form. These base languages have already been successfully employed in similar projects ABIZprd () and arx2 (), so we are highly confident in their reliability.
In the future, we plan to address the remaining twoloop electroweak corrections to match the planned precision of the MOLLER experiment and the possible future experiments at ILC. Clearly, for the electroweak lowenergy experiments briefly outlined in this paper and for other potential future measurements, it is absolutely essential for an excellent control of oneloop and good understanding of twoloop effects KK2013 ().
In this paper, we outlined motivation and summarized some of our work on the twoloop electroweak radiative corrections involving the SM particles. Even if the LHC continues to agree with the Standard Model up to 14 TeV, the MOLLER experiment will continue to look for new physics scenarios that could escape LHC detection, like various hidden weak scale scenarios. If the LHC does observe an anomaly, then MOLLER will have enough sensitivity to provide sufficient constraints to distinguish between the possible new physics scenarios (new massive or supermassive bosons, for example). To have that kind of sensitivity, the MOLLER aims to measure the PV asymmetry predicted within SM to be about 33 ppb with an overall precision of 0.7 ppb. The advantage of trying to access new physics via such lowenergy ee scattering asymmetry is that a purely leptonic PV asymmetry is one of the few observables whose theoretical uncertainties are well under control. There is no significant contribution from the hadronic sector, the SM Higgs mass, one of our input parameters, is known well enough for our needs, and the full set of NLO (oneloop) electroweak radiative corrections, although large, is now known to better than 0.1%. Just a decade ago, such precision would not be feasible. Now, with the recent development in computer algebra and increased accessibility of computing facilities, we can aim to further improve the SM prediction for PV asymmetry by calculating the radiative corrections at the NNLO (twoloop) level. Since the EWC corrections depend quite significantly on the energy and scattering angles, they would need to be evaluated for each experiment specifically. For example, at the MOLLER kinematic conditions, the part of EWC induced by quadratic oneloop amplitudes () will increase the asymmetry up to 4%, but increases dramatically in the higherenergy region Qpart (). This by itself is not a problem, since the part is now know. However, we still far from making the final conclusion on behavior of the part corresponding to the interference of the Born and twoloop diagrams (). So far, dominant twoloops contributions to the PV asymmetry are at the order of 1% and they are coming from (Ver + BSE) and double boxes. As far as we know at the moment, the newphysics particles are not going to contribute significantly enough at two loops to warrant fullscale calculations, but they may contribute quite noticeably at the oneloop level, depending on the SM extension employed.
Vi Acknowledgments
The authors gratefully acknowledge Yu. Bystritskiy and E. Kuraev for their help with this project, and the the Joint Institute for Nuclear Research for hospitality in 2013. This work is supported by the Natural Sciences and Engineering Research Council of Canada and Belarus scientific program "Convergence". AA and SB thank JLab Theory Group for hospitality during their stay in 2014.
References
 (1) C. Møller, Ann. der Physik 14, 531 (1932).
 (2) K. S. Kumar et al., Mod. Phys. Lett. A 10, 2979 (1995); Eur. Phys. J. A 32, 531 (2007); SLAC E158 Collab. P. L. Anthony et al., Phys. Rev. Lett. 92, 181602 (2004) [hepex/0312035].
 (3) P. L. Anthony et al., Phys. Rev. Lett. 95, 081601 (2005).
 (4) R. D. Carlini et al., Qweak Proposal to Jefferson Lab Program Advisory Committee, arXiv:1202.1255 (2007)
 (5) Qweak Collaboration, Phys.Rev.Lett. 111 (2013) 14, 141803
 (6) The MOLLER Collaboration, "The MOLLER Experiment Measurement Of a Lepton Lepton Electroweak Reaction An Ultraprecise Measurement of the Weak Mixing Angle using Møller Scattering", MIE proposal to US DOE (2014), http://http://hallaweb.jlab.org/12GeV/Moller/

(7)
J. Benesch et al.,
www.jlab.org/~armd/moller_proposal.pdf (2008)
 (8) N. L. Hall et al., ADP1326/T846, JLABTHY131821, arXiv:1311.3389
 (9) A. Czarnecki and W. J. Marciano, Phys. Rev. D 53, 1066 (1996) [hepph/9507420].
 (10) F. J. Petriello, Phys. Rev. D 67 (2003) 033006, [hepph/0210259].
 (11) J. Erler, M. J. RamseyMusolf, Phys.Rev. D72:073003 (2005)
 (12) A. Aleksejevs et al., Phys. Rev. D 82, 093013 (2010).
 (13) A. Aleksejevs et al., Phys.Part.Nucl. 44 (2013) 161174
 (14) K.S. Kumar et al., Ann.Rev.Nucl.Part.Sci. 63 (2013) 237267
 (15) A. Aleksejevs et al., Phys. Rev. D 85, 013007 (2012), arXiv:1110.1750 [hepph].
 (16) M. Böhm, H. Spiesberger, and W. Hollik, Fortschr. Phys. 34, 687 (1986).
 (17) A. Denner, Fortsch. Phys. 41, 307 (1993).
 (18) E. A. Kuraev and V. S. Fadin, Yad. Fiz. 41, 733 (1985).
 (19) G. ’t Hooft and M. Veltman, Nucl. Phys. B. 153, 365 (1979).
 (20) E. A. Kuraev, N. R. Merenkov, and V. S. Fadin, Yad. Fiz. 45, 782 (1987).
 (21) C. Amsler et al., Phys. Lett. B 667, 1 (2008).
 (22) F. Jegerlehner, J. Phys. G. 29 101 (2003) [hepph/0104304].
 (23) A. Denner and S. Pozzorini, Eur. Phys. J. C 7, 185 (1999).
 (24) A. Aleksejevs et al., Eur.Phys.J. C72 (2012) 2249
 (25) A. Aleksejevs et al., Phys.Atom.Nucl. 76 (2013) 888900, Yad.Fiz. 76 (2013) 942954