A new parton model for the soft interactions at high energies

TL;DR

This paper introduces a novel parton model based on Pomeron calculus that effectively describes high-energy soft interactions, ensuring unitarity and accommodating various scattering scenarios, bridging gaps in current QCD approaches.

Contribution

The paper presents a new parton model derived from Pomeron calculus that incorporates impact parameter dependence and unitarity, applicable to multiple scattering processes at high energies.

Findings

Model describes experimental data at high energies.

Ensures both t-channel and s-channel unitarity.

Applicable to hadron-hadron, hadron-nucleus, and nucleus-nucleus scattering.

Abstract

We propose a new parton model and demonstrate that the model describes the relevant experimental data at high energies. The model is based on Pomeron calculus in 1+1 space-time dimensions, as suggested in Ref. [18] and on simple assumptions regarding the hadron structure, related to the impact parameter dependence of the scattering amplitude. This parton model evolves from QCD, assuming that the unknown non-perturbative corrections lead to fixing size of the interacting dipoles. The advantage of this approach is that it satisfies both t-channel and s-channel unitarity, and can be used for summing all diagrams of Pomeron interactions, including Pomeron loops. We can use this approach for all reactions: dilute-dilute (hadron-hadron), dilute-dense (hadron-nucleus) and dense-dense (nucleus-nucleus)for the scattering of parton systems. Unfortunately, we are still far from being able to tackle this problem in the effective QCD theory at high energy (i.e. in the CGC /saturation approach).

Tables

Table 1. Table 1: Fitted parameters.

$\Delta$	$p_0$	m (GeV)	${\chi }^2$/d.o.f.
0.6488 $\pm$ 0.030	0.489 $\pm$ 0.030	0.867 $\pm$ 0.005	1.3

Equations

$$Y\le \frac{2}{{\Delta }_{\text{mbox}BFKL}}\ln \left(\frac{1}{{\Delta }_{\text{mbox}BFKL}^2}\right)$$

$$\mathcal{H}$$

$$K(x,y|z)=\frac{(\mathbf{x}-\mathbf{y}{)}^2}{(\mathbf{x}-\mathbf{z}{)}^2(\mathbf{y}-\mathbf{z}{)}^2},$$

$$[\bar{P}(x,y),P(u,v)]=\gamma (x,y;u,v)$$

$$\gamma (x,y;u,v)=\frac{{\alpha }_s^2}{32{\pi }^2}{\ln }^2\frac{(\mathbf{x}-\mathbf{u}{)}^2(\mathbf{y}-\mathbf{v}{)}^2}{(\mathbf{x}-\mathbf{v}{)}^2(\mathbf{y}-\mathbf{u}{)}^2}$$

$$\frac{d\sigma \left(Y\right)}{dY}=\Delta \sigma \left(Y\right),$$

$$\frac{dN\left(Y\right)}{dY}\,\,=\,\,\Delta\,\Bigg{(}N\left(Y\right)\,\,-\,\,N^{2}\left(Y\right)\Bigg{)}$$

$$⟨m|\bar{n}⟩=\int d\bar{P}\delta (\bar{P})(1-P{)}^m(1-\bar{P}{)}^{\bar{n}}$$

$$⟨m|\bar{n}{⟩}_Y=\int d\bar{P}\delta (\bar{P})(1-P{)}^m{e}^{HY}(1-\bar{P}{)}^{\bar{n}}$$

$${\cal H}_{\rm BK}=-\frac{1}{\gamma}\Big{(}\bar{P}P-\bar{P}P^{2}\Big{)}$$

$$P=-\gamma \frac{d}{d\bar{P}};\gamma \sim {\bar{\alpha }}_S^2>0\leftarrow \text{mbox}Eq.(\text{ref}DDAM)in1+1RFT$$

$$⟨m|e^{{\mathcal{H}}_{\mathrm{BK}}\delta Y}\propto (1-m\delta Y)⟨m|+m\delta Y⟨m+1|$$

$$e^{{\mathcal{H}}_{\mathrm{BK}}\delta Y}|\bar{n}⟩=(1+\bar{n}\delta Y)|\bar{n}⟩-\bar{n}[1+\gamma (\bar{n}-1)]\delta Y|\bar{n}-1⟩+\gamma \bar{n}(\bar{n}-1)\delta Y|\bar{n}-2⟩$$

