2 2 A pr 2 01 9 New feature in the di ff erential cross sections at 13 TeV be measured at the LHC

Analysis of dσ/dt of the TOTEM Collaboration data, carried out without model assumptions, showed the existence of a new effect in the behavior of the hadron scattering amplitude at a small momentum transfer at a high confidence level. The quantitatively description of the data in the framework of the HEGS model supports the existence of such phenomenon which can be connect with quark potentials at large distances.


Introduction
Research into the structure of the elastic hadron scattering amplitude at superhigh energies and small momentum transfert can give a connection between experimental knowledge and the basic asymptotic theorems, which are based on first principles [1,2]. It gives information about hadron interaction at large distances where the perturbative QCD does not work, and a new theory, as, for example, instantons or string theories, must be developed. There were many works in which the consequences of breaking the Pomeranchuk theorem were investigated [3]. It was shown [4] that if the Pomeranchuk theorem was broken and the scattering amplitude grew to a maximal possible extent, but not break the Froissart boundary, many zeros in the scattering amplitude should be available in the nearest domain of t → 0 in the limit s → ∞. Hence, with increasing energy of colliding beams, some new effects [5] in differential cross sections can be discovered at small t [6]. Such posibility were explored in some works [7,8]. They had represent an explanation of some separate ( individual) effects in the definite domain of momentum transfer and with a period proportional to ∆t. In the previous paper [9], it was shown that "AKM oscillations" could exist in high-precision experimental data of the UA4/4 Collaboration. Now we examined this effect assuming the existence of the potential of hadron-hadron interactions at large distance and made a new quantitative treatment of experimental data at 13 TeV at a high confidence level.

Hadron potential at large distance
The standard fast falling potentials of the Gaussian type lead to the exponential dependence of the scattering amplitude in the range of small momentum transfer. During a long time there have been different attempts to find unusual behavior of the amplitude of the elastic hadron scattering. In [10], it was shown that peripheral contributions of the inelastic diffraction processes lead to the appearance in the elastic cross sections of large and small periodical structures over t. Note also, that in [11], was shown that the potential of the rigid string led to oscillations of the pion elastic form factors at large distances (q = 10 − 20 f m −1 ). So let us take an additional potential in our case in the impact parameter space in the following form: here ω is the frequency and ϕ(s) is the phase of the additional potential. This potential leads to the scattering amplitude in the form . In the case of cutting the potential, for example by the θ-function, we shall get a strongly oscillating expression as t → 0 with the period, which is approximately proportional t, and with the amplitude, which is suppressed as t grows.
In the case of cutting the potential, for example by the θ-function, we shall get a strongly oscillating expression as t → 0 with the period, which is approximately proportional t, and with the amplitude, which is suppressed as t grows. As, example, note the potentials with peculiarities on the determined distance, considered in [12], for example lead to similar behavior of the additional term. In contrast with the potential of the multi-gluon exchange there is no necessity to enter the additional cut of the potential. But on the other hand, generally speaking, we have to get eikonalization of the Born terms. At first, we need to calculate the eikonal and then the amplitude in the t-representation together with the leading eikonal. Numerical integration allows one to calculate an additional term in both the direct and eikonal approach. The results of calculations are shown in Fig.1.

Analysis of the TOTEM 13 TeV data by the selection method
The usual method of minimization χ 2 in this situation often works poorly. On the one hand, we should define a certain model for part of the scattering amplitude having zeros in the domain of small t. However, this model may slightly differ from a real physical picture. On the other hand, the effect is rather small and gives an insignificant change in the sum of χ 2 . Therefore, in this work let us apply another method in first, namely, the method of comparison of two statistically independent choices, for example [13]. If we have two statistically independent choices x ′ n 1 and x " n 2 of values of the quantity X distributed around a definite value of A with the standard error equal to 1, we can try to find the difference between x ′ n 1 and x " n 2 . For that we can compare the arithmetic mean of these choices The standard deviation for that case will be δ x = [1/n 1 + 1/n 2 ] 1/2 And if ∆X/δ x is larger than 3, we can say that the difference between these two choices has the 99% probability .
The deviations ∆R i of experimental data from these theoretical cross sections we will be measured in units of experimental error for an appropriate point where δ exp i is an experimental error. To take this effect into account, we break the whole studied interval of momentum transfer into k equal pieces kδt = (t max − t min ), where sum ∆R i separately over even and odd pieces. Thus, we receive two sums S up and S dn for n 1 even and n 2 odd interval. At this n 1 + n 2 = k and |n 1 − n 2 | = 0 or 1 In the case of some difference of experimental data from the theoretical behavior expected by us or incorrectly determined parameters, these two sums will deviate from zero; however, their sizes should coincide within experimental errors. However, this will be so in the case if experimental data have no any periodic structure or a sharp effect coincides with one interval.
We assume that such a periodic structure is available and its period coincides with the chosen interval 2δt. In this case, the sum S up will contain, say, all positive half-cycles; and the sum S dn , all negative half-cycles. The difference between S up and S dn will show the magnitude of an additional effect summed over the whole researched domain. Our method does not require exact representation of the oscillatory part of the scattering amplitude, and now let us apply it to new LHC data of the TOTEM Collaboration at 13 TeV [14,15]. It gives us two sets of the data: one 0.000879 < |t| < 0.201041 GeV 2 includes the Coulomb-hadron interference range and the other 0.0384 < |t| < 3.822 GeV 2 for the large t. Both sets have the overlapping region. We find that three first points of the second set hardly differ from the data of the first set and we removed them. For the first analysis we do not include the region of the diffraction minimum as we try, using this new method, to examine some additional oscillation behaviour with minimum model assumption, but the region of the dip requires some model for the description of the elastic scattering in a wide region of t, Here and below we use only statistical errors, and systematic errors are taken into account as an additional normalization coefficient for both the sets. For the basic (non-oscillating ) amplitude we use the standard exponential form with three slopes multiplied by t, t 2 , t 3 . In Fig.2, such sums in eq.(3) are represented for the period that is proportional to √ −t. Then we move the segments by one-half the segment and calculate these sums again. The results are shown in Fig.2 by the long-dashed and dotted lines. The large difference between the first and second cases clearly shows the existence of some oscillation contributions in the scattering amplitude especially at large t. Now let us calculate such sums with the period that is proportional t. The results are shown in Fig.3 Obviously, it is show the existence of the oscillatory contributions on high confidence level.

