Heavy Hadrons Production by Coalescence Plus Fragmentation in AA Collisions at RHIC and LHC †

The hadronization process of heavy hadrons with bottom and charm quarks, especially for baryons Λc, in a dense QGP medium is largely not understood. We present within a coalescence plus fragmentation model the predictions for D0 and Λc spectra and the related baryon to meson ratios at RHIC and LHC. We will discuss how our model can predict values for Λc/D of the order of O(1), which is much larger than the expectations from fragmentation, and in agreement whit early data from STAR collaboration. Furthermore in the same scheme can be predicted a baryon to meson ratio Λc/D in pp collisions assuming that at the LHC top energies there can be the formation of QGP matter. The results show a considerable volume effects that significantly reduce the ratios, but still predict quite larger values with respect to fragmentation, in agreement with recent data from ALICE in pp collisions.


Introduction
Ultra-relativistic heavy ion collision at Large Hadron Collider (LHC) and at Relativistic Heavy-Ion Collider (RHIC) have been designed to reach a new state of matter composed of a strongly interacting plasma of deconfined quark and gluons the so called Quark-Gluon Plasma (QGP).The bulk properties of the matter created are governed by the light quarks and gluons while heavy quarks like charm or bottom quarks are useful probes of the QGP properties [1][2][3][4][5][6][7][8][9][10][11][12].In their final state the charm quarks appear as constituent of charmed hadrons mainly D mesons and Λ c , Σ c baryons.Recent experimental results from STAR collaboration have shown an enhancement of the baryon/meson ratio in the heavy flavor sector like the one observed for light and strange hadrons compared to the one for p-p collision [13][14][15].In particular the experimental data in 10%-60% central Au + Au collisions have shown a Λ c /D 0 ∼ 0.8 ÷ 1.5 for 3 < p T < 6 GeV which is a very large enhancement compared to the value predicted by the charm hadron fragmentation ratio for p+p collisions [16].The idea of the coalescence model comes from the fact that comoving partons in the QGP combine their transverse momentum to produce a final-state meson or baryon with higher transverse momentum [17][18][19][20].Few studies have investigated the modification of the relative abundance of the different heavy hadron species produced.In particular this can manifests in a baryon-to-meson enhancement for charmed hadrons [21,22].

Coalescence plus Fragmentation Model
The coalescence approach is based on the Wigner formalism, the momentum spectrum of hadrons formed by coalescence of quarks can be written as: where dσ i denotes an element of a space-like hypersurface, g H is the statistical factor to form a colorless hadron while f q i are the quark (anti-quark) phase-space distribution functions for i-th quark (anti-quark).f H (x 1 ...x n , p 1 ...p n ) is the Wigner function and describes the spatial and momentum distribution of quarks in a hadron and can be related to the hadron wave function.The Wigner distribution function used has a Gaussian shape in space and momentum, f M (x 1 , x 2 ; p 1 , p 2 ) = ) where x r1 and p r1 are the 4-vectors for the relative coordinates.A W is a normalization constant fixed to guarantee that in the limit p → 0 we have all the charm hadronizing.While σ r is the covariant width parameter and it can be related to the oscillator frequency ω by is the reduced mass.The width of f M can be related to the size of the hadron and in particular to the root mean square charge radius of the meson.For D + meson r 2 ch = 0.184 f m 2 corresponding to a σ p = σ −1 r = 0.283 GeV; for Λ + c the widths are fixed by the mean square charge radius of Λ + c which is given by r 2 ch = 0.15 f m 2 .We compute the coalescence probability P coal for each charm quark then we can assign a probability of fragmentation as P f rag (p T ) = 1 − P coal (p T ).Therefore the hadron momentum spectra from the charm spectrum dN f ragm /d 2 p T dy that do not undergo to coalescence is given by the convolution with the fragmentation function,for D and Λ + c we employ the Peterson fragmentation function [23] , where c is a free parameter to fix the shape of the fragmentation function and is determined assuring that the experimental data on D and Λ c production in p + p collisions are well described by a fragmentation hadronization mechanism.The value it has been fixed to c = 0.06 and c = 0.12 as discussed in [5].The relative ratios between different hadron channels are properly calculated and normalized according to the ratio of fragmentation fraction in [16].

