On higher order radiative corrections to elastic electron–proton scattering

QED radiative corrections to elastic electron–proton scattering at low energies are discussed. Corrections to the electron line and effects due to vacuum polarization are computed. Higher order effects are estimated for the conditions of the experiment on the electric and magnetic proton form factors by the A1 Collaboration. Calculations are performed within the next-to-leading approximation. The inclusion of the higher order effects can affect the value of the proton charge radius extracted from the experimental data.


Introduction
First of all, our paper is motivated by recent very accurate experimental measurements of the electron-proton (ep) elastic scattering at the Mainz Microtron (MAMI) [1]. The average point-to-point errors in the cross sections measurement was of the order of a few permille.
Besides extraction of the proton electromagnetic form factors, the experiment managed to define the value of the proton electric charge radius with high precision. It is worth to note that the result for the charge radius extracted from the electron-proton scattering data was found to be inconsistent with the one obtained from muonic hydrogen [2]. The disagreement stimulates theoretical studies aimed at its resolution. In the present paper we are going to discuss several effects which can affect the data analysis of low-energy elastic electron-proton scattering.
Higher order corrections to elastic and inelastic ep scattering were discussed in the literature for a long time. In particular, higher order leading logarithmic corrections were implemented for several event selection cases in the computer code HECTOR [3]. The complete set of the next-to-leading QED corrections to deep inelastic ep scattering was first computed in Ref. [4]. Universal higher order QED corrections to a e-mail: arbuzov@theor.jinr.ru polarized and unpolarized lepton scattering are known up to O α 5 in the leading logarithmic approximation [5,6]. Here we are not going to repeat or extend the existing results, but rather reduce and scrutinize them for a concrete special case.
The high precision of the experimental measurement of the differential cross section provides the clear requirement on the magnitude of effects which should be taken into account. We assume that aiming at the one-permille accuracy, we have to treat systematically all relative corrections being at least of 10 −4 size.

Preliminaries and notation
Let us consider the process e( p 1 ) + p(P 1 ) −→ e( p 2 ) + p(P 2 ) + (nγ, e + e − ). (1) The initial electron energy E 1 = p 0 1 ≡ E is of the order 1 GeV, E m e . The momentum transfer squared Q 2 = −( p 2 − p 1 ) 2 will be taken in the range 0.003 < Q 2 < 1 GeV 2 , which was explored in the experiment. Note that the condition Q 2 m 2 e holds for the whole range. The differential cross section dσ /d e of process (1) will be considered.
One-loop QED corrections to the process under consideration are well known. They are naturally separated into the following parts: -real and virtual corrections to the electron line, -real and virtual corrections to the proton line, -interference of amplitudes of the first two types, -the effect due to vacuum polarization.
The corresponding analytic results were reproduced in [1] 1 Among one-loop corrections, there is still an open discussion as regards the proper treatment of double photon exchange contributions; see e.g. Refs. [7,8] and the references therein. We agree with the importance of this point, but it goes beyond the scope of our present study.
To estimate the numerical effect of radiative corrections one has to take into account concrete experimental conditions. Of course, to get the final answer one should include the corrections into the whole program of the data analysis. But our task here will be just to present analytic results with simple estimates of their impact. So we will simplify the set-up (still following the main features of the experiment): -we assume that the measurement is based on the detection of the final electron energy and momentum, -the electron is detected "bare", i.e. without possible accompanying photons, -there is just a simple cut on the lost energy: p 0 1 − p 0 2 ≥ E where is a dimensionless parameter, 1 and E m e . The typical magnitude of the O (α) corrections to the differential cross section is defined by three major factors: The enhancement by the so-called large logarithm L ≡ ln(Q 2 /m 2 e ) and by the logarithm of the cut-off parameter make the size of the one-loop correction to be of the order of a few percent. Since the experimental uncertainties are well below this order, the one-loop corrections were treated in the data analysis with care; see the details in [1].
The purpose of our paper is to estimate the leading and next-to-leading higher order corrections. We will consider one by one the following higher order contributions: 1. higher order effects in vacuum polarization; 2. cut-off dependence of the photonic corrections; 3. light pair corrections in the leading logarithmic approximation; 4. complete next-to-leading O α 2 L 1 corrections to the lepton line.
As can be seen from the first order, higher order corrections only to the electron line and to vacuum polarization can be numerically important.