HEGS model analysis
Now let us try to find the form of such additional oscillation contribution to the basic elastic scattering amplitude. As a basis, take our high energy generalized structure (HEGS) model [16,17] which quantitatively describes, with only a few parameters, the differential cross section of pp and pp from √ s = 9 GeV up to 13 TeV, includes the Coulomb-hadron interference region and the high-|t| region up to |t| = 15 GeV 2 and quantitatively well describes the energy dependence of the form of the diffraction minimum [18]. However, to avoid possible problems connected with the low-energy region, we consider here only the asymptotic variant of the model [19]. The total elastic amplitude in general receives five helicity contributions, but at high energy it is enough to write it as F(s, t) = F h (s, t) + F em (s, t)e ϕ(s,t) , where F h (s, t) comes from the strong interactions, F em (s, t) from the electromagnetic interactions and ϕ(s, t) is the interference phase factor between the electromagnetic and strong interactions [20,21,22,23]. The Born term of the elastic hadron amplitude at large energy can be written as a sum of two pomeron and odderon contributions, All terms are supposed to have the same intercept α 0 = 1+ǫ 0 = 1.11, and the pomeron slope is fixed at α ′ = 0.24 GeV −2 . The model takes into account the two hadron form factors F 1 (t) and A(t), which correspond to the charge and matter distributions [24]. Both form factors are calculated as the first and second moments of the same Generalized Parton Distributions (GPDs). As a probe for the oscillatory function take here J 1 (τ) is the Bessel function of the first order. This form has only a few additional fitting parameters and allows one to represent a wide range of possible oscillation functions.
After the fitting procedure, with the modern version of FUMILY [26], we obtain χ 2 /do f = 1.24 (remember that we used only statistical errors). One should note that the last points of the second set above −t = 2.8 GeV 2 show an essentially different slope, and we removed them. The total number of experimental points of both sets equals 415. If we remove the oscillatory function, then χ 2 /do f = 2.7, so an increase is more than two times. If we make a new fit without f osc , then χ 2 /do f = 2.5 decreases but remains large. Our model calculations are represented in Fig.4(a,b,c). It can be seen that the model gives a beautiful description of the differential cross section in both the Coulomb-hadron interference region and the region of the diffraction dip.
To see the oscillations in the differential cross sections, let us determine two values -one is pure theoretical and the other with experimental data The corresponding values calculated from the fit of two sets of the TOTEM data at 13 TeV are presented in Fig.5. At small t there is a large noise; however, the oscillation contributions can be seen. This corresponds to the small size of the S uS d values in Fig.3; however, at large −t > 0.1 GeV 2 we can see that R∆ th is similar to the value R∆ Exp . The oscillation contribution is small; however, the noise of the background decreases at this t and does not dump the oscillation part.

Conclusion
We show that HEGS model describe on quantitatively level the new experimental data at 13 TeV with taking into account only statistical errors. However, only adding the small oscillatory term allow obtain χ 2 do f = 1.25. The phenomenon of oscillations of the elastic scattering amplitude will give us the important information about the behavior of the hadron interaction potential at large distances. We have shown the existence of such oscillations at the statistical level by three methods: a) the method of statistically independent selection; b) the comparison of the χ 2 without oscillation ( χ 2 = 1140) and with oscillation ( χ 2 = 515); c) the comparison of R∆ th and R∆ exp , Fig.5). All three methods show the presence of the oscilatory behavior. Very likely that such effects exist also in experimental data at essentially smaller energies [27] but, maybe, they have a more complicated form (with two different periods, for example). The phenomena of oscillations are also related to the asymptotic properties of the scattering amplitude. They can impact the determination of the sizes of the total cross sections, the ratio of the elastic to the total cross sections and the size of the ρ(s, t) -the ratio of the real to imaginary part of the elastic scattering amplitude.