Forward two-photon exchange in elastic lepton–proton scattering and hyperfine-splitting correction

We relate the forward two-photon exchange (TPE) amplitudes to integrals of the inclusive lepton–proton scattering cross sections. These relations yield an alternative way for the evaluation of the TPE correction to hyperfine-splitting (HFS) in the hydrogen-like atoms with an equivalent to the standard approach (Iddings, Drell and Sullivan) result implying the Burkhardt–Cottingham sum rule. For evaluation of the individual effects (e.g., elastic contribution) our approach yields a distinct result. We compare both methods numerically on examples of the elastic contribution and the full TPE correction to HFS in electronic and muonic hydrogen.


Introduction
The two-photon exchange (TPE) corrections could be responsible for the discrepancy in the ratio of the proton electric to magnetic form factors between the traditional Rosenbluth separation [1] and the polarization transfer [2] methods; cf. [3][4][5][6] for reviews. These corrections introduce the largest hadronic uncertainty in the elastic electron-proton scattering experiments and should be precisely accounted for in analysis of the modern data. Dedicated experiments to measure the TPE correction from the ratio of the elastic electronproton to positron-proton scattering cross sections, where TPE enters with different signs, have recently been carried out at VEPP-3 (Novosibirsk) [7], OLYMPUS (DESY) [8], and CLAS (JLab) [9]. These measurements show evidence of a significant TPE effect.
The TPE correction also play a prominent role in evaluation of the proton structure (finite-size) contributions to the spectrum of hydrogen (H) and muonic hydrogen (μH). The uncertainty of TPE is the dominant uncertainty in the precise spectroscopy measurements with the muonic hydrogen. These corrections received a renewed attention in light of the so-called proton radius puzzle, the discrepancy in the extracted proton charge radius from the muonic hydrogen Lamb shift [10,11] and measurements with electrons [12][13][14]; see [11,15] for recent reviews. One of the recent achievements in this field is the measurement of the 2S hyperfine splitting (HFS) in muonic hydrogen by the CREMA Collaboration at PSI [11]. Also a high-precision measurement of 1S HFS in muonic hydrogen with an unprecedented ppm accuracy level is planned by this [16,17] and other collaborations [18,19]. Such accurate measurements will strictly constrain the elastic and inelastic proton structure.
To evaluate the contribution from the diagrams with two exchanged photons to atomic spectrum one notices that the lower part of the TPE graph, see Fig. 1, corresponds with the process of forward doubly virtual Compton scat- tering (VVCS). The latter in the forward kinematics can be expressed through dispersion relations (DRs) in terms of the proton structure functions (SFs) measured in elastic and inelastic electron-proton (ep) scattering. The forward TPE is thus evaluated, performing the Wick rotation, as an integral over the photon energy, ν γ , and virtuality, Q 2 > 0. For the Lamb shift correction evaluation [20][21][22][23][24][25] the subtraction function is needed in addition. However, the latter can be estimated at low virtualities from the chiral perturbation theory [26][27][28] or non-relativistic quantum electrodynamics (NRQED) [29] and at high virtualities from the operator product expansion [30]. This function, in principle, can also be determined with account of the high-energy SFs data [31][32][33]. The leading in the electromagnetic coupling constant TPE effects of the proton structure in HFS can be entirely expressed in terms of SFs exploiting the formalism of the projection operators on the singlet and triplet states [34][35][36][37][38][39][40][41][42][43]. For recent numerical evaluations of the TPE correction to HFS see Refs. [41][42][43], for the model-independent evaluation within frameworks of NRQED and ChPT exploiting the hydrogen HFS measurement see Ref. [44] and for results in chiral EFT see Ref. [45].
In this work we provide a different way to express the forward TPE contributions in terms of the proton SFs. We are working on the level of the forward lepton-proton scattering amplitudes and account for the forward double spinflip amplitude (i.e., the one where the helicities of both lepton and proton are flipped) for the first time. The method we are using to derive these relations is akin to deriving sum rules for Compton or light-by-light scattering [46,47]. For the scattering of two charged particles the derivation changes quite a bit. The crossing relates the amplitude of the particle scattering with the amplitude of the antiparticle scattering [48,49]. Nevertheless, for the TPE amplitudes we exploit the crossing in one channel due to the charge independence of the TPE contributions. In this way our dispersion relations (DRs) for the elastic lepton-proton scattering do not involve the crossed channels, and hence are different from generic DRs for charged particles, such as the DRs for the nucleon-nucleon scattering [48][49][50]. Advantageously, the method based on DRs for lp amplitudes [51][52][53] does not require one to evaluate the poles contributions arising after the Wick rotation for amplitudes above the threshold.
We express the TPE corrections to the Lamb shift and HFS of S energy levels through the forward TPE amplitudes at threshold. With DRs for the lepton-proton amplitudes we are not able to express the correction to the Lamb shift through the experimental information, but the correction to the hyperfine splitting is entirely expressed through the proton spin SFs. The resulting HFS correction agrees with the standard approach of Iddings, Drell and Sullivan et al. [35][36][37][38][39][40] only after account for the Burkhardt-Cottingham (BC) sum rule [54]. However, the contribution of each individual channel to the TPE correction in this work differs from the literature result, which can be also obtained exploiting the DRs for the VVCS amplitudes. Afterwards, we compare the proton state and total TPE correction to HFS in H and μH exploiting the DRs for the lepton-proton amplitudes and for the VVCS amplitudes. We connect the region with small photons virtualities in HFS integrand that we express in terms of moments of the spin SFs to the region with large photons virtualities, where the two methods give the same integrands.
The paper is organized as follows. We write down DRs for the forward elastic lp scattering TPE amplitudes in terms of the inclusive lp cross sections at leading α order in Sect. 2 and verify them in QED in Appendix B. We relate the inclusive lp cross sections and the forward TPE amplitudes to proton SFs in Sect. 3. Subsequently, we derive the leading in the electromagnetic coupling constant proton structure corrections to S energy levels coming from TPE and compare the unitarity-based method with the standard approach in Sect. 4. The numerical comparison of the proton and inelastic intermediate states TPE contributions to HFS within two methods is given in Sect. 5. We give our conclusions in Sect. 6. We also provide the derivation of the TPE correction to HFS using the forward Compton scattering tensor in Appendix D. In Appendix E we classify all possible elastic scattering polarization observables in the forward kinematics.

