Photoproduction of doubly heavy baryon at the ILC

In the present paper, we make a detailed study on the doubly heavy baryon photoproduction in the future e+e− International Linear Collider (ILC). The baryons Ξcc, Ξbc, and Ξbb are produced via the channel γγ→ΞQQ′+Q¯′+Q¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \gamma \gamma \to {\Xi}_{QQ\prime }+\overline{Q}^{\prime }+\overline{Q} $$\end{document}, where Q and Q′ stand for heavy c or b quark, respectively. As for the ΞQQ′-baryon production, it shall first generate a (QQ′)[n]-diquark and then form the final baryon via fragmentation, where [n] stands for the color- and spin- configurations for the (QQ′)-diquark states. According to the non-relativistic QCD theory, four diquark configurations shall provide sizable contributions to the baryon production, e.g., [n] equals 3S13¯,1S06,3S16,or1S03¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\left[{}^3{S}_1\right]}_{\overline{3}},\;{\left[{}^1{S}_0\right]}_6,\;{\left[{}^3{S}_1\right]}_6,\;\mathrm{or}\kern0.37em {\left[{}^1{S}_0\right]}_{\overline{3}} $$\end{document}, respectively. We adopt the improved helicity amplitude approach for the hard scattering amplitude to improve the calculation efficiency. Total and differential cross sections of those channels, as well as the theoretical uncertainties, are presented. We show that sizable amounts of baryon events can be generated at the ILC, i.e., about 2.0 × 106 Ξcc, 2.2 × 105 Ξbc, as well as 3.0 × 103 Ξbb events are to be generated in one operation year for S=500\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sqrt{S}=500 $$\end{document} GeV and ℒ ≃ 1036 cm−2s−1.


JHEP12(2014)018
According to NRQCD, the Ξ QQ ′ baryon can be expanded over the Fock states, where v is the relative velocity of the constituent heavy quarks in the baryon rest frame. Usually, it is stated that all the baryons are dominated by the first Fock state |(QQ ′ )q , then the emitted gluon from the heavy quark for (QQ ′ ) in [ 1 S 0 ] 6 state must change the spin of the heavy quark; thus, the probability coefficient c 1 (v) shall dominant over other coefficients, or equivalently, h 6 shall be at least v 2 -suppressed to h3 and can be neglected. Here h3 stands for the probability of transforming the color antitriplet diquark into the baryon and h 6 stands for the probability of transforming the color sextuplet diquark into the baryon.
A different power counting rule over v-expansion has also been suggested in the literature. Ref. [5] suggests that the second Fock state |(QQ ′ )qg can be of the same importance as |(QQ ′ )q [5]. Its main idea lies in that one of the heavy quarks can emit a gluon, which does not need to change the spin of the heavy quark, and this gluon can further split into a light qq pair; the light quarks can also emit gluons, and finally, the baryon components can be formed with a light quark q plus one or more soft gluons. Since the light quark can emit gluons easily, we have c 1 (v) ∼ c 2 (v) ∼ c 3 (v). As a rough order estimation, we take the transition probabilities for those diquark states to form the corresponding baryon to be the same, i.e., h 6 ≃ h3. Following this approximation, we shall find that the color sextuplet diquark component can also provide sizable contributions to the baryon production. It is found that those matrix elements are overall parameters, and their uncertainties can be conveniently discussed when we know their values well. For convenience, we will adopt the assumption h 6 ≃ h3 to do our discussions throughout the paper.
The remaining parts of the paper are organized as follows. In section 2, we present the formulation for dealing with the photoproduction channel γγ → Ξ QQ ′ +Q ′ +Q at the leading-order level. In section 3, we give the numerical results. Section 4 is reserved for a summary.

