# TRANSVERSE MOMENTUM DEPENDENCE OF CUMULATIVE PIONS This work is supported by the Russian Foundation for Fundamental Research, Grant No. 97-02-18123

## 1 Introduction

A few years ago we have proposed a quark-parton model of cumulative phenomena in the interactions with nuclei  based on perturbative QCD calculations of the corresponding quark diagrams near the thresholds, at which other quarks (”donors”) in the nuclear flucton transfer all their longitudinal momenta to the distinguished active quark and become soft.

Consider the scattering of a hadronic projectile off a nucleus with the momentum $`P`$ in the c.m. system. At high energies the momentum $`K`$ of the produced pion belongs to the cumulative region if $`K_z>P_z/A`$. As a reasonable first approximation, we treat the nucleus as a collection of $`N=3A`$ valence quarks, which, on the average, carry each longitudinal momentum $`x_0P_z/A`$ with $`x_0=1/3`$. In our approach the cumulative pion production proceeds in two steps. First a valence quark with a scaling variable $`x>1`$ is created. Afterwards it decays into the observed hadron with its scaling variable $`x`$ close to the initial cumulative quark’s one. This second step is described by the well-known quark fragmentation functions  and will not be discussed here.

The produced cumulative (”active”) quark acquires the momentum much greater than $`x_0P_z/A`$ only if this quark has interacted by means of gluon exchanges with other $`p`$ quarks of flucton (”donors”) and has taken some of their longitudinal momenta (see Fig.1). If this active quark accumulates all longitudinal momentum of these $`p`$ quarks then $`K_z=(p+1)x_0P_z/A`$ and the donors become soft. It is well-known that interactions which make the longitudinal momentum of one of the quark equal to zero may be treated by perturbation theory . This allows to calculate the part of Fig. 1 responsible for the creation of a cumulative quark explicitly. This was done in , to which papers we refer the reader for all the details. As a result we were able to explain the exponential fall-off of the production rate in the cumulative region.

Since with the rise of $`x`$ the active quark has to interact with a greater number of donors, one expects that its average transverse momentum also grows with $`x`$. Roughly one expects that $`K_{}^2`$ is proportional to the number of interactions, that is, to $`x`$. In  this point was not studied: we have limited ourselves with the inclusive cross-section integrated over the transverse momenta, which lead to some simplifications. The aim of the present paper is to find the pion production rate dependence on the transverse momentum and the mean value of the latter as a function of $`x`$ in the cumulative region. This dependence and also the magnitude of $`K_{}^2`$ have been studied experimentally. The comparison of our predictions with the data allows to obtain further support for our model and fix one of the two its parameters (the infrared cutoff).

## 2 The $`K_{}`$ dependence

Repeating the calculations of the diagram in Fig.1 described in  but not limiting ourselves with the inclusive cross-section integrated over the transverse momentum, we readily find that all dependence upon the transverse momentum $`K_{}`$ of the produced particle is concentrated in a factor:

$$J(K_{})=\rho _A(\underset{p+1}{\underset{}{r,\mathrm{},r}}|\underset{p+1}{\underset{}{\overline{r},\mathrm{},\overline{r}}})G(c_1,\mathrm{},c_p)\underset{i=1}{\overset{p}{}}\lambda (c_ir)\lambda (c_i\overline{r})d^2c_ie^{i(\overline{r}r)K_{}}d^2rd^2\overline{r}$$
(1)   
Here $`\rho _A`$ is the (translationally invariant) quark density matrix of the nucleus:

$$\rho _A(r_i|\overline{r}_i)\psi _A(r_i,r_m)\psi _A^{}(\overline{r}_i,r_m)\underset{m=p+2}{\overset{N}{}}d^2r_m$$
(2)   
where $`\psi _A`$ is the transverse part of the nuclear quark wave function. The propgation of soft donor quarks is decribed by

$$\lambda (c)=\frac{K_0(m|c|)}{2\pi }$$
(3)   
where $`m`$ is the constituent quark mass and $`K_0`$ is the modified Bessel function (the Mac-Donald function). The interaction with the projectile contributes a factor

