application/xmlThe Fermi motion contribution to [formula omitted] production at the hadron collidersM.A. Gomshi NobaryB. NikoobakhtQuarkoniaThe [formula omitted]Hard QCDFragmentationPhysics Letters B 639 (2006) 488-493. doi:10.1016/j.physletb.2006.06.047journalPhysics Letters BCopyright © 2006 Elsevier B.V. All rights reserved.Elsevier B.V.0370-2693639517 August 20062006-08-17488-49348849310.1016/j.physletb.2006.06.047http://dx.doi.org/10.1016/j.physletb.2006.06.047doi:10.1016/j.physletb.2006.06.047http://vtw.elsevier.com/data/voc/oa/OpenAccessStatus#Full2014-01-01T00:14:32ZSCOAP3 - Sponsoring Consortium for Open Access Publishing in Particle Physicshttp://vtw.elsevier.com/data/voc/oa/SponsorType#FundingBodyhttp://creativecommons.org/licenses/by/3.0/

Journals

S300.2

PLB23152S0370-2693(06)00776-310.1016/j.physletb.2006.06.047Elsevier B.V.Phenomenology

Fig. 1

Fragmentation function for J/ψ production. While the dashed curve represents the function when the Fermi motion is off, the solid ones show the behavior of this function when the Fermi motion is on. The β values are indicated. The function picks up as β increases. This gives raise to increase in the fragmentation probability.

Fig. 2

The differential cross section for direct fragmentation production of J/ψ and its subsequent decay J/ψ→μ+μ− at the Tevatron Run I energies. While the dashed curve is obtained using (5) or equally (7) with β=0, the solid one is due to (7) with β=0.6 GeV. The result is compared with the CDF Run I data. Other contributions are not included. The scale is chosen to be 2μR.

Fig. 3

The differential cross section for direct fragmentation production of J/ψ at Tevatron Run II. The two curves are obtained using (7) with β=0 and 0.6 GeV respectively. The scale has been set to 2μR.

Fig. 4

The differential cross section for direct fragmentation production of J/ψ at the RHIC and the CERN LHC. The shift due to the Fermi motion is significantly increased in the case of LHC but shows to be less important at the RHIC. In both cases the fragmentation function (7) is employed with β=0 and 0.6 GeV, respectively.

The Fermi motion contribution to J/ψ production at the hadron colliders

M.A.

Gomshi Nobary

⁎

mnobary@razi.ac.ir

Nikoobakht

aDepartment of Physics, Faculty of Science, Razi University, Kermanshah, IranbThe Center for Theoretical Physics and Mathematics, AEOI, Roosbeh Building, P.O. Box 11365-8486 Tehran, Iran⁎Corresponding author.

Editor: N. Glover

Abstract

We investigate the relativistic Fermi motion effect in the case of J/ψ production in various hadron colliders. A light-cone wave function is adopted to represent the J/ψ final state. The change in the confinement parameter which sets a scale for the size of the final state, allows one to see the effect in an explicit manner. While the effect has considerable influence on the fragmentation probabilities and the differential cross sections, the total cross sections essentially are left unchanged. Such a feature is in agreement with the momentum sum rule which the fragmentation functions should satisfy.

Keywords

QuarkoniaThe J/ψHard QCDFragmentation

Introduction

Evaluation of the J/ψ cross section at Tevatron energies has been one of the interesting problems of QCD in theory [1–3] and in experiment [4]. Predictions using the QCD calculations where the production of J/ψ is assumed to occur in color singlet form, fails to agree with the experimental results. To bring about the agreement the mechanism of color octet is introduced in which the bound state is produced originally in the color octet state at the production point [5]. Then by emitting a soft gluon the colored object is transformed into a color singlet state [6].Generally in the formation of quarkonia bound state it is assumed that the constituents are not relativistic and therefore they are let to fly together within the bound state ignoring their respective motion. Here we propose that the constituents within such states specially the charm quark is not that heavy to let one to ignore its relativistic effects particularly in relation to the energies at which such a particle is produced. Therefore the effect should be accounted for in the production process in a more explicit and accurate form. The results of such a study may shed light on the problem of J/ψ production cross section.In this work, we introduce the Fermi motion into the J/ψ production in direct fragmentation process using a light-cone wave function. Since our aim is to show the size of the effect, we have not included other contributions. We demonstrate the enhancement of the fragmentation function due to this effect. With such a significant enhancement, we evaluate the differential cross section times the branching ratio and the total integrated cross section for the process p¯p→c→J/ψ→μ+μ− at the Tevatron Run I energies and compare them with the CDF data. We also present similar results for the case of the RHIC, the Tevatron Run II, and the RHIC and the CERN LHC pp collisions.2

