Local one-dimensional reggeon model of the interaction of pomerons and odderons

We propose the one-dimensional reggeon theory describing local pomerons and odderons. It generalizes the well-known one-dimensional theory of pomerons (the Gribov model) and includes only triple interaction vertices. The proposed theory is studied by numerical methods: the one-particle pomeron and odderon propagators and the pA amplitude are found as functions of rapidity by integrating the evolution equation.


Introduction
In the framework of the Quantum Chromodynamics in the kinematic region where energy is much greater than the transferred momenta ("the Regge kinematics") the strong interactions can be described by the exchange of pomerons, which can be interpreted as bound states of pairs of the so-called reggeized gluons. In the quasiclassical approximation (which neglects pomeron loops) and in the approximation of a large number of colors, for the scattering of a small projectile off a large target ("dilute-dense scattering") it leads to the well-known Balitsky-Kovchegov (BK) equation widely used for the description of DIS and particle-nucleus (pA) scattering. The BK equation corresponds to summing fan diagrams going from the projectile to the nuclear target with the propagator given by the well-known BFKL equation and the triple pomeron vertices responsible for the splitting of a pomeron in two. Going beyond the quasiclassical approximation and taking account of loops presents a hardly surmountable problem, which has been not solved until now.
In connection with this difficulty much attention was given to the previous attempts to study the strong interactions before the era of the QCD and using the old reggeon theory [1,2,3] introduced by V.N.Gribov and based on the phenomenological local pomeron and its interaction vertices. In fact the pA interaction in this framework was considered in [4] where the sum of all fan diagrams similar to the BK equation in the QCD was found. Unlike the QCD, in the local pomeron model both the pomeron intercept and coupling constant for the triple pomeron vertex are taken as phenomenological parameters adjusted to the experimental data. Needless to say the local pomeron is much poorer in his physical content as compared to its QCD counterpart, which makes it unfit to describe processes with hard momentum transfer like DIS. However, the local pomeron theory is much simpler than the QCD one and admits various methods which make it possible to go beyond the quasiclassical approximation. In particular renormalization group methods have been widely applied to establish some basic properties of the Gribov model. In [3,5,6,7] a reggeon field theory with zero renormalized reggeon mass has been investigated by means of the Callan-Symanzik equation and the ǫ-expansion was studied in the vicinity (D = 4 − ǫ) of the critical dimension of transverse space. The key result was the existence of an infrared fixed point which leads to scaling laws for the pomeron Green function. More recently a general reggeon field theory as investigated in the framework of the functional renormalization group formalism with the aim of finding fixed points for the coupling constants flows, which may shed light on comparing it with the QCD and relevance for the description of the region of low momenta [8].
Still in the realistic three-dimensional world even the local pomeron model does not allow to find the full quantum-mechanical solution of the problem with the contribution of pomeron loops fully taken into account. In view of this trouble a still simpler model ("toy" model) was considered and studied in some detail. The local pomeron model was taken in the onedimensional world, that is depending only on rapidity y [9,10,11,12,13]. Such a model essentially was equivalent to the standard Quantum Mechanics and studied by methods of the quantum mechanics. Later also numerical methods were used both for the evaluation of the eigenvalues of the quantum Hamiltonian [14,15] and for integrating the partial differential equation which governs the evolution of this system in rapidity [16,17]. The important message which follows from these studies is that the quantum effects, that is the loops, change cardinally the high-energy behaviour of the amplitudes and so the neglect of them is at most a very crude approximation.
It is remarkable that in the QCD, apart from the pomeron with the positive C-parity and signature, a compound state of three reggeized gluons with the negative C-parity and signature, the odderon, appears. Its possible experimental manifestations has not been found with certainty up to now, which may be explained both by its behaviour with energy and its small coupling to participant hadrons. On the theoretical level two species of the odderon were found, the Bartels-Lipatov-Vacca (BLV) odderon [18] with the intercept exactly equal to unity, in which the three reggeized gluons are pairwise located at the same spatial point, and the more complicated Janik-Wosiek odderon [19,20] with all three reggeized gluons at different points, the intercept somewhat below unity and so probably subdominant at high energies. It was noted in [22] that the BLV odderon is in a certain sense an imaginary part of the full S-matrix with both C = ±1 exchanges whose real part is the pomeron. Having this in mind, in analogy with the BK equation a system of equations summing fan diagrams made of both pomerons P and odderons O has been derived [21,22]. It takes into account transitions P ↔ P P , P ↔ OO and O ↔ P O with the coupling constants related as 1:-1:2, respectively. This system is more complicated than the BK equation but can be solved numerically as the latter [22,23,24]. However, going beyond the quasiclassical approximation remains unrealizable. In the framework of the local reggeon model the odderon has been introduced in [25] again within the functional renormalization group formalism and the same goal of finding possible fixed points for the coupling constants flows. However, even forgetting strong limitations on the number of terms in the effective action, the found fixed points tell nothing about the actual values for amplitudes with given phenomenological parameters.
In view of this and learning the lesson from the studies of the theory with only pomerons, in this paper we propose a generalization of the one-dimensional "toy" model to include interactions of both the pomerons and odderons. Note that earlier results found in the functional renormalization group formalism in [25] do not help much to construct such a model due to the essential use of infrared regulators inherent in this formalism. The Lagrangian function of the one-dimensional model with only the pomerons is where Φ, Φ + are the complex pomeron field and its conjugate, the mass parameter µ P = α(0)−1 is defined by the intercept of the pomeron Regge trajectory and λ is the effective coupling constant.
To generalize (1) we are guided by the properties of the transition vertices P ↔ P P , P ↔ OO and O ↔ P O found both in the analysis of the loop diagrams in the QCD [26] and in the quasiclassical equations [21,22]. As derived in [26] the one loop corrections should be negative when a pomeron is present in the loop and positive for the odderon loop in the pomeron propagator. From [22] we deduce that the intercept of the odderon is unity and relation of the coupling constants P → P P , P → OO and O → P O should be as in the quasiclassical equations. We have also to take into account the negative signature of the odderon, which makes its contribution to the amplitude real in contrast to the positive imaginary one for the pomeron. As will be discussed in the next Section, these properties dictate the form of the Lagrangian in our new local reggeon model as Here Φ is the complex pomeron field, Ψ is the complex odderon field, Φ + and Ψ + are their conjugates, µ P and µ O are the pomeron and odderon mass parameters equal to the intercepts minus 1, respectively, the effective coupling constant λ is the same as in the pure pomeron theory (1). To take the odderon signature into account its free Lagrangian is multiplied by i, which multiplies its propagator by −i. Apart from this we have taken real the coupling constants for Ψ + ΨΦ and its conjugate in contrast with all the rest coupling constants kept positive imaginary. This is necessary to conform to the QCD loop contributions and quasiclassical equations [21,22,26]. This theory is formally different with the one introduced earlier in [25], however, we shall see that in fact they are fully equivalent after the appropriate transformation of field variables. This theory describes both creation and absorption of pomerons and odderons and contains interactions only with even powers of the odderon field, hereby conserving the C-symmetry. It also has the relation between the coupling constants of different transitions in correspondence with the QCD. In our work we will start with the arbitrary µ O , but in the numerical studies we adopt the zero mass parameter for the odderon in accordance with property of leading BLV odderon.
In the next Section from the Lagrangian function (2) we construct the quantum Hamiltonian of the model, study the equations of motion and compare with the quasiclassical (fan) equations in the QCD [21,22]. In Section 3 we consider evolution in rapidity of the pomeron and odderon propagators and pA amplitude in our model and suggest two calculational schemes for its numerical study: development in powers and point-like evolution. In Section 4 the numerical results found by these methods are presented. In Section 5 the effect of introducing of the quartic interaction in the theory is studied. Section 6 contains some conclusions.

