Hadronic production of the doubly charmed baryon via the proton–nucleus and the nucleus–nucleus collisions at the RHIC and LHC

The European Physical Journal C, Oct 2018

We present a detailed discussion on the doubly charmed baryon \(\Xi _{cc}\) production at the RHIC and LHC via the proton–nucleus (p–N) and nucleus–nucleus (N–N) collision modes. The extrinsic charm mechanism via the subprocesses \(g+c\rightarrow (cc)[n]+\bar{c}\) and \(c+c\rightarrow (cc)[n]+g\) together with the gluon-gluon fusion mechanism via the subprocess \(g+g\rightarrow (cc)[n]+\bar{c}+\bar{c}\) have been taken into consideration, where the intermediate diquark is in \([n]=[^1S_0]_\mathbf{6}\)-state or \([^3S_1]_{\bar{\mathbf{3}}}\)-state, respectively. Total and differential cross sections have been discussed under various collision energies. To compare with the \(\Xi _{cc}\) production via proton-proton collision mode at the LHC, we observe that sizable \(\Xi _{cc}\) events can also be generated via p–N and N–N collision modes at the RHIC and LHC. For examples, about \(8.1\times 10^7\) and \(6.7\times 10^7\) \(\Xi _{cc}\) events can be accumulated in p-Pb and Pb-Pb collision modes at the LHC within one operation year.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-018-6283-1.pdf

Hadronic production of the doubly charmed baryon via the proton–nucleus and the nucleus–nucleus collisions at the RHIC and LHC