$${\mathcal{H}}_{\mathrm{B}}=-\frac{1}{\gamma }\left[\bar{P}P-\bar{P}{P}^2-{\bar{P}}^{2}P\right]$$

$$e^{{\mathcal{H}}_{\mathrm{B}}\delta Y}|\bar{n}⟩\approx (1+{\mathcal{H}}_{\mathrm{B}}\delta Y)|\bar{n}⟩=(1-\bar{n}\delta Y)|\bar{n}⟩+\bar{n}\delta Y|\bar{n}+1⟩-\gamma \bar{n}(\bar{n}-1)\delta Y|\bar{n}-1⟩+\gamma \bar{n}(\bar{n}-1)\delta Y|\bar{n}-2⟩$$

$$⟨m|\bar{n}⟩=\sum _{l=0}^m\frac{m!\bar{n}!}{(m-l)!(\bar{n}-l)!l!}(-\gamma {)}^l$$

$$⟨1|\bar{n}⟩=1-\bar{n}\gamma$$

$$⟨1|\bar{n}⟩=\sum _{k=0}^{\bar{n}}\frac{\bar{n}!}{(\bar{n}-k)!k!}(-\gamma {)}^k$$

$$\Big{(}1\,\,-\,\,P\Big{)}\Big{(}1\,\,-\,\,\bar{P}\Big{)}\,\,=\,\,(1-\gamma)\Big{(}1\,\,-\,\,\bar{P}\Big{)}\Big{(}1\,\,-\,\,P\Big{)}$$

$$[P,\bar{P}]=-\gamma +...$$

$$1-\bar{P}=e^{-\ln (1-\gamma )\frac{d}{d\Phi }},;1-P=e^{-\Phi }$$

$$⟨m|\bar{n}⟩=(1-\gamma {)}^{m\bar{n}}$$

$${\cal H}_{\rm B}\,\,=\,\,-\frac{1}{\gamma}\Big{[}(1-\bar{P})P-(1-\bar{P})^{2}P+(1-\bar{P})P^{2}-P^{2}\Big{]}$$

$$e^{{\mathcal{H}}_{\mathrm{B}}\delta Y}|\bar{n}⟩=\left[1+\frac{\delta Y}{\gamma }{\left[1-(1-\gamma {)}^{\bar{n}}\right]}^2\right]|\bar{n}⟩-\frac{\delta Y}{\gamma }\left[1-(1-\gamma {)}^{\bar{n}}\right]\left[2-(1-\gamma {)}^{\bar{n}}\right]|\bar{n}+1⟩+\frac{\delta Y}{\gamma }\left[1-(1-\gamma {)}^{\bar{n}}\right]|\bar{n}+2⟩$$

$$⟨m|e^{{\mathcal{H}}_{\mathrm{B}}\delta Y}=\left[1+\frac{\delta Y}{\gamma }{\left[1-(1-\gamma {)}^m\right]}^2\right]⟨m|-\frac{\delta Y}{\gamma }\left[1-(1-\gamma {)}^m\right]\left[2-(1-\gamma {)}^m\right]⟨m+1|+\frac{\delta Y}{\gamma }\left[1-(1-\gamma {)}^m\right]⟨m+2|$$

$${\mathcal{H}}_{\mathrm{NPM}}=-\frac{1}{\gamma }\bar{P}P$$

$$e^{{\mathcal{H}}_{\mathrm{NPM}}\delta Y}|\bar{n}⟩=\left[1-\frac{\delta Y}{\gamma }[1-(1-\gamma {)}^{\bar{n}}]\right]|\bar{n}⟩+\frac{\delta Y}{\gamma }[1-(1-\gamma {)}^{\bar{n}}]|\bar{n}+1⟩$$

$$P^{†}=\frac{d}{d\Phi }{e}^{\Phi }=\frac{1}{\gamma }\ln (1-\bar{P})\frac{1}{1-P};{\bar{P}}^{†}=-\frac{1}{\gamma }{e}^{-\gamma \frac{d}{d\Phi }}\Phi =\frac{1}{\gamma }\frac{1}{1-\bar{P}}\ln (1-P)$$

