Multi-charmed and singled charmed hadrons from coalescence: yields and ratios in different collision systems at LHC

We study the production of charmed and multi-charmed hadrons in ultra-relativistic Heavy Ion Collisions coupling the transport approach for charm dynamics in the medium to an hybrid hadronization model of coalescence plus fragmentation. In this paper, we mainly discuss the particle yields for single charmed and multi-charmed baryons focusing mainly on the production of Ξ cc and Ω ccc . We provide first predictions for PbPb collision in 0 − 10% centrality class and then we explore the system size dependence through KrKr , to ArAr and OO collisions, planned within the ALICE3 experiment. In these cases, a monotonic behavior for the yields emerges which can be tested in future experimental data. We found about three order of magnitude increase in the production of Ω ccc in PbPb collisions compared with the yield in small collision systems like OO collisions. Furthermore, we investigate the effects on the Ω ccc particle production and spectra coming from the modification of the charm quark distribution due to the different size of the collision systems comparing also to the case of thermalized charm distributions. These results suggest that observation on the Ω ccc spectra and their evolution across system size can give information about the partial thermalization of the charm quark distribution as well as to its wave function width.


I. INTRODUCTION
Many probes have been proposed to investigate the properties of the matter created in ultra-Relativistic Heavy Ion collisions (uRHICs).A preminent one is the heavy quark hadron production; all the observables related to them have been considered one of the most useful probes to characterize the quark-gluon plasma.Due to their large masses Heavy quarks, namely charm and bottom, are considered as a solid probe to characterize the QGP phase [1][2][3].They are produced by pQCD processes with a formation time τ 0 < 0.08f m/c ≪ τ QGP that permits to probe also the strong electromagnetic and vortical fields expected in the initial stage of the collision [4][5][6].On the other hand, the large mass implies a larger thermalization time w.r.t.light counterpart and appears currently to be comparable to the one of the QGP itself [1,7] and HQs can probe the whole evolution of the plasma.Furthermore, they are expected to conserve memory of the history of the system evolution and the final hadrons can keep information about the out-of-equilibrium initial conditions.Two main observables have been studied in uRHICs for HF hadrons: the heavy mesons nuclear modification factor R AA (p T ) [8][9][10] , and the elliptic flow, v 2 (p T ) [11,12].These observables have been studied in order to extract the heavy flavour transport coefficients and understand the HQs dynamics in QGP from a theoretical point of view [13][14][15][16][17][18][19][20][21][22][23].Recently, further efforts have been done to extend the analysis to higher order anisotropic flows v n [24][25][26][27][28] that can provide more constraints on the extraction of the transport coefficients.The coalescence mechanism is one of the possible description for the hadronization process present in the Quark-Gluon Plasma, it is able to explain the p T baryon to meson spectra and the splitting of elliptic flow of light mesons and baryons produced in heavy ion collisions at top RHIC energies [29][30][31][32][33] as well as at LHC [34].In the last decades coalescence models were extended to include finite width to take into account for off-shell effects [35][36][37] which however preserve the main features of the initial modeling like the enhancement of baryon production and an approximate quark number scaling of the elliptic flow v 2 .In particular, for open heavy flavor, the hadronization by coalescence play an important role to determine the R AA (p T ) and the elliptic flow v 2 (p T ) = ⟨cos(2ϕ)⟩ of D meson affecting the evolution of D s [2,7,18,20,23,38,39]. Furthermore, coalescence approaches have predicted an unexpected large Λ c /D 0 ∼ 0.5 − 1 in AA collisions [40,41], that has been recently observed at RHIC energies [42] and in pp, pA, AA collisions at LHC [43][44][45].Such a large ratio appears to be a strong violation of the universality of the fragmentation function seen in elementary collision systems, where the observed Λ c /D 0 is of O(10 −1 ) [46].A large baryon over meson ratio in the charm sector, as observed in AA collision, is compatible with hadronization by coalescence [40,47,48].Furthermore, very recently, the Catania coalescence model have correctly predicted the ratio of Λ c /D 0 , Ξ c /D 0 and Ω c /D 0 even in pp collisions [49][50][51].

