arxiv:hep-ph/ v1 5 Mar 2003

Transcription

1 NEXT-to-LEADING ORDER CONSTITUENT QUARK arxiv:he-h/333v 5 Mar 23 STRUCTURE AND HADRONIC STRUCTURE FUNCTIONS Firooz Arash (a,b,c) and Ali N. Khorramian (c,d) (a) Center for Theoretical Physics and Mathematics, AEOI, P.O. Box , Tehran, Iran (b) Physics Deartment, AmirKabir University (Tafresh Camus), Hafez Ave.,Tehran Iran 9834 (c) Institute for Studies in Theoretical Physics and Mathematics P.O.Box , Tehran, Iran (d) Physics Deartment, Semnan University, Semnan, Iran, January 23, 28 Abstract We calculate the artonic structure of a constituent quark in the Next-to-Leading Order framework. The structure of any hadron can be obtained thereafter using a convolution method. Such a rocedure is used to generate the structure function of roton and ion in NLO, neglecting certain corrections to Λ QCD. It is shown that while the constituent quark structure is generated urely erturbatively and accounts for the most art of the hadronic structure, there is a few ercent contributions coming farash@cic.aut.ac.ir khorramiana@theory.im.ac.ir

2 from the nonerturbative sector in the hadronic structure. This contribution lays the key role in exlaining the SU(2) symmetry breaking of the nucleon sea and the observed violation of Gottfried sum rule. These effects are calculated. We obtained an Excellent agreement with the exerimental data in a wide range of x = [ 6,] and Q 2 = [.5,5]GeV 2 for the roton structure function. We have also calculated Pion structure and comared it with the existing data. Again, the model calculations agree rather well with the data from exeriment. PACS Numbers: 3.6 Hb, x, e, 2.2.Fv INTRODUCTION Our knowledge of hadronic structure is based on the hadron sectroscoy and the Dee Inelastic Scattering (DIS) data. In the former icture quarks are massive articles and their bound states describe the static roerties of hadrons. On the other hand, the interretation of DIS data relies uon the quarks of QCD Lagrangian with very small mass. The hadronic structure in this icture is intimately connected with the resence of a large number of artons (quarks and gluons). Mass is not the only difference between these two tyes of quarks; they also differ in other imortant roerties. For examle, the color charge of quark field in QCD Lagrangian is ill-defined and is not gauge invariant, reflecting the color of gluons in an interacting theory; whereas, color associated with a constituent quark (CQ) is a well defined entity. It is shown that [][2] one can erturbatively dress a QCD Lagrangian field to all orders and construct a CQ in conformity with the color confinement. From this oint of view a CQ is defined as a quasi-article emerging from the dressing of valence quark with gluons and q q airs in QCD. Of course, the concet of CQ as an intermediate ste between the quarks of QCD Lagrangian and hadrons is not new. In fact, Altarelli and Cabibo [3] used this concet in the context of SU(6)XO(3). R.C. Hwa in his work used the term valon,for them, elaborated on the concet and showed its alications to many hysical rocesses [4]. In Ref.[2] it 2

3 is suggested that the concet of dressed quark and gluon might be useful in the area of jet hysics and heavy quark effective theory. Desite the ever resence of CQ no one has calculated its content and artonic structure without resorting to hadronic data and the rocess of deconvolution. Thus, the urose of this aer is threefold: (a) to evaluate the structure of a CQ in the Next-to-Leading Order (NLO) framework of QCD, (b) to verify its conformity with the structure function data of nucleon and ion for which amle data are available, (c) to include additional refinments that are needed to account for the violation of Gottfried Sum Rule (GSR) and the binding effect of constituent quarks to form a hysical hadron. 2 FORMALISM 2. Moments of Parton Distribution in a CQ In this subsection we will utilize the extended work done on the develoment of the NLO calculations of moments, to calculate the structure of a constituent quark. What will follow in the rest of this subsection is not a new next-to-leading order calculation. However, we find it interesting to exlore the existing calculations in the valon framework. By definition, a CQ is a universal building block for every hadron; that is, the structure of a CQ is common to all hadrons and is generated erturbatively. Once its structure is calculated, it would be ossible, in rincile, to calculate the structure of any hadron. In doing so, we will follow the hilosohy that in a DIS exeriment at high enough Q 2, it is the structure of a CQ which is being robed, and that at sufficiently low value of Q 2 this structure cannot be resolved. Thus, a CQ behaves as a valence quark, and a hadron is viewed as the bound state of its constituent quarks. Under these criteria, artons of DIS exeriments are comonents of CQ. The structure function of a U-tye CQ at high Q 2 3

