# Reduction of the proton radius discrepancy by 3

###### Abstract

We show that in previous analyses of electron-proton scattering, the uncertainties in the statistical procedure to extract the proton charge radius are underestimated. Using a fit function based on a conformal mapping, we can describe the scattering data with high precision and extract a radius value in agreement with the one obtained from muonic hydrogen.

1. Two principally different methods are commonly used to determine the proton charge radius . On the one hand, it enters the QED calculations of atomic energy splittings (electronic and muonic Pohl:2010zza () hydrogen) and can thus be obtained from measurements of these. On the other hand, can be obtained from elastic electron-proton scattering. The corresponding cross sections can be parameterized in terms of the electric and magnetic Sachs form factors and , respectively, that depend on the invariant momentum transfer squared . Positive -values refer to the scattering process, negative to annihilation/creation. The reduced cross section, here in the one-photon approximation, describes the deviation from the scattering off a point-like particle:

(1) |

where is the virtual
photon polarization, is the electron scattering angle in the
laboratory frame and , with the nucleon mass.

Both methods refer to the same quantity, the slope of the proton form
factor at the origin:

(2) |

The form factor obtained from the cross sections has to be extrapolated from
the data at lowest momentum transfer to the origin.
The most precise electron-proton scattering data from Ref. Bernauer:2010wm ()
analyzed using spline and polynomial fit functions lead to a proton charge
radius that differs by from the muonic hydrogen
radius of Ref. Pohl:2010zza (), when averaged with measurements in electronic
hydrogen Beyer13 ().

The purpose of this letter is to illustrate that such extrapolations
lack precision in purely statistical analyses with arbi-
trary fit functions. For example, the fit functions quoted in the final results of
Ref. Bernauer:2010wm () are polynomials and splines.
In this letter, we construct a simple function, that
describes the data equally well and
corresponds to a small radius in agreement with the one obtained from
muonic hydrogen spectroscopy.
This function is based on a conformal mapping and thus obeys the analytic structure
of the form factors. The following function maps the cut in the -plane onto
the unit circle in a new variable :

(3) |

where is the lowest singularity of the form factors with the charged pion mass. The Sachs form factors can then be expanded in the new variable :

(4) |

Here, the form factors are normalized to the charge and anomalous magnetic
moment of the proton, respectively.

Conformal mapping techniques are a standard tool in hadron physics. So far,
they have not been applied to the electron-proton scattering data by the
A1-collaboration, the data for this process with the highest quoted
precision. A previous elaborate analysis of world form factor data in a similar approach
was carried out by Hill and Paz Hill:2010yb (). In contrast to their analysis,
we do not constrain the parameters any further to have a most flexible fit
function, which is needed for the statistical reasoning here. Moreover, the results by Hill and Paz refer to older form factor data, that are extracted from cross sections mainly via the Rosenbluth method. We avoid the systematical uncertainties related to this procedure by directly fitting the cross sections. Also, the results by Hill and Paz show a strong ambiguity due to the included fit range. As we have shown before Lorenz12 (), this can be avoided in a full dispersion relation approach, since this makes use of the complete available information on the spectral function. Loose constraints on the coefficients in a -expansion neglect the mass-related information on the spectral function.

We emphasize the mainly illustrative purposes of this work. This means that
the significance of these fits lies in the comparison to the data analysis by
the A1-collaboration Bernauer:2010wm (). To allow for a direct comparison
to that work, we use exactly the same data without further radiative
corrections and with fixed normalization parameters (see the next section for
details). Physically, the main advantage of the
function used here, compared to the polynomials and splines used by the
experimenters, is the correct inclusion of the lowest singularity of the form
factors. Besides the basic analytic structure of the form factors, the
imaginary part of the form factors can be constrained further due to
unitarity, see e.g. Lorenz12 (), but this goes beyond the aim of this work.
However, arbitrary coefficients correspond to an unconstrained spectral function.
Therefore, the fits shown here have to be treated on the same footing as polynomial or spline fits.

