Convergence properties of L\'evy expansions: implications for Odderon and proton structure

TL;DR

This paper investigates the convergence properties of the Levy expansion method for analyzing elastic hadron-hadron scattering, demonstrating its stability and effectiveness in revealing proton sub-structures and Odderon effects.

Contribution

It introduces a new model-independent Levy expansion technique and shows its convergence and stability in analyzing differential cross sections in hadron scattering.

Findings

Levy expansion converges and is stable for analyzing scattering data.

Reveals sub-structures inside protons.

Detects Odderon effects in proton scattering.

Abstract

We detail here the convergence properties of a new model-independent imaging method, the L\'evy expansion, that seems to play an important role in the analysis of the differential cross section of elastic hadron-hadron scattering. We demonstrate, how our earlier results concerning the Odderon effects in the differential cross-section of elastic proton-proton and proton-antiproton scattering as well as those related to apparent sub-structures inside the protons were obtained in a convergent and stable manner.

Equations

$$\frac{d\sigma }{dt}=\frac{1}{4\pi }|T_{\mathrm{el}}(\Delta ){|}^2,T_{\mathrm{el}}(\Delta )$$

$${\int }_0^{\infty }dz\exp (-z^{\alpha })l_n(z|\alpha )l_m(z|\alpha )$$

$$l_j(z|\alpha )$$

$$L_0(z|\alpha )$$

$$L_1(z|\alpha )$$

$$L_2(z|\alpha )$$

$$D_0(\alpha )$$

$$D_1(\alpha )$$

$$D_2(\alpha )$$

$$D_3(\alpha )$$

$${\mu }_{n,\alpha }={\int }_0^{\infty }dz{z}^n\exp (-z^{\alpha })=\frac{1}{\alpha }\Gamma \left(\frac{n+1}{\alpha }\right),\Gamma (x)={\int }_0^{\infty }dz{z}^{x-1}{e}^{-z},$$

$${\sigma }_{\mathrm{tot}}\equiv 2\mathrm{Im}{T}_{\mathrm{el}}(\Delta =0)=2\sqrt{4\pi A}\left(1+\sum _{i=1}^{\infty }{a}_i{l}_i(0|\alpha )\right).$$

$${\sigma }_{\mathrm{el}}={\int }_{-t=0}^{\infty }dt\frac{d\sigma }{dt}=\frac{A}{R^2}\left[\frac{1}{\alpha }\Gamma \left(\frac{1}{\alpha }\right)+\sum _{i=1}^{\infty }(a_i^2+b_i^2)\right].$$

$$B(t)=\frac{d}{dt}\left(\ln \frac{d\sigma }{dt}\right).$$

$$\rho (t)\equiv \frac{\mathrm{Re}{T}_{\mathrm{el}}(t)}{\mathrm{Im}{T}_{\mathrm{el}}(t)}=-{\frac{{\sum }_{i=1}^{\infty }{b}_i{l}_i(z|\alpha )}{1+{\sum }_{i=1}^{\infty }{a}_i{l}_i(z|\alpha )}|}_{z=t{R}^2}.$$

$$T_{\mathrm{el}}(t)=|T_{\mathrm{el}}(t)|\exp \left[i{\varphi }_1(t)\right].$$

$${\varphi }_2(t)=\frac{\pi }{2}-\arctan \rho (t).$$

$$P(b)\equiv 1-{\left|e^{-\Omega (b)}\right|}^2=[2-\mathrm{Im}{t}_{\mathrm{el}}(b)]\mathrm{Im}{t}_{\mathrm{el}}(b)-[\mathrm{Re}{t}_{\mathrm{el}}(b){]}^2,$$

$$t_{\mathrm{el}}(b)=\int \frac{d^2\Delta }{(2\pi {)}^2}{e}^{-i\mathbf{\Delta }\mathbf{b}}{T}_{\mathrm{el}}(\Delta )=i\left[1-e^{-\Omega (b)}\right],\Delta \equiv |\mathbf{\Delta }|,b\equiv |\mathbf{b}|.$$

Full text

1 Introduction

