Elastic amplitudes studied with the LHC measurements at 7 and 8 TeV

Recent measurements of the differential cross sections in the forward region of pp elastic scattering at 7 and 8 TeV show the precise form of the t dependence. We present a detailed analysis of these measurements including the structures of the real and imaginary parts of the scattering amplitude. A good description is achieved, confirming in all experiments the existence of a zero in the real part in the forward region close to the origin, in agreement with the prediction of a theorem by Martin, with an important role in the observed form of dσ/dt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{d}\sigma /\mathrm{d}t$$\end{document}. A universal value for the position of this zero and regularity in other features of the amplitudes are found, leading to quantitative predictions for the forward elastic scattering at 13 TeV.


LHC experiments in elastic pp scattering
With an enormous gap in the center-of-mass energy with respect to previous data in pp and pp scattering, the Totem and Atlas experimental groups at LHC have recently measured dσ/dt in the forward t ranges at √ s = 7 and 8 TeV [1][2][3][4][5]. These measurements offer a unique opportunity to investigate the behavior of p-p collisions at the highest energies reached in laboratory. A detailed and precise analysis of these data can establish a precious milestone for the understanding of the high energy behavior of p-p interactions. The datasets and their t ranges are listed in Table 1, where we use the obvious notation T7, T8, A7, A8 to specify Totem (T) and Atlas (A) Collaborations and center-of-mass energies 7 and 8 TeV.
In order to build a bridge towards theoretical models aiming at the understanding of the dynamics, it is important that the analysis of these LHC data be made with identification of the structures of the individual parts of the complex scattering amplitude. The disentanglement of the two terms in the a e-mail: kendi@if.ufrj.br observed modulus dσ/dt is the crucial task. At each energy, parameterizations must search to exhibit clearly the properties of magnitudes, signs, slopes and zeros of the real and imaginary parts. External support, as dispersion relations and connections with analyses at other energies, gives important clues. The intervention of the electromagnetic interactions must be treated coherently with a proposed analytical form for the nuclear part, and account must be taken of phase of the Coulomb-Nuclear Interference (CNI), which is calculated in Appendix A.
In the present work we perform a detailed examination of the data trying to satisfy these requirements. Each part of the amplitude is written with an exponential factor with a slope, multiplying a linear term in t, thus with three parameters. These analytical forms are sufficient to describe the properties of the nuclear parts. The six parameters for each dataset are studied using fits to data with appropriate statistical control. Correlations are studied, and resulting values are proposed for each dataset. Good description of the measurements is obtained, with details in the shape of the forward diffractive peak, exhibiting the zero of the real part predicted in the theorem by Martin [6], and with observation in the forward range of the ingredients that construct the imaginary zero responsible for the dip in dσ/dt observed when data exist at higher |t|.
According to analyses and models of pp and pp scattering at high energies [7][8][9][10][11][12][13][14][15][16], including full-t ranges, the imaginary part has a zero near the marked dip observed in dσ/dt at about 0.4-0.6 GeV 2 , while the real part starts positive at t = 0, has a zero for small |t| (Martin's theorem), and a second zero after the dip. The descriptions of the models differ for higher |t| about the position of this second real zero and about the existence of further imaginary or real zeros. Models describing large |t| ranges have different motivations and dynamical structures, and may be analytically very sophisticated, trying to represent the observed shapes of dip, bump, tail, in Table 1 Measurements at √ s =7 and 8 TeV from Totem and Atlas Collaborations at LHC [1][2][3][4][5]. The ρ values for a b are taken from COM-PETE Collaboration [19]. For c the ρ value was obtained by the authors with a forward SET-I and kept fixed in a complete SET-II as explained in the section with the analysis of T8. For d the ρ value is taken from [ the angular dependence. However, in the forward range the analytical structure required to describe the data may be very simple. In the present work we show that real and imaginary parts including exponential slopes and linear factors, combined with the electromagnetic interference, contain the essential ingredients for a precise representation. The second real zero occurs outside the studied range and does not influence the analysis. Particularly due to the small magnitude of the ρ parameter, the disentanglement of the two parts is not trivial, requiring careful analysis, and still leaving room for some subjective but physically reasonable choice. Once the amplitudes are identified, the results provide a necessary connection between data and theoretical models of microscopic nature for the strong interaction dynamics. The deviation from the pure exponential form in dσ/dt is obvious beforehand, since dσ/dt is a sum of two independent squares. More clarity can be obtained in the analysis with the identification of the two parts of the complex amplitude and their control by fundamental constraints (dispersion relations, Martin's theorem for the real part, zero of the imaginary part anticipating the dip).
Of the four datasets, T8 is the only one reaching very small |t|, allowing for a more complete investigation of some details, such as the influence of the Coulomb phase. However, a comparative analysis of the four cases is extremely important, since on observing coherence in some characteristics we may believe that there is reliability in the descriptions. The energy dependence of pp elastic scattering is very smooth, and features of 7 and 8 TeV data must support each other in a unified treatment.
A theorem by Martin proves that the real part has a zero close to the origin [6]. The abstract of the paper says, literally We show that if for fixed negative (physical) square of the momentum transfer t, the differential cross section dσ/dt tends to zero and if the total cross section tends to infinity, when the energy goes to infinity, the real part of the even signature amplitude cannot have a constant sign near t = 0. @ 1997 Elsevier Science B.V.
Thus the real part in pp and pp scattering has a zero close to the origin, with location approaching t = 0 as the energy increases. This constraint has been confirmed in previous analyses at LHC and lower energies [7,8], with the conclusion that the position of the first zero of the real part behaves like |t R | = A + 1/(c 0 + c 1 log s + c 2 log 2 s).
Although the analytical properties of the amplitude are defined in t-space, the insights for the construction of theoretical models are natural in the geometrical space, where the physical intuition to build amplitudes may be represented. The Fourier transformed space is appropriate to study asymptotic properties of the cross sections such as questioning whether the proton behaves as a black or a gray disk. The profile functions are also convenient tools to study the unitarity constraint. Of course, we recognize that the amplitudes written for the short forward t-range cannot lead to a sufficient understanding of b-space properties. Even so, we believe that the relationship is important, and the Fourier transformation of our amplitudes to b-space is analytically performed and its properties are discussed in Appendix B.
This paper is organized as follows. Section 2 describes the formalism of the proposed model; Sect. 3 presents the results of the model fits to the four LHC measurements; Sect. 4 summarizes the numerical analysis; Sect. 5 concludes the paper.