arXiv:2305.03687v3 [hep-ph] 13 Mar 2024
At the same time, in the heavy quark sector, theoretical studies have shown that the investigation on the bottom quarks can provide further insight on the HQ thermalization looking at R AA , anisotropic flows and hadronization mechanism of B mesons and Λ b [24,52,53].In this paper we are concentrating our focus on the production of the multi-charmed hadrons, which have two or three charm as constituent quarks.The study of these probes is, first of all, a natural extension of the investigation done for single-charmed hadrons, moreover it can provide further information about the hadronization mechanism, that should be much more sensitive to the charm quark features in the QGP medium.The production of multicharmed hadrons are part of the physical motivation of the ALICE3 proposal for HI-LHC [54].First observations of multi-charmed baryons was reported in 2002 by SELEX collaboration at Fermilab for Ξ + cc [55].Recently, LHCb collaboration has observed Ξ ++ cc in pp collisions at top LHC energies in two different decay channels [56,57].First theoretical studies on the multi-charmed production were made with statistical hadronization model in [58], and more recently an increasing interest has grown about the topic, with different calculations present in literature [59][60][61][62][63].The multi-charmed hadrons, as well as the exotic states, e.g.T + cc , and pentaquarks, has been also indicated as interesting for future experimental developments and investigations [54].The paper is organized in the following way: in Section II, we describe our hybrid hadronization model by coalescence plus fragmentation.In Section III, we discuss the collision systems characteristics and parton distribution setup for our calculations.In Section IV, we present our results for the singlecharmed and multi-charmed yields in PbPb collision at LHC energies compared with Statistical Hadronization Model (SHM), and the effect on the production coming from microscopical details of our model.In Section V, we study the multi-charm production varying the collision systems, i.e.PbPb , KrKr , ArAr and OO , and the role of non-equilibrium behaviour of charm quark distribution.Finally, in Section VI, we give our conclusions.