4 can be written as[4]: F2 U (z,q 2 ) = 4 9 z(g u +G U ū )+ U 9 z(g d U +G d U +Gs U +G s )+... () U where all the functions on the right-hand side are the robability functions for quarks havingmomentumfraction z ofau-tyecqat Q 2. Asimilarexressioncanbewrittenfor a D-tye CQ. Following Ref. [4], we define the singlet (S) and nonsinglet (NS) comonents of the CQ structure functions as: G S = f (G q i CQ i= +G q i ) = G f +(2f )G uf (2) CQ G NS = f (G q i CQ i= G q i ) = G f G uf (3) CQ where G f is the favored distribution describing the structure function of a quark within a CQ of the same flavor, while the unfavored distribution, G uf, describes the structure function of any quark of different flavor within the CQ. We have used, f, for the number of active flavors. Solving (2) and (3) for G f and G uf, we get, G f = 2f (GS +(2f )G NS ), (4) G uf = 2f (GS G NS ). (5) Having exressed all the structure functions of a CQ in terms of G f and G uf, we now go to the N- moment sace and define the moments of these distributions as: M 2 (N,Q 2 ) = M i (N,Q 2 ) = x N 2 F 2 (x,q 2 )dx (6) x N G i (x,q 2 )dx (7) where, the subscrition i stands for S or NS. Due to charge symmetry, in the following, we only write the CQ distribution for roton. In the NLO aroximation the deendence of the running couling constant, α s, on Q 2 is given by: α s (Q 2 ) 4π β ln( Q2 ) ( β lnln( Q2 ) Λ 2 ), (8) β Λ 2 2 ln(q2 ) Λ 2 4

