it is easy to find products of the colour factors Ci,j = â color. CiCâ ... Î¨(0). 2âmcc. , amplitudes Mi, with corresponding colour factors Ci, ...

0 downloads 0 Views 125KB Size

arXiv:hep-ph/0104173v1 18 Apr 2001

Abstract The production of baryons containing two charmed quarks (Ξ∗cc or Ξcc ) in hadronic interactions at high energies and large transverse momenta is considered. It is supposed, that Ξcc -baryon is formed during a non-perturbative fragmentation of the (cc)diquark, which was produced in the hard process of c-quark scattering from the colliding protons: c + c → (cc) + g. It is shown that such mechanism enhances the expected doubly charmed baryon production cross section on Tevatron and LHC colliders approximately 2 times in contrast to predictions, obtained in the model of gluon - gluon production of (cc)-diquarks in the leading order of perturbative QCD.

1

Introduction

Doubly heavy baryons take the special place among baryons which contain heavy quarks. The existence of two heavy quarks causes the brightly expressed quark - diquark structure of the Ξcc baryon, in a wave function which one’s the configuration with compact heavy (QQ)diquark dominates. Regularity in a spectrum of mass of doubly heavy baryons appear in many respects to a similar case of mesons containing one heavy quark [1, 2, 3, 4]. Production mechanisms for (QQq)-baryons and (Q¯ q )-mesons also have common features. At the first stage compact heavy (QQ)-diquark is formed, than it fragments in a final (QQq)-baryon, picking up a light quark. The calculations of production cross sections for doubly heavy baryons in ep- and pp-interactions were made recently as in the model of a hard fragmentation of a heavy quark in doubly heavy diquark [5, 6, 7] as within the framework of the model of precise calculation of cross section of a gluon – gluon fusion into doubly heavy diquarks and two heavy antiquarks in the leading order of the perturbation theory of QCD [8, 9, 10]. Mechanism of production of hadrons containing charmed quarks, based on consideration of hard parton subprocesses with one c-quark in an initial state, was discussed earlier in papers [11, 12, 13]. It was shown, that in the region of a large transferred momentum (Q2 >> m2c , where mc is charmed quark mass) the concept of a charm excitation in a hadron does not contradict parton model and allows to effectively take into account the contribution of the high orders of the perturbative QCD theory to the Born approximation. However, there is open problem of the ”double score”, which is determined by the fact that the part of the Born diagrams of birth of two heavy quarks in a gluon - gluon fusion can be interpreted, as the diagrams with charm excitation in one of initial protons. These diagrams give leading in αs contribution in the c-quark perturbative, so-called point-like structure function (SF) of a proton. As to the non-perturbative contribution in c−quark SF of a proton [14] it does not depend from Q2 and becomes very small at Q2 >> m2c . For example, fig. 1 shows one of 36 Born diagrams, which have order αs4 , describing production of the (cc)-diquark in the gluon-gluon fusion subprocess. The experience in 1

[email protected]

1

calculation of heavy quark production cross sections in a gluon - gluon fusion demonstrates that the contribution of the next order of the perturbative QCD in αs can be comparable with the contribution of the Born diagrams. In the case of gluon – gluon production of the two pairs of heavy quarks there will be more than three hundred diagrams with additional gluon in the final state, which have order αs5 , and their direct calculation is now considered difficultly feasible.

2

Subprocess c + c → (cc) + g

In this paper the model of (cc)-diquark production in proton - proton interactions, based on the mechanism of the charm excitation in a proton is considered. It is supposed, that the (cc)-diquark is formed during scattering of c-quarks from colliding protons with radiation of a hard gluon, i.e. in the parton subprocess: c + c → (cc) + g

(1)

The Feynman diagrams of the parton subprocess (1) are shown in fig. 2, where q1 and q2 are 4-momenta of the initial c-quarks, k is 4-momentum of the final gluon, p is 4-momentum of the diquark, which one is divided equally between the final c-quarks. The doubly heavy diquark is considered as bound state of two c-quarks in the antitriplet colour state and in the vector spin state. If i and j are colour indexes of initial quarks, and m is colour index of a final diquark, the amplitude of production of the (cc)-diquark Mijm (c + c → (cc) + g) is p connected with the amplitude of production of two c-quarks with 4-momenta p1 = p2 = 2 as follows: εnmk p Mijm (c + c → (cc) + g, p) = K0 √ Mijnk (c + c → c + c + g, p1 = p2 = ), 2 2