$$e^{{\mathcal{H}}_{\mathrm{NPM}}\delta Y}|\bar{n}⟩=\left[1-\frac{\delta Y}{\gamma }[1-(1-\gamma {)}^{\bar{n}}]\right]|\bar{n}⟩+\frac{\delta Y}{\gamma }[1-(1-\gamma {)}^{\bar{n}}]|\bar{n}+1⟩$$

$$\frac{dP}{d\eta}\,=\,\Big{[}{\cal H},P\Big{]}\,;~{}~{}~{}~{}{\rm and}~{}~{}~{}~{}\frac{d\bar{P}}{d\eta}\,=\,\Big{[}{\cal H},\bar{P}\Big{]}$$

Full text

BFKL Pomeron, soft interaction, CGC/saturation approach, correlations

††preprint: TAUP-

1. I Introduction

2. II A new parton model

3. II.1 Violation of the s-channel unitarity in 1 + 1 RFT

4. II.2 Commutators

5. II.3 The Hamiltonian of the new parton model

6. II.4 Equations of motion and the scattering amplitude

7. III The model

8. III.1 Explicit solutions

9. III.2 Non-perturbative structure of hadrons

10. IV Comparison with experimental data

11. IV.1 Fitting $\sigma_{\rm tot},\sigma_{\rm el}$ and $B_{\rm el}$

12. IV.2 $A_{el}(Y,b)$

13. V Conclusions

I Introduction

In our previous papers GLP1 ; GLP2 ; GLMNI ; GLM2CH ; GLMINCL ; GLMCOR ; GLMSP ; GLMACOR we demonstrated that it is possible to build a model, based on the effective QCD theory at high energies: i.e. Colour Glass Condensate (CGC) approach (see Ref.KOLEB for the review). The success of the model emanates from two principle ideas, which are supported by experimental data: (i) typical distances in the soft processes at high energies, turn out to be rather short; and (ii) the CGC approach can be re-written in an equivalent form, as the interaction of BFKL PomeronsAKLL in a limited range of rapidities ( $Y\leq Y_{\rm max}$):

[TABLE]

where $\Delta_{\mbox{\tiny BFKL}}$ denotes the intercept of the BFKL PomeronBFKL . In our model $\Delta_{\mbox{\tiny BFKL}}\,\approx\,0.2-0.25$ leading to $Y_{max}=20-30$, which covers all collider energies.

The equivalence between the CGC approach and the BFKL Pomeron calculus is very important, since (i) it shows that the CGC approach satisfies t-channel unitarity, which is a prerequisite for any effective theory at high energies111We wish to remind the reader that t-channel unitarity follows from two sources: the t-channel unitarity of the BFKL Pomeron, proven in Ref.BFKL , and from the Gribov Reggeon Diagram TechniqueGRIBRC . and (ii) it allows us to consider elastic and diffraction processes, and processes of multiparticle generation, on the same footing. The last statement follows from the AGK cutting rulesAGK , which has been proven for inclusive productionKOTU , but not for correlations KOJA .

In Ref.AKLL it was shown than in the rapidity range of Eq. (1) we can use the MPSI approximationMPSI , which sums the large Pomeron loops for the following Hamiltonian

[TABLE]

where

[TABLE]

$P$ and $\bar{P}$ denote the BFKL Pomeron fields.

The commutation relations are of the form

[TABLE]

where $\gamma(x,y;u,v)$ is the scattering amplitude of a dipole $(x,y)$ on the dipole $(u,v)$

[TABLE]

The above form of the Hamiltonian was suggested in Ref.BRAUN as the natural way for summing all BFKL Pomeron diagrams with only triple Pomeron interactions. Below, we will denote this Hamiltonian by ${\cal H}_{\rm B}$.

In Ref.KLL it is shown that the Hamiltonian of Eq. (2) cannot be the correct one, since it violates s-channel unitarity. In other words, we can use ${\cal H}_{\rm B}$ only for summing the large Pomeron loops in the MPSI approximation, but cannot use it for a general description of the high energy interaction in QCD, nor even for the description of the DIS, which is given by Balitsky - Kovchegov equationBK , and was considered to be the most reliable equation in the framework of the CGC approach. It should be stressed, that in the CGC approach, we do not have a general Hamiltonian that describes the interaction of two dense, or of two dilute systems of partons, although some work in this direction has been doneFOAM . Unfortunately, no progress in formulating the BFKL Pomeron calculus for these systems has been achieved in Ref.KLL .

