Finite nuclei in the reggeon “toy model”

Hadron–nucleus amplitudes at high energies are studied in the “toy” Regge model in zero transverse dimension for finite nuclei, when the standard series of fan diagrams is converted into a finite sum and loses physical sense at quite low energies. Taking into account all the loop contributions by numerical methods we find a physically meaningful amplitudes at all energies. They practically coincide with the amplitudes for infinite nuclei. A surprising result is that for finite nuclei and small enough triple pomeron coupling the infinite series of fan diagrams describes the amplitude quite well in spite of the fact that in reality the series should be cut and as such deprived of any physical sense at high energies.


Motivation
At high energies in the framework of the perturbative quantum chromodynamics strong interactions are mediated by the exchenge of hard pomerons, which are non-local entities propagating according to the BFKL equation and splitting into two or merging from two to one with the known triple pomeron vertex. Neglecting pomeron loops and choosing the projectile to have a short range and the target to be a heavy nucleus one comes to the well-known Balitski-Kovchegov (BK) equation, which sums pomeron fan diagrams going from the projectile towards the target [1,2]. This equation with certain degrees of sophistication, including higher orders and running coupling, is widely used in applications with very positive results. However, from the start it is clear that summing all fans this equation neglects the obvious limitation that the nucleus in fact is finite, so that the number of splittings in fans is restricted. It is not clear how this restriction affects the resulting amplitude. This fact cannot be too small. For instance, taking instead of a heavy nucleus a light one, say the deuteron, one sees that the results become drastically different, since then the amplitude essena e-mails: braun1@pobox.spbu.ru; mibraun@yandex.ru tially reduces to the exchange of two pomerons and grows correspondingly at high energies.
This problem is difficult to study in the framework of the full-fledged QCD pomeron theory. However, very long ago a description of the hA interaction by means of the sum of fan diagrams was proposed by Schwimmer in the reggeon field theory with a local supercritical pomeron [3]. There one easily obtains a solution of the equation, which sums all fans and is basically similar to the BK equation. Taking into account that the momenta transferred to the nucleus are small one gets for the hA amplitude at given rapidity y and fixed impact parameter b Here μ is the pomeron intercept minus unity, assumed to be positive, T (b) is the nuclear profile function normalised to unity, λ > 0 is the triple pomeron coupling with the opposite sign and g is the pomeron-nucleon coupling. This old formula possesses some nice features. In particular the Schwimmer amplitude goes to a constant value at very high energies, implying that at such energies the nucleus behaves like a (grey) disk, which more or less agrees with the later more sophisticated treatments. We shall be interested in not so much of its physical applications but rather as a tool to study our problem: how the limitations on the number of splitting coming from the actual finiteness of A change the result.
Having the explicit solution (1) this is quite trivial to see. Let so that the Schwimmer amplitude is just For a finite nucleus with atomic number A obviously we get just A first terms of the expansion of (3) in powers of z Here and in the following we denote with subindex A the amplitudes which refer to finite nuclei of atomic number A, leaving the amplitude for infinite nucleus without any subindex. One immediately sees that if z < 1 then the series in (4) is convergent and the error in using the infinite nuclei instead of the finite (and physical) is exponentially small. On the other hand, if z > 1 the series (4) is divergent and the Schwimmer formula for infinite nucleus has nothing in common with the real amplitude for finite nucleus. Condition of validity of the fan diagram amplitude z < 1 translates into the restriction on the highest rapidity where this amplitude has the physical meaning If the internucleon distance in the nucleus is R 0 then crudely estimating we have and (5) gives With fixed λ independent of A and large A >> 1 this degenerates into which has little physical meaning. The only possibility to have some sense for the Schwimmer amplitude for realistic nuclei is to choose λ extremely small to compensate factor A −1/3 in (6). Otherwise this formula is just an analytic continuation of the physical amplitude having no relation to reality. This problem is of course not a new one and is not restricted to only fan diagrams. In a simpler case of the scattering of a hadron on a nuclear target in the Glauber approximation one gets the amplitude where a is the forward proton-proton scattering amplitude.
In the limit A → ∞ and T (b) ∝ A −2/3 one gets the standard expression which is quite attractive, since it is explicitly unitary. However, Eq. (8) is unitary only while |a(y)T (b)| < 1. This is always so when the proton-proton scattering amplitude is unitary itself. However, if one takes for a(y) the amplitude corresponding to the exchange of a supercritical pomeron and so rising with y as exp(μy) then the expression for finite nuclei loses sense for high enough y. The limiting expression for infinite nuclei preserves its unitary character but its relation to the physical amplitude becomes lost.
Observing these examples we may conclude that in both cases, fans and Glauber, the origin of the difficulty lies in the wrong behaviour of the elementary proton-proton scattering amplitude (pomeron propagator) growing exponentially with rapidity and violating unitarity. This wrong behaviour is possibly cured by inclusion of contributions from pomeron loops. We cannot be sure that this happens in the reggeon field theory, which lies at the basis of the Schwimmer equation (1). Still less is known about the behaviour of the propagator of the non-local pomeron in the perturbative QCD. However, there is a simple model where calculation of all loop contributions is possible. This is a reggeon "toy model" in the zero-dimensional transverse space. Having rapidity as the only variable it actually reduces to a sort of quantum mechanics with a non-Hermithean interaction. This theory has been extensively studied in the past [4][5][6][7][8][9] and in the limit λ → 0 it was shown analytically that inclusion of loops makes the pomeron propagator vanish in the high-energy limit. More recently it was considered in [10,11] and in [11] a calculational technique was elaborated which allowed one to numerically sum all contributions for arbitrary values of parameters. In this paper we apply this technique to study the behaviour of the hA amplitude in the model for finite nuclei and establish the relation between the sum of fan diagrams and the full amplitude in this case.
Note that the fan amplitude in the toy model is identical to the more physical Schwimmer amplitude (1). The simplification of neglecting the transverse space is felt only in the loop diagrams, which in the toy model are certainly different from the ones in the reggeon field theory with transverse dimensions. Still we expect that lessons known in the study of the toy model will be instructive to consider the situation in more physical theories including the perturbative QCD.

