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

TL;DR

This paper develops a new method to evaluate the two-photon exchange correction to hyperfine splitting in hydrogen-like atoms by relating TPE amplitudes to inclusive scattering cross sections, offering an alternative to traditional approaches.

Contribution

It introduces a novel approach linking TPE amplitudes to scattering cross sections, providing a different evaluation method for hyperfine splitting corrections.

Findings

The new method reproduces the standard sum rule results.

Numerical comparisons show differences in elastic contributions.

The approach offers a distinct evaluation of TPE effects in hydrogen.

Abstract

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.

Figures

Tables

Table 1. Table 1: Finite-size TPE correction to the hyperfine splitting of S energy levels in e 𝑒 e H and μ 𝜇 \mu H. The Fermi energy HFS is 5.86785 ​ μ ​ eV 5.86785 𝜇 eV 5.86785~{}\mathrm{\mu eV} for the 1S level in e 𝑒 e H and 182.4432 ​ meV 182.4432 meV 182.4432~{}\mathrm{meV} for the 1S level in μ 𝜇 \mu H.

	${10}^6\Delta$, $eH$	${10}^3\Delta$, $\mu H$
Zemach, ${\Delta }_Z$	$-$39.59(75)	$-$7.36(14)
Non-pole term, ${\Delta }^{F_P^2}$	22.53(7)	1.11(1)
Recoil, ${\Delta }_R^p$	5.31(7)	0.8476(84)
Traditional elastic, ${\Delta }_0^{\mathrm{el}}$	$-$34.29(75)	$-$6.51(14)
Total elastic, ${\Delta }_{\mathrm{HFS}}^{\mathrm{el}}$	$-$43.40(74)	$-$7.03(14)
Total ${\Delta }_{\mathrm{HFS}}$, within VVCS DRs	$-$32.60(2.56)	$-$6.18(49)
Total ${\Delta }_{\mathrm{HFS}}$ within $lp$ DRs	$-$32.80(1.56)	$-$6.22(29)

Equations

$$T_{++++},T_{+-+-},T_{--++},$$

$$f_{\pm }\left(\omega \right)$$

$$g\left(\omega \right)$$

$$T_{h^{\prime }{\lambda }^{\prime },h\lambda }\left(\omega \right)$$

$$\bar{u}(k,h^{\prime })u(k,h)=2m{\delta }_{h^{\prime },h},\bar{N}(p,{\lambda }^{\prime })N(p,\lambda )=2M{\delta }_{{\lambda }^{\prime },\lambda }.$$

$$f_{+}^{l^{+}p}\left(\omega \right)$$

$$f_{-}^{l^{+}p}\left(\omega \right)$$

$$g^{l^{+}p}\left(\omega \right)$$

$$f_{-}^{l^{\pm }p}=\sum _{n=1}^{\infty }\left(f_{-}^{(2n)\gamma }\pm f_{-}^{(2n-1)\gamma }\right),$$

$$f_{+}^{n\gamma }\left(\omega \right)$$

$$f_{-}^{n\gamma }\left(\omega \right)$$

$$g^{n\gamma }\left(\omega \right)$$

$$\Im f_{\pm }\left(\omega \right)$$

$$\Im g\left(\omega \right)$$

$$\displaystyle\Re\left\{\begin{array}[]{l}f_{\pm}\left(\omega\right)\\ g(\omega)\end{array}\right\}$$

$$\Re f_{+}^{2\gamma }\left(\omega \right)$$

$$\Re f_{-}^{2\gamma }\left(\omega \right)$$

$$\Re g^{2\gamma }\left(\omega \right)$$

$$\displaystyle\Re f^{2\gamma}_{+}\left(\omega\right)-\Re f^{2\gamma}_{+}\left(\omega_{s}\right)=\frac{4M\left(\omega^{2}-\omega_{s}^{2}\right)}{\pi}\fint\limits^{\infty}_{m}\frac{\omega^{\prime}|\vec{k}^{\prime}|\sigma\left(\omega^{\prime}\right)}{\left(\omega^{\prime 2}-\omega^{2}\right)\left(\omega^{\prime 2}-\omega_{s}^{2}\right)}\mathrm{d}\omega^{\prime}.$$

$$f_{-}^{1\gamma }$$

$$g^{1\gamma }$$

$$d\sigma$$

$$x_{\mathrm{Bj}}$$

$$L^{\mu \nu }$$

$$W_{\mu \nu }$$

$$\frac{{\mathrm{d}}^2\sigma }{\mathrm{d}{\nu }_{\gamma }d{Q}^2}$$

$$\frac{{\mathrm{d}}^2{\sigma }_{++}-{\mathrm{d}}^2{\sigma }_{+-}}{\mathrm{d}{\nu }_{\gamma }\mathrm{d}{Q}^2}$$

$$\frac{{\mathrm{d}}^2{\sigma }_{\perp }-{\mathrm{d}}^2{\sigma }_{\parallel }}{\mathrm{d}{\nu }_{\gamma }\mathrm{d}{Q}^2}$$

$$F_1\left(x_{\mathrm{Bj}},Q^2\right)$$

$$F_2\left(x_{\mathrm{Bj}},Q^2\right)$$

