pp Elastic Scattering at LHC Energies

Using a unified analytic representation for the elastic scattering amplitudes of pp scattering valid for all energy region, the behavior of observables in the LHC collisions in the range $\sqrt{s}$= 2.76 - 14 TeV is discussed. Similarly to the case of 7 TeV data, the proposed amplitudes give excellent description of the preliminary 8 TeV data. We discuss the expected energy dependence of the observable quantities, and present predictions for the experiments at 2.76, 13 and 14 TeV.


I. INTRODUCTION
Elastic scattering is described by one single complex function depending on two kinetic variables: the incident center of mass energy √ s and momentum transfer q. In high energy pp(p) scattering, the scattering amplitude is usually represented as T (s, t), where t is the four momentum transfer squared. More than a decade ago, Ferreira and Pereira analyzed all available elastic scattering data for energies above 20 GeV [1] and all |t|, identifying properties of the amplitudes (zeros, signs, magnitudes), with proper attention given to the real part, which plays a critical role in differential cross sections for mid and large |t| ranges.
Recently, this analysis was extended [2] to the LHC-TOTEM elastic scattering 7 TeV data [3], and also the behavior of proposed amplitudes was re-examined in the whole energy region from 20 GeV to 14 TeV to determine the precise energy dependence of the model parameters [4], and applied to the cosmic energy domain with calculation of p-air cross sections [5]. From this analysis, an analytic representation of scattering amplitudes as function of √ s and t was established. In the present work we apply these analytical forms to investigate in detail the LHC energy region from 2.76 to 14 TeV.
We stress that we establish explicitly disentangled real and imaginary amplitudes based on a QCD motivated model, and not just fit pure phenomenological expressions to observables. Besides, since the so-called impact parameter representation (s, b) and its Fourier transform in (s, q) space are both represented by simple analytical forms, we are able to control unitarity and dispersion relation constraints, and provide geometric interpretation of the interaction range. The regularity that we obtain in our treatment of the data and associated reasonable physical interpretation of the consequences give reliability to our proposal of disentanglement of the amplitudes [4].
The present work is organized as follows. In the next section, we describe briefly the amplitudes and their energy dependences in both t-and b-representations. In Sec. III we apply these amplitude to describe observables and discuss their energy dependence, and also we investigate consequences for very high energies of the form of the amplitudes in b-space. In Sec. IV we make use of the properties of our amplitudes and observables for the LHC range, particularly for √ s = 8 TeV, where preliminary information on dσ/dt starts to become available.
The last section is devoted to discuss further our results and perspectives, together with geometric interpretation in the b-space representation.

II. ANALYTIC REPRESENTATION OF THE AMPLITUDES
A.

