Spin-dependent Pomeron and Odderon in elastic proton–proton scattering

We introduce a new model of near-forward elastic proton–(anti)proton scattering at high energy based on the modern formulation of Pomeron and Odderon in terms of Wilson lines and generalized transverse momentum dependent distributions. We compute the helicity-dependent elastic amplitudes ϕ1,2,3,4,5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi _{1,2,3,4,5}$$\end{document} in this model and study their energy dependence from the nonlinear small-x evolution equations. While both Pomeron and Odderon contribute to helicity-flip processes in general, in the forward limit t=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t=0$$\end{document} only the double helicity-flip amplitude ϕ2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi _2$$\end{document}, dominated by the spin-dependent Odderon, survives. This may affect the extraction of the ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document} parameter as well as the total cross section in the LHC energy domain and beyond.


Introduction
The elastic proton-(anti)proton scattering at high energies becomes an important source of information about the multilayer proton structure [1]. While the gluon-driven exchanges are dominant at asymptotically high energies (small-x), an elastic scattering implies, at least, a (colour-singlet) pair of correlated gluons propagating in the t-channel known as the QCD Pomeron (see e.g. Ref. [2] and references therein), in analogy to the leading pole exchange with Regge trajectory of the highest intercept [3] (for more detailed on the Regge theory, see Ref. [4]). Even larger numbers of interacting gluons can be exchanged in an elastic scattering process, but the role of such multi-gluon interactions in elastic scattering yet remains uncertain, particularly, from the QCD point of view.
An odd-number gluon exchange starting from the leading triple-gluon one corresponds to the crossing-odd Odderon contribution in the Regge picture [5,6] (see also Ref. [7]). It a e-mail: yhatta@bnl.gov (corresponding author) was proposed back in the 70s in Ref. [8] that the Odderon contribution may be non-negligible compared to that of the Pomeron in the high-energy limit. However, while an experimental observation of the Odderon is yet unavailable, an exact magnitude and characteristics of such an elusive effect from theoretical viewpoint remain largely unknown and are the subjects of an intense debate and even controversial statements in the literature.
The recent outbreak of Odderon activity (see e.g. Refs. [9][10][11]) is largely triggered by the precision TOTEM data at the highest energy of the LHC, √ s = 13 TeV, on total σ tot [12] and differential dσ/dt [13] pp cross sections, as well as on the real-to-imaginary ratio of the elastic nuclear amplitude at the optical point, the so-called ρ-parameter [14]. Introducing the total helicity non-flip elastic amplitude T (s, t) as a function of the total c.m. energy squared s and four-momentum transfer squared t, the basic measurable quantities of the elastic scattering read so that the ρ-parameter is related to the total and differential (at vanishing momentum transfer) cross sections as follows where M is the proton mass. The ρ-parameter is small at TeV energies, ρ ∼ 0.1, and has been extracted by the TOTEM collaboration in Ref. [14] from the experimental data on dσ/dt near t ≈ 0 using the Coulomb-Nuclear Interference (CNI). As long as ρ(s) is known with sufficiently high precision, Eq. (2) is used to determine σ tot (s). The dominating claim is that a growth of the total cross section with energy, together with a decreasing ρ-parameter, as well as a qualitative difference of differential cross sections of pp and pp collisions [15][16][17], all are associated with the Odderon effect. There are some concerns in the literature, however, about the validity of the experimental procedure of ρ extraction (see e.g. [18,19]) and to the Odderon interpretation of its decrease with energy (see e.g. Refs. [11,20]), and hence more care is needed to justify the magnitude and the significance of the Odderon effect in the ρ measurement. In off-forward kinematics, a substantiated claim about the Odderon effect and its significance is made recently from the shape analysis of the elastic differential pp and pp cross sections based upon their scaling properties in Ref. [21]. In the current study, we instead consider a possible Odderon effect and its energy dependence at the optical point of vanishing t ≈ 0 only and leave the analysis of t-dependence for a future work. The usual rationale, similarly to the Pomeron, is that the Odderon is assumed to not flip the helicities of the scattered hadrons. Indeed, the existing theoretical formulations and the procedure of ρ extraction from the experimental data itself heavily rely on the presumption about an absence or a large suppression of helicity-flip processes at high energies. In this work, we question this convention and, in particular, explore a viable possibility that the helicity-flip elastic amplitude may be non-negligible at high energies. To our knowledge, this has neither been confirmed nor disproved by direct experimental measurements in the TeV region. On the other hand, it has been suggested in the literature that the Odderon can contribute to helicity-flip amplitudes [22][23][24][25]. Yet, the exact treatment of the problem has been difficult due to the lack of a systematic way to connect the Pomeron and Odderon with the spin degrees of freedom of the scattering (composite) particles such as protons in QCD. 1 Recently, however, there has been a significant progress in our understanding of the interplay between the Odderon and the proton spin [29][30][31][32][33][34]. In the Deep Inelastic Scattering (DIS) at small-x, the Color Glass Condensate (CGC) framework [35] provides a consistent description of the Pomeron and Odderon in terms of Wilson line correlators. Their couplings with various proton polarization states can be completely parametrized by the generalized transverse momentum dependent distributions (GTMDs) [34]. Indeed, the gluon Sivers function [36] at small-x is connected to the Odderon in the forward limit [29][30][31] and participates in the proton helicity-flip reactions including the unpolarised elastic scattering processes. In particular, it has been observed that the so-called spin-dependent Odderon [29] can flip the proton helicity even in the forward limit, and this effect can 1 For an earlier discussion of the Pomeron helicity flip observables for phenomenological scalar, vector and tensor Pomeron-proton couplings, see Ref. [26], and also, Refs. [27,28]. survive at high energies since the Odderon intercept is exactly equal to unity [7].
Motivated by these developments, in this paper we introduce a new model of near-forward elastic proton-proton scattering designed for the TeV region and beyond. By treating one of the protons within the quark-diquark model, we can devise a setup analogous to DIS in the so-called dipole frame. In this frame, helicity-flip amplitudes can be calculated by exchanging the spin-dependent Pomeron and Odderon. We then study their energy dependence at t = 0 by numerically solving the small-x evolution equations for Pomeron and Odderon. Of course, in near-forward pp scattering there is no apparent hard scale (like the photon virtuality Q 2 in DIS) which guarantees the use of perturbative approaches. However, in the TeV region one can consider the saturation momentum Q s as a dynamically generated hard scale.
The paper is organised as follows. In Sect. 2, we introduce the basic helicity amplitudes of the elastic pp scattering and discuss their role at high energies. In Sect. 3, we derive the helicity amplitudes in the quark-diquark dipole model and discuss their main properties. In Sect. 4, energy dependence of the helicity amplitudes and their ratios is numerically studied from the nonlinear small-x evolution equations. Finally, a brief summary and concluding remarks are given in Sect. 5.

