Three-quark force in p–p elastic scattering

In the framework of Glauber optical limit approximation, considering the proton has an outer pion cloud of radius ∼0.87fm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 0.87\;fm$$\end{document} and an inner core of radius ∼0.44fm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 0.44\;fm$$\end{document} where the valence three quarks are confined, and including two-gluon exchange three-quark force, a good fit with the experimental data of p–p elastic scattering differential cross section up to q2≈3(GeV/c)2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q^{2}\approx 3\;({\mathrm{GeV}/\mathrm{c})}^{\mathrm {2}}$$\end{document}, total cross section and the ratio of real to imaginary parts of elastic scattering amplitude in the forward direction is obtained at laboratory momenta 200, 290, 500, 1070 and 1500 GeV/c. The radii of two-quark interaction rt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_t$$\end{document} and three-quark force rth\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{th}$$\end{document} are calculated. The quant energy representing the gluon Eg\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_g$$\end{document} is evaluated.


Introduction
Proton-proton scattering at very high energy taking into account the quark structure of proton located is becoming in the centre of the attention of researchers, especially with many measurements of proton-proton scattering cross sections. Different theoretical models are suggested to study different properties of proton structure and of proton-proton scattering mechanism. In most cases the phenomenological models with some physical picture of proton structure are used. Despite the availability of LHC data at very high energy [1][2][3][4], we will focus attention on the relatively low energy of FNAL and CERN-ISR data, where the centre of mass energy √ s = 19.42 − 53 GeV ( p L = 200 − 1500 GeV/c) [5][6][7][8][9]. One of the most prevalent models of proton-proton elastic scattering calculations at high energy is the Regge-pole theory with the Pomeron exchange. In general, in the framework of this theory, a good agreement with the proton-proton elastic scattering differential cross section, total cross section σ t and the ratio ρ of real to imaginary parts of protonproton elastic scattering amplitude in the forward direction was obtained at FNAL and CERN-ISR energy [10][11][12][13][14]. A a e-mail: el_hamam@sci.asu.edu.eg similar fit with the experimental data of the same protonproton quantities at the same energy was obtained by using many different other approaches, for example, the eikonalization approach [15], Bialas-Bzdak model [16], axiomatic quantum field theory [17] and others [18][19][20][21][22].
At higher energies of CDF, HERA and LHC the picture of the theoretical results of proton-proton elastic scattering may be a little different and we may need a new physics [23]. With Regge-pole theory, considering a parametrization of the total cross section σ t and the ratio ρ with two Reggeons and four Pomeron contributions the authors in [24] cannot obtain a good agreement with the data at √ s = 13 TeV. However, in the framework of constituent quark model with pomeron exchange [25] a fit with the CDF and HERA data of protonproton total cross section σ t and elastic cross section σ e was improved including the triple pomeron vertex and with the double pomeron exchange and the authors in [25] show that the radius of the constituent quark R q ≈ 0.0624−0.0882 fm. With the same model including all possible quark interactions and taking quark-quark scattering amplitude in eikonal approximation a good fit with experimental data of protonproton total cross section σ t and elastic cross section σ e at √ s = 23−1855 GeV was obtained [23]. The used quark radius to obtain the fit is R q ≈ 0.0789 fm. At the same time, at √ s = 18000 GeV the authors found R q ≈ 0.3−0.4 fm. In additive quark model with pomeron exchange approach the total cross section σ t , differential cross section dσ dt [26] and the ratio ρ [27] for proton-proton elastic scattering at LHC energies are calculated. In this model with only four order of interaction a good fit with data of dσ dt at √ s = 7 TeV was obtained up to q 2 = 2.5 (GeV/c) 2 with the quark radius R q ≈ 0.44 fm. However, with the all order of quark-quark interactions the agreement with the data was observed up to q 2 = 0.4 (GeV/c) 2 only. Something must be discussed with increasing of q 2 . At LHC energy, for σ t a good fit was obtained [26], while the values of the ratio ρ were little below the data [27].
One of the successful models in the study of proton-proton scattering as composite systems is the multiple scattering theory of Glauber [28]. Using Glauber theory [29][30][31] the authors show that the multiple scattering mechanism with the quark degree of freedom describes the diffraction pattern of particle-particle elastic scattering and quarks having spatial dimensions which are small compared to the corresponding hadrons and a good agreement with the experimental data of high energy was obtained. With simple Gaussian form of proton wave function, by introducing the time-ordering effect in quark-quark multi-scattering the good agreement with the proton-proton data was obtained at CERN-ISR energies [35]. Different wave functions were used in the framework of Glauber model to describe the quark distribution in the proton and obtained a good fit with the experimental data at the CERN-ISR energies [31,32]. In [32], by introducing quark-quark short rang correlations in the wave function obtained the quark radius R q ≈ 0.17 fm. Using the geometrical impact parameter representation [33] obtained R q ≈ 0.15.5 fm. On the other hand, evaluations of the quark radius in the framework of Glauber theory at CERN-ISR energies are of order R q ≈ 0.4 fm [34].
Thus, approximately same value of quark radius R q at different energies, (CERN-ISR and LHC) may be related to different models and then different parametrization of quarkquark interaction. This interpret, also, the similar results of proton-proton cross sections in the framework of different models.
The optical limit approximation in the framework of Glauber theory is not tested for proton-proton elastic scattering. Therefore, in this work, using this approximation, we consider the physical picture of proton as used in Ref. [36], where the proton has an outer pion cloud of radius ∼ 0.87 f m and an inner core of radius ∼ 0.44 f m where the valence three quarks are confined. With this picture of the proton, the proton-proton elastic scattering amplitude can be written as [36] T where T 0 (q) represents the cloud effect, i.e., the interaction of incident proton pion cloud with the target proton pion cloud (soft Pomeron), and F(q) represents the contributions of quark-quark interaction in terms of quark-quark elastic scattering amplitudes (hard Pomeron). The expressions between brackets may explain the common idea between the Pomeron approach and Glauber theory. The pion cloud has an effect in the forward direction of proton-proton scattering. In fact, since q ∼ 1/r , the outer cloud affects at small values of momentum transfer q. At the same time, for the quarks, we have two types of interactions. The first is quarkquark interaction inside the core of the proton and the other is the interaction of a quark of incident proton with a quark of target proton in the case of proton-proton scattering.
With small radius of proton core, r cor = 0.44 fm, the probability of three-quark force in proton-proton scattering cannot be neglected. However, the three-quark force was not studied. We can, simply, introduce the two-gluon exchange three-quark force effect through the profile function of quarkquark interaction in the Glauber theory in similar fashion to the used approaches in the case of 2π -exchange threenucleon force [37,38]. The optical limit approximation [38] is used to study three-quark force in proton-proton elastic scattering at FNAL and CERN-ISR energy. We will compare our results with the experimental data at laboratory momentum in the range 200-1500 GeV/c for the differential cross section of [5][6][7][8][9] and total cross section of [39,40]. We aim to obtain an information about radius of quark-quark interaction and three-quark force radius. Also, using the obtained value of two-quark force radius we aim to evaluate the quant energy representing the mediator of the force between two-quarks.