$$g_1\left(x_{\mathrm{Bj}},Q^2\right)$$

$$g_2\left(x_{\mathrm{Bj}},Q^2\right)$$

$$\Re g^{2\gamma }\left(\omega \right)$$

$$\Re f_{-}^{2\gamma }\left(\omega \right)$$

$$|\vec{k}|$$

$${\omega }_0$$

$${\tau }_l$$

$$f_{-}^{2\gamma }\left(m\right)$$

Full text

1. I Introduction
2. II Forward lepton-proton scattering
3. III Relation of the forward two-photon amplitudes to the proton structure
4. IV Lamb shift and hyperfine splitting (HFS)
5. V HFS correction in both approaches
6. VI Conclusions and outlook
7. VII Acknowledgments
8. A Crossing in lepton-proton scattering
9. B Dispersion relations verification in QED
10. C Forward invariant amplitudes in terms of the proton structure functions
11. D HFS through the forward double virtual Compton scattering amplitudes
12. E Forward scattering observables

I 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 Rosenbluth:1950yq and the polarization transfer Jones:1999rz methods, cf. Guichon:2003qm ; Blunden:2003sp and Carlson:2007sp ; Arrington:2011dn 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 electron-proton to positron-proton scattering cross sections, where TPE enters with different signs, have recently been carried out at VEPP-3 (Novosibirsk) Rachek:2014fam , OLYMPUS (DESY) Henderson:2016dea , and CLAS (JLab) Adikaram:2014ykv . 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 ($\mu$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 Pohl:2010zza ; Antognini:1900ns and measurements with electrons Bernauer:2010wm ; Bernauer:2013tpr ; Mohr:2012tt , see Antognini:1900ns ; Carlson:2015jba 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 Antognini:1900ns . Also a high-precision measurement of 1S HFS in muonic hydrogen with an unprecendented ppm accuracy level is planned by this Pohl:2016tqq ; Dupays:2003zz and other collaborations Ma:2016etb ; Adamczak:2016pdb . 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 scattering (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, $\nu_{\gamma}$, and virtuality, $Q^{2}>0$. For the Lamb shift correction evaluation Pachucki:1996zza ; Faustov:1999ga ; Pineda:2002as ; Pineda:2004mx ; Nevado:2007dd ; Carlson:2011zd the subtraction function is needed in addition. However, the latter can be estimated at low virtualities from the chiral perturbation theory Birse:2012eb ; Alarcon:2013cba ; Peset:2014jxa or non-relativistic quantum electrodynamics (NRQED) Hill:2012rh and at high virtualities from the operator product expansion Hill:2016bjv . This function, in principle, can also be determined with account of the high-energy SFs data Gorchtein:2013yga ; Tomalak:2015hva ; Caprini:2016wvy . 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 Zemach:1956zz ; Iddings:1959zz ; Iddings:1965zz ; Drell:1966kk ; Faustov:1966 ; Faustov:1970 ; Bodwin:1987mj ; Faustov:2001pn ; Carlson:2008ke ; Carlson:2011af . For recent numerical evaluations of the TPE correction to HFS see Refs. Faustov:2001pn ; Carlson:2008ke ; Carlson:2011af , for the model-independent evaluation within the frameworks of NRQED and ChPT exploiting the hydrogen HFS measurement see Ref. Peset:2016wjq and for results in chiral EFT see Ref. Hagelstein:2015egb .

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 spin-flip 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 GellMann:1954db ; Pascalutsa:2010sj . 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 Grein:1977mn ; Kroll:1973re . 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 Grein:1977mn ; Kroll:1973re ; Guichon:1982eb . Advantageously, the method based on DRs for $lp$ amplitudes Gorchtein:2006mq ; Borisyuk:2008es ; Tomalak_PhD 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. Iddings:1959zz ; Iddings:1965zz ; Drell:1966kk ; Faustov:1966 ; Faustov:1970 ; Bodwin:1987mj only after account for the Burkhardt-Cottingham (BC) sum rule Burkhardt:1970ti . 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 $\mathrm{H}$ and $\mu\mathrm{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 $\alpha$ order in Sec. II and verify them in QED in App. B. We relate the inclusive $lp$ cross sections and the forward TPE amplitudes to proton SFs in Sec. III. 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 Sec. IV. The numerical comparison of the proton and inelastic intermediate states TPE contributions to HFS within two methods is given in Sec. V. We give our conclusions in Sec. VI. We also provide the derivation of the TPE correction to HFS using the forward Compton scattering tensor in App. D. In App. E we classify all possible elastic scattering polarization observables in the forward kinematics.

II Forward lepton-proton scattering

The forward elastic $lp$ scattering is described by three Gribov:1963gx non-vanishing independent helicity amplitudes $T_{h^{\prime}\lambda^{\prime}h\lambda}$, with $h$ ($h^{\prime}$) and $\lambda$ ($\lambda^{\prime}$) the helicities of the inital (final) lepton and proton, see Fig. 1 for the notations of kinematics and helicities:

[TABLE]

which are functions of the lepton energy $\omega$ 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 even-odd properties with respect to the crossing $\omega\to-\omega$:

[TABLE]

The Lorentz structure of the forward amplitude is then given by

[TABLE]

with $m$ and $u$ the lepton mass and spinor, $\vec{k}$ the lepton momentum in the lab frame, $M$ and $N$ the proton mass and spinor, $\gamma^{\mu\nu}=\frac{1}{2}[\gamma^{\mu},\gamma^{\nu}]$, where $\gamma^{\mu}$ are the Dirac matrices, the spinors are normalized as

[TABLE]

In order to establish the even-odd properties for the invariant amplitudes under $\omega\to-\omega$, we first perform the crossing on the lepton line and relate amplitudes of the lepton-proton scattering $f^{l^{-}p}(\omega)$ in the physical region ($\omega>0$) to amplitudes of the antilepton-proton scattering $f^{l^{+}p}(-\omega)$ in the unphysical region ($\omega<0$):

[TABLE]

where $\omega$ 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)\gamma}\left(\omega\right)$ and odd $f^{(2n-1)\gamma}\left(\omega\right)$ number of photons connected to both lepton and proton lines, e.g.:

[TABLE]

and obtain the following crossing relations for the contribution of graphs with $n$ exchanged photons on the real $\omega$ axis:

[TABLE]

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:

[TABLE]

where $\sigma_{h\lambda}$ is the inclusive cross section with the incoming lepton helicity $h$ and the incoming proton helicity $\lambda$; $\sigma_{\perp}$ ($\sigma_{\parallel}$) is the inclusive cross section with lepton and proton polarized transversely and perpendicular (parallel) to each other. 111The TPE correction from the inelastic intermediate states was related to the unpolarized photoabsorption cross section in Ref. Brown:1970te .

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 Section IV.

We express the latter cross sections in terms of the proton SFs up to the order $\alpha^{2}$ in Section III and obtain the imaginary parts of the TPE amplitudes at the leading $\alpha$ order with Eqs. (12-13).

Assuming analyticity of these amplitudes in the entire complex $\omega$-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:

[TABLE]

where $\fint$ stands for the principal-value integration. Using properties of the TPE amplitudes under the crossing $\omega\to-\omega$, see Eqs. (9-11), and assuming the sufficient high-energy 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 $\alpha$ order and thereby account for both direct and crossed graphs:

[TABLE]

with the lepton momentum in the lab frame $|\vec{k}^{\prime}|=\sqrt{\omega^{\prime 2}-m^{2}}$. These DRs are written for amplitudes in one channel contrary to the DRs for the forward proton-proton scattering Grein:1977mn .

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. (19). One should perform a subtraction in the DR for the amplitude $f_{+}^{2\gamma}$, e.g., at the point $\omega_{s}$: 222In 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. Miller:2012ne .

[TABLE]

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 $\sigma$ divergent due to the virtual photons with high energy in the lab frame. Consequently the DRs of Eqs. (19) and (22) cannot be written for the inelastic TPE contribution to the unpolarized amplitude $f_{+}$.

The DRs of Eqs. (20-22) have the same form as the DRs for the light-by-light scattering Pascalutsa:2010sj .

We decompose the forward scattering amplitudes into a sum of one-photon exchange (OPE) and TPE contributions. The TPE amplitudes, except for $\Re f^{2\gamma}_{+}$, are obtained with DRs and unitarity relations as described above. The OPE amplitudes are real and given by

[TABLE]

with the proton magnetic moment $\mu_{P}\approx 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}$: $\delta 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 Maximon:2000hm . 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, 20, 21), at $O(\alpha^{2})$ up to the one unknown spin-independent amplitude $\Re f_{+}$.