Helicity amplitudes
Consider near-forward proton-proton elastic scattering P 1 P 2 → P 1 P 2 at high energies schematically shown in Fig. 1, with 4-momenta satisfying where i = 1, 2 denotes the transverse momentum components. We introduce the spin-dependent elastic amplitudes λ 1 λ 2 |T |λ 1 λ 2 [37][38][39][40] where λ = 2h = ±1 represents the helicity of each proton (multiplied by two, for convenience). These helicity amplitudes depend on s ≈ 2P as well as the azimuthal angle ϕ = Arg( 1 ⊥ + i 2 ⊥ ). We factor out the ϕ-dependence as 2 and switch to the commonly used notation while the total cross section is A general argument shows that φ 4 ∝ t, φ 5 ∝ √ −t as t → 0, whereas φ 1,2,3 go to a constant in this limit [40,41]. (Here we focus on the QCD part of the amplitude. The QED part behaves differently, see Appendix A.) Given these limiting behaviors, it is convenient to rescale the spin-dependent amplitudes as [24] where The complex functions r i=2,4,5 have a finite limit as t → 0. They can be experimentally accessed by measuring various spin asymmetries [24]. For example, r 5 is closely related to single spin asymmetry A N , and r 2 is related to double spin asymmetry A N N . The results from fixed-target experiments at RHIC at √ s = 13.76 GeV and 21.92 GeV [42,43] indicate that the parameters r 2 , r 5 are small, of order 10 −3 in this lowenergy region. There are also RHIC data in the collider mode at √ s = 200 GeV [44,45]. The analysis mostly focused on A N and a rather small value of r 5 has been reported.
At higher energies, however, nothing is known about the fate of the helicity-flip amplitudes since there is no polarized proton collider beyond the RHIC energies. They are rarely discussed in connection with the ongoing measurements at the LHC, or with the earlier measurements at the Tevatron. It is usually assumed, often without even mentioning it, that φ 1 ≈ φ 3 and φ 2,4,5 ≈ 0 for all values of t. There is then only one (complex) amplitude T = 8πφ 1 , and (6) and (7) reduce to the formulas mentioned in the introduction. Yet, even in unpolarized scattering, the helicity-flip amplitudes affect the observables. In the presence of nonvanishing φ 2 , Eq. (2) should be modified as Also, in the non-forward scattering with |t| > 0, φ 2,4,5 amplitudes could affect the shape of dσ/dt, especially, in the dip region where |φ 1,3 | become small. In the next section, we compute all the φ's in a model which incorporates Pomeron and Odderon in the Wilson line formulation of small-x QCD. We do not make the usual assumption that the helicity-flip amplitudes φ 2,4,5 are negligibly small. As we shall see very clearly below, the spindependent Pomeron and Odderon exchanges naturally generate non-negligible helicity-flip amplitudes. The latter are a priori not suppressed at high energies since they share the same energy dependence ('Regge intercept') as for the helicity-conserving ones.