Amplitudes and observables in forward scattering
In the analysis limited to the forward ranges shown in Table  1, the expectations are satisfied writing the differential cross section in the form where t ≡ −|t|, α is the fine structure constant and (hc) 2 = 0.3894 mb GeV 2 . This expression is applied for pp an pp, and the parameters are specific for each case. F C (t) and αφ(t) represent the form factor and phase of the Coulomb interaction. The phase αφ has opposite signs for pp and pp scattering. The real electromagnetic amplitude is given in terms of the proton form factor for the pp/pp collisions. The proton form factor is taken as where 2 = 0.71 GeV 2 . The phase of the Coulomb-Nuclear interference is discussed in Appendix A.
We have thus assumed the imaginary amplitude with an exponential that accounts for the forward diffractive peak and a linear factor that accounts for the zero that occurs near the dip in dσ/dt, writing the simple form and The influence of the parameter μ I depends on the range of the data analyzed. To include the influence of the first real zero, we write and In order to check the influence of a second zero in the real part, we could add in the amplitude a term ξ R t 2 , but actually it has no effect in the present analysis. The normalization is defined by and for the pure nuclear interaction At t = 0, we have the usual definition of the ρ parameter With positive ρ and negative μ R (this is what we have in pp at high energies, as our analysis shows), there is a zero in the real amplitude, namely Martin's zero, located at The position of this zero and the magnitudes of the real and imaginary amplitudes in its neighborhood are responsible for details in the deviation of the differential cross section from a pure exponential behavior. The derivatives of the nuclear amplitudes at t = 0 are and The average slope measured directly in dσ/dt is the quantity We remark that parameters are determined fitting data in limited |t| ranges, at finite distance from the origin, so that the values obtained depend on the analytical forms (4), (5), (6), and (7) of the amplitudes. In particular, the slope parameter usually written in the form dσ/dt = σ 2 (ρ 2 +1) exp(Bt) does not agree with the expression for the differential cross section as sum of two independent squared magnitudes, each with its own slope. The assumption that B R and B I are equal is not justified. The average quantity B alone gives rough and unsatisfactory information. The importance of the different slopes in the analysis of pp elastic scattering has been investigated in the framework of the so-called dispersion relations for slopes [17]. It is important to note that also the Coulomb-Nuclear phase φ(t) depends essentially on the form of the nuclear amplitudes [20,21]. In Appendix A, generalizing previous work [7,[22][23][24], we derive the expression for the phase to be used with the assumed amplitudes written above.
Of course the six parameters are correlated, and in the present work we investigate the bounds of the correlations. We attempt to identify the values of parameters that may be considered as common representatives for different measurements. We show that the differences between the two experimental collaborations may be restricted to quantities characterizing normalization. The question of normalization is essential, and our inputs are the values of dσ/dt given in the experimental literature [1][2][3][4][5].
The extraction of forward parameters in the pp scattering has difficulties due to the small value of the ρ parameter, and consequently has suffered in many analyses from neglect of the properties of the real part. In our view the values of σ , ρ, B appearing in universal databases [18,19] as if they were direct experimental measurements should give room for critically controlled phenomenological determinations. A proper consideration for the properties of the complex amplitude is necessary. We observe that the properties that B R = B I and of the presence of zeros are common to several models [7][8][9][10][11][12][13][14][15][16]. The determination of the amplitudes for all |t| in several models is coherent.
We observe that the polynomial factors written in the exponent in some parameterizations of data [3,4] correspond to the linear and quadratic factors mentioned above, if the assumption is made that they are much smaller than 1 and can be converted into exponentials. However, this substitution is not convenient, because it does not show explicitly the essential zeros, and it also gives an unsatisfactory parameterization that cannot be extended even to nearby |t| values.
We thus have the framework necessary for the analysis of the data, with clear identification of the role of the free parameters. The quantities to be determined for each dataset are σ , ρ, B I , B R , μ I , μ R .

