Two-photon exchange correction to the hyperfine splitting in muonic hydrogen

We reevaluate the Zemach, recoil and polarizability corrections to the hyperfine splitting in muonic hydrogen expressing them through the low-energy proton structure constants and obtain the precise values of the Zemach radius and two-photon exchange (TPE) contribution. The uncertainty of TPE correction to S energy levels in muonic hydrogen of 105 ppm exceeds the ppm accuracy level of the forthcoming 1S hyperfine splitting measurements at PSI, J-PARC and RIKEN-RAL.


I. INTRODUCTION
The first spectroscopy measurements with muonic atoms by the CREMA Collaboration at PSI [1] allowed us to study the proton electromagnetic structure with unprecedented precision. The accurate extraction of the proton charge radius from the muonic hydrogen Lamb shift [1,2] gave the discrepancy to measurements with electrons [3][4][5]; see [2,6] for recent reviews. This problem is known as the proton radius puzzle.
The precise spectroscopy measurements require an improvement in the theoretical knowledge of the radiative corrections. The dominant theoretical uncertainty in the proton size extractions from the Lamb shift is coming from the graph with two exchanged photons.
The new highly precise insights on the proton electromagnetic structure will be obtained by the forthcoming measurements of 1S hyperfine splitting (HFS) in muonic hydrogen with an unprecedented ppm precision by the CREMA [24] and FAMU [25,26] Collaborations as well as at J-PARC [27]. In these experiments, the expected accuracy level is two orders of magnitude smaller than the theoretical knowledge of the TPE correction with 213 ppm uncertainty in the dispersive estimate [28] and 109 ppm in the effective field theory approach [29]. The leading TPE effects of the proton structure in HFS are expressed in terms of the proton spin structure functions and form factors [28][29][30][31][32][33][34][35][36][37][38][39][40][41]. Consequently, the dominant uncertainty from the TPE correction can be reduced by the precise measurements of the proton electric and magnetic form factors in the low-Q 2 region [42] and studies of the proton spin structure functions g 1 and g 2 by EG4, SANE and g2p experiments at JLab [43][44][45].
In Ref. [41], we proved the standard expressions for the TPE correction [28,39] expressing it in terms of the forward lepton-proton scattering amplitudes. With the aim to decrease the uncertainty of the α 5 TPE contribution to HFS, we reevaluate the Zemach, recoil and polarizability corrections expressing the region with small photon virtuality in terms of proton radii [46,47], which was introduced by Karshenboim to constrain the values of the electric and magnetic radii from the atomic spectroscopy measurements, and moments of the spin structure functions. We exploit the elastic proton form factors fit, which is based on the unpolarized and polarization transfer world data [3,4], and the latest parametrization of the proton spin structure functions [48][49][50][51]. Additionally, we express the polarizability correction in terms of the measurable spin asymmetry, which provides a direct relation to the experimental observables.
The paper is organized as follows. We describe the standard framework of the TPE correction to S-level HFS and evaluate the proton and inelastic intermediate states contributions in Sect. II. Afterwards, we present the comparison with previous computations. We give our conclusions with an outlook of the forthcoming 1S HFS measurements in Sect. III.
where M and m are the proton and the lepton masses, m r = M m/(M + m) is the reduced mass, µ P ≈ 2.793 is the proton magnetic moment and α ≈ 1/137 is the electromagnetic coupling constant.
The TPE correction is given by a sum of diagrams with proton and with inelastic intermediate states. Conventionally, it is expressed as a sum of the Zemach correction ∆ Z , the recoil correction ∆ p R and the polarizability correction ∆ pol [52]: with the photon energy ν γ and the photon virtuality are the Dirac, Pauli, Sachs electric and magnetic proton form factors (FFs), g 1 (ν γ , Q 2 ) and g 2 (ν γ , Q 2 ) are the spin-dependent inelastic proton structure functions. The following definitions were introduced: The inelastic threshold is given by ν inel thr = m π + (m 2 π + Q 2 ) / (2M ), with the pion mass m π . In the following sections, we evaluate the contributions of Eqs. (4)-(6) separately performing the low-energy expansion in the region of low photon virtuality.