However, in Ref.KLL such a Hamiltonian is constructed for the 1+1 Reggeon Field Theory, which corresponds to QCD in which the size of the interacting dipoles is fixedMUDI ; LELU . In this model the BFKL equation for the dipole scattering cross section $\sigma$ at a rapidity $Y$ is reduced to

[TABLE]

where $\Delta$ denotes the BFKL intercept. The Eq. (5) reproduces the power-like increase of the cross section with energy, $\exp(\Delta Y)=(1/x)^{\Delta}$. For the scattering amplitude $N$ we obtain the non-linear equation:

[TABLE]

this form is similar to that of the Balitsky-Kovchegov non-linear equationBK . Note that Eq. (6) does not depend on the impact parameter . This fact is also in agreement with the BFKL approach, in which the average momentum of the produced colourless dipoles increases with energy since $p_{T}\propto\,s^{\lambda}$, and Gribov diffusionGRIBDIF in impact parameter with $\Delta b\propto\frac{1}{p_{T}}n\,=\,\frac{1}{p_{T}}\Delta Y$, leads to $\Delta b\to 0$ at high energies. We can consider this simple model as the QCD approach in which the typical size of the colourless dipoles are fixed, and do not depend on energy. In other words, we can view this simple approach as a new parton model for the high energy interactionPARTMOD .

In this model we encountered the same problems with $s$-channel unitarity as in QCD at high energy, however we have identified the Pomeron Hamiltonian which cures all these problemsKLL . In the next section we will briefly review our finding. In section 3 we propose a model based on this Hamiltonian, while in section 4 we compare the predictions of this model with the relevant experimental data. We summarize our results in the Conclusions.

II A new parton model

Below we give a short review of sections 4 and 5 of Ref.KLL , where 1+1 Reggeon Field Theory(RFT) is discussed. We emphasis section 5, which contains the results that we are going to use in building our model.

II.1 Violation of the s-channel unitarity in 1 + 1 RFT

As a direct generalization of the original QCD to 1+1 RFT, the scattering matrix of the projectile consisting of $m$ dipoles on a target consisting of $\bar{n}$ dipoles, is given as by

[TABLE]

The evolution in rapidity takes the form

[TABLE]

The 1+1 analog of the BK evolution (see Eq. (6)) is given by the Hamltonian

[TABLE]

As previously, we take $P$ and $\bar{P}$ to satisfy the dilute limit algebra with the commutator Eq. (3), such that

[TABLE]

To comprehend the origin of $s$-channel unitarity, we consider the infinitesimal evolution ($\delta Y$) of the projectile and target wave functions with the Hamiltonian ${\cal H}_{\rm BK}$:

[TABLE]

[TABLE]

Note the difference between these two equations. The coefficient in front of an $n$-dipole state has the meaning of, the probability to find this number of dipoles in the wave function. The projectile evolution, given by Eq. (11) is unitary: all the probabilities in the evolved state remain positive and smaller than unity, and the sum of the probabilities add up to unity.

On the other hand, the target evolution is non-unitary. Indeed, we face two difficulties with the unitarity: (i) the probability to find the initial state $|\bar{n}\rangle$ after a short interval of evolution exceeds unity; and (ii) the probability to find a state $|\bar{n}-1\rangle$ is negative. The coefficients still sum to unity as for the projectile, but clearly the target evolution violates unitarity.

The Braun Hamiltonian of Eq. (2) takes the following form in 1+1 RFT:

[TABLE]

Exploring the same question about $s$-channel unitarity as in Eq. (11) and Eq. (12) we obtain for the evolution of the target wave function:

[TABLE]

One can see that Eq. (14) violates unitarity, but this violation is smaller than in Eq. (12), being $O(\gamma)$, and it is small for small $\bar{n}$ . However, the coefficient of the term $|\bar{n}\rangle$ is still negative, and becomes large parametrically, long before the saturation limit is reached. Since Eq. (13) is symmetric between the target and the projectile, the projectile evolution now is also non-unitary, and involves negative probabilities. The fact that the violation of the unitarity occurs at large $\bar{n}$ suggets that one reconsiders the commutation relations of Eq. (3) and Eq. (10), which we discuss in the next subsection.

II.2 Commutators