The model-independent Lévy imaging technique  Csorgo:2018uyp ; Csorgo:2018ruk ; Csorgo:2019rsr has recently become a useful, simple and unambiguous tool for extracting the physics information from the elastic hadron-hadron scattering data in a statistically acceptable manner. The power of this technique has been first demonstrated in our earlier analysis of the most recent data sets from the total, elastic and differential cross-section measurements in elastic $pp$ collisions $\sqrt{s}=13$ TeV performed by the TOTEM Collaboration at the Large Hadron Collider (LHC) (for a few most recent TOTEM publications, see Refs. Antchev:2017dia ; Antchev:2017yns ; Antchev:2018edk ; Antchev:2018rec ). These TOTEM results, and in particular the comparison of the differential cross-section of elastic proton-proton scattering at $\sqrt{s}=2.76$ TeV with D0 results on elastic proton-antiproton scattering at 1.96 TeV Abazov:2012qb indicate several Odderon effects, as discussed recently in Refs. Csorgo:2018uyp ; Antchev:2018rec

In particular, indirect signatures of the Odderon exchange in differential elastic $pp$ and $p\bar{p}$ cross-sections have been identified by using the Lévy imaging technique, also known as the model-independent Lévy expansion method. Another important implication of this technique is that it enables to probe the internal structure of the proton by identifying its smaller substructures imprinted in the behaviour of the $t$-dependent elastic slope $B(t)$. In particular, the proton substructures of two distinct sizes in the low (a few tens of GeV) and high (a few TeV) energy regimes, respectively, have been found and discussed in Refs. (Csorgo:2018uyp, ; Csorgo:2018ruk, ; Csorgo:2019rsr, ). A remarkable feature of the Lévy expansion of the elastic amplitude is that the diffractive cone is described fairly well by the Lévy-stable (or stretched exponential) distribution in terms of two free parameters only, the Lévy scale parameter $R$ characterising the length-scale of the scattered systems, and the exponent $\alpha=\alpha_{L}/2$. One of the conventions of elastic scattering is to parameterize the diffractive cone with an exponential behaviour, corresponding to $\alpha=1$, i.e. to a Gaussian scattering amplitude in the impact parameter space. The exponent $\alpha_{L}=2\alpha$ is known as the Lévy index of stability providing a small but significant deviation of the cone shape from the classically expected Gaussian shape. A detailed presentation of this sophisticated new imaging technique is detailed in Ref. Csorgo:2018uyp , while for a brief summary of the main results, see Refs. Csorgo:2018ruk ; Csorgo:2019rsr .

Given a power of the Lévy imaging technique operating with very few initial assumption for description of a large amount of data, a natural question arises about the stability and convergence of the associated Lévy series that are used for mapping the elastic amplitude with this method. In this contribution, we attempt to give a short description of the main results of the method, as well as demonstrate its convergence properties for a large variety of data sets at different scattering energies. The main results were summarized recently in two short contributions Csorgo:2018ruk ; Csorgo:2019rsr , however the convergence properties of this model-independent Lévy expansion method were not yet detailed in the literature before.

2 Model-independent Lévy expansion of the scattering amplitude

In order to study systematically the deviations of a given data set on the elastic cross section $d\sigma/dt$, differential in four-momentum transfer $t=(p_{1}-p_{3})^{2}<0$, from an approximate Lévy-stable shape apparent at larger $t$ beyond the diffractive cone, we adopt the Lévy series expansion for the elastic amplitude $T_{\rm el}(\Delta)$, $\Delta=\sqrt{|t|}$, represented in terms of a complete orthonormal set of Lévy polynomials (with complex coefficients) and a Lévy weight function $w(z|\alpha)=\exp(-z^{\alpha})$ Csorgo:2018uyp ; Csorgo:2018ruk . While the Lévy polynomials exhibit an oscillatory behaviour in dimensionless scaling variable $z=|t|R^{2}$, the differential cross-section of elastic scattering is proportional to a hit distribution which is, by construction, a positively-definite function of $t$, i.e.

[TABLE]

where $c_{j}=a_{j}+ib_{j}$ are the complex coefficients of the Lévy expansion ($j=0,1,...,)$ and $a_{0}=1$ fixed so that the overall normalization can be absorbed into the coefficient $A$. This simple expansion form indicates three underlying physical assumptions, with further details explained below.