Quantization
The free part of the Lagrangian (2) has the standard form and admit the canonical quantization. Treating y as the imaginary time one finds that the canonical conjugates to the variables Φ and Ψ are respectively. The canonical quantization in the imaginary time leads to the following commutation relations The Lagrangian can be presented in the form Here V (Φ, Ψ) is the interaction term which does not contain derivatives in y. Then the Euclidean Hamiltonian is where the interaction V is Comparing with the pomeron-odderon interaction introduced in [25] we observe that in our case all triple interactions are imaginary but the relative sign of P → OO and OO → P transitions is opposite. This guarantees that the OO loops are have a different sign as compared to P P loops. If one does the canonical phase transformation Ψ * → e −iπ/4 Ψ * , Ψ → e iπ/4 Ψ, then the two last terms in (7) will appear with the same real coupling constant λ as in [25]. So our pomeron-odderon interaction is fully equivalent to one in the latter, provided the three coupling constants there have the same magnitude.
The relations (4) allow one to interpret Φ * , Ψ * as the creation operators and Φ, Ψ as the annihilation operators. The vacuum satisfies the condition Φ|0 = Ψ|0 = 0 and the Fock space is considered as a full space of states created from the vacuum by action of any number of Φ * and Ψ * . The Schrödinger equation for the wave function of the system F (y, Φ * , Ψ * ) is Note that the original field Ψ + = −iΨ * has a meaning of the creation operator of the odderon state.