II. HYBRID HADRONIZATION BY COALESCENCE AND FRAGMENTATION
Coalescence models have been widely applied as a mechanism of hadronization in HICs.These models were successful for the explanation of the constituent-quark number scaling of the elliptic flow v 2 (p T ) in AA collisions for light hadrons and the large baryon-to-meson ratios for light hadron production at "intermediate" p T [64].In recent years, for the HF hadron chemistry in AA collisions, the charmed hadrons production has been investigated within these models predicting a large Λ c /D 0 [40,41,65,66].This approach is suitable also to provide results in different collision systems, and it was proposed as an alternative approach of hadronization in very small collisions systems like pp collisions, assuming that also in this systems a small QGP droplet can be formed where the charm quarks, produced in perturbative processes, can hadronize via recombination with light thermal partons present in the medium [50].It was shown that this hadronization mechanism can explain in a natural way the p T dependence of different charmed hadrons ratio like Λ c /D 0 , Ξ c /D 0 , Ω c /D 0 at top LHC energies.
In this section, we recall the basic elements of the coalescence model developed in [30][31][32]38] and based on the Wigner formalism.The momentum spectrum of hadrons formed by coalescence of quarks can be written as: with g H we indicate the statistical factor to form a colorless hadron from quarks and antiquarks with spin 1/2.For mesons with spin-0 the statistical factors g = 1/36 gives the probability that two random quarks have the right colour, spin, isospin to match the quantum number of the considered mesons.For baryons with spin-1/2 the statistical factors is g = 1/108.The dσ i denotes an element of a space-like hypersurface, while f qi are the quark (anti-quark) phase-space distribution functions for i-th quark (anti-quark).C H (...) = N f H (...) with N a normalization factor; while f H (x 1 ...x Nq , p 1 ...p Nq ) is the Wigner function which describes the spatial and momentum distribution of quarks in a hadron.Following the Refs.[38,40,41] we adopt for the Wigner distribution function a Gaussian shape in space and momentum, where N q is the number of constituent quarks.Notice that N has been fixed to guarantee that in the limit p → 0 all the charm hadronize by coalescence in a heavy hadron.This is imposed by requiring that the total coalescence probability for charm quarks gives lim p→0 P tot coal = 1.It has been shown, by other studies, that the inclusion of missing charm-baryon states [66] or the variation of the width of the D meson wave function [48,65], can permit that all the zero momentum charm quarks can be converted to charmed hadrons.The relative coordinate are evaluated going into the CM frame of the particles involved in the process and are defined as follows.For mesons the relative coordinates (r 1 , p r1 ) are given by, while for baryons are defined as and r 2 , p r2 are given by r 2 = 2 3 The σ ri are the covariant widths, that can be related to the oscillator frequency ω by σ ri = 1/ √ µ i ω where µ i are the reduced masses In our calculations the masses of light and heavy quarks have been fixed to m u,d = 300 MeV, m s = 380 MeV, m c = 1.5 GeV.The Wigner function for the heavy meson has only one parameter which is the width σ r1 , while for baryons there are two parameters σ r1 and σ r2 that are related by the oscillatory frequency ω through the reduced masses by The widths of the Wigner function f H are fixed by using the relation with the size of the hadron and in particular to the root-mean-square charge-radius of the hadron, with N = 2, 3 for mesons and baryons respectively; see ref. [41,50] for single-charmed mesons and baryons.The mean square charge radius of mesons and single charmed baryons used as reference come from quark model [67,68].The widths for heavy hadron used in this work are shown in Table I.The corresponding mean square charge radii evaluated from these widths have values within the uncertainties coming from the quark models calculation mentioned above.In table II we report the widths and radii for multi-charm hadrons.To fix the widths values of the Wigner function in the case of multi-charmed hadrons (Ξ ++ cc and Ω ccc ) we start from the frequency ω of the single charmed hadrons, Ξ c and Ω c respectively; and then we calculate the widths with a scaling from the frequency and the new reduced masses according to Eq.7.As known from previous works on coalescence [19,34,40,41,69] the coalescence probability decreases at increasing p T , this behaviour let the the standard independent fragmentation to be the dominant hadronization process for the production at high p T .Hence, the inclusion of the hadronization by fragmentation is necessary to describe correctly the transition to the high momentum regime but does not affect significantly the yield.In our approach the smooth transition from low to high p T regime is given by introducing a fragmentation probability P f rag (p T ).As done in Ref. [41,50] we start from the probability that one charm quark can hadronize by coalescence and we assume that charm quarks that do not hadronize via coalescence are converted to hadrons by fragmentation.The fragmentation probability is given by P f rag (p T ) = 1 − P tot coal (p T ), where P tot coal is the total coalescence probability.Notice that with the same approach used in [41,50] the coalescence probability for a charm quark to hadronize via coalescence is forced to be 1 at p T ≈ 0. The hadron momentum spectra from the charm parton fragmentation is given by: is the fragmentation function and z = p had /p c is the momentum fraction of heavy quarks transferred to the final heavy hadron while Q 2 = (p had /2z) 2 is the momentum scale for the fragmentation process.In our calculations for charm quarks we have used the Peterson fragmentation function [70] where ϵ c is a free parameter that is determined assuring that the shape of the fragmentation function agrees with the experimental data on p T distributions for D 0 and Λ c at p T > 10 GeV The ϵ c parameters used here are the same as in [41].We notice that at high p T the fragmentation becomes to be the dominant charm hadronization mechanism and a quark will hadronize according to the different fragmentation fractions into specific final charmed hadron channels.The fragmentation fraction is evaluated according to PYTHIA8 ratios at high p T > 10 GeV that are similar to the e + + e − [46] apart from an increase of the fraction for Λ c and moderate decrease of the fraction going to D 0 , as already done in [41,50].In our calculation the multi-dimensional integrals in the coalescence formula are evaluated by using a Monte-Carlo method, see [29,41] for more details.