The leading order behaviour corresponds to a nearly exponential (in other words, a non-exponential) distribution. This leading order behaviour is $\frac{d\sigma}{dt}=A\exp\left[-(R^{2}|t|)^{\alpha}\right]$, which is consistent with the TOTEM observation of a non-exponential low-$|t|$ behaviour of the differential cross-section of elastic $pp$ scattering at $\sqrt{s}=8$ TeV  Antchev:2015zza . In the limit of $\alpha$ $\rightarrow$ 1, the exponential cone behaviour is recovered, with $\lim_{\alpha\rightarrow 1}\frac{d\sigma}{dt}=A\exp\left[-B|t|\right]$ with $B=R^{2}$. 2. 2.

The leading order scattering amplitude is assumed to be imaginary, with vanishing real part, corresponding to $b_{0}=0$ and to all the $c_{j}$ expansion coefficients vanish for $j\geq 1$. Such an assumption can be relaxed by assuming that $b_{0}$ can be different from zero, but so far all the fits that we have performed, we have found $b_{0}$ to be consistent within errors with zero, so we have fixed this possible expansion parameter to zero, accordingly. The possibility of relaxing $b_{0}$ to be different from zero is investigated in Fig. 1 and it is briefly discussed there in the context of the $\sqrt{s}=13$ TeV elastic $pp$ scattering data. 3. 3.

For predominantly imaginary scattering amplitudes with $b_{0}=0$, the existence of a non-vanishing real to imaginary ratio $\rho(t)$, in particular, $\rho_{0}=\rho(t=0)$, is mapped one-to-one to the existence of a diffractive interference structure, i.e. a diffractive minimum in elastic scattering cross section. Given that for $t\neq 0$ the Lévy expansion is an analytic function, and that at the first diffractive minimum the imaginary part is expected to vanish, the phase of the scattering amplitude can be uniquely determined with the Lévy expansion method up to an overall sign, that can be fixed from measurements at the Coulomb-nuclear interference region. This implies that the $\rho_{0}$ parameter can be determined from the Lévy expansions if the fit range includes a diffractive minimum.

The set of orthonormal polynomials are denoted above by $\left\{l_{j}(z|\alpha)\right\}_{j=0}^{\infty}$, while the set of orthogonal, but unnormalized Lévy polynomials are denoted by $\left\{L_{j}(z|\alpha)\right\}_{j=0}^{\infty}$. These unnormalized Lévy polynomials were introduced in Refs. DeKock:2012gp ; Novak:2016cyc , to analyze nearly Lévy shaped Bose-Einstein correlations in two-particle Bose-Einstein correlations or particle interferometry. Their orthogonality is defined with respect to the weight function $w(z|\alpha)=\exp\left(-|z|^{\alpha}\right)$ where the dimensionless variable $z=|t|R^{2}$. The weight function acts simultaneously as the leading order approximate form of the measured distribution, and at the same time also as a measure in an abstract Hilbert space, where the convergence properties of the Lévy series can be investigated, following the general ideas proposed in Ref. Csorgo:1999wx . Thus, this expansion method allows to keep the number of expansion coefficients at the minimal necessary level.

The orthonormality of the Lévy polynomials $l_{n}(z\,|\,\alpha)$ is expressed as

[TABLE]

These orthonormalized polynomials $l_{j}(z\,|\alpha)$ are proportional to the orthogonal – but unnormalized — Lévy polynomials, $L_{j}(z\,|\,\alpha)$,

[TABLE]

which in turn are given by a Gram-Schmidt orthogonalization procedure as

[TABLE]

introduced previously in Ref. DeKock:2012gp . The Gram-determinants of order $j$, $D_{j}\equiv D_{j}(\alpha)$ defined as

[TABLE]

In the above expressions,

[TABLE]

where $\Gamma(x)$ the Euler’s gamma function.

3 Observables

Observables can be calculated from Eq. (1). In what follows, we present the results for the fixed $b_{0}=0$.

The total cross-section is obtained, using the optical theorem, as follows

[TABLE]

The elastic cross-section, $\sigma_{\rm el}$ can be obtained using the orthonormality, Eq. (2) as

[TABLE]