Elastic scattering in the quark-diquark model
In this section, we calculate the helicity amplitudes φ 1,...,5 in the dipole model of high-energy pp (and pp) scattering illustrated in Fig. 1. Our setup is similar to the description of DIS at small-x in the so-called 'dipole frame' where the virtual photon fluctuates into a quark-antiquark pair long before interacting with the target proton. Specifically, we work in an asymmetric frame in which P + 1 P − 2 . The 'slow', leftmoving proton 2 is treated in the quark-diquark model [46]. It fluctuates into a quark and a scalar diquark, and the pair interacts with the shockwave created by the 'fast' proton 1 in the eikonal approximation. The corresponding scattering amplitude is given by Here, is the light-front wave function of the proton 2 fluctuation into a q −qq pair to be specified shortly, λ 1,2 and λ 1,2 denote helicities of protons 1,2 in the initial and final states, respectively, r ⊥ is the transverse distance between the quark and the diquark, z is the longitudinal momentum fraction of the proton 2 carried by the quark, and N is the so-called dipole scattering amplitude defined by in terms of a lightlike Wilson line in the fundamental representation which describes the quark scattering off the target color field, and that for the diquark, U † , which has the same color representation as an antiquark. As usual, g and N c = 3 denote the QCD coupling and number of colors. Following [34], we parametrize the dipole amplitude as where the gluon GTMDS f 1,n and g 1,n (n = 1, 2, 3) are functions of k 2 ⊥ , 2 ⊥ and |k ⊥ · ⊥ | as well as the Bjorken-x variable. At small-x, they come from the real and imaginary parts of the operator TrUU † and represent the Pomeron [47] and Odderon [48] exchanges, respectively. The apparent pole at k 2 ⊥ = 2 ⊥ /4 is innocuous because f and g are proportional to k 2 ⊥ − 2 ⊥ /4, see for example (36) below. Equation (14) describes the most general coupling between the Pomeron/Odderon and the proton, consistent with the symmetries of the GTMDs. By using the Gordon identity, one can write (14) as the linear combination of scalar and vector couplings. The tensor coupling [26] is absent in our framework.
Let us now work out the product of spinors explicitly. Up to corrections of order M/P + 1 and ⊥ /P + 1 , we getū where we introduced the 'polarization vector' In the r ⊥ -space, the parametrization takes the form, with and for n = 1, 2, 3. All functions defined above are real-valued. The wave function of proton 2 in the quark-diquark model is given by (see "Appendix B" for the relevant Feynman rules) where c s is a constant normalisation. The constituent quark has momentum fraction z, transverse momentum l ⊥ , mass m q and helicity λ q . The scalar diquark has momentum fraction z = 1 − z and mass m s . Due to a finite binding energy, m q + m s ≥ M, and this ensures thatM 2 ≡zm 2 q + zm 2 s − zz M 2 ≥ 0. After carrying out the integration over l ⊥ , one can rewrite the wave function in Eq. (23) as, This leads to the following expression for the wave function squared in the forward limit, in terms of the helicity flip and helicity non-flip parts of the wave function respectively. In the non-forward case, a nontrivial phase e i z− 1 2 ⊥ ·r ⊥ emerges 3 [49]. We keep the subleading terms up to quadratic order in ⊥ , so in practice we use 3 For the reader's convenience, here we briefly recapitulate the discussion in Ref. [49]. The non-forward amplitude in dipole models typically has the following structure in impact parameter space b ⊥ , where N is the dipole scattering amplitude. The shift b ⊥ → b ⊥ − zr ⊥ is caused by the phase factor e izr ⊥ · ⊥ which generically appears in non-forward impact factors [50], see Eq. (36) below for example. This implies that b ⊥ + (1 − z)r ⊥ and b ⊥ − zr ⊥ can be interpreted as the coordinate of the quark and antiquark (or diquark), respectively. We can thus identify where we used the abbreviation z * = z − 1/2. Assembling the above pieces together, we get where we defined and Footnote 3 continued This gives Above, f 1,1 and E are nothing but the GTMD version of the GPDs of the fast proton (c.f., Eq. (4.48) of Ref. [51]) normalized as where H and E are the standard gluon GPDs. Note that, since we are colliding identical particles, by symmetry the coefficients of δ λ 2 ,−λ 2 δ λ 1 ,λ 1 and δ λ 2 ,λ 2 δ λ 1 ,−λ 1 have to be equal (up to a sign and trivial relabeling). However, in the asymmetric frame in which we are working, this is not obvious at first sight. While we do not have an explicit proof, we nevertheless argue that the two expressions are indeed equivalent. The functions n and f introduced in Eqs. (26) and (27) are related to the helicity non-flip and helicity flip parts of the gluon GTMD of the slow proton, respectively, The t-channel gluon propagators in Eq. (36) (as well as the small-x evolution) are absorbed intoẼ andH . Thus, the terms proportional to nẼ and fH in the last two lines of Eq. (32) are both the convolution of the H -type GTMD of one proton and the E-type GTMD of the other proton, and are thus equal. A similar argument applies to the imaginary parts proportional to ng1,2 and fg1,1 . Although there is in general no relation betweeng 1,2 and g 1,1 , they satisfy the same evolution equation. The only difference is the way the t-channel Odderon amplitude T O couples to the proton, and this coupling is proportional to f and n , respectively, cf., Ref. [30]. Thus, the imaginary terms in the last two lines of Eq. (32) both have the structure n ⊗ T O ⊗ f , and are thus equivalent. After removing the phase according to Eq. (4), we arrive at The sign in front ofg 1,2 in Eq. (40) has been fixed using the relation + + |T | + − = − + + |T | − + [52] 4 . We immediately notice that φ 1,3 are purely imaginary and φ 2 is purely real. Therefore, the usual ρ-parameter (1) vanishes at t = 0 in this model. Away from t = 0, the ρ-parameter is dominated by the spin-independent Odderong 1,1 . We also see that the Pomeron (H ,Ẽ) and Odderon (g 1,2,3 ) contributions are always relatively imaginary. This means that there is no interference when squaring the amplitudes |φ i | 2 , and dσ/dt is insensitive to the sign ofg 1,2,3 . In other words, dσ/dt is identical for pp and pp scatterings in this model. Recently, there are indications that the difference dσ pp /dt − dσ pp /dt is nonvanishing from an analysis of the LHC and Tevatron data [21,53]. In order to explain this, the Odderon has to have a small imaginary part (and the Pomeron has a small real part). It may be possible to generalize our model to accommodate this effect, for example, by using the dispersion relation or invoking Regge theory or the AdS/CFT correspondence [54]. This is however beyond the scope of this work.
As for the ratios (8), we get Since φ 4,5 vanish at t = 0, r 4,5 are not well-defined at t = 0, and of course measurements are always performed at t = 0. On the other hand, r 2 has a well-defined limit t → 0, and there we need to subtract (2π) 2 δ 2 ( ⊥ = 0) ≡ A, the transverse area of the proton, fromH in the denominator. This converts the S-matrix (H ) into the T -matrix, and is crucial for the r ⊥ integral to converge at small r ⊥ . [Note that n, f (r ⊥ ) ∼ 1/r 2 ⊥ at small-r ⊥ .] When t is nonzero, the subtraction is absent but there is no convergence problem sinceH vanishes at r ⊥ = 0 if ⊥ = 0. However, we can keep this subtraction in the denominator of r 4,5 and evaluate it at t = 0 thanks to the fact that the limit Im φ 1 (t → 0) is smooth. Also,H in the numerator of r 5 can be safely evaluated at t = 0 since the factor r 2 ⊥ kills the divergence at r ⊥ = 0. In the present work, these tricks are crucial for the numerical study in the next section since we do not have a numerical solution ofH at t = 0.
We see that the real parts R 2,4,5 entirely come from the Odderon. In particular, R 2 at t = 0 is nonvanishing due to the spin-dependent Odderon g 1,2 , and this can contribute to the differential and total cross section according to Eq. (10). The imaginary parts I 4,5 come from the Pomeron and I 2 vanishes in this model. It is interesting to notice that I 5 is parametrically of order unity if the typical value is r ⊥ ∼ 1/M. However, at high energy the integrand is more localized at small-r ⊥ , and then the factor r 2 ⊥ leads to a suppression of I 5 (see below). We also expect |I 5 | |I 4 | assuming |H | |E|.