$$G(c_1,\mathrm{},c_p)=\underset{i=1}{\overset{p}{}}\sigma _{qq}(c_ib_i)\eta _H(b_1,\mathrm{},b_p)d^2b_i$$
(4)   
where $`\sigma _{qq}(c)`$ is the quark-quark cross-section at a given value of impact parameter $`c`$ and

$$\eta _H(b_1,\mathrm{},b_p)=\underset{Lp}{}\frac{L!}{(Lp)!}|\psi _H(b_i)|^2\delta ^{(2)}(\frac{1}{L}\underset{i=1}{\overset{L}{}}b_i)d^2b_{p+1}\mathrm{}d^2b_L$$
(5)   
is a multiparton distribution in the projectile, expressed via the transverse part of its partonic wave function $`\psi _H`$ . If we integrate $`J(K_{})`$ over $`K_{}`$ we come back to our old result (Eq. (33) in ):

$$J(K_{})\frac{d^2K_{}}{(2\pi )^2}=\rho _A(\underset{p+1}{\underset{}{0,\mathrm{},0}}|\underset{p+1}{\underset{}{0,\mathrm{},0}})G(c_1,\mathrm{},c_p)\underset{i=1}{\overset{p}{}}\lambda ^2(c_ir)d^2c_id^2r$$

If one assumes factorization of the multiparton distribution $`\eta _H(b_1,\mathrm{},b_p)`$ then   
$`G(c_1,\mathrm{},c_p)`$ also factorizes:

$$G(c_1,\mathrm{},c_p)=\underset{i=1}{\overset{p}{}}G_0(c_i)$$
(6)   
Following  we use the quasi-eikonal approximation for $`\eta _H`$:

$$\eta _H(b_1,\mathrm{},b_p)=\xi ^{(p1)/2}\nu _H^p\underset{i=1}{\overset{p}{}}\eta _H(b_i)$$

where $`\xi `$ is the quasi-eikonal diffraction factor, $`\nu _H`$ is the mean number of partons in the projectile hadron and the single parton distribution $`\eta _H(b)`$ is normalized to unity. In a Gaussian approximation for $`\sigma (c)`$ and $`\eta _H(b)`$ we find:

$$G_0(c)=\xi ^{\frac{1}{2}\frac{1}{2p}}\frac{\nu _H\sigma _{qq}}{\pi r_{0H}^2}e^{\frac{c^2}{r_{0H}^2}}$$

where $`\sigma _{qq}`$ is the total quark-quark cross-section, $`r_{0H}^2=r_0^2+r_H^2`$, $`r_0`$ and $`r_H`$ are the widths of $`\sigma (c)`$ and $`\eta _H(b)`$ respectively.

With the factorised $`G(c_1,\mathrm{},c_p)`$ (6) we have

$$J(K_{})=\rho _A(\underset{p+1}{\underset{}{0,\mathrm{},0}}|\underset{p+1}{\underset{}{\overline{r}r,\mathrm{},\overline{r}r}})j^p(r,\overline{r})e^{i(\overline{r}r)K_{}}d^2rd^2\overline{r}$$

where

$$j(r,\overline{r})=d^2cG_0(c)\lambda (cr)\lambda (c\overline{r})$$

We also have used the translational invariance of the $`\rho `$-matrix. Note that near the real threshold we have no spectators and

$$\rho _A(\underset{p+1}{\underset{}{0,\mathrm{},0}}|\underset{p+1}{\underset{}{\overline{r}r,\mathrm{},\overline{r}r}})=\rho _A(\underset{p+1}{\underset{}{0,\mathrm{},0}}|\underset{p+1}{\underset{}{0,\mathrm{},0}})$$

In any case large $`K_{}`$ corresponds to small $`\overline{r}r`$ so we factor $`\rho _A`$ out of the integral in zero point. In the rest integral we pass to the variables

$$B=\frac{r+\overline{r}}{2},b=\overline{r}r$$

and shift the integration variable $`c`$, then

$$J(K_{})=\rho _A(\underset{p+1}{\underset{}{0,\mathrm{},0}}|\underset{p+1}{\underset{}{0,\mathrm{},0}})j^p(B,b)e^{ibK_{}}d^2bd^2B$$
(7)   
where