The following $t$-dependent functions have also been analyzed:

• •

The four-momentum transfer dependent elastic slope $B(t)$, defined as

[TABLE]

• •

The four-momentum dependent $\rho(t)$, defined as the ratio of the real to imaginary parts of $T_{\rm el}$

[TABLE]

The value of $\rho(t)$ at $t=0$ can be measured in the Coulomb-nuclear interference region and in this work we refer to this value as $\rho_{0}=\rho(t=0)$.

• •

The first definition for the nuclear phase $\phi_{1}(t)$, that can be introduced as

[TABLE]

An alternative definition was used recently by the TOTEM collaboration, corresponding to the principal value of the nuclear phase, that reads as

[TABLE]

If the nuclear phase $\phi_{1}(t)$ satisfies $0\leq\phi_{1}(t)\leq\pi$, then the above two definitions are equivalent. However, for complex arguments, $\arctan(z)$ has branch cut discontinuities on the complex plane hence in general the two definitions $\phi_{1}(t)$ and $\phi_{2}(t)$ are inequivalent, as detailed recently in Ref. Csorgo:2018uyp . In the present work we plot the principal value of the nuclear phase, that corresponds to the definition with Eq. (27), which by definition satisfies $0\leq\phi_{2}(t)\leq\pi$.

• •

We shall also study the convergence properties of the shadow profile function,

[TABLE]

where $\Omega(b)$ is the complex-valued opacity function. These functions can be represented in terms of the impact-parameter dependent elastic amplitude $t_{\rm el}(b)$, see for example Ref. Csorgo:2018uyp :

[TABLE]

4 Convergence properties

The Lévy series (1) enables one to precisely and model-independently characterise the $t$-dependence of the elastic cross-section $\frac{d\sigma}{dt}$ not only in the diffractive cone, but also significantly away from it, namely, in a vicinity of the diffractive minimum as well as beyond it in the large $t$ regime (diffractive tail). This is achieved by a fit of Eq. (1) to the existing data controlling the fit quality and the convergence of the Lévy ansatz to the data points in each particular region of 4-momentum transfers.

We adopted the following procedure: a fourth order Lévy polynomial fit was performed to the measured $\frac{d\sigma}{dt}$ data for elastic $pp$ scattering. The expansion coefficients $(a_{i},b_{i})$ are determined for $i=1,...,i_{\rm max}=4$. On the subsequent plots, the contributions are summed only up to a given order $j\leq i_{\rm max}$. Thus a procedure not unlike to a Taylor series expansion is defined and the inclusion of each subsequent order enables us to see how the series converges to the measurements on a more and more extended domain. As described in Ref. Csorgo:2018uyp , good quality fits of elastic $pp$ scatterings were achieved for almost all of the published data sets, using $i_{\rm max}=4$. To fit the $p\bar{p}$ datasets, that were measured in a more limited $t$-range and with reduced statistics, third order Lévy expansions with $i_{\rm max}=3$ were sufficient.

As we demonstrate below, these fits converge rather fast to the data when increasing the order of the Lévy expansion from zeroth to the fourth order.

As an example of the convergence property, in Fig. 2, we illustrate partial sums contributions to the Lévy series (1) for the differential elastic cross-section against the TOTEM 7 TeV (left) and D0 1.96 TeV (right) data sets. As Figure 2 indicates, the zeroth order term fixes the total cross-section, i. e. the contribution at $t=0$ only, it apparently insufficient even at small $t$ to approximate the shape of the distribution. The first order Lévy expansions is required to describe $B(t=0)$ and the beginining of the diffractive cone. Although the first order contribution exhibits a structure with a diffractive dip, the position and the size of the dip is not yet picked up correctly. The second order partial sum is needed for getting the position of the diffractive minimum correctly, while the third order term fixes the position and the magnitude of the diffractive maximum. The fourth order terms are necessary to describe the data also well beyond the diffractive minimum and maximum.

While at zeroth order the real part of the elastic amplitude is found to be vanishing in our current approach, it may get (re)generated starting from the first order partial sums, and appears to be necessarily small. Remarkably, the considered quality fits shown for example in Fig. 1 for the 13 TeV elastic $pp$ scattering data of TOTEM enable us to reproduce the measured real-to-imaginary parts ratio $\rho(t)$ at $t=0$ with an excellent precision.