Formalism
In the optical limit approximation the elastic scattering amplitude F(q) of two colliding protons with three-quark constitutes for each [34], including three-quark force [38] and taking into account the geometrical scaling property of protonproton reaction at CERN-ISR energies [41] is given by, where q 0 = R 0 q, b 0 = b/R 0 are dimensionless variable, q and b are the momentum transfer vector and the two-dimensional impact vector between the mass centres of the two pro- is the inelastic radius of proton and σ inel (s) is the total inelastic cross section for proton-proton collision. In fact, the geometric scaling property is related to quark-quark collision. According to this property, the quark-quark profile function jk (b jk ) of scattering process is replaced by jk (b jk /R 0 ), see Eqs. (4) and (7). The total optical phase shift is [38], where χ t jk (b 0 ) is the two-quark force phase shift, χ th j (kl) (b 0 ) is the three-quark force correction where the quark k and quark l are target quarks and χ th ( jm)k (b 0 ) is the three-quark force correction where quark j and quark m are incident quarks. In the used approach, j(kl) means that the incident quark j interacts with the target quark k and, at the same Table 1 The dimensionless values of α 0 = R 2 0 / r 2 p , two-quark force parameters A 0 , ε and β 0 , three-quark force parameters A, B and γ 0 and cloud parameters B 0 and c 0 time, the target quark k interacts with another target quark l. Also, (jm)k means that the incident quark j interacts with the target quark k and, at the same time, the incident quark j interacts with another quark m in the incident proton. In terms of two-quark and three-quark force profile functions where ϕ I (ρ 0 , r 0 ) and ϕ T (ρ 0 , r 0 ) are the wave functions of the relative motion of the incident and target proton, respectively, with dimensionless variable (ρ 0 , r 0 ) and (ρ 0 , r 0 ), respectively. ρ = R 0 ρ 0 = r 3 − (r 1 + r 2 )/2, r = R 0 r 0 = r 1 − r 2 are the relative coordinates of quarks where r 1 , r 2 and r 3 are the position vectors of the quarks inside the proton. Also, where the sign ⊥ means that the projection is on the impact plane.
The two-quark force profile function Γ t jk (b 0 ) [31] and the three-quark force profile functions th j (kl) (b 0 jk , t 0kl ) and th ( jm)k (b 0 jk , t 0 jm ) [38] are taken as follows where (A 0 jk , ε jk , β 0 jk ) are two-quark force dimensionless parameters. The quantity ε jk is the ratio of the real to imaginary parts of the quark-quark amplitude in the forward direction and β 0 jk is the slope parameter and is related to the radius of two-quark interaction. Note that, from the profile function Γ t jk (b 0 jk ), we can take r t = 2β jk as a radius of quark-quark interaction, where β jk = β 0 jk R 2 0 . The parameters (A j (kl) , B j (kl) , A ( jm)k , B ( jm)k , γ 0 ) are the three-quark force dimensionless parameters. Also, r th = √ γ , where γ = γ 0 R 2 0 , can be considered as the three-quark force radius [37,38]. From the uncertainty relation of energy and time, we know that r th = 1 2 r t . Therefore, we can write √ γ = 1 2 2β jk = β jk 2 . Thus, we can take γ = β jk 2 . The wave functions of incident and target protons are taken in the same Gauss form, where α = 1 r 2 p [35]. The value of r 2 1/2 p is taken as given in the Refs. [31,35], r 2 1/2 p ≈ 0.59 fm. This value is less than the charge radius of proton, ≈ 0.87 f m. In fact, the wave functions (10) describe the motion of the three quarks inside the core of the proton. Therefore, the root mean square radius r 2 1/2 p in this case has a value less than the value of electromagnetic radius of the proton where the outer pion cloud is included in the measurements. In Ref. [25] the value r 2 1/2 p = 0, 68 fm was used and the difference with the electromagnetic radius interpreted by the difference in the distribution of charged and strong interacting matter inside the proton.
Thus, the total proton-proton optical phase shift χ opt with three-quark force correction, from Eqs. (4)-(10), is Fig. 1 The p-p elastic scattering differential cross section. The results for 200, 290, 500, 1070 and 1500 GeV/c are given in figures (a-e), respectively. The dashed, dotdashed and solid curves represent the results of (two-quark force) only, (two-quark force + pion cloud effect), and (two-+three-quark force+ pion-cloud effect), respectively. The experimental data are taken from [5][6][7][8][9] where, Since the experimental evaluation of quark-quark scattering amplitude parameters is not available and for simplicity we will assume that A 0 jk = A 0 , ε jk = ε and β 0 jk = β 0 for all j and k. Also, we assume that A j(kl) = A ( jm)k = A and B j(kl) = B ( jm)k = B for all j, k, l and m. Thus, with these approximations the total proton-proton optical phase shift χ opt becomes, To obtain a good fit in the forward direction of proton-proton elastic scattering differential cross section we will take into account the pion-cloud effect. In Ref. [31] this pion-cloud effect is represented by an additional phenomenological term which takes the form where c 0 and B 0 are dimensionless constants. Therefore, in the calculations of p-p elastic scattering differential cross section we will use the equation [35] dσ For the total cross section we will use The ratio ρ pp of real part to imaginary part of the elastic scattering amplitude in the forward direction is given by