Data analysis
The range of -t covered in this analysis corresponds to a forward region, with |t| from |t| min = 0.00074 GeV 2 for T8 to |t| max = 0.3709 GeV 2 for T7. In this range the Coulomb effects play an important role and the relative Coulomb phase is taken into account. We compare results for the relative phase φ = φ ff , calculated in Appendix A with proton form factor (ff) and the reference case of phase zero φ = φ 0 = 0. This alternative is examined to set reference values because of a possible lack of understanding of electromagnetic effects, as in the calculation of the phase and the possible influence of the proton radius at high energies.
The statistical methods used in the analyses are performed with CERN-ROOT software [27], accounting for statistical and systematics uncertainties. However, since the values of χ 2 do not change much compared with the statistical uncertainties only, we understand that the statistical errors are sufficient for the analyses with our amplitudes.
We also study the correlation between the parameters, which is an useful tool to control possible instabilities of the fits. The correlations between the parameters are defined as where α and β are any two parameters, the brackets are the expectation values of the fitted parameters and σ α and σ β the variances associated with them. This correlation criterion is known as the Pearson coefficient and it varies from − 1 to 1, where − 1 is a complete anti-correlation, meaning that if one parameter is increased the other one decreases, and 1 is a complete correlation which means that if one parameter is increased the other one also increases. If the correlation coefficient is zero the parameters are said to be non-correlated, or independent. The correlation factors help us to understand the relation between the determination of the parameter and the range of |t| where this determination is performed.
With the purpose of identifying generic or universal values for parameters, for all measurements we study four different conditions in the fit: • (I) all six parameters are free; • (II) fixing ρ at 0.14, as suggested by dispersion relations; • (III) fixing μ I from the expected positions of imaginary zero [7,8] and dip in dσ/dt; • (IV) fixing simultaneously ρ and μ I at the above values.
We present our analysis for the four experiments separately in the next subsections. Since T8 has more precise data in the very forward region and the experimental paper has provided a detailed description of the observed structure, we investigate it in more detail. For this purpose we introduce a new diagram to represent the structure in the data at low |t|, plotting the ratio T 2 R /T 2 I against |t|, with This ratio does not depend on the total cross section, and therefore normalization uncertainties are canceled, allowing for identification of the zero in the real amplitude.
The results of the analysis are presented in Table 2. The headings of the table indicate the quantities determined in fits, namely the six parameters σ ,ρ,B I ,B R , μ I and μ R . The other columns give the derived quantities t R and B eff I , and the estimated χ 2 values. The first three columns specify the measurement, with the number N of data points, and the phase option (either the true phase φ ff or the reference option of zero phase).