In Figs. 3, 4 and 5 we show the convergence properties of the shadow profile (left panels) and the elastic slope (right panels) for Lévy fits of the TOTEM data at 13 TeV, 7 TeV ($pp$ collisions) and D0 data at 1.96 TeV ($p\bar{p}$ collisions) at different consecutive orders as previously illustrated in Fig. 2, respectively.

The convergence of partial sums contributions to the shadow profile turns out to be faster than that for the elastic slope. Namely, already starting from the second Lévy order the result for $P(b)$ remains fairly stable at large $b$ for all the considered data sets, but higher order terms are needed for a precise result at small values of $b$. For $pp$ collisions TOTEM data at 13 and 7 TeV, the behavior of the elastic slope stabilizes in a vicinity of the diffractive dip/bump structure only at third Lévy order, in consistency with the observations made from Fig. 2. So, for all TOTEM data sets the minimal preferred order enabling to extend the data description significantly beyond the dip/bump structure is the fourth. Due to a lack of data in the large-$t$ region for $p\bar{p}$ collisions at 1.96 TeV, we limit ourselves to the third order that provides the results for $B(t)$ sufficiently stable in the “shoulder” region.

Finally, in Fig. 6 we demonstrate the corresponding convergence property of $\phi_{2}(t)$, the principal value of the nuclear phase for $\sqrt{s}=13$ TeV $pp$ elastic scattering. We observe that there is no significant difference between the third- and the fourth-order Lévy expansion results indicating a good stability property of this function in the range of the diffractive minimum, where $\phi_{2}(t_{\rm dip})=\pi$.

5 Odderon

The concept of Odderon exchange corresponds to a crossing-odd term in the $pp$ elastic scattering amplitude. This concept was introduced by Lukaszuk and Nicolescu in 1973 Lukaszuk:1973nt .

The recent TOTEM results Antchev:2017dia ; Antchev:2017yns generated a burst of high level and intense theoretical debate about the correct interpretation of these data, see Refs. Samokhin:2017kde ; Khoze:2017swe ; Petrov:2018xma ; Khoze:2018bus ; Goncalves:2018yxc ; Shabelski:2018jfq ; Broilo:2018els ; Lebiedowicz:2018eui ; Martynov:2018nyb ; Troshin:2018ihb ; Dremin:2018uwt ; Broniowski:2018xbg ; Khoze:2018kna . All possible extreme views were present among these first interpretative papers, including claims for a maximal Odderon effect Martynov:2018nyb and claims of lack of any significant Odderon effects, Refs. Shabelski:2018jfq ; Khoze:2018kna . A review of recent theoretical developments on possible Odderon effects is given in Refs. Martynov:2018pye ; Martynov:2018yas .

At sufficiently high energies, the relative contribution from secondary Regge trajectories is suppressed, as they decay as negative powers of the colliding energy $\sqrt{s}$. The vanishing nature of these Reggeon contributions offers a direct way of extracting the Odderon as well as the Pomeron contributions, $T_{el}^{O}(s,t)$ and $T_{el}^{P}(s,t)$, respectively, from elastic scattering data at sufficiently high colliding energies. Thus, the Odderon effects are expected to be detected clearly with measurements in the TeV energy range Avila:2006wy ; Martynov:2017zjz .

In Ref. Ster:2015esa , the authors also argued that the LHC energy scale is already sufficiently large to suppress the Reggeon contributions, and they presented the $(s,t)$ dependent contributions of Odderon exchange to the differential and total cross-sections at LHC energies. That analysis relied on a model-dependent, phenomenological extension of the Phillips-Barger model Phillips:1974vt , that focused on fitting the dip region of elastic proton-proton scattering, but did not analyze in detail the tails and the cone region. The fitted model parameters of proton-proton and proton-antiproton reactions were extrapolated to exactly the same energies, and the results were recently confirmed and extended in Ref. Goncalves:2018nsp . Similarly, Ref. Lebiedowicz:2018eui also argued that the currently highest LHC energy of $\sqrt{s}=$ 13 TeV is sufficiently high to see various Odderon contributions. In particular, the Pomeron and the Odderon contributions can be extracted from the forward scattering amplitudes at sufficiently high energies as discussed, for example, in Refs. Ster:2015esa ; Csorgo:2018uyp . Elastic proton-proton and proton-antiproton scattering data were not measured at the same energy in the TeV region so far. However, we have identified two robust-looking features of the already performed measurements, that provided not only an Odderon signal, but they also clearly indicated the existence of two different sizes for some sub-structures inside the protons, as imaged by elastic proton-proton scattering.