Equations of motion and comparison with the QCD equations for amplitudes
The corresponding to (6) differential equations of motion for the field operators are Note that the connection (3) with the Hermitian conjugates with respect to the standard scalar product is not conserved during the evolution in y since the Hamiltonian is not Hermitian and the dynamics of the system is not unitary.
We are interested in these equations to compare them in the fan approximation with the QCD. The fan equations are obtained when all terms which contain two creation operators and one annihilation operator are dropped from the Hamiltonian. In our case one has to consider each of the equations (9) without the last term in the parentheses. Here we can introduce Φ = −iχ, Ψ = −iω to get from two first equations One can see that these fan equations coincide in their form with the fan equations in the QCD. Indeed, the coupled equations for evolution of the odderon together with the pomeron, derived in [22] in the transverse coordinate space are where N(x, y; y) and O(x, y; y) are the C-even and C-odd dipole amplitudes, respectively, symmetric and antisymmetric in x, y. Hereᾱ s = N c α s /π. Following [23] consider the translational invariant situation when the amplitudes depend only on the difference r = x − y. Then passing to the momentum space we get the equations [23] ∂χ(k, y) ∂y where H BF KL is the standard LL BFKL Hamiltonian. For the local pomeron and odderon both χ and ω do not depend on k and depend only on rapidity y. Then the equations simplify to ∂χ(y) ∂y = µ P χ(y) −ᾱ s χ 2 (y) +ᾱ s ω 2 (y), where the mass parameters µ P and µ O have a sense of the pomeron and odderon intercepts minus 1, respectively. If one chooses λ =ᾱ s ≡ N c α s /π then the equations (15) completely coincide with our equations (10).
It is to be noted that the choice of real vertices for transitions O → O + P and O + P → O in our original Lagrangian (2) was dictated by the form of the equation for ∂ω/∂y. With the imaginary vertices one would obtain factor i in the right-hand part, which could not be eliminated by rescaling of ω, since the equation is linear in it.
We acknowledge that the use of the one-dimensional description without the transverse space dependence for the odderon field, which is antisymmetric in k, is problematic. The comparison of terms of the equations connected with local interactions can, however, have a sense. It shows, at least, that coefficients before different terms and their signs coincide in (9) and (15).
Finally, the comparison of the fan equations fixes the coupling constant for Φ * ΨΨ interaction but not for the last term in (7). Here we choose the constant for the last term in (2) negative imaginary (or, equivalently, the opposite signs before the two last terms in (7)) to take into account the signature properties of odderon. With the given choice of constants the simple loop constructed from two odderon lines gives the positive correction to the pomeron propagator, whereas the simple loop from one odderon line and one pomeron line gives the negative correction both to the pomeron and odderon propagators.

Passage to the real Hamiltonian
The Lagrangian function (2) is to be accompanied by the form of the coupling of the pomerons and odderons to the external particles. We assume the eikonal form for the operators of creation of the initial state and annihilation of the final state (analogous to ones taken in [16,17]): with possibly different coupling constants of interaction with the projectile and the target. We choose g p with a negative sign to later deal with well behaved expressions. The Hamiltonian (6) with interaction (7) is complex, which makes practical use of it inconvenient. So one can pass to a real Hamiltonian. One possibility (case A) to do it is to introduce the Fock-Bargmann representation of operators In terms of these operators The initial and final states become So the coupling of the odderon to the projectile and the target becomes pure imaginary. Accordingly, our Hamiltonian is now real but the initial and final states are not. To separate real and imaginary amplitudes one has to separate terms with even and odd numbers of odderons, that is having C = +1 and C = −1. This can be done separating terms proportional to even or odd powers m = m i + m f of the product g The alternative possibility (case B) is to only pass to imaginary Φ and Φ * and do not change Ψ and Ψ * : and to retain the form of the odderon coupling to the initial and final states. The alternative Hamiltonian is As compared to the previous case the signs of the two last terms will be opposite. However, it is trivial to find out that this possibility is equivalent to the canonical transformation Ψ → iΨ, Ψ * → −iΨ * and in the end gives the same result for the amplitudes.
3 Evolution in rapidity

