Coulomb-Nuclear Interference in polarized pA scattering

We made the first attempt to understand the observed unusual t dependence of single-spin asymmetry observed in the HJET experiment at RHIC. Usually, the interaction of hadrons is presented as a long-range Coulomb interaction and a short-range strong interaction with Coulomb corrections. Such a division gives rise to a Coulomb phase of the hadronic term. Conversely, here we consider short-range hadronic interaction as a correction to the long-range electromagnetic term, i.e., we treat it as an absorptive correction. This significantly affects the Coulomb-nuclear interference, which is a source of single-spin azimuthal asymmetry at small angles.


Introduction
Elastic scattering is usually characterized by spin non-flip f n f and spin-flip f s f amplitudes, which determine the differential elastic cross-sections and single-spin azimuthal asymmetry A N (t), where t is 4-momentum transfer squared.The latter contains interference of the helicity amplitudes, which makes polarization effects sensitive to the hadron interaction dynamics [1,2].According to Equation ( 2), an important condition for azimuthal asymmetry is the existence of a phase shift between the spin amplitudes.However, to the best of our knowledge, the phases of f n f and f s f at high energies are similar.For example, in the Regge pole model, f h s f and f h n f have the same phase given by the signature factor.To maximize A N (t), one should combine hadronic and electromagnetic amplitudes [3].While the former is predominantly imaginary, the latter is nearly real.Apparently, a sizable effect is expected at small momentum transfer squared t, where the Coulomb and hadronic amplitudes are of the same order.The t dependence of asymmetry in pp scattering is described in a simplified approximation [3] by where y = −t/t pp p and Here, α em = 1/137 is the fine structure constant; µ p = 2.79 is the proton magnetic moment.The asymmetry A pp N (t) reaches a maximum (5) of about 4-5%, at t = t p ≈ 2 × 10 −3 GeV 2 .For the sake of simplicity, we assume here, like in [3], a pure non-flip and imaginary hadronic amplitude, no Coulomb phase, etc.In what follows, we present most of these simplifications.
The CNI asymmetry Equation (3) predicted in [3] was confirmed by measurements in [4,5].It is worth mentioning at least two important practical applications of the CNI effect.First, it is a unique opportunity to measure the spin-flip component of the hadronic amplitude at high energies [6,7].
Similar relations can be applied to proton-nucleus elastic scattering, with Correspondingly, the maximal value of A N at t = t pA p reads, A simple estimate σ pA tot ∼A 2/3 σ pp tot leads to non-dramatic modification of the CNI asymmetry, with about a 20% increase in A pA N (t pA p ) even for heavy nuclei, e.g., gold.The energy of 100 GeV in the Lab frame is not sufficiently high to suppress iso-vector Reggeons, which have quite a large (dominating) spin-flip component.This is why it is difficult to disentangle large Reggeon and small Pomeron spin-flip terms.On iso-scalar nuclei, e.g., carbon, copper, etc., iso-vector Reggeons are completely excluded, otherwise they are suppressed by a small factor (A − 2Z)/A.

Born Approximation
The elastic pA amplitude is fully described by two helicity amplitudes f n f ,s f defined in [2], each having hadronic and electromagnetic parts, f i (q T ) = f h i (q T ) + f em i (q T ), where for small-angle elastic scattering t ≡ −q 2 ≈ −q 2 T , they are comparable.Hadronic Born amplitudes can be represented as where r pA 5 is a nuclear analog of r pp For pp elastic scattering, r pp 5 was fitted to data in [5,7].The hadronic nuclear form factor F h A (q 2 T ) in Equations ( 9)-( 10) can be evaluated within the Glauber approximation [22], as well as the total nuclear cross-section, Correspondingly, the Born electromagnetic amplitudes read where the small fictitious photon mass λ is introduced to avoid infrared divergence.The final results are checked for stability at λ → 0. The nuclear electromagnetic form factor in Equations ( 14) and ( 15) has the form The nuclear thickness functions are defined as where ρ(b, z) is the nuclear density distribution function, and Following [22,23], next we replace the nuclear thickness function with a more accurate effective thickness function convoluted with NN elastic amplitude where where B hN is the slope of the differential hN elastic cross-section.This effective nuclear thickness function can be simplified to