Forward lepton-proton scattering
The forward elastic lp scattering is described by three [55,56] non-vanishing independent helicity amplitudes T h λ hλ , with h (h ) and λ (λ ) the helicities of the initial (final) lepton and proton, see Fig. 1 for the notations of kinematics and helicities: which are functions of the lepton energy ω in the lab frame.
We denote the positive helicity as + and the negative helicity as −. The contribution with exchange of a fixed number of photons to the following three amplitudes has definite evenodd properties with respect to the crossing ω → −ω: The Lorentz structure of the forward amplitude is then given by with m and u the lepton mass and spinor, k the lepton momentum in the lab frame, M and N the proton mass and spinor, where γ μ are the Dirac matrices, the spinors are normalized as In order to establish the even-odd properties for the invariant amplitudes under ω → −ω, we first perform the crossing on the lepton line and relate amplitudes of the lepton-proton scattering f l − p (ω) in the physical region (ω > 0) to amplitudes of the antilepton-proton scattering f l + p (−ω) in the unphysical region (ω < 0): where ω is treated as a complex variable; see Appendix A for details of this derivation. The perturbative contributions with odd number of photons connected to the lepton (antilepton) line have different sign in the amplitudes of the lepton-proton and antilepton-proton scattering as compared to the contributions with an even number of photons, which have the same sign. We express the scattering amplitudes in terms of the contributions with even f (2n)γ (ω) and odd f (2n−1)γ (ω) number of photons connected to both lepton and proton lines, e.g.: and we obtain the following crossing relations for the contribution of graphs with n exchanged photons on the real ω axis: As usual, the optical theorem establishes the relations between the imaginary parts of the forward amplitudes and the total inclusive cross sections of lp collisions: where σ hλ is the inclusive cross section with the incoming lepton helicity h and the incoming proton helicity λ; σ ⊥ (σ ) is the inclusive cross section with lepton and proton polarized transversely and perpendicular (parallel) to each other. 1 The elastic (proton) contribution to the inclusive cross section is infrared divergent. This divergence should be subtracted in a proper way in all three amplitudes. We realize this subtraction only for the case of amplitudes at threshold in Sect. 4.
We express the latter cross sections in terms of the proton SFs up to the order α 2 in Sect. 3 and obtain the imaginary parts of the TPE amplitudes at the leading α order with Eqs. (12)- (13).
Assuming analyticity of these amplitudes in the entire complex ω-plane, except for the branch cuts along the real axis extending from threshold to infinity, see Fig. 2, we can write down the standard DRs: where − stands for the principal-value integration. Using properties of the TPE amplitudes under the crossing ω → −ω, see Eqs. (9)- (11), and assuming the sufficient highenergy behavior of the TPE amplitudes in order to neglect the contribution of the contour at infinity in the Cauchy integral formula, we write down DRs valid at the leading α order and thereby account for both direct and crossed graphs: with the lepton momentum in the lab frame | k | = √ ω 2 − m 2 . These DRs are written for amplitudes in one channel contrary to the DRs for the forward proton-proton scattering [48].
The high-energy behavior of the total unpolarized inclusive cross section does not allow to neglect the contribution of the contour at infinity and to write down Eq. (15). One should perform a subtraction in the DR for the amplitude f 2γ + , e.g., at the point ω s : 2 The lepton-proton DRs were checked at the leading QED order, see Appendix B for details. However, the Regge behavior makes the inclusive cross section σ divergent due to the virtual photons with high energy in the lab frame. Consequently the DRs of Eqs. (15) and (18) [47].
We decompose the forward scattering amplitudes into a sum of one-photon exchange (OPE) and TPE contributions. The TPE amplitudes, except for f 2γ + , are obtained with 2 In the language of effective field theories this subtraction corresponds to the counter term, which was studied in context of the atomic physics and muon-proton scattering in Ref. [58].
DRs and unitarity relations as described above. The OPE amplitudes are real and given by with the proton magnetic moment μ P ≈ 2.793. The vacuum polarization correction is zero in the forward scattering. The lepton vertex correction does not change the amplitude g, modifies the amplitude f − by the lepton anomalous magnetic moment a l : δ f − = a l e 2 and contributes to the amplitude f + . The proton vertex correction with the proton intermediate state contributes to f − and f + amplitudes [59]. The contribution of inelastic intermediate states is expected to be small. However, this correction requires an additional theoretical investigation. Therefore, the forward amplitudes are completely expressed in terms of the total inclusive cross sections, see Eqs. (12), (13), (16), (17), at O(α 2 ) up to the one unknown spin-independent amplitude f + .