Higher order effects in vacuum polarization
Running of the QED coupling constant can be naturally represented as where α(0) ≡ α ≈ 1/137.036. A discussion of the relative size of different contributions to (Q 2 ) for low Q 2 values can be found in Ref. [9]. The magnitude of (Q 2 ) for the range of momentum transfer under consideration is about 0.01. The bulk of the vacuum polarization effect comes from one-loop e + e − pair insertion into the photon propagator, One can note that the O α 2 contribution is of the next-toleading order, since it contains only the first power of the large logarithm L. So it makes only a ∼10 −5 effect well below the precision tag. The resummation of the vacuum polarization effect gives The polarization of the vacuum by virtual μ + μ − pairs is not as large as by the e + e − ones. But in the bulk of the kinematical domain the suppression is only logarithmic. So, has to be taken into account at least in the first order in α.
For Q 2 = 1 GeV 2 it reaches 2 · 10 −3 . Instead of the resummed geometrical series of Eq. (5), the A1 Collaboration in Ref. [1] used exponentiation of the effect of the vacuum polarization by leptons, which is close numerically for the given Q 2 range; see Fig. 2 below.
The hadronic contribution hadr (Q 2 ) is rather small at Q 2 ≤ 1 GeV 2 , but at the right edge it is rising steeply and reaches a few permille. Contributions of tau leptons and electroweak bosons are obviously numerically negligible in our case. More detailed numerical estimates of the vacuum polarization effects will be presented below in Sect. 3.

Cut-off dependence of the photonic corrections
The Yennie-Frautschi-Suura theorem [10] proves that emission of each soft photon can be treated as an independent process. As the result, multiple emission of soft photons can be resummed into an exponent. By construction in the case of independent emission of soft photons, the maximal energy of each photon is limited independently. But in the given experimental set-up, we have a cut-off on the total lost energy. The corresponding effect was considered e.g. in Ref. [11]. For double soft photon emission in gives the following shift: At Q 2 = 1 GeV 2 this leads to a visible relative shift of the cross section of about −3.5 · 10 −3 .
To have the theoretical precision under control we can estimate the effect also for the leading logarithmic photonic correction in the third order. The relative correction reads So, the treatment of the cut-off results in the relative shift of the order which is not small and reaches about 2 · 10 −3 .
In the same way one can verify that the naïve exponentiation leads to a considerable off-set in the fourth order leading logarithmic correction: can be found in Ref. [12]. The proper exponentiation of radiative corrections in the leading logarithmic approximation is based on the exact solution of the renormalization group equation; see [13]. But for the practical application under consideration it is sufficient to compute the effect order by order and keep the theoretical uncertainty under control in this way.

Light pair corrections; the leading logarithmic approximation
The contribution of e + e − pairs can easily be estimated with the help of the leading logarithmic approximation (LLA) in QED [13][14][15]: where the so-called -parts of splitting functions (see e.g. Refs. [12,14]) read Note that in the third order in α we have an effect due to simultaneous (either virtual or soft) radiation of a pair and a photon.
To have a better control on the precision level, we can include also the next-to-leading pair corrections in the order O α 2 L where some enhancement due to the experimental cut-off takes place. The corresponding effect will be estimated below.