T8
The T8 dataset contains measurements with two different optics. The first set of data (30 points) covers a very forward region, 0.000741 < −t < 0.191 GeV 2 , and the second set (30 points) starts a bit later, 0.0285 < −t < 0.195 GeV 2 , overlapping partially with the first, but with better statistical precision. They analyzed separately the first set, here called In all cases comparison is made with results obtained with Coulomb-Nuclear interference phase φ put equal to zero (φ 0 = 0) and calculated with electromagnetic proton form factor (φ = φ ff ) as  Table 2 continued SET I, with 30 points, and the combination of the first and the second measurements, with N = 60 points, here called SET II. Analyzing SET I they obtained ρ = 0.12 ± 0.03. In the analysis with SET II, ρ was not independently determined, but it was rather fixed at the same value of SET I. Repeating this analysis in two steps, using our nuclear amplitudes and testing phase values φ = φ 0 = 0 and φ = φ ff , we find that the determination of the ρ parameter is critical, strongly influenced by points in the most forward region. Using Condition (I), with all six parameters free, the results are the following. With SET I we obtain ρ = 0.106 ± 0.021 for φ ff case, which is a little bit below the expectation 0.14 from dispersion relations With the zero reference phase φ 0 we obtain ρ = 0.133 ± 0.021. With the complete SET II the ρ values come out considerably smaller than in SET-I for both φ 0 and φ ff phases, deviating still more strongly from 0.14. These results show the difficulty of T8 for free independent determination of the parameter ρ.
SET II is used in the description that follows, and the resulting parameters are registered in Table 2.
Using Condition (I) we obtain equivalent μ R = −2.78 GeV −2 for both φ ff and 0 phases, but since the values of ρ are different the position of the zero t R = ρ/μ R changes from − 0.032 ± 0.004 to − 0.040 ± 0.004 GeV 2 for the φ ff and φ 0 phases, respectively. We remark that precise data in the vicinity of the zero is important for the determination of ρ.
In Fig. 1 we present the fit results of the T8 data with all parameters free, as under Condition I, and phase φ ff . The parametrization is able to describe the T8 measurements. With plots showing local displacements, shown in Fig. 2, we are able to exhibit details of the structure of the amplitudes. In the LHS we plot a diagram similar to that presented by the experimental group [3], showing the valley structure appearing when we subtract the simple exponential form ref = A e Bt from the best fit solution for the differential cross section. The structure appears neatly, but the band of normalization errors is very large, and the roles played by the amplitudes are not clear. In the right hand side of Fig. 2 shows the turning point due to the zero of the real part of the amplitude, with a much narrower band of the systematic errors due to the cancellation of σ dependence in the ratio. The interplay of the magnitudes of the real and imaginary amplitudes influences the structure of dσ/dt. The deviation of a pure exponential behavior is inherent to the sum of two independent squared quantities.
From the fits performed in T8 we obtain from the correlation coefficients that ρ is weakly correlated with all other parameters. This is the reason why, when we fix its value according to Conditions (II) and (IV), we obtain a value of  Fig. 2 The left plot shows the non-exponential behavior of the differential cross section for T8. The figure is obtained subtracting from the best fit of the differential cross section a reference function which is dσ/dt written with a pure exponential form and dividing the subtraction by this reference function. The dashed lines show the normalization error band in dσ/dt, which is quite large. The plot in the RHS shows the ratio T 2 R /T 2 I , which exhibits information of a non-exponential behavior with advantages compared with the first plot, since σ is canceled, and with it most of normalization systematic error, and the second advantage is the clear indication of the structure of the real amplitude showing the position of its zero |t| R χ 2 larger than with Condition (I): the other parameters are not able to compensate for the change in ρ.
Under both Conditions (I) and (II) the parameter of the imaginary part μ I is very small, which means that the trange in the data is not far enough to feel the zero of the imaginary amplitude. On the other hand, fixing μ I under Condition (III), ρ is seen as strongly dependent on the phase.
It is important that, although the parameter μ I has low Pearson coefficients, the statistical error associated with this parameter is large, which means that it cannot be well determined with Condition (I) or (II). The χ 2 values do not change considerably from Condition (I) to Condition (III), but the parameters B I , B R and μ R change. The presence of non-zero negative μ I forces the imaginary amplitude to point towards zero. The real amplitude must compensate for the decrease of the imaginary part, reducing the magnitude of the real slope, and the value of μ R is also affected. The imaginary slope compensates for the increase in the magnitude of μ I introduced in (III), thus preserving the value of B eff I . Figure 3 shows the correlation maps between pairs of parameters. These figures correspond to Pearson coefficients corr(B R , ρ) = 0.26, corr(μ R , ρ) = −0.39 and corr(B R , μ R ) = −0.98 showing weak correlations for the first two cases, and a strong anti-correlation for the latter. The lines represent the allowed regions at different standard deviations. Since ρ is weakly correlated with the parameters of the real part according to Condition (III), under Condition (IV) we expect small deviation in B R and μ R . Under this condition we obtain very similar results for the free parameters for both choices of phase.

A8
In A8 [5] Fig. 4 and we can clearly note that the position of the Martin zero is in agreement between the measurements. Important differences between T8 and A8 are in the parameters of the imaginary part. Thus μ I = −1.65 GeV −2 for A8, while for T8 it is compatible with the zero due to the short |t| range of the data.
The values of σ in A8 and T8 differ by ∼ 5%. This difference may be due relative normalization, and we may wonder whether a unification of dσ/dt through a constant factor could unify the solutions for σ , while leading also ρ and μ R to common values.
Since the fits for A8 show a strong anti-correlation factor between σ and ρ, and σ is a stable parameter, we use Condition (II) to fix ρ at 0.14 and obtain a value of χ 2 very similar to the result obtained under Condition (I) with free ρ, which is natural since in (I) ρ in not very different from 0.14. Under Condition (II), with fixed ρ, the results for T8 and A8 give compatible μ R values.
In A8, for both (I) and (II), the μ I parameter predicts |t I | = 0.59 GeV 2 , which is far on the right of the region where the dip in dσ/dt is expected to occur. By fixing μ I in Condition (III), B I for A8 agrees with the value found in T8, as expected, since values of B eff I are in agreement for all conditions.
The use of Condition (IV) does not change considerably the parameters of A8 when compared with Condition (III), because the central values of ρ in (III) are close to 0.14.