Impact parameter Representation
The Fourier transform of the momentum transfer q amplitudes to the b-space defines the impact parameter (or simply b-space) representation. Since the impact parameter variable b is not observable, the treatments of data are made usually in (s, t) space, except for integrated cross sections. However the b-space description gives insight in geometric aspects of the collision, since in the classical limit the variable b reduces to the physical impact parameter. Besides, it plays important role in the eikonal representation, where unitarity constraints are more simply formulated. On the other hand, the dispersion relation (causality) constraint is properly dealt in t-space. In the following discussion, we do not consider effects of spin or polarization.
The amplitudes for the description of pp scattering in the Stochastic Vacuum Model [6] are originally constructed through prof ile functions in b-space and here we describe this formalism first. The dimensionless (s, b) amplitudes due to the nuclear interaction are written as with the characteristic shape function The label K = R, I indicates either the real or the imaginary part of the complex amplitude. The fixed quantity a 0 = 1.39 GeV −2 is related to the square of the correlation length a of the correlation function of the gluon condensate, with a = (0.2 ∼ 0.3) fm , as measured in hadronic interactions and in lattice QCD, with our best choice 0.27 fm. In the large b behaviour of the profile function of the Stochastic Vacuum Model there appears the dimensionless combination b 2 /a 0 where a 0 = [a/(3π/8)] 2 , that fixes the value of a 0 appearing in Eq. (2). The quantity 3π/8 is a feature of the correlation function [6] .
The Gaussian form of the first term in Eq. (1) is similar to the usual formalism of reggeon exchanges [7]. The second term, referred to as shape function, represents contributions from the perturbed vacuum structure around the protons at larger b values. It is zero at b = 0 and is normalized as In Eq. (1) we have introduced four energy dependent parameters for each amplitude, α K , β K , γ K , λ K , with γ K dimensionless, while α K , γ K and β K are like GeV −2 .
In the small and mid b ranges there is superposition of the contributions of the two parts, that, for convenience of language, we may call respectively Regge phenomenology and loop-loop interaction. The resulting parameter values are determined describing with accuracy the imaginary and real amplitudes as a whole, there is no case of double counting effects, and each part is duly represented, if one thinks of each one separately. Actually Eq.
(1) represents an extension of parametrization of results of the Stochastic Vaccum Model, opening possibilities of introducing proper s and t dependences.
Although b is not exactly the physical impact parameter, neither observable, the b-space representation permits a geometrical interpretation of the behavior of the amplitude. For large b, which corresponds to peripheral collisions, the amplitudes fall down with a Yukawa-like tail, that reflects the effects of virtual partons (the modified gluon field) at large distance in the Stochastic Vacuum Model. A feature of the b-space representation is that it can be directly related the eikonal formalism, as shown below. We introduce the eikonal function χ (s, b) through with Separating real and imaginary parts, we have sin From Eq. (8) we have immediately and thus the general unitarity constraint is witten as Our solution, at all energies, satisfy this bound condition. Satisfying a monotonic behavior of the scattering amplitudes, our solutions are restricted to the branch where χ R ≥ 0, and thus in turn, we have Under these conditions, our analysis shows that for a fixed √ s, the function T I (s, b) is monotonically decreasing in b. The maximum of the imaginary amplitude, T I (s, b = 0) tends to its limiting value √ π for asymptotic large energies [4].
In terms of the T K (s, b) amplitudes, the elastic, total and inelastic cross sections are written respectively and In terms of the eikonal function, we write d σ el (s, b) B.

t-space representation
The comparison with dσ/dt data and determination of parameters are made with the amplitudes in t-space. The quantities Ψ K (γ K (s), t = − q 2 T ) obtained by Fourier transform of Eq. (1) are written (18) with K = R, I, and the shape functions in t− space take the form = 2 e γK e −γK √ 1+a0|t| 1 + a 0 |t| − e γK e −γK √ 4+a0|t| 4 + a 0 |t| , with the property that corresponds to Eq. (3). The expression (18) represents the nuclear amplitude due to the non-perturbative QCD interactions that dominate the low and mid |t| regions. To describe elastic dσ/dt data for all |t|, we should account for contributions from perturbative processes. We thus add a term representing the perturbative three-gluon exchange amplitude [8] that may appear in the large |t| region, and the complete nuclear amplitudes are then written where the Kronecker delta symbol δ K,R is introduced since we define R ggg (t) as the real contribution from the perturbative three-gluon exchange amplitude. The effect of the tail term R ggg (t), producing a universal (not energy dependent) |t| −8 form for large |t| in dσ/dt, was studied in the analysis of the experiments at CERN-ISR, CERN-SPS [1], 1.8 TeV [9] and 7 TeV [2]. We write where the last two factors cut-off this term smoothly in the non-perturbative domain, and the signs ± refer to the pp and pp amplitudes respectively. Although the cut-off factors written in Eq. (22) have been adequate for all cases that were examined, their detailed forms in the transition range (2.5 < |t| < 4) GeV 2 must be examined with data. For a complete analysis of elastic scattering, we must also take into account the contribution from the Coulomb interaction.
The complete amplitudes T R (s, t) and T I (s, t), with dimensions GeV −2 , contain the nuclear and the Coulomb parts as and where α is the fine-structure constant, Φ(s, t) is the Coulomb phase and F C (t) is related with the proton form factor for the pp/pp collisions. The proton form factor is taken as where t 0 = 0.71 GeV 2 . Note that the strong interaction part of the amplitudes are smooth and regular functions of s and t, while the Coulomb amplitude is relevant in the very forward range |t| < 10 −2 GeV 2 . In our normalization the elastic differential cross section is written and the total pp cross section is given by the optical theorem The analysis of all pp elastic scattering data for √ s from 20 GeV to 7 TeV leads to a separate identification of the real and imaginary parts contributing to Eq. (27). The energy dependence of the eight parameters is given below, with √ s in TeV, and GeV −2 in the units of the parameters that are not dimensionless ( γ I and γ R are dimensionless).
α I (s) = 11.0935 + 1.35479 log √ s, λ I (s) = 14.02008 + 3.23842 log √ s + 0.444594 log 2 √ s , The peculiar (not so simple) expression for β I (s) is constructed in order to satisfy both the low-energy phenomenology and unitarity constraints at all energies, as given in Eq. (11), and leads to the asymptotic behavior For very high energy and considerations of asymptotic behaviour, it is useful to use the simpler form for β I (s) The first term in Eq. (18) can be written in the usual notation of Regge phenomenology, with the dimensionless scattering amplitude A(s, t) where |t| is in GeV 2 and √ s in TeV, with a t-dependent residue and a trajectory with intercept 1 and angular coefficient 0.19 GeV −2 . The log term corresponds to a double pole, arising from derivative with respect to the trajectory [7]. These expressions are able to give high precision representation for all data [1,2,4,9], with coherent and reliable identification of the real and imaginary amplitudes. Properties and consequences for the energy range above 1 TeV are discussed in the present paper, with particular attention to the experimental LHC energies. Cosmic ray energies up to √ s =100 TeV and asymptotic behaviour have been discussed elsewhere [5].