Complete next-to-leading logarithmic corrections to the lepton line
In order to control the precision of theoretical estimates we can compute the complete set of next-to-leading order (NLO) corrections to the given process by means of the renormalization group approach to QED [13]. The NLO QED structure functions were first introduced in [16]. The corresponding fragmentation functions were used in [17,18] to evaluate NLO corrections to the muon decay spectrum. Here we can follow the paper [19], where NLO QED corrections were computed in a similar set-up for the case of Bhabha scattering. The relevant photonic and e + e − pair contributions to the NLO electron structure (str) and fragmentation (frg) functions have the form 2 D str,frg ee Explicit expressions for the splitting functions P (n) ee and d 1 can be found in [19]. The master formula for NLO photonic corrections to the elastic electron-proton scattering reads dσ = 1 z dzD str ee (z)(dσ (0) (z) + dσ (1) (z) where dσ (1) is the O (α) correction to the ep scattering with a "massless electron", calculated using the MS scheme to subtract the lepton mass singularities. The energy fraction of the incoming parton is z, and Y is the energy fraction of the outgoing (observed) electron. As concerning the factorization scale, it is natural to choose it to be equal to the momentum transfer: L ≡ ln(Q 2 /m 2 e ). Here we are interested in the contributions due to virtual and soft photons, so both integrals have the same lower limit, being equal to 1 − . First we can perform a convolution of the structure and fragmentation functions entering Eq. (14) with each other, D str ee ⊗ D frg ee (z). If z = 1 − and 1, the result of the convolution gives the probability density to find such a situation where one loses in total, due to photon emission, E beam from the total energy of the process under consideration.
A convolution of the function found above with the Born part of the kernel cross section gives us the corresponding part to the cross section (with the upper limit on the lost energy): where the relevant -parts read (2)) The values of the Riemann zeta function are ζ(2) = π 2 /6 and ζ(3) ≈ 1.202. The -parts of the structure and fragmentation splitting functions (P (1,γ (pair))str ee ) and (P (1,γ (pair))frg ee ) coincide, so the notation is simplified.
Note that by construction in the MS scheme, the complete first order correction is reproduced since The factor 2 before the subtracted term on the right hand side reflects the presence of mass singularities in both the initial and the final state corrections.
3 Numerical results Figure 1 shows different contributions to the vacuum polarization correction This figure was obtained with the help of the Fortran package alphaQED by Jegerlehner [20]. One can see that vacuum polarization by muons and hadrons contributes by up to one percent. That is a rather large effect for the given precision tag. Moreover, the momentum dependence of the total vacuum polarization correction is different from the pure electron one. That can affect the extrapolation procedure which is applied for extraction of the proton charge radius.
As for the hadronic contribution to vacuum polarization: it can be either treated as a part of radiative corrections or as a part of the proton form factor. To our mind, the former treatment has two advantages. First, this contribution is always there as for point-like as well as for non-point-like particles. Second, in higher order corrections it is not factorized out as can be seen already in Eq. (5). At first glance the hadronic contribution should not affect the value of the proton charge radius since it is defined at the zero momentum transfer, where this effect is vanishing. Nevertheless, the effect has a pronounced Q 2 dependence in the explored domain and it certainly affects the extrapolation to the zero  Figure 2 shows the difference between the corrections due to vacuum polarization by electron and muons between the result obtained with the help of the Fortran package alphaQED (taking into account also known 2-loop contributions) and the exponentiated treatment of the effect described in Ref. [1]. One can see that the difference is of the order of 2 · 10 4 , which might be relevant for a better control of systematic errors.
The relative QED corrections to the electron line are presented in Fig. 3. The index i runs over: LLA,pair + δ (4)

Conclusions
We presented results for higher order corrections to elastic electron-proton scattering which can be relevant for modern high-accuracy experiments. The corrections are presented in an analytic form. Numerical results are given for a simplified The quantity (20) plotted in Fig. 3 is an estimate of the effect due to an advanced treatment of higher order corrections to the electron line in the process of ep scattering, which is presented here. We have shown also that an accurate treatment of the vacuum polarization effects is also important for getting a high precision. An adequate treatment of all other relevant effects (double photon exchange, radiative corrections to the proton line, details of the experimental set-up, etc.) is also required.