A. Zemach and recoil correction evaluation
The Zemach correction can be evaluated accounting for the measured values of the proton charge and magnetic radii. We split the Q-integration in the Zemach contribution at the small enough scale Q 0 and exploit the radii expansion at low Q 2 [47], thus with the approximate value Q 0 ∼ 0.1 − 0.2 GeV and the definition of the proton radii: For the numerical evaluations, we exploit the elastic proton form factors fit of Ref. [4], which is based on a global analysis of the electron-proton scattering data at Q 2 < 10 GeV 2 with an account of TPE corrections. The resulting uncertainty is evaluated as a sum of the form factors uncertainties [3,4] for Q 2 > Q 2 0 and radii uncertainties for Q 2 < Q 2 0 in quadrature. We add the point Q 2 = 0 with a zero uncertainty to the fit of form factors.
We select Q 0 = 0.15 GeV in the following and estimate the error due to this choice as a difference between our results with Q 0 = 0.15 GeV and Q 0 = 0.2 GeV. We account for the Q 4 and Q 6 terms, exploiting the chiral perturbation theory expansion coefficients [53], and we add the uncertainty of the higher-order contributions as the difference between the calculation with higher-order terms in expansion and result based on Eq. (8), which contributes 13 ppm to the Zemach correction. We substitute the values of the electric charge radius r e E = 0.879 ± 0.008 fm from the electron-proton scattering data [4] as well as r µH E = 0.84087±0.00039 fm from the muonic hydrogen spectroscopy experiments [2]. For the proton magnetic radius, we choose the extraction of the A1 Collaboration r e M = 0.799 ± 0.017 fm [4] and the later more conservative analysis of Ref. [54] r W M = 0.844±0.038 fm. In the following, we study the systematic uncertainty due to the pure knowledge of the proton radii performing the calculation for all possible combinations of chosen electric and magnetic radii.
As a consistency check, we show the dependence of the Zemach contribution on the splitting parameter Q 0 in Fig. 1. The upper plots with the magnetic radius value of A1 Collaboration are closer to the plateau behavior at small Q 0 , which has to appear for the consistent experimental input. However, neither µH nor the electron-proton scattering charge radius passes this check.
In Table I we provide results for different contributions to Zemach term with the corresponding uncertainties. In the calculation with the magnetic radius r e M of Ref. [4], the main uncertainty comes from the error in the proton magnetic radius and form factors in the region of A1/MAMI data Q 2 (0.6 − 1) GeV 2 . The dependence on the splitting parameter is larger in the case of the electron-proton scattering charge radius corresponding to a better consistency of plots in the right panel of Fig. 1. In the calculation with the magnetic radius r W M of Ref. [54], the uncertainty is dominated by the conservative error estimate of the proton magnetic radius extraction. The error due to the choice of the splitting parameter is also enhanced in the case of larger r M , which can be read off from the lower panel of Fig. 1 representing pure consistency of r W M with other electromagnetic proton properties.  The evaluation with the form factor parametrizations at high-Q 2 of Refs. [55,56] gives the same 36 ppm as the fit of A1 Collaboration from the region Q 2 > 0.6 GeV 2 . However, the earlier parametrization of Ref. [57] results in 39 ppm. The possible 3 ppm error is negligible in the evaluation of the resulting uncertainty in quadrature.
The difference between two results based on µH spectroscopy and electron data in 73 ppm can give a hint on the correct radius value in new HFS measurements with ppm accuracy level. Only an improved precision of the magnetic form factor and radius as well as the reduction of the uncertainty in the polarizability correction will make it possible. We also evaluate the recoil correction ∆ p R and the sum ∆ Z +∆ p R performing the similar radii expansion and present the results in Table IV of Sect. II C. We evaluate the error of ∆ Z + ∆ p R adding uncertainties from the form factors under the integral in quadrature.
For completeness, we provide a detailed study of the recoil correction in Table II. The main contribution and uncertainty come from form factors in the region of A1/MAMI data Q 2 (0.6 − 1) GeV 2 . Other form factor parametrizations at high-Q 2 of Refs. [55][56][57] contribute only around 1 ppm from the region Q 2 > 0.6 GeV 2 within the uncertainty of the estimate in Table II  Additionally, we obtain the precise value for the Zemach radius r Z which is defined as Substituting the electric charge radius from the scattering data and the magnetic radius r e M (r W M ), we evaluate the Zemach radius as r Z = 1.0544 ± 0.0079 fm (r Z = 1.0660 ± 0.0135 fm). With the substitution of the charge radius from the muonic hydrogen spectroscopy and the magnetic radius r e M (r W M ), the Zemach radius is given by r Z = 1.0440 ± 0.0068 fm (r Z = 1.0557 ± 0.0120 fm). The results are in reasonable agreement between each other and with the extractions from atomic spectroscopy of Refs.