T7
At √ s = 7 TeV the Totem Collaboration measured elastic p-p cross sections with two sets [1], the first in the range (0.005149 < −t < 0.3709 GeV 2 ) with 87 points, and the second in the range 0.377 < −t < 2.443 GeV 2 with 78 points. Using our expressions for the forward amplitudes, we analyze the forward set, obtaining the results shown in Table  2. With all parameters free, for both φ ff and φ 0 phases the ρ values are less than 0.14. The values for μ I are similar to those in A8. The total cross section is higher but compatible with the original paper [1], as given in Table 1.
The correlation factor corr(σ, ρ) = −0.95 obtained under Condition (I), shows a strong anti-correlation. Since σ is a stable parameter and the central ρ value is larger than the expected 0.14, Condition (II) fixing ρ, leads for all other parameters and for χ 2 results similar to Condition (I).
With Condition (III) we observe that the central ρ values are smaller than obtained with (I), and the statistical errors are about 55% (φ 0 ) and 88% (φ ff ) of the central value.

Analysis of an extended set
The special availability in T7 of data beyond the |t| ≈ 0.4 GeV 2 limit of the measurements in Table 1, may be used to study a range where the μ I parameter becomes more effective, pointing towards the zero of the imaginary amplitude, and being determined with higher accuracy. We thus add to the forward set the first 17 points of the second dataset, reaching |t max | = 0.52 GeV 2 . The results of the analysis with the combined set of 104 points are presented in the Table 3.
The study of the extended set with Condition (I) leads to the correlation map for the quantities B I and μ I shown in Fig. 5. The Pearson coefficient depends on the range on t where the fit is made and of course the fit conditions used, and for this extended range it shows a slight anti-correlation between μ I and B I .
In Fig. 6 we compare the squared magnitudes of the real and imaginary amplitudes, and also the simplified single exponential amplitude assuming the effective slope B eff I , for the extended T7 with Condition (I). We observe that the imaginary part starts to deviate from the simple exponential near |t| ∼ 0.2 GeV 2 , and at |t| ∼ 0.47 GeV 2 it passes through a zero. Beyond this range the real amplitude would be modified to incorporate other terms (say a quadratic ξ R t 2 ) that may play an important role in the construction of the dip structure.
As with the forward set, here the ρ values come out larger than 0.14, but within error bars. Fixing ρ = 0.14 under Condition (II), the changes in parameters and in χ 2 are very small.
It is very important that the parameter μ I is here determined with more precision and is compatible with the value fixed to establish Condition (III) in Table 2.
The use of Condition (IV) for this extended set, shown in Table 3, leads to a decrease in the magnitudes of B R and μ R when we compare with the results of the forward set under the same condition. Since ρ is fixed, the decrease in μ R implies an increase in the magnitude of |t R |.

A7
In A7 the Atlas Collaboration measured [2] 40 points in the region 0.0062 < −t < 0.3636 ( GeV 2 ). This experiment is challenging, since the experimental authors recognize that ρ cannot be determined from the data with usual forms of amplitudes. The natural result for ρ as a free parameter is a negative quantity, and we run Condition (I) imposing a lower bound 0 for ρ, and of course the expected value is zero, with large error bars. Under this condition the total cross section for φ ff phase is σ = 95.77 mb, which is well below, while the parameters B I , B R , μ R and μ I are compatible with T7. In these measurements the correlation factors between ρ and all other parameters are very small. Thus with Condition (II) Table 3 Fits for the extended T7 data (N = 104 points in the range 0.00515 < −t < 0.52 GeV 2 ) according to Conditions (I), (II) and (IV)   we expect an obvious increase of χ 2 , but the parameters B I , B R and μ R are compatible with T7. With Condition (I), μ I is compatible with zero, with large error. With Condition (III), χ 2 is not much changed, as well as the total cross section. The ρ value is still undetermined but B I , B R and μ R are compatible with the other experiments. Thus, in spite of the smaller value of total cross section and instabilities in the determination of ρ, Condition (IV) shows a similar t dependence in A7 with respect to the other datasets, with a similar behavior of the amplitudes.