The toy model
In this section we briefly recapitulate the definition and properties of the toy model indispensable for our study, referring the reader either to old papers [4][5][6][7][8][9] or to the comparatively recent paper [11]. The toy model is the Regge-Gribov theory of a pomeron field φ(y) depending only on rapidity y (zero-dimensional transverse space), which may be defined by the functional integral where Here μ is the pomeron intercept (α(0) − 1). For the supercritical pomeron μ > 0. Triple pomeron coupling constant λ is also positive. The functional integral (10) converges for μ < 0 (subcritical pomeron). But in the physically interesting case that μ > 0 the integral does not exist. Then in fact it only serves to introduce perturbative diagrams in the Regge-Gribov approach.
One can pass to an alternative, Hamiltonian formalism, which reproduces the perturbative diagrams but is free from the restriction μ < 0. It is based on a quasi-Schroedinger equation in rapidity for the wave function (y), with the Hamiltonian H which can be chosen to be real and is a function of two operators u and v, which are anti-Hermitean to each other, and satisfy the commutation relation The operators u and v have the meaning of creation and annihilation operators of the pomeron, respectively. The vacuum state 0 , normalised to unity, satisfies v 0 = 0. All other states are built from 0 by application of some number of operators u. The transition amplitude from the initial state i at rapidity y = 0 to the final state f at rapidity y is given by The amplitude A f i is imaginary positive so that the matrix element on the right-hand side of (16) is negative. Some care should be taken to express the initial and final scattering states i (u) and f (u) via creation operators. We take them also to be real. Assuming that the initial state representing a heavy nucleus with A → ∞ has an eikonal structure we take where g i is a positive coupling constant with the initial nucleus. It is important that the final state should be taken not as an immediate copy of (17) (with a maybe different coupling constant) but with an additional change u → −u As we shall see this immediately follows from the form of the amplitude at y = 0. So for the scattering of two nuclei we get the amplitude in terms of purely real quantities In the last formula the vacuum matrix element is implied.
Since H 0 = 0 the term independent of g i and g f vanishes, so that we can also write where F i (y, u) is the operator which creates the evolved initial state. It satisfies the equation with the initial condition The commutation relation (15) allows one to represent and then (20) implies that to find the amplitude one has to substitute u by g f in F i (y, u) At y = 0 this gives the initial amplitude in clear correspondence with the nucleus-nucleus amplitude in the so-called optical approximation. Should we take the final state without reversing the sign of u we would get the sign plus in the exponent in obvious contradiction with the optical amplitude.
Taking the complex conjugate of (20) we find Having in mind that the amplitude is pure imaginary, we see that interchanging the target and projectile leads to the overall change of sign u → −u. However, this will not change the amplitude. Indeed after evolution we shall get function F i (y, −u). But the change u → −u in the final state requires that now we have to substitute u by −g f so that the result will be the same F i (y, g f ) as for the direct transition. So the interchange of the target and projectile does not change the amplitude.

