Initial conditions from the shadowed Glauber model for Pb+Pb at $\sqrt{s_{\rm NN}}=2.76$ TeV
Abstract
We study the initial conditions for Pb+Pb collisions at TeV using the two component Monte-Carlo Glauber model with shadowing of the nucleons in the interior by the leading ones. The model parameters are fixed by comparing to the multiplicity data of p+Pb and Pb+Pb at and TeV respectively. We then compute the centrality dependence of the eccentricities upto the fourth order as well as their event by event distributions. The inclusion of shadowing brings the Monte-Carlo Glauber model predictions in agreement with data as well as with results from other dynamical models of initial conditions based on gluon saturation at high energy nuclear collisions. Further, we find that the shadowed Glauber model provides the desired relative magnitude between the ellipticity and triangularity of the initial energy distribution required to explain the data on the even and odd flow harmonics and respectively at the LHC.
pacs:
25.75.-q,24.10.HtI Introduction
One of the most important ingredients in understanding the evolution of matter formed in heavy ion collisions (HIC) is its initial condition (IC). Currently there are several models of IC available with varying degrees of success in explaining the data Bialas et al. (1976); Eskola et al. (2000); Kharzeev and Nardi (2001); Kharzeev et al. (2005); Hirano et al. (2006); Miller et al. (2007); Broniowski et al. (2009); Schenke et al. (2012); Gale et al. (2013); Rybczynski et al. (2013); Goldschmidt et al. (2015); Niemi et al. (2015). The Monte-Carlo based IC models generate the event by event (E/E) fluctuations in observables which can be compared to those measured in experiments. Most of these models share the first step- sampling the positions of the constituent nucleons of the two colliding nuclei from their nuclear density distribution which is usually taken to be a Woods-Saxon profile Woods and Saxon (1954). In the second step, they all differ in the energy deposition scheme corresponding to a specific configuration of the nucleon positions. This finally results in different predictions of centrality dependence of multiplicity, eccentricities and their event by event distributions. The largest source of uncertainties on the extracted values of the medium properties obtained by comparing the predictions of the theoretical models to data is known to stem from the choice of the IC Song et al. (2011).
Monte-Carlo Glauber models (MCGMs) have been reasonably successful in describing the qualitative features of various observables Miller et al. (2007); Broniowski et al. (2009). The energy deposition scheme is largely geometrical with the only dynamical input being a constant nucleon-nucleon cross-section . The recent data on correlation for the top ZDC events in U+U collisions at GeV and Au+Au interactions at GeV could not be reproduced within the ambit of the standard MCGM Adamczyk et al. (2015); Goldschmidt et al. (2015). The MCGM predictions are also in disagreement with the E/E distribution of the second flow harmonic for Pb+Pb collisions at TeV. However, dynamical models based on gluon saturation physics like IP-Glasma Gale et al. (2013) and EKRT Niemi et al. (2015) are in agreement with data.
MCGMs provide a simple and intuitive description of the IC and hence there have been considerable efforts to address the above issues within the geometric approach of the MCGM Rybczynski et al. (2013); Moreland et al. (2015). Recently, we have shown that the inclusion of shadowing effect due to the leading nucleons on those located deep inside provides a simple and physical picture that brings the predictions of the shadowed MCGM (shMCGM) in agreement with that of data as well as dynamical models like IP-Glasma at the top RHIC energy for Au+Au as well as U+U collisions Chatterjee et al. (2015). In this paper we present the results of the shMCGM for Pb+Pb collisions at TeV and compare with IP-Glasma model predictions Gale et al. (2013) as well as the LHC data Aamodt et al. (2011); Timmins (2013); Aad et al. (2013, 2014).
In the next section II, we provide the details of our shMCGM and the values of the parameters of the model estimated by comparing with data. In section III we present the results obtained in the shMCGM. We provide estimates for the centrality dependence of various eccentricities and their E/E distributions. We find good quantitative agreement with data as well as IP-Glasma results. Finally, in section IV we summarise.
Ii The Model
The details of the shMCGM are given in Ref. Chatterjee et al. (2015). Here we summarise its main features. The shMCGM is an extension of the two component MCGM. In the latter, the energy deposited at on the plane transverse to the beam axis (which is along the z axis) is assumed to be a linear superposition of two terms- and where counts the number of participant nucleons and is the number of possible binary collisions between them. The total charged multiplicity is also assumed to have a similar linear relation with the total and
(1) | |||||
(2) |
where the s are sampled from a negative binomial distribution () with variance Bozek and Broniowski (2013)
(3) |
and are the overall normalization parameters for the energy deposited and multiplicity produced respectively. is usually called the hardness factor which is fixed by comparing with data. The criterion for a binary collision between nucleon from nucleus at and nucleon from nucleus at is given by where is the squared distance in the transverse plane between the two nucleons
(4) |
All nucleons that suffer at least a single binary collision is treated as a participant. Thus in the standard two component MCGM approach all the participants are treated democratically irrespective of their positions along -axis. In the shMCGM we introduce the effect of shadowing due to leading participants on the others through the following ansatz
(5) |
where is the shadowing effect on a participant due to other nucleons from the same nucleus which are in front and shadow it. Thus all the participants are no more treated on equal footing- the leading nucleons contribute to energy deposition more than those located deep inside. Thus overall we have the following four parameters in the shMCGM- , , and which are constrained by data.
Iii Results
We will now present the results of the shMCGM. In Fig. 1 (a) we have plotted the probability distribution of charged particle multiplicity observed in p+Pb collision at 5.02 TeV for Chatrchyan et al. (). We have compared the data with results from both MCGM as well as shMCGM. The nucleons are sampled from a Woods-Saxon distribution of the Pb nuclear density. The Woods-Saxon parameters of the Pb nucleus are: the radius fm obtained from the parametrisation and the surface diffusion fm Kolb and Heinz (2003). We find good description of the long tail in the distribution only after including NBD fluctuation Bozek and Broniowski (2013) with . In Fig. 1 (b) we display the centrality dependence of the charged particle multiplicity in Pb+Pb collisions at TeV and contrast the results with the available data Aamodt et al. (2011). It is clear that this plot does not discriminate between MCGM and shMCGM as both the models describe the data well for suitably adjusted values of the model parameters. Thus the shadow parameter can not be fixed from this plot. We have fixed from the E/E distribution plot of scaled assuming linear hydrodynamic response, i.e. . The values of the parameters so obtained is tabulated in Table 1. Apart from the parameters listed in Table 1, we need the nucleon-nucleon cross-section which is taken here as mb (the corresponding value at RHIC is mb). We take fm, which is used in the Gaussian ansatz below to smear the energy deposited by a participant located at ,
(6) |
Model | ||||
---|---|---|---|---|
MCGM | 4.05 | 0.11 | 1 | - |
shMCGM | 4.05 | 0.32 | 1 | 0.08 |
We notice that in both MCGM as well as shMCGM, nearly doubles from GeV to TeV. On the other hand, drops by as we go from RHIC to LHC in MCGM while it stays almost constant in shMCGM. However, the shadow parameter drops by at LHC compared to highest RHIC energy. This suggests that increases with decrease in and hence the distinction between the two Glauber approaches will become even more significant at lower energies (e.g. GSI-FAIR energies).
Thus, having fixed the parameters of the Glauber models, we now look at the model predictions for other observables and compare them with data as well as results from other dynamical models of IC like IP-Glasma. We first study the centrality dependence of the mean eccentricity harmonics of the initial energy deposited in the overlap region
(7) |
where . The coordinate system is chosen such that . represents averaging over the transverse plane with the initial energy deposited on the transverse plane as the weight function. In Fig. 2, we have plotted the ensemble average of , and vs their centrality which is determined from the final state (FS) charged hadron multiplicity. We have shown the results for MCGM, shMCGM and IP-Glasma. All the models show a rising trend for with centrality which is also expected from geometrical arguments- events with lower multiplicity occur with larger impact parameter which results in larger ellipticity of the overlap region. For all centralities, in shMCGM is enhanced as compared to MCGM and is also in good agreement with IP-Glasma results. The higher value of in shMCGM as compared to MCGM was also found at RHIC energy and is a typical effect due to nucleon shadowing Chatterjee et al. (2015). The ends of the minor axis of the collision zone have on an average larger number of participants as compared to the ends of the major axis of the overlap region which calls for larger shadowing effect at the ends of the minor axis. This effectively reduces the minor axis more than the major axis, thus increasing the of the overlap region.
The higher harmonics arise due to granularity in the IC which is controlled by : smaller corresponds to larger granularity Bhalerao et al. (2011); Schenke et al. (2014). This explains the rising trend in and with centrality. We note that both MCGM and shMCGM yield almost similar values for these higher eccentricities even though in shMCGM the effective number of participants is smaller than in MCGM due to the shadowing effect. This is so because the shadowing effect systematically weakens only those sources which have other sources in front. Thus, weakening of these sources in the bulk do not increase the granularity in the transverse plane.
The correlation between the even-odd harmonics largely stem from the correlation of the initial state (IS). Starting from the observed correlation in the data of at the LHC, an allowed band for the ratio of r.m.s values of to was obtained in Ref. Retinskaya et al. (2014) within the realm of linear response. In Fig. 3 we have shown this band. We also show the values obtained for the same quantity in MCGM and shMCGM. The enhancement of in shMCGM as compared to MCGM as noted earlier in Fig. 2 also helps here- it pushes the prediction for the ratio of r.m.s. of to into the band that is favored by data unlike the case of MCGM which underpredicts as compared to the band.
We now turn our attention from the mean geometric properties in the IC to their fluctuations. We first analyse the E/E distributions of the scaled by their ensemble average values and compare with that of IP-Glasma Gale et al. (2013) as well as ATLAS data of Timmins (2013); Aad et al. (2013). As long as the hydrodynamic response is linear ( where is a constant), we expect the E/E distributions of to be a good representative of . In Figs. 4, 5 and 6 we have plotted the E/E distribution plots for , and for the following centrality classes: , and . Overall, there is good quantitative agreement between shMCGM and data as well as IP-Glasma. It is well known that the standard MCGM produces a broader E/E distribution as compared to data as well as IP-Glasma results Schenke et al. (2014); Niemi et al. (2015). However as already argued earlier in Ref. Chatterjee et al. (2015), the shadowing effect shadows the participants as well as their E/E fluctuations in position which eventually results in narrower E/E distribution that is in good agreement with data and IP-Glasma predictions.
For further peripheral centralities, the agreement in case of the E/E distributions of worsen. On the other hand, the distributions of the higher harmonics continue to be in good agreement. Recently, this has been shown as evidence of cubic hydrodynamic response to ellipticity in peripheral collisions Noronha-Hostler et al. (2015). A modified predictor of , in terms of was shown to do a much better job in the EKRT model for such peripheral bins Niemi et al. (2015) where it is defined as:
(8) |
In Fig. 7 we have plotted the E/E distribution for and compared with the case of as well as from data for the centrality bins of and . We indeed find better agreement between data of and shMCGM prediction for than .
The standard deviation is a good measure of fluctuation about the mean value. ATLAS data is available on the scaled variance of the flow coefficients which can be compared to that of initial eccentricity assuming linear response Aad et al. (2013). These quantities have been recently shown to agree well with MCGM computations Rybczynski and Broniowski (2015). In Fig. 8, we have plotted the centrality dependence of in MCGM and shMCGM. We find almost equally good agreement between data and the different Glauber approaches. Two particle cumulants and four particle cumulants are other quantities of interest that are often measured in experiments and compared to theory,
(9) | |||||
(10) |
From these cumulants it is possible to define the following scaled moments which within linear response will have one to one correspondence between the initial state(IS) and the final state (FS)
(11) |
We have plotted the centrality dependence of in Fig. 9 and find good qualitative agreement between data Aad et al. (2014), MCGM and shMCGM. Thus it is clear from Figs. 3, 4, 8 and 9 that the relative magnitudes of and E/E distribution plots of scaled can discriminate clearly between MCGM and shMCGM. The rest of the observables have weaker discriminatory power.
Iv Summary and Discussions
Glauber models provide a simple and intuitive picture of the initial conditions for the system produced in heavy ion collisions. This model provides the centrality dependence of various ensemble average observables like charged particle multiplicity, anisotropies of the initial energy deposited etc. The MCGM generates event by event distributions of these quantities. The above observables including their event by event distribution can be measured in experiments and compared to such models of initial condition. Overall, Monte-Carlo Glauber models manage to provide a good qualitative description of the data. However, such geometrical models can not describe the recent data on correlation from U+U at GeV and event by event flow data for Pb+Pb at TeV. This has called for a lot of ongoing effort to address such issues within the geometric approach of the Glauber model. We have earlier successfully addressed these issues at the top RHIC energy by introducing the effect of shadowing due to leading nucleons on the other participants located in the bulk Chatterjee et al. (2015). Here we have extended our study to Pb+Pb collisions at TeV. This idea of eclipse of nucleons in presence of other nucleons inside a nucleus has been proposed about sixty years back Glauber (1955). However, the significance of this effect in the phenomenology of heavy ion collisions has been hardly explored before. The current as well as our earlier work Chatterjee et al. (2015) suggest that the nucleon shadowing plays an important role in the ICs of HICs.
We find that for all centralities is enhanced while and do not change much due to the inclusion of shadowing. This brings the predictions of the shadowed Glauber model for the mean values as well as event by event distributions of the eccentricities in agreement with IP-Glasma results. The relative magnitudes of and event by event distributions of clearly demonstrate the superior performance of the shadowed Glauber model as compared to its conventional version, MCGM. The shadow parameter drops by as we go from top RHIC to LHC energies. This implies that at lower energies, the effect of shadowing in Glauber models will be even more important. Currently for low energies where the fireball is expected to carry a non-zero net baryon number, dynamical models are yet to be formulated and predictions made. On the other hand, within the ambit of the shadowed Monte-Carlo Glauber model, it is straightforward to make predictions at all energies as long as there is a good understanding of the dependence of the model parameters.
Acknowledgement: We would like to thank Prithwish Tribedy for providing the IP-Glasma data. SC acknowledges him for many fruitful discussions on the initial condition and thanks “Centre for Nuclear Theory” [PIC XII-RD-VEC-5.02.0500], Variable Energy Cyclotron Centre for support. SG acknowledges Department of Atomic Energy, Govt. of India for support.
References
- A. Bialas, M. Bleszynski, and W. Czyz, Nucl. Phys. B111, 461 (1976).
- K. J. Eskola, K. Kajantie, P. V. Ruuskanen, and K. Tuominen, Nucl. Phys. B570, 379 (2000).
- D. Kharzeev and M. Nardi, Phys. Lett. B507, 121 (2001).
- D. Kharzeev, E. Levin, and M. Nardi, Phys. Rev. C71, 054903 (2005).
- T. Hirano, U. W. Heinz, D. Kharzeev, R. Lacey, and Y. Nara, Phys. Lett. B636, 299 (2006).
- M. L. Miller, K. Reygers, S. J. Sanders, and P. Steinberg, Ann. Rev. Nucl. Part. Sci. 57, 205 (2007).
- W. Broniowski, M. Rybczynski, and P. Bozek, Comput. Phys. Commun. 180, 69 (2009).
- B. Schenke, P. Tribedy, and R. Venugopalan, Phys. Rev. Lett. 108, 252301 (2012).
- C. Gale, S. Jeon, B. Schenke, P. Tribedy, and R. Venugopalan, Phys. Rev. Lett. 110, 012302 (2013).
- M. Rybczynski, W. Broniowski, and G. Stefanek, Phys. Rev. C87, 044908 (2013).
- A. Goldschmidt, Z. Qiu, C. Shen, and U. Heinz, Phys. Rev. C92, 044903 (2015).
- H. Niemi, K. J. Eskola, and R. Paatelainen, (2015), arXiv:1505.02677 [hep-ph] .
- R. D. Woods and D. S. Saxon, Phys. Rev. 95, 577 (1954).
- H. Song, S. A. Bass, U. Heinz, T. Hirano, and C. Shen, Phys. Rev. Lett. 106, 192301 (2011), [Erratum: Phys. Rev. Lett.109,139904(2012)].
- L. Adamczyk et al. (STAR), Phys. Rev. Lett. 115, 222301 (2015).
- J. S. Moreland, J. E. Bernhard, and S. A. Bass, Phys. Rev. C92, 011901 (2015).
- S. Chatterjee, S. K. Singh, S. Ghosh, M. Hasanujjaman, J. Alam, and S. Sarkar, (2015), arXiv:1510.01311 [nucl-th] .
- K. Aamodt et al. (ALICE), Phys. Rev. Lett. 106, 032301 (2011).
- A. R. Timmins (ALICE), J. Phys. Conf. Ser. 446, 012031 (2013).
- G. Aad et al. (ATLAS), JHEP 11, 183 (2013).
- G. Aad et al. (ATLAS), Eur. Phys. J. C74, 3157 (2014).
- P. Bozek and W. Broniowski, Phys. Rev. C88, 014903 (2013).
- S. Chatrchyan et al. (CMS), PhysicsResultsHIN12015 .
- P. F. Kolb and U. W. Heinz, (2003), arXiv:nucl-th/0305084 [nucl-th] .
- R. S. Bhalerao, M. Luzum, and J.-Y. Ollitrault, Phys. Rev. C84, 054901 (2011).
- B. Schenke, P. Tribedy, and R. Venugopalan, Phys. Rev. C89, 064908 (2014).
- E. Retinskaya, M. Luzum, and J.-Y. Ollitrault, Phys. Rev. C89, 014902 (2014).
- J. Noronha-Hostler, L. Yan, F. G. Gardim, and J.-Y. Ollitrault, (2015), arXiv:1511.03896 [nucl-th] .
- M. Rybczynski and W. Broniowski, (2015), arXiv:1510.08242 [nucl-th] .
- R. J. Glauber, Phys. Rev. 100, 242 (1955).