Summary of the analysis
The LHC measurements at √ s = 7 and 8 TeV shown in Table 2 are analyzed under the assumption of the analytical forms for the scattering amplitudes given in Eqs. (4) and (6). These forms are considered to be as simple as possible under theoretical conditions to describe the scattering amplitudes in the forward region. The analysis aims at the determination of six intervening parameters, three for the imaginary part (σ , B I , μ I ) and three (ρ, B R , μ R ) for the real part, with expected smooth energy dependence. Only the T8 measurements cover the very small |t| range, so that the analysis is made with nonhomogeneous inputs, and four specific conditions, named (I), (II), (III), (IV) in Table 2, are studied separately.
The Coulomb-Nuclear interference is a crucial ingredient, and its phase φ ff is treated compatibly with the forms of the amplitudes, as presented in Appendix A. In order to have a reference (although it is not realistic) we also give results of fits with phase φ 0 put equal to zero.
As expected, the direct results of the fits with all parameters kept free are rather dispersive in some aspects, as shown in the sub-table with title Condition (I). The fitted values of ρ do not agree among the measurements within the statistical uncertainties, but once the normalization uncertainty is considered, the values are in agreement. We observe that the other parameters of the real part (B R and μ R ), related to the shape of the amplitudes but not so strongly to its normalization, appear to be more regular.
Under Condition (II) ρ is fixed at a reference value 0.14 suggested by dispersion relations for the 7-8 TeV range. Compared to (I) there is a loss in χ 2 in the cases with free ρ far from 0.14, namely T8 and A7, but not in the other cases, where we observe only a rather slight adaptation in the other real parameters: B R and μ R appear as regular quantities. Also B I remains the same (except in A7), and it is particularly important to remark that the effective slope B eff I = B I − 2μ I in Eq. (12), which compensates for the influence of μ I on the |t| dependence of the imaginary part, appears as a very regular common quantity. In spite of the differences, at this point we may say that we are lead accept the value ρ ≈ 0.14 for all measurements. On the other hand, μ I does not remain regular among the experiments under this Condition (II). Since μ I is responsible for the presence of a zero (t I = 1/μ I ) in the imaginary amplitude (which occurs near the dip in dσ/dt at about 0.4 GeV 2 ), it is natural that T8 (limited to 0.2 GeV 2 ) is not sensitive to μ I , and puts it at zero under Conditions (I) and (II).
It is interesting to observe the effects of correlations. For instance, in T8 and A7, ρ is weakly correlated with the other parameters. Then fixing ρ under Condition (II) worsens strongly χ 2 for these datasets, more than in sets A8 and T7 where the strong correlation between ρ and the other parameters absorb the effects of the fixing condition. For T8 in the φ ff case the worsening in χ 2 is more dramatic, corroborating our concern about proton form factor and Coulomb phase.
With Condition (III) we fix μ I according with the expected positions of the imaginary zero and dip in dσ/dt [7,8], and let ρ free. This is successful, as ρ results about the same or improved with respect to Condition (I) (all parameters free), except for A7, as expected. The parameters B R and μ R of the real part remain the same, and it is remarkable that B I changes, becoming very regular, absorbing the influence of μ I (now fixed) and keeping the constant and regular effective |t| dependence of the imaginary amplitude, represented by B eff I . The μ I parameter determines a zero at t I = 1/μ I in the imaginary amplitude and is related with the position of the dip in the differential cross section. A precise determination of μ I depends on the existence of data in a region near the dip. Thus in the T8 dataset, without points for small |t|, the central values of μ I are near zero in (I) and (II), while in A8, T7 and A7 the values have larger magnitudes. We study this question analyzing the T7 experiment with inclusion of a second set of points [1]. Forming a larger dataset in the range 0.005149-0.52 GeV 2 , the best solution with fixed ρ = 0.14 gives μ I = −2.15 GeV −2 shown in Table 3, in good agreement with the prediction [8]. We are thus lead to Conditions (III) and (IV), which fix μ I .
Fixing both ρ and μ I at their expected values under Condition (IV) we obtain good modeling for all measurements, except for the total cross sections, which separate Atlas from Totem. It is particularly significant that the position of the zero of the real amplitude is nearly the same for all cases.
The deviation from pure exponential form in the differential cross section is interpreted as due to the shape difference between the real and imaginary amplitudes. The T8 experiment presents very precise data at low |t|, showing a valley behavior in the differential cross section, while the T7, A7 and A8 data also indicate a structure for the data at low |t|, but with large uncertainties. The shape is constructed when the real amplitude crosses zero, passing to negative values. After the zero, the action of the real slope pushes this negative value back to zero, and the structure is formed.
This mechanism suggests that the determination of ρ depends not only on the extrapolation to the limit at |t| = 0, but that it also depends on the form of the real amplitude around its zero. A precise extraction of ρ depends both on the specific analytical model used and on the data in the whole range 0 < −t < 0.2 GeV 2 where the valley occurs.
The position of the real zero is obtained from the equation t R = ρ/μ R (values are given in Table 2). Important quantities are the derivatives of the amplitudes and their effective slopes determined in each measurement. The exponent B I written in Eq. (4) is not the logarithmic derivative of the imaginary amplitude, which is given by Eq. (12). The effective slope at small |t|, is seen in the approximation of the linear factor put in exponential form 1 − μ I t ≈ exp(−μ I t). Thus the determinations of μ I and B I depend on data in the very forward region, and also in the region near the dip. Table 2 shows the interesting regularity of the quantity B eff I , as opposed to B I . The average slope B measured in the differential cross section is given in Eq. (14). In our analysis B appears to be stable, with a value B 20.5 GeV −2 . Comparing our result with the values in Table 1 we see deviations of about 1 GeV −2 , and we thus remark that the measured average slope B depends on the parameters μ I and μ R , which are influenced by data in the large |t| region.
The quality of the representations of the data can be read from the pull plots in Fig. 7. The y-axis represents the standard deviations at each |t| and is defined as (dσ/dt i expdσ/dt the (t i ))/er i , where dσ/dt i exp is the experimental value at some t i with error er i , and dσ/dt the (t i ) is the theoretical value calculated at t i . Assuming that the statistical errors follow Gaussian distributions, the most probable solution should contain about 68% of the points within ±1σ (deviation) and about 95% of the points within ±2σ . For T8 we see that about 65% of the points are within ±1σ following this criterion, and for A8 about 74% of the points are within ±1σ . Of course, care should be taken in this analysis because a large number of experimental points are needed for good statistics.
For T7, about 65% of the experimental points are within ±1σ and about 94% within ±2σ of deviation with respect to the theoretical curve (fitted curve). Since these are the measurements with larger number of points (N=87), the maximum likelihood criterion for Gaussian statistical errors shows that our curve is a good representation of the data.
The regularity on the values of μ R is remarkable. The zero of the real part t R determined by the parameters μ R and ρ is associated with the predicted zero of the theorem by Martin [6]. We see that the position of the zero is stable in all experiments and in Table 4 we observe agreement at t R −0.037 GeV 2 within the statistical errors. In terms of amplitudes we observe that the position of the zero together with the magnitude of B R determines the structure shown in Fig. 2. The existence of this zero is very important for the superposition of the real and imaginary parts that controls the detailed structure of dσ/dt in the low |t| region. Finally, as a general remark we observe that total cross sections maintain that Totem values are higher than Atlas by 4-6 mb (about 5%), and the numbers remain stable under all conditions studied. The difference is not due to statistical fluctuation, but rather it is related to normalization procedures, both experiments informing a relatively high systematic uncertainty. Using Condition (I), we observe at 8 TeV Therefore, in our proposed model the discrepancy between the two measurements is less than 2 standard deviations. As another check, we observe that under Condition IV the systematic differences in σ do not affect the other parameters, which are all coherent.
Another general remark is that differences in the results obtained with the two assumptions for the Coulomb phase are relevant only for T8 that has experimental points for very small |t|. It may be fortuitous, but we observe that χ 2 is smaller when the phase is put at φ 0 = 0.
The difficulty in the determination of ρ may be due to Coulomb-Nuclear interference. The distribution of electric charge in the proton determined by electromagnetic scattering at low energies ( 2 = 0.71 GeV 2 ) may be unrealistic for high energy hadronic scattering. In Appendix A we show that with expanded proton size (as may be the case at high energies) the Coulomb phase decreases. We see that the problem of the Coulomb interference in elastic pp scattering still has open questions.

