1. Introduction
In 1973, Lukaszuk and Nicolescu [
1] proposed that a noticeable crossing-odd contribution called odderon may be present in the amplitude of elastic proton–proton (
) and proton–antiproton (
) scattering at asymptotically high energies. In the field theory of strong interactions, quantum chromodynamics (QCD), odderon exchange corresponds to the
t-channel exchange of a colour-neutral gluonic compound state consisting of an odd number of gluons [
2]. For more than 20 years since the theoretical prediction of the odderon in the framework of QCD, and for more than 47 years since the odderon concept was introduced in Regge phenomenology, the odderon has remained elusive until the time of 2019–2021, due to the lack of definitive and statistically significant experimental evidence of odderon exchange.
A direct way to probe odderon exchange in elastic scattering is by comparing the differential cross-section of particle–particle and particle–antiparticle scattering at the same and sufficiently high energy [
3,
4]. The TeV energy scale is high enough to render the amplitudes for the massive (Reggeon) exchange to negligibly small values. The differential cross-sections for elastic proton–proton and proton–antiproton scatterings are proportional to the modulus squares of the corresponding amplitudes. If the modulus squared amplitudes are different, then the two amplitudes cannot be equal. As the amplitude for odderon exchange is proportional to the difference between the amplitudes of the proton–antiproton and the proton–proton elastic scattering, it follows that the amplitude for odderon exchange cannot vanish if the differential cross-sections of elastic proton–proton and proton–antiproton scatterings are significantly different at the TeV energy scale, as detailed in ref. [
5]. A search performed at the ISR energy of
GeV in 1985 [
6] resulted in an indication of the odderon at the 3.35
significance level. That analysis, however, did not utilize all the available data in the overlapping acceptance of the
and
measurements. Furthermore, at
GeV, Reggeon exchanges may play a significant role, rendering the odderon search at the ISR energies rather inconclusive.
As far as we know, the first anonymously peer-reviewed publication of a statistically significant, at least 6.26
signal of odderon exchange was published in June 2020 by the authors of the present manuscript, in the proceedings of the ISMD 2019 (Santa Fe, NM, USA) [
7]. This peer-reviewed conference contribution was backed up in our February 2021 paper, with the same statistical significance of an at least 6.26
signal for odderon exchange [
5]. Our Hungarian–Swedish team introduced a new scaling function that turned out to be energy independent in the LHC energies of
= 2.76–7.0 TeV in elastic proton–proton (
) collisions, based on a model-independent, direct data-to-data comparison [
5].
Our idea was based on the existence of a diffraction cone at small values of the squared four-momentum transfer, the Mandelstam variable t, characterized by an exponential fall. By plotting the and data as a function of a dimensionless variable , where B is the characteristic scale of this exponential fall, and by dividing the differential cross-section by , where is the integral of the differential cross-section (dominated by the contribution of the diffractive cone), a dimensionless scaling function is obtained. In the region of the diffractive cone, this scaling function is approximately a structureless exponential, . As the diffractive cone is located at small values of t, this scaling removes the trivial s-dependences of the differential cross-section, where the Mandelstam variable s is the squared center of mass energy.
Thus, such an
scaling removes the trivial energy-dependent terms, due to the known
s-dependence of the elastic slope
, the elastic and total cross-sections
and
, and the real-to-imaginary ratio
[
8].
It follows that the amplitude for odderon exchange cannot vanish if the dimensionless
scaling functions of elastic proton–proton and proton–antiproton scatterings are significantly different at the TeV energy scale [
5].
These results as well as the model-
dependent investigation of the domain of validity of the
scaling have been seconded by a theoretical paper of T. Csörgő and I. Szanyi [
9], increasing the statistical significance of the observation of odderon exchange to at least 7.08
. At the same time, this model-dependent investigation found that the domain of validity in
s of the
scaling (where
stands for a fixed energy scale in the TeV region, e.g.,
TeV) extends also to the top Tevatron energies of
= 1.96 TeV [
9]. This theoretical work utilized a validated model, proposed in its first form by A. Bialas and A. Bzdak [
10]; however, the original model lacked a real part hence the possibility of odderon exchange. However, the so-called Real-extended Bialas–Bzdak (ReBB) model of ref. [
11] fixed these shortcomings and has been utilized in ref. [
9] to extrapolate not only the elastic proton–proton scattering data from the LHC energies of
and 7 TeV to the D0 energy of
TeV but also to extrapolate the elastic proton–antiproton scattering data from
and 1.96 TeV of the UA4 [
12,
13], and the D0 Collaborations [
14] to the LHC energies of 2.76 TeV and 7 TeV. Evaluating the proton–proton data with a model increased the uncertainty and decreased the odderon signal from proton–proton scattering data alone, but this decrease was well over-compensated with the ability of the model to evaluate theoretically the proton–antiproton scattering at all the LHC energies. Overall, this procedure resulted in a model-dependent increase in the statistical significance from odderon exchange from 6.26 to 7.08
[
9] as published in final form in July 2021, but limited the comparison to the diffractive minimum and maximum region in the four-momentum transfer range, to the verified domain of validity of the ReBB model. More recently, these results were extended to the new TOTEM data on elastic
scattering at
TeV, published in March 2022 [
15] in ref. [
9]. When TOTEM data on elastic
collisions at
, 7, and 2.76 TeV are analysed simultaneously with D0 data at
1.96 TeV in the framework of the ReBB model, a combined statistical significance greater than 32.4
can be achieved, rendering the statistical significance of odderon observation, in any practical terms, to a certainty [
9].
A series of important papers has been published by the TOTEM Collaboration investigating the properties of elastic
scattering in the LHC energy range between
and 13 TeV in refs. [
16,
17,
18,
19]. Most recently, the latest measurement performed by TOTEM at
TeV [
15] extended the earlier analysis up to
-values of 1.9 GeV
2 [
20]. An increase in the total cross-section,
, associated with a decrease in the real-to-imaginary ratio,
, with energy, first identified at
TeV [
16,
17] indicated a possible odderon effect triggering an intense discussion in the literature (see, e.g., refs. [
21,
22,
23,
24,
25,
26,
27,
28,
29,
30,
31,
32,
33,
34,
35,
36,
37,
38]). The persistent diffractive minimum–maximum structure in the
t-dependent profile of
in elastic
collisions observed by TOTEM at
= 2.76, 7, 8 and 13 TeV, and the lack of such structure in elastic
collisions measured by D0 [
14] at
TeV indicate a qualitatively clear odderon effect. The possibility of utilizing experimental data at the TeV energy scale for a search for a statistically significant odderon exchange has been proposed by Jenkovszky et al. in refs. [
4,
23]. In 2020, the TOTEM Collaboration made a qualitative conclusion about odderon exchange in ref. [
19] as follows:
“Under the condition that the effects due to the energy difference between TOTEM and D0 can be neglected, the result provides evidence for a colourless 3-gluon bound state exchange in the t-channel of the elastic scattering”. However, no statistical significance for this observation has been evaluated in ref. [
19].
More recently, in August 2021, a properly quantified statistical significance of the odderon signal has been published by the TOTEM and D0 Collaborations [
39] employing different methods and techniques, obtaining an at least 5.2
combined statistical significance for an almost model-independent observation of odderon exchange, a first statistically significant result obtained by two experimental collaborations. This result was based on the extrapolation of TOTEM experimental data of the differential cross-section of elastic
scattering from
13, 8, 7 and 2.76 TeV down to 1.96 TeV using an almost model-independent technique and comparing the results with D0 data in a limited four-momentum transfer range, resulting in an at least 3.4
signal for odderon exchange. TOTEM has also measured the pair of the total cross-section and the real-to-imaginary ratio
, and compared it with a set of models without odderon exchange. When a partial combination of the TOTEM
measurements is carried out at
TeV, the obtained partial significances range between 3.4 and 4.6
for the considered models. The full combination of the signal of odderon exchange from TOTEM
measurements at
TeV with the signal of the comparison of extrapolated TOTEM
data to
TeV with a subset of 8 out of 17 D0 data points [
14] on elastic
scattering leads to total significances ranging from 5.2 to 5.7
for
t-channel odderon exchange for each of the considered models [
39].
The validity of this D0–TOTEM proof of odderon exchange has been questioned in several published papers by now. Most importantly, the ATLAS collaboration [
40] published a significantly different total cross-section hence a significantly different pair of
at
TeV, questioning the significance of the signal of odderon exchange from these low
observations. Such a significant incompatibility between the ATLAS and TOTEM measurements of total cross-section and the ratio of the real to imaginary parts of the scattering amplitude, that is, between the pairs of
at
TeV as published by the ATLAS and by the TOTEM experiments, has recently been confirmed by the recent results of Petrov and Tkachenko in refs. [
41,
42,
43]. Furthermore, Donnachie and Landshoff [
44,
45] stressed the point that the phase of an elastic scattering amplitude is related to its energy variation and, as a consequence, they have questioned the validity of the D0–TOTEM signal of odderon exchange at
. Petrov and Tkachenko obtained results similar to those of Donnachie and Landshoff in ref. [
41], suggesting that the systematic error on the determination of the
parameter at
TeV might have been slightly but significantly underestimated by TOTEM in ref. [
17]. Let us mention that refs. [
5,
7] scale out the
observables from their analysis of the
scaling functions, while the low-
t domain has been explicitly excluded from finding a statistically significant signal of odderon exchange in refs. [
9,
46]. Hence, these odderon discovery papers are not affected by the above-mentioned criticisms of the D0–TOTEM publication of a statistically significant, at least 5.2
experimental observation of odderon exchange.
Let us note that the ATLAS data that were measured in the diffraction cone at low
at
7, 8 and 13 TeV in refs. [
40,
47,
48], respectively, collapse to the same
scaling function, that is, close to a structureless
. Hence, the method discussed in this manuscript cannot be utilized for a search of odderon effects in these ATLAS datasets.
In addition to the criticism of the D0–TOTEM method of using
data at
= 13 TeV for the observation of
t-channel odderon exchange, Cui and collaborators [
49] utilized a mathematical approach based on interpolation via continued fractions enhanced by statistical sampling and suggested that a model-independent extrapolation of TOTEM experimental data of the differential cross-section of elastic
scattering from
13, 8, 7 and 2.76 TeV down to 1.96 TeV, comparing the results with D0 data in a limited four-momentum transfer range, results in only an at least 2.2
signal for odderon exchange. This result alone decreases the significance of the D0–TOTEM combined result for odderon exchange from an at least 5.2
to an at least 4.0
signal for odderon exchange [
49], suggesting that the D0–TOTEM method of proving the significance of odderon exchange is only on the level of an indication (defined as a significance between 3.0
and 5.0
), but falls a little bit too short of being an experimentally conclusive, definitive discovery as the corrected value falls short of the discovery threshold of 5
. Such an at least 2.2
signal for odderon exchange from extrapolating the TOTEM measured differential cross-sections of
8, 7 and 2.76 TeV down to 1.96 TeV confirms the model-dependent results of refs. [
9,
46] as well.
A response to these published criticisms was given by the talk of K. Österberg, the physics coordinator of the TOTEM experiment at the ISMD 2023 conference in Gyöngyös, Hungary extending an earlier D0–TOTEM response also given by him at the 2021 Low-x Workshop at La Biodola, Italy [
50,
51]. We have good reasons to expect that a detailed D0–TOTEM response to the above criticisms will be submitted for a publication as soon as reasonably possible. Furthermore, we also second the suggestion of Petrov and Tkachenko in ref. [
43], calling for a joint ATLAS–TOTEM analysis to sort out the differences between their low-
t measurements at
TeV, also proposing a comparison of ATLAS and TOTEM data at low
at
TeV [
47,
52,
53] and also the same comparison at 8 TeV [
48,
54].
Let us note that this ongoing debate in the literature focuses on questioning the validity of certain D0–TOTEM proofs of a statistically significant observation of odderon exchange, but this debate does not question the existence and statistical significance of odderon exchange in all the four published papers on this topic at this energy scale. The statistically significant, well above the 5.0
observations of a
t-channel odderon exchange, as published in refs. [
5,
7] as well as in refs. [
9,
46], have not been affected by the above criticism and have not been challenged so far in other publications either, as far as we know. Furthermore, the statistical significance of odderon exchange as determined from the ReBB model analysis has been increased by taking into account the new 8 TeV data points of the TOTEM experiment by I. Szanyi and T. Csörgő: in any practical terms, within the framework of the ReBB model, the signal for odderon exchange in the limited
GeV
2 and
TeV kinematic region is so large that it amounts to not a probability but a certainty [
9].
In the present manuscript, we summarize our
model-independent analysis of the statistical significance of the odderon observation using the recently published [
15] and extended [
39]
TeV dataset of the TOTEM Collaboration in elastic
collisions, together with earlier data from the D0 [
14] and TOTEM [
18,
19,
52] Collaborations, extending our earlier scaling studies of the differential
scattering cross-section at TeV energies up to 8 TeV.
Our approach is model-independent in the sense that it does not rely on any fitting function or any theoretical input; it uses only linear and log-linear interpolation techniques between neighbouring data points, to allow for data-to-data comparison at the same values of the horizontal axes (using the scaling variable
). As we compare pairwise the
scaling functions constructed at different energies based only upon the available data and look for statistically significant differences within any pair of TeV-scale
and
datasets depending on the collision energy, we need to utilize rebinning that also includes model-independent linear interpolation methods, as detailed in ref. [
5], where the basic concepts and methodology have been explained in detail.
Note that we have determined the domain of validity of the applicability of the
scaling at
TeV so far model-dependently only, in ref. [
5], based on the ReBB model of refs. [
9,
11,
46]. The domain of validity of this
scaling has been found to include
,
, 7 and 8 TeV in a model-independent way as well, as presented at various conferences, e.g., [
55], but these results go way beyond the scope of the present manuscript and will be detailed elsewhere.
2. Scaling of 2.76, 7 and 8 TeV Data of TOTEM
Our analysis is based on our recent discovery of a novel scaling law of elastic
scattering at TeV energies, referred to as
scaling, as described in ref. [
5]. The scaling function
is defined for elastic proton–proton (
) collisions as follows:
A similar function
can be introduced for elastic proton–antiproton (
collisions. In ref. [
5], we have shown that
is independent of the energy in elastic
collisions, utilizing a direct data-to-data comparison of the elastic
data of the TOTEM Collaboration at
and 7 TeV, as described in refs. [
19,
52,
56]. In ref. [
5] we have also shown, in a model-dependent way, that this scaling extends down to
TeV, and we have also found that the scaling violations are already significant at
TeV. In the present study, we investigate whether the energy independence of this scaling is, within experimental errors, valid for the recently published TOTEM data at
TeV or not, within the experimental errors as determined by the TOTEM Collaboration in ref. [
15]. This question is particularly interesting as the elastic-to-total cross-section ratio increases with increasing energies, and at the LHC energies between 7 and 13 TeV it crosses significantly the important limit of
, see Figure 4 of ref. [
8]. This ratio reaches such a critical value in the region of
–
TeV and it clearly exceeds it at
TeV, where we find statistically significant violations of the
scaling law. The search for the onset of the violation of
scaling motivates our current investigation.
The
scaling function was found to increase the statistical significance of the odderon signal, based on a direct comparison of
and
at
TeV, due to at least two reasons. First of all, the overall normalization uncertainty, the largest source of the systematic error in differential cross-section measurements, cancels from
(as well as from
). Secondly, this
scaling law reduces the collision-energy-related uncertainties in the data which consequently increases the precision of the extracted quantities. We may also mention that, in the diffractive cone, expected to be valid up to at least
, the
scaling function is expected to start as
, with a well-defined normalization at the optical point:
as follows from the definitions of this scaling function. This property of the scaling functions removes the uncertainty related to possible small differences between the optical points of the differential cross-sections in elastic
and
collisions that may exist between
and
elastic collisions even if they are measured at the same centre of mass energies.
In the left panels of
Figure 1 and
Figure 2, we demonstrate the energy independence of the
scaling within statistical errors, where the TOTEM datasets at
= 2.76 and 8 TeV as well as at
and 8 TeV are compared pairwise, respectively. The agreement between these datasets corresponds to a confidence level (CL) of at least 98%.
However, one should note that in a local low
x (low
) interval one can see larger deviations than the reported data errors suggest: in the case of the comparison of the
scaling functions at
TeV and 8 TeV, we observe a partial disagreement beyond the published experimental errors in a small
x range of about
x = [3.48–7.45] on
Figure 2 corresponding to
= [0.175–0.375] GeV
2, which points toward a possible adjustment problem of the two parts of the datasets measured with different optics (
m and
m and published separately in Tables 4 and 5 of ref. [
52], respectively) and suggests the need of a further, more detailed investigation of the published systematic errors in the low
= [0.175–0.375] GeV
2 range of the dataset measured by TOTEM at
TeV with the
m LHC optics [
52].
Difficulties in describing this part of the TOTEM data at
TeV have been reported in several earlier analyses, but these problems were typically attributed to the insufficiency of the applied analysis methods; for examples, see Figures 1 and 2 of ref. [
11] or Figure 2 of ref. [
4]. So, we suggest caution when using the
TeV TOTEM data in the
= [0.175–0.375] GeV
2 four-momentum transfer region.
Fortunately, this
range, when mapped to
, overlaps only marginally with the D0 acceptance and data, at three points only, not impacting our conclusions about a model-independent observation of odderon exchange [
5] in a significant manner: the exclusion of these first three D0 data points at
TeV minimally modifies our earlier result of the odderon significance given in ref. [
5]. We have repeated the analysis of ref. [
5] by leaving out the overlapping acceptance, by removing the first 3 D0 points, but in an otherwise unchanged manner as compared to ref. [
5], and we have found that the statistical significance of odderon exchange decreased only slightly, from 6.26 to 6.10
, remaining safely well above the 5
discovery threshold. This is due to the fact that most of the signal of odderon exchange arises from the kinematic range of the diffractive minimum and maximum. A more detailed investigation about this
range dependence of odderon exchange has been presented by T. Csörgő in his invited talk at the XXXVth IWHEP conference [
55].
3. Quantification of Significance of Odderon Exchange
The
TeV and
TeV projections correspond to keeping the measured
x values at
1.96 and 8 TeV, respectively, and determining by interpolation the
scaling functions at these
x values, but at the other energy,
8 and 1.96 TeV, respectively. This procedure is described in full detail in ref. [
5]. The left panel of
Figure 3 compares the
scaling function of elastic
collisions at
TeV to that of
collisions at
TeV. In this case, adopting the method of ref. [
57], the confidence level of the agreement of the
scaling functions is found to be lower than
%, with a minimum of
. Hence, the difference between these scaling functions is statistically significant, an at least 3.79
effect, and it represents our result for the odderon observation from the comparison of the
TeV
and
TeV
datasets in the
acceptance domain. This seems to be a conservative and robust result as we find that this value can only increase for variations both in the procedure itself and in the
definition.
The quantification of the significance of the odderon is based on a method developed by the PHENIX Collaboration in ref. [
57] using a specific
definition that effectively diagonalizes the covariance matrix. In the case considered by the PHENIX Collaboration in Appendix A of ref. [
57], the experimental data are compared to a theoretical calculation. In our analysis, we adapt the PHENIX method [
57] for a comparison of one set of data directly to another set of data, without using any theory or fitting functions. Following the PHENIX method, we classify the experimental errors of a given dataset into three different types: (i) type A, point-to-point fluctuating (uncorrelated) systematic and statistical errors, (ii) type B errors that are point-to-point dependent, but 100% correlated systematic errors, and (iii) type C errors, that are point-to-point independent, but fully correlated systematic errors to evaluate the significance of correlated data [
57], when the full covariance matrix is not publicly available. Since the
t-dependent systematic errors in TOTEM measurements are almost 100% correlated, we classified them as type B errors, while the
t-independent overall normalization errors are type C errors, and the statistical errors are type A errors.
The covariance matrix has been published together with the differential cross-section of elastic scattering at
TeV by TOTEM [
18] and, on this dataset, we have checked, with a Levy series expansion method [
23], that a fit with the full covariance matrix and another fit with our adopted PHENIX method gave the same minimum, the same central values for the fit parameters and the same errors of these parameters, within one standard deviation. This suggests that the two methods are indeed equivalent in our case too, and we utilized and adapted the PHENIX method for the comparison of two datasets, where the covariance matrix of at least one of these datasets was not publicly available.
At
TeV, in ref. [
19], the TOTEM Collaboration published the
differential cross-section data with separated type A and type B errors. The source of the TOTEM
differential cross-section data, measured at
TeV, is ref. [
52]. In addition, the values of
were determined together with their errors of type A and B as given in Table 5 of ref. [
52] and Table 3 of ref. [
56]. The
t-independent, type C errors cancel from the
scaling functions, as they multiply both the numerator and the denominator of
.
The D0 collaboration did not publish type B errors for its differential cross-section data at
TeV [
14]. We have thus fixed the correlation coefficient of these D0 type B errors to zero. The input values of the nuclear slope parameters
B and the elastic cross-sections
are summarized in
Table 1, together with the appropriate references. The sources of the TOTEM
differential cross-section data, measured at
TeV, are refs. [
15,
20] with errors of type A and B. The second reference specifies data without horizontal errors; only the bin widths are given that significantly overestimate those errors. Therefore, we have performed calculations with more realistic but likely still well-overestimated horizontal errors of half bin widths, as well, resulting in significances much higher as compared to the significances presented in
Table 2.
We define the significance of the agreement between the dataset
and the projection
of dataset
to
in their overlapping acceptance, with the following
definition [
5]:
where
is the number of data points
in
indexed by
j, the same as in
but remaining in the overlapping acceptance of
sets,
,
, are the type
errors found in terms of the type M vertical errors on data point
j,
, added in quadrature with the corresponding type M vertical errors that were evaluated from the corresponding errors on the horizontal axis
x with the scaled variance method,
, where
stands for the numerical derivative of the measured quantity in dataset
at point
j in the common acceptance and
is the
j-dependent type M horizontal error. The overall correlation coefficients of the type B errors
of the
datasets are denoted by
. These coefficients are usually unknown; therefore, we perform a scan of them to find the minimum of the
, whose corresponding values are indicated on the plots.
In the right panel of
Figure 3, we present visible and statistically significant deviation from unity in the ratio of the scaling functions of
and
elastic scattering. The ratio of the
scaling functions is shown for elastic
collisions at
TeV over that of
collisions at
TeV. As a cross-check, we show the results of two different projection procedures: direct
TeV and inverse
TeV denoted by blue and red circle points, respectively. No significant variation with respect to the direction of projection has been found. In both methods, we observe an odderon effect as a peak in the
region, followed by a factor of two suppression or decrease from unity in a broad range of
. The statistical significance of the observed difference between the
and
scaling functions has been found to be at least 3.79
, which corresponds to an indication but not alone a statistically significant proof of odderon exchange from the comparison of these two datasets only. It is important to emphasize that this result is obtained by utilizing the scaling properties of
scattering without any reference to modelling and without removing (or adding to) any of the published D0 or TOTEM data points, that is, all the available data (in particular, all the 17 D0 points) were utilized in this analysis. This is in contrast to the D0–TOTEM publication on observation of odderon exchange [
39], where the analysis was limited to only 8 of all the 17 available D0 points [
14].