C. Forward Amplitudes and Associated Observables
In the very forward direction, where the elastic pp and pp scattering amplitudes can be approximated by pure exponential forms, the differential cross section is written where t ≡ −|t| and we must allow different values for the slopes B I and B R of the imaginary and real amplitudes. With σ in milibarns and |t| in GeV 2 , we have (hc) 2 = 0.3894. Since we work with B R = B I , treatment of the Coulomb interference requires a more general expression for the Coulomb phase, which has been developed before [2]. The limits of the amplitudes for small |t| give the total cross section σ, the ratio ρ of the real to imaginary amplitudes, and the slopes B R,I at t = 0 through where √ s is in TeV, σ in milibarns, B I and B R are in GeV −2 ; ρ is dimensionless, passes through a maximum at about 1.8 TeV, and decreases at higher energies, with asymptotic value zero. The ratio B R /B I is always larger than one, as expected from dispersion relations [10]. The ratio B R /B I as function of the energy is shown in Fig. (1). There is a finite asymptotic value B R /B I → 1.887.
This treatment of pp forward scattering has been applied to the calculation of p-air cross sections measured in Extend Air Showers studies in cosmic ray experiments. Covering the range from 1 to 100 TeV in pp energies, our input amplitudes are used as basis of Glauber calculations, giving good description [5] of all cosmic ray data. . At all energies it is BR > BI , as predicted by dispersion relations [10]. In the figure, the ratio BR/BI is plotted as function of the energy, indicating the finite asymptotic limit.