The light-cone wave function and the bound state formation via fragmentation

There have been different approaches to introduce the Fermi motion into the production and decay processes of various meson states and quarkonia using different forms of wave functions [7]. The light-cone wave functions have been employed in the case of charmonia to study their production in B-decay [8] and in photoproduction [9]. Here motivated by harmonic oscillator model, we have picked up a wave function in the light-cone quantization to represent the J/ψ bound state. It has the following form [7](1)ψc→J/ψ(x1,x2,qT)=Ac→J/ψexp[−18β2m2+qT2x1x2], where Ac→J/ψ is the normalization coefficient, m is the quark (anti-quark) mass and qT is the transverse momentum of the constituents. The x's are the energy momentum ratios and finally the parameter β is known as the confinement parameter which controls the width of the wave packet representing the bound state. The normalization condition is(2)∑n,λi∫[dx][d2qT]|ψn(xi,qTi,λi)|2=1, where(3)[dx]≡∏i=1ndxiδ[1−∑i=1nxi], and(4)[d2qT]≡∏i=1nd2qTi16π3δ2[∑i=1nqTi]. The sum is over all Fock states and helicities.With the choice of the wave function (1), we fix all degrees of freedom (transverse and longitudinal) of the constituents within the J/ψ bound state by matching them with those in the matrix elements relevant to the fragmentation function and vary the parameter β to change the size of the wave packet representing the bound state. Benefiting such a method, in the leading order perturbative regime, the fragmentation functions for J/ψ production without and with the Fermi motion are obtained as follows [10].2.1

Fermi motion off

In this case the transverse momenta of the constituents are set equal to zero and that the longitudinal components are chosen to be equal. The confinement parameter is β=0 in this case. The fragmentation function is obtained as(5)Dc→J/ψ(z,μ○,β=0)=αs2CF2〈kT2〉1/216mF(z){z(1−z)2[ξ2z4+2ξz2(4−4z+5z2)+(16−32z+24z2−8z3+9z4)]}, where αs is the strong interaction coupling constant and CF is the color factor. The quantity 〈kT2〉 is the average transverse momentum squared of the initial state heavy quark, the parameter ξ is defined as ξ=〈kT2〉/m2 and finally F(x) is given by(6)F(z)=[ξ2z4−(z−2)2(3z−4)+ξz2(8−7z+z2)]2.2.2

Fermi motion on

Here the different momentum components of the constituents, i.e. the quantities q=|qT| and x=x1=1−x2 in the matrix elements and the wave function (1) are integrated over. The only remaining parameter is the confinement parameter β. The fragmentation function in this case is obtained as(7)Dc→J/ψ(z,μ○,β)=π2αs2CF2〈kT2〉1/22m∫dqdx|ψM|2x2(1−z)2zqG(z){1−4(1−x)z+2(4−10x+7x2)z2+4(−1+x3−5x2+4x)z3+(1−4x+8x2−4x3+x4)z4+ηξz2[1−2x+z2+x2(2−2z+z2)]+η[2+(−6+4x)z+(9−8x+2x2)z2−2(2−x+x2)z3+(1+x2)z4]+ξz2[1+2x3(2−3z)z+z2+2x4z2+2x(−1+z−2z2)+x2(2−8z+9z2)]+η2(1−z)2+ξ2(1−x)2x2z4}. The function G(z) reads(8)G(z)={[η(1−z)2+ξx2z2+(1−(1−x)z)2][η(−1+z)+ξ(−1+x)xz2−1+(1−x+x2)z]}2. Here we have defined η=q2/m2. In our study we have set the charm quark mass equal to 1.25 GeV.The fragmentation functions (5) and (7) provide a comparison of the cases of the Fermi motion off and on. The two functions coincide at sufficiently low β. As β increases, (7) raises considerably at the peak region which leads to the respective increase in fragmentation probability. To choose the β value, we note that different values within the range of β=0.250–0.750 GeV have been employed for different meson states [7]. Even the case of β=mq is adopted by Tao Huang et al. in [7]. On the other hand in [10], we have argued that the size of the wave packet representing the bound state is important here. Therefore we have raised the value of β up to 0.6 GeV which we believe is quite safe. The behavior of the fragmentation function (7) is shown in Fig. 1