Hadronic vs. Electromagnetic Amplitudes
The long-range Coulomb forces affect the strong-interaction amplitude by generating a phase shift, known as the Coulomb phase [24][25][26].The interplay of Coulomb and hadronic interaction mechanisms is illustrated in Figure 1, following the consideration of this problem in [24].In [24] the last two graphs, (N) and (CN) were combined and treated as a Coulombmodified strong-interaction amplitude.The modification was approximated by giving an extra phase factor to the hadronic amplitude.In the literature, this factor is called the Coulomb phase [24][25][26].Alternatively, one can combine (C) and (CN) and obtain a hadronic correction to the Coulomb amplitude.In what follows, we call this absorptive correction [27,28].
As usual, multiple interaction amplitudes, depicted in Figure 1, are easily calculated in impact parameter representation, where the result is just a product of multiple amplitudes.Thus, we switch from q T to b dependent amplitudes, and simultaneously, from the Born approximation to the eikonal optical model.
These hadronic eikonal phases were used in the Glauber model expressions ( 12) and (13).
Adding higher order terms to the Born Coulomb amplitudes ( 22) and ( 23), as is explained in [28], we arrive at the eikonal form for the electromagnetic amplitudes with Coulomb eikonal phases, given by the Born amplitudes in impact parameter representation.
Combining the Coulomb (C) with Coulomb-nuclear (CN) mechanisms depicted in Figure 1, one obtains the Coulomb terms, Equations ( 24) or (25), which acquire an absorption factor, given by the standard Glauber eikonal approximation [29], As usual, the complicated multi-loop integrations in momentum representation are essentially simplified to a multiplicative combination in impact parameters.This is why the correction S(b) of Equation ( 28) enters the final expressions as a factor.Thus, the absorptioncorrected amplitudes f take the form Eventually, we arrive at the full amplitudes.
First of all, we performed calculations within the Born approximation, Equations ( 9) and ( 10) and ( 14) and (15).The hadronic spin-flip component was set to zero.The results are depicted with black dotted curves in Figures 2 and 5.The magnitude of A N substantially exceeds the data.
Then, we relied on the eikonal form of higher order terms, Equations ( 24) and ( 25), keeping r 5 = 0.The results are depicted with blue double-dot-dashed curves.The effect of eikonalization turns out to be rather mild, and the discrepancy with the data remains significant.
The next step is the introduction of absorptive corrections, which significantly reduce the values of A N (t) as is demonstrated by the red solid curves, calculated with Equations ( 29) and (30) (r 5 is still zero).
Eventually, we can adjust the single unknown parameter, r pA 5 , for each nuclear target.The results of the fit are presented in Table 1.
The found values of r pA 5 are pretty close to the values for the Pomeron spin component found by fitting to the pp data in [7].
With values of r 5 in Table 1 and amplitude Equations ( 29)-( 30), we plot the green dashed curves, which we treat as our final results.The same as in Figure 2, but for proton-gold elastic scattering.Data at 100 GeV are from [31,32].Notice that the nuclear data are quite sensitive to the value of r pA 5 ; this is why the CNI method was proposed [6] as a unique way for measuring the hadronic spin-flip amplitudes at high energies.

Conclusions
Concluding, we performed the first calculation of single-spin asymmetry in polarized proton-nucleus elastic scattering in the CNI region.We achieved a reasonable agreement with the data, in spite of the rather complicated theoretical construction.The remarkable feature of the nuclear form factor ( 12) is a change in the sign of the imaginary part of the elastic hadronic non-flip amplitude corresponding to the first zero of the Bessel function J 1 (t).Since the non-flip hadronic amplitude significantly exceeds the spin-flip part, they become of the same order right before and after the former changes sign.The asymmetry A N (t) Equation ( 2) reaches maximum at f n f = f s f , so A N (t) should abruptly vary between positive and negative maxima in the vicinity of the Bessel zero.Comparison with the data in Figures 4 and 5 confirms such behavior; indeed, the measured A N (t) exposes two maxima with opposite signs with positions close to predicted.
However, the magnitude of these maxima is exaggerated in our calculations, so there is still room for improvements.Here, nuclear effects were evaluated within the Glauber approximation, which is subject to Gribov inelastic shadowing corrections [33,34].Their calculation requires substantial modeling, including knowledge of the proton wave function, interaction mechanism, etc.This needs a detailed study, as was performed in [23].We leave this issue for future development.

Figure 1 .
Figure 1.Three types of interaction: pure electromagnetic (a), pure strong interaction (b), and combined strong and electromagnetic interactions (c). .

Figure 3 .
Figure3.The same as in Figure2, for proton-aluminum elastic scattering.Data at 100 GeV are from[31].