Relation of the forward two-photon amplitudes to the proton structure
In this section, we first express the total inclusive cross sections in terms of the experimentally measured proton SFs. Exploiting these relations, we express the real parts of the forward TPE amplitudes as integrals over the photon energy ν γ in the lab frame and photon virtuality Q 2 .
A common way to express the differential inelastic e − p scattering cross section in terms of the proton structure assumes the exchange of one photon. The cross section is given by the contraction of the leptonic tensor L μν and the hadronic tensor W μν [60]. It is proportional to the phase space of the final lepton (with 4-momentum k = ω , k in the lab frame) and given by with the unit of electric charge e. The kinematics are traditionally described by the kinematical Bjorken variable x Bj , the variable y related to the energy transferred by the virtual photon relative to the beam energy and the momentum transfer Q 2 : The leptonic tensor is evaluated in QED. It is given by where s μ is the lepton spin vector: The general Lorentz and gauge-invariant structure of the hadronic tensor W μν , which preserves parity and charge conjugation invariance is given by with the virtual photon energy in the laboratory frame ν γ = ( p · q)/M and the proton SFs , which are extracted from the experimental data. The proton spin 4-vector satisfies S 2 = −1, (S· p) = 0. The total unpolarized cross section is given by This expression reduces to the well-known expression [60] in the massless limit. Consider a scattering of longitudinally polarized leptons on the proton polarized in the lepton momentum direction σ hλ = σ +− and the scattering on the proton polarized in the opposite direction σ hλ = σ ++ with the proton (lepton) spin vector in the laboratory frame S μ = (0, −λˆ k) (s μ = (| k|, ωˆ k) /m) andˆ k = k/| k|. For the cross sections difference we obtain This expression reduces to the well-known expression [60,61] in the massless limit.
Consider the scattering of transversely polarized leptons on transversely polarized protons. Denoting the averaged over the azimuthal angle cross section σ ⊥ (σ ) for scattering with perpendicular (parallel) spin vectors of lepton (s μ = (0, cos φ l , sin φ l , 0)) and proton (S μ = (0, cos φ P , sin φ P , 0)), 3 i.e. φ l − φ P = ±π/2 (φ l − φ P = 0) for the perpendicular (parallel) configuration, we obtain Consequently, a measurement of the inclusive e − p cross sections accesses the proton spin SFs g 1 and g 2 . The elastic scattering cross sections l − p → l − p are obtained by substitution of the inelastic SFs in Eqs. (25)- (27) by the elastic contribution to them: where  17), changing the integration order, as detailed in Appendix C, and expressing the spin-dependent forward TPE amplitudes in terms of the proton SFs, we obtain where we have introduced the notation In Eqs. (32) and (33) the elastic threshold ν thr and the inelastic threshold ν inel thr are given by ν thr = 0 and ν inel thr = m π + (m 2 π + Q 2 )/ (2M) respectively, where m π denotes the pion mass.
The leading TPE correction to the atomic energy levels is given by the values of the amplitudes at threshold ω = m. The TPE amplitudes f 2γ − , g 2γ at threshold can then be expressed in terms of the proton spin SFs as with Evaluating the sum of the spin-dependent lepton-proton TPE amplitudes in the DR approach, see Eqs. (37) and (38), we obtain which is a trivial relation due to the Burkhardt-Cottingham (BC) sum rule [54]:

Lamb shift and hyperfine splitting (HFS)
In this section, we express the TPE correction to S energy levels in the hydrogen-like atom in terms of the forward TPE amplitudes at threshold and compare our result with the standard approach.
In the lp center-of-mass (c.m.) reference frame the TPE forward scattering amplitude T 2γ is expressed in terms of the forward amplitudes f 2γ with s (S) andˆ k (ˆ p) the lepton (proton) spin and momentum direction vectors. This decomposition often arises in the analysis of the non-relativistic forward neutron-proton scattering; see e.g. [50]. It is then easy to see that, considering T 2γ as correction to the Coulomb potential, its effect on the nS-state energy level is given by with |ψ nS (0)| 2 = α 3 m 3 r /(π n 3 )-the non-relativistic squared wave function of the hydrogen atom, where m r = Mm/(M + m) is the reduced mass of the lepton and proton bound state. Using the DR for f 2γ + with only the elastic part of the unpolarized cross section and subtracting the accounted TPE contribution in the hydrogen wave functions, as well as the OPE finite-size correction, we reproduce the non-relativistic limit of the TPE contribution: This correction yields the third Zemach moment [45].
Considering the spin part of T 2γ as correction to the Hamiltonian of the lepton-proton spin-spin interaction, we express the leading TPE proton structure correction to the S-level HFS in terms of the amplitudes f 2γ − , g 2γ at threshold (ω = m): The TPE correction to the S-level HFS 0 of Refs. [35][36][37][38][39][40][41][42][43][44], which was derived from the difference between the spinaveraged forward matrix elements in the singlet and triplet states with a subsequent loop integration, can be also obtained adding g 2γ (m) + f 2γ − (m) = 0 to Eq. (47): Consequently, we have verified the TPE correction to HFS of S energy levels: 0 = HFS . The HFS correction of Refs. [35][36][37][38][39][40][41][42][43] can be reproduced performing the loop integration for the amplitudes g 2γ (m) and f 2γ − (m) expressions in terms of the forward Compton scattering as a lower part of the TPE graph; see Appendix D for details.
In the following, we study the difference in the individual channel contribution to HFS correction between the traditional HFS expressions and the DR approach based on the unitarity and analyticity of the forward lepton-proton scattering amplitudes. Such a difference between the direct box graph evaluation and unsubtracted dispersion relation result was studied in Refs. [52,53,63] in the case of the TPE correction to the non-forward elastic lepton-proton scattering.
Traditionally, the TPE correction to HFS is expressed by a sum of the proton state contribution el 0 and the polarizability correction pol 0 : The proton intermediate state TPE correction to HFS el 0 is expressed as a sum of the Zemach correction Z with subtraction of the TPE contribution, which is already accounted for in the hydrogen wave functions, and the recoil correction p R : with β 1 (τ ) = −3τ + 2τ 2 + 2(2 − τ ) √ τ (1 + τ ) [42]. The substitution of the photon-proton-proton vertex in the TPE box graph by the on-shell vertex with the Dirac and Pauli couplings [4,35,36], which reproduces the proton state contribution to the VVCS tensor, with a subsequent loop integration, leads to the result for p 0 : which differs from el 0 by the non-pole contribution F 2 P : This correction comes from the proton non-pole term in the spin-dependent forward double virtual Compton scattering tensor, see Appendix D for a description of the VVCS tensor, and should not contribute to the resulting TPE correction as the non-pole terms are not reproduced by the dispersion relations [43]. We express the elastic TPE contribution to the S-level HFS el HFS in the lepton-proton amplitudes DR framework with subtraction of the TPE contribution, which is already accounted for in the hydrogen wave functions, in terms of the proton electric and magnetic form factors as We reproduce the Zemach correction [34] as the nonrelativistic limit of the elastic TPE correction el HFS . In our approach, the non-pole term of Eq. (54) does not appear. The difference between the proton state contribution in the box graph model with an on-shell vertex [4,35,36]  Traditionally, the polarizability correction pol 0 is given by [42,43] The non-pole term F 2 P allows one to expand the HFS integrand near Q 2 = 0 in terms of the proton spin polarizabilities [45] and to be consistent with the massless lepton limit of this correction [43]. In order to avoid the non-pole contributions in the resulting TPE correction to HFS, it was subtracted from the recoil correction in Refs. [42,43]; see Eq. (52).
We express the inelastic α 5 -correction to the S-level HFS in the lepton-proton amplitudes DR approach in terms of the proton inelastic spin SFs g 1 and g 2 as The remaining difference in the inelastic correction pol 0 − F 2 P − inel HFS is given by the contribution from the spin SF g 2 to the amplitude g 2γ (m) + f 2γ − (m). Consequently, the resulting TPE correction to HFS of S energy levels can be equivalently expressed as

