Production of doubly charmed baryons nearly at rest

TL;DR

This paper studies the production characteristics of doubly charmed baryons, focusing on their near-rest production via the intrinsic heavy quark mechanism, with implications for fixed-target experiments.

Contribution

It introduces the production of doubly charmed baryons near rest via the intrinsic heavy quark mechanism, providing predictions for experimental measurements.

Findings

Predicted production cross sections and distributions for doubly charmed baryons.

Near-rest production makes these baryons accessible in fixed-target experiments.

Potential observability at STAR@RHIC and AFTER@LHC.

Abstract

We investigate the production cross sections, momentum distributions and rapidity distributions for doubly charmed baryons which according to the intrinsic heavy quark mechanism are produced nearly at rest. These events should be measurable at fixed-target experiments like STAR@RHIC and AFTER@LHC.

Equations

$$R=\frac{\sigma ({\Xi }_{cc}^{++})\cdot Br({\Xi }_{cc}^{++}\to {\Lambda }_c^{+}{K}^{-}{\pi }^{+}{\pi }^{+})}{\sigma ({\Lambda }_c^{+})}=\frac{N_{{\Xi }_{cc}^{++}}}{{ϵ}_{++}}\cdot \frac{{ϵ}_{{\Lambda }_c^{+}}}{N_{{\Lambda }_c^{+}}},$$

$$R\sim (2-3)\times 10^{-3}.$$

$$\displaystyle\frac{dP_{iQ_{1}Q_{2}}}{\prod_{i=1}^{n}dx_{i}}\propto\alpha_{s}^{4}(M_{Q_{1}\bar{Q}_{1}})\alpha_{s}^{4}(M_{Q_{2}\bar{Q}_{2}})\frac{\delta\big{(}1-\sum_{i=1}^{n}x_{i}\big{)}}{\big{(}\sum_{i=1}^{n}\hat{m}_{T,i}^{2}/x_{i}\big{)}^{2}},$$

$${\sigma }_{icc}=\frac{P_{icc}}{P_{ic}}\cdot {\sigma }_{ic},{\sigma }_{ic}=P_{ic}\cdot {\sigma }^{inel}\frac{{\mu }^2}{4{\hat{m}}_c^2},$$

$$\sigma (cc)=\frac{1}{3}{f}_{cc}{\sigma }_{icc},$$

$$f^{icc}_{cc}=\int_{4m_{c}^{2}}^{4m_{D}^{2}}dM_{cc}^{2}\frac{dP_{icc}}{dM_{cc}^{2}}\,\,\bigg{/}\int_{4m_{c}^{2}}^{s}dM_{cc}^{2}\frac{dP_{icc}}{dM_{cc}^{2}}.$$

$$\sigma ({\Xi }_{cc})\sim (0.2-0.3)\cdot \sigma (cc).$$

$$\sigma ({\Xi }_{cc})\approx (0.2-0.3)\times 75\text{nb}.$$

$$E_{lab}=\frac{1}{2m_{tar}}(m_{cc}^2+m_{tar}^2),p_{lab}=\frac{1}{2m_{tar}}(m_{cc}^2-m_{tar}^2).$$

$$\sigma ({\Xi }_{cc})\approx (0.2-0.3)\times 65\text{mbox}nb.$$

Full text

1 Introduction

Doubly heavy baryons are a rigorous prediction of quantum chromodynamics. In the 2000’s, evidence for the existence of the doubly charmed baryons was reported by the SELEX Collaboration [1, 2, 3, 4, 5, 6]. However, the discrepancy of the production properties derived from the SELEX data with predictions of perturbative QCD induced wrong expectations for the production rates of the doubly charmed baryons in the different production environments, leading to reports stating the non-evidence of doubly heavy baryons [7, 8, 9, 10].

In Ref. [11] it is shown that the intrinsic charm mechanism can fill the gap between the rates motivated by the SELEX data and the theoretically predicted production rates. Moreover, the ratio $\sigma(c\bar{c}c\bar{c})/\sigma(c\bar{c})$ calculated from the SELEX data is comparable or even smaller than the same rate derived from the double $J/\psi$ production data at the NA3 experiment [12, 13]. In addition, the intrinsic charm mechanism can perfectly explain the kinematics of the doubly charmed baryons with $\langle x_{F}\rangle\sim 0.33$ and the relatively small mean transverse momentum of approx. $1\mbox{\,GeV}/c$.