Numerical studies
Calculation of the scattering amplitude reduces to the solution of the differential equation in two variables, y and u which determines evolution in y of the function F(y, u) initially given at y = 0: F(y, u)| y=0 = F 0 (u). Note that apart from the chosen A the amplitude depends on y and two parameters μ and λ. From the form of Eq. (27) it follows that this last three variables are combined in two: the scaled rapidityȳ = μy and ratio ρ = μ/λ. So one can explore the whole domain of rapidities y and values of μ and λ by limitingȳ ≤ȳ max and changing values of ρ appropriately. Equation (27) can be solved analytically only in the case when one drops the term with the second derivative, which describes fusing of pomerons. The remaining equation with only the first derivatives describes propagating pomerons and their consecutive splittings, that is, fan diagrams. Its solution can easily be obtained [11]: For the fan amplitude one chooses the initial state to be a single pomeron F 0 (u) = g i u to find F f an (y, u) = g i ue μy 1 + u ρ (e μy − 1) .

(28)
The amplitude itself is obtained from (28) by putting u = g f and multiplication by i (it is essentially identical to the Schwimmer amplitude (1)).
With the second derivative term included, solution of Eq. (27) gives the complete amplitude with all tree diagrams and loop diagrams taken into account. Note that for nucleusnucleus scattering the set of tree diagrams is much wider than the set of fan diagrams. Unfortunately Eq. (27) in this case cannot be solved analytically (Note, however, some important estimates at very small λ in older papers [4][5][6][7][8][9]). So one is compelled to recur to numerical methods for the solution of Eq. (27). In [11] it was found that the most straightforward approach of evolving the initial function in rapidity by the Runge-Kutta method proved to be quite feasible, provided the step in y is small enough and correlated with the step in u. In our present calculation we find good convergence with y = 5.10 −7 , u = 1.10 −2 and the interval in u taken as 0 < u < 20. Further diminishing of y or u or raising the maxinum value of u have been found to produce no change whatsoever.
Calculations in [11] pursued a somewhat restricted goal to only illustrate the feasibility of the numerical approach and see the limiting behaviour of the propagator and hA amplitude in the limit of very high rapidities. Here we study the A dependence of the hA amplitude having in mind finite nuclei of different atomic numbers and comparison with the results for infinite nuclei. In the standard eikonal picture the effective coupling to the nucleus grows as A 1/3 , which comes from the product AT (b). Accordingly we take for the nucleus g f = A 1/3 and for the nucleon g i = 1 for simplicity. The hA amplitude for the infinite nucleus is then found as explained earlier. At y = 0 we start from F 0 (u) = u, evolve this function according to Eq. (27) and take the final function F(y, u) at u = g f . However, the main purpose of our calculations is to find what will happen when we consider realistic nuclei with finite atomic numbers A. To pass to finite nuclei we change the eikonal amplitude for the infinite nucleus to its standard Glauber form for the finite nucleus with atomic number A: If we change correspondingly the final state in the matrix element for the amplitude then after we evolve the initial function F 0 (u) = u to the desired rapidity the amplitude will be given by For not very small values of A numerical calculations of this expression are hardly feasible due to necessity to find highorder derivatives. So instead we use the discussed symmetry under the interchange of the projectile and target and calculate the inverse amplitude with the initial state represented by the finite nucleus and the final one by the proton, that is, Our calculations were performed for the interval 0 ≤ȳ ≤ 5 of the scaled rapidity for three values of ρ = 10, 2 and 0.5. We recall that greater values of ρ correspond to smaller values of the triple pomeron coupling λ. Taking μ = 0.1 more or less in correspondence with the soft pomeron properties our amplitudes are found at rapidities up to 50. The values of λ studied are then λ = 0.01 for ρ = 10, λ = 0.05 for ρ = 2 and λ = 0.2 for ρ = 0.5.
Our results for A = 8, 27, 64 and 125 are shown in Figs. 1, 2 and 3 for the three mentioned values of ρ, respectively. In each figure we compare the fan amplitude for the finite nucleus with a given A (1), the full amplitude with loops for the finite nucleus (2), the fan amplitude for infinite nucleus Eq. (28) with g f given by A 1/3 (3) and finally the full amplitude with loops calculated for the infinite nucleus (with eikonal initial function) and also with g f = A 1/3 (4).

