arxiv: v1 [hep-ph] 21 Nov 2018

Transcription

1 Model-independent femtoscopic Lévy imaging for elastic proton-proton scattering T. Csörgő, 1,, 3, R. Pasechnik, 4, 5, 6, and A. Ster 3, 1 CERN, CH-111 Geneva 3, Switzerland EKE KRC, H-300 Gyöngyös, Mátrai út 36, Hungary arxiv: v1 [hep-ph] 1 Nov MTA Wigner FK, H-155 Budapest 114, P.O.Box 49, Hungary 4 Department of Astronomy and Theoretical Physics, Lund University, SE-3 6 Lund, Sweden 5 Nuclear Physics Institute ASCR, 5068 Řež, Czech Republic 6 Departamento de Física, CFM, Universidade Federal de Santa Catarina, C.P. 476, CEP , Florianópolis, SC, Brazil Abstract A model independent Lévy expansion method is introduced to describe nearly Lévy shaped squared moduli of Fourier transforms. We apply this method to precisely characterize the most recent elastic scattering data of proton-proton collisions at s = 13 TeV. The results reveal a substructure of protons, which is found to be significantly larger and darker in high energy proton-proton collisions at the TeV scale, as compared to a rather faint and apparently overlooked substructure, that we also identify in elastic proton-proton scattering at the ISR energy range of s = 3 6 GeV. PACS numbers: Dz, Lg, Pj, 4.30.Wb tcsorgo@cern.ch Roman.Pasechnik@thep.lu.se Ster.Andras@wigner.mta.hu 1

2 Introduction Lévy stable distributions are ubiquitous in Nature. They correspond to the applicability of generalized central limit theorems, when the final distribution is a convolution of several apparently random elementary processes [1, ]. In high energy particle physics, squared modulus of a Fourier-transformed source distribution may have an approximately Lévy (or, stretched exponential) shape for example in the analysis of the two-particle Bose-Einstein correlation functions [3], as well as in the case of the differential cross-section of elastic scattering at low four-momentum transfers [4]. Recently, the TOTEM collaboration presented preliminary data on elastic proton-proton scattering [5] at the Large Hadron Collider (LHC) energy of s = 13 TeV, the currently largest colliding energy ever achieved in a laboratory experiment. The precision and the range of this data set provides a challenge for a statistically acceptable, quantitative description, which is a precondition for a reliable imaging process that aims to determine also the statistical significance of the analysis results. Here we present a precise imaging of this TOTEM data set with the help of a novel and model-independent Lévy expansion method, that has its roots in precision analysis of Bose-Einstein correlation functions in high energy particle and heavy-ion physics [6 8]. The results reveal not only the structure (the impact parameter dependent shadow or inelastic profile function) of the protons on the femtometer scale, but for the first time the shadow profiles of a substructure inside the protons is also reconstructed and its energy evolution is investigated from the top TeV LHC colliding energy of s = 13 TeV down to the lowest ISR energy of 3.5 GeV. Formalism The differential cross-section of elastic scattering is proportional to the modulus square of the forward scattering amplitude. The TOTEM Collaboration established [5, 9], that this differential cross-section deviates significantly from a simple exponential shape, dσ/dt = A exp( B t ) at low values of the four-momentum transfer t = (p 1 p 3 ). We introduce the exponent α to characterize this deviation by a single parameter, so we assume that to leading order, dσ/dt = A exp [ (R t ) α ]. This form is frequently called as the Fourier-transform of a symmetric Lévy-stable source distribution, sometimes also called as a stretched exponential distribution. We have recently shown that this approximation with fixed α = 0.9 gives a statistically acceptable description of proton-proton elastic scattering at low-t from s = 3.5 GeV to 13 TeV [4]. At larger absolute values of the four-momentum transfer t, an interference or dip-bump