Using the commutators of Eq. (3) and Eq. (10), the scattering matrix can be calculated explicitly (for $m+1<\bar{n}$)

[TABLE]

In particular

[TABLE]

From the above one can see that the amplitudes become negative for $\bar{n}\,>\,1/\gamma$, and when a single dipole of the projectile scatters on several dipoles of the target, our commutation relation does not allow us to account for multiple scattering corrections.

Therefore, as well as curing the problems with unitarity, we need to change Eq. (3)(Eq. (10)) in a such way that the scattering matrix will have the form:

[TABLE]

so as to correctly account for multiple rescatterings.

In Ref.KLL the following commutator which reproduce Eq. (17) is proposed.

[TABLE]

Eq. (18) gives the correct factor $(1-\gamma)^{\bar{n}}$ that includes all multiple scattering corrections, while all the dipoles remain intact, and can subsequently scatter on additional projectile or target dipoles. For small $\gamma$, and in the regime where $P$ and $\bar{P}$ are also small, we obtain

[TABLE]

consistent with our original expression.

Note, that the algebra of Eq. (18) is equivalent to the following representation

[TABLE]

where $\Phi$ and $\bar{\Phi}$ have the same meaning as in Eq. (2).

In the calculation of an amplitude of the type of Eq. (17), once all the factors of $1-\bar{P}$ are commuted through to the left, then in the remaining matrix element $\bar{P}$ operates on the $\delta$-function and thus vanishes. The remaining factors of $(1-P)$ also become unity, since a factor of $\Phi$ is equivalent to a derivative acting on the $\delta$-function, and vanishes when integrated over $\bar{P}$.

With the new algebra we have

[TABLE]

which is a simple and intuitive result: the S-matrix of dipole-dipole scattering raised to the power of the number of dipole pairs that scatter.

It should be stressed, that the modification of the Pomeron algebra is not a matter of choice, but is necessary to obtain the amplitude of Eq. (17), which is unitary for arbitrary numbers of colliding dipoles. However, the question of the unitarity of the evolution is a completely separate one. We will examine the s-channel unitarity for the Braun Hamiltonian, since we plan to use this Hamiltonian for the description of hadron-hadron scattering at high energy.

We write the Braun Hamiltonian in a more convenient form:

[TABLE]

The action on the projectile and the target is obviously symmetric, as the Hamiltonian is self dual under the transformation $P\rightarrow\bar{P}$.

[TABLE]

[TABLE]

This result is a disaster: the evolution of both, projectile and target are non-unitary. In fact the lack of unitarity occurs for an arbitrary number of dipoles $m$ and $\bar{n}$. It is shown in Ref.KLL that the terms that include the four Pomeron interaction, do not help.

II.3 The Hamiltonian of the new parton model

In Ref.KLL a new Hamiltonian is proposed

[TABLE]

where NPM stands for “new parton model”. The fact that it is self dual is obvious. In the limit of small $\bar{P}$ this Hamiltonian reproduces ${\cal H}_{\rm BK}$ ( see Ref.KLL for details). This condition is important for fixing the form of ${\cal H}_{\rm NPM}$.

To check unitarity we consider:

[TABLE]

Performing the operation in Eq. (26), we first express the $P^{\dagger}$ in terms of $\bar{P}$. To do this, recall that $P^{\dagger}$ should annihilate a dipole when acting on the wave function. Using eq.(20) we can write

[TABLE]

For simplicity, in the above equations we have used $\ln(1-\gamma)\approx-\gamma$, since $\gamma\sim\alpha_{s}^{2}\ll 1$

To check for unitarity, we see that

[TABLE]

showing that this evolution is clearly unitary. Due to self duality, it is clear that the evolution of the projectile wave function is unitary as well.

It is interesting that Eq. (28) displays saturation behavior very similar to that expected from QCD evolution, namely at large $\bar{n}$, the change in the wave function is independent of the number of dipoles $\bar{n}$. In the BK approach in QCD, the wave function never saturates(see Eq. (11)), and saturation of the scattering amplitudes is due to multiple scattering effects.

II.4 Equations of motion and the scattering amplitude

The general form of the equation of motion follows from

[TABLE]

For the Hamiltonian ${\cal H}_{\rm NPM}$ we get

[TABLE]