2.

To illustrate the comparison to the original analysis Bernauer:2010wm (),
we follow the same procedure to choose the number of parameters. This means,
we increase the number of terms in the expansion Eq. (4) until the
of the fits reach a plateau and the fits stabilize.
We have performed the calculations in python and checked the results with mathematica. The fits are carried out using the lmfit package with several optimization methods as the Levenberg-Marquardt and simulated annealing algorithms fitting ().

The extracted
radii and -values are shown in Fig. 1.
For , the
absolute -value of 1563 is reached. This is exactly the value found
in Bernauer10 () for the best polynomial fit. For , both
electric and magnetic radii start to stabilize. The proton charge radius is
found to be fm, consistent with the muonic hydrogen value
of Ref. Pohl:2010zza () and also the one obtained from a dispersion
theoretical analysis of the Mainz and older data, including the ones for the
neutron Lorenz12 (). The level of variation in the magnetic radii is
larger than in the electric
case, as we expect, since is suppressed by a factor of the momentum
transfer squared in the reduced cross section, cf. Eq. (1).
The proton magnetic radius comes out as
fm, somewhat larger than the value found by the
A1-collaboration and within large variations compatible with the one obtained
in Ref. Lorenz12 ().
The level of precision of the original analysis Bernauer:2010wm () is reproduced here
only with unconstrained parameters. As an example, we give the parameters for the fit: for the electric form factor and
for the magnetic form factor. Clearly, such large coefficients generate an unphysical spectral function and thus are far from realistic, just like polynomial or spline fits. However, the crucial point is that constraints on the coefficients as suggested by Hill and Paz increase the in our fits. This is also the case in a dispersive framework. Currently, there exist no fits in the literature with an equally good to these data that obey all known physics constraints. In order to unambiguously distinguish between purely statistical and theoretically motivated data analyses, here we treat only statistical issues. The consideration of more physically motivated and constrained spectral functions will be covered in a later publication LDHM () using the statistical reasoning and unphysical fits shown here as a prerequisite. In contrast to here, the constraints from a realistic spectral function and the asymptotic behaviour naively expected from quark counting rules will be
treated. The possible impact of data at larger on the radius term is of course an additional source of uncertainty. The main point of this paper is to clearly show that the combined uncertainties of fits without physics input do not allow to distinguish between the ‘small’ and ‘large’ radius.

3.

The cross sections corresponding to this fit are shown in Fig. 2 for the conformal mapping function with . The six data sets for different energy settings of the incoming electron beam are separated by an offset. Note that each of these data sets contains measurements from three different spectrometers. The different experimental settings give rise to a normalization uncertainty between the individual data sets. One can take this into account by 31 floating normalization parameters, as described in the original analysis. According to Ref. Bernauer10 (), the data contains the normalizations determined by the spline fit. In order to show the underestimation of uncertainties in previous analyses, it is sufficient to keep the normalization parameters fixed. However, the additional normalization uncertainty covers an even larger range of radius values than given here when considered in floating normalizations. To be more precise, we exactly fit to the same data as was done by the A1-collaboration just using an alternative fit function. Other issues like an improved treatment of radiative and two-photon corrections are not of relevance for this letter but will be taken up in a later publication LDHM (), as well as the rigorous inclusion of physics constraints. Here, this procedure is necessary for a proper comparison to the Mainz analysis.

4.

As a further check on our fits, we now consider the form factor ratio that has been measured precisely using recoil polarization techniques. The form factor ratio from the illustrative fit with compares well to the recent measurements at Jefferson Laboratory, see Fig. 3. The displayed ratio is very similar to the one obtained in the spline fit in Ref. Bernauer10 (). The ‘wiggle’ they found below GeV from the magnetic form factor is also reproduced. This could be interpreted as a result from overfitting, which would be consistent with the fact that the wiggle vanishes when including further physical constraints, see e.g. Ref. Lorenz12 (). We do not want to enter this issue here in more detail but refer the reader to Ref. Meissner:2007tp () for a discussion. However, under the assumption that all statistical and systematic errors are sufficiently under control, a good data description in terms of a low is required. In this case, one has to consider polynomial, spline and unconstrained conformal mapping fits on the same footing. In principle, the latter is to be preferred due to the requirements from analyticity. Even if one were to neglect this fact, one can see from this work that a disagreement between the proton charge radius extracted from electron-proton scattering data and muonic hydrogen cannot be inferred from polynomial or spline fits, as one neglects a sizeable source of uncertainty.