3 structure is also apparent in these differential cross-sections, that indicates a rather significant deviation from the symmetric Lévy source distributions and from the approximate applicability of generalized central limit theorems. To quantify the deviations from a symmetric Lévy source distribution, let us utilize an orthonormal series expansion, with complex coefficients c i, where the Lévy polynomials are introduced with the help of a Gram-Schmidt orthonormalization as follows: dσ 1 dt = A w(z α) + c j l j (z α), (1) j=1 w(z α) = e zα, z = t R 0, c j = a j + i b j, () as detailed in Ref. [4]. Here w(z α) stands for the Lévy stable weight function, the dimensionless variable z is t multiplied by the Lévy scale parameter R squared. The real and the imaginary parts of the expansion coefficients are denoted by a i and b i, respectively, and l j (z α) stands for the normalized Lévy polynomial of order j. The n-th moment of w(z α), µ α n is µ α n = 0 dz z n e zα = 1 ( ) n + 1 α Γ, (3) α where Γ(x) = 0 dz z x 1 e z stands for Euler s gamma function. The Gram-determinants D j D j (α) are D 0 (α) = 1, D 1 (α) = µ α 0, (4) D (α) = det µα 0 µ α 1, (5) µ α 1 µ α µ α 0 µ α 1 µ α D 3 (α) = det µ α 1 µ α µ α 3,... etc. (6) µ α µ α 3 µ α 4 The normalized Lévy polynomials l j (z α) are introduced in terms of the Gram-determinants and the unnormalized Lévy polynomials [4] as follows: l j (z α) = D 1 j D 1 j+1 L j(z α), for j 0. (7) In this equation the unnormalized Lévy polynomials are denoted as L i (z α). These poly- 3