A. Differential Cross Sections and Amplitudes
In Fig. 2 we show the predictions for dσ/dt for the LHC energies 2.76 , 8 , 13 and 14 TeV. We first observe that the dip and the bump peak displace to the left as the energy increases and in this figure these displacements follow almost straight lines, as indicated by marks with black circles and open squares. For the sake of convenience, we list the values of parameters for these energies in Table I, where γ I , λ I , α R and γ R are substituted by more commonly used quantities σ, ρ together with the slope parameters B I and B R . In Table II we show the values of several quantities obtained in the numerical calculation of the amplitudes and of observables in the elastic process. Some characteristic features are exhibited below in plots.
In Fig. 3 we use the energy √ s = 8 TeV as an example to show the imaginary and real amplitudes T N I (s, t), T N R (s, t) as functions of |t| as predicted by Eq. (22). For all energies the characteristic features are the two zeros of the real part, and the single zero of the imaginary part appearing in the plotted range (a second zero of T N I would appear in a much larger |t|, outside experimental visibility). The interplay of the imaginary and real amplitudes at mid values of |t| is responsible for the dip-bump structure of the differential cross section, that was shown before [2] for √ s = 7 TeV, and is exemplified for 8 TeV in the next section. For |t| ≥ 1.5 GeV 2 the real part becomes dominant, with positive sign. The inset shows the small |t| range, in log scale, normalized to one at |t| = 0. The straight exponential slopes are shown in dashed lines, with the dramatic difference between the real and imaginary amplitudes. Soon the exact amplitudes leave the straight line and curve down, searching for their respective zeros. As shown in the next section, the consequences for the behavior of dσ/dt at 8 TeV will be visible for |t| larger than about 0.2 GeV 2 .
The difference in slopes B R and B I that is required by dispersion relations [10], is often neglected. The real part is small for small |t|, due to the small value of ρ, but becomes influential or dominant for mid and large |t|. The amplitudes must be treated as functions for the whole |t| range. Our unique analytical form connects all regions and controls the behavior both at small and large |t|. Thus, for example, the value of ρ is very important for the shape of the dip-bump structure.
The regular energy dependence of the positions of the zeros and of dips and peaks of bumps is shown in Fig. 4. We see that all these characteristic quantities move towards smaller |t| with increasing energy, following forms like The general behaviour is the same for all energies, with one and two zeros respectively for the imaginary and real parts. The behaviour for small |t| is shown in the inset, indicating the difference of slopes BR and BI at the origin, and the deviations of the exponential forms that occur as |t| increases, each amplitude going towards its zero. A second zero of the imaginary part occurs at much higher |t|.  (2) |t| dip dσ/dt| dip |t| peak dσ/dt| peak ratio |t| infl dσ/dt| infl σ inel σ el σ I el σ R el σ el /σ TeV GeV 2 GeV 2 GeV 2 GeV 2 mb/GeV 2 GeV 2 mb/GeV 2 R GeV 2 mb/GeV 2 mb mb mb mb 1.
R , that at very high energies behaves as above, with A = 0 and c = 0, according to a theorem by A. Martin [11]. This behaviour is obviously connected with a fast increase of the slope B R .
It is interesting to observe the relative positions of the dip and the peak of the bump in dσ/dt and the zeros of the imaginary and real parts, shown in Fig. 4. This question has been discussed a long time ago [1]. The figure shows that Z I and the dip position tend to the (apparently) common finite limit. Dips and peaks are always located between Z I and Z (2) R . All energy dependences are simple and can be easily parameterized.
It is interesting to note that the ratio between the maximum of the mid-|t| bump (called peak) and the dip minimum increases with energy rather rapidly (see Fig. 4-b), like ∼ ln 2 √ s, while the distance |t| peak − |t| dip remains practically constant (Fig. 4-a).
In Fig. 5 we plot dσ/dt for 2.76 and 8 TeV, showing that the characteristic dip/bump structure of dσ/dt occurs in the interval between the imaginary zero and the second real zero.

B. Integrated Quantities, Ratios and Asymptotic Limits
The integrated elastic cross section due to the imaginary amplitude can be represented by with √ s in TeV and σ I el (s) in mb. The accuracy of this representation is very good, particularly for energies equal and above 7 TeV . The ratio with the total cross section has a finite asymptotic limit at high energies σ I el /σ → 0.354. This result is very important for a geometrical description of pp scattering, as it means that pp collision does not follows a black disk form at high energies (see below).
For the contribution of the real part to the elastic cross section the quantity that is related to the exponential behaviour in the forward direction and that presents a finite asymptotic ratio with σ requires an extra factor 1/ρ 2 . We have the representation = 10.2037 + 2.47691 log √ s + 0.23108 log 2 √ s .
The dimensionless ratios are related to σ I el /σ and (1/ρ 2 )σ R el /σ when the amplitudes are of pure exponential forms with B I and B R slopes. The imaginary part is studied to investigate the occurrence of black disk behaviour (assuming zero real part), where the ratios σ I el /σ and σ/(16πB I ) are both equal to 1/2. As shown in Fig. 6 our solutions lead to values about 1/3 for the imaginary part case, which is a more realistic expectation [12] than the black disk hypothesis.