for β=0, 0.2, 0.4 and 0.6 GeV. We take the maximum value of β=0.6 GeV and use it for our further considerations.

Inclusive production cross section

We have employed the idea of factorization to evaluate the J/ψ production cross section at hadron colliders. For p¯p collisions we may write(9)dσdpT(p¯p→c→J/ψ(pT)X)=∑i,j∫dx1dx2dzfi/p¯(x1,μ)fj/p(x2,μ)[σˆ(ij→c(pT/z)X,μ)Dc→J/ψ(z,μ,β)].Where fi,j are parton distribution functions with momentum fractions of x1 and x2 (different from x1 and x2 which appear in (1)), σˆ is the charm quark production cross section and D(z,μ,β) represents the fragmentation of the produced heavy quark into c¯c state with confinement parameter β at the scale μ. We have used the parameterization due to Martin–Roberts–Stiriling (MRS) [11] for parton distribution functions and have included the heavy quark production cross section up to the order of αs3[12]. The dependence on μ is estimated by choosing the transverse mass of the heavy quark as our central choice of scale defined by(10)μR=pT2(parton)+mc2, and vary it appropriate to the fragmentation scale of the particle state to be considered. This choice of scale, which is of the order of pT (parton), avoids the large logarithms in the process of the form ln(mQ/μ) or ln(pT/μ). However, we have to sum up the logarithms of order of μR/mQ in the fragmentation functions. But this can be implemented by evolving them by the Altarelli–Parisi equation [13]. The following form of this equation is used here(11)μ∂∂μDc→J/ψ(z,μ,β)=∫z1dyyPQ→Q(z/y,μ)Dc→J/ψ(y,μ,β). Here PQ→Q(x=z/y,μ) is the Altarelli–Parisi splitting function. The boundary condition on the evolution equation (11) is the initial fragmentation function Dc→J/ψ(z,μ,β) at some scale μ=μ○. In principle this function may be calculated perturbatively as a series in αs at this scale.Detection of final state requires kinematical cuts of the transverse momentum, pT, and the rapidity, y. We have imposed the required pTcut and ycut in our simulations for different colliders as required and have used the following definition of rapidity(12)y=12log[E−pLE+pL].4

Results and discussion

We have used a light-cone wave function to introduce the Fermi motion in J/ψ production in direct fragmentation channel and obtained its fragmentation function in leading order perturbative regime. In this function the confinement parameter switches the effect of Fermi motion and within its physically acceptable values, the function demonstrates the effect of Fermi motion in an interesting manner. The motivation of introducing this effect is to see its role in improving the QCD versus experimental results for the J/ψ cross section at the Tevatron energies.In the case of the J/ψ state such a fragmentation function, with a reasonable choice of the confinement parameter, gives raise to a significant increase in the fragmentation probability due to the Fermi motion. Such a behavior is illustrated in Fig. 1.To see the effect in the production rates, we have used the usual procedure of factorization for J/ψ production in hadron colliders. We present the pT distribution of BR(J/ψ→μ+μ−)dσ/dpT for the cases of the Fermi motion off and on along with the CDF data at Run I in the Fig. 2

. The branching ratio BR(J/ψ→μ+μ−)=0.0597 is taken from [14]. The poor agreement with data is due to the fact that here we have only considered the contribution of p¯p→c→J/ψ→μ+μ−. Similar behavior at Tevatron Run II energies is shown in Fig. 3

. We have also extended our study to the cases of the RHIC and the CERN LHC pp colliders. Here we provide the pT distributions of the differential cross sections for c→J/ψ and compare the two cases of the Fermi motion off and on in the Fig. 4

. In all cases we have used β=0.6 GeV for the confinement parameter in the fragmentation functions. Naturally, the results for β in the range of 0–0.6 GeV fall between the above results.