(2)

s

2 Ψcc (0), mcc = 2mc is the diquark mass, Ψcc (0) is the diquark wave function mcc εnmk in zero point, √ is the colour part of a diquark wave function. Considering spin degrees 2 of freedom of c−quarks and (cc)-diquark, we have following conformity between amplitudes of birth of free quarks and a diquark with fixed spin projections (without colour indexes and common factor K0 ): where K0 =

M(c + c → (cc) + g, sz = +1) ∼ M(c + c → c + c + g, s1z = +1/2, s2z = +1/2) (3) 1 1 M(c + c → (cc) + g, sz = −1) ∼ M(c + c → c + c + g, s1z = − , s2z = − ) (4) 2 2 " 1 1 1 M(c + c → (cc) + g, sz = 0) ∼ √ M(c + c → c + c + g, s1 z = + , s2z = − ) + 2 2 2 # 1 1 (5) + M(c + c → c + c + g, s1z = − , s2z = + ) . 2 2

2

Because the wave function of the (cc)-diquark is antisymmetric on colour index and symmetric on remaining indexes, the production of the scalar (cc)-diquark is forbidden, i.e. 1 1 1 1 M(c+c → c+c+g, s1z = + , s2z = − )−M(c+c → c+c+g, s1z = − , s2z = + ) = 0. (6) 2 2 2 2 Amplitudes adequate to the diagrams in fig. 2, where the final c-quarks are in the arbitrary spin states, are written out below, without the colour factors and the common factor K0 : ¯ 1 )γ µ (ˆ M1 = gs3εµ (k)U(p p1 + kˆ + mc )γ ν U(q1 )U¯ (p2 )γν U(q2 )/((p1 + k)2 − m2c )(p2 − q2 )2(7) ¯ 1 )γ ν U(q1 )U¯ (p2 )γ µ (ˆ M2 = g 3εµ (k)U(p p2 + kˆ + mc )γν U(q2 )/((p2 + k)2 − m2 )(q1 − p1 )2(8) M3 =

M4 = M5 =

s c 3 ν µ 2 2 ˆ ¯ 1 )γ (ˆ gs εµ (k)U(p q1 − k + mc )γ U(q1 )U¯ (p2 )γν U(q2 )/((q1 − k) − mc )(p2 − q2 )2(9) 2 ¯ 1 )γ ν U(q1 )U¯ (p2 )γ ν (ˆ gs3εµ (k)U(p q2 − kˆ + mc )γµ U(q2 )/((q2 − k)2 − m2c )(q1 − p1 )(10) ¯ 1 )γ ν U(q1 )U¯ (p2 )γ λ U(q2 )Gλµν (p2 − q2 , k, p1 − q1 )/(q1 − p1 )2 (p2 − q2(11) gs3εµ (k)U(p )2

√ where gs = 4παs , αs is strong coupling constant, Gλµν (p, k, q) = (p − k)ν gλµ + (k − q)λgνµ + (q − p)µ gνλ . Let’s remark, that the amplitudes M6 − M10 are received by replacement of the initial quarks momenta q1 ↔ q2 in the amplitudes M1 − M5 ans are taken with a minus sign, that allows for the antisymmetrization of the initial state of two identical c-quarks. The corresponding colour factors are presented by the following expressions: εnmk c b εnmk c b b C1 = √ (Tnl Tli )(Tkj ), C6 = √ (Tnl Tlj )(Tkib ), 2 2 εnmk c εnmk b )(Tklc Tljb ), C7 = √ (Tnj )(Tklc Tlib ), C2 = √ (Tni 2 2 nmk nmk ε ε b c b b c C3 = √ (Tnl Tli )(Tkj ), C8 = √ (Tnl Tlj )(Tkib ), 2 2 nmk nmk ε ε b b )(Tklb Tljc ), C8 = √ (Tnj )(Tklb Tlic ), C4 = √ (Tni 2 2 iεnmk b iεnmk b a )(Tkj )f bac , C10 = √ (Tnj )(Tkia )f bac . C5 = √ (Tni 2 2

(12)

Using known property of a completely antisymmetric tensor of the third rank εn′mk′ εnmk = δnn′ δkk′ − δnk′ δkn′ , it is easy to find products of the colour factors Ci,j =

X

Ci Cj∗ , which ones are presented in

color

the Appendix A. The method of the calculation of a production amplitude of a bound nonrelativistic state of quarks in the fixed spin state is based on a formalism of the projection operator [15]. Using properties of the charge conjugation matrix C = iγ2 γ0 , we can link a scattering amplitude of a quark on a quark with a scattering amplitude of an antiquark on a quark, for example: M1 = gs3 εµ (k)U¯ (p1 )γ µ (ˆ p1 + kˆ + mc )γ ν U(q1 )U¯ (p2 )γν U(q2 )/((p1 + k)2 − m2c )(p2 − q2 )2 = = gs3 εµ (k)V¯ (q1 )γ ν (−ˆ p1 − kˆ + mc )γ µ V (p1 )U¯ (p2 )γν U(q2 )/((p1 + k)2 − m2c )(p2 − q2 )2 . (13) 3

As it may be shown, at p1 = p2 =

p one has: 2

1 1 V (p1 , s1z = − )U¯ (p2 , s2z = + ) ∼ εˆ(p, sz = +1)(ˆ p + mcc ), 2 2 1 1 p + mcc ), (14) V (p1 , s1z = + )U¯ (p2 , s2z = − ) ∼ εˆ(p, sz = −1)(ˆ 2 2 1 1 h 1 ¯ √ V (p1 , s1z = − )U(p 2 , s2z = + ) + 2 2 2 i 1 1 +V (p1 , s1z = + )U¯ (p2 , s2z = − ) ∼ εˆ(p, sz = 0)(ˆ p + mcc ), 2 2 where εµ (p) – is polarization 4-vector of a spin-1 particle. After following effective replacements V (p1 )U¯ (p2 ) → εˆ(p)(ˆ p + mcc ) and K0 → K, (15) Ψ(0) where p1 = p2 = p/2 and K = √ , amplitudes Mi , with corresponding colour factors Ci , 2 mcc describe production of the (cc)-diquark with fixed polarization. The square of the module of amplitude of (cc)-diquark production after average on spin and colour degrees of freedom is given by the following expression: |M|2 =

10 X 1 2 X Mi (p1 , p2 )Mj∗ (p1 , p2 )Ci,j , K 36 i,j=1 spin

(16)

p . The summation on vector diquark 2 polarizations in the square of the amplitude of the process (1) was done using the standard formula: X pµ pν µ ∗ν µν ε (p)ε (p) = −g + 2 . (17) mcc spin where in the amplitudes Mi we have put p1 = p2 =

The calculation of the value F =

10 X X

Mi Mj∗ Ci,j have been executed using the package

i,j spin

of an analytical calculations FeynCalc [16]. The answer is shown in the Appendix B, as a function of standard Mandelstam variables sˆ and tˆ.

3

RESULTS OF CALCULATIONS

In the parton model the cross section of a (cc)-diquark production in pp-interactions is represented as follows:

where

Z ymax Z 1 dσ (pp → (cc) + X) = p⊥ dy dx1 Cp (x1 , Q2 )Cp (x2 , Q2 ) × dp⊥ ymin x1min |M|2 1 √ × , × 16π(s(s − m2cc ))1/2 x1 s − sm⊥ ey

√ x1 sm⊥ e−y − 23 m2cc √ , x2 = x1 s − sm⊥ ey 4

(18)

x1min =

√

sm⊥ ey − 23 m2cc √ , s − sm⊥ e−y

Cp (x, Q2 ) is the c-quark distribution function in a proton at Q2 = m2⊥ = m2cc + p2⊥ , p⊥ is the diquark transverse momentum, y is the rapidity of the diquark in c.m.f. of colliding protons, m2cc , 2 √ 3 tˆ = (q1 − p)2 = m2cc − x1 sm⊥ e−y , 2 √ 3 uˆ = (q2 − p)2 = m2cc − x2 sm⊥ ey . 2 sˆ = (q1 + q2 )2 = x1 x2 s +

(19)

It is supposed that spin- 12 and spin- 23 Ξcc −baryons relative yield is 1 : 2 as it is predicted by the simple counting rule for the spin states. The production cross section of Ξcc -baryons plus Ξ∗cc −baryons in our approach is connected with the production cross section of (cc)diquark within the framework of a model of a non-perturbative fragmentation as follows: dσ (pp → Ξcc + X) = dp⊥

Z

0

1

p⊥ dz dσ (pp → (cc)X, p′⊥ = )D(cc)→Ξcc (z, Q2 ), ′ z dp⊥ z

(20)

where D(cc)→Ξcc (z, Q2 ) is the phenomenological function of a fragmentation, normalized approximately on unity, as a total probability of transition (cc)-diquark in final doubly charmed baryon. At Q2 = m2cc the fragmentation function is selected in the standard form [17]: D(cc)→Ξcc (z, Q2 ) =

D0 , m2q 2 m2Ξ 2 − ) z(mcc − z 1−z

(21)

where mΞ = mcc + mq is the Ξcc -baryon mass, mq is the light quark mass, D0 the is ratefixing constant. The fragmentation function for Q2 > Q20 can be determined by the solving the DGLAP evolution equation [18]. Following to paper [10], at numerical calculations we have used following values of parameters: mcc = 3.4 GeV, αs = 0.2, |Ψcc(0)|2 = 0.03 GeV3 , 2 mq = 0.3 GeV. For a c−quark distribution function in a proton Cp (x, √ √ Q ) the parametrization CTEQ5 [19] was used. In figures 3 and 4 at s = 1.8 TeV and s = 14 TeV, accordingly, the curves show results of our calculations of p⊥ -spectra (|y| < 1) of Ξcc -baryons, the stars show results of the calculations from paper [10], adequate to the contribution of the gluongluon fusion production of Ξcc -baryons in a Born approximation. Thus, our calculations demonstrate, that the observed production cross section of Ξcc -baryons on colliders Tevatron and LHC can be approximately 2 times more at the expense of the contribution of the parton subprocess c + c → (cc) + g, than it was predicted earlier in the papers [8, 10]. The authors thank S.P. Baranov, V.V. Kiselev and A.K. Likhoded for useful discussions. The work is executed at support of the Program ”Universities of Russia – Basic Researches” (Project 02.01.03).

5

References [1] S.Fleck and J.Richard, Prog.Theor.Phys.82, 760 (1989). [2] E.Bagan et al., Z.Phys. C64, 57 (1994). [3] V.V.Kiselev et al., Phys.Lett. B332, 411 (1994). [4] D.Ebert et al., Z.Phys.C76, 111 (1997). [5] A.F.Falk et al., Phys.Rev. D49, 555 (1994). [6] A.P.Martynenko and V.A.Saleev, Phys.Lett. B385, 297 (1996). [7] V.A. Saleev, Phys.Lett., B426, 384 (1998). [8] A.V.Berezhnoy et al., Yad.Fiz. 59, 909 (1996). [9] S.P.Baranov, Phys.Rev. D54, 3228 (1996). [10] A.V.Berezhnoy et al., Phys.Rev. D57, 4385 (1998). [11] V.A.Saleev, Mod.Phys.Lett. A9, 1083 (1994). [12] A.P.Martynenko and V.A.Saleev, Phys.Lett. B343, 381 (1995). [13] S.P.Baranov, Phys.Rev. D56, 3046 (1997). [14] S.J.Brodsky and R.Vogt, Nucl.Phys. 478, 311 (1996). [15] B. Guberina et al., Nucl. Phys., B174, 317 (1980). [16] R.Mertig, The FeynCalc Book, Mertig Research and Consulting, 1999. [17] C. Peterson, Phys.Rev. D27, 105 (1983). [18] V.N. Gribov and L.N. Lipatov, Sov.J.Nucl.Phys. 15, 438 (1972); Yu.A. Dokshitser, Sov.Phys.JETP, 46, 641 (1977); G. Altarelli and G.Parisi, Nucl.Phys., B126, 298 (1977). [19] H.L. Lai et al., (CTEQ Coll.), Preprint hep-ph 9903282.

6

Appendix A C1,1 = 97 , C1,2 = 19 , C1,3 = − 92 , C1,4 = 10 , 9 C1,5 = −1, C1,6 = − 19 , C1,7 = − 97 , C1,8 = − 10 , 9 2 C1,9 = 9 , C1,10 = −1, C2,2 = 79 ,

C2,3 = 10 , 9 2 C2,4 = − 9 , C2,5 = 1, C2,6 = − 97 , C2,7 = − 91 , C2,8 = 29 , C2,9 = − 10 , 9 C2,10 = 1, C3,3 = 16 , 9 8 C3,4 = − 9 , C3,5 = 2,

C3,6 = − 10 , 9 2 C3,7 = 9 , C3,8 = 89 , C3,9 = − 16 , 9 C3,10 = 2, C4,4 = 16 , 9 C4,5 = −2, C4,6 = 92 , C4,7 = − 10 , 9 16 C4,8 = − 9 , C4,9 = 98 ,

7

C4,10 = −2, C5,5 = 3, C5,6 = −1, C5,7 = 1, C5,8 = 2, C5,9 = −2, C5,10 = 3, C6,6 = 97 , C6,7 = 91 , C6,8 = − 92 , C6,9 = 10 , 9

C6,10 = 1 C7,7 = 97 C7,8 = 10 9 C7,9 = − 29 C7,10 = 1 C8,8 = 16 9 C8,9 = − 89 C8,10 = 2 C9,9 = 16 9 C9,10 = −2 C10,10 = 3

Appendix B F = −(4παs )3

512FN 9FD

FN = 26361 M 18 − 6 M 16 20513 s + 67472 t +

(22)

+16 M 14 14621 s2 + 100076 s t − 86020 t2 −

−16 M 12 14873 s3 + 122408 s2 t − 657280 s t2 − 382560 t3 +

+64 M 10 2101 s4 − 658 s3 t − 509652 s2 t2 − 468736 s t3 − 170408 t4 + 2

+65536 s t2 s + t

9 s4 + 11 s3 t + 13 s2 t2 + 4 s t3 + 2 t4

−256 M 8 120 s5 − 8749 s4 t − 201737 s3 t2 −

−

−255896 s2 t3 − 149332 s t4 − 44640 t5 −

−1024 M 6 7 s6 + 2180 s5 t + 44390 s4 t2 +

+74060 s3 t3 + 57876 s2 t4 + 28176 s t5 + 7184 t6

−16384 M 2 t 10 s7 + 353 s6 t + 924 s5 t2 +

−

+1151 s4 t3 + 898 s3 t4 + 460 s2 t5 + 160 s t6 + 28 t7 +

+4096 M 4 s7 + 235 s6 t + 5484 s5 t2 + 11610 s4 t3 + +11609 s3 t4 + 7368 s2 t5 + 3056 s t6 + 672 t7

FD = M 2 − s

2

4

M2 − 4 t

8

(23)

4

5 M2 − 4 s + t

(24)

Figure 1: One of the Born diagrams used for description subprocess g + g → (cc) + c¯ + c¯.

q1

k p

q2

+

( q1

!q

2)

Figure 2: Diagrams used for description subprocess c + c → (cc) + g.

9

Figure √ 3: Cross section of Ξcc -baryon production at s = 1.8 TeV and |y| < 1. Stars (*) show the results of calculation from paper [10], curve is our result obtained in the model of a charm excitation in colliding protons.

10

Figure √ 4: Cross section of Ξcc -baryon production at s = 14 TeV and |y| < 1. Stars (*) show the results of calculation from paper [10], curve is our result obtained in the model of a charm excitation in colliding protons.

11