III 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 $\nu_{\gamma}$ 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^{\mu\nu}$ and the hadronic tensor $W^{\mu\nu}$ Agashe:2014kda . It is proportional to the phase space of the final lepton (with 4-momentum $k^{\prime}=\left(\omega^{\prime},\vec{k}^{\prime}\right)$ in the lab frame) and given by

[TABLE]

with the unit of electric charge $e$. The kinematics are traditionally described by the kinematical Bjorken variable $x_{\mathrm{Bj}}$, the variable $y$ related to the energy transferred by the virtual photon relative to the beam energy and the momentum transfer $Q^{2}$:

[TABLE]

The leptonic tensor is evaluated in QED. It is given by

[TABLE]

where $s^{\mu}$ is the lepton spin vector: $s^{\mu}s_{\mu}=-1,~{}(s\cdot k)=0$. The general Lorentz and gauge-invariant structure of the hadronic tensor $W^{\mu\nu}$ which preserves parity and charge conjugation invariance is given by

[TABLE]

with the virtual photon energy in the laboratory frame $\nu_{\gamma}=\left(p\cdot q\right)/M$ and the proton SFs $F_{1}\left(\nu_{\gamma},Q^{2}\right),~{}F_{2}\left(\nu_{\gamma},Q^{2}\right),~{}g_{1}\left(\nu_{\gamma},Q^{2}\right),~{}g_{2}\left(\nu_{\gamma},Q^{2}\right)$, which are extracted from the experimental data. The proton spin 4-vector satisfies: $S^{2}=-1,~{}(S\cdot p)=0$.