We have also calculated the total integrated cross sections for each case. We found that the total cross sections for with and without Fermi motion essentially remain unchanged within the uncertainties of Mote Carlo simulations. The reason is first due to the momentum sum rule which the fragmentation functions should satisfy. In other words although the modification of fragmentation functions by the Fermi motion redistributes the final states, the integrated cross sections are left unchanged. Alternatively although the Fermi motion increases the fragmentation probability for the state, i.e., introduces a state with overall higher mass, the cross section is lowered by just the same amount when we introduce the effect in calculation of the total integrated cross section.First of all, we note that the effect of Fermi motion is indeed significant and that the implementation of such a study in all production channels of J/ψ may give raise to considerable enhancement of the color singlet differential cross section.It is evident from Figs. 2, 3 and 4 that the effect increases with increasing s. The kinematical cuts play important role apart from s. The large cross section at the RHIC compared with the LHC in the Fig. 4 is an example.Finally, we consider the uncertainties involved in our study. There are two main sources of uncertainties.The first is about the simulation of J/ψ production at hadron colliders such as the uncertainties along with the fragmentation functions and parton distribution functions. These kind of uncertainties are well discussed in the literature. The second source of uncertainty is due to the choice of the confinement parameter in the fragmentation function. Relying on our discussion in Section 2, our choice of β=0.6 GeV seems to be justified. Future determination of this parameter will shed more light on this situation. It is worth mentioning that our choice of charm quark mass, i.e., 1.25 GeV, have put our results in their upper side and that the change of the charm quark mass in its acceptable range does not have significant impact on the Fermi motion effect in J/ψ production.

References

[1]

E.L.

Berger

JonesPhys. Rev. D

1981

1521

Baier

RucklZ. Phys. C

1983

251

[2]

S.D.

Ellis

Einhorn

QuiggPhys. Rev. Lett.

1976

1263

[3]

Cacciari

GrecoPhys. Rev. Lett.

1994

1586

Sridhar

A.D.

Martin

W.J.

StirlingPhys. Lett. B

438

1998

211

PetrellNucl. Phys. B (Proc. Supp.)

2000

533

Cooper

hep-ph/0402219

[4]

Albajar

Phys. Lett. B

256

1991

112

Abe

Phys. Rev. Lett.

1992

3704

Abe

Phys. Rev. Lett.

1993

2537

Granier de CassagnacPHENIX CollaborationJ. Phys. G: Nucl. Part. Phys.

2004

S1341

Acosta

Phys. Rev. D

2005

032001

[5]

G.T.

Bodwin

Braaten

G.P.

LepagePhys. Rev. D

1995

1125

Braaten

M.A.

Doncheski

Fleming

ManganoPhys. Lett. B

333

1994

548

G.T.

Bodwin

Braaten

G.P.

LepagePhys. Rev. D

1992

R1914

Cacciari

GrecoPhys. Rev. Lett.

1994

1586

[6]

V.A.

Khoze

A.D.

Martin

M.G.

Ryskin

W.J.

Stirling

hep-ph/0410020

[7]

G.P.

Lepage

S.J.

BrodskyPhys. Rev. D

1980

2157

X.-h.

Guo

HuangPhys. Rev. D

1991

2931

M.G.

RyskinZ. Phys. C

1993

Szczepaniak

hep-ph/9309284

S.J.

Brodsky

Phys. Rev. D

1994

3134

Huang

B.-Q.

Q.-X.

ShenPhys. Rev. D

1994

1490

J.P.

Wang

hep-ph/0412282

J.P.

Z.G.

hep-ph/0405111

Frankfurt

Koepf

StrikmanPhys. Rev. D

1996

3194

[8]

H.Y.

Cheng

K.C.

YangPhys. Rev. D

2001

074011

H.Y.

Cheng

Y.Y.

Keum

K.C.

YangPhys. Rev. D

2002

094023

Z.Z.

Song

K.T.

ChaoPhys. Lett. B

568

2003

127

Z.Z.

Song

Phys. Rev. D

2004

054009

[9]

Frankfurt

Mc Derott

StrikmanJHEP

0103

2001

045

L. Frankfurt et al., in [7]

Hayashigaki

Suzuki

TanakaPhys. Rev. D

2003

093002

[10]

M.A.

Gomshi Nobary

JvadiEur. Phys. J. C

2005