5 with β = 3 (33 2f) and β = 2 38f 3. The evolution of arton distributions is controlled by the anomalous dimensions: γns N = αs 4π γ()n qq +( αs 4π )2 γ ()N NS γij N = αs 4π γ()n ij +( αs 4π )2 γ ()N ij where i,j = q,g; and γ (,)N are one and two loo anomalous dimension matrices. The nonsinglet moments in the NLO are as follows: M NS (N,Q 2 ) = [+ α s(q 2 ) α s (Q 2 ) 4π β γ qq ()N )]( α s(q 2 ) 2β ( γ()n NS 2β 2 ()N γqq α s (Q 2 )) 2β. (9) The anomalous dimension matrices, γ (,)N,are given by, γ (,)N = γ (,)N qq γ (,)N gq γ (,)N qg γ (,)N gg () These matrices, govern the moments of the singlet quark and sector as given by: M S (N,Q 2 ) = {( α s(q 2 λ N ) M G (N,Q 2 α ) s (Q 2 )) 2β [ N α s (Q 2 ) α s(q 2 ) N 2β 4π γn N ( α s(q 2 ) 4π α s(q 2 ) 4π (α s(q 2 λ N ) + λn α s (Q 2 )) 2β ) N γn N + 2β +λ N + λ N ]}+{(Q + )} () where γ N = γ ()N β β γ ()N and γ ()N, γ ()N are one and two-loo anomalous dimension matrices, resectively. is theunitmatrix. Notice that thelast term ineq. (), Q +, means that all subscrits are exchanged. The leading order behavior is obtained from the first term in the square brackets. λ N ± denote the eigenvalues of the one-loo anomalous dimension matrix, γ ()N : λ N ± = 2 [γ()n qq +γ gg ()N ± (γ gg ()N γ qq ()N ) 2 +4γ qg ()N γ gq ()N ] (2) These anomalous dimensions and γ (,)N kl are calculated u to second order in α s and can be found in [5] [6]. N ± are as follows: N ± = ±(γ()n λ N )/(λn + λn ) (3) InRef. [4], d () ± is used to denote the leading order anomalous dimensions; these are related to λ ± by: d ± = λn ± 2β. 5

6 In evaluating these moments, we have taken Q 2 =.283 GeV 2 and Λ =.22 GeV, as our initial scales. It seems that evolution of arton distributions from such a low value of Q 2 is not justified theoretically; we shall take u this issue shortly in subsection Parton Densities In CQ In the revious subsection we outlined how to calculate the moments. In this subsection we will give the corresonding arton distributions in a CQ. Themoments of thecqstructurefunction, F CQ 2 (z,q 2 ), areexressed comletely interms of Q 2 through the evolution arameter t: The moments of valence and sea quarks in a CQ are, Q2 ln Λ t = ln 2 (4) ln Q2 Λ 2 Mvalence CQ = M NS (N,Q 2 ) (5) Msea CQ = 2f (MS M NS ) (6) wherem S,NS are given by eqs. (9) and (). It is now straightforward, to evaluate Mvalence CQ and Msea CQ at any Q2 or t. These are calculated numerically for a range of N at fixed Q 2 or t. In Figure () we resent these moments both in the leading and next-to-leading orders. For every value of t we fit the moments by a beta function, which in turn is the moment of the arton distribution as given below, zq val. (z,q 2 ) = az b ( z) c (7) CQ zq sea CQ (z,q2 ) = αz β ( z) γ [+ηz +ξz.5 ] (8) Thearametersa,b, c, α, β, γ, η, andξ arefunctionsofq 2 throughtheevolutionarameter t. The same form as in Eq. (8) is obtained for the gluon distribution in a CQ, only with different arameter values. The functional forms of these arameters are olynomials of order three in t and are given in the aendix. The goodness of these we fits are checked 6

7 by χ 2 minimization rocedure. We find that, for all arameters a, b, c, and d, of the valence quark sector χ 2 /DOF =.99, with the standard error of order 3. For the sea quark and gluon sector, χ 2 /DOF carries between.9 and.98. Of course, the forms (7) and (8) are not unique; rather, they are the most simle and commonly used forms. We note that the sum rule reflecting the fact that each CQ contains only one valence quark is satisfied at all values of Q 2 : qval. (z,q 2 )dz =. (9) CQ Substituting these results into Eq. () comletes the evaluation of CQ structure function in NLO. In Figure (2), we lot various arton distributions inside a CQ. 2.3 discussion In subsection 2. we indicated that our initial scales are Q 2 =.283 GeV 2 and Λ =.22 GeV. The above value of Q corresonds to a distance of.36 fm which is roughly equal to or slightly less than the radius of a CQ. It may be objected that such distances are robablytoolargeforameaningfulureerturbativetreatment. WenotethatF CQ 2 (z,q 2 ) has the roerty that it becomes δ(z ) as Q 2 is extraolated to Q 2, which is beyond the region of validity. This mathematical boundary condition signifies that the internal structure of a CQ cannot be resolved at Q in the NLO aroximation. Consequently, when this roerty is alied to Eq.(2), below, the structure function of the nucleon becomes directly related to xg CQ(x) at those values of Q ; that is, Q is the leading P order effective value at which the hadron can be considered as consisting of only three (two) CQ, for baryons (mesons). We have checked that when Q 2 aroaches to Q 2 the quark moments aroach to unity and the moments of gluon aroach to zero. In fact for Q 2 =.2839 gluon moments are at the order of 4 and singlet and nonsinglet moments are In fact our results are only meaningful for Q 2.4 GeV 2. From the theoretical standoint, both Λ and Q deend on the order of the moments, N; 7

8 buthere, wehaveassumedthattheyareindeendentofn. Inthisway, wehaveintroduced some degree of aroximation to the Q 2 evolution of the valence and sea quarks. However, on one hand there are other contributions like target-mass effects, which add uncertainties to the theoretical redictions of erturbative QCD, while on the other hand since we are dealing with the CQ, there is no exerimental data to invalidate an N indeendent Λ assumtion. 3 HADRONIC STRUCTURE 3. Constituent Quark Distribution in hadron So far, the structure of a CQ is calculated in the NLO framework. In this section we will use the convolution method to calculate the structure functions of roton, F 2 (x,q2 ), and ion. Let us denote the structure function of a CQ by F CQ 2 (z,q 2 ) and the robability of finding a CQ carrying momentum fraction y of the hadron by GCQ(y). The corresonding h structure function of a hadron is obtained by convolution of F CQ 2 (z,q 2 ) and GCQ(y), F h 2 (x,q 2 ) = CQ x dygcq(y)f CQ 2 ( x h y,q2 ) (2) where summation runs over the number of CQ s in a articular hadron. We note that GCQ h (y) is indeendent of the nature of the robe, and its Q 2 value. In effect GCQ(y) h describes the wave function of hadron in CQ reresentation, containing all the comlications due to confinement. From the theoretical oint of view, this function cannot be evaluated accurately. To facilitate henomenological analysis, following [4], we assume a simle form for the exclusive CQ distribution in roton and ion as follows, h G UUD/ (y,y 2,y 3 ) = l(y y 2 ) m y n 3 δ(y +y 2 +y 3 ) (2) G ŪD/π (y,y 2 ) = qy µ yν 2δ(y +y 2 ) (22) 8

9 where l and q are normalization arameters. After integrating over unwanted momenta, we can arrive at the inclusive distribution of individual CQ: G U/ (y) = G D/ (y) = B(m+,n+m+2) ym ( y) ( m+n+), (23) B(n+,2m+2) yn ( y) 2m+, (24) G Ū/π (y) = B(µ+,ν +) yµ ( y) ν. (25) Similar exression for G D/π is obtained with the interchange of µ ν. In the above equations, B(i, j) is the Euler Beta function. The arguments of this function, as well as, l and q of Eqs.(2) and (22), are fixed by the number and momentum sum rules: GCQ(y)dy =, h CQ ygcq(y)dy =, (26) h where CQ = U,D,Ū and h =,π. Numerical values are: µ =., ν =.6, m =.65 and n =.35. These arameters are indeendent of Q 2 and are given in [4] and [7] for roton and ion, resectively. In [8], a new set of values for m and n are suggested which differ significantly from the values quoted above. Reference [8] uses the CTEQ arton distributions in the Next-to-Leading order to fix these arameters. The authors of Ref. [8] have found yet another set of values in [9], which differs from those given in [4] and [8]. They attribute these differences to the different theoretical assumtions. We have tried a range of values for m, n, ν and µ, and checked that the sum rules (26) are satisfied. Because of the Beta function in the denominator of Eqs.(23-25) the number sum rule does not change by changing m and n and we found little sensitivity on the momentum sum rule. In Fact, we have varied both m and n, randomly in the interval [.2,.95], and checked it against F 2 data. We found very little sensitivity to the choices made. For the data that we shall analyze, these two arameter formulas will rove to be adequate. The differences between arameters of Ref. [8] and Ref. [4] (which we have adoted) can be attributed to the different methods of extracting the arameters. We have calculated 9

10 the arton distributions from the erturbative QCD, whereas Ref.[8] have used CTEQ s global fit to extract arton content of nucleon. The considerations on the form of CQ distribution that are used here, does not exclude other ossibilities. However, in our model the quark distribution, say, in a roton is described as the convolution of the CQ distribution, G CQ, and the quark distribution in the CQ. In the moment form the latter two become the roduct of their moments. Mathematically, these two moments can be modified without changing the roduct. That would, of course, modify the original G CQ and arton distribution in a CQ. But the arton distribution has a definite meaning. It is the evolution of a quark as Q increases. At Q = Q, arton distribution is a delta function, imlying that the CQ has no internal structure that can be robed. Our model is to say that the three valence quark at high Q becomes three CQ s at low Q without anything else. Such a hysical statement of the model does not ermit to redefine GCQ criticism still remains. or the arton distribution in a CQ, although mathematical Our choice of GCQ is sufficient to fit all the data of dee inelastic scattering in the range of interest. It does not mean that a more comlicated GCQ is ruled out. It is a subject that can be revisited in the future when roblem arise. At this oint our formula for the CQ distribution is satisfactory. A comment regarding the moments is in order: The definition of the moments are not arbitrary if we are to attach hysical interretation to the constituent quark distribution. There are only two U and one D in a roton, which imose constrains on the first moment. Also, the sum of the CQ momenta should be one, which imoses a constraint on the second moment. If one writes the valence quark moments as a roduct of G n and M NS then M NS must have connection with QCD evolution to ensure that there is only one valence quark in every CQ and hence, the constraint becomes M NS =. With these constrains, one cannot freely redefine moments of CQ distribution, M NS and M S. yet, there are infinite number of CQ moments and we have cited only two constraints. Does this leave an infinite horde number of undefined quantities? The higher

11 moments of the CQ distribution are affected mainly by the large x behavior of roton wave function, or will affect only the even larger x behavior of the arton distribution. Until one focuses on that region and finds inconsistency, one chooses the simlest arameterization. There are similar situations both in classical and contemorary hysics. The question can be turned around and asked why did Newton sto at r 2 law for gravitational force? Could he have excluded the ossible existence of r n terms with n = 3,4,5,...? By only citing the existence of two constraints, Newton avoided confronting the infinite horde of undefined quantities. Our model is inductive hysics. Valons or CQ s are not derived by deductive logic from the first rinciles. In modeling, one finds the simlest way to cature the essence of the hysics involved. The situation is also resent in contemorary hysics: we could have asked why the Cornell otential for the bound-state roblem of quark and antiquark does not contain all ossible ower law terms beside the r and the r terms. The answer is simle: it is what fits the article mass sectrum that counts, again in the simlest way ossible that catures the essence of the roblem. 3.2 Parton distribution in hadrons Having secified the arton distributions in the CQ and the constituent quark distribution in roton and ion, it is ossible to determine various arton distributions in any hadron. For roton, we get: dy q val./ (x,q 2 ) = 2 x y G U/(y)q val./u ( x y,q2 ) dy + x y G D/(y)q val./d ( x y,q2 ) = u val./ (x,q 2 )+d val./ (x,q 2 ) (27) q sea/ (x,q 2 dy ) = 2 x y G U/(y)q sea/u ( x dy y,q2 )+ x y G D/(y)q sea/d ( x y,q2 ) (28) The above equations reresent the contribution of constituent quarks to the sea and the valence quark distributions in roton. In Figure (3), the CQ contribution to arton

12 distributions in and π is shown. Comaring with roton structure function data, shows that the results fall short of reresenting the exerimental data by a few ercent (3%-8%, deending on Q 2 and x). This is not surrising, for it is a well established fact that there are soft gluons in the nucleon. We attribute their resence to the fact that constituent quarks are not free in a hadron, but interact with each other to form the bound states. At each Q 2 value we have comuted the soft gluon contribution using direct comarison with exerimental data. The final result is arameterized as: xg soft = R(Q 2 )x.29 ( x) 2.79 ; (29) where, R(Q 2 ) = t t 2.82t 3. In figure (4), the shae of xg soft is resented for a range of Q 2 values. Slitting of these gluons could in turn roduce q q airs which can combine with the original constituent quarks and fluctuate the baryonic state to a meson-baryon state. This also will break the SU(2) symmetry of the sea as we shall see next. 3.3 The Role of Soft Gluon In our model there is no room in the CQ structure for the breaking of SU(2) symmetry of the nucleon sea. But, after creation of q q air by the soft gluons, these quarks can recombine with CQ to fluctuate into meson-nucleon state which breaks the symmetry of the nucleon sea. In what follows we will comute this comonent, and its contribution to the nucleon structure function, and its relevance to the violation of the Gottfried sum rule. We will follow the rescrition of [] and []. The above roosed model, is similar to the effective chiral quark theory, in which the relevant degrees of freedom are constituent quarks, gluons and Goldstone bosons. The couling of ions roduces ion cloud of the constituent quark, in which we are interested. The oint we would like to make is that, there are two tyes of sea quark and gluon that contribute to the nucleon sea measured 2

13 in DIS. One tye is generated from QCD hard bremsstrahlung and gluon slitting. This comonent is discussed in subsection 2.2 and is associated with the internal comosition of the constituent quark rather than the roton itself. The other contribution is generated nonerturbatively via meson-baryon fluctuation. As we will show, according to our results, this comonent rovides a consistent framework to understand the flavor asymmetry of the nucleon sea. It also comensates for the shortcoming of erturbatively generated arton densities in reroducing the exerimental data on F 2. In order to distinguish these artons from those confined inside the CQ, we will term them as inherent artons. This comonent is intimately related to the bound state roblem, and hence it has a nonerturbative origin. However, for the rocess of CQ CQ+gluon q q, at an initial value of Q 2 =.65GeV 2 where α s is still small enough, we will calculate it erturbatively. The corresonding slitting functions are as follows, P gq (z) = 4 +( z) 2, (3) 3 z P qg (z) = 2 (z2 +(+z) 2 ). (3) For the joint robability distribution of the rocess at hand, we get, q inh. (x,q 2 ) = q inh (x,q 2 ) = N α s dy (2π) x y P qg( x y )g soft(y)dy (32) where q inh. is the inherent arton density. The slitting functions, and the q inh (x,q 2 ) = ( q inh (x,q 2 )), are those of the leading order, rather than NLO. We do not exect that this aroximation will affect the whole structure function, as can be seen from Figure (5). In the above equation N, is a factor deending on Q 2, and G CQ is the constituent quark distribution in the roton as given reviously. The same rocess, can also be a source of SU(2) symmetry breaking of nucleon sea; namely, u sea d sea, and hence, the violation of the Gottfried sum rule (GSR). We will take u this issue in the next subsection. The rocess is deicted in Figure (6). Probability to form a meson-baryon state can be 3

14 written as in Ref. [], P MB (x) = dy y dz F(y,z)R(y,z;x) (33) z where F(y,z) is the joint robability of finding a CQ with momentum fraction y, and an inherent quark or anti-quark of momentum fraction z in the roton. R(y,z;x) is the robability of recombining a CQ of momentum y, with an inherent quark of momentum z, to form a meson of momentum fraction x in the roton. For a more general case, the evaluation of both of these robability functions are discussed in [2]. An earlier and ioneering version was roosed in [3]. In the resent case, these functions are much simler. Guided by the works of Ref. [, 2, 3], we can write, F(y,z) = ΩyGCQ(y)z q inh. (z)( y z) δ (34) R(y,z;x) = ρy a z b δ(y +z ) (35) Here, we take a = b = reflecting that the two CQ s in meson share its momentum almost equally. The exonent δ is fixed for the n and ++ states, using the data from E866 exeriment [4], and the mass ratio of to n. They turn out to be aroximately 8 and 3, resectively. Ω and ρ are the normalization constants, also fixed by data. We recognize that the discussion leading to Eqs. (34, 35) and fixing the arameters are based on a henomenological ansatz. The ansatz and the arameters of Eqs. (34, 35) seems working remarkably well, and further effects that might be of hadronic origin need not be larger than a few ercent. Now it is ossible to evaluate ū M and d M quarks, associated with the formation of meson states[]: d M (x,q 2 ) = x ū M (x,q 2 ) = 2 dy y [P π + n(y)+ 6 P π + (y)] d π ( x y,q2 ) (36) x dy y P π ++(y)ū π( x y,q2 ) (37) where ū π and d π are the valence quark densities in the ion at scale Q 2. The coefficients 2 and 6 are due to isosin consideration. Using Eqs. (7), (8), and (25), we can calculate 4

15 various arton distributions in a ion. Those ertinent to Eqs. (35) and (36) are: ū π val.(x,q 2 ) = G Ū/π (y)ūval./ū(x x y,q2 ) dy y, (38) d π val. (x,q2 ) = G D/π (y)d val./d ( x x y,q2 ) dy y, (39) ū val./ Ū = u val./u. (4) There are some data on the valence structure function of π [5], which can serve as a check for the roosed model. Defining valence structure function of π as F π val. = xūπ val. = xdπ val. (4) Weresenttheresultsofourcalculation forfval. π infigure(7), alongwiththeexerimental data at Q 2 around 6 Gev Asymmetry of the Nucleon Sea and the Gottfried Sum Rule There are exerimental evidences [4],[6] that the Gottfried integral, S G = [F 2 (x) Fn 2 (x)]dx x = dx[ d(x) ū(x)], (42) is less than 3 which is the value exected in the simle quark model. The NMC collaboration [6], result at Q 2 = 4 GeV 2 is, S G =.235±.26, which is significantly smaller than 3. There are several exlanations for this observed violation of the Gottfreid sum rule (GSR). Among them are: flavor asymmetry of the nucleon sea, that is, ū sea d sea,[7] [8], isosin symmetry breaking between roton and neutron, and Pauli blocking, among others. One of these exlanations fits well within our model. It was roosed by Eichten, Inchliffe, and Quigg [9], that valence quark fluctuates into a quark and a ion. This exlanation fits very well in our model. In other words, a nucleon can fluctuate into a meson-nucleon state. The idea is aealing to us and in our model it can be calculated rather easily. According to our model, after a air of inherent q q, is created, a ū can coule to a D-tye CQ to form an intermediate π = Dū, while the 5

16 u quark combines with the other two U-tye CQ s to form a ++. This is the lowest uū fluctuation. Similarly, a d d can fluctuate into the π + n state. Since the ++ state is more massive than the n state, the robability of d d fluctuation will dominate over uū fluctuation, which naturally leads to an excess of d d airs over uū in the roton sea. This rocess is deicted in Figure (6). To summarize, there are three sources that contribute to the sea artons in the roton: Constituent quarks of the roton, the mesonic cloud of the CQ, and the slitting of soft gluon. The combined contribution to the ratio of d ū is as follows: ( d ū ) roton = d M + d inh.+cq ū M +ū inh.+cq. (43) NUSea collaboration at FermiLab, E866 exeriment, has ublished their results [4], for the integral of d ū and d ū at Q = 7.35 GeV. With the rocedure described, we have calculated these values at the same Q and for the same range of x, that is, x = [.2,.35]. The result of the model is, dx( d ū) =.85, (44) to be comared with the exerimental value of.68 ±.6. For the entire range in x, we get: dx( d ū) =.3. The corresonding exerimentally extraolated value is.±.8, which is in excellent agreement with our calculations. At Q = 7.35 GeV, this gives S G =.264. In Figure (8), d(x) ū(x) and d ū of x at Q = 7.35 GeV along with the measured results. in the roton are shown as a function 3.5 Proton Structure Function F 2 We are now in a osition to resent the results for the roton structure function, F 2. In Eqs.(7,8) we resented the form of arton distributions in a CQ. Using those relations, along with the numerical values given in the Aendix, the structure function of the CQ is obtained via Eq. (). The shae of CQ distributions in the roton is rovided by Eqs. 6

17 (23) and (24). Now, all the ingredients are in lace to calculate F 2 from Eq. (2). In Figure (9) the results are shown at many values of Q 2. As it is evident from Figure (4), sole use of CQ structure to reconstruct the F 2 data is not enough. The CQ contribution fall a few ercent short of reresenting the actual data. However, as mentioned earlier, there is an additional contribution from the inherent artons to F 2, which is calculated in Eq. (3). Adding this comonent reresents the data rather well as can be seen from Figure (9). The data oints are from [2]. For the urose of comaring our results with other calculations, we have also included in Figure (9), the GRV s NLO results [2], as well as the rediction of CTEQ4M [22]. Notice that we have taken the number of active flavors to be three for Q 2 5GeV 2 and four, otherwise. In Figure (), the gluon distribution redicted by the model is resented along with those from Ref.[6, 2]. 4 Summary and Conclusion In this aer we have used the notion of constituent quark as a well defined entity being common to all hadrons. Its structure can be calculated erturbatively to all orders in QCD. A CQ receives its own structure by dressing a valence quark with gluon and q q airs in QCD. We have calculated its structure function in the Next-to-Leading order framework. Considering a hadron as the bound states of these constituent quarks, we have used the convolution theorem to extract the hadronic structure functions for roton and ion. Besides the CQ structure contribution to the hadrons, there is also a nonerturbative contribution. This comonent is small, being only a few ercent. But it is crucial in roviding a framework to understand and exlain the violation of the Gottfried sum rule and the excess of d over ū in the nucleon sea. A mechanism is devised for this urose and necessary calculations are erformed. We have resented the results and comared them with all the available and relevant data, as well as with the work of others. We found that our results are in good agreement with the exerimental data. 7

18 Acknowledgment We are grateful to rofessor Rudolh C. Hwa for reading the manuscrit and his helful comments and suggestions. A. N. K. is also grateful to Semnan University artially suorting his research. 5 APPENDIX In this aendix we will give the functional form of arameters of Eqs. (7, 8) in terms of the evolution arameter, t. This will comletely determines artonic structure of CQ and their evolution. The results are valid for three and four flavors, although the flavor number is not exlicitly resent but they have entered in through the calculation of moments. As we exlained in the text, we have taken the number of flavors to be three for Q 2 5GeV 2 and four for higher Q 2 values. I) Valence quark in CQ (Eq. 7): a = t.45t t 3 b =.46.37t +.47t 2.89t 3 c =.3+.37t.23t t 3 II) Sea quark in CQ (Eq. 8): α =.7.23t+.247t 2.8t 3 β = t+.495t 2.455t 3 γ = t t t 3 η = t+.t t 3 ξ = t 8.23t t 3 N = t 2.36t t 3 8