Fireball parameters and quark distribution
We consider the systems created at RHIC in Au+Au collisions at √ s NN = 200 GeV and at LHC in Pb+Pb collisions at √ s NN = 2.76 TeV.Our approach is based on a fireball where the bulk of particles is a thermalized system of gluons and u, d, s quarks and anti-quarks.The fireball is considered at τ = 7.8 f m/c, for LHC, and τ = 4.5 f m/c, for RHIC, and the system has a temperature of T C = 165 MeV.To take into account for the collective flow, we assume a radial flow profile as β T (r T ) = R , where R is the transverse radius of the fireball.For partons at low transverse momentum, p T < 2 GeV, hence we consider a thermal distribution, instead for p T > 2.5 GeV, we consider the minijets that have undergone the jet quenching mechanism.For heavy quarks we use the transverse momentum distribution obtained by solving the relativistic Boltzmann equation [5] giving a good description of R AA and v 2 of D mesons.The heavy quark numbers are estimated to be dN c /dy 2 at RHIC and dN c /dy 15 at LHC in agreement with the energy dependence of charm production cross section [24].In the following calculation the charm quark mass used is m c = 1.3 GeV.

Results
The coalescence probability is a decreasing function with p T , and at low p T having a coalescence probability for Λ c even larger than for D 0 is a quite peculiar feature of the coalescence mechanism that we expect to lead to large values of the Λ c /D 0 ratio [22].
In Figure 1a are shown the transverse momentum spectra at midrapidity for Au + Au collisions at √ s = 200 GeV and for (0%-10%) centrality for D 0 meson (left panel) and for Λ + c baryon (right panel) [25], we can see that for D 0 the contribution of both mechanism is about similar for p T < 3 GeV and at higher p T the fragmentation becomes the dominant.For Λ + c and D 0 we have included the main hadronic channels that comes from D * 0 , D * + ,Σ * c (2520) and Σ c (2455).The coalescence mechanism is the dominant mechanism for the Λ + c production for p T < 7 GeV and it is mainly related to the fragmentation fraction from the analysis in Ref. [16], where this fraction is about the 6% of the total produced heavy hadrons.In Figure 1b we show the results for the Λ + c /D 0 ratio.Coalescence by itself predicts a rise and fall of the baryon/meson ratio, the inclusion of fragmentation reduces the ratio, and we can see a quite good agreement with the experimental data in the peak region (orange solid line) in comparison with the STAR experimental data shown by circle [15,26].In Figure 2a is shown the comparison between RHIC and LHC for the Λ + c /D 0 ratio.Coalescence predicts a similar ratio for both energies, and the same for fragmentation, because the ratio established from the experimental measured fragmentation fraction remains the same changing the collision energy.Even if the only coalescence and the only fragmentation ratio remain similar, the combined ratio is different because, for each species, the production ratio between coalescence and fragmentation is smaller at LHC than at RHIC.Therefore, at LHC the larger contribution in particle production from fragmentation [22,27] leads to a final ratio that is smaller than at RHIC.A baryon over meson ratio that is so large at low momenta, can lead also to a smaller D 0 R AA in this region.It is consequence of the charm quark number conservation and the dominance of D mesons in the total particle production, in pp collisions.
In recent years there has been a broadly discussed idea about the possible formation of QGP also in systems smaller than the one formed in heavy ion collision.We have applied our model in the case of pp collisions, assuming that a medium is formed also in this small system like the one simulated in hydrodynamics calculations [28].In Figure 2b is shown with the blue dashed line the Λ + c /D 0 ratio obtained for this kind of system.Our calculations predict the disappearance of the peak, but an enhancement at low momenta that is significantly different from the ratio obtained with the only fragmentation.Moreover, the presence of a coalescence mechanism can have a deep impact on the pp baseline used to evaluate the R AA , in particular in the case of Λ c , where the presence of coalescence implies a different behavior especially at low momenta.This point is still completely open, because of the not yet available experimental data.

Figure 2 .
Figure 2. (Color online) (a) (left) Λ + c to D 0 ratio as a function of p T and at mid-rapidity for Au + Au collisions at √ s = 200 GeV (left panel) and for Pb + Pb collisions at √ s = 2.76 TeV (right panel).(b) (right) Λ + c to D 0 ratio as a function of p T and at mid-rapidity for Pb + Pb collisions at √ s = 5.02 TeV