III. FIREBALL AND PARTON DISTRIBUTION
In our calculation the bulk of particles that we assume is a thermalized system of gluons and u, d, s quarks and anti-quarks.The partons are distributed uniformly in the transverse plane and in the rapidity range |y z | < 0.5.The longitudinal momentum distribution is assumed to be boost-invariant in the range y ∈ (−0.5, +0.5), and is included a radial flow with the following radial profile R , where R is the transverse radius of the fireball.Partons at low transverse momentum, p T < 2 GeV, are considered thermally distributed where m T = p 2 T + m 2 q,q is the transverse mass.The factors g q = g q = 6 are the spin-color degeneracy.The presence of gluons in the quark-gluon plasma is taken into account by converting them to quarks and anti-quark pairs according to the flavour compositions, as assumed in [31,71].
The volume of the fireball in one unit of rapidity is given by V = πR 2 ⊥ τ where R ⊥ is the radius of the fireball taking into account the radial expansion.We fix the radial flow and the volume by imposing the total multiplicity dN/dy and the total transverse energy dE T /dy to be equal to the experimental data.This constraints lead to values for radial flow, radius and τ as shown in Table III, in quite good agreement also with simulations in kinetic transport approaches, and with Statistical Hadronization Model estimate for the active hadronization Volume [60].Such values correspond in one unity of rapidity to fireball volume for P bP b, KrKr, ArAr and OO that are, respectively, 5000f m 3 , 2000f m 3 , 920f m 3 and 345f m 3 .For the initial p T distribution of partons at high transverse momentum, p T > 2.5GeV , we have considered the mini-jets that have undergone the jet quenching mechanism.As done in Ref. [41], we have considered the initial p T distribution from pQCD calculation and the thickness function of the Glauber model to go from pp collisions to AA ones.The charm pair production is described by hard process and it is described by perturbative QCD (pQCD) at NNLO.Therefore, the starting point to compute the initial heavy quarks spectra at LHC collision energy of √ s = 5.02 T eV is by pQCD calculation.In our calculation the initial charm quark spectrum have been taken in accordance to the Fixed Order + Next-to-Leading Log (FONLL), as given in Refs.[72,73].The number of charm quarks has been chosen in accordance with a charm quark cross section of dσ cc /dy ∼ 0.500 that scales from P bP b to ArAr, KrKr and OO with the mean nuclear thickness function T AA .[74].This cross section and T AA scaling gives, for these calculations, a number of charmed quarks that are N Finally, the charm quark distribution evolution is obtained starting from the perturbative distribution (FONLL) and solving the relativistic Boltzmann transport equations for charm quarks scattering in a bulk medium of quarks and gluons.The scattering cross section giving a drag and diffusion corresponding to a space transport coefficient D s (p → 0) in agreement to lQCD [1,7,53].In PbPb such an approach is able to provide a good prediction for R AA (p t ),v 2 (p t ) and v 3 (p t ) [7,20,27,28,53].In Fig. 1 are shown the final charm quark transverse momentum distribution obtained after the evolution in QGP, in different collision systems.The shown spectra are normalized to the same number of charm quarks in order to compare them as a function of p T .In particular we compare these spectra with the initial FONLL distribution, the fireball lifetime decreases with the system size, therefore in smaller systems charm quarks have a thermalization time that is larger w.r.t. the fireball lifetime, suggesting that the high p T distribution remains close to the non-equilibrium pQCD initial distribution.In fact, in Fig. 1, moving from larger systems, like PbPb (blue solid line), to small systems, i.e.OO (gold solid line), we observe a flattening of the charm quark spectra.Furthermore, we have explored also the extreme scenario where the charm quark are fully thermalized as in Eq.9 (black solid line in Fig. 1); from one hand this case allow us to have a more direct comparison with SHM model, and from the other hand it give us the opportunity to explore the sensitivity to the charm distribution function on the multi-charmed hadrons production.