19 III) Gluon in CQ (Eq. 8) α = t+.856t 2.564t 3 β = t+.95t 2.242t 3 γ = t.223t 2.23t 3 η = t+4.66t 2.44t 3 ξ = t 7.458t t 3 N = t+6.879t 2.876t 3 References [] M. Lavelle and D. McMullan, Phys. Lett. B 37, 83 (996). [2] M. Lavelle and D. McMullan, Phys. Re. 279, (997). [3] G. Altarelli, N. Cabibbo, L. Maiani, and R. Petronzio, Nucl. Phys. B 69, 53 (974); N. Cabbibo and R. Petronzio, ibid. B37, 395 (978). [4] R. C. Hwa and M. S. Zahir, Phys. Rev. D 23, 2539 (98); [5] W. Furmanski, R. Petronzio, Z. Phys.C, 293 (982). [6] E. G. Floratos, C. Kounnas, and R. Lacaze, Nucl. Phys. B (98). [7] H. W. Kua, L. C. Kwek, and C. H. Oh, Phys. Rev. D 59, 7425 (999). [8] Rudolh Hwa and C. B. Yang, Phys. Rev. C66, 2524 (22). [9] Rudolh Hwa and C. B. Yang, Phys. Rev. C (22). 9

20 [] J. Magnin and H. R. Christiansen, Phys.Rev. D6, 546 (2). [] H. R. Christiansen and J. Magnin, Phys. Lett. B 445 (998) 8. [2] F. Arash and L. Tomio, Phys. Lett. B 4 (997) 27. [3] K. P. Das and R. C. Hwa, Phys. Lett. B 68 (977) [4] E. A. Hawker Phys. Rev. Lett. 8 (998) 375; J. C. Peng et al., Phys. ReV. D 58 (998) 924. [5] J. S. Comway et al., Phys. Rev. D 39, 92 (989); P. J. Sutton, A. D. Martin, W. J. Stirling, and R. G. Roberts, Phys. Rev D 45 (992) [6] NM Collab., P. Amaudruz et al., Phys. Rev. Lett.66 (99) 272; M. Arneodo et al., Phys. Rev. D5 (994) R. [7] G. Prearata, P. G. Ratcliffe, and J. Soffer, Phys. Rev. Lett. 66 (99) 687. [8] E. M. Henley and G. A. Miller, Phys. lett. B 25 (99) 453. [9] E. J. Eichten, I. Inchliffe and C. Quigg, Phys. Rev. D45 (992) [2] J. Breitweg, et al., The Euro. Phys. Jour C 7, 69 (999); DESY-98-2 (998); M. Adamus, itet al., ZEUS Coll., Phys.Lett. B 47, 432 (997); M. Derrick, et al., Z. Phys. C 72, 399 (996); ibid, Z. Phys. C 65, 293 (995); C. Adloff, H Coll. et al., Nucl. Phys. B 497, 3 (997); T. Ahmed, H Coll. et al., Nucl. Phys. B 439, 47 (995), ibid, Nucl. Phys. B 47, 55 (993). [2] M. Gluck, E. Reya, and A. Vogt, Z. Phys. C 67, 433 (995). [22] CTEQ4, I. L. Lai, et al., Phys. Rev. D 55 (997) 28. 2