The total unpolarized cross section is given by

[TABLE]

This expression reduces to the well-known expression Agashe:2014kda in the massless limit.

Consider a scattering of longitudinally polarized leptons on the proton polarized in the lepton momentum direction $\sigma_{h\lambda}=\sigma_{{+}{-}}$ and the scattering on the proton polarized in the opposite direction $\sigma_{h\lambda}=\sigma_{{+}{+}}$ with the proton (lepton) spin vector in the laboratory frame $S^{\mu}=(0,-\lambda\hat{\vec{k}})$ (s^{\mu}=(|\vec{k}|,\omega\hat{\vec{k}})$$/m) and $\hat{\vec{k}}=\vec{k}/|\vec{k}|$. For the cross sections difference we obtain:

[TABLE]

This expression reduces to the known expression Agashe:2014kda , Blumlein:1996vs in the massless limit.

Consider the scattering of transversely polarized leptons on transversely polarized protons. Denoting the averaged over the azimuthal angle cross section $\sigma_{\perp}$ ($\sigma_{\parallel}$) for scattering with perpendicular (parallel) spin vectors of lepton ($s^{\mu}=(0,\cos\phi_{l},\sin\phi_{l},0)$) and proton ($S^{\mu}=(0,\cos\phi_{p},\sin\phi_{p},0)$), 333The nontrivial relation $\sigma_{\perp}=\sigma$ holds for the lepton-proton scattering. We have obtained this relation in the OPE approximation of Eq. (25) and have proved it by the direct cross sections evaluation for the case of elastic scattering and exploiting the symmetry properties of the helicity amplitudes for the scattering to arbitrary channel $lp\to lX$. i. e. $\phi_{l}-\phi_{p}=\pm\pi/2$ ($\phi_{l}-\phi_{p}=0$) for the perpendicular (parallel) configuration, we obtain:

[TABLE]

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\to l^{-}p$ are obtained by substitution of the inelastic SFs in Eqs. (III-31) by the elastic contribution to them:

[TABLE]

where $F_{D}(Q^{2})$, $F_{P}(Q^{2})$, $G_{E}(Q^{2})$, $G_{M}(Q^{2})$ are the Dirac, Pauli, Sachs electric and magnetic proton form factors and $\tau_{P}=Q^{2}/(4M^{2})$, following by the integration over the virtual photon energy.

Substituting the expressions for the inclusive cross sections of Eqs. (III-31) into the DRs, see Eqs. (20-21), 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:

[TABLE]

where we have introduced the notations:

[TABLE]

In Eqs. (37, 36) the elastic threshold $\nu_{\mathrm{thr}}$ and the inelastic threshold $\nu^{\mathrm{inel}}_{\mathrm{thr}}$ are given by $\nu_{\mathrm{thr}}=0$ and $\nu^{\mathrm{inel}}_{\mathrm{thr}}=m_{\pi}+\left(m_{\pi}^{2}+Q^{2}\right)/\left(2M\right)$ respectively, where $m_{\pi}$ denotes the pion mass.

The leading TPE correction to the atomic energy levels is given by the values of the amplitudes at threshold $\omega=m$. The TPE amplitudes $f_{-}^{2\gamma},~{}g^{2\gamma}$ at threshold can then be expressed in terms of the proton spin SFs as

[TABLE]

with

[TABLE]

Evaluating the sum of the spin-dependent lepton-proton TPE amplitudes in the DR approach, see Eqs. (LABEL:fminus) and (42), we obtain:

[TABLE]

which is a trivial relation due to the Burkhardt-Cottingham (BC) sum rule Burkhardt:1970ti :

[TABLE]

IV 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\gamma}$ is expressed in terms of the forward amplitudes $f^{2\gamma}_{+},~{}f^{2\gamma}_{-},~{}g^{2\gamma}$ as

[TABLE]