IV. PRODUCTION IN PB-PB COLLISIONS
In this section, we discuss the results for the total yields of charmed and multi-charmed hadrons using the model described in the previous section in the case of PbPb collisions at √ s = 5.02 TeV.The presence of the resonance has a relevant impact because it supply a substantial contribution in addition to the ground state production.In this work we include multiple states for the different species, the complete set of considered states is listed in Table IV for single charmed hadrons and in Table V for multi-charmed hadrons.The single-charmed resonances considered are the ones present and confirmed, at the moment, by the Particle Data Group [75].Following the same approach used in [41] we consider a statistical factor given by [m TABLE IV: Ground states of charmed mesons and baryons as well as their first excited states including their decay modes with their corresponding branching ratios as given in Particle Data Group [76,77].
TABLE V: Ground states and first excited states for multicharmed baryons.
shown the yields for single charmed hadrons in P bP b collisions for 0−10% centrality at mid-rapidity obtained with our model of coalescence plus fragmentation (redyellow square points).Some parameters (Temperature and Volume) in our model have been set in order to compare the results with the particle production predicted by the Statistical Hadronization Model (blue open circle points) [60].In Fig. 2 are also shown the results for SHM considering an enhanced set of charmed baryons (green open triangle points) with respect to the ones listed by the PDG, as suggested in studies on charm hadron production with statistical models [66,78].This result seems in accord with former indication for the production of charmed baryons in pp and P bP b collisions; where our coalescence plus fragmentation model predicts a baryon enhancement with respect to SHM.Our results are in very good agreement with recent measurements in P bP b collisions at 5.02 TeV with a Λ c production dN Λc /dy = 3.28±0.42(stat)±0.44(syst)±0.16(BR)[43].
In pp collisions the results of our coalescence plus fragmentation model are similar to the SHM ones, within the uncertainties, in the case of D 0 ,D + ,D * ,D s , but exhibit a significant difference for Λ c and Ξ c considering statistical models with the same number of resonances [50,78]; this behaviour seems to be confirmed also in PbPb collisions.In Fig. 3 are shown the results for all charmed hadrons considered, adding with respect to Fig. 2 the J/ψ, the Ω c and the multi-charmed baryons Ξ ++ cc , Ω cc and Ω ccc .The Ω c production in our model is about one order of magnitude smaller with respect to the Λ c yield, a results that is in line with what has been obtained in our previous work in pp collision [50], that turns out to be larger than the statistical yield of about a factor ten.The multicharmed baryons that we consider, i.e.Ξ ++ cc ,Ω cc and Ω ccc , have yields of about 3 and 5 order of magnitude smaller than the total charm quarks available in the system formed during the collision, i.e. the dN/dy for Ξ ++ cc is about (8± 2) •10 −3 , for Ω cc is about (0.5± 0.02) •10 −3 and for Ω ccc is between (0.12−1.01)•10 −4 , see also Table VI.For multi-charmed baryons we have found a quite large sensitivity to the underlying charm p T distribution function.If we assume a fully thermalized charm distribution we obtain an enhancement of the yields w.r.t. the realistic distribution shown in Fig. 1, this can be expected because when a thermal distribution is considered there is a large presence of charm quarks concentrated in a small momentum region at very low p T w.r.t. the realistic distributions.This feature, facilitates the recombination mechanism, because of the larger probability to find charm quarks close in phase space in the region where is present the peak of the charm distribution function.For the multi-charmed hadrons, this property results in an enhancement in the final total yields, that is more sensitive with respect to the single-charmed hadrons because of their quark content.In fact the range shown for this particles in the plot, corresponds to the yields obtained with realistic distribution for the lower limit and the one obtained with thermal distribution for the upper limit.In Fig. 2 and Fig. 3 we compare our results for all the charmed hadrons with the production obtained with SHM (blue open circles) [60].From the comparison we observe that our model gives an enhancement for the single-charmed baryons, i.e.Λ c , Ξ c and Ω c of about a factor ∼ 2−3.However the baryons production becomes similar when in the SHM an enhanced set of baryon resonances is considered, as shown in figures by the green open triangles.For the multi-charmed hadrons production shown in Fig. 3, especially for Ω ccc SHM predicts a larger production w.r.t. the coalescence results with realistic distribution (lower limit of the band), but a similar production is obtained when the thermal distribution is used also in the coalescence model.In Fig. 3 we have reported also the J/ψ production in both dynamical and thermal charm cases, with a resulting yields that is in the range (1−6)•10 −2 , about a factor 2 smaller than the SHM result and experimental measurements [79,80].However, with our modeling we are considering that the production for the J/ψ happens at the freeze-out hypersurface with temperature T c , while in dynamical dissociation and recombination approaches consider the survival probability of J/Ψ above T c [81,82] that we have not included in our modeling.Considering a recombination model based on a dynamical approach goes quite beyond the scope of this paper that is focused in providing first predictions and system size dependence of charmed baryons.

A. Sensitivity to hadron size
In this section, we explore how the microscopic details of our hadronization model have a role on the Ω ccc production.In fact while for D, D s , Λ c and Ξ c we employ the results from the quarks model, for Ω ccc the mean square radius and in general the wave function is quite unknown.As discussed in the previous section the Wigner function depends only on the widths, that are directly related to the mean charge radius.In order to obtain a quantitative estimate of the effect produced by the variation of the Wigner function widths we have performed a calculation in PbPb varying the widths in such a way to have a resulting variation in the square root of the charged squared radius in a range between -50% and 50% of the established value discussed in Section II.It is also interesting to compare the effects on the Ω ccc production with those for D 0 and Λ c .In Fig. 4 are shown the variation of the particle yield for D 0 (black open diamonds with line), Λ c (blue triangle with line) and Ω ccc (orange circles with line) with respect to the yield obtained with the previous widths, used in this way as a baseline reference.This variation is shown as a function of the variation of ⟨r 2 ch ⟩ 1/2 .We notice that the D 0 (and so D + , D * ,D s ) and Λ c yields are almost unchanged in the two extreme cases, this behaviour can be explained recalling that this two species provide the majority of the charm hadron production; an enlargement or a shrinking of the spatial widths would normally lead to a larger or a smaller production respectively, but the imposed charm quark conservation in conjunction with the constrained coalescence probability at zero momentum engender a compensation of the size change effect.The just mentioned concurrence of constraints has not a big effect on the Ω ccc production; in particular because, in this case, the particle production is five orders of magnitude smaller than the two aforementioned particles.As a consequence, the production of this multi-charmed particle have an impact that is negligible on the condition about the charm quark conservation.In this way, the outcome of changing the hadron radius, in Ω ccc case, is a larger variation of the total production w.r.t.D 0 and Λ c ; it can be quantitative described by an increase of about 60% when the radius is increased of the 50%, and a reduction of about 80% when the radius is decreased by 50%.Furthermore we have found for the Ξ ++ cc a very similar behaviour like the one shown for the Ω ccc , we don't show it in Fig. 4.  II .Finally assuming that the Ω ccc production via fragmentation is marginal, due to the very large mass of the baryons, this result suggests that the observations about this multi-charm particle production is very sensitive to the microscopic characteristics of the hadronization process and in particular to the wave function.This is particularly interesting because for charm quark it could be employed a potential model to compute the Ω ccc wave function using the heavy quark free energy from the lQCD.This would be similar to the J/Ψ, but Ω ccc production, having 3 charms, should be even more sensitive to the wave function.A seminal work in this direction can be found in [61].Moreover the production of this multicharmed hadrons can be a clearer probe with respect to other charmed hadrons production, because it is partially disentangled from the effects of charm number conservation constraint.

V. PRODUCTION EVOLUTION WITH COLLISION SYSTEM:PB-PB, KR-KR, AR-AR, O-O
In this section, we want to examine the effect on the production caused by the specifics of parton distributions in different collision systems.We have employed the model set for PbPb collision system in some other collision system, in particular in KrKr , ArAr and OO .As discussed in Sec.III the fireball parameters are summarised in in Table III and the charm dN c /dy scales with the T AA thickness function.As in the case of PbPb collisions, we start from FONLL p T distribution that evolves in a QGP medium described by a relativistic Boltzmann approach.The differences in the final p T distribution are shown in Fig. 1.As expected, in smaller systems, the final transverse momentum charm quark distributions are flatter.The dN/dy obtained by coalescence plus fragmentation for each species, as a function of A 1/3 , are shown in Fig. 5 in comparison with SHM results and summarized in Table VI.As shown, our productions scales accordingly to SHM for the single charmed particle, albeit an overall larger production of Λ c (blue circles with line) as already shown in Fig. 2.
For Ω ccc (orange circles with line) and Ξ ++ cc (green circles with line), we evaluate the yields scaling the Wigner function widths assuming the same oscillator frequency of Ξ c and Ω c as discussed in Section II (see Eq.7).We observe that our model predictions for the particle yields, have a similar but still different scaling with the decrease of the system size in comparison with SHM.Notice that for the absolute yields we show a band where the upper limit comes from the assumption of fully thermalized charm quarks while the lower limit comes from the realistic simulation of partial thermalization evaluated through the Boltzmann transport approach.In the realistic case, the Ξ ++ cc absolute production shown in Figure 5 is close to SHM values (full square) for all the systems.Instead for the Ω ccc absolute production we observe a reduction for the realistic case in all systems compared to SHM results; for the fully thermalized case there is an agreement in PbPb , as discussed before and showed in Fig. 3, and an enhancement w.r.t.SHM in smaller systems.Moreover, going back to the Ω ccc yield, its absolute value depends also on the assumption for the wave function width (as discussed in Sec.IV, see Fig. 4).An increase of 50% of the Ω ccc charged radius ⟨r 2 ⟩ 1/2 ch can give an increment of about a 60% for the lower limit of the band in Fig. 5.As already mentioned a study of Ω ccc from lQCD potential can asses this aspect in a quite solid way.However, in Fig. 5 it is clear that the degree of charm thermalization play a major role.We discuss here the full thermal case for the sake of a more direct comparison to SHM.Upcoming data on p T distribution of D 0 , D + should put more stringent constraints on the p T distribution, but we see from Fig. 5 that Ω ccc will be quite sensitive to the degree of charm thermalization.Furthermore, as we will discuss later in this Section, the Ω ccc momentum spectra will be a strong meter of the charm thermalization.In general the charmed hadrons production depends on the charm quark number (N c ) and the system size at the freeze-out, i.e. fireball volume (V ).In order to understand how the production of multi-charmed hadrons changes with the system size, we can relate N c and V to the mass number A. In first approximation, neglecting the difference in the radial flow effect that are not so different from one system to another, the volume is proportional to the mass number V ∝ A. The charm quark number comes from the perturbative hard processes in the initial stage of the collision, so the scaling of this quantity is expected to be proportional to N coll .From Glauber model, in a central collision of two iden-tical nuclei AA, the mean number of collisions scales as N coll ∝ A 4/3 , then charm quark number should scales as N c ∝ A 4/3 .From naive considerations, the coalescence mechanism has a production that is proportional to the product of the volume and the densities of the constituent quarks involved in the hadron formation.Therefore for a charmed hadron, this scaling can be expressed as where C is the number of charm quarks contained by the hadron.Finally, considering the above relation between these quantities and the mass number, the hadron production scaling roughly results in N H ∝ A C+3 3 .Such a scaling is the one expected also in SHM if canonical suppression is discarded.We now concentrate our focus on the yield dependence as a function of A in the Ω ccc case, which is particularly sensitive to this scaling because of the presence of three charm quarks.In order to disentangle the effects coming from the distribution function and the ones coming from the system size change, we have performed calculation where we consider fireball size and parameters as in Table III, but employing a charm momentum distribution with the same p T dependence of PbPb even in KrKr , ArAr and OO .In Fig. 6 the Ω ccc productions are shown, and for clarity's sake we have scaled all the curves in such a way to obtain that different cases has the same yields when PbPb collisions are considered.In the above figure different cases are shown: the yields obtained with our model in the case of realistic distribution (red-yellow square with line), the expected scaling with A (black dashed line), the SHM production (blue circles with line), and finally what we obtain considering for KrKr , ArAr and OO the charm distribution used in PbPb in a fireball consistent with the reduced dimension, as said before (green open squares with line).The last case has been realized aiming to see an hypothetical effect of the increasingly non-equilibrium of the charm quark distribution with the decreasing of the system size.In the inset of Fig. 6 we show the cube of the normalized charm distribution in all the systems.One can easily realize that the integration in p T of this quantity provides a rough estimate of the effect to the Ω ccc yield due to the difference in momentum dependence.For example going from PbPb to ArAr there is a reduction of the production by a factor ∼ 1.7.Hence, in a coalescence approach, one should see a similar reduction of the Ω ccc yield due only to the change of the p T charm distribution function between PbPb and other collision systems.Therefore the results shown in Fig. 6 give the indication that the simple scaling ∝ A C+3 3 lacked information about the p T distribution impact on the final total production.On the other hand, this scaling is compatible with the production obtained when the charm distribution is maintained fixed, has can be seen comparing black dotted line and green solid line.It is relevant to underline that the SHM model uses a factor that consider the canonical suppression, this factor is close to the unity for single charmed hadrons in PbPb and has a decreas- FIG. 7: Ωccc spectra showed in different collision system, P bP b (blue lines), KrKr (red lines), ArAr (green lines) and OO (yellow lines).Two cases are shown, results with realistic distributions (solid lines) and thermal distributions (dashed lines).
ing behaviour going to systems with smaller A, and it becomes larger increasing the charm quark content C of the hadron [60].In our model this suppression factor is not present; the impact of taking into account this factor in the multi-charmed hadrons yield turn out to be an underestimation when the system size decreases.In Fig. 7 the Ω ccc momentum distribution are shown in different collision systems in both realistic (solid lines) and thermal (dashed lines) cases.In Sec.IV was pointed out that the shape of thermal distribution falls down quickly w.r.t. the realistic distributions, as one should expect such a behaviour is reflected also in the Ω ccc spectra.Moreover, as shown for the Ω ccc yield in PbPb collisions in Fig. 3, the production results enhanced in the thermal case.This behaviour should be expected looking at the parton distribution in Fig. 1, considering that the hadron p T is the sum of quark momenta and the fact that thermal distribution is higher than the dynamical one below to ∼ 2 GeV.Hence at p T above 6 GeV the slope of the Ω ccc spectra are very different and at high momentum the tail of realistic distributions result in a p T spectrum various order of magnitude larger w.r.t. the thermal case.
Finally, in Fig. 8 the ratios between the Ω ccc transverse momentum spectra are shown in different collision systems.Each spectrum is normalized with a factor N c (N c /V ) 2 in order to have ratios that are comparable along the different system size.The low-p T part, up to 4 GeV, of the PbPb ratios shows an increasing trend that reflects the presence of a stronger flow with respect to KrKr , ArAr and OO ; on the other hand the decreasing ratio at high p T is determined by the large nonthermalization in KrKr, ArAr and OO distributions.Notice that the PbPb /OO ratio (orange line) with the normalization N c (N c /V ) 2 has a value that is approximately 2 at low p T .This effect is related to the nonthermalization in OO , as pointed out in the discussion about Fig. 6, where the yields obtained in smaller systems with realistic distributions manifest a reduction with respect to the A (C+3)/3 scaling.In the inset of Fig. 8 we show the corresponding ratios of charm distribution normalized to N c in various systems, with the same colours as for the Ω ccc ratios.We can observe that the Ω ccc ratios reflect the charm ratios; the role of the hadronization is to stretch in momentum the ratio obtained at the charm quark level, so that the behaviour at charm very low p T (∼ 1GeV) is transposed to the Ω ccc inter-  mediate momentum region, i.e. 4-5 GeV.The ratios at charm quark level reflect the non-equilibrium effect in the evolution through the QGP as discussed for Fig. 1.This behaviour, given the particularity of the Ω ccc , is much evident and sensitive with respect to similar ratios evaluated for D 0 and Λ c .In these cases besides the effect due to the presence of light quarks there is also a strong contribution that comes from the fragmentation, that can hide the direct effect coming from the charm distribution.Hence the Ω ccc distribution might unveil direct information about the charm quark distribution with a larger sensitivity w.r.t.D 0 or Λ c .Summarizing, the consequence of the spectrum flattening going from larger to smaller systems is observed as a reduction of the ratio at higher transverse momentum; producing, on the other hand, an enhancement in the very low momentum region.This result suggests that a systematic study of the multi-charm production can provide information about the charm distribution in momentum region hardly accessible trough observables derived from other charmed hadrons, i.e.D 0 and Λ c ; providing in this way a further observable able to infer charm quark interaction in the hot QCD matter.

VI. CONCLUSIONS
In this paper, we have studied the single-charmed and multi-charmed hadron production by using an hadronization model by coalescence manly focusing on Ξ ++ cc and Ω ccc that are new heavy baryons likely to be detected in the next ALICE3 experiment at LHC.In particular, we have discussed the HF production in different collision systems, P bP b, KrKr, ArAr and OO.We have considered in our study both realistic ("dynamical") charm distribution function coming from the evolution in QGP described by Boltzmann Transport Equation and thermal distribution with flow coming from the assumption of a thermalized source for charm quarks in order to have a better comparison to SHM.For D 0 , D + , D * and Λ c we have found a scaling with A of the colliding system that is quite similar to SHM, but with absolute yield of Λ c , Ξ c and Ω c that are quite larger for the coalescence model.A result for PbPb that is in line with what found in pp collisions [50].We have found that the charmed hadron yields in PbPb show a mass ordering with an enhancement for single-charmed baryons.We have found that for D mesons and Λ c our model has compatible results with SHM in the case where the latter includes an enhanced set for baryons resonances.We predict yields for multi-charmed hadrons that is of the order of 10 −2 −10 −3 for doubled charmed hadrons, while is nearly two order of magnitude smaller in the case of Ω ccc .When we consider thermal distribution form charm quark we obtain a production that is comparable with the one from SHM.Moreover, we have studied the production of all hadrons changing the system size.Going from large to smaller collision systems we obtain roughly a decrease of the production of single-charmed and multi-charmed hadrons scaling with volume and number of charm as N c (N c /V ) C−1 .But, focusing on the Ω ccc , this study shows a breaking of this simple scaling due to the change in shape of the charm quark distribution when a realistic simulation of the QGP medium is performed; suggesting that, looking at the yields, the Ω ccc production is an observable sensitive to the non-equilibrium features of charm quarks.In particular, we have seen that the ratio of Ω ccc p T -spectra between different system size can provide a solid signature of a lack of thermalization when going to smaller system like O-O.Furthermore, in PbPb we have seen a strong sensitivity of Ω ccc yield to its mean sqaure radius, at variance with D mesons and Λ c whose yields are mainly constrained by charm conservation.

FIG. 1 :
FIG.1: Normalized charm distributions: from FONLL (black dashed line) before evolution, after evolution in P bP b (blue solid line), KrKr (red solid line), ArAr (green solid line), OO (gold solid line) collisions and the thermal production (black solid line) where m H * is the mass of the resonance.This statistical factor is given by the Boltzmann probability to populate an excited state of mass m + ∆m, at a temperature T .

TABLE III :
Fireball radii, lifetime, flow and volume considered in the different collision systems studied.

TABLE VI :
D 0 , Λc , Ξ ++ cc and Ωccc hadron dN/dy in OO, ArAr, KrKr and P bP b.Upper values corresponds to thermal charm pT distribution, lower cases to realistic ones.See text for details.