Discussion
Inspection of our numerical results in Figs. 1, 2 and 3 leads to the following conclusions.
1. As expected, the finite sum of fan diagrams corresponding to the given nucleus with a finite atomic number correctly describes the hA amplitude for very small values of the triple pomeron coupling (ρ = 10) up to a certain value of scaled rapidityȳ <ȳ max after which the result sharply blows up practically to infinity. The limiting rapidityȳ depends on A very weakly, diminishing from 1.5 for A = 8 to 1.1 at A = 125. Of course these features can be immediately read from Eq. (3). Note that at ρ = 0.5 the interval ofȳ where finite fans make some sense is close to zero, so that the corresponding curve (1) is not visible in Fig. 3. 2. Remarkably with loop taken into account the results for finite nuclei and infinite ones practically coincide in all cases except for A = 8 at ρ = 10. Moreover, at comparatively high rapidities they are very weakly dependent on A (but strongly dependent on ρ). For instance for ρ = 2 atȳ = 5 the hA amplitudes for A = 8, 27, 64 and 125 are found to be 1.23, 1.28, 1.31 and 1.34 with the total difference less than 10%. This means that the structure of the hA amplitude with loops taken into account is due mostly to formation of loops during evolution so that the amplitude quickly forgets the initial state. In the exceptional case A = 8 and ρ = 10 the small values of λ and A evidently do not allow one to form enough loops to strongly influence the evolving amplitude. 3. At small values of the triple pomeron coupling λ (ρ = 10) fan diagrams with infinite number of splitting (the Schwimmer formula (1)) describe the total amplitude with loops very well. This is in spite of the fact that with a finite nucleus infinite fans seem to have nothing to do with the physical amplitude. So the analytic continuation involved in extending the validity of the finite fans to infinite ones seems to effectively take into account contributions from loops at small enough λ. This surprising result might have some bearing on the validity of the fan amplitude in general and in the pertubative QCD (BK eqiation) in particular. 4. With the growth of λ also infinite fans cease to describe the amplitude, which can be seen from Figs. 2 and 3.

Conclusion
For realistic nuclei with finite atomic numbers A the standard infinite series of fan diagrams in the Regge theory for hA amplitude converts into a finite sum. With the elementary hadron cross-section rising at large energies this finite sum preserves physical meaning only up to a certain maximal energy, which in fact is not very large and goes down with the rise of A. Inclusion of loop diagrams may cure this situation.
To study this problem we considered the toy Regge model existing in zero transverse dimensions where the loops can be taken into account by numerical methods.
Our results first show that for finite nuclei the model with contribution from loops included gives reasonable results up to very high energies. Second we discover that the found hA amplitude is practically identical to the one which corresponds to the infinite nucleus. In fact at large energies the found amplitude is weakly dependent on the initial amplitude at zero energies, so that it is formed completely from loop contributions.
Finally we found that at small enough triple pomeron coupling the infinite series of fan diagrams gives a good description of the amplitude for finite nuclei in spite of the fact that in reality the series should be cut and the cut series has no physical sense at large energies. The analytic continuation in parameters involved in making the cut series to converge seems to somehow take the loop contribution into account. This conclusion, as mentioned, could explain the success of using infinite fan diagrams in the description of hA amplitudes both in the local Regge theory and the perturbative QCD.