The European Physical Journal C October 2018, 78:801 | Cite as Hadronic production of the doubly charmed baryon via the proton–nucleus and the nucleus–nucleus collisions at the RHIC and LHC AuthorsAuthors and affiliations Gu ChenChao-Hsi ChangXing-Gang Wu Open Access Regular Article - Theoretical Physics First Online: 01 October 2018 101 Downloads Abstract We present a detailed discussion on the doubly charmed baryon \(\Xi _{cc}\) production at the RHIC and LHC via the proton–nucleus (p–N) and nucleus–nucleus (N–N) collision modes. The extrinsic charm mechanism via the subprocesses \(g+c\rightarrow (cc)[n]+\bar{c}\) and \(c+c\rightarrow (cc)[n]+g\) together with the gluon-gluon fusion mechanism via the subprocess \(g+g\rightarrow (cc)[n]+\bar{c}+\bar{c}\) have been taken into consideration, where the intermediate diquark is in \([n]=[^1S_0]_\mathbf{6}\)-state or \([^3S_1]_{\bar{\mathbf{3}}}\)-state, respectively. Total and differential cross sections have been discussed under various collision energies. To compare with the \(\Xi _{cc}\) production via proton-proton collision mode at the LHC, we observe that sizable \(\Xi _{cc}\) events can also be generated via p–N and N–N collision modes at the RHIC and LHC. For examples, about \(8.1\times 10^7\) and \(6.7\times 10^7\) \(\Xi _{cc}\) events can be accumulated in p-Pb and Pb-Pb collision modes at the LHC within one operation year. 1 Introduction Within the quark model, the doubly heavy baryon is regarded as a three quark state with two heavy quarks (c or b) and a light quark q (\(q=u, d, s\)) [1, 2, 3, 4, 5]. The doubly heavy baryons are important for the understanding of Quantum Chromodynamics (QCD) theory. For convenience, throughout the paper, we adopt \(\Xi _{QQ'}\) as short notation for the baryon \(\Xi _{QQ'q}\), where Q and \(Q'\) stand for the heavy c or b quark, respectively. In year 2000, the SELEX collaboration [6, 7] reported the observation of \(\Xi _{cc}^+\) via its decay channels \(\Xi _{cc}^+ \rightarrow \Lambda _c^+K^-\pi ^+\) and \(\Xi _{cc}^+ \rightarrow pD^+K^-\). Later one, more experimental measurements have been carried out to confirm this observation by the FOCUS [8], the BABAR [9], the Belle [10, 11], and the LHCb collaboration [12]. However, all of those experiments were fail to reproduce the SELEX observation. In year 2017 the LHCb collaboration released their first observation of \(\Xi _{cc}^{++}\) [13] via its weak decay channel \(\Xi _{cc}^{++} \rightarrow \Lambda _c^+K^-\pi ^+\pi ^+\). This stimulates many new works, either experimentally or theoretically, for the doubly heavy baryons. The LHCb observation is based on the simulation by using a dedicated generator GENXICC [14, 15, 16], which is designed to simulate the doubly heavy baryon production via proton-proton (pp) collision at the Large Hadron Collider (LHC). Theoretically, the production of doubly heavy baryons at different types of high-energy colliders, such as the \(e^+e^-\), the electron-proton (ep), and the pp (or \(p\bar{p}\)) colliders, has been studied in detail in many works, cf. Refs.[17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41]. At the hadronic colliders, those works culminated in the generator GENXICC. This generator not only contains gluon-gluon fusion and the extrinsic charm mechanisms via the \(g+g\rightarrow \Xi _{cc}+X\), \(g+c\rightarrow \Xi _{cc}+X\) and \(c+c\rightarrow \Xi _{cc}+X\) subprocesses, but also contains various contributions from high Fock states in the baryons.1 Similar to the case of \(B_c\)-meson production [47, 48, 49], it is interesting to show that whether sizable number of \(\Xi _{cc}^{++}\) events can also be produced at the heavy ion colliders, such as the STAR experiment at the Relativistic Heavy Ion Collider (RHIC) and a large ion collider experiment (ALICE) at the LHC. In the heavy ion collisions, the \(\Xi _{cc}^{++}\) baryons can be produced via the proton–nucleus (p–N) and the nucleus–nucleus (N–N) collision modes, respectively. The production of doubly heavy baryon via p–N or N–N collision should also provide an alternative candidate of studying on the quark-gluon plasma (QGP), similar to the case of doubly heavy meson production via heavy ion collision [49]. Through high-energy p–N and N–N collisions, the generated heavy quarks will either combine and evolve into doubly heavy baryon before QGP formation or hadronize into colorless baryons via the transition from the QGP phase to the hadronic phase. This makes the production of \(\Xi _{cc}^{++}\) baryons via the p–N and N–N modes quite different from the usual pp-collision mode. Especially, the nuclear effects, e.g., the shadowing effect and the modifications of nuclear parton density functions (nPDFs) etc., shall play a significant role. Consequently, by measuring the properties of the doubly heavy baryon, one may achieve important information about QGP and nuclear properties. The remaining parts of the paper are organized as follows. In Sect. 2, we present the calculation technology. In Sect. 3, we present the numerical results for the \(\Xi _{cc}^{++}\) production via p–N or N–N collision model, and the results for the pp collision mode are also presented for comparison. Sect. 4 is reserved for a summary. 2 Explanation of the calculation technology The production of \(\Xi _{cc}\) baryon can be factorized into two steps: The first step is to produce two \(c\bar{c}\) pairs. This step is pQCD calculable, since the intermediate gluon should be hard enough to generate a heavy \(c\bar{c}\) pair. The second step is to make the heavy c-quarks into a bounding (cc)-diquark in the spin-triplet (\([^3S_1]\)) or spin-singlet (\([^1S_0]\)), and the \({\bar{\mathbf {3}}}\) (or \(\mathbf {6}\)) color state, accordingly. More explicitly, the intermediate diquark in \(\Xi _{cc}\) has two spin-and-color configurations, \([^3S_1]_{\bar{\mathbf{3}}}\) and \([^1S_0]_\mathbf{6}\). It will then be hadronized into a \(\Xi _{cc}\) baryon via fragmentation, whose probability is characterized by the non-perturbative matrix element: 1) We adopt the usual assumption that the diquark shall evolve into the baryon with \(100\%\) probability2; 2) The intermediate \((c\bar{c})\)-diquark shall grab a light quark with possible soft gluons from the hadron to form the final colorless baryon with a relative possibility for different light quarks as \(u :d :s \simeq 1:1:0.3\) [51]. If the diquark \((cc)[^3S_1]_{\bar{3}}\) is produced, then it will fragment into \(\Xi _{cc}^{++}\) with \(43\%\) probability, \(\Xi _{cc}^{+}\) with \(43\%\) probability, and \(\Omega _{cc}^{+}\) with \(14\%\) probability. Consequently, the production of \(\Xi _{cc}\) baryon can be expressed as [46], $$\begin{aligned}&d\sigma (\mathrm{A}\mathrm{B} \rightarrow {\Xi _{cc}}+X) = \sum _{i,j=g,c}\sum _{n} \int dx_1 dx_2 (N_A f^{h/\mathrm{A}}_{i}\nonumber \\&\quad \times (x_1, \mu _{f})) (N_B f^{h/\mathrm{B}}_{j}(x_2, \mu _{f})) d{\hat{\sigma }}_{ij \rightarrow (cc)[n]+ X} \langle \mathcal{O}^\mathcal{H}[n] \rangle ,\nonumber \\ \end{aligned}$$ (1) where the symbol [n] stands for the spin-and-color state of the (cc)-diquark, and the symbols \(\mathrm A\) and \(\mathrm B\) stand for p or N for the incident hadron to be proton or nucleus, respectively. The functions \((N_A f_{i, j}^{h/A})\) and \((N_B f_{i, j}^{h/B})\) with \((i,j=g,c)\) are effective nucleus parton distribution functions (PDFs) for the nucleus \(\mathrm A\) or \(\mathrm B\) accordingly, which stands for the parton density of bound-nucleon h in nucleus A and carries the fraction \(x_n\,(n=1,2)\) of the hadron momentum at the factorization scale \(\mu _f\). h stands for the nucleon, proton or neutron, respectively. Here \(N_\mathrm{A}\) or \(N_\mathrm{B}\) is the nuclear number in the incident nucleus. For examples, \(N_\mathrm{Au}=197\) for the gold nucleus (\(^{197}_{79}\mathrm Au\)), and \(N_\mathrm{Pb}=208\) for the lead nucleus (\(^{208}_{82}\mathrm Pb\)). Many PDF models have been suggested to study the heavy-ion collisions, such as the Heavy-Ion Jet INteraction Generator (HIJING) model [52], A Multiphase Transport (AMPT) model [53], the Monte Carlo Glauber Model [54, 55, 56, 57], and etc.. Here, following the same idea of CTEQ group [58], we adopt the PDF of bound-nucleon in nucleus as the heavy ion PDF. For the pp collision mode, we have \(N_\mathrm{A}=N_\mathrm{B}=1\), and the PDFs are reduced to the usual PDFs inside the free proton. For the case of \(N_\mathrm{A}\ne 1\) and \(N_\mathrm{B}\ne 1\), we need to consider the collision geometry and the spatial dependence of the shadowing parameterization effect [59, 60] to the nucleus PDF. More explicitly, for the case of p–N and N–N collisions, we need to know the nuclear gluon/charm distribution functions. For the nuclear A, its overall gluon/charm distribution function can be expressed as $$\begin{aligned} f^\mathrm{A}_{g(c)}(x_n, \mu _f) = \int d^2 \vec {r} dz \,f^\mathrm{A}_{g(c)}(x_n, \mu _{f}, \vec {r}, z), \end{aligned}$$ (2) where \(\vec {r}\) and z are transverse and longitudinal location of the parton in the coordinate space. Considering the nuclear effects and collision geometry, the nuclear densities \(f^\mathrm{A}_{g(c)}(x_n, \mu _f, \vec {r}, z)\) can be factorized as the product of the nucleon density in the nucleus \(\rho _\mathrm{A}(s)\), the free-nucleon density \(f^{h}_{g(c)}(x_n, \mu _f)\), and the shadowing ratio \(S^{g(c)}_\mathrm{P, S}(N_\mathrm{A}, x_n, \mu _f, \vec {r}, z)\), i.e. $$\begin{aligned}&f^\mathrm{A}_{g(c)}(x_n, \mu _f, \vec {r}, z) \nonumber \\&\quad = \rho _\mathrm{A}(s) S^{g(c)}_\mathrm{P, S}(N_\mathrm{A}, x_n, \mu _f, \vec {r}, z) f^{h}_{g(c)}(x_n, \mu _f), \end{aligned}$$ (3) where \(s=\sqrt{r^2+z^2}\), \(\rho _\mathrm{A}(s)\) is assumed to be Woods-Saxon distribution [61], which satisfies the normalization condition $$\begin{aligned} \int \rho _\mathrm{A}(s) d^2\vec {r} dz= & {} N_\mathrm{A}. \end{aligned}$$ (4) Then, we obtain $$\begin{aligned}&\int d^2 \vec {r} dz \rho _\mathrm{A}(s) S^{g(c)}_\mathrm{P, S}(N_\mathrm{A}, x_n, \mu _f, \vec {r}, z) \nonumber \\&\quad = N_\mathrm{A}\,S_\mathrm{P}^{g(c)}(N_\mathrm{A}, x_n, \mu _f). \end{aligned}$$ (5) Consequently, we have [62, 63] $$\begin{aligned} f^\mathrm{A}_{g(c)}(x_n, \mu _f)= & {} N_\mathrm{A} \,S_\mathrm{P}^{g(c)}(N_\mathrm{A}, x_n, \mu _f) f^h_{g(c)}(x_n, \mu _f) \nonumber \\= & {} N_\mathrm{A}\,f_{g(c)}^{h/\mathrm A}(x_n, \mu _f). \end{aligned}$$ (6) Here \(f_{g(c)}^{h/\mathrm A}(x_n, \mu _f)\) is the wanted effective bound PDF for nucleon, which describes the gluon/charm density of the bound-nucleon h in nucleus A. We adopt the nCTEQ15 version to calculate \(f_{g(c)}^{h/\mathrm A}(x_n, \mu _f)\), which is fixed via a global fit by using the experimental data on nuclei all the way up to \(^{208}\mathrm Pb\) [58]. Furthermore, the \(\Xi _{cc}\) baryon can be expanded as a series of Fock states over the relative velocity (v) of the constituent heavy quarks in the baryon rest frame, $$\begin{aligned} \vert \Xi _{cc} \rangle = c_1 \vert (cc) q\rangle +c_2 \vert (cc) qg \rangle +c_3 \vert (cc) q gg \rangle +\cdots , \end{aligned}$$ where the expansion coefficients \(c_i(i=1,2,\cdots )\) are functions of v. \(\langle \mathcal{O}^\mathcal{H}[n] \rangle \) is the long-distance matrix element, which is proportional to the inclusive transition probability of the perturbative diquark state (cc)[n] into the heavy baryon \(\Xi _{cc}\). For convenience, we adopt the assumption that the transition probability for the color anti-triplet or color sextuplet of the (cc)-diquark are the same [23]. The non-perturbative long-distance matrix element can be related to the Schr\(\ddot{\mathrm{o}}\)dinger wavefunction at the origin as [38], \(\langle \mathcal{O}^\mathcal{H}[n] \rangle \simeq |\psi _{cc}(0)|^2\). \(d\hat{\sigma }_{ij \rightarrow (cc)[n]+ X}\) is the differential cross-section of the hard subprocess, which are different for different channels and shall be dealt with by using the generator GENXICC. 3 Numerical results and discussions To do the numerical calculation, we take \(|\Psi _{cc}(0)|^2 =0.039\) GeV\(^3\) [19], \(M_{\Xi _{cc}}=3.50\) GeV with \(m_c=M_{\Xi _{cc}}/2\), and the nCTEQ15 [58] as the nucleon PDF. The renormalization scale and the factorization scale are set to be the transverse mass of \(\Xi _{cc}\), \(M_t=\sqrt{M^2_{\Xi _{QQ'}}+p_t^2}\). For the collision energies, we adopt [49, 64]: \(\sqrt{S_{p\mathrm Pb}}= 8.16\;\mathrm TeV\) and \(\sqrt{S_\mathrm{PbPb}}= 5.02\;\mathrm TeV\) at the LHC, and \(\sqrt{S_{p\mathrm Au}}= 0.2\;\mathrm TeV\) and \(\sqrt{S_\mathrm{AuAu}}= 0.2\;\mathrm TeV\) at the RHIC. 3.1 Basic results for \(\Xi _{cc}\) production Table 1 Total cross sections (in unit \(\mu \)b) for \(\Xi _{cc}\) in pp, p–N, and N–N collision modes at the RHIC and the LHC   RHIC LHC \(\sqrt{S_\mathrm{NN}}\) (TeV) p-Au (0.2) Au–Au (0.2) pp (13) p-Pb (8.16) Pb–Pb (5.02) \(\sigma (gg\rightarrow (cc)_\mathbf{6}[^1S_0])\) \(1.05\times 10^{-1}\) \(2.30\times 10^{1}\) \(7.96\times 10^{-2}\) 9.20 \(1.04\times 10^{3}\) \(\sigma (gg\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\) \(5.79\times 10^{-1}\) \(1.27\times 10^{2}\) \(4.19\times 10^{-1}\) \(4.79\times 10^{1}\) \(5.38\times 10^{3}\) \(\sigma (gc\rightarrow (cc)_\mathbf{6}[^1S_0])\) \(2.83\times 10^{-1}\) \(6.26\times 10^{1}\) \(9.40\times 10^{-2}\) \(1.17\times 10^{1}\) \(1.19\times 10^{3}\) \(\sigma (gc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\) 2.58 \(5.73\times 10^{2}\) \(8.64\times 10^{-1}\) \(9.31\times 10^{1}\) \(1.09\times 10^{4}\) \(\sigma (cc\rightarrow (cc)_\mathbf{6}[^1S_0])\) \(4.06\times 10^{-5}\) \(9.00\times 10^{-3}\) \(8.80\times 10^{-5}\) \(1.05\times 10^{-2}\) 1.25 \(\sigma (cc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\) \(1.01\times 10^{-3}\) \(2.23\times 10^{-1}\) \(2.33\times 10^{-3}\) \(2.94\times 10^{-1}\) \(3.33\times 10^{1}\) Table 2 R value which shows the importance of extrinsic charm mechanisms   RHIC LHC \(\sqrt{S_\mathrm{NN}}\) (TeV) p-Au(0.2) Au–Au(0.2) p-Pb(8.16) Pb–Pb(5.02) R   6.13 6.19 3.39 3.44 In Table 1, we present the total cross sections for the production of \(\Xi _{cc}\) via the p–N and N–N collision modes at the RHIC and the LHC. By summing up contributions from \(g+g\), \(g+c\) and \(c+c\) channels and different spin-and-color configurations of the intermediate (cc)-diquark together, we obtain $$\begin{aligned}&\left. \sigma ^\mathrm{tot}_{p\mathrm Au}(\Xi _{cc})\right| _\mathrm{RHIC} = 3.55\;\mathrm{\mu b}, \end{aligned}$$ (7) $$\begin{aligned}&\left. \sigma ^\mathrm{tot}_\mathrm{AuAu}(\Xi _{cc}) \right| _\mathrm{RHIC} = 7.85\times 10^{2}\;\mathrm{\mu b}, \end{aligned}$$ (8) $$\begin{aligned}&\left. \sigma ^\mathrm{tot}_{pp}(\Xi _{cc})\right| _\mathrm{LHC} = 1.46\;\mathrm{\mu b}, \end{aligned}$$ (9) $$\begin{aligned}&\left. \sigma ^\mathrm{tot}_{p\mathrm Pb}(\Xi _{cc})\right| _\mathrm{LHC} = 1.62\times 10^{2}\;\mathrm{\mu b}, \end{aligned}$$ (10) $$\begin{aligned}&\left. \sigma ^\mathrm{tot}_\mathrm{PbPb}(\Xi _{cc})\right| _\mathrm{LHC} = 1.85\times 10^4\;\mathrm{\mu b}. \end{aligned}$$ (11) To compare with \(\Xi _{cc}\) production via the pp collision, we observe that the total cross sections of \(\Xi _{cc}\) are enhanced by about 2–4 orders of magnitude in p–N and N–N collision modes at the RHIC and the LHC. At the RHIC, the designed luminosities are \(4.5\times 10^{29}\,\mathrm cm^{-2}s^{-1}\) and \(8.0\times 10^{27}\,\mathrm cm^{-2}s^{-1}\) for the p-Au and the Au-Au collisions.3 Thus, at the RHIC, we shall have \(1.6 \times 10^7\) \(\Xi _{cc}\) events to be generated for the p-Au collision in one operation year, \(6.3 \times 10^7\) \(\Xi _{cc}\) events for the Au-Au collision. At the LHC, the designed luminosities are \(5.0\times 10^{33}\,\mathrm cm^{-2}s^{-1}\), \(5.0\times 10^{29}\,\mathrm cm^{-2}s^{-1}\) and \(3.6\times 10^{27}\,\mathrm cm^{-2}s^{-1}\) for the pp, p-Pb and the Pb-Pb collisions. Thus, at the LHC, we shall have \(7.3 \times 10^{10}\) \(\Xi _{cc}\) events to be generated in pp collision in one operation year, \(8.1 \times 10^7\) \(\Xi _{cc}\) events to be generated in p-Pb collision, \(6.7 \times 10^7\) \(\Xi _{cc}\) events for the Pb-Pb collision. Total cross-sections at the RHIC and LHC are in the following sequential order, \(\sigma _{pp}<\sigma _{p\mathrm{A}}<\sigma _\mathrm{AA}\); due to the shadowing effect, the relative ratio among those cross-sections is smaller than the naive ratio, \(1:N_\mathrm{A}:N^2_\mathrm{A}\), with A=Au for RHIC and A=Pb for LHC. At the LHC, the \(\Xi _{cc}\) events generated via p-Pb and Pb-Pb collisions are about three-orders lower than those via pp collision, which are still sizable; thus the STAR and ALICE experiments at the RHIC and LHC can also be potential platform for studying the properties of \(\Xi _{cc}\). For the p-Pb (p-Au) collision at the LHC (RHIC), Table 1 shows that the contribution from \((cc)_ {\bar{\mathbf{3}}}[^3S_1]\) is larger than that of \((cc)_\mathbf{6}[^1S_0]\) by about five (six) times for \(g+g\)-mechanism, eight (nine) times for \(g+c\)-mechanism, and twenty-eight (twenty-five) times for \(c+c\)-mechanism. The conditions for the N–N collision are similar. Even though the contribution from \((cc)_\mathbf{6}[^1S_0]\) is small, one should take it into consideration for an accurate prediction. By summing up the contributions from \((cc)_\mathbf{6}[^1S_0]\) and \((cc)_{\bar{\mathbf{3}}}[^3S_1]\) configurations together, we obtain the relative importance of different production mechanisms $$\begin{aligned}&\left. \sigma _{g+g}^{p\mathrm Au}: \sigma _{g+c}^{p\mathrm Au}: \sigma _{c+c}^{p\mathrm Au}\right| _\mathrm{RHIC} \simeq 647 : 2720 : 1, \end{aligned}$$ (12) $$\begin{aligned}&\left. \sigma _{g+g}^\mathrm{AuAu}: \sigma _{g+c}^\mathrm{AuAu}: \sigma _{c+c}^\mathrm{AuAu}\right| _\mathrm{RHIC} \simeq 645 : 2735 : 1, \end{aligned}$$ (13) $$\begin{aligned}&\left. \sigma _{g+g}^{p\mathrm Pb}: \sigma _{g+c}^{p\mathrm Pb}: \sigma _{c+c}^{p\mathrm Pb}\right| _\mathrm{LHC} \simeq 187 : 344 : 1, \end{aligned}$$ (14) $$\begin{aligned}&\left. \sigma _{g+g}^\mathrm{PbPb}: \sigma _{g+c}^\mathrm{PbPb}: \sigma _{c+c}^\mathrm{PbPb}\right| _\mathrm{LHC} \simeq 187 : 351 : 1. \end{aligned}$$ (15) Contributions from the \(g+c\) mechanism are larger than the usually considered gluon-gluon fusion mechanism. As shall be shown below, this is caused by the fact that the cross-section for \(g+c\) mechanism is larger than that of \(g+g\) mechanism in small \(p_t\) region. To see how the extrinsic charm mechanism affects the \(\Xi _{cc}\) production, we define a ratio $$\begin{aligned} R = \frac{\sigma _{tot}}{\sigma _{gg \rightarrow \Xi _{cc}(cc)_ {\bar{\mathbf{3}}}[^3S_1]}}, \end{aligned}$$ (16) where \(\sigma _{tot}\) stands for the summation of total cross sections for all the considered production mechanisms and diquark configurations, and \(\sigma _{gg \rightarrow \Xi _{cc}(cc)_{ \bar{\mathbf{3}}}[^3S_1]}\) is the cross section for usually considered gluon-gluon fusion via \(gg \rightarrow \Xi _{cc}(cc)_{\bar{\mathbf{3}}}[^3S_1]\). The R values are put in Table 2. Table 2 shows that the extrinsic charm mechanism plays significant role via the p–N and the N–N collisions at the RHIC and the LHC, respectively. 3.2 Differential distributions of the \(\Xi _{cc}\) production via p–N and N–N collisions Table 3 Total cross sections (in unit: \(\mu \)b) for the \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms via p-Au and Au–Au collisions under various transverse momentum (\(p_t\)) cuts at the RHIC   p-Au (0.2 TeV) Au–Au (0.2 TeV)   \(p_t>2\,\mathrm GeV\) \(p_t>4\,\mathrm GeV\) \(p_t>6\,\mathrm GeV\) \(p_t>2\,\mathrm GeV\) \(p_t>4\,\mathrm GeV\) \(p_t>6\,\mathrm GeV\) \(\sigma (gg\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(4.79\times 10^{-2}\) \(8.58\times 10^{-3}\) \(1.36\times 10^{-3}\) \(1.05\times 10^{1}\) 1.88 \(2.89\times 10^{-1}\) \(\sigma (gg\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(2.21\times 10^{-1}\) \(3.08\times 10^{-2}\) \(4.12 \times 10^{-3}\) \(4.91\times 10^{1}\) 6.74 \(8.68\times 10^{-1}\) \(\sigma (gc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(7.17\times 10^{-2}\) \(7.51\times 10^{-3}\) \(8.93 \times 10^{-4}\) \(1.70 \times 10^{1}\) 1.76 \(2.00\times 10^{-1}\) \(\sigma (gc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(6.12\times 10^{-1}\) \(4.41\times 10^{-2}\) \(4.33\times 10^{-3}\) \(1.46\times 10^{2}\) \(1.04\times 10^{1}\) \(9.53\times 10^{-1}\) \(\sigma (cc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(4.05\times 10^{-5}\) \(4.05\times 10^{-5}\) \(3.30\times 10^{-5}\) \(8.96\times 10^{-3}\) \(8.96\times 10^{-3}\) \(7.24\times 10^{-3}\) \(\sigma (cc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(1.01\times 10^{-3}\) \(1.01 \times 10^{-3}\) \(8.64\times 10^{-4}\) \(2.22\times 10^{-1}\) \(2.22\times 10^{-1}\) \(1.89\times 10^{-1}\) Table 4 Total cross sections (in unit: \(\mu \)b) for the \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms via p-Pb and Pb–Pb collisions under various transverse momentum (\(p_t\)) cuts at the LHC   p-Pb (8.16 TeV) Pb–Pb (5.02 TeV)   \(p_t>2\,\mathrm GeV\) \(p_t>4\,\mathrm GeV\) \(p_t>6\,\mathrm GeV\) \(p_t>2\,\mathrm GeV\) \(p_t>4\,\mathrm GeV\) \(p_t>6\,\mathrm GeV\) \(\sigma (gg\rightarrow (cc)_\mathbf{6}[^1S_0])\)   5.64 1.95 \(6.38\times 10^{-1}\) \(6.39\times 10^{2}\) \(2.20\times 10^{2}\) \(7.08\times 10^{1}\) \(\sigma (gg\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(2.64\times 10^{1}\) 7.81 2.27 \(2.99\times 10^{3}\) \(8.86\times 10^{2}\) \(2.53\times 10^{2}\) \(\sigma (gc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   3.81 \(7.51\times 10^{-1}\) \(1.72\times 10^{-1}\) \(4.33\times 10^{2}\) \(8.48\times 10^{1}\) \(1.98\times 10^{1}\) \(\sigma (gc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(3.42\times 10^{1}\) 5.08 1.00 \(3.85\times 10^{3}\) \(5.73\times 10^{2}\) \(1.14\times 10^{2}\) \(\sigma (cc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(1.05\times 10^{-2}\) \(1.05\times 10^{-2}\) \(9.46\times 10^{-3}\) 1.25 1.25 1.16 \(\sigma (cc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(2.94\times 10^{-1}\) \(2.94\times 10^{-1}\) \(2.67\times 10^{-1}\) \(3.32\times 10^{1}\) \(3.32\times 10^{1}\) \(3.13\times 10^{1}\) Table 5 Total cross sections (in unit: \(\mu \)b) for the \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms via p-Au and Au–Au collisions under various rapidity (y) cuts at the RHIC   p-Au (0.2 TeV) Au–Au (0.2 TeV) \(y_\mathrm{cut}\) \(|y|<1\) \(|y|<2\) \(|y|<3\) \(|y|<1\) \(|y|<2\) \(|y|<3\) \(\sigma (gg\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(6.72\times 10^{-2}\) \(9.94\times 10^{-2}\) \(1.05\times 10^{-1}\) \(1.54\times 10^{1}\) \(2.22\times 10^{1}\) \(2.30\times 10^{1}\) \(\sigma (gg\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(3.72\times 10^{-1}\) \(5.47\times 10^{-1}\) \(5.75\times 10^{-1}\) \(8.56\times 10^{1}\) \(1.22\times 10^{2}\) \(1.26\times 10^{2}\) \(\sigma (gc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(1.42\times 10^{-1}\) \(2.44\times 10^{-1}\) \(2.79\times 10^{-1}\) \(4.66\times 10^{1}\) \(6.15\times 10^{1}\) \(6.30\times 10^{1}\) \(\sigma (gc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   1.23 2.21 2.54 \(4.28\times 10^{2}\) \(5.61\times 10^{2}\) \(5.75\times 10^{2}\) \(\sigma (cc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(3.16\times 10^{-5}\) \(4.02\times 10^{-5}\) \(4.05\times 10^{-5}\) \(7.45\times 10^{-3}\) \(8.95\times 10^{-3}\) \(8.96\times 10^{-3}\) \(\sigma (cc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(7.88\times 10^{-4}\) \(1.00\times 10^{-3}\) \(1.01\times 10^{-3}\) \(1.84\times 10^{-1}\) \(2.22\times 10^{-1}\) \(2.22\times 10^{-1}\) Table 6 Total cross sections (in unit: \(\mu \)b) for the \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms via p-Pb and Pb–Pb collisions under various rapidity (y) cuts at the LHC   p-Pb (8.16 TeV) Pb–Pb (5.02 TeV) \(y_\mathrm{cut}\) \(|y|<1\) \(|y|<2\) \(|y|<3\) \(|y|<1\) \(|y|<2\) \(|y|<3\) \(\sigma (gg\rightarrow (cc)_\mathbf{6}[^1S_0])\)   2.53 4.96 6.96 \(2.99\times 10^{2}\) \(5.83\times 10^{2}\) \(8.20\times 10^{2}\) \(\sigma (gg\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(1.33\times 10^{1}\) \(2.59\times 10^{1}\) \(3.62\times 10^{1}\) \(1.54\times 10^{3}\) \(3.02\times 10^{3}\) \(4.28\times 10^{3}\) \(\sigma (gc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   5.17 7.51 9.36 \(4.49\times 10^{2}\) \(7.25\times 10^{2}\) \(9.42\times 10^{2}\) \(\sigma (gc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(3.30\times 10^{1}\) \(5.36\times 10^{1}\) \(7.06\times 10^{1}\) \(4.10\times 10^{3}\) \(6.53\times 10^{3}\) \(8.52\times 10^{3}\) \(\sigma (cc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(4.54\times 10^{-3}\) \(8.39\times 10^{-3}\) \(9.95\times 10^{-3}\) \(5.86\times 10^{-1}\) 1.07 1.22 \(\sigma (cc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(1.21\times 10^{-1}\) \(2.35\times 10^{-1}\) \(2.83\times 10^{-1}\) \(1.56\times 10^{1}\) \(2.86\times 10^{1}\) \(3.25\times 10^{1}\) At the hadronic colliders, the small \(p_t\) and/or large rapidity y (the produced baryons move very close to the beam direction) cannot be detected by the detectors directly, and such kind of events cannot be utilized for experimental studies in common cases. In this subsection, we perform detailed calculations and discussions under different \(p_t\) cuts and y cuts. We adopt three typical \(p_t\) cuts, \(p_t>2\) GeV, \(p_t>4\) GeV and \(p_t>6\) GeV, to show how the production cross sections change with the \(\Xi _{cc}\) transverse momentum. The results are presented in Tables 3 and 4. Table 3 shows that the total cross sections for the \(\Xi _{cc}\) via p-Au and Au-Au collisions at the RHIC shall be reduced by about \(54\% \sim 99\%\) when the \(p_t\)-cut varying from 2 GeV to 6 GeV. Similarly, the production rates via p-Pb and Pb-Pb collisions shall be reduced by about \(38\% \sim 98\%\). This shows that small \(p_t\) region provides significant contributions to the \(\Xi _{cc}\) production via p–N and N–N collisions at the RHIC and the LHC. We adopt three typical rapidity cuts, \(|y|<1\), \(|y|<2\), and \(|y|<3\), to show how the production cross sections change with the \(\Xi _{cc}\) rapidity. The results are presented in Tables 5 and 6. Table 5 shows that at the RHIC, the \(\Xi _{cc}\) events mainly distribute in the rapidity region of \(y \in [-3,3]\); At the LHC, the \(\Xi _{cc}\) events distribute within a broader rapidity range. Open image in new window Fig. 1 The \(p_t\)-distributions of the \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms via p-Au and Au-Au collisions with \(\sqrt{S_{p\mathrm Au}}=0.2\) TeV and \(\sqrt{S_\mathrm{AuAu}}=0.2\) TeV at the RHIC Open image in new window Fig. 2 The \(p_t\)-distribution of the \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms via p-Pb and Pb-Pb collisions with \(\sqrt{S_{p\mathrm Pb}}=8.16\) TeV and \(\sqrt{S_\mathrm{PbPb}}=5.02\) TeV at the LHC We present the \(\Xi _{cc}\) \(p_t\) distributions via p–N and N–N collisions at the RHIC and the LHC in Figs. 1, 2. For each production mechanism, contributions from \((cc)_\mathbf{6}[^1S_0]\) and \((cc)_{\bar{\mathbf{3}}}[^3S_1]\) configurations have been summed for convenience. The \(p_t\) distributions at the RHIC and LHC have similar shapes, which shall first increase with the increment of \(p_t\) in small \(p_t\) region and then decrease quickly in large \(p_t\) region. In small \(p_t\) region, the \(g+c\) mechanism is larger than \(g+g\) mechanism via p–N and N–N collisions. If the experimental measurements can be extended to small \(p_t\) region, then one may study extrinsic charm mechanism by measuring the \(\Xi _{cc}\) events. Open image in new window Fig. 3 The y-distributions of the \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms via p-Au and Au-Au collisions with \(\sqrt{S_{p\mathrm Au}}=0.2\) TeV and \(\sqrt{S_\mathrm{AuAu}}=0.2\) TeV at the RHIC Open image in new window Fig. 4 The y-distributions of the \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms via p-Pb and Pb-Pb collisions with \(\sqrt{S_{p\mathrm Pb}}=8.16\) TeV and \(\sqrt{S_\mathrm{PbPb}}=5.02\) TeV at the LHC Open image in new window Fig. 5 The pseudo-rapidity (\(y_p\))-distributions of the \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms via p-Au and Au-Au collisions with \(\sqrt{S_{p\mathrm Au}}=0.2\) TeV and \(\sqrt{S_\mathrm{AuAu}}=0.2\) TeV at the RHIC Open image in new window Fig. 6 The pseudo-rapidity (\(y_p\))-distributions of the \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms via p-Pb and Pb-Pb collisions with \(\sqrt{S_{p\mathrm Pb}}=8.16\) TeV and \(\sqrt{S_\mathrm{PbPb}}=5.02\) TeV at the LHC We present the rapidity (y) and pseudo-rapidity (\(y_p\)) distributions of the \(\Xi _{cc}\) production at the RHIC and LHC in Figs. 3, 4, 5, and 6, respectively. For the dominant \(g+g\) and \(g+c\) mechanisms, there are plateaus for \(|y| \le 2\) or \(|y_p| \le 2\) at the RHIC, which become broader at the LHC, i.e. \(|y| \le 5\) and \(|y_p| \le 5\). 3.3 Theoretical uncertainties for \(\Xi _{cc}\) production The non-perturbative matrix elements are overall parameters, thus their uncertainties can be conveniently suppressed when we know their exact values. In this subsection, we discuss the other two important uncertainties, which are from the choices of charm quark mass and renormalization scale. For the purpose, we take \(m_c=1.75 \pm 0.1\,\mathrm GeV\) and three frequently used scales, e.g. the transverse mass \(M_t\) of \(\Xi _{cc}\), \(\sqrt{\hat{s}}\) and \(\sqrt{\hat{s}}/2\), where \(\sqrt{\hat{s}}\) is the center-of-mass energy of the corresponding subprocess, to do the calculation. For convenience, when considering one parameter’s uncertainty, the other parameters shall be kept to their central values. Total cross sections of \(\Xi _{cc}\) production for \(m_c=1.75 \pm 0.1\,\mathrm GeV\) are presented in Tables 7 and 8, which are for RHIC and LHC, respectively. Tables 7 and 8 indicate the total cross sections decrease with the increment of \(m_c\). The uncertainties at the RHIC via p-Au and Au-Au collisions with \(\sqrt{S_{p\mathrm{Au}}}=\sqrt{S_{\mathrm{Au Au}}}=0.2\mathrm{TeV}\) are $$\begin{aligned} \sigma _{gg\rightarrow (cc)_\mathbf{6}[^1S_0]}^{p\mathrm Au}= & {} \left( 1.05^{+0.67}_{-0.397}\right) \times 10^{-1}\;\mathrm{\mu b}, \end{aligned}$$ (17) $$\begin{aligned} \sigma _{gg\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1]}^{p\mathrm Au}= & {} \left( 5.79^{+3.73}_{-2.19}\right) \times 10^{-1}\;\mathrm{\mu b}, \end{aligned}$$ (18) $$\begin{aligned} \sigma _{gc\rightarrow (cc)_\mathbf{6}[^1S_0]}^{p\mathrm Au}= & {} \left( 2.83^{+1.24}_{-0.85}\right) \times 10^{-1}\;\mathrm{\mu b}, \end{aligned}$$ (19) $$\begin{aligned} \sigma _{gc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1]}^{p\mathrm Au}= & {} 2.58^{+1.12}_{-0.78}\;\mathrm{\mu b}, \end{aligned}$$ (20) $$\begin{aligned} \sigma _{cc\rightarrow (cc)_\mathbf{6}[^1S_0]}^{p\mathrm Au}= & {} \left( 4.06^{+3.57}_{-1.18}\right) \times 10^{-5}\;\mathrm{\mu b}, \end{aligned}$$ (21) $$\begin{aligned} \sigma _{cc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1]}^{p\mathrm Au}= & {} \left( 1.01^{+0.85}_{-0.309}\right) \times 10^{-3}\;\mathrm{\mu b}, \end{aligned}$$ (22) $$\begin{aligned} \sigma _{gg\rightarrow (cc)_\mathbf{6}[^1S_0]}^\mathrm{AuAu}= & {} \left( 2.30^{+1.46}_{-0.86}\right) \times 10^{1}\;\mathrm{\mu b}, \end{aligned}$$ (23) $$\begin{aligned} \sigma _{gg\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1]}^\mathrm{AuAu}= & {} \left( 1.27^{+0.800}_{-0.479}\right) \times 10^{2}\;\mathrm{\mu b}, \end{aligned}$$ (24) $$\begin{aligned} \sigma _{gc\rightarrow (cc)_\mathbf{6}[^1S_0]}^\mathrm{AuAu}= & {} \left( 6.26^{+2.44}_{-1.68}\right) \times 10^1\;\mathrm{\mu b}, \end{aligned}$$ (25) $$\begin{aligned} \sigma _{gc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1]}^\mathrm{AuAu}= & {} \left( 5.73^{+2.16}_{-1.58}\right) \times 10^{2}\;\mathrm{\mu b}, \end{aligned}$$ (26) $$\begin{aligned} \sigma _{cc\rightarrow (cc)_\mathbf{6}[^1S_0]}^\mathrm{AuAu}= & {} \left( 9.00^{+8.20}_{-2.65}\right) \times 10^{-3}\;\mathrm{\mu b}, \end{aligned}$$ (27) $$\begin{aligned} \sigma _{cc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1]}^\mathrm{AuAu}= & {} \left( 2.23^{+1.92}_{-0.69}\right) \times 10^{-1}\;\mathrm{\mu b}. \end{aligned}$$ (28) The uncertainties at the LHC via p-Pb and Pb-Pb collisions with \(\sqrt{S_{p\mathrm{Pb}}}=8.16\) TeV and \(\sqrt{S_{\mathrm{Pb Pb}}}=5.02\) TeV are $$\begin{aligned} \sigma _{gg\rightarrow (cc)_\mathbf{6}[^1S_0]}^{p\mathrm Pb}= & {} 9.20^{+3.70}_{-2.62}\;\mathrm{\mu b}, \end{aligned}$$ (29) $$\begin{aligned} \sigma _{gg\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1]}^{p\mathrm Pb}= & {} \left( 4.79^{+1.92}_{-1.33}\right) \times 10^1\;\mathrm{\mu b}, \end{aligned}$$ (30) $$\begin{aligned} \sigma _{gc\rightarrow (cc)_\mathbf{6}[^1S_0]}^{p\mathrm Pb}= & {} \left( 1.17^{+0.310}_{-0.230}\right) \times 10^1\;\mathrm{\mu b}, \end{aligned}$$ (31) $$\begin{aligned} \sigma _{gc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1]}^{p\mathrm Pb}= & {} \left( 9.31^{+2.49}_{-1.87}\right) \times 10^1\;\mathrm{\mu b}, \end{aligned}$$ (32) $$\begin{aligned} \sigma _{cc\rightarrow (cc)_\mathbf{6}[^1S_0]}^{p\mathrm Pb}= & {} \left( 1.05^{+0.550}_{-0.226}\right) \times 10^{-2}\;\mathrm{\mu b}, \end{aligned}$$ (33) $$\begin{aligned} \sigma _{cc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1]}^{p\mathrm Pb}= & {} \left( 2.94^{+1.45}_{-0.66}\right) \times 10^{-1}\;\mathrm{\mu b}, \end{aligned}$$ (34) $$\begin{aligned} \sigma _{gg\rightarrow (cc)_\mathbf{6}[^1S_0]}^\mathrm{PbPb}= & {} \left( 1.04^{+0.420}_{-0.288}\right) \times 10^{3}\;\mathrm{\mu b}, \end{aligned}$$ (35) $$\begin{aligned} \sigma _{gg\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1]}^\mathrm{PbPb}= & {} \left( 5.38^{+2.23}_{-1.46}\right) \times 10^{3}\;\mathrm{\mu b}, \end{aligned}$$ (36) $$\begin{aligned} \sigma _{gc\rightarrow (cc)_\mathbf{6}[^1S_0]}^\mathrm{PbPb}= & {} \left( 1.19^{+0.310}_{-0.238}\right) \times 10^{3}\;\mathrm{\mu b}, \end{aligned}$$ (37) $$\begin{aligned} \sigma _{gc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1]}^\mathrm{PbPb}= & {} \left( 1.09^{+0.290}_{-0.223}\right) \times 10^{4}\;\mathrm{\mu b}, \end{aligned}$$ (38) $$\begin{aligned} \sigma _{cc\rightarrow (cc)_\mathbf{6}[^1S_0]}^\mathrm{PbPb}= & {} 1.25^{+0.650}_{-0.269}\;\mathrm{\mu b}, \end{aligned}$$ (39) $$\begin{aligned} \sigma _{cc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1]}^\mathrm{PbPb}= & {} \left( 3.33^{+1.64}_{-0.76}\right) \times 10^1\;\mathrm{\mu b}. \end{aligned}$$ (40) Table 7 Total cross sections (in unit: \(\mu \)b) of \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms for three typical c-quark masses via p-Au and Au–Au collisions at the RHIC   p-Au (0.2 TeV) Au–Au (0.2 TeV) \(m_c\) 1.65 GeV 1.75 GeV 1.85 GeV 1.65 GeV 1.75 GeV 1.85 GeV \(\sigma (gg\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(1.72\times 10^{-1}\) \(1.05\times 10^{-1}\) \(6.53\times 10^{-2}\) \(3.76\times 10^{1}\) \(2.30\times 10^{1}\) \(1.44\times 10^{1}\) \(\sigma (gg\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(9.52\times 10^{-1}\) \(5.79\times 10^{-1}\) \(3.60\times 10^{-1}\) \(2.07\times 10^{2}\) \(1.27\times 10^{2}\) \(7.91\times 10^{1}\) \(\sigma (gc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(4.07\times 10^{-1}\) \(2.83\times 10^{-1}\) \(1.98\times 10^{-1}\) \(8.70\times 10^{1}\) \(6.26\times 10^{1}\) \(4.58\times 10^{1}\) \(\sigma (gc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   3.70 2.58 1.80 \(7.89\times 10^{2}\) \(5.73\times 10^{2}\) \(4.15\times 10^{2}\) \(\sigma (cc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(7.63\times 10^{-5}\) \(4.06\times 10^{-5}\) \(2.88\times 10^{-5}\) \(1.72\times 10^{-2}\) \(9.00\times 10^{-3}\) \(6.35\times 10^{-3}\) \(\sigma (cc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(1.86\times 10^{-3}\) \(1.01\times 10^{-3}\) \(7.01\times 10^{-4}\) \(4.15\times 10^{-1}\) \(2.23\times 10^{-1}\) \(1.54\times 10^{-1}\) Table 8 Total cross sections (in unit: \(\mu \)b) of \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms for three typical c-quark masses via p-Pb and Pb–Pb collisions at the LHC   p-Pb (8.16 TeV) Pb–Pb (5.02 TeV) \(m_c\) 1.65 GeV 1.75 GeV 1.85 GeV 1.65 GeV 1.75 GeV 1.85 GeV \(\sigma (gg\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(1.29\times 10^{1}\) 9.20 6.58 \(1.46\times 10^{3}\) \(1.04\times 10^{3}\) \(7.52\times 10^{2}\) \(\sigma (gg\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(6.71\times 10^{1}\) \(4.79\times 10^{1}\) \(3.46\times 10^{1}\) \(7.61\times 10^{3}\) \(5.38\times 10^{3}\) \(3.92\times 10^{3}\) \(\sigma (gc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(1.48\times 10^{1}\) \(1.17\times 10^{1}\) 9.40 \(1.50\times 10^{3}\) \(1.19\times 10^{3}\) \(9.52\times 10^{2}\) \(\sigma (gc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(1.18\times 10^{2}\) \(9.31\times 10^{1}\) \(7.44\times 10^{1}\) \(1.38\times 10^{4}\) \(1.09\times 10^{4}\) \(8.67\times 10^{3}\) \(\sigma (cc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(1.60\times 10^{-2}\) \(1.05\times 10^{-2}\) \(8.24\times 10^{-3}\) 1.90 1.25 \(9.81\times 10^{-1}\) \(\sigma (cc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(4.39\times 10^{-1}\) \(2.94\times 10^{-1}\) \(2.28\times 10^{-1}\) \(4.96\times 10^{1}\) \(3.32\times 10^{1}\) \(2.56\times 10^{1}\) Table 9 Total cross sections (in unit: \(\mu \)b) of \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms for three typical renormalization scale \(\mu _R\) via p-Au and Au–Au collisions at the RHIC   p-Au (0.2 TeV) Au–Au (0.2 TeV) \(\mu _R\)  \(\sqrt{\hat{s}}\) \(\sqrt{\hat{s}}/2\) \(M_t\) \(\sqrt{\hat{s}}\) \(\sqrt{\hat{s}}/2\) \(M_t\) \(\sigma (gg\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(2.71\times 10^{-2}\) \(5.98\times 10^{-2}\) \(1.05\times 10^{-1}\) 5.85 \(1.30\times 10^{1}\) \(2.30\times 10^{1}\) \(\sigma (gg\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(1.42\times 10^{-1}\) \(3.15\times 10^{-1}\) \(5.79\times 10^{-1}\) \(3.09\times 10^{1}\) \(6.85\times 10^{1}\) \(1.27\times 10^{2}\) \(\sigma (gc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(2.09\times 10^{-1}\) \(2.86\times 10^{-1}\) \(2.83\times 10^{-1}\) \(4.21\times 10^{1}\) \(6.11\times 10^{1}\) \(6.26\times 10^{1}\) \(\sigma (gc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   1.89 2.58 2.58 \(3.83\times 10^{2}\) \(5.56\times 10^{2}\) \(5.73\times 10^{2}\) \(\sigma (cc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(3.11\times 10^{-5}\) \(4.00\times 10^{-5}\) \(4.06\times 10^{-5}\) \(6.81\times 10^{-3}\) \(8.91\times 10^{-3}\) \(9.00\times 10^{-3}\) \(\sigma (cc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(7.70\times 10^{-4}\) \(9.98\times 10^{-4}\) \(1.01\times 10^{-3}\) \(1.68\times 10^{-1}\) \(2.21\times 10^{-1}\) \(2.23\times 10^{-1}\) Table 10 Total cross sections (in unit: \(\mu \)b) of \(\Xi _{cc}\) production via \(g+g\), \(g+c\), and \(c+c\) mechanisms for three typical renormalization scale \(\mu _R\) via p-Pb and Pb–Pb collisions at the LHC   p-Pb (8.16 TeV) Pb–Pb (5.02 TeV) \(\mu _R\)  \(\sqrt{\hat{s}}\) \(\sqrt{\hat{s}}/2\) \(M_t\) \(\sqrt{\hat{s}}\) \(\sqrt{\hat{s}}/2\) \(M_t\) \(\sigma (gg\rightarrow (cc)_\mathbf{6}[^1S_0])\)   6.07 7.50 9.20 \(6.70\times 10^{2}\) \(8.37\times 10^{2}\) \(1.04\times 10^{3}\) \(\sigma (gg\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(3.17\times 10^{1}\) \(3.91\times 10^{1}\) \(4.79\times 10^{1}\) \(3.48\times 10^{3}\) \(4.38\times 10^{3}\) \(5.38\times 10^{3}\) \(\sigma (gc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(1.12\times 10^{1}\) \(1.10\times 10^{1}\) \(1.17\times 10^{1}\) \(1.16\times 10^{3}\) \(1.11\times 10^{3}\) \(1.19\times 10^{3}\) \(\sigma (gc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(9.18\times 10^{1}\) \(8.76\times 10^{1}\) \(9.31\times 10^{1}\) \(1.06\times 10^{4}\) \(1.02\times 10^{4}\) \(1.09\times 10^{4}\) \(\sigma (cc\rightarrow (cc)_\mathbf{6}[^1S_0])\)   \(1.11\times 10^{-2}\) \(1.05\times 10^{-2}\) \(1.05\times 10^{-2}\) 1.31 1.25 1.25 \(\sigma (cc\rightarrow (cc)_{\bar{\mathbf{3}}}[^3S_1])\)   \(3.13\times 10^{-1}\) \(2.96\times 10^{-1}\) \(2.94\times 10^{-1}\) \(3.48\times 10^{1}\) \(3.33\times 10^{1}\) \(3.33\times 10^{1}\) Table 11 Total cross sections (in unit \(\mu \)b) for doubly heavy baryons \(\Xi _{bc}\) and \(\Xi _{bb}\) via the dominant gluon-gluon fusion mechanism via p–N, and N–N collisions at the RHIC and the LHC   RHIC LHC \(\sqrt{S_\mathrm{NN}}\) (TeV) p-Au (0.2) Au–Au (0.2) p-Pb (8.16) Pb–Pb (5.02) \(\sigma (gg\rightarrow (bc)_\mathbf{6}[^1S_0])\) \(1.53\times 10^{-3}\) \(3.27\times 10^{-1}\) \(7.16\times 10^{-1}\) \(7.97\times 10^{1}\) \(\sigma (gg\rightarrow (bc)_{\bar{\mathbf{3}}}[^3S_1])\) \(5.98\times 10^{-3}\) 1.28 2.92 \(3.28\times 10^{2}\) \(\sigma (gg\rightarrow (bc)_\mathbf{6}[^3S_1])\) \(1.06\times 10^{-2}\) 2.28 4.32 \(4.83\times 10^{2}\) \(\sigma (gg\rightarrow (bc)_{\bar{\mathbf{3}}}[^1S_0])\) \(1.65\times 10^{-3}\) \(3.52\times 10^{-1}\) \(7.70\times 10^{-1}\) \(8.51\times 10^{1}\) \(\sigma ^{tot}(\Xi _{bc})\) \(1.98\times 10^{-2}\) 4.24 8.73 \(9.75\times 10^{2}\) \(\sigma (gg\rightarrow (bb)_\mathbf{6}[^1S_0])\) \(1.94\times 10^{-5}\) \(3.71\times 10^{-3}\) \(4.07\times 10^{-2}\) 4.40 \(\sigma (gg\rightarrow (bb)_{\bar{\mathbf{3}}}[^3S_1])\) \(1.05\times 10^{-4}\) \(2.01\times 10^{-2}\) \(2.15\times 10^{-1}\) \(2.35\times 10^{1}\) \(\sigma ^{tot}(\Xi _{bb})\) \(1.25\times 10^{-4}\) \(2.38\times 10^{-2}\) \(2.55\times 10^{-1}\) \(2.79\times 10^{1}\) The renormalization scale-setting problem is an important problem of fixed-order pQCD predictions [67]. As an quantitative estimation of renormalization scale dependence, we choose three usually adopted values as the renormalization scale, i.e. the transverse mass \(M_t\) of \(\Xi _{cc}\), \(\sqrt{\hat{s}}\) and \(\sqrt{\hat{s}}/2\), where \(\sqrt{\hat{s}}\) is the center-of-mass energy of the subprocess. Numerical results are presented in Tables 9 and 10. The scale uncertainties at the RHIC are large, which varies from \(23\%\) to \(76\%\) for various mechanisms via p-Au and Au-Au collisions, accordingly; while the scale uncertainty at the LHC is smaller, which varies from \(4\%\) to \(36\%\) for various mechanisms via p-Pb and Pb-Pb collisions. Thus we need a next-to-leading order calculation to achieve more accurate predictions, especially, by applying the principle of maximum conformality scale-setting approach [68, 69], the renormalization scale uncertainties can be eliminated. 3.4 A simple discussion of \(\Xi _{bc}\) and \(\Xi _{bb}\) production at the RHIC and LHC In this subsection, we present a simple discussion of \(\Xi _{bc}\) and \(\Xi _{bb}\) production properties via p–N and N–N collisions at the RHIC and LHC. Their production mechanisms can be treated via the same way as those of \(\Xi _{cc}\) production, and we adopt the generator GENXICC to do the calculation. As for the input parameters, we take: \(|\Psi _{bc}(0)|^2=0.065\) GeV\(^3\) and \(|\Psi _{bb}(0)|^2 = 0.152\) GeV\(^3\) [19], and \(M_{\Xi _{bc}}=6.9\) with \(m_c=1.8\) GeV and \(m_b=5.1\) GeV, \(M_{\Xi _{bb}}=10.2\) GeV with \(m_b=M_{\Xi _{bb}}/2\). And we set the renormalization scale as \(M_t\). In different to the \(\Xi _{cc}\) production, for the present case of \(\Xi _{bc}\) and \(\Xi _{bb}\), the extrinsic mechanisms shall be highly suppressed by the much smaller bottom-quark PDF, thus we shall only consider the dominant gluon-gluon fusion mechanism. We present the total cross sections for \(\Xi _{bc}\) and \(\Xi _{bb}\) produced via p–N and N–N collisions at the RHIC and the LHC in Table 11. By summing up different spin-color diquark configurations, we obtain $$\begin{aligned}&\left. \sigma ^\mathrm{tot}_{p\mathrm Au}(\Xi _{bc})\right| _\mathrm{RHIC} = 1.98\times 10^{-2}\;\mathrm{\mu b}, \end{aligned}$$ (41) $$\begin{aligned}&\left. \sigma ^\mathrm{tot}_{p\mathrm Au}(\Xi _{bb})\right| _\mathrm{RHIC} = 1.25\times 10^{-4}\;\mathrm{\mu b}, \end{aligned}$$ (42) $$\begin{aligned}&\left. \sigma ^\mathrm{tot}_\mathrm{AuAu}(\Xi _{bc})\right| _\mathrm{RHIC} = 4.24\;\mathrm{\mu b}, \end{aligned}$$ (43) $$\begin{aligned}&\left. \sigma ^\mathrm{tot}_\mathrm{AuAu}(\Xi _{bb})\right| _\mathrm{RHIC} = 2.38\times 10^{-2}\;\mathrm{\mu b}, \end{aligned}$$ (44) $$\begin{aligned}&\left. \sigma ^\mathrm{tot}_{p\mathrm Pb}(\Xi _{bc})\right| _\mathrm{LHC} = 8.73\;\mathrm{\mu b}, \end{aligned}$$ (45) $$\begin{aligned}&\left. \sigma ^\mathrm{tot}_{p\mathrm Pb}(\Xi _{bb})\right| _\mathrm{LHC} = 2.55\times 10^{-1}\;\mathrm{\mu b}, \end{aligned}$$ (46) $$\begin{aligned}&\left. \sigma ^\mathrm{tot}_\mathrm{PbPb}(\Xi _{bc})\right| _\mathrm{LHC} = 9.75\times 10^{2}\;\mathrm{\mu b}, \end{aligned}$$ (47) $$\begin{aligned}&\left. \sigma ^\mathrm{tot}_\mathrm{PbPb}(\Xi _{bb})\right| _\mathrm{LHC} = 2.79\times 10^{1}\;\mathrm{\mu b}. \end{aligned}$$ (48) To estimate the event numbers for \(\Xi _{bc}\) and \(\Xi _{bb}\) porduction, we adopt the same luminosities, as mentioned above, for the p–N and N–N collisions at the RHIC and LHC. Our results show that at the RHIC, \(8.9 \times 10^4\) \(\Xi _{bc}\) and \(5.6 \times 10^2\) \(\Xi _{bb}\) events can be generated via p-Au collision, \(3.4 \times 10^5\) \(\Xi _{bc}\) and \(1.9 \times 10^3\) \(\Xi _{bb}\) events via Au-Au collision; At the LHC, \(4.4 \times 10^6\) \(\Xi _{bc}\) and \(1.3 \times 10^5\) \(\Xi _{bb}\) events can be generated via p-Pb collision, \(3.5 \times 10^6\) \(\Xi _{bc}\) and \(1.0 \times 10^5\) \(\Xi _{bb}\) events can be generated via Pb-Pb collision. The number of \(\Xi _{bc}\) events to be generated at the RHIC or LHC are smaller than that of the \(\Xi _{cc}\) events under the same collision by about one order; and the number of \(\Xi _{bb}\) events to be generated at the RHIC or LHC are smaller than that of the \(\Xi _{cc}\) events under the same collision by about two order. Those results show that if more experimental data have been accumulated at the RHIC and LHC, one may also have the chance to study the properties of the other two doubly heavy baryons, \(\Xi _{bc}\) and \(\Xi _{bb}\). 4 Summary We have studied the \(\Xi _{cc}\) production via p–N and N–N collisions at the RHIC and LHC. The generator GENXICC with suitable changes for the use of nuclear PDF has been adopted for the calculation. Our results show that in addition to the gluon-gluon fusion mechanism, the extrinsic charm mechanisms via \(g+c\) and \(c+c\) subprocesses, are important to achieve a sound prediction for the \(\Xi _{cc}\) production. By summing up contributions from \(g+g\), \(g+c\) and \(c+c\) mechanisms and contributions from different spin-and-color configurations of the intermediate (cc)-diquark together, we observe that sizable number of \(\Xi _{cc}\) events can be produced via p–N and N–N collisions at the RHIC and LHC. More explicitly, we have shown that \(1.6 \times 10^6\) and \(6.3 \times 10^6\) \(\Xi _{cc}\) events can be produced via p-Au and Au-Au collisions at the RHIC, respectively; \(8.1 \times 10^7\) and \(6.7 \times 10^7\) \(\Xi _{cc}\) events can be produced via p-Pb and Pb-Pu collisions at the LHC, respectively. Sizable number of \(\Xi _{cc}\) events can be accumulated via p–N and N–N collisions at the RHIC and LHC. Thus, in addition to the pp collision as now been performed by LHCb experiment, the p–N and N–N collisions at the hadron colliders can be a good platform for investigating the properties of the doubly heavy baryons. As shown by Figs. 1, 2, the \(g+c\) mechanism is larger than \(g+g\) mechanism in small \(p_t\) region via p–N and N–N collisions. If the experimental measurements can be extended to small \(p_t\) region, then one may study extrinsic charm mechanism by carefully measuring the \(\Xi _{cc}\) events. Footnotes 1. The general-mass variable-flavor-number scheme [42, 43, 44, 45] has been adopted to deal with the double coupling problem among different channels and the QCD factorization is done within the framework of nonrelativistic QCD (NRQCD) [46]. 2. By using a concrete fragmentation function, as suggested by Ref.[50], to deal with the production shall give quite small difference to this simple assumption [31]. 3. In predicting the event numbers, at the RHIC, one operation year means \(10^7\)s for p-Au and Au-Au collisions; At the LHC, one operation year means \(10^7\)s for pp-collision and \(10^6\)s for p-Pb and Pb-Pb collisions [65, 66]. Notes Acknowledgements This work was supported in part by the Natural Science Foundation of China under Grant No.11605029, No.11675239, No.11535002 and No.11625520, the science project of colleges and universities directly under the Guangzhou Education Bureau No.1201630158, the Foundation for Fostering the Scientific and Technical Innovation of Guangzhou University, and by the Fundamental Research Funds for the Central Universities under the Grant No.2018CDPTCG0001/3. References 1. M. Gell-Mann, A schematic model of Baryons and Mesons. Phys. Lett. 8, 214 (1964)ADSCrossRefGoogle Scholar 2. G. Zweig, An SU(3) model for strong interaction symmetry and its breaking. Version 1, CERN-TH-401Google Scholar 3. D. Ebert, R.N. Faustov, V.O. Galkin, A.P. Martynenko, V.A. Saleev, Heavy baryons in the relativistic quark model. Z. Phys. C 76, 111 (1997)CrossRefGoogle Scholar 4. S.M. Gerasyuta, D.V. Ivanov, Charmed baryons in bootstrap quark model. Nuovo Cim. A 112, 261 (1999)ADSCrossRefGoogle Scholar 5. C. Itoh, T. Minamikawa, K. Miura, T. Watanabe, Doubly charmed baryon masses and quark wave functions in baryons. Phys. Rev. D 61, 057502 (2000)ADSCrossRefGoogle Scholar 6. M. Mattson, SELEX Collaboration, et al., First observation of the doubly charmed baryon \(\Xi _{cc}^{+}\). Phys. Rev. Lett. 89, 112001 (2002)Google Scholar 7. A. Ocherashvili, SELEX Collaboration, et al., Confirmation of the double charm baryon \(\Xi _{cc}^{+}\)(3520) via its decay to \(p D^+ K^-\). Phys. Lett. B 628, 18 (2005)Google Scholar 8. S.P. Ratti, New results on c-baryons and a search for cc-baryons in FOCUS. Nucl. Phys. Proc. Suppl. 115, 33 (2003)ADSCrossRefGoogle Scholar 9. B. Aubert, BaBar Collaboration, et al., Search for doubly charmed baryons \(\Xi _{cc}^{+}\) and \(\Xi _{cc}^{++}\) in BABAR. Phys. Rev. D 74, 011103 (2006)Google Scholar 10. R. Chistov, Belle Collaboration, et al., Observation of new states decaying into \(\Lambda _c^+ K^- \pi ^+\) and \(\Lambda _c^+ K^0_S \pi ^-\). Phys. Rev. Lett. 97, 162001 (2006)Google Scholar 11. Y. Kato, Belle Collaboration, et al., Search for doubly charmed baryons and study of charmed strange baryons at Belle. Phys. Rev. D 89, 052003 (2014)Google Scholar 12. R. Aaij, LHCb Collaboration, et al., Search for the doubly charmed baryon \(\Xi _{cc}^+\). J. High Energy Phys. 12, 090 (2013)Google Scholar 13. R. Aaij, LHCb Collaboration, et al., Observation of the doubly charmed baryon \(\Xi _{cc}^{++}\). Phys. Rev. Lett. 119, 112001 (2017)Google Scholar 14. C.H. Chang, J.X. Wang, X.G. Wu, GENXICC: A generator for hadronic production of the double heavy baryons \(\Xi _{cc}\), \(\Xi _{bc}\) and \(\Xi _{bb}\). Comput. Phys. Commun. 177, 467 (2007)ADSCrossRefGoogle Scholar 15. C.H. Chang, J.X. Wang, X.G. Wu, GENXICC2.0: An upgraded version of the generator for hadronic production of double heavy baryons \(\Xi _{cc}\), \(\Xi _{bc}\) and \(\Xi _{bb}\). Comput. Phys. Commun. 181, 1144 (2010)ADSCrossRefGoogle Scholar 16. X.Y. Wang, X.G. Wu, GENXICC2.1: An improved version of GENXICC for hadronic production of doubly heavy baryons. Comput. Phys. Commun. 184, 1070 (2013)ADSCrossRefGoogle Scholar 17. A.F. Falk, M.E. Luke, M.J. Savage, M.B. Wise, Heavy quark fragmentation to baryons containing two heavy quarks. Phys. Rev. D 49, 555 (1994)ADSCrossRefGoogle Scholar 18. V.V. Kiselev, A.K. Likhoded, M.V. Shevlyagin, Double charmed baryon production at B factory. Phys. Lett. B 332, 411 (1994)ADSCrossRefGoogle Scholar 19. S.P. Baranov, On the production of doubly flavored baryons in \(pp\), \(ep\), and \(\gamma \gamma \) collisions. Phys. Rev. D 54, 3228 (1996)ADSCrossRefGoogle Scholar 20. A.V. Berezhnoy, V.V. Kiselev, A.K. Likhoded, A.I. Onishchenko, Doubly charmed baryon production in hadronic experiments. Phys. Rev. D 57, 4385 (1998)ADSCrossRefGoogle Scholar 21. D.A. Gunter, V.A. Saleev, Hadronic production of doubly charmed baryons via charm excitation in proton. Phys. Rev. D 64, 034006 (2001)ADSCrossRefGoogle Scholar 22. V.V. Braguta, V.V. Kiselev, A.E. Chalov, Pair production of doubly heavy diquarks. Phys. Atom. Nucl. 65, 1537 (2002)ADSCrossRefGoogle Scholar 23. J.P. Ma, Z.G. Si, Factorization approach for inclusive production of doubly heavy baryon. Phys. Lett. B 568, 135 (2003)ADSCrossRefGoogle Scholar 24. E. Braaten, M. Kusunoki, Y. Jia, T. Mehen, \(\Lambda _{c}^{+}/\Lambda _{c}^{-}\) asymmetry in hadroproduction from heavy quark recombination. Phys. Rev. D 70, 054021 (2004)ADSCrossRefGoogle Scholar 25. S.Y. Li, Z.G. Si, Z.J. Yang, Doubly heavy baryon production at gamma gamma collider. Phys. Lett. B 648, 284 (2007)ADSCrossRefGoogle Scholar 26. Z.J. Yang, T. Yao, Doubly heavy baryon production at polarized photon collider. Chin. Phys. Lett. 24, 3378 (2007)ADSCrossRefGoogle Scholar 27. J.W. Zhang, X.G. Wu, T. Zhong, Y. Yu, Z.Y. Fang, Hadronic production of the doubly heavy baryon \(\Xi _{bc}\) at LHC. Phys. Rev. D 83, 034026 (2011)ADSCrossRefGoogle Scholar 28. J. Jiang, X.G. Wu, Q.L. Liao, X.C. Zheng, Z.Y. Fang, Doubly heavy baryon production at a high luminosity \(e^+ e^-\) collider. Phys. Rev. D 86, 054021 (2012)ADSCrossRefGoogle Scholar 29. J. Jiang, X.G. Wu, S.M. Wang, J.W. Zhang, Z.Y. Fang, A further study on the doubly heavy baryon production around the \(Z^0\) peak at a high luminosity \(e^+ e^-\) collider. Phys. Rev. D 87, 054027 (2013)ADSCrossRefGoogle Scholar 30. A.P. Martynenko, A.M. Trunin, Relativistic corrections to the pair double heavy diquark production in \(e^+e^- \rightarrow \) annihilation. Phys. Rev. D 89, 014004 (2014)ADSCrossRefGoogle Scholar 31. G. Chen, X.G. Wu, Z. Sun, Y. Ma, H.B. Fu, Photoproduction of doubly heavy baryon at the ILC. J. High Energy Phys. 12, 018 (2014)ADSCrossRefGoogle Scholar 32. Z.J. Yang, X.X. Zhao, The production of \(\Xi _{bb}\) at photon collider. Chin. Phys. Lett. 31, 091301 (2014)ADSCrossRefGoogle Scholar 33. Z.J. Yang, P.F. Zhang, Y.J. Zheng, Doubly heavy baryon production in \(e^{+}e^{-}\) annihilation. Chin. Phys. Lett. 31, 051301 (2014)ADSCrossRefGoogle Scholar 34. A.P. Martynenko, A.M. Trunin, Pair double heavy diquark production in high energy proton-proton collisions. Eur. Phys. J. C 75, 138 (2015)ADSCrossRefGoogle Scholar 35. W.K. Lai, A.K. Leibovich, \(\Lambda _{c}^{+}/\Lambda _{c}^{-}\) and \(\Lambda _{b}^{0}/\bar{\Lambda }_{b}^0\) production asymmetry at the LHC from heavy quark recombination. Phys. Rev. D 91, 054022 (2015)ADSCrossRefGoogle Scholar 36. S.J. Brodsky, S. Groote, S. Koshkarev, Resolving the SELEXCLHCb double-charm baryon conflict: the impact of intrinsic heavy-quark hadroproduction and supersymmetric light-front holographic QCD. Eur. Phys. J. C 78, 483 (2018)ADSCrossRefGoogle Scholar 37. X. Yao, B. Muller, Doubly charmed baryon production in heavy ion collisions. Phys. Rev. D 97, 074003 (2018)ADSCrossRefGoogle Scholar 38. C.H. Chang, C.F. Qiao, J.X. Wang, X.G. Wu, Estimate of the hadronic production of the doubly charmed baryon \(\Xi _{cc}\) under GM-VFN scheme. Phys. Rev. D 73, 094022 (2006)ADSCrossRefGoogle Scholar 39. C.H. Chang, J.P. Ma, C.F. Qiao, X.G. Wu, Hadronic production of the doubly charmed baryon \(\Xi _{cc}\) with intrinsic charm. J. Phys. G 34, 845 (2007)CrossRefGoogle Scholar 40. G. Chen, X.G. Wu, J.W. Zhang, H.Y. Han, H.B. Fu, Hadronic production of \(\Xi _{cc}\) at a fixed-target experiment at the LHC. Phys. Rev. D 89, 074020 (2014)ADSCrossRefGoogle Scholar 41. X.C. Zheng, C.H. Chang, Z. Pan, Production of doubly heavy-flavored hadrons at \(e^+e^-\) colliders. Phys. Rev. D 93, 034019 (2016)ADSCrossRefGoogle Scholar 42. M.A.G. Aivazis, F.I. Olness, W.K. Tung, Leptoproduction of heavy quarks. 1. General formalism and kinematics of charged current and neutral current production processes. Phys. Rev. D 50, 3085 (1994)ADSCrossRefGoogle Scholar 43. M.A.G. Aivazis, J.C. Collins, F.I. Olness, W.K. Tung, Leptoproduction of heavy quarks. 2. A unified QCD formulation of charged and neutral current processes from fixed target to collider energies. Phys. Rev. D 50, 3102 (1994)ADSCrossRefGoogle Scholar 44. F.I. Olness, R.J. Scalise, W.K. Tung, Heavy quark hadroproduction in perturbative QCD. Phys. Rev. D 59, 014506 (1999)ADSCrossRefGoogle Scholar 45. J. Amundson, C. Schmidt, W.K. Tung, X. Wang, Charm production in deep inelastic scattering from threshold to high \(Q^{2}\). J. High Energy Phys. 10, 031 (2000)ADSCrossRefGoogle Scholar 46. G.T. Bodwin, E. Braaten, G.P. Lepage, Rigorous QCD analysis of inclusive annihilation and production of heavy quarkonium. Phys. Rev. D 51, 1125 (1995)ADSCrossRefGoogle Scholar 47. N.B. Chang et al., Physics Perspectives of Heavy-Ion Collisions at Very High Energy. Sci. China Phys. Mech. Astron. 59, 621001 (2016)CrossRefGoogle Scholar 48. G. Aarts et al., Heavy-flavor production and medium properties in high-energy nuclear collisions - What next? Eur. Phys. J. A 53, 93 (2017)ADSCrossRefGoogle Scholar 49. G. Chen, C.H. Chang, X.G. Wu, \(B_c (B_c^*)\) meson production via the proton-nucleus and the nucleus-nucleus collision modes at the colliders RHIC and LHC. Phys. Rev. D 97, 114022 (2018)ADSCrossRefGoogle Scholar 50. C. Peterson, D. Schlatter, I. Schmitt, P.M. Zerwas, Scaling Violations in Inclusive e+ e- Annihilation Spectra. Phys. Rev. D 27, 105 (1983)ADSCrossRefGoogle Scholar 51. T. Sjostrand, S. Mrenna, P.Z. Skands, PYTHIA 6.4 physics and manual. J. High Energy Phys. 05, 026 (2006)ADSCrossRefGoogle Scholar 52. X.N. Wang, M. Gyulassy, HIJING: A Monte Carlo model for multiple jet production in \(pp\), \(p\)A and AA collisions. Phys. Rev. D 44, 3501 (1991)ADSCrossRefGoogle Scholar 53. Z.W. Lin, C.M. Ko, B.A. Li, B. Zhang, S. Pal, A Multi-phase transport model for relativistic heavy ion collisions. Phys. Rev. C 72, 064901 (2005)ADSCrossRefGoogle Scholar 54. B. Alver, M. Baker, C. Loizides, P. Steinberg, The PHOBOS Glauber Monte Carlo. arXiv:0805.4411 [nucl-ex] 55. W. Broniowski, M. Rybczynski, P. Bozek, GLISSANDO: Glauber initial-state simulation and more. Comput. Phys. Commun. 180, 69 (2009)ADSCrossRefGoogle Scholar 56. M. Rybczynski, G. Stefanek, W. Broniowski, P. Bozek, GLISSANDO 2: GLauber Initial-State Simulation AND mOre ver. 2. Comput. Phys. Commun. 185, 1759 (2014)ADSCrossRefGoogle Scholar 57. C. Loizides, J. Kamin, D. d’Enterria, Improved Monte Carlo Glauber predictions at present and future nuclear colliders. Phys. Rev. C 97(5), 054910 (2018)ADSCrossRefGoogle Scholar 58. K. Kovarik et al., nCTEQ15 - Global analysis of nuclear parton distributions with uncertainties in the CTEQ framework. Phys. Rev. D 93, 085037 (2016)ADSCrossRefGoogle Scholar 59. J.J. Aubert, European Muon Collaboration, et al., Measurements of the nucleon structure functions \(F2_n\) in deep inelastic muon scattering from deuterium and comparison with those from hydrogen and iron. Nucl. Phys. B 293, 740 (1987)Google Scholar 60. S. Gavin, R. Vogt, \(J/\psi \) Suppression From Hadron - Nucleus to Nucleus-nucleus Collisions. Nucl. Phys. B 345, 104 (1990)ADSCrossRefGoogle Scholar 61. C.W. De Jager, H. De Vries, C. De Vries, Nuclear charge and magnetization density distribution parameters from elastic electron scattering. Atom. Data Nucl. Data Tabl. 14, 479 (1974)ADSCrossRefGoogle Scholar 62. R. Vogt, Cold Nuclear Matter Effects on \(J/\psi \) and \(\Upsilon \) Production at the LHC. Phys. Rev. C 81, 044903 (2010)ADSCrossRefGoogle Scholar 63. N. Brambilla, Quarkonium Working Group, et al., Heavy quarkonium physics. arXiv:hep-ph/0412158 64. M. Tanabashi et al., Particle Data Group. Phys. Rev. D 98, 030001 (2018)ADSCrossRefGoogle Scholar 65. F. Carminati, ALICE Collaboration, et al., ALICE: Physics performance report, volume I. J. Phys. G 30, 1517 (2004)Google Scholar 66. B. Alessandro, ALICE Collaboration, et al., ALICE: Physics performance report, volume II. J. Phys. G 32, 1295 (2006)Google Scholar 67. X.G. Wu, S.J. Brodsky, M. Mojaza, The Renormalization Scale-Setting Problem in QCD. Prog. Part. Nucl. Phys. 72, 44 (2013)ADSCrossRefGoogle Scholar 68. S.J. Brodsky, X.G. Wu, Eliminating the Renormalization Scale Ambiguity for Top-Pair Production Using the Principle of Maximum Conformality. Phys. Rev. Lett. 109, 042002 (2012)ADSCrossRefGoogle Scholar 69. M. Mojaza, S.J. Brodsky, X.G. Wu, Systematic All-Orders Method to Eliminate Renormalization-Scale and Scheme Ambiguities in Perturbative QCD. Phys. Rev. Lett. 110, 192001 (2013)ADSCrossRefGoogle Scholar Copyright information © The Author(s) 2018 Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP3 Authors and Affiliations Gu Chen1Chao-Hsi Chang23Xing-Gang Wu4Email authorView author's OrcID profile1.School of Physics and Electronic EngineeringGuangzhou UniversityGuangzhouPeople’s Republic of China2.Institute of Theoretical PhysicsChinese Academy of SciencesBeijingPeople’s Republic of China3.School of Physical SciencesUniversity of Chinese Academy of SciencesBeijingChina4.Department of PhysicsChongqing UniversityChongqingPeople’s Republic of China


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-018-6283-1.pdf

Gu Chen, Chao-Hsi Chang, Xing-Gang Wu. Hadronic production of the doubly charmed baryon via the proton–nucleus and the nucleus–nucleus collisions at the RHIC and LHC, The European Physical Journal C, 2018, 801, DOI: 10.1140/epjc/s10052-018-6283-1