HFS correction in both approaches
For the numerical evaluation of the TPE corrections to HFS from the proton intermediate state and the F 2 P part of the polarizability correction pol 0 we exploit the elastic form factor parametrizations from Refs. [12,13]. For the Zemach correction we make two evaluations for the 1 − σ band curves coming from the elastic proton form factor uncertainties of Refs. [12,13], where a global analysis of the electronproton scattering data with account of TPE corrections for Q 2 < 10 GeV 2 was performed, and estimate the uncertainty as a half of a difference between these two curves. For the numerical evaluation of the inelastic correction we exploit the spin SFs data parametrization from Refs. [64][65][66] in the region of large Q 2 . In the region of low Q 2 , we expand the Q 2 -integrand from the proton spin SFs in terms of small x Bj and account for the leading non-vanishing moments: with the moments of the proton spin SFs: with x inel thr = Q 2 /(2Mν inel thr ) and the low-energy constants values [67][68][69][70][71][72]: In Fig. 3 we show the integrand I HFS (Q) entering the TPE correction: in the case of eH and μH. The low-Q behavior based on the moments of the proton spin SFs of Eqs. (61), (62)-(64) I LE HFS and the high-Q behavior based on the data I d HFS are almost independent of the way to evaluate the TPE correction (for Q > 0.5 GeV the two methods agree within 2.5%). While in the region 0.2 GeV Q 0.5 GeV the HFS evaluation with the DRs for the lepton-proton amplitudes ( HFS = el HFS + inel HFS ) and the traditional HFS evaluation ( HFS = 0 = el 0 + pol 0 ) slightly differ. New data in this kinematical region will be very useful for such evaluation. In order to avoid any model dependence, we connect the two with f (Q) given by Furthermore, Q LE = 0.2 GeV, Q d = 0.5 GeV, a 0 = 0.1 GeV, and the constants c 1 , c 2 , Q LE , a 0 were chosen as those that preserve the regularity and smoothness of the integrand I HFS (Q).
In the low-Q region, we make two evaluations for the 1−σ bands of the elastic proton form factors from Refs. [12,13]. We add the combined uncertainty from γ 0 , δ LT , I 1 (0) lin- For the larger Q > 0.5 GeV region, we add the uncertainty of the proton spin structure function parametrization in quadrature to the uncertainties coming from the proton elastic form factors. We connect the high-Q integrands I d by two curves to the 1 − σ boundaries in the low-Q region I LE . We estimate the uncertainty from the difference between the integral of Eq. (68) for these two curves, which are shown in Fig. 4 for eH (μH), and take the averaged central value. In the region Q 2 > 10 GeV 2 , the sizable contribution comes only from the μ P term in Eqs. (50)-(55) and does not introduce any sizable additional uncertainty. We evaluate the proton TPE correction to HFS in the DR approach either for the forward double virtual Compton scattering (TPE correction of Refs. [42,43]) or for the lepton-  (29) proton amplitudes. Also we present results for the Zemach correction, the recoil correction and the F 2 P contribution in Table 1. The evaluation of the resulting TPE correction is performed for the sum of elastic and inelastic contributions using the traditional expressions and the expressions based on the DRs for lp amplitudes. The latter leads to twice smaller uncertainties. However, the two evaluations agree within errors, which is a good test of the proton spin structure function g 2 parametrization for such calculation. The leading Zemach correction is a bit smaller than the result based on the typical proton form factors parametrization in Ref. [42] due to the suppressed low-Q 2 behavior of the magnetic proton form factor measured by the A1 Collaboration at MAMI [12,13].