Interestingly, although it is not obvious from the form of the Hamiltonian (see Eq.(22)), the evolution has the same fixed points as in two transverse dimensions, $(0,0),\ (1,0),\ (0,1),\ (1,1)$.

Since the Hamiltonian is conserved, we have

[TABLE]

Using this conservation relation we obtain the following equations of motion:

[TABLE]

It is instructive to note that fixed points $(0,0),\ (0,1)$ and $(1,0)$ are not present in these equations, which means that they are not reachable at $\alpha\neq 0$. The point $(1,1)$ is also not reachable by evolution for $\alpha\neq 1$. Eq. (32) has only two interesting fixed points: $\left(1,\alpha\right)$ and $\left(\alpha,1\right)$. Since for any physical initial condition $P(0)>\alpha$, the asymptotics at $\eta\rightarrow\infty$ is always dominated by the fixed point $(P=1,\bar{P}=\alpha)$, while for $\eta\rightarrow 0$ the point $(P=\alpha,\bar{P}=1)$ is approached.

The general solution to Eq. (30) takes the form:

[TABLE]

where the parameters $\beta$ and $\alpha$ should be determined from the boundary conditions:

[TABLE]

One can see that for $p_{0}\,>\,\bar{p}_{0}$ and $e^{(1-\alpha)Y}\,\gg\,1$ Eq. (34) leads to

[TABLE]

For a symmetric boundary condition $p_{0}=\bar{p}_{0}$, Eq. (34) gives $P(0)=\bar{P}(Y)$, and the solution takes the form

[TABLE]

This solution has a distinct BFKL limit. Indeed, for $\alpha\ll 1$ and $e^{-Y/2}=a\sqrt{\alpha}$ with $1/\alpha\gg a\gg 1$.

We now have the BFKL-like contribution

[TABLE]

Comparing Eq. (37) with Eq. (5) and Eq. (6) one can see that the variable $Y$ in all our formulae is actually equal to $\Delta Y$. The exponential “BFKL-like” growth continues until the Pomeron reaches the value $P(Y)=1$.

The classical solutions determine the scattering amplitude in the classical approximation. The scattering amplitude has a path integral representation. The Pomeron Lagrangian that generates the equations of motion eq.(30) is

[TABLE]

The scattering amplitude is then given by

[TABLE]

In the classical approximation

[TABLE]

where $p(\eta)$ and $\bar{p}(\eta)$ denote the solutions of the classical equations of motion, with the boundary conditions specified in Eq. (40)).

It is interesting to compare the scattering amplitude given by this expression, to the one obtained from the BK equation, which in QCD describes deep inelastic scattering with nuclei. For the latter we have

[TABLE]

In the classical approximation

[TABLE]

Note that the solution for $\bar{P}$ is irrelevant for the BK amplitude, which is determined entirely by $P(Y)$. On the other hand, the scattering amplitude in NPM depends on $\bar{P}$. Nevertheless the two models should be approximately the same in the regime where the BK evolution applies. The results of the estimates in Ref.KLL shows that in the region close to saturation the differences between BK and NPM are quite significant. We will continue the comparison of the two approaches in the following sections, where we construct a realistic model based on the NPM approach.

III The model

To build a model we need to solve two problems:(i) to express parameters $\alpha$ and $\beta$ in Eq. (33) in terms of $p_{0}$ and $\bar{p}_{0}$, and to take the integral over $\eta$ in Eq. (40); (ii) to introduce the non-perturbative structure of hadrons.

III.1 Explicit solutions

We start with solving the first problem, which although it is technical, simplifies the fitting procedure. In the general case we need to replace Eq. (34) and Eq. (35) by the following expressions:

[TABLE]

After doing the integration in Eq. (40) we obtain :

[TABLE]

with

[TABLE]

where $z_{m}=e^{\Delta(1-p_{0})Y}$.

III.2 Non-perturbative structure of hadrons

In this paper we only discuss the non-perturbative corrections, which are related to the impact parameter dependence of the scattering amplitudes. Our assumptions regarding the hadron structure are based on the following features of the Pomeron interactions which stem from the Colour Glass Condensate (CGC/saturation) effective QCD theory at high energies (see Ref.KOLEB for a review):