21 6 Figure Cation Figure-. Moments of artons in a CQ at Q 2 = 8.5,25,2 GeV 2. The dashed curve reresents the variation of leading order moments and the solid curve is that of next to leading order. Figure-2. Parton distributions in a CQ at a tyical value of Q 2 = 2 GeV 2 Figure-3. Parton distributions in roton and π at Q 2 = 2 GeV 2 as a function of x. Here only erturbatively generated artons are resented. Figure-4 Soft gluon distribution, xg soft, as a function of x at various Q 2 values. Figure-5. The effect of inherent comonent to a sea arton distribution in roton is resented as a function of x. The dashed-dotted line is the erturbatively generated inconstituent quark contribution. The solid line reresents the sum of the two comonents. Figure-6. Processes resonsible to SU(2) symmetry breaking in the nucleon sea and violation of Gottfried sum rule. Figure-7. Pion valence structure function as a function of x at Q 2 = 5.5 (Gev/c) 2. Solid curve is the result of the model calculations and the data oints are from Ref.[7]. Figure-8. The ratio d ū and the difference d ū as a function of x. The solid line in the model calculation and the dotted line is the rediction of CTEQ4M [22]. Data are from Ref. [8]. Figure-9. Proton structure function F 2 as a function of x calculated using the model and comared with the data from Ref. [2] for different Q 2 values. The dot-dot-dashed line is the rediction of GRV Ref.[2] and the dashed line is that of CTEQ4M Ref.[22]. Figure-. The gluon distribution in roton as a function of x at Q 2 = 2 GeV 2. We have also shown the rediction of GRV (dashed-dotted line) and CTEQ4M (dashed line). The data oints are from H collaboration. 2