The most recent Belle and LHCb results cannot argue against the SELEX data. The Belle experiment [9] presented an upper limit on the cross section $\sigma(e^{+}e^{-}\to\Xi_{cc}^{+}+X)$ of $82-500\mbox{\,fb}$ at 95% confidence level for the decay mode with $\Lambda_{c}^{+}$ at $\sqrt{s}=10.58\mbox{\,GeV}$, using a luminosity of $980\mbox{\,fb}^{-1}$. This turns out to be at least twice as much as the theoretical upper limit given by $\sigma(\Xi_{cc}^{+})\approx 35\pm 10\mbox{\,fb}$ [14, 15].

The LHCb Collaboration [10] published upper limits at 95% confidence level on the ratio $\sigma(\Xi_{cc}^{+})\cdot Br(\Xi_{cc}^{+}\to\Lambda_{c}^{+}K^{-}\pi^{+})/\sigma(\Lambda_{c}^{+})$ of $1.5\times 10^{-2}$ and $3.9\times 10^{-4}$ for the lifetimes of 100 fs and 400 fs, respectively, and for an integrated luminosity of $0.65\mbox{\,fb}^{-1}$. This is comparable with results from Refs. [14, 16, 17, 18] of about $10^{-4}-10^{-3}$. However, the minimum lifetime reached by the LHCb is about three times lager than the one measured by the SELEX experiment which is $\tau(\Xi_{cc}^{+})<33$ fs at 90% confidence level [1], and almost two times larger than the theoretical prediction of $\tau(\Xi_{cc}^{+})\approx 53$ fs [19]. In other words, the LHCb Collaboration provided an analysis outside the signal region.

The production of doubly heavy baryons at high Feynman-$x$ at the scheduled future fixed-target experiment at the LHC (AFTER@LHC) via the intrinsic heavy quark mechanism is already discussed in Refs. [11, 20]. In addition, some existing experiments have fixed-target programs [21, 22, 23]. However, as these experiments were built as collider detectors, they have very limited access to high Feynman-$x$.

In a recent talk [24], Stanley Brodsky proposed the production of charmed hadrons via the intrinsic charm mechanism from the target (see also chapter 6.1 in Ref [25]). In this paper we investigate the opportunity to produce doubly charmed baryons via the intrinsic heavy quark mechanism from the target at the fixed-target experiment at STAR [21], the fixed-target detector at the LHCb experiment (SMOG@LHCb) [22, 23] and at the scheduled future fixed-target experiment at the LHC (AFTER@LHC) [26, 27].

2 Revisiting the SELEX data

The SELEX experiment was a fixed-target experiment utilizing the Fermilab negative and positive charged beams at $600\mbox{\,GeV}/c$ to produce charm particles in a set of thin foil of Cu or in a diamond. It was operated in the kinematic region $x_{F}>0.1$. The negative beam composition was about 50% $\Sigma^{-}$ and 50% $\pi^{-}$ while the positive beam was composed of 90% protons. The experimental data recorded used both positive and negative beams. 67% of the events were induced by $\Sigma^{-}$, 13% by $\pi^{-}$, and 18% by protons.

The production cross section was not provided by the SELEX collaboration. However, the production properties of the doubly charmed baryons can be compared to that of the $\Lambda_{c}^{+}$ baryon. To simplify the analysis we will take a look only at the production of the 20 signal events for the decay mode $\Xi_{cc}^{++}\to\Lambda_{c}^{+}K^{-}\pi^{+}\pi^{+}$ at a mass of $3.76\mbox{\,GeV}$ over a sample of 1656 events for $\Lambda_{c}^{+}\to pK^{-}\pi^{+}$ [2] with $x_{F}(\Lambda_{c}^{+})>0.15$. This sample was previously used for a precision measurement of the lifetime of $\Lambda_{c}^{+}$ [28, 29]. The measured production ratio $R$ is defined by

[TABLE]