Calculation technology
At the leading order O(α 2 α 2 s ), within the NRQCD factorization approach, the differential cross section for the channel γγ → Ξ QQ ′ +Q ′ +Q can be formulated as where O H (n) stands for the long-distance matrix element, which is proportional to the inclusive transition probability of the perturbative state, (QQ ′ )[n] pair into the heavy baryon Ξ QQ ′ . M is the hard scattering amplitude, which is calculable since the intermediate gluon should be hard enough to generate a heavy QQ or Q ′Q′ pair. means we need to average over the spin states of the electron and positron and sum over the color and spin of all final particles. dΦ 3 is the conventional three-body phase space. f γ (x) is the density function of the incident photon [25,26]. Figure 1. Typical Feynman diagrams for the production channel γ(k 1 ) + γ(k 2 ) → Ξ QQ ′ (p 3 ) + Q ′ (p 4 ) +Q(p 5 ) at the tree level. Other ten diagrams can be obtained by exchanging the positions of the initial photons attached to the quark lines.
The hard scattering amplitude M = 20 k=1 M k for the process can be written in a general form as where the subindices m and n are color indices of the constituent heavy quarks, and l is the color of the diquark. a = 1, . . . , 8 is the color index of the gluon propagator.
the normalization factor. The function G mnl equals the antisymmetric ε mnl (the symmetric f mnl ) for the color antitriplet3 (the color sextuplet 6) of (QQ ′ ) diquark. The sum of the anti-symmetric and the symmetric functions satisfy the following equations and With the help of the above relations, we obtain C 2 ijl = 4 3 for the color-antitriplet diquark state and C 2 ijl = 2 3 for the color-sextuplet diquark state, respectively. All the amplitudes M k with k = (1, . . . , 20) contain massive quark lines, so it is too complicated and lengthy by using the conventional trace technique to deal with the amplitude square. To shorten the calculations and to make the results more compact, we adopt the improved helicity amplitude approach [22] to deal with the difficulty of calculating the expressions for the yields when the quark masses cannot be neglected. It is found that we can connect the doubly heavy baryon production with those of doubly heavy quarkonium production. We have made a detailed discussion on the heavy quarkonium production at the ILC under the improved helicity amplitude approach via the channel γγ → |[QQ ′ ](n) + Q ′ +Q in ref. [46]. To compare with the quarkonium case, by applying the charge conjugation matrix C = −iγ 2 γ 0 and the transverse of the matrix element to the amplitude M k , we can transform eq. (2.3) as where ρ stands for the number of the γ-matrixes appearing in the amplitude M k . The second line is the matrix element for the heavy quarkonium production, which indicates that the amplitudes for the diquark production are merely different from those of the heavy quarkonium case with an overall factor (−1) ρ+1 . Thus, inversely, we can conveniently derive the hard scattering amplitudes M k from ref. [46] after proper transformation. To shorten the paper, we will not put the detailed calculation technology for the baryon production here, the interesting readers may turn to ref. [46] for details of the improved helicity amplitude approach. Here, to derive eq. (2.7), we have implicitly applied the relations:

Numerical results and discussions
As discussed in the Introduction, we adopt h 6 ≃ h3 to do our discussion. The nonperturbative matrix element with color antitriplet diquark, h3, can be related to the Schrödinger wave functions at the origin |ψ (QQ ′ ) (0)| as [5]: Since the spin-splitting effect is small, we do not distinguish the bound state parameters for the spin-singlet and the spin-triplet states; i.e., those parameters, such as the constituent quark masses, the bound state mass, and the wave function, are taken to be the same for the spin-singlet and spin-triplet states. We take the wavefunctions at the origin as [3]: |Ψ cc (0)| 2 = 0.039 GeV 3 , |Ψ bc (0)| 2 = 0.065 GeV 3 , and |Ψ bb (0)| 2 = 0.152 GeV 3 . The heavy quark masses are taken as: m c = 1.5 GeV and m b = 4.9 GeV. The doubly heavy baryon mass is taken as M Ξ QQ ′ = m Q + m Q ′ . The other parameters are taken as the same as those of ref. [46], e.g., the renormalization scale is taken as the transverse mass, µ r = M t = M 2 QQ ′ + p 2 t . As a cross check of our calculation, we obtain same numerical results as those derived from the conventional squared amplitude approach.
Total cross sections for the doubly heavy baryon photoproduction with three collision energies, i.e., √ S = 250 GeV, 500 GeV, and 1 TeV, are put in table 1. Summing up the contributions from different color-and spin-configurations, we find that the total cross sections decrease with the increment of √ S, i.e., σ Ξcc | 250GeV : σ Ξcc | 500GeV : σ Ξcc | 1TeV ≃ 6 : 3 : 1, It is noted that the relative importance among different color-and spin-configurations for the total cross sections and the differential distributions are similar under different collision energies. In the following, we take √ S = 500 GeV as the e + e− collision energy. Under the condition of √ S = 500 GeV, we obtain It indicates that the [ 3 S 1 ]3 diquark state provides the dominant contribution, while other configurations may also provide significant contributions. By summing up all the possible diquark configurations, we obtain σ Ξcc = 202.26 fb, σ Ξ bc = 22.47 fb, and σ Ξ bb = 0.3 fb. If the integrated luminosity is as high as 10 4 fb −1 , we shall have about 2.0 × 10 6 Ξ cc , 2.2 × 10 5 Ξ bc , and 3.0 × 10 3 Ξ bb events to be generated through the direct photon collision at the ILC in an operation year. The Ξ cc production rate is larger than those of Ξ bc and Ξ bb , i.e., σ Ξcc : σ Ξ bc : σ Ξ bb = 647 : 75 : 1. Thus, in the following, we shall focus on the photoproduction of Ξ cc and Ξ bc . Figure 2 shows the baryon transverse momentum (p t ) distributions for the photoproduction of Ξ cc and Ξ bc . Similar to the above conclusion, the [ 3 S 1 ]3 configuration for both Ξ cc and Ξ bc production provides dominant contributions over the other configurations in the whole p t region. We present the rapidity (y) and pseudorapidity (y p ) distributions in figures 3 and 4. There is a plateau within |y| < 4 or |y p | < 4. We present the differential cross sections dσ/dz in figure 5, where z = 2 s (k 1 +k 2 )·p 3 withŝ = x 1 x 2 S being the invariant mass of the initial photons of the subprocess. In the subprocess center-of-mass frame, z is simply twice the fraction of the total energy carried by the baryon and is experimentally observable. To be useful references, we present the total cross sections under various p t or y cuts in tables 2 and 3.
In the literature, people usually takes a simple assumption by treating the evolution from the diquark to doubly heavy baryon with 100% probability and with equal importance for all phase-space point; we call it the "direct evolution". In the present paper, we have adopted the fragmentation approach with the help of the fragmentation function (2.2) to deal with such evolution; we call it the "evolution via fragmentation". In table 4, we present a comparison of the total cross sections for the baryon photoproduction at the ILC under   p t-cut 1 GeV 2 GeV 3 GeV (cc) 6 Table 3. Total cross sections (in units fb) for the photoproduction of Ξ cc and Ξ bc with √ S = 500 GeV under various color-and spin-configurations and various rapidity cuts. those two treatments. The subscript "d" stands for the "direct evolution", the subscript "f" stands for the "evolution via fragmentation". Table 4 shows the discrepancies of total cross sections for those two treatments are quite small, i.e., less than ∼ 1%. On the other hand, the differences for the p t distributions are also very small in the whole p t region. For example, we put a comparison of the p t distributions under those two treatments in figure 6. We also present a comparison of the z distributions under those two treatments in figure 7. In those two figures, we have summed up the contributions from the mentioned color-and spin-diquark-configurations for convenience. As for the Ξ cc production, it is found that the differential cross sections for the "evolution via fragmentation" are slightly larger in small z region, while slightly smaller in large z region. This result is consistent with the small differences for the p t distributions under those two treatments. Thus the conventional treatment is viable and provides a good approximation to deal with the heavy baryon production. As a sound estimation, we take the fragmentation approach to do the discussion.

JHEP12(2014)018
σ d σ f (cc) 6 Table 4. Comparison of the total cross sections (in units fb) for the baryon photoproduction at the ILC with √ S = 500 GeV. The subscript "d" stands for the "direct evolution", the subscript "f" stands for the "evolution via fragmentation".  Figure 6. Comparison of the p t distributions for the photoproduction of the Ξ cc and Ξ bc baryons at the ILC with √ S = 500 GeV. The superscript "d" stands for the "direct evolution", the superscript "f" stands for the "evolution via fragmentation".  Figure 7. Comparison of the z distributions for the photoproduction of the Ξ cc and Ξ bc baryons at the ILC with √ S = 500 GeV. The superscript "d" stands for the "direct evolution", the superscript "f" stands for the "evolution via fragmentation".  Table 6. Uncertainties for the total cross sections (in units fb) by taking m b = 4.9 ± 0.2 GeV. m c = 1.5 GeV and µ r = M t .
As a final remark, we make a discussion on the theoretical uncertainties from the heavy quark masses. For the purpose, we set m c = 1.50 ± 0.10 GeV and m b = 4.9 ± 0.20 GeV. As shown in table 5, the uncertainties for m c = 1.50 ± 0.10 GeV are σ (cc) 6

Summary
We have investigated the photoproduction of the doubly heavy baryons at the ILC within NRQCD. The improved helicity amplitude approach has been adopted to improve the calculation efficiency. By taking the assumption, h 6 ≃ h3, we observe that the channel via the intermediate [ 3 S 1 ]3 diquark state provides the dominant contribution, while other configurations may also provide significant contributions. Total and differential cross sections, together with their theoretical uncertainties, have been presented. By taking the errors from the heavy quark masses into consideration, we shall have 2.0 +0.68 −0.47 × 10 6 Ξ cc and 2.2 +0. 37 −0.29 × 10 5 Ξ bc events to be produced in one operation year at the ILC with √ S = 500 GeV and L ≃ 10 36 cm −2 s −1 . Thus, the ILC would provide another good platform for studying Ξ QQ ′ -baryon properties.