Conclusions and outlook
In this work we have derived dispersion relations for the forward lepton-proton scattering TPE amplitudes. We have fully expressed the spin-dependent TPE amplitudes in terms of the inclusive total cross sections at the leading α order or equivalently in terms of the proton spin SFs. With these relations we have found a new way to determine the O(α 5 ) TPE proton structure correction to the S-level HFS in the hydrogen-like atoms explicitly accounting for the helicity double spin-flip amplitude. The result for individual channel HFS contribution is distinct to the standard approach, e.g., the elastic TPE contribution to HFS differs by the finite-sizable correction. Only accounting for contributions from all possible channels through the Burkard-Cottingham sum rule the two methods agree. We have reevaluated the TPE correction to HFS in electronic and muonic hydrogen connecting the region with small photons virtualities, which is expressed through the first and third moments of the proton spin SFs, to the region with large photons virtualities, where the data (parametrizations) on the proton form factors and SFs exists. The resulting TPE correction is similar to the literature result with slightly larger uncertainties in the standard approach to HFS and slightly smaller uncertainties within our approach. In a subsequent paper, we plan to decrease the errors of the TPE correction to HFS in the muonic hydrogen.
We introduce the finite photon mass μ with the aim to have no deal with IR divergences. Such process is completely The high-energy behavior of the relevant cross sections in the OPE approximation is as follows: The unsubtracted DR for the amplitude f 2γ + of Eq. (15) is divergent as ω, therefore we use the subtracted DR of Eq. (18) for this amplitude.
The helicity amplitude corresponding with the TPE direct box graph T 2γ h λ hλ is given by with γ.a ≡ γ μ a μ . We find the contribution of the direct box graph to the forward amplitudes f 2γ,dir + , f 2γ,dir − with exchange of two photons (see Fig. 1) multiplying the fermion spinors by the spin projection operators. Then we sum over all possible polarizations evaluating traces of the Dirac matrices.
The direct amplitudes f with the lepton and proton spin vectors in the laboratory frame s = (| k|, 0, 0, ω)/m and S = (0, 0, 0, −1). For the double spin-flip amplitude g 2γ,dir we use the decomposition in terms of the scalar integrals, since the evaluation of traces cannot be exploited here due to the different spin directions of the initial and final fermions. Furthermore, we repeat the same steps for the crossed box graph contribution. The optical theorem of Eqs. (12) and (13), the oncesubtracted DR for f 2γ + amplitude of Eq. (18), the unsubtracted DRs for f 2γ − , g 2γ amplitudes of Eqs. (16,17) and the amplitudes properties under the crossing ω → −ω of Eqs. (9)-(11) were checked comparing with the sum of the direct and crossed box graphs.
The reasonable result for the nucleon or narrow contribution, when we are allowed to interchange the ω and Q 2 , ν γ integration order, is given by the once-subtracted dispersion relation of Eq. (18). It can be formally obtained from the subtracted at the point ω s Eq. (C1) by f 2γ The ν γ integrals for the spin-dependent amplitudes of Eqs. (33) and (32) are convergent. Due to the Regge behavior of the F 1 proton structure function given by the Pomeron exchange, the ν γ integrals are divergent, and the DR for the amplitude f 2γ + is not applicable for the inelastic intermediate states TPE contribution.

Appendix D: HFS through the forward double virtual Compton scattering amplitudes
It is instructive to evaluate the HFS correction cutting only the lower blob of the TPE graph [42,43], see right panel in Fig. 6.