22 M valence/ CQ GeV 2 (x 4) 25 GeV 2 (x 2) 2 GeV N M sea/ CQ GeV 2 25 GeV 2 (x 2) 2 GeV 2 (x 4) N M gluon/ CQ.. NLO LO 8.5 GeV 2 25 GeV 2 (x 2) 2 GeV 2 (x 4) N Fig.

23 .8 Q 2 = 2 GeV 2.6 zq gluon/cq zq i /CQ.4 zq val. /CQ.2 zq sea/cq z Fig. 2

24 .8.6 xg (a) Q 2 = 2 GeV 2 Model Ref. [22].2.8 xg π Q 2 = 2 GeV 2 (b) xq i xq i π.6.4 xd v xu v.4 π (xu+xd) v.2 xq s.2 xq s π x x Fig. 3

25 .2 xg soft GeV 2 5 GeV 2 2 GeV x Fig. 4

26 xu Q 2 = 35 GeV 2 xu CQ + inh. xu CQ x Fig. 5

27 U U u d d u d π + n D U d u u u u π ++ D U Fig. 6

28 .5 Q 2 = 5.5 GeV 2 π F val..4 Ref. [7] x.3.2. Model Fig. 7

29 d - u d / u Q= 7.35 GeV Model CTEQ4M[22] E866[8] NA 5, X x X x Fig. 8