The BFKL Pomeron exchange occurs at fixed impact parameters. In other words the Green function of the Pomeron $\propto\,\delta^{(2)}\left(\boldsymbol{b}\right)$. The triple Pomeron vertex (see Fig. 1) does not change the impact parameters. We have discussed this property in the introduction. 2. 2.

For nucleus-nucleus collisions the main contributions stems from the ’net’ diagrams of Fig. 1-a BRAUN . In these diagrams the dependence on the impact parameters are concentrated in the vertices of the Pomeron interaction with the nuclei, this dependence has a clear meaning: i.e. the number of nucleons in a nucleus at fixed impact parameter. 3. 3.

For DIS we have the following formula for the total cross section:

[TABLE]

where $Y\,=\,\ln\left(1/x_{Bj}\right)$ and $x_{Bj}$ is the Bjorken $x$. $z$ is the fraction of energy carried by quark. $Q$ denotes the photon virtuality. $N\left(r,Y;b\right)$ is the scattering amplitude for a dipole of size $r$ at impact parameter $b$ which is the solution to the Balitsky-Kovchegov equationBK . $\Psi_{\gamma^{*}}\left(Q,r,z\right)$ is the wave function of the dipole in the virtual photon.

Eq. (47) splits the calculation of the scattering amplitude into two steps: (i) calculation of the wave functions, and (ii) estimates of the dipole scattering amplitude. 4. 4.

The initial condition for the amplitude $N\left(r,Y;b\right)$ can be written as follows GLA ; MUGLA ; MV

[TABLE]

where $N^{\rm BA}\left(r,Q_{T}\right)$ is the amplitude for dipole scattering in the Borm approximation . $T\left(b\right)$ denotes the number of the nucleons at fixed impact parameter $b$.

In Eq. (48), $N^{\rm BA}\left(r,Q_{T}=0\right)\,=\,\frac{\alpha_{S}\pi}{2N_{c}}\pi r^{2}\ln\left(\frac{r^{2}}{R_{N}^{2}}\right)$ while $T\left(b\right)\propto 1/(\pi R^{2}_{A})$ ($R_{N}$ and $R_{A}$ are the radii of nucleon and nucleus, respectively). Therefore, we can re-write $N^{\rm BA}\left(r,Q_{T}=0\right)\,T_{A}\left(b\right)$ in Eq. (48) as $\gamma\,S\left(b\right)$ where $\gamma$ is the dipole scattering amplitude which is equal to $\frac{\alpha_{S}\pi}{2N_{c}}\ln\left(\frac{r^{2}}{R_{N}^{2}}\right)$, and $S_{A}\left(b\right)$ denotes the number of the dipoles with size $r$ that we have in the nucleus at impact parameter $b$. Eq. (48) takes the form

[TABLE]

Comparing Eq. (49) with Eq. (40) we see that we can interpret the initial condition $p(0)=1-\exp\left(-\gamma n\right)$ as $\gamma=\gamma^{\rm BA}\left(r\right)$ and $n=S_{A}\left(b\right)$.

Bearing this in mind, we introduce the following initial conditions for $p(y)$ and $\bar{p}(y)$, considering $\gamma S_{A}(b)\,\ll\,1$:

[TABLE]

choosing

[TABLE]

where $K_{1}$ is the McDonald function (see Ref.RY formula 8.43). Comparing Eq. (50) with Eq. (49) we see that $\gamma=p_{0}$.

The natural generalization of Eq. (47) for diagram of Fig. 1-b is

[TABLE]

which at $Y\to 0$ has the form

[TABLE]

where $n=\,\frac{m^{2}}{4\pi}\int d^{2}b\,S\left(\boldsymbol{b}-\boldsymbol{b}^{\prime},m\right)\,S\left(\boldsymbol{b}^{\prime},m\right)$ is the number of dipoles of size $1/m$ in the two interacting nucleons. It should be stressed that in Eq. (53) we assumed that the typical size of the dipole that is described by the new parton model is $2/m$. In principle, we can introduce this size as a new parameter whose value we will need to determine from comparison with the experimental data.

The expression for $z_{m}=e^{\Delta(1-p_{0})Y}$ sums diagrams of the type shown in Fig. 1-c . Since the triple Pomeron interaction does not induce any impact parameter dependence, we consider $p_{0}=\gamma={\rm Const}$ as a function of $b$.