5. In this Letter, we have reanalyzed the recent elastic electron-proton scattering data from Mainz with a fit function that is sufficiently flexible to describe the data with a given precision, ie. with the same precision as achieved by the experimenters using spline and polynomial fit functions. The results for the proton charge radius are in perfect agreement with the values obtained via a dispersion relation approach Lorenz12 () and the recent muonic hydrogen measurements. The remaining -discrepancy is the deviation between the average of the spectroscopic measurements in electronic hydrogen and those in muonic hydrogen, see e.g. Ref. Beyer13 (). To solve this, further measurements in ordinary hydrogen are under way Beyer13 (). The planned muon-proton scattering experiment MUSE MUSE () might also shed further light on these issues.

###### Acknowledgements.

We thank Jan Bernauer for useful discussions. We also thank Hans-Werner Hammer for collaboration during the initial stage of this project. One of the authors (IL) would like to express a special thanks to the Mainz Institute for Theoretical Physics (MITP) for its hospitality and support. This work is supported in part by the DFG and the NSFC through funds provided to the Sino-German CRC 110 “Symmetries and the Emergence of Structure in QCD”, and the EU I3HP “Study of Strongly Interacting Matter” under the Seventh Framework Program of the EU.## References

- (1) R. Pohl, A. Antognini, F. Nez, F. D. Amaro, F. Biraben, J. M. R. Cardoso, D. S. Covita and A. Dax et al., Nature 466 (2010) 213.
- (2) J. C. Bernauer et al. [A1 Collaboration], Phys. Rev. Lett. 105 (2010) 242001. [arXiv:1007.5076 [nucl-ex]].
- (3) A. Beyer et al., J. Phys. Conf. Ser. 467 (2013) 012003.
- (4) R. J. Hill and G. Paz, Phys. Rev. D 82 (2010) 113005 [arXiv:1008.4619 [hep-ph]].
- (5) I. T. Lorenz, H. -W. Hammer and U.-G. Meißner, Eur. Phys. J. A 48 (2012) 151 [arXiv:1205.6628 [hep-ph]].
- (6) J. C. Bernauer, “Measurement of the elastic electron-proton cross section and separation of the electric and magnetic form factor in the Q range from 0.004 to 1 (GeV/c)”, PhD thesis, Johannes Gutenberg-Universität Mainz, 2010, http://wwwa1.kph.uni-mainz.de/A1/publications/doctor/bernauer.pdf.
- (7) E. Jones, T. Oliphant and P. Peterson, SciPy: Open source scientific tools for Python, 2001, http://www.scipy.org/. M. Newville et al., lmfit-py-package, https://github.com/lmfit/lmfit-py/.
- (8) I. T. Lorenz, Y.-B. Dong, H.-W. Hammer and U.-G. Meißner, in preparation.
- (9) G. Ron et al. [JLab Hall A Collaboration], Phys. Rev. C 84 (2011) 055204. [arXiv:1103.5784v1 [nucl-ex]]
- (10) X. Zhan et al. [JLab Hall A Collaboration], Phys. Lett. B 705 (2011) 59-64. [arXiv:1102.0318v1 [nucl-ex]]
- (11) U.-G. Meißner, AIP Conf. Proc. 904 (2007) 142 [arXiv:nucl-th/0701094].
- (12) MUSE Collaboration (R. Gilman et al.), [arXiv:1303.2160 [nucl-ex]]