B. Polarizability correction evaluation
For the numerical evaluation of the polarizability correction, we subtract the leading moment of the spin structure function g 1 and separate contributions from the structure functions g 1 and g 2 [40]: with the corresponding integrands: where the first moment I 1 (Q 2 ) of the structure function g 1 is given by In order to evaluate the contribution from 4I 1 + F 2 P , we approximate I 1 (Q 2 ) = I 1 (0) + I 1 (0) Q 2 up to Q I 1 = 0.25 GeV with the low-energy constant I 1 (0) = 7.6 ± 2.5 GeV −2 [63].
For larger Q 2 , we exploit the spin structure functions data parametrization of Refs. [48][49][50][51] (JLab parametrization). We show the corresponding Q 2 -dependence of the integrand I I 1 in The low-Q 2 and the larger-Q 2 integrands have an intersection point slightly above the lowest data at Q 2 ∼ 0.05 GeV 2 , which was used in the structure functions parametrization [48][49][50][51]. At smaller values of Q 2 , the uncertainty of the data parametrization rapidly increases, and the integration can give an overestimated value. The recent JLAB data [45] confirms the smaller effective value of I 1 (0) for the parametrization. Consequently, the criterion of the same integrand values and similar uncertainties in both regions helps us to choose We estimate the inaccuracy due to this choice as a difference in the HFS correction between calculations with two splitting parameters: Q I 1 = 0.2 GeV and Q I 1 = 0.25 GeV, which evaluates to 8.6 ppm, and add it in quadrature.
For the remaining polarizability corrections ∆ pol 1 and ∆ pol 2 from the proton spin structure functions we use the JLab parametrization only, which is in fair agreement with the MAID model [64,65] in the region of low Q 2 , see Fig. 3 for details.
We add the uncertainties coming from the Pauli form factor F P [3,4], the spin structure functions g 1 , g 2 and the parameter I 1 (0) in quadrature under the HFS integrand and treat the uncertainties from the two Q-integration regions in ∆ 1 and I I 1 contributions as uncorrelated uncertainties. One of important error sources in the resulting polarizability correction is 47 ppm uncertainty from the error of I 1 (0) .
Additionally, we add the 23 ppm error of the higher-order terms in the low-energy expansion as a difference of our evaluation and the calculation [41] with the replacement of the following leading moments of the spin structure functions at Q 2 < 0.25 GeV 2 : by the low-energy constants [63][64][65][66][67][68]: We present the results for different contributions of Eqs. (11)- (13) to the S-level HFS in µH and compare them to Refs. [37,39,59] in Table III, where for results of Ref. [37] we have accounted for the convention conversion correction of Ref. [28]. Though the contributions from the structure functions g 1 and g 2 are slightly different to previous dispersive evaluations of Refs. [28,37,39], the resulting polarizability correction is in good agreement with the results of Ref. [   Within the dispersion relation approach of Ref. [41], we express the polarizability correction ∆ pol 0 directly in terms of the measurable inclusive inelastic lp cross sections as where σ inel hλ denotes the inclusive inelastic cross section with the incoming lepton (proton) helicity h(λ). The integration starts from the inelastic threshold, i.e. the pion production threshold ω thr = m + m π (2M + 2m + m π )/(2M ). Such an expression can be used for the independent direct evaluation of the polarizability correction from the data input.

C. Comparison with literature
In  [7,28,29,[37][38][39][40]59], where we have subtracted the recoil correction of order α 2 [28,39], the radiative correction to the Zemach contribution [28,29,39] and the convention conversion correction of Ref. [28] when it is needed [37,38]. The absolute value of the Zemach contribution is smaller than results of previous estimates [28,39,61] based on the existing form factors parametrizations before the A1/MAMI data, which has a larger value of the magnetic form factor at low-Q 2 region. The recoil correction is in reasonable agreement with other estimates [28,59]. The polarizability correction is in good agreement with dispersive calculations, though all dispersive results are in contradiction to the ChPT prediction.   We finish the comparison to previous results for the total HFS correction in Fig. 4. The difference from Refs. [28,39] is mainly due to the smaller value of the Zemach radius in our evaluation. The smaller value of the polarizability contribution in Refs. [40,59]

III. CONCLUSIONS AND OUTLOOK
In view of the forthcoming high-precision measurements of the 1S hyperfine splitting in muonic hydrogen with ppm precision level [24,26,27], we provide the corresponding best estimates of the TPE correction in Table V. The uncertainty of our result is 100 times larger than the expected experimental accuracy. The proton state contribution allows one to determine the precise value of the Zemach radius, which shows a consistent within 1-1.5σ dependence on the proton radii. The error of the polarizability contribution is almost two times larger than the uncertainty of the Zemach term in the evaluation with the proton magnetic radius of Ref. [4] and is of the similar size exploiting the radius of Ref. [54]. It is dominated by the pure knowledge of I 1 (0) and spin structure functions g 1 , g 2 . The forthcoming data from EG4, SANE and g2p experiments at JLab on the proton spin structure functions g 1 , g 2 [43][44][45] will improve the knowledge of the polarizability correction. The precise measurements of the proton magnetic form factors at low Q 2 [42] and the reextraction of the magnetic radius [69] will allow us to decrease the uncertainty of the Zemach contribution.
Consequently, after accounting for all corrections at the 1−10 ppm level, the forthcoming measurements can constrain the low-Q 2 proton structure contribution to HFS ∆ structure with the following combination of the radii and I 1 (0) : where h(τ ) = (9 − 4τ ) τ 2 + 15 2 ln