C. Geometric Scaling and Ratio of Cross Sections
In Fig. 7-a, we plot d 2 σ inel /d b 2 defined in Eq. (14) as function of b for √ s = 2.76 and 14 TeV. The behavior at very high energies ( √ s = 10 4 , 10 5 and 10 6 TeV) is also shown. We clearly see the increase of effective radius of the interaction range with increasing energy. In Fig. 7-b, we plot the same quantities with respect to the variable This figure shows clearly that there exist a universal function ξ (x) such that for √ s ≫ 10 4 TeV. An important point is that ξ (x) is far from the Heaviside step function, rather possesing a considerably diffused surface. In this asymptotic limit, we can safely set cos χ R → 1 so that the total cross section is [4] for all x where the equality holds if and only if ξ = 0 or ξ = 1. Therefore, whenever the function ξ is different from a sharp-cut Heaviside theta For our amplitudes, as shown in Fig. 7, ξ clearly does not converge to a sharp-cut θ function, preserving an appreciable diffused surface for asymptotic energies. Therefore, we have or This means that our amplitudes do not show the black disk behavior at very large energies, deviating from the well known result for a black disk σ el /σ (s) → 1/2. From the above discussion, we can also easily see that the more diffused surface ξ (x) has, the less the ratio σ el (s) /σ (s) becomes. In fact, for our case this ratio is close to 1/3 (see also [12]). Note that this is somewhat different scenario compared to [13], where ξ (x) would not have surface diffuseness.  4. a) Positions of the zeros of the amplitudes, and of the dip and peak at the bump of dσ/dt. There appears one zero in the imaginary and two in the real amplitude. A second imaginary zero occurring at very large |t| is outside the physically accessible range. All quantities move towards small values with increasing energies. The dips tend to coincide with the imaginary zero at high energies. The remarkable dip/bump structure in pp scattering occurs in the interval between the imaginary zero and the second real zero. The first real and the imaginary zero move towards smaller |t|, indicating the log 2 √ s increase of the real and imaginary slopes. The dots are put to help the connection of values of the quantities for different energies. b) There is a regular and fast increase of the ratio R = [dσ/dt]peak/[dσ/dt]dip, with increasing sharpness of the dip/bump structure although the distance |t| peak − |t| dip between them varies very little. These symptoms come from the increasing proximity of |t| dip and ZI , and to the convergence to finite asymptotic limits of both |t| peak and |t| dip . R . This behavior is shown in this figure for the energies 2.76 and 8 TeV. As the energy increases |t| dip approaches ZI from the right to the left. Fig. 4 illustrates these properties again, in another way.

IV. COMPARISON WITH DATA AND PREDICTIONS
Our description [2] of the elastic scattering data at 7 TeV from the TOTEM Collaboration [3] reproduces N=165 points in dσ/dt with an impressive squared average relative deviation < χ 2 >= 0.31. Characteristic quantities at this energy, shown in Tables I and II  After the successful description of the 7 TeV data [2], we now present comparison and predictions for other LHC energies.