Conclusions
In this work we study the properties of the amplitudes in pp elastic scattering analyzing experimental data at the LHC center-of-mass energies 7 and 8 TeV, based on a model for the complex amplitude, with explicit real and imaginary parts, each containing an exponential slope and a linear factor to account for the existence of a zero. The zero of the real part, close to the origin, corresponds to Martin's theorem, and the zero of the imaginary part anticipates the dip in the differential cross section that occurs beyond the range of the available data under study.
Our study shows that the real amplitude plays crucial role in the description of the differential cross section in the forward region. Interference with the Coulomb interaction is properly accounted for, and use is made of information from external sources, such as dispersion relations and predictions for the imaginary zero obtained in studies of full-t behavior of the differential cross section [7,8]. We organize the analysis under four conditions, according to the specifications of the parameters with values fixed in each case. Comparison is made of the results obtained for the four experimental measurements. We obtain the results shown in Table 4, which we believe to be a good representation of the experimental data of Table 1.
Assuming σ ∼ log 2 s at high energies, dispersion relations give ρ ∼ π/ log s, and we know that the position of the first real zero [8] behaves like |t R | A + 1/(c 0 + c 1 log s + c 2 log 2 s). From Eq. (11) it then follows that μ R increases like −λ log s at high energies. With s in GeV 2 , from Table 4 we obtain λ 0.2 GeV −2 . Future precise measurements in LHC at 13 TeV may be investigated with this purpose. Our present work predicts that the zero at √ s = 13 TeV is at t R = −0.037 GeV 2 .
The quantity μ R is related with the scaling variable τ = t log 2 s introduced by Dias de Deus [28,29] connecting s and t dependences in the amplitudes at high energies and small |t|. Martin [30] uses the same idea of a scaling variable, writing an equation for the real part ρ(s, t) using crossing symmetric scattering amplitudes of a complex s variable, valid in a forward range. The proposed ratio is where f (τ ) is a damping function, with the implicit existence of a real zero. The form of f (τ ) determines the properties of the real zero [31], found in the analysis of the data. This may be a clue for the introduction of explicit crossing symmetry and analyticity in our phenomenological treatment of the data. Other models [9][10][11][12][13][14][15][16] also deal with the position of the real zero, discussing different analytical forms for the amplitudes, and it would be interesting to have their predictions for the amplitudes in the forward range.
In Appendix B we investigate the behavior the amplitudes in b-space for a geometric insight of the scattering process. We find that although our amplitudes never violate Eq. (B11) up to √ s = 8 TeV at LHC, the differential cross sections in b-space are all monotonically decreasing functions in b. This means that our amplitudes do not indicate any signal for peripheral dominant elastic scattering process. However, in an analysis of T8 the dominance of the square of the real amplitude in b-space in Ref. [32] was interpreted as periph- Table 4 Proposed values of parameters for the four datasets. The phase of the Coulomb interference is calculated with electromagnetic proton form factor as described in Appendix A. The T7 data are also shown with inclusion of points at higher |t| that are important for confirmation of the value of μ I Fixed quantities: ρ = 0.14, μ I = −2.16 GeV −2 (8 TeV) [8], In non-perturbative QCD, in several instances, the proton appears as a structure with expanding size as the energy increases [33][34][35][36][37][38][39][40][41][42][43][44][45][46][47], with varied mechanisms, as distribution of valence quarks in a cloud around a core, modifications in QCD vacuum in the region of the colliding particles, and so on. Together with the evolution of the proton hadronic size, its electromagnetic properties, as they appear in high energy collisions, may change also. A linear increase in log s is a usual assumption for the effective proton radius, and the form-factor parameter 2 would then be reduced by about 1/2, corresponding to an increase of about 40 % in the proton radius. In Appendix A we calculate the interference phase with this example.
We expect that future data in pp elastic scattering at 13 TeV and higher energies will have a high quality, covering a wide t range to allow for a determination of the properties of the real and imaginary amplitudes in pp elastic scattering, including studies of the amplitudes up to the perturbative tail of large |t|. Hopefully the experimental groups will receive the necessary support and encouragement for this effort.
We start from the exact expression for the phase the signs ∓ corresponding to the pp/pp systems. Cahn has considered the approximation justifying that the momentum transfer scale q 2 is much smaller than the inverse size of the proton. This is a very strong assumption, which should be tested in specific cases. Following the suggestion, Eq. (A8) becomes where for simplicity we use the t variable instead of q 2 . As a generalization with respect to Cahn's calculation, we take for the nuclear amplitude the same expression in Eq. (A1). Then we need to evaluate the integrals where we have used the definition The results of the integrations (N = 0, 1, 2) are where γ = 0.5772 is the Euler gamma constant. The phase is then written with real and imaginary parts, respectively, and Equations (A15) and (A16) are our final results for the phase calculated with form factors, in a generalization of the work by Cahn [24], assuming more complete structures for the real and imaginary parts of the scattering amplitude. It may be of practical usefulness to define and and then write φ R (t) = ∓ 1 1 + ρ 2 C I + ρ C R + γ + log(−t) (A19) and φ I (t) = ∓ 1 1 + ρ 2 ρ C I − C R . (A20) It must be observed that in these expressions b R , b I and −t have compatible units, as GeV −2 and GeV 2 . The result is simple: in the real part the t dependence is purely linear in log(−t), the imaginary part is a very small constant, and there is no explicit energy dependence.
In the simplified case μ R = μ I = 0, ξ R = ξ I = 0, we obtain Cahn's original form. The real and imaginary parts of the phase are plotted in Fig. 8 with an example of values for the parameters.
In the figure we plot also a (dot-dashed) line representing the real part of the phase calculated for a larger proton, with 2 = 0.305 GeV 2 .
of the amplitude T (s, t) in Eqs. (4) and (6) is given bỹ and Figure 9, withT at √ s= 8 TeV as a function of b for T8 and A8, shows the behavior of the amplitudes, which differ by only a few % in the two measurements. We note that the magnitude of real parts become comparable and even larger than the imaginary parts for large b. This behavior, which occurs also for the 7 TeV datasets, is not observed at lower energies.
To investigate in more detail the significance of the dominance of the real part in the peripheral region, we introduce where χ R = Re (χ ) and χ I = Im (χ ). These b-space representations of the differential cross sections offer a geometric view of the pp interactions, although such interpretation should be taken with caution, since they are not physical observables. From the unitarity condition of the scattering amplitude, we must have χ I ≥ 0, and dσ inel /d 2 b ≤ 1. For b > 10 GeV −1 we have T 1, which implies |χ | 1. Thus, in this region, taking up to the leading orders in χ R and χ I , respectively, we obtain dσ el (s, b) From the above, it is clear that in the domain where T 1 (that is, where the Born approximation is valid), the inelastic part comes totally fromT I , while the elastic contribution may exceed the inelastic one if Thus the appearance of the peripheral domain whereT R > T I may indicate that the contribution of elastic scattering is going to be significant in the peripheral region, and even can be dominant if Eq. (B11) is satisfied. If such a situation occurs, the scattering becomes basically elastic. Figure 10 shows various b-space differential cross sections using the amplitudes obtained in the present analysis. Note that the corresponding b-space amplitude satisfies the unitarity condition mentioned before, and 0 ≤ dσ el /db, dσ inel /db ≤ 1.
As seen from the figure, in all cases, although the real part amplitude becomes dominant over the imaginary part in the peripheral region, Eq. (B11) is far from being satisfied, and the inelastic contribution is totally dominant over the elastic contribution at these energies (7 and 8 TeV).