Energy dependence of the helicity amplitudes
In this section, we study the center-of-mass energy √ s dependence of the helicity amplitudes φ i and their ratios obtained in the previous section. f 1,n and g 1,n are the real and imaginary parts of the dipole scattering amplitude (12), respectively. The latter satisfies the Balitsky-Kovchegov (BK) equation [47,55] which is an evolution equation in ln s including the gluon saturation effect. Thus, the Pomeron and Odderon amplitudes can be obtained from the real and imaginary parts of the BK equation with appropriate initial conditions [48,56]. We restrict ourselves to the forward limit ⊥ = 0, which means that we concentrate on f 1,1 and g 1,2 . Solving the BK equation with finite ⊥ is numerically more involved, and to our knowledge this has not been done for the Odderon.
Admittedly, the use of the BK equation for our problem must be legitimately criticized. Being an equation originally derived in perturbation theory, in principle the BK equation can only apply to processes which involve a hard scale. However, in near-forward elastic pp scattering, apparently there is no such hard scale. Yet, the idea of gluon saturation and the Color Glass Condensate [35] is that at asymptotically high energies, the gluon distribution in the colliding particles is characterized by a dynamically generated hard scale, called the saturation momentum Q s (s) which is an increasing function of √ s. It has been demonstrated [57] that the theory of saturation with Q 2 s ∝ ( √ s) 0.23 > 1 GeV 2 can successfully describe the multiplicity and the mean p T distribution in pp collisions across the whole energy region of the LHC, 0.9 GeV < √ s < 7 TeV. This partly justifies our approach at least for the Pomeron, and allows us to calculate the perturbative part of the growth of the total cross section with energy. Of course there are also nonperturbative contributions to the total cross section, but in our model these are absorbed into the parameter A. As a matter of fact, the same argument does not quite hold for the Odderon. It has been noticed that the characteristic momentum scale of the Odderon amplitude does not grow like Q s [33,58]. Therefore, the results involving Odderon below are at best a crude estimate of the possible high energy behavior suggested by perturbation theory. In reality the dominance of the nonperturbative effects may be overwhelming.
We basically follow Ref. [33] for the numerical evaluation of f 1,1 and g 1,2 , except that we now include the running coupling effect. Reference [33] considered a transversely polarized proton and studied the gluon Sivers function which is the forward limit of g 1,2 . On the other hand, in our problem the proton is longitudinally polarized. We thus need a little spinor algebra to connect the two works. Let us return to Eq. (14) and take the forward limit ⊥ = 0, Here we assume that the proton is transversely polarised, with the transverse spin vector S ⊥ normalised as | S ⊥ | = 1.
In the r ⊥ -space, This can be written as (compare with Eq. (5) of [33]) where for the Pomeron and the spin-dependent Odderon components of the dipole S-matrix, respectively. We compute P(s, r ⊥ ) and Q(s, r ⊥ ) as functions of the center-of-mass energy squared s from the solution of the BK equation with running coupling as prescribed in Ref. [59]. Then, using Eq. (47) we accessH (r ⊥ ) andg 1,2 (r ⊥ ) that are further employed in computing r 2,5 through Eqs. (41) and (43). We adopt the following form for the coupling constant  Fig. 2 The energy dependence of the total cross section. We take Q s0 = 1.0 GeV to fit the experimental data with b 0 = 9 4π (corresponding to n f = 3), = 0.241 GeV and a = e 8π 9 . The initial conditions are given at the starting energy scale s 0 as follows: The initial saturation scale Q s0 is expected to be around 1 GeV in the TeV region, while the strength of Odderon κ is an unknown parameter including its sign (see, however, [29]) which should be fitted to the data [34]. The other parameters in this model are m q , m s and c 2 s A (only this product enters our observables). We fix m q = 0.3 GeV and m s = M − m q , while c 2 s A is fitted to the total cross section. The energy dependence of the total cross section computed in our approach is shown in Fig. 2 with Q s0 = 1.0 GeV. Here, the green and pink lines denote two different values of the starting energy scale, √ s 0 = 0.1 and 0.5 TeV, respectively. The result is in reasonable agreement with the corresponding measurements in pp collisions performed at several distinct energies, such as those by the TOTEM LHC Collaboration at 13 TeV [12] [16], as well as in pp collisions by D0 Tevatron Collaboration at 1.96 TeV [63] and by UA4 CERN SPS Collaboration at 546 GeV [64] and 630 GeV [65]. Since the measured values for σ tot (s) are sometimes not available in the experimental articles, in those cases the σ tot values have been taken from the global Lévy analysis of the corresponding elastic pp and pp cross section data performed recently in Ref. [17]. Incidentally, we have also tried Q s0 = 0.5 GeV, but the quality of the fit is noticeably worse in this case.
The results for R 2 , I 5 and R 5 are plotted in Fig. 3 as functions of √ s in upper-left, upper-right and bottom panels, respectively. Note that the normalization and sign of R 2,5 are arbitrary, as it is proportional to the unknown parameter κ, and we have chosen R 5 to be negative following the recent suggestion in [66]. Irrespective of this, we can predict that R 2 and R 5 have the same sign and that |R 2 | is roughly two times larger than |R 5 |. We also see a clear tendency that the magnitude of R 2,5 decreases with increasing energy. This is because, although the Odderon intercept is unity in the dilute (BFKL) regime, the nonlinear saturation effect tends to suppress the Odderon amplitude [33,48,58]. On the other hand, the value of I 5 is a prediction of this model, since both the denominator and numerator of (43) come from the Pomeron. It is negative and the magnitude decreases with energy because of the factor r 2 ⊥ in the numerator of (43): The r ⊥ -integral is dominated by r ⊥ ∼ 1/Q s (s), and Q s (s) is an increasing function of energy.
The data on single spin asymmetry A N in small-angle elastic pp collisions have recently become available from the fixed-target measurement HJET at BNL [67] as well as earlier from the STAR measurements of polarized elastic pp collisions at √ s = 200 GeV [44]. These data have enabled to extract the real and imaginary parts of r 5 ratio in a wide energy domain. The values of R 5 published by the experimental collaborations were found (by STAR measurement and by an extrapolation from the lower HJET energies) to be either small positive or consistent with zero at √ s = 200 GeV (see also [26]), while the hadronic contribution predicted in Fig. 3 (upper-right panel) is found to be larger than the ballpark of experimental values.
We note, however, that the CNI contribution has to be taken into consideration as its impact on r 5 can be rather important, whereas the current analysis only focuses on the hadronic contribution to φ 5 . Indeed, as was recently advocated in Ref. [66] relying on a Regge analysis and a dominance of the Pomeron spin-flip contribution, the absorptive corrections to the Coulomb spin-flip amplitude significantly modify the CNI mechanism. As a result, this modification affects the extracted values of r 5 , in particular making the spin-flip Pomeron I 5 rather large and negative, at the level of −5 to −10 % at √ s = 200 GeV, in consistency with expectations [68]. The fact that our QCD-based approach predicts non-vanishing and negative I 5 is encouraging, although as we explained above it falls with energy, in contrast to the behavior predicted by the Regge fit of Ref. [66]. These results are not inconsistent and rather suggest that the gluon saturation regime has not been reached at RHIC energies. We leave a thorough analysis of the CNI effects in the current framework for a future work.