30 F Q 2 = <Q 2 >=.25 GeV 2 % % *** % * % % % * %.2... X F % % * *% * % % %* %* % % % Q 2 = <Q 2 >=2.2 GeV X X x x x F * % * " % ' ** % % * % % * % % Q 2 = 3.5 GeV 2 Model GRV[2] CTEQ4M[22] ZEUS [2] H 97 [2] ZEUS 96 [2] H [2] ZEUS [2] F Q 2 = 2 GeV 2 " *' " *" ' ' " * " ' " ' X F " Q 2 = 5 GeV 2 " ' "' " " ' " ' ' " " Q 2 = 8-22 " <Q 2 >= 2 GeV 2 " """ " " " "" "" " " X X x x x Fig.9 F 2.5.5

31 F 2 F ' " ' " " ' " ' " '" Q 2 = 2-25 <Q 2 >=22.5 GeV 2 ' " ' " '' " " '"" ' " " " ' Q "' 2 = 35 GeV 2 " "' " ' " ' " ' " " " Q 2 = <Q 2 >=26 GeV 2.3 " F 2 F " " " ' ' Q 2 = GeV 2 " ' ".. " F ' ' "' " "'" ' " " ' Q 2 = 2 GeV 2 ' Q 2 = 65-7 <Q 2 >=67.5 GeV 2 " ' ".5.5 " x x x F 2 x x x Fig.9(continued)

32 ,,,,, Q 2 = 2 GeV 2 xg(x,q 2 )..., Model CTEQ4M [22] GRV [2] H x Fig.