Basic equations and methods for solution
If the initial state in the Schrödinger picture at rapidity zero is given by then at rapidity y it will be given by with the evolution to rapidity y given by the equation The amplitude for the transition from the state F y=0 |0 to a state 0|F (f in) at rapidity y will be given by the matrix element We take the Hamiltonian in the real form (18) with the initial and final states in the form (19). The amplitude becomes where is the result of the evolution of the initial state to rapidity y.
Applying the derivatives we find that is the amplitude is obtained by just substituting u and w in the evolved F (u, w) with the coupling constant for the target. To find the propagator one evidently has to take the terms linear in these coupling constants (or in u or w). So the pomeron and odderon propagators are obtained as with the odderon propagator actually carrying extra factor i. Passing to the possible methods to solve the evolution equation (24) with the Hamiltonian given by (18) or (21) we can find two basic alternatives. The first one is to study the Hamiltonian and find its eigenvalues. This is a difficult problem already without odderon, as demonstrated by earlier studies. With only the pomeron the Hamiltonian is non-Hermitian and acting on a complex field variable. Fortunately, in this case it could be transformed to a Hermitian Hamiltonian, which allowed to approximately find its ground state at very small coupling [11]. Later a formalism has been developed to find all eigenvalues from the initial Hamiltonian in the complex variable [15]. Inclusion of the odderon substantially aggravates the situation. Now the initial non-Hermitian Hamiltonian depends on two complex variables. Its analysis requires much more efforts and transformation to a Hermitian Hamiltonian becomes hardly possible.
So we are left with the second alternative, which is to directly integrate the evolution equation numerically. This procedure turned out quite efficient for the pure pomeron model at physically reachable energies [16], although of course it cannot be stretched to infinite energies when the groundstate eigenvalue becomes essential. So we turn to direct integration of (24) starting from some initial condition at y = 0.
One may consider two possibilities for the choice of variables for the wave function. One may develop F y (u, w) in power expansion in u and w This form has a nice physical interpretation. Term with u n w m corresponds to a state with n pomerons and m odderons. So expansion (30) immediately gives the pomeron-odderon content of the wave function. However, convergence of this expansion is not guaranteed and the cut series obviously has a very bad behaviour at large u and w. So, as we shall see, the applicability of this method is severely restricted to small values of y and the coupling constant λ.
Another possibility actually employed in [16] is to take the initial function F y=0 (u, w) on the lattice (u i , w j ) and evolve F on this lattice. In this case the evolution equation itself is responsible for the high y and λ behaviour. As a result this method can be applied in a wide region of y and λ.
In the next subsections we consider these numerical methods in more detail.