Results and discussion
The p-p elastic scattering differential cross section at the laboratory momenta 200, 290, 500, 1070 and 1500 GeV/c is calculated including two-gluon three-quark force and pion cloud effect. The experimental data are taken from [5] for 200 GeV/c and from [6][7][8][9] for 290, 500, 1070 and 1500 GeV/c. The values of the quark-quark parameters which are used in the calculations are given in Table 1. The results of the differential cross section are given in the Fig. 1. It is clear, from the Fig. 1, that the optical limit approximation with two-quark force only cannot describe the differential cross section at all the used energies, dashed curves in the figure.
The values of theoretical results are clearly below the experimental data. This may be due, partially, to the absence of the multi-scattering terms in this approximation. The inclusion of pion-cloud effect leads to an improvement of the results in the forward direction. We obtained a good fit with the experimental data up to q 2 ≈ 3 (GeV/c) 2 by including the three-quark force effect. These results are obtained with quark-quark radius of interaction (in average) is r t = √ 2β = 0.738 fm and the three-quark force radius (in average) is r th = √ γ = 0.364 fm = (1/2)r t , see Table 2. The obtained value for strong force radius between two quarks inside the proton is consistent with the measured values of the proton radius, r p = 0.84 − 0.87 fm [42].  Table 1, respectively Table 3 Total cross section with (two-quark force, σ t tpp ), with (twoquark force +pion cloud, σ t+c tpp ) and with (two-quark force +pion cloud+three-quark force, σ t+c+th tpp ). The ratio of real to imaginary parts of p-p elastic scattering amplitude in the forward direc-tion with(two-quark force, ρ t pp ), with (two-quark force +pion cloud, ρ t+c pp ) and with (two-quark force +pion cloud+three-quark force, ρ t+c+th pp ). The experimental data are taken from [39,40]  From the uncertainty relation of time and energy t E ≥ h and using the obtained radius of quark-quark interaction we get the quant energy representing the mediator of the force between two quarks (gluon) E g ≈ 271.695 MeV in average. E 2g = 2E g ≈ 543.39 MeV, in the Table 2, represent twogluon exchange in the case of three-quark force. We must note that the values of each of quantities r t , r th ,E g have the same order at different considered energies. Therefore, we take the average of all cases.
To confirm the obtained radius of the two-and three quark forces and other parameter values, the p-p total cross section σ t pp and the ratio ρ pp are calculated. The results are given in the Table 3. The contributions of pion-cloud and three-quark force are very clear and are important for the two quantities. Inclusion of pion cloud improves the results, but with threequark force the results match with experimental data for the total cross section and for the ratio the results are in good agreement with the experimental data.
Since the rest mass of three quarks inside the proton is considered a small part of proton mass, and the other part of the proton mass is coming from the field particles inside the proton (gluons) and the field particle of the external strong field of proton (π − meson), we will define the residual energy of proton as follows Since E π ≈ 1 2 E g , this interprets the absence of the quark outside of proton, especially if we note that the three-quark force, with small radius r th = 0.364 fm, takes an important place inside the core of the proton of radius 0.4 fm and each quark is tied with the other two quarks by two gluons and E π ≈ 1 4 E 2g . By using the average values of the parameters in the Table 1 to calculate the elastic scattering differential cross section, Fig. 2, we have approximately the same results as in the Fig. 1 with a little difference at the first minimum. Therefore, we can consider that two-quark force, three-quark force and pion-cloud parameters are slightly dependent on energy.

Conclusions
In conclusion, taking into account the pion-cloud and threequark force effects, by using the optical limit approximation a good fit with the experimental data of p-p elastic scattering differential cross section up to q 2 ≈ 3 (GeV/c) 2 , total cross section and the ratio of real to imaginary parts of elastic scattering amplitude in the forward direction is obtained at the laboratory momenta 200, 290, 500, 1070 and 1500 GeV/c. The radii of two-quark interaction r t and three-quark force r th are calculated. The quant energy representing gluon E g and the mass of the field particle of strong interaction outside the proton m π are evaluated.