1. Introduction
High-energy hadron diffraction has been studied experimentally and theoretically by many groups of authors. Apart from details, there is consensus concerning global features such as the rise with energy or the shrinkage of the diffraction cone, etc. At the same time, many important details are still either unknown or disputable, such as the production of resonances in missing mass and possible structure in the differential cross section of diffractive dissociation, similar to the well-known dip-bump in high-energy proton–proton scattering. In the present paper we suggest a model predicting such structures in single dissociation. We give numerical predictions to be tested in future experiments. As far as we know, this is the first prediction of this kind.
The majority of papers on high-energy particle physics start and end with a phrase stating that the right theory is the standard model, with quantum chromodynamics (QCD) being its foundation. However attractive, this is not true until the problem of confinement is solved. Examples are hadron scattering in the “soft” (low and moderate t) region, considered in this paper, or hadron spectroscopy or the nature of confining forces producing heavy glueballs lying on the Pomeron or Odderon trajectory from massless gluons.
Our goal is to study the elastic and inelastic diffractive scattering of hadrons (protons and antiprotons) including resonance production in these reactions. We rely on methods of the analytic S matrix theory, including the Regge-pole model appended and extended by analyticity, duality and unitarity, fixing the available freedom and flexibility of the Regge-pole model.
Most of the existing criticism of the Regge-pole approach uses the epithet “old fashioned” in contrast with the “advanced” standard model (“antique” vs. “modern”), implying the uniqueness of the QCD Lagrangian compared to the flexibility of the Regge approach. Let us remind readers that the available freedom and flexibility of the Regge approach can be used for its perfection. We prefer to refer to local field theory and analytic S matrix theory in comparing these two possible approaches. We use maximal analyticity, crossing, duality and unitarity in constructing the S matrix (scattering amplitude) to restrict the available freedom, e.g., in the Regge-pole approach.
Specifically, we use double poles [
1,
2] instead of simple poles with linear Regge trajectories in the
-channel at small
as input [
3,
4]. Otherwise, the trajectories are non-linear complex functions, their form being strongly constrained, as shown in [
5] and papers cited therein. Without going into details here, we only mention that
-channel unitarity requires the presence of trajectory threshold singularities, with the lowest (in the Pomeron trajectory) at
, while at large
, the trajectories flatten, resulting in the “hardening” of dynamics and a transition from an exponential to a power decrease in the amplitude (or cross section). This is important for the direct channel trajectory
A typical trajectory featuring the transition from nearly forward to large
behaviour is
The trajectories in this duality approach play the role of somewhat dynamical variables.
In the present paper, we use a dipole Pomeron as the input, extended by a dip-bump mechanism and non-linear trajectories whose form is constrained [
5] by unitarity and analyticity. Apart from the Pomeron, the Odderon trajectory, required by the data, is included. Lower-lying, non-leading trajectories, importantly below the LHC energy range, are ignored here for simplicity when analyzing elastic scattering.
The non-linearity of Regge trajectories is an essential premise of our approach. The shape of the complex, non-linear trajectories and their observable consequences will be discussed.
High-energy hadron scattering is characterized by a forward peak followed by possible dips and bumps. The forward peak shrinks with energy (its slope increases), and this is related to the radii of the scattering particles. Multiple dips and bumps may appear in scattering of nuclei, but only a single structure was seen in
scattering [
6].
Proton elastic scattering and diffractive dissociation, single (SD) and double (DD), are closely related reactions. They were studied at the ISR, SPS and FNAL and are being studied at the LHC. Until now, no structures were seen in the differential cross sections of SD or DD.
For SD and DD, we use the model developed in a number of papers [
7,
8]. The scattering amplitude in those papers is similar to that of elastic scattering, just as the elastic vertices are replaced by DIS structure functions.
The position of the structures depends on missing masses as well as on the slopes of the SD and DD cones. The general trend is that the decreasing slope pushes the structures towards larger . Most of the present measurements at the LHC are in the region of large missing masses in a wide range of t. With the present paper, we encourage experimentalists to measure the dependence of SD and DD beyond several for different diffrative excitation values.
2. Two Components
Any hadron–hadron total cross section (or scattering amplitude) is a sum of two contributions: diffractive and non-diffractive (see
Table 1). According to the concept of two-component duality [
3], the diffractive component, the smooth background at low energies (here: missing masses) is dual to a Pomeron
-channel exchange at high energies, while the non-diffractive component contains direct channel resonances, dual to high-energy
-channel (sub-leading) Reggeon exchanges.
According to our present knowledge about two-body hadron reactions, two distinct classes of reaction mechanisms exist.
The first one includes the formation of resonances in the
-channel and the exchange of particles, resonances or Regge trajectories in the
-channel. The low-energy resonance behaviour and the high-energy Regge asymptotics are related by duality, which, at the Born level, or, alternatively, for tree diagrams can be mathematically formalized in the Veneziano model, which is a combination of Euler Beta-functions [
9].
The second class does not exhibit resonances at low energies and its high-energy behaviour is governed by the exchange of a vacuum Regge trajectory, the Pomeron, with an intercept equal to or slightly greater than one. Harari and Rosner hypothesized that the low-energy non-resonating background is dual to the high-energy Pomeron exchange, or diffraction. In other words, the low energy background should extrapolate to high-energy diffraction in the same way that the sum of narrow resonances sum up to produce Regge behaviour. However, contrary to the case of narrow resonances, the Veneziano amplitude, by construction, cannot be applied to (infinitely) broad resonances. This becomes possible in a generalization of narrow resonance dual models called Dual Amplitudes with Mandelstam Analyticity (DAMA), allowing for (infinitely) broad resonances; see [
10] and the references therein.
In the resonance region, roughly GeV, the non-diffractive component of the amplitude is adequately described by a “reggeized Breit-Wigner” term, following from the low-energy decomposition of a dual amplitude, with a direct-channel baryon (nucleon, in our case) trajectory with relevant nucleon resonances lying on it.
By duality, a proper sum of direct channel resonances produces smooth Regge behaviour, and, to avoid “double counting”, one should not add the two. Actually, this is true only for an infinite number of resonance poles. For technical reasons, we include only a finite number of resonance poles; moreover, apart from the “regular” contribution of the nucleon resonances, we only include poles lying on the
trajectory (see [
8]).
The second, diffractive component is essentially the contribution from a Pomeron pole exchange
, with
[
4], and as shown by Donnachie and Landshoff [
4], this term also can give some contribution to the low-energy (here, missing mass) flat background.
3. Analyticity, Unitarity, Crossing and Duality
The present model is inspired by the main ingredients of the analytic S matrix theory, accumulated in dual amplitudes with Mandelstam analyticity (DAMA). DAMA is a flexible representation valid in a wide kinematic range, with “soft” and “hard” modes of elastic and inelastic hadron scattering. It is flexible enough to enable various applications, in particular elastic and inelastic diffraction, the subject of the present paper. To make this point clear, we start by introducing the main ingredients of the approach, namely Regge-poles and non-linear complex Regge trajectories, playing in dual models the role of somewhat dynamical variables. Note that a single-variable Regge trajectory, connects the scattering region, with particle spectroscopy, , thus anticipating duality.
In this paper, we use the standard Mandelstam relativistic invariant kinematical variables: is the squared incoming energy, and is the squared momentum transfer, where and are the four-momenta of the incoming hadrons and is the four-momentum of one of the outgoing hadrons.
3.1. Regge Poles and Trajectories
The formal quantum-mechanical derivation of Regge-poles (Tullio Regge) was preceded by the important empirical observation by G. Chew and S. Frauchi that the hadron scattering amplitude does not factorize in
s and
t. This is supported by the behaviour of the slope of the forward peak, defined as
where
and
is the scattering amplitude.
In most of the papers on the subject, the trajectories are assumed to be real, linear functions. The origins of this prejudice are as follows: (1) the spectrum of most of the observed resonances is nearly compatible with linear real parts of the trajectories; (2) the Veneziano model works only for linear trajectories; (3) simplicity.
To complete this brief introduction to Regge trajectories, we remind that they differ by quantum numbers and by their parameters.
There are boson and fermion trajectories associated with known resonances. Typical (simplified, approximate) examples are meson
or baryonic
The above are simple examples. We use more advanced non-linear trajectories and much higher precision.
The above trajectories are non-leading, also called “secondary”, meaning that the relevant cross sections decrease with energy due to their intercepts being lower than one, Apart from the above non-leading trajectories, there are two more, the Pomeron and the Odderon, with vacuum quantum numbers introduced in the theory to explain the rise of the total cross sections (the Pomeron) and create different cross sections involving particles (e.g., protons) and antiparticles (antiprotons). The Pomeron and the Odderon differ through their -parity, positive for the Pomeron and negative for the Odderon.
The scattering amplitude with a single Regge-pole exchange is
where
is the residue and
is the signature factor,
. For small
, the signature factor simplifies to
In most of the papers on Regge-poles, one assumes that the input is a simple pole on a linear trajectory, subject to a subsequent iteration procedure based on unitarity, also called a Regge-eikonal. We instead start from a dipole and non-linear, complex trajectories. The choice of the Regge dipole is motivated by unitarity. A Regge dipole in the impact parameter representation of the form scales with the free parameters being fitted to the data. The dipole (derivative of a simple Regge-pole) produces logarithmically rising cross sections at the unit Pomeron intercept.
3.2. Duality
3.2.1. The Veneziano Model
The notion “duality” has many meanings and applications. We use it as it was introduced on the basis of finite energy sum rules (FESR) and realized in the Veneziano model in the late 1960s. It implies the equality of the properly summed contribution from direct-channel resonance to the smooth Regge asymptotic. This observation was due to numerical calculations by R. Dolen, D. Horn and C. Schmid and subsequent analytic ones by Ademollo, Veneziano and Virasoro, culminating in a simple explicit model of the amplitude in the form of the Euler Beta function by G. Veneziano [
9].
The model raised much interest, enthusiasm and hopes in the high-energy community, generating a flood of papers on and around the subject. Promising was (and remains) the generalization of the elastic amplitude to multi-particle reactions. At the same time, the limitations of the Veneziano Beta function were also evident: it is valid only for real and linear trajectories. Thus, infinitely narrow resonances, excluding scattering, result in the production of real, finite-width resonances. While immediate practical applications seemed to be impossible, the interest in duality and dual models continued in more formal directions, namely in the string theory, supersymmetry, etc., away from practical application.
Searches for applications in physics continued resulting at least in two modifications. One is the logarithmic dual model, and the other one is a class of dual amplitudes with Mandelstam Analyticity (DAMA) [
10]. Below, we briefly discuss the main ideas behind duality and present the main features of DAMA.
An elastic scattering amplitude with an infinite number of poles in the
s channel can be written in the form of the following integral representation:
where
is regular at
. By expanding
in a series near
and integrating term-by-term, one obtains
To cause the equation to become cross symmetric, we set
where
is regular at
and
and is symmetric with respect to
In this way, instead of Equation (
6), we obtain
In the simplest case of
, expression Equation (
8) reduces to the familiar Euler
B function, suggested by G. Veneziano as the scattering amplitude:
Its properties are as follows:
1. Crossing symmetry: ;
2. Meromorphic: i.e., the only singularities are simple poles in s (and in t) at (same in t).
3. Pole decomposition:
can be expanded in a pole series both in the
s or
t channel:
The Veneziano amplitude is valid only in the narrow resonance approximation, with linear trajectories.
3.2.2. DAMA
DA with MA (DAMA) is defined in [
10] as
In DAMA, trajectories are replaced by homotopies
, where
are physical trajectories and
The Veneziano limit is:
For more details, see [
10].
3.3. Super Broad vs. Narrow Resonance Approximation
Iteration/unitarization procedures are based on two premises: the choices of the input and of the iteration procedures. The procedures that were borrowed form quantum electrodynamics failed in the strong interaction theory. Below, we sketch an alternative [
10].
For simplicity, let us start with a square-root trajectory both in the
s and
t channels:
with a corresponding equation for
.
While the iteration/unitarization procedure for the Veneziano amplitude aimed to shift the infinite number of poles from the real axis to the non-physical sheet (see
Figure 1), in DAMA, it will move some poles from the real negative axis to the fourth quadrant of the non-physical
s plane, as shown in
Figure 2, producing a finite number of resonances.
Note that while in the Veneziano model, it it impossible to shift the poles from the real axis, e.g., by adding a non-linear part to the trajectories without damaging it, this is possible in DAMA, which not only allows but even requires complex trajectories. This procedure is similar to the transition from a Born series in coupling constant g to the “anti-Born” expansion in the inverse coupling .
Duality in our approach is essential in the inelastic structure functions, where the dependence on the missing mass via the direct-channel trajectory connects the low-mass, resonance region (resonances) with its Regge asymptotic.
4. Regge Dipole Model of Elastic Scattering
The basic object of this approach is a double Regge–Pomeron pole, alternatively called a dipole Pomeron (DP) [
1,
2,
11]. This object is the unique alternative to a simple pole. It has a number of advantages over the simple pole; in particular, it produces logarithmically rising cross sections at an intercept of one. It is also unique since higher-order poles are forbidden by unitarity, since they would violate unitarity. Furthermore, the dipole Pomeron scales and is invariant against
-channel unitarization, leaving its form invariant, merely rescaling the parameters that are already adjusted to the data. These are the main properties of the DP used as input.
In early 1970s, a dip-bump structure was observed in proton–proton scattering at the ISR [
6]. A single dip followed by a maximum (bump) was observed at other proton accelerators. The phenomenon was a puzzle since it was not clear whether it results from the composite structure of the nucleon (rescattering of quarks) or from geometrical considerations (shape of the nucleon) or just from unitarity, as calculated, e.g., by means of eikonal corrections, known also as a Regge eikonal [
12]. While most of these approaches predicted many dips and bumps, experimentally, a single minimum (dip) followed by maximum (bump) was observed. At the new energy frontier, the LHC confirmed this trend, leaving the puzzle unresolved.
A single Regge–Pomeron pole model does not produce any structure; rather, it produces a smooth diffraction cone, the exponential shape being slightly modified by non-linear trajectories or
-channel unitarity corrections. In Ref. [
11], the DP model was appended by a dip-bump mechanism based on a non-trivial interference between the simple and double poles moving with energy. Remarkably, the extended DP model does not generate extra dips and bumps, and it fits the data.
In view of the above considerations, we find it important to investigate the existence of similar dip-bump structures in proton diffraction dissociation. We start with the better known single diffraction dissociation (SD). This reaction was studied in detail in Ref. [
7] and reviewed in [
3,
13]. We start by reviewing the dip-bump mechanism in the DP model applied to elastic scattering to be extended to SD in the subsequent section.
The dipole Pomeron (DP) amplitude extended by the dip-bump mechanism, mentioned above, is defined as [
11]
where
, and
s and
t are the standard Mandelstam variables. Since
in the second squared brackets determines the shape of the cone, we fix
where
is the intercept of
.
is recovered via integration:
By introducing the parameter
, the Pomeron amplitude Equation (
12) can be rewritten in a “geometrical“ form:
where
and
. The Odderon contribution is
Typical fits to the data by the Regge dipole Pomeron and Odderon model are shown in
Figure 3 and
Figure 4. The values of the fitted parameters are given in
Table 2.
Figure 5 shows the contribution of particular components in the differential cross section at
TeV.
5. Dip and Bump in Elastic Scattering
A dipole amplitude generates a dip-bump structure in the differential cross section. The positions of the minimum (dip) and the maximum (bump) of the elastic differential cross section are
and
Figure 6 shows the position of the dip and the bump as a function of the slope
b. This parameter is correlated with the slope of the differential cross section. The smaller the value of the
b parameter, the higher the
value where the dip-bump structure appears. The slope of the
distribution decreases by about a factor of two when we go from elastic scattering to single diffractive dissociation. Thus, in the picture provided by the dipole Regge framework, it is natural to expect that in single diffractive dissociation, a possible dip-bump structure appears at higher
values than it does in elastic scattering.
It is known that in elastic scattering, the slope of the diffraction cone is energy dependent and it rises with increasing energy. As the energy rises, the position of the dip-bump structure moves to smaller values. The energy-dependent slope is given by the derivative of the logarithm of the differential cross section at . In the case of single diffraction dissociation, the slope of the differential cross section depends not only on the energy but also on the mass squared of the produced hadronic system, as detailed below.
Note that the above dip-bump mechanism does not depend on the form of the exchanged trajectory. For small , a linear trajectory can be used, while in the missing mass channel, the non-linearity of the complex trajectory is essential in providing finite widths of resonances.
6. Single Diffractive Dissociation (SD)
Proton’s diffractive dissociations, both simple (SD) and double (DD), were studied in many papers [
13,
14,
15,
16,
17,
18]. In a recent paper [
19], SD was considered within the Good–Walker formalism. Calculations of a particular case of central
production in pp collisions with single proton diffractive dissociation at the LHC can be found in a recent paper [
20]. To our best knowledge, possible structures in the
t distribution of SD, DD or CED have not been studied until now.
To generate a dip-bump structure in single diffraction dissociation, we introduce a dipole Pomeron (DP) and Odderon exchange to the differential cross section.
In the triple Regge approach (see [
13,
14,
15,
16,
17,
18]), the triple PPP exchange contribution to the SD differential cross section is
is nearly
t-independent [
3]. Thus, for the
t-dependent part of the amplitude of SD, we have
where the
t-dependence resulting from
is accounted for by
. Hence the
t-dependent part of the dipole Pomeron amplitude is
where
The resulting double differential cross section for the SD process is
where
By using the proton relative momentum loss variable
, we obtain
where
Using Equations (
13) and (
14), one finds the position of the dip and the bump in the SD differential cross section:
We also include the contribution of the Odderon in form of the Odderon–Odderon–Pomeron (OOP) triple exchange. In standard triple Regge formalism, it takes the following form:
We use its dipole form (the derivation is the same as in case of the dipole PPP contribution):
where
is given by Equation (
26), and the possible difference between
and
is accounted for by
. The dip and bump positions in the dipole OOP contribution are given by Equations (
29) and (
30).
The complete model for double differential cross section of SD is a sum of three triple Regge contributions appended by a pion exchange contribution:
The values of the parameters are summarized in
Table 3. The parameters
a and
b were fitted to the data. The parameter
is fixed at the value obtained in the analysis of the elastic scattering data. The intercepts of the Pomeron and the Odderon trajectories are 1, i.e.,
. The slopes of the Pomeron and the Odderon trajectories are fixed at the values obtained in the analysis of the elastic scattering data.
Interestingly, the DP Pomeron model of proton–proton SD produces properly rising total integrated cross sections with a unit Pomeron intercept, i.e., with .
7. Dip-Bump in SD
In this section, we show our predictions concerning the dip-bump structure in single diffractive dissociation at SPS and LHC.
The predicted structure in the
t distribution of the differential cross section at
546 GeV integrated in
is shown in
Figure 7.
Figure 8 shows the dip-bump structure at fixed
value. One can see that with increasing energy and fixed
value or with decreasing
at a fixed energy, the position of the dip-bump structure goes to smaller
values. The latter property is related to the
-dependence of the slope of the
t distribution of the differential cross section of the single diffractive dissociation. As
increases, the slope of the
t distribution decreases and the dip-bump structure moves to higher
values (see below).
Figure 9 shows the prediction at
and
630 GeV. Different contributions to the differential cross section are also shown. One can see that both the dipole PPP triple exchange and the dipole OOP triple exchange generate a dip-bump structure. The experimentally observable effect is expected from the OOP contribution.
The predicted dip-bump at
TeV for selected missing masses
is shown in
Figure 10.
Our descriptions for double differential cross section of proton–antiproton SD at fixed
and
GeV as a function of
are shown in
Figure 11.
Predictions for the integrated differential cross section of proton–proton SD at
TeV as a function of
are shown in
Figure 12.
Our prediction for the dip-bump structure of the
integrated differential cross section of SD at
8 TeV is shown in
Figure 13.
Figure 14 shows the motion of the dip and bump as functions of
at
8 TeV. As
increases, the slope of the
t distribution decreases and, consequently, the dip and bump move to higher
values.
The necessary data sample for measuring the dip-bump structure in SD at LHC energies in future experiments can be estimated based on the following assumptions. The differential cross section at the dip position is on the order of . Taking bins in missing mass of 5000 and in t of = 0.1 , the cross section per such a phase space cell is about 0.05 pb.
A measurement with statistical uncertainty better than 10 requires a data sample of a minimum of 100 events in the phase space cell at the dip position. Assuming a factor of 0.2 representing geometrical acceptance, pile-up rejection and reconstruction efficiency results in a necessary integrated luminosity of . If the SD events can be measured in both beam directions, then this integrated luminosity requirement is reduced by a factor of two.
8. Summary
In this paper, predictions concerning structures in the differential cross section of single diffraction in the LHC energy region are presented. Our model was tested on elastic proton–proton data, including the dip-bump structure in the differential cross section, then it was generalized to proton single diffractive dissociation. We have elaborated the model to account for the expected structures in the differential cross section of diffractive dissociation. This is the first prediction of these structures—the existence of a dip and bump. The predictions concern the position of these structures, their dependence on the incoming energy and missing mass. Our predictions (numbers) may slightly change (by a maximum of ) depending on the results of future measurements, for example, those on a slope, without affecting our main message: a single dip and bump are expected, and their position and motion are predicted.
While structures in elastic scattering still remain an important part in studies of high-energy hadron diffraction, a subject of controversy, related studies in SD and DD are missing in the existing literature, probably because of the complexity of the problem. Its understanding is important for further progress in the field.
Our quantitative predictions concerning the structure are based on the relevant values of the SD cross sections and slopes, measured only indirectly or fragmentarily, the missing values/intervals being filled in either by extrapolations or model predictions that make these studies exciting, especially in view of future measurements that will verify the predictions.
To summarize, we conclude the following:
(1) The Regge-pole model constrained and extended by analyticity, unitarity and duality is the basic tool to handle hadron scattering in the ”soft” (low and moderate t) region.
(2) In constructing realistic models, one should not necessarily follow the “standard” procedure based on a simple Regge-pole with a linear trajectory as input. We have presented an alternative approach based on a Regge dipole which can be extended by a non-linear, complex trajectory as input.
(3) The present approach enables a smooth, continuous transition between “soft” and “hard” dynamics, i.e., from the “soft” exponential to the “hard”, power-like decrease in the cross sections, obeying quark-counting rules and/or perturbative QCD. In terms of the Regge-pole theory, this transition is mimicked by logarithmic Regge trajectories.
To conclude, a dip followed by a bump is predicted in the t distribution of single diffractive dissociation at in the range of 3 7 . Apart from the dependence on s (total energy squared) and t (squared momentum transfer), the dependence on missing masses, both at large and low masses in the resonance region, is predicted. We are looking forward to future SD measurements at the upgraded LHC and other future experimental facilities.