Combining all the above, we obtain the following equation for the scattering amplitude:

[TABLE]

IV Comparison with experimental data

IV.1 Fitting $\sigma_{\rm tot},\sigma_{\rm el}$ and $B_{\rm el}$

As we have seen in the previous section, we introduce two dimensionless parameters: $\Delta$ - the intercept of the BFKL Pomeron, and $\gamma$ - the amplitude for dipole-dipole scattering at low energies. For $b$-dependence we suggested a specific form for $b$-dependence (see Eq. (50)) which is characterized by the dimensional factor $m$. All three parameters were determined by fitting to the experimental data. We choose to describe three observables: total and elastic cross section and the elastic slope. They have the following expressions through the partial amplitudes:

[TABLE]

From Fig. 2 one can see that we can describe the data for $W\geq 1\,TeV$. The values of parameters are shown in Table 1. Comparing these parameters with the resulting curves in Fig. 2 we note that the shadowing corrections play an essential role. First, the corrections to the Green function of the Pomeron reduce the Pomeron intercept from $\Delta=0.644$ to $\Delta_{\rm dresssed}\,\,=\,\,0.35$. The other shadowing corrections lead to the effective intercept $\Delta_{\rm eff}\approx 0.07$.

From Fig. 2 we see that we fail to describe the experimental data for $W<1\,TeV$. However, we would like to stress that we used a very naive model for the hadron structure. Our previous experienceGLMNI ; GLM2CH shows that we need to take into account the processes of diffraction production, which have been neglected in this model. Our main goal is to demonstrate that the suggested model is able to describe the experimental data at high energies.

IV.2 $A_{el}(Y,b)$

In Fig. 3 we plot the elastic scattering amplitude $A_{\rm el}\left(Y,b\right)$, as a function of the impact parameter. Note that this amplitude has reached the unitary limit 1 at $W=13\,TeV$ and shows the increasing with energy radius, of the interaction. For a comparison, we include in this picture the scattering amplitude, calculated in Balitsky-Kovchegov non-linear equation (see $A_{\rm el}(X)$). This amplitude is far from the unitarity limit and shows only a small increment with increasing energy.

Such behaviour of the scattering amplitude $A_{\rm el}(X)$ reflects the fact that in our approach there is no fixed point (1,0) (or/ and (0,1)), as occurs in the Braun Hamiltonian, and the scattering of hadrons or nucleus cannot be reduced to BK evolution at high energies.

V Conclusions

In this paper we showed that the experimental data at high energies, can be described in the framework of the new parton model. The model is based on the Pomeron calculus in 1+1 space-time, suggested in Ref. AKLL , and on simple assumptions on the hadron structure, related to the impact parameter dependence of the scattering amplitude. This parton model stems from QCD, assuming that the unknown non-perturbative corrections lead to fixing the size of the interacting dipoles. The advantage of this approach is that it satisfies both the t-channel and s-channel unitarity, and can be used for summing all diagrams of the Pomeron interaction, including Pomeron loops. In other words, we can use this approach for all possible reactions: dilute-dilute (hadron-hadron), dilute-dense (hadron - nucleus) and dense-dense (nucleus-nucleus) parton systems scattering. Unfortunately, we are still far from tackling this problem in the framework of QCD effective theory at high energy (CGC /saturation approach).

We achieved quite good descriptions of the three experimental observables: $\sigma_{\rm tot}$,$\sigma_{\rm el}$ and $B_{\rm el}$, especially regarding the energy dependence of these observables. We consider this paper as the first attempt to show that the new parton model can be relevant to the discussion of the experimental data. In spite of the embryonic state of the theory of the quark-gluon confinement, we hope that our model can be viewed as the first step in the right direction regarding the theoretical description of the dilute-dilute parton system scattering at high energy at which, we believe, that such systems become dense.

We are aware that our model is very naive in the description of the hadron structure. We are planning to include diffraction production in our formalism, and to develop a theoretical approach in the framework of the new parton model, so as to be able to treat processes of multiparticle generation.

Acknowledgements.

We thank our colleagues at Tel Aviv University and UTFSM for encouraging discussions. This research was supported by the BSF grant 2012124, by Proyecto Basal FB 0821(Chile) , Fondecyt (Chile) grants 1140842 and 1180118 and by CONICYT grant PIA ACT1406.