with $\boldsymbol{s}$ ($\boldsymbol{S}$) and $\hat{\vec{k}}$ ($\hat{\vec{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. Guichon:1982eb .

It is then easy to see that, considering $T^{2\gamma}$ as correction to the Coulomb potential, its effect on the nS-state energy level is given by

[TABLE]

with $|\psi_{\mathrm{nS}}(0)|^{2}=\alpha^{3}m_{r}^{3}/(\pi 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\gamma}$ 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:

[TABLE]

This correction yields the third Zemach moment Hagelstein:2015egb .

The TPE contribution to the nS-level HFS $\delta E^{\mathrm{HFS}}_{\mathrm{nS}}$ is expressed in terms of the relative correction $\Delta_{\mathrm{HFS}}$ and the leading order nS-level HFS $E^{\mathrm{HFS},0}_{\mathrm{nS}}$ (Fermi energy) as 444Note the absence of the muon anomalous magnetic moment in Eq. (50) Eides:2000xc .

[TABLE]

Considering the spin part of $T^{2\gamma}$ 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\gamma},~{}g^{2\gamma}$ at threshold ($\omega=m$):

[TABLE]

The TPE correction to the S-level HFS $\Delta_{\mathrm{0}}$ of Refs. Iddings:1959zz ; Iddings:1965zz ; Drell:1966kk ; Faustov:1966 ; Faustov:1970 ; Bodwin:1987mj ; Faustov:2001pn ; Carlson:2008ke ; Carlson:2011af ; Peset:2016wjq , which was derived from the difference between the spin-averaged forward matrix elements in the singlet and triplet states with a subsequent loop integration, can be also obtained adding $g^{2\gamma}(m)+f_{-}^{2\gamma}(m)=0$ to Eq. (51):

[TABLE]

Consequently, we have verified the TPE correction to HFS of S energy levels: $\Delta_{\mathrm{0}}=\Delta_{\mathrm{HFS}}$. The HFS correction of Refs. Iddings:1959zz ; Iddings:1965zz ; Drell:1966kk ; Faustov:1966 ; Faustov:1970 ; Bodwin:1987mj ; Faustov:2001pn ; Carlson:2008ke ; Carlson:2011af can be reproduced performing the loop integration for the amplitudes $g^{2\gamma}(m)$ and $f_{-}^{2\gamma}(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. Borisyuk:2008es ; Tomalak_PhD ; Blunden:2017nby 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 $\Delta^{\mathrm{el}}_{0}$ and the polarizability correction $\Delta^{\mathrm{pol}}_{\mathrm{0}}$:

[TABLE]

The proton intermediate state TPE correction to HFS $\Delta^{\mathrm{el}}_{0}$ is expressed as a sum of the Zemach correction $\Delta_{\mathrm{Z}}$ with subtraction of the TPE contribution, which is already accounted for in the hydrogen wave functions, and the recoil correction $\Delta^{\mathrm{p}}_{\mathrm{R}}$:

[TABLE]

with $\beta_{1}(\tau)=-3\tau+2\tau^{2}+2(2-\tau)\sqrt{\tau(1+\tau)}$ Carlson:2008ke . The substitution of the photon-proton-proton vertex in the TPE box graph by the on-shell vertex with the Dirac and Pauli couplings Iddings:1959zz ; Iddings:1965zz ; Blunden:2003sp , which reproduces the proton state contribution to the VVCS tensor, with a subsequent loop integration, leads to the correction $\Delta^{\mathrm{el}}_{0}$, while the result based on the dispersion relations for the VVCS tensor $\Delta^{\mathrm{p}}_{0}$ is given by

[TABLE]

which differs from $\Delta^{\mathrm{el}}_{0}$ by the non-pole contribution $\Delta^{\mathrm{F^{2}_{P}}}$:

[TABLE]

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 Carlson:2011af .

We express the elastic TPE contribution to the S-level HFS $\Delta^{\mathrm{el}}_{\mathrm{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

[TABLE]

We reproduce the Zemach correction Zemach:1956zz as the non-relativistic limit of the elastic TPE correction $\Delta^{\mathrm{el}}_{\mathrm{HFS}}$. In our approach, the non-pole term of Eq. (58) does not appear. The difference between the proton state contribution in the box graph model with an on-shell vertex Iddings:1959zz ; Iddings:1965zz ; Blunden:2003sp $\Delta^{\mathrm{p}}_{0}$ and within dispersion relation approach $\Delta^{\mathrm{el}}_{\mathrm{HFS}}$, see Eqs. (54-58) and Eq. (59), $\Delta^{\mathrm{p}}_{0}-\Delta^{\mathrm{el}}_{\mathrm{HFS}}$ is given by the elastic contribution to the amplitude $g^{2\gamma}(m)+f_{-}^{2\gamma}(m)$.

Traditionally, the polarizability correction $\Delta^{\mathrm{pol}}_{\mathrm{0}}$ is given by Carlson:2008ke ; Carlson:2011af

[TABLE]

The non-pole term $\Delta^{\mathrm{F^{2}_{P}}}$ allows one to expand the HFS integrand near $Q^{2}=0$ in terms of the proton spin polarizabilities Hagelstein:2015egb and to be consistent with the massless lepton limit of this correction Carlson:2011af . In order to avoid the non-pole contributions in the resulting TPE correction to HFS, it was subtracted from the recoil correction in Refs. Carlson:2008ke ; Carlson:2011af , see Eq. (56).

We express the inelastic $\alpha^{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

[TABLE]

The remaining difference in the inelastic correction $\Delta^{\mathrm{pol}}_{\mathrm{0}}-\Delta^{\mathrm{F^{2}_{P}}}-\Delta^{\mathrm{inel}}_{\mathrm{HFS}}$ is given by the contribution from the spin SF $g_{2}$ to the amplitude $g^{2\gamma}(m)+f_{-}^{2\gamma}(m)$.

Consequently, the resulting TPE correction to HFS of S energy levels can be equivalently expressed as

[TABLE]

V HFS correction in both approaches

For the numerical evaluation of the TPE corrections to HFS from the proton intermediate state and the $\Delta^{\mathrm{F^{2}_{P}}}$ part of the polarizability correction $\Delta^{\mathrm{pol}}_{\mathrm{0}}$ we exploit the elastic form factor parametrizations from Refs. Bernauer:2010wm ; Bernauer:2013tpr . For the Zemach correction we make two evaluations for the 1-$\sigma$ band curves coming from the elastic proton form factor uncertainties of Refs. Bernauer:2010wm ; Bernauer:2013tpr , where a global analysis of the electron-proton scattering data with account of TPE corrections for $Q^{2}<10~{}\mathrm{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. Kuhn:2008sy ; griffioen ; Sato:2016tuz 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_{\mathrm{Bj}}$ and account for the leading non-vanishing moments:

[TABLE]

with the moments of the proton spin SFs:

[TABLE]

with $x^{\mathrm{inel}}_{\mathrm{thr}}=Q^{2}/(2M\nu^{\mathrm{inel}}_{\mathrm{thr}})$ and the low-energy constants values Drechsel:2000ct ; Drechsel:2007if ; Drechsel:2002ar ; Prok:2008ev ; Pascalutsa:2014zna ; Lensky:2017dlc :

[TABLE]

In Fig. 3 we show the integrand $I_{\mathrm{HFS}}(Q)$ entering the TPE correction:

[TABLE]

in the case of $e\mathrm{H}$ and $\mu\mathrm{H}$. The low-$Q$ behavior based on the moments of the proton spin SFs of Eqs. (65, 66-68) $I^{\mathrm{LE}}_{\mathrm{HFS}}$ and the high-$Q$ behavior based on the data $I^{\mathrm{d}}_{\mathrm{HFS}}$ are almost independent of the way to evaluate the TPE correction (for $Q>0.5~{}\mathrm{GeV}$ the two methods agree within $2.5\%$). While in the region $0.2~{}\mathrm{GeV}\lesssim Q\lesssim 0.5~{}\mathrm{GeV}$ the HFS evaluation with the DRs for the lepton-proton amplitudes ($\Delta_{\mathrm{HFS}}=\Delta_{\mathrm{HFS}}^{\mathrm{el}}+\Delta_{\mathrm{HFS}}^{\mathrm{inel}}$) and the traditional HFS evaluation ($\Delta_{\mathrm{HFS}}=\Delta_{0}=\Delta^{\mathrm{el}}_{0}+\Delta_{0}^{\mathrm{pol}}$) 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 model-independent regions by the function of the Fermi-Dirac distribution type:

[TABLE]

with $f(Q)$ given by

[TABLE]

Furthermore, $Q_{\mathrm{LE}}=0.2~{}\mathrm{GeV}$, $Q_{\mathrm{d}}=0.5~{}\mathrm{GeV}$, $a_{0}=0.1~{}\mathrm{GeV}$, and the constants $c_{1},~{}c_{2},~{}Q_{\mathrm{LE}},~{}a_{0}$ were chosen as those that preserve the regularity and smoothness of the integrand $I_{\mathrm{HFS}}(Q)$.

In the low-$Q$ region, we make two evaluations for the 1-$\sigma$ bands of the elastic proton form factors from Refs. Bernauer:2010wm ; Bernauer:2013tpr . We add the combined uncertainty from $\gamma_{0},~{}\delta_{LT},~{}{I^{\prime}}_{1}\left(0\right)$ linearly to the form factors uncertainty. For the larger $Q>(0.013-0.017)~{}\mathrm{GeV}$, we make evaluation for the central values of the proton elastic form factors and add the uncertainty of the proton elastic form factors to the uncertainties from $\gamma_{0},~{}\delta_{LT},~{}{I^{\prime}}_{1}\left(0\right)$ in quadrature. The boundary $Q$ value is chosen as the value that leads to the same uncertainties in the proton intermediate state HFS contribution in both ways of the error estimate described in this paragraph. For the larger $Q>0.5~{}\mathrm{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^{\mathrm{d}}$ by two curves to the 1-$\sigma$ boundaries in the low-$Q$ region $I^{\mathrm{LE}}$. We estimate the uncertainty from the difference between the integral of Eq. (72) for these two curves, which are shown in Fig. 4 for $e\mathrm{H}$ ($\mu\mathrm{H}$), and take the averaged central value. In the region $Q^{2}>10~{}\mathrm{GeV}^{2}$, the sizable contribution comes only from the $\mu_{P}$ term in Eqs. (54-59) 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. Carlson:2008ke ; Carlson:2011af ) or for the lepton-proton amplitudes. Also we present results for the Zemach correction, the recoil correction and the $\Delta^{\mathrm{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. Carlson:2008ke due to the enhanced low-$Q^{2}$ behavior of the magnetic proton form factor measured by the A1 Collaboration at MAMI Bernauer:2010wm ; Bernauer:2013tpr .

VI 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 $\alpha$ order or equivalently in terms of the proton spin SFs. With these relations we have found a new way to determine the $O(\alpha^{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 based on the unitarity and analyticity of the lepton-proton scattering amplitudes 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.

VII Acknowledgments

We thank Vladimir Pascalutsa for useful discussions about general formalism and advices given during this work and manuscript preparation, Marc Vanderhaeghen for useful discussions and advises concerning the numerical evaluations, Carl Carlson for the discussions concerning literature. We thank Keith Griffioen, Sebastian Kuhn, Nevzat Guler and Jacob Ethier for providing us with results on the proton spin SFs. This work was supported by the Deutsche Forschungsgemeinschaft (DFG) through Collaborative Research Center “The Low-Energy Frontier of the Standard Model” (SFB 1044), and Graduate School “Symmetry Breaking in Fundamental Interactions” (DFG/GRK 1581).

Appendix A Crossing in lepton-proton scattering

In order to establish the even-odd properties for the invariant amplitudes under the crossing $\omega\to-\omega$, we first perform the crossing on the lepton line and relate the amplitudes of the lepton-proton scattering $f^{l^{-}p}(\omega)$ in the physical region ($\omega>0$) to the amplitudes of the antilepton-proton scattering $f^{l^{+}p}(-\omega)$ in the unphysical region ($\omega<0$). We write the general form of the amplitude as

[TABLE]

with $O=\left(1,\,\gamma^{\mu\nu},\,\gamma^{\mu}\gamma_{5}\right)$. We observe that after the replacement in the lepton line $k\to-k$, the amplitude transforms to

[TABLE]

We can rewrite the lepton spinor $u$ in terms of the antilepton spinor $v$ as $u\left(-k,-h\right)=\gamma^{2}v^{*}\left(k,h\right)$, where we exploit the same form $u$ for the antilepton spinor as only particles or antiparticles participate in interaction. The expression for the helicity amplitude is given by

[TABLE]

Transposing the lepton line we obtain:

[TABLE]

The tensor structure of Eq. (3) transforms to

[TABLE]

which corresponds to the scattering of the antilepton off the proton.

According to the crossing properties we can write amplitudes for the scattering of the antilepton $f^{l^{+}p}$ in terms of the lepton scattering amplitudes $f^{l^{-}p}$ as

[TABLE]

Appendix B Dispersion relations verification in QED

In this Appendix, we verify the lepton-proton forward dispersion relations in QED. We reconstruct the real parts of the TPE amplitudes with the relations of Eqs. (20-22) and compare them with the sum of the direct and crossed box graphs. The OPE helicity amplitude $T^{1\gamma}_{h^{\prime}\lambda^{\prime}h\lambda}$ for the lepton scattering off the charged point proton: $l(k,h)+p(p,\lambda)\to l(k^{\prime},h^{\prime})+p(p^{\prime},\lambda^{\prime})$, where $h(h^{\prime})$ denote the incoming (outgoing) lepton helicities and $\lambda(\lambda^{\prime})$ the corresponding proton helicities respectively, see Fig. 5, is given by

[TABLE]

We introduce the finite photon mass $\mu$ with the aim to have no deal with IR divergences. Such process is completely described by 2 Mandelstam variables, e.g., $Q^{2}=-(k-k^{\prime})^{2}$ - the squared momentum transfer, and $s=(p+k)^{2}=M^{2}+2M\omega+m^{2}$ - the squared energy in the lepton-proton c.m. reference frame.

The relevant OPE cross sections are given by

[TABLE]

with $\Sigma_{s}=\left(s-(M+m)^{2}\right)\left(s-(M-m)^{2}\right)=4M^{2}\vec{k}^{2}$. The high-energy behavior of the relevant cross sections in the OPE approximation is as follows:

[TABLE]

The unsubtracted DR for the amplitude $f_{+}^{2\gamma}$ of Eq. (19) is divergent as $\omega$, therefore we use the subtracted DR of Eq. (22) for this amplitude.

The helicity amplitude corresponding with the TPE direct box graph $T^{\mathrm{2\gamma}}_{h^{\prime}\lambda^{\prime}h\lambda}$ is given by

[TABLE]

with $\gamma.a\equiv\gamma^{\mu}a_{\mu}$. We find the contribution of the direct box graph to the forward amplitudes $f_{+}^{2\gamma,\mathrm{dir}},~{}f_{-}^{2\gamma,\mathrm{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_{+}^{2\gamma,\mathrm{dir}},~{}f_{-}^{2\gamma,\mathrm{dir}}$ are given by

[TABLE]

with the lepton and proton spin vectors in the laboratory frame $s=(|\vec{k}|,~{}0,~{}0,~{}\omega)/m$ and $S=(0,~{}0,~{}0,-1)$. For the double spin-flip amplitude $g^{2\gamma,\mathrm{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 once-subtracted DR for $f_{+}^{2\gamma}$ amplitude of Eq. (22), the unsubtracted DRs for $f_{-}^{2\gamma},~{}g^{2\gamma}$ amplitudes of Eqs. (20, 21) and the amplitudes properties under the crossing $\omega\to-\omega$ of Eqs. (9-11) were checked comparing with the sum of the direct and crossed box graphs.

Appendix C Forward invariant amplitudes in terms of the proton structure functions

Substituting expressions for the inclusive cross sections of Eqs. (III-31) into DRs, see Eqs. (19-21), we change the integration order and express the forward TPE amplitudes in terms of the proton SFs:

[TABLE]

with the DR master integrals (we introduce the cut-off $\Lambda$ in the divergent integrals):

[TABLE]

[TABLE]

The reasonable result for the nucleon or narrow $\Delta$ contribution, when we are allowed to interchange the $\omega$ and $Q^{2},~{}\nu_{\gamma}$ integration order, is given by the once-subtracted dispersion relation of Eq. (22). It can be formally obtained from the subtracted at the point $\omega_{s}$ Eq. (91) by $\Re f^{2\gamma}_{+}\left(\omega\right)-\Re f^{2\gamma}_{+}\left(\omega_{s}\right)$:

[TABLE]

The $\nu_{\gamma}$ integrals for the spin-dependent amplitudes of Eqs. (36) and (37) are convergent.

Due to the Regge behavior of the $F_{1}$ proton structure function given by the Pomeron exchange, the $\nu_{\gamma}$ integrals are divergent, and the DR for the amplitude $f^{2\gamma}_{+}$ 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 Carlson:2008ke ; Carlson:2011af , see right panel in Fig. 6.

The forward TPE amplitude can be evaluated considering the lower blob of the TPE graph in Fig. 1 as a forward VVCS process on a proton: $\gamma^{\ast}\left(q,\lambda_{1}\right)+N\left(p,\lambda\right)\to\gamma^{\ast}\left(q,\lambda_{2}\right)+N\left(p,\lambda^{\prime}\right)$, which is shown in Fig. 7.

The forward VVCS amplitude $T_{\lambda_{2}\lambda^{\prime},\lambda_{1}\lambda}$ can be written in terms of the forward VVCS tensor $M^{\mu\nu}$ as

[TABLE]

where $\varepsilon_{\nu},~{}\varepsilon^{\ast}_{\mu}$ denote the virtual photon polarization vectors, $N,\bar{N}$ the proton spinors, and $\lambda_{1},\lambda_{2}~{}(\lambda,\lambda^{\prime})$ the photon (proton) helicities.

The forward VVCS tensor $M^{\mu\nu}$ can be expressed as a sum of the symmetric (spin-independent) $M^{\mu\nu}_{S}$ and the antisymmetric (spin-dependent) $M^{\mu\nu}_{A}$ parts:

[TABLE]

[TABLE]

with the forward Compton amplitudes $\mathrm{T}_{1},~{}\mathrm{T}_{2},~{}\mathrm{S}_{1},~{}\mathrm{S}_{2}$, which are functions of $\nu_{\gamma},~{}Q^{2}$ and enter Eq. (102) in a gauge-invariant way, e.g., $q_{\mu}M^{\mu\nu}=q_{\nu}M^{\mu\nu}=0$. The optical theorem relates the imaginary parts of the forward Compton amplitudes to the proton SFs by

[TABLE]

The forward lepton-proton scattering TPE amplitude can be expressed in terms of the forward VVCS amplitude $M^{\mu\nu}$ as

[TABLE]

with the leptonic tensor $\tilde{L}^{\mu\nu}_{h^{\prime}h}$:

[TABLE]

The expression for the spin-independent forward TPE amplitude $f_{+}^{2\gamma}$ is given by

[TABLE]

This result at threshold is in agreement with Ref. Carlson:2011zd .

The expressions for the spin-dependent forward TPE amplitudes $f_{-}^{2\gamma},~{}g^{2\gamma}$ are given by

[TABLE]

Evaluating these expressions in the nucleon rest frame at threshold we obtain:

[TABLE]

in agreement with Eq. (44) and the standard HFS derivation of Eq. (52) Iddings:1959zz ; Iddings:1965zz ; Drell:1966kk ; Faustov:1966 ; Faustov:1970 ; Bodwin:1987mj .

Exploiting DRs for the spin-dependent Compton amplitudes $\mathrm{S}_{1}$ and $\mathrm{S}_{2}$:

[TABLE]

we obtain the same expression for the amplitude $f^{2\gamma}_{-}$, see Eq. (LABEL:fminus), as with the DRs for the forward $lp$ amplitudes. The difference in the amplitude $g^{2\gamma}$, see Eq. (42), is given by the $\omega$-independent term, which is just the constant real part in $g^{2\gamma}$ amplitude. It vanishes with account of the BC sum rule. However, if one uses the DR for the amplitude $\nu_{\gamma}\mathrm{S}_{2}$ Hagelstein:2015egb , the result for $f^{2\gamma}_{-}(m)$ coincides with the expression $-g^{2\gamma}(m)$.

The forward TPE amplitudes evaluated within the DR approach coincide with the amplitudes evaluated with a help of the DRs for the forward VVCS amplitudes. However, the contribution of an individual TPE intermediate state differs in these two approaches.

Appendix E Forward scattering observables

The forward unpolarized elastic scattering cross section in the c.m. reference frame is given by

[TABLE]

with the electron scattering angles $\Omega$ and the azimuthal angle $\theta$.

All possible single-spin asymmetries are zero for the scattering in the forward direction. We denote the lepton spin asymmetry for the scattering on the polarized proton as $A$ and for the scattering on the unpolarized proton with the polarization transfer to the final proton as $P$. The asymmetries for the longitudinally polarized lepton and the longitudinally polarized proton in the forward scattering are expressed as

[TABLE]

The asymmetries for the transversely polarized lepton and the transversely polarized proton in the forward scattering are expressed as

[TABLE]

The asymmetries in the case of one transverse and one longitudinal polarizations vanish.