where $N$ is the number of events in the respective sample and the reconstruction efficiency of $\Xi_{cc}^{++}$ is given by $1/\epsilon_{++}\simeq 3.7$ [2]. Central values for the number $N_{\Lambda_{c}^{+}}/\epsilon_{\Lambda_{c}^{+}}$ of corrected events can be found in Ref. [30] to lie between $13326$ and $10010$ according whether the lowest bin with $x_{F}\in[0.125,0.175]$ is taken into account or not. However, using the intrinsic charm mechanism as the production mechanism, $\epsilon_{++}$ will be at least $2.3$ times bigger (see Figs. 1 and 2 for the acceptance for perturbative QCD and intrinsic charm, respectively). Therefore, we obtain

[TABLE]

This value is at least an order of magnitude smaller than the production ratio in the sample $N_{\Xi_{cc}^{++}}/\epsilon_{++}\cdot 1/N_{\Lambda_{c}^{+}}\approx 0.045$ which is usually erroneously used as the production ratio.

3 Theoretical background

Origin, properties and possibilities for the detection of a non-perturbative intrinsic heavy flavor component in the nucleon are widely discussed in the literature [31, 32, 33, 34, 35, 26, 27, 36, 37]. QCD predicts such components from the outset. Intrinsic charm and bottom quarks are contained in the wavefunction of a light hadron and originate from diagrams where the heavy quarks are multiply attached via gluons to the valence quarks. Intrinsic heavy flavor components are contributed by the twist-six contribution of the operator product expansion proportional to $1/m_{Q}^{2}$ [31, 32]. In this case, the frame-independent light-front wavefunction of the light hadron has maximum probability if the Fock state is minimally off-shell. This means that all the constituents are at rest in the hadron rest frame and thus have the same rapidity $y$ if the hadron is boosted. Equal rapidity occurs if the light-front momentum fractions $x=k^{+}/P^{+}$ of the Fock state constituents are proportional to their transverse masses, $x_{i}\propto m_{T,i}=(m^{2}_{i}+k^{2}_{T,i})^{1/2}$, i.e. if the heavy constituents have the largest momentum fractions. This features the BHPS model given by Brodsky, Hoyer, Peterson and Sakai for the distribution of intrinsic heavy quarks [33, 34].

In the BHPS model the wavefunction of a hadron in QCD can be represented as a superposition of Fock state fluctuations, e.g. $|h\rangle\sim|h_{l}\rangle+|h_{l}g\rangle+|h_{l}Q\bar{Q}\rangle\ldots$, where $h_{l}$ is the light quark content, and $Q=c,b$. If the projectile interacts with the target, the coherence of the Fock components is broken and the fluctuation can hadronize. The intrinsic heavy flavor Fock components are generated by virtual interactions such as $gg\to Q\bar{Q}$ where the gluons couple to two or more valence quarks of the projectile. The probability to produce such $Q\bar{Q}$ fluctuations scales as $\alpha_{s}^{2}(m_{Q}^{2})/m_{Q}^{2}$ relative to the leading-twist production.

Following Refs. [33, 34, 35], the general formula for the probability distribution of two heavy quark pairs by intrinsic heavy flavor Fock state as a function of the momentum fractions $x_{i}$ is given by

[TABLE]

where $\hat{m}_{i}=(m_{i}^{2}+\langle k^{2}_{T,i}\rangle)^{1/2}$ is the effective mass, $k^{2}_{T,i}$ is the mean transverse momenta, and the masses of the light quarks are neglected.

The normalization of the production cross section of two charm pairs is given by

[TABLE]

where $\mu^{2}\approx 0.2\mbox{\,GeV}^{2}$ denotes the squared soft interaction scale parameter, the effective transverse $c$-quark mass is given by $\hat{m}_{c}=1.5\mbox{\,GeV}$, and $\sigma^{\it inel}$ is the inelastic proton–proton cross section.

The nuclear dependence scaling from the manifestation of intrinsic charm is expected to be $\sigma_{A}\approx\sigma_{icc}\cdot A^{2/3}$ for production from the beam and $\sigma_{A}\approx\sigma_{icc}\cdot A$ for production from the target [24].

4 Production of the doubly charmed baryons

The production cross section of the doubly charmed baryon can be obtained as an application of the principle of quark–hadron duality. According to this principle, the cross section of the baryon is obtained by calculating the production of a $QQ$ pair in the small invariant mass interval between $2m_{Q}$ and the threshold to produce open heavy quark hadrons, $2m_{H}$. The $QQ$ pair has $3\times 3=(\bar{3}+6)$ color components, consisting of a color antitriplet and a color sextet. The probability that a $QQ$ pair forms an antitriplet state is $\bar{3}/(\bar{3}+6)=1/3$. Therefore, in case of the doubly charmed baryon the cross section will be