A. Inelastic and Total Cross Sections
For the inelastic cross section we assume the difference σ inel = σ − σ el and then we have 73.26 mb at 7 TeV.. Published values of the TOTEM Coll. using different methods are 73.15 ± 1.26 [3], 73.7 ± 3.4 [14] and 72.9 ± 1.5 [15]. ALICE Coll. [16] gives σ inel = 73.2 ± 5.3 mb , and ATLAS Coll. σ inel = 69.4 ± 2.4 ± 6.9 mb [17]. We are not able to understand the CMS results [18] in terms of pure σ inel due to non-informed missing contributions. In these measurements there are extrapolations to using Monte Carlo models to include diffractive events of low mass. Of course all these results are compatible with our calculations.
A measurement to be compared with our predictions  is the √ s = 2.76 TeV value of ALICE Coll., that gives σ inel = 62.8 ± 4.2 mb , while our tables give the compatible value 63.11 mb.
The analysis of compatibility for the 1.8 TeV measurements of σ inel by CDF and E811 in Fermilab [19] suggests the value (1 + ρ 2 )σ inel = (60.3 ± 2.3 mb, that with with our ρ value gives σ inel = (59.1 ± 2.3 mb. Our table gives 58.89 mb for 1.8 TeV, once more in very good agreement. Finally, at 57 TeV the Auger Cosmic Ray experiment [20], using other models for the pp input, evaluates σ inel = 92 ± 14.8 mb , while our extrapolation gives 101 mb. We have discussed this measurement [21] together with other CR Extended Air Showers (EAS) experiments, using our amplitudes as inputs and a basic Glauber method to connect pp and p-air processes. Our calculation reproduces well all CR data for p-air cross sections with √ s (in the pp system) up to 100 TeV. For 8 TeV we have predictions σ = 101.00 mb , σ el = 26.18 mb , σ inel = 74.82 mb , σ el /σ = 0.26 shown in the tables. The measurements by TOTEM [22] give for the same quantities σ = 101.7 ± 2.9 mb , σ el = 27.1 ± 1.4 mb, σ inel = 74.7 ± 1.7 mb , σ el /σ = 0.266 ± 0.006. Of course these numbers are very encouraging, indicating also good expectations for dσ/dt at this energy.
The data and our curve for σ inel (s) are shown in Fig.  8. All this information shows that our formulae for the energy dependence of σ(s) and σ inel (s) in pp scattering work very well.

B. Expected data for dσ/dt at 8 TeV
The preliminary data for dσ/dt at 8 TeV, shown in talks by members of the TOTEM Collaboration [23], are encouraging for the application of our method of analysis. We recall that in the treatment of the 7 TeV data, we obtained precise description, with average χ 2 = 0.34 for 165 data points in the whole |t| interval of measurements. If Fig. 9 we shown our calculation for dσ/dt covering the whole |t| range of the preliminary information, using the amplitudes defined in Sec. III. The characteristic features of the forward peak and of the dip/bump structure are expected to represent accurately the angular dependence. Numerical values for characteristic features are given in Tables I and II. This is the description of the global dσ/dt data at 8 TeV, that promises to be more complete and regular than the 7 TeV data, except for not reaching larger |t| values. In the following we discuss the forward region in more detail.
In Fig. 10 we plot the calculations in the small |t| range , including the influence of the Coulomb phase [2]. The calculation with Coulomb phase put equal to zero is represented by the dashed line, showing that its influence is small. Our specific calculation of the Coulomb phase takes into account the difference in values of the B R and B I slopes. Other calculations for the interference phase [24] also show that its influence is small, reducing dσ/dt by a few percent.
Our values for B I and B R given in Table I lead to the dσ/dt effective slope at 8 TeV Our predictions seem to be in accordance with the eyeguided reading of the preliminary data of dσ/dt that appear in presentations of the TOTEM group in workshops, at least at the qualitative level. At 7 TeV our expressions perform extremely well when compared to the published experimental information, and we expect that the same will happen at 8, 13 and 14 TeV .

C. Other models
The complete identification of the real and imaginary parts of the complex pp elastic amplitude is fundamental for the knowledge of the dynamics of the collision, being an essential bridge between the observed quantities and fundamental QCD dynamical processes. Our determination, though consistent and complete, depends on the analytical forms used for the representation. It is thus important to compare our predictions with the results obtained with other input assumptions.
The important Yukawa-like behaviour of the amplitudes in b-space based on the behaviour of the loop-loop interaction for large b in the Stochastic Vacuum Model, that is incorporated in our input amplitudes in Eq. (2), is confirmed in a recent treatment of the pp interaction through Wilson correlation functions [25]. The representation of amplitudes in b-space from the ISR to the SPS energies shows at b = 0 a slow increase with the energy [26], remaining below saturation, which seems to be approached asymptotically, as can be seen in the present work for the LHC energies and also in studies at higher cosmic ray energies [5]. With parameters adjusted to describe the energy dependence in the 23 GeV -546 GeV range [27] the model predicts characteristic quantities of pp forward scattering for the range 7 -14 TeV, with results for the total cross section and the slope parameter that agree very well with the numbers given in Table I, and in particular with the TOTEM values at 7 TeV.
The model proposed by Bourrely, Soffer and Wu (hereafter called BSW model) [28] gives explicitly the full s,t dependence of the elastic scattering amplitudes and is appropriate for the comparison with our results. Important similarities and differences were discussed in detail in the 7 TeV case [2], and we now compare the predictions for 14 TeV. Fig. 11 shows that the dip-bump structure occurs in similar |t| regions, but there is a difference in dσ/dt by a factor larger than 2. This difference results from the larger magnitudes (with negative signs) of both real and imaginary parts in the BSW model, as can be observed in the second part of the figure. The second real zero occurs for a larger |t| in the BSW calculation.
To raise interest on measurements at higher |t|, in Fig.  12 the 14 TeV plot is extended to very large |t| pointing out the possible smooth connection with the supposedly universal tail at 27.4 GeV [29]. We recall the situation with similar plot drawn in the 7 TeV case [2], where the measurements reached larger |t| ≈ 2.5 GeV 2 and the conjecture of the universality of the tail at such high energy encounters motivation. The broad dip in the region of 6 GeV 2 in the BSW calculation is due to a zero in its imaginary amplitude, as was also indicated at 7 TeV .
The structure of the pp and pp interactions studied by O. Selyugin [30], based on the analysis of different sets of Parton Distribution Functions and introducing t-dependence in the Generalized Parton distributions, gives good representation of dσ/dt data in large energy range, up to the LHC Totem experiment at 7 TeV. We include in Fig. 11 the t dependences of dσ/dt and amplitudes in this model for 14 TeV. It is important to observe the similarity with our results in the forms of the amplitudes, that differ essentially only in the real part for large |t|, with different locations of the second zero. The similarity, that must be investigated at all energies, reinforces the expectation of the present work, that is to find a realistic and accurate disentanglement of the elastic amplitudes.
Recently, T. Csörgö et al. applied the multiple diffraction calculation of the Glauber-Valesco Model [31] to the TOTEM results for 7 TeV and fitted the parameters of partonic charge form factor, obtaining good representation for the scattering cross section, including the dip-bump structure [32]. Parametrizing the partonic distribution function for each energy, the Glauber-Valesco model also describes the lower energy data from ISR and FNAL. The model leads to the conclusion that at LHC energies the pp interaction is not that of sharp-edged black-discs, but presents a considerable tail in the profile function [33]. These observations are essentially the same as in predictions of our work [2,4,5].
In addition, more detailed measurements of inelastic pp cross section in LHC experiments [34], especially of diffractive-dissociation processes (DD), add very important information on the dynamics of the interaction, as discussed by Lipari and Lusignoli [35]. These studies and ours are complementary in the sense that while our approach is based on the field theoretical model of interaction between two geometrical objects, their approach applies the multiple diffractive model and introduces parametrizations of the partonic form factor of the proton. A detailed comparison of the two approaches, together with the identification of the presence of diffraction in the inelastic part in our formalism, will give more precise insight on the physical roles behind the features of our description. Studies in this direction are under investigation.

V. FINAL REMARKS AND COMMENTS
In this paper we present predictions for observables of pp elastic scattering above √ s = 1 TeV up to coming LHC energies in terms of analytic forms for the real and imaginary parts of the complex scattering amplitude. The representation proposes a separate identification of the two parts, which are both constructed respecting unitarity and dispersion relation constraints, and precisely determine their influences in the observed quantities.
The amplitudes have simple analytical forms, that can be directly evaluated with few operations with elementary functions. The shape of the dip-bump behavior results from a delicate interplay of the imaginary and real amplitudes. All intervening quantities and derived properties are connected by smooth energy dependences.
The zeros of the real and imaginary parts have very regular displacements, converging to finite limits as the energy increases. There is remarkable connection between positions of zeros and positions and heights of dips and bumps and inflections in dσ/dt.
The slopes B I and B R at the origin, with their characteristic difference in values, together with the ratio ρ, are essential quantities that participate in the definition, through the unique analytical forms of the amplitudes, of the properties of the observed dσ/dt in the whole t range. Their values are thus fixed with high accuracy. It is very important that the slopes show quadratic dependence in log s, instead of the linear dependence suggested by Regge phenomenology.
The integrated elastic cross sections are evaluated in their separate parts, obtained from the real and imaginary amplitudes, and are also represented by simple parabolic forms in log s.
The properties of ratios (with respect to the total cross section) of slopes and of integrated elastic cross sections, that tend to finite asymptotic limits, are studied, showing that the hypothesis of a black disk limit in the behaviour of the pp interaction seems to be excluded by phenomenology.
In Sec. IV we give predictions, presented in tables I and II, in equations and in figures. Taking into account previous publications at 1.8 and 7 TeV, the present paper give explicit predictions of cross sections at 2.76, 8 , 13 and 14 TeV, with no free numbers. More precise future data may confirm our predictions more firmly.
We also discuss the geometrical interpretation of our amplitudes, showing that the effective interaction radius in b-space increases with the energy. Our amplitudes obey a geometric scaling in asymptotic energies, and indicate that the profile function d 2 σ inel /d 2 b tends to a universal (energy independent) function with respect to a scaling variable, x ∼ b/ √ σ. This universal function exhibits a considerable diffused surface, indicating a scenario different from the commonly accepted black disk. At LHC energies, the saturation seems to start (the central value of d 2 σ inel /d 2 b is almost unity), but the asymptotic profile is still far and only can be reached for √ s > 10 4 TeV. The connection between the diffused surface of long range and inelastic diffractive processes will be an interesting line of investigation.
We believe that our analytic representation of the scattering amplitudes will serve as important guidance, not only for the future measurements in LHC, but also for a theoretical understanding of the intermediate region of partonic saturation phenomena.