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 looses 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 surprizing 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 wellknown 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 essentially reduces to the exchange of two pomerons ang grows correspondingly at high energies.
This problem is difficult to study in the framework of the fully-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 A.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 normalized 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 phyiscal) 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 Schwmimmer 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 atttractive, since it is explicitly unitary. However expression (8) is unitary only while |a(y)T (b)| < 1. This is always so when the proton-proton scattering amplituide 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 looses 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 formula (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 zerodimensional 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]- [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 to numerically sum all contributioins 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 propreties of the toy model indispensable for our study, referring the reader either to old papers [4]- [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 interestng case when µ > 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-Hermithean to each other and satisfy the commutation relation Operators u and v have the meaning of creation and annihilation operators of the pomeron respectively. The vacuum state Ψ 0 , normalized 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 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 function F (y, u) initally 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 the Eq. (27) it folows 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 anaytically 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 be easily obtained to be [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 .
The amplitude itself is obtained from (28) by putting u = g f and multiplcation by i (it is essentially identical with 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 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]- [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 fuction in rapidity by the Runge=Kutta method proved to be quite feasible, provided the step in y is small enough and correlated wth 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 realisitic 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 high-order 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 in more or less corresondence 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-3 leads to the following conclusions. 1) As expected, the finite sum of fan diagrams corresponding to the given nucleus with a finite atomic numver 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 have 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  Figure 2: Same as in Fig. 1 for ρ = 2 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 exeptional case A = 8 and ρ = 10 the small values of λ and A evidently do not allow 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 surprizing 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. With ρ = 0.5 atȳ = 5 they overestimate the amplitude by 3 orders of magnitude. 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 with 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.

Conclusion
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.

Acknowledgments
The author thanks the Bologna Physics department and INFN for hospitality. He also thanks G.P.Vacca for very interesting and constructive discussions. The author acknowledges Saint-