Evolution by points. Composite fields
The method similar to one used in the previous papers [16,17] is to use the evolution equation as it is, choosing some initial function at a set of points u 0 , . . . , u N , w 0 , . . . , w N and evolving it with rapidity y by the Runge-Kutta method. As a result of this numerical exercise one gets F y (u 0 , ...u N |w 0 , ...w N ). This method worked for the pomeron evolution without the odderon in a very wide choice of the model parameters.
Note, however, that this method requires a very particular choice of variables u and w. In the Gribov model studied in [16,17] the wave function is an analytic functions in u. Obviously, any analytical function can be completely determined by its values on some ray in the complex plane. However, one experience drawn from those papers was that the initial condition for the evolution has to be imposed on the positive real axis iΦ * = u > 0, otherwise the numerical evolution process breaks down. The reason may be in that in practice one has to limit the values of Φ * (or u). Convergence requires the wave function to fall at large values of the argument and this is possible only for a very particular direction in the complex Φ * -plane. As a result one has obligatory to pass from Φ * to u = iΦ * and start the evolution from positive u.
With the odderon included we have the same situation with two variables Φ * and Ψ * and from the start it is not clear which values for these variables in their complex plane allow convergence. We have immediately discovered that with the choice u = iΦ * and w = iΨ * evolution from the positive axes for both variables is impossible, leading to divergence at quite small values of y < 1.
So we are bound to seek for different variables in which evolution becomes realizable. As a hint one may consider a model (which will be the subject of our subsequent paper) in which the fake odderon does not possess the negative signature and moreover has the same intercept as the pomeron, their difference reduced only to the C-parity. In such a model the choice of variables for evolution is trivial, since it follows this choice for the model with only the pomeron.
Based on this hint we introduce composite fields as linear combinations of the pomeron and odderon fields. We define Their commutation relations are which shows that ϕ andφ are independent dynamical variables. If one takes the interaction Hamiltonian as a sum of two interactions of the Gribov model (1) depending on different composite variables and with an appropriate coupling constants one finds or in terms of elementary fields Comparing with the true interaction V with the odderon signature taken into account we find where So the full Hamiltonian in terms of composite fields is where the free part (in the physical case µ O = 0) is Following the experience with only the pomeron field we pass to real variables introducing We stress that these u and w are not the old u and w in (17). We denote them so only to economize in notations. Then the final Hamiltonian in variables u and w is found as This Hamiltonian is symmetric in variables u and w. So once the initial wave function F (u, w) is (anti)symmetric in u and w it will preserve this property in the evolution in rapidity.
It is remarkable that if one takes the initial function F (u, v) on the positive axis of both u and v then its evolution in rapidity become possible in a rather wide region of parameters µ P and λ and also to rather high values of the rapidity. The results of this evolution for different sets of (µ P , λ) are reported in the next Section.
Note that the pomeron and odderon fields are expressed as or in terms of u and w and the evolution equation for the wave function takes the form (24) with the Hamiltonian (47). To start evolution in terms of ϕ,φ we choose (49), dropping the common factor (−i), as the initial state F y=0 (u, w) for the pomeron and odderon, respectively. The expressions for the propagators become instead of (29). Since the exchange of u and w is the C-transformation, the C-symmetry of the Hamiltonian assures that the wave function remains symmetrical under u ↔ w for the pomeron and antisymmetrical for the odderon at any y. Furthermore, C-symmetry guarantees that the symmetrical and antisymmetrical parts of the wave function evolve independently that allows one to choose the sum u √ 2 of (49) as the initial state and to find both propagators from (50) in one numerical simulation.

The propagators
The pomeron propagator P (y) calculated by power expansion for µ = 0.1 is shown in Fig. 1 by upper curves in the left upper panel for λ = 0.01 and the left lower panel for λ = 0.03. To compare and see the precision the lower curves present our far more precise results obtained by point evolution. In the right panels we similarly present the odderon propagator for the same values of µ and λ. For the odderon the lower curves were obtained by power expansion and the upper ones by point evolution. One observes that the power expansion gives satisfactory results at smaller values of λ and rapidity. With the growth of both its precision worsens notably and at λ = 0.03 its convergence breaks down already at y ∼ 15. The influence of loops and the odderon in the whole region 0 < y < 30 will be illustrated later, in Fig. 5 based on the precise calculations by the point-like evolution.

The pA amplitude
The pA amplitude may be C-even and C-odd and depend on the couplings of the pomeron and odderon to the proton in the projectile and target. We take g The C-even and C-odd amplitudes are presented in Figs. 2 and Fig. 3. In Fig. 4 we show the even amplitude with the odderon coupling equal to zero, when the whole influence of the odderon reduces to its loop contribution. In all these pictures the pA amplitude is shown together with the predictions from fan diagrams. As one observes in the C-even amplitude the combined effect of the pomeron and odderon loops is the same as in the propagators: they diminish the amplitudes, just as without odderon. The situation with the C-odd amplitude is not clear, since it is quite small and the loop effects become visible at rapidities probably outside the region of convergence, which is observed at y > 5 for λ = 0.03.

Evolution by points
As mentioned, unlike evolution by powers, evolution of Eq. (24) on the (u, w) lattice admits wide areas in y, µ and λ using different steps in y. The latter depend on the values of λ. For λ ≤ 0.05 satisfactory results were obtained with N = 400 points in the interval (0, 20) for both u and w and 2000 points for ∆y = 1. However, for λ = 0.1 we had to raise the precision in y to take 20000 points for ∆y = 1. With λ = 1 one has to take 100000 points for ∆y = 1 with the corresponding rise in the processor time. Raising the number of points in u and w to N = 500 changes both propagators only by few percents. Below we report on the calculated pomeron and odderon propagators, for which we took symmetric and antisymmetric initial conditions (u ± w)/ √ 2, respectively. We do not report our results for pA amplitudes, since their behaviour more or less follows the one found above via the evolution by powers.