In particular, we found that clear, but indirect signals of Odderon effects are present in the difference between the $t$-dependent nuclear slopes of elastic proton-proton and proton-antiproton scattering Csorgo:2018uyp ; Csorgo:2019rsr . Our results for the existence of a well-defined and negative minimum of the $B(t)$ functions for $pp$ reactions and a lack of significantly negative values of $B(t)$ in $p\overline{p}$ reactions at the TeV region has recently been confirmed using the maximal Odderon model of Martynov and Nicolescu in Ref. Martynov:2018sga . From the experimental point of view, $B(t)$ is straightforward to measure and can be used for an experimental search for Odderon effects independently of Lévy expansion and imaging results Csorgo:2019fbf .

In addition, we have also identified a clear difference between the principal values of the nuclear phase, $\phi_{2}(t)$ of proton-proton and proton-antiproton collisions in the TeV energy range Csorgo:2018uyp ; Csorgo:2019rsr . However, the $t$-dependence of the nuclear phase is rather difficult, close to impossible to access experimentally, in particular independently of the Lévy imaging methods.

Both of these Odderon effects are obtained with a convergent series expansion, and are stable for higher order Lévy expansion coefficients, as indicated on the right panels of Figs. 3, 4 and 5 in case of the nuclear slope $B(t)$ and on Fig. 6 in case of the principal value of the nuclear phase $\phi_{2}(t)$.

6 Summary and conclusions

Our analysis in Refs. Csorgo:2018uyp ; Csorgo:2018ruk ; Csorgo:2019rsr has been primarily motivated by the search for Odderon effects. We have identified two independent Odderon effects in TOTEM differential cross-section measurements. The comparision of $B(t)$ for $pp$ and $p\overline{p}$ reactions at exactly the same $\sqrt{s}$ in the TeV region is one of the most promising channel for the experimental observation of Odderon effects.

One of the most obvious but nevertheless striking feature of the elastic $pp$ scattering at TeV energies is that the differential cross-section has a unique, single minimum. In multiple diffractive scattering theory, single diffractive minimum may appear in symmetric collisions of composite objects if and only if the colliding systems have two internal substructures Czyz:1969jg . This suggests that the quark-diquark picture of elastic proton-proton collisions, where a diquark that acts as a single unit in elastic scattering even at the LHC energies, formulated in terms of the real extended Bialas-Bzdak model of elastic proton-proton scattering Bialas:2006qf ; Nemes:2015iia , may indeed capture correctly some of the most fundamental properties of elastic proton-proton collisions at the LHC energies.

One of our most surprising result was a clear-cut evidence for two different sub-structures inside the protons, as detailed in Ref. Csorgo:2018uyp . We determined the significance of these substructure effects and estimated the sizes of these sub-structures and their contributions to the total and elastic proton-proton cross-section in Ref. Csorgo:2018ruk .

Acknowledgments

We would like to thank S. Giani, G. Gustafson, V. Khoze, F. Nemes, B. Nicolescu and K. Österberg for inspiring and clarifying discussions. T. Cs. expresses his gratitude to the Organizers of ISMD 2018 for the opportunity to present these results, for partial support and for an outstanding, inspiring and useful meeting. T. Cs. was partially supported by the Hungarian NKIFH grants FK-123842 and FK-123959 and the EFOP-3.6.1-16-2016-00001 grants (Hungary). R.P. was partially supported by the Swedish Research Council grants, contract numbers 621-2013-4287 and 2016-05996, CONICYT grant MEC80170112, and by the Ministry of Education, Youth and Sports of the Czech Republic project LT17018. This research was partially supported by the THOR project, COST Action CA15213 of the European Union.