Investigating the Timelike Compton scattering in photon-proton interactions at the EIC

Information about the 3-dimensional description of the quark and gluon content of hadrons, described by the generalized parton distributions (GPDs), can be probed by exclusive processes in electron - proton ($ep$) collisions. In this letter, we investigate the timelike Compton scattering (TCS) in $ep$ collisions at the future electron - ion collider (EIC). Such process is characterized by the exclusive dilepton production through the subprocess $\gamma p \rightarrow \gamma^*p \rightarrow l^+ l^- p$, with the real photon in the initial state being emitted by the incoming electron. Assuming a given model for the GPDs, the TCS differential cross-section is estimated, as well the contribution associated to the interference between the TCS and Bethe - Heitler (BH) amplitudes. Predictions for the TCS, BH and interference contributions are presented considering the kinematical range expected to be covered by the EIC detectors. Moreover, the polarized photon asymmetry is also studied. Our results indicated that a future experimental analysis, considering photon circular polarizations, can be useful to probe the interference contribution and constrain the description of the GPDs for the proton.

The tomography picture from the hadrons can be revealed in exclusive photon -hadron interactions, which are characterized by the fact that the hadron remains intact after scattering.Such processes can be studied in electron-hadron colliders as well as ultraperipheral collisions at RHIC and LHC (For recent reviews see, e.g.Refs.[1][2][3][4][5][6][7][8]) and has a scattering amplitude that can be expressed in terms of generalized parton distributions (GPDs) [9][10][11][12].As a consequence, they provide information about the 3-dimensional description of the quark and gluon content of hadrons [13,14].Important constraints about the GPDs have been derived from the experimental results obtained in ep collisions at HERA, JLab and COM-PASS for the exclusive vector meson production, deeply virtual Compton scattering (DVCS) and for the timelike Compton scattering (TCS), with the promising expectation that future electron -ion colliders at BNL (EIC) [4], China (EicC) [6] and CERN (LHeC) [7,8] will allow us to improve our understanding of the hadronic structure in a larger kinematical range.
In this letter we will focus on the timelike Compton scattering (TCS), which is the exclusive photoproduction of a lepton pair with large invariant mass, considering the kinematical range that will be probed in ep collisions at the EIC.In previous investigations, TCS have been studied in LHC and JLab kinematics [15,16].Our goal is to complement the studies performed in Refs.[17,18], where we have demonstrated that the TCS process can be investigated in ultraperipheral pA collisions at RHIC and LHC.In particular, we will consider the contribution induced by real photons, represented in Fig. 1, which is expected to be dominant in ep collisions, and present predictions for the dependence of the differential cross-section on the squared transferred momentum t, azimuthal angle ϕ, center of mass photon -proton energy W and invariant mass of the dilepton system Q 2 .Moreover, predictions for the photon beam circular polarization asymmetry will also be presented.We will show the TCS results, as well as those associated to the Bethe -Heitler (BH) process, which contributes to the same final state, and for the TCS -BH interference.
Initially, let's present a brief review of the formalism need to describe the TCS process in terms of generalized parton distributions (For a detailed discussion see, e.g., Refs. Dilepton production in photon-proton interactions at the EIC.[13,14,[19][20][21]).As in these previous studies, we will assume the validity of the factorization theorem, which allow us to separate the hard scattering process, calculated using perturbation theory, from the non-perturbative dynamics encoded in GPDs.As in Ref. [17], in this exploratory study, the TCS process will be estimated at leading order of the strong running coupling constant α s , which implies that the Compton amplitude is dominated by the quark handbag diagrams described in terms of the quark GPDs1 .The impact of the contributions proportional to gluon GPDs will be presented in a forthcoming study.As the Bethe -Heitler (BH) process generates the same final state than the time-like Compton scattering, both processes contribute at the amplitude level, which implies that the differential cross-section for the dilepton production in ep collisions will be expressed as a sum of three contributions: , where INT denotes the term associated to the TCS -BH interference.The BH contribution had been estimated in Ref. [19], where the differential cross-section is explicitly presented.In addition, as a detailed derivation of the main formulas associated to the TCS and INT contributions was presented in Refs.[19,22,23], here we will only present the final expressions, which have been used in phenomenological studies of the TCS process at the leading order in α s .One has that the TCS contribution is given by [19] dσ where η = Q 2 /(2s γp − Q 2 ), s γp = W 2 , M p is the proton mass and the Compton form factors H 1 , H 1 , E 1 , and E 1 are expressed in terms of the hard -scattering kernels T q,g and the GPDs H, H, E and Ẽ, defined in Ref. [13] (For details, see Refs.[15,19]).At leading order, one has that T g = 0 and the expressions for the quark sector can be found in Refs.[15,19].The NLO kernels are presented in Refs.[22,23].As statement in Ref. [19], the θ in above equation is the angle between virtual photon and lepton in the lepton pair cms frame while the ϕ angle is the angle between lepton pair plane and photon-proton plane.On the other hand, the contribution associated to the interference between the TCS and BH processes, considering unpolarized protons and photons, can be expressed as follows [19,24] ), and we have neglected the lepton mass and assumed that s γp , Q 2 ≫ t, M 2 p .Moreover, F 1 (t) and F 2 (t) are the usual Dirac and Pauli form factors, with F 2 (0) normalized to the anomalous magnetic moment of the proton.One has that differently from the BH and TCS contributions, which are even under the transformation ϕ → ϕ + π for integration limits symmetric about θ = π/2, the interference term is odd due to charge conjugation.Such aspect allow us to probe the Compton amplitude through a study of 2π 0 dϕ cos ϕdσ/dϕ.One has that for the case of a scattering between unpolarized protons and photons, one probes the real part of the Compton helicity amplitudes.On the other hand, the imaginary part can be accessed with polarized photon beams.If the photons have a circular polarization ν, as is the case for a bremsstrahlung beam emitted from longitudinally polarized leptons, one has that the interference term will be given by [19] One has that this interference contribution can be probed by investigating the photon beam circular polarization asymmetry, A LU , defined by The main ingredient to estimate the TCS and INT cross-sections is a realistic model for the GPDs.In this exploratory study, we will consider the GPDs proposed and detailed in Refs.[25][26][27], usually called Goloskokov -Kroll (GK) model, which is based on fits of meson electroproduction data.In this model, the GPDs are expressed by ( where F = H, H, E, Ẽ and i denotes the quark flavor of the double distribution f i , which are given by the following expressions One has that n i is set to 1 for valence quarks and 2 for sea quarks.Moreover, h q sea (β, 0) = q sea (|β|)sign(β) and h q val (β, 0) = q val (|β|)Θ(β), where q sea and q val are the usual unpolarized PDFs.Finally, the t-dependence of the double distributions are described by g i (β, t), which is assumed to have a Regge behavior with linear trajectories, being given by The parameters considered in our analysis are detailed in Refs.[25][26][27].It is important to emphasize that the GK model provides a satisfactory description of the recent TCS data measured by the CLAS Collaboration [28].
In what follows, we will present our predictions for the exclusive dilepton production considering the TCS, BH and interference contributions.We will initially consider unpolarized protons and photons and assume W = 60 GeV, which is a typical value for the photonproton center -of -mass energy, expected to be probed in future ep collisions at the EIC.The results for the |t| -distribution, obtained considering two values of Q 2 , fixed azimuthal angle ϕ and integrating θ from π/4 to 3π/4, are presented in Fig. 2. In agreement with  the predictions derived in previous studies [15,17,19], the BH contribution is larger than the TCS and INT contributions.However, our results indicate that the INT contribution is not negligible in the kinematical range considered.Moreover, one has that the cross-section decreases with the increasing of |t| and Q 2 .
In Fig. 3 we present our predictions for the dependence on ϕ of the differential crosssection, derived considering W = 60 GeV, Q 2 = 5 and 10 GeV 2 , θ integrated from π/4 to 3π/4 and |t| over the full kinematical range.One has that TCS contribution is a constant, since it is independent of ϕ.On the other hand, as the unpolarized INT contribution is proportional cos ϕ, one has that it is larger for ϕ ≈ π.We have verified that the INT contribution is about 4 percent of the total cross-section in the kinematical range considered.Thus, this contribution is not negligible in the total cross-section.Moreover, one has that the peak for ϕ ≈ π in larger for smaller values of the dilepton invariant mass Q 2 .
Let's now analyze the behavior of the observable A LU , which can accessed making use of photon circular polarizations.The ϕ dependence of the photon polarized asymmetry A LU (W, Q 2 , ϕ) is presented in Fig. 4, derived by integrating θ in the range [π/4, 3π/4] and |t| in the full kinematical range, for two values of Q 2 (left panel) and two values of W (right panel).The proportionality of A LU with sin ϕ is clearly observed.One has that it reaches ≈ 10% for Q 2 = 5 GeV 2 , decreasing for larger values of Q 2 .On the other hand, one has that A LU is almost energy independent in the range considered.
Finally, in Fig. 5, we present our predictions for the |t| (left panel) and W (right panel) dependencies of A LU (W, Q 2 , |t|, ϕ).One has that the largest value of A LU occurs for |t| ≈ 0.2 GeV 2 .Moreover, the dependence on W is similar for the two values of Q 2 considered, with the normalization decreasing with Q 2 .These results indicate that a future measurement of the photon beam circular polarization asymmetry for small values of the invariant mass of the dilepton system and |t| ≈ 0.2 GeV 2 can be useful to probe the INT contribution and probe the description of the GPDs.
As a summary, the improvement of our understanding about the quantum 3D imaging of the partons inside the protons and nuclei is one of the main goals of Particle Physics.The quantum information of how partons are distributed inside hadrons is encoded in the quantum phase space Wigner distributions, which include information on both generalized parton distributions (GPDs) and transverse momentum parton distributions (TMDs).In this letter, we have investigated the timelike Compton scattering, which is one of the promising process to probe the quark and gluon GPDs of the proton.In particular, we focus on the kinematical range that will be probed by the future Electron -Ion Collider at the BNL/USA.The contributions of the Bethe -Heitler, TCS and interference were estimated at leading order in α s and predictions for the kinematical range expected to be covered by the future EIC detectors were presented assuming the GK model for the quark GPDs.Our results indicate that the exclusive dilepton production in ep collisions is dominated by the BH contribution.However, if this contribution is subtracted, our results indicate that the TCS and INT contributions are not negligible and can be probed in future experimental analysis.We also provided predictions for the photon beam circular polarization asymmetry, A LU , which indicated that it is larger for smaller values of |t| and Q 2 .Such results indicate that a future experimental determination of the INT contribution is, in principle, feasible.

FIG. 2 .
FIG. 2. Predictions for |t| dependence of the differential cross-section dσ/dQ 2 d|t|dϕ associated with the exclusive dilepton production in ep collisions.The results for the total (BH + TCS + INT), TCS and INT contributions are presented separately.

FIG. 3 .
FIG. 3. Predictions for azimuthal angle distribution of the differential cross-section dσ/dQ 2 dϕ associated with the exclusive dilepton production in ep collisions.The results for the total (BH + TCS + INT), TCS and INT contributions are presented separately.

FIG. 4 .FIG. 5 .
FIG. 4. Predictions for the dependence of A LU on the azimuthal angle ϕ for different values of Q 2 (left panel) and distinct photon -proton center -of -mass energies (right panel).