Abstract
We report on a comprehensive reinterpretation of the existing cross-section data for elastic electron-proton scattering obtained by the initial-state radiation technique, resulting in a significantly improved accuracy of the extracted proton charge radius. By refining the external energy corrections we have achieved an outstanding description of the radiative tail, essential for a detailed investigation of the proton finite-size effects on the measured cross sections. This development, together with a novel framework for determining the radius, based on a regression analysis of the cross sections employing a polynomial model for the form factor, led us to a new value for the charge radius, which is \((0.878 \pm 0.011_\mathrm {stat.}\pm 0.031_\mathrm {sys.}\pm 0.002_\mathrm {mod.})\,\mathrm {fm}\)
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
The problem of the proton charge radius arose from the significant deviation of the very precise Lamb shift measurements in muonic hydrogen [1, 2], which gave a value of \(0.84087(39)\,\mathrm {fm}\), from the CODATA [3] value of \(0.8751(61)\,\mathrm {fm}\), compiled from electron scattering and the old atomic Lamb shift measurements. This discrepancy has motivated several follow-up measurements, but despite the efforts remains unresolved. While the measurement of the 2S–4P transition in hydrogen [4] and the new scattering experiment at Jefferson Lab [5] yield values of \(0.8335(95)\,\mathrm {fm}\) and \(0.831(14)\,\mathrm {fm}\), which are in agreement with the smaller radius, the measurement of the 1S–3S transition [6] gives a value of \(0.877(13)\,\mathrm {fm}\), and supports previous scattering experiments [7, 8] and the hypothesis of a large proton radius. Therefore, additional experiments, both scattering and spectroscopic, still have the potential to make valuable contributions to the proton size problem [9, 10].
In scattering experiments the charge radius of the proton is traditionally determined by measuring the cross section for elastic scattering of electrons from hydrogen, which depends on \(G_E^p\) and carries information about the charge distribution of the proton. The proton charge radius, \(r_p\), is given by
where \(Q^2\) is the absolute value of the square of the four-momentum transfer to the proton. The accuracy of the radius obtained in this manner is limited by the extent of available data sets, which dictates the approach to the extrapolation of \(G_E^p\) needed to determine the slope at \(Q^2=0\). Hence, to ensure a reliable extraction of the radius, measurements of \(G_E^p\) are needed in the region of \(Q^2\lesssim 0.01\,\mathrm {GeV}^2/c^2\).
Efforts to perform such measurements with the standard approaches are limited by the minimum \(Q^2\) accessible with the experimental apparatus at hand, predominantly due to the restrictions in the available electron beam energy and the minimum scattering angle. Therefore, a new experimental approach based on initial-state radiation has been introduced [11] that allows for cross-section measurements down to \(0.001\,\mathrm {GeV^2}/c^2\) with sub-percent precision by using information about the charge form factor that is implicit in the radiative tail of the elastic peak.
2 Initial-state radiation experiment
The radiative tail of an elastic peak is dominated by the coherent sum of two Bethe-Heitler diagrams [12] shown in Fig. 1. The initial-state radiation diagram (BH-i) describes the process where the incident electron emits a real photon before interacting with the proton. Since the emitted photon carries away a fraction of the incident energy, the momentum transfer to the proton is decreased. Hence, this process probes the proton structure at values of \(Q^2\) smaller than the value fixed by the experimental kinematics and is thus sensitive to the form factors at \(Q^2\) smaller than those corresponding to the elastic setting. On the other hand, the final state radiation diagram (BH-f) corresponds to the reaction where the real photon is emitted after the interaction with the nucleon. Consequently, \(Q^2\) at the vertex remains constant, while the detected four-momentum transfer changes.
In an inclusive experiment \(Q^2\) can not be measured directly, which means that looking only at the data the initial-state radiation processes cannot be distinguished from the final state radiation. Hence, in order to get information on \(G_E^{p}\) at \(Q^2\) smaller than the elastic setting, the data must be studied in conjunction with a Monte-Carlo simulation, which includes a detailed description of internal radiative corrections and considers \(G_E^{p}\) as its free parameter. This is the basic idea of the MAMI experiment, which opened the door of obtaining \(G_{E}^{p}\) down to \(Q^2 \simeq 10^{-4}\,\mathrm {GeV^2}/c^2 \) [11].
The measurement of the radiative tail has been performed at the Mainz Microtron (MAMI) in 2013 using the spectrometer setup of the A1-Collaboration [13]. A rastered electron beam with energies of \(E_0=195\), 330 and \(495\,\mathrm {MeV}\) was used in combination with a hydrogen target, which consisted of a \(5\,\mathrm {cm}\) long cigar-shaped Havar cell filled with liquid hydrogen and placed in an evacuated scattering chamber. For the cross-section measurements the single-dipole magnetic spectrometer B was operated at a fixed angle of \(15.21^\circ \), while its momentum settings were adjusted to scan the complete radiative tail for each beam energy. The central momentum of each setting was measured with an NMR probe to a relative accuracy of \(8\times 10^{-5}\). The spectrometer was equipped with the standard detector package consisting of two layers of vertical drift chambers (VDCs) for tracking, two layers of scintillation detectors for triggering, and a threshold Cherenkov detector for particle identification. The kinematic settings of the experiment were chosen such that the radiative tails scanned at three beam energies overlap.
The beam current was between \(10\,\mathrm {nA}\) and \(1\,\mu \mathrm {A}\) and was limited by the maximum rate allowed in the VDCs (\(\approx 1\,\mathrm {kHz/wire}\)), resulting in raw rates up to \(20\,\mathrm {kHz}\). The beam current was determined by a non-invasive fluxgate-magnetometer and from the collected charge of the beam stopped in a Faraday cup. At low beam currents and low beam energies the accuracy of both approaches is not better than \( 2\,\mathrm {\%}\), which is insufficient for precise cross-section measurements. Hence spectrometer A was used at a fixed momentum and angular setting for precise monitoring of the relative luminosity.
The first analysis of the data, presented in [11], revealed inconsistencies between data and simulation on the order of \(10\,{\%}\) at the top of the elastic peak, see Fig. 2, which led to the omission of the most statistically relevant, elastic data points in the analysed sample. The inconsistency arose due to the incomplete correction for the electron energy losses in the target material. To be able to incorporate the elastic data in the analysis and ensure a more precise extraction of the proton charge radius, the estimations of the external corrections had been investigated in detail.
The external radiative corrections were considered using the formalism of Mo and Tsai [14], while the collisional corrections were approximated by the Landau distribution [15]. The uncertainty of the applied energy corrections was estimated to be smaller than \(1\,\%\) [14]. This was confirmed by the dedicated calibration data collected using plastic (\([\mathrm {CH_2}]_n\)) targets with different thicknesses, which created a perfect testbed for validating the applied corrections, since the spectra differ only in the size of the external energy loss correction. The comparison of the elastic peak shapes for different target thicknesses with the simulations presented the correct scaling of the corrections with the thickness of the target.
The impact of the external radiative and collisional corrections on the shape of the radiative tail is observable only in the first few \(\mathrm {MeV}\) of the radiative tail, see Fig. 2. Hence, an inconsistency between the data and simulation in this region is an indication of an unaccounted for material traversed by the electrons. These inconsistencies are related to the traces of cryogenic deposits on the end caps and side walls of the target cell. They consist mostly of residual nitrogen and oxygen still present in the scattering chamber in spite of the good vacuum conditions (\(10^{-6}\,\mathrm {mbar}\)) [16]. The extra material affects the measured spectra and thus needs to be included in the simulation. This requires knowing the amount of cryogenic deposits at the target walls. The thickness of the depositions on the target entrance window was determined using the nitrogen/oxygen elastic data. Since this wall is exposed to the electron beam, the electrons scattered from the cryogens enter the physics spectra and can be monitored. For this purpose spectrometer A was positioned such that the nitrogen/oxygen elastic lines were inside its acceptance. The collected spectra, together with the known elastic cross sections for these elements were used to determine the thickness of the deposited layer. On the other hand, incident electrons do not scatter from the material on the side walls, hence the cryogens there are not directly detectable by the spectrometers. Furthermore, depositions on the side walls can be much thicker than those on the end caps, because the former are not heated by the electron beam. The amount of cryogens there was estimated by matching the functional dependence of the simulation to the measured elastic spectra, which is uniquely correlated to the thickness of the traversed material, see Fig. 2. The analysis has shown that the layer of cryogens on the side can be as much as 200 times thicker (roughly \(4\cdot 10^{-3}\,\mathrm {g/cm^2}\)) than at the end caps.
With this advancement the agreement between the data and simulation improved significantly and allowed us to include the elastic data in the new analysis presented in this paper. Following the approach described in [11] the full set of 25 data points could then be used to extract the proton-charge form factors for \(0.001 \le Q^2 \le 0.017\,\mathrm {GeV}^2/c^2\). The values and the details of the extraction are presented in Appendix B of this paper.
3 Experimental uncertainties
Although the ISR experiment provides remarkable control over the systematic uncertainties [11], a few ambiguities remain and limit the precision of the measurements. The contributions relevant for the extraction of the proton charge radius include the uncertainty in the relative luminosity (\(0.2\,{\%}\)), the uncertainty in the detector efficiencies \((0.2\, {\%})\) and the contamination coming from the target support frame and the spectrometer entrance flange \((\lesssim 0.7\, {\%})\). The uncertainty of the data associated by the contamination of the spectra with the cryogenic depositions was assessed using a dedicated simulation, normalized to the intensity of the nitrogen, oxygen and Havar elastic lines in the measured spectra. The relative uncertainty of the simulation was estimated to be \(5\,\%\), which leads to a \(0.1\,\%-1.1\,\%\) uncertainty of the measured cross-section ratios. For further details see Appendix A.
4 Parameterisation of the form factor
For \(Q^2<0.02\,\mathrm {GeV}^2/c^2\) it suffices for all practical purposes to parameterise the measured proton charge form factor by using a polynomial of the form
where \(n_{E_0}\) represents the normalisation of the data, \(r_p\) represents the radius and a and b are the higher moments of the model. Although the \(Q^4\)- and \(Q^6\)-terms at \(Q^2<0.02\,\mathrm {GeV}^2/c^2\) account for only about a percent of the form-factor value, they need to be considered in the fit. Due to a strong correlation between radius and a (b), which was estimated to be 0.97 (0.92), a wrong choice of these parameters could shift the value of the radius at the level of \(0.01 (0.003)\,\mathrm {fm}\), see Fig. 3. The data in the available \(Q^2\) range (even with a superior experimental precision) do not permit simultaneous determination of all three parameters [17]. Therefore, a and b need to be taken from the literature. To minimise the bias of the extracted radius, we considered values obtained from the analysis of the available world data [18] which are consistent with the latest experimental results from Jefferson Lab [5]:
5 Original determination of the radius
The prevailing way of determining the proton charge radius is by comparing the measured sets of proton-charge form factors with a selected model. Following the steps in [11] the data were fit by a polynomial (2), by using a common parameter for the radius, \(r_p\), and different normalisation factors, \(n_{E_0}\), one for each energy, see Appendix B for details. In terms of this fit with 21 degrees of freedom, the radius was determined to be \(r_p = (0.873 \pm 0.017_{\mathrm {stat.}} \pm 0.059_{\mathrm {syst.}} \pm 0.003_\mathrm {mod.})\,\mathrm {fm}\,\). The value of \(\chi ^2\), when both statistical and systematical uncertainties are considered, equals 17.7 and remains comparable to the value of 15.8 obtained with the previous dataset [11]. However, even by improving the analysis in the manner described above, the extracted radius still depends critically on the available \(Q^2\) range and the number of fitting parameters. It also remains burdened by large systematic uncertainties of the measured cross-section ratios. The systematic uncertainty is a combination of the experimental uncertainties, presented in Sect. 3, and the uncertainties of the simulation employed to calculate the radiative tail of the elastic peak, which are presented in Appendix A.
6 New extraction of the radius
To further improve the result within the scope of available data, an alternative approach was considered, applicable at the level of measured cross sections. In particular, we studied the ratios of measured and simulated cross sections. First, the cross-section ratios were normalised to the elastic point. To first order the \(E'\) evolution of the rescaled ratios between the data and the simulation at each energy setting depends linearly on the proton charge radius. Furthermore, since all points for a single energy configuration are strongly correlated due to the nature of the experimental approach, the effect of changing the radius appears as a change of the slope of the ratio, \(k(r_p)\). Relying on the chosen model (2), the simulation was performed for different values of \(r_p\) between \(0.76\,\mathrm {fm}\) and \(1.05\,\mathrm {fm}\), and compared to the data; see Fig. 4. The results of the comparison for the two highest energy settings (\(330\,\mathrm {MeV}\) and \(495\,\mathrm {MeV}\)) exhibit a clear dependence and sensitivity to small changes in the proton radius. On the other hand, the \(Q^2\) values of the data at \(195\,\mathrm {MeV}\) are so small that within the measured uncertainties these data alone demonstrate no detectable dependence on the radius and were thus excluded from the analysis.
The best estimate for the proton charge radius should reveal a constant ratio between the data and the simulation. The ratios presented in Fig. 4 indicate that the \(495\,\mathrm {MeV}\) setting favours a radius of \(\approx 0.84\,\mathrm {fm}\), while the \(330\,\mathrm {MeV}\) data suggest \(\approx 1.0\,\mathrm {fm}\). Finding the \(r_p\) at which the simulation matches the data corresponds to finding a point where \(k(r_p)=0\), see Fig. 5. The extracted radii for the two energy settings exhibit a \(2.5\sigma \) tension. This could be a consequence of a statistical fluctuation, but it could also hint at underestimated or missing sources of the systematic uncertainty. These have been studied in detail, see Appendix A, but the discrepancy remained. Consequently, the proton charge radius that agrees with both data sets corresponds to the weighted average of the results for the two beam energies, which is:
Following this approach, the radius is the only free parameter. By investigating the slopes, the normalisations \(n_{E_0}\) disappear from the analysis, resulting in a more robust extraction of the radius, which is less sensitive to the systematic uncertainties than the original result, presented in Sect. 5. The uncertainty of the radius directly follows from the uncertainties of \(k(r_p)\) shown in Figs. 4 and 5. The statistical uncertainty combines contributions of data and simulation, added in quadrature. The number of simulated events was chosen such that the simulated uncertainty is always smaller than the experimental one. The systematic uncertainty is dominated by the pointwise contributions presented in Sect. 3. The model dependent uncertainty is dominated by the uncertainty of the parameter a considered in Eq. 2 to model \(G_E^p\) in the simulation, see Fig. 3.
7 Conclusions
The initial-state radiation experiment at MAMI [11] established a new method for precise investigations of the electromagnetic structure of the nucleon and underlying electromagnetic processes at extremely small \(Q^2\). In this paper we present our findings on the improved data analysis, which revealed the necessity of a complete consideration of cryogens deposited on the liquid hydrogen cell and their influence on the e-p scattering results. The analysis also demonstrated the precision with which these effects could be studied and offered new, improved values of the \(G_E^p\) form factor not accessible in the original work. Furthermore, by studying the slopes of the measured radiative tails relative to the simulated ones, an alternative approach for the extraction of the proton charge radius was developed, which yielded a new best value for the proton charge radius that could be extracted from the ISR data. The new result, see “Cross-section study” point in Fig. 6, is consistent with the value obtained from the original form factor analysis, but is almost three times as precise as the first result [11].
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: The full data set is presented in Appendix B.]
References
R. Pohl et al., Nature 466, 213 (2010)
A. Antognini et al., Science 339, 417 (2013). https://doi.org/10.1126/science.1230016
P.J. Mohr, D.B. Newell, B.N. Taylor, Rev. Mod. Phys. 88, 035009 (2016). https://doi.org/10.1103/RevModPhys.88.035009
A. Beyer, L. Maisenbacher, A. Matveev, R. Pohl, K. Khabarova, A. Grinin, T. Lamour, D.C. Yost, T.W. Hänsch, N. Kolachevsky, T. Udem, Science 358, 79 (2017) https://doi.org/10.1126/science.aah6677. http://science.sciencemag.org/content/358/6359/79
W. Xiong et al., Nature 575(7781), 147 (2019). https://doi.org/10.1038/s41586-019-1721-2
H. Fleurbaey, S. Galtier, S. Thomas, M. Bonnaud, L. Julien, F.m.c. Biraben, F.m.c. Nez, M. Abgrall, J. Guéna, Phys. Rev. Lett. 120, 183001 (2018). https://doi.org/10.1103/PhysRevLett.120.183001.https://link.aps.org/doi/10.1103/PhysRevLett.120.183001
J.C. Bernauer et al., Phys. Rev. Lett. 105, 242001 (2010). https://doi.org/10.1103/PhysRevLett.105.242001
X. Zhan et al., Phys. Lett. B 705, 59 (2011). https://doi.org/10.1016/j.physletb.2011.10.002
R. Pohl, R. Gilman, G.A. Miller, K. Pachucki, Ann. Rev. Nucl. Part. Sci. 63, 175 (2013). https://doi.org/10.1146/annurev-nucl-102212-170627
C.E. Carlson, Prog. Part. Nucl. Phys. 82, 59 (2015). https://doi.org/10.1016/j.ppnp.2015.01.002
M. Mihovilovič et al., Phys. Lett. B 771, 194 (2017). https://doi.org/10.1016/j.physletb.2017.05.031
M. Vanderhaeghen, J.M. Friedrich, D. Lhuillier, D. Marchand, L. Van Hoorebeke, J. Van de Wiele, Phys. Rev. C 62, 025501 (2000). https://doi.org/10.1103/PhysRevC.62.025501
K. Blomqvist et al., Nucl. Instr. Meth. A 403, 263 (1998). https://doi.org/10.1016/S0168-9002(97)01133-9
L.W. Mo, Y.S. Tsai, Rev. Mod. Phys. 41, 205 (1969). https://doi.org/10.1103/RevModPhys.41.205
L.D. Landau, in Collected Papers of L.D. Landau, ed. by D.T. HAAR (Pergamon, 1965), pp. 417 – 424. https://doi.org/10.1016/B978-0-08-010586-4.50061-4
M. Mihovilovič et al., EPJ Web Conf. 72, 00017 (2014). https://doi.org/10.1051/epjconf/20147200017
I. Sick, D. Trautmann, Phys. Rev. C 95, 012501 (2017) https://doi.org/10.1103/PhysRevC.95.012501. https://link.aps.org/doi/10.1103/PhysRevC.95.012501
M.O. Distler, J.C. Bernauer, T. Walcher, Phys. Lett. B 696, 343 (2011). https://doi.org/10.1016/j.physletb.2010.12.067
L.N. Hand, D.G. Miller, R. Wilson, Rev. Mod. Phys. 35, 335 (1963) https://doi.org/10.1103/RevModPhys.35.335. https://link.aps.org/doi/10.1103/RevModPhys.35.335
D. Frèrejacque, D. Benaksas, D. Drickey, Phys. Rev. 141, 1308 (1966) https://doi.org/10.1103/PhysRev.141.1308. https://link.aps.org/doi/10.1103/PhysRev.141.1308
Yu .K. Akimov, et al., Sov. Phys. JETP 35, 651 (1972). [Zh. Eksp. Teor. Fiz. 62, 1231 (1972)]
J.J. Murphy, Y.M. Shin, D.M. Skopik, Phys. Rev. C 9, 2125 (1974) https://doi.org/10.1103/PhysRevC.9.2125. https://link.aps.org/doi/10.1103/PhysRevC.9.2125
F. Borkowski, G.G. Simon, V.H. Walther, R.D. Wendling, Z. Phys. A 275, 29 (1975). https://doi.org/10.1007/BF01409496
G. Simon, C. Schmitt, F. Borkowski, V. Walther, Nucl. Phys. A 333, 381 (1980). https://doi.org/10.1016/0375-9474(80)90104-9
M. McCord, H. Crannell, L.W. Fagg, J.T. O’Brien, D.I. Sober, J.W. Lightbody, X.K. Maruyama, P.A. Treado, Nucl. Instrum. Meth. B 56–57, 496 (1991). https://doi.org/10.1016/0168-583X(91)96079-Z
I. Eschrich et al., Physics Letters B 522, 233 (2001) https://doi.org/10.1016/S0370-2693(01)01285-0. http://www.sciencedirect.com/science/article/pii/S0370269301012850
L.C. Maximon, J.A. Tjon, Phys. Rev. C 62, 054320 (2000). https://doi.org/10.1103/PhysRevC.62.054320
H. Arenhövel, D. Drechsel, Nucl. Phys. A 233(1), 153 (1974). https://doi.org/10.1016/0375-9474(74)90248-6
D. Drechsel, S. Kamalov, L. Tiator, Eur. Phys. J. A 34(1), 69 (2007). https://doi.org/10.1140/epja/i2007-10490-6
F. Borkowski, P. Peuser, G. Simon, V. Walther, R. Wendling, Nucl. Phys. A 222, 269 (1974). https://doi.org/10.1016/0375-9474(74)90392-3
Acknowledgements
The authors would like to thank the MAMI accelerator group for the excellent beam quality which made this experiment possible. This work is supported by the Federal State of Rhineland-Palatinate, by the Deutsche Forschungsgemeinschaft with the Collaborative Research Center 1044, by the Slovenian Research Agency under Grant Z1-7305, by Croatian Science Foundation under the project IP-2018-01-8570 and U. S. Department of Energy under Award Numbers DE-FG02-96ER41003 and DE-FG02-94ER40818.
Funding
Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Klaus Peters.
Now at Stony Brook University, Stony Brook, NY 11794 and Riken BNL Research Center, Upton, NY 11793.
Appendices
Appendix A: Systematic uncertainties
The uncertainties of the cross-section ratios and deduced form factors are dominated by the combination of the experimental uncertainties, presented in Sect. 3, and the uncertainties of the simulation employed to calculate the radiative tail of the elastic scattering process. These need to be studied carefully in order to extract credible results, since some parts depend strongly on the energy of the scattered electron, see Fig. 7. The obtained uncertainties are gathered in Table 1 and discussed in this section.
The simulation exactly calculates the coherent sum of the amplitudes for the leading diagrams shown in Fig. 1. The next-order vacuum polarization diagrams (with electrons inside the fermion loop) are exactly calculable and represent as a multiplicative factor to the cross section. The uncertainty of this correction is given by the size of the omitted vacuum polarization due to muon loops and is smaller than \(0.2\,\%\). The virtual corrections to the Bethe-Heitler diagrams (self-energy corrections and various vertex corrections) require integration of the loop diagrams and are considered as effective corrections to the cross section, using calculations of Ref. [12]. These are accurate to about \(1\,\%\), which in turn corresponds to an error of the cross-section ratios smaller than \(0.016\,\mathrm {\%}\).
The non-virtual second-order correction (emission of two real photons) is approximated by the corrections to the elastic cross section [12, 27]. Relying on the soft-photon peaking approximation the contribution of this correction to the final uncertainty was estimated from the size of the next-order correction and determined to be smaller than \(0.5\,\%\). Hadronic corrections are also considered in the elastic limit using the calculations of Ref. [27]. They contribute only up to \(0.5\,\mathrm {\%}\) to the cross section at lowest energies. In the simulation the proton is always on-shell. Hence, the uncertainty of hadronic corrections arises from the part that is related to the internal structure of the proton, described by the generallized polarisabilities [28]. These effects are tiny at our values of \(Q^2\), and hadronic corrections contribute less than \(\le 5\cdot 10^{-4}\,\%\) to the measured cross-sections. Far away from the elastic peak, one also needs to consider \(\mathrm {H}(e,e')n\pi ^{+}\) and \(\mathrm {H}(e,e')p\pi ^{0}\) reactions, which contribute up to \(10\,\mathrm {\%}\) of all events. These processes were simulated using the MAID model [29] and subtracted from the data. The relative uncertainty of the MAID calculations was estimated to be \(5\,\mathrm {\%}\), and is significant only for the \(495\,\mathrm {MeV}\) setting, where it contributes \(0.5\,{\%}\) to the uncertainty of the measured ratios.
The last important contribution to the systematic uncertainty can be traced back to the formalism of Mo and Tsai [14], employed to describe external corrections to the cross-section, see Sect. 2. These corrections add up to \(0.24\,\%\) to the total uncertainty of the measured cross-section ratios.
The obtained uncertainties determine the errors of the deduced form factors discussed in Sect. 5. In the approximation of small \(Q^2\) these are linearly proportional to the uncertainties of the cross-section ratios, where the proportionality factor, \((1+ \sigma _{\mathrm {BH-f}}/\sigma _{\mathrm {BH-i}})/2\), depends only on the ratio of the initial-state and final-state radiation parts of the cross-section, \(\sigma _{\mathrm {BH-i}}/ \sigma _{\mathrm {BH-f}}\), and approximately equals 5/4.
Appendix B: Proton charge form factors
Table 2 collects the proton charge form factors determined by the Initial-State Radiation experiment for \(0.001 \le Q^2 \le 0.017\,\mathrm {GeV}^2/c^2\). The results were obtained by studying the reaction \(\mathrm {H}(\mathrm {e},\mathrm {e'})\mathrm {p}\gamma \) and comparing the shape of the measured radiative tail with the simulation, using the approach explained in Ref. [11]. The extracted values, presented in the third and fourth column, were then compared to the polynomial Eq. 2. The three data sets were fit with a common parameter for the radius, \(r_p\), but with different renormalisation factors, \(n_{E_0}\), for each energy, disregarding the original normalisations of the data, determined from the analysis of cross-sections. The analysis considered both statistical and systematic uncertainty of the extracted values of \(G_E^p\). In terms of this fit with 21 degrees of freedom and \(\chi ^2\) of 18.3, the normalisations and the radius were determined to be:
The obtained normalisation parameters were applied to renormalise the experimentally determined form factors in order to match the values of the three energy settings and to ensure consistency of the results with the statical limit \(G_E^p(0) = 1\). The corrected values of \(G_E^p\) with the corresponding statistical uncertainties are model dependent and are gathered in the fifth and sixth columns of the table. They are also presented in Fig. 8.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Mihovilovič, M., Achenbach, P., Beranek, T. et al. The proton charge radius extracted from the initial-state radiation experiment at MAMI. Eur. Phys. J. A 57, 107 (2021). https://doi.org/10.1140/epja/s10050-021-00414-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epja/s10050-021-00414-x