Pomeron and odderon propagators
We start with our results for the pomeron at small values of µ and λ, where also power expansion gives converging results previously shown in Fig. 1.
In Fig. 5 the solid (middle) curves show the pomeron propagator with µ = 0.1 and λ = 0.01 (left panel) and 0.03 (right panel). The bottom curves show the results obtained without the odderon (only pomeron loops) and the upper curves show the propagator without any loops. One concludes that the odderon loops somewhat enhance the propagator, which is natural as they bear the opposite sign. At λ ≥ 0.04 evolution by power expansion breaks down. So the rest plots for the pomeron and odderon propagators shown in Figs. 6-10 use evolution by points. For the pomeron we compared our results including odderon loops with our old calculations without odderon [15,16]. We consider the cases with (µ, λ) = (1, 0.1), (−1, 0.1), (1, 1/3), (1, 1), (0.1, 1), which fully enough illustrate the dependence on the intercept 1 + µ and triple coupling constant λ > 0. The dashed curves in all plots of the pomeron propagator correspond the results without odderon.
Observing our results we see that the effect of the odderon loops is universally constructive. As expected they act in the opposite direction compared to the pomeron loops making the propagator somewhat larger. This effect is enhanced with the growth of rapidity. At large values of λ the effect of odderon loops starts at quite small y and becomes very strong, see Figs. 8-10. In all cases at large enough y the pomeron propagator goes to zero.

Quadruple interaction
The conclusion that only the triple interactions are essential at large energies is based on the renormalization group analysis [2,3] and is not applicable for the one-dimensional model. To introduce the quadruple interaction in our theory we shall exploit the idea that the model can be expressed via the composite fields (38) and add a term to the Hamiltonian In terms of the pomeron and odderon fields one finds Then the quadruple interaction (51) can be written as Passing to variables u, w (17) we have One can try to use these form of the interaction to perform the evolution by points in the space (u, w). Unfortunately the additional quadruple terms drastically spoil convergence and do not allow to move above y = 1.
In contrary, with a non-zero λ ′ evolution by powers allows to calculate the pomeron and odderon propagators at rather high values of λ and λ ′ at least up to rapidity y = 10.
In the evolution in powers of u and w this interaction generates an additional term f In Fig. 11 we show the propagators at µ = 0.1, λ = 0.1 and values of λ ′ = 0, 0.5, 1 and 2. Convergence becomes better with the growth of λ ′ . In particular at λ ′ = 0 that is without the quadruple interaction, the power series method diverges and the corresponding curves could only be obtained by the point evolution.
Calculation shows that the quadruple interaction somewhat raises both the pomeron and odderon propagators which are illustrated in Fig. 11. The curves show the propagators successively enhanced as the values of λ ′ are raised: λ ′ = 0, 0.5, 1, 2. So the quadruple interaction (51) acts in the opposite direction to the loops coming from the triple interaction providing a constructive contribution.

Conclusions
We have proposed an one-dimensional reggeon model for the interaction of local pomerons and odderons at different energies. The model is a natural generalization of the well-known "toy" model for the pomerons only and follows the structure of the pomeron-odderon evolution in the QCD. It allows to study evolution of interacting pomerons and odderons with quantum effects (loops) fully taken into account and so estimate these effects and validity of the quasiclassical approximation.
Practical methods to study evolution in rapidity were proposed. In the general case they can be pursued by the same numerical procedure that turned out to be successful for the model without odderons, provided the initial wave function is adequately chosen. In particular, naive independent initial choice of the pomeron and odderon contribution is not allowed for evolution, since, probably, it contradicts the asymptotic properties of the wave function considered as a function of two complex field describing the pomeron-odderon system.
The found numerical results show that the inclusion of odderon loops enhances the pomeron propagators but somewhat damps the pA amplitudes. In particular, the odderon propagator strongly falls with energy in agreement with conclusions from the QCD. Remarkably the influence of odderon loops is clearly felt both in the propagators and amplitudes. also leads to overall damping with energies.
We also studied a certain quadruple interaction choosing it in a form symmetric in pomerons and odderons. This interaction acts constructively and raises both the pomeron and odderon propagators.