Conclusions
In this work, we have presented a new QCD-inspired model for small-angle elastic proton-(anti)proton scattering in terms of spin-dependent Pomeron and Odderon helicity amplitudes in the dipole picture based upon the Wilson line approach. The elastic amplitudes φ 1,...,5 are effectively described in near-forward kinematics by means of a scattering of the lowest Fock state p → q + (qq) of projectile proton (i.e. the quark-diquark dipole) off the proton target, i.e. in a similar fashion as DIS. The corresponding dipole S-matrix receives contributions from non-flip and spin-flip Pomeron and Odderon exchanges that are represented in terms of GTMDs of different types.
Connecting to the numerical analysis of the small-x Odderon evolution equation performed earlier in Ref. [33] and incorporating in addition the QCD running coupling effect, we explore the relative importance of spin-flip contributions to the elastic pp scattering at high energies. In particular, we analyse the energy dependence of the spin-flip Pomeron (I 5 ) and spin-flip Odderon (R 5 ) amplitudes, as well as doublespin-flip Odderon (R 2 ) amplitude relative to the non-flip one. At variance with an earlier Regge-based calculation of Ref. [66] incorporating for the first time the absorptive corrections in the CNI mechanism, we do not assume that the exchanged spin-independent and spin-dependent Regge trajectories have different intercepts and do not neglect the Odderon contributions. Yet, we have reached a qualitatively similar conclusion about a significant and negative contribu-tion to the single helicity-flip amplitude I 5 . Moreover, the measured value of R 5 can be used to determine the Odderon coupling κ, which in turn determines the value of R 2 . The energy dependence of r 5 in our approach is decaying and hence is strictly opposite to the steeply rising behavior from the Regge analysis [66] obtained in the lower energy region. This suggests that once the gluon saturation effect kicks in, the behavior of r 5 changes. A further analysis of this issue is certainly needed.
The experimentally probed energies in the existing measurements of the spin-flip contributions may not be high enough to make a conclusive statement about the energy dependence of spin-dependent Pomeron and, especially, Odderon effects. Indeed, at such low energies as √ s = 200 GeV the C-odd effects may come mostly from secondary Reggeon exchanges, not due to spin-dependent Odderon studied in our analysis here. It is therefore of high importance to perform a new measurement of r 2 and r 5 in a TeV energy range to make a definite conclusion about the energy dependence of spin-dependent Pomeron and Odderon in the future. Note that this does not necessarily require polarized proton beams which are not available at the LHC. The differential cross section (6) gets contributions from the helicity-flip amplitudes, but they are usually ignored in the CNI analysis. It would be very interesting to test more flexible parametrizations of the CNI effect including the hadronic and electromagnetic contributions to φ 2,4,5 . This could eventually affect the value of the ρ-parameter, and also the total cross section via (10).
Finally, it is of course necessary to extend the present calculation to finite momentum transfer t, in particular up to the 'dip' region of dσ/dt. The basic formulas are given in (37)- (40), but we are missing models of the spin-independent and spin-dependent Pomeron and Odderon amplitudes at finite t (see for example [69,70] for a model of g 1,1 at finite impact parameter). They can also serve as an initial condition for the impact-parameter dependent BK equation to determine the energy dependence. It is also interesting to consider different models for the 'slow' proton such as a bound state of three quarks. We hope to address these issues elsewhere.