$$j(B,b)=G_0(B+c)\lambda (\frac{b}{2}c)\lambda (\frac{b}{2}+c)d^2c$$
(8)   
## 3 The calculation of $`|K_{}|`$

Now we would like to find the width of the distribution on $`K_{}`$ as a function of $`p`$ or what is the same of the cumulative number $`x=(p+1)/3`$. From the mathematical point of view it is simpler to calculate the mean squared width of the distribution $`K_{}^2`$. Unfortunately in our case this quantity is logarithmically divergent at large $`K_{}`$. This divergency results from the behavior of $`j(B,b)`$ at small $`b`$. This behavior is determined by the behavior of the $`\lambda (b)=K_0(m|b|)/(2\pi )`$ (3), which has a logarithmical singularity at $`|b|=0`$. Smooth $`G_0(B+c)`$ does not affect this behavior.

For this reason we shall rather calculate $`|K_{}|`$:

$$|K_{}|=\frac{1}{J_N}j^p(B,b)|K_{}|e^{ibK_{}}d^2bd^2B\frac{d^2K_{}}{(2\pi )^2}$$
(9)   
where $`J_N`$ is the same integral as in the numerator but without $`|K_{}|`$. Presenting $`|K_{}|`$ as $`K_{}^2/|K_{}|`$ and $`K_{}^2`$ as the Laplacian $`\mathrm{\Delta }_b`$ applied to the exponent we find

$$|K_{}|=\frac{1}{J_N}j^p(B,b)\mathrm{\Delta }_be^{ibK_{}}d^2bd^2B\frac{d^2K_{}}{|K_{}|(2\pi )^2}$$

Twice integrating by parts and using the formula

$$\frac{d^2K_{}}{|K_{}|}e^{ibK_{}}=\frac{2\pi }{|b|}$$

we find

$$|K_{}|=\frac{1}{2\pi J_N}\frac{1}{|b|}\mathrm{\Delta }_bj^p(B,b)d^2bd^2B$$

Now we again integrate by parts once to find

$$|K_{}|=\frac{1}{2\pi J_N}d^2B\frac{d^2b}{|b|^2}(n_b_b)j^p(B,b)$$

where $`n_b=b/|b|`$. This leads to our final formula

$$|K_{}|=\frac{p}{2\pi J_N}d^2B\frac{d^2b}{|b|^2}j^{p1}(B,b)(n_b_b)j(B,b)$$
(10)   
where $`j(B,b)`$ is given by (8), $`\lambda (b)`$ is given by (3) and

$$J_N=d^2Bj^p(B,b=0)$$

## 4 Approximations

To simplify numerical calculations we make some additional approximations, which are not very essential but are well supported by the comparison with exact calulations at a few sample points.

As follows from the the asymptotics of $`K_0(z)`$ at large $`z`$

$$K_0(z)\sqrt{\frac{\pi }{2z}}e^z$$

the width of $`\lambda (b)`$ (3) is of the order $`m^1`$. The function $`G_0`$ is smooth in the vicinity of the origin and its width $`r_{0H}=\sqrt{r_0^2+r_H^2}`$ is substancially larger than the width of $`\lambda `$. For this reason we factor $`G_0(B+c)`$ out of the integral (8) over $`c`$ at the point $`B`$:

$$j(B,b)=G_0(B)\mathrm{\Lambda }(b),\mathrm{\Lambda }(b)\lambda (c)\lambda (cb)d^2c=\frac{|b|}{4\pi m}K_1(m|b|)$$
(11)   
Then we find that the integrals over $`B`$ and $`b`$ decouple

$$J(K_{})=\rho _A(\underset{p+1}{\underset{}{0,\mathrm{},0}}|\underset{p+1}{\underset{}{0,\mathrm{},0}})G_0^p(B)d^2B\mathrm{\Lambda }^p(b)e^{ibK_{}}d^2b$$
(12)   
In this approximation we find that $`|K_{}|`$ depends only on one parameter - the constituent quark mass $`m`$, which in our approach plays the role of an infrared cutoff:

$$|K_{}|=pm_0^{\mathrm{}}𝑑zK_0(z)(zK_1(z))^{p1}$$
(13)   
This allows to relate $`m`$ directly to the experimental data on the transverse momentum dependence.

## 5 Comparison with the data and discussion

The integral in (13) can be easy calculated numerically. For values of $`p=1,\mathrm{},12`$ it is very well approximated by a power dependence (see Fig.2), so that we obtain

$$|K_{}|/m=1.594p^{0.625}$$
(14)   
As we observe, the rise of $`|K_{}|`$ turns out to be even faster than expected on naive physical grounds mentioned in the Introduction ($`\sqrt{p}`$). The resulting plots for $`|K_{}|^2`$ as a function of the cumulative number $`x=(p+1)/3`$ at different values of parameter $`m`$ are shown in Fig.3 together with avaiable experimental data from  on $`K_{}^2`$ for pion production obtained in experiments - with 10 $`GeV`$ protons and  with 8.94 $`GeV`$ protons.

Note that earlier publications of the first group  reported a much stronger increase of $`K_{}^2`$ with $`x`$, up to value 2 $`(GeV/c)^2`$ at $`x=3`$ for pion production. In our approach such an increase would require the quark mass to be as high as $`m225MeV`$. In a more recent publication  the rise of $`K_{}^2`$ is substancially weaker (it corresponds to $`m175MeV`$ in our approach). The authors of  explain this by new experimental data obtained and by a cutoff $`K_{max}`$ introduced in calculations of $`K_{}^2`$ in . The introduction of this cutoff considerably (approximatly two times) decreases the experimental value of $`K_{}^2`$ at $`x=3`$. In our opinion this is a confirmation that the cumulative pion production rate only weakly decreases with $`K_{}`$ in the cumulative region so that the the integral over $`K_{}^2`$ which enters the definition of $`K_{}^2`$ is weakly convergent or even divergent, as in our approach. Undoubtedly presentation of the experimental data in terms of the mean value $`|K_{}|^2`$, rather than $`K_{}^2`$ should reduce the dependence on the cutoff $`K_{max}`$ and make the results more informative.

One of the ideas behind the investigations of the cumulative phenomena is that they may be a manifestation of a cold quark-gluon plasma formed when several nucleons overlap in the nuclear matter. In  we pointed out that our model does not correspond to this picture. It implies coherent interactions of the active quark with donors and, as a result, strong correlations between the longitudinal and transverse motion. Predictions for the dependence of $`|K_{}|`$ on $`x`$ are also different. From the cold quark-gluon plasma model one expects $`|K_{}|`$ to behave as $`x^{1/3}`$, since the Fermi momentum of the quarks inside the overlap volume is proportional to the cubic root of the quark density. Our model predicts a much faster increase, with a power twice larger. The experimental data seem to support our predictions.

## 6 Acknowledgments

The authors are greatly thankful to Prof. P.Hoyer who attracted their attention to the problem.

This work is supported by the Russian Foundation for Fundamental Research, Grant No. 97-02-18123.

Figure captions

The diagram for the production of a cumulative quark with the momentum $`K`$ in the scattering of a projectile hadron with the momentum $`H`$ off a nucleus $`A`$ with the momentum $`P`$. Dashed and chain lines show gluon and pomeron exchanges, respectively.

The $`|K_{}|/m`$ as a function of $`p`$. The points are the results of calculations on (13). The line is the best power fit.

The $`|K_{}|^2`$ as a function of the cumulative number $`x=(p+1)/3`$. The lines are the results of calculations on (13) at different values of parameter $`m`$. The points ($``$) are the experimental data from  on $`K_{}^2`$ for pion production with a cutoff (see the text) obtained in experiments - on a bombardment of nuclei by 10 $`GeV`$ and 9 $`GeV`$ protons. The points ($``$) are the data from the earlier publications of the group  without a cutoff.

M.Braun and V.Vechernin

Fig. 1

M.Braun and V.Vechernin

Fig. 2

M.Braun and V.Vechernin

Fig. 3

M.Braun and V.Vechernin