4 nomials were introduced in ref. [8] and are defined as follows: Observables L 0 (z α) = 1, L 1 (z α) = det µα 0 µ α 1, (8) 1 z µ α 0 µ α 1 µ α L (z α) = det µ α 1 µ α µ α 3,... etc. (9) 1 z z The elastic differential cross-section is related to the modulus squared Fourier transformed forward scattering amplitude T el ( ) as dσ dt = 1 4π T el( ), = t. () We may assume, that the forward scattering amplitude is predominantly imaginary at high energies. In this case, the zeroth order normalization coefficient can be absorbed to the overall normalization coefficient A, and the following Lévy expansion is obtained for the complex-valued forward scattering amplitude T el ( ): T el ( ) = i ( ) 4πA e zα / 1 + c i l i (z α) The total and elastic cross-sections are given as σ tot Im T el (0) = 4πA σ el 0 d t dσ dt = A R [ 1 α Γ i=1 ( 1 + ( ) 1 + α z= R. (11) ) a i l i (0 α), (1) i=1 ] (a i + b i ). (13) The ratio of the real to imaginary parts of T el at t = 0 is ρ Re T el(0) Im T el (0 = i=1 b il i (0 α) 1 + i=1 a il i (0 α). (14) Imaging on the femtometer scale The angular distribution of elastic scattering provides pictures of the elastically scattered particles in impact parameter or b space. independent processes, t el (b) = 1 π i=1 For spin J 0 ( b) T el ( ) d, (15) is the impact parameter dependent forward scattering amplitude, where, b b and J 0 (x) is the zeroth-order Bessel function of the first kind. represented in an eikonal form as This t el (b) can also be t el (b) = i [ 1 e Ω(b)], (16) 4

5 ) dσ/dt (mb/gev pp: TOTEM - 13 TeV preliminary data Levy expansion (pos.def) σ tot = 115. ± 0.8 mb σ el = 31.4 ± 0. mb ρ = ± χ /NDF = 330 / 79, CL =.009e-0 R = 0.7 ± fm - A = ± 5.mbGeV α = ± a 1 = ± b 1 = ± a = ± 0.00 b = ± 0.00 a 3 = ± b 3 = 0.07 ± a 4 = ± b 4 = ± t (GeV ) FIG. 1. (Color online) Eq. (1) fitted to TOTEM dσ/dt preliminary data at s = 13 TeV in elastic proton-proton collisions. where Ω(b) is the so-called opacity function, which is in general complex. Images of the protons are given by the P (b) shadow profile or inelastic profile functions: P (b) = 1 e Ω(b). (17) where Ω(b) stands for the opacity function (known also as the eikonal). This way, the rather fundamental results of multiple diffraction theory [] when powered with the help of the newly found Lévy expansion of the forward scattering amplitude, together make it possible now to perform a high precision optical imaging of elastic hadron-hadron scattering and obtain pictures of the (sub)structure(s) of the protons on the femtometer scale. Results The model-independent Lévy expansion method is shown to describe the recently presented TOTEM preliminary data [5] of elastic proton-proton scattering at s = 5

6 ) dσ/dt (mb/gev pp data 7 TeV 13 TeV α A exp( -(R t ) α = t (GeV ) ) 3.5 GeV 30.7 GeV 44.7 GeV 5.8 GeV 6.5 GeV FIG.. (Color online) Eq. (1) with (a i, b i ) (0, 0) for i 1 describes the nearly exponential behaviour of dσ/dt of elastic proton-proton collisions, with α = 0.9 fixed, in the secondary cone after the dip and the bump, corresponding to.5 t 5 GeV at the ISR and 1.5 t 3 GeV at LHC energies, with parameters of Tables I and II. 13 TeV on Figure 1. These data span orders of magnitude on the vertical scale, indicate a nearly exponential cone region at low values of the four-momentum transfer t, that is followed by a dip and bump structure and a secondary tail. The Lévy expansion has a confidence level CL = %. This indicates that the Lévy expansion represents these data in 6

7 0.6 The shadow profile for the proton substructure TeV 7 TeV 6 GeV 3 GeV 0.0 b, fm FIG. 3. (Color online) Images and the imaging errors of the internal substructures inside the protons at various colliding energies. For s = GeV at ISR, the parameters of these substructures do not depend significantly on the center of mass energy of the collisions, see Table I. However, images of proton-substructures at the LHC energies of 7 and 13 TeV are significantly darker and larger, see also Table II. an acceptable manner. As far as we know, the precision of this data set is unprecedented in the field of elastic scattering in high energy physics and its description on a quantitative, statistically acceptable manner remains a challenge at present for all other, model dependent attempts. At s = 13 TeV, the maximal measured four-momentum transfer values are t max 4 GeV. This allows for a minimal spatial resolution of / t max 0.1 fm, and makes it possible to resolve some of the substructures inside the colliding protons at LHC energies. The model independent Lévy expansion technique also well suited to describe earlier measurements of elastic proton-proton scattering with s ranging from 3.5 GeV to 7 TeV, as detailed in Appendix A of ref. [4]. Each of the differential cross sections, summarized on Figure, lack secondary dip and bump structure. As far as we know, the significance of this observation was noted first in ref. [11]. As we detailed in Ref. [4], the data at t larger than the dip and bump region are described in a statistically acceptable manner with a zeroth order Lévy fit, corresponding 7

8 FIG. 4. (Color online) Images of protons and their substructures at various colliding energies. The shadow profile of the protons is growing with increasing colliding energies. At s = GeV, a shadow of a substructure is significantly smaller and fainter as compared to that of LHC energies of s = 7 and 13 TeV. to dσ/dt = A exp( ( t R ) α, with α = 0.9 fixed value. This value of α indicates a small but systematic deviation from the exponential shape. The results, summarized on Figure, reveal another striking and unexpected detail: at the ISR energy range of 3.5 s 6.5 GeV the fits with the zeroth order Lévy expansion yield nearly parallel lines. This indicates that neither the Lévy exponent α nor the Lévy scale R of this tail region evolves significantly with changing the colliding energy in this lower, ISR energy range. A corresponding substructure is reconstructed inside the proton with the help of eq. (17) as indicated on Figure 3. At ISR energies, this substructure is very faint, with maximum value of P (b = 0) of about Apparently, there are no statistically significant differences among the faint shadows of these proton substructures in this ISR energy region of s = 3.5 to 6.5 GeV. A significantly different, larger and darker substructure is however found in elastic proton-proton collisions at 7 TeV. Increasing the center of mass energy from 7 to 13 8

9 s (GeV) A (mb GeV ) ± ± ± ± ± ± R (fm) 0.30 ± ± ± ± ± 0.01 χ /NDF.8 / / 1 19 / / 1 1. / 1 CL [%] σ tot sub σ el sub (mb) 0.4 ± ± ± ± ± 0.03 (µb) 1.9 ± ± ± 0.3. ± ± 0.4 TABLE I. Substructure parameters from Lévy fits, dσ/dt = A exp( ( t R ) α ) at various ISR colliding energies, for α = 0.9 fixed. The values of the Lévy scales R are rounded up to two decimal digits, corresponding to the estimated resolution. s (TeV) 7 13 A (mb GeV ) 1.5 ± ± 0.57 R (fm) 0.48 ± ± 0.01 χ /NDF 4.93 / 4. / 8 CL [%] σ tot sub σ el sub (mb) 5.46 ± ± 0.60 (mb) 0.5 ± ± 0.14 TABLE II. Substructure parameters from Lévy fits, dσ/dt = A exp( ( t R ) α ) at LHC colliding energies, for α = 0.9 fixed. The fit range is the secondary cone after the dip and the bump, corresponding to 1.5 t 3 GeV at LHC energies. The values of the Lévy scale R and the corresponding contribution of these substructures to the total cross-section are rounded up to two decimal digits, corresponding to the estimated resolution of the method. TeV, the same kind of larger and darker substructure is observed again, with a size that does not grow significantly from 7 TeV to 13 TeV, but has a significantly darker shadow, as shown on Figure 3. The zeroth order Lévy fit parameters that characterize the secondary cone region and the corresponding substructures at various colliding energies are summarized in Tables I,II. The calculated contributions of these substructures to the total and to the elastic cross-section, σ sub tot and σ sub el are also shown, as obtained from eqs. (1,13) with (a i, b i ) = (0, 0). Figure 4 indicates the shadow profiles of the protons at s = 13 TeV, 7 TeV and 6 GeV, 9

10 and compares them with the shadow profiles of the substructures of the protons at the same energies. The shadow profiles of these substructures are significantly smaller and fainter than the shadow of the protons at each investigated energies. With increasing colliding energy, the size and the darkness of the proton is increasing. In the LHC energy range, a nearly constant maximally dark region develops in the central region of the proton shadow, that increases with increasing energies. This structure is qualitatively different from the nearly Gaussian proton structure apparent in the collisions at lower energies. Surprisingly, a substructure is also resolved in the ISR energy range: it is found to be rather faint and its size within the errors of resolution is a constant. At the LHC energies, the imaged substructure is significantly larger and darker as compared to that at the ISR. In summary, we have imaged the overall space-time structure of elastic proton-proton collisions at the largest currently available colliding energy of s = 13 TeV and CERN LHC, using recent preliminary TOTEM data [5]. Uniting R. J. Glauber s multiple diffraction theory with a new orthonormal Lévy series expansion of the forward scattering amplitude, have achieved an unprecendented precision in describing the differential cross section of elastic proton-proton scattering, over orders of magnitude. Surprisingly, we have identified and quantified a substructure of the protons, that contributes to the total cross-section with about 0.3 mb in the energy range of s = ISR energy range and about 5.5 and 8.9 mb at s = 7 and 13 TeV LHC energies. The imaging method presented here can only report on the observation of these substructures but, being model independent, it cannot be overused to determine the detailed physical meaning of these newly identified substructures. Nevertheless, let us mention that these are strikingly similar in magnitude to the size of the dressed quarks and diquarks that have been determined from a model dependent analysis of proton-proton elastic scattering at the ISR and LHC energies [1 16]. The similarity of the model independent Lévy imaging results to the results of this Bialas-Bzdak quark-diquark picture of elastic proton-proton scattering may also provide a quantitative support to the quark-diquark picture of baryons, that were recently obtained by Brodsky and collaborators using AdS/QCD strong coupling calculations [17]. Our model independent imaging results are consistent with this QCD based interpretation as well, but they are open for other possible quantitative models and descriptions to explain the origin of the observed substructures. Lévy stable distributions are ubiquitous in Nature, with applications in probabilistic mod-

11 els, correlated systems and fractals, anomalous diffusion and chaos, physics, radiophysics, astrophysics, stochastic algorithms, quantitative finance, biology and geology [1]. Hence, we may reasonably hope that the model independent method that we have presented here can also be successfully applied in broad areas of quantitative sciences to search for and to quantify the deviations from Lévy distributions for positive definite observables. Acknowledgments T. Cs. thanks to R.J. Glauber for inspiring discussions and for proposing a series expansion to describe elastic scattering. We acknowledge inspiring discussions with S. Giani, W. Guryn, G. Gustafson, L. Lönnblad, K. Österberg and M. Šumbera and thank the TOTEM Collaboration for making their preliminary data set available for us. R. P. is partially supported by the Swedish Research Council grants No and , by CONICYT grants PIA ACT1406 and MEC , by an ERC H00 grant No , as well as by the Ministry of Education, Youth and Sports of the Czech Republic project LT T. Cs. and A. S. were partially supported by the NKIFH grants No. FK-1384 and FK-13959, and by the EFOP grant (Hungary) and the NKM-9/017 projects of the Hungarian and the Ukrainian Academy of Sciences. Our collaboration was supported by THOR, the EU COST Action CA1513. [1] V. V. Uchaikin and V. M. Zolotarev, Chance and Stability, Stable Distributions and Their Applications, (VSP Science, 1999) ISBN: , 596 pp. [] C. Tsallis, S. V. F. Lévy, A. M. C. Souza and R. Maynard, Phys. Rev. Lett. 75, 3589 (1995). [3] T. Csörgő, S. Hegyi and W. A. Zajc, Eur. Phys. J. C 36 (004) 67. [4] T. Csörgő, R. Pasechnik and A. Ster, arxiv: [hep-ph], submitted for a publication. [5] F. Nemes for the TOTEM Collaboration, Proc. of the 4th Elba workshop on Forward LHC energy, 4-6 May 018, Elba, Italy, [6] T. Csörgő and S. Hegyi, Phys. Lett. B 489, 15 (000) [7] M. B. De Kock, H. C. Eggers and T. Csörgő, PoS WPCF 011 (011) 033. [8] T. Novák, T. Csörgő, H. C. Eggers and M. de Kock, Acta Phys. Polon. Supp. 9, 89 (016) [9] G. Antchev et al. [TOTEM Collaboration], Nucl. Phys. B 899, 57 (015) [] R. J. Glauber and J. Velasco, Phys. Lett. 147B, 380 (1984). [11] I. M. Dremin, Universe 4, (5), 65 (018). 11

12 [1] A. Bialas and A. Bzdak, Acta Phys. Polon. B 38, 159 (007). [13] F. Nemes and T. Csörgő, Int. J. Mod. Phys. A 7, (01). [14] T. Csörgő and F. Nemes, Int. J. Mod. Phys. A 9, (014). [15] T. Csörgő, R. J. Glauber and F. Nemes, arxiv: [hep-ph], Proceedings of Low-X 013, (Rehovot, Israel). [16] F. Nemes, T. Csörgő and M. Csanád, Int. J. Mod. Phys. A 30, no. 14, (015). [17] S. J. Brodsky, Adv. High Energy Phys. 018, (018). 1