[TABLE]

where $f_{cc}$ is the fragmentation ratio of the $cc$ pair written as

[TABLE]

Not all $cc$ pairs form a doubly charmed baryon. Unfortunately, the fragmentation rate is unknown. However, if we consider the diquark $cc$ as heavy antiquark and analyze the fragmentation of $c$ and $b$ quarks into mesons, for example $f(c\to D^{+})=0.217\pm 0.043$ [38] or $f(b\to B^{+},B^{0})=0.344\pm 0.021$ [39], we can assume the doubly charmed baryonic cross section to be

[TABLE]

The STAR fixed-target program is a fixed-target experiment using the proton beam of the Relativistic Heavy Ion Collider (RHIC) up to $250\mbox{\,GeV}/c$ and the Au beam up to $100\mbox{\,GeV}/c$ colliding with a wired target. Combining Eqs. (2), (3) and (4) and using $\sigma^{\it inel}_{pp}(p_{\rm beam}=200\mbox{\,GeV})\approx 32\text{\,mb}$ [40], we may expect the production cross section of the $\Xi_{cc}$ to be

[TABLE]

The kinematic limits on the energy and the momentum of the doubly charmed baryon formed by the intrinsic charm from the target are given by

[TABLE]

These last expressions depend solely on the two masses $m_{cc}$ and $m_{\it tar}$ and no longer on the beam energy. Upon combining Eqs. (1) and (6) we can find the momentum distribution (Fig. 3) and the distribution of the rapidity difference $\Delta y=y-y_{\it tar}$ (Fig. 4) in the laboratory frame. It is obvious that experiments at the STAR detector, typical at rapidities $|y|<1$ for track selection, have the potential to observe doubly charmed baryons.

The SMOG@LHCb is a fixed-target experiment using the LHC beam at $6500\mbox{\,GeV}/c$ dumped in the Helium gas target. Following the logic from above and using $\sigma^{\it inel}_{pp}(\sqrt{s}=110\,\text{GeV})\approx 37\mbox{\,mb}$ [41] we obtain the production cross section of the $\Xi_{cc}$ to be

[TABLE]

The kinematic distributions have the same shape with slightly different mean values. Unfortunately, SMOG@LHCb has acceptance only for $\Delta y>2$ which will make the detection of the doubly charmed baryons problematic.

AFTER@LHC is the scheduled future fixed-target experiment at the LHC operating at $\sqrt{s}=115\mbox{\,GeV}$. Therefore, the estimation provided for SMOG@LHCb can also be used for AFTER@LHC.

5 Conclusion

In this paper we investigated the beautiful prediction of the intrinsic heavy quark mechanism. The doubly charmed baryons are produced from the target with an approximate mean value for the momentum of about $1.5\mbox{\,GeV}/c$. Such “soft” final states can be observed at the current and future fixed-target experiments. The production cross sections are presented.

From the calculation we see that the intrinsic heavy quark mechanism does not contribute to the region of negative $x_{F}$. The double intrinsic heavy quark mechanism is not the leading production mechanism at the modern accelerators [42]. However, it still can aim to searching exotic states like double charmed baryons. In addition, we reinterpreted the SELEX data and obtained a realistic production ratio for doubly charmed baryons $(ccu)$ over $\Lambda_{c}^{+}$ in the kinematic region of SELEX (see the discussion in Sec. 2).

In the end it is interesting to note that using the fragmentation ratios of $c$ and $b$ quarks into baryons, $f(c\to\Lambda_{c}^{+})=0.071\pm 0.21$ [43] and $f(b\to b\text{-baryons})=0.197\pm 0.046$ [39] we can also roughly calculate the doubly, hidden and open charm, tetraquark production cross section as $0.1\cdot\sigma(cc)$.

Acknowledgments

We would like to thank S.J. Brodsky for very detailed and productive discussions on the production of heavy quark states from the target via the intrinsic heavy quark mechanism. In addition, we would like to thank Guannan Xie for comments on acceptance of the STAR detector and Yu. Shcheglov for comments on the acceptance of the LHCb detector. This work was